Solvability of the p-Adic Analogue of Navier–Stokes Equation via the Wavelet Theory

P-adic numbers serve as the simplest ultrametric model for the tree-like structures arising in various physical and biological phenomena. Recently p-adic dynamical equations started to be applied to geophysics, to model propagation of fluids (oil, water, and oil-in-water and water-in-oil emulsion) in capillary networks in porous random media. In particular, a p-adic analog of the Navier–Stokes equation was derived starting with a system of differential equations respecting the hierarchic structure of a capillary tree. In this paper, using the Schauder fixed point theorem together with the wavelet functions, we extend the study of the solvability of a p-adic field analog of the Navier–Stokes equation derived from a system of hierarchic equations for fluid flow in a capillary network in porous medium. This equation describes propagation of fluid’s flow through Geo-conduits, consisting of the mixture of fractures (as well as fracture’s corridors) and capillary networks, detected by seismic as joint wave/mass conducts. Furthermore, applying the Adomian decomposition method we formulate the solution of the p-adic analog of the Navier–Stokes equation in term of series in general form. This solution may help researchers to come closer and find more facts, taking into consideration the scaling, hierarchies, and formal derivations, imprinted from the analogous aspects of the real world phenomena.


Introduction
The last decades have witnessed great use of Fourier and more generally wavelet analysis over the p-adic fields, and its various physical applications in physics, biology and cognitive science, and recently in geophysics. The keyword of these applications is "hierarchy". These applications are based on representation of hierarchies by tree-like geometry. Hierarchy is also a natural attribute of ultrametric spaces which mathematically can be represented as a duality between ultrametric (non-Archimedean) spaces and trees of balls in these spaces where mathematical tools such as integral and series are frequently used (see also [1,2]). Thus, ultrametric (non-Archimedean) spaces play the  In 2017, Oleschko et al. [38] focused the p-adic dynamics described by fractional differential operators (Vladimirov operators) starting with discrete dynamics based on hierarchically-structured interactions between the fluids' volumes concentrated at different levels of the percolation tree and coming to the multiscale universal topology of the percolating nets. They presented a system of dynamical equations reflecting the tree structure of a capillary network in porous media and then derived the following nonlinear p-adic pseudo-differential equation for fluid's velocity u(t.x) along capillaries: x)Du(t, x) − θD 2 u(t, x) + G(t, x), x ∈ Q p , t, u ∈ R, t ≥ 0, where θ is the viscosity parameter with the initial condition u(0, x) = ϕ(x). This equation can be considered as the p-adic analogue of the Navier-Stokes equation. We stress that this is just an analog of the Navier-Stokes equation. The tree-like configuration space differs crucially from the real space of hydrodynamics and this difference is reflected in dynamical equations. In particular, u is a real scalar and not a real vector with three coordinates as in the Navier-Stokes equation in the Euclidean space, see Section 3 for details. Alternative nonlinear term in p-adic analogue of the Navier-Stokes equation was considered in paper of Kozyrev [39]. Moreover, the mathematical theory of such equations has not yet been developed. Very recently, Khrennikov and Kochubei [43] investigated the local solvability of Equation (1) using the von Wahl's theorem for the case G(t, x) as the source term vanishes.
In the present paper, inspired by [43], our attention will be turned exclusively to study the solvability of Equation (1) using the technique of wavelet basis and a well-known fixed point theorem. Theory of p-adic wavelets was initiated by Kozyrev [11] and it found numerous applications (see, e.g., [2,12,13,35]), including modeling of fluid propagation in capillary networks in random disordered medium [36][37][38]. In Section 4, under certain conditions we study the solvability of an infinite system which is derived from (1) (with wavelet expansion of solution) in the sequence space c 0 . In Section 5, we employ a numerical method, the so-called Adomian decomposition method (ADM), to formulate the solution of Equation (1) represented by series.
This paper demonstrates that, for fluids' propagation through capillary networks in porous disordered media, p-adic linear models developed and investigated in our previous works [35,36] can be successfully generalized (at the mathematical level of rigorousness) to nonlinear phenomena.

Geophysics: From Fractal to Tree-Like Models
The starting point of our research was fractal/multifractal modeling in geophysics [44,45]. The detailed presentation on such an approach can be found in Section 1.1 of our paper [36]. Here we briefly point to the most important moments. The fractal/multifractal scaling features of capillary networks was studied (both theoretically and experimentally) since the early 1980s and 1990s [46,47]: for invasion percolation, diffusion-limited aggregation (DLA), anti-DLA processes [48].
Later fractal modeling of fluids' flows in porous random media [49], including transport through tree-like networks and diffusion on fractals [50] was widely used in oil recovery studies [51]. This modeling was supported by theoretical and experimental studies demonstrating that fluid's flow through tree-like networks is faster [52,53].
Stanley and Meakin [46] have discussed the important thermodynamics aspects of multifractality in physics and chemistry, founding the formal analogy among the probability distribution function Z(q) and partition function Z(β). Therefore, the analogy between the Legendre transform f (α) and entropy H, as well as between the function α and energy E was found. These analogies are the key points for physics of fractal capillaries patterns treelike morphology. The difference in the medium heterogeneity can be quantified by several multifractal indicators (for instance, the degree of the graph symmetry or strength of singularity).
Thus, fractal/multifractal studies led to an understanding of the importance of tree-like structures in mathematical modeling of fluids' flows through capillary networks in porous disordered media.

Navier-Stokes Equation on Tree-Like Configuration Space and Its Generalizations
Here we briefly repeat the basics of the physical model leading to the p-adic analog of the Navier-Stokes equation, see [38] for details. The main motivation for derivation of this equation is application to modeling propagation of fluid's flow through Geo-conduits, consisting of the mixture of fractures (as well as fracture's corridors) and capillary networks, detected by seismic as joint wave/mass conducts (see Figures 1 and 2).
In such modeling [35,36], the tree-like structure of capillaries in real porous disordered media is represented by trees endowed with the root distance-ultrametric spaces. The rock environment of capillary networks is ignored. We explore only geometry of the network of capillaries. So, instead of the Euclidean configuration space, we use the tree-like configuration space. The rock environment is encoded in the coefficients of the dynamical equations describing fluids' propagation through capillary networks in porous media. Tree-like geometry and ultrametric spaces give the proper mathematical model for such networks.
In our model, the "spatial-variable" x belongs to an ultrametric space denoted by symbol X and the time variable t is real. We operate with functions f (t, x) depending on real and ultrametric variables. Here x is the "pore network coordinate", each pathway of pore capillaries is encoded by a point x of the ultrametric space (or in the tree-like representation-by a branch of the tree, see [36] for details). Time is usual real time. Thus, by assigning the ultrametric coordinate x to a system (e.g., oil or water droplet) we know in which pathway composed of capillaries this system is located. The ultrametric model provides a fuzzy description of system's location in a pore network. The simplest trees are homogeneous trees, the p-adic trees. Furthermore, in this paper, we concentrate our study on such configuration spaces. The general ultrametric case was considered in [36]. Of course, real capillary structures are described by non-homogeneous trees. However, the general case is essentially more complex mathematically. Here we were able to study (at the mathematical level of rigorousness) only linear pseudo-differential equations for fluid's propagation through capillary networks in porous media. Now, we turn to derivation of the p-adic analog of the Navier-Stokes equation, for simplicity we consider the case p = 2. Consider the homogeneous tree with two branches leaving each vertex. The root of the tree is denoted by I 0 ; the levels of the tree are enumerated, n = 0, 1, 2, ...., where n = 0 corresponds to the root I 0 . Vertexes at the nth level of the tree are enumerated as I n,j , j = 1, ..., 2 n . We are interested in capillaries connecting successive vertexes, i.e., from I n−1,j to the corresponding two vertexes at level n. They can be labeled in the same way as corresponding vertices, i.e., E 1,1 = I 0 I 1,1 and E 1,2 = I 0 I 1,2 , the capillaries going from root I 0 to vertexes I 1 and I 2 . There are 2n edges connecting vertexes of the (n − 1)-th level with vertices of the n-th level: E n,1 = I n−1,1 I n,1 , E n,2 = I n−1,1 I n,2 , E n,3 = I n−1,2 I n,3 , E n,4 = I n−1,2 I n,4 , ..., or E n,2j−1 = In − 1, jI n,2j−1 , E n,2j = I n−1,j I n,2j .
In the present model the diameters and lengths of edges are coupled by scaling 1/2. Thus at each vertex I n,j fluid's flow is split into two capillaries that are two times shorter and thinner than the capillary incoming to this vertex. Denote by u n,j (t) the average velocity of fluid along capillary E n,j . Thus, in our model we are not interested in fluid's velocity in each point of a capillary, but only in its average with respect to capillary's volume. The velocity is a real number. Its sign encodes the (average) direction of the velocity along the capillary, towards and backwards with respect to the root of the tree. In the mathematical model, instead of the discrete variables n, j we can operate with the continuous 2-adic variable x ∈ Q 2 and velocity u(t, x), where t, u ∈ R, x ∈ Q 2 . In [38], we derived nonlinear pseudo-differential dynamical Equation (1) for velocity u(t, x). Its form is analogous to the form of the standard Navier-Stokes equation. However, it is the scalar-function equation, since fluid's propagation is only along capillary's axis.
Derivation of the p-adic analog of the Navier-Stokes equation opens the door to consideration of a bunch of interesting problems on dynamics of fluids in capillary networks. One of such problems is coupling of the dynamics of fluid with electromagnetic field, derivation of the p-adic (capillary network) analog of the system of magneto-fluid dynamics equations (cf. [54]). The main difficulty is derivation of "restriction" of the Maxwell equations onto the capillary network and then their coupling with the p-adic Navier-Stokes equation (cf. [54]). Such a theory can find applications in geophysics (cf. [55]).
Finally, we remark that the tree-like structure of capillary networks is not only very common in nature, but recently they started to be artificially manufactured-to speed up fluid's propagation (see [56] and references herein). This applied research is theoretically justified by the works of Shou et al. [52,53] on the tree-like basis of acceleration of fluid's flows (see also [57] on p-adic modeling). Such industrial applications stimulate p-adic (and more general ultrametric) modeling of fluids' propagation in tree-like capillary networks.

