Wenjun Kou

A fully resolved active musculo-mechanical model for esophageal transport

Esophageal transport is a physiological process that mechanically transports an ingested food bolus from the pharynx to the stomach via the esophagus, a multilayered muscular tube. This process involves interactions between the bolus, the esophagus, and the neurally coordinated activation of the esophageal muscles. In this work, we use an immersed boundary (IB) approach to simulate peristaltic transport in the esophagus. The bolus is treated as a viscous fluid that is actively transported by the muscular esophagus, and the esophagus is modeled as an actively contracting, fiber-reinforced tube. Before considering the full model of the esophagus, however, we first consider a standard benchmark problem of flow past a cylinder. Next a simplified version of our model is verified by comparison to an analytic solution to the tube dilation problem. Finally, three different complex models of the multi-layered esophagus, which differ in their activation patterns and the layouts of the mucosal layers, are extensively tested. To our knowledge, these simulations are the first of their kind to incorporate the bolus, the multi-layered esophagus tube, and muscle activation into an integrated model. Consistent with experimental observations, our simulations capture the pressure peak generated by the muscle activation pulse that travels along the bolus tail. These fully resolved simulations provide new insights into roles of the mucosal layers during bolus transport. In addition, the information on pressure and the kinematics of the esophageal wall resulting from the coordination of muscle activation is provided, which may help relate clinical data from manometry and ultrasound images to the underlying esophageal motor function.

Simulation studies of circular muscle contraction, longitudinal muscle shortening, and their coordination in esophageal transport

On the basis of a fully coupled active musculomechanical model for esophageal transport, we aimed to find the roles of circular muscle (CM) contraction and longitudinal muscle (LM) shortening in esophageal transport, and the influence of their coordination. Two groups of studies were conducted using a computational model. In the first group, bolus transport with only CM contraction, only LM shortening, or both was simulated. Overall features and detailed information on pressure and the cross-sectional area (CSA) of mucosal and the two muscle layers were analyzed. In the second group, bolus transport with varying delay in CM contraction or LM shortening was simulated. The effect of delay on esophageal transport was studied. For cases showing abnormal transport, pressure and CSA were further analyzed. CM contraction by itself was sufficient to transport bolus, but LM shortening by itself was not. CM contraction decreased the CSA and the radius of the muscle layer locally, but LM shortening increased the CSA. Synchronized CM contraction and LM shortening led to overlapping of muscle CSA and pressure peaks. Advancing LM shortening adversely influenced bolus transport, whereas lagging LM shortening was irrelevant to bolus transport. In conclusion, CM contraction generates high squeezing pressure, which plays a primary role in esophageal transport. LM shortening increases muscle CSA, which helps to strengthen CM contraction. Advancing LM shortening decreases esophageal distensibility in the bolus region. Lagging LM shortening no longer helps esophageal transport. Synchronized CM contraction and LM shortening seems to be most effective for esophageal transport.

A continuum mechanics-based musculo-mechanical model for esophageal transport

In this work, we extend our previous esophageal transport model using an immersed boundary (IB) method with discrete fiber-based structural model, to one using a continuum mechanics-based model that is approximated based on finite elements (IB-FE). To deal with the leakage of flow when the Lagrangian mesh becomes coarser than the fluid mesh, we employ adaptive interaction quadrature points to deal with Lagrangian-Eulerian interaction equations based on a previous work (Griffith and Luo [1]). In particular, we introduce a new anisotropic adaptive interaction quadrature rule. The new rule permits us to vary the interaction quadrature points not only at each time-step and element but also at different orientations per element. This helps to avoid the leakage issue without sacrificing the computational efficiency and accuracy in dealing with the interaction equations. For the material model, we extend our previous fiber-based model to a continuum-based model. We present formulations for general fiber-reinforced material models in the IB-FE framework. The new material model can handle non-linear elasticity and fiber-matrix interactions, and thus permits us to consider more realistic material behavior of biological tissues. To validate our method, we first study a case in which a three-dimensional short tube is dilated. Results on the pressure-displacement relationship and the stress distribution matches very well with those obtained from the implicit FE method. We remark that in our IB-FE case, the three-dimensional tube undergoes a very large deformation and the Lagrangian mesh-size becomes about 6 times of Eulerian mesh-size in the circumferential orientation. To validate the performance of the method in handling fiber-matrix material models, we perform a second study on dilating a long fiber-reinforced tube. Errors are small when we compare numerical solutions with analytical solutions. The technique is then applied to the problem of esophageal transport. We use two fiber-reinforced models for the esophageal tissue: a bi-linear model and an exponential model. We present three cases on esophageal transport that differ in the material model and the muscle fiber architecture. The overall transport features are consistent with those observed from the previous model. We remark that the continuum-based model can handle more realistic and complicated material behavior. This is demonstrated in our third case where a spatially varying fiber architecture is included based on experimental study. We find that this unique muscle fiber architecture could generate a so-called pressure transition zone, which is a luminal pressure pattern that is of clinical interest. This suggests an important role of muscle fiber architecture in esophageal transport.

Inter-rater agreement of novel high-resolution impedance manometry metrics: Bolus flow time and esophageal impedance integral ratio

Novel high‐resolution impedance manometry (HRIM) metrics of bolus flow time (BFT) and esophageal impedance integral (EII) ratio have demonstrated clinical utility, though the reliability of their analysis has not been assessed. We aimed to evaluate the inter‐rater agreement of the BFT and EII ratio.

