Loading [MathJax]/extensions/TeX/boldsymbol.js
Research article Special Issues

On slip of a viscous fluid through proximal renal tubule with linear reabsorption

  • The hydrodynamical problem of flow in proximal renal tubule is investigated. Axisymmetric flow of viscous, incompressible fluid through the proximal renal tubule that undergoes linear reabsorption with slip at the wall is considered. The stream function is used to transform the governing equations to system of ordinary differential equations. The analytical solutions for velocity components, pressure distribution, fractional reabsorption and the shear stress are found. The effect of slip parameter and reabsorption rate on the flow have been investigated. The points of extreme values for the axial and radial velocity components are identified. The solution is applied to physiological data from human and rat kidney, and the results are presented in tables and graphs.

    Citation: Abdul M. Siddiqui, Getinet A. Gawo, Khadija Maqbool. On slip of a viscous fluid through proximal renal tubule with linear reabsorption[J]. AIMS Biophysics, 2021, 8(1): 80-102. doi: 10.3934/biophy.2021006

    Related Papers:

    [1] JinJun Yong, Changlun Ye, Xianbing Luo . A fully discrete HDG ensemble Monte Carlo algorithm for a heat equation under uncertainty. Networks and Heterogeneous Media, 2025, 20(1): 65-88. doi: 10.3934/nhm.2025005
    [2] Serge Nicaise, Cristina Pignotti . Asymptotic analysis of a simple model of fluid-structure interaction. Networks and Heterogeneous Media, 2008, 3(4): 787-813. doi: 10.3934/nhm.2008.3.787
    [3] Seung-Yeal Ha, Shi Jin, Jinwook Jung . A local sensitivity analysis for the kinetic Kuramoto equation with random inputs. Networks and Heterogeneous Media, 2019, 14(2): 317-340. doi: 10.3934/nhm.2019013
    [4] Xiangdong Du, Martin Ostoja-Starzewski . On the scaling from statistical to representative volume element in thermoelasticity of random materials. Networks and Heterogeneous Media, 2006, 1(2): 259-274. doi: 10.3934/nhm.2006.1.259
    [5] Thomas Abballe, Grégoire Allaire, Éli Laucoin, Philippe Montarnal . Application of a coupled FV/FE multiscale method to cement media. Networks and Heterogeneous Media, 2010, 5(3): 603-615. doi: 10.3934/nhm.2010.5.603
    [6] Tasnim Fatima, Ekeoma Ijioma, Toshiyuki Ogawa, Adrian Muntean . Homogenization and dimension reduction of filtration combustion in heterogeneous thin layers. Networks and Heterogeneous Media, 2014, 9(4): 709-737. doi: 10.3934/nhm.2014.9.709
    [7] Martin Heida, Benedikt Jahnel, Anh Duc Vu . Regularized homogenization on irregularly perforated domains. Networks and Heterogeneous Media, 2025, 20(1): 165-212. doi: 10.3934/nhm.2025010
    [8] Tong Yan . The numerical solutions for the nonhomogeneous Burgers' equation with the generalized Hopf-Cole transformation. Networks and Heterogeneous Media, 2023, 18(1): 359-379. doi: 10.3934/nhm.2023014
    [9] Nimisha Pandey, Pramod Kumar Mishra . Detection of DDoS attack in IoT traffic using ensemble machine learning techniques. Networks and Heterogeneous Media, 2023, 18(4): 1393-1409. doi: 10.3934/nhm.2023061
    [10] Elisabeth Logak, Isabelle Passat . An epidemic model with nonlocal diffusion on networks. Networks and Heterogeneous Media, 2016, 11(4): 693-719. doi: 10.3934/nhm.2016014
  • The hydrodynamical problem of flow in proximal renal tubule is investigated. Axisymmetric flow of viscous, incompressible fluid through the proximal renal tubule that undergoes linear reabsorption with slip at the wall is considered. The stream function is used to transform the governing equations to system of ordinary differential equations. The analytical solutions for velocity components, pressure distribution, fractional reabsorption and the shear stress are found. The effect of slip parameter and reabsorption rate on the flow have been investigated. The points of extreme values for the axial and radial velocity components are identified. The solution is applied to physiological data from human and rat kidney, and the results are presented in tables and graphs.



    Let (Ω,F,P) be a complete probability space, where Ω is a sample space, F2Ω is the σ-algebra of events, and P:F[0,1] is a probability measure. [0,T] is the temporal interval. DR2 is an open, bounded, Lipschitz physical space. D=D0D1 is the union of disjoint portions for the boundary. The unit outward normal vector to D is denoted by n. We carry out the numerical simulations of a transient heat model with random Robin coefficients and diffusion coefficients. That's to find a random function y:[0,T]×D×ΩR conforming to equations almost surely ωΩ,

    {yt[a(t,x,ω)y]=f(t,x,ω), in [0,T]×D×Ω,y(0,x,ω)=y0(x,ω), in D×Ω,ay(t,x,ω)n=0, on [0,T]×D0×Ω,ay(t,x,ω)n=α(t,x,ω)(u(t,x,ω)y(t,x,ω)), on [0,T]×D1×Ω, (1.1)

    where Robin coefficients α, diffusion coefficients a, boundary condition u, source term f, and initial condition y0 are random fields with continuous and bounded covariance functions.

    The model (1.1) describes the process of heat as well as mass transfer with random Robin factors and diffusion factors, which can be found in [3,4] and the references therein. These uncertain factors are involved in thermal environments. The temperature is quantified stochastically. In the review [3], this model can describe random high-cycle temperature fluctuations observed at the upper core structure of fast-breeder reactors or random variations in heat transfer coefficients around the stator vanes of gas turbines. The model (1.1) also appears in robust optimal boundary control problems as state equations, such as [25,26]. The effective implementation of these applications needs to solve (1.1) efficiently.

    To solve PDEs with random coefficients, there are many numerical algorithms that include polynomial chaos method, stochastic collocation method, stochastic finite element method, MC method, and so on (see, e.g., [1,9,11,21,22,27,29,31]). Where the MC method (see, e.g., [8,10,12,24]) is easy to implement and the convergence property of the MC method is not dependent on the dimension of the parameters in a random model. When the MC method is adopted, we first need independent sampling and then deal with independent numerical simulation per sample. These simulations are affected by independent diffusion factors, Robin factors, forces, initial and boundary conditions. To improve calculating efficiency, the ensemble algorithms are widely applied (see, e.g., [13,17,20,23,24]). In [23,24], the authors studied the parabolic problem with random coefficients by using the ensemble method, and obtained an error estimate. But the error estimate therein is not optimal in physical space. The authors of [18] combined the ensemble MC with the HDG method to obtain an optimal error estimate in space. For the heat conduction model with random Robin factors, we have not found other relevant results about ensemble algorithms besides our recent work [30].

    Look back at the numerical algorithm in the article [30]. ˉα(t,x) and ˉa(t,x) stand for the ensemble mean for the Robin factors as well as the diffusion factors, respectively, and we denote the fluctuations, αj:=ˉααj, aj:=ˉaaj. We adopt a uniform partition on the interval [0,T] as well as tn=nΔt,n=0,1,2,...,N. Let ynj,fnj,unj,anj, αnj, ˉan, ˉαn, anj and αnj be the values of functions yj,fj,uj,aj, αj, ˉa, ˉα, aj and αj at t=tn, respectively. yj=y(t,x;ωj) and the others are similar. In [30], based on the implicit Euler method, we developed a variational average ensemble method for the random transient heat problem (1.1). Conveniently, we called the variational average ensemble Monte Carlo (VAEMC) method: find yn+1j(j=1,2,,J;n=0,1,,N1) such that for all vH1(ˉD),

    (yn+1jynjΔt,v)+(ˉan+1yn+1j,v)+((an+1jˉan+1)ynj,v)+(ˉαn+1yn+1j,v)D1+((αn+1jˉαn+1)ynj,v)D1=(fn+1j,v)+(αn+1jun+1j,v)D1,

    After arrangement, that is

    (1Δtyn+1j,v)+(ˉan+1yn+1j,v)+(ˉαn+1yn+1j,v)D1=(fn+1j,v)+(αn+1jun+1j,v)D1+(1Δtynj,v)+((an+1j)ynj,v)+((αn+1j)ynj,v)D1. (1.2)

    Denote a=max1jJajˉa and α=max1jJαjˉα. From [30], we base on the following numerical stability condition a<1 and α<1 for the first-order ensemble time-stepping methods (1.2). In the ensemble algorithm in the literatures [17,18,19,20,23,24,30], both stability and convergence were conditionally dependent on the ratio between the fluctuating and average values of the relevant factors.

    Now, instead of employing ensemble means of Robin coefficients as well as diffusion coefficients in scheme (1.2), we propose a maximum ensemble of these coefficients at each time step. The process gets a shared coefficient matrix for the ensemble group to better improve computational efficiency. This is a new variational ensemble Monte Carlo numerical method, which we called the variational MAX ensemble Monte Carlo (VMEMC) method to conveniently distinguish it from other methods. In contrast with related methodologies, the novelty of this algorithm is that the numerical stability and convergence are unconditional.

    The following is the outline of this work. Some mathematical preliminaries and notations are listed in Section 2. The full discretization VMEMC scheme for the random transient heat model (1.1) is given in Section 3. In Section 4, the unconditional numerical stability and error estimate are analyzed. Several numerical tests are shown in Section 5. Related discussions of this algorithm and corresponding summaries are presented in Section 6.

    In this section, we will give some notations and mathematical preliminaries. For simplicity, dx,ds, and dt in some expressions will be omitted when there is no confusion. The boundaries D0 and D1 concern the experimentally accessible and inaccessible parts, respectively. Throughout this work C is a positive constant; it has different values in different places and does not rely on time step Δt, sample size J, or mesh size h.

    Let and (,) be the L2(D) norm as well as inner product, respectively. Simultaneously, D and (,)D stand for the L2(D) norm and the inner product. The Sobolev space Ws,q(D) with the norm vWs,q(D), here sN (natural number set) and q1. We denote Hs(D)=Ws,2(D). Particularly, H1(D) is equipped with the norm 1=1,D, which is defined by

    y1,D=(y2+y2)12.

    Let Hs(D) be the dual space of the bounded linear functions on Hs(D), with norm fs=sup0vHs(D)(f,v)/vs. Denote

    L(D)={v:v is a measurable function and |v|<+},

    where |v|=esssupxD|v|.

    (Ω,F,P) is a complete probability space. ZL1P(Ω) is a random variable. Denote

    E[Z]=ΩZ(ω)dP(ω).

    Let δ=(δ1,,δd) is a d-tuple with the length |δ|=di=1δi,δiN+. The stochastic Sobolev spaces ˜Ws,q(D)= LqP(Ω,Ws,q(D)) contains stochastic function v:Ω×DR, which is measurable with respect to the product σ-algebra FB(D). The norm of ˜Ws,q(D) is defined by

    v˜Ws,q(D)=(E[vqWs,q(D)])1/q=(E[|δs|D|δv|q])1/q,1q<+.

    Let ˜Hs(D)=˜Ws,2(D) L2P(Ω)Hs(D).

    We will use the tensor product Hilbert space

    X=˜L2(0,T;H1(D))L2P(0,T;H1(D);Ω)

    with its inner product

    (v,u)XE[T0D(vu+vu)].

    The induced norm is given by

    vX=(E[T0D(|v|2+v2)])1/2.

    Lemma 2.1. The norm 1,D1 is defined by

    y21,D1=D|y|2+D1|y|2.

    This norm and the standard norm 1 in H1(D) are equivalent. Particularly, there exist constants CP,CP>0 such that vH1(D)

    CPv1v1,D1,CPv1,D1v1.

    Proof. cf. [14,15].

    The weak solution of the problem (1.1) is defined as follows: a function yX is a weak solution of (1.1) if it satisfies the initial condition y(0,x,ω)=y0(x,ω)L2P(H1(D);Ω) and for T>0,

    T0(yt,v)dt+T0ΩDayvdxdP(ω)dt+T0ΩD1αyvdsdP(ω)dt=T0ΩDfvdxdP(ω)dt+T0ΩD1αuvdsdP(ω)dt (2.1)

    for all vX.

    We emphasize that the assumed regularity only requires the random fields to be square integrable. Assumptions on the finite-dimensional noise, that is, all the involved random input data depend on a finite number of real-valued random variables. Assumptions of input data for coefficients are considered:

    (A1) a=a(t,x,ω)L2P(0,T;L2(D);Ω) is uniformly bounded and coercive, i.e., amin, amax>0 such that Pro{ωΩ:a(t,x,ω)[amin,amax],(t,x)[0,T]×D}=1.

    (A2) α=α(t,x,ω)L2P(0,T;L2(D1);Ω) is uniformly bounded and coercive too, i.e., αmin, αmax>0 such that Pro{ωΩ:α(t,x,ω)[αmin,αmax],(t,x)[0,T]×D1}=1.

    Denote the ensemble group of solution variables for the equation (1.1) by yj(t,x)=y(t,x;ωj), j=1,2,...,J, with corresponding boundary conditions. [0,T]= N1n=0[tn,tn+1] is an isometric partition of the interval. These norms are utilized in the numerical error discussion: 1s<,

    We use the MC method for random sampling because it is non-intrusive, easy to implement, and its convergence is independent of the dimension of the uncertain model parameters. For , we let and , where , . Without causing confusion, we abuse the fluctuating component symbols and . Suppressing the spatial discretization, utilizing an implicit-explicit for time-direct to the model (1.1), we propose the VMEMC scheme: find such that

    (2.2)

    Using a standard finite element (FE) method to discretize the space, we obtain the following block linear system for each ensemble member :

    (2.3)

    where is the mass matrix, is the boundary mass matrix, and is the diffusion matrix.

    The above linear system is equivalent to the following:

    (2.4)

    where is the resulting coefficient matrix (independent of ). The matrix is symmetric positive definite (SPD) since both , , and are SPD. The system (2.4) can be solved with efficient block solvers [7,16,28]. Further, since only one coefficient matrix is required for computation each time step, the storage requirement is thereby reduced.

    Suppose is a quasi-uniform triangulation of the domain , such that . Let be the diameter of the element . Define . Let be the set of polynomials of degree . Denote the finite elements (FEs) space of the domain by

    the FEs space of the trace space corresponding to by here the space is defined as the restriction of on , is the number of nodes on , is the number of nodes on . The spaces above satisfy the following approximation properties [2,5] : ,

    (2.5)

    or denotes the interpolation operator related to the spaces or .

    We first present a full discretization VMEMC method with finite elements for the VMEMC scheme (2.2) of the random transient heat model (1.1). And then the VMEMC algorithm is proposed. The full-discrete VMEMC scheme: find such that

    (3.1)

    here the initial value which .

    As a comparison, here let's list the full-discrete VAEMC scheme [30] : find such that,

    (3.2)

    here the initial value which .

    By the way, the full-discrete variational Monte Carlo (VMC) scheme: find such that,

    (3.3)

    here the initial value which .

    Remark 3.1. (1) The symbols and both appear in VMEMC and VAEMC schemes. Although they all represent fluctuations in corresponding random coefficients, they all mean different things.

    (2) Although the VMEMC and VAEMC schemes are formally similar, the novelty of the VMEMC scheme is its unconditional numerical stability and convergence. And the VMEMC scheme is more efficient in the numerical simulation.

    (3) If the random thermal conductivity and convective coefficient are independent of time, such as, , , the accordingly analyzed analysis is analogous to that presented below. It (VMEMC) is unconditionally stable as well as the first-order accurate approach.

    (4) Assume that the random thermal conductivity and convective coefficient are provided with uncertain temperature-dependent values as [6], e.g., , . Let be the fully discrete approximate solution at time , , , , . Then we can replace them accordingly in the VMEMC scheme. Similarly, unconditionally stable and first-order accurate methods under assumptions that two coefficients are continuously differentiable can also be obtained.

    Applying the variational MAX ensemble scheme to the stochastic model (1.1), we first take the MC approach for random sampling. After sampling with independent identically distributed (i.i.d.), we solve the corresponding deterministic PDE for each sample. Furthermore, the solution for model (1.1) is made from the expectation of the solutions of the deterministic PDEs. A variational MAX ensemble-based Monte-Carlo algorithm (i.e., VMEMC algorithm) is proposed to quantifies uncertainty and improve its computational efficiency. The VMEMC algorithm includes the following three steps (S):

    VMEMC algorithm:

    Choose a group of random samples for the stochastic Robin coefficient, diffusion coefficient, source term, boundary, and initial conditions, , , , , and for per . Therefore, the solutions are i.i.d..

    Denote , and . For and the j-th member, we find to satisfy the VMEMC scheme. One finds the FE solution after choosing an approximate FE space.

    Take the VMEMC sample average as an approximation to the expectation . Given an object of interest , we are going to analyze these outputs for the variational ensemble systems, , , , to obtain the stochastic information.

    The unconditional stability and error estimates are presented in this section.

    The unconditional stability for the VMEMC method is below.

    Theorem 4.1. Assume , , and hypothesises and are satisfied. Then the VMEMC method is unconditionally stable. Furthermore, the numerical solution of (3.1) holds

    (4.1)

    Proof. Taking in (3.1), we obtain

    (4.2)

    Noting that

    (4.3)

    and

    (4.4)

    applying the polarization identity, integrating over the probability space, and rearranging terms in equation (4.2), we obtain

    (4.5)

    Combing Cauchy-Schwarz, Young's inequalities with the trace theorem and the Sobolev embedding theorem to the right-hand side (RHS) of (4.5), we have

    (4.6)
    (4.7)

    Under Lemma 2.1, we obtain . Multiplying by and integrating over the probability space, we can yield

    (4.8)

    For Eq (4.5), dropping two positive terms and , useing estimate Eqs (4.6) and (4.7) with , we can obtain

    (4.9)

    Summing from to , one obtains the stability result Eq (4.1).

    Remark 4.1. (1) The stability result above shows that we have controlled over the temperature approximation in both and , unconditionally. Specifically, this numerical stability is independent of both the time step size and the grid size constraints, as well as the fluctuations in the random coefficients. This feature is different from the VAEMC approach ((see, e.g., [23,30]).

    (2)Compared to the standard Backward Euler method, we see that the numerical dissipation is enhanced with the additional term

    (3) If we consider the assumption of minimum regularity, that is,

    hypotheses and are satisfied, then the VMEMC method also possesses unconditional stability.

    In this subsection, we first decompose the overall error of the VMEMC algorithm into the statistical error and the discretization error. Next, we sequentially estimate the statistical error (Theorem 4.2) and the discretization error (Theorem 4.3). Finally, we combine these two error estimates to obtain the overall error estimate of the VMEMC algorithm (Theorem 4.4).

    The full-discrete VMEMC approximate solution of (1.1) is defined as . Noting that , we can divide the overall error as

    Here the statistical error, , is estimated using the sample number . The second term, , that is corresponding to the total error of the time discretization and the FE discretization, is controlled by the time step as well as mesh size. Next, we examine the bounds of and .

    Theorem 4.2. Assume , , and hypotheses and are held, for the numerical solution to the VMEMC (3.1), then

    Proof. The process of this proof is similar to the Theorem 5 in [30].

    The is examined in the following. Recall, the true solution satisfies

    (4.10)

    The error at the time is denoted by , where is a Lagrange interpolation [2,5] of . Taking and subtracting Eq (3.1) from (4.10), respectively, we have the error equation:

    (4.11)

    here is defined as

    (4.12)

    The following regularity assumptions are needed:

    (4.13)

    As a preparatory step for estimating the discretization error, we first estimate each term of the in the following lemma.

    Lemma 4.1. Suppose every fulfills the equation (1.1) and the regularity assumptions (4.13), then,

    (4.14)

    Proof. Applying Lemma 2.1, the Cauchy-Schwarz and Young's inequalities, the Taylor's Theorem with integral remainder, we obtain the following results for is defined by Eq (4.12). The first term of is estimated as

    (4.15)

    The second and third terms of are estimated as

    (4.16)
    (4.17)
    (4.18)

    The last two terms of are estimated as

    (4.19)
    (4.20)
    (4.21)

    Combining the above estimates Eqs (4.15)–(4.21), using and , and integrating over the probability space, we can yield Eq (4.13).

    With the results of Lemma 4.1, the major convergence result of this work can be proved next.

    Theorem 4.3. Suppose satisfies the assumptions of Lemma 4.1. Moreover, suppose is an approximation of with the accuracy of the interpolation. Then, the VMEMC solution of (3.1) satisfies

    (4.22)

    Proof. Noting , rearranging the equation (4.11) with , we have

    (4.23)

    Integrating over the probability space and applying the polarization identity, we arrive at

    (4.24)

    From Lemma 4.1 with , the first term on the RHS of Eq (4.24) can be expressed as

    (4.25)

    Applying Lemma 2.1, the Taylor's theorem with integral remainder, the Young's inequality, and the Cauchy-Schwarz inequality on the second term of RHS of (4.24), we have

    (4.26)

    Using the Young's inequality and the Cauchy-Schwarz inequality to the third and fourth terms, we obtain

    (4.27)
    (4.28)

    Using the Young's inequality and the Cauchy-Schwarz inequality as well as the trace theorem to the two boundary terms, we have

    (4.29)
    (4.30)

    Under Lemma 2.1, proceeding in similar fashion with equation (4.8), we obtain

    (4.31)

    Dropping the positive terms and , using the estimate (4.25)-(4.31), noting that , , selecting for , in equation (4.24), then we can obtain,

    (4.32)

    Summing to , collecting the constants, and reorganizing the items, we obtain

    (4.33)

    Then,

    (4.34)

    Inequalities and imply the result (4.22).

    Combining the MC sampling error and the discretization error, one can acquire the total error of the VMEMC FE approximation.

    Theorem 4.4. Assume that , , satisfies the assumptions of Lemma 4.1. Moreover, suppose is an approximation of with the accuracy of the interpolation. Then,

    (4.35)

    Proof. In view of the first term on the left-hand side (LHS) of (4.35), the triangle inequality and Young inequality suggest

    Using Jensen inequality to the first term on the RHS of the above inequality, one gets

    Thus, the result can be obtained by Theorem 4.2 and Theorem 4.3. The remaining terms of the LHS of (4.35) can be handled similarly.

    Remark 4.2. Collecting these constants, the RHS of (4.35) can be expressed as

    For comparison purposes, here we use the numerical examples in our recent work [30] to test the various performances of the VMEMC algorithm. Example One is a deterministic heat transfer model with a known exact solution, which is intended to show the convergence rate of time and space, respectively. A random transient heat equation is in Example Two, which is used to test Theorem 4.4 and reveal the effectiveness of the VMEMC algorithm.

    We take as the manufactured solution; here , is a random variable, , and . The upper and lower boundary of the domain are zero Neumann boundary; the left and right boundary are Robin boundaries. The Robin coefficient and thermal conductivity diffusion factor are specified by . The source term, initial condition, and Robin boundary function should be adjusted appropriately.

    In this experiment, the ensemble schemes (3.1) and (3.2) are applied to simulate the ensemble set that includes three samples where , , . We consider . We use isometric time partition and linear FEs on uniform space partition. Moreover, we take from to with to verify the convergence order in the temporal direction. The numerical results on the VMEMC and the VAEMC schemes are shown in Table 1. Meanwhile, we select from to with to consider the convergence order in space. The resulted numerical errors on the VMEMC and the VAEMC schemes are listed in Table 2.

    Table 1.  Errors and convergence orders on temporal step.
    ,
    0.04729 0.02254 0.01091 0.00537
    Rate 1.068 1.046 1.022
    0.07742 0.03724 0.01812 0.00893
    Rate 1.055 1.039 1.020
    ,
    0.02952 0.01384 0.00659 0.00317
    Rate 1.093 1.069 1.055
    0.01504 0.00711 0.00344 0.00171
    Rate 1.079 1.048 1.009
    ,
    0.05090 0.02377 0.01132 0.00546
    Rate 1.098 1.070 1.049
    0.00600 0.00286 0.00134 0.00060
    Rate 1.068 1.093 1.150

     | Show Table
    DownLoad: CSV
    Table 2.  Errors and convergence orders on spatial step.
    ,
    0.14645 0.03808 0.00965 0.00246
    Rate 1.943 1.979 1.969
    0.14648 0.03812 0.00970 0.00251
    Rate 1.942 1.974 1.946
    ,
    0.20886 0.05419 0.01364 0.00338
    Rate 1.946 1.989 2.011
    0.20891 0.05424 0.01370 0.00344
    Rate 1.945 1.985 1.993
    ,
    0.22145 0.05743 0.01444 0.00356
    Rate 1.946 1.991 2.018
    0.22150 0.05749 0.01450 0.00362
    Rate 1.945 1.987 2.001

     | Show Table
    DownLoad: CSV

    Both the VMEMC and the VAEMC algorithms, the convergence orders on temporal direction are about one from Table 1, and coincide with theoretical findings. From Table 2, the rates of convergence on the VMEMC scheme in spatial direction are close to 2, which is higher than the theoretical result (Theorem 4.4). This may be the cause of the higher regularity of source term and boundary function in this numerical test.

    The transient heat model (1.1) with random factors is considered on . The upper and lower boundary of the domain are zero Neumann boundary; the left and right boundary are Robin boundaries. Inspired by numerical examples in literature [25,30], we choose the exact solution, the diffusion factor and Robin factor are as follow:

    The random variables are independent. is uniformly distributed on . are uniformly distributed in the interval . The Robin boundary condition, the source term, and initial conditions, are chosen to match the exact solution. Let

    (1) Convergence orders for time and physic space. Analogy to Example One, we use linear FEs on uniform partition physical space and isometric partition in time direction. The corresponding average numerical errors on VAEMC and VMEMC schemes are shown in Table 3 () and Table 4 ().

    Table 3.  Average errors and convergence orders on temporal step.
    0.01472 0.00727 0.00358 0.00180
    Rate 1.017 1.020 0.990
    0.02233 0.01134 0.00574 0.00294
    Rate 0.977 0.982 0.962

     | Show Table
    DownLoad: CSV
    Table 4.  Average errors and convergence orders on spatial step.
    0.16484 0.04296 0.01084 0.00270
    Rate 1.939 1.985 2.001
    0.18079 0.04713 0.01191 0.00299
    Rate 1.939 1.984 1.993

     | Show Table
    DownLoad: CSV

    From Table 3, the convergence order on temporal direction is according to Theorem 4.4. Looking at Table 4, the orders of convergence on the VMEMC scheme on spatial direction are close to , which is higher than the theoretical result (Theorem 4.4). This may be the cause of the higher regularity of source term and boundary function .

    (2) Convergence order for MC sampling. We take the VMEMC solution choosing samples as our benchmark, vary the values of in the VMEMC simulations, and then evaluate the approximation errors based on the benchmark. We repeat such error analysis for independent replicas and compute the average of the output errors. Denote the VMEMC solution at time in the th independent replica by where is result of the VMEMC (3.1) in the th experiment. Denote

    We run times and take the logarithm for the numerical results of at are depicted in the left subgraph of Figure 1. One group of values of can be seen in Table 5 and the corresponding linear regression model is shown in the right subgraph of Figure 1.

    Figure 1.  Convergence order of the VMEMC algorithm for sampling number . Left: the numerical results of . Right: the linear regression for one group of numerical results of .
    Table 5.  One group of values of ().
    10 20 40 80 160
    0.03431 0.01763 0.02839 0.01077 0.00661

     | Show Table
    DownLoad: CSV

    We apply linear regression analysis on the numerical errors, which show

    These results shows that the rate of convergence with respect to is close to , which coincides with our theoretical results in Theorem 4.4.

    (3) Comparison. Using the mesh size and the sample size , we obtain the mean and the variance of the solutions at . The results are plotted in Figure 2 and Figure 3 for different algorithms.

    Figure 2.  Three means. Left: mean of the VMEMC solution. Middle: mean of the VAEMC solution. Right: mean of the IFE-MC solution.
    Figure 3.  Variances. Left: variance of the VMEMC solution. Middle: variance of the VAEMC solution. Right: variance of the IFE-MC solution.

    From Figure 2 and Figure 3, it is seen that both VMEMC and VAEMC approximations are excellent.

    To check out the performance of the VMEMC algorithm, we compare the result with that of individual finite element MC (IFE-MC) simulations applying the same groups of sample values. The difference from the mean to the IFE-MC solution is plotted in Figure 4.

    Figure 4.  Difference of the VMEMC mean and the IFE-MC solution.

    The accurate approximation to the IFE-MC implemented for the VMEMC scheme is very nice. From Figure 4, the difference on the order of signifies that the VMEMC approximation can also achieve basically accurate approximation as the IFE-MC implements.

    (4) Efficiency. To test the efficiency, setting the same time as well as space size, we record the CPU running times of the VMEMC, VAEMC, and MC algorithms in Table 6. Here we apply LU factorization because the discrete system's size is small.

    Table 6.  CPU time (s) ().
    10 50 100 200
    CPU time for VMEMC 37.65 150.21 295.59 569.66
    CPU time for VAEMC 38.37 152.47 285.91 607.01
    CPU time for MC 65.27 288.10 559.29 1100.71

     | Show Table
    DownLoad: CSV

    From Table 6, when , we can show that both the VMEMC and the VAEMC algorithms take fewer CPU times than the MC algorithm. In contrast with the non-ensemble scheme, these ensemble methods improve the calculational efficiency by about 70%–80%. As the size of the ensemble increases, the efficiency of these ensemble algorithms improves, particularly for the VMEMC algorithm.

    The core idea of the ensemble algorithm is to develop a semi-implicit and semi-explicit temporal discretization scheme for PDEs. It will select a representative of the random coefficient as the implicit term, and take the fluctuation between random coefficients and this representation as the explicit term. This results in a shared coefficient matrix for each ensemble member at every time step. Consequently, storage requirements are reduced, and solution speed is increased. In general, as long as the relative fluctuation is not very large, such as , the ensemble numerical scheme can be stable. A large amount of literature [17,18,19,20,23,24,30] discussed the numerical stability conditions and error analysis of the ensemble algorithm based on taking the mean as representation .

    Based on taking the maximum value of the random coefficient as representation, an unconditionally stable ensemble scheme, VMEMC, is established in this work. First, we apply the VMEMC algorithm to deal with the transient heat equation with the stochastic diffusion factors and the Robin factors. Without any stability conditions, we prove that the algorithm is numerically stable. Additionally, we also provide the overall error estimate for the algorithm. Moreover, numerical experiments are conducted to demonstrate the application of the ensemble schemes and to validate the established properties of the algorithm.

    Important next steps include accommodating phase changes in the solid material (e.g., transitions to a liquid phase) and incorporating additional physical processes into the boundary conditions (e.g., surface-to-ambient radiation). The stochastic Robin boundary control problem is worth investigating. The VMEMC method can also be developed for nonlinear stochastic parabolic equations.

    Tingfu Yao: Methodology, writing original draft, software; Changlun Ye: Software and supervision; Xianbing Luo: Review, editing, and supervision; Shuwen Xiang: Supervision and validation.

    The authors declare they have not used Artificial Intelligence (AI) tools in the creation of this article.

    This work is supported by the National Natural Science Foundation of China (Granted No. 11961008, 12361033) and The Natural Science Research Foundation of Education Department of Guizhou Province (Granted NO. QJJ[2024]190).

    The authors declare there is no conflict of interest.


    Acknowledgments



    Getinet Gawo would like to thank Pennsylvania State University, York Campus for the financial support provided for this work through York Campus Advisory Board Scholarly Activity Grant.

    Conflict of interest



    All authors have no conflict of interest.

    Authors contribution



    A.M Siddiqui provided the formulation of the problem, defining governing equations and boundary conditions and supervised whole paper. G. Gawo solved the governing equations and results are compared with previous works, and results are also implemented on rat kidney data. K.Maqbool revised the literature review, implemented the result on human kidney and draft the result and handled the paper submission process.

    [1] Macey RI (1963) Pressure flow patterns in a cylinder with reabsorbing walls. B Math Biophys 25: 1-9. doi: 10.1007/BF02477766
    [2] Berman AS (1953) Laminar flow in channels with porous walls. J Appl Phys 24: 1232-1235. doi: 10.1063/1.1721476
    [3] Yuan SW (1955) Laminar pipe flow with injection and suction through a porous wall. Princeton Univ nj James Forrestal Research Center 78: 719-724.
    [4] Yuan SW (1956) Further investigation of laminar flow in channels with porous walls. J Appl Phys 27: 267-269. doi: 10.1063/1.1722355
    [5] Terrill RM (1982) An exact solution for flow in a porous pipe. Z Angew Math Phys 33: 547-552. doi: 10.1007/BF00955703
    [6] Macey RI (1965) Hydrodynamics in the renal tubule. B Math Biophys 27: 117. doi: 10.1007/BF02498766
    [7] Gilmer GG, Deshpande VG, Chou CL, et al. (2018) Flow resistance along the rat renal tubule. Am J Physiol-Renal 315: F1398-F1405. doi: 10.1152/ajprenal.00219.2018
    [8] Achala LN, Shreenivas KR (2011) Two dimensional flow in renal tubules with linear model. Adv Appl Math Biosci 2: 47-59.
    [9] Siddiqui AM, Haroon T, Shahzad A (2016) Hydrodynamics of viscous fluid through porous slit with linear absorption. Appl Math Mech 37: 361-378. doi: 10.1007/s10483-016-2032-6
    [10] Navier C (1823) Thesis on the Laws of Motion of Fluids. Memoirs of the Royal Academy of Sciences of the Institut de France 6: 389-440.
    [11] Darrigol O (2002) Between hydrodynamics and elasticity theory: the first five births of the Navier-Stokes equation. Arch Hist Exact Sci 56: 95-150. doi: 10.1007/s004070200000
    [12] Mooney M (1931) Explicit formulas for slip and fluidity. J Rheol (1929–1932) 2: 210-222. doi: 10.1122/1.2116364
    [13] Chauffoureaux JC, Dehennau C, Van Rijckevorsel J (1979) Flow and thermal stability of rigid PVC. J Rheol 23: 1-24. doi: 10.1122/1.549513
    [14] Lau HC, Schowalter WR (1986) A model for adhesive failure of viscoelastic fluids during flow. J Rheol 30: 193-206. doi: 10.1122/1.549888
    [15] Cohen Y, Metzner AB (1985) Apparent slip flow of polymer solutions. J Rheol 29: 67-102. doi: 10.1122/1.549811
    [16] Hatzikiriakos SG, Dealy JM (1991) Wall slip of molten high density polyethylene. I. Sliding plate rheometer studies. J Rheol 35: 497-523. doi: 10.1122/1.550178
    [17] Hatzikiriakos SG, Dealy JM (1992) Wall slip of molten high density polyethylenes. II. Capillary rheometer studies. J Rheol 36: 703-741. doi: 10.1122/1.550313
    [18] Rao IJ, Rajagopal KR (1999) The effect of the slip boundary condition on the flow of fluids in a channel. Acta Mech 135: 113-126. doi: 10.1007/BF01305747
    [19] Elshahed M (2004) Blood flow in capillary under starling hypothesis. Appl Math Comput 149: 431-439.
    [20] Singh R, Laurence RL (1979) Influence of slip velocity at a membrane surface on ultrafiltration performance—I. Channel flow system. Int J Heat Mass Tran 22: 721-729. doi: 10.1016/0017-9310(79)90119-4
    [21] Chu ZKH (2000) Slip flow in an annulus with corrugated walls. J Phys D: Appl Phys 33: 627. doi: 10.1088/0022-3727/33/6/307
    [22] Beavers GS, Joseph DD (1967) Boundary conditions at a naturally permeable wall. J Fluid Mech 30: 197-207. doi: 10.1017/S0022112067001375
    [23] Priezjev NV, Darhuber AA, Troian SM (2005) Slip behavior in liquid films on surfaces of patterned wettability: Comparison between continuum and molecular dynamics simulations. Phys Rev E 71: 041608. doi: 10.1103/PhysRevE.71.041608
    [24] Palatt PJ, Sackin H, Tanner RI (1974) A hydrodynamic model of a permeable tubule. J Theor Biol 44: 287-303. doi: 10.1016/0022-5193(74)90161-1
    [25] Siddiqui AM, Haroon T, Kahshan M, et al. (2015) Slip effects on the flow of Newtonian fluid in renal tubule. J Comput Theor Nanos 12: 4319-4328. doi: 10.1166/jctn.2015.4358
    [26] Kapur JN (1988)  Mathematical Modelling New Age International.
    [27] Curthoys NP, Moe OW (2014) Proximal tubule function and response to acidosis. Clin J Am Soc Nephro 9: 1627-1638. doi: 10.2215/CJN.10391012
  • Reader Comments
  • © 2021 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(2851) PDF downloads(165) Cited by(3)

Figures and Tables

Figures(8)  /  Tables(4)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog