Infiltration in Porous Media with Dynamic Capillary Pressure: Travelling Waves
|
|
- Ella Charles
- 5 years ago
- Views:
Transcription
1 Infiltration in Porous Media with Dynamic Capillary Pressure: Travelling Waves C. Cuesta, C.J. van Duijn, J. Hulshof Abstract We consider a model for non-static groundwater flow where the saturation-pressure relation is extended by a dynamic term. This approach together with a convective term due to gravity, results in a pseudo-parabolic Burgers type equation. We give a rigorous study of global travelling wave solutions, with emphasis on the role played by the dynamic term and the appearance of fronts. Introduction Capillary pressure is an essential characteristic of two-phase flow in porous media. It is the empirical macroscopic description of the pressure differences, due to surface tension, between the phases in the pores of an elementary representation volume. In the standard approach, capillary pressure is expressed as a monotone function of one of the phase saturations: decreasing in terms of the wetting phase saturation and increasing in terms of the non-wetting phase saturation. Such expressions are based on static conditions within an elementary representative volume, see for instance Bear [3, 4] or Bedrikovetsky [5]. Recently, see Gray and Hassanisadeh [], new and more realistic models have been proposed to include non-static conditions as well. In its simplest form this leads to a capillary pressure, which is now a function of the wetting phase saturation and its time derivative. The formulation, with this improved capillary pressure, results in a transport equation containing higher order mixed derivatives. This equation is subject of study in this paper. We confine ourselves to the particular case of unsaturated groundwater flow, where imbibition takes places under influence of gravity. The two phases are water (wetting phase) and air (non-wetting phase). Throughout we assume relative small values of the water saturation, so that regions where the porous medium is fully saturated, as described by Hulshof [] and van Duijn and Peletier [8] for example, do not occur. According to Bear [3, 4] we have for water in a homogeneous and isotropic porous medium the momentum balance equation and the mass balance q = K(S)( p w + ρge z ) (Darcy s law), (.) φ t (ρs) + div(ρq) = 0. (.) CWI, P.O Box 94079, 090 GB Amsterdam, The Netherlands, carlota@cwi.nl. CWI, P.O Box 94079, 090 GB Amsterdam, The Netherlands, hansd@cwi.nl. Mathematical Department of Leiden University, Niels Bohrweg, 333 CA Leiden, The Netherlands, hulshof@wi.leidenuniv.nl.
2 Here q denotes volumetric water flux, S water saturation, K(S) hydraulic conductivity, p w water pressure, ρ water density, φ porosity and g gravity constant. Further e z is the unit vector in positive z-direction, i.e. against the direction of gravity. To solve equation (.) and (.) an additional relation between p w and S is needed. This relation is based on the assumption that the air pressure, p a, is constant in space and time. This assumption is common practice in unsaturated groundwater flow and valid for small air viscosity. Under static conditions, see Bear [3, 4] and Bedrikovetsky [5], one uses p a p w = p c (S), (.3) where p c denotes the monotone capillary pressure. Gray and Hassanisadeh [] proposed to use Under non-static conditions, p a p w = p c (S) φl(s) t S, (.4) to capture the difference between drainage and imbibition. Here L(S) is a nonlinear damping term. Equations (.), (.) and (.4) can be combined into a single equation for S, which reads φ t (ρs) = div{ρk(s)ρge z + ρk(s) ( p c (S) + φl(s) t S)}. (.5) Assuming now a one-dimensional flow in vertical z-direction, with ρ and φ constant, and applying a straightforward scaling, equation (.5) reduces to t S = z {K(S) + K(S) z ( p c (S) + L(S) t S)}. (.6) In a previous work, see [] and [], G.I. Barenblatt proposed a different model to describe the non-static situation. He modified (.) and (.3) by replacing S by S + τ t S (τ > 0), into p c (S) and K(S). The resulting equation then is of the form t S = zz {Φ (S + τ t S)} + z {Φ (S + τ t S)}, where Φ and Φ are nonlinear positive functions and degenerate at S + τ t S = 0. This equation admits a splitting into two equations which decouple the space and time derivatives. Namely, { w = S + τst τ(φ (w) zz + Φ (w) z ) + w = S. Thus it can be studied as an ODE problem in Banach spaces. Moreover maximum principles are applicable, see [7]. There is a similar splitting for (.6), which reads { w = St p c (S) (K(S)w z ) z + w = K(S) z p c (S) The appearance of K(S) as a coefficient in the w-equation, however does not allow the maximum principle to be applied as in [7].
3 We note here that if the Barenblatt ansatz (S S + τ S) is applied only in (.3), then (.6) would result with L(S) = τ dp c φ ds, implying that L(S) becomes unbounded as S 0. However, experiments carried out by Smiles et al [4], see also Hassanizadeh [0] for an overview, show that L(S) vanishes as S 0. To investigate the role of the nonlinear terms (i.e. K(S), p c (S) and L(S)), we replace them by power-law relations. Note that this is consistent with the assumption of small water saturation. We approximate K(S) = S α (α > ), p c (S) = + S ( > 0), L(S) = εs γ (γ > 0), where ε > 0 is introduced as a parameter to investigate the consequence of the third order mixed term. The parameter ranges are chosen to capture the relevant physical properties of unsaturated flow. In particular, we want K and L to be non-negative, with L(0) = K(0) = K (0) = 0, and p c (0 + ) =. Using these power law relations, equation (.6) reduces to t S = z {S α + S α z S + εs α z (S γ t S)}. (.7) The static capillary pressure relation (.3), would have resulted in the convectiondiffusion equation t S = z {S α + S α z S}. (.8) It is well-known, see e.g. Gilding [9], that this equation has finite speed of propagation if and only if δ 0 D(S) S ds < for some δ > 0, with D(S) = Sα. (.9) This requires α >. Because occurrence of fronts has our special interest we analyse equation (.7) in the parameter range, γ > 0 and < α. (.0) In this paper we analyse travelling waves solutions of (.7). They are conjectured to describe the large time behaviour of solutions resulting from a certain class of initial conditions. The stability properties of travelling waves are subject of a separated future study. Thus we consider subject to the boundary conditions S(z, t) = f(η) with η = z + ct, (.) f( ) = A > 0, f( ) = δ 0 (δ < A). (.) 3
4 Hence the fluid moves downwards whenever the wave speed c is positive. For f we obtain the equation cf = {f α + f α f + cεf α (f γ f ) } on R, (.3) Integration and application of boundary conditions give the equation c(f δ) = f α δ α + f α f + cεf α (f γ f ), (.4) and the second order boundary value problem c(f δ) = f α δ α + f α f + cεf α (f γ f ) on R (T W ) f( ) = δ, f(+ ) = A, where c is given by c = Aα δ α A δ, (.5) which is the Rankine-Hugoniot wave speed if we interpret (.7) as a regularization of the hyperbolic equation t S = x S α. Note that c A α as δ 0, c αa α as δ A, with dc > 0 for 0 < δ < A. dδ In Section we show existence of travelling waves for fixed positive values of ε and δ. They are unique up to translations in η. This analysis also shows an oscillatory, but non-periodic, behaviour of the profiles. Here the value of ε is crucial; for ε sufficiently small (depending on α,, γ, δ and A) we obtain monotone profiles, this is shown in Section 3. In Section 4 we study the limit case ε 0, while δ > 0 is kept fixed. Using essentially monotonicity for small ε, we obtain convergence to travelling waves of equation (.8), i.e. the standard model based on the static pressure saturation relation (.3). In Section 5 we analyse existence of front solutions to Problem T W with δ = 0. It turns out that there are two relevant ranges of powers α,, γ for which fronts appear. In the range > α γ there exists a family of solutions which degenerate at a finite value of η. When = α γ existence of fronts is shown provided ε. Uniqueness does not hold. Nevertheless we have discerned 4A in each of the previous (α ) cases a unique (up to translations in η) exceptional profile, which is the limit profile to (T W ) when letting δ 0. This is shown in Section 6. In the other cases, = α γ, ε, and < α γ, bounded travelling 4A wave solutions do not exist. These results (α ) correspond to the formal asymptotic analysis made in [3]. It is worthy to observe that the limits ε 0 and δ 0 do not commute. We can always take the limit ε 0 followed by δ 0. However, the converse order is only possible if α γ. Remark. Fixing δ (0, A), the existence is demonstrated for α > and, γ > 0. All other results require in addition < α. 4
5 Existence and Uniqueness The main result of this section is: Theorem. Let α >,, γ, ε > 0 and 0 < δ < A. Further, let c be given by (.5). Then there exists a C solution of Problem T W, unique up to translations in η. Proof. We transform equation (.4) into a planar system and apply a phase plane analysis. First we set u = f, which gives cu α = (cδ δ α )u α + + γ u( with boundary conditions ) u + εc u( ) = δ, u( ) = A. + γ u on R (.) Next we put this equation in the Liénard phase-plane, by considering u and v := εc u u, as independent variables. This results in the system ) εu = F (u, v) = (v + u c (P ) v = G(u) = + cu α (cδ δ α )u α. A solution f of (T W ) is an integral curve of (P ) connecting the equilibria p = (δ, δ ) and p + = (A, A ). The phase plane, with the isoclines Γ u = {(u, v) : F (u, v) = 0} and Γ v = {(u, v) : G(u) = 0}, is drawn in Figure. v δ S A S u p + Γ u p - S 4 S 3 Figure : Phase plane for (P ) 5
6 The matrix of the linearised system is X γ εc εc X α γ (c( α)x + (cx X α )α) 0 where X = δ corresponds to the equilibrium p, and X = A to p +. The eigenvalues in p and p + are, expressed in X, λ = a ± a + b with a = c αxα and b = εc X+γ+ X α+γ εc. (.) It is clear that a > 0 for positive X, and b has the sign of c X α α. Introducing its primitive i(x) = cx X α, (.3) we note that A and δ are related by i(a) = i(δ) > 0, so that b > 0 in X = δ and b < 0 in X = A. Therefore at p the two eigenvalues are real and of opposite sign, whence p is a saddle point. The point p + is a sink and, 0 a + b < a implies two real eigenvalues, λ λ < 0, and a + b < 0 implies two complex eigenvalues, with negative real part. (.4) To prove Theorem. we have to show that an orbit leaving p connects to p +. Because p is a saddle point there exist locally two orbits containing solutions (u(η), v(η)) of (P ), satisfying lim (u(η), v(η)) = p. η Let C be the orbit for which u, v > 0. Inspection of Figure shows that the other orbit will never reach p +. The only possibility giving existence of a travelling wave is for C to end up in p +. The corresponding solutions will satisfy (T W ), and uniqueness up to translations in η will hold. To show that C connects to p + we consider the sets, see also Figure, S = {(u, v) R : δ γ+ < u < A γ+, v > u }, S = {(u, v) R : u > A γ+, v > u }, S 3 = {(u, v) R : u > A γ+, v < u }, S 4 = {(u, v) R : δ γ+ < u < A γ+, v < u }. Note that the boundaries of S, S, S 3, S 4 are the isoclines of system (P ). Furthermore C enters S from p. We have the following two possibilities. Lemma. The orbit C rotates around p +, going from S i to S i+ mod 4, or it enters p + from S or S 3. 6
7 Proof. Let (u(η), v(η)) be a solution of (P ) parametrising C. Further, let η max be the maximum η-value for which the solution can be extended to (, η max ). Near points where u (η) 0, we can express any solution of (P ) locally as v = v(u), satisfying dv du = v = ε G(u) u F (u, v). (.5) Below we exhaust all possibilities.. Suppose (u(η), v(η)) S for all 0 < η < η max. Then two cases need to be considered: (i): u ū and v v, with ū, v <. This implies η max = + and (ū, v) = p +. (ii): u ū and v. Since F (u, v) becomes unbounded as v, while G(u) remains bounded, (.5) directly gives a contradiction. Thus either (u(η), v(η)) p + as η or C crosses u = A at η = η (0, η max ) and enters S.. Suppose (u(η), v(η)) remains in S for all η (η, η max ). Again two cases are possible. (i): u ū and v v with ū < and, of course, v >. Then η max = + and (ū, v) is an equilibrium point, which is impossible because there are no equilibria to the right of p +. (ii): u + and v v, with v 0. Then the v-equation implies v contradicting v 0. Thus, there exists η (η, η max ) such that (u(η), v(η)) crosses Γ u and enters S Suppose (u(η), v(η)) S 3 for all η (η, η max ). Completely analogous to., we find that either (u(η), v(η)) p + or C crosses A at some η 3 (η 0, η max ) and enters S Suppose (u(η), v(η)) S 4 for all η (η 3, η max ). Again we distinguish (i): u ū and v v. As before η max = + and consequently (ū, v) = p. Thus C is a homoclinic orbit. Next we consider the domain D enclosed by C. Since its boundary is smooth, except at p where it is Lipschitz, we may apply Gauss theorem. Using the fact that C is an orbit of (P ) we get ( F (u, v) 0 = C ( ) F (u, v)dv G(u)du = div ε D contradicting ( ) div F (u, v), G(u) ε ε ), G(u), = +γ+ u γ+ < 0 for all u > δ γ+. cε (ii): C crosses u = δ γ+ below p at (δ, ṽ). Consider the closed curve C C, where C is the straight line segment parametrised by { u(s) = δ γ+, v(s) = s s [ṽ, δ ]. 7
8 As before we call D the interior of C C, and apply Gauss theorem. This gives ( ) ( ) F (u, v) 0 > div, G(u) = F (u, v)dv G(u)du = D ε C C ε ( ) δ F (u, v)dv G(u)du = ε ε F (δγ+, s)ds > 0, C a contradiction. Hence there exists η 4 (η 3, η max ) at which (u(η), v(η)) crosses Γ u and enters S. To complete the proof of the theorem we use Lemma. and the Poincaré-Bendixon theorem. This leave us with two possibilities. Either (u(η), v(η)) p + as η or C approaches a periodic orbit. The latter is impossible by the same argument as in 4(i), now applied to the periodic orbit instead of C. The C regularity for f follows from a bootstrap argument. 3 Monotonicity In this section we derive a sufficient condition for the travelling wave solution to have a monotone profile. This condition is related to the value of ε. Therefore we write C = C ε and P = P ε whenever appropriate. The eigenvalues at p + are (see expressions (.)) λ, = εc A +γ+ ± ṽ (εc) A + 4(c Aα α), (+γ+) A α+γ εc with corresponding eigenvectors ( ± ) (u,, v, ) =, + 4(cA A α α)εca α+γ+. ( + γ) A + Thus for all ε such that ε < ε := 4c(A α α c) Aα γ (ε > 0), (3.) λ, λ are real and strictly negative, which is a necessary condition for the travelling wave profile to be monotone. Henceforth we suppose being in this situation. With λ < λ < 0 we call (u, v ) the slow eigenvector and (u, v ) the fast eigenvector at p +. By standard local analysis, see e.g. [6], there exist exactly two orbits entering p + tangent to the (u, v )-direction: one along (u, v ), the other along ( u, v ). The connecting orbit goes around p + at most a finite number of times. Lemma 3. Let (.0) be satisfied. For δ > 0 fixed, there exists 0 < ε < ε, such that for every 0 < ε < ε the travelling wave obtained in Theorem. is strictly increasing on R. 8
9 Lemma 3. is a direct consequence of Proposition 3. below, in which we construct an invariant region which contains C and which itself is contained in S. More specifically, for fixed δ (0, A) and µ (0, ), let S δ µ denote the set enclosed by the curves u = δ, v = u and v = g µ (u) := ( µu ( µ)a ). We will show in Proposition 3. that for ε sufficiently small (i.e. ε < ε δ µ), Sµ δ is invariant for Problem P ε or, equivalently, C ε Sµ δ. Proposition 3. will also be helpful in the study of the limits ε 0 and δ 0. In particular we want to bound ε δ µ away from zero as δ 0. It will appear that this is only possible if + γ + α 0. Proposition 3. For any fixed δ in (0, A) and µ (0, ) there exists ε δ µ (0, ε ), such that for every ε (0, ε δ µ ), C ε Sµ δ. Further, if + γ + α 0, and Proof. lim δ 0 εδ µ = 4µ( µ) lim δ 0 ε = µ( µ) 4(α ) A γ, lim δ 0 Sδ µ is invariant for Problem T W with δ = 0. Observe that the eigenvectors at p + satisfy (u, v ) (, ( + γ)a + ), and (u, v ) (, 0) as ε 0, where (, ) is the tangent vector at p ()A + + to Γ u. The invariant region Sµ δ is below the horizontal line v = A, and only contains orbits entering p + along the slow eigenvector (u, v ). Observe that g µ (u) > u for all 0 < u < A, and g µ (A ) = A. Obviously the vector field is pointing inwards at boundary points of Sµ δ on u = δ and Γ u. It remains to examine the vector field on v = g µ (u). We clearly must have dv du = εc + γ G(u) g µ (u) + u for S δ µ to be invariant. This is equivalent to where H δ (u) = c g µ (u) for all u (δ, A ), (3.) εh δ (u) µ( µ) (3.3) + cu α (cδ δ α )u α (u A ) Note that H δ (δ ) = 0 and, by L Hôspital s rule, H δ (A ) = lim H δ(u) = c αaα c u A u +. (3.4) A +γ+ α (α ) A +γ+α as δ 0. 9
10 For u > 0 fixed, lim H δ (u) behaves as u +γ+ α near 0, which suggests to write δ 0 H δ (u) as where H δ (u) = c u +γ+ α h δ (u) (3.5) h δ (u) = u α + cu (cδ δ α ). u ( A u ) Observe that for every u (0, A ) and every δ > 0, recalling (.5), So that h 0 (u) h δ (u) = (cδ δα ) u (c A α ) > 0. (3.6) u ( A u ) h δ (u) < h 0 (u) = u α + A α u u ( A u α = A α A α u ) A u which is increasing for u (0, A ). Here we used < α from (.0). Setting (α )Aα M := sup h 0 (u) =, u (0,A ) an upper bound for H δ (u) is given by H δ (u) < M c A+γ+ α if + γ + α 0,, (3.7) and H δ (u) < M c δ+γ+ α if + γ + α 0. Thus a sufficient condition for (3.3) to hold is µ( µ) ε < ε δ c(α ) A γ if + γ + α 0 µ := µ( µ) δ α γ if + γ + α 0. c(α )A α (3.8) This completes the proof of the first statement. The statements about the δ 0 limit follow immediately from (3.8) with α γ, (3.) and (.5). 4 The ε 0 limit case Let δ (0, A) be fixed and (.0) be satisfied. In this section we examine the behaviour of the connecting orbit C ε and that of the corresponding travelling wave f = f ε as ε 0. For ε < ε we denote C ε by As a first convergence result we have v = ϕ ε (u), δ γ+ u A γ+. (4.) Proposition 4. ϕ ε (u) u uniformly on [δ, A ] as ε 0. 0
11 (a) monotone profile: ε = ε (b) oscillatory profile: ε = Figure : Travelling wave solutions for different values of ε, where f(0) = A+δ, α = 4 3, = 4, γ =, A =, δ = 6 Proof. Proposition 3. implies u < ϕε (u) < g µ (u) (4.) for all u (δ, A ) and for all ε (0, ε δ µ ). Since g µ(u) u as µ, the result is immediate. For the travelling waves f ε we have Theorem 4. Translate f ε so that f ε (0) = δ+a for all ε > 0. Then f ε f C (R) as ε 0, uniformly on R, where f satisfies Problem T W with ε = 0. Proof. First we employ the scaling η = ετ, so that in the τ-variable (P ε ) reads ( P ε ) { u = F (u, v), v = εg(u). (4.3) Unlike (P 0 ) the limit system ( P 0 ), is well defined. The one-dimensional manifold of critical points M 0 = {F (v, u) = 0} = {v = u } is invariant and normally hyperbolic in the sense of geometric singular perturbation theory, see [5], because for ( P 0 ) the only pure imaginary eigenvalue is zero, and has a one-dimensional eigenspace tangential to M 0. Let K be a neighbourhood of {(u, v) M 0 : δ u A }, and choose 0 < δ < δ < A < A such that {(u, v) M 0 : δ u A } K. By Fenichel s invariant manifold theorem [5], there exists, for given k N, a number ε 0 > 0 and a function h C k ([δ, A ] [0, ε 0 ]) with h(u, 0) = u, such that for every 0 < ε < ε0 M ε = {(u, v) K : v = h(u, ε), δ u A }
12 is locally invariant. The manifold M ε is not uniquely determined. However between u = δ and u = A it must coincide with the connecting orbit v = ϕ ε (u), because this is the only orbit which remains close to {v = u : δ < u < A }. Using (4.) and the v-equation in (P ε ), we note that u ε = fε satisfies, u = G(u) ϕ ε(u), (4.4) and connects the two zeros of G. Since h C k and h(u, ε) = ϕ ε (u) for δ u A, we have as a result of Proposition 4., that ϕ ε(u) G(u) ϕ ε(u) u uniformly on [δ, A ] and thus G(u) u uniformly on [δ, A ] as ε 0. In this limit the differential equation (4.4) is identical to equation (.) with ε = 0. Using the fact that u ε (0) is fixed for all ε > 0, standard arguments imply that u ε converges uniformly on R to the corresponding solution of the limit equation. 5 δ = 0 system In this section we consider the limit case δ = 0 directly. Thus we study the system { ) εu (Pε 0 = F 0 (u, v) = (v + u, c ) v = G 0 (u) = + cu α, where c = A α, and we look for orbits connecting u = 0 to u = A. The critical point corresponding to the latter now has real eigenvalues, see also (3.), for 0 < ε ε = 4(α ) A α γ. The phase plane, see Figure, clearly implies, that the desired orbit has to be originated from the segment {(u, v) : u = 0, v 0} where the equations are singular. Since we are interested in (Pε 0 ) as limit of (Pε δ ), and in particular of a possible limit orbit of the connecting orbit C, we expect such a limit orbit, if it exists, to behave as v d u, 0 < d, as u 0. Thus a convenient new dependent variable is Z = u q v, where q, for later purposes, is not fixed yet. Whenever u 0, Z satisfies the equation u dz du = qz + εc + γ u+q+ α c u α Z + u q. (5.) Below we investigate the solvability of (5.) for 0 < u < A. The analysis and results critically depend on the value of the parameters α, and γ. In particular the value of α + γ + plays a crucial role, which is to be expected considering the results of the formal analysis by Hulshof and King [3]. With q appropriately chosen, we consider the cases:
13 5. > α γ. Here we take q = +γ+ α, and set W = u. Then (5.) becomes ( + γ + α) dz dw = Z W + εc(c W Z + α +γ+ α ). (5.) We look for solutions of (5.) with Z > as W 0 (i.e. u > 0 as u 0). In Figure 3 we sketch the (W, Z)-phase plane. Equation (5.) and the phase plane imply that Z Z 0 {0, } as W 0, where orbits with Z 0 = 0 have v = o( u ) while orbits with Z0 = have v u Proposition 5. For > α γ there is a unique orbit C 0 with u 0 and v u as η decreases. This orbit reaches (u, v) = (0, ) at some finite η-value, implying the existence of a travelling wave with a front. The local behavior of the front is determined by the relation. f c f + α as f 0. Proof. We first prove existence. Choose W 0 small and denote the solution of (5.) with Z = ξ in W = W 0 by Z = Z(W, ξ). Let S + = {ξ (, 0) : Z(W, ξ) 0 as W 0} S = {ξ (, 0) : W (0, W 0 ) Z(W, ξ) as W W > 0} S 0 = {ξ (, 0) : Z(W, ξ) as W 0}. By standard arguments we have for W 0 sufficiently small that (0, ) = S S 0 S + and S and S + are nonempty and open. Hence S 0 is nonempty, which gives existence. We observe that for such an orbit, see (5.), Z + aw as W 0 where a = εc. (5.3) Next we prove uniqueness. Suppose there are two solutions Z = Z (W ) and Z = Z (W ) with Z as W 0. Since Z + Z + = ( Z + ) (Z Z ), where Z lies between Z and Z, we have for Y = Z Z > 0, say, where b =, and +γ+ α b = W 0, contradicting (5.3). dy dw b Y W Y b as W 0, W. Here we used (5.3). Hence Y as εc (+γ+ α) 3
14 Expressing W and Z in terms of f we observe that (5.3) implies the behaviour f c f + α as f 0. Since + α < we find that f reaches zero at some finite η-value; i.e. the travelling wave has a finite front. Z 0 W - Figure 3: Phase plane realted to (5.), with q =. Note that we did not put any restriction on ε. Thus the conclusion about the behaviour as f 0, provided > α γ, is valid for any value of ε 0. This ε-uniformity is lost in the next case. Indeed, the behaviour of the travelling wave near the front when it vanishes corresponds to criterion (.9) in the case where the capillary damping is absent. 5. = α γ. Now we take q =, and set W = u α. Then (5.) becomes (α ) dz dw = ( Z + εc(c W ) ) (Z + ) W. (5.4) Again we look for solutions satisfying Z > and Z Z 0 as W 0. It follows that Z 0 ± = ± ) ( 4εc, implying ε 4c. With reference to Figure 4 we have Proposition 5. (i) If ε <, there exists a family of orbits satisfying v 4c Z 0 + u as u 0, and a unique orbit, denoted by C 0, which satisfies v Z0 u as u 0. All cases give travelling waves with finite fronts. In particular the orbit C 0 implies ( f ( ) ) εc 4εc f + α, as f 0. 4
15 (ii) If ε = 4c there exists a family of orbits having v u as u 0. In particular there is a unique orbit, again denoted by C 0, satisfying u v as u 0. The orbit C 0 implies again a traveling wave with a finite front, such that f c f + α, as f 0. (iii) If ε > 4c, there is no orbit with v > u and u 0. Remark 5. (i) Comparing (3.) and (3.8) we observe that ε 0 = ε = 4(α )c < 4c, implying that, depending on the ε-value, monotone and oscillatory waves with finite fronts occur. (ii) Since now + α = γ < 0, finite front waves have f as f 0. Proof. (i) It is immediate from the phase-plane in Figure 4 that (0, Z 0 + ) is a source, and (0, Z0 ) is a saddle, with one( unique orbit Z = ) Z(W ) leaving in the direction W > 0. It behaves as Z Z0 4εc ( ) as W 0. (ii) Now Z 0 = at local analysis shows that (0, Z 0) is a saddle-node, with a unique orbit in the direction of W > 0. This orbit does not cross the isocline, and behaves as Z Z 0 W. (iii) Now the segment {(W, Z) : W = 0, Z (, 0)} is disconnected from the isocline and hence no connecting orbit exists. 5.3 < α γ. Here we choose q = α γ () > and W = uq, yielding the equation (α ) α γ ) (α γ ) dz dw = (α γ ) Z W + εc(c W (W + Z)W (5.5) Proposition 5.3 There exists no orbits with v > u and u 0. Proof. Suppose such an orbit exists. Then we would have (α γ ) dz dw ) ((α γ )Z + εc, Z > 0, W Z as W 0. Since (α γ )Z + εc Z gives a contradiction. is negative and bounded away from zero, this 5
16 Z 0 + Z0 W Z 0 W Z 0 - Z0 - - ε < - ε = 4 c - 4 c Z 0 W - ε > - 4c Figure 4: q = 6 δ 0 limit, = α γ phase-plane for different ε-values. For the purpose of this section we denote the connecting orbit of (P δ ε ) by v = ϕ δ (u) and, in the cases for which C 0 from Proposition 5. and 5. exists, we call its graph v = ϕ 0 (u). Lemma 6. There exists δ > 0 such that, with c = c(δ) = Aα δ α F u (δ) = c(δ)g δ (u) = c(δ) A δ, ( + c(δ)u α c(δ)a A α u α is decreasing in 0 δ < δ for any fixed 0 < u ( A ). Proof. Since df u dδ = dc ) ( + cu α + (A α ca)u α, dδ and dc > 0, we need to show that the term between brackets is negative for small dδ δ. At δ = 0 it becomes + A α u α (u A) < 0, ) 6
17 for all 0 < u < ( A ). Proposition 6. For α,, γ such that > α γ, or = α γ and ε (0, ) fixed, translate f 4c δ such that f δ (0) = A for all δ (0, δ ). Then f δ f C (R) uniformly on R. Hence f satisfies Problem T W with δ = 0. Proof. It will be sufficient to show that ϕ δ (u) ϕ 0 (u) locally uniformly. By Lemma 6. we have for any 0 < δ δ u < ϕδ (u) < ϕ δ (u) < ϕ 0 (u) for 0 < δ < δ < δ and u ( δ, ( A ) ). Also εc(δ) + γ G δ (u) F δ (u, v) εaα + γ uniformly on [ δ, ( A ) ]. Therefore ϕ δ (u) ϕ(u) ϕ 0 (u) G 0 (u) F 0 (u, v) as δ 0 where v = ϕ 0 (u) is a solution of (Pε 0 ). The reasoning above holds for every 0 < δ δ, which implies that ϕ 0 (u) exists for all u (0, ( ) A ). In view of Section 5 and ϕ ϕ 0 this implies that ϕ 0 (u) = ϕ 0 (u). Using that f δ (0) = A is fixed for all 0 < δ < δ, standard arguments imply that f δ converges uniformly on (, A ) to the corresponding solution of the limit equation. Existence of global travelling waves and uniqueness of the initial value problem for all 0 δ < δ, implies the uniform convergence on R. 7 Concluding Remarks In this paper we present a model for unsaturated groundwater flow which includes an expression for the non-static phase pressure difference, see (.4). Replacing the nonlinearities in the transport equation by power-law expressions we arrive at (.7). We study travelling wave solutions representing moisture profiles moving downwards due to gravity. For positive initial saturation (δ > 0) we demonstrate existence and uniqueness (up to translations). Small values of the damping coefficient ε result in monotone saturation profiles. Large values of ε result in profiles which exhibit oscillatory behaviour near the injection saturation A. When initially no moisture is present (δ = 0) existence of bounded travelling waves depends critically on the exponents of the power-law expressions. This is related to the occurrence of finite fronts in the moisture profiles: i.e. descending planes (in the direction of gravity) below which the water saturation remains zero. Related to equation (.7) we have shown the following. 7
18 If > α γ, then travelling wave solutions with fronts exist for all ε > 0. In other words, for S(z, t) = f(η), with η = z + ct, there exists η 0 R such that f(η) = 0 for all η η 0. Moreover near η = η 0 the profile satisfies f Aα f + α. This corresponds to the front behaviour of solutions of the convection diffusion equation under static conditions (equation (.8)); i.e. ε and γ are absent in this asymptotic expression. If = α γ we obtain a similar result provided the damping coefficient ε is 4A (α ) sufficiently small: i.e. ε. For larger values of ε, no waves exist satisfying f( ) = 0. Finally, if < α γ, again no such waves exist, regardless the value of ε > 0. We also investigate the limit ε 0 (for δ > 0, fixed) and δ 0 (for ε > 0 fixed). In particular the latter provides a uniqueness criterion for the degenerate case when δ = 0. We also note that the limits ε 0 and δ 0 do not commute: ε 0 followed by δ 0 is always possible, while δ 0 followed by ε 0 is only possible when α γ. References [] G.I. Barenblatt, J. García-Azorero, A. De Pablo, J.L. Vázquez, Mathematical Model of the Non-Equilibrium Water-Oil Displacement in Porous Strata, Applicable Analysis, 65: 9-45, 997. [] G.I. Barenblatt, V.M. Entov and V.M. Ryzhik, Theory of fluid flow through natural rocks, Kluwer, Dordrecht-Boston-London, 990. [3] J. Bear, Dynamics of fluids in porous media, Dover (nd Publisher), New York, 988. [4] J. Bear, Hydraulics of groundwater, McGraw-Hill: Series in Water Resources and Enviromental Engineering, New York, 979. [5] P. Bedrikovetsky, Mathematical theory of oil and gas recovery, Kluwer, Dordrecht-Boston-London 993. [6] E. Coddington and N. Levinson, Theory of ordinary differential equations, Mc- Graw Hill Book Company, Inc., New York, 955. [7] E. Dibenedetto and M. Pierre, On the maximum principle for pseudoparabolic equations, Indiana Univ. Math. J., 30: 8-854, 98. [8] C.J. van Duijn and L. A. Peletier, Nonstationary filtration in partially saturated porous media, Archs. ratio. Mech. Analysis, 78: 73-98, 98. [9] B.H. Gilding, The occurrence of interfaces in nonlinear diffusion-advection processes, Arch. Rational Mech. Anal. 00, 43-63,
19 [0] S.M. Hassanizadeh, Dynamic Effects in the capillary pressure saturation relationship, Proceedings of 4th international conference on Civil Engineering, Volume 4: Water Resources and envirom. eng., 4-49, Sharif U. of Technology, Tehran, Iran, May 997. [] S.M. Hassanizadeh and W.G. Gray, Thermodynamic basis of capillary pressure in porous media, Water Resources Reseach, 9: , 993. [] J. Hulshof, An elliptic-parabolic free boundary problem: continuity of the free boundary, Proc. Roy. Soc., Edinb. 06A: , 987. [3] J. Hulshof and J.R. King, Analysis of a Darcy flow model with a dynamic pressure saturation relation, SIAM, Journal of Applied Math. 59(), , 998. [4] D.E. Smiles, G. Vachaud, and M. Vauclin, A test of the uniqueness of the soil moisture characteristic during transient, non-hysteretic flow of water in rigid soil, Soil Sci. Soc. Amer. Proc., 35: , 97. [5] J. Smoller, Shock waves and reaction-diffusion equations, Springer, nd Edition,
Infiltration in porous media with dynamic capillary pressure: travelling waves
Euro. Jnl of Applied Mathematics (000), vol. 11, pp. 381 397. c 000 Cambridge University Press Printed in the United Kingdom 381 Infiltration in porous media with dynamic capillary pressure: travelling
More informationIntroduction Multiphase ows in porous media occur in various practical situations including oil recovery and unsaturated groundwater ow, see for insta
Analysis of a Darcy ow model with a dynamic pressure saturation relation Josephus Hulshof Mathematical Department of the Leiden University Niels Bohrweg, 333 CA Leiden, The Netherlands E-mail: hulshof@wi.leidenuniv.nl
More informationTravelling wave solutions for the Richards equation incorporating non-equilibrium effects in the capillarity pressure
UHasselt Computational Mathematics Preprint eries Travelling wave solutions for the Richards equation incorporating non-equilibrium effects in the capillarity pressure Hans van Duijn, Koondanibha Mitra
More informationTravelling waves. Chapter 8. 1 Introduction
Chapter 8 Travelling waves 1 Introduction One of the cornerstones in the study of both linear and nonlinear PDEs is the wave propagation. A wave is a recognizable signal which is transferred from one part
More informationA new class of entropy solutions of the Buckley-Leverett equation
A new class of entropy solutions of the Buckley-Leverett equation C.J. van Duijn, L.A. Peletier & I.S. Pop Abstract We discuss an extension of the Buckley-Leverett (BL) equation describing twophase flow
More informationTRAVELING WAVES IN 2D REACTIVE BOUSSINESQ SYSTEMS WITH NO-SLIP BOUNDARY CONDITIONS
TRAVELING WAVES IN 2D REACTIVE BOUSSINESQ SYSTEMS WITH NO-SLIP BOUNDARY CONDITIONS PETER CONSTANTIN, MARTA LEWICKA AND LENYA RYZHIK Abstract. We consider systems of reactive Boussinesq equations in two
More informationAuthor(s) Huang, Feimin; Matsumura, Akitaka; Citation Osaka Journal of Mathematics. 41(1)
Title On the stability of contact Navier-Stokes equations with discont free b Authors Huang, Feimin; Matsumura, Akitaka; Citation Osaka Journal of Mathematics. 4 Issue 4-3 Date Text Version publisher URL
More informationTHE CRITICAL WAVE SPEED FOR THE FISHER-KOLMOGOROV-PETROWSKII-PISCOUNOV EQUATION WITH CUT-OFF
THE CRITICAL WAVE SPEED FOR THE FISHER-KOLMOGOROV-PETROWSKII-PISCOUNOV EQUATION WITH CUT-OFF FREDDY DUMORTIER, NIKOLA POPOVIĆ, AND TASSO J. KAPER Abstract. The Fisher-Kolmogorov-Petrowskii-Piscounov (FKPP)
More informationTraveling stripes in the Klausmeier model of vegetation pattern formation
Traveling stripes in the Klausmeier model of vegetation pattern formation Paul Carter Arjen Doelman Abstract The Klausmeier equation is a widely studied reaction-diffusion-advection model of vegetation
More informationUNIQUENESS OF SELF-SIMILAR VERY SINGULAR SOLUTION FOR NON-NEWTONIAN POLYTROPIC FILTRATION EQUATIONS WITH GRADIENT ABSORPTION
Electronic Journal of Differential Equations, Vol. 2015 2015), No. 83, pp. 1 9. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu UNIQUENESS OF SELF-SIMILAR
More information2 The second case, in which Problem (P 1 ) reduces to the \one-phase" problem (P 2 ) 8 >< >: u t = u xx + uu x t > 0, x < (t) ; u((t); t) = q t > 0 ;
1 ON A FREE BOUNDARY PROBLEM ARISING IN DETONATION THEORY: CONVERGENCE TO TRAVELLING WAVES 1. INTRODUCTION. by M.Bertsch Dipartimento di Matematica Universita di Torino Via Principe Amedeo 8 10123 Torino,
More informationPiecewise Smooth Solutions to the Burgers-Hilbert Equation
Piecewise Smooth Solutions to the Burgers-Hilbert Equation Alberto Bressan and Tianyou Zhang Department of Mathematics, Penn State University, University Park, Pa 68, USA e-mails: bressan@mathpsuedu, zhang
More informationLecture No 2 Degenerate Diffusion Free boundary problems
Lecture No 2 Degenerate Diffusion Free boundary problems Columbia University IAS summer program June, 2009 Outline We will discuss non-linear parabolic equations of slow diffusion. Our model is the porous
More informationHomogenization and error estimates of free boundary velocities in periodic media
Homogenization and error estimates of free boundary velocities in periodic media Inwon C. Kim October 7, 2011 Abstract In this note I describe a recent result ([14]-[15]) on homogenization and error estimates
More informationSHOCK WAVES FOR RADIATIVE HYPERBOLIC ELLIPTIC SYSTEMS
SHOCK WAVES FOR RADIATIVE HYPERBOLIC ELLIPTIC SYSTEMS CORRADO LATTANZIO, CORRADO MASCIA, AND DENIS SERRE Abstract. The present paper deals with the following hyperbolic elliptic coupled system, modelling
More informationdynamics of f luids in porous media
dynamics of f luids in porous media Jacob Bear Department of Civil Engineering Technion Israel Institute of Technology, Haifa DOVER PUBLICATIONS, INC. New York Contents Preface xvii CHAPTER 1 Introduction
More informationIntroduction LECTURE 1
LECTURE 1 Introduction The source of all great mathematics is the special case, the concrete example. It is frequent in mathematics that every instance of a concept of seemingly great generality is in
More informationOn the Cauchy Problems for Polymer Flooding with Gravitation
On the Cauchy Problems for Polymer Flooding with Gravitation Wen Shen Mathematics Department, Penn State University Abstract We study two systems of conservation laws for polymer flooding in secondary
More informationOn the Cauchy Problems for Polymer Flooding with Gravitation
On the Cauchy Problems for Polymer Flooding with Gravitation Wen Shen Mathematics Department, Penn State University. Email: wxs27@psu.edu November 5, 2015 Abstract We study two systems of conservation
More informationThe stability of travelling fronts for general scalar viscous balance law
J. Math. Anal. Appl. 35 25) 698 711 www.elsevier.com/locate/jmaa The stability of travelling fronts for general scalar viscous balance law Yaping Wu, Xiuxia Xing Department of Mathematics, Capital Normal
More informationA REMARK ON THE GLOBAL DYNAMICS OF COMPETITIVE SYSTEMS ON ORDERED BANACH SPACES
PROCEEDINGS OF THE AMERICAN MATHEMATICAL SOCIETY Volume 00, Number 0, Pages 000 000 S 0002-9939(XX)0000-0 A REMARK ON THE GLOBAL DYNAMICS OF COMPETITIVE SYSTEMS ON ORDERED BANACH SPACES KING-YEUNG LAM
More informationTopic 6: Projected Dynamical Systems
Topic 6: Projected Dynamical Systems John F. Smith Memorial Professor and Director Virtual Center for Supernetworks Isenberg School of Management University of Massachusetts Amherst, Massachusetts 01003
More informationStability of Feedback Solutions for Infinite Horizon Noncooperative Differential Games
Stability of Feedback Solutions for Infinite Horizon Noncooperative Differential Games Alberto Bressan ) and Khai T. Nguyen ) *) Department of Mathematics, Penn State University **) Department of Mathematics,
More informationON THE SOLVABILITY OF A NONLINEAR PSEUDOPARABOLIC PROBLEM
Indian J. Pure Appl. Math., 44(3): 343-354, June 2013 c Indian National Science Academy ON THE SOLVABILITY OF A NONLINEAR PSEUDOPARABOLIC PROBLEM S. Mesloub and T. Hayat Mathematics Department, College
More informationẋ = f(x, y), ẏ = g(x, y), (x, y) D, can only have periodic solutions if (f,g) changes sign in D or if (f,g)=0in D.
4 Periodic Solutions We have shown that in the case of an autonomous equation the periodic solutions correspond with closed orbits in phase-space. Autonomous two-dimensional systems with phase-space R
More informationSome aspects of vanishing properties of solutions to nonlinear elliptic equations
RIMS Kôkyûroku, 2014, pp. 1 9 Some aspects of vanishing properties of solutions to nonlinear elliptic equations By Seppo Granlund and Niko Marola Abstract We discuss some aspects of vanishing properties
More informationMEAN CURVATURE FLOW OF ENTIRE GRAPHS EVOLVING AWAY FROM THE HEAT FLOW
MEAN CURVATURE FLOW OF ENTIRE GRAPHS EVOLVING AWAY FROM THE HEAT FLOW GREGORY DRUGAN AND XUAN HIEN NGUYEN Abstract. We present two initial graphs over the entire R n, n 2 for which the mean curvature flow
More informationVarious lecture notes for
Various lecture notes for 18311. R. R. Rosales (MIT, Math. Dept., 2-337) April 12, 2013 Abstract Notes, both complete and/or incomplete, for MIT s 18.311 (Principles of Applied Mathematics). These notes
More informationNon-linear Wave Propagation and Non-Equilibrium Thermodynamics - Part 3
Non-linear Wave Propagation and Non-Equilibrium Thermodynamics - Part 3 Tommaso Ruggeri Department of Mathematics and Research Center of Applied Mathematics University of Bologna January 21, 2017 ommaso
More informationPropagating terraces and the dynamics of front-like solutions of reaction-diffusion equations on R
Propagating terraces and the dynamics of front-like solutions of reaction-diffusion equations on R P. Poláčik School of Mathematics, University of Minnesota Minneapolis, MN 55455 Abstract We consider semilinear
More informationNoether s Theorem. 4.1 Ignorable Coordinates
4 Noether s Theorem 4.1 Ignorable Coordinates A central recurring theme in mathematical physics is the connection between symmetries and conservation laws, in particular the connection between the symmetries
More informationNONCLASSICAL SHOCK WAVES OF CONSERVATION LAWS: FLUX FUNCTION HAVING TWO INFLECTION POINTS
Electronic Journal of Differential Equations, Vol. 2006(2006), No. 149, pp. 1 18. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu (login: ftp) NONCLASSICAL
More informationNon-degeneracy of perturbed solutions of semilinear partial differential equations
Non-degeneracy of perturbed solutions of semilinear partial differential equations Robert Magnus, Olivier Moschetta Abstract The equation u + F(V (εx, u = 0 is considered in R n. For small ε > 0 it is
More informationPart 1 Introduction Degenerate Diffusion and Free-boundaries
Part 1 Introduction Degenerate Diffusion and Free-boundaries Columbia University De Giorgi Center - Pisa June 2012 Introduction We will discuss, in these lectures, certain geometric and analytical aspects
More informationSynchronization, Chaos, and the Dynamics of Coupled Oscillators. Supplemental 1. Winter Zachary Adams Undergraduate in Mathematics and Biology
Synchronization, Chaos, and the Dynamics of Coupled Oscillators Supplemental 1 Winter 2017 Zachary Adams Undergraduate in Mathematics and Biology Outline: The shift map is discussed, and a rigorous proof
More informationDYNAMIC BIFURCATION THEORY OF RAYLEIGH-BÉNARD CONVECTION WITH INFINITE PRANDTL NUMBER
DYNAMIC BIFURCATION THEORY OF RAYLEIGH-BÉNARD CONVECTION WITH INFINITE PRANDTL NUMBER JUNGHO PARK Abstract. We study in this paper the bifurcation and stability of the solutions of the Rayleigh-Bénard
More informationTRANSPORT IN POROUS MEDIA
1 TRANSPORT IN POROUS MEDIA G. ALLAIRE CMAP, Ecole Polytechnique 1. Introduction 2. Main result in an unbounded domain 3. Asymptotic expansions with drift 4. Two-scale convergence with drift 5. The case
More informationψ ae is equal to the height of the capillary rise in the soil. Ranges from about 10mm for gravel to 1.5m for silt to several meters for clay.
Contents 1 Infiltration 1 1a Hydrologic soil horizons...................... 1 1b Infiltration Process......................... 2 1c Measurement............................ 2 1d Richard s Equation.........................
More informationTwo Phase Transport in Porous Media
Two Phase Transport in Porous Media Lloyd Bridge Iain Moyles Brian Wetton Mathematics Department University of British Columbia www.math.ubc.ca/ wetton CRM CAMP Seminar, October 19, 2011 Overview Two phase
More informationIntegrodifferential Hyperbolic Equations and its Application for 2-D Rotational Fluid Flows
Integrodifferential Hyperbolic Equations and its Application for 2-D Rotational Fluid Flows Alexander Chesnokov Lavrentyev Institute of Hydrodynamics Novosibirsk, Russia chesnokov@hydro.nsc.ru July 14,
More informationThuong Nguyen. SADCO Internal Review Metting
Asymptotic behavior of singularly perturbed control system: non-periodic setting Thuong Nguyen (Joint work with A. Siconolfi) SADCO Internal Review Metting Rome, Nov 10-12, 2014 Thuong Nguyen (Roma Sapienza)
More informationRionero s critical perturbations method via weighted energy for stability of convective motions in anisotropic porous media
1 Rionero s critical perturbations method via weighted energy for stability of convective motions in anisotropic porous media F. Capone,a, M. Gentile,a, A. A. Hill,b a Department of Mathematics and Applications
More informationINVARIANT CURVES AND FOCAL POINTS IN A LYNESS ITERATIVE PROCESS
International Journal of Bifurcation and Chaos, Vol. 13, No. 7 (2003) 1841 1852 c World Scientific Publishing Company INVARIANT CURVES AND FOCAL POINTS IN A LYNESS ITERATIVE PROCESS L. GARDINI and G. I.
More informationarxiv: v2 [math.ap] 1 Jul 2011
A Godunov-type method for the shallow water equations with discontinuous topography in the resonant regime arxiv:1105.3074v2 [math.ap] 1 Jul 2011 Abstract Philippe G. efloch 1 and Mai Duc Thanh 2 1 aboratoire
More informationAutonomous systems. Ordinary differential equations which do not contain the independent variable explicitly are said to be autonomous.
Autonomous equations Autonomous systems Ordinary differential equations which do not contain the independent variable explicitly are said to be autonomous. i f i(x 1, x 2,..., x n ) for i 1,..., n As you
More informationEVANESCENT SOLUTIONS FOR LINEAR ORDINARY DIFFERENTIAL EQUATIONS
EVANESCENT SOLUTIONS FOR LINEAR ORDINARY DIFFERENTIAL EQUATIONS Cezar Avramescu Abstract The problem of existence of the solutions for ordinary differential equations vanishing at ± is considered. AMS
More informationUniversität des Saarlandes. Fachrichtung 6.1 Mathematik
Universität des Saarlandes U N I V E R S I T A S S A R A V I E N I S S Fachrichtung 6.1 Mathematik Preprint Nr. 155 A posteriori error estimates for stationary slow flows of power-law fluids Michael Bildhauer,
More informationLINEAR FLOW IN POROUS MEDIA WITH DOUBLE PERIODICITY
PORTUGALIAE MATHEMATICA Vol. 56 Fasc. 2 1999 LINEAR FLOW IN POROUS MEDIA WITH DOUBLE PERIODICITY R. Bunoiu and J. Saint Jean Paulin Abstract: We study the classical steady Stokes equations with homogeneous
More informationNon-degeneracy of perturbed solutions of semilinear partial differential equations
Non-degeneracy of perturbed solutions of semilinear partial differential equations Robert Magnus, Olivier Moschetta Abstract The equation u + FV εx, u = 0 is considered in R n. For small ε > 0 it is shown
More informationPERIODIC SOLUTIONS WITH NONCONSTANT SIGN IN ABEL EQUATIONS OF THE SECOND KIND
PERIODIC SOLUTIONS WITH NONCONSTANT SIGN IN ABEL EQUATIONS OF THE SECOND KIND JOSEP M. OLM, XAVIER ROS-OTON, AND TERE M. SEARA Abstract. The study of periodic solutions with constant sign in the Abel equation
More informationEffects of Viscous Dissipation on Unsteady Free Convection in a Fluid past a Vertical Plate Immersed in a Porous Medium
Transport in Porous Media (2006) 64: 1 14 Springer 2006 DOI 10.1007/s11242-005-1126-6 Effects of Viscous Dissipation on Unsteady Free Convection in a Fluid past a Vertical Plate Immersed in a Porous Medium
More informationUniform estimates for Stokes equations in domains with small holes and applications in homogenization problems
Uniform estimates for Stokes equations in domains with small holes and applications in homogenization problems Yong Lu Abstract We consider the Dirichlet problem for the Stokes equations in a domain with
More informationS chauder Theory. x 2. = log( x 1 + x 2 ) + 1 ( x 1 + x 2 ) 2. ( 5) x 1 + x 2 x 1 + x 2. 2 = 2 x 1. x 1 x 2. 1 x 1.
Sep. 1 9 Intuitively, the solution u to the Poisson equation S chauder Theory u = f 1 should have better regularity than the right hand side f. In particular one expects u to be twice more differentiable
More information1.4 The Jacobian of a map
1.4 The Jacobian of a map Derivative of a differentiable map Let F : M n N m be a differentiable map between two C 1 manifolds. Given a point p M we define the derivative of F at p by df p df (p) : T p
More information7 Planar systems of linear ODE
7 Planar systems of linear ODE Here I restrict my attention to a very special class of autonomous ODE: linear ODE with constant coefficients This is arguably the only class of ODE for which explicit solution
More informationRadon measure solutions for scalar. conservation laws. A. Terracina. A.Terracina La Sapienza, Università di Roma 06/09/2017
Radon measure A.Terracina La Sapienza, Università di Roma 06/09/2017 Collaboration Michiel Bertsch Flavia Smarrazzo Alberto Tesei Introduction Consider the following Cauchy problem { ut + ϕ(u) x = 0 in
More informationUnsaturated Flow (brief lecture)
Physical Hydrogeology Unsaturated Flow (brief lecture) Why study the unsaturated zone? Evapotranspiration Infiltration Toxic Waste Leak Irrigation UNSATURATAED ZONE Aquifer Important to: Agriculture (most
More informationLipschitz shadowing implies structural stability
Lipschitz shadowing implies structural stability Sergei Yu. Pilyugin Sergei B. Tihomirov Abstract We show that the Lipschitz shadowing property of a diffeomorphism is equivalent to structural stability.
More informationn v molecules will pass per unit time through the area from left to
3 iscosity and Heat Conduction in Gas Dynamics Equations of One-Dimensional Gas Flow The dissipative processes - viscosity (internal friction) and heat conduction - are connected with existence of molecular
More informationMildly degenerate Kirchhoff equations with weak dissipation: global existence and time decay
arxiv:93.273v [math.ap] 6 Mar 29 Mildly degenerate Kirchhoff equations with weak dissipation: global existence and time decay Marina Ghisi Università degli Studi di Pisa Dipartimento di Matematica Leonida
More informationDELFT UNIVERSITY OF TECHNOLOGY
DELF UNIVERSIY OF ECHNOLOGY REPOR -2 A model for a viscous preflush prior to gelation in a porous medium F.J. Vermolen, P.L.J. Zitha and J. Bruining ISSN 389-652 Reports of the Department of Applied Mathematical
More informationResearch Scholar, PACIFIC University, Udaipur, Rajasthan, India. Department of Mathematics, UGC Visiting Professor, HNGU Patan, Gujarat, India
2018 IJSRSET Volume 4 Issue 4 Print ISSN: 2395-1990 Online ISSN : 2394-4099 Themed Section : Engineering and Technology Adomian Decomposition Study of Unidimensional Flow in Unsaturated Porous Media Khatri
More informationIntroduction to Partial Differential Equations
Introduction to Partial Differential Equations Partial differential equations arise in a number of physical problems, such as fluid flow, heat transfer, solid mechanics and biological processes. These
More informationTakens embedding theorem for infinite-dimensional dynamical systems
Takens embedding theorem for infinite-dimensional dynamical systems James C. Robinson Mathematics Institute, University of Warwick, Coventry, CV4 7AL, U.K. E-mail: jcr@maths.warwick.ac.uk Abstract. Takens
More informationResearch interests: Nonlinear partial differential equations and their applications
Research interests: Nonlinear partial differential equations and their applications Degenerate diffusion, higher order equations, reaction diffusion equations, heat flow of harmonic maps, pseudo-parabolic
More informationParameter Dependent Quasi-Linear Parabolic Equations
CADERNOS DE MATEMÁTICA 4, 39 33 October (23) ARTIGO NÚMERO SMA#79 Parameter Dependent Quasi-Linear Parabolic Equations Cláudia Buttarello Gentile Departamento de Matemática, Universidade Federal de São
More informationThe Riemann problem. The Riemann problem Rarefaction waves and shock waves
The Riemann problem Rarefaction waves and shock waves 1. An illuminating example A Heaviside function as initial datum Solving the Riemann problem for the Hopf equation consists in describing the solutions
More informationHeat kernels of some Schrödinger operators
Heat kernels of some Schrödinger operators Alexander Grigor yan Tsinghua University 28 September 2016 Consider an elliptic Schrödinger operator H = Δ + Φ, where Δ = n 2 i=1 is the Laplace operator in R
More informationAn introduction to Mathematical Theory of Control
An introduction to Mathematical Theory of Control Vasile Staicu University of Aveiro UNICA, May 2018 Vasile Staicu (University of Aveiro) An introduction to Mathematical Theory of Control UNICA, May 2018
More informationNonlinear Autonomous Systems of Differential
Chapter 4 Nonlinear Autonomous Systems of Differential Equations 4.0 The Phase Plane: Linear Systems 4.0.1 Introduction Consider a system of the form x = A(x), (4.0.1) where A is independent of t. Such
More informationFlow and Transport. c(s, t)s ds,
Flow and Transport 1. The Transport Equation We shall describe the transport of a dissolved chemical by water that is traveling with uniform velocity ν through a long thin tube G with uniform cross section
More informationProblem set 7 Math 207A, Fall 2011 Solutions
Problem set 7 Math 207A, Fall 2011 s 1. Classify the equilibrium (x, y) = (0, 0) of the system x t = x, y t = y + x 2. Is the equilibrium hyperbolic? Find an equation for the trajectories in (x, y)- phase
More information1 Lyapunov theory of stability
M.Kawski, APM 581 Diff Equns Intro to Lyapunov theory. November 15, 29 1 1 Lyapunov theory of stability Introduction. Lyapunov s second (or direct) method provides tools for studying (asymptotic) stability
More informationWe denote the space of distributions on Ω by D ( Ω) 2.
Sep. 1 0, 008 Distributions Distributions are generalized functions. Some familiarity with the theory of distributions helps understanding of various function spaces which play important roles in the study
More informationApplications of the compensated compactness method on hyperbolic conservation systems
Applications of the compensated compactness method on hyperbolic conservation systems Yunguang Lu Department of Mathematics National University of Colombia e-mail:ylu@unal.edu.co ALAMMI 2009 In this talk,
More informationProblem: A class of dynamical systems characterized by a fast divergence of the orbits. A paradigmatic example: the Arnold cat.
À È Ê ÇÄÁ Ë ËÌ ÅË Problem: A class of dynamical systems characterized by a fast divergence of the orbits A paradigmatic example: the Arnold cat. The closure of a homoclinic orbit. The shadowing lemma.
More informationFirst Boundary Value Problem for the Diffusion. Equation. I. Iterated Logarithm Test for the Boundary. Regularity and Solvability
First Boundary Value Problem for the Diffusion Equation I. Iterated Logarithm Test for the Boundary Regularity and Solvability Ugur G. Abdulla Max-Planck Institute for Mathematics in the Sciences, Leipzig
More informationTurning points and traveling waves in FitzHugh-Nagumo type equations
Turning points and traveling waves in FitzHugh-Nagumo type equations Weishi Liu and Erik Van Vleck Department of Mathematics University of Kansas, Lawrence, KS 66045 E-mail: wliu@math.ku.edu, evanvleck@math.ku.edu
More informationMath 273 (51) - Final
Name: Id #: Math 273 (5) - Final Autumn Quarter 26 Thursday, December 8, 26-6: to 8: Instructions: Prob. Points Score possible 25 2 25 3 25 TOTAL 75 Read each problem carefully. Write legibly. Show all
More informationSingular perturbation methods for slow fast dynamics
Nonlinear Dyn (2007) 50:747 753 DOI 10.1007/s11071-007-9236-z ORIGINAL ARTICLE Singular perturbation methods for slow fast dynamics Ferdinand Verhulst Received: 2 February 2006 / Accepted: 4 April 2006
More informationExponential Stability of the Traveling Fronts for a Pseudo-Para. Pseudo-Parabolic Fisher-KPP Equation
Exponential Stability of the Traveling Fronts for a Pseudo-Parabolic Fisher-KPP Equation Based on joint work with Xueli Bai (Center for PDE, East China Normal Univ.) and Yang Cao (Dalian University of
More informationA Ginzburg-Landau approach to dislocations. Marcello Ponsiglione Sapienza Università di Roma
Marcello Ponsiglione Sapienza Università di Roma Description of a dislocation line A deformed crystal C can be described by The displacement function u : C R 3. The strain function β = u. A dislocation
More informationFundamentals of Fluid Dynamics: Ideal Flow Theory & Basic Aerodynamics
Fundamentals of Fluid Dynamics: Ideal Flow Theory & Basic Aerodynamics Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI (after: D.J. ACHESON s Elementary Fluid Dynamics ) bluebox.ippt.pan.pl/
More informationA symmetry analysis of Richard s equation describing flow in porous media. Ron Wiltshire
BULLETIN OF THE GREEK MATHEMATICAL SOCIETY Volume 51, 005 93 103) A symmetry analysis of Richard s equation describing flow in porous media Ron Wiltshire Abstract A summary of the classical Lie method
More informationEquivariant self-similar wave maps from Minkowski spacetime into 3-sphere
Equivariant self-similar wave maps from Minkowski spacetime into 3-sphere arxiv:math-ph/99126v1 17 Oct 1999 Piotr Bizoń Institute of Physics, Jagellonian University, Kraków, Poland March 26, 28 Abstract
More informationRegularity for Poisson Equation
Regularity for Poisson Equation OcMountain Daylight Time. 4, 20 Intuitively, the solution u to the Poisson equation u= f () should have better regularity than the right hand side f. In particular one expects
More informationGlobal unbounded solutions of the Fujita equation in the intermediate range
Global unbounded solutions of the Fujita equation in the intermediate range Peter Poláčik School of Mathematics, University of Minnesota, Minneapolis, MN 55455, USA Eiji Yanagida Department of Mathematics,
More informationUniform estimates for Stokes equations in domains with small holes and applications in homogenization problems
Uniform estimates for Stokes equations in domains with small holes and applications in homogenization problems Yong Lu Abstract We consider the Dirichlet problem for the Stokes equations in a domain with
More informationComputation of homoclinic and heteroclinic orbits for flows
Computation of homoclinic and heteroclinic orbits for flows Jean-Philippe Lessard BCAM BCAM Mini-symposium on Computational Math February 1st, 2011 Rigorous Computations Connecting Orbits Compute a set
More informationClassical solutions for the quasi-stationary Stefan problem with surface tension
Classical solutions for the quasi-stationary Stefan problem with surface tension Joachim Escher, Gieri Simonett We show that the quasi-stationary two-phase Stefan problem with surface tension has a unique
More informationPROBLEM OF CRACK UNDER QUASIBRITTLE FRACTURE V.A. KOVTUNENKO
PROBLEM OF CRACK UNDER QUASIBRITTLE FRACTURE V.A. KOVTUNENKO Overview: 1. Motivation 1.1. Evolutionary problem of crack propagation 1.2. Stationary problem of crack equilibrium 1.3. Interaction (contact+cohesion)
More informationHyperbolic Systems of Conservation Laws. I - Basic Concepts
Hyperbolic Systems of Conservation Laws I - Basic Concepts Alberto Bressan Mathematics Department, Penn State University Alberto Bressan (Penn State) Hyperbolic Systems of Conservation Laws 1 / 27 The
More informationFundamentals of Fluid Dynamics: Elementary Viscous Flow
Fundamentals of Fluid Dynamics: Elementary Viscous Flow Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI bluebox.ippt.pan.pl/ tzielins/ Institute of Fundamental Technological Research
More informationInvariant Sets for non Classical Reaction-Diffusion Systems
Global Journal of Pure and Applied Mathematics. ISSN 0973-1768 Volume 1, Number 6 016, pp. 5105 5117 Research India Publications http://www.ripublication.com/gjpam.htm Invariant Sets for non Classical
More informationFlux Invariants for Shape
Flux Invariants for Shape avel Dimitrov James N. Damon Kaleem Siddiqi School of Computer Science Department of Mathematics McGill University University of North Carolina, Chapel Hill Montreal, Canada Chapel
More informationTHE NEUMANN PROBLEM FOR THE -LAPLACIAN AND THE MONGE-KANTOROVICH MASS TRANSFER PROBLEM
THE NEUMANN PROBLEM FOR THE -LAPLACIAN AND THE MONGE-KANTOROVICH MASS TRANSFER PROBLEM J. GARCÍA-AZORERO, J. J. MANFREDI, I. PERAL AND J. D. ROSSI Abstract. We consider the natural Neumann boundary condition
More informationAPPROXIMATE SIMILARITY SOLUTION TO A NONLINEAR DIFFUSION EQUATION WITH SPHERICAL SYMMETRY
INTERNATIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING Volume 9, Number, Pages 4 c 22 Institute for Scientific Computing and Information APPROXIMATE SIMILARITY SOLUTION TO A NONLINEAR DIFFUSION EQUATION
More informationAA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS
AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 1 / 31 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS Linearization and Characteristic Relations 1 / 31 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS
More informationCombustion fronts in porous media with two layers
Combustion fronts in porous media with two layers layer 1 layer 2 Steve Schecter North Carolina State University (USA) Jesus Carlos da Mota Universidade Federal de Goiás (Brazil) 1 Overview Subject: Propagation
More informationTHE ROSENZWEIG-MACARTHUR PREDATOR-PREY MODEL
THE ROSENZWEIG-MACARTHUR PREDATOR-PREY MODEL HAL L. SMITH* SCHOOL OF MATHEMATICAL AND STATISTICAL SCIENCES ARIZONA STATE UNIVERSITY TEMPE, AZ, USA 8587 Abstract. This is intended as lecture notes for nd
More information