Convergence of Finite Volume Schemes for Hamilton-Jacobi Equations with Dirichlet Boundary Conditions

We study numerical methods for time-dependent Hamilton-Jacobi equations with weak Dirichlet boundary conditions. We first propose a new class of abstract monotone approximation schemes and get a convergence rate of 1/2 . Then, according to the abstract convergence results, by newly constructing monotone finite volume approximations on interior and boundary points, we obtain convergent finite volume schemes for time-dependent Hamilton-Jacobi equations with weak Dirichlet boundary conditions. Finally give some numerical results.     

A new method for the connection of non-conforming NURBS patches in isogeometric analysis

The main objective of isogeometric analysis is the use of one common data set for design and analysis. Geometrical models usually consist of multiple NURBS surface patches with non-matching parametrizations along common edges. The ability to handle non-conforming meshes is essential for isogeometric finite element systems to avoid manipulations of the geometry model. In general two possible strategies exist. The first one is to enhance the weak form of the potential with additional terms to constrain the deformations of adjacent edges to be equal, e.g. with the Lagrange Multiplier method. The second one is to establish a relation between the deformations of adjacent patches. Thus, superfluous degrees of freedom are condensated out of the system and the patches are connected naturally by shared degrees of freedom. No alteration of the weak form is necessary which simplifies the implementation in existing finite element systems. For adjacent edges with hierarchic knot vectors an analytical solution exists. In the general case of non-hierarchic knot vectors it is not possible to find an analytical solution. This contribution presents a new numerical method to establish this relation. A collocation along the connection line links the deformations of adjacent patches. A theoretical derivation of the method is given and the choice of the collocation points is investigated. Numerical studies compare the new method with the Lagrange Multiplier method and show its applicability. (© 2013 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

A boundary integral equation procedure for the mixed boundary value problem of the vector helmholtz equation

A boundary integral method is developed for the mixed boundary value problem for the vector Helmholtz equation in R³. The obtained boundary integral equations for the unknown Cauchy data build a strong elliptic system of pseudodifferential equations which can therefore be used for numerical computations using Galerkin's procedure. We show existence, uniqueness and regularity of the solution of the integral equations. Especially we give the local "edge" behavior of the solution near the submanifold which divides the Dirichlet boundary from the Neumann boundary     

Nonuniformly nonlinear elliptic equations of p-biharmonic type

This paper deals with the existence and multiplicity of weak solutions to nonlinear differential equations involving a general p-biharmonic operator (in particular, p-biharmonic operator) under Dirichlet boundary conditions or Navier boundary conditions. Our method is mainly based on variational arguments.     

A posteriori error estimation and adaptive mesh-refinement techniques

We analyse three different a posteriori error estimators for elliptic partial differential equations. They are based on the evaluation of local residuals with respect to the strong form of the differential equation, on the solution of local problems with Neumann boundary conditions, and on the solution of local problems with Dirichlet boundary conditions. We prove that all three are equivalent and yield global upper and local lower bounds for the true error. Thus adaptive mesh-refinement techniques based on these estimators are capable to detect local singularities of the solution and to appropriately refine the grid near these singularities. Some numerical examples prove the efficiency of the error estimators and the mesh-refinement techniques.     

Free vibration analysis of cracked composite beams subjected to coupled bending–torsion loads based on a first order shear deformation theory

In this paper, free vibration analysis of cracked composite beam subjected to coupled bending–torsion loading is presented. The composite beam is assumed to have an open edge crack of length a. A first order shear deformation theory is applied to count for the effect of shear deformations on natural frequencies as well as the effect of coupling in torsion and bending modes of vibration. Governing equations and boundary conditions are derived using Hamilton principle. Local flexibility matrix is used to obtain the additional boundary conditions of the beam in cracked area. After obtaining the governing equations and boundary conditions, generalized differential quadrature (GDQ) method is applied to solve the obtained eigenvalue problem. Finally, some numerical results of beams with various boundary conditions and different fiber orientations are given to show the efficiency of the method. In addition, to study the effect of shear deformations, numerical results of the current model are compared with previously given results in which shear deformations were neglected.     

Generalized impedance boundary conditions for the maxwell equations as singular perturbations problems

In this paper the Maxwell equations in an exterior domain with generalaized impedance boundary conditions of Engquist-Nédélec are considered. The particular form of the assumed boundary conditions can be considered to be a singular perturbation of the Dirichlet boundary conditions. The convergence of the solution of the Maxwell equations with these generalized impedance boundary conditions to that of the corresponding Dirichlet problem is proven. The proof uses a new integral equations method combined with results dealing with singular perturbation problems of a class of pseudo-differential operators.     

Multiple solutions for a class of nonlinear elliptic equations on the Sierpinski gasket

This paper investigates a class of nonlinear elliptic equations on a fractal domain. We establish a strong Sobolev-type inequality which leads to the existence of multiple non-trivial solutions of △u+ c(x)u = f(x, u), with zero Dirichlet boundary conditions on the Sierpihski gasket. Our existence results do not require any growth conditions of f(x,t) in t, in contrast to the classical theory of elliptic equations on smooth domains.     

An initial boundary value problem for a pseudoparabolic equation with a nonlinear boundary condition

An initial boundary value problem for a quasilinear equation of pseudoparabolic type with a nonlinear boundary condition of the Neumann–Dirichlet type is investigated in this work. From a physical point of view, the initial boundary value problem considered here is a mathematical model of quasistationary processes in semiconductors and magnets, which takes into account a wide variety of physical factors. Many approximate methods are suitable for finding eigenvalues and eigenfunctions in problems where the boundary conditions are linear with respect to the desired function and its derivatives. Among these methods, the Galerkin method leads to the simplest calculations. On the basis of a priori estimates, we prove a local existence theorem and uniqueness for a weak generalized solution of the initial boundary value problem for the quasilinear pseudoparabolic equation. A special place in the theory of nonlinear equations is occupied by the study of unbounded solutions, or, as they are called in another way, blow-up regimes. Nonlinear evolutionary problems admitting unbounded solutions are globally unsolvable. In the article, sufficient conditions for the blow-up of a solution in a finite time in a limited area with a nonlinear Neumann–Dirichlet boundary condition are obtained.     

Wavelet-Galerkin method for solving parabolic equations in finite domains

A novel wavelet-Galerkin method tailored to solve parabolic equations in finite domains is presented. The emphasis of the paper is on the development of the discretization formulations that are specific to finite domain parabolic equations with arbitrary boundary conditions based on weak form functionals. The proposed method also deals with the development of algorithms for computing the associated connection coefficients at arbitrary points. Here the Lagrange multiplier method is used to enforce the essential boundary conditions. The numerical results on a two-dimensional transient heat conducting problem are used to validate the proposed wavelet-Galerkin algorithm as an effective numerical method to solve finite domain parabolic equations.     

Strong solutions for semilinear wave equations with damping and source terms

This article deals with local existence of strong solutions for semilinear wave equations with power-like interior damping and source terms. A long-standing restriction on the range of exponents for the two nonlinearities governs the literature on wellposedness of weak solutions of finite energy. We show that this restriction may be eliminated for the existence of higher regularity solutions by employing natural methods that use the physics of the problem. This approach applies to the Cauchy problem posed on the entire ? ⁿ as well as for initial boundary problems with homogeneous Dirichlet boundary conditions.     

The Levenberg-Marquardt method and its modified versions for solving nonlinear equations with application to the inverse gravimetry problem

The Levenberg-Marquardt method and its modified versions are studied. Under some local conditions on the operator (in a neighborhood of a solution), strong and weak convergence of iterations is established with the solution error monotonically decreasing. The conditions are shown to be true for one class of nonlinear integral equations, in particular, for the structural gravimetry problem. Results of model numerical experiments for the inverse nonlinear gravimetry problem are presented.     

Stabilization of a thermoelastic Mindlin–Timoshenko plate model revisited

This paper is a continuation of our work in Grobbelaar-Van Dalsen (Appl Anal 90:1419–1449, 2011) where we showed the strong stability of models involving the thermoelastic Mindlin–Timoshenko plate equations with second sound. For the case of a plate configuration consisting of a single plate, this was accomplished in radially symmetric domains without applying any mechanical damping mechanism. Further to this result, we establish in this paper the non-exponential stability of the model for a particular configuration under mixed boundary conditions on the shear angle variables and Dirichlet boundary conditions on the displacement and thermal variables when the heat flux is described by Fourier’s law of heat conduction. We also determine the rate of polynomial decay of weak solutions of the model in a radially symmetric region under Dirichlet boundary conditions on the displacement and thermal variables and free boundary conditions on the shear angle variables.     

A STUDY ON GRADIENT BLOW UP FOR VISCOSITY SOLUTIONS OF FULLY NONLINEAR, UNIFORMLY ELLIPTIC EQUATIONS

We investigate sharp conditions for boundary and interior gradient estimates of continuous viscosity solutions to fully nonlinear,uniformly elliptic equations under Dirichlet boundary conditions. When ...     

Application and Implementation of Incorporating Local Boundary Conditions into Nonlocal Problems

We study nonlocal equations from the area of peridynamics, an instance of nonlocal wave equation, and nonlocal diffusion on bounded domains whose governing equations contain a convolution operator based on integrals. We generalize the notion of convolution to accommodate local boundary conditions. On a bounded domain, the classical operator with local boundary conditions has a purely discrete spectrum, and hence, provides a Hilbert basis. We define an abstract convolution operator using this Hilbert basis, thereby automatically satisfying local boundary conditions. The main goal in this paper is twofold: apply the concept of abstract convolution operator to nonlocal problems and carry out a numerical study of the resulting operators. We study the corresponding initial value problems with prominent boundary conditions such as periodic, antiperiodic, Neumann, and Dirichlet. To connect to the standard convolution, we give an integral representation of the abstract convolution operator. For discretization, we use a weak formulation based on a Galerkin projection and use piecewise polynomials on each element which allows discontinuities of the approximate solution at the element borders. We study convergence order of solutions with respect to polynomial order and observe optimal convergence. We depict the solutions for each boundary condition.     

New higher-order compact finite difference schemes for 1D heat conduction equations

In this paper, we present two higher-order compact finite difference schemes for solving one-dimensional (1D) heat conduction equations with Dirichlet and Neumann boundary conditions, respectively. In particular, we delicately adjust the location of the interior grid point that is next to the boundary so that the Dirichlet or Neumann boundary condition can be applied directly without discretization, and at the same time, the fifth or sixth-order compact finite difference approximations at the grid point can be obtained. On the other hand, an eighth-order compact finite difference approximation is employed for the spatial derivative at other interior grid points. Combined with the Crank–Nicholson finite difference method and Richardson extrapolation, the overall scheme can be unconditionally stable and provides much more accurate numerical solutions. Numerical errors and convergence rates of these two schemes are tested by two examples.     

On the Solution of the Exterior Dirichlet Problem for an Axisymmetric Potential

For an unbounded domain of the meridian plane with bounded smooth boundary that satisfies certain additional conditions, we develop a method for the reduction of the Dirichlet problem for an axisymmetric potential to Fredholm integral equations. In the case where the boundary of the domain is a unit circle, we obtain a solution of the exterior Dirichlet problem in explicit form.     

On the effects of the Bohm potential on a macroscopic system of self-interacting particles

We consider a nonstationary macroscopic system of self-interacting particles with an additional potential, the so called Bohm potential. We study the existence of nonnegative global solutions to the system of equations and allude to the differences to results obtained for classical models. The problem is considered on a bounded domain up to three spatial dimensions, subject to initial and Neumann boundary conditions for the particle density, and the Dirichlet boundary condition for the self-interacting potential. Moreover, the initial datum is only assumed to be nonnegative and to satisfy a weak integrability condition.     

Approximate solutions and error estimates for a Stokes boundary value problem

The aim of this paper is the numerical treatment of a boundary value problem for the system of Stokes’ equations. For this we extend the method of approximate approximations to boundary value problems. This method was introduced by Maz’ya (DFG-Kolloquium des DFG-Forschungsschwerpunktes Randelementmethoden, 1991) and has been used until now for the approximation of smooth functions defined on the whole space and for the approximation of volume potentials. In the present paper we develop an approximation procedure for the solution of the interior Dirichlet problem for the system of Stokes’ equations in two dimensions. The procedure is based on potential theoretical considerations in connection with a boundary integral equations method and consists of three approximation steps as follows. In a first step the unknown source density in the potential representation of the solution is replaced by approximate approximations. In a second step the decay behavior of the generating functions is used to gain a suitable approximation for the potential kernel, and in a third step Nyström’s method leads to a linear algebraic system for the approximate source density. For every step a convergence analysis is established and corresponding error estimates are given.     

Hydrostatic Stokes equations with non-smooth data for mixed boundary conditions

The main subject of this work is to study the concept of very weak solution for the hydrostatic Stokes system with mixed boundary conditions (non-smooth Neumann conditions on the rigid surface and homogeneous Dirichlet conditions elsewhere on the boundary). In the Stokes framework, this concept has been studied by Conca [Rev. Mat. Apl. 10 (1989)] imposing non-smooth Dirichlet boundary conditions.In this paper, we introduce the dual problem that turns out to be a hydrostatic Stokes system with non-free divergence condition. First, we obtain strong regularity for this dual problem (which can be viewed as a generalisation of the regularity results for the hydrostatic Stokes system with free divergence condition obtained by Ziane [Appl. Anal. 58 (1995)]). Afterwards, we prove existence and uniqueness of very weak solution for the (primal) problem.As a consequence of this result, the existence of strong solution for the non-stationary hydrostatic Navier-Stokes equations is proved, weakening the hypothesis over the time derivative of the wind stress tensor imposed by Guillén-González, Masmoudi and Rodríguez-Bellido [Differential Integral Equations 50 (2001)].