Computational Tools for the Investigation of the Male Lower Urinary Tract Functionality in Health and Disease

This paper aims to show the potentialities of computational bioengineering in the field of lower urinary tract pathophysiology. Engineering methods allow the investigation of urine flow in healthy and pathologic conditions and the analysis of urethral occlusion by means of artificial urinary sphincters. Computational models of bladder and urethra were developed and exploited to investigate the lower urinary tract physiology in health and in disease. Average male morphometric configurations were assumed, together with typical properties of both biological tissues and fluids. The reliability of the models was assessed by the mutual comparison of results and the investigation of data from experimental and clinical activities. The developed models allowed to analyze typical situations, such as the micturition in health and in disease, and the lumen occlusion by external devices. The models provided information that clinical and experimental tests barely provide, as the occurrence of turbulent phenomena within urine flow, the shear stresses at the lumen wall, the external pressure that is strictly required to occlude the lumen. The methods of bioengineering allow broadening and deepening the knowledge of the lower urinary tract functionality. More in detail, modeling techniques provide information that contributes to explain the occurrence of pathological situations, and allows to design and to optimize clinical-surgical procedures and devices.


Introduction
The comprehension of lower urinary tract physiology in health and in disease requires the investigation of urethral mechanical functionality [1][2][3][4]. The analysis pertains to the framework of biomechanics and provides considerable complexity, because it couples both solid and fluid mechanics problems [5,6]. Bladder and urethra are among the principal components of the lower urinary tract. The bladder is a hollow and distensible viscus, and its wall is a complex structure made up of mucosa, submucosa, muscolaris and serosa layers. The male urethra is a hollow tube, which is composed by different regions, such as prostatic, membranous, bulbar and penile urethra [7]. The cross-sectional shape of the urethral lumen ranges between slit-and starlike ones in the absence of urine flow, while it assumes an almost circular shape during micturition [8,9]. The urethral wall is composed of an inner lamina propria, such as a mucosa-submucosa layer, while an outer muscular layer characterizes the prostatic and the membranous urethra, and a thick spongy layer defines the outer shell of both the bulbar and the penile urethra [10]. The natural urethral sphincter is located all around the prostatic and the membranous urethra.
The analyses of urine flow and pressure along the urethra are typical tests that are performed to investigate how the bladder and the urethral sphincter are working, to discover 1 3 obstructions or to diagnose patients who have urinary incontinence or other symptoms and diseases, such as prostatic hyperplasia, reflux of urine in the prostatic ducts, prostatitis and dysuria [7,11,12]. The accurate knowledge of the pattern of urine fluid dynamics from the bladder up to the external meatus should provide a reliable tool for a better comprehension and explanation of results from uroflowmetry and urodynamic tests. Furthermore, such clinical investigations are usually invasive and painful. Computational Fluid Dynamics (CFD) [13] provides quantitative data about the urine flow regime depending on many different influencing factors, such as the specific morphometry of the lower urinary tract and obstruction phenomena [14][15][16][17]. The development of CFD models requires the acquisition of morphometric data by means of non-invasive techniques, such as echography or MRI. In detail, fluid dynamics analyses provide data about urine pressure, velocity and turbulence. The latter frequently entails inflammation of urethral mucosa, as well as stagnation and proliferation of bacteria, which are among the principal causes of dysuria [18,19].
Moreover, urethra is not a rigid tube, but soft biological tissues build up the urethral wall providing for a considerably compliant duct. Urethral wall dilatation because of urine pressure should be analyzed [20,21]. Computational Solid Mechanics (CSM) [22] provides tools for the investigation of such dilatation phenomena, together with data about stress and strain fields within urethral wall [23]. The deformation of the urethra, in turn, affects the urine flow regime, entailing the coupling of fluid flow and solid mechanics problems, which is the framework of the so-called Fluid Structure Interaction (FSI) phenomena [5,24]. Nowadays, many different applications in the cardiovascular field have proven the strength and the reliability of the FSI approach [25].
A further field of investigation for CSM and FSI pertains to urinary incontinence, which is one of the principal health concerns worldwide. Despite the Artificial Urinary Sphincter (AUS) is the gold standard for the treatment of male incontinence [26], up to 50% of patients receiving AUS requires surgical revision after initial placement [27][28][29][30]. Device malfunction, urethral atrophy, urethral erosion or infection processes mainly determine AUS failure. AUS application usually occurs all around the bulbar urethra, which is a non sphincteric region, where the local tissues are not used to support occluding actions. AUS has to ensure continence up to the highest bladder pressure conditions, and consequently it constantly applies a tiring occluding action. Such mechanical stimulation of urethral tissues is responsible for damage phenomena that may evolve to atrophy or erosion [27]. CSM and FSI allow analyzing mechanical interaction phenomena between AUS and urethra and provide quantitative information about AUS mechanical reliability. More in detail, the computational methods provide information about stress and strain within ureteral tissues, which are responsible for many different mechano-biological phenomena, such as tissue damage and vasoconstriction. In this sense, the results from computational investigations are mandatory for the optimal design of novel and more reliably AUS devices [31][32][33][34].
Previous studies reported about the computational investigation of lower urinary tract functionality, by means of a specific technique only, with particular regard to CFD methods [14][15][16][17]. The present work envisages the integration and the coupling of different computational techniques, showing the potentialities and the limitations of the different approaches. The activities aim to highlight the potentialities of computational methods of CFD, CSM and FSI for the investigation of lower urinary tract functionality in health and in disease. The different simulation scenarios are summarized in Table 1.

