Solvability and asymptotics of the heat equation with mixed variable lateral conditions and applications in the opening of the exocytotic fusion pore in cells

We investigate a mixed problem with variable lateral conditions for the heat equation that arises in modelling exocytosis, i.e. the opening of a cell boundary in specific biological species for the release of certain molecules to the exterior of the cell. The Dirichlet condition is imposed on a surface patch of the boundary and this patch is occupying a larger part of the boundary as time increases modelling where the cell is opening (the fusion pore), and on the remaining part, a zero Neumann condition is imposed (no molecules can cross this boundary). Uniform concentration is assumed at the initial time. We introduce a weak formulation of this problem and show that there is a unique weak solution. Moreover, we give an asymptotic expansion for the behaviour of the solution near the opening point and for small values in time. We also give an integral equation for the numerical construction of the leading term in this expansion.

An alternating boundary integral based method for a Cauchy problem for the Laplace equation in semi-infinite regions

We consider a Cauchy problem for the Laplace equation in a two-dimensional semi-infinite region with a bounded inclusion, i.e. the region is the intersection between a half-plane and the exterior of a bounded closed curve contained in the half-plane. The Cauchy data are given on the unbounded part of the boundary of the region and the aim is to construct the solution on the boundary of the inclusion. In 1989, Kozlov and Maz'ya [10] proposed an alternating iterative method for solving Cauchy problems for general strongly elliptic and formally self-adjoint systems in bounded domains. We extend their approach to our setting and in each iteration step mixed boundary value problems for the Laplace equation in the semi-infinite region are solved. Well-posedness of these mixed problems are investigated and convergence of the alternating procedure is examined. For the numerical implementation an efficient boundary integral equation method is proposed, based on the indirect variant of the boundary integral equation approach. The mixed problems are reduced to integral equations over the (bounded) boundary of the inclusion. Numerical examples are included showing the feasibility of the proposed method.

An integral equation method for a mixed initial boundary value problem for unsteady Stokes system in a doubly-connected domain

We present a novel numerical method for a mixed initial boundary value problem for the unsteady Stokes system in a planar doubly-connected domain. Using a Laguerre transformation the unsteady problem is reduced to a system of boundary value problems for the Stokes resolvent equations. Employing a modied potential approach we obtain a system of boundary integral equations with various singularities and we use a trigonometric quadrature method for their numerical solution. Numerical examples are presented showing that accurate approximations can be obtained with low computational cost.

On an indirect integral equation approach for stationary heat transfer in semi-infinite layered domains in R3 with cavities

Regions containing internal boundaries such as composite materials arise in many applications.We consider a situation of a layered domain in IR3 containing a nite number of bounded cavities. The model is stationary heat transfer given by the Laplace equation with piecewise constant conductivity. The heat ux (a Neumann condition) is imposed on the bottom of the layered region and various boundary conditions are imposed on the cavities. The usual transmission (interface) conditions are satised at the interface layer, that is continuity of the solution and its normal derivative. To eciently calculate the stationary temperature eld in the semi-innite region, we employ a Green's matrix technique and reduce the problem to boundary integral equations (weakly singular) over the bounded surfaces of the cavities. For the numerical solution of these integral equations, we use Wienert's approach [20]. Assuming that each cavity is homeomorphic with the unit sphere, a fully discrete projection method with super-algebraic convergence order is proposed. A proof of an error estimate for the approximation is given as well. Numerical examples are presented that further highlights the eciency and accuracy of the proposed method.

An alternating potential based approach for a Cauchy problem for the Laplace equation in a planar domain with a cut

We consider a Cauchy problem for the Laplace equation in a bounded region containing a cut, where the region is formed by removing a sufficiently smooth arc (the cut) from a bounded simply connected domain D. The aim is to reconstruct the solution on the cut from the values of the solution and its normal derivative on the boundary of the domain D. We propose an alternating iterative method which involves solving direct mixed problems for the Laplace operator in the same region. These mixed problems have either a Dirichlet or a Neumann boundary condition imposed on the cut and are solved by a potential approach. Each of these mixed problems is reduced to a system of integral equations of the first kind with logarithmic and hypersingular kernels and at most a square root singularity in the densities at the endpoints of the cut. The full discretization of the direct problems is realized by a trigonometric quadrature method which has super-algebraic convergence. The numerical examples presented illustrate the feasibility of the proposed method.

