Regularity for the stationary Navier-Stokes equations in bounded domain

Let u, p be a weak solution of the stationary Navier-Stokes equations in a bounded domain ^N, 5N . If u, p satisfy the additional conditions

Singular perturbations of two-point boundary problems for systems of ordinary differential equations

Asymptotic solutions of linear systems of ordinary differential equations are employed to discuss the relationship of the solution of a certain “complete” boundary problem.

$$\begin{gathered} \left\{ \begin{gathered} {\text{ }}\frac{{d{\text{ }}x_1 }}{{d{\text{ }}t}} = A_{11} (t,\varepsilon ){\text{ }}x_1 (t,\varepsilon ){\text{ }} + \cdots + A_{1p} (t,\varepsilon ){\text{ }}x_p (t,\varepsilon ) \hfill \\ \varepsilon ^{h_2 } \frac{{d{\text{ }}x_2 }}{{d{\text{ }}t}} = A_{21} (t,\varepsilon ){\text{ }}x_1 (t,\varepsilon ){\text{ }} + \cdots + A_{2p} (t,\varepsilon ){\text{ }}x_p (t,\varepsilon ) \hfill \\ {\text{ }} \vdots {\text{ }} \vdots {\text{ }} \vdots \hfill \\ \varepsilon ^{h_p } \frac{{d{\text{ }}x_2 }}{{d{\text{ }}t}} = A_{p1} (t,\varepsilon ){\text{ }}x_1 (t,\varepsilon ){\text{ }} + \cdots + A_{pp} (t,\varepsilon ){\text{ }}x_p (t,\varepsilon ) \hfill \\ \end{gathered} \right\} \hfill \\ {\text{ }}R(\varepsilon ){\text{ }}x(a,{\text{ }}\varepsilon ){\text{ }} + {\text{ }}S(\varepsilon ){\text{ }}x(b,{\text{ }}\varepsilon ) = c(\varepsilon ){\text{ }} \hfill \\ \end{gathered}$$     

Remainder estimates for the asymptotics of elliptic eigenvalue problems with indefinite weights

Let A be a positive self-adjoint elliptic operator of order 2m on a bounded open set Ω ?? ^k . We consider the variational eigenvalue problem (P) $$\mathcal{A}u = \lambda r{\text{(}}x{\text{)}}u,{\text{ }}x \in \Omega ,$$ , with Dirichlet or Neumann boundary conditions; here the “weight” r is a real-valued function on Ω which is allowed to change sign in Ω or to be discontinuous. Such problems occur naturally in the study of many nonlinear elliptic equations. In an earlier work [Trans. Amer. Math. Soc. 295 (1986), pp. 305–324], we have determined the leading term for the asymptotics of the eigenvalues λ of (P). In the present paper, we obtain, under more stringent assumptions, the corresponding remainder estimates. More precisely, let N ^±(λ) be the number of positive (respectively, negative) eigenvalues of (P) less than λ>0 (respectively, greater than λ<0); set r _± = max (±r, 0) and \(\Omega _ \pm = {\text{\{ }}x \in \Omega :r{\text{(}}x{\text{)}} \gtrless {\text{0\} }}\) . We show that $$N^ \pm {\text{(}}\lambda {\text{) = }}\mathop \smallint \limits_{\Omega _ \pm } {\text{(}}\lambda r{\text{(}}x{\text{))}}^{\frac{k}{{{\text{2}}m}}} {\text{ }}\mu \prime _\mathcal{A} {\text{(}}x{\text{) }}dx + 0{\text{(}}\left| \lambda \right|^{\frac{{k - 1}}{{{\text{2}}m}} + \delta } {\text{) as }}\lambda \to \pm \infty {\text{,}}$$ , where δ>0 and μ^′ _A (x) is the Browder-Gårding density associated with the principal part of A. How small δ can be chosen depends on the “regularity” of the leading coefficients of A, r _±, and of the boundary of Ω _±. These results seem to be new even for positive weights.     

Degenerate solutions of electron-beam equations