Materials and Methods
In the last decades, the methods of engineering burst into the field of medicine and biology and provided innovative tools for the investigation of functionality of biological systems, the reliability assessment, the planning and the design of procedures and devices. The engineering approach is usually based on the definition of a physical model of the investigated system by means of a coupled experimental and computational approach. Experimental activities are mandatory for model definition, identification and validation. With particular regard to the Actual evaluation of urethral behavior, considering the influence of both wall mechanics and urine flow High computational effort fields of solid mechanics and fluid-dynamics, experimentations must provide data about the geometrical configuration of the biological district, the mechanical behavior of the building tissues and of the involved fluids. In this sense, notes about lower urinary tract histo-morphometry and mechanical experimentations are reported in Sect. 2.1. Computational methods allow to broad experimental results to an extremely wider scenario, considering many different configurations of the system, as many different patient situations, and many different loading conditions, considering also surgical procedures and devices. Furthermore, computational models provide information that barely can be achieved by experimental activities, as the stress and the strain fields, the fluid velocity and turbulence, etc. The potentialities of different computational methods, as CFD, CSM and FSI are consequently reported in sections from 2.2 to 2.5. Stress and strain provide relevant information in the field of tissue mechanics, because they are responsible for tissue damage and failure, tissue adaptation and mechano-transduction phenomena. On the other hand, the precise prediction of the flow within biological ducts is mandatory to evaluate the effects of pathologies, or to investigate the consequences of fluid actions on the duct wall.

Experimental Characterization of the Male Lower Urinary Tract
The geometrical configuration of the male lower urinary tract was obtained by the investigation of bladder and urethra histoanatomical configurations. Average data obtained by bladder and urethra morphometry and histology [7][8][9][10] made it possible to develop 3D CAD models of the lower urinary tract (UGS NX, Siemens PLM Software, Plano, Texas, USA). The mechanical behavior of both bladder and urethra walls is usually investigated by means of experimentations performed on tissue samples from animal models, human cadavers and surgical residuals. The tissues mechanical response has been mostly investigated via uni-axial [20,35] and bi-axial [36] tensile tests. The mechanical behavior of the entire biological structure is usually analyzed by means of insufflation tests [20,37,38]. The tissues microstructural changes occurring during stretching (i.e.: collagen fibers uncramping and alignment) induce a progressive increase of material stiffness, which entails the non-linearity of the stress-strain response. This is a typical phenomenon of soft biological tissues mechanics, which is usually interpreted by means of hyperelastic or viscohyperelastic constitutive formulations [23].

