Solving fractional boundary value problems with Dirichlet boundary conditions using a new iterative method

Linear/nonlinear fractional diffusion-wave equations on finite domains with Dirichlet boundary conditions have been solved using a new iterative method proposed by Daftardar-Gejji and Jafari [V. Daftardar-Gejji, H. Jafari, An iterative method for solving nonlinear functional equations, J. Math. Anal. Appl. 316 (2006) 753–763].     

Quadrature imposition of compatibility conditions in Chebyshev methods

Often, in solving an elliptic equation with Neumann boundary conditions, a compatibility condition has to be imposed for well-posedness. This condition involves integrals of the forcing function. When pseudospectral Chebyshev methods are used to discretize the partial differential equation, these integrals have to be approximated by an appropriate quadrature formula. The Gauss-Chebyshev (or any variant of it, like the Gauss-Lobatto) formula cannot be used here since the integrals under consideration do not include the weight function. A natural candidate to be used in approximating the integrals is the Clenshaw-Curtis formula; however, we show in this article that this is the wrong choice and it may lead to divergence if time-dependent methods are used to march the solution to steady state. We develop, in this paper, the correct quadrature formula for these problems. This formula takes into account the degree of the polynomials involved. We show that this formula leads to a well-conditioned Chebyshev approximation to the differential equations and that the compatibility condition is automatically satisfied.     

Deterministic learning of completely resonant nonlinear wave systems with Dirichlet boundary conditions

In this paper, we investigate the approximation of completely resonant nonlinear wave systems via deterministic learning. The plants are distributed parameter systems (DPS) describing homogeneous and isotropic elastic vibrating strings with fixed endpoints. The purpose of the paper is to approximate the infinite-dimensional dynamics, rather than the parameters of the wave systems. To solve the problem, the wave systems are first transformed into finite-dimensional dynamical systems described by ordinary differential equation (ODE). The properties of the finite-dimensional systems, including the convergence of the solution, as well as the dominance of partial system dynamics according to point-wise measurements, are analyzed. Based on the properties, second, by using the deterministic learning algorithm, an approximately accurate neural network (NN) approximation of the the finite-dimensional system dynamics is achieved in a local region along the recurrent trajectories. Simulation studies are included to demonstrate the effectiveness of the proposed approach.     

Solving PDEs in irregular geometries with multiresolution methods I: Embedded Dirichlet boundary conditions

In this work, we develop and analyze a formalism for solving boundary value problems in arbitrarily-shaped domains using the MADNESS (multiresolution adaptive numerical environment for scientific simulation) package for adaptive computation with multiresolution algorithms. We begin by implementing a previously-reported diffuse domain approximation for embedding the domain of interest into a larger domain (Li et al., 2009 [1]). Numerical and analytical tests both demonstrate that this approximation yields non-physical solutions with zero first and second derivatives at the boundary. This excessive smoothness leads to large numerical cancellation and confounds the dynamically-adaptive, multiresolution algorithms inside MADNESS. We thus generalize the diffuse domain approximation, producing a formalism that demonstrates first-order convergence in both near- and far-field errors. We finally apply our formalism to an electrostatics problem from nanoscience with characteristic length scales ranging from 0.0001 to 300 nm.     

Discrete radiation boundary conditions for a semi-infinite layer of fluid

The paper deals with the discrete formulation of radiation boundary conditions for a layer of fluid. The problem is examined with the help of the finite difference method. The proposed radiation boundary enables us to replace an infinite layer by a finite domain. The conditions ensure near equivalence between the infinite layer and the proposed finite model. The method is consistent itself and operates on a finite number of points. The results of numerical solutions are in good agreement with the results of analytical solutions of the problem.     

Application of Bernoulli matrix method for solving two-dimensional hyperbolic telegraph equations with Dirichlet boundary conditions

The present article is devoted to develop a new approach and methodology to find the approximate solution of second order two-dimensional telegraph equations with the Dirichlet boundary conditions. We first transform the telegraph equations into equivalent partial integro-differential equations (PIDEs) which contain both initial and boundary conditions and therefore can be solved numerically in a more appropriate manner. Operational matrices of integration and differentiation of Bernoulli polynomials together with the completeness of these polynomials are used to reduce the PIDEs into the associated algebraic generalized Sylvester equations which can be solved by an efficient Krylov subspace iterative (i.e., BICGSTAB) method. The efficiency of the proposed method has been confirmed with several test examples and it is clear that the results are acceptable and found to be in good agreement with exact solutions. We have compared the numerical results of the proposed method with radial basis function method and differential quadrature method. Also, the method is simple, efficient and produces very accurate numerical results in considerably small number of basis functions and hence reduces computational effort. Moreover, the technique is easy to apply for multidimensional problems.     

Robust stability of delayed reaction-diffusion recurrent neural networks with Dirichlet boundary conditions on time scales

In this paper, the global robust exponential stability of equilibrium solution to delayed reaction-diffusion recurrent neural networks with Dirichlet boundary conditions on time scales is studied. Using topological degree theory, M-matrix method, Lyapunov functional and inequality skills, we establish some sufficient conditions for the existence, uniqueness and global robust exponential stability of equilibrium solution to delayed reaction-diffusion recurrent neural networks with Dirichlet boundary conditions on time scales. One example is given to illustrate the effectiveness of our results.     