In this paper, exact solutions are constructed for stationary election beams that are degenerate in the Cartesian (x,y,z), axisymmetric (r,θ,z), and spiral (in the planes y=const (u,y,v)) coordinate systems. The degeneracy is determined by the fact that at least two coordinates in such a solution are cyclic or are integrals of motion. Mainly, rotational beams are considered. Invariant solutions for beams in which the presence of vorticity resulted in a linear dependence of the electric-field potential ? on the above coordinates were considered in [1], In degenerate solutions, the presence of vorticity results in a quadratic or more complex dependence of the potential on the coordinates that are integrals of motion. In [2] and in a number of papers referred to in [2], the degenerate states of irrotational beams are described. The known degenerate solutions for rotational beams apply to an axisymmetric one-dimensional (r) beam with an azimuthal velocity component [3] and to relativistic conical flow [1]. The equations used below follow from the system of electron hydrodynamic equations for a stationary relativistic beam $$\begin{array}{*{20}c} {\sum\limits_{\beta = 1}^3 {\frac{\partial }{{\partial q^\beta }}\left[ {\sqrt \gamma g^{\beta \beta } g^{\alpha \alpha } \left( {\frac{{\partial A_\alpha }}{{\partial q^\beta }} - \frac{{\partial A_\beta }}{{\partial q^\alpha }}} \right)} \right]} = 4\pi \rho \sqrt \gamma g^{\alpha \alpha } u_\alpha ,} \\ {\sum\limits_{\beta = 1}^3 {\frac{\partial }{{\partial q^\beta }}\left( {\sqrt \gamma g^{\beta \beta } \frac{{\partial \varphi }}{{\partial q^\beta }}} \right)} = 4\pi \rho \sqrt {\gamma u} ,\sum\limits_{\beta = 1}^3 {g^{\beta \beta } u_\beta ^2 + 1 = u^2 } } \\ \begin{gathered} \frac{\eta }{c}u\frac{{\partial \mathcal{E}}}{{\partial q^\alpha }} = \sum\limits_{\beta = 1}^3 {g^{\beta \beta } u_\beta } \left( {\frac{{\partial p_\beta }}{{\partial q^\alpha }} - \frac{{\partial p_\alpha }}{{\partial q^\beta }}} \right), \hfill \\ \begin{array}{*{20}c} {\sum\limits_{\beta = 1}^3 {\frac{\partial }{{\partial q^\beta }}(\sqrt \gamma g^{\beta \beta } \rho u_\beta ) = 0,u \equiv \frac{\eta }{{c^2 }}(\varphi + \mathcal{E}) + 1,} } \\ {cu_\alpha \equiv \frac{\eta }{c}A_\alpha + p_\alpha ,\alpha ,\beta = 1,2,3,\gamma \equiv g_{11} g_{22} g_{33} } \\ \end{array} \hfill \\ \end{gathered} \\ \end{array} $$ where q^β denotes orthogonal coordinates with the metric tensor g^ββ (β=1,2,3); A_α is the magnetic potential; A_α = (u_α/u)c is the electron velocity; ρ is the scalar space-charge density (ρ > 0); is the energy in eV; p_α is the generalized momentum of an electron per unit mass; η is the electron charge-mass ratio.     

Phragmén-Lindelöf theorems for nonlinear elliptic equations

