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.     

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.     

Asymptotic analysis of linearly elastic shells. III. Justification of Koiter's shell equations

We consider as in Parts I and II a family of linearly elastic shells of thickness 2?, all having the same middle surfaceS=?(?)?R ³, whereω?R ² is a bounded and connected open set with a Lipschitz-continuous boundary, and? ∈ ?³ (?;R ³). The shells are clamped on a portion of their lateral face, whose middle line is?(γ ₀), whereγ ₀ is a portion of?ω withlength γ ₀>0. For all?>0, let $\zeta _i^\varepsilon$ denote the covariant components of the displacement $u_i^\varepsilon g^{i,\varepsilon }$ of the points of the shell, obtained by solving the three-dimensional problem; let $\zeta _i^\varepsilon$ denote the covariant components of the displacement $\zeta _i^\varepsilon$ a ⁱ of the points of the middle surfaceS, obtained by solving the two-dimensional model ofW.T. Koiter, which consists in finding $$\zeta ^\varepsilon = \left( {\zeta _i^\varepsilon } \right) \in V_K (\omega ) = \left\{ {\eta = (\eta _\iota ) \in {\rm H}^1 (\omega ) \times H^1 (\omega ) \times H^2 (\omega ); \eta _i = \partial _v \eta _3 = 0 on \gamma _0 } \right\}$$ such that $$\begin{gathered} \varepsilon \mathop \smallint \limits_\omega a^{\alpha \beta \sigma \tau } \gamma _{\sigma \tau } (\zeta ^\varepsilon )\gamma _{\alpha \beta } (\eta )\sqrt a dy + \frac{{\varepsilon ^3 }}{3} \mathop \smallint \limits_\omega a^{\alpha \beta \sigma \tau } \rho _{\sigma \tau } (\zeta ^\varepsilon )\rho _{\alpha \beta } (\eta )\sqrt a dy \hfill \\ = \mathop \smallint \limits_\omega p^{i,\varepsilon } \eta _i \sqrt a dy for all \eta = (\eta _i ) \in V_K (\omega ), \hfill \\ \end{gathered}$$ where $a^{\alpha \beta \sigma \tau }$ are the components of the two-dimensional elasticity tensor ofS, $\gamma _{\alpha \beta }$ (η) and $\rho _{\alpha \beta }$ (η) are the components of the linearized change of metric and change of curvature tensors ofS, and $p^{i,\varepsilon }$ are the components of the resultant of the applied forces. Under the same assumptions as in Part I, we show that the fields $\frac{1}{{2_\varepsilon }}\smallint _{ - \varepsilon }^\varepsilon u_i^\varepsilon g^{i,\varepsilon } dx_3^\varepsilon$ and $\zeta _i^\varepsilon$ a ⁱ , both defined on the surfaceS, have the same principal part as? → 0, inH ¹ (ω) for the tangential components, and inL ²(ω) for the normal component; under the same assumptions as in Part II, we show that the same fields again have the same principal part as? → 0, inH ¹ (ω) for all their components. For “membrane” and “flexural” shells, the two-dimensional model ofW.T. Koiter is therefore justified.     

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} }}$$ .     

Linear control semigroups acting on projective space

Linear control semigroupsL?=Gl(d,R) are associated with semilinear control systems of the form whereA:R ^m → gl(d,R) is continuous in some open set containingU. The semigroupL then corresponds to the solutions with piecewise constant controls, i.e., L acts in a natural way onR ^d {0}, on the sphereS ^d?1, and on the projective spaceP ^d?1. Under the assumption that the group generated byL in Gl(d,R) acts transitively onP ^d?1, we analyze the control structure of the action ofL onP ^d?1: We characterize the sets inP ^d?1, where the system is controllable (the control sets) using perturbation theory of eigenvalues and (generalized) eigenspaces of the matrices g εL For nonlinear control systems on finitedimensional manifoldsM, we study the linearization on the tangent bundleTM and the projective bundleP M via the theory of Morse decompositions, to obtain a characterization of the chain-recurrent components of the control flow onU×PM. These components correspond uniquely to the chain control sets onP M, and they induce a subbundle decomposition ofU×TM. These results are used to characterize the chain control sets ofL acting onP ^d?1 and to compare the control sets and chain control sets.     

