Dynamics of the upper urinary tract: I. Peristaltic flow through a distensible tube of limited length

Abnormal urodynamic changes in post-upper urinary tract dysfunction in ureteral obstruction rat models

Objects: This study investigated changes in upper urinary tract urodynamics (UUTU) after upper urinary tract dysfunction (UUTD). Methods: The UUTD model was induced through unilateral ureteral obstruction. To measure the renal pelvis volume, and resting pressure. Ureteral electromyography (EMG) and in situ ureteral constriction experiments were performed. Ureteral tissue was obtained for HE and masson staining, IF staining and IHC research to explore the distribution of Piezo1, and the expression of Piezo1 was studied using Western blotting. Results: The study showed that the renal pelvis volumes and the renal pelvis resting pressures gradually increased post surgery in the experimental group. The degree of ureteral tissue edema, cell necrosis and fibrosis gradually increased. The maximum contraction force and frequency of ureter in the experimental group post surgery were significantly higher than in the sham group. Western blotting showed that the expression intensity of Piezo1 gradually increased and was significantly higher than in the sham group. Further analysis of each sub-layer of the ureter revealed that Piezo1 was highly expressed in the urothelium layer, followed by the suburothelium layer, and had low expression in the smooth muscle cell layer. Conclusion: The study observed that morphological and electrophysiological changes in the upper urinary tract may be important mechanisms of abnormal UUTU. Increased expression of the Piezo1 may be a new molecular mechanism of abnormal urodynamics after UUTD.

CFD investigation of multiple peristaltic waves in a 3D unobstructed ureter

Ureters are essential components of the urinary system and play a crucial role in the transportation of urine from the kidneys to the bladder. In the current study, a three-dimensional ureter is modelled. A series of peristaltic waves are made to travel on the ureter wall to analyse and measure parameter effects such as pressure, velocity, gradient pressure, and wall shear at different time steps. The flow dynamics in the ureters are thoroughly analysed using the commercially available ANSYS-CFX software. The maximum pressure is found in the triple wave at the ureteropelvic junction and maximum velocity is observed in the single and double wave motion due to the contraction produced by the peristalsis motion. The pressure gradient is maximum at the inlet of the ureter during the single bolus motion. The contraction produces a high jet of velocity due to neck formation and also helps in urine trapping in the form of a bolus, which leads to the formation of reverse flow. Due to the reduction in area, shear stress builds on the ureter wall. The high shear stress may rupture the junctions in the ureter.

Effect of ureteral stent length and implantation position on migration after implantation

BACKGROUND

Ureteral obstruction is a urinary system disease that causes urinary retention, renal injury, renal colic, and infection. Ureteral stents are often used for conservative treatment in clinics, and their migration usually results in ureteral stent failure. The migrations include proximal migration to the kidney side and distal migration to the bladder side, but the biomechanism of stent migration is still unknown.

METHOD

Finite element models of stents with lengths from 6-30 cm were developed. The stents were implanted into the middle of the ureter to analyze the effect of stent length on its migration, and the effect of stent implantation position on 6-cm-long stent migration was also observed. The stents' maximum axial displacement was used to assess the ease of stent migration. A time-varying pressure was applied to the ureter outer wall to simulate peristalsis. The stent and ureter adopted friction contact conditions. The two ends of the ureter were fixed. The radial displacement of the ureter was used to evaluate the effect of the stent on peristalsis.

RESULTS AND DISCUSSION

The maximum migration occurs in the positive direction for a 6-cm-long stent implanted at the proximal ureter (CD and DE), but in the negative direction at the distal ureter (FG and GH). The 6-cm-long stent demonstrated almost no effect on ureteral peristalsis. The 12-cm-long stent diminished the radial displacement of the ureter from 3-5 s. The 18-cm stent diminished the radial displacement of the ureter from 0-8 s, and the radial displacement within 2-6 s was weaker than other time. The 24-cm stent diminished the radial displacement of the ureter from 0-8 s, and the radial displacement within 1-7 s was weaker than other time.

CONCLUSION

