Kidney tubules are lined with flow-sensing structures, yet information about the flow itself is not easily obtained. We aim to generate a multiscale biomechanical model for analyzing fluid flow and fluid-structure interactions within an elastic kidney tubule when the driving pressure is pulsatile. We developed a two-dimensional macroscopic mathematical model of a single fluid-filled tubule corresponding to a distal nephron segment and determined both flow dynamics and wall strains over a range of driving frequencies and wall compliances using finite-element analysis. The results presented here demonstrate good agreement with available analytical solutions and form a foundation for future inclusion of elastohydrodynamic coupling by neighboring tubules. Overall, we are interested in exploring the idea of dynamic pathology to better understand the progression of chronic kidney diseases such as Polycystic Kidney Disease.
Citation: Niksa Praljak, Shawn D. Ryan, Andrew Resnick. Pulsatile flow through idealized renal tubules: Fluid-structure interaction and dynamic pathologies[J]. Mathematical Biosciences and Engineering, 2020, 17(2): 1787-1807. doi: 10.3934/mbe.2020094
Related Papers:
[1]
Ray-Ming Chen .
Extracted features of national and continental daily biweekly growth rates of confirmed COVID-19 cases and deaths via Fourier analysis. Mathematical Biosciences and Engineering, 2021, 18(5): 6216-6238.
doi: 10.3934/mbe.2021311
[2]
Aili Wang, Xueying Zhang, Rong Yan, Duo Bai, Jingmin He .
Evaluating the impact of multiple factors on the control of COVID-19 epidemic: A modelling analysis using India as a case study. Mathematical Biosciences and Engineering, 2023, 20(4): 6237-6272.
doi: 10.3934/mbe.2023269
[3]
Peng An, Xiumei Li, Ping Qin, YingJian Ye, Junyan Zhang, Hongyan Guo, Peng Duan, Zhibing He, Ping Song, Mingqun Li, Jinsong Wang, Yan Hu, Guoyan Feng, Yong Lin .
Predicting model of mild and severe types of COVID-19 patients using Thymus CT radiomics model: A preliminary study. Mathematical Biosciences and Engineering, 2023, 20(4): 6612-6629.
doi: 10.3934/mbe.2023284
[4]
Sarafa A. Iyaniwura, Rabiu Musa, Jude D. Kong .
A generalized distributed delay model of COVID-19: An endemic model with immunity waning. Mathematical Biosciences and Engineering, 2023, 20(3): 5379-5412.
doi: 10.3934/mbe.2023249
[5]
Rahat Zarin, Usa Wannasingha Humphries, Amir Khan, Aeshah A. Raezah .
Computational modeling of fractional COVID-19 model by Haar wavelet collocation Methods with real data. Mathematical Biosciences and Engineering, 2023, 20(6): 11281-11312.
doi: 10.3934/mbe.2023500
[6]
Amir Yasseen Mahdi, Siti Sophiayati Yuhaniz .
Optimal feature selection using novel flamingo search algorithm for classification of COVID-19 patients from clinical text. Mathematical Biosciences and Engineering, 2023, 20(3): 5268-5297.
doi: 10.3934/mbe.2023244
[7]
Zhenwu Xiang, Qi Mao, Jintao Wang, Yi Tian, Yan Zhang, Wenfeng Wang .
Dmbg-Net: Dilated multiresidual boundary guidance network for COVID-19 infection segmentation. Mathematical Biosciences and Engineering, 2023, 20(11): 20135-20154.
doi: 10.3934/mbe.2023892
[8]
Fang Wang, Lianying Cao, Xiaoji Song .
Mathematical modeling of mutated COVID-19 transmission with quarantine, isolation and vaccination. Mathematical Biosciences and Engineering, 2022, 19(8): 8035-8056.
doi: 10.3934/mbe.2022376
[9]
Xianfang Wang, Qimeng Li, Yifeng Liu, Zhiyong Du, Ruixia Jin .
Drug repositioning of COVID-19 based on mixed graph network and ion channel. Mathematical Biosciences and Engineering, 2022, 19(4): 3269-3284.
doi: 10.3934/mbe.2022151
[10]
Javad Hassannataj Joloudari, Faezeh Azizi, Issa Nodehi, Mohammad Ali Nematollahi, Fateme Kamrannejhad, Edris Hassannatajjeloudari, Roohallah Alizadehsani, Sheikh Mohammed Shariful Islam .
Developing a Deep Neural Network model for COVID-19 diagnosis based on CT scan images. Mathematical Biosciences and Engineering, 2023, 20(9): 16236-16258.
doi: 10.3934/mbe.2023725
Abstract
Kidney tubules are lined with flow-sensing structures, yet information about the flow itself is not easily obtained. We aim to generate a multiscale biomechanical model for analyzing fluid flow and fluid-structure interactions within an elastic kidney tubule when the driving pressure is pulsatile. We developed a two-dimensional macroscopic mathematical model of a single fluid-filled tubule corresponding to a distal nephron segment and determined both flow dynamics and wall strains over a range of driving frequencies and wall compliances using finite-element analysis. The results presented here demonstrate good agreement with available analytical solutions and form a foundation for future inclusion of elastohydrodynamic coupling by neighboring tubules. Overall, we are interested in exploring the idea of dynamic pathology to better understand the progression of chronic kidney diseases such as Polycystic Kidney Disease.
1.
Introduction
In 14th December 2020, the UK has found a new variant of SARS-CoV-2 or SARS-CoV-2 VOC 202012/01 (Variant of Concern, year 2020, month 12, variant 01) [1,2,3] or B.1.1.7. There are some researches on the emergence of this SARS-CoV-2 variant [4,5]. The variant [6] has caused a huge concern across the globe [7] and has spread across other continents. Based on known epidemiologic, scientific modellings, and clinical findings, it suggests this variant has indeed increased the overall transmissibility.
On the other hand, it also indicates that there is no change in severity [8] - measured by length of hospitalization and 28-day case fatality - or occurrence of reinfection in the UK [9], despite some scientists' dire warnings [10]. Even though the evaluation of variant is still going on [11], we need to find a systematic method to delve into the fundamental features of this mutant variant. This might concern whether 202012/01 diminishes the potency of the developed vaccines [12].
To trace the fundamental features of this new variant, we sample 9 regional daily COVID-19 data in the UK from 27th March, 2020 to 8th, January, 2021 [13] as the preliminary study and sample data from 28th July 2020 to 11th May 2021 as the follow-up study. Then we create a set of batches of data vectors whose dimensions are decided by the days chosen. In our research, we choose 7 days as the datum dimension. Furthermore, to find the fundamental features of B1.1.7, we introduce the 3-valued featured vectors whose elements consist of only $ -, 0 and 1. These features are easy to apply and would serve as our fundamental database for further feature recognition. We then utilise the concept of vector rejection to measure the similarities between the set of data vectors and the features based on the total norms for the vector rejections. Afterwards, we obtain distance tensor products which record such similarities. By choosing the minimal total norms regarding the features, we retrieve the minimal total norms and their corresponding features. Since these extracted features vary, we trace their paths via cosine values for the dynamical features. This path would reveal the fundamental trend and dynamical structures of the features for B1.1.7 in the UK. The preliminary focuses on comparing the dynamical structures of COVID-19 with respect to data with London (9 regions), where it is supposed to be the onset of B1.1.7., and without London (8 regions).
Our results show there is no obvious shift for the main features before and after the variant - this confirms there is no fundamental change in the behaviour of the variant. These findings indicate B1.1.7 contain no leading features, and, in principle, shall bolster the efficacy of the vaccines.
2.
Methodology
We devise rejection distance which derives from the concept of vector rejections and three-valued featured database. Both shall serve as our fundamental methods in tracking the features of B1.1.7.
2.1. Distance tensor product
Claim 1.(rejection distance) The norm of →c, a vector rejection of →a on →b, is ||→c||=√(|→a||⋅||→b||)2−(→a⋅→b)2||→b|| or specifically the rejection distance between →a and →b
RejDist(→a,→b):=√(|→a||⋅||→b||)2−(→a⋅→b)2||→b||.
Proof. By the property that →c=→a−||→a||⋅cosθ⋅→b||→b||=→a−→a⋅→b||→b||2⋅→b and the operation of the inner product ⋅, the result follows immediately.
Rejection distance acts as a role for measuring the distances between features, including empirical features or theoretical features (3-valued features, for example). This measurement will serve our optimal feature findings.
Definition 2.1.(Distance Tensor Product, DTP) For any ordered set of vectors, whose dimension are the same, V=(→v1,→v2,⋯,→vm) and W=(→w1,→w2,⋯,→wn), we define their distance tensor product by a m-by-n matrix DTP(V,W) whose (i,j) element is RejDist(→vi,→wj).
DTP shows the norms of the vector rejections, which reveal the distances between a set of given datum vectors and the featured vectors or the distance between feature vectors.
2.2. 3-valued features
Unlike Principal Component Analysis, or PCA [14], which heavily relies on the statistical property, the 3-valued features are applied directly to trace the trend of mutant variants of COVID-19. 3-valued features consist of a set of value {−1,0,1}, which is also used in wavelet analysis [15]. In total, there are 3n−1 features (or vectors) - the 0 vector is excluded from the whole potential features, where n denotes the dimension of the datum vector.
3.
Implementation procedures
3.1. Preliminary settings
In this section, we specify the implementation processes mentioned in section 3. From the database, we select three metrics: newCasesBySpecimenDate (M1, or 1), newDeaths28DaysByDeathDate (M2, or 2), and newDeaths60DaysByDeathDate (M3, or 3) - there indexes 1, 2, 3 will be used later on; and Area Type "Region" (there are 9 regions: Yorkshire and The Humber (R1), East of England (R2), North East(R3), North West(R4), South East(R5), East Midlands(R6), South West(R7), London(R8), West Midlands(R9)). The data are sorted by date, from 27th, March, 2020 to 8th, January, 2021. The data are stored in three matrices, according to M1, M2, and M3. The we calculate their first difference (change) of the new cases, 28-day deaths, and 60-day deaths, as shown in the following (the length for each vector is 287, since we sample 288 points of time). The value in interval i is the original value in i+1 minus the one in i. When compared with the original data, these values are the second-order difference.
3.2. Data description
The data come from two main parts. One is for the preliminary study for the onset of the variant spread in the UK. The periods cover from 27th, March, 2020 to 8th, January, 2021. It lists the daily new cases, 28-day deaths and 60-day deaths of COVID-19 for 9 regions in the UK. Specifically, the metrics (or indicators) from the database are M1, M2, and M3. M2 is a variable for the new increase of people who had a positive test result for COVID-19 and died within 28 days of the first positive test, up to the date of death. M3 is a variable for the new increase of people who had a positive test result for COVID-19 and died within 60 days of the first positive test, up to the date of death.
This imported data are further separated into two sets of data - the preliminary study - one covers only 8 regions by excluding London and the other covers 9 regions which contain London. The second part is a follow-up study covering a period from 28th, July, 2020 to 11th, May, 2021 [16]. The data are not separated and are treated as a whole, i.e., the 9 regions.
3.3. Procedures
Then we set every 7 days as the dimension of one datum vector and 9 regions as the number of data vectors. In other words, there are 287/7=41 ordered batches of data sets - each set contains 9 points (or datum sub-vectors, or simply datum vectors), and each point is a 7-dimensional datum vector. We form these data with respect to cases, 28-day deaths and 60-day deaths for the 287 intervals. There are several steps in the implementation. These procedures are designed for preliminary study, but the follow-up study would also follow suit.
1. Download the data related to daily COVID-19 cases and deaths.
2. Sort the data by the 288 dates, 9 regions and 3 metrics, according to our analytical purposes. The results are presented in Table 1. Let us still call them raw data, though they are actually processed to fit our purpose.
Table 1.
Time-series sorted data for 9 regions and 3 metrics.
3. Compute the 287 first-order difference of the 288 time-series data, according to the 9 regions and 3 metrics. This will reveal the accelerated spreading trend of B.1.1.7. Let us call them difference data. Then the data are further extracted by regions and we name them difference datum vectors. The results are presented in Appendix B: difference datum vectors.
4. Split each 287 difference datum vector into 41 sub-vectors - each of which contains 7 elements (or daily data). Let us use {→Rkij} to denote the j'th 7-day sub-vector for region i, where k∈{1,2,3} (k=1 to denote "cases", k=2 to denote "28-day deaths", and k=3 to denote "60-day deaths"), 1≤i≤9 and 1≤j≤41. For example →R111=(54,21,57,34,−22,−55,51).
5. Specify the 3-valued features (7-element vectors), which is described in section 2.2, as our potential feature database. For 7-day period, there are 37−1=2186 three-valued features. The one being excluded is the zero vector. The resulting feature database goes as follows (→fi stands for i'th feature):
→f3=(1,−1,−1,−1,−1,−1,−1);→f4=(−1,0,−1,−1,−1,−1,−1); ⋯; →f2183=(1,0,1,1,1,1,1); →f2184=(−1,1,1,1,1,1,1); →f2185=(0,1,1,1,1,1,1); →f2186=(1,1,1,1,1,1,1). Let us use FDB to denote the 3-valued feature database, i.e., FDB={→f1,→f2,⋯,→f2186}.
6. Calculate the distance tensor product (DTP) which is defined in Definition 2.1, where each element is defined by
DTPki,j,h:=RejDist(→Rkij,→fh),
where 1≤i≤9;1≤j≤41;1≤h≤2186, and k∈{1,2,3}. The generated DTP could be regarded as three (regarding k) three-dimensional matrices (regarding i,j,h). The results are presented in Table 2.
8. Compute maxkj=max{DTPkj,h:1≤h≤2186},minkj=min{DTPkj,h:1≤h≤2186},argmaxkj:=argmax{DTPkj,h:1≤h≤2186},argminkj:=argmin{DTPkj,h:1≤h≤2186} for all 1≤j≤41. The results are shown in Table 4. Meanwhile, plot {minkj:1≤j≤41} and {maxkj:1≤j≤41} for k∈{1,2,3} as presented in Figure 1.
Table 4.
Min/max rejection distances and features: cases (top block), 28-day deaths (middle block) and 60-day deaths (bottom block).
Figure 1.
Min/max rejection distances between datum vectors and featured vectors for new cases (left), 28-day deaths (middle) and 60-day deaths (right).
9. Compute {cos(argminkv,argminkv+1):1≤v≤40} and {cos(argmaxkv,argmaxkv+1):1≤v≤40} for k∈{1,2,3}, i.e., the similarity indexes (cosine values) between the extracted minimal (or maximal) features to yield a dynamical trend of the representative features. The results are presented in Figures 2–4.
Figure 2.
Minimal and maximal raw and difference trends without London.
In this article, we explore the trend of underlying features of the COVID-19 new cases, 28-day deaths and 60-day deaths to reveal whether the new mutant variants of COVID-19 in the UK have caused fundamental features, which are based on 3-valued featured vectors. The results are presented in two parts: preliminary study and follow-up study.
4.1. Preliminary results
It covers the period from 27th, March, 2020 to 8th, January, 2021. The onset of the mutant variant spreads mainly in London. By comparing the dynamical trend of time-series of COVID-19 with London and without London, we shall see the role of B1.1.7 in the preliminary stage.
1. From Figure 1, we observe that the maximal and minimal rejection distance go almost together – an indication that the features chosen are qualified to capture their B.1.1.7's dynamical behaviours. Furthermore, no matter in the new cases, 28-day death, or 60-day deaths, when the time mutant variant takes place (around the tail of the three graphs), the diluting features are reverted in comparison to previous intervals. This indicates the new mutant variant changes the usual course of the collective viruses. In some sense, at this stage, it has not yet created a lead role in the pandemic, but it did affect the course of the development.
2. From the left-hand sides of Figures 2–4, the raw data for the minimal similarities almost reach 1, while the maximal similarities are pretty random. This indicates the optimal features are representative of the behaviour of B.1.1.7.
3. From the right-hand side of graphs in Figures 2 and 3, there is no obvious shift in the features when comparing the eight regions (without London) and the 9 regions (with London). This preliminary indication shows there is no clear feature from B.1.1.7 in terms of leading the pandemic. Since in the very beginning, the mutant variant is located mainly in London [18]. If there is no clear feature change between the two situations, it shows the variant is still gaining momentum and not yet to produce a leading feature. This might also provide a good indication that the vaccines for other COVID-19 viruses should be working for this mutant variant.
4.2. Follow-up results
From Figure 4, we observe that feature similarities fluctuates around 0. This indicates the leading features are pretty independent of each other. In some sense, the underlying behaviour of B.1.1.7 is versatile and less predictable. Combining the preliminary results and the follow-up ones, we shall reach a conclusion that B.1.1.7 has no clear feature and might be very easy in adopting new environments.
5.
Conclusions and future work
In this study, we devise 3-valued features to track the dynamical behaviours of SARS-CoV-2 VOC 202012/01. There features serve as the feature database for matching. By preliminary study and follow-up research, our study shows there is no clear leading features for B1.1.7. It also indicates that the virus is hard to pin down its properties and is versatile and hard to predict. The advantages for this study on the features are
a. It is a complete set features and serve as a matching database;
b. It could couple with machine learning and big data analysis to train and locate the features;
c. It could server as the pilot study to understand the potential features of the virus.
There are also some disadvantages:
a. As the length of each datum vector increases, the complexity of computation also exponentially increases;
b. It is hard to practically link the 3-valued features with real natural or scientific properties.
As for the future research, we could apply machine learning techniques on searching and selecting the optimal datum dimension. We could even lift the 3-valued features by adding more features into the candidate featured vectors. We could also associate the 3-valued features with other variables to yield a meaning interpretation of the extracted optimal features or associate them with randomness of genetic codes [17]. Another issue is whether a virus contains a leading feature would increase the severity of death rate of COVID-19 [18].
Acknowledgements
This work is supported by the Humanities and Social Science Research Planning Fund Project under the Ministry of Education of China (Grant No. 20XJAGAT001).
Conflict of interest
No potential conflict of interest was reported by the authors.
Appendix
Appendix A: Time-series raw data for 9 regions and 3 metrics
Appendix B: Difference datum vectors
(1) Difference datum vector for cases: (the length of each Rki is 287)
Appendix E: Optimal rejection distances and features
References
[1]
I. Mnassri, A. El Baroudi, Vibrational frequency analysis of finite elastic tube filled with compressible viscous fluid, Acta Mech. Solida Sin., 30 (2017), 435-444.
[2]
R. M. Terrill, An exact solution for flow in a porous pipe, Z. Angew. Math. Phys., 33 (1982), 547-552.
[3]
O. San, A. E. Staples, Dynamics of pulsatile flows through elastic microtubes, Int. J. Appl. Mech., 4 (2012).
[4]
S. Nag, A. Resnick, Biophysics and biofluid dynamics of primary cilia: Evidence for and against the flow-sensing function, Am. J. Physiol. Renal Physiol., 313 (2017), 706-720.
[5]
M. J. Lighthill, Mathematical Biofluiddynamics, Society for Industrial and Applied Mathematics, Philadelphia, 1975.
[6]
A. El Baroudi, F. Razafimahery, L. Rakotomanana, Fluid-structure interaction within three-dimensional models of an idealized arterial wall, Int. J. Eng. Sci., 84 (2014), 113-126.
[7]
M. Zamir, The Physics of Pulsatile Flow, Springer-Verlag, New York, 2000.
[8]
P. D. Cabral, J. L. Garvin, Luminal flow regulates no and o2(-) along the nephron, Am. J. Physiol. Renal Physiol., 300 (2011), 1047-1053.
[9]
L. M. Satlin, S. Sheng, C. B. Woda, T. R. Kleyman, Epithelial Na(+) channels are regulated by flow, Am. J. Physiol. Renal Physiol., 280 (2001), 1010-1018.
[10]
M. Essig, G. Friedlander, Tubular shear stress and phenotype of renal proximal tubular cells, J. Am. Soc. Nephrol, 14 (2003), S33-35.
[11]
J. B. Freund, J. G. Goetz, K. L. Hill, J. Vermot, Fluid flows and forces in development: Functions, features and biophysical principles, Development, 139 (2012), 1229-1245.
[12]
F. Kotsis, R. Nitschke, M. Doerken, G. WalzE, W, Kuehn, Flow modulates centriole movements in tubular epithelial cells, Pflugers Arch, 456 (2008), 1025-1035.
[13]
A. B. Maunsbach, G. H. Giebisch, B. A. Stanton, Effects of flow rate on proximal tubule ultrastructure, Am. J. Physiol., 253 (1987), 582-587.
[14]
H. A. Praetorius, K. R. Spring, Removal of the mdck cell primary cilium abolishes flow sensing, J. Membr. Biol., 191 (2002), 69-76.
[15]
H. A. Praetorius, K. R. Spring, The renal cell primary cilium functions as a flow sensor, Curr. Opin. Nephrol Hypertens, 12 (2003), 517-520.
[16]
D. J. Furley, J. S. Wilkie, Galen on Respiration and the Arteries: an Edition with English Translation and Commentary of De usu Respirationis, An in Arteriis Natura Sanguis Contineatur, De usu Pulsum, and De Causis Respirationis, Princeton University Press, Guildford, 1984.
[17]
W. Harvey, On the Motion of the Heart and Blood in Animals; and on the Circulation of the Blood; and on the Generation of Animals, William Benton, Chicago, 1952.
[18]
J. R. Womersley, Method for the calculation of velocity, rate of flow and viscous drag in arteries when the pressure gradient is known, J. Physiol., 127 (1955), 553-563.
[19]
J. R. Womersley, Xxiv. oscillatory motion of a viscous liquid in a thin-walled elastic tubei: The linear approximation for long waves, London, Edinburgh, Dublin Philos. Mag. J. Sci., 46 (1955), 199-221.
[20]
J. Womersley, An elastic tube theory of pulse transmission and oscillatory flow in mammalian arteries, 1957.
[21]
R. I. Macey, Pressure flow patterns in a cylinder with reabsorbing walls, Bull. Math. Biophy., 25 (1963), 1-9.
[22]
E. A. Marshall, E. A. Trowbridge, Flow of a newtonian fluid through a permeable tube - application to proximal renal tubule, Bull. Math. Bio., 36 (1974), 457-476.
[23]
C. Pozrikidis, Stokes flow through a permeable tube, Arch. Appl. Mech., 80 (2010), 323-333.
[24]
A. M. Weinstein, Nonequilibrium thermodynamic model of the rat proximal tubule epithelium, Biophys. J., 44 (1983), 153-70.
[25]
A. M. Weinstein, A mathematical model of the rat proximal tubule, Am. J. Physiol., 250 (1986), 860-873.
[26]
A. M. Weinstein, A mathematical model of rat ascending henle limb. i. cotransporter function, Am. J. Physiol. Renal Physiol., 298 (2010), 512-524.
[27]
A. M. Weinstein, A mathematical model of rat ascending henle limb. iii. tubular function, Am. J. Physiol. Renal Physiol., 298 (2010), 543-556.
[28]
A. M. Weinstein, T. A. Krahn, A mathematical model of rat ascending henle limb. ii. epithelial function, Am. J. Physiol. Renal Physiol., 298 (2010), 525-542.
[29]
A. M. Weinstein, A mathematical model of rat collecting duct. i. flow effects on transport and urinary acidification, Am. J. Physiol. Renal Physiol., 283 (2002), 1237-1251.
[30]
M. E. Downs, A. M. Nguyen, F. A. Herzog, D. A. Hoey, C. R. Jacobs, An experimental and computational analysis of primary cilia deflection under fluid flow, Comput. Methods Biomech. Biomed. Eng., 17 (2014), 2-10.
[31]
W. Liu, S. Xu, C. Woda, P. Kim, S. Weinbaum, L. M. Satlin, Effect of flow and stretch on the [Ca2+] i response of principal and intercalated cells in cortical collecting duct, Am. J. Physiol. Renal Physiol., 285 (2003), 998-1012.
[32]
A. K. O'Connor, E. B. Malarkey, N. F. Berbari, M. J. Croyle, C. J. Haycraft, P. D. Bell, et al., An inducible ciliagfp mouse model for in vivo visualization and analysis of cilia in live tissue, Cilia, 2 (2013), 8.
[33]
A. T. Layton, L. C. Moore, H. E. Layton, Signal transduction in a compliant thick ascending limb, Am. J. Physiol. Renal Physiol., 302 (2012), 1188-1202.
[34]
D. J. Marsh, O. V. Sosnovtseva, E. Mosekilde, N. H. Holstein-Rathlou, Vascular coupling induces synchronization, quasiperiodicity, and chaos in a nephron tree, Chaos, 17 (2007), 015114.
[35]
A. T. Layton, H. E. Layton, A computational model of epithelial solute and water transport along a human nephron, PLoS Comput. Biol., 15 (2019), e1006108.
[36]
COMSOL Multiphysics® v. 5.4. www.comsol.com. COMSOL AB, Stockholm, Sweden.
[37]
M. Dejam, Dispersion in non-newtonian fluid flows in a conduit with porous walls, Chem. Eng. Sci., 189 (2018), 296-310.
[38]
M. Dejam, H. Hassanzadeh, Z. X. Chen, Shear dispersion in combined pressure-driven and electro-osmotic flows in a channel with porous walls, Chem. Eng. Sci., 137 (2015), 205-215.
[39]
M. Zamir, Hemo-Dynamics, Springer International Publishing, 2015.
[40]
A. T. Layton, Modeling transport and flow regulatory mechanisms of the kidney, ISRN Biomath., 2012 (2012).
[41]
A. T. Layton, L. C. Moore, H. E. Layton, Multistable dynamics mediated by tubuloglomerular feedback in a model of coupled nephrons, Bull. Math. Bio., 71 (2009).
[42]
H. B. Atabek, Wave propagation through a viscous fluid contained in a tethered, initially stresses, orthotropic elastic tube, Biophys. J., 8 (1968), 626-649.
[43]
M. A. Day, The no-slip condition of fluid dynamics, Erkenntnis, 33 (1990), 285-296.
[44]
B. J. Cox, J. M. Hill, Flow through a circular tube with a permeable navier slip boundary, Nanoscale Res. Lett., 6 (2011), 389.
[45]
L. N. Reinking, B. Schmidt-Nielsen, Peristaltic flow of urine in the renal capillary collecting ducts of hamsters, Kidney Int., 20 (1981), 55-60.
[46]
M. Pradella, R. M. Dorizzi, F. Rigolin, Relative density of urine: Methods and clinical significance, Crit. Rev. Clin. Lab. Sci., 26 (1988), 195-242.
[47]
B. Vahidi, N. Fatouraee, A. Imanparast, A. N. Moghadam, A mathematical simulation of the ureter: Effects of the model parameters on ureteral pressure/flow relations, J. Biomech. Eng., 133 (2011), 031004.
[48]
S. Cortell, F. J. Gennari, M. Davidman, W. H. Bossert, W. B. Schwartz, A definition of proximal and distal tubular compliance. Practical and theoretical implications, J. Clin. Invest., 52 (1973), 2330-2339.
[49]
N. H. Holstein-Rathlou, D. J. Marsh, Oscillations of tubular pressure, flow, and distal chloride concentration in rats, Am. J. Physiol., 256 (1989), 1007-1014.
[50]
E. Gonzalez, P. Carpi-Medina, G. Whittembury, Cell osmotic water permeability of isolated rabbit proximal straight tubules, Am. J. Physiol., 242 (1982), 321-330.
[51]
E. Frömter, C. Mller, T. Wick, Permeability properties of the proximal tubular epithelium of the rat kidney studied with electrophysiological methods, Electrophysiol. Epithelial Cells, (1971), 119-146.
[52]
J. A. Schafer, Transepithelial osmolality differences, hydraulic conductivities, and volume absorption in the proximal tubule, Annu. Rev. Physiol., 52 (1990), 709-726.
[53]
J. Howard, Mechanics of Motor Proteins and the Cytoskeleton, Sinauer, 2001.
[54]
T. Sakai, D. A. Craig, A. S. Wexler, D. J. Marsh, Fluid waves in renal tubules, Biophys. J., 50 (1986), 805-813.
[55]
M. Mahran, A. ELsabbagh, H. Negm, A comparison between different finite elements for elastic and aero-elastic analyses, J. Adv. Res., 8 (2017), 635-648.
[56]
R. Carrisoza-Gaytan, Y. Liu, D. Flores, C. Else, H. G. Lee, G. Rhodes, et al., Effects of biomechanical forces on signaling in the cortical collecting duct (ccd), Am. J. Physiol. Renal Physiol., 307 (2014), 195-204.
[57]
J. J. Kang, I. Toma, A. Sipos, F. McCulloch, J. Peti-Peterdi, Quantitative imaging of basic functions in renal (patho) physiology, Am. J. Physiol. Renal Physiol., 291 (2006), 495-502.
[58]
A. T. Layton, L. C. Moore, H. E. Layton, Multistability in tubuloglomerular feedback and spectral complexity in spontaneously hypertensive rats, Am. J. Physiol. Renal Physiol., 291 (2006) 79-97.
[59]
E. B. Pitman, R. M. Zaritski, K. J. Kesseler, L. C. Moore, H. E. Layton, Feedback-mediated dynamics in two coupled nephrons, Bul. Math. Bio., 66 (2004), 1463-1492.
[60]
J. L. Laugesen, E. Mosekilde, N.-H. Holstein-Rathlou, Synchronization of period-doubling oscillations in vascular coupled nephrons, Chaos: Interdiscip. J. Nonlinear Sci., 21 (2011), 033128.
[61]
D. J. Marsh, O. V. Sosnovtseva, E. Mosekilde, N. H. Holstein-Rathlou, Vascular coupling induces synchronization, quasiperiodicity, and chaos in a nephron tree, Chaos: Interdiscip. J. Nonlinear Sci., 17 (2007), 015114.
This article has been cited by:
1.
Ray-Ming Chen,
The Similarities and Distances of Growth Rates Related to COVID-19 Between Different Countries Based on Spectral Analysis,
2021,
9,
2296-2565,
10.3389/fpubh.2021.695141
2.
Ray-Ming Chen,
Analysing deaths and confirmed cases of COVID-19 pandemic by analytical approaches,
2022,
231,
1951-6355,
3603,
10.1140/epjs/s11734-022-00535-4
3.
Ray-Ming Chen,
Extracted features of national and continental daily biweekly growth rates of confirmed COVID-19 cases and deaths via Fourier analysis,
2021,
18,
1551-0018,
6216,
10.3934/mbe.2021311
4.
Juan Ma, Shaojuan Ma, Xinyu Bai, Jinhua Ran,
Ergodic stationary distribution of stochastic virus mutation model with time delay,
2023,
8,
2473-6988,
21371,
10.3934/math.20231089
Figure 1. The nephron is the functional unit of the mammalian kidney, composed of a malpighian body and renal tubule. The malpighian body, located in the renal cortex, consists of a capillary tuft (glomerulus) that delivers blood from the cardiovascular system and is surrounded by a double-walled capsule, Bowman's capsule. The glomerular filter separates blood and capsular space. Primary urine flows through the Proximal Tubule (PT) into the Loop of Henle which consists of a descending limb and thick ascending limb (TAL). The ultrafiltrate then moves through the distal tubule and empties via the collecting duct (CD) which opens into the renal pelvis, ureter and bladder. Salt and water balance is achieved by the selective transfer of fluid and solutes from the tubule to the renal vein
Figure 2. Schematic plan of axisymmetric geometry with fluid-wall interplay
Figure 3. An example pressure boundary condition applied at the inlet, expressed in Eq (2.9). We set our initial condition to be Po (i.e., the steady pressure component). Here, the angular frequency ω=πrad/s and amplitude α=0.1
Figure 4. Geometry of force balance between applied fluid pressure Pf as well as shear τf, radial σrr, and axial σzz wall stresses
Figure 5. Grid geometry at the interface between the fluid and structural domains, labeled with F and S respectively for each mesh
Figure 6. Left: Axial inlet velocity profiles for the four prior mesh geometries at t=1.5s. Right: A detailed close-up at the center of tube where the fluid velocity is maximum. The peak axial velocity converges to approximately 6.34mms as we refine the mesh
Figure 7. Left: Radial percentage change for four frequencies at fixed α = 0.1 and Po = 13mmHg. Right: Analysis of the radial behavior as the amplitudes of the pulsatile flow vary over time at fixed ω = 2πrad/s. For both, Young's Modulus was fixed at E = 3×104Pa
Figure 8. Plots of time-averaged R% as a function of frequency. Resonance-type behavior is observed in the neighborhood of f = 1Hz when assigning E = 1×105Pa
Figure 9. Top and bottom velocity profiles present the corresponding time (t = 1.25, 1.5s). Peak uz and radial expansion can be seen to differ with compliance
Figure 10. Animation based on Figure 9 at time (t=1.25s) representing fluid-wall interplay. A, B, C and D present the corresponding Young's Modulus (E = 3×104, 5×104, 1×105, 5×105Pa). Total velocity and radial expansion can be seen to differ with compliance
Figure 11. Plot of the peak uz and Q at the midline of the tubular structure during time (t=1.25,1.5s) as tubular stiffening occurs. Axial velocity and volume flow rate are in units of centimeters per second and microliters per minute. Parameters α, Po, ω was fixed at 0.1, 13mmHg, and 2πradss