Mathematical Preliminaries
In this section, we recall some auxiliary facts concerned with p-adic fields and wavelet theory. In view of the Ostrovski theorem (see [6], Ch. I, § 1.1), there exists, in some sense, only two "universes" of equal status: the real universe and the p-adic one. The real "universe" is structured by the field of real numbers R, which is introduced by the completion of the field Q of rational numbers with relevance to the usual Euclidean norm, and the p-adic "universe" is based on the field Q p of p-adic numbers, which is given as the completion of the field Q with respect to the p-adic norm | · | p . This norm is defined as below. | · | p : |0| p = 0; if an arbitrary rational number x = 0 is represented as ∈ Z and m, n are not divisible by p then |x| p = p −γ . This norm satisfies the following properties: |x + y| p ≤ max{|x| p , |y| p }, for every x, y ∈ Q p , and when |x| p = |y| p , we have |x + y| p = max{|x| p , |y| p }.
The condition (iii) as the strong triangle inequality makes the norm | · | p non-Archimedean and hence the space (Q p , | · | p ) is an ultrametric space.
We shall systematically utilize the notation and results from [6]. Denote by N, Z, C the sets of positive integers, integers, and complex numbers, respectively.
Any p-adic number x ∈ Q p , x = 0, is represented in the canonical form as follows where γ = γ(x) ∈ Z, and x k = 0, 1, . . . , p − 1, x 0 = 0, k = 0, 1, . . . The series converges in the p-adic norm | · | p to p −γ , that is, |x| p = p −γ . Hence, the absolute value | · | p takes the discrete set of nonzero values p γ , for γ ∈ Z. The fractional part of a p-adic number x ∈ Q p given by (2) is defined as The additive character χ p of the field Q p is given by The topology equipped with | · | p in Q p are known by as balls and spheres of radius p γ with center at a, respectively. It is worth mentioning that any point of the ball is its center, besides, any two balls in Q p are either disjoint or one is included in the other. Furthermore, all balls and spheres are simultaneously open and closed sets in Q p . For the certain case, the unit ball Z p = B 0 (0) is the ring of p-adic integers consisting of the elements represented by the sum of p mutually disjoint balls.
The topological group (Q p , +) is locally compact commutative and thus there is a additive Haar measure dx, which is positive and invariant under the translation, i.e., d(x + a) = dx, a ∈ Q p . This measure is unique by normalizing dx so that Regarding with the additive normalized character where Ω(t) is the characteristic function of the segment [0, 1] ⊂ R. A complex-valued function f in Q p is said to be a locally constant function if for any x ∈ Q p , there exists an integer l(x) ∈ Z such that f (x + y) = f (x), for every y ∈ B l(x) . The largest of these numbers, l = l( f ), is called the parameter of constancy of the function f . We denote the space of locally constant functions on Q p by E (Q p ). Indicate by D(Q p ) the space of Bruhat-Schwartz test functions, i.e., the subspace of E (Q p ) including compactly supported functions. Moreover, denote by D (Q p ) the set of all linear functionals on D(Q p ) (see also ([6], VI.3)).
The Fourier transform of test function ϕ ∈ D(Q p ) is given by the formulâ which is evidently a Hilbert space with the inner product This guarantees a linear isomorphism taking D(Q p ) onto D(Q p ). It can be uniquely extended to a linear isomorphism of L 2 (Q p ). Moreover, the Plancherel equality holds

p-Adic Wavelet Theory
Throughout this section, we gather some facts related with the theory of p-adic wavelets which is widely employed in so many applications. It is now hard to find an area of engineering where wavelets are not applied. In 1910, Haar [58] initially presented the wavelet basis by an orthonormal basis in L 2 (Q p ) including dyadic translations and dilations of a single function; since then various generalizations of it have been revealed in several results. It is interesting to know that it took almost a century to create another wavelet function whose shifts and dilations would bring an orthogonal basis. The intensive progression in wavelet theory initiated only in the 1990s. At that moment Meyer [59] and Mallat [60,61] improved a scheme of structure for wavelet functions based on the concept of multiresolution analysis (MRA); see, for instance, [62], Ch. 5, [63], § 2.1.
Regarding with the theory of p-adic wavelets, it has a short background in comparison with that in the real status. In 2002, Kozyrev [11] found a compactly supported p-adic wavelet basis, similar to the real Haar basis, for L 2 (Q p ). Kozyrev's wavelet functions have the following structure: where χ p and Ω are the standard additive character of Q p and characteristic function of [0, 1], respectively. This wavelet basis (created by the shifts and dilations of the wavelet functions (4)) contains of the wavelet functions where k ∈ J p = {1, 2, . . . , p − 1}, j ∈ Z, and n is taken as an element of the m-direct product of factor group That is, n belongs to {x ∈ Q p : {x} p = x}.

Vladimirov's Operator and p-Adic Lizorkin Spaces
Introduced by V.S. Vladimirov [6], pseudo-differential operator A (on the field of p-adic numbers) in an open set O ⊂ Q p is given by which acts on C-valued functions ϕ(x) of p-adic arguments x ∈ O. Here we assume that functions ϕ(x) are extended by zero from the set O on whole space Q p , andφ(ξ) are their Fourier transforms recalled previously. The function A(ξ, x), ξ ∈ Q p , x ∈ O is called symbol of the operator A.
In [64,65] Lizorkin presented spaces invariant under the real actions of fractional operators. These spaces can be defined in p-adic case. In view of [25,66], the p-adic Lizorkin space of test functions is described as follows: Clearly, Ψ, Φ = ∅. Regarding the fact that Fourier transform is a linear isomorphism D(Q p ) into D(Q p ), one can see that Ψ, Φ ∈ D(Q p ). The space Φ(Q p ) can be decorated with the topology of the space D(Q p ), which turns it into a complete space. The space Φ can be determined by the following characterization: In addition, the space Φ = Φ (Q p ) as the topological dual of Φ is said to be the Lizorkin space of p-adic distributions (see also [25]). The Vladimirov operator D α (initially introduced by Taibleson) also is included in the class (6) with symbol A(ξ) = |ξ| α p , i.e., where α ∈ C. The formula (7) can be rewritten as a convolution of the following functions: where the distribution κ α ∈ Φ (Q p ) is called the Riesz kernel given by x ∈ Q p and Γ p (α) = 1−p α−1 1−p −α is the Γ-function (for more details see [6]). The domain of D α is given by We remark that all the concepts as above can be reconsidered in multidimensional p-adic field Q n p which is not in our considerations in the current paper.

Solvability of the p-Adic Navier-Stokes Equation
This section deals with the main result of our paper. Namely, we will apply the wavelet and fixed point theories to show the existence of solutions for our problem. Indeed, we study the existence of solution of the p-adic pseudo-differential equation in [0, ∞) given as (1).
Let us assume there exists u ∈ U I := C 1 (I, R) ∩ C(I, M(D 1 )) ∩ C(I, M(D 2 )), for any interval I in terms of wavelet functions ψ k;jn (x) with coefficients u k;jn (t), that is, u(x, t) = ∑ u k;jn (t)ψ k;jn (x), ϕ k;jn := u k;jn (0) = ϕ(x), ψ k;jn (x) . Moreover, On the other hand, taking into account that it implies that Replacing (8), (9) and (10) in Equation (1) and regarding the action of Vladimirov operator on u we derive the following infinite system: where G ı (t) = G(t, x), ψ ı (x) . Removing the index ı and fixing j yields the following equation Generally, the solution of this nonlinear differential equation cannot be presented explicitly but for the special cases we have: and using the initial condition we see that c = H ı ((ϕ ı ) ı ), that is, the solution of Equation (1) takes the following form • If F(û) = Aû and G(t) = Bt then the solution is represented in the parametric form: To find more cases of Equation (13) we refer the reader to see ( [67], [Section 1.6.3]).
In the sequel, we focus on the solvability of Equation (13) in the general form.