The biomechanism of stent migration and ureteral peristalsis weakening after stent implantation was explored. Shorter stents were more likely to migrate. The implantation position had less influence on ureteral peristalsis compared with the stent length, which provided a reference for stent design aimed at reducing stent migration. Stent length was the main factor affecting ureteral peristalsis. This study provides a reference for the study of ureteral peristalsis.

Fluid mechanical modeling of the upper urinary tract

The upper urinary tract (UUT) consists of kidneys and ureters, and is an integral part of the human urogenital system. Yet malfunctioning and complications of the UUT can happen at all stages of life, attributed to reasons such as congenital anomalies, urinary tract infections, urolithiasis and urothelial cancers, all of which require urological interventions and significantly compromise patients' quality of life. Therefore, many models have been developed to address the relevant scientific and clinical challenges of the UUT. Of all approaches, fluid mechanical modeling serves a pivotal role and various methods have been employed to develop physiologically meaningful models. In this article, we provide an overview on the historical evolution of fluid mechanical models of UUT that utilize theoretical, computational, and experimental approaches. Descriptions of the physiological functionality of each component are also given and the mechanical characterizations associated with the UUT are provided. As such, it is our aim to offer a brief summary of the current knowledge of the subject, and provide a comprehensive introduction for engineers, scientists, and clinicians who are interested in the field of fluid mechanical modeling of UUT. This article is categorized under: Cancer > Biomedical Engineering Infectious Diseases > Biomedical Engineering Reproductive System Diseases > Biomedical Engineering.

Pumping Patterns and Work Done During Peristalsis in Finite-Length Elastic Tubes

Balloon dilation catheters are often used to quantify the physiological state of peristaltic activity in tubular organs and comment on their ability to propel fluid which is important for healthy human function. To fully understand this system's behavior, we analyzed the effect of a solitary peristaltic wave on a fluid-filled elastic tube with closed ends. A reduced order model that predicts the resulting tube wall deformations, flow velocities, and pressure variations is presented. This simplified model is compared with detailed fluid-structure three-dimensional (3D) immersed boundary (IB) simulations of peristaltic pumping in tube walls made of hyperelastic material. The major dynamics observed in the 3D simulations were also displayed by our one-dimensional (1D) model under laminar flow conditions. Using the 1D model, several pumping regimes were investigated and presented in the form of a regime map that summarizes the system's response for a range of physiological conditions. Finally, the amount of work done during a peristaltic event in this configuration was defined and quantified. The variation of elastic energy and work done during pumping was found to have a unique signature for each regime. An extension of the 1D model is applied to enhance patient data collected by the device and find the work done for a typical esophageal peristaltic wave. This detailed characterization of the system's behavior aids in better interpreting the clinical data obtained from dilation catheters. Additionally, the pumping capacity of the esophagus can be quantified for comparative studies between disease groups.

In Situ and Intraoperative Detection of the Ureter Injury Using a Highly Sensitive Piezoresistive Sensor with a Tunable Porous Structure

Iatrogenic ureteral injury, as a commonly encountered problem in gynecologic, colorectal, and pelvic surgeries, is known to be difficult to detect in situ and in real-time. Consequently, this injury may be left untreated, thereby leading to serious complications such as infections, renal failure, or even death. Here, high-performance tubular porous pressure sensors were proposed to identify the ureter in situ intraoperatively. The electrical conductivity, mechanical compressibility, and sensor sensitivity can be tuned by changing the pore structure of porous conductive composites. A low percolation threshold of 0.33 vol % was achieved due to the segregated conductive network by pores. Pores also lead to a low effective Young's modulus and high compressibility of the composites and thus result in a high sensitivity of 448.2 kPa^-1 of sensors, which is consistent with the results of COMSOL simulation. Self-mounted on the tip of forceps, the sensors can monitor tube pressures with different frequencies and amplitudes, as demonstrated using an artificial pump system. The sensors can also differentiate ureter pulses from aorta pulses of a Bama minipig in situ and in real-time. This work provides a facile, cost-effective, and nondestructive method to identify the ureter intraoperatively, which cannot be effectively achieved by traditional methods.