CFD Analysis of Micturition
Computational Fluid Dynamics allows to simulate and to analyze the flow of urine from the bladder up to the external meatus. The bladder and the urethra inner surfaces define the mathematical domain of urine flow and, in the framework of CFD, they are assumed to be not deformable. The average three-dimensional model of the urine domain was defined by the analysis of the bladder and urethra morphometry (Fig. 1a): an ellipsoidal conformation was assumed for the bladder, while the urethra was defined by means of a circular tube with variable diameter (prostatic region: 25 mm length, 1.9 mm radius; membranous region: 15 mm length, 1.9 mm radius; bulbar region: 40 mm length, 2.0 mm radius; penile region: 110 mm length, 1.5 mm radius; fossa navicularis: 1.8 mm radius; meatus: 1.2 mm radius). The urine domain was discretized by means of hexahedral elements within the boundary layer and tetrahedral elements elsewhere. The final mesh consisted of about 2,840,000 elements, considering different mesh sizes in order to provide a reliable discretization and a reasonable computational effort. No slip condition characterized the interaction between urine and wall domain. The classical Navier-Stokes equations defined the mathematical problem, together with the k-ε model for the description of turbulent flow mean conditions [13]. The material characterization of urine assumed incompressible and Newtonian behavior. Density and viscosity at 37 °C were identified on the basis of data from Inman et al. (2013) [39], as 1.02 kg/m 3 and 0.83 cSt, respectively. The boundary conditions of the CFD problem aimed at analyzing the typical micturition, characterized by a bladder wall pressure (inlet pressure condition) and atmospheric pressure at the end of the urinary meatus (outlet pressure condition) (Fig. 1a). The bladder pressure includes the contributions from both detrusor and abdominal pressure and ranges between 0 and 12 kPa [40][41][42][43]. Transient CFD simulations were performed: the inlet bladder pressure condition was increased from 0 up to the target value in 2 s, subsequently the bladder pressure was kept constant for 5 s aiming to achieve an almost steady state flow situation.
Aiming to further exploit the CFD approach to the investigation of the lower urinary tract functionality, modifications of the urethra model geometry allowed to analyze different stricture scenarios. Urethral stricture was applied in the bulbar region to simulate typical clinical situations. Different stricture seizes were investigated, such as 40%, 55% and 70% obstruction conditions. CFD simulations have been performed by means of the general purpose code Comsol Multiphysics 5.4 (Comsol Inc., Burlington, MA, USA).

FSI Analysis of Micturition
The CFD approach to the investigation of micturition assumes rigid body behavior of both bladder and urethra walls. However, urine flow depends on the geometrical configuration of bladder and urethra domains, while the urine pressure induces dilatation phenomena of bladder and urethra walls, leading to modifications of the fluid domain. The interaction between the fluid dynamics and solid mechanics has been investigated by the methods of FSI [24]. Due to the computational effort of FSI techniques, the investigation was performed on a 30 mm length region from bulbar urethra (Fig. 1b). The urethra wall was the solid domain and it was assumed to be composed by a 0.15 mm thick inner layer of dense connective tissue, as the lamina propria, and a 5.00 mm thick outer layer of spongy tissue [20]. A fully incompressible Ogden hyperelastic formulation characterized the mechanical behavior of biological tissues [22]. The identification and validation of constitutive parameters of urethral tissues was performed on the basis of data from tensile tests performed on urethral wall samples  [20]. The inverse analysis of experimental data led to the constitutive parameters of lamina propria and spongy tissue, such as μ = 20 kPa and α = 6.178, μ = 1.1 kPa and α = 6.178, respectively [23]. A cylinder with a radius of 2 mm defined the fluid domain in the undeformed configuration. Urine properties were assumed as reported above for CFD analyses. The geometrical configuration of the system suggested to develop an axis-symmetric model. The solid domain was discretized by 4400 and 9000 rectangular elements for lamina propria and spongy layers. The discretization of the fluid domain assumed rectangular elements within the boundary layer and triangular elements elsewhere (about 30,000 elements). On the lower and the upper boundaries of the solid domain, the displacement along the axial direction was constrained. No slip condition characterized the coupling of solid and fluid domains.
The assumption allowed to approximate the continuity of bulbar urethra with the adjacent membranous and penile regions. The fluid boundary conditions aimed to interpret the typical pressure conditions that occur during micturition. Processing results from the previously described CFD analyses provided information about the pressure values along the bulbar urethra. With specific regard to a bladder pressure between 4.0 and 5.0 kPa, which is a typical condition that brings to urinate [44], the pressure within the bulbar region is about 3.0 kPa, and a variation of at most 0.1 kPa along the assumed 30 mm length region (Fig. 2a, d) should be expected. Transient FSI simulations were performed. A preliminary step (2 s) aimed to increase the urine pressure from 0 to 2.9 kPa within the entire fluid domain. During the next step, the pressure at the outlet section was kept constant, while at the inlet section the pressure was increased in 2 s up to a target value, which ranged between 2.9 and 3.0 kPa. Afterwards, both inlet and outlet conditions were kept constant for 5 s aiming to reach an almost steady state flow situation.
FSI analyses have been performed by means of the general purpose code Comsol Multiphysics 5.4 (Comsol Inc., Burlington, MA, USA).

