Under arid or semiarid conditions, evaporation demands usually exceed the ability of an exposed porous medium to provide liquid-phase water. The water content of subsurface zones of coarse-grained rocks or sandy soils is usually far below its residual value, and consequently only the gas-phase water can flow through. These, mostly up to a few centimeters thick, zones are referred to as vapor zones, dry surface layers or evaporation zones e.g.,.

The extent and development of the dry surface layer significantly affect the material's decay: directly by changes in its wetness, by frost and particularly by salt weathering, since dissolved salts are transported by the capillary water and form crystals at places of evaporation e.g.,.

The phenomenon is preconditioned by the fact that the porous medium becomes impervious to the liquid-phase water if the water content becomes sufficiently small. A limit arises inside the porous medium behind which the transport of water is only possible in the form of vapor. Such a soil profile can be divided into two parts that can be referred to as the zone of water flow and the zone of vapor diffusion. This formulation, however, is too vague and leaves the intermediate zone, its extent and its nature unclear.

Several studies were published, giving a detailed description of the evaporation process and the development of the transition zone e.g.,. These papers are mostly focused on special problems, unlike the present paper, which studies the three-dimensional problem under general transport conditions.

studied the problem, considering both hydraulic and thermal processes, and detected a narrow transition layer at the bottom of the dry surface layer. Another approach connects the interface between the region of water flow and the region of vapor flow, denoted as the evaporation front, with a critical value of the water content that can be determined directly from such porous-medium characteristics as hydraulic conductivity and vapor diffusivity.

studied the problem of water vapor transport through a region of dry material from a receding evaporation front. In the paper, the heat balance equation was involved in the final system of equations, the evaporation front was considered to be a sharp interface between the saturated zone and the dry (without liquid water) zone, the front was fixed and given a priori, and the liquid water was unmovable.

studied the time–space development of discontinuities in a one-dimensional porous medium. Liquid water, vapor and mixture of liquid water and vapor were assumed in the void space, and two governing equations, water and heat flow, were considered. No particular interfaces were defined, and discontinuities, in general, were studied in time–space. started their study with similar assumptions concerning the governing laws and investigated the resulting transition surfaces and conditions of loss of their stability.

Unlike these studies, the present paper aims to define the evaporation front by means of porous-medium characteristics and to formulate the law of its motion generally not involving any particular law governing the water transport. This approach makes it possible to use any set of flow and transport laws when formulating a problem of the evaporation front motion.

Several methods using dyes were developed to visualize the dry and wet regions within soil or rock profiles e.g.,. These methods proved their efficiency in laboratory conditions when utilized to visualize the dry surface layers and, in particular, the evaporation front positions. The applied water–dye solutions increase their concentration at places of evaporation and indicate these places by changing their color.

A special method was developed that minimizes the medium destruction and is usable under the field conditions. A very thin rod covered by a layer of color is inserted into a narrow hole drilled to the investigated material where there is the sought evaporation front. The present liquid-phase water colors the corresponding part of the rod showing its extent.

A number of experiments aiming at seeking and visualizing the evaporation front, see and , show that its position can be detected as a sharp line. The present paper tries to respect this experimental result in the definition presented below.

The goal of this paper is to give an exact definition of the evaporation front and to formulate the law of its motion.

We study such processes of water transport in porous media, where both the fluid phases, gaseous and liquid, are present and evaporation is taken into account. A porous-medium domain is considered that is in contact with a wet neighborhood at a part of its boundary, while at the other part of the boundary, it is in contact with a dry neighborhood. Here wet and dry are understood as containing liquid water and without liquid water, respectively. Such a part of the domain's boundary which is open to the atmosphere is considered to be the dry contact.

Under these conditions, there necessarily exists a set of points inside the studied domain or upon its boundary that makes an interface between the wet medium (porous medium) and the dry medium (porous medium or air). In view of the above-introduced terms, these points can be considered to be points of the evaporation front.