A two way fully coupled fluid structure simulation of human ureter peristalsis

Numerical simulations of ureter peristalsis have been carried out in the past to understand both the flow field and ureter wall mechanics. The main objective of the current investigations is to have a better understanding of the urine transport due to the peristalsis in the ureter, thus making the information helpful for a better treatment and diagnosis of ureteral complications like urine reflux. In the current study, a numerical simulation is performed using a finite-element-based solver with a two-way fully coupled fluid structure interaction approach between the ureter wall and urine. For the first time, the ureter wall is modeled as an anisotropic hyper-elastic material based on experiments performed in previous literature on the human ureter. Peristalsis in the ureter is modeled as a series of isolated boluses. By observing the flow field it is clear that the peristalsis mechanism has a natural tendency to create a backflow as the isolated bolus moves forward. As a result, the urine can flow back from the bladder to the ureter at the ureterovesical (ureter-bladder) junctions, if the one-way valve starts to malfunction.

Fluid-structure interaction simulation of ureter with vesicoureteral reflux and primary obstructed megaureter

Two common abnormalities in ureters include primary refluxing megaureter (PRM) and primary obstructed megaureter (POM). The aim of this study was to represent the numerical simulation of the urine flow at the end of the ureter with vesicoureteral reflux (VUR) and POM during peristalsis. Methodologically, the peristalsis in the ureter wall was created using Gaussian distribution. Fluid-structure interaction (FSI) was applied to simulate urine-elastic wall interactions; and governing equations were solved using the arbitrary Lagrangian-Eulerian method. Theories such as wall elasticity, Newtonian fluid, and incompressible Navier-Stokes equations were used. Velocity fields, viscous stresses and volumetric outflow rate profiles were obtained through the simulation of the ureter with VUR and POM during peristalsis. In addition, the effect of urine viscosity on flow rate was investigated. When the bladder pressure increased, VUR occurred because of the ureterovesical junction (UVJ) dysfunction, leading to high stresses on the wall. In the POM, the outflow rate was ultimately zero, and stresses on the wall were severe in the obstructed section. Comparing the results demonstrated that the peristalsis leads to even further dilation of the prestenosis portion. It was also observed that the reflux occurs in the ureter with VUR when the bladder pressure is high. Additionally, the urine velocity during the peristalsis was higher than the non-peristaltic ureter.

A fluid-structure interaction (FSI)-based numerical investigation of peristalsis in an obstructed human ureter

Urine moves from the kidney to the bladder through the ureter. A series of compression waves facilitates this transport. Due to the highly concentrated mineral deposits in urine, stones are formed in the kidney and travel down through the urinary tract. While passing, a larger stone can get stuck and cause severe damage to ureter wall. Also, stones in the ureter obstructing the urine flow can cause pain and backflow of urine which in turn might require surgical intervention. The current study develops a 2D axisymmetric numerical model to gain an understanding of the ureter obstruction and its effects on the flow, which are critical in assessing the different treatment options. Transient computational analysis involving a two-way fully coupled fluid-structure interaction with the arbitrary Lagrangian-Eulerian method between the ureteral wall and urine flow is conducted with an obstruction in the ureter. The ureter wall is modeled as an anisotropic hyperelastic material, data of which, is based on biaxial tests on human ureter from previous literature, while the incompressible Navier-Stokes equations are solved to calculate urine flow. A finite element-based monolithic solver is used for the simulations here. The obstruction is placed in the fluid domain as a circular stone at the proximal part of the ureter. One of the objectives of this study is to quantify the effect of the ureteral obstruction. A sharp jump in pressure gradient and wall shear stress, as well as retrograde urine flow, is observed as a result of the obstruction.

A biomechanical simulation of ureteral flow during peristalsis using intraluminal morphometric data