Computational Structural Mechanics Analysis of Lumen Occlusion
The methods of Computational Solid Mechanics allow investigating a further main topic of urology, as the mechanical actions that are required to provide for lumen occlusion. The relevance of the investigation pertains to the design and the setting of AUS devices. A fully three-dimensional model of a 30 mm length portion of the bulbar urethra was developed (Fig. 1c). The morphometrical configuration and the mechanical behavior of urethral tissues were defined as previously reported in paragraph 2.3. The finite element mesh was developed by means of about 6000 and 20,000 hexahedral elements for lamina propria and spongy layer, respectively. Constraints were imposed on the bottom and the upper surfaces of the cylindrical model aiming to prevent rigid body motions. The occlusive action was defined by means of a radial pressure field that was uniformly applied over a 10 mm length central region of the outer urethra surface. The action of urine was simulated by means of an intraluminal pressure field that was applied on the lumen surface (Fig. 1c). The developed model was exploited to investigate the relationship between the occlusive action and the intraluminal pressure that is required to open the lumen. The first step of the analysis pertained to lumen occlusion. The intraluminal pressure was null, while the occlusive pressure increased from 0 to the target value. During the next step, the occlusive pressure was kept constant at the target value, while the intraluminal pressure increased from 0 up to 10 kPa (Fig. 1d). The analysis was performed assuming different target values of the occlusive pressure, which ranged between 1 and 12 kPa. The computational analysis was performed by means of the explicit dynamic solver [22] Abaqus Explicit 2018 (Dassault Systèmes, Simulia Corp., Providence, RI, USA).

FSI Analysis of Lumen Occlusion
The CSM investigation of lumen occlusion does not contemplate all the involved physical phenomena. In fact, the action of urine was simulated by means of the pressure it applies on the lumen surface. Aiming to a more reliable investigation, FSI techniques allow to account for the actual fluid behavior of urine. Again, an axis-symmetric model formulation was assumed. The FSI analysis of lumen occlusion and of the subsequent opening because of urine inflow is a real challenge. In fact, the complete occlusion of the lumen should require an unattainable distortion of the fluid domain. Therefore, the FSI investigation was performed by assuming an inverse approach, as the occluding action was imposed after the establishment of the urine flow. The steady state flow condition was provided as described in paragraph 2.3 by assuming a 3 kPa inlet pressure condition. Subsequently, an occlusive pressure field was applied over a 10 mm length central region of the urethra external surface (Fig. 1b). The occlusive pressure was progressively increased up to the target value in 2 s, and then it was kept constant up to a new steady state condition. Different analyses were performed by assuming different target values of the occlusive pressure, which values were assumed as for CSM analyses. FSI analyses have been performed by means of the general purpose code Comsol Multiphysics 5.4 (Comsol Inc., Burlington, MA, USA).