Recovering boundary data in planar heat conduction using boundary integral equation method

We consider a Cauchy problem for the heat equation, where the temperature field is to be reconstructed from the temperature and heat flux given on a part of the boundary of the solution domain. We employ a Landweber type method proposed in [2], where a sequence of mixed well-posed problems are solved at each iteration step to obtain a stable approximation to the original Cauchy problem. We develop an efficient boundary integral equation method for the numerical solution of these mixed problems, based on the method of Rothe. Numerical examples are presented both with exact and noisy data, showing the efficiency and stability of the proposed procedure and approximations.

Computer methods for the heat conduction equation

The merits of various numerical methods for the solution of the one and two dimensional heat conduction equation with a radiation boundary condition have been examined from a practical standpoint in order to determine accuracies and efficiencies. It is found that the use of five increments to approximate the space derivatives gives sufficiently accurate results provided the time step is not too large; further, the implicit backward difference method of Liebmann (27) is found to be the most accurate method. On this basis, a new implicit method is proposed for the solution of the three-dimensional heat conduction equation with radiation boundary conditions. The accuracies of the integral and analogue computer methods are also investigated.

Breakup of a multisoliton state of the linearly damped nonlinear Schrödinger equation

We address the breakup (splitting) of multisoliton solutions of the nonlinear Schrödinger equation (NLSE), occurring due to linear loss. Two different approaches are used for the study of the splitting process. The first one is based on the direct numerical solution of the linearly damped NLSE and the subsequent analysis of the eigenvalue drift for the associated Zakharov-Shabat spectral problem. The second one involves the multisoliton adiabatic perturbation theory applied for studying the evolution of the solution parameters, with the linear loss taken as a small perturbation. We demonstrate that in the case of strong nonadiabatic loss the evolution of the Zakharov-Shabat eigenvalues can be quite nontrivial. We also demonstrate that the multisoliton breakup can be correctly described within the framework of the adiabatic perturbation theory and can take place even due to small linear loss. Eventually we elucidate the occurrence of the splitting and its dependence on the phase mismatch between the solitons forming a two-soliton bound state. © 2007 The American Physical Society.

Nonlinear integral equation methods for the reconstruction of an acoustically sound-soft obstacle

The problem considered is that of determining the shape of a planar acoustically sound-soft obstacle from knowledge of the far-field pattern for one time-harmonic incident field. Two methods, which are based on the solution of a pair of integral equations representing the incoming wave and the far-field pattern, respectively, are proposed and investigated for finding the unknown boundary. Numerical resultsare included which show that the methods give accurate numerical approximations in relatively few iterations.

An iterative method for a Cauchy problem for the heat equation

An iterative method for reconstruction of the solution to a parabolic initial boundary value problem of second order from Cauchy data is presented. The data are given on a part of the boundary. At each iteration step, a series of well-posed mixed boundary value problems are solved for the parabolic operator and its adjoint. The convergence proof of this method in a weighted L2-space is included.

Fokker-Planck equation approach to the description of soliton statistics in optical fiber transmission systems

We derive rigorously the Fokker-Planck equation that governs the statistics of soliton parameters in optical transmission lines in the presence of additive amplifier spontaneous emission. We demonstrate that these statistics are generally non-Gaussian. We present exact marginal probability-density functions for soliton parameters for some cases. A WKB approach is applied to describe the tails of the probability-density functions. © 2005 Optical Society of America.

Strichartz Type Estimates for Oscillatory Problems for Semilinear Wave Equation

The author is partially supported by: M. U. R. S. T. Prog. Nazionale “Problemi e Metodi nella Teoria delle Equazioni Iperboliche”.