Dynamics in dumbbell domains II. The limiting problem

In this work we continue the analysis of the asymptotic dynamics of reaction-diffusion problems in a dumbbell domain started in [J.M. Arrieta, A.N. Carvalho, G. Lozada-Cruz, Dynamics in dumbbell domains I. Continuity of the set of equilibria, J. Differential Equations 231 (2) (2006) 551-597]. Here we study the limiting problem, that is, an evolution problem in a “domain” which consists of an open, bounded and smooth set Ω⊂RN with a curve R₀ attached to it. The evolution in both parts of the domain is governed by a parabolic equation. In Ω the evolution is independent of the evolution in R₀ whereas in R₀ the evolution depends on the evolution in Ω through the continuity condition of the solution at the junction points. We analyze in detail the linear elliptic and parabolic problem, the generation of linear and nonlinear semigroups, the existence and structure of attractors.     

The Hele-Shaw problem with surface tension in a half-plane: A model problem

Consider the Hele-Shaw problem with surface tension in the half-plane {y₁>0} when at time t=0 the domain Ω(t) lies partly on the line y₁=0, and partly in {y₁>0}. In order to establish existence of a solution to this free boundary problem we need to study the (linear) model problem when the Ω(t) is a fixed angular domain. In this paper we consider this model problem and establish existence of a solution satisfying sharp weighted Hölder estimates. These estimates will be used in subsequent work to solve the full Hele-Shaw problem.     

On a method for computing waveguide scattering matrices in the presence of point spectrum

A waveguide occupies a domain G in ? ⁿ⁺¹, n ? 1, having several cylindrical outlets to infinity. The waveguide is described by a general elliptic boundary value problem that is self-adjoint with respect to the Green formula and contains a spectral parameter µ. As an approximation to a row of the scattering matrix S(µ) we suggest a minimizer of a quadratic functional J ^R (·, µ). To construct such a functional, we solve an auxiliary boundary value problem in the bounded domain obtained by cutting off, at a distance R, the waveguide outlets to infinity. It is proved that, if a finite interval [µ₁, µ₂] of the continuous spectrum contains no thresholds, then, as R → ∞, the minimizer tends to the row of the scattering matrix at an exponential rate uniformly with respect to µ ∈ [µ₁, µ₂]. The interval may contain some waveguide eigenvalues whose eigenfunctions exponentially decay at infinity.     

A mixed boundary value problem for a singularly perturbed reaction-diffusion equation in an L-shaped domain

A mixed boundary value problem for a singularly perturbed reaction-diffusion equation in an L-shaped domain is considered for when the solution has singularities at the corners of the domain. The densification of the Shishkin mesh near the inner corner where different boundary conditions meet is such that the solution obtained by the classical five-point difference scheme converges to the solution of the initial problem in the mesh norm L _?? ^h uniformly with respect to the small parameter with almost second order, i.e., as a smooth solution. Numerical analysis confirms the theoretical result.     

Evaluation of a boundary integral representation for the conformal mapping of the unit disk onto a simply-connected domain

A representation of the conformal mapping g of the interior or exterior of the unit circle onto a simply-connected domain Ω as a boundary integral in terms of?|_?Ω is obtained, where? :=g ^-l. A product integration scheme for the approximation of the boundary integral is described and analysed. An ill-conditioning problem related to the domain geometry is discussed. Numerical examples confirm the conclusions of this discussion and support the analysis of the quadrature scheme.     

Large time and space size behaviour of the heat equation in non-cylindrical domains

In this note we consider a linear parabolic problem defined on a non-cylindrical unbounded domain Q. If Ω _t denotes the section of Q above t, the Ω _t size goes to +∞, when t → +∞, i.e. the sections Ω _t become unbounded in some directions when the time t becomes large. Here a model problem is studied, but the technique used can be applied for a wide class of problems, as nonlinear ones, defined on more general domains Q as those introduced by Lions [11]. An asymptotically exponential convergence of the solutions of such problems towards the solution of an elliptic problem defined on a lower dimensional domain is established.     

On Parameter Depending Differential Systems Under Carathéodory Conditions

For every λ in a complex domain G, consider on some interval I the initial value problem y′(λ,x) = A(λ,x)y(λ,x) + b(λ,x), y(λ,x₀) - y₀. If this problem satisfies the Carathéodory conditions for every A, then there exist locally absolutely continuous and almost everywhere differentiable solutions y(λ,· ) of the initial value problem. In general, the union N of the exceptional sets N _λ ? I where y(λ, ·) is not differentiate or does not fulfill the differential equation, is not of Lebesgue measure zero. It will be shown that N is of Lebesgue measure zero provided that A and b are holomorphic with respect to λ and their integrals with respect to x are locally bounded on G × I.     

Nonlinear stability analysis of a regular vortex pentagon outside a circle

A nonlinear stability analysis of the stationary rotation of a system of five identical point vortices lying uniformly on a circle of radius R ₀ outside a circular domain of radius R is performed. The problem is reduced to the problem of stability of an equilibrium position of a Hamiltonian system with a cyclic variable. The stability of stationary motion is interpreted as Routh stability. Conditions for stability, formal stability and instability are obtained depending on the values of the parameter q = R ²/R ₀ ² .     

Parabolic maximal functions and potentials of distributions in Hp

A linear homogeneous ODE is constructed, among whose solutions are all products of solutions of two given linear homogeneous ODE's L_m[u]=0, M_n[v]=0, in some classes. Its order is the minimum and its coefficients can be obtained by a finite number of rational operations and differentiations on the coefficients of L_m, M_n. The problem is considered (locally) both in the real and in the complex domain, around an isolated singularity. Examples are also given.     

