Solving high-order partial differential equations with indirect radial basis function networks

This paper reports a new numerical method based on radial basis function networks (RBFNs) for solving high-order partial differential equations (PDEs). The variables and their derivatives in the governing equations are represented by integrated RBFNs. The use of integration in constructing neural networks allows the straightforward implementation of multiple boundary conditions and the accurate approximation of high-order derivatives. The proposed RBFN method is verified successfully through the solution of thin-plate bending and viscous flow problems which are governed by biharmonic equations. For thermally driven cavity flows, the solutions are obtained up to a high Rayleigh number

Force dipoles and stable local defects on fluid vesicles

An exact description is provided of an almost spherical fluid vesicle with a fixed area and a fixed enclosed volume locally deformed by external normal forces bringing two nearby points on the surface together symmetrically. The conformal invariance of the two-dimensional bending energy is used to identify the distribution of energy as well as the stress established in the vesicle. While these states are local minima of the energy, this energy is degenerate; there is a zero mode in the energy fluctuation spectrum, associated with area and volume preserving conformal transformations, which breaks the symmetry between the two points. The volume constraint fixes the distance $S$, measured along the surface, between the two points; if it is relaxed, a second zero mode appears, reflecting the independence of the energy on $S$; in the absence of this constraint a pathway opens for the membrane to slip out of the defect. Logarithmic curvature singularities in the surface geometry at the points of contact signal the presence of external forces. The magnitude of these forces varies inversely with $S$ and so diverges as the points merge; the corresponding torques vanish in these defects. The geometry behaves near each of the singularities as a biharmonic monopole, in the region between them as a surface of constant mean curvature, and in distant regions as a biharmonic quadrupole. Comparison of the distribution of stress with the quadratic approximation in the height functions points to shortcomings of the latter representation. Radial tension is accompanied by lateral compression, both near the singularities and far away, with a crossover from tension to compression occurring in the region between them.Comment: 26 pages, 10 figure

Optimising Spatial and Tonal Data for PDE-based Inpainting

Some recent methods for lossy signal and image compression store only a few selected pixels and fill in the missing structures by inpainting with a partial differential equation (PDE). Suitable operators include the Laplacian, the biharmonic operator, and edge-enhancing anisotropic diffusion (EED). The quality of such approaches depends substantially on the selection of the data that is kept. Optimising this data in the domain and codomain gives rise to challenging mathematical problems that shall be addressed in our work. In the 1D case, we prove results that provide insights into the difficulty of this problem, and we give evidence that a splitting into spatial and tonal (i.e. function value) optimisation does hardly deteriorate the results. In the 2D setting, we present generic algorithms that achieve a high reconstruction quality even if the specified data is very sparse. To optimise the spatial data, we use a probabilistic sparsification, followed by a nonlocal pixel exchange that avoids getting trapped in bad local optima. After this spatial optimisation we perform a tonal optimisation that modifies the function values in order to reduce the global reconstruction error. For homogeneous diffusion inpainting, this comes down to a least squares problem for which we prove that it has a unique solution. We demonstrate that it can be found efficiently with a gradient descent approach that is accelerated with fast explicit diffusion (FED) cycles. Our framework allows to specify the desired density of the inpainting mask a priori. Moreover, is more generic than other data optimisation approaches for the sparse inpainting problem, since it can also be extended to nonlinear inpainting operators such as EED. This is exploited to achieve reconstructions with state-of-the-art quality. We also give an extensive literature survey on PDE-based image compression methods

Boundary knot method: A meshless, exponential convergence, integration-free, and boundary-only RBF technique

Based on the radial basis function (RBF), non-singular general solution and dual reciprocity principle (DRM), this paper presents an inheretnly meshless, exponential convergence, integration-free, boundary-only collocation techniques for numerical solution of general partial differential equation systems. The basic ideas behind this methodology are very mathematically simple and generally effective. The RBFs are used in this study to approximate the inhomogeneous terms of system equations in terms of the DRM, while non-singular general solution leads to a boundary-only RBF formulation. The present method is named as the boundary knot method (BKM) to differentiate it from the other numerical techniques. In particular, due to the use of non-singular general solutions rather than singular fundamental solutions, the BKM is different from the method of fundamental solution in that the former does no need to introduce the artificial boundary and results in the symmetric system equations under certain conditions. It is also found that the BKM can solve nonlinear partial differential equations one-step without iteration if only boundary knots are used. The efficiency and utility of this new technique are validated through some typical numerical examples. Some promising developments of the BKM are also discussed.Comment: 36 pages, 2 figures, Welcome to contact me on this paper: Email: [email protected] or [email protected]

A meshless, integration-free, and boundary-only RBF technique

Based on the radial basis function (RBF), non-singular general solution and dual reciprocity method (DRM), this paper presents an inherently meshless, integration-free, boundary-only RBF collocation techniques for numerical solution of various partial differential equation systems. The basic ideas behind this methodology are very mathematically simple. In this study, the RBFs are employed to approximate the inhomogeneous terms via the DRM, while non-singular general solution leads to a boundary-only RBF formulation for homogenous solution. The present scheme is named as the boundary knot method (BKM) to differentiate it from the other numerical techniques. In particular, due to the use of nonsingular general solutions rather than singular fundamental solutions, the BKM is different from the method of fundamental solution in that the former does no require the artificial boundary and results in the symmetric system equations under certain conditions. The efficiency and utility of this new technique are validated through a number of typical numerical examples. Completeness concern of the BKM due to the only use of non-singular part of complete fundamental solution is also discussed

Finite elements modelling of scattering problems for flexural waves in thin plates: Application to elliptic invisibility cloaks, rotators and the mirage effect

We propose a finite elements algorithm to solve a fourth order partial differential equation governing the propagation of time-harmonic bending waves in thin elastic plates. Specially designed perfectly matched layers are implemented to deal with the infinite extent of the plates. These are deduced from a geometric transform in the biharmonic equation. To numerically illustrate the power of elastodynamic transformations, we analyse the elastic response of an elliptic invisibility cloak surrounding a clamped obstacle in the presence of a cylindrical excitation i.e. a concentrated point force. Elliptic cloaking for flexural waves involves a density and an orthotropic Young's modulus which depend on the radial and azimuthal positions, as deduced from a coordinates transformation for circular cloaks in the spirit of Pendry et al. [Science {\bf 312}, 1780 (2006)], but with a further stretch of a coordinate axis. We find that a wave radiated by a concentrated point force located a couple of wavelengths away from the cloak is almost unperturbed in magnitude and in phase. However, when the point force lies within the coating, it seems to radiate from a shifted location. Finally, we emphasize the versatility of transformation elastodynamics with the design of an elliptic cloak which rotates the polarization of a flexural wave within its core.Comment: 14 pages, 5 figure