Generally, the porous-medium profile can be divided into three parts: (a) the dry zone, where just two phases, solid and gaseous (air), are present and water exists in the form of vapor as a component of the gaseous phase, (b) the wet zone, where the movable liquid water exists, and (c) the intermediate zone, where the liquid water is present but only in such a contact with the solid phase that makes it unmovable. Here, such liquid-phase water is understood as movable that moves due to the hydraulic head gradient.

The evaporation front does not exist in itself; it is a matter of definition. It seems natural to place the evaporation front in the intermediate zone or in an interface between the intermediate zone and one of the neighboring zones. It can be expected that during the process of evaporation, the depth of the intermediate zone will become small. The present water evaporates quickly due to its immobility, its small amount and contact with the solid phase. In view of this and the fact that experimentally the evaporation front can be indicated as a sharp interface between two neighboring zones, we assume that the extent of zone (c) can be neglected, and the evaporation front is defined as the common boundary of the zone without liquid water and the zone with movable liquid water. The concept evidently enables existence of a jump in water content values.

We do not consider the temperature distribution and heat flow and balance, since, by virtue of its definition, the evaporation front results from the water transport data. Though unknown, the heat flow within the profile provides the latent heat of vaporization that is necessary for the evaporation resulting from the actual process of water transport.

The evaporation front changes its position with time according to the outer conditions. Its shape and motion result from mutual relations (water transfer) between the wet zone and the dry zone. The front moves towards the wet region if the evaporation exceeds the flow of the liquid water towards the interface through the wet zone and vice versa. Since the evaporation front inside porous media, e.g., in a rock massif, is difficult to detect, mathematical modeling becomes an important tool always if the knowledge of its position and motion is required.

In what follows, all the introduced characteristics are macroscale porous-medium characteristics; e.g., a domain is a macroscale domain, a surface is a macroscale surface, etc.

Denote by Ω, Ω⊂R3 the domain in space and by (0,T) the time interval in which we study the transport process and suppose that the movable liquid-phase water occupies an open part Gw of the time–space domain G (i.e., the water content is positive and sufficient to enable the water flow in Gw), where 1G=(0,T)×Ω. We further define 2Gd=G\Gw‾; by virtue of our assumptions, θ=0 in Gd and θ>0 in Gw, where θ denotes the water content.

For any time t∈(0,T) we define the wet zone Ωtw and the dry zone Ωtd by putting 3Ωtw={x∈R3;(t,x)∈Gw}andΩtd={x∈R3;(t,x)∈Gd}.

It holds that 4Ωtw∩Ωtd=∅,Ωtw‾∪Ω‾td=Ω‾fort∈(0,T).

We define the evaporation front γt at time t∈(0,T) as 5γt=Ωtw‾∩Ω‾td.

This definition covers the wet–dry interfaces inside Ω. In order to enable the location of the evaporation front in the domain's surface, we define the wet and dry contacts from outside of Ω. The sets Ωtw and Ωtd divide the boundary ∂Ω into wet and dry parts with respect to wet–dry conditions inside Ω at t∈(0,T). Denote by Btw and Btd such two parts of ∂Ω that are at time t in contact with wet and dry conditions outside Ω‾, respectively, and define a boundary point x∈∂Ω as belonging to the evaporation front at time t if it satisfies 6x∈Ω‾tw∩Btdorx∈Ω‾td∩Btw. The complete evaporation front γt at time t∈(0,T) is then 7γt=Ω‾tw∩Ωtd‾∪Ω‾tw∩Btd∪Ω‾td∩Btw.

The image of the evaporation front in time–space is 8Γ={(t,x)∈R4;x∈γt,t∈(0,T)}. We assume that Γ is a smooth hypersurface in R4 and denote by νΓ∈R4 the unit vector normal to Γ that points out of the wet part or into the dry part of G. Since we defined the wet and dry boundaries Btw and Btd of Ω at time t with respect to the outer conditions, this orientation has sense everywhere on Γ.