We obtain theorems of Phragmén-Lindelöf type for the following classes of elliptic partial differential inequalities in an arbitrary unbounded domain \(\Omega \subseteq \mathbb{R}^n ,{\text{ }}n \geqq 2\) (A.1) $$\sum\limits_{i,j = 1}^n {\frac{\partial }{{\partial x_i }}\left( {a_{ij} 9(x)\frac{{\partial u}}{{\partial xj}}} \right)} + \sum\limits_{i = 1}^n {b_i (x,{\text{ }}u,{\text{ }}\nabla u)\frac{\partial }{{\partial x_i }}} \geqq f(x,{\text{ }}u)$$ where a _ij are elliptic in Ω and b _i ε L^∞(Ω) and where also a _ij are uniformly elliptic and Holder continuous at infinity and b _i = O(|x|⁺¹) as x → ∞; (A.2) $${\text{(A}}{\text{.2) }}\sum\limits_{i,j = 1}^n {a_{ij} (x,{\text{ }}u,{\text{ }}\nabla u)\frac{{\partial ^2 u}}{{\partial x_i \partial x_j }}} + \sum\limits_{i = 1}^n {b_i (x,{\text{ }}u,{\text{ }}\nabla u)\frac{\partial }{{\partial x_i }}} \geqq f(x,{\text{ }}u)$$ where a_ijare uniformly elliptic in Ω and b _iε L^∞(Ω); and finally (A.3) $${\text{div(}}\nabla u^p \nabla u {\text{)}} \geqq f{\text{(}}u{\text{), }}p > - 1,$$ where the operator on the left is the so-called P-Laplacian. The function f is always supposed positive and continuous. Moreover u is assumed throughout to be in the natural weak Sobolev space corresponding to the particular inequality under consideration, namely u ε. W _loc ^1,2 (Ω) ∩L _loc ^t8 (Ω) for (A.I), W _loc ^2,n(Ω) for (A.2), and W _loc ^1,p+2 (Ω) ∩ L _loc ^t8 (Ω) for (A.3). As a consequence of our results we obtain both non-existence and Liouville theorems, as well as existence theorems for (A.1).     

Forced convection heat transfer and friction characteristics in multi-exit -rectangular package with inclined tube bundles

An experiment was carried out to investigate the characteristics of the heat transfer and pressure drop for forced convection airflow over tube bundles that are inclined relative to the on-coming flow in a rectangular package with one outlet and two inlets. The experiments included a wide range of angles of attack and were extended over a Reynolds number range from about 250 to 12,500. Correlations for the Nusselt number and pressure drop factor are reported and discussed. As a result, it was found that at a fixed Re, for the tube bundles with attack angle of 45 ° has the best heat transfer coefficient, followed by 60, 75 and 90 °, respectively. This investigation also introduces the factors which can be used for finding the heat transfer and the pressure drop factor on the tube bundles positioned at different angles to the flow direction. Moreover, no perceptible dependence of Cand C on Re was detected. In addition, flow visualizations were explored to broaden our fundamental understanding of the heat transfer for the present study.     

The principle of Saint-Venant in linear and non-linear plane elasticity

A two-dimensional elastic isotropic body

$$R:{\text{ }}[0 < x{\text{ < 2}}l;{\text{0 < }}y < h]{\text{ with }}\frac{h}{l} \ll {\text{1}}$$     

Pressure drop in alternating curved tubes