Adjoint Methods for the Infinity Laplacian Partial Differential Equation

To study fine properties of certain smooth approximations ${u^\varepsilon}$ to a viscosity solution u of the infinity Laplacian partial differential equation (PDE), we introduce Green??s function ${\sigma^\varepsilon}$ for the linearization. We can then integrate by parts with respect to ${\sigma^\varepsilon}$ and derive various useful integral estimates. We are, in particular, able to use these estimates (i) to prove the everywhere differentiability of u and (ii) to rigorously justify interpreting the infinity Laplacian equation as a parabolic PDE.     

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}$$     

On Saint-Venant’s Problem with Concentrated Loads

The classical Saint-Venant problem is to find a solution of the traction problem of elastostatics in a finite cylinder ?? loaded over its bases. We prove that the problem has a unique solution for equilibrated surface forces $\hat{ \boldsymbol { s}}\in W^{-1,q}(\partial\Omega)$ , with q??(2?? ₀,+??) for some positive ? ₀ depending on ??. Hence $\hat{ \boldsymbol { s}}$ can model force acting on ???, concentrated on sets of zero Lebesgue surface measure of ???. Moreover, if $\hat{ \boldsymbol { s}}$ is equilibrated on each basis, we give a simple proof of the Toupin estimate expressing Saint-Venant??s principle.     

An experimental study on turbulence modification in the near-wall boundary layer of a dilute gas-particle channel flow

Turbulence modifications of a dilute gas-particle flow are experimentally investigated in the lower boundary layer of a horizontal channel by means of a simultaneous two-phase PIV measurement technique. The measurements are conducted in the near-wall region with y ⁺?<?250 at Re _τ (based on the wall friction velocity u _τ and half channel height h)?=?430. High spatial resolution and small interrogation window are used to minimize the PIV measurement uncertainty due to the velocity gradient near the wall. Polythene beads with the diameter of 60?μm (d _p ⁺ ?=?1.71, normalized by the fluid kinematic viscosity ν and u _τ) are used as dispersed phase, and three low mass loading ratios (Φ _m ) ranging from 10^?4 to 10^?3 are tested. It is found that the addition of the particles noticeably modifies the mean velocity and turbulent intensities of the gas-phase, as well as the turbulence coherent structures, even at Φ _m ?=?0.025?%. Particle inertia changes the viscous sublayer of the gas turbulence with a smaller thickness and a larger streamwise velocity gradient, which increases the peak value of the streamwise fluctuation velocity ( $ u_{\text{rms}}^{ + } $ ) of the gas-phase with its location shifting to the wall. Particle sedimentation increases the roughness of the bottom wall, which significantly increases the wall-normal fluctuation velocity ( $ v_{\text{rms}}^{ + } $ ) and Reynolds shear stress ( $ - \langle u^{ \prime } v^{\prime } \rangle^{ + } $ ) of the gas-phase in the inner region of the boundary layer (y ⁺?<?10). Under effect of particle–wall collision, the Q2 events (ejections) of the gas-phase are slightly increased by particles, while the Q4 events (sweeps) are obviously decreased. The spatial scale of the coherent structures near the wall shrinks remarkably with the presence of the particles, which may be attributed to the intensified crossing-trajectory effects due to particle saltation near the bottom wall. Meanwhile, the $ v_{\text{rms}}^{ + } $ and $ - \langle u^{ \prime } v^{\prime } \rangle^{ + } $ of the gas-phase are significantly reduced in the outer region of the boundary layer (y ⁺?>?20).     

Universal Moduli of Continuity for Solutions to Fully Nonlinear Elliptic Equations