Reflux nephropathy and vesicoureteral reflux are two of the most important abnormalities in the upper urinary system in which toxins and bacteria from the bladder infect the ureter and the kidney and initiate renal scar formation. A quantitative analysis that characterizes urine flow will further help our understanding of the ureter and also assist in the design of flow aided devices such as valves and stents to correct reflux situations. Here, A numerical simulation with fluid-structure interactions (FSI) using arbitrary Lagrangian-Eulerian (ALE) formulation and adaptive mesh procedure was introduced and solved to perform ureteral flow analysis. Incompressible Navier-Stokes equations were utilized as the governing equations of fluid domain. Ureteral in-vivo morphometric data during peristalsis were used to construct the presented model. A nonlinear material model was used to exhibit ureteral wall mechanical properties. Direct coupling method was used to solve the solid, fluid and interface equations simultaneously. Results showed that recirculation regions formed against the jet flow, neighboring the bolus peak. Through wave propagation, separation occurred behind the moving bolus on the wall and ureteropelvic reflux began from that location and extended upstream to the ureteral inlet. The maximum luminal pressure consistently occurred behind the urine bolus during peristalsis. The measured magnitude of maximum volumetric flow rate resulted from isolated bolus transportation was 0.92 ml/min. In conclusion; due to presence of fluid inertial forces during peristalsis, the function of ureteropelvic junction in prevention of reflux is significant, especially at the beginning of peristaltic wave propagation. Moreover, modeling of ureteral function using imaging data will be valuable and it may help physicians to diagnose and cure the abnormalities.

A mathematical simulation of the ureter: effects of the model parameters on ureteral pressure/flow relations

Ureteral peristaltic mechanism facilitates urine transport from the kidney to the bladder. Numerical analysis of the peristaltic flow in the ureter aims to further our understanding of the reflux phenomenon and other ureteral abnormalities. Fluid-structure interaction (FSI) plays an important role in accuracy of this approach and the arbitrary Lagrangian-Eulerian (ALE) formulation is a strong method to analyze the coupled fluid-structure interaction between the compliant wall and the surrounding fluid. This formulation, however, was not used in previous studies of peristalsis in living organisms. In the present investigation, a numerical simulation is introduced and solved through ALE formulation to perform the ureteral flow and stress analysis. The incompressible Navier-Stokes equations are used as the governing equations for the fluid, and a linear elastic model is utilized for the compliant wall. The wall stimulation is modeled by nonlinear contact analysis using a rigid contact surface since an appropriate model for simulation of ureteral peristalsis needs to contain cell-to-cell wall stimulation. In contrast to previous studies, the wall displacements are not predetermined in the presented model of this finite-length compliant tube, neither the peristalsis needs to be periodic. Moreover, the temporal changes of ureteral wall intraluminal shear stress during peristalsis are included in our study. Iterative computing of two-way coupling is used to solve the governing equations. Two phases of nonperistaltic and peristaltic transport of urine in the ureter are discussed. Results are obtained following an analysis of the effects of the ureteral wall compliance, the pressure difference between the ureteral inlet and outlet, the maximum height of the contraction wave, the contraction wave velocity, and the number of contraction waves on the ureteral outlet flow. The results indicate that the proximal part of the ureter is prone to a higher shear stress during peristalsis compared with its middle and distal parts. It is also shown that the peristalsis is more efficient as the maximum height of the contraction wave increases. Finally, it is concluded that improper function of ureteropelvic junction results in the passage of part of urine back flow even in the case of slow start-up of the peristaltic contraction wave.

Fast uncoiling kinetics of F1C pili expressed by uropathogenic Escherichia coli are revealed on a single pilus level using force-measuring optical tweezers

Uropathogenic Escherichia coli (UPEC) express various kinds of organelles, so-called pili or fimbriae, that mediate adhesion to host tissue in the urinary tract through specific receptor-adhesin interactions. The biomechanical properties of these pili have been considered important for the ability of bacteria to withstand shear forces from rinsing urine flows. Force-measuring optical tweezers have been used to characterize individual organelles of F1C type expressed by UPEC bacteria with respect to such properties. Qualitatively, the force-versus-elongation response was found to be similar to that of other types of helix-like pili expressed by UPEC, i.e., type 1, P, and S, with force-induced elongation in three regions, one of which represents the important uncoiling mechanism of the helix-like quaternary structure. Quantitatively, the steady-state uncoiling force was assessed as 26.4 ±1.4 pN, which is similar to those of other pili (which range from 21 pN for S(I) to 30 pN for type 1). The corner velocity for dynamic response (1,400 nm/s) was found to be larger than those of the other pili (400-700 nm/s for S and P pili, and 6 nm/s for type 1). The kinetics were found to be faster, with a thermal opening rate of 17 Hz, a few times higher than S and P pili, and three orders of magnitude higher than type 1. These data suggest that F1C pili are, like P and S pili, evolutionarily selected to primarily withstand the conditions expressed in the upper urinary tract.