Suppose that (t,ξ)∈Γ and 9νtΓ(t,ξ)≠0, where νΓ=(νtΓ,ν1Γ,ν2Γ,ν3Γ). Then the hypersurface Γ can be in a certain neighborhood of (t,ξ) expressed by a function τ in the form of the equation 10t=τ(x1,x2,x3). Then νtΓ<0 implies the existence of a positive value ϵ such that (t+ϑ,ξ)∈Gw for ϑ∈(0,ϵ). Consequently, the evaporation front γt moves at its point ξ towards the dry zone if νtΓ is negative at (t,ξ) and vice versa.

The position of the evaporation front results from the mutual relations between the water transport in the wet zone and in the dry zone. Denote by n the porosity, by w the volumetric flux density of liquid water in the wet zone, by vd and vw the volumetric flux density of the gaseous phase in the dry zone and in the wet zone, and by bd and bw the water vapor flux density by diffusion in the gaseous phase within the dry zone and the wet zone. Let furthermore cd and cw denote the water vapor concentration in the gaseous phase within the dry zone and the wet zone. We suppose that functions bd,cd,vd are continuous in Gd‾ and functions bw,cw,vw,w are continuous in Gw‾.

The evaporation front is not connected with a fixed set of mass points, and the problem of its motion is not a problem of the particle tracking. The evaporation front moves in such a direction and with such a velocity that are given by the balance of mass of water. Since the tangential motion of the evaporation front at its point does not change, the front's position and the evaporation front move at each point in the direction of its normal.

Let ξ=ξ(t) be a point upon the evaporation front at time t: ξ∈γt. Let δt be an elementary time step and δS an elementary surface surrounding the point ξ, δS⊂γt. Denote by ν(t,ξ) the unit normal vector to γt at point ξ, oriented out of the wet zone, and by δs the distance between the positions γt and γt+δt at ξ. Then the flow and transport of water coming to the elementary surface δS from the wet zone, ρw+bw+cwvv, push the front towards the dry zone, and the transport of water out of δS into the dry zone, bd+cdvd, pushes the front towards the wet zone. The excess of water coming to δS during the time interval δt, δtδS(ρw+bw+cwvw-bd-cdvd,ν), is compensated for by the mass of water δsδS(ρθ+cw(n-θ)-ncd) shifting the surface δS to its new position distant by δs in the direction ν. Here and in the sequel, (u,v), where u,v are two vectors, denotes the scalar product. This account gives the balance equation 11δSδt(ρw(t,ξ)+bw(t,ξ)+cw(t,ξ)vw(t,ξ)-bd(t,ξ)-cd(t,ξ)vd(t,ξ),ν(t,ξ))=δSδs(θ(t,ξ)ρ+(n-θ(t,ξ))cw(t,ξ)-ncd(t,ξ),) where values of higher orders of δt and δs are neglected. The idea of this calculation can be seen in Fig. . It shows two positions of the front, at times t and t+δt, and two possible vectors, W and B, where, for the sake of simplicity, B denotes the vector sum of fluxes by diffusion, 12B=bd(t,ξ)-bw(t,ξ), and W denotes the sum fluxes due to water flow and advection by the gaseous phase, 13W=ρw(t,ξ)+cw(t,ξ)vw(t,ξ)-cd(t,ξ)vd(t,ξ). The depicted directions of these vectors suggest simultaneous flooding of the dry zone and drying of the wet zone. These two processes act against each other; the figure shows that the vapor transport predominates and the front moves towards the wet zone.

The limit form of Eq. () for δt approaching zero gives 14∂s∂t(t,ξ)=(ρw(t,ξ)+bw(t,ξ)+cw(t,ξ)vw(t,ξ)-bd(t,ξ)-cd(t,ξ)vd(t,ξ),ν(t,ξ))cw(t,ξ)(n-θ(t,ξ))-ncd(t,ξ)+ρθ(t,ξ). By this equation and by the assumption of the direction of the front motion, we obtain the following ordinary first-order differential equation 15dξdt(t)=(ρw(t,ξ)+cw(t,ξ)vw(t,ξ)+bw(t,ξ)-bd(t,ξ)-cd(t,ξ)vd(t,ξ),ν(t,ξ))cw(t,ξ)(n-θ(t,ξ))-ncd(t,ξ)+ρθ(t,ξ)ν(t,ξ). Equation () is the governing equation of the evaporation front motion in the interval (0,T). The water density and porosity are presented as constants in Eqs. () and (). Such an assumption is evidently not necessary, and the equations remain unchanged for n=n(t,ξ) and ρ=ρ(t,ξ).