This paper provides universal, optimal moduli of continuity for viscosity solutions to fully nonlinear elliptic equations F(X, D ² u) =  f(X), based on the weakest and borderline integrability properties of the source function f in different scenarios. The primary result established in this work is a sharp Log-Lipschitz estimate on u based on the L ⁿ norm of f, which corresponds to optimal regularity bounds for the critical threshold case. Optimal C ^1,α regularity estimates are also delivered when ${f\in L^{n+\varepsilon}}$ . The limiting upper borderline case, ${f \in L^\infty}$ , also has transcendental importance to elliptic regularity theory and its applications. In this paper we show, under the convexity assumption on F, that ${u \in C^{1,{\rm Log-Lip}}}$ , provided f has bounded mean oscillation. Once more, such an estimate is optimal. For the lower borderline integrability condition allowed by the theory, we establish interior a priori estimates on the ${C^{0,\frac{n-2\varepsilon}{n-\varepsilon}}}$ norm of u based on the L ^n-ε norm of f, where ? is the Escauriaza universal constant. The exponent ${\frac{n-2\varepsilon}{n-\varepsilon}}$ is optimal. When the source function f lies in L ^q , n > q > n?ε, we also obtain the exact, improved sharp Hölder exponent of continuity.     

Regularity for Energy-Minimizing Area-Preserving Deformations

In this paper we establish the square integrability of the nonnegative hydrostatic pressure p, that emerges in the minimization problem $$\inf_{\mathcal{K}}\int_{\varOmega}|\nabla \textbf {v}|^2, \quad\varOmega\subset \mathbb {R}^2 $$ as the Lagrange multiplier corresponding to the incompressibility constraint det?v=1 a.e. in Ω. Our method employs the Euler-Lagrange equation for the mollified Cauchy stress C satisfied in the image domain Ω ^?=u(Ω). This allows to construct a convex function ψ, defined in the image domain, such that the measure of the normal mapping of ψ controls the L ² norm of the pressure. As a by-product we conclude that $\textbf {u}\in C^{\frac{1}{2}}_{\textrm {loc}}(\varOmega)$ if the dual pressure (introduced in Karakhanyan, Manuscr. Math. 138:463, 2012) is nonnegative.     

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.     

Extensions of Serrin’s Uniqueness and Regularity Conditions for the Navier–Stokes Equations

Consider a smooth bounded domain ${\Omega \subseteq {\mathbb{R}}^3}$ , a time interval [0, T), 0?<?T?≤?∞, and a weak solution u of the Navier–Stokes system. Our aim is to develop several new sufficient conditions on u yielding uniqueness and/or regularity. Based on semigroup properties of the Stokes operator we obtain that the local left-hand Serrin condition for each ${t\in (0,T)}$ is sufficient for the regularity of u. Somehow optimal conditions are obtained in terms of Besov spaces. In particular we obtain such properties under the limiting Serrin condition ${u \in L_{\rm loc}^\infty([0,T);L^3(\Omega))}$ . The complete regularity under this condition has been shown recently for bounded domains using some additional assumptions in particular on the pressure. Our result avoids such assumptions but yields global uniqueness and the right-hand regularity at each time when ${u \in L_{\rm loc}^\infty([0,T);L^3(\Omega))}$ or when ${u(t)\in L^3(\Omega)}$ pointwise and u satisfies the energy equality. In the last section we obtain uniqueness and right-hand regularity for completely general domains.     

Well-Posedness of Nematic Liquid Crystal Flow in {L^{3}_{\rm uloc}(\mathbb{R}^3)}

In this paper, we establish the local well-posedness for the Cauchy problem of a simplified version of hydrodynamic flow of nematic liquid crystals in ${\mathbb{R}^3}$ for any initial data (u ₀, d ₀) having small ${L^{3}_{\rm uloc}}$ -norm of ${(u_{0}, \nabla d_{0})}$ . Here ${L^{3}_{\rm uloc}(\mathbb{R}^3)}$ is the space of uniformly locally L ³-integrable functions. For any initial data (u ₀, d ₀) with small ${\|(u_0, \nabla d_0)\|_{L^{3}(\mathbb{R}^3)}}$ , we show that there exists a unique, global solution to the problem under consideration which is smooth for t > 0 and has monotone deceasing L ³-energy for ${t \geqq 0}$ .     

Hadamard Variational Formula for the Green’s Function of the Boundary Value Problem on the Stokes Equations

For every ${\varepsilon > 0}$ , we consider the Green’s matrix ${G_{\varepsilon}(x, y)}$ of the Stokes equations describing the motion of incompressible fluids in a bounded domain ${\Omega_{\varepsilon} \subset \mathbb{R}^d}$ , which is a family of perturbation of domains from ${\Omega\equiv \Omega_0}$ with the smooth boundary ${\partial\Omega}$ . Assuming the volume preserving property, that is, ${\mbox{vol.}\Omega_{\varepsilon} = \mbox{vol.}\Omega}$ for all ${\varepsilon > 0}$ , we give an explicit representation formula for ${\delta G(x, y) \equiv \lim_{\varepsilon\to +0}\varepsilon^{-1}(G_{\varepsilon}(x, y) - G_0(x, y))}$ in terms of the boundary integral on ${\partial \Omega}$ of ${G_0(x, y)}$ . Our result may be regarded as a classical Hadamard variational formula for the Green’s functions of the elliptic boundary value problems.     

Decay estimates for some semilinear damped hyperbolic problems

Let Ω be a bounded open domain in R ⁿ , g ∶ R → R a non-decreasing continuous function such that g(0)=0 and h ε L _loc ¹ (R⁺; L ²(Ω)). Under suitable assumptions on g and h, the rate of decay of the difference of two solutions is studied for some abstract evolution equations of the general form u ^′′ + Lu + g(u ^′) = h(t,x) as t → + ∞. The results, obtained by use of differential inequalities, can be applied to the case of the semilinear wave equation $$u_u - \Delta u + g{\text{(}}u_t {\text{) = }}h{\text{ in }}R^ + \times \Omega ,{\text{ }}u = {\text{0 on }}R^ + \times \partial \Omega$$ in R ⁺×Ω, u=0 on R ⁺×?Ω. For instance if \(g(s) = c\left| s \right|^{p - 1} s + d\left| s \right|^{q - 1} s\) with c, d>0 and 1 < p≦q, (n?2)q≦n+2, then if \(h \in L^\infty (R + ;L^2 (\Omega ))\) , all solutions are bounded in the energy space for t≧0 and if u, v are two such solutions, the energy norm of u(t) ? v(t) decays like t ^?1/p?1 as t → + ∞.     

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).     

Stochastic Variational Inequalities and Applications to the Total Variation Flow Perturbed by Linear Multiplicative Noise

In this work, we introduce a new method to prove the existence and uniqueness of a variational solution to the stochastic nonlinear diffusion equation ${{\rm d}X(t) = {\rm div} \left[\frac{\nabla X(t)}{|\nabla X(t)|}\right]{\rm d}t + X(t){\rm d}W(t) {\rm in} (0, \infty) \times \mathcal{O},}$ where ${\mathcal{O}}$ is a bounded and open domain in ${\mathbb{R}^N, N \geqq 1}$ and W(t) is a Wiener process of the form ${W(t) = \sum^{\infty}_{k = 1}\mu_{k}e_{k}\beta_{k}(t), e_{k} \in C^{2}(\overline{\mathcal{O}}) \cap H^{1}_{0}(\mathcal{O}),}$ and ${\beta_{k}, k \in \mathbb{N}}$ are independent Brownian motions. This is a stochastic diffusion equation with a highly singular diffusivity term. One main result established here is that for all initial conditions in ${L^2(\mathcal{O})}$ , it is well posed in a class of continuous solutions to the corresponding stochastic variational inequality. Thus, one obtains a stochastic version of the (minimal) total variation flow. The new approach developed here also allows us to prove the finite time extinction of solutions in dimensions ${1\leqq N \leqq3}$ , which is another main result of this work.     

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}$ .