Results
Different computational techniques have been exploited aiming to analyze the functionality of the lower urinary tract, with regard to both micturition and lumen occluding actions.
The CFD investigation of urine flow during micturition was performed by assuming a fully 3D model of both bladder and urethra. The analyses reveled the actual complexity of urine flow. The non-uniformity of lumen dimension and the three-dimensional shape of the urethral duct determine the irregular pressure drop from the bladder to the external meatus (Fig. 2a, d) and the discontinuous changes of urine velocity (Fig. 2b, e). The CFD analysis further allowed to appreciate the occurrence of turbulent phenomena, as described by the turbulent kinetic energy (Fig. 2c, f), which specifies the energy dissipation because of turbulent eddies. Finally, the relationship between bladder pressure and urine flow was analyzed and compared with average data from experimental measurements [43]. Typical urodynamic test usually entails evaluation of detrusor pressure at the maximum urine flow. The ranges of such values were compared with the pressure-flow curve from CFD analysis (Fig. 2g). Bladder pressure was evaluated by considering the additive contributions of both detrusor pressure and abdominal pressure. Further CFD analyses made it possible to evaluate the influence of urethral strictures on the lower urinary tract functionality, considering both the effects on the pressure-flow relationship (Fig. 3a) and the urine velocity field (Fig. 3b).
Aiming to a more comprehensive investigation of urethral functionality during micturition, the FSI investigation of a urethral segment from the bulbar region was performed. The analyses reveled the actual deformability of the urethra during urine flow, as shown by the displacement field within the urethral wall (Fig. 4a), the relationship between urine flow and the local pressure variation (Fig. 4b) and the profile of urine velocity across the lumen section (Fig. 4c). The investigation further allowed to evaluate the influence of urine flow on the shear stresses that act on the lumen epithelium during micturition (Fig. 4d).
The methods of CSM were exploited aiming to a preliminary investigation of lumen occlusion because of external actions (Fig. 5). The simulation approach allowed to interpret the instability phenomena that the occluding action determines within the urethral wall, and the complex shapes that the lumen assumes (Fig. 5a). The subsequent application of an intraluminal pressure allowed analyzing the effects of urine pressure on the progressive opening of the lumen (Fig. 5b). Results from computational simulations were Fig. 2 CFD investigation of micturition: pressure field at a 6 kPa bladder pressure condition (a); urine velocity field (6 kPa bladder pressure) (b); turbulent kinetic energy field (6 kPa bladder pressure) (c); distribution of urine pressure (d), velocity (e) and turbulent kinetic energy (f) along the urethra at different bladder pressure conditions (ranging between 0.5 and 12 kPa). Relationship between bladder pressure and urine flow: computational results (black line) and average data from clinical measurements (pink rectangle)-bladder pressure accounts for the contributions of both detrusor pressure and abdominal pressure (g) post-processed aiming to identify the lumen status depending on both occlusive and intraluminal pressure conditions (Fig. 5c). More in details, the progressive increase of occlusive pressure determines, in stage #α (Figs. 1d, 5a), the continuous reduction of lumen radius. When the occlusive pressure reaches a critical value, the lumen collapses (lumen collapse condition). Instability phenomena implicate a sudden reduction of lumen area and the lumen shapes becomes irregular [45,46]. Subsequently, in stage #β, the lumen area continuously decreases with the occlusive pressure up the complete lumen occlusion (lumen "zero" area condition). In stage #χ, the occlusive pressure continues to increase up to the target value and, successively, the application of intraluminal pressure begins. During this stage the lumen remains occluded. When the intraluminal pressure reaches a critical value, the lumen starts to open and assumes a lenticular shape. The corresponding value of the intraluminal pressure defines the instable opening pressure. The lumen area increases with the intraluminal pressure in stage #δ. When the lumen area returns to the value of the lumen collapse condition, the lumen gets its original circular shape. The corresponding value of the intraluminal pressure defines the stable opening pressure. Subsequently, in stage #ε, the radius increases with the intraluminal pressure. The specific values of both the stable and the instable opening pressure depend on the applied occlusive pressure. Different analyses were performed to investigate such relationships (Fig. 5d). Finally, FSI techniques made it possible to investigate the comprehensive fluid-mechanical functionality of the urethra accounting for the influence of external occluding actions. Given a specific bladder pressure condition, the developed analyses provided data about the actual influence of the occlusive pressure on lumen dimension (Fig. 6a), urine flow ( Fig. 6b) and urine velocity profile (Fig. 6c). The investigation further allowed to evaluate the effects of occluding actions on flow turbulence (Fig. 6d), actual urine pressure in the occluded region (Fig. 6e) and shear stress acting on the lumen wall (Fig. 6f).  Fig. 1d), instable contraction of the lumen (stage #β), lumen complete occlusion (stage #χ) and lumen instable opening (stage #δ) (a); evolution of the lumen area throughout the time history-during stage #α (before condition A) the lumen shape is regular and its area continuously decreases with the occluding pressure up to a critical condition (lumen collapse condition), when the area suddenly falls (the derivative of lumen area reaches a minimum); during stage #β (between points A and B) the lumen shape seems indented and its area continues to decrease up to lumen closure (lumen "zero" area condition); during stage #χ (between point B and C) the lumen is occluded, the occlusive pressure reaches the target value and the intraluminal pressure increases; stage #δ (between points C and D) starts when the lumen begins to open according to a lenticular shape; stage #δ lasts until the lumen returns to a stable shape; during stage #ε (after point D) the lumen shape is regular again and its area continuously increases with intraluminal pressure (b); evaluation of the relationship between lumen area and intraluminal pressure for different target values of the occlusive pressure (c); plot of the lumen status depending on both occlusive pressure and intraluminal pressure conditionsthe lumen is stably opened within the upper-left region of the graph, while it is occluded within the bottom-right region; the two portions of the pressure map are separated by the instability region, where the lumen shape appears indented or lenticular (d)