In the cases we commonly meet in connection with problems of evaporation, the flow of the gaseous phase is restricted to balancing the changing volume of liquid water, i.e., 16‖w‖≈‖vd‖≈‖vw‖, and since 17ρ‖w‖≫cd‖vd‖andρ‖w‖≫cw‖vw‖, the advective transport of water vapor can be neglected. The governing equation becomes 18dξdt(t)=ρw(t,ξ)+bw(t,ξ)-bd(t,ξ),ν(t,ξ)cw(t,ξ)(n-θ(t,ξ))-ncd(t,ξ)+ρθ(t,ξ)ν(t,ξ).

The front's motion reflects the proportions between the water flow and transport out of the front and towards the front which are given by the laws of flow and transport in the wet zone and in the dry zone. In order to evaluate flow and transport in porous media, Darcy's law and Fick's law can be utilized: 19w=-k(h)gradx3+h,bd=-Ddgradcd,andbw=-Dw(θ)gradcw, where the third coordinate x3 is oriented vertically upwards, h is the pressure head, k is the hydraulic conductivity, and Dd and Dw are the coefficients of water diffusion in air within the porous medium.

In the dry zone, water is present in the form of water vapor, and its motion is governed by the continuity equation with the use of Fick's law: 20∂∂t(ncd)-∂∂xiDd∂cd∂xi+∂∂xi(cdvid)=0. The vapor motion in the wet zone is governed by the same laws, 21∂∂t(n-θ)cw-∂∂xiDw(θ)∂cw∂xi+∂∂xi(cwviw)=0, and the motion of liquid-phase water in the wet zone is governed by Richards' equation: 22∂θ∂t-∂∂xik(h)∂x3∂xi+∂h∂xi=0.

By virtue of the introduced theory, the evaporation front motion is governed by Eqs. (), (), () and () or (). The unknown functions are θ and h, connected by the retention curve, cw, cd and ξ defined in Gw‾, and Gd‾ and [0,T). The functions n,k,D,v and ρ are supposed to be known or given by additional equations.

In this way, Eqs. () and () in Gw and Eq. () in Gd stand for a coupled moving boundary problem, and Eq. () is a condition imposed upon the movable common part of the boundaries of the domains. The unknown function ξ defined in [0,T) is then given as the position of the moving boundary.

Another possible formulation of the problem is to solve the ordinary differential Eq. () in the interval [0,T), where the right-hand side of the equation is given as the solution of the Eqs. (), () and () defined in Gw and Gd.

The latter approach was utilized when solving the problem presented in Sect. .

Let now the studied domain be an interval Ω⊂R1: 23Ω=(0,L),Ωtw=(0,ξ(t)),Ωtd=(ξ(t),L), where the function ξ(t) that denotes the position of the evaporation front at time t is a scalar function defined in (0,T). The sets Gw and Gd are 24Gw={(t,x)∈R2;t∈(0,T),x∈(0,ξ(t))}andGd={(t,x)∈R2;t∈(0,T),x∈(ξ(t),L)}. The image of the evaporation front in time–space is 25Γ={(t,x)∈R2;ξ(t)-x=0,t∈(0,T)}, and νΓ∈R2, the unit normal to Γ, is 26νΓ(t,ξ(t))=-dξdt(t),1/dξdt(t)2+1. The above-presented assertion concerning the direction of the evaporation front motion with respect to sgn(νtΓ) is now evident.

