A Viscous Fluid Flow through a Thin Channel with Mixed Rigid-Elastic Boundary: Variational and Asymptotic Analysis
R. Fares,1G. P. Panasenko,1and R. Stavre2
Academic Editor: D. O'Regan
Received23 Jan 2012
Accepted29 Feb 2012
Published21 Jun 2012
Abstract
We study the nonsteady Stokes flow in a thin tube structure composed by two thin rectangles with lateral elastic boundaries which are connected by a domain with rigid boundaries. After a variational approach of the problem which gives us existence, uniqueness, regularity results, and some a priori estimates, we construct an asymptotic solution. The existence of a junction region between the two rectangles imposes to consider, as part of the asymptotic solution, some boundary layer correctors that correspond to this region. We present and solve the problems for all the terms of the asymptotic expansion. For two different cases, we describe the order of steps of the algorithm of solving the problem and we construct the main term of the asymptotic expansion. By means of the a priori estimates, we justify our asymptotic construction, by obtaining a small error between the exact and the asymptotic solutions.
1. Introduction
A few years ago we began to publish several papers dealing with an asymptotic analysis for an interaction problem between a viscous fluid and an elastic structure. Problems involving the interaction between a fluid and a deformable structure have been studied extensively in the last years due to their applications in many areas such as: engineering, biomechanics, biology, and hydroelasticity.
In order to model the blood flow through a vessel we considered in [1] a non steady viscous flow in a thin rectangle with elastic walls, when at the ends of the flow domain periodicity conditions are imposed. The asymptotic approach was then extended for the nonperiodic case in [2], when the inflow and outflow velocities are given functions. Generalizations of the previous cited papers were obtained in [3โ5]. In the first two papers we extended the results of [1] to a-three dimensional case, by considering a periodic, axisymmetric flow in a thin cylinder with elastic lateral boundary, while in [4] we studied two different cases: when the inertial term of the equation for the wall displacement is much smaller than the stress term and when these two terms are of the same order. Moreover, [5] deals with the asymptotic analysis for a fluid with variable viscosity.
In all these papers, the fluid flow was described by the Stokes equations and the behavior of the elastic structure was simulated either by the Sophie Germainโs or by the Koiterโs equation. The fluid-structure interaction was mathematically expressed by the equality of the fluid velocity at the boundary and the time derivative of the wall displacement. We constructed an asymptotic solution and we proved that it represents a good approximation for the exact solution, by obtaining a small error between them.
In this paper we consider the nonsteady Stokes flow in a thin tube structure. In two-dimensional case, a tube structure is some connected union of thin rectangles, having a dimension much smaller than the other one. The main difference between this paper and our previous works [2, 5] is that now the boundary of the flow domain contains elastic parts and rigid parts, as well.
The flow domain consists of two thin rectangles with elastic boundary. The junction between the rectangles is realized by means of rigid boundaries. The interaction between the viscous fluid and the elastic boundaries produces normal displacements. The elastic boundaries behavior is described by the Sophie Germainโs equation.
This domain models a vessel structure where a stent was placed. In fact, for treating arterial stenoses or occlusions, percutaneous angioplasty is indicated. It uses small inflatable balloons, single-use, that help dilate the artery at the site of narrowing. During angioplasty, a small wire mesh tube called a stent may be permanently placed in the newly opened artery or vein to help it remain open.
There are two types of stents: bare stents (wire mesh) and covered stents (also commonly called stent grafts). The first are in use as part of the expansion of retrecissements vessels, the second impermeable and can be used to prevent or treat hemorrhage. Stents are used for arteries in the heart, the kidney, the arm, or the leg. They are also used for the aorta in the abdomen or chest (see Figures 1 and 2).
We assimilate then the part of the vessel where the stent is placed to a domain with rigid boundaries.
We suppose that the viscous fluid has a variable viscosity depending on a longitudinal variable for each rectangle. This situation models a blood flow in a vessel structure where the viscosity depends on the concentration of some substances diluted in blood or some blood cells. Indeed, the asymptotic analysis of the convection-diffusion equation set in such domains [6, 7] shows that in the case of the Neumann (impermeability) condition at the lateral boundary and small Reynolds numbers, the concentration is asymptotically close to the one-dimensional description, that is, the convection-diffusion equation set on the graph. The solution of the problem on the graph is the leading term of the asymptotic expansion, and it evidently depends on the longitudinal variable. On the other hand, the viscosity often depends on the concentration of the diluted substances or distributed cells, and so, it depends on the longitudinal variable. Of course, the fluid motion equation is coupled with the diffusion-convection equation in this case. However, if the velocity is small (in our case, it is of order ), then neglecting the convection, in comparison with the diffusion term or iterating with respect to the small term,we get the steady state diffusion equation; in absence of the source term in the right-hand side, it has a piecewise-linear asymptotic solution on the graph for the concentration. So, in this simplified situation, the diffusion equation can be solved before the fluid motion equation, and we obtain for the flow, the Stokes or Navier-Stokes equation with a variable viscosity depending (via concentration) on the longitudinal variable.
There are, of course, many other practical problems involving fluids with variable viscosity. For example, the presence of bacteria in suspension (see [8]) may change locally the viscosity.
The outline of the paper is as follows. In Section 2 we give a description of the flow domain and of the coupled system which models our problem. The next section presents the variational formulation of the problem. The literature contains an important number of papers dealing with a variational approach of fluid-structure interaction problems. For instance, results concerning the existence of weak or strong solutions when the fluid domain is either fixed or depending on time can be found in [9โ12]. We establish in this section results such as: existence, uniqueness, regularity and a priori estimates. In Section 4 we construct the asymptotic solution. As in our previous papers [1โ5], the problem depends on two small parameters. The first small parameter, , is defined as the ratio of the dimensions of the thin rectangles; the second one, , corresponds to the softness of the wall. For various values of the small parameters and , an asymptotic expansion of the solution is constructed; the parameter is taken of the form , with , . The asymptotic expansion is different for the cases: and . The asymptotic solution contains three types of terms: the regular part, defined as in [1], the boundary layer correctors corresponding to the boundary conditions and the boundary layer correctors which realize the junction between the motion in the two rectangles. The first two types of correctors have already been introduced in [1, 2]. The third type is characteristic for structures with junction regions. The asymptotic solution is more complicated in this case since it contains also some truncation functions introduced in order to restrict the influence of the boundary layer correctors to the regions to which they correspond. We present and solve the problems for all the components of the asymptotic solution. For the two cases and the order of solving the problems is presented and the leading term of the asymptotic expansion is described. In the last section, we establish the error between the exact solution and the asymptotic one, by means of the a priori estimates obtained in Section 3. The small error between the two solutions justifies our asymptotic expansion.
2. The Physical Problem
We consider an incompressible, viscous fluid, with variable viscosity, occupying a thin domain, . The flow domain is a thin tube structure, composed by two thin rectangles with lateral elastic boundaries, connected by a region with rigid boundaries. We introduce the first small parameter of our problem, , , in connection with the ratio of the two dimensions of the rectangles, as below. The thin rectangles are given by:
and the junction region is
The flow domain is given by , as shown in Figure 3.
The main difference between this paper and our previous works [2, 5] is that now the boundary of the flow domain contains elastic parts and rigid parts as well. Let and be the rigid parts of the boundary of defined as follows:
the elastic parts of are given by:
Let and be the inflow and outflow parts of the boundary of . Denoting , we can write .
We study the nonsteady, slow flow of the viscous fluid in the domain previously described, when the fluid interacts with the elastic boundaries, . The interaction between the fluid and the elastic boundaries produces the normal displacements . We neglect the tangential displacements and we consider that the elastic boundaries are clamped. We study the problem for , with an arbitrary positive constant independent on and we assume that the membranes are not very elastic so that the displacement of the boundaries is small enough. Consequently, at each time , we approximate the position of the elastic membranes by its initial position and, hence, the fluid flow equations are considered in the initial configuration. We suppose that the displacements have the following form:
The problem described above, with nonhomogeneous boundary conditions for the velocity, is modeled by the following coupled system:
with the outer unit normal on the boundary of and the unit tangent vector to . The given data contained by the previous system are: some material constants and some given functions. are positive given constants in connection with the properties of materials representing the density of the fluid, the density of the elastic walls, a viscosity coefficient, and the two-dimensional Youngโs modulus (The two dimensional Youngโs modulus is defined as where is the common three dimensional Young modulus and is the Poisson ratio.), respectively, and the positive constant stands for the thickness of the elastic walls.
The given functions are: , the variable viscosity of the fluid, which satisfies for all , , the exterior force applied to the fluid, , the exterior forces applied on the elastic walls, with ,
and a small inflow-outflow velocity defined by means of the function defined as follows:
with . The function is the trace of a function denoted also by with the following properties:
where . The unknowns of the system (2.6) are the velocity of the fluid, , the pressure of the fluid, , and the normal displacements of the elastic walls, .
The fluid flow is described by the nonsteady Stokes equations. For the normal displacements we consider the Sophie-Germainโs equation. A โviscousโ type term, , is added to the usual forth-order equation for the normal displacements. It corresponds to the viscoelastic behavior of the wall (the so called Kelvin-Voigt model). Usually, the Youngโs modulus, , has a value of โPa. On the other hand, we assume that the characteristic longitudinal space scale for vessels is of order of cm and the time scale is of order of seconds. Let us use the S.I. system of units. This leads us to the necessity of scaling of every derivative is by the factor ; that is, the fourth derivative will contain the additional factor . If is of order โm or โm, then the coefficient can be taken in the further analysis as a value of order of 1. The coefficient in (2.6) will be replaced (after scaling in ) by a great coefficient with the value of order of to . If the ratio of thickness and the length of the vessel are of order , then is of order from to . We assume that the โviscousโ term is much smaller than the term with the coefficient and hence the new coefficient denoted also by , obtained after scaling in , is . More details concerning (2.6) can be found, for instance in [1].
Due to the properties of the function , the compatibility condition for the coupled system which describes the physical problem is
Using next the initial condition for the displacements, condition for the above coupled system becomes
This condition states that the global area of the flow domain is preserved.
3. Variational Formulation: Existence, Uniqueness, Regularity, and A Priori Estimates
In order to obtain the above properties for the solution of the physical problem, we introduce the variational framework.
To simplify the computations, we consider first (2.6) with homogeneous boundary conditions on and , that is, the problem for . Then, the same properties for the solution of (2.6) follow with the usual technique for nonhomogeneous problems.
Taking into account the conditions for the velocity and for the displacements and the condition (2.11) we introduce the following spaces:
Choosing for the data the regularity and we consider the following variational problem:
Here and below and .
For the nonhomogeneous boundary conditions we still obtain the variational formulation (3.2) with replaced by and replaced by .
Theorem 3.1. The variational problem (3.2) has the unique solution with .
Proof. Let us start with the proof of the uniqueness for the solution of (3.2). Consider and two solutions of (3.2) and define . Subtracting the two relations (3.2) and taking as test function we get:
Integrating from 0 to this equality and taking into account the initial conditions, we obtain: a.e. in and a.e. in . Hence the problem (3.2) has a unique solution. For proving the existence and the regularity of the functions and , we will use the Galerkinโs method. We begin with the construction of a basis for the space . Let be a basis of chosen by considering the eigenfunctions of the following problem:
where is the fourth derivative of and , for all . We define as follows: where
It is easy to check that is a basis for . We choose the functions of the basis such that
As a consequence of the previous relation we also get:
We consider now a basis of the space , constructed with the eigenfunctions of the following Stokes problem:
with for all . The functions are uniquely determined from the condition
which implies
Next, for any we consider the following problem:
Following the classical results of [13] for nonhomogeneous Stokes problems, we obtain a unique and a function unique up to an additive constant. Moreover, for any
By means of the functions , we are now in a position to define, for each , an approximate solution of (3.2) as follows:
with scalar unknown functions. These functions are determined below from the condition that is the solution for the problem:
We introduce the notations:
Taking into account the previous notations and the relations (3.6), (3.7), (3.9), (3.10), (3.12) we get from (3.14) the following linear differential system for the unknown functions :
The previous system uniquely determines the unknown functions . For more details see [1, Sectionโโ3]. In the sequel we obtain some a priori estimates which give the regularity of the solution for (3.2). Computing (3.14)(3.14) and using (3.13) we get:
Integrating from 0 to , using the property of and the initial conditions we obtain, as in [1], the first estimates:
with and is the variable of the parametrisation for or on . The second estimates are obtained computing (3.14)(3.14):
From (3.18) and (3.19) we get the boundedness of in . Consequently, we can pass to the limit in (3.14) which yields the existence result of the theorem. To achieve the proof, we notice that the regularity stated in (3.2) follows from the estimates (3.18) and (3.19).
As a consequence of Theorem 3.1 we introduce the pressure which appears in (2.6).
Corollary 3.2. There exists a unique function such that satisfies (2.6) a.e. in and on , respectively.
To obtain this result we follow the ideas of [1, Sectionโโ3].
The last result of this section presents the estimates for , the unique solution of (2.6).
Corollary 3.3. Let be the solution of the problem (2.6) corresponding to the data , with and . Then the following estimates hold:
Remark 3.4. If we consider the nonhomogeneous problem (which corresponds to ) we obtain for the corresponding solution the same estimates (3.20) with a different , but which depends in the same way on .
4. Asymptotic Analysis
In order to approximate the solution of (2.6), we consider more regular data than in Section 3. We suppose that:
and , and are chosen as follows:
with constant.
4.1. Construction of the Asymptotic Solution
In the sequel we introduce the second small parameter and we take , with , . The asymptotic solution approximating the periodic flow in an infinite rectangle (see [1, Sectionโโ5]) is modified by using two types of correctors: the first type corresponds to the boundary conditions on and and the second one represents the boundary layer functions in a neighborhood of . So, the asymptotic solution contains the regular part of the solution (which has two terms, due to the form of the flow domain), two boundary layer functions corresponding to and to , respectively, and the correctors in . Since the terms of the asymptotic solution should vanish in different parts of the flow domain, we multiply them with the truncation functions defined as follows: , , ,
We construct the asymptotic solution as below:
Due to the definition of the truncation functions, we notice that in , that is, in a neighborhood of the region of variable viscosity, the asymptotic solution reduces to its regular part, , in a neighborhood of it reduces to , , while in the asymptotic solution is equal to the corrector in . This means that the regular part of the asymptotic solution (which can be computed and has a very simple expression) represents an approximation for the exact solution in all the flow domain except some neighborhoods of , and .
We give next the expressions of the three components of the asymptotic solution and we explain their role in this construction.
4.1.1. The Regular Part of the Asymptotic Solution
For each rectangle and we define the corresponding regular part of the asymptotic expansion as in [1]. The regular part corresponding to has the expression:
It represents the solution of the problem (2.6) set in the infinite rectangle in direction, ; the functions which appear in (4.5) can be computed explicitly and represent a good approximation of the exact solution of (2.6) in a neighborhood of the region with variable viscosity of (as we will prove in the last section).
In a similar way we introduce the regular part of the asymptotic solution corresponding to .
Hence, the sum of the first two terms of the asymptotic solution represents a good approximation of the exact solution of (2.6) in some neighborhood of the region of variable viscosity of the flow domain, but it is not close to it on the inflow/outflow boundaries and in . Since the purpose of the asymptotic construction is to approximate the exact solution with a small error between the exact and asymptotic solution, we modify and complete the regular part by two types of boundary layer correctors.
4.1.2. The Boundary Layer Correctors for and
These boundary layer functions are introduced in order to repair the traces of the regular part of the asymptotic solution on and . They are given by:
The corrector with corresponds to the end and that with corresponds to .
From the definition of the truncation functions and of the asymptotic solution we notice that the influence of each corrector defined above is significant only near the corresponding end of the flow domain.
4.1.3. The Boundary Layer Corrector in
This corrector is necessary in order to realize the junction between the two parts of the asymptotic solution, corresponding to the two branches of the flow domain and to obtain the conditions on the rigid boundaries of the domain. The expressions of the correctors in are given by:
From the definition of the truncation function it follows that the corrector in appears in the expression of the asymptotic solution (4.4) only in , which represents a neighborhood of with constant viscosity and where the forces are equal to zero.
Remark 4.1. The boundary layer method is close to the well-known method of matching of asymptotic expansions (see [14]). However there is a difference between these two methods. In the boundary layer method the regular expansion is reexpanded with respect to the fast variable (as in the matching method), but then the new expansion is multiplied by a cutting function. This product being inserted into the equation gives a discrepancy in the right-hand side which is then compensated by an appropriate boundary layer corrector. So, we get finally a unique expression for the asymptotic approximation in each point of the domain. In the matching approach the reexpanded regular ansatz (outer expansion) coexists with the so-called inner expansion defined near the boundary, so that there are several overlapping expressions for an asymptotic approximation in different parts of the domain. We emphasize that this difference is not too important and corresponds more to the form of presentation of the result.
4.2. The Determination of the Asymptotic Solution
This subsection is devoted to the resolution of the problems satisfied by the three different components of the asymptotic solution presented in the previous subsection. Since some computations are different with respect to the values of , we will analyse the problems and the order of solving them for and for . Moreover, in each case we will specify the leading term of the asymptotic expansion.
We begin the approach with the problems for the boundary layer correctors corresponding to and since the study of these problems is the same both for and for . To fix the ideas, we obtain in the sequel the problems for the corrector corresponding to the end . As we noticed before, the term containing this corrector is not equal to zero only in a neighborhood of the boundary . So, the problems and the other relations corresponding to this corrector are obtained substituting the asymptotic solution in (2.6). We obtain two separate problems: one for the velocity-pressure correctors and the other one for the displacements correctors.
Since in this neighborhood the viscosity is constant, the problem for has constant coefficients. Denoting by the semi-infinite rectangle and imposing for the velocity and for the pressure the condition of decay at , we obtain for the problem:
The compatibility condition for (4.8) reads:
For any the right-hand side of (4.8) is known; so the boundary layer correctors for the velocity and for the pressure corresponding to are uniquely determined from (4.8) (see [2, Sectionโโ4]). The condition (4.9) represents a relation for the regular part of the asymptotic solution.
The boundary layer correctors for the displacements exponentially stabilizing to zero at are obtained as the unique solution of the following problems:
Since at the step the problem (4.10) gives both and , introducing the asymptotic solution into (2.6) we obtain two boundary conditions for the regular part of the asymptotic solution for the displacements:
In a similar way, we obtain the boundary layer correctors corresponding to the end . The boundary layers for the velocity-pressure are defined on , with , and the boundary layers for the displacements are defined also on .
We study next the problems for the regular parts of the asymptotic solution. The results are obtained for the regular part corresponding to ; the regular part corresponding to may be obtained from the previous with some obvious changes.
Introducing (4.5) into (2.6) and collecting together the terms of the same order with respect to we are leaded to consider the following problem for :
The two cases, and , appear because of the last relation of the previous system. We can see, indeed, that for the unknown of this relation is , while for , (4.12) contains two unknowns. From this point, the computations are different with respect to the values of .
We introduce the functions:
with the properties: and .
We also use the notations:
4.2.1. The Order of Solving the Problems for
The regular part corresponding to is computed by integrating (4.12), as stated below.
Proposition 4.2. The unknowns are determined from (4.12), up to nine functions of .
Proof. Integrating twice (4.12) from โ1 to and using the boundary conditions (4.12) we get:
which contains as unknowns and . The other functions contained by this relation are either known from previous computations or equal to zero. We integrate next the incompressibility condition (4.12) with respect to with the boundary condition (4.12) for and we obtain
The previous two relations give the components of the velocity with respect to . The pressure approximations are determined from (4.12), supposing that the integration functions, depending on , are equal to zero, since we consider that any function depending only on , could be contained in .
Taking in (4.16) and using the boundary condition (4.12) for we obtain the following second order differential equation for the function :
Integrating (4.18) from to 1, we express by means of , which represents the only unknown of this expression. This function of is obtained as follows: we take in (4.15) and we introduce the result into (4.9). Hence, we determined the expression of in , which is:
Introducing (4.19) into (4.15) and (4.18) into (4.16) we determine in . We integrate next โยทโ(4.19) from 0 to and we get determined up to the function of , :
The functions satisfy the fourth-order differential equations:
with given by (4.20). Writting as and integrating four times (4.21) with respect to we obtain the following expressions for :
Hence, the regular part of the asymptotic solution corresponding to is determined up to the functions , , , , , which achieves the proof.
In a similar way we express the regular part of the asymptotic solution corresponding to depending on 9 unknown functions of .
We continue the construction with the problems for the corrector in . As we previously said, the term of the asymptotic solution containing this corrector is not equal to zero only in a neighborhood of . In this neighborhood, the expression of the asymptotic solution reduces to . For obtaining the problems satisfied by the correctors in , we introduce the previous expression of the asymptotic solution in (2.6), with and in (2.6); for derivating the terms which contain the two types of variables and we proceed as follows: we replace by , we replace by and we expand the functions as a Taylor expansion with respect to . We introduce the notations and we obtain for the following nondivergence free problem:
for and
Unlike the problem (4.8) which give the correctors for , the problems (4.23), (4.24) have unknown right hand sides. The functions are known, but and are elements from the next approximation, so they are unknown. The function is given by:
and it has to satisfy a compatibility condition.
The problem (4.24) has the unique solution (with the pressure unique up to an additive function of )
We obtain next the problems for the correctors in corresponding to the displacements. We first notice that
Let us denote by the unique solution of
and by the unique solution of
Introducing the asymptotic solution into (2.6) and identifying the coefficients of the power of we obtain:
The relation (4.30) gives the corrector in for the displacements corresponding to one branch of the flow domain.
We notice that the right-hand side of (4.31) is known, since the problems (4.28) and (4.29) can be solved also if we replace by .
We are now in a position to return to the regular part of the asymptotic solution. For each regular part of the asymptotic solution, we need 9 relations, in order to determine the 9 unknown functions of integration. For the regular part corresponding to , (4.11) and (4.31) give 8 functions: , , , . It remains undetermined . We show next that at the jth step this function is already determined at the previous step. For this purpose, we return to the problem (4.23). As usual, we should construct boundary layer functions stabilizing to zero at infinity. As we imposed this condition for the velocity in (4.23), it is known that the pressure stabilizes to some functions of . At the jth step, we determine and from the condition of the exponential decay of at infinity. To this end, we consider instead of (4.23) and (4.24), the following problem with known right-hand side for :
with unknown. Supposing for the time being that the function satisfies the compatibility condition, it follows that (4.32) has a unique solution (the uniqueness for the pressure being understood up to an additive function of ). This means that the function , with is fixed. We define for
Standard computations show that is solution for (4.32) iff satisfies (4.23) or (4.24) for with uniformly when . To obtain the desired behavior for when , it suffices to take .
A first equation for the unknowns and is
with . The other equation is obtained introducing the asymptotic solution into the compatibility condition (2.11):
The right-hand side of (4.35) is known from the jth approximation. Replacing in the left hand side of (4.35) given by (4.22) (for replaced by ) and given by a similar relation, we are leaded to the second equation for the unknowns and :
with a known function for , with
From (4.34) and (4.36) we determine and , which means that the right-hand side of (4.23) is now known, that allows to solve (4.23), if we show that the problem (4.23) is well posed, that is, the compatibility condition is fulfilled. Moreover, this justifies the previous assertion that are computed at the ()th step. We can say that the entire th approximation was obtained for if we prove.
Proposition 4.3. For any , the function satisfies the compatibility condition:
Proof. From (4.23) and (4.35) it follows:
We express next the three terms of the right-hand side of the previous equality with respect to the components of the velocity corresponding to the regular parts of the asymptotic solution. From (4.9), the corresponding relation for and (2.9) the first term of the right-hand side is equal to
The third term is expressed by means of (4.12) and its analogous for and it is equal to
Computing now and using (4.12) and the similar one for we get
The second term is given by
Finally, computing , we obtain the following expression of the right-hand side of (4.38):
We compute next the left-hand side of (4.38). Since on , we can write:
We express the first term of the right-hand side of the previous relation taking into account that for , we have :
For obtaining the above equality we also used the property:
A similar result is obtained for the second term of the right-hand side and, comparing the expression of with the right-hand side of (4.44), the proof is complete.
The last result of this subsection gives the leading term of the regular part of the asymptotic solution corresponding to .
Proposition 4.4. For the regular part of the asymptotic solution corresponding to is given by
with defined by (4.37) and
4.2.2. The Solution of the Problems for
Unlike in the previous case, for we cannot solve the problems corresponding to and to separately. We begin with the problems corresponding to the functions with the index 1. The relations (4.15)โ(4.17) are still true. Hence, is also given by (4.17). The difference with respect to the previous case is that now we cannot express the other unknowns only in function of , since in the right-hand side of (4.16) is also unknown.
We introduce the notations:
with given by (4.14). The function is obtained as the unique solution of a fourth-order differential problem:
the boundary conditions being obtained by means of (4.11) and (4.31). This means that can also be expressed via . We show that all the unknowns of the regular part of the asymptotic solution can be expressed in function of ; hence, the regular part of the asymptotic solution corresponding to is determined if we obtain and solve a problem for . As we cannot obtain all the necessary conditions for , we have to consider the problem for the couple of unknowns, , .
Theorem 4.5. The approximations of the displacement, , , are obtained as the unique solution of the following system of two sixth-order parabolic equations:
where , , , , and are known functions, given below.
Proof. The first relation between the unknowns and is given by (4.12), that is,
The second one is obtained from (4.16) for and , (4.12) for and (4.50):
Eliminating from (4.53) and (4.54) we get (4.52) with
As we can see, depends only on the data and depends on some functions determined at the previous approximations. Similar expressions are obtained for and . Relations (4.52) are in fact (4.11) and (4.31), respectively, and the same remark holds for (4.52). We obtain next (4.52) as follows: we derivate (4.53) with respect to and we take ; we also take in (4.15) and we eliminate from these two relations, by means of (4.9). The right-hand side of (4.52) is
and a similar expression we obtain for . Relation (4.52) is a consequence of (4.34) for replaced by , of (4.53) for and of the corresponding relation for the index 2. The right hand side of (4.52) is
Finally, (4.52) is a consequence of the compatibility condition for (4.32), obtained as follows: (4.25) can be written as
where
from (4.32) we get
Next, we derivate (4.53), we take and we add with the same computation for the indexโโ2. Replacing the previous relation, we obtain (4.52), with
which completes the proof.
The other components of the regular part of the asymptotic solution corresponding to have the following expressions:
As in the previous subsection, we discuss the expression of the leading term of the asymptotic solution.
Proposition 4.6. For the regular part of the asymptotic solution corresponding to is given by
with satisfying the problems presented in the proof, (4.64) and (4.65).
Proof. The relations (4.63) are consequences of (4.15)โ(4.17) and of (4.53) for . It remains to show that can be determined as unique solutions of some problems. For the problem for becomes:
Finally, the problem for is the following:
At the end of this subsection we notice that the system (4.52) can be solved also for replaced by which gives, by means of (4.53), , . Hence, the system (4.23) can be solved, and all the th approximation is determined for .
4.2.3. Conclusions
In the previous two subsections we determined the th approximation of the asymptotic solution for two different cases: and . For we can separate the problems for the two branches of the elastic tube and we obtain the regular part of the asymptotic solution (for ) as follows:(i)the macroscopic variable is the unique solution of the problem:
with known functions;(ii)the macroscopic variables are the unique solutions of the fourth-order differential problems:
and a similar one for ;(iii) given by (4.15)โ(4.17). The main term of the regular part of the asymptotic solution is presented in (4.48)-(4.49). Unlike in the case , for the problems must be solved for the two branches together, since they cannot be uncoupled;(iv)we first introduce the auxiliary functions , with the unique solution of (4.51); then the regular part of the asymptotic solution is obtained as follows;(v)the macroscopic variables are obtained as the unique solution of the sixth-order parabolic system (4.52);(vi)the macroscopic variables are given by (4.53) and the similar equation for the index 2;(vii) given also by (4.15)โ(4.17), but we obtain different expressions.
The main term of the regular part of the asymptotic solution is given by (4.63) and (4.65).
5. Error Estimates
In the last section, we estimate the error between the exact solution and the asymptotic one, in order to justify the asymptotic expansion. Introducing the asymptotic solution of order , , into (2.6) we want to obtain for it a problem of the same type as (2.6).
The next computations hold for any integer . We begin with the continuity equation satisfied by the asymptotic velocity of order . Standard computations give
with
We notice that the asymptotic velocity of order is not a divergence free function.
We present next the equations for the asymptotic solution, which correspond to (2.6) and (2.6):
with residual known functions. Since the expressions of these residuals are rather complicated, but not too technical to obtain, let us give below only the estimates for their norms in the corresponding spaces: ,
These estimates will be used to derive finally the error estimate between the exact and the asymptotic solution. More details concerning these estimates can be found in [2, Theoremโโ5.1].
We derive in the sequel the conditions on the elastic boundaries. Introducing the asymptotic solution of order into (2.6) we get
where the expression of on is given by:
and an analogous expression can be written on . We notice that for the first terms in the above definition are zero.
We complete the problem (5.1), (5.3), (5.5) with the following boundary conditions:
and initial conditions:
The compatibility condition for the problem (5.1), (5.3) with the boundary and initial conditions (5.5), (5.7), (5.8) is given by:
We notice that the boundary conditions for and are of different type from those for and , respectively; so we cannot apply directly the estimates (3.20) which hold for different right-hand sides in the equations but for the same boundary conditions. Moreover, the asymptotic velocity of order is not a divergence free function. To overcome this difficulty, we define new functions and in order to obtain for a problem of the same type as (2.6).
Construction of For this construction we use some ideas from [5]. We first define , satisfying:
Let us consider
with four functions which are determined by imposing some conditions, as we show below. It is obvious that the functions satisfy (5.10) without any conditions on .
We introduce next the following notations:
We notice that from (5.2) we obtain . A relation between the functions and follows from , with the initial conditions (5.10). In a similar way we establish a relation between and . So we proved the possibility of construction of functions which satisfy (5.10) and (5.11).
We define next
and we verify that it satisfies the same boundary conditions as .
Construction of We look for a function satisfying the following problem, for :
Remark 5.1. Since vanishes in , (5.15) represents in fact two separate problems, one in and the second one on . The compatibility condition for each problem is fulfilled due to the construction of .
We can prove that.
Proposition 5.2. The problem (5.15) has at least one solution, with the property
Proof. Due to (5.15), we can consider the problem (5.15) as two separate problems: one in and the second in , each problem with its compatibility condition. So, we can obtain the estimate (5.16) only for . For any we introduce the new variable ; obviously . We define a new function , by . Obvious computations lead to the following problem for :
The right-hand sides of (5.17) contain either correctors (which are of order , with independent on ), or terms of order . Applying a result of [15], Chapter III, page 127 we get:
Expressing next the norm of with respect to the norm of we obtain , that is, . With similar computations we get the same estimate for , which completes the proof.
We notice that from (5.16) and from the regularity of the right-hand sides of (5.17) it follows that .
The function
satisfies the same type of boundary conditions as on . Moreover , due to (4.1). As a consequence of the previous constructions, we obtain for the following problem:
with .
We are now in a position to apply (3.20) in order to establish the error between the exact solution of (2.6) and its approximation given by (4.4).
Theorem 5.3. Let be the asymptotic solution of order given by (4.4) and the exact solution of the physical problem (2.6). Then the following estimates hold:
Proof. Let us obtain the estimate (5.21). From (5.16) and (5.19) it follows that
It remains to estimate the first term of the right-hand side of the previous inequality. For this we first estimate and then we apply (3.20). This estimate is given either by or by , the other terms having greater order than these. Following the ideas of [2, Theoremโโ5.1] we obtain (5.21); the other estimates of (5.21) are proved with the same technique and the proof is achieved.
In order to improve the estimates (5.21) we analyze the order of the leading term of the asymptotic solution. We prove that.
Proposition 5.4. For the leading term of the asymptotic solution one has the following estimates:
Proof. Let us establish (5.23). From (4.48) (for ) and from (4.63) (for ) we obtain for the regular velocity corresponding to
This means that the regular part of the asymptotic velocity in the norm of (5.24) is . It remains to estimate in the same norm the boundary layer functions. For the corrector corresponding to we have
which yields
Finally, for the corrector in we have
that leads to an estimate of the same type as (5.26). Hence, we obtain (5.23) as a consequence of (5.24) and (5.26). More details concerning this type of computations can be found in [5, Theoremโโ5.1].
The previous Proposition allows us to improve the estimates given by Theorem 5.3 in the following sense.
Theorem 5.5. Let be the asymptotic solution of order and the exact solution of the physical problem. Then the following estimates hold:
Proof. Let be a fixed integer and . Then from (5.21) and (5.23) we get
The other estimates of (5.28) are obtained in a similar way, and the proof is achieved.
Acknowledgments
The authors were partially supported by the following grants: SFR MODMAD of the University of Saint-Etienne, ENISE and the School of Mines of Saint-Etienne, the joint French-Russian PICS CNRS grant โMathematical modeling of blood diseasesโ and by the Grant no. 14.740.11.0875 โMultiscale problems: analysis and methodsโ of the Ministry of Education and Research of Russian Federation.
References
G. P. Panasenko and R. Stavre, โAsymptotic analysis of a periodic flow in a thin channel with visco-elastic wall,โ Journal de Mathématiques Pures et Appliquées, vol. 85, no. 4, pp. 558โ579, 2006.
G. P. Panasenko and R. Stavre, โAsymptotic analysis of a non-periodic flow in a thin channel with visco-elastic wall,โ Networks and Heterogeneous Media, vol. 3, no. 3, pp. 651โ673, 2008.
G. P. Panasenko and R. Stavre, โAsymptotic expansion of the solution to the Stokes flow problem in a thin cylindrical elastic tube,โ Applied Analysis. In press.
G. P. Panasenko and R. Stavre, โWell posedness and asymptotic expansion of solution of Stokes equation set in a thin cylindrical elastic tube,โ in Around the Research of Vladimir Maz'ya. II, A. Laptev, Ed., vol. 12 of International Mathematical Series, Springer, New York, NY, USA, 2010.
G. P. Panasenko and R. Stavre, โAsymptotic analysis of the Stokes flow with variable viscosity in a thin elastic channel,โ Networks and Heterogeneous Media, vol. 5, no. 4, pp. 783โ812, 2010.
G. Cardone, G. P. Panasenko, and Y. Sirakov, โAsymptotic analysis and numerical modeling of mass transport in tubular structures,โ Mathematical Models & Methods in Applied Sciences, vol. 20, no. 3, pp. 397โ421, 2010.
G. P. Panasenko, Multi-Scale Modelling for Structures and Composites, Springer, Dordrecht, The Netherlands, 2005.
B. M. Haines, I. S. Aranson, L. Berlyand, and D. A. Karpeev, โEffective viscosity of dilute bacterial suspensions: a two-dimensional model,โ Physical Biology, vol. 5, no. 4, Article ID 046003, pp. 1โ9, 2008.
B. Desjardins, M. J. Esteban, C. Grandmont, and P. le Tallec, โWeak solutions for a fluid-elastic structure interaction model,โ Revista Matemática Complutense, vol. 14, no. 2, pp. 523โ538, 2001.
B. Desjardins and M. J. Esteban, โExistence of weak solutions for the motion of rigid bodies in a viscous fluid,โ Archive for Rational Mechanics and Analysis, vol. 146, no. 1, pp. 59โ71, 1999.
C. Grandmont and Y. Maday, โExistence for an unsteady fluid-structure interaction problem,โ Mathematical Modelling and Numerical Analysis, vol. 34, no. 3, pp. 609โ636, 2000.
D. Serre, โChute libre d'un solide dans un fluide visqueux incompressible. Existence,โ Japan Journal of Applied Mathematics, vol. 4, no. 1, pp. 99โ110, 1987.
R. Temam, Navier-Stokes Equations. Theory and Numerical Analysis, North-Holland, Amsterdam, The Netherlands, 1977.
A. M. Ilyin, Matching of Asymptotic Expansions of Solutions of Boundary Value Problems, Translations of Mathematical Monographs, American Mathematical Society, Washington, DC, USA, 1992.
G. P. Galdi, An Introduction to the Mathematical Theory of the Navier-Stokes Equations, Springer, New York, NY, USA, 1994.