Multipili attachment of bacteria with helixlike pili exposed to stress

A number of biomechanical properties of various types of pili expressed by Escherichia coli, predominantly their force-versus-elongation behavior, have previously been assessed in detail on a single pilus level. In vivo, however, bacteria bind in general to host cells by a multitude of pili, which presumably provides them with adhesion properties that differs from those of single pili. Based upon the previously assessed biomechanical properties of individual pili, this work presents a theoretical analysis of the adhesion properties of multipili-attaching bacteria expressing helixlike pili exposed to an external force. Expressions for the adhesion lifetime of dual- and multipili-attaching bacteria are derived and their validity is verified by Monte Carlo simulations. It is demonstrated that the adhesion lifetime of a multipili-binding bacterium depends to a large degree on the cooperativity of the attaching pili, which, in turn, depends strongly on their internal biomechanical properties, in particular their helixlike structure and its ability to elongate, which, in turn, depends on the intrinsic properties of the bonds, e.g., their lengths and activation energies. It is shown, for example, that a decrease in the length of a layer-to-layer bond in the rod of P pili, expressed by E. coli, by 50% leads to a decrease in the adhesion lifetime of a bacterium attaching by ten pili and exposed to a force of 500 pN by three orders of magnitude. The results indicate moreover that the intrinsic properties of the rod for this particular type of pili are optimized for multipili attachment under a broad range of external forces and presumably also to its in vivo environment. For example, P pili seems to be optimized to withstand a force exposure during approximately 3 s, which correspond to the time it takes for a bolus to pass a bacterium attached to the ureteral wall. Even though the results presented in this work apply quantitatively to one type of pilus, they are assumed to apply qualitatively to all helixlike pili systems expressing slip bonds.

A numerical simulation of peristaltic motion in the ureter using fluid structure interactions

An axisymmetric model with fluid-structure interactions (FSI) is introduced and solved to perform ureter flow and stress analysis. The Navier-Stokes equations are solved for the fluid and a linear elastic model for ureter is used. The finite element equations for both the structure and the fluid were solved by the Newton-Raphson iterative method. Our results indicated that shear stresses were high around the throat of moving contracted wall. The pressure gradient magnitude along the ureter wall and the symmetry line had the maximum value around the throat of moving contracted wall which decreased as the peristalsis propagates toward the bladder. The flow rate at the ureter outlet at the end of the peristaltic motion was about 650 mm3/s. During propagation of the peristalsis toward the bladder, the inlet backward flow region was limited to the areas near symmetry line but the inner ureter backward flow regions extended to the whole ureter contraction part. The backward flow was vanished after 1.5 seconds of peristalsis propagation start up and after that time the urine flow was forward in the whole ureter length, so reflux is more probable to be present at the beginning of the wall peristaltic motion.

Measurement of flow speed in the channels of novel threadlike structures on the surfaces of mammalian organs