By virtue of the introduced theory, the law of the evaporation front motion is given by Eq. (). In one space dimension, since ν(t,ξ) is either 1 or -1, the equation reads as 27dξdt(t)=ρw(t,ξ)+cw(t,ξ)vw(t,ξ)+bw(t,ξ)-bd(t,ξ)-cd(t,ξ)vd(t,ξ)cw(t,ξ)(n-θ(t,ξ))-ncd(t,ξ)+ρθ(t,ξ) and its simplified form () is 28dξdt(t)=ρw(t,ξ)+bw(t,ξ)-bd(t,ξ)cw(t,ξ)(n-θ(t,ξ))-ncd(t,ξ)+ρθ(t,ξ).

In order to complete the problem formulation, to determine the right-hand-side function, the one-dimensional form of Eqs. () to () can be utilized.

In the framework of wider research that concerns evaporation from rock surfaces and the time dependence of evaporation front positions, an experiment was carried out. It was an unpublished auxiliary experiment performed by several of the author's colleagues, and since it was not sufficiently documented from the point of view of this study, it cannot be simulated. Nevertheless, part of its results can be utilized here in order to present an example of possible use of the achieved theoretical results. The missing data were simply chosen, not optimized. The following description should be understood as a problem formulation, not as a documentation of measurements.

A cylinder-shaped sample of the studied rock was put into position with the horizontal axis. The jacket of the cylinder was insulated so that no water (of any phase) could penetrate, and the motion of water through the sample was possible only in the horizontal direction along the cylinder's axis.

The length of the sample was L=49 mm, and the length of the time interval was T=63 d. One open end of the sample, say x=0, was equipped so that it was possible to measure the pressure head and the rate of water inflow into the sample at this point. The obtained discrete data were approximated by smooth functions h0(t) and w0(t) (pressure head and volumetric flux density), h0,w0∈C1(0,T), that were utilized as the imposed boundary conditions.

Figure  shows the pressure head values; the squares are the measured data, and the solid curve is their approximation h0.

The volumetric flux density was measured indirectly in the form of discrete values of the cumulative flux to the sample. In Fig. , the step function represents the measured data and the smooth function is the boundary condition w0. The integral values over (0,T) (the total inflow to the sample per unit surface) of both functions are equal. The other end of the sample, x=L, was left open to the outer (atmospheric) conditions which were not measured.

Soil moisture retention data and the hydraulic conductivity (at saturation) Ks were obtained elsewhere using samples of similar material. Making use of these characteristics, the Mualem and van Genuchten parameters θr,θs,α,m and n were determined and, hence, the hydraulic conductivity and the soil moisture retention curve as functions of either pressure head or water content were defined.

Similarly, the value of the diffusion coefficient of water in the gaseous phase within the porous medium, Dd, was obtained measured on samples of the utilized material. The characteristic surface layer ε was added from outside to the domain Ω, and its value was taken from the paper of , where the layer is referred to as “calibrated L” and also as “boundary layer”.

The fluorescein visualization method was utilized to detect the front's position during the experiment. Fluorescein was applied at the water input side of the sample, and a set of couples (t,ξ), time and the front's position was registered. The data were calibrated (using a simple linear transformation) to agree with the value measured after finishing the experiment.

The following problem was formulated. In the interval [0,T], the solution t↦ξ(t), ξ∈C1(0,T)∩C[0,T], of Eq. () is sought that satisfies the initial condition 29ξ(0)=L, where ρ and n (density and porosity) are known constants.

Since neither functions cd and cw nor their boundary values at x=0,ξ(t) and L were measured during the experiment, the following two assumptions were introduced.

(*) The vapor concentration has a constant value cf at the evaporation front in the gaseous phase and at points of positive value of water content θ, i.e., at points of contact with liquid-phase water.

(**) There are steady-state outer atmospheric conditions during the process giving a constant value ce of vapor concentration in the gaseous phase at x=L.

The assumptions do not contradict each other at the point ξ(0), since the characteristic surface layer ε was accepted in the model. As presented above, the right-hand side of Eq. () can be determined by solving the related initial boundary value problems with Eqs. (), (), and ().

In view of assumption (*) and Eqs. () and (), it holds that 30cw(t,x)=cfandbw(t,x)=0inGw.