A multigrid fluid pressure solver handling separating solid boundary conditions

We present a multigrid method for solving the linear complementarity problem (LCP) resulting from discretizing the Poisson equation subject to separating solid boundary conditions in an Eulerian liquid simulation’s pressure projection step. The method requires only a few small changes to a multigrid solver for linear systems. Our generalized solver is fast enough to handle 3D liquid simulations with separating boundary conditions in practical domain sizes. Previous methods could only handle relatively small 2D domains in reasonable time, because they used expensive quadratic programming (QP) solvers. We demonstrate our technique in several practical scenarios, including nonaxis-aligned containers and moving solids in which the omission of separating boundary conditions results in disturbing artifacts of liquid sticking to solids. Our measurements show, that the convergence rate of our LCP solver is close to that of a standard multigrid solver.     

On boundary conditions for the numerical solution of fluid dynamic problems

The stability and accuracy of various boundary treatments are analyzed for a finite difference scheme proposed by the author for the numerical solution of problems in fluid dynamics. The theory of Gustafsson, Kreiss and Sundstrom is used to establish stability and the theory of Skollermo is used to compare the accuracy of the various methods. The accuracy preductions are compared with computed results.     

A posteriori error analysis of an augmented mixed formulation in linear elasticity with mixed and Dirichlet boundary conditions

We develop a residual-based a posteriori error analysis for the augmented mixed methods introduced in [13], [14] for the problem of linear elasticity in the plane. We prove that the proposed a posteriori error estimators are both reliable and efficient. Numerical experiments confirm these theoretical properties and illustrate the ability of the corresponding adaptive algorithms to localize the singularities and large stress regions of the solutions.     

High-accuracy quadrature methods for solving boundary integral equations of steady-state anisotropic heat conduction problems with Dirichlet conditions

This paper presents the mechanical quadrature methods (MQMs) for solving the boundary integral equations of steady-state anisotropic heat conduction equation on the smooth domains and polygons, respectively. The costless and high-accurate Sidi–Israeli quadrature formula are applied to deal with the integrals in which the kernels have a logarithmic singularity. Especially, the Sidi transformation is used for the polygon cases in order to obtain a rapid convergence by degrading the singularity at the corners on the boundary. The convergence and stability of the MQMs solution are proved based on Anselone's collective compact theory. In addition, asymptotic error expansion of the MQMs shows that the approximation order is of O(h³), where h is the partition size of the boundary. Finally, numerical examples are tested and results verify the theoretical analysis.     

Wall functions in computational fluid mechanics

Wall functions are a powerful tool in complex flows calculations. Unfortunately, often blind and inappropriate use of them have given the impression that they should be avoided. Our aim though this paper is, without being exhaustive, to correct some of these prejudges and to point some advantages of wall functions, not only from the complexity point of view but also from their greater modeling capacity.     

Lagrangian finite element modelling of dam–fluid interaction: Accurate absorbing boundary conditions

The dynamic dam–fluid interaction is considered via a Lagrangian approach, based on a fluid finite element (FE) model under the assumption of small displacement and inviscid fluid. The fluid domain is discretized by enhanced displacement-based finite elements, which can be considered an evolution of those derived from the pioneering works of Bathe and Hahn [Bathe KJ, Hahn WF. On transient analysis of fluid–structure system. Comp Struct 1979;10:383–93] and of Wilson and Khalvati [Wilson EL, Khalvati M. Finite element for the dynamic analysis of fluid–solid system. Int J Numer Methods Eng 1983;19:1657–68]. The irrotational condition for inviscid fluids is imposed by the penalty method and consequentially leads to a type of micropolar media. The model is implemented using a FE code, and the numerical results of a rectangular bidimensional basin (subjected to horizontal sinusoidal acceleration) are compared with the analytical solution. It is demonstrated that the Lagrangian model is able to perform pressure and gravity wave propagation analysis, even if the gravity (or surface) waves are dispersive. The dispersion nature of surface waves indicates that the wave propagation velocity is dependent on the wave frequency.For the practical analysis of the coupled dam–fluid problem the analysed region of the basin must be reduced and the use of suitable asymptotic boundary conditions must be investigated. The classical Sommerfeld condition is implemented by means of a boundary layer of dampers and the analysis results are shown for the cases of sinusoidal forcing.The classical Sommerfeld condition is highly efficient for pressure-based FE modelling, but may not be considered fully adequate for the displacement-based FE approach. In the present paper a high-order boundary condition proposed by Higdom [Higdom RL. Radiation boundary condition for dispersive waves. SIAM J Numer Anal 1994;31:64–100] is considered. Its implementation requires the resolution of a multifreedom constraint problem, defined in terms of incremental displacements, in the ambit of dynamic time integration problems. The first- and second-order Higdon conditions are developed and implemented. The results are compared with the Sommerfeld condition results, and with the analytical unbounded problem results.Finally, a number of finite element results are presented and their related features are discussed and critically compared.