There have been several reports on novel threadlike structures (NTSs) on the surfaces of the internal organs of rats and rabbits since their first observation by Bonghan Kim in 1963. To confirm this novel circulatory function, it is necessary to observe the flow of liquid through the NTS as well as the structurally corroborating channels in the NTS. In this article, we report on the measurement of the flow speed of Alcian blue solution in the NTSs on the organ surfaces of rabbits, and we present electron microscopic images depicting the cribrous cross-section with channels. The speed was measured as 0.3 +/- 0.1 mm/s, and the flow distance was up to 12 cm. The flow was unidirectional, and the phase contrast microscopic images showed that the NTSs were strongly stained with Alcian blue. The ultrastructure of the NTSs revealed by cryo-scanning electron microscopy and high-voltage electron microscopy showed that (1) there were cell-like bodies and globular clumps of matter inside the sinus of the channel with thin strands of segregated zones which is a microscopic evidence of the liquid flow, (2) the sinuses have wall structures surrounded with extracellular matrices of collagenous fibers, and (3) there exists a cribriform structure of sinuses. To understand the mechanism for the circulation, a quantitative analysis of the flow speed has been undertaken applying a simplified windkessel model. In this analysis, it was shown that the liquid flow through the NTSs could be due to peristaltic motion of the NTS itself.

Modeling physiology of the urinary tract

We review mathematical and physical models of physiology of the organs of the urinary tract and their functions of producing, excreting, and voiding urine. Models for urine concentration in the kidney, urine flow in the ureters, bladder filling and emptying, urethral function during micturition, pelvic floor muscles, and neural control are reviewed in the context of their application to the development of new diagnostic and therapeutic techniques. The focus of this review is on modeling of physiology and function at the tissue and organ level, as almost all research to date has been done in those areas. Although physiological models of the lower urinary tract are in their infancy, they have the long-term potential to improve our understanding of physiological mechanisms, as well as to provide environments for simulation or testing in silico of new therapies and techniques.

Progress in urodynamic research on the upper urinary tract: implications for practical urology

The development of new surgical techniques for bladder substitution and continent urinary diversion has extended interest in urodynamics of the upper urinary tract. From a subdiscipline attracting mainly scientists and bioengineers, renal pelvic kinetics and ureteral peristalsis have evolved as important factors in routine clinical urology. The observed changes in peristaltic pattern during high diuresis, obstruction and urinary reflux have influenced management of stone disease and neurogenic bladder. The demonstration that high intravesical pressure is reflected to the kidney not only when the ureteric orifice is incompetent, but also during high diuresis, established the necessity for low pressures in neobladders. Much further clarification of urinary transport from the renal tubules to the bladder should be achievable by refined techniques of fluoroscopy, isotopic renography and manometry.

Urinary tract reconstruction in children

The successful introduction of clean intermittent catheterization and increased awareness of urinary tract physiology and urodynamics have been the basis for recent major advances in urinary tract reconstructive surgery. Surgical techniques are now available to manage anatomical and functional deficiencies of any isolated or combined components of the urinary tract. The high incidence of unsatisfactory long-term results with ileal conduit diversion has led to increased popularity in urinary tract undiversion and greater utilization of reconstructive principles. As with any new surgical field of endeavour, new operative techniques are appearing at a rapid rate. In particular, there has been a recent proliferation of surgical procedures that provide a continent, low pressure, catheterizable reservoir for urine storage. Most children with major urinary tract deficiencies can now be offered socially and cosmetically unobtrusive surgical solutions without jeopardizing renal function.

Dynamics of the upper urinary tract: II. The effect of variations of peristaltic frequency and bladder pressure on pyeloureteral pressure/flow relations

The dynamics of pyeloureteral flow is described when there is no peristalsis and for peristalsis of high and intermediate frequencies, on the assumption that the ureter is uniform except in the mid-ureter and at the outlet. The possibility of upstream transmission of bladder pressure variations to the renal pelvis is considered. The overall behaviour depends on three principal variables, the maximum tube pressure in the contraction waves, the intrinsic peristaltic carrying capacity and the peristaltic frequency f, expressed in the form fT where T is the time for a peristaltic contraction wave to sweep through the ureter. At intermediate peristaltic frequencies (fT less than but comparable with one) oscillatory flow patterns can occur, in which periods of peristaltically driven flow alternate with extraperistaltic periods of flow through the open ureter. The kidney is better isolated from bladder pressure variations when the peristaltic frequency is high, but high peristaltic frequency can by itself lead to elevated renal pelvic pressure if the flow rate is high. Experimental observations in pigs are presented to support these conclusions.