Simulation studies of the role of esophageal mucosa in bolus transport

Based on a fully coupled computational model for esophageal transport, we analyzed the role of the mucosa (including the submucosa) in esophageal bolus transport and how bolus transport is affected by mucosal stiffness. Two groups of studies were conducted using a computational model. In the first group, a base case that represents normal esophageal transport and two hypothetical cases were simulated: (1) esophageal mucosa replaced by muscle and (2) esophagus without mucosa. For the base case, the geometric configuration of the esophageal wall was examined and the mechanical role of mucosa was analyzed. For the hypothetical cases, the pressure field and transport features were examined. In the second group of studies, cases with mucosa of varying stiffness were simulated. Overall transport characteristics were examined, and both pressure and geometry were analyzed. Results show that a compliant mucosa helped accommodate the incoming bolus and lubricate the moving bolus. Bolus transport was marginally achieved without mucosa or with mucosa replaced by muscle. A stiff mucosa greatly impaired bolus transport due to the lowered esophageal distensibility and increased luminal pressure. We conclude that mucosa is essential for normal esophageal transport function. Mechanically stiffened mucosa reduces the distensibility of the esophagus by obstructing luminal opening and bolus transport. Mucosal stiffening may be relevant in diseases characterized by reduced esophageal distensibility, elevated intrabolus pressure, and/or hypertensive muscle contraction such as eosinophilic esophagitis and jackhammer esophagus.

Normal Values of Esophageal Distensibility and Distension-Induced Contractility Measured by Functional Luminal Imaging Probe Panometry

Background & Aims: Functional luminal imaging probe (FLIP) panometry provides a comprehensive evaluation of esophageal functional at the time of endoscopy, including assessment of esophageal distensibility and distension‐induced esophageal contractility. However, the few and inconsistent findings from healthy individuals pose challenges to the application of FLIP to research and clinical practice. We performed FLIP panometry in asymptomatic volunteers. Methods: We performed a prospective study of 20 asymptomatic volunteers (ages, 23–44; 14 women) who were evaluated with 16‐cm FLIP positioned across the esophagogastric junction (EGJ) and distal esophagus (and in 8 subjects also repositioned at the proximal esophagus) during sedated upper endoscopy. FLIP data were analyzed with a customized program that generated FLIP panometry plots and calculated the EGJ‐distensibility index (DI) and distensibility plateaus (DP) of distal and proximal esophageal body. Distension‐induced esophageal contractility was also assessed. Results: The median EGJ‐DI was 5.8 mm2/mm Hg (interquartile range [IQR], 4.9–6.7 mm2/mm Hg); all 20 subjects had an EGJ‐DI greater than 2.8 mm2/mm Hg. The median DP values from all subjects tested were 20.2 mm (IQR, 19.8–20.8 mm) at the distal body, 21.1 mm (IQR, 20.3–22.9 mm) at the proximal body, and greater than 18 mm at both locations. Repetitive antegrade contractions (RACs) were observed in all 20 subjects; in 19 of 20 (95%) subjects, the RAC pattern persisted for 10 or more consecutive antegrade contractions. Conclusions: Normal parameters of FLIP panometry are EGJ‐DI greater than 2.8 mm2/mm Hg, DP greater than 18 mm, and antegrade contractions that occur in a repetitive pattern (RACs)—these can be used as normal findings for esophageal distensibility and distension‐induced contractility. These values can be used in comparative studies of esophageal diseases, such as achalasia and eosinophilic esophagitis, and will facilitate application of FLIP panometry to clinical practice.

Studies of abnormalities of the lower esophageal sphincter during esophageal emptying based on a fully coupled bolus–esophageal–gastric model

The aim of this work was to develop a fully coupled bolus–esophageal–gastric model based on the immersed boundary–finite element method to study the process of esophageal emptying across the esophagogastric junction (EGJ). The model included an esophageal segment, an ellipsoid-shaped stomach, a bolus, and a simple model of the passive and active sphincteric functions of the lower esophageal sphincter (LES). We conducted three sets of case studies: (1) the effect of a non-relaxing LES; (2) the influence of the tissue anisotropy in the form of asymmetrical right- and left-sided compliance of the LES segment; and (3) the influence of LES and gastric wall stiffness on bulge formation of the distal esophageal wall. We found that a non-relaxing LES caused sustained high wall stress along the LES segment and obstruction of bolus emptying. From the simulations of tissue anisotropy, we found that the weaker side (i.e., more compliant) of the LES segment sustained greater deformation, greater wall shear stress, and a greater high-pressure load during bolus transit. In the third set of studies, we found that a right-sided bulge in the esophageal wall tends to develop during esophageal emptying when LES stiffness was decreased or gastric wall stiffness was increased. Hence, the bulge may be partly due to the asymmetric configuration of the gastric wall with respect to the esophageal tube. Together, the observations from these simulations provide insight into the genesis of epiphrenic diverticula, a complication observed with esophageal motility disorders. Future work, with additional layers of complexity to the model, will delve into the mechanics of gastroesophageal reflux and the effects of hiatus hernia on EGJ function.

Could the peristaltic transition zone be caused by non-uniform esophageal muscle fiber architecture? A simulation study.

Based on a fully coupled computational model of esophageal transport, we analyzed how varied esophageal muscle fiber architecture and/or dual contraction waves (CWs) affect bolus transport. Specifically, we studied the luminal pressure profile in those cases to better understand possible origins of the peristaltic transition zone.