Solvability of Infinite System (12) over the Sequence Space c 0
Let us first convert the infinite system of differential equations (12) into the following infinite system of integral equations The existence theory concerning the infinite systems of integral equations is satisfactorily developed up to now and we are interested in study the system (14) in the Banach sequence space c 0 containing sequences of real numbers converging to zero. We recall that c 0 is a closed subspace of c ⊂ ∞ as the space of convergent sequences. It turns out this space is very convenient and natural for investigations of infinite systems both differential and integral equations. In what follows, we proceed our study in the Banach space c 0 including of real sequences converging to zero with the standard norm û c 0 = {|u ı | : ı ∈ J} forû = (u ı ) ı . Note that the index has been altered following our notation and the subjected problem.
In the following we intend to apply the generalized theorem of Arzéla (see also [68]) which describes a criterion of compactness in the space C(I, E) for the arbitrary interval I and the Banach space E. Suppose that G I := max t∈I |G(t)| < ∞ and the interval I = [0, T] for T > 0 is given. To investigate the solvability of nonlinear pseudo-differential equation (1) it only needs to focus on the existence of u(t) from the system (14). To do this, let us first present the following well-known fixed point result.
In what follows, system (14) will be investigated under the following hypotheses.
(i) The functions F ı are given on the set I × R ∞ and take real values (ı ∈ J). Further, the operator F is defined on the space I × c 0 in the following way: There exist nonnegative functions α ı (t) and β ı (t) defined, integrable and uniformly bounded on I and such that T 0 β ı (s)ds < 1. Furthermore, the function sequence ( t 0 α ı (s)ds) converges monotonically to zero at each point t ∈ I while the function sequence ( t 0 β ı (s)ds) is non-increasing at each point t ∈ I and the following estimate is satisfied: The functions G ı + ϕ ı : I → R are continuous on I and the sequence (|G ı (t) + ϕ ı |) ı converges monotonically to zero at each point t ∈ I.
In the following we recall the well-known Schauder fixed point theorem which is crucial to present our result.  1])). Let U be a nonempty and convex subset of a normed space E. Let T be a continuous mapping of U into a compact set K ⊂ U. Then T has a fixed point.
Proof. Indicate by S the subset of space C(I, c 0 ) including all functionsû(t) = (u ı (t)) ı so that for ı ∈ J and t ∈ I, where ≤ in index is the usual partial order in J, and a ı (t), b ı (t) are defined in the following way: Remark that the functions a ı (t) and b ı (t) are nondecreasing functions on the interval I and non-increasing sequences. Besides, from the assumptions, it follows that the functional sequences (a ı (t)) and (b ı (t)) converge uniformly on I to the function vanishing identically on I.
Let us assume the operator Γ defined on the space C(I, c 0 ) as follows: Notice that the operator Γ maps the set S into itself. Indeed, fix arbitrarily ı andû(t) ∈ S. Then for  ≥ ı, we get Now we prove that the operator Γ is continuous on the set S. Consider > 0 arbitrarily fixed andû 0 ∈ S. Then, taking into account the equicontinuity of the family of functions revealed in assumption (i) let us take δ = δ(û 0 , ), i.e., forû ∈ S such that û −û 0 c 0 ≤ δ we derive ( Fû)(t) − ( Fû 0 )(t) c 0 ≤ for each t ∈ I. Moving forward, which implies the desired claim. Now, let us take the set S 1 = ΓS. Recall that this set contains equicontinuous functions on I. In fact, taking an arbitraryû = (u ı ) ı ∈ S, and bringing in mind our hypotheses, we conclude Since the function sequences (α ı (t)), (β ı (t)), (a ı (t)), and (b ı (t)) are uniformly bounded on I and the function sequence (G ı (t) + ϕ ı ) is equicontinuous on I, from the above estimation, we conclude that the set S 1 = ΓS is equicontinuous on I.
Suppose that S 2 stands for ConvS 1 (i.e., the closed convex hull of the set S 1 ). Clearly, S 2 is closed, bounded, and equicontinuous on I. Furthermore, ΓS 2 ⊂ S 2 ⊂ S.
From another point of view, forû ∈ S, we get Since the sequence (a ı (t) + b ı (t)) converges uniformly on I to the function vanishing identically on I, we conclude that for each > 0, there exists an index ı 0 such that |(Γû) ı (t)| ≤ , for each ı ≥ ı 0 and for any t ∈ I. Therefore, by virtue of the criterion of compactness in the space c 0 as mentioned before, we infer that for each t ∈ I, the set S 1 (t) is relatively compact in the space c 0 . The above arguments allow us to deduce that the set S 2 is relatively compact in the space C(I, c 0 ). Besides, the closedness of S 2 yields that it is compact. Hence, keeping in mind that Γ transforms continuously the set S 2 into itself, we result (by the Schauder fixed-point principle) that the operator Γ has a fixed point in the set S 2 being a solution of our problem. This completes the proof.
then one can easily utilize the Banach contraction principle and find the unique solution for the subjected system.
Theorem 3 implies the following result immediately.