Applying assumption (*), Eq. () reads as 31∂θ∂t-∂vw∂x=0. The comparison with Richards' Eq. () gives vw=-w; the gaseous phase continuously replaces the leaving water. A similar result has already been expected even for more general cases; see the estimations utilized when replacing Eq. () by Eq. (). Making use of these results, Eq. () becomes 32dξdt(t)=ρw(t,ξ)-b(t,ξ)θ(t,ξ)(ρ-cf), where the unknown parameters on the right-hand side of Eq. () are solutions of the following two initial boundary value problems and, since bw, cw and Dw do not appear in what follows, b, c and D stand for bd, cd and Dd.

To determine functions θ and w, function h defined in Gw is sought that satisfies Eq. (), now in the form 33∂θ∂t-∂∂xk(h)∂h∂x=0 and the conditions 34h(0,x)=hi(x)in(0,L) and 35h(t,0)=h0(t)and-k(h(t,0))∂h∂x(t,0)=w0(t)in(0,T), where θ(h) is the retention curve. Functions h0 and w0 are known from the experiment, and function hi, the initial pressure head distribution, was not measured and has to be determined.

Incorporating the characteristic surface layer into the problem formulation, we change the set Gd to 36Gd={(t,x)∈R2;x∈(ξ(t),L+ε),t∈(0,T)}. Now, c solves the equation 37∂∂t(nc)-∂∂xD∂c∂x=0 in Gd and satisfies the conditions 38c(0,x)=ci(x)in(L,L+ε) and 39c(t,ξ(t))=cfandc(t,L+ε)=cein(0,T), where ci is the initial distribution of the water concentration in the gaseous phase. The values ce and cf were chosen with respect to these requirements: the relative humidity at the front was 100 %, the outer relative humidity was between 30 % and 50 % and the temperature was between 18 and 23 ∘C.

In order to define the initial state of the sample, functions ci and hi, it was assumed that the process started from the steady-state water vapor diffusion determined in (L,L+ε) by the boundary values cf and ce and from the steady-state water flow determined in (0,L) by the initial conditions h0(0) and w0(0). The error involved in functions ci and hi vanishes soon, since the sample is small and the imposed boundary conditions take over the dominant role.

Our requirement for the initial conditions is easy to satisfy in the case of function ci in the domain Ωt=0d, since, in view of Eq. (), the governing equation is 40ddxDdcidx=0in(L,L+ε). With respect to the boundary conditions, the solution reads as 41ci(x)=ce-cfε(x-L)+cf.

In the case of function hi and domain Ωt=0w, the governing equation is 42ddxρk(hi)dhidx=0in(0,L), with the initial conditions 43hi(0)=h0(0)and-k(hi(0))dhidx(0)=w0(0). The solution to the problem (), (), is equivalent to the solution to the problem 44dhidx=-w0(0)k(hi)in(0,L), with the initial condition hi(0)=h0(0). Let the maximum solution to this initial value problem be defined in [0,λ). It can be shown that, in our case, λ>L and the wet zone reach the point x=L. Hence, our choice of the initial condition hi does not contradict the initial condition ().

Now the problem (), (), can be solved numerically. The utilized values are the Mualem and van Genuchten parameters α=0.0223 cm-1, n=1.99, and m=1-1/n, the saturated hydraulic conductivity Ks=0.0071 cm s-1 and the diffusion coefficient of water in the gaseous phase within the porous medium Dd=0.045 cm2 s-1. Figure  shows three solutions obtained for three different couples of chosen parameters and the boundary values of water concentration in the gaseous phase. The chosen values were (1)

cf=2.27×10-5 g/cm3, ce=9.35×10-6 g/cm3,

The evaporation front has been defined as an interface separating two different zones, wet and dry, inside or upon the boundary of a porous-medium domain. The exact definition of these zones, presented in this paper, is based on the form of water they contain. Subsequently, the law of the evaporation front motion was formulated in the form of the vector Eq. (). Since the law is based on the complete mass balance of water, i.e., liquid water and water vapor, it holds generally and does not need any additional account of energy. The laws of heat transfer and heat balance do not affect the presented equations which define the evaporation front motion. On the contrary, solving problems that are fully determined by water transport data, the equations of the evaporation front motion can give certain insight into the energy requirements of such processes, e.g., the final part of Sect. .