Pressure drop measurements in the laminar and turbulent regions for water flowing through an alternating curved circular tube (x=h sin 2πz/λ) are presented. Using the minimum radius of curvature of this curved tube in place of that of the toroidally curved one in calculating the Dean number (ND=Re(D/2R _c )², it is found that the resulting Dean number can help in characterizing this flow. Also, the ratio between the height and length of the tube waves which represents the degree of waveness affects significantly the pressure drop and the transition Dean number. The following correlations have been found:

1. For laminar flow: $$F_w \left( {\frac{{2R_c }}{D}} \right)^{{1 \mathord{\left/ {\vphantom {1 2}} \right. \kern-\nulldelimiterspace} 2}} = F_s \left( {\frac{{2R_c }}{D}} \right)^{{1 \mathord{\left/ {\vphantom {1 2}} \right. \kern-\nulldelimiterspace} 2}} + 0.03,\operatorname{Re}< 2000.$$
2. For turbulent flow: $$F_w \left( {\frac{{2R_c }}{D}} \right)^{{1 \mathord{\left/ {\vphantom {1 2}} \right. \kern-\nulldelimiterspace} 2}} = F_s \left( {\frac{{2R_c }}{D}} \right)^{{1 \mathord{\left/ {\vphantom {1 2}} \right. \kern-\nulldelimiterspace} 2}} + 0.005,2000< \operatorname{Re}< 15000.$$
3. The transition Dean number: $$ND_{crit} = 5.012 \times 10^3 \left( {\frac{D}{{2R}}} \right)^{2.1} ,0.0111< {D \mathord{\left/ {\vphantom {D {2R_c }}} \right. \kern-\nulldelimiterspace} {2R_c }}< 0.71.$$

     

Two-phase flow in heterogeneous porous media I: The influence of large spatial and temporal gradients

In order to capture the complexities of two-phase flow in heterogeneous porous media, we have used the method of large-scale averaging and spatially periodic models of the local heterogeneities. The analysis leads to the large-scale form of the momentum equations for the two immiscible fluids, a theoretical representation for the large-scale permeability tensor, and a dynamic, large-scale capillary pressure. The prediction of the permeability tensor and the dynamic capillary pressure requires the solution of a large-scale closure problem. In our initial study (Quintard and Whitaker, 1988), the solution to the closure problem was restricted to the quasi-steady condition and small spatial gradients. In this work, we have relaxed the constraint of small spatial gradients and developed a dynamic solution to the closure problem that takes into account some, but not all, of the transient effects that occur at the closure level. The analysis leads to continuity and momentum equations for the-phase that are given by

A non-linear seales method for strongly non-linear oscillators

A non-linear seales method is presented for the analysis of strongly non-linear oseillators of the form % MathType!MTEF!2!1!+-% feaafiart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXafv3ySLgzGmvETj2BSbqefm0B1jxALjhiov2D% aebbfv3ySLgzGueE0jxyaibaiGc9yrFr0xXdbba91rFfpec8Eeeu0x% Xdbba9frFj0-OqFfea0dXdd9vqaq-JfrVkFHe9pgea0dXdar-Jb9hs% 0dXdbPYxe9vr0-vr0-vqpWqaaeaabiGaciaacaqabeaadaqaaqGaaO% qaamXvP5wqonvsaeHbfv3ySLgzaGqbdiqb-Hha4zaadaGaey4kaSIa% am4zaiaacIcacqWF4baEcaGGPaGae8xpa0JaeqyTduMaamOzaiaacI% cacqWF4baEcqWFSaalcuWF4baEgaGaaiaabMcaaaa!4FEC!\[\ddot x + g(x) = \varepsilon f(x,\dot x{\text{)}}\], where g(x) is an arbitrary non-linear function of the displacement x. We assumed that % MathType!MTEF!2!1!+-% feaafiart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXafv3ySLgzGmvETj2BSbqefm0B1jxALjhiov2D% aebbfv3ySLgzGueE0jxyaibaiGc9yrFr0xXdbba91rFfpec8Eeeu0x% Xdbba9frFj0-OqFfea0dXdd9vqaq-JfrVkFHe9pgea0dXdar-Jb9hs% 0dXdbPYxe9vr0-vr0-vqpWqaaeaabiGaciaacaqabeaadaqaaqGaaO% qaamXvP5wqonvsaeHbfv3ySLgzaGqbdiab-Hha4jaacIcacqWF0baD% cqWFSaalcqaH1oqzcaGGPaGaeyypa0Jae8hEaG3aaSbaaSqaaiaaic% daaeqaaOGaaiikaiabe67a4jaacYcacqaH3oaAcaGGPaGaey4kaSYa% aabmaeaacqaH1oqzdaahaaWcbeqaaiaad6gaaaaabaGaamOBaiabg2% da9iaaigdaaeaacaWGTbGaeyOeI0IaaGymaaqdcqGHris5aOGae8hE% aG3aaSbaaSqaaiab-5gaUbqabaGccaGGOaGaeqOVdGNaaiykaiabgU% caRiaad+eacaGGOaGaeqyTdu2aaWbaaSqabeaacaWGTbaaaOGaaiyk% aaaa!67B9!\[x(t,\varepsilon ) = x_0 (\xi ,\eta ) + \sum\nolimits_{n = 1}^{m - 1} {\varepsilon ^n } x_n (\xi ) + O(\varepsilon ^m )\], where % MathType!MTEF!2!1!+-% feaafiart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXafv3ySLgzGmvETj2BSbqefm0B1jxALjhiov2D% aebbfv3ySLgzGueE0jxyaibaiGc9yrFr0xXdbba91rFfpec8Eeeu0x% Xdbba9frFj0-OqFfea0dXdd9vqaq-JfrVkFHe9pgea0dXdar-Jb9hs% 0dXdbPYxe9vr0-vr0-vqpWqaaeaabiGaciaacaqabeaadaqaaqGaaO% qaaiaabsgacqaH+oaEcaGGVaGaaeizaiaadshacqGH9aqpdaaeWaqa% aiabew7aLnaaCaaaleqabaGaamOBaaaaaeaacaWGUbGaeyypa0JaaG% ymaaqaaiaad2gaa0GaeyyeIuoakiaadkfadaWgaaWcbaGaamOBaaqa% baGccaGGOaGaeqOVdGNaaiykaaaa!4FFC!\[{\text{d}}\xi /{\text{d}}t = \sum\nolimits_{n = 1}^m {\varepsilon ^n } R_n (\xi )\], % MathType!MTEF!2!1!+-% feaafiart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXafv3ySLgzGmvETj2BSbqefm0B1jxALjhiov2D% aebbfv3ySLgzGueE0jxyaibaiGc9yrFr0xXdbba91rFfpec8Eeeu0x% Xdbba9frFj0-OqFfea0dXdd9vqaq-JfrVkFHe9pgea0dXdar-Jb9hs% 0dXdbPYxe9vr0-vr0-vqpWqaaeaabiGaciaacaqabeaadaqaaqGaaO% qaaiaabsgacqaH3oaAcaGGVaGaaeizaiaadshacqGH9aqpdaaeWaqa% aiabew7aLnaaCaaaleqabaGaamOBaaaaaeaacaWGUbGaeyypa0JaaG% imaaqaaiaad2gaa0GaeyyeIuoakiaadofadaWgaaWcbaGaamOBaaqa% baGccaGGOaGaeqOVdGNaaiilaiabeE7aOjaacMcaaaa!5241!\[{\text{d}}\eta /{\text{d}}t = \sum\nolimits_{n = 0}^m {\varepsilon ^n } S_n (\xi ,\eta )\], and R _n,S _nare to be determined in the course of the analysis. This method is suitable for the systems with even non-linearities as well as with odd non-linearities. It can be viewed as a generalization of the two-variable expansion procedure. Using the present method we obtained a modified Krylov-Bogoliubov method. Four numerical examples are presented which served to demonstrate the effectiveness of the present method.     

On the stability or instability of the singular solution of the semilinear heat equation with exponential reaction term

We study questions of existence, uniqueness and asymptotic behaviour for the solutions of u(x, t) of the problem $$\begin{gathered} {\text{ }}u_t - \Delta u = \lambda e^u ,{\text{ }}\lambda {\text{ > 0, }}t > 0,{\text{ }}x{\text{ }}\varepsilon B, \hfill \\ (P){\text{ }}u(x,0) = u_0 (x),{\text{ }}x{\text{ }}\varepsilon B, \hfill \\ {\text{ }}u(x,t) = 0{\text{ }}on{\text{ }}\partial B \times (0,\infty ), \hfill \\ \end{gathered} $$ where B is the unit ball $\{ x\varepsilon R^N :|x|{\text{ }} \leqq {\text{ }}1\} {\text{ and }}N \geqq 3$ . Our interest is focused on the parameter λ ₀=2(N?2) for which (P) admits a singular stationary solution of the form $$S(x) = - 2log|x|$$ . We study the dynamical stability or instability of S, which depends on the dimension. In particular, there exists a minimal bounded stationary solution u which is stable if $3 \leqq N \leqq 9$ , while S is unstable. For $N \geqq 10$ there is no bounded minimal solution and S is an attractor from below but not from above. In fact, solutions larger than S cannot exist in any time interval (there is instantaneous blow-up), and this happens for all dimensions.     

Fluid flow and heat transfer in microchannels with rectangular cross section

Forced convective heat transfer coefficients and friction factors for flow of water in microchannels with a rectangular cross section were measured. An integrated microsystem consisting of five microchannels on one side and a localized heater and seven polysilicon temperature sensors along the selected channels on the other side was fabricated using a double-polished-prime silicon wafer. For the microchannels tested, the friction factor constant obtained are values between 53.7 and 60.4, which are close to the theoretical value from a correlation for macroscopic dimension, 56.9 for D _h  = 100 μm. The heat transfer coefficients obtained by measuring the wall temperature along the micro channels were linearly dependent on the wall temperature, in turn, the heat transfer mechanism is strongly dependent on the fluid properties such as viscosity. The measured Nusselt number in the laminar flow regime tested could be correlated by which is quite different from the constant value obtained in macrochannels.     

A boundary value problem associated with the second painlevé transcendent and the Korteweg-de Vries equation

The differential equation considered is \(y'' - xy = y|y|^\alpha \) . For general positive α this equation arises in plasma physics, in work of De Boer & Ludford. For α=2, it yields similarity solutions to the well-known Korteweg-de Vries equation. Solutions are sought which satisfy the boundary conditions (1) y(∞)=0 (2) (1) $$y{\text{(}}\infty {\text{)}} = {\text{0}}$$ (2) $$y{\text{(}}x{\text{) \~( - }}\tfrac{{\text{1}}}{{\text{2}}}x{\text{)}}^{{{\text{1}} \mathord{\left/ {\vphantom {{\text{1}} \alpha }} \right. \kern-\nulldelimiterspace} \alpha }} {\text{ as }}x \to - \infty $$ It is shown that there is a unique such solution, and that it is, in a certain sense, the boundary between solutions which exist on the whole real line and solutions which, while tending to zero at plus infinity, blow up at a finite x. More precisely, any solution satisfying (1) is asymptotic at plus infinity to some multiple kA i(x) of Airy's function. We show that there is a unique k^*(α) such that when k=k^*(α) the condition (2) is also satisfied. If 0 ^*, the solution exists for all x and tends to zero as x→-∞, while if k>k ^* then the solution blows up at a finite x. For the special case α=2 the differential equation is classical, having been studied by Painlevé around the turn of the century. In this case, using an integral equation derived by inverse scattering techniques by Ablowitz & Segur, we are able to show that k^*=1, confirming previous numerical estimates.     

Uniqueness of Positive Ground State Solutions of the Logarithmic Schrödinger Equation

We prove the uniqueness of positive ground state solutions of the problem \({ {\frac {d^{2}u}{dr^{2}}} + {\frac {n-1}{r}}{\frac {du}{dr}} + u \ln(|u|) = 0}\), \({u(r) > 0~\forall r \ge 0}\), and \({(u(r),u'(r)) \to (0, 0)}\) as \({r \to \infty}\). This equation is derived from the logarithmic Schrödinger equation \({{\rm i}\psi_{t} = {\Delta} \psi + u \ln \left(|u|^{2}\right)}\), and also from the classical equation \({{\frac {\partial u}{\partial t}} = {\Delta} u +u \left(|u|^{p-1}\right) -u}\). For each \({n \ge 1}\), a positive ground state solution is \({ u_{0}(r) = \exp \left(-{\frac{r^2}{4}} + {\frac{n}{2}}\right),~0 \le r < \infty}\). We combine \({u_{0}(r)}\) with energy estimates and associated Ricatti equation estimates to prove that, for each \({n \in \left[1, 9 \right]}\), \({u_{0}(r)}\) is the only positive ground state. We also investigate the stability of \({u_{0}(r)}\). Several open problems are stated.     

Global Well-posedness of Incompressible Inhomogeneous Fluid Systems with Bounded Density or Non-Lipschitz Velocity

In this paper, we first prove the global existence of weak solutions to the d-dimensional incompressible inhomogeneous Navier–Stokes equations with initial data ${a_0 \in L^\infty (\mathbb{R}^d), u_0 = (u_0^h, u_0^d) \in \dot{B}^{-1+\frac{d}{p}}_{p, r} (\mathbb{R}^d)}$ , which satisfy ${(\mu \| a_0 \|_{L^\infty} + \|u_0^h\|_{\dot{B}^{-1+\frac{d}{p}}_{p, r}}) {\rm exp}(C_r{\mu^{-2r}}\|u_0^d\|_{\dot{B}^{-1+\frac{d}{p}}_{p,r}}^{2r}) \leqq c_0\mu}$ for some positive constants c ₀, C _r and 1 < p < d, 1 < r < ∞. The regularity of the initial velocity is critical to the scaling of this system and is general enough to generate non-Lipschitz velocity fields. Furthermore, with additional regularity assumptions on the initial velocity or on the initial density, we can also prove the uniqueness of such a solution. We should mention that the classical maximal L ^p (L ^q ) regularity theorem for the heat kernel plays an essential role in this context.     

Regularization of Point Vortices Pairs for the Euler Equation in Dimension Two

In this paper, we construct stationary classical solutions of the incompressible Euler equation approximating singular stationary solutions of this equation. This procedure is carried out by constructing solutions to the following elliptic problem $$\left\{\begin{array}{l@{\quad}l} -\varepsilon^2 \Delta u = \sum\limits_{i=1}^m \chi_{\Omega_i^{+}} \left(u - q - \frac{\kappa_i^{+}}{2\pi} {\rm ln} \frac{1}{\varepsilon}\right)_+^p\\ \quad - \sum_{j=1}^n \chi_{\Omega_j^{-}} \left(q - \frac{\kappa_j^{-}}{2\pi} {\rm \ln} \frac{1}{\varepsilon} - u\right)_+^p , \quad \quad x \in \Omega,\\ u = 0, \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad x \in \partial \Omega,\end{array}\right.$$ where p > 1, ${\Omega \subset \mathbb{R}^2}$ is a bounded domain, ${\Omega_i^{+}}$ and ${\Omega_j^{-}}$ are mutually disjoint subdomains of Ω and ${\chi_{\Omega_i^{+}} ({\rm resp}.\; \chi_{\Omega_j^{-}})}$ are characteristic functions of ${\Omega_i^{+}({\rm resp}. \;\Omega_j^{-}})$ , q is a harmonic function. We show that if Ω is a simply-connected smooth domain, then for any given C ¹-stable critical point of Kirchhoff–Routh function ${\mathcal{W}\;(x_1^{+},\ldots, x_m^{+}, x_1^{-}, \ldots, x_n^{-})}$ with ${\kappa^{+}_i > 0\,(i = 1,\ldots, m)}$ and ${\kappa^{-}_j > 0\,(j = 1,\ldots,n)}$ , there is a stationary classical solution approximating stationary m + n points vortex solution of incompressible Euler equations with total vorticity ${\sum_{i=1}^m \kappa^{+}_i -\sum_{j=1}^n \kappa_j^{-}}$ . The case that n = 0 can be dealt with in the same way as well by taking each ${\Omega_j^{-}}$ as an empty set and set ${\chi_{\Omega_j^{-}} \equiv 0,\,\kappa^{-}_j=0}$ .     

A note on the soil-water conductivity of a fractal soil

We show that for a fractal soil the soil-water conductivity, K, is given by $$\frac{K}{{K_\varepsilon }} = (\Theta /\varepsilon )^{2D/3 + 2/(3 - D)}$$ where $K_\varepsilon$ is the saturated conductivity, θ the water content, ? its saturated value and D is the fractal dimension obtained from reinterpreting Millington and Quirk's equation for practical values of the porosity ?, as $$D = 2 + 3\frac{{\varepsilon ^{4/3} + (1 - \varepsilon )^{2/3} - 1}}{{2\varepsilon ^{4/3} \ln ,{\text{ }}\varepsilon ^{ - 1} + (1 - \varepsilon )^{2/3} \ln (1 - \varepsilon )^{ - 1} }}$$ .