Adomian Decomposition Method (ADM)
In this section we give standard description of the ADM to find the solution of Equation (1) in which we proved its existence in the last section. Consider the general equation where u is the function subjected to be found, L is the linear differential operator of higher order which is simply invertible. Suppose that its inverse is L −1 and it will be an integral operator, N is the nonlinear operator, R is the remaining linear part and g is a given function (source). Taking L −1 to both sides of (16) we get: hence, where φ is chosen from the initial conditions or from the boundary conditions or both, it depends on how we select differential operator that solve the given problem. The ADM considers that solution u of the functional equation can be decomposed into infinite series u = ∞ ∑ n=0 u n and the nonlinear term N(u) can be expressed as infinite series Nu = ∑ ∞ n=0 A n where the A n 's are the Adomian polynomials, which depend upon u 0 , u 1 , . . . , u n . We recall that the Adomian polynomials A n 's are first constructed by Adomian in 1992, he gave a general formula to determine the values of A n 's: Therefore, Equation (16) takes the following form: Now from Equation (17), we can derive the solution algorithm as follows: Given u 0 , the other terms of u can be defined, respectively. Hence, the existing solution u of Equation (16) can be determined by the series of recursive sequence u n . To do this, let us rewrite the nonlinear differential equation (13) as beloŵ Assumingû = ∑ ∞ n=0û n , and applying ADM we obtain , n = 0, 1, 2, . . . . (20) where φ would be specified by the initial condition in (8). Moreover, the convergence of Adomian's decomposition method is discussed in the Appendix A.
Since we have no information about the form of G(t, x) and the initial condition of the problem (1), we present a model to see how ADM works. (1), let us suppose the non-homogeneity term G(t, x) = t s ln |x| p , s = −1, and the initial condition u(0, x) = ϕ(x) := Ω(|x| p ) which is the refinable function.

Example 1. In Equation
To formulate the function G in terms of basis ψ k;jn , using the symbol ξ := p j x − n we arrive at G k;jn (t) = G(t, x), ψ k;jn (x) = Q p t s ln |x| p · ψ k;jn (x)dx For the case k = 0 we get G 0;jn (t) = t s p −3j 2 ln p · (j − ord p (n)), and, otherwise, for the case 1 ≤ k ≤ p − 1, we have G k;jn (t) = 0.
Otherwise, for the case 1 ≤ k ≤ p − 1, using the fact that B r we see that On the other hand, looking at (11) one can see that Hence, Now, we are ready to apply ADM as follows.

Concluding Remarks
The use of tree-like (ultrametric) geometry is the promising direction in modeling of fluids' transport through capillary networks in porous disordered media. Such geometry approximates fractal (and multi-fractal) structures in Geo-conduits.
Theory of linear dynamical equations (especially, p-adic) is well developed and its application to geophysics (see [36]) did not demand essential mathematical efforts. However, as well as in Euclidean geophysical modeling, the basic ultrametric dynamical equations are nonlinear. One of such equations, an analog of the Navier-Stokes equation, was derived in recent paper [38]. Its study posed a variety of new problems. This study is especially complicated in the absence of the general theory of nonlinear (pseudo-)differential equations on p-adic spaces; just the first steps in this direction were done in articles [37,43]. The present paper is the important step towards establishing theory p-adic Navier-Stokes equation.
Fractal and multifractal mathematical models are widely used for diagnostic of hydrocarbon-reservoirs stratigraphic patterns anisotropy (see, e.g., [44,45]). The concrete "on field applications" are based on software; one of the promising complexes of diagnostic programs was developed by the research group of K. Oleschko. This complex was actively used in realization of the projects for Mexican oil-industry, e.g., the project SENER-CONACYT-Hidrocarburos, Yacimiento petrolero como un reactor fractal. Development of ultrametric models for oil transport plays the important role in mathematical justification of application of (multi-)fractal models for Mexican Petroleum Industry.