A mesh-free numerical method for solution of the family of Kuramoto-Sivashinsky equations

In this paper, radial basis function (RBFs) based mesh-free method is implemented to find numerical solution of the Kuramoto-Sivashinsky equations. This approach has an edge over traditional methods such as finite-difference and finite element methods because it does not require a mesh to discretize the problem domain, and a set of scattered nodes in the domain of influence provided by initial data is required for the realization of the method. The accuracy of the method is assessed in terms of the error norms L₂,L_∞, number of nodes in the domain of influence, free parameter, dependent parameter RBFs and time step length. Numerical experiments demonstrate accuracy and robustness of the method for solving a class of nonlinear partial differential equations.     

Venttsel’ problems in fractal domains

A nonsteady Venttsel’ problem in a fractal domain Ω or in the corresponding prefractal domain Ω _h is studied. Existence, uniqueness, and regularity results for the strict solution, in both cases, are established as well as convergence results of the solutions of the approximating problems in varying Hilbert spaces.     

A note on the regularity criteria for the Navier-Stokes equations

We consider the regularity problem for 3D Navier-Stokes equations in a bounded domain with smooth boundary. A new sufficient condition which guarantees the regularity of weak solutions on the quotient ∇p/(1+|u|^δ₁+|∇u|^δ₂) for the Navier-Stokes equations is established.     

Sobolev Trace Theorem and the Dirichlet Problem in the Unit Disk

We give a direct and elementary proof for the trace theorem in L _p -based Sobolev spaces, when the domain is the unit disk. We also consider the Dirichlet boundary problem for the Laplace equation, where the boundary value is a function in the Besov space. The Poisson kernel enables us to solve this problem in the unit disk more easily than in a general domain.     

Global existence of solutions for semilinear damped wave equation in 2-D exterior domain

We consider a mixed problem of a damped wave equation u_tt−Δu+u_t=|u|^p in the two dimensional exterior domain case. Small global in time solutions can be constructed in the case when the power p on the nonlinear term |u|^p satisfies p∗=2<p<+∞. For this purpose we shall deal with a radially symmetric solution in the exterior domain. A new device developed in Ikehata-Matsuyama (Sci. Math. Japon. 55 (2002) 33) plays an effective role.     

Sequential and parallel domain decomposition methods for a singularly perturbed parabolic convection-diffusion equation

A boundary value problem for a singularly perturbed parabolic convection-diffusion equation is considered in a rectangular domain in x and t; the perturbation parameter ? multiplying the highest derivative takes arbitrary values in the half-open interval (0,1]. For the boundary value problem, we construct a scheme based on the method of lines in x passing through N ₀+1 points of the mesh with respect to t. To solve the problem on a set of intervals, we apply a domain decomposition method (on overlapping subdomains with the overlap width δ), which is a modification of the Schwarz method. For the continual schemes of the decomposition method, we study how sequential and parallel computations, the order of priority in which the subproblems are sequentially solved on the subdomains, and the value of the parameter ? (as well as the values of N ₀, δ) influence the convergence rate of the decomposition scheme (as N ₀ → ∞), and also computational costs for solving the scheme and time required for its solution (unless a prescribed tolerance is achieved). For convection-diffusion equations, in contrast to reaction-diffusion ones, the sequential scheme turns out to be more efficient than the parallel scheme.     

A method for finding coefficients of a quasilinear hyperbolic equation

The inverse problem of finding the coefficients q(s) and p(s) in the equation u _tt = a ² u _xx + q(u)u _t ? p(u)u _x is investigated. As overdetermination required in the inverse setting, two additional conditions are set: a boundary condition and a condition with a fixed value of the timelike variable. An iteration method for solving the inverse problem is proposed based on an equivalent system of integral equations of the second kind. A uniqueness theorem and an existence theorem in a small domain are proved for the inverse problem to substantiate the convergence of the algorithm.     

Existence of weak solutions of a nonlinear parabolic equation in a noncylindrical domain

Consider the mixed boundary value problem ?_tu + L[u] = f with a squareintegrable initial value and with zero boundary values in a domain Q. L[u] is a nonlinear elliptic operator in divergence form, defined on a domain with timedependent boundary. Weak solutions in cylindrical domains are used to construct a weak solution in Q by approximating Q by a system of cylinders. It is shown that this solution is continuously dependent on the initial value.     

Control problems for parabolic equations on controlled domains

Given the one-dimensional heat equation v_t = v_xx on the controlled domain Q(y) = {(t, x); 0 < x < y(t), 0 < t < T} subject to some initial-boundary conditions, we study the problem of optimally selecting y(·) from some admissible class so as to maximize a given payoff of fixed duration. Q(y) is thus a controlled domain. We also study the problem in which the heat equation holds in Q(y, z) = {z(t) < x < y(t), 0 < t < T}; z minimizing, y maximizing, i.e., the differential game. The principle techniques involved are (i) transforming the controlled domain to an uncontrolled domain and then (ii) using the method of lines for parabolic equations to enable us to use known results for control systems governed by ordinary differential equations. Sufficient conditions for existence in an admissible class is given and the method of lines allows numerical techniques to be applied to determine the optimal control in our class.     

Quantization of periodic motions on compact surfaces of constant negative curvature in a magnetic field

We use the semiclassical approach to study the spectral problem for the Schrödinger operator of a charged particle confined to a two-dimensional compact surface of constant negative curvature. We classify modes of classical motion in the integrable domain E < E _cr and obtain a classification of semiclassical solutions as a consequence. We construct a spectral series (spectrum part approximated by semiclassical eigenvalues) corresponding to energies not exceeding the threshold value E _cr; the degeneration multiplicity is computed for each eigenvalue.