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

Desingularization of Vortex Rings and Shallow Water Vortices by a Semilinear Elliptic Problem

Steady vortices for the three-dimensional Euler equation for inviscid incompressible flows and for the shallow water equation are constructed and shown to tend asymptotically to singular vortex filaments. The construction is based on a study of solutions to the semilinear elliptic problem $$ \left\{ \begin{aligned} -{\rm div} \left(\frac{\nabla u_{\varepsilon}}{b}\right) & = \frac{1}{\varepsilon^2} b f \left(u_{\varepsilon} - \log \tfrac{1}{\varepsilon} q \right) & & \text{ in } \; \Omega, \\u_\varepsilon & = 0 & & \text{ on } \; \partial \Omega, \end{aligned}\right.$$ for small values of ${\varepsilon > 0}$ .     

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.     

Monotonicity and 1-Dimensional Symmetry for Solutions of an Elliptic System Arising in Bose–Einstein Condensation

We study monotonicity and 1-dimensional symmetry for positive solutions with algebraic growth of the following elliptic system: $$\left\{\begin{array}{ll} -\Delta u = -u \upsilon^2 &\quad {\rm in}\, \mathbb{R}^N\\ -\Delta \upsilon= -u^2 \upsilon &\quad {{\rm in}\, \mathbb{R}^N},\end{array}\right.$$ for every dimension ${N \geqq 2}$ . In particular, we prove a Gibbons-type conjecture proposed by Berestycki et al.     

On the effect of non-smooth Coulomb damping on flutter-type self-excitation in a non-gyroscopic circulatory 2-DoF-system

This article deals with self-excited vibrations, attractivity of stationary solutions, and the corresponding bifurcation behavior of two-dimensional differential inclusions of the type $\mathbf{M}\mathbf{q}'' + \mathbf{D}\mathbf{q}' + (\mathbf{K} + \bar{\mu}\mathbf{N})\mathbf{q} \in-\mathbf{R}\operatorname{Sign}(\mathbf{q}')$ . For the smooth case R=0, the equilibrium may become unstable due to non-conservative positional forces stemming from the circulatory matrix N. This type of instability is usually referred to as flutter instability and the loss of stability is related to a Hopf bifurcation of the steady state, which occurs for a critical parameter $\bar{\mu}= \bar{\mu}_{\mathrm{crit}}$ . For R≠0, the steady state is a set of equilibria, which turns out to be attractive for all values of the bifurcation parameter $\bar{\mu}$ . Depending on $\bar{\mu}$ , the basin of attraction of the equilibrium set can be infinite or finite. The transition from an infinite to a finite basin of attraction occurs at the stability threshold $\bar{\mu}_{\mathrm{crit}}$ of the underlying smooth problem. For the finite basin of attraction, its size is proportional to the Coulomb friction and inverse-proportional to $(\bar{\mu}- \bar{\mu}_{\mathrm{crit}})$ . By adding Coulomb damping the notion of steady state stability for the smooth problem is replaced by the question whether the basin of attraction of the steady state is infinite or finite. Simultaneously, the local Hopf-bifurcation is replaced by a global bifurcation. This implies that in the presence of Coulomb damping the occurrence of self-excited vibrations can only be investigated with regard to the perturbation level.     

Stabilization in a two-dimensional chemotaxis-Navier–Stokes system

This paper deals with an initial-boundary value problem for the system $$\left\{ \begin{array}{llll} n_t + u\cdot\nabla n &=& \Delta n -\nabla \cdot (n\chi(c)\nabla c), \quad\quad & x\in\Omega, \, t > 0,\\ c_t + u\cdot\nabla c &=& \Delta c-nf(c), \quad\quad & x\in\Omega, \, t > 0,\\ u_t + \kappa (u\cdot \nabla) u &=& \Delta u + \nabla P + n \nabla\phi, \qquad & x\in\Omega, \, t > 0,\\ \nabla \cdot u &=& 0, \qquad & x\in\Omega, \, t > 0,\end{array} \right.$$ which has been proposed as a model for the spatio-temporal evolution of populations of swimming aerobic bacteria. It is known that in bounded convex domains ${\Omega \subset \mathbb{R}^2}$ and under appropriate assumptions on the parameter functions χ, f and ?, for each ${\kappa\in\mathbb{R}}$ and all sufficiently smooth initial data this problem possesses a unique global-in-time classical solution. The present work asserts that this solution stabilizes to the spatially uniform equilibrium ${(\overline{n_0},0,0)}$ , where ${\overline{n_0}:=\frac{1}{|\Omega|} \int_\Omega n(x,0)\,{\rm d}x}$ , in the sense that as t→∞, $$n(\cdot,t) \to \overline{n_0}, \qquad c(\cdot,t) \to 0 \qquad \text{and}\qquad u(\cdot,t) \to 0$$ hold with respect to the norm in ${L^\infty(\Omega)}$ .     

Global Compensated Compactness Theorem for General Differential Operators of First Order

Let A ₁(x, D) and A ₂(x, D) be differential operators of the first order acting on l-vector functions ${u= (u_1, \ldots, u_l)}$ in a bounded domain ${\Omega \subset \mathbb{R}^{n}}$ with the smooth boundary ${\partial\Omega}$ . We assume that the H ¹-norm ${\|u\|_{H^{1}(\Omega)}}$ is equivalent to ${\sum_{i=1}^2\|A_iu\|_{L^2(\Omega)} + \|B_1u\|_{H^{\frac{1}{2}}(\partial\Omega)}}$ and ${\sum_{i=1}^2\|A_iu\|_{L^2(\Omega)} + \|B_2u\|_{H^{\frac{1}{2}}(\partial\Omega)}}$ , where B _i  = B _i (x, ν) is the trace operator onto ${\partial\Omega}$ associated with A _i (x, D) for i = 1, 2 which is determined by the Stokes integral formula (ν: unit outer normal to ${\partial\Omega}$ ). Furthermore, we impose on A ₁ and A ₂ a cancellation property such as ${A_1A_2^{\prime}=0}$ and ${A_2A_1^{\prime}=0}$ , where ${A^{\prime}_i}$ is the formal adjoint differential operator of A _i (i = 1, 2). Suppose that ${\{u_m\}_{m=1}^{\infty}}$ and ${\{v_m\}_{m=1}^{\infty}}$ converge to u and v weakly in ${L^2(\Omega)}$ , respectively. Assume also that ${\{A_{1}u_m\}_{m=1}^{\infty}}$ and ${\{A_{2}v_{m}\}_{m=1}^{\infty}}$ are bounded in ${L^{2}(\Omega)}$ . If either ${\{B_{1}u_m\}_{m=1}^{\infty}}$ or ${\{B_{2}v_m\}_{m=1}^{\infty}}$ is bounded in ${H^{\frac{1}{2}}(\partial\Omega)}$ , then it holds that ${\int_{\Omega}u_m\cdot v_m \,{\rm d}x \to \int_{\Omega}u\cdot v \,{\rm d}x}$ . We also discuss a corresponding result on compact Riemannian manifolds with boundary.     

Segregated and Synchronized Vector Solutions for Nonlinear Schrödinger Systems

We consider the following nonlinear Schrödinger system in ${\mathbb{R}^3}$ $$\left\{\begin{array}{ll}-\Delta u + P(|x|)u = \mu u^{2}u + \beta v^2u,\quad x \in \mathbb{R}^3,\\-\Delta v + Q(|x|)v = \nu v^{2}v + \beta u^2v,\quad x \in \mathbb{R}^3,\end{array}\right.$$ where P(r) and Q(r) are positive radial potentials, ${\mu > 0, \nu > 0}$ and ${\beta \in \mathbb{R}}$ is a coupling constant. This type of system arises, in particular, in models in Bose–Einstein condensates theory. We examine the effect of nonlinear coupling on the solution structure. In the repulsive case, we construct an unbounded sequence of non-radial positive vector solutions of segregated type, and in the attractive case we construct an unbounded sequence of non-radial positive vector solutions of synchronized type. Depending upon the system being repulsive or attractive, our results exhibit distinct characteristic features of vector solutions.     

Generalized Flows Satisfying Spatial Boundary Conditions

In a region D in ${\mathbb{R}^2}$ or ${\mathbb{R}^3}$ , the classical Euler equation for the regular motion of an inviscid and incompressible fluid of constant density is given by $$\partial_t v+(v\cdot \nabla_x)v=-\nabla_x p, {\rm div}_x v=0,$$ where v(t, x) is the velocity of the particle located at ${x\in D}$ at time t and ${p(t,x)\in\mathbb{R}}$ is the pressure. Solutions v and p to the Euler equation can be obtained by solving $$\left\{\begin{array}{l} \nabla_x\left\{\partial_t\phi(t,x,a) + p(t,x)+(1/2)|\nabla_x\phi(t,x,a)|^2 \right\}=0\,{\rm at}\,a=\kappa(t,x),\\ v(t,x)=\nabla_x \phi(t,x,a)\,{\rm at}\,a=\kappa(t,x), \\ \partial_t\kappa(t,x)+(v\cdot\nabla_x)\kappa(t,x)=0, \\ {\rm div}_x v(t,x)=0, \end{array}\right. \quad\quad\quad\quad\quad(0.1)$$ where $$\phi:\mathbb{R}\times D\times \mathbb{R}^l\rightarrow\mathbb{R}\,{\rm and}\, \kappa:\mathbb{R}\times D \rightarrow \mathbb{R}^l$$ are additional unknown mappings (l?≥ 1 is prescribed). The third equation in the system says that ${\kappa\in\mathbb{R}^l}$ is convected by the flow and the second one that ${\phi}$ can be interpreted as some kind of velocity potential. However vorticity is not precluded thanks to the dependence on a. With the additional condition κ(0, x)?=?x on D (and thus l?=?2 or 3), this formulation was developed by Brenier (Commun Pure Appl Math 52:411–452, 1999) in his Eulerian–Lagrangian variational approach to the Euler equation. He considered generalized flows that do not cross ${\partial D}$ and that carry each “particle” at time t?=?0 at a prescribed location at time t?=?T?>?0, that is, κ(T, x) is prescribed in D for all ${x\in D}$ . We are concerned with flows that are periodic in time and with prescribed flux through each point of the boundary ${\partial D}$ of the bounded region D (a two- or three-dimensional straight pipe). More precisely, the boundary condition is on the flux through ${\partial D}$ of particles labelled by each value of κ at each point of ${\partial D}$ . One of the main novelties is the introduction of a prescribed “generalized” Bernoulli’s function ${H:\mathbb{R}^l\rightarrow \mathbb{R}}$ , namely, we add to (0.1) the requirement that $$\partial_t\phi(t,x,a) +p(t,x)+(1/2)|\nabla_x\phi(t,x,a)|^2=H(a)\,{\rm at}\,a=\kappa(t,x)\quad\quad\quad\quad\quad(0.2)$$ with ${\phi,p,\kappa}$ periodic in time of prescribed period T?>?0. Equations (0.1) and (0.2) have a geometrical interpretation that is related to the notions of “Lamb’s surfaces” and “isotropic manifolds” in symplectic geometry. They may lead to flows with vorticity. An important advantage of Brenier’s formulation and its present adaptation consists in the fact that, under natural hypotheses, a solution in some weak sense always exists (if the boundary conditions are not contradictory). It is found by considering the functional $$(\kappa,v)\rightarrow \int\limits_{0}^T \int\limits_D\left\{\frac 1 2 |v(t,x)|^2+H(\kappa(t,x))\right\}dt\, dx$$ defined for κ and v that are T-periodic in t, such that $$\partial_t\kappa(t,x)+(v\cdot\nabla_x)\kappa(t,x)=0, {\rm div}_x v(t,x)=0,$$ and such that they satisfy the boundary conditions. The domain of this functional is enlarged to some set of vector measures and then a minimizer can be obtained. For stationary planar flows, the approach is compared with the following standard minimization method: to minimize $$\int\limits_{]0,L[\times]0,1[} \{(1/2)|\nabla \psi|^2+H(\psi)\}dx\,{\rm for}\,\psi\in W^{1,2}(]0,L[\times]0,1[)$$ under appropriate boundary conditions, where ψ is the stream function. For a minimizer, corresponding functions ${\phi}$ and κ are given in terms of the stream function ψ.     

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.     

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.     

Regularity and Asymptotic Behavior of Nonlinear Stefan Problems

We study the following nonlinear Stefan problem $$\left\{\begin{aligned}\!\!&u_t\,-\,d\Delta u = g(u) & &\quad{\rm for}\,x\,\in\,\Omega(t), t > 0, \\ & u = 0 \, {\rm and} u_t = \mu|\nabla_{x} u|^{2} &&\quad {\rm for}\,x\,\in\,\Gamma(t), t > 0, \\ &u(0, x) = u_{0}(x) &&\quad {\rm for}\,x\,\in\,\Omega_0,\end{aligned} \right.$$ where ${\Omega(t) \subset \mathbb{R}^{n}}$ ( ${n \geqq 2}$ ) is bounded by the free boundary ${\Gamma(t)}$ , with ${\Omega(0) = \Omega_0}$ , μ and d are given positive constants. The initial function u ₀ is positive in ${\Omega_0}$ and vanishes on ${\partial \Omega_0}$ . The class of nonlinear functions g(u) includes the standard monostable, bistable and combustion type nonlinearities. We show that the free boundary ${\Gamma(t)}$ is smooth outside the closed convex hull of ${\Omega_0}$ , and as ${t \to \infty}$ , either ${\Omega(t)}$ expands to the entire ${\mathbb{R}^n}$ , or it stays bounded. Moreover, in the former case, ${\Gamma(t)}$ converges to the unit sphere when normalized, and in the latter case, ${u \to 0}$ uniformly. When ${g(u) = au - bu^2}$ , we further prove that in the case ${\Omega(t)}$ expands to ${{\mathbb R}^n}$ , ${u \to a/b}$ as ${t \to \infty}$ , and the spreading speed of the free boundary converges to a positive constant; moreover, there exists ${\mu^* \geqq 0}$ such that ${\Omega(t)}$ expands to ${{\mathbb{R}}^n}$ exactly when ${\mu > \mu^*}$ .     

C 1 Solutions for Semi-Implicit Systems of Differential Equations

The present note is a continuation of the author??s effort to study the existence of continuously differentiable solutions to the semi-implicit system of differential equations (1) $$f(x^{\prime}(t)) = g(t, x(t))$$ (2) $$\quad x(0) = x_0,$$ where

${\quad\Omega_g \subseteq \mathbb{R} \times\mathbb{R}^n}$ is an open set containing (0, x ₀) and ${g:\Omega_g \rightarrow\mathbb{R}^n}$ is a continuous function,

${\quad\Omega_f \subseteq \mathbb{R}^n}$ is an open set and ${f:\Omega_f\rightarrow\mathbb{R}^n}$ is a continuous function.

The transformation of (1)?C(2) into a solvable explicit system of differential equations is trivial if f is locally injective around an element ${\gamma\in \Omega_f\cap f^{-1}(g(0,x_0))}$ . In this paper, we study (1)?C(2) when such a translation is not possible because of the inherent multivalued nature of f ^?1.     

KdV Limit of the Euler–Poisson System

Consider the scaling ${\varepsilon^{1/2}(x-Vt) \to x, \varepsilon^{3/2}t \to t}$ in the Euler–Poisson system for ion-acoustic waves (1). We establish that as ${\varepsilon \to 0}$ , the solutions to such Euler–Poisson systems converge globally in time to the solutions of the Korteweg–de Vries equation.     

Real Solutions to the Nonlinear Helmholtz Equation with Local Nonlinearity

In this paper, we study real solutions of the nonlinear Helmholtz equation $$- \Delta u - k^2 u = f(x,u),\quad x\in \mathbb{R}^N$$ satisfying the asymptotic conditions $$u(x)=O\left(|x|^{\frac{1-N}{2}}\right) \quad {\rm and} \quad \frac{\partial^2 u}{\partial r^2}(x)+k^2u(x)=o\left(|x|^{\frac{1-N}{2}}\right) \quad {\rm as}\, r=|x| \to\infty.$$ We develop the variational framework to prove the existence of nontrivial solutions for compactly supported nonlinearities without any symmetry assumptions. In addition, we consider the radial case in which, for a larger class of nonlinearities, infinitely many solutions are shown to exist. Our results give rise to the existence of standing wave solutions of corresponding nonlinear Klein–Gordon equations with arbitrarily large frequency.     

Multiple Blow-Up Phenomena for the Sinh-Poisson Equation

We consider the sinh-Poisson equation $$(P) _ \lambda - \Delta{u} = \lambda \, {\rm sinh} \, u \quad {\rm in} \, \Omega, \quad u = 0 \quad {\rm on} \, \partial\Omega$$ , where Ω is a smooth bounded domain in ${\mathbb{R}^2}$ and λ is a small positive parameter. If ${0 \in \Omega}$ and Ω is symmetric with respect to the origin, for any integer k if λ is small enough, we construct a family of solutions to (P) _λ , which blows up at the origin, whose positive mass is 4πk(k?1) and negative mass is 4πk(k + 1). This gives a complete answer to an open problem formulated by Jost et al. (Calc Var PDE 31(2):263–276, 2008).     

On the Stokes and Oseen Problems with Singular Data

We consider the steady Stokes and Oseen problems in bounded and exterior domains of ${\mathbb{R}^n}$ of class C ^k-1,1 (n = 2, 3; k ≥ 2). We prove existence and uniqueness of a very weak solution for boundary data a in ${W^{2-k-1/q,q} (\partial\Omega)}$ . If ${\Omega}$ is of class ${C^\infty}$ , we can assume a to be a distribution on ${\partial\Omega}$ .