and studied the process of evaporation from soils with particular attention to the phase change and found that nonequilibrium models yield better agreement with experimental data than equilibrium models. The nature of the phase change process does not affect the results presented here directly, since the equation of the evaporation front motion requires other kinds of data. The process of phase change enters Eq. () through its actual effect on the transport of water. On the other hand, the constitution laws like Darcy's law or the retention curve that may be utilized when solving problems with Eq. () are equilibrium laws. In the example presented in Sect. , equilibrium laws were utilized. However, the governing Eq. (), being general, makes it possible to use nonequilibrium laws as well. presented a general nonequilibrium approach to two-phase systems that keeps Darcy's law valid.

and investigated the process of evaporation from the top of an initially saturated vertical column. They introduced the term “characteristic length” as the distance between the surface and the receding drying front (interface between the saturated zone and the unsaturated zone) and described different stages of the evaporation process. No evaporation front was introduced. By virtue of the present theory, the evaporation front cannot move in the positive direction of ν(t,ξ) if 49ξ(t)∈Ωtw‾∩Btd, i.e., if ξ(t) is a point of a “dry-from-outside” part of the domain boundary. Hence, see Eq. (): the front does not move until the condition 50((ρw(t,ξ)+cw(t,ξ)vw(t,ξ)+bw(t,ξ)-bd(t,ξ)-cd(t,ξ)vd(t,ξ)),ν(t,ξ))<0 is satisfied. In the simplified one-dimensional case, Eq. (), this condition reads as 51ρw(t,ξ)+bw(t,ξ)-bd(t,ξ)ν(t,ξ)<0, where ν(t,ξ)=1 if x∈Ωtw⇒x<ξ(t) and vice versa. Under conditions of sufficiently small values of |bd|, condition () is not satisfied for a period, and the evaporation front does not move. Consequently, the evaporation rate does not change significantly. In the solved problem, the chosen initial conditions and the size of ε make Eq. () valid even at t=0, and the front moves from the very beginning of the process. If the flux density of liquid water from inside of the wet zone exceeds the flux density b(t,L), either θ(t,L) increases or, being θ(t,L)=θsat, flux density ρw(t,L)-b(t,L) of liquid water discharges out of the porous-medium domain.

The characteristic surface layer ε was found experimentally, see , and accepted in this paper as a part of the measured data. Note that studied similar problems and also introduced a special diffusion layer outside the porous medium. From the viewpoint of moving front equations, the characteristic surface layer prevents infinite values of function bd at x=0 and t=0 which may be obtained when solving a problem with Eq. (). This possibility originates in the fact that the equation is a balance equation that contains an equilibrium law – Fick's law; for more on this problem and an alternate approach, see .

The process of evaporation alone does not determine the direction of the evaporation front motion. Since the denominator of the right-hand side of Eq. () is positive, the direction of the evaporation front motion is determined by the sign of the scalar product in the numerator. Consequently, both the processes of wetting or drying (increasing or decreasing the wet zone) can take place while evaporating water out of the profile; compare also Eqs. () and ().

The presented problem example and its numerical solutions were aimed at showing the ability of the theory to simulate real processes, not at getting an optimized agreement. Most of the measured parameters of the solved problem were obtained independently of the experiment. Only the pressure head and water flow data shown in Figs.  and were measured during the experiment and utilized as the imposed boundary conditions of the problem. The concentration of water in the gaseous phase, the functions cf and ce, and the initial values of functions h and w were not measured but chosen. For the sake of their simple interpretation by means of acceptable values of temperatures and relative humidities, cf and ce were kept constant. No method of fitting was applied, and a different choice of functions hi and wi and constants ce and cf can give a better agreement between the measured and computed values.

The presented theory is now prepared to prove its reliability in such problems that are fully documented and to be used when solving a wide range of problems of evaporation from a rock or soil profile.