Discussion and Conclusions
The reported models, analyses and results proved the reliability and effectiveness of computational methods for the investigation of the lower urinary tract functionality and physiology in health and in disease. CFD investigations of micturition highlight the actual complexity of urine flow throughout the urethra. The anatomical configuration of the urethra entails the nonlinear pressure drop from bladder to the urinary meatus (Fig. 2a, d) and the non-homogenous distribution of urine velocity (Figs. 2b, e). Both inlet pressure conditions and complex geometrical conformation entail turbulent phenomena, with particular regard to the bulbar region and the urinary meatus. The CFD approach further allows analyzing the effects of urethral obstructions, such in the case of urethral strictures (Fig. 3). Such information are complementary with data from urodynamic clinical investigations, providing detailed descriptions of urine flow parameters, such as the pressure-flow trend (Fig. 3a) and the urine velocity field (Fig. 3b), which can be correlated with pathologies and diseases. Computational results highlight the significant influence of urethral strictures on flow resistance, as it is reported by the clinical experience [47].
Previous CFD investigations reported about the influence of urethral strictures on the urine flow field assuming laminar flow conditions [17]. On the other side, turbulent phenomena may characterize urine flow, as previously reported by means of 2D CFD models [15]. The here proposed CFD approach accounts for both the actual 3D conformation of the lower urinary tract and turbulent phenomena. FSI investigations further prove the arising of urine turbulences within the bulbar region, because Reynolds numbers greater than 2000 characterize the local flow [13] (Fig. 4b). The investigation and the prediction of turbulences is a relevant topic, because turbulent urinary flow in the urethra may be a causal factor for different pathologies and diseases, as prostatic hyperplasia [48], reflux of urine in the prostatic ducts, prostatitis [49,50] and urinary tract infections, such as dysuria [51]. Furthermore CFD and FSI techniques allow the quantitative identification of the shear stresses that urine flow plays on the lumen wall (Figs. 4d and 6f). Shear stresses stimulate the epithelium and the local urethral-neuroendocrine cells. Such cells release mediators and, via activation of adjacent sensory nerves, stimulate urethral reflexes [52] and influence continence and sensation [53]. Furthermore, many investigations have been reported in the field of vascular fluid-dynamics about the role of shear stresses on endothelial cells physiology [54]. Shear stresses induce morphological changes of endothelium and vessel wall, but also trigger biochemical and biological events. There is considerable evidence that physiologic stresses exert vasoprotective roles and a perturbation of the stresses can lead to remodeling and possible dysfunction [55]. Similar phenomena may characterize the urethral epithelium. It further follows the relevance of investigating the influence of occlusive pressure on wall shear stresses. With regard to AUS patients, nevertheless the device is deactivated during micturition, it anyway exerts mechanical actions on the urethra, which lead to a magnification of shear stresses (Fig. 6f). With specific regard to the lower urinary tract, previous FSI investigations were performed to analyze the ureteral pressure-flow relationship, considering also the influence of occlusive actions [5]. The here proposed analyses highlight the relevance of the FSI approach for the investigation of urethra functionality considering both the ordinary micturition condition and the occlusion by AUS devices.
The reliability of the developed models is assessed by the mutual comparison of results and the evaluation of data from experimental activities and clinical tests. As an example, the CFD model prediction of the relationship between bladder pressure and urine flow well interpret results from urological tests (Figs. 2g, 3a) [47]. The urine flow predictions from both CFD and FSI analyses provide comparable results at corresponding local pressure conditions (Figs. 2d, g and 4b). The CSM predictions of occlusive pressure values that ensure lumen occlusion are comparable with clinical data from AUS surgery [26,56].
The here proposed models just provide a preliminary approach to the computational investigation of the lower urinary tract. A more refined investigation of urine flow during micturition should require a FSI model of both bladder and urethra. Furthermore, the modeling of the lumen shape should be performed accounting for its actual complex morphology. Such a model should be exploited aiming to analyze many different diseases, as urethra strictures or prostatic obstructions [17]. With regard to the investigation of lumen occlusion mechanisms by AUS devices, the modeling of the prosthetic device and of its interaction with the urethra is mandatory. This is part of further analyses that are under development at the University of Padova [32].  6 Fluid Structure Interaction analysis of urethral lumen occlusion during micturition at a 0.1 kPa gradient pressure condition: relationship between lumen radius and occlusive pressure (a); time history of urine flow at different occlusive pressure conditions-the urine flow progressively increases with the pressure gradient (0-2 s); the pressure gradient is kept constant and the urine flow reaches a steady state condition (2-7 s); the occlusive pressure is progressively applied up to the target value and the urine flow decreases (7-9 s); the occlusive pressure is kept constant and the urine flow reaches a new steady state condition (9-13 s) (b); profile of urine velocity along the lumen diameter at different occlusive pressure conditions (c); relationship between occlusive pressure energy dissipation because of turbulent phenomena (d), urine pressure in the occluded region (e), and maximum shear stress at the lumen wall (f) ◂ Funding Open Access funding provided by Università degli Studi di Padova. This study has been partially supported by University of Padova, Project FONT_BIRD2020_01 "Characterization of Artificial Urinary Sphincters for the identification of new DEVices (AUS-DEV)".

Data Availability
The models, which have been described in the paper and exploited to achieve the provided results, have been deposited within the following dataset: https ://datad ryad.org/stash /share /SqRj2 waxWH tI3bo O0moM ww6e3 5ZoAt EvrYD CagqW C4E

Compliance with Ethical Standards
Conflict of interest The author declares that he has no conflict of interest.

Ethical Approval
The author has seen and approved the final version of the manuscript. The author warrants that the article is the authors' original work, hasn't received prior publication and isn't under consideration for publication elsewhere.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creat iveco mmons .org/licen ses/by/4.0/.