WorldWideScience

Sample records for boundary integral equations

  1. A Boundary Integral Equation Approach for Boundary Problem of Laplace Equation

    Institute of Scientific and Technical Information of China (English)

    SUNJian-she; YELiu-qing

    2003-01-01

    Using the second Green formula, the boundary problem of Laplace equation satisfied by potential function of static electric field is transformed to the problem of the boundary integral equation,and then a boundary integral equation approach is established by partitioning boundary using linear boundary element.

  2. Iterative solution of Hermite boundary integral equations

    Energy Technology Data Exchange (ETDEWEB)

    Gray, Leonard J [ORNL; Nintcheu Fata, Sylvain [ORNL; Ma, Ding [ORNL

    2008-01-01

    An efficient iterative method for the solution of the linear equations arising from a Hermite boundary integral approximation has been developed. Along with equations for the boundary unknowns, the Hermite system incorporates equations for the first-order surface derivatives (gradient) of the potential, and is therefore substantially larger than the matrix for a corresponding linear approximation. However, by exploiting the structure of the Hermite matrix, a two-level iterative algorithm has been shown to provide a very efficient solution algorithm. In this approach, the boundary function unknowns are treated separately from the gradient, taking advantage of the sparsity and near-positive definiteness of the gradient equations. In test problems, the new algorithm significantly reduced computation time compared to iterative solution applied to the full matrix. This approach should prove to be even more effective for the larger systems encountered in three-dimensional analysis, and increased efficiency should come from pre-conditioning of the non-sparse matrix component.

  3. A NOVEL BOUNDARY INTEGRAL EQUATION METHOD FOR LINEAR ELASTICITY--NATURAL BOUNDARY INTEGRAL EQUATION

    Institute of Scientific and Technical Information of China (English)

    Niu Zhongrong; Wang Xiuxi; Zhou Huanlin; Zhang Chenli

    2001-01-01

    The boundary integral equation (BIE) of displacement derivatives is put at a disadvantage for the difficulty involved in the evaluation of the hypersingular integrals. In this paper, the operators δij and εij are used to act on the derivative BIE. The boundary displacements, tractions and displacement derivatives are transformed into a set of new boundary tensors as boundary variables. A new BIE formulation termed natural boundary integral equation (NBIE) is obtained. The NBIE is applied to solving two-dimensional elasticity problems. In the NBIE only the strongly singular integrals are contained. The Cauchy principal value integrals occurring in the NBIE are evaluated. A combination of the NBIE and displacement BIE can be used to directly calculate the boundary stresses. The numerical results of several examples demonstrate the accuracy of the NBIE.

  4. Computing the Casimir force using regularized boundary integral equations

    Science.gov (United States)

    Kilen, Isak; Jakobsen, Per Kristen

    2014-11-01

    In this paper we use a novel regularization procedure to reduce the calculation of the Casimir force for 2D scalar fields between compact objects to the solution of a classical integral equation defined on the boundaries of the objects. The scalar fields are subject to Dirichlet boundary conditions on the object boundaries. We test the integral equation by comparing with what we get for parallel plates, concentric circles and adjacent circles using mode summation and the functional integral method. We show how symmetries in the shapes and configuration of boundaries can easily be incorporated into our method and that it leads to fast evaluation of the Casimir force for symmetric situations.

  5. APPLICATION OF BOUNDARY INTEGRAL EQUATION METHOD FOR THERMOELASTICITY PROBLEMS

    OpenAIRE

    Vorona Yu.V.; Kara I.D.

    2015-01-01

    Boundary Integral Equation Method is used for solving analytically the problems of coupled thermoelastic spherical wave propagation. The resulting mathematical expressions coincide with the solutions obtained in a conventional manner.

  6. APPLICATION OF BOUNDARY INTEGRAL EQUATION METHOD FOR THERMOELASTICITY PROBLEMS

    Directory of Open Access Journals (Sweden)

    Vorona Yu.V.

    2015-12-01

    Full Text Available Boundary Integral Equation Method is used for solving analytically the problems of coupled thermoelastic spherical wave propagation. The resulting mathematical expressions coincide with the solutions obtained in a conventional manner.

  7. On the Implementation of 3D Galerkin Boundary Integral Equations

    Energy Technology Data Exchange (ETDEWEB)

    Nintcheu Fata, Sylvain [ORNL; Gray, Leonard J [ORNL

    2010-01-01

    In this article, a reverse contribution technique is proposed to accelerate the construction of the dense influence matrices associated with a Galerkin approximation of singular and hypersingular boundary integral equations of mixed-type in potential theory. In addition, a general-purpose sparse preconditioner for boundary element methods has also been developed to successfully deal with ill-conditioned linear systems arising from the discretization of mixed boundary-value problems on non-smooth surfaces. The proposed preconditioner, which originates from the precorrected-FFT method, is sparse, easy to generate and apply in a Krylov subspace iterative solution of discretized boundary integral equations. Moreover, an approximate inverse of the preconditioner is implicitly built by employing an incomplete LU factorization. Numerical experiments involving mixed boundary-value problems for the Laplace equation are included to illustrate the performance and validity of the proposed techniques.

  8. Boundary Integral Equations and A Posteriori Error Estimates

    Institute of Scientific and Technical Information of China (English)

    YU Dehao; ZHAO Longhua

    2005-01-01

    Adaptive methods have been rapidly developed and applied in many fields of scientific and engineering computing. Reliable and efficient a posteriori error estimates play key roles for both adaptive finite element and boundary element methods. The aim of this paper is to develop a posteriori error estimates for boundary element methods. The standard a posteriori error estimates for boundary element methods are obtained from the classical boundary integral equations. This paper presents hyper-singular a posteriori error estimates based on the hyper-singular integral equations. Three kinds of residuals are used as the estimates for boundary element errors. The theoretical analysis and numerical examples show that the hyper-singular residuals are good a posteriori error indicators in many adaptive boundary element computations.

  9. Galerkin Boundary Integral Analysis for the 3D Helmholtz Equation

    Energy Technology Data Exchange (ETDEWEB)

    Swager, Melissa [Emporia State University; Gray, Leonard J [ORNL; Nintcheu Fata, Sylvain [ORNL

    2010-01-01

    A linear element Galerkin boundary integral analysis for the three-dimensional Helmholtz equation is presented. The emphasis is on solving acoustic scattering by an open (crack) surface, and to this end both a dual equation formulation and a symmetric hypersingular formulation have been developed. All singular integrals are defined and evaluated via a boundary limit process, facilitating the evaluation of the (finite) hypersingular Galerkin integral. This limit process is also the basis for the algorithm for post-processing of the surface gradient. The analytic integrations required by the limit process are carried out by employing a Taylor series expansion for the exponential factor in the Helmholtz fundamental solutions. For the open surface, the implementations are validated by comparing the numerical results obtained by using the two different methods.

  10. Advanced applications of boundary-integral equation methods

    International Nuclear Information System (INIS)

    The BIE (boundary integral equation) method is based on the numerical solution of a set of integral constraint equations which couple boundary tractions (stresses) to boundary displacements. Thus the dimensionality of the problem is reduced by one; only boundary geometry and data are discretized. Stresses at any set of selected interior points are computed following the boundary solution without any further numerical approximations. Thus, the BIE method has inherently greater resolution capability for stress gradients than does the finite element method. Conversely, the BIE method is not efficient for problems involving significant inhomogeneity such as in multi-thin-layered materials, or in elastoplasticity. Some progress in applyiing the BIE method to the latter problem has been made but much more work remains. Further, the BIE method is only optional for problems with significant stress risers, and only when boundary stresses are most important. Interior stress calculations are expensive, per point, and can drive the solution costs up rapidly. The current report summarizes some of the advanced elastic applications of fracture mechanics and three-dimensional stress analysis, while referencing some of the much broader developmental effort. Future emphasis is needed to exploit the BIE method in conjunction with other techniques such as the finite element method through the creation of hybrid stress analysis methods

  11. Advanced applications of boundary-integral equation methods

    International Nuclear Information System (INIS)

    Numerical analysis has become the basic tool for both design and research problems in solid mechanics. The boundary-integral equation (BIE) method is based on classical mathematical techniques but is finding new life as a basic stress analysis tool for engineering applications. The BIE method is based on the numerical solution of a set of integral constraint equations which couple boundary tractions (stresses) to boundary displacements. Thus the dimensionality of the problem is reduced by one; only boundary geometry and data are discretized. Stresses at any set of selected interior points are computed following the boundary solution without any further numerical approximations. Thus, the BIE method has inherently greater resolution capability for stress gradients than does the finite element method. Conversely, the BIE method is not efficient for problems involving significant inhomogeneity such as in multi-thin-layered materials, or in elastoplasticity. Some progress in applying the BIE method to the latter problem has been made but much more work remains. Further, the BIE method is only optional for problems with significant stress risers, and only when boundary stresses are more important. Interior stress calculations are expensive, per point, and can drive the solution costs up rapidly. The current report summarizes some of the advanced elastic applications of fracture mechanics and three-dimensional stress analysis, while referring some of the much broader developmental effort. (Auth.)

  12. Tokamak plasma shape identification based on the boundary integral equations

    International Nuclear Information System (INIS)

    A necessary condition for tokamak plasma shape identification is discussed and a new identification method is proposed in this article. This method is based on the boundary integral equations governing a vacuum region around a plasma with only the measurement of either magnetic fluxes or magnetic flux intensities. It can identify various plasmas with low to high ellipticities with the precision determined by the number of the magnetic sensors. This method is applicable to real-time control and visualization using a 'table-look-up' procedure. (author)

  13. A spectral boundary integral equation method for the 2D Helmholtz equation

    International Nuclear Information System (INIS)

    In this paper, we present a new numerical formulation of solving the boundary integral equations reformulated from the Helmholtz equation. The boundaries of the problems are assumed to be smooth closed contours. The solution on the boundary is treated as a periodic function, which is in turn approximated by a truncated Fourier series. A Fourier collocation method is followed in which the boundary integral equation is transformed into a system of algebraic equations. It is shown that in order to achieve spectral accuracy for the numerical formulation, the non-smoothness of the integral kernels, associated with the Helmholtz equation, must be carefully removed. The emphasis of the paper is on investigating the essential elements of removing the non-smoothness of the integral kernels in the spectral implementation. The present method is robust for a general smooth boundary contour. Aspects of efficient implementation of the method using FFT are also discussed. Numerical examples of wave scattering are given in which the exponential accuracy of the present numerical method is demonstrated. 15 refs., 3 figs., 4 tabs

  14. Advanced applications of boundary-integral equation methods

    International Nuclear Information System (INIS)

    Numerical analysis has become the basic tool for both design and research problems in solid mechanics. The need for accuracy and detail, plus the availablity of the high speed computer has led to the development of many new modeling methods ranging from general purpose structural analysis finite element programs to special purpose research programs. The boundary-integral equation (BIE) method is based on classical mathematical techniques but is finding new life as a basic stress analysis tool for engineering applications. The paper summarizes some advanced elastic applications of fracture mechanics and three-dimensional stress analysis, while referencing some of the much broader developmental effort. Future emphasis is needed to exploit the BIE method in conjunction with other techniques such as the finite element method through the creation of hybrid stress analysis methods. (Auth.)

  15. Supersymmetric extension of the reflection equation algebra and integrable boundary conditions in doped spin-1 model

    International Nuclear Information System (INIS)

    The reflection equation algebra of Sklyanin is extended to the supersymmetric case. A graded reflection equation algebra is proposed and the corresponding graded (supersymmetric) boundary quantum inverse scattering method (QISM) is formulated. As an application, integrable open-boundary conditions for the doped spin-1 chain of the supersymmetric t-J model are studied in the framework of the boundary QISM. Diagonal boundary K-matrices are found and four classes of integrable boundary terms are determined. (author)

  16. RESTRICTED NONLINEAR APPROXIMATION AND SINGULAR SOLUTIONS OF BOUNDARY INTEGRAL EQUATIONS

    Institute of Scientific and Technical Information of China (English)

    Reinhard Hochmuth

    2002-01-01

    This paper studies several problems, which are potentially relevant for the construction of adaptive numerical schemes. First, biorthogonal spline wavelets on [0,1 ] are chosen as a starting point for characterizations of functions in Besov spaces B , (0,1) with 0<σ<∞ and (1+σ)-1<τ<∞. Such function spaces are known to be related to nonlinear approximation. Then so called restricted nonlinear approximation procedures with respect to Sobolev space norms are considered. Besides characterization results Jackson type estimates for various tree-type and tresholding algorithms are investigated. Finally known approximation results for geometry induced singularity functions of boundary integeral equations are combined with the characterization results for restricted nonlinear approximation to show Besov space regularity results.

  17. Integral equations

    CERN Document Server

    Moiseiwitsch, B L

    2005-01-01

    Two distinct but related approaches hold the solutions to many mathematical problems--the forms of expression known as differential and integral equations. The method employed by the integral equation approach specifically includes the boundary conditions, which confers a valuable advantage. In addition, the integral equation approach leads naturally to the solution of the problem--under suitable conditions--in the form of an infinite series.Geared toward upper-level undergraduate students, this text focuses chiefly upon linear integral equations. It begins with a straightforward account, acco

  18. New Existence Results for Nonlinear Fractional Differential Equations with Three-Point Integral Boundary Conditions

    Directory of Open Access Journals (Sweden)

    Ntouyas SotirisK

    2011-01-01

    Full Text Available This paper studies a boundary value problem of nonlinear fractional differential equations of order with three-point integral boundary conditions. Some new existence and uniqueness results are obtained by using standard fixed point theorems and Leray-Schauder degree theory. Our results are new in the sense that the nonlocal parameter in three-point integral boundary conditions appears in the integral part of the conditions in contrast to the available literature on three-point boundary value problems which deals with the three-point boundary conditions restrictions on the solution or gradient of the solution of the problem. Some illustrative examples are also discussed.

  19. Weak solutions for nonlinear fractional differential equations with integral boundary conditions in Banach spaces

    OpenAIRE

    Mouffak Benchohra; Fatima-Zohra Mostefai

    2012-01-01

    The aim of this paper is to investigate a class of boundary value problems for fractional differential equations involving nonlinear integral conditions. The main tool used in our considerations is the technique associated with measures of weak noncompactness.

  20. A well-conditioned boundary integral equation for transmission problems of electromagnetism

    OpenAIRE

    Levadoux, David; Millot, Florence; Pernet, Sébastien

    2015-01-01

    We propose a new well-conditioned boundary integral equation to solve transmission problems of electromagnetism. This equation is well posed and appears as a compact perturbation of the identity leading to fast iterative solutions without the help of any preconditioner. Some numerical experiments confirm this result.

  1. Boundary integral equation Neumann-to-Dirichlet map method for gratings in conical diffraction.

    Science.gov (United States)

    Wu, Yumao; Lu, Ya Yan

    2011-06-01

    Boundary integral equation methods for diffraction gratings are particularly suitable for gratings with complicated material interfaces but are difficult to implement due to the quasi-periodic Green's function and the singular integrals at the corners. In this paper, the boundary integral equation Neumann-to-Dirichlet map method for in-plane diffraction problems of gratings [Y. Wu and Y. Y. Lu, J. Opt. Soc. Am. A26, 2444 (2009)] is extended to conical diffraction problems. The method uses boundary integral equations to calculate the so-called Neumann-to-Dirichlet maps for homogeneous subdomains of the grating, so that the quasi-periodic Green's functions can be avoided. Since wave field components are coupled on material interfaces with the involvement of tangential derivatives, a least squares polynomial approximation technique is developed to evaluate tangential derivatives along these interfaces for conical diffraction problems. Numerical examples indicate that the method performs equally well for dielectric or metallic gratings. PMID:21643404

  2. An integral equation formulation for the boundary-finite element model in eccentrically stiffened plate bending

    Energy Technology Data Exchange (ETDEWEB)

    Shih, H.R. [Jackson State Univ., MS (United States); Duffield, R.C.; Lin, J. [Univ. of Missouri, Columbia, MO (United States)

    1996-10-01

    An integral equation formulation and a numerical procedure for a boundary-finite element technique are developed for the static analysis of a stiffened plate with eccentric stiffeners. This formulation employs the fundamental solution associated with unstiffened plate bending and plane stress problems. With this approach, the resulting integral equation not only contained integrals along the perimeter of the stiffened but additional integrals along the stiffeners and the interface between the plate and its stiffeners. Thus the domain of the plate has to be divided into zones between the stiffeners. Each zone is modeled by boundary elements and stiffeners by finite elements. In this paper, the boundary element solution procedures for plate bending and in-plane problems are presented. The zone technique which permits coupling of unstiffened plate boundary element with stiffener finite elements is presented as well. Numerical example is given to demonstrate the effectiveness of this approach.

  3. Analytical Solution of Boundary Integral Equations for 2-D Steady Linear Wave Problems

    Institute of Scientific and Technical Information of China (English)

    J.M. Chuang

    2005-01-01

    Based on the Fourier transform, the analytical solution of boundary integral equations formulated for the complex velocity of a 2-D steady linear surface flow is derived. It has been found that before the radiation condition is imposed,free waves appear both far upstream and downstream. In order to cancel the free waves in far upstream regions, the eigensolution of a specific eigenvalue, which satisfies the homogeneous boundary integral equation, is found and superposed to the analytical solution. An example, a submerged vortex, is used to demonstrate the derived analytical solution. Furthermore,an analytical approach to imposing the radiation condition in the numerical solution of boundary integral equations for 2-D steady linear wave problems is proposed.

  4. Retarded potentials and time domain boundary integral equations a road map

    CERN Document Server

    Sayas, Francisco-Javier

    2016-01-01

    This book offers a thorough and self-contained exposition of the mathematics of time-domain boundary integral equations associated to the wave equation, including applications to scattering of acoustic and elastic waves. The book offers two different approaches for the analysis of these integral equations, including a systematic treatment of their numerical discretization using Galerkin (Boundary Element) methods in the space variables and Convolution Quadrature in the time variable. The first approach follows classical work started in the late eighties, based on Laplace transforms estimates. This approach has been refined and made more accessible by tailoring the necessary mathematical tools, avoiding an excess of generality. A second approach contains a novel point of view that the author and some of his collaborators have been developing in recent years, using the semigroup theory of evolution equations to obtain improved results. The extension to electromagnetic waves is explained in one of the appendices...

  5. Fractal boundary value problems for integral and differential equations with local fractional operators

    Directory of Open Access Journals (Sweden)

    Yang Xiao-Jun

    2015-01-01

    Full Text Available In the present paper we investigate the fractal boundary value problems for the Fredholm\\Volterra integral equations, heat conduction and wave equations by using the local fractional decomposition method. The operator is described by the local fractional operators. The four illustrative examples are given to elaborate the accuracy and reliability of the obtained results. [Projekat Ministarstva nauke Republike Srbije, br. OI 174001, III41006 i br. TI 35006

  6. Weak solutions for nonlinear fractional differential equations with integral boundary conditions in Banach spaces

    Directory of Open Access Journals (Sweden)

    Mouffak Benchohra

    2012-01-01

    Full Text Available The aim of this paper is to investigate a class of boundary value problems for fractional differential equations involving nonlinear integral conditions. The main tool used in our considerations is the technique associated with measures of weak noncompactness.

  7. Existence and uniqueness of solutions for nonlinear hyperbolic fractional differential equation with integral boundary conditions

    Directory of Open Access Journals (Sweden)

    Brahim Tellab

    2016-01-01

    Full Text Available In this paper, we investigate the existence and uniqueness of solutions for second order nonlinear fractional differential equation with integral boundary conditions. Our result is an application of the Banach contraction principle and the Krasnoselskii fixed point theorem.

  8. Boundary integral equation formulations for free-surface flow problems in two and three dimensions

    OpenAIRE

    Romate, J.E.; Zandbergen, P.J.

    1989-01-01

    To compute the transient solution of free-surface flow problems in two and three dimensions boundary integral equation formulations are considered. Consistent lower and higher order approximations based on small curvature expansions are compared and applied to a time-dependent, linear free-surface wave problem.

  9. On preconditioning techniques for dense linear systems arising from singular boundary integral equations

    Energy Technology Data Exchange (ETDEWEB)

    Chen, Ke [Univ. of Liverpool (United Kingdom)

    1996-12-31

    We study various preconditioning techniques for the iterative solution of boundary integral equations, and aim to provide a theory for a class of sparse preconditioners. Two related ideas are explored here: singularity separation and inverse approximation. Our preliminary conclusion is that singularity separation based preconditioners perform better than approximate inverse based while it is desirable to have both features.

  10. A boundary integral equation method for the calculation of the effective permittivity of periodic composites

    OpenAIRE

    Sareni, Bruno; Krähenbühl, Laurent; Beroual, Abderrahmane; Nicolas, Alain; Brosseau, C.

    1997-01-01

    We present a numerical method based upon the resolution of boundary integral equations for the calculation of the effective permittivity of a lossless composite structure consisting of a two component mixture, each with its own dielectric anti shape characteristics. The topological arrangements considered are periodic lattices inhomogeneities. Our numerical simulations are compared to the effective medium approach and with results of previous works.

  11. Open spin chains with generic integrable boundaries: Baxter equation and Bethe ansatz completeness from SOV

    CERN Document Server

    Kitanine, N; Niccoli, G

    2014-01-01

    We solve the longstanding problem to define a functional characterization of the spectrum of the transfer matrix associated to the most general spin-1/2 representations of the 6-vertex reflection algebra for general inhomogeneous chains. The corresponding homogeneous limit reproduces the spectrum of the Hamiltonian of the spin-1/2 open XXZ and XXX quantum chains with the most general integrable boundaries. The spectrum is characterized by a second order finite difference functional equation of Baxter type with an inhomogeneous term which vanishes only for some special but yet interesting non-diagonal boundary conditions. This functional equation is shown to be equivalent to the known separation of variable (SOV) representation hence proving that it defines a complete characterization of the transfer matrix spectrum. The polynomial character of the Q-function allows us then to show that a finite system of equations of generalized Bethe type can be similarly used to describe the complete transfer matrix spectru...

  12. Boundary integral equation methods in eigenvalue problems of elastodynamics and thin plates

    CERN Document Server

    Kitahara, M

    1985-01-01

    The boundary integral equation (BIE) method has been used more and more in the last 20 years for solving various engineering problems. It has important advantages over other techniques for numerical treatment of a wide class of boundary value problems and is now regarded as an indispensable tool for potential problems, electromagnetism problems, heat transfer, fluid flow, elastostatics, stress concentration and fracture problems, geomechanical problems, and steady-state and transient electrodynamics.In this book, the author gives a complete, thorough and detailed survey of the method. It pro

  13. Inverse Coefficient Problem of the Parabolic Equation with Periodic Boundary and Integral Overdetermination Conditions

    OpenAIRE

    Fatma Kanca

    2013-01-01

    This paper investigates the inverse problem of finding a time-dependent diffusion coefficient in a parabolic equation with the periodic boundary and integral overdetermination conditions. Under some assumption on the data, the existence, uniqueness, and continuous dependence on the data of the solution are shown by using the generalized Fourier method. The accuracy and computational efficiency of the proposed method are verified with the help of the numerical examples.

  14. Solution of the Stokes system by boundary integral equations and fixed point iterative schemes

    International Nuclear Information System (INIS)

    The solution to the exterior three dimensional Stokes problem is sought in the form of a single layer potential of unknown density. This reduces the problem to a boundary integral equation of the first kind whose operator is the velocity component of the single layer potential. It is shown that this component is an isomorphism between two appropriate Sobolev spaces containing the unknown densities and the data respectively. The isomorphism corresponds to a variational problem with coercive bilinear form. The latter property allows us to consider various fixed point iterative schemes that converge to the unique solution of the integral equation. Explicit error estimates are also obtained. The successive approximations are also considered in a more computable form by using the product integration method of Atkinson. (author). 47 refs

  15. A general integral form of the boundary-layer equation for incompressible flow with an application to the calculation of the separation point of turbulent boundary layers

    Science.gov (United States)

    Tetervin, Neal; Lin, Chia Chiao

    1951-01-01

    A general integral form of the boundary-layer equation, valid for either laminar or turbulent incompressible boundary-layer flow, is derived. By using the experimental finding that all velocity profiles of the turbulent boundary layer form essentially a single-parameter family, the general equation is changed to an equation for the space rate of change of the velocity-profile shape parameter. The lack of precise knowledge concerning the surface shear and the distribution of the shearing stress across turbulent boundary layers prevented the attainment of a reliable method for calculating the behavior of turbulent boundary layers.

  16. On the Assessment of Acoustic Scattering and Shielding by Time Domain Boundary Integral Equation Solutions

    Science.gov (United States)

    Hu, Fang Q.; Pizzo, Michelle E.; Nark, Douglas M.

    2016-01-01

    Based on the time domain boundary integral equation formulation of the linear convective wave equation, a computational tool dubbed Time Domain Fast Acoustic Scattering Toolkit (TD-FAST) has recently been under development. The time domain approach has a distinct advantage that the solutions at all frequencies are obtained in a single computation. In this paper, the formulation of the integral equation, as well as its stabilization by the Burton-Miller type reformulation, is extended to cases of a constant mean flow in an arbitrary direction. In addition, a "Source Surface" is also introduced in the formulation that can be employed to encapsulate regions of noise sources and to facilitate coupling with CFD simulations. This is particularly useful for applications where the noise sources are not easily described by analytical source terms. Numerical examples are presented to assess the accuracy of the formulation, including a computation of noise shielding by a thin barrier motivated by recent Historical Baseline F31A31 open rotor noise shielding experiments. Furthermore, spatial resolution requirements of the time domain boundary element method are also assessed using point per wavelength metrics. It is found that, using only constant basis functions and high-order quadrature for surface integration, relative errors of less than 2% may be obtained when the surface spatial resolution is 5 points-per-wavelength (PPW) or 25 points-per-wavelength squared (PPW2).

  17. Boundary-integral equation in engineering stress and fracture mechanics. Technical report

    International Nuclear Information System (INIS)

    A boundary-integral equation (BIE) method, as employed by a computer program called PESTIE (Plane Elastic Solution Technique of Integral Equations), is described for application to engineering stress and fracture mechanics problems of interest to the electric power generating industry. The paper describes a unique combination of several numerical analysis capabilities and user-oriented features that produce significant advantages over finite element programs and earlier BIE programs. These advantages are demonstrated by comparing numerical results and performance of PESTIE with those of other programs for a series of stress concentration and stress intensity factor problems. The examples include the effect of closely spaced notches (e.g. as encountered in the rim area of a turbine rotor spindle), and the stress intensity factor of a shallow wide crack (e.g. a longitudinal crack in a large diameter pipe)

  18. Application of the boundary-integral-equation method to the three-dimensional thermoelastic problem

    International Nuclear Information System (INIS)

    It is possible to solve three-dimensional steady state thermoelastic problems from data of temperature, fluxes, tractions, and displacements on the boundary alone. The numerical discretisation of these integral equations is based on a technique similar to that of finite elements. The boundary alone needs to be discretised. Two-dimensional examples are given to prove the accuracy and the feasibility. In the case of three-dimensional problems, the surface is represented by eight nodes quadrilateral elements and six nodes triangular elements. The unknown (temperature, flux, displacement, traction) may be considered to vary linearly, quadratically, or cubically, with respect to he intrinsic coordinates of each element. The integration is performed numerically using Gaussian quadrature formulas, for which the number of integration points is chosen automatically by the program, so that the upper bound of error in integration is minimized. In order to obtain a banded form matrix and also to be able to study elongated structures, the body is divided into subregions, for each of which the integral equations are written. The thermoelastic stress and deformation field are obtained from two successive calculations on the same mesh. The first gives the thermal field. The second, taking account of the thermal field, calculates the thermoelastic displacement and stress field. This type of approach is especially suited to complicated three-dimensional thick structures for which finite element procedures are very expensive. The data are simple to generate, for data on the boundary alone have to be given. Much time can be saved by the user of the program in data generations and data checks

  19. A Fast Spectral Galerkin Method for Hypersingular Boundary Integral Equations in Potential Theory

    Energy Technology Data Exchange (ETDEWEB)

    Nintcheu Fata, Sylvain [ORNL; Gray, Leonard J [ORNL

    2009-01-01

    This research is focused on the development of a fast spectral method to accelerate the solution of three-dimensional hypersingular boundary integral equations of potential theory. Based on a Galerkin approximation, the Fast Fourier Transform and local interpolation operators, the proposed method is a generalization of the Precorrected-FFT technique to deal with double-layer potential kernels, hypersingular kernels and higher-order basis functions. Numerical examples utilizing piecewise linear shape functions are included to illustrate the performance of the method.

  20. TIME–HARMONIC BEHAVIOUR OF CRACKED PIEZOELECTRIC SOLID BY BOUNDARY INTEGRAL EQUATION METHOD

    Directory of Open Access Journals (Sweden)

    Rangelov Tsviatko

    2014-03-01

    Full Text Available Anti-plane cracked functionally graded finite piezoelectric solid under time-harmonic elecromechanical load is studied by a non-hypersingular traction boundary integral equation method (BIEM. Exponentially varying material properties are considered. Numerical solutions are obtained by using Mathematica. The dependance of the intensity factors (IF - mechanical stress intensity factor (SIF and electrical field intensity factor (FIF on the inhomogeneous material parameters, on the type and frequency of the dynamic load and on the crack position are analyzed by numerical illustrative examples

  1. An exterior Poisson solver using fast direct methods and boundary integral equations with applications to nonlinear potential flow

    Science.gov (United States)

    Young, D. P.; Woo, A. C.; Bussoletti, J. E.; Johnson, F. T.

    1986-01-01

    A general method is developed combining fast direct methods and boundary integral equation methods to solve Poisson's equation on irregular exterior regions. The method requires O(N log N) operations where N is the number of grid points. Error estimates are given that hold for regions with corners and other boundary irregularities. Computational results are given in the context of computational aerodynamics for a two-dimensional lifting airfoil. Solutions of boundary integral equations for lifting and nonlifting aerodynamic configurations using preconditioned conjugate gradient are examined for varying degrees of thinness.

  2. Boundary integral equation method for added mass in arrays of cylinders

    International Nuclear Information System (INIS)

    The dynamic behavior of a group of cylinders in fluid, such as heat exchanger tubes and nuclear fuel assemblies, is strongly influenced by the surrounding fluid. Although, the added mass of such clusters of cylinders has been studied by many researchers with various analytical methods and numerical methods, no attempt has been made so far to analyze these problems by the Boundary Integral Equation Method (BIEM). This paper presents a BIEM model simulating the added mass arising when clusters of cylinders vibrate in an inviscid and incompressible fluid. In this model, perturbed fluid pressure is described by a two- dimensional Laplace equation. The primary advantage of this approach compared with other numerical methods, e.g., the finite element method (FEM), is that the integration and discretization of the model are only needed on boundary rather than in whole domain. Therefore, the proposed approach is much more economical than the finite element method. Various numerical examples are subsequently presented in this paper to illustrate the methodology and to demonstrate its accuracy

  3. MULTIPLICITY OF NONNEGATIVE SOLUTIONS TO SECOND-ORDER SINGULAR DIFFERENTIAL EQUATION WITH INTEGRAL BOUNDARY CONDITIONS IN BANACH SPACES

    Institute of Scientific and Technical Information of China (English)

    2010-01-01

    This paper investigates the existence and multiplicity of nonnegative solutions to a singular nonlinear boundary value problem of second order differential equations with integral boundary conditions in a Banach space. The arguments are based on the construction of a nonempty bounded open convex set and fixed point index theory. Our nonlinearity possesses singularity and first derivative which makes it different with that in [10].

  4. Boundary integral equation methods and numerical solutions thin plates on an elastic foundation

    CERN Document Server

    Constanda, Christian; Hamill, William

    2016-01-01

    This book presents and explains a general, efficient, and elegant method for solving the Dirichlet, Neumann, and Robin boundary value problems for the extensional deformation of a thin plate on an elastic foundation. The solutions of these problems are obtained both analytically—by means of direct and indirect boundary integral equation methods (BIEMs)—and numerically, through the application of a boundary element technique. The text discusses the methodology for constructing a BIEM, deriving all the attending mathematical properties with full rigor. The model investigated in the book can serve as a template for the study of any linear elliptic two-dimensional problem with constant coefficients. The representation of the solution in terms of single-layer and double-layer potentials is pivotal in the development of a BIEM, which, in turn, forms the basis for the second part of the book, where approximate solutions are computed with a high degree of accuracy. The book is intended for graduate students and r...

  5. Initial-boundary-value problems for discrete evolution equations: discrete linear Schrödinger and integrable discrete nonlinear Schrödinger equations

    International Nuclear Information System (INIS)

    We present a method to solve initial-boundary-value problems for linear and integrable nonlinear differential-difference evolution equations. The method is the discrete version of the one developed by A S Fokas to solve initial-boundary-value problems for linear and integrable nonlinear partial differential equations via an extension of the inverse scattering transform. The method takes advantage of the Lax pair formulation for both linear and nonlinear equations, and is based on the simultaneous spectral analysis of both parts of the Lax pair. A key role is also played by the global algebraic relation that couples all known and unknown boundary values. Even though additional technical complications arise in discrete problems compared to continuum ones, we show that a similar approach can also solve initial-boundary-value problems for linear and integrable nonlinear differential-difference equations. We demonstrate the method by solving initial-boundary-value problems for the discrete analogue of both the linear and the nonlinear Schrödinger equations, comparing the solution to those of the corresponding continuum problems. In the linear case we also explicitly discuss Robin-type boundary conditions not solvable by Fourier series. In the nonlinear case, we also identify the linearizable boundary conditions, we discuss the elimination of the unknown boundary datum, we obtain explicitly the linear and continuum limit of the solution, and we write the soliton solutions

  6. Boundary integral equation method calculations of surface regression effects in flame spreading

    Science.gov (United States)

    Altenkirch, R. A.; Rezayat, M.; Eichhorn, R.; Rizzo, F. J.

    1982-01-01

    A solid-phase conduction problem that is a modified version of one that has been treated previously in the literature and is applicable to flame spreading over a pyrolyzing fuel is solved using a boundary integral equation (BIE) method. Results are compared to surface temperature measurements that can be found in the literature. In addition, the heat conducted through the solid forward of the flame, the heat transfer responsible for sustaining the flame, is also computed in terms of the Peclet number based on a heated layer depth using the BIE method and approximate methods based on asymptotic expansions. Agreement between computed and experimental results is quite good as is agreement between the BIE and the approximate results.

  7. INTERFACIAL CRACK ANALYSIS IN THREE-DIMENSIONAL TRANSVERSELY ISOTROPIC BI-MATERIALS BY BOUNDARY INTEGRAL EQUATION METHOD

    Institute of Scientific and Technical Information of China (English)

    ZHAO Ming-hao; LI Dong-xia; SHEN Ya-peng

    2005-01-01

    The integral-differential equations for three-dimensional planar interfacial cracks of arbitrary shape in transversely isotropic bimaterials were derived by virtue of the Somigliana identity and the fundamental solutions, in which the displacement discontinuities across the crack faces are the unknowns to be determined. The interface is parallel to both the planes of isotropy. The singular behaviors of displacement and stress near the crack border were analyzed and the stress singularity indexes were obtained by integral equation method. The stress intensity factors were expressed in terms of the displacement discontinuities. In the non-oscillatory case, the hyper-singular boundary integral-differential equations were reduced to hyper-singular boundary integral equations similar to those of homogeneously isotropie materials.

  8. Existence and Analytic Approximation of Solutions of Duffing Type Nonlinear Integro-Differential Equation with Integral Boundary Conditions

    Directory of Open Access Journals (Sweden)

    Alsaedi Ahmed

    2009-01-01

    Full Text Available A generalized quasilinearization technique is developed to obtain a sequence of approximate solutions converging monotonically and quadratically to a unique solution of a boundary value problem involving Duffing type nonlinear integro-differential equation with integral boundary conditions. The convergence of order for the sequence of iterates is also established. It is found that the work presented in this paper not only produces new results but also yields several old results in certain limits.

  9. Boundary transfer matrices and boundary quantum KZ equations

    Energy Technology Data Exchange (ETDEWEB)

    Vlaar, Bart, E-mail: Bart.Vlaar@nottingham.ac.uk [School of Mathematical Sciences, University of Nottingham, Nottingham NG7 2RD (United Kingdom)

    2015-07-15

    A simple relation between inhomogeneous transfer matrices and boundary quantum Knizhnik-Zamolodchikov (KZ) equations is exhibited for quantum integrable systems with reflecting boundary conditions, analogous to an observation by Gaudin for periodic systems. Thus, the boundary quantum KZ equations receive a new motivation. We also derive the commutativity of Sklyanin’s boundary transfer matrices by merely imposing appropriate reflection equations, in particular without using the conditions of crossing symmetry and unitarity of the R-matrix.

  10. Periodic Boundary Value Problems for First-Order Impulsive Functional Integrodifferential Equations with Integral-Jump Conditions

    Directory of Open Access Journals (Sweden)

    Chatthai Thaiprayoon

    2014-01-01

    Full Text Available By developing a new comparison result and using the monotone iterative technique, we are able to obtain existence of minimal and maximal solutions of periodic boundary value problems for first-order impulsive functional integrodifferential equations with integral-jump conditions. An example is also given to illustrate our results.

  11. Sequential Derivatives of Nonlinear q-Difference Equations with Three-Point q-Integral Boundary Conditions

    Directory of Open Access Journals (Sweden)

    Nittaya Pongarm

    2013-01-01

    Full Text Available This paper studies sufficient conditions for the existence of solutions to the problem of sequential derivatives of nonlinear q-difference equations with three-point q-integral boundary conditions. Our results are concerned with several quantum numbers of derivatives and integrals. By using Banach's contraction mapping, Krasnoselskii's fixed-point theorem, and Leray-Schauder degree theory, some new existence results are obtained. Two examples illustrate our results.

  12. A Collocation Method for Volterra Integral Equations with Diagonal and Boundary Singularities

    Science.gov (United States)

    Kolk, Marek; Pedas, Arvet; Vainikko, Gennadi

    2009-08-01

    We propose a smoothing technique associated with piecewise polynomial collocation methods for solving linear weakly singular Volterra integral equations of the second kind with kernels which, in addition to a diagonal singularity, may have a singularity at the initial point of the interval of integration.

  13. The open boundary equation

    Directory of Open Access Journals (Sweden)

    D. Diederen

    2015-06-01

    Full Text Available We present a new equation describing the hydrodynamics in infinitely long tidal channels (i.e., no reflection under the influence of oceanic forcing. The proposed equation is a simple relationship between partial derivatives of water level and velocity. It is formally derived for a progressive wave in a frictionless, prismatic, tidal channel with a horizontal bed. Assessment of a large number of numerical simulations, where an open boundary condition is posed at a certain distance landward, suggests that it can also be considered accurate in the more natural case of converging estuaries with nonlinear friction and a bed slope. The equation follows from the open boundary condition and is therefore a part of the problem formulation for an infinite tidal channel. This finding provides a practical tool for evaluating tidal wave dynamics, by reconstructing the temporal variation of the velocity based on local observations of the water level, providing a fully local open boundary condition and allowing for local friction calibration.

  14. Integral equations

    CERN Document Server

    Tricomi, Francesco Giacomo

    1957-01-01

    This classic text on integral equations by the late Professor F. G. Tricomi, of the Mathematics Faculty of the University of Turin, Italy, presents an authoritative, well-written treatment of the subject at the graduate or advanced undergraduate level. To render the book accessible to as wide an audience as possible, the author has kept the mathematical knowledge required on the part of the reader to a minimum; a solid foundation in differential and integral calculus, together with some knowledge of the theory of functions is sufficient. The book is divided into four chapters, with two useful

  15. An application of the K-integrals for solving the radiative transfer equation with Fresnel boundary conditions

    International Nuclear Information System (INIS)

    We introduce the reader to an approximate method of solving the transport equation which was developed in the context of neutron thermalisation by Kladnik and Kuscer in 1962 [Kladnik R, Kuscer I. Velocity dependent Milne's problem. Nucl Sci Eng 1962;13:149]. Essentially the method is based upon two special weighted integrals of the one-dimensional transport equation which are valid regardless of the boundary conditions, and any solution must satisfy these integral relationships which are called the K-integrals. To obtain an approximate solution to the transport equation we turn the argument around and insist that any approximate solution must also satisfy the K-integrals. These integrals are particularly useful when the problem under consideration cannot be solved easily by analytic methods. It also has the marked advantage of being applicable to problems where there is energy exchange in a collision and anisotropy of scattering. To establish the feasibility of the method we obtain a number of approximate solutions using the K-integral method for problems to which we have exact analytical solutions. This enables us to validate the method. It is then applied to a new problem that has not yet been solved; namely the calculation of the discontinuity in the scalar intensity at the boundary between two optically dissimilar materials

  16. Integral equations and boundary-element solution for static potential in a general piece-wise homogeneous volume conductor

    CERN Document Server

    Stenroos, Matti

    2016-01-01

    Boundary element methods (BEM) are used for forward computation of bioelectromagnetic fields in multi-compartment volume conductor models. Most BEM approaches assume that each compartment is in contact with at most one external compartment. In this work, I present a general surface integral equation and BEM discretization that remove this limitation and allow BEM modeling of general piecewise-homogeneous medium. The new integral equation allows positioning of field points at junctioned boundary of more than two compartments, enabling the use of linear collocation BEM in such a complex geometry. A modular BEM implementation is presented for linear collocation and Galerkin approaches, starting from standard formulation. The approach and resulting solver are verified in three ways, including comparison to finite element method (FEM). In a two-compartment split-sphere model with two spaced monopoles, the results obtained with high-resolution FEM and the BEMs were almost identical (relative difference < 0.003).

  17. On the removal of boundary errors caused by Runge-Kutta integration of non-linear partial differential equations

    Science.gov (United States)

    Abarbanel, Saul; Gottlieb, David; Carpenter, Mark H.

    1994-01-01

    It has been previously shown that the temporal integration of hyperbolic partial differential equations (PDE's) may, because of boundary conditions, lead to deterioration of accuracy of the solution. A procedure for removal of this error in the linear case has been established previously. In the present paper we consider hyperbolic (PDE's) (linear and non-linear) whose boundary treatment is done via the SAT-procedure. A methodology is present for recovery of the full order of accuracy, and has been applied to the case of a 4th order explicit finite difference scheme.

  18. Numerical experiments using CHIEF to treat the nonuniqueness in solving acoustic axisymmetric exterior problems via boundary integral equations

    Directory of Open Access Journals (Sweden)

    Adel A.K. Mohsen

    2010-07-01

    Full Text Available The problem of nonuniqueness (NU of the solution of exterior acoustic problems via boundary integral equations is discussed in this article. The efficient implementation of the CHIEF (Combined Helmholtz Integral Equations Formulation method to axisymmetric problems is studied. Interior axial fields are used to indicate the solution error and to select proper CHIEF points. The procedure makes full use of LU-decomposition as well as the forward solution derived in the solution. Implementations of the procedure for hard spheres are presented. Accurate results are obtained up to a normalised radius of ka = 20.983, using only one CHIEF point. The radiation from a uniformly vibrating sphere is also considered. Accurate results for ka up to 16.927 are obtained using two CHIEF points.

  19. Accurate simulation of two-dimensional optical microcavities with uniquely solvable boundary integral equations and trigonometric Galerkin discretization.

    Science.gov (United States)

    Boriskina, Svetlana V; Sewell, Phillip; Benson, Trevor M; Nosich, Alexander I

    2004-03-01

    A fast and accurate method is developed to compute the natural frequencies and scattering characteristics of arbitrary-shape two-dimensional dielectric resonators. The problem is formulated in terms of a uniquely solvable set of second-kind boundary integral equations and discretized by the Galerkin method with angular exponents as global test and trial functions. The log-singular term is extracted from one of the kernels, and closed-form expressions are derived for the main parts of all the integral operators. The resulting discrete scheme has a very high convergence rate. The method is used in the simulation of several optical microcavities for modern dense wavelength-division-multiplexed systems. PMID:15005404

  20. 3-D nonlinear force-free field reconstruction of solar active region 11158 by direct boundary integral equation

    CERN Document Server

    Wang, Rui; Tan, Baolin

    2013-01-01

    A 3-D coronal magnetic field is reconstructed for NOAA 11158 on Feb 14, 2011. A GPU-accelerated direct boundary integral equation (DBIE) method is implemented. This is about 1000 times faster than the original DBIE used on solar NLFFF modeling. Using the SDO/HMI vector magnetogram as the bottom boundary condition, the reconstructed magnetic field lines are compared with the projected EUV loop structures from different views three-dimensionally by SDO/AIA and STEREO A/B spacecraft simultaneously for the first time. They show very good agreement so that the topological configurations of the magnetic fields can be analyzed, thus its role in the flare process of the active region can be better understood. A quantitative comparison with some stereoscopically reconstructed coronal loops shows that the present averaged misalignment angles are at the same order as the state-of-the-art results obtained with reconstructed coronal loops as prescribed conditions and better than other NLFFF methods. It is found that the o...

  1. Rigorous electromagnetic analysis of two dimensional micro-axicon by boundary integral equations.

    Science.gov (United States)

    Lin, Jie; Tan, Jiubin; Liu, Jian; Liu, Shutian

    2009-02-01

    The focal performance of the micro-axicon and the Fresnel axicon (fraxicon) are investigated, for the first time, by the rigorous electromagnetic theory and boundary element method. The micro-axicon with different angle of apex and the fraxicon with various period and angle of apex are investigated. The dark segments of the fraxicon are explored numerically. Rigorous results of focal performance of the micro-axicon and the fraxicon are different from the results given by the approximation of geometrical optics and the scalar diffraction theory. The scattering effects are dominant in the fraxicon with small size of feature. It is expected that our study can provides very useful information in analyzing the axicon in optical trapping systems. PMID:19188975

  2. Three-dimensional layerwise modeling of layered media with boundary integral equations

    OpenAIRE

    Kokkinos, Filis-Triantaphyllos T.

    1995-01-01

    A hybrid method is presented for the analysis of layers, plates, and multi-layered systems consisting of isotropic and linear elastic materials. The problem is formulated for the general case of a multi-layered system using a total potential energy formulation and employing the layerwise laminate theory of Reddy. A one-dimensional finite element model is used for the analysis of the multi-layered system through its thickness, and integral Fourier transforms are used to obtai...

  3. A STUDY ON THE WEIGHT FUNCTION OF THE MOVING LEAST SQUARE APPROXIMATION IN THE LOCAL BOUNDARY INTEGRAL EQUATION METHOD

    Institute of Scientific and Technical Information of China (English)

    LongShuyao; HuDe'an

    2003-01-01

    The meshless method is a new numerical technique presented in recent years .It uses the moving least square (MLS) approximation as a shape function . The smoothness of the MLS approximation is determined by that of the basic function and of the weight function, and is mainly determined by that of the weight function. Therefore, the weight function greatly affects the accuracy of results obtained. Different kinds of weight functions, such as the spline function, the Gauss function and so on, are proposed recently by many researchers. In the present work, the features of various weight functions are illustrated through solving elasto-static problems using the local boundary integral equation method. The effect of various weight functions on the accuracy, convergence and stability of results obtained is also discussed. Examples show that the weight function proposed by Zhou Weiyuan and Gauss and the quartic spline weight function are better than the others if parameters c and a in Gauss and exponential weight functions are in the range of reasonable values, respectively, and the higher the smoothness of the weight function, the better the features of the solutions.

  4. Geophysical interpretation using integral equations

    CERN Document Server

    Eskola, L

    1992-01-01

    Along with the general development of numerical methods in pure and applied to apply integral equations to geophysical modelling has sciences, the ability improved considerably within the last thirty years or so. This is due to the successful derivation of integral equations that are applicable to the modelling of complex structures, and efficient numerical algorithms for their solution. A significant stimulus for this development has been the advent of fast digital computers. The purpose of this book is to give an idea of the principles by which boundary-value problems describing geophysical models can be converted into integral equations. The end results are the integral formulas and integral equations that form the theoretical framework for practical applications. The details of mathematical analysis have been kept to a minimum. Numerical algorithms are discussed only in connection with some illustrative examples involving well-documented numerical modelling results. The reader is assu­ med to have a back...

  5. Applied singular integral equations

    CERN Document Server

    Mandal, B N

    2011-01-01

    The book is devoted to varieties of linear singular integral equations, with special emphasis on their methods of solution. It introduces the singular integral equations and their applications to researchers as well as graduate students of this fascinating and growing branch of applied mathematics.

  6. Nonlinear streak computation using boundary region equations

    International Nuclear Information System (INIS)

    The boundary region equations (BREs) are applied for the simulation of the nonlinear evolution of a spanwise periodic array of streaks in a flat plate boundary layer. The well-known BRE formulation is obtained from the complete Navier–Stokes equations in the high Reynolds number limit, and provides the correct asymptotic description of three-dimensional boundary layer streaks. In this paper, a fast and robust streamwise marching scheme is introduced to perform their numerical integration. Typical streak computations present in the literature correspond to linear streaks or to small-amplitude nonlinear streaks computed using direct numerical simulation (DNS) or the nonlinear parabolized stability equations (PSEs). We use the BREs to numerically compute high-amplitude streaks, a method which requires much lower computational effort than DNS and does not have the consistency and convergence problems of the PSE. It is found that the flow configuration changes substantially as the amplitude of the streaks grows and the nonlinear effects come into play. The transversal motion (in the wall normal-streamwise plane) becomes more important and strongly distorts the streamwise velocity profiles, which end up being quite different from those of the linear case. We analyze in detail the resulting flow patterns for the nonlinearly saturated streaks and compare them with available experimental results. (paper)

  7. Reaction diffusion equations with boundary degeneracy

    Directory of Open Access Journals (Sweden)

    Huashui Zhan

    2016-03-01

    Full Text Available In this article, we consider the reaction diffusion equation $$ \\frac{\\partial u}{\\partial t} = \\Delta A(u,\\quad (x,t\\in \\Omega \\times (0,T, $$ with the homogeneous boundary condition. Inspired by the Fichera-Oleinik theory, if the equation is not only strongly degenerate in the interior of $\\Omega$, but also degenerate on the boundary, we show that the solution of the equation is free from any limitation of the boundary condition.

  8. Solution of the advection-diffusion equation for a nonhomogeneous and nonstationary Planetary Boundary Layer by GILTT (Generalized Integral Laplace Transform Technique)

    International Nuclear Information System (INIS)

    In this work is shown the solution of the advection-diffusion equation to simulate a pollutant dispersion in the Planetary Boundary Layer. The solution is obtained through of the GILTT (Generalized Integral Laplace Transform Technique) analytic method and of the numerical inversion Gauss Quadrature. The validity of the solution is proved using concentration obtained from the model with concentration obtained for Copenhagen experiment. In this comparison was utilized potential and logarithmic wind profile and eddy diffusivity derived by Degrazia et al (1997) [17] and (2002) [19]. The best results was using the potential wind profile and the eddy diffusivity derived by Degrazia et al (1997). The vertical velocity influence is shown in the plume behavior of the pollutant concentration. Moreover, the vertical and longitudinal velocity provided by Large Eddy Simulation (LES) was stood in the model to simulate the turbulent boundary layer more realistic, the result was satisfactory when compared with contained in the literature. (author)

  9. Soliton equations solved by the boundary CFT

    OpenAIRE

    SAITO, Satoru; Sato, Ryuichi

    2003-01-01

    Soliton equations are derived which characterize the boundary CFT a la Callan et al. Soliton fields of classical soliton equations are shown to appear as a neutral bound state of a pair of soliton fields of BCFT. One soliton amplitude under the influence of the boundary is calculated explicitly and is shown that it is frozen at the Dirichlet limit.

  10. Integral equations and computation problems

    International Nuclear Information System (INIS)

    Volterra's Integral Equations and Fredholm's Integral Equations of the second kind are discussed. Computational problems are found in the derivations and the computations. The theorem of the solution of the Fredholm's Integral Equation is discussed in detail. (author)

  11. BOUNDARY CONTROL OF MKDV-BURGERS EQUATION

    Institute of Scientific and Technical Information of China (English)

    TIAN Li-xin; ZHAO Zhi-feng; WANG Jing-feng

    2006-01-01

    The boundary control of MKdV-Burgers equation was considered by feedback control on the domain [0,1]. The existence of the solution of MKdV-Burgers equation with the feedback control law was proved. On the base, priori estimates for the solution was given. At last, the existence of the weak solution of MKdV-Burgers equation was proved and the global-exponential and asymptotic stability of the solution of MKdV-Burgers equation was given.

  12. On Certain Dual Integral Equations

    Directory of Open Access Journals (Sweden)

    R. S. Pathak

    1974-01-01

    Full Text Available Dual integral equations involving H-Functions have been solved by using the theory of Mellin transforms. The proof is analogous to that of Busbridge on solutions of dual integral equations involving Bessel functions.

  13. Functional equations for Feynman integrals

    International Nuclear Information System (INIS)

    New types of equations for Feynman integrals are found. It is shown that Feynman integrals satisfy functional equations connecting integrals with different kinematics. A regular method is proposed for obtaining such relations. The derivation of functional equations for one-loop two-, three- and four-point functions with arbitrary masses and external momenta is given. It is demonstrated that functional equations can be used for the analytic continuation of Feynman integrals to different kinematic domains

  14. Boundary value problems and partial differential equations

    CERN Document Server

    Powers, David L

    2005-01-01

    Boundary Value Problems is the leading text on boundary value problems and Fourier series. The author, David Powers, (Clarkson) has written a thorough, theoretical overview of solving boundary value problems involving partial differential equations by the methods of separation of variables. Professors and students agree that the author is a master at creating linear problems that adroitly illustrate the techniques of separation of variables used to solve science and engineering.* CD with animations and graphics of solutions, additional exercises and chapter review questions* Nearly 900 exercises ranging in difficulty* Many fully worked examples

  15. Boundary stabilization of wave equations with variable coefficients

    Institute of Scientific and Technical Information of China (English)

    FENG; Shaoji

    2001-01-01

    [1]Chen, G., Energy decay estimates and exact boundary value controllability for the wave equation in a bounded domain, J. Math. Pures. & Appl., 1979, 58: 249.[2]Komornik, V., Exact controllability and stabilization, Research in Applied Mathematics (Series Editors: Ciarlet, P. G., Lions, J.), New York: Masson/John Wiley copublication, 1994.[3]Komornik, V., Zuazua, E., A direct method for the boundary stabilization of the wave equation, J. Math. Pures. & Appl., 1990, 69: 33.[4]Lagnese, J., Decay of solutions of wave equations in a bounded region with boundary dissipation, J. Differential Equations, 1983, 50: 163.[5]Lasiecka, I., Triggiani, R., Uniform stabilization of the wave equation with Dirichlet or Neumann feedback control without geometrical conditions, Appl. Math. & Optim., 1992, 25: 189.[6]Wyler, A., Stability of wave equations with dissipative boundary conditions in a bounded domain, Differential and Integral Equations,1994, 7: 345.[7]Yao, P. F., On the observability inequality for exact controllability of wave equations with variable coefficients, SIAM J. Control & Optimization, 1999, 37, 5: 1568.[8]Wu, H., Shen, C. L., Yu, Y. L., Introduction to Riemannian Geometry (in Chinese), Beijing: Peking University Press, 1989.

  16. Integral equations and their applications

    CERN Document Server

    Rahman, M

    2007-01-01

    For many years, the subject of functional equations has held a prominent place in the attention of mathematicians. In more recent years this attention has been directed to a particular kind of functional equation, an integral equation, wherein the unknown function occurs under the integral sign. The study of this kind of equation is sometimes referred to as the inversion of a definite integral. While scientists and engineers can already choose from a number of books on integral equations, this new book encompasses recent developments including some preliminary backgrounds of formulations of in

  17. Radioisotope diffusion in grain textures by boundary integral method

    International Nuclear Information System (INIS)

    Aim of this contribution is to deal with radioisotope diffusion in grain texture by Boundary Integral method (BIM). Governing partial integral equation is transformed to an equivalent boundary integral equation, which is written in a discrete form and a system of linear algebraic equations is thus obtained. Advantage of BIM is that the system of equations is solved only for unknown values on the boundary. values in the domain are calculated explicitly in a down stream procedure. A given example indicates a good agreement with analytical results. (author)

  18. New Efficient Dynamic 3-D Boundary Integral Equation Method and Application to Non-Planar Fault Geometry Dipping in Elastic Half Space

    Science.gov (United States)

    Ando, R.

    2014-12-01

    The boundary integral equation method formulated in the real space and time domain (BIEM-ST) has been used as a powerful tool to analyze the earthquake rupture dynamics on non-planar faults. Generally, BIEM is more accurate than volumetric methods such as the finite difference method and the finite difference method. With the recent development of the high performance computing environment, the earthquake rupture simulation studies have been conducted considering three dimensional realistic fault geometry models. However, the utility of BIEM-ST has been limited due to its heavy computational demanding increased depending on square of time steps (N2), which was needed to evaluate the historic integration. While BIEM can be efficient with the spectral domain formulation, the applications of such a method are limited to planar fault cases. In this study, we propose a new method to reduce the calculation time of BIEM-ST to linear of time step (N) without degrading the accuracy in the 3 dimensional modeling space. We extends the method proposed earlier for the case of the 2 dimensional framework, applying the asymptotic expressions of the elasto-dynamic Green's functions. This method uses the physical nature of the stress Green's function as dividing the causality cone according to the distances from the wave-fronts. The scalability of this method is shown on the parallel computing environment of the distributed memory. We demonstrate the applicability to analyses of subduction earthquake cases, suffering long time from the numerical limitations of previously available BIEMs. We analyze the dynamic rupture processes on dipping reverse faults embed in a three dimensional elastic half space.

  19. Integral form of the radiation transfer equation

    International Nuclear Information System (INIS)

    The integral form of the radiation transfer equation is given in a non-scattering medium for which the source and absorption terms are known explicitly. The problem is solved for an one-dimensional, inhomogeneous, non stationary, non isotropic configuration, in cartesian and spherical coordinates for arbitrary initial and boundary conditions. The same problem is solved for a boundary condition that is given on a moving surface, then the three-dimensional problem is examined in cartesian coordinates

  20. Enclosing Solutions of Integral Equations

    DEFF Research Database (Denmark)

    Madsen, Kaj; NA NA NA Caprani, Ole; Stauning, Ole

    1996-01-01

    We present a method for enclosing the solution of an integral equation. It is assumed that a solution exists and that the corresponding integral operator T is a contraction near y. When solving the integral equation by iteration we obtain a result which is normally different from y because of...

  1. The Dirichlet problem for the two-dimensional Helmholtz equation for an open boundary

    Science.gov (United States)

    Hayashi, Y.

    1973-01-01

    Development of a complete theory of the two-dimensional Dirichlet problem for an open boundary. It is shown that the solution of the Dirichlet problem for an open boundary requires the solution of a Fredholm integral equation of the first kind. Although a Fredholm integral equation of the first kind usually has no solution if the kernel is continuous, owing to the logarithmic singularity of the kernel, the equation in this case is converted to a singular integral equation with a Cauchy kernel. It is proven that the homogeneous adjoint equation of the singular integral equation has no nonzero solution. By virtue of this result, and with the aid of an existence theorem known in the theory of singular integral equations, the existence of solutions of the singular integral equation, and then of the unique solution of the Fredholm integral equation of the first kind is proved.

  2. 弹性平面问题超奇异边界积分方程的近似列式与数值求解%Approximate Formulation and Numerical Solution for Hypersingular Boundary Integral Equations in Plane Elasticity

    Institute of Scientific and Technical Information of China (English)

    马杭; 黄兴

    2003-01-01

    Based on the fact that the singular boundary integrals in the sense of Cauchy principal value can be represented approxi-mately by the mean values of two companion nearly singular boundary integrals, a vary general approach was developed in the paper.In the approach, the approximate formulation before discretization was constructed to cope with the difficulties encountered in the cor-ner treatment in the formulations of hypersingular boundary integral equations. This makes it possible to solve the hypersingularboundary integral equation numerically in a non-regularized form and in a local manner by using conforming C0 quadratic boundary ele-ments and standard Gaussian quadratures similar to those employed in the conventional displacement-BIE formulations. The approxi-mate formulation is very convenient to use because the corner information is comprised naturally in the representations of those ap-proximate integrals. Numerical examples in plane elasticity show that with the present approach, the compatible or better results canbe achieved in comparison with those of the conventional BIE formulations.

  3. Singular-boundary reductions of type-Q ABS equations

    OpenAIRE

    Atkinson, James; Joshi, Nalini

    2011-01-01

    We study the fully discrete elliptic integrable model Q4 and its immediate trigonometric and rational counterparts (Q3, Q2 and Q1). Singular boundary problems for these equations are systematised in the framework of global singularity analysis. We introduce a technique to obtain solutions of such problems, in particular constructing the exact solution on a regular singularity-bounded strip. The solution technique is based on the multidimensional consistency and uses new insights into these eq...

  4. 非线性积分-微分方程周期边值问题的单调迭代方法%Monotonic Iteration of Nonlinear Integral-differential Equation with Periodic Boundary Condition

    Institute of Scientific and Technical Information of China (English)

    王家玉; 陈金兰

    2001-01-01

    By using the concept of abstract cone,improved comparable theoremand fixed point theorem of nonlinear operator and by discussing the periodic boundary condition of nonlinear integral-differential equation in Banach spaces,this paper gave two monotonic iterative sequences,also verified the minimum and maximum theorem for nonlinear integral-differential equation with periodic boundary condition in Banach spaces.%本文考虑Banach空间中非线性积分-微分方程的周期边值问题,利用抽象锥、推广了的比较定理及非线性算子的不动点定理,构造出两个单调迭代序列,证明了Banach空间中非线性积分-微分方程具有周期边值的最小解、最大解存在定理。

  5. Boundary differentiability for inhomogeneous infinity Laplace equations

    Directory of Open Access Journals (Sweden)

    Guanghao Hong

    2014-03-01

    Full Text Available We study the boundary regularity of the solutions to inhomogeneous infinity Laplace equations. We prove that if $u\\in C(\\bar{\\Omega}$ is a viscosity solution to $\\Delta_{\\infty}u:=\\sum_{i,j=1}^n u_{x_i}u_{x_j}u_{x_ix_j}=f$ with $f\\in C(\\Omega\\cap L^{\\infty}(\\Omega$ and for $x_0\\in \\partial\\Omega$ both $\\partial\\Omega$ and $g:=u|_{\\partial\\Omega}$ are differentiable at $x_0$, then u is differentiable at $x_0$.

  6. Boundary conditions for hyperbolic formulations of the Einstein equations

    OpenAIRE

    Frittelli, Simonetta; Gomez, Roberto

    2003-01-01

    In regards to the initial-boundary value problem of the Einstein equations, we argue that the projection of the Einstein equations along the normal to the boundary yields necessary and appropriate boundary conditions for a wide class of equivalent formulations. We explicitly show that this is so for the Einstein-Christoffel formulation of the Einstein equations in the case of spherical symmetry.

  7. Solvability of boundary value problems with Riemann-Stieltjes Δ-integral conditions for second-order dynamic equations on time scales at resonance

    OpenAIRE

    Li Yongkun; Shu Jiangye

    2011-01-01

    Abstract In this paper, by making use of the coincidence degree theory of Mawhin, the existence of the nontrivial solution for the boundary value problem with Riemann-Stieltjes Δ-integral conditions on time scales at resonance x Δ Δ ( t ) = f ( t , x ( t ) , x Δ ( t ) ) + e ( t ) , a . e . t ∈ [ 0 , T ] T , x Δ ( 0 ) = 0 , x ( T ) = ∫ 0 T x σ ( s ) Δ g ...

  8. Picard iterations of boundary-layer equations. [in singular-perturbation analysis of flightpath optimization problems

    Science.gov (United States)

    Ardema, M. D.; Yang, L.

    1985-01-01

    A method of solving the boundary-layer equations that arise in singular-perturbation analysis of flightpath optimization problems is presented. The method is based on Picard iterations of the integrated form of the equations and does not require iteration to find unknown boundary conditions. As an example, the method is used to develop a solution algorithm for the zero-order boundary-layer equations of the aircraft minimum-time-to-climb problem.

  9. Integral equation methods for electromagnetics

    CERN Document Server

    Volakis, John

    2012-01-01

    This text/reference is a detailed look at the development and use of integral equation methods for electromagnetic analysis, specifically for antennas and radar scattering. Developers and practitioners will appreciate the broad-based approach to understanding and utilizing integral equation methods and the unique coverage of historical developments that led to the current state-of-the-art. In contrast to existing books, Integral Equation Methods for Electromagnetics lays the groundwork in the initial chapters so students and basic users can solve simple problems and work their way up to the mo

  10. Recovering an obstacle using integral equations

    KAUST Repository

    Rundell, William

    2009-05-01

    We consider the inverse problem of recovering the shape, location and surface properties of an object where the surrounding medium is both conductive and homogeneous and we measure Cauchy data on an accessible part of the exterior boundary. It is assumed that the physical situation is modelled by harmonic functions and the boundary condition on the obstacle is one of Dirichlet type. The purpose of this paper is to answer some of the questions raised in a recent paper that introduced a nonlinear integral equation approach for the solution of this type of problem.

  11. Calculation of Turbulent Boundary Layers Using the Dissipation Integral Method

    Institute of Scientific and Technical Information of China (English)

    MatthiasBuschmann

    1999-01-01

    This paper gives an introduction into the dissipation integral method.The general integral equations for the three-dimensional case are derved.It is found that for a practical calculation algorithm the integral monentum equation and the integral energy equation are msot useful.Using Two different sets of mean velocity profiles the hyperbolical character of a dissipation integral method is shown.Test cases for two-and three-dimensional boundary layers are analysed and discussed.The paper concludes with a discussion of the advantages and limits of dissipation integral methods.

  12. LINEAR SINGULAR INTEGRAL EQUATION ON DOMAINS COMPOSED BY BALLS

    Institute of Scientific and Technical Information of China (English)

    2006-01-01

    For domains composed by balls in Cn, this paper studies the boundary behaviour of Cauchy type integrals with discrete holomorphic kernels and the corresponding linear singular integral equation on each piece of smooth lower dimensional edges on the boundary of the domain.

  13. Free Boundary Value Problems for Abstract Elliptic Equations and Applications

    Institute of Scientific and Technical Information of China (English)

    Veli SHAKHMUROV

    2011-01-01

    The free boundary value problems for elliptic differential-operator equations are studied.Several conditions for the uniform maximal regularity with respect to boundary parameters and the Fredholmness in abstract Lp-spaces are given.In application,the nonlocal free boundary problems for finite or infinite systems of elliptic and anisotropic type equations are studied.

  14. The free boundary Euler equations with large surface tension

    OpenAIRE

    Disconzi, Marcelo M.; Ebin, David G.

    2015-01-01

    We study the free boundary Euler equations with surface tension in three spatial dimensions, showing that the equations are well-posed if the coefficient of surface tension is positive. Then we prove that under natural assumptions, the solutions of the free boundary motion converge to solutions of the Euler equations in a domain with fixed boundary when the coefficient of surface tension tends to infinity.

  15. The free boundary Euler equations with large surface tension

    Science.gov (United States)

    Disconzi, Marcelo M.; Ebin, David G.

    2016-07-01

    We study the free boundary Euler equations with surface tension in three spatial dimensions, showing that the equations are well-posed if the coefficient of surface tension is positive. Then we prove that under natural assumptions, the solutions of the free boundary motion converge to solutions of the Euler equations in a domain with fixed boundary when the coefficient of surface tension tends to infinity.

  16. SOLVABILITY FOR NONLINEAR ELLIPTIC EQUATION WITH BOUNDARY PERTURBATION

    Institute of Scientific and Technical Information of China (English)

    2007-01-01

    The solvability of nonlinear elliptic equation with boundary perturbation is considered. The perturbed solution of original problem is obtained and the uniformly valid expansion of solution is proved.

  17. Treatment of domain integrals in boundary element methods

    Energy Technology Data Exchange (ETDEWEB)

    Nintcheu Fata, Sylvain [ORNL

    2012-01-01

    A systematic and rigorous technique to calculate domain integrals without a volume-fitted mesh has been developed and validated in the context of a boundary element approximation. In the proposed approach, a domain integral involving a continuous or weakly-singular integrand is first converted into a surface integral by means of straight-path integrals that intersect the underlying domain. Then, the resulting surface integral is carried out either via analytic integration over boundary elements or by use of standard quadrature rules. This domain-to-boundary integral transformation is derived from an extension of the fundamental theorem of calculus to higher dimension, and the divergence theorem. In establishing the method, it is shown that the higher-dimensional version of the first fundamental theorem of calculus corresponds to the well-known Poincare lemma. The proposed technique can be employed to evaluate integrals defined over simply- or multiply-connected domains with Lipschitz boundaries which are embedded in an Euclidean space of arbitrary but finite dimension. Combined with the singular treatment of surface integrals that is widely available in the literature, this approach can also be utilized to effectively deal with boundary-value problems involving non-homogeneous source terms by way of a collocation or a Galerkin boundary integral equation method using only the prescribed surface discretization. Sample problems associated with the three-dimensional Poisson equation and featuring the Newton potential are successfully solved by a constant element collocation method to validate this study.

  18. Boundary Integral Solutions to Three-Dimensional Unconfined Darcy's Flow

    Science.gov (United States)

    Lennon, Gerard P.; Liu, Philip L.-F.; Liggett, James A.

    1980-08-01

    The boundary integral equation method (BIEM) is used to solve three-dimensional potential flow problems in porous media. The problems considered here are time dependent and have a nonlinear boundary condition on the free surface. The entire boundary, including the moving free surface, discretized into linear finite elements for the purpose of evaluating the boundary integrals. The technique allows transient, three-dimensional problems to be solved with reasonable computational costs. Numerical examples include recharge through rectangular and circular areas and seepage flow from a surface pond. The examples are used to illustrate the method and show the nonlinear effects.

  19. Explicit Expressions for 3D Boundary Integrals in Potential Theory

    Energy Technology Data Exchange (ETDEWEB)

    Nintcheu Fata, Sylvain [ORNL

    2009-01-01

    On employing isoparametric, piecewise linear shape functions over a flat triangular domain, exact expressions are derived for all surface potentials involved in the numerical solution of three-dimensional singular and hyper-singular boundary integral equations of potential theory. These formulae, which are valid for an arbitrary source point in space, are represented as analytic expressions over the edges of the integration triangle. They can be used to solve integral equations defined on polygonal boundaries via the collocation method or may be utilized as analytic expressions for the inner integrals in the Galerkin technique. Also, the constant element approximation can be directly obtained with no extra effort. Sample problems solved by the collocation boundary element method for the Laplace equation are included to validate the proposed formulae.

  20. Compatibility between shape equation and boundary conditions of lipid membranes with free edges.

    Science.gov (United States)

    Tu, Z C

    2010-02-28

    Only some special open surfaces satisfying the shape equation of lipid membranes can be compatible with the boundary conditions. As a result of this compatibility, the first integral of the shape equation should vanish for axisymmetric lipid membranes, from which two theorems of nonexistence are verified: (i) there is no axisymmetric open membrane being a part of torus satisfying the shape equation; (ii) there is no axisymmetric open membrane being a part of a biconcave discodal surface satisfying the shape equation. Additionally, the shape equation is reduced to a second-order differential equation while the boundary conditions are reduced to two equations due to this compatibility. Numerical solutions to the reduced shape equation and boundary conditions agree well with the experimental data [A. Saitoh et al., Proc. Natl. Acad. Sci. U.S.A. 95, 1026 (1998)]. PMID:20192294

  1. The Pauli equation with complex boundary conditions

    CERN Document Server

    Kochan, D; Novak, R; Siegl, P

    2012-01-01

    We consider one-dimensional Pauli Hamiltonians in a bounded interval with possibly non-self-adjoint Robin-type boundary conditions. We study the influence of the spin-magnetic interaction on the interplay between the type of boundary conditions and the spectrum. A special attention is paid to PT-symmetric boundary conditions with the physical choice of the time-reversal operator T.

  2. Absorbing boundary conditions for second-order hyperbolic equations

    Science.gov (United States)

    Jiang, Hong; Wong, Yau Shu

    1990-01-01

    A uniform approach to construct absorbing artificial boundary conditions for second-order linear hyperbolic equations is proposed. The nonlocal boundary condition is given by a pseudodifferential operator that annihilates travelling waves. It is obtained through the dispersion relation of the differential equation by requiring that the initial-boundary value problem admits the wave solutions travelling in one direction only. Local approximation of this global boundary condition yields an nth-order differential operator. It is shown that the best approximations must be in the canonical forms which can be factorized into first-order operators. These boundary conditions are perfectly absorbing for wave packets propagating at certain group velocities. A hierarchy of absorbing boundary conditions is derived for transonic small perturbation equations of unsteady flows. These examples illustrate that the absorbing boundary conditions are easy to derive, and the effectiveness is demonstrated by the numerical experiments.

  3. Numerical integration of asymptotic solutions of ordinary differential equations

    Science.gov (United States)

    Thurston, Gaylen A.

    1989-01-01

    Classical asymptotic analysis of ordinary differential equations derives approximate solutions that are numerically stable. However, the analysis also leads to tedious expansions in powers of the relevant parameter for a particular problem. The expansions are replaced with integrals that can be evaluated by numerical integration. The resulting numerical solutions retain the linear independence that is the main advantage of asymptotic solutions. Examples, including the Falkner-Skan equation from laminar boundary layer theory, illustrate the method of asymptotic analysis with numerical integration.

  4. A New Algorithm for System of Integral Equations

    Directory of Open Access Journals (Sweden)

    Abdujabar Rasulov

    2014-01-01

    Full Text Available We develop a new algorithm to solve the system of integral equations. In this new method no need to use matrix weights. Beacause of it, we reduce computational complexity considerable. Using the new algorithm it is also possible to solve an initial boundary value problem for system of parabolic equations. To verify the efficiency, the results of computational experiments are given.

  5. Transmission problem for the Laplace equation and the integral equation method

    Czech Academy of Sciences Publication Activity Database

    Medková, Dagmar

    2012-01-01

    Roč. 387, č. 2 (2012), s. 837-843. ISSN 0022-247X Institutional research plan: CEZ:AV0Z10190503 Keywords : transmission problem * Laplace equation * boundary integral equation Subject RIV: BA - General Mathematics Impact factor: 1.050, year: 2012 http://www.sciencedirect.com/science/article/pii/S0022247X11008985

  6. The Pauli equation with complex boundary conditions

    International Nuclear Information System (INIS)

    We consider one-dimensional Pauli Hamiltonians in a bounded interval with possibly non-self-adjoint Robin-type boundary conditions. We study the influence of the spin–magnetic interaction on the interplay between the type of boundary conditions and the spectrum. Special attention is paid to PT-symmetric boundary conditions with the physical choice of the time-reversal operator T. This article is part of a special issue of Journal of Physics A: Mathematical and Theoretical devoted to ‘Quantum physics with non-Hermitian operators’. (paper)

  7. The Pauli equation with complex boundary conditions

    Science.gov (United States)

    Kochan, D.; Krejčiřík, D.; Novák, R.; Siegl, P.

    2012-11-01

    We consider one-dimensional Pauli Hamiltonians in a bounded interval with possibly non-self-adjoint Robin-type boundary conditions. We study the influence of the spin-magnetic interaction on the interplay between the type of boundary conditions and the spectrum. Special attention is paid to {PT}-symmetric boundary conditions with the physical choice of the time-reversal operator {T}. This article is part of a special issue of Journal of Physics A: Mathematical and Theoretical devoted to ‘Quantum physics with non-Hermitian operators’.

  8. Integrability of equations for soliton's eigenfunctions

    International Nuclear Information System (INIS)

    Eigenfunctions of the auxiliary linear problems for the soliton equations obey the nonlinear evolution equations. It is shown that these eigenfunction equations are integrable by the inverse spectral transform method. Eigenfunction equations are also the generating equations. Several (1+1) and (2+1) dimensional eigenfunction equations and their properties are considered. 11 refs

  9. The open boundary equation (discussion paper)

    NARCIS (Netherlands)

    Diederen, D.; Savenije, H.H.G.; Toffolon, M.

    2015-01-01

    We present a new equation describing the hydrodynamics in infinitely long tidal channels (i.e., no reflection) under the influence of oceanic forcing. The proposed equation is a simple relationship between partial derivatives of water level and velocity. It is formally derived for a progressive wave

  10. A LOCAL BOUNDARY INTEGRAL EQUATION[1mm]METHOD FOR THE ELASTICITY PROBLEM%弹性力学问题的局部边界积分方程方法

    Institute of Scientific and Technical Information of China (English)

    龙述尧; 许敬晓

    2000-01-01

    The basic concept and numerical implementation of a local boundaryintegral equation formulation for solving the elasticity problem havebeen presented in the present paper. It is a new truly meshless method,because the numerical implementation of the method leads to an efficientmeshless discrete model. The concept of a companion solution isintroduced, such that the traction terms would not appear in theintegrals over the local boundary after the modified integral kernel isused for all nodes whose local boundary s falls within theglobal boundary of the given problem; it uses the moving leastsquare approximations, and involves only boundary integration over alocal boundary centered at the node in question. It poses nodifficulties in satisfying essential boundary conditions, and leads to abanded and sparse system matrix. The undependence of the solution onthe size of the integral local boundary provides a great flexibilityin dealing with the numerical model of the elastic plane problems undervarious boundary conditions with arbitrary shapes. Convergence studiesin the numerical examples show that the present method possesses anexcellent rate of convergence and reasonably accurate results for boththe unknown displacement and strain energy, as the originalapproximated trial solutions have nice continuity and smoothness. Thenumerical results also show that using both linear and quadratic basesas well as spline and Gaussian weight functions in approximationfunctions can give quite accurate numerical results.   Compared with the conventional boundary element method based on the globalboundary integral equations, the present method is advantageous in thefollowing aspects:   i) No boundary and domain element needed to be constructed in the presentmethod, while it is necessary to discretize both the entire domain andits boundary for the conventional boundary element method in general. Thevolume and boundary integrals in the present method are evaluated onlyover small

  11. Integrable (2k)-Dimensional Hitchin Equations

    CERN Document Server

    Ward, R S

    2016-01-01

    This letter describes a completely-integrable system of Yang-Mills-Higgs equations which generalizes the Hitchin equations on a Riemann surface to arbitrary k-dimensional complex manifolds. The system arises as a dimensional reduction of a set of integrable Yang-Mills equations in 4k real dimensions. Our integrable system implies other generalizations such as the Simpson equations and the non-abelian Seiberg-Witten equations. Some simple solutions in the k=2 case are described.

  12. From the Boltzmann Equation to the Euler Equations in the Presence of Boundaries

    CERN Document Server

    Golse, François

    2011-01-01

    The fluid dynamic limit of the Boltzmann equation leading to the Euler equations for an incompressible fluid with constant density in the presence of material boundaries shares some important features with the better known inviscid limit of the Navier-Stokes equations. The present paper slightly extends recent results from [C. Bardos, F. Golse, L. Paillard, Comm. Math. Sci., 10 (2012), 159--190] to the case of boundary conditions for the Boltzmann equation more general than Maxwell's accomodation condition.

  13. From the Boltzmann Equation to the Euler Equations in the Presence of Boundaries

    OpenAIRE

    Golse, François

    2011-01-01

    The fluid dynamic limit of the Boltzmann equation leading to the Euler equations for an incompressible fluid with constant density in the presence of material boundaries shares some important features with the better known inviscid limit of the Navier-Stokes equations. The present paper slightly extends recent results from [C. Bardos, F. Golse, L. Paillard, Comm. Math. Sci., 10 (2012), 159--190] to the case of boundary conditions for the Boltzmann equation more general than Maxwell's accomoda...

  14. Boundary value problemfor multidimensional fractional advection-dispersion equation

    Directory of Open Access Journals (Sweden)

    Khasambiev Mokhammad Vakhaevich

    2015-05-01

    Full Text Available In recent time there is a very great interest in the study of differential equations of fractional order, in which the unknown function is under the symbol of fractional derivative. It is due to the development of the theory of fractional integro-differential theory and application of it in different fields.The fractional integrals and derivatives of fractional integro-differential equations are widely used in modern investigations of theoretical physics, mechanics, and applied mathematics. The fractional calculus is a very powerful tool for describing physical systems, which have a memory and are non-local. Many processes in complex systems have nonlocality and long-time memory. Fractional integral operators and fractional differential operators allow describing some of these properties. The use of the fractional calculus will be helpful for obtaining the dynamical models, in which integro-differential operators describe power long-time memory by time and coordinates, and three-dimensional nonlocality for complex medium and processes.Differential equations of fractional order appear when we use fractal conception in physics of the condensed medium. The transfer, described by the operator with fractional derivatives at a long distance from the sources, leads to other behavior of relatively small concentrations as compared with classic diffusion. This fact redefines the existing ideas about safety, based on the ideas on exponential velocity of damping. Fractional calculus in the fractal theory and the systems with memory have the same importance as the classic analysis in mechanics of continuous medium.In recent years, the application of fractional derivatives for describing and studying the physical processes of stochastic transfer is very popular too. Many problems of filtration of liquids in fractal (high porous medium lead to the need to study boundary value problems for partial differential equations in fractional order.In this paper the

  15. Stochastic integration and differential equations

    CERN Document Server

    Protter, Philip E

    2003-01-01

    It has been 15 years since the first edition of Stochastic Integration and Differential Equations, A New Approach appeared, and in those years many other texts on the same subject have been published, often with connections to applications, especially mathematical finance. Yet in spite of the apparent simplicity of approach, none of these books has used the functional analytic method of presenting semimartingales and stochastic integration. Thus a 2nd edition seems worthwhile and timely, though it is no longer appropriate to call it "a new approach". The new edition has several significant changes, most prominently the addition of exercises for solution. These are intended to supplement the text, but lemmas needed in a proof are never relegated to the exercises. Many of the exercises have been tested by graduate students at Purdue and Cornell Universities. Chapter 3 has been completely redone, with a new, more intuitive and simultaneously elementary proof of the fundamental Doob-Meyer decomposition theorem, t...

  16. Boundary Element Method with Non—overlapping Domain Decomposition for Diffusion Equation

    Institute of Scientific and Technical Information of China (English)

    ZHUJialin; ZHANGTaiping

    2002-01-01

    A boundary element method based on non-overlapping domain decomposition method to solve the time-dependent diffusion equations is presented.The time-dependent fundamental solution is used in the formulation of boundary integrals and the time integratioin process always restarts from the initial time condition.The process of replacing the interface values,which needs a summation of boundary integrals related to the boundary values at previous time steps can be treated in parallel parallel iterative procedure,Numerical experiments demonstrate that the implementation of the present alogrithm is efficient.

  17. Analysis of Laminar Boundary Layer Equations

    Directory of Open Access Journals (Sweden)

    R. Yesman

    2012-01-01

    Full Text Available The paper proposes methodology for analysis and calculation of laminar fluid flow processes in a boundary layer.The presented dependences can be used for practical calculations while power carriers of various application are moving in the channels of heat and power devices. 

  18. Homogenization of the evolution Stokes equation in a perforated domain with a stochastic Fourier boundary condition

    KAUST Repository

    Bessaih, Hakima

    2015-04-01

    The evolution Stokes equation in a domain containing periodically distributed obstacles subject to Fourier boundary condition on the boundaries is considered. We assume that the dynamic is driven by a stochastic perturbation on the interior of the domain and another stochastic perturbation on the boundaries of the obstacles. We represent the solid obstacles by holes in the fluid domain. The macroscopic (homogenized) equation is derived as another stochastic partial differential equation, defined in the whole non perforated domain. Here, the initial stochastic perturbation on the boundary becomes part of the homogenized equation as another stochastic force. We use the twoscale convergence method after extending the solution with 0 in the holes to pass to the limit. By Itô stochastic calculus, we get uniform estimates on the solution in appropriate spaces. In order to pass to the limit on the boundary integrals, we rewrite them in terms of integrals in the whole domain. In particular, for the stochastic integral on the boundary, we combine the previous idea of rewriting it on the whole domain with the assumption that the Brownian motion is of trace class. Due to the particular boundary condition dealt with, we get that the solution of the stochastic homogenized equation is not divergence free. However, it is coupled with the cell problem that has a divergence free solution. This paper represents an extension of the results of Duan and Wang (Comm. Math. Phys. 275:1508-1527, 2007), where a reaction diffusion equation with a dynamical boundary condition with a noise source term on both the interior of the domain and on the boundary was studied, and through a tightness argument and a pointwise two scale convergence method the homogenized equation was derived. © American Institute of Mathematical Sciences.

  19. Boundary value problems of discrete generalized Emden-Fowler equation

    Institute of Scientific and Technical Information of China (English)

    YU; Jianshe; GUO; Zhiming

    2006-01-01

    By using the critical point theory, some sufficient conditions for the existence of the solutions to the boundary value problems of a discrete generalized Emden-Fowler equation are obtained. In a special case, a sharp condition is obtained for the existence of the boundary value problems of the above equation. For a linear case, by the discrete variational theory, a necessary and sufficient condition for the existence, uniqueness and multiplicity of the solutions is also established.

  20. Integral Transform Approach to Generalized Tricomi Equations

    OpenAIRE

    Yagdjian, Karen

    2014-01-01

    We present some integral transform that allows to obtain solutions of the generalized Tricomi equation from solutions of a simpler equation. We used in [13,14],[41]-[46] the particular version of this transform in order to investigate in a unified way several equations such as the linear and semilinear Tricomi equations, Gellerstedt equation, the wave equation in Einstein-de Sitter spacetime, the wave and the Klein-Gordon equations in the de Sitter and anti-de Sitter spacetimes.

  1. Linear integral equations and soliton systems

    International Nuclear Information System (INIS)

    A study is presented of classical integrable dynamical systems in one temporal and one spatial dimension. The direct linearizations are given of several nonlinear partial differential equations, for example the Korteweg-de Vries equation, the modified Korteweg-de Vries equation, the sine-Gordon equation, the nonlinear Schroedinger equation, and the equation of motion for the isotropic Heisenberg spin chain; the author also discusses several relations between these equations. The Baecklund transformations of these partial differential equations are treated on the basis of a singular transformation of the measure (or equivalently of the plane-wave factor) occurring in the corresponding linear integral equations, and the Baecklund transformations are used to derive the direct linearization of a chain of so-called modified partial differential equations. Finally it is shown that the singular linear integral equations lead in a natural way to the direct linearizations of various nonlinear difference-difference equations. (Auth.)

  2. Multiple integral representation for the trigonometric SOS model with domain wall boundaries

    International Nuclear Information System (INIS)

    Using the dynamical Yang-Baxter algebra we derive a functional equation for the partition function of the trigonometric SOS model with domain wall boundary conditions. The solution of the equation is given in terms of a multiple contour integral.

  3. Multiple integral representation for the trigonometric SOS model with domain wall boundaries

    Science.gov (United States)

    Galleas, W.

    2012-05-01

    Using the dynamical Yang-Baxter algebra we derive a functional equation for the partition function of the trigonometric SOS model with domain wall boundary conditions. The solution of the equation is given in terms of a multiple contour integral.

  4. On a stochastic Burgers equation with Dirichlet boundary conditions

    Directory of Open Access Journals (Sweden)

    Ekaterina T. Kolkovska

    2003-01-01

    Full Text Available We consider the one-dimensional Burgers equation perturbed by a white noise term with Dirichlet boundary conditions and a non-Lipschitz coefficient. We obtain existence of a weak solution proving tightness for a sequence of polygonal approximations for the equation and solving a martingale problem for the weak limit.

  5. Thermodynamically admissible boundary conditions for the regularized 13 moment equations

    International Nuclear Information System (INIS)

    A phenomenological approach to the boundary conditions for linearized R13 equations is derived using the second law of thermodynamics. The phenomenological coefficients appearing in the boundary conditions are calculated by comparing the slip, jump, and thermal creep coefficients with linearized Boltzmann solutions for Maxwell’s accommodation model for different values of the accommodation coefficient. For this, the linearized R13 equations are solved for viscous slip, thermal creep, and temperature jump problems and the results are compared to the solutions of the linearized Boltzmann equation. The influence of different collision models (hard-sphere, Bhatnagar–Gross–Krook, and Maxwell molecules) and accommodation coefficients on the phenomenological coefficients is studied

  6. Thermodynamically admissible boundary conditions for the regularized 13 moment equations

    Science.gov (United States)

    Rana, Anirudh Singh; Struchtrup, Henning

    2016-02-01

    A phenomenological approach to the boundary conditions for linearized R13 equations is derived using the second law of thermodynamics. The phenomenological coefficients appearing in the boundary conditions are calculated by comparing the slip, jump, and thermal creep coefficients with linearized Boltzmann solutions for Maxwell's accommodation model for different values of the accommodation coefficient. For this, the linearized R13 equations are solved for viscous slip, thermal creep, and temperature jump problems and the results are compared to the solutions of the linearized Boltzmann equation. The influence of different collision models (hard-sphere, Bhatnagar-Gross-Krook, and Maxwell molecules) and accommodation coefficients on the phenomenological coefficients is studied.

  7. Thermodynamically admissible boundary conditions for the regularized 13 moment equations

    Energy Technology Data Exchange (ETDEWEB)

    Rana, Anirudh Singh, E-mail: anirudh@uvic.ca [Department of Mechanical and Aerospace Engineering, Gyeongsang National University, Jinju, Gyeongnam 52828 (Korea, Republic of); Struchtrup, Henning, E-mail: struchtr@uvic.ca [Department of Mechanical Engineering, University of Victoria, Victoria, British Columbia V8W 2Y2 (Canada)

    2016-02-15

    A phenomenological approach to the boundary conditions for linearized R13 equations is derived using the second law of thermodynamics. The phenomenological coefficients appearing in the boundary conditions are calculated by comparing the slip, jump, and thermal creep coefficients with linearized Boltzmann solutions for Maxwell’s accommodation model for different values of the accommodation coefficient. For this, the linearized R13 equations are solved for viscous slip, thermal creep, and temperature jump problems and the results are compared to the solutions of the linearized Boltzmann equation. The influence of different collision models (hard-sphere, Bhatnagar–Gross–Krook, and Maxwell molecules) and accommodation coefficients on the phenomenological coefficients is studied.

  8. Path Integral Solution of Linear Second Order Partial Differential Equations I. The General Construction

    OpenAIRE

    LaChapelle, J.

    2004-01-01

    A path integral is presented that solves a general class of linear second order partial differential equations with Dirichlet/Neumann boundary conditions. Elementary kernels are constructed for both Dirichlet and Neumann boundary conditions. The general solution can be specialized to solve elliptic, parabolic, and hyperbolic partial differential equations with boundary conditions. This extends the well-known path integral solution of the Schr\\"{o}dinger/diffusion equation in unbounded space. ...

  9. First integrals of difference Hamiltonian equations

    International Nuclear Information System (INIS)

    In the present paper, the well-known Noether's identity, which represents the connection between symmetries and first integrals of Euler-Lagrange equations, is rewritten in terms of the Hamiltonian function. This approach, based on the Hamiltonian identity, provides a simple and clear way to find first integrals of canonical Hamiltonian equations without integration. A discrete analog of the Hamiltonian identity is developed. It leads to a connection between symmetries and first integrals of difference Hamiltonian equations that can be used to conserve the structural properties of Hamiltonian equations under discretizaton. The results are illustrated by a number of examples for both continuous and difference Hamiltonian equations.

  10. On Impulsive Boundary Value Problems of Fractional Differential Equations with Irregular Boundary Conditions

    Directory of Open Access Journals (Sweden)

    Guotao Wang

    2012-01-01

    Full Text Available We study nonlinear impulsive differential equations of fractional order with irregular boundary conditions. Some existence and uniqueness results are obtained by applying standard fixed-point theorems. For illustration of the results, some examples are discussed.

  11. Solitons induced by boundary conditions from the Boussinesq equation

    Science.gov (United States)

    Chou, Ru Ling; Chu, C. K.

    1990-01-01

    The behavior of solitons induced by boundary excitation is investigated at various time-dependent conditions and different unperturbed water depths, using the Korteweg-de Vries (KdV) equation. Then, solitons induced from Boussinesq equations under similar conditions were studied, making it possible to remove the restriction in the KdV equation and to treat soliton head-on collisions (as well as overtaking collisions) and reflections. It is found that the results obtained from the KdV and the Boussinesq equations are in good agreement.

  12. SINGULARLY PERTURBED BOUNDARY VALUE PROBLEMS FOR SEMI-LINEAR RETARDED DIFFERENTIAL EQUATIONS WITH NONLINEAR BOUNDARY CONDITIONS

    Institute of Scientific and Technical Information of China (English)

    任景莉; 葛渭高

    2003-01-01

    A boundary value problems f or functional differenatial equations, with nonlinear boundary condition, is studied by the theorem of differential inequality. Using new method to construct the upper solution and lower solution, sufficient conditions for the existence of the problems' solution are established. A uniformly valid asymptotic expansions of the solution is also given.

  13. Periodic solutions to nonlinear equations with oblique boundary conditions

    OpenAIRE

    Allergretto, Walter; Papini, Duccio

    2012-01-01

    We study the existence of positive periodic solutions to nonlinear elliptic and parabolic equations with oblique and dynamical boundary conditions and non-local terms. The results are obtained through fixed point theory, topological degree methods and properties of related linear elliptic problems with natural boundary conditions and possibly non-symmetric principal part. As immediate consequences, we also obtain estimates on the principal eigenvalue for non-symmetric elliptic ...

  14. Distribution theory for Schrödinger’s integral equation

    Energy Technology Data Exchange (ETDEWEB)

    Lange, Rutger-Jan, E-mail: rutger-jan.lange@cantab.net [VU University Amsterdam, 1081 HV Amsterdam (Netherlands)

    2015-12-15

    Much of the literature on point interactions in quantum mechanics has focused on the differential form of Schrödinger’s equation. This paper, in contrast, investigates the integral form of Schrödinger’s equation. While both forms are known to be equivalent for smooth potentials, this is not true for distributional potentials. Here, we assume that the potential is given by a distribution defined on the space of discontinuous test functions. First, by using Schrödinger’s integral equation, we confirm a seminal result by Kurasov, which was originally obtained in the context of Schrödinger’s differential equation. This hints at a possible deeper connection between both forms of the equation. We also sketch a generalisation of Kurasov’s [J. Math. Anal. Appl. 201(1), 297–323 (1996)] result to hypersurfaces. Second, we derive a new closed-form solution to Schrödinger’s integral equation with a delta prime potential. This potential has attracted considerable attention, including some controversy. Interestingly, the derived propagator satisfies boundary conditions that were previously derived using Schrödinger’s differential equation. Third, we derive boundary conditions for “super-singular” potentials given by higher-order derivatives of the delta potential. These boundary conditions cannot be incorporated into the normal framework of self-adjoint extensions. We show that the boundary conditions depend on the energy of the solution and that probability is conserved. This paper thereby confirms several seminal results and derives some new ones. In sum, it shows that Schrödinger’s integral equation is a viable tool for studying singular interactions in quantum mechanics.

  15. Distribution theory for Schrödinger’s integral equation

    International Nuclear Information System (INIS)

    Much of the literature on point interactions in quantum mechanics has focused on the differential form of Schrödinger’s equation. This paper, in contrast, investigates the integral form of Schrödinger’s equation. While both forms are known to be equivalent for smooth potentials, this is not true for distributional potentials. Here, we assume that the potential is given by a distribution defined on the space of discontinuous test functions. First, by using Schrödinger’s integral equation, we confirm a seminal result by Kurasov, which was originally obtained in the context of Schrödinger’s differential equation. This hints at a possible deeper connection between both forms of the equation. We also sketch a generalisation of Kurasov’s [J. Math. Anal. Appl. 201(1), 297–323 (1996)] result to hypersurfaces. Second, we derive a new closed-form solution to Schrödinger’s integral equation with a delta prime potential. This potential has attracted considerable attention, including some controversy. Interestingly, the derived propagator satisfies boundary conditions that were previously derived using Schrödinger’s differential equation. Third, we derive boundary conditions for “super-singular” potentials given by higher-order derivatives of the delta potential. These boundary conditions cannot be incorporated into the normal framework of self-adjoint extensions. We show that the boundary conditions depend on the energy of the solution and that probability is conserved. This paper thereby confirms several seminal results and derives some new ones. In sum, it shows that Schrödinger’s integral equation is a viable tool for studying singular interactions in quantum mechanics

  16. Numerical Quadrature of Periodic Singular Integral Equations

    DEFF Research Database (Denmark)

    Krenk, Steen

    1978-01-01

    This paper presents quadrature formulae for the numerical integration of a singular integral equation with Hilbert kernel. The formulae are based on trigonometric interpolation. By integration a quadrature formula for an integral with a logarithmic singularity is obtained. Finally it is demonstra......This paper presents quadrature formulae for the numerical integration of a singular integral equation with Hilbert kernel. The formulae are based on trigonometric interpolation. By integration a quadrature formula for an integral with a logarithmic singularity is obtained. Finally it is...

  17. Traveling waves for a boundary reaction-diffusion equation

    CERN Document Server

    Caffarelli, L; Sire, Y

    2011-01-01

    We prove the existence of a traveling wave solution for a boundary reaction diffusion equation when the reaction term is the combustion nonlinearity with ignition temperature. A key role in the proof is plaid by an explicit formula for traveling wave solutions of a free boundary problem obtained as singular limit for the reaction-diffusion equation (the so-called high energy activation energy limit). This explicit formula, which is interesting in itself, also allows us to get an estimate on the decay at infinity of the traveling wave (which turns out to be faster than the usual exponential decay).

  18. Linear boundary value problems for differential algebraic equations

    OpenAIRE

    Balla, Katalin; März, Roswitha

    2003-01-01

    By the use of the corresponding shift matrix, the paper gives a criterion for the unique solvability of linear boundary value problems posed for linear differential algebraic equations up to index 2 with well-matched leading coefficients. The solution is constructed by a proper Green function. Another characterization of the solutions is based upon the description of arbitrary affine linear subspaces of solutions to linear differential algebraic equations in terms of solutions to the adjoint ...

  19. The Einstein Constraint Equations on Compact Manifolds with Boundary

    CERN Document Server

    Dilts, James

    2013-01-01

    We continue the study of the Einstein constraint equations on compact manifolds with boundary initiated by Holst and Tsogtgerel in their 2013 paper. In particular, we consider the full system and prove existence of solutions in both the near-CMC and far-from-CMC (for Yamabe positive metrics) cases. We also prove analogues many of the useful inequalities and results in previous "limit equation" papers by Dahl, Gicquaud, Humbert and others.

  20. Numerical solutions of telegraph equations with the Dirichlet boundary condition

    Science.gov (United States)

    Ashyralyev, Allaberen; Turkcan, Kadriye Tuba; Koksal, Mehmet Emir

    2016-08-01

    In this study, the Cauchy problem for telegraph equations in a Hilbert space is considered. Stability estimates for the solution of this problem are presented. The third order of accuracy difference scheme is constructed for approximate solutions of the problem. Stability estimates for the solution of this difference scheme are established. As a test problem to support theoretical results, one-dimensional telegraph equation with the Dirichlet boundary condition is considered. Numerical solutions of this equation are obtained by first, second and third order of accuracy difference schemes.

  1. Navier-Stokes equation with slip-like boundary condition

    OpenAIRE

    Ken'ichi Hashizume; Tetsuya Koyama; Mitsuharu Otani

    2009-01-01

    The aim of this note is to investigate a time-discretized 2-dimensional Navier-Stokes equation with a slip-like boundary condition, which arises in the melting ice problem. We prove the existence and uniqueness of a weak solution.

  2. Positive Solutions for Nonlinear Differential Equations with Periodic Boundary Condition

    Directory of Open Access Journals (Sweden)

    Shengjun Li

    2012-01-01

    Full Text Available We study the existence of positive solutions for second-order nonlinear differential equations with nonseparated boundary conditions. Our nonlinearity may be singular in its dependent variable. The proof of the main result relies on a nonlinear alternative principle of Leray-Schauder. Recent results in the literature are generalized and significantly improved.

  3. QUENCHING ON BOUNDARY TO THE NEWTON FILTRATION EQUATION (Ⅰ)

    Institute of Scientific and Technical Information of China (English)

    段志文; 谢春红; 卢伟明

    2003-01-01

    This paper discusses the global existence and quenching of the solution to the Newton filtration equation with the nonlinear boundary condition.The authors also discuss the profile of the quenching solution in the quenching time and obtain the quenching rate of the quenching solution.

  4. Algebraic Bethe Ansatz for O(2N) sigma models with integrable diagonal boundaries

    CERN Document Server

    Gombor, Tamas

    2015-01-01

    The finite volume problem of O(2N) sigma models with integrable diagonal boundaries on a finite interval is investigated. The double row transfer matrix is diagonalized by Algebraic Bethe Ansatz. The boundary Bethe Yang equations for the particle rapidities and the accompanying Bethe Ansatz equations are derived.

  5. An integral equation representation approach for valuing Russian options with a finite time horizon

    Science.gov (United States)

    Jeon, Junkee; Han, Heejae; Kim, Hyeonuk; Kang, Myungjoo

    2016-07-01

    In this paper, we first describe a general solution for the inhomogeneous Black-Scholes partial differential equation with mixed boundary conditions using Mellin transform techniques. Since Russian options with a finite time horizon are usually formulated into the inhomogeneous free-boundary Black-Scholes partial differential equation with a mixed boundary condition, we apply our method to Russian options and derive an integral equation satisfied by Russian options with a finite time horizon. Furthermore, we present some numerical solutions and plots of the integral equation using recursive integration methods and demonstrate the computational accuracy and efficiency of our method compared to other competing approaches.

  6. A CLASS OF NONLINEAR SINGULARLY PERTURBED INITIAL BOUNDARY VALUE PROBLEMS FOR REACTION DIFFUSION EQUATIONS WITH BOUNDARY PERTURBATION

    Institute of Scientific and Technical Information of China (English)

    Mo Jiaqi

    2007-01-01

    A class of nonlinear initial boundary value problems for reaction diffusion equations with boundary perturbation is considered. Under suitable conditions and using the theory of differential inequalities the asymptotic solution of the initial boundary value problems is studied.

  7. On the wave equation with semilinear porous acoustic boundary conditions

    KAUST Repository

    Graber, Philip Jameson

    2012-05-01

    The goal of this work is to study a model of the wave equation with semilinear porous acoustic boundary conditions with nonlinear boundary/interior sources and a nonlinear boundary/interior damping. First, applying the nonlinear semigroup theory, we show the existence and uniqueness of local in time solutions. The main difficulty in proving the local existence result is that the Neumann boundary conditions experience loss of regularity due to boundary sources. Using an approximation method involving truncated sources and adapting the ideas in Lasiecka and Tataru (1993) [28], we show that the existence of solutions can still be obtained. Second, we prove that under some restrictions on the source terms, then the local solution can be extended to be global in time. In addition, it has been shown that the decay rates of the solution are given implicitly as solutions to a first order ODE and depends on the behavior of the damping terms. In several situations, the obtained ODE can be easily solved and the decay rates can be given explicitly. Third, we show that under some restrictions on the initial data and if the interior source dominates the interior damping term and if the boundary source dominates the boundary damping, then the solution ceases to exists and blows up in finite time. Moreover, in either the absence of the interior source or the boundary source, then we prove that the solution is unbounded and grows as an exponential function. © 2012 Elsevier Inc.

  8. Gravity actions, boundary terms and second-order field equations

    International Nuclear Information System (INIS)

    The Lovelock action is a natural generalization of the Einstein-Hilbert gravity action in higher (more than four) dimensions. Certain boundary terms have to be added to this action because of various reasons which we briefly review. Dimensional reduction of the Lovelock action leads to a generalized Einstein-Yang-Mills (+ scalar field or σ-models) action in lower dimensions. Since the field equations derived from the Lovelock action are of second order only, the same holds for the dimensionally reduced action. The necessary boundary terms for the latter action are then obtained by dimensional reduction of the boundary terms in the Lovelock action. We consider reduction from n to n-1 dimensions in some detail and obtain generalizations of Horndeski's non-minimally coupled Einstein-Maxwell action with second-order field equations in more than four dimensions. (orig./HSI)

  9. Solving wave equation with spectral methods and nonreflecting boundary conditions

    CERN Document Server

    Novák, J; Novak, Jerome; Bonazzola, Silvano

    2002-01-01

    A multidomain spectral method for solving wave equations is presented. This method relies on the expansion of functions on basis of spherical harmonics $(Y_l^m(\\theta, \\phi))$ for the angular dependence and of Chebyshev polynomials $T_n(x)$ for the radial part. The spherical domains consist of shells surrounding a nucleus and cover the space up to a finite radius $R$ at which boundary conditions are imposed. Time derivatives are estimated using standard finite-differences second order schemes, which are chosen to be implicit to allow for (almost) any size of time-step. Emphasis is put on the implementation of absorbing boundary conditions that allow for the numerical boundary to be completely transparent to the physical wave. This is done using a multipolar expansion of an exact boundary condition for outgoing waves, which is truncated at some point. Using an auxiliary function, which is solution of a wave equation on the sphere defining the outer boundary of the numerical grid, the absorbing boundary conditi...

  10. Partial differential equations & boundary value problems with Maple

    CERN Document Server

    Articolo, George A

    2009-01-01

    Partial Differential Equations and Boundary Value Problems with Maple presents all of the material normally covered in a standard course on partial differential equations, while focusing on the natural union between this material and the powerful computational software, Maple. The Maple commands are so intuitive and easy to learn, students can learn what they need to know about the software in a matter of hours- an investment that provides substantial returns. Maple''s animation capabilities allow students and practitioners to see real-time displays of the solutions of partial differential equations.  Maple files can be found on the books website. Ancillary list: Maple files- http://www.elsevierdirect.com/companion.jsp?ISBN=9780123747327  Provides a quick overview of the software w/simple commands needed to get startedIncludes review material on linear algebra and Ordinary Differential equations, and their contribution in solving partial differential equationsIncorporates an early introduction to Sturm-L...

  11. Integral equations for the electromagnetic field in dielectrics

    Science.gov (United States)

    Mostowski, Jan; Załuska-Kotur, Magdalena A.

    2016-09-01

    We study static the electric field and electromagnetic waves in dielectric media. In contrast to the standard approach, we use, formulate and solve integral equations for the field. We discuss the case of an electrostatic field of a point charge placed inside a dielectric; the integral equation approach allows us to find and interpret the dielectric constant in terms of molecular polarizability. Next we discuss propagation of electromagnetic waves using the same integral equation approach. We derive the dispersion relation and find the reflection and transmission coefficients at the boundary between the vacuum and the dielectric. The present approach supplements the standard approach based on macroscopic Maxwell equations and contributes to better a understanding of some electromagnetic effects.

  12. A kernel-free boundary integral method for elliptic boundary value problems

    Science.gov (United States)

    Ying, Wenjun; Henriquez, Craig S.

    2007-12-01

    This paper presents a class of kernel-free boundary integral (KFBI) methods for general elliptic boundary value problems (BVPs). The boundary integral equations reformulated from the BVPs are solved iteratively with the GMRES method. During the iteration, the boundary and volume integrals involving Green's functions are approximated by structured grid-based numerical solutions, which avoids the need to know the analytical expressions of Green's functions. The KFBI method assumes that the larger regular domain, which embeds the original complex domain, can be easily partitioned into a hierarchy of structured grids so that fast elliptic solvers such as the fast Fourier transform (FFT) based Poisson/Helmholtz solvers or those based on geometric multigrid iterations are applicable. The structured grid-based solutions are obtained with standard finite difference method (FDM) or finite element method (FEM), where the right hand side of the resulting linear system is appropriately modified at irregular grid nodes to recover the formal accuracy of the underlying numerical scheme. Numerical results demonstrating the efficiency and accuracy of the KFBI methods are presented. It is observed that the number of GMRES iterations used by the method for solving isotropic and moderately anisotropic BVPs is independent of the sizes of the grids that are employed to approximate the boundary and volume integrals. With the standard second-order FEMs and FDMs, the KFBI method shows a second-order convergence rate in accuracy for all of the tested Dirichlet/Neumann BVPs when the anisotropy of the diffusion tensor is not too strong.

  13. The determination of an unknown boundary condition in a fractional diffusion equation

    KAUST Repository

    Rundell, William

    2013-07-01

    In this article we consider an inverse boundary problem, in which the unknown boundary function ∂u/∂v = f(u) is to be determined from overposed data in a time-fractional diffusion equation. Based upon the free space fundamental solution, we derive a representation for the solution f as a nonlinear Volterra integral equation of second kind with a weakly singular kernel. Uniqueness and reconstructibility by iteration is an immediate result of a priori assumption on f and applying the fixed point theorem. Numerical examples are presented to illustrate the validity and effectiveness of the proposed method. © 2013 Copyright Taylor and Francis Group, LLC.

  14. Optimal control problems for impulsive systems with integral boundary conditions

    Directory of Open Access Journals (Sweden)

    Allaberen Ashyralyev

    2013-03-01

    Full Text Available In this article, the optimal control problem is considered when the state of the system is described by the impulsive differential equations with integral boundary conditions. Applying the Banach contraction principle the existence and uniqueness of the solution is proved for the corresponding boundary problem by the fixed admissible control. The first and second variation of the functional is calculated. Various necessary conditions of optimality of the first and second order are obtained by the help of the variation of the controls.

  15. Dual integral equations with Fox's H-function kernel

    Directory of Open Access Journals (Sweden)

    R. N. Kalia

    1984-01-01

    Full Text Available The dual integral equations involving Bessel function kernels were first considered by Weber in 1873. The problem comprised of finding potential of an electrified disc which belongs to a general category of mixed boundary value problems. Titchmarsh gave the formal solution using Wiener-Hopf procedure. We use this direct method as improvised by Busbridge to solve a class of dual integral equations which can be reduced to other known kernels by particularizing the parameters in the Fox's H-function.

  16. Maxwell boundary conditions impose non-Lindblad master equation

    CERN Document Server

    Bamba, Motoaki

    2016-01-01

    From the Hamiltonian connecting the inside and outside of an Fabry-Perot cavity, which is derived from the Maxwell boundary conditions at a mirror of the cavity, a master equation of a non-Lindblad form is derived when the cavity embeds matters, although we can transform it to the Lindblad form by performing the rotating-wave approximation to that Hamiltonian. We calculate absorption spectra by these Lindblad and non-Lindblad master equations and also by the Maxwell boundary conditions in framework of the classical electrodynamics, which we consider the most reliable approach. We found that, compared to the Lindblad master equation, the absorption spectra by the non-Lindblad one agree better with those by the Maxwell boundary conditions. Although the discrepancy is highlighted only in the ultra-strong light-matter interaction regime with a relatively large broadening, the master equation of the non-Lindblad form is preferable rather than of the Lindblad one for pursuing the consistency with the classical elec...

  17. ON A PARABOLIC FREE BOUNDARY EQUATION MODELING PRICE FORMATION

    KAUST Repository

    MARKOWICH, P. A.

    2009-10-01

    We discuss existence and uniqueness of solutions for a one-dimensional parabolic evolution equation with a free boundary. This problem was introduced by Lasry and Lions as description of the dynamical formation of the price of a trading good. Short time existence and uniqueness is established by a contraction argument. Then we discuss the issue of global-in-time-extension of the local solution which is closely related to the regularity of the free boundary. We also present numerical results. © 2009 World Scientific Publishing Company.

  18. Counting master integrals. Integration by parts vs. functional equations

    Energy Technology Data Exchange (ETDEWEB)

    Kniehl, Bernd A.; Tarasov, Oleg V. [Hamburg Univ. (Germany). II. Inst. fuer Theoretische Physik

    2016-01-15

    We illustrate the usefulness of functional equations in establishing relationships between master integrals under the integration-by-parts reduction procedure by considering a certain two-loop propagator-type diagram as an example.

  19. Counting master integrals. Integration by parts vs. functional equations

    International Nuclear Information System (INIS)

    We illustrate the usefulness of functional equations in establishing relationships between master integrals under the integration-by-parts reduction procedure by considering a certain two-loop propagator-type diagram as an example.

  20. Counting master integrals: Integration by parts vs. functional equations

    CERN Document Server

    Kniehl, Bernd A

    2016-01-01

    We illustrate the usefulness of functional equations in establishing relationships between master integrals under the integration-by-parts reduction procedure by considering a certain two-loop propagator-type diagram as an example.

  1. Optimal Control of a Parabolic Equation with Dynamic Boundary Condition

    International Nuclear Information System (INIS)

    We investigate a control problem for the heat equation. The goal is to find an optimal heat transfer coefficient in the dynamic boundary condition such that a desired temperature distribution at the boundary is adhered. To this end we consider a function space setting in which the heat flux across the boundary is forced to be an Lp function with respect to the surface measure, which in turn implies higher regularity for the time derivative of temperature. We show that the corresponding elliptic operator generates a strongly continuous semigroup of contractions and apply the concept of maximal parabolic regularity. This allows to show the existence of an optimal control and the derivation of necessary and sufficient optimality conditions.

  2. Engineering integrable nonautonomous nonlinear Schroedinger equations

    International Nuclear Information System (INIS)

    We investigate Painleve integrability of a generalized nonautonomous one-dimensional nonlinear Schroedinger (NLS) equation with time- and space-dependent dispersion, nonlinearity, and external potentials. Through the Painleve analysis some explicit requirements on the dispersion, nonlinearity, dissipation/gain, and the external potential as well as the constraint conditions are identified. It provides an explicit way to engineer integrable nonautonomous NLS equations at least in the sense of Painleve integrability. Furthermore analytical solutions of this class of integrable nonautonomous NLS equations can be obtained explicitly from the solutions of the standard NLS equation by a general transformation. The result provides a significant way to control coherently the soliton dynamics in the corresponding nonlinear systems, as that in Bose-Einstein condensate experiments. We analyze explicitly the soliton dynamics under the nonlinearity management and the external potentials and discuss its application in the matter-wave dynamics. Some comparisons with the previous works have also been discussed.

  3. Coupling Integrable Couplings of an Equation Hierarchy

    International Nuclear Information System (INIS)

    Based on a kind of Lie algebra G proposed by Zhang, one isospectral problem is designed. Under the framework of zero curvature equation, a new kind of integrable coupling of an equation hierarchy is generated using the methods proposed by Ma and Gao. With the help of variational identity, we get the Hamiltonian structure of the hierarchy. (general)

  4. Integrability of doubly-periodic Riccati equation

    Directory of Open Access Journals (Sweden)

    Guan Ke-Ying

    1998-12-01

    Full Text Available By the structure of solvable subgroup of SL(2,ℂ (see [1], the integrability and properties of solutions of a Riccati equation with an elliptic function coefficient, which is related to a Fuchsian equation on the torus T2,is studied.

  5. THE COUPLING OF NATURAL BOUNDARY ELEMENT AND FINITE ELEMENT METHOD FOR 2D HYPERBOLIC EQUATIONS

    Institute of Scientific and Technical Information of China (English)

    De-hao Yu; Qi-kui Du

    2003-01-01

    In this paper, we investigate the coupling of natural boundary element and finite ele-ment methods of exterior initial boundary value problems for hyperbolic equations. Thegoverning equation is first discretized in time, leading to a time-step scheme, where anexterior elliptic problem has to be solved in each time step. Second, a circular artifi-in an unbounded domain is transformed into the nonlocal boundary value problem in abounded subdomain. And the natural integral equation and the Poisson integral formulaare obtained in the infinite domain Ω2 outside circle of radius R. The coupled variationalformulation is given. Only the function itself, not its normal derivative at artificial bound-and the boundary element stiffness matrix has a few different elements. Such a coupledmethod is superior to the one based on direct boundary element method. This paper dis-cusses finite element discretization for variational problem and its corresponding numericaltechnique, and the convergence for the numerical solutions. Finally, the numerical exampleis presented to illustrate feasibility and efficiency of this method.

  6. Solution of a Problem Linear Plane Elasticity with Mixed Boundary Conditions by the Method of Boundary Integrals

    Directory of Open Access Journals (Sweden)

    Nahed S. Hussein

    2014-01-01

    Full Text Available A numerical boundary integral scheme is proposed for the solution to the system of …eld equations of plane. The stresses are prescribed on one-half of the circle, while the displacements are given. The considered problem with mixed boundary conditions in the circle is replaced by two problems with homogeneous boundary conditions, one of each type, having a common solution. The equations are reduced to a system of boundary integral equations, which is then discretized in the usual way, and the problem at this stage is reduced to the solution to a rectangular linear system of algebraic equations. The unknowns in this system of equations are the boundary values of four harmonic functions which define the full elastic solution and the unknown boundary values of stresses or displacements on proper parts of the boundary. On the basis of the obtained results, it is inferred that a stress component has a singularity at each of the two separation points, thought to be of logarithmic type. The results are discussed and boundary plots are given. We have also calculated the unknown functions in the bulk directly from the given boundary conditions using the boundary collocation method. The obtained results in the bulk are discussed and three-dimensional plots are given. A tentative form for the singular solution is proposed and the corresponding singular stresses and displacements are plotted in the bulk. The form of the singular tangential stress is seen to be compatible with the boundary values obtained earlier. The efficiency of the used numerical schemes is discussed.

  7. Boundary value problemfor multidimensional fractional advection-dispersion equation

    OpenAIRE

    Khasambiev Mokhammad Vakhaevich

    2015-01-01

    In recent time there is a very great interest in the study of differential equations of fractional order, in which the unknown function is under the symbol of fractional derivative. It is due to the development of the theory of fractional integro-differential theory and application of it in different fields.The fractional integrals and derivatives of fractional integro-differential equations are widely used in modern investigations of theoretical physics, mechanics, and applied mathematics. T...

  8. High Accuracy Method for Integral Equations with Discontinuous Kernels

    OpenAIRE

    Kang, Sheon-Young; Koltracht, Israel; Rawitscher, George

    1999-01-01

    A new highly accurate numerical approximation scheme based on a Gauss type Clenshaw-Curtis Quadrature for Fredholm integral equations of the second kind, whose kernel is either discontinuous or not smooth along the main diagonal, is presented. This scheme is of spectral accuracy when the kernel is infinitely differentiable away from the main diagonal, and is also applicable when the kernel is singular along the boundary, and at isolated points on the main diagonal. The corresponding composite...

  9. On Approximate Asymptotic Solution of Integral Equations

    CERN Document Server

    Jikia, Vagner

    2013-01-01

    It is well known that multi-particle integral equations of collision theory, in general, are not compact. At the same time it has been shown that the motion of three and four particles is described with consistent integral equations. In particular, by using identical transformations of the kernel of the Lipman-Schwinger equation for certain classes of potentials Faddeev obtained Fredholm type integral equations for three-particle problems $[1]$. The motion of for bodies is described by equations of Yakubovsky and Alt-Grassberger-Sandhas-Khelashvili $[2.3]$, which are obtained as a result of two subsequent transpormations of the kernel of Lipman-Schwinger equation. in the case of $N>4$ the compactness of multi-particle equations has not been proven yet. In turn out that for sufficiently high energies the $N$-particle $\\left( {N \\ge 3} \\right)$ dynamic equations have correct asymptotic solutions satisfying unitary condition $[4]$. In present paper by using the Heitler formalism we obtain the results briefly sum...

  10. Path Integral for the Dirac Equation

    OpenAIRE

    Polonyi, Janos

    1998-01-01

    A c-number path integral representation is constructed for the solution of the Dirac equation. The integration is over the real trajectories in the continuous three-space and other two canonical pairs of compact variables controlling the spin and the chirality flips.

  11. Path Integral for Relativistic Equations of Motion

    OpenAIRE

    Gosselin, Pierre; Polonyi, Janos

    1997-01-01

    A non-Grassmanian path integral representation is given for the solution of the Klein-Gordon and the Dirac equations. The trajectories of the path integral are rendered differentiable by the relativistic corrections. The nonrelativistic limit is briefly discussed from the point of view of the renormalization group.

  12. A Collocation Method for Volterra Integral Equations

    Science.gov (United States)

    Kolk, Marek

    2010-09-01

    We propose a piecewise polynomial collocation method for solving linear Volterra integral equations of the second kind with logarithmic kernels which, in addition to a diagonal singularity, may have a singularity at the initial point of the interval of integration. An attainable order of the convergence of the method is studied. We illustrate our results with a numerical example.

  13. A Nonlinear Singularly Perturbed Problem for Reaction Diffusion Equations with Boundary Perturbation

    Institute of Scientific and Technical Information of China (English)

    Jia-qi Mo; Wan-tao Lin

    2005-01-01

    A nonlinear singularly perturbed problems for reaction diffusion equation with boundary perturbation is considered. Under suitable conditions, the asymptotic behavior of solution for the initial boundary value problems of reaction diffusion equations is studied using the theory of differential inequalities.

  14. A Numerical Method for 1-D Parabolic Equation with Nonlocal Boundary Conditions

    Directory of Open Access Journals (Sweden)

    M. Tadi

    2014-01-01

    Full Text Available This paper is concerned with a local method for the solution of one-dimensional parabolic equation with nonlocal boundary conditions. The method uses a coordinate transformation. After the coordinate transformation, it is then possible to obtain exact solutions for the resulting equations in terms of the local variables. These exact solutions are in terms of constants of integration that are unknown. By imposing the given boundary conditions and smoothness requirements for the solution, it is possible to furnish a set of linearly independent conditions that can be used to solve for the constants of integration. A number of examples are used to study the applicability of the method. In particular, three nonlinear problems are used to show the novelty of the method.

  15. Spectral integration of linear boundary value problems

    CERN Document Server

    Viswanath, Divakar

    2012-01-01

    Spectral integration is a method for solving linear boundary value problems which uses the Chebyshev series representation of functions to avoid the numerical discretization of derivatives. It is occasionally attributed to Zebib (J. of Computational Physics vol. 53 (1984), p. 443-455) and more often to Greengard (SIAM J. on Numerical Analysis vol. 28 (1991), p. 1071-1080). Its advantage is believed to be its relative immunity to errors that arise when nearby grid points are used to approximate derivatives. In this paper, we reformulate the method of spectral integration by changing it in four different ways. The changes consist of a more convenient integral formulation, a different way to treat and interpret boundary conditions, treatment of higher order problems in factored form, and the use of piecewise Chebyshev grid points. Our formulation of spectral integration is more flexible and powerful as show by its ability to solve a problem that would otherwise take 8192 grid points using only 96 grid points. So...

  16. Universal equations of unsteady two-dimensional MHD boundary layer whose temperature varies with time

    Directory of Open Access Journals (Sweden)

    Boričić Zoran

    2009-01-01

    Full Text Available This paper concerns with unsteady two-dimensional temperature laminar magnetohydrodynamic (MHD boundary layer of incompressible fluid. It is assumed that induction of outer magnetic field is function of longitudinal coordinate with force lines perpendicular to the body surface on which boundary layer forms. Outer electric filed is neglected and magnetic Reynolds number is significantly lower then one i.e. considered problem is in inductionless approximation. Characteristic properties of fluid are constant because velocity of flow is much lower than speed of light and temperature difference is small enough (under 50ºC . Introduced assumptions simplify considered problem in sake of mathematical solving, but adopted physical model is interesting from practical point of view, because its relation with large number of technically significant MHD flows. Obtained partial differential equations can be solved with modern numerical methods for every particular problem. Conclusions based on these solutions are related only with specific temperature MHD boundary layer problem. In this paper, quite different approach is used. First new variables are introduced and then sets of similarity parameters which transform equations on the form which don't contain inside and in corresponding boundary conditions characteristics of particular problems and in that sense equations are considered as universal. Obtained universal equations in appropriate approximation can be solved numerically once for all. So-called universal solutions of equations can be used to carry out general conclusions about temperature MHD boundary layer and for calculation of arbitrary particular problems. To calculate any particular problem it is necessary also to solve corresponding momentum integral equation.

  17. About potential of double layer and boundary value problems for Laplace equation

    International Nuclear Information System (INIS)

    An integral operator raisen by a kernel of the double layer's potential is investigated. The kernel is defined on S (S - two-digit variety of C2 class presented by a boundary of the finite domain in R3). The operator is considered on C(S). Following results are received: the operator's spectrum belongs to [-1,1]; it's eigenvalues and eigenfunctions may be found by Kellog's method; knowledge of the operator's spectrum is enough to construct it's resolvent. These properties permit to point out the determined interation processes, solving boundary value problems for Laplace equation. One of such processes - solving of Roben problem - is generalized on electrostatic problems. 6 refs

  18. The Navier-Stokes Equations with the Kinematic and Vorticity Boundary Conditions on Non-Flat Boundaries

    OpenAIRE

    Chen, Gui-Qiang; Osborne, Dan; Qian, Zhongmin

    2008-01-01

    We study the initial-boundary value problem of the Navier-Stokes equations for incompressible fluids in a general domain in $\\R^n$ with compact and smooth boundary, subject to the kinematic and vorticity boundary conditions on the non-flat boundary. We observe that, under the nonhomogeneous boundary conditions, the pressure $p$ can be still recovered by solving the Neumann problem for the Poisson equation. Then we establish the well-posedness of the unsteady Stokes equations and employ the so...

  19. Asymptotic behavior of solutions to wave equations with a memory condition at the boundary

    Directory of Open Access Journals (Sweden)

    Mauro de Lima Santos

    2001-11-01

    Full Text Available In this paper, we study the stability of solutions for wave equations whose boundary condition includes a integral that represents the memory effect. We show that the dissipation is strong enough to produce exponential decay of the solution, provided the relaxation function also decays exponentially. When the relaxation function decays polynomially, we show that the solution decays polynomially and with the same rate.

  20. Periodic and Boundary Value Problems for Second Order Differential Equations

    Indian Academy of Sciences (India)

    Nikolaos S Papageorgiou; Francesca Papalini

    2001-02-01

    In this paper we study second order scalar differential equations with Sturm–Liouville and periodic boundary conditions. The vector field (, , ) is Caratheodory and in some instances the continuity condition on or is replaced by a monotonicity type hypothesis. Using the method of upper and lower solutions as well as truncation and penalization techniques, we show the existence of solutions and extremal solutions in the order interval determined by the upper and lower solutions. Also we establish some properties of the solutions and of the set they form.

  1. Adaptive integral equation methods in transport theory

    International Nuclear Information System (INIS)

    In this paper, an adaptive multilevel algorithm for integral equations is described that has been developed with the Chandrasekhar H equation and its generalizations in mind. The algorithm maintains good performance when the Frechet derivative of the nonlinear map is singular at the solution, as happens in radiative transfer with conservative scattering and in critical neutron transport. Numerical examples that demonstrate the algorithm's effectiveness are presented

  2. Inequalities for differential and integral equations

    CERN Document Server

    Ames, William F

    1997-01-01

    Inequalities for Differential and Integral Equations has long been needed; it contains material which is hard to find in other books. Written by a major contributor to the field, this comprehensive resource contains many inequalities which have only recently appeared in the literature and which can be used as powerful tools in the development of applications in the theory of new classes of differential and integral equations. For researchers working in this area, it will be a valuable source of reference and inspiration. It could also be used as the text for an advanced graduate course.Key Features* Covers a variety of linear and nonlinear inequalities which find widespread applications in the theory of various classes of differential and integral equations* Contains many inequalities which have only recently appeared in literature and cannot yet be found in other books* Provides a valuable reference to engineers and graduate students

  3. Quadratic integral solutions to double Pell equations

    CERN Document Server

    Veneziano, Francesco

    2011-01-01

    We study the quadratic integral points-that is, (S-)integral points defined over any extension of degree two of the base field-on a curve defined in P_3 by a system of two Pell equations. Such points belong to three families explicitly described, or belong to a finite set whose cardinality may be explicitly bounded in terms of the base field, the equations defining the curve and the set S. We exploit the peculiar geometry of the curve to adapt the proof of a theorem of Vojta, which in this case does not apply.

  4. Algorithms For Integrating Nonlinear Differential Equations

    Science.gov (United States)

    Freed, A. D.; Walker, K. P.

    1994-01-01

    Improved algorithms developed for use in numerical integration of systems of nonhomogenous, nonlinear, first-order, ordinary differential equations. In comparison with integration algorithms, these algorithms offer greater stability and accuracy. Several asymptotically correct, thereby enabling retention of stability and accuracy when large increments of independent variable used. Accuracies attainable demonstrated by applying them to systems of nonlinear, first-order, differential equations that arise in study of viscoplastic behavior, spread of acquired immune-deficiency syndrome (AIDS) virus and predator/prey populations.

  5. Boundary integral method applied in chaotic quantum billiards

    CERN Document Server

    Li, B; Li, Baowen; Robnik, Marko

    1995-01-01

    The boundary integral method (BIM) is a formulation of Helmholtz equation in the form of an integral equation suitable for numerical discretization to solve the quantum billiard. This paper is an extensive numerical survey of BIM in a variety of quantum billiards, integrable (circle, rectangle), KAM systems (Robnik billiard) and fully chaotic (ergodic, such as stadium, Sinai billiard and cardioid billiard). On the theoretical side we point out some serious flaws in the derivation of BIM in the literature and show how the final formula (which nevertheless was correct) should be derived in a sound way and we also argue that a simple minded application of BIM in nonconvex geometries presents serious difficulties or even fails. On the numerical side we have analyzed the scaling of the averaged absolute value of the systematic error \\Delta E of the eigenenergy in units of mean level spacing with the density of discretization (b = number of numerical nodes on the boundary within one de Broglie wavelength), and we f...

  6. PREFACE: Symmetries and Integrability of Difference Equations

    Science.gov (United States)

    Doliwa, Adam; Korhonen, Risto; Lafortune, Stéphane

    2007-10-01

    The notion of integrability was first introduced in the 19th century in the context of classical mechanics with the definition of Liouville integrability for Hamiltonian flows. Since then, several notions of integrability have been introduced for partial and ordinary differential equations. Closely related to integrability theory is the symmetry analysis of nonlinear evolution equations. Symmetry analysis takes advantage of the Lie group structure of a given equation to study its properties. Together, integrability theory and symmetry analysis provide the main method by which nonlinear evolution equations can be solved explicitly. Difference equations (DE), like differential equations, are important in numerous fields of science and have a wide variety of applications in such areas as mathematical physics, computer visualization, numerical analysis, mathematical biology, economics, combinatorics, and quantum field theory. It is thus crucial to develop tools to study and solve DEs. While the theory of symmetry and integrability for differential equations is now largely well-established, this is not yet the case for discrete equations. Although over recent years there has been significant progress in the development of a complete analytic theory of difference equations, further tools are still needed to fully understand, for instance, the symmetries, asymptotics and the singularity structure of difference equations. The series of SIDE meetings on Symmetries and Integrability of Difference Equations started in 1994. Its goal is to provide a platform for an international and interdisciplinary communication for researchers working in areas associated with integrable discrete systems, such as classical and quantum physics, computer science and numerical analysis, mathematical biology and economics, discrete geometry and combinatorics, theory of special functions, etc. The previous SIDE meetings took place in Estérel near Montréal, Canada (1994), at the University of

  7. THE NAVIER-STOKES EQUATIONS WITH THE KINEMATIC AND VORTICITY BOUNDARY CONDITIONS ON NON-FLAT BOUNDARIES

    Institute of Scientific and Technical Information of China (English)

    Gui-Qiang Chen; Dan Osborne; Zhongmin Qian

    2009-01-01

    We study the initial-boundary value problem of the Navier-Stokes equations for incompressible fluids in a general domain in RN with compact and smooth boundary, subject to the kinematic and vorticity boundary conditions on the non-fiat boundary. We observe that, under the nonhomogeneons boundary conditions, the pressure p can be still recovered by solving the Neumann problem for the Poisson equation. Then we establish the well-posedness of the unsteady Stokes equations and employ the solution to reduce our initial-boundary value problem into an initial-boundary value problem with absolute boundary conditions. Based on this, we first establish the well-posedness for an appropriate local linearized problem with the absolute boundary conditions and the initial condition (without the incompressibility condition), which establishes a velocity mapping. Then we develop apriori estimates for the velocity mapping, especially involving the Sobolev norm for the time-derivative of the mapping to deal with the complicated boundary conditions, which leads to the existence of the fixed point of the mapping and the existence of solutions to our initial-boundary value problem. Finally, we establish that, when the viscosity coefficient tends zero, the strong solutions of the initial-boundary value problem in RN(n≥3) with nonhomogeneous vorticity boundary condition converge in L2 to the corresponding Euler equations satisfying the kinematic condition.

  8. Path Integral Methods for Stochastic Differential Equations

    OpenAIRE

    Chow, Carson C.; Buice, Michael A.

    2015-01-01

    Stochastic differential equations (SDEs) have multiple applications in mathematical neuroscience and are notoriously difficult. Here, we give a self-contained pedagogical review of perturbative field theoretic and path integral methods to calculate moments of the probability density function of SDEs. The methods can be extended to high dimensional systems such as networks of coupled neurons and even deterministic systems with quenched disorder.

  9. Nonparametric (smoothed) likelihood and integral equations

    CERN Document Server

    Groeneboom, Piet

    2012-01-01

    We show that there is an intimate connection between the theory of nonparametric (smoothed) maximum likelihood estimators for certain inverse problems and integral equations. This is illustrated by estimators for interval censoring and deconvolution problems. We also discuss the asymptotic efficiency of the MLE for smooth functionals in these models.

  10. Simplified differential equations approach for Master Integrals

    OpenAIRE

    Papadopoulos, Costas G.

    2014-01-01

    A simplified differential equations approach for Master Integrals is presented. It allows to express them, straightforwardly, in terms of Goncharov Polylogarithms. As a proof-of-concept of the proposed method, results at one and two loops are presented, including the massless one-loop pentagon with up to one off-shell leg at order epsilon.

  11. Bethe ansatz for the Temperley–Lieb spin chain with integrable open boundaries

    International Nuclear Information System (INIS)

    In this paper we study the spectrum of the spin-1 Temperley–Lieb spin chain with integrable open boundary conditions. We obtain the eigenvalue expressions as well as its associated Bethe ansatz equations by means of the coordinate Bethe ansatz. These equations provide the complete description of the spectrum of the model. (paper)

  12. Three New Integrable Hierarchies of Equations

    International Nuclear Information System (INIS)

    A general Lie algebra Vs and the corresponding loop algebra V-tildes are constructed, from which the linear isospectral Lax pairs are established, whose compatibility presents the zero curvature equation. As its application, a new Lax integrable hierarchy containing two parameters is worked out. It is not Liouville-integrable, however, its two reduced systems are Liouville-integrable, whose Hamiltonian structures are derived by making use of the quadratic-form identity and the γ formula (i.e. the computational formula on the constant γ appeared in the trace identity and the quadratic-form identity).

  13. Three New Integrable Hierarchies of Equations

    Institute of Scientific and Technical Information of China (English)

    2007-01-01

    A general Lie algebra Vs and the corresponding loop algebra ~Vs are constructed, from which the linear isospectral Lax pairs are established, whose compatibility presents the zero curvature equation. As its application, a new Lax integrable hierarchy containing two parameters is worked out. It is not Liouville-integrable, however, its two reduced systems are Liouville-integrable, whose Hamiltonian structures are derived by making use of the quadratic-form identity and the γ formula (i.e. the computational formula on the constant γ appeared in the trace identity and the quadratic-form identity).

  14. Lectures on differential equations for Feynman integrals

    International Nuclear Information System (INIS)

    Over the last year significant progress was made in the understanding of the computation of Feynman integrals using differential equations (DE). These lectures give a review of these developments, while not assuming any prior knowledge of the subject. After an introduction to DE for Feynman integrals, we point out how they can be simplified using algorithms available in the mathematical literature. We discuss how this is related to a recent conjecture for a canonical form of the equations. We also discuss a complementary approach that is based on properties of the space–time loop integrands, and explain how the ideas of leading singularities and d-log representations can be used to find an optimal basis for the DE. Finally, as an application of these ideas we show how single-scale integrals can be bootstrapped using the Drinfeld associator of a DE. (topical review)

  15. Non-local problems with integral gluing condition for loaded mixed type equations involving the Caputo fractional derivative

    Directory of Open Access Journals (Sweden)

    Obidjon Kh. Abdullaev

    2016-06-01

    Full Text Available In this work, we study the existence and uniqueness of solutions to non-local boundary value problems with integral gluing condition. Mixed type equations (parabolic-hyperbolic involving the Caputo fractional derivative have loaded parts in Riemann-Liouville integrals. Thus we use the method of integral energy to prove uniqueness, and the method of integral equations to prove existence.

  16. A direct solver with O(N) complexity for integral equations on one-dimensional domains

    OpenAIRE

    Gillman, Adrianna; Young, Patrick; Martinsson, Per-Gunnar

    2011-01-01

    An algorithm for the direct inversion of the linear systems arising from Nystrom discretization of integral equations on one-dimensional domains is described. The method typically has O(N) complexity when applied to boundary integral equations (BIEs) in the plane with non-oscillatory kernels such as those associated with the Laplace and Stokes' equations. The scaling coefficient suppressed by the "big-O" notation depends logarithmically on the requested accuracy. The method can also be applie...

  17. Nonlocal singular problem with integral condition for a second-order parabolic equation

    Directory of Open Access Journals (Sweden)

    Ahmed Lakhdar Marhoune

    2015-03-01

    Full Text Available We prove the existence and uniqueness of a strong solution for a parabolic singular equation in which we combine Dirichlet with integral boundary conditions given only on parts of the boundary. The proof uses a priori estimate and the density of the range of the operator generated by the problem considered.

  18. Reducing differential equations for multiloop master integrals

    OpenAIRE

    Lee, Roman

    2014-01-01

    We present an algorithm of the reduction of the differential equations for master integrals the Fuchsian form with the right-hand side matrix linearly depending on dimensional regularization parameter $\\epsilon$. We consider linear transformations of the functions column which are rational in the variable and in $\\epsilon$. Apart from some degenerate cases described below, the algorithm allows one to obtain the required transformation or to ascertain irreducibility to the form required. Degen...

  19. Boundary conditions for the Pauli equation: application to photodetachment of Cs-.

    Science.gov (United States)

    Bahrim, C; Fabrikant, I I; Thumm, U

    2001-09-17

    We formulate the boundary conditions near the atomic nucleus for solving the Pauli equation, based on the analytic solution of the Dirac equation for a Coulomb potential. We then integrate the Pauli equation using an effective potential that is adjusted to reproduce Dirac R-matrix scattering phase shifts, and find the (3)P(o)(1) resonance contribution to the photodetachment cross section of Cs-. Our photodetachment cross sections agree with recent experiments by Scheer et al. [Phys. Rev. Lett. 80, 684 (1998)] after tuning the resonance position by 2.4 meV. We also provide angle-differential photodetachment cross sections and the corresponding asymmetry parameter beta near the Cs(6s) threshold. PMID:11580504

  20. Recursion operators, conservation laws, and integrability conditions for difference equations

    Science.gov (United States)

    Mikhailov, A. V.; Wang, Jing Ping; Xenitidis, P.

    2011-04-01

    We attempt to propose an algebraic approach to the theory of integrable difference equations. We define the concept of a recursion operator for difference equations and show that it generates an infinite sequence of symmetries and canonical conservation laws for a difference equation. As in the case of partial differential equations, these canonical densities can serve as integrability conditions for difference equations. We obtain the recursion operators for the Viallet equation and all the Adler-Bobenko-Suris equations.

  1. SOLVING INTEGRAL EQUATIONS WITH LOGARITHMIC KERNEL BY USING PERIODIC QUASI-WAVELET

    Institute of Scientific and Technical Information of China (English)

    Han-lin Chen; Si-long Peng

    2000-01-01

    In solving integral equations with logarithmic kernel which arises from the boundary integral equation reformulation of some boundary value problems for the two dimensional Helmholtz equation, we combine the Galerkin method with Beylkin's ([2]) approach, series of dense and nonsymmetric matrices may appear if we use traditional method. By appealing the so-called periodic quasi-wavelet (PQW in abbr.) ([5]), some of these matrices become diagonal, therefore we can find a algorithm with only O(K(m)2) arithmetic operations where m is the highest level. The Galerkin approximation has a polynomial rate of convergence.

  2. A Class of Nonlinear Singularly Perturbed Problems for Reaction Diffusion Equations with Boundary Perturbation

    Institute of Scientific and Technical Information of China (English)

    Jia-qi Mo; Wan-tao Lin

    2006-01-01

    A class of nonlinear singularly perturbed problems for reaction diffusion equations with boundary perturbation are considered. Under suitable conditions, the asymptotic behavior of solution for the initial boundary value problems is studied using the theory of differential inequalities.

  3. THE NONLINEAR SINGULARLY PERTURBED PROBLEMS FOR REACTION DIFFUSION EQUATIONS WITH BOUNDARY PERTURBATION

    Institute of Scientific and Technical Information of China (English)

    OuyangCheng; MoJiaqi

    2005-01-01

    The nonlinear singularly perturbed problems for reaction diffusion equations with boundary perturbation are considered. Under suitable conditions, using the theory of differential inequalities the asymptotic behavior of solution for the initial boundary value problems is studied.

  4. A Priori Estimates for Solutions of Boundary Value Problems for Fractional-Order Equations

    CERN Document Server

    Alikhanov, A A

    2011-01-01

    We consider boundary value problems of the first and third kind for the diffusionwave equation. By using the method of energy inequalities, we find a priori estimates for the solutions of these boundary value problems.

  5. An initial-boundary value problem for three-dimensional Zakharov-Kuznetsov equation

    Science.gov (United States)

    Faminskii, Andrei V.

    2016-02-01

    An initial-boundary value problem with homogeneous Dirichlet boundary conditions for three-dimensional Zakharov-Kuznetsov equation is considered. Results on global existence, uniqueness and large-time decay of weak solutions in certain weighted spaces are established.

  6. LIMIT BEHAVIOUR OF SOLUTIONS TO EQUIVALUED SURFACE BOUNDARY VALUE PROBLEM FOR PARABOLIC EQUATIONS

    Institute of Scientific and Technical Information of China (English)

    LI Fengquan

    2002-01-01

    In this paper, we discuss the limit behaviour of solutions to equivalued surface boundary value problem for parabolic equations when the equivalued surface boundary shrinks to a point and the space dimension of the domain is two or more.

  7. SINGULARLY PERTURBED BOUNDARY VALUE PROBLEMS FOR ELLIPTIC EQUATION WITH A CURVE OF TURNING POINT

    Institute of Scientific and Technical Information of China (English)

    MoJiaqi

    2002-01-01

    The singularly perturbed elliptic equation boundary value problem with a curve of turning point is considered. Using the method of multiple scales and the comparison theorem,the asymptotic behavior of solution for the boundary value problem is studied.

  8. Integrated watershed planning across jurisdictional boundaries

    Science.gov (United States)

    Watts, A. W.; Roseen, R.; Stacey, P.; Bourdeau, R.

    2014-12-01

    We will present the foundation for an Coastal Watershed Integrated Plan for three communities in southern New Hampshire. Small communities are often challenged by complex regulatory requirements and limited resources, but are wary of perceived risks in engaging in collaborative projects with other communities. Potential concerns include loss of control, lack of resources to engage in collaboration, technical complexity, and unclear benefits. This project explores a multi-town subwatershed application of integrated planning across jurisdictional boundaries that addresses some of today's highest priority water quality issues: wastewater treatment plant upgrades for nutrient removal; green infrastructure stormwater management for developing and re-developing areas; and regional monitoring of ecosystem indicators in support of adaptive management to achieve nutrient reduction and other water quality goals in local and downstream waters. The project outcome is a collaboratively-developed inter-municipal integrated plan, and a monitoring framework to support cross jurisdictional planning and assess attainment of water quality management goals. This research project has several primary components: 1) assessment of initial conditions, including both the pollutant load inputs and the political, economic and regulatory status within each community, 2) a pollutant load model for point and non-point sources, 3) multi-criteria evaluation of load reduction alternatives 4) a watershed management plan optimized for each community, and for Subwatersheds combining multiple communities. The final plan will quantify the financial and other benefits/drawbacks to each community for both inter municipal and individual pollution control approaches. We will discuss both the technical and collaborative aspects of the work, with lessons learned regarding science to action, incorporation of social, economic and water quality assessment parameters, and stakeholder/researcher interaction.

  9. A Hierarchy of Lax Integrable Lattice Equations,Liouville Integrability and a New Integrable Symplectic Map

    Institute of Scientific and Technical Information of China (English)

    XU Xi-Xiang; ZHANG Yu-Feng

    2004-01-01

    A discrete matrix spectral problem and the associated hierarchy of Lax integrable lattice equations are presented, and it is shown that the resulting Lax integrable lattice equations are all Liouville integrable discrete Hamiltonian systems. A new integrable symplectic map is given by binary Bargmann .constraint of the resulting hierarchy.Finally, an infinite set of conservation laws is given for the resulting hierarchy.

  10. Monotone iterative technique for fractional differential equations with periodic boundary conditions

    Directory of Open Access Journals (Sweden)

    J. D. Ramírez

    2009-01-01

    Full Text Available In this paper we develop Monotone Method using upper and lower solutions for fractional differential equations with periodic boundary conditions. Initially we develop a comparison result and prove that the solution of the linear fractional differential equation with periodic boundary condition exists and is unique. Using this we develop iterates which converge uniformly monotonically to minimal and maximal solutions of the nonlinear fractional differential equations with periodic boundary conditions in the weighted norm.

  11. Monotone iterative technique for fractional differential equations with periodic boundary conditions

    OpenAIRE

    J. D. Ramírez; A. S. Vatsala

    2009-01-01

    In this paper we develop Monotone Method using upper and lower solutions for fractional differential equations with periodic boundary conditions. Initially we develop a comparison result and prove that the solution of the linear fractional differential equation with periodic boundary condition exists and is unique. Using this we develop iterates which converge uniformly monotonically to minimal and maximal solutions of the nonlinear fractional differential equations with periodic boundary con...

  12. Difference equations for correlation functions of $A^{(1)}_{n-1}$-face model with boundary reflection

    CERN Document Server

    Quano, Y H

    2001-01-01

    The $A^{(1)}_{n-1}$-face model with boundary reflection is considered on the basis of the boundary CTM bootstrap. We construct the fused boundary Boltzmann weights to determine the normalization factor. We derive difference equations of the quantum Knizhnik-Zamolodchikov type for correlation functions of the present model. The simplest difference equations are solved the in the case of the free boundary condition.

  13. Boundary value problems of the linearized non-homogeneous Navier-Stokes equations for the axisymmetrical slow motion

    International Nuclear Information System (INIS)

    The boundary value problem of axisymmetrical slow motion of viscous incompressible fluid in the upper half-space or in the whole space surrounding a fixed infinite circular cylinder is considered. The solution of the linearized non-homogeneous Navier-Stokes equations is obtained in each case in quadratures by using Abel integral equations, which transform the axisymmetrical problem to solvable two dimensional plane problem. (author)

  14. SECOND-ORDER ANALYSIS OF A BOUNDARY CONTROL PROBLEM FOR THE VISCOUS CAHN–HILLIARD EQUATION WITH DYNAMIC BOUNDARY CONDITION

    Directory of Open Access Journals (Sweden)

    Pierluigi Colli

    2015-07-01

    Full Text Available In this paper we establish second-order sufficient optimality conditions for a boundary control problem that has been introduced and studied by three of the authors in the preprint arXiv:1407.3916. This control problem regards the viscous Cahn–Hilliard equation with possibly singular potentials and dynamic boundary conditions.

  15. ON SOLVABILITY OF A BOUNDARY VALUE PROBLEM FOR A NONHOMOGENEOUS BIHARMONIC EQUATION WITH A BOUNDARY OPERATOR OF A FRACTIONAL ORDER

    Institute of Scientific and Technical Information of China (English)

    A.S. BERDYSHEV; A. CABADA; B.Kh. TURMETOV

    2014-01-01

    This paper is concerned with the solvability of a boundary value problem for a nonhomogeneous biharmonic equation. The boundary data is determined by a differential operator of fractional order in the Riemann-Liouville sense. The considered problem is a generalization of the known Dirichlet and Neumann problems.

  16. Bounded solutions for fuzzy differential and integral equations

    Energy Technology Data Exchange (ETDEWEB)

    Nieto, Juan J. [Departamento de Analisis Matematico Facultad de Matematicas Universidad de Santiago de Compostela, 15782 (Spain)] e-mail: amnieto@usc.es; Rodriguez-Lopez, Rosana [Departamento de Analisis Matematico Facultad de Matematicas Universidad de Santiago de Compostela, 15782 (Spain)] e-mail: amrosana@usc.es

    2006-03-01

    We find sufficient conditions for the boundness of every solution of first-order fuzzy differential equations as well as certain fuzzy integral equations. Our results are based on several theorems concerning crisp differential and integral inequalities.

  17. A Study of the Navier-Stokes Equations with the Kinematic and Navier Boundary Conditions

    OpenAIRE

    Chen, Gui-Qiang; Qian, Zhongmin

    2008-01-01

    We study the initial-boundary value problem of the Navier-Stokes equations for incompressible fluids in a domain in $\\R^3$ with compact and smooth boundary, subject to the kinematic and Navier boundary conditions. We first reformulate the Navier boundary condition in terms of the vorticity, which is motivated by the Hodge theory on manifolds with boundary from the viewpoint of differential geometry, and establish basic elliptic estimates for vector fields subject to the kinematic and Navier b...

  18. Path Integral and Solutions of the Constraint Equations: The Case of Reducible Gauge Theories

    OpenAIRE

    Ferraro, R.; Henneaux, M.; Puchin, M.

    1994-01-01

    It is shown that the BRST path integral for reducible gauge theories, with appropriate boundary conditions on the ghosts, is a solution of the constraint equations. This is done by relating the BRST path integral to the kernel of the evolution operator projected on the physical subspace.

  19. Explicit expressions for three-dimensional boundary integrals in linear elasticity

    Energy Technology Data Exchange (ETDEWEB)

    Nintcheu Fata, Sylvain [ORNL

    2011-01-01

    On employing isoparametric, piecewise linear shape functions over a flat triangle, exact formulae are derived for all surface potentials involved in the numerical treatment of three-dimensional singular and hyper-singular boundary integral equations in linear elasticity. These formulae are valid for an arbitrary source point in space and are represented as analytical expressions along the edges of the integration triangle. They can be employed to solve integral equations defined on triangulated surfaces via a collocation method or may be utilized as analytical expressions for the inner integrals in a Galerkin technique. A numerical example involving a unit triangle and a source point located at various distances above it, as well as sample problems solved by a collocation boundary element method for the Lame equation are included to validate the proposed formulae.

  20. Funnel Equations and Multivalued Integration Problems for Control Synthesis

    OpenAIRE

    Kurzhanski, A.B.; Nikonov, O.I.

    1989-01-01

    This paper indicates a sequence of evolutionary "funnel equations" with set-valued solutions which are crucial for the construction of respective feed-back control strategies along the schemes introduced by N.N. Krasovskii. The integration of these funnel equations leads to a sequence of multivalued integrals that generalize some of those introduced earlier (Auman's integral, the Convolution integral, Pontryagin's Alternated integral).

  1. Evolution equations: Frobenius integrability, conservation laws and travelling waves

    International Nuclear Information System (INIS)

    We give new results concerning the Frobenius integrability and solution of evolution equations admitting travelling wave solutions. In particular, we give a powerful result which explains the extraordinary integrability of some of these equations. We also discuss ‘local’ conservations laws for evolution equations in general and demonstrate all the results for the Korteweg–de Vries equation. (paper)

  2. An accurate numerical method for systems of differentio-integral equations associated with multiphase flow

    CERN Document Server

    Shu, Jian-Jun

    2014-01-01

    A very simple and accurate numerical method which is applicable to systems of differentio-integral equations with quite general boundary conditions has been devised. Although the basic idea of this method stems from the Keller Box method, it solves the problem of systems of differential equations involving integral operators not previously considered by the Keller Box method. Two main preparatory stages are required: (i) a merging procedure for differential equations and conditions without integral operators and; (ii) a reduction procedure for differential equations and conditions with integral operators. The differencing processes are effectively simplified by means of the unit-step function. The nonlinear difference equations are solved by Newton method using an efficient block arrow-like matrix factorization technique. As an example of the application of this method, the systems of equations for combined gravity body force and forced convection in laminar film condensation can be solved for prescribed valu...

  3. A boundary condition to the Khokhlov-Zabolotskaya equation for modeling strongly focused nonlinear ultrasound fields

    Science.gov (United States)

    Rosnitskiy, P.; Yuldashev, P.; Khokhlova, V.

    2015-10-01

    An equivalent source model was proposed as a boundary condition to the nonlinear parabolic Khokhlov-Zabolotskaya (KZ) equation to simulate high intensity focused ultrasound (HIFU) fields generated by medical ultrasound transducers with the shape of a spherical shell. The boundary condition was set in the initial plane; the aperture, the focal distance, and the initial pressure of the source were chosen based on the best match of the axial pressure amplitude and phase distributions in the Rayleigh integral analytic solution for a spherical transducer and the linear parabolic approximation solution for the equivalent source. Analytic expressions for the equivalent source parameters were derived. It was shown that the proposed approach allowed us to transfer the boundary condition from the spherical surface to the plane and to achieve a very good match between the linear field solutions of the parabolic and full diffraction models even for highly focused sources with F-number less than unity. The proposed method can be further used to expand the capabilities of the KZ nonlinear parabolic equation for efficient modeling of HIFU fields generated by strongly focused sources.

  4. Adomian Method for Solving Fuzzy Fredholm-Volterra Integral Equations

    Directory of Open Access Journals (Sweden)

    M. Barkhordari Ahmadi

    2013-09-01

    Full Text Available In this paper, Adomian method has been applied to approximate the solution of fuzzy volterra-fredholm integral equation. That, by using parametric form of fuzzy numbers, a fuzzy volterra-fredholm integral equation has been converted to a system of volterra-fredholm integral equation in crisp case. Finally, the method is explained with illustrative examples.

  5. BOUNDARY INTEGRAL FORMULA OF ELASTIC PROBLEMS IN CIRCLE PLANE

    Institute of Scientific and Technical Information of China (English)

    DONG Zheng-zhu; LI Shun-cai; YU De-hao

    2005-01-01

    By bianalytic functions, the boundary integral formula of the stress function for the elastic problem in a circle plane is developed. But this integral formula includes a strongly singular integral and can not be directly calculated. After the stress function is expounded to Fourier series, making use of some formulas in generalized functions to the convolutions, the boundary integral formula which does not include strongly singular integral is derived further. Then the stress function can be got simply by the integration of the values of the stress function and its derivative on the boundary. Some examples are given. It shows that the boundary integral formula of the stress function for the elastic problem is convenient.

  6. One-way spatial integration of hyperbolic equations

    Science.gov (United States)

    Towne, Aaron; Colonius, Tim

    2015-11-01

    In this paper, we develop and demonstrate a method for constructing well-posed one-way approximations of linear hyperbolic systems. We use a semi-discrete approach that allows the method to be applied to a wider class of problems than existing methods based on analytical factorization of idealized dispersion relations. After establishing the existence of an exact one-way equation for systems whose coefficients do not vary along the axis of integration, efficient approximations of the one-way operator are constructed by generalizing techniques previously used to create nonreflecting boundary conditions. When physically justified, the method can be applied to systems with slowly varying coefficients in the direction of integration. To demonstrate the accuracy and computational efficiency of the approach, the method is applied to model problems in acoustics and fluid dynamics via the linearized Euler equations; in particular we consider the scattering of sound waves from a vortex and the evolution of hydrodynamic wavepackets in a spatially evolving jet. The latter problem shows the potential of the method to offer a systematic, convergent alternative to ad hoc regularizations such as the parabolized stability equations.

  7. Well-posed initial-boundary value problems for the Zakharov-Kuznetsov equation

    Directory of Open Access Journals (Sweden)

    Andrei V. Faminskii

    2008-09-01

    Full Text Available This paper deals with non-homogeneous initial-boundary value problems for the Zakharov-Kuznetsov equation, which is one of the variants of multidimensional generalizations of the Korteweg-de Vries equation. Results on local and global well-posedness are established in a scale of Sobolev-type spaces under natural assumptions on initial and boundary data.

  8. REGULARITY THEORY FOR SYSTEMS OF PARTIAL DIFFERENTIAL EQUATIONS WITH NEUMANN BOUNDARY CONDITIONS

    Institute of Scientific and Technical Information of China (English)

    2002-01-01

    The objective of this paper is to consider the theory of regularity of systems of partial differential equations with Neumann boundary conditions. It complements previous works of the authors for the Dirichlet case. This type of problem is motivated by stochastic differential games. The Neumann case corresponds to stochastic differential equations with reflection on boundary of the domain.

  9. A Cartesian Embedded Boundary Method for the Compressible Navier-Stokes Equations

    Energy Technology Data Exchange (ETDEWEB)

    Kupiainen, M; Sjogreen, B

    2008-03-21

    We here generalize the embedded boundary method that was developed for boundary discretizations of the wave equation in second order formulation in [6] and for the Euler equations of compressible fluid flow in [11], to the compressible Navier-Stokes equations. We describe the method and we implement it on a parallel computer. The implementation is tested for accuracy and correctness. The ability of the embedded boundary technique to resolve boundary layers is investigated by computing skin-friction profiles along the surfaces of the embedded objects. The accuracy is assessed by comparing the computed skin-friction profiles with those obtained by a body fitted discretization.

  10. Material derivatives of boundary integral operators in electromagnetism and application to inverse scattering problems

    Science.gov (United States)

    Ivanyshyn Yaman, Olha; Le Louër, Frédérique

    2016-09-01

    This paper deals with the material derivative analysis of the boundary integral operators arising from the scattering theory of time-harmonic electromagnetic waves and its application to inverse problems. We present new results using the Piola transform of the boundary parametrisation to transport the integral operators on a fixed reference boundary. The transported integral operators are infinitely differentiable with respect to the parametrisations and simplified expressions of the material derivatives are obtained. Using these results, we extend a nonlinear integral equations approach developed for solving acoustic inverse obstacle scattering problems to electromagnetism. The inverse problem is formulated as a pair of nonlinear and ill-posed integral equations for the unknown boundary representing the boundary condition and the measurements, for which the iteratively regularized Gauss-Newton method can be applied. The algorithm has the interesting feature that it avoids the numerous numerical solution of boundary value problems at each iteration step. Numerical experiments are presented in the special case of star-shaped obstacles.

  11. Partial differential equations IX elliptic boundary value problems

    CERN Document Server

    Egorov, Yu; Shubin, M

    1997-01-01

    This EMS volume gives an overview of the modern theory of elliptic boundary value problems. The contribution by M.S. Agranovich is devoted to differential elliptic boundary problems, mainly in smooth bounded domains, and their spectral properties. This article continues his contribution to EMS 63. The contribution by A. Brenner and E. Shargorodsky concerns the theory of boundary value problems for elliptic pseudodifferential operators. Problems both with and without the transmission property, as well as parameter-dependent problems are considered. The article by B. Plamenevskij deals with general differential elliptic boundary value problems in domains with singularities.

  12. Local existence of solution to free boundary value problem for compressible Navier-Stokes equations

    OpenAIRE

    Liu, Jian

    2015-01-01

    This paper is concerned with the free boundary value problem for multi-dimensional Navier-Stokes equations with density-dependent viscosity where the flow density vanishes continuously across the free boundary. A local (in time) existence of weak solution is established, in particular, the density is positive and the solution is regular away from the free boundary.

  13. Boundary conditions for the Einstein-Christoffel formulation of Einstein's equations

    OpenAIRE

    Arnold, Douglas N.; Nicolae Tarfulea

    2007-01-01

    Specifying boundary conditions continues to be a challenge in numerical relativity in order to obtain a long time convergent numerical simulation of Einstein's equations in domains with artificial boundaries. In this paper, we address this problem for the Einstein--Christoffel (EC) symmetric hyperbolic formulation of Einstein's equations linearized around flat spacetime. First, we prescribe simple boundary conditions that make the problem well posed and preserve the constraints. Next, we indi...

  14. A Kind of Discrete Non-Reflecting Boundary Conditions for Varieties of Wave Equations

    Institute of Scientific and Technical Information of China (English)

    Xiu-min Shao; Zhi-ling Lan

    2002-01-01

    In this paper, a new kind of discrete non-reflecting boundary conditions is developed. It can be used for a variety of wave equations such as the acoustic wave equation, the isotropic and anisotropic elastic wave equations and the equations for wave propagation in multi-phase media and so on. In this kind of boundary conditions, the composition of all artificial reflected waves, but not the individual reflected ones, is considered and eliminated. Thus, it has a uniform formula for different wave equations. The velocity CA of the composed reflected wave is determined in the way to make the reflection coefficients minimal, the value of which depends on equations. In this paper, the construction of the boundary conditions is illustrated and CA is found, numerical results are presented to illustrate the effectiveness of the boundary conditions.

  15. Boundary singularities and characteristics of Hamilton–Jacobi equation

    NARCIS (Netherlands)

    Melikyan, A.; Olsder, G.J.

    2010-01-01

    Boundary-value problems for first order PDEs are locally considered, when the classical sufficient condition for the solution existence does not hold, but a solution still exists, possibly defined on one or both sides of the boundary surface. We note three situations when such a surface (locally) ar

  16. The Neumann problem in toroidal domains solved by an integral equation technique for the surface current

    International Nuclear Information System (INIS)

    The 3D MHD free boundary equilibrium NEMEC code combines the fixed boundary VMEC code with the vacuum NESTOR code. In NESTOR a Neumann problem is solved for a single valued potential Φ by an integral equation technique satisfying the condition rvec BV· rvec n≡(rvec B0+∇Φ)· rvec n=0 on the plasma-vacuum boundary (rvec B0=field produced by external currents and the total toroidal plasma current It, rvec n=exterior normal to the boundary). The only quantity required for the free boundary equilibrium computation is the value of BV2/2 at the boundary. Here, an equivalent method for computing the vacuum field rvec BV is considered which is particularly suited for the case It≠0 and makes use of the fact that this field can be generated by the external current jext and by a surface current rvec j on the plasma vacuum boundary. The difficulty inherent to the Green's function method is the numerical integration of the singular Green's function. This is solved by regularizing the kernel: an analytically integrable function with the same singularity is subtracted and its integral is added to the integral equation

  17. Numerical strategies to solve coupled, singular, nonlinear integral equations

    International Nuclear Information System (INIS)

    Full text: In QFT, the Dyson-Schwinger equations constitute a set of nonlinear integral equations, whose solutions are the Green's functions of the theory. Implementing consistent pictures of confinement, these integral equations become highly infrared singular. The latter property especially slows down convergence in any iterative algorithm. Different numerical strategies to solve such integral equations are compared and benchmarked results for the efficiency of different algorithms are presented. (author)

  18. Airfoil design using a coupled euler and integral boundary layer method with adjoint based sensitivities

    Science.gov (United States)

    Edwards, S.; Reuther, J.; Chattot, J. J.

    The objective of this paper is to present a control theory approach for the design of airfoils in the presence of viscous compressible flows. A coupled system of the integral boundary layer and the Euler equations is solved to provide rapid flow simulations. An adjoint approach consistent with the complete coupled state equations is employed to obtain the sensitivities needed to drive a numerical optimization algorithm. Design to a target pressure distribution is demonstrated on an RAE 2822 airfoil at transonic speeds.

  19. Solution of Moving Boundary Space-Time Fractional Burger’s Equation

    Directory of Open Access Journals (Sweden)

    E. A.-B. Abdel-Salam

    2014-01-01

    Full Text Available The fractional Riccati expansion method is used to solve fractional differential equations with variable coefficients. To illustrate the effectiveness of the method, the moving boundary space-time fractional Burger’s equation is studied. The obtained solutions include generalized trigonometric and hyperbolic function solutions. Among these solutions, some are found for the first time. The linear and periodic moving boundaries for the kink solution of the Burger’s equation are presented graphically and discussed.

  20. Porosity of Free Boundaries in the Obstacle Problem for Quasilinear Elliptic Equations

    Indian Academy of Sciences (India)

    Jun Zheng; Zhihua Zhang; Peihao Zhao

    2013-08-01

    In this paper, we establish growth rate of solutions near free boundaries in the identical zero obstacle problem for quasilinear elliptic equations. As a result, we obtain porosity of free boundaries, which is naturally an extension of the previous works by Karp et al. (J. Diff. Equ. 164 (2000) 110–117) for -Laplacian equations, and by Zheng and Zhang (J. Shaanxi Normal Univ. 40(2) (2012) 11–13, 18) for -Laplacian type equations.

  1. Hyers-Ulam stability for second-order linear differential equations with boundary conditions

    Directory of Open Access Journals (Sweden)

    Pasc Gavruta

    2011-06-01

    Full Text Available We prove the Hyers-Ulam stability of linear differential equations of second-order with boundary conditions or with initial conditions. That is, if y is an approximate solution of the differential equation $y''+ eta (x y = 0$ with $y(a = y(b =0$, then there exists an exact solution of the differential equation, near y.

  2. KdV equation under periodic boundary conditions and its perturbations

    OpenAIRE

    Huang, Guan; Kuksin, Sergei

    2013-01-01

    In this paper we discuss properties of the KdV equation under periodic boundary conditions, especially those which are important to study perturbations of the equation. Next we review what is known now about long-time behaviour of solutions for perturbed KdV equations.

  3. The KdV equation under periodic boundary conditions and its perturbations

    Science.gov (United States)

    Guan, Huang; Kuksin, Sergei

    2014-09-01

    In this paper we discuss properties of the Korteweg-de Vries (KdV) equation under periodic boundary conditions, especially those which are important for studying perturbations of the equation. We then review what is known about the long-time behaviour of solutions for perturbed KdV equations.

  4. Boundary Layer Equations and Lie Group Analysis of a Sisko Fluid

    Directory of Open Access Journals (Sweden)

    Gözde Sarı

    2012-01-01

    Full Text Available Boundary layer equations are derived for the Sisko fluid. Using Lie group theory, a symmetry analysis of the equations is performed. A partial differential system is transferred to an ordinary differential system via symmetries. Resulting equations are numerically solved. Effects of non-Newtonian parameters on the solutions are discussed.

  5. Dimension reduction for periodic boundary value problems of functional differential equations

    CERN Document Server

    Sieber, Jan

    2010-01-01

    Periodic boundary-value problems for functional differential equations can be reduced to finite-dimensional algebraic systems of equations. The smoothness assumptions on the right-hand side follow those of the review by Hartung et al. (2006) and are set up such that the result can be applied to differential equations with state-dependent delays.

  6. The adjoint equation method for constructing first integrals of difference equations

    International Nuclear Information System (INIS)

    A new method for finding first integrals of discrete equations is presented. It can be used for discrete equations which do not possess a variational (Lagrangian or Hamiltonian) formulation. The method is based on a newly established identity which links symmetries of the underlying discrete equations, solutions of the discrete adjoint equations and first integrals. The method is applied to an invariant mapping and to discretizations of second order and third order ordinary differential equations. In examples the set of independent first integrals makes it possible to find the general solution of the discrete equations. (paper)

  7. Integrability of two coupled Kadomtsev–Petviashvili equations

    Indian Academy of Sciences (India)

    Abdul-Majid Wazwaz

    2011-08-01

    The integrability of two coupled KP equations is studied. The simplified Hereman form of Hirota’s bilinear method is used to examine the integrability of each coupled equation. Multiplesoliton solutions and multiple singular soliton solutions are formally derived for each coupled KdV equation.

  8. BOUNDARY VALUE PROBLEM TO DYNAMIC EQUATION ON TIME SCALE

    Institute of Scientific and Technical Information of China (English)

    2011-01-01

    In this paper we consider a nonlinear first-order boundary value problem on a time scale. The existence results of three positive solutions are obtained using fixed point theorems. Finally,examples are presented to illustrate the main results.

  9. PREFACE: Symmetries and integrability of difference equations Symmetries and integrability of difference equations

    Science.gov (United States)

    Levi, Decio; Olver, Peter; Thomova, Zora; Winternitz, Pavel

    2009-11-01

    The concept of integrability was introduced in classical mechanics in the 19th century for finite dimensional continuous Hamiltonian systems. It was extended to certain classes of nonlinear differential equations in the second half of the 20th century with the discovery of the inverse scattering transform and the birth of soliton theory. Also at the end of the 19th century Lie group theory was invented as a powerful tool for obtaining exact analytical solutions of large classes of differential equations. Together, Lie group theory and integrability theory in its most general sense provide the main tools for solving nonlinear differential equations. Like differential equations, difference equations play an important role in physics and other sciences. They occur very naturally in the description of phenomena that are genuinely discrete. Indeed, they may actually be more fundamental than differential equations if space-time is actually discrete at very short distances. On the other hand, even when treating continuous phenomena described by differential equations it is very often necessary to resort to numerical methods. This involves a discretization of the differential equation, i.e. a replacement of the differential equation by a difference one. Given the well developed and understood techniques of symmetry and integrability for differential equations a natural question to ask is whether it is possible to develop similar techniques for difference equations. The aim is, on one hand, to obtain powerful methods for solving `integrable' difference equations and to establish practical integrability criteria, telling us when the methods are applicable. On the other hand, Lie group methods can be adapted to solve difference equations analytically. Finally, integrability and symmetry methods can be combined with numerical methods to obtain improved numerical solutions of differential equations. The origin of the SIDE meetings goes back to the early 1990s and the first

  10. Toward a theory of integrable hyperbolic equations of third order

    International Nuclear Information System (INIS)

    Examples are considered of integrable third order hyperbolic equations with two independent variables. In particular, an equation is found which admits as evolutionary symmetries the Krichever–Novikov equation and the modified Landau–Lifshitz system. The problem of the choice of dynamical variables for the hyperbolic equations is discussed. (paper)

  11. Towards the theory of integrable hyperbolic equations of third order

    OpenAIRE

    Adler, V. E.; Shabat, A. B.

    2012-01-01

    The examples are considered of integrable hyperbolic equations of third order with two independent variables. In particular, an equation is found which admits as evolutionary symmetries the Krichever--Novikov equation and the modified Landau--Lifshitz system. The problem of choice of dynamical variables for the hyperbolic equations is discussed.

  12. Difference equations and cluster algebras I: Poisson bracket for integrable difference equations

    CERN Document Server

    Inoue, Rei

    2010-01-01

    We introduce the cluster algebraic formulation of the integrable difference equations, the discrete Lotka-Volterra equation and the discrete Liouville equation, from the view point of the general T-system and Y-system. We also study the Poisson structure for the cluster algebra, and give the associated Poisson bracket for the two difference equations.

  13. Variation-difference method for solving boundary value problems for linear elliptic complex equations

    Institute of Scientific and Technical Information of China (English)

    2001-01-01

    This paper deals with boundary value problems for linear uniformly elliptic systems. First the general linear uniformly elliptic system of the first order equations is reduced to complex form, and then the compound boundary value problem for the complex equations of the first order is discussed. The approximate solutions of the boundary value problem are found by the variation-difference method, and the error estimates for the approximate solutions are derived.Finally the approximate method of the oblique derivative problem for linear uniformly elliptic equations of the second or der is introduced.

  14. On the weak solution of a three-point boundary value problem for a class of parabolic equations with energy specification

    Directory of Open Access Journals (Sweden)

    Abdelfatah Bouziani

    2003-01-01

    Full Text Available This paper deals with weak solution in weighted Sobolev spaces, of three-point boundary value problems which combine Dirichlet and integral conditions, for linear and quasilinear parabolic equations in a domain with curved lateral boundaries. We, firstly, prove the existence, uniqueness, and continuous dependence of the solution for the linear equation. Next, analogous results are established for the quasilinear problem, using an iterative process based on results obtained for the linear problem.

  15. Functional derivatives, Schr\\"{o}dinger equations, and Feynman integration

    OpenAIRE

    Dynin, Alexander

    2007-01-01

    Schr\\"{o}dinger equations in \\emph{functional derivatives} are solved via quantized Galerkin limit of antinormal functional Feynman integrals for Schr\\"{o}dinger equations in \\emph{partial derivatives

  16. Solving Abel integral equations of first kind via fractional calculus

    Directory of Open Access Journals (Sweden)

    Salman Jahanshahi

    2015-04-01

    Full Text Available We give a new method for numerically solving Abel integral equations of first kind. An estimation for the error is obtained. The method is based on approximations of fractional integrals and Caputo derivatives. Using trapezoidal rule and Computer Algebra System Maple, the exact and approximation values of three Abel integral equations are found, illustrating the effectiveness of the proposed approach.

  17. The Exact General Solution of Painlev\\'e's Sixth Equation (PVI) and The Exact General Solution of the Navier Stokes Equations with Applications to Boundary Layer Problems

    CERN Document Server

    Roman-Miller, Lance Arthur

    2011-01-01

    This paper provides the first known exact general solutions of Painlev\\'e's sixth equation (PVI) and the exact general solutions of the Navier Stokes equations and Prandtl's boundary layer equations.

  18. Non-singular field-only surface integral equations for electromagnetic scattering

    CERN Document Server

    Klaseboer, Evert; Chan, Derek Y C

    2016-01-01

    A boundary integral formulation of electromagnetics that involves only the components of E and H is derived without the use of surface currents that appear in the classical PMCHWT formulation. The kernels of the boundary integral equations for E and H are non-singular so that all field quantities at the surface can be determined to high precision and also geometries with closely spaced surfaces present no numerical difficulties. Quadratic elements can readily be used to represent the surfaces so that the surface integrals can be calculated to higher numerical precision than using planar elements for the same numbers of degrees of freedom.

  19. Exponential integrators for the incompressible Navier-Stokes equations.

    Energy Technology Data Exchange (ETDEWEB)

    Newman, Christopher K.

    2004-07-01

    We provide an algorithm and analysis of a high order projection scheme for time integration of the incompressible Navier-Stokes equations (NSE). The method is based on a projection onto the subspace of divergence-free (incompressible) functions interleaved with a Krylov-based exponential time integration (KBEI). These time integration methods provide a high order accurate, stable approach with many of the advantages of explicit methods, and can reduce the computational resources over conventional methods. The method is scalable in the sense that the computational costs grow linearly with problem size. Exponential integrators, used typically to solve systems of ODEs, utilize matrix vector products of the exponential of the Jacobian on a vector. For large systems, this product can be approximated efficiently by Krylov subspace methods. However, in contrast to explicit methods, KBEIs are not restricted by the time step. While implicit methods require a solution of a linear system with the Jacobian, KBEIs only require matrix vector products of the Jacobian. Furthermore, these methods are based on linearization, so there is no non-linear system solve at each time step. Differential-algebraic equations (DAEs) are ordinary differential equations (ODEs) subject to algebraic constraints. The discretized NSE constitute a system of DAEs, where the incompressibility condition is the algebraic constraint. Exponential integrators can be extended to DAEs with linear constraints imposed via a projection onto the constraint manifold. This results in a projected ODE that is integrated by a KBEI. In this approach, the Krylov subspace satisfies the constraint, hence the solution at the advanced time step automatically satisfies the constraint as well. For the NSE, the projection onto the constraint is typically achieved by a projection induced by the L{sup 2} inner product. We examine this L{sup 2} projection and an H{sup 1} projection induced by the H{sup 1} semi-inner product. The H

  20. Calculation of transonic flows using an extended integral equation method

    Science.gov (United States)

    Nixon, D.

    1976-01-01

    An extended integral equation method for transonic flows is developed. In the extended integral equation method velocities in the flow field are calculated in addition to values on the aerofoil surface, in contrast with the less accurate 'standard' integral equation method in which only surface velocities are calculated. The results obtained for aerofoils in subcritical flow and in supercritical flow when shock waves are present compare satisfactorily with the results of recent finite difference methods.

  1. Integrable dissipative nonlinear second order differential equations via factorizations and Abel equations

    International Nuclear Information System (INIS)

    We emphasize two connections, one well known and another less known, between the dissipative nonlinear second order differential equations and the Abel equations which in their first-kind form have only cubic and quadratic terms. Then, employing an old integrability criterion due to Chiellini, we introduce the corresponding integrable dissipative equations. For illustration, we present the cases of some integrable dissipative Fisher, nonlinear pendulum, and Burgers–Huxley type equations which are obtained in this way and can be of interest in applications. We also show how to obtain Abel solutions directly from the factorization of second order nonlinear equations.

  2. ICM: an Integrated Compartment Method for numerically solving partial differential equations

    Energy Technology Data Exchange (ETDEWEB)

    Yeh, G.T.

    1981-05-01

    An integrated compartment method (ICM) is proposed to construct a set of algebraic equations from a system of partial differential equations. The ICM combines the utility of integral formulation of finite element approach, the simplicity of interpolation of finite difference approximation, and the flexibility of compartment analyses. The integral formulation eases the treatment of boundary conditions, in particular, the Neumann-type boundary conditions. The simplicity of interpolation provides great economy in computation. The flexibility of discretization with irregular compartments of various shapes and sizes offers advantages in resolving complex boundaries enclosing compound regions of interest. The basic procedures of ICM are first to discretize the region of interest into compartments, then to apply three integral theorems of vectors to transform the volume integral to the surface integral, and finally to use interpolation to relate the interfacial values in terms of compartment values to close the system. The Navier-Stokes equations are used as an example of how to derive the corresponding ICM alogrithm for a given set of partial differential equations. Because of the structure of the algorithm, the basic computer program remains the same for cases in one-, two-, or three-dimensional problems.

  3. ICM: an Integrated Compartment Method for numerically solving partial differential equations

    International Nuclear Information System (INIS)

    An integrated compartment method (ICM) is proposed to construct a set of algebraic equations from a system of partial differential equations. The ICM combines the utility of integral formulation of finite element approach, the simplicity of interpolation of finite difference approximation, and the flexibility of compartment analyses. The integral formulation eases the treatment of boundary conditions, in particular, the Neumann-type boundary conditions. The simplicity of interpolation provides great economy in computation. The flexibility of discretization with irregular compartments of various shapes and sizes offers advantages in resolving complex boundaries enclosing compound regions of interest. The basic procedures of ICM are first to discretize the region of interest into compartments, then to apply three integral theorems of vectors to transform the volume integral to the surface integral, and finally to use interpolation to relate the interfacial values in terms of compartment values to close the system. The Navier-Stokes equations are used as an example of how to derive the corresponding ICM alogrithm for a given set of partial differential equations. Because of the structure of the algorithm, the basic computer program remains the same for cases in one-, two-, or three-dimensional problems

  4. Picard-Fuchs Equations, Hauptmoduls and Integrable Systems

    OpenAIRE

    Harnad, J.

    1999-01-01

    The Schwarzian equations satisfied by certain Hauptmoduls (i.e., uniformizing functions for Riemann surfaces of genus zero) are derived from the Picard-Fuchs equations for families of elliptic curves and associated surfaces. The inhomogeneous Picard-Fuchs equations associated to elliptic integrals with varying endpoints are derived and used to determine solutions of equations that are algebraically related to a class of Painlev\\'e VI equations.

  5. The method of duel integral equations for analysis of heat transfer processes

    OpenAIRE

    Mandrik, P. A.

    2001-01-01

    The paper is devoted to the study of a heat conduction equations with mixed boundary conditions on a surface of a researched solid. The method based on the Laplace an Hankel transforms is suggested for the first time to reduce such differential problems to dual integral equations and to obtain the solution in the close form. Method allows to determine analytically regularities of time-space development of appropriate temperature fields.

  6. High Order Numerical Solution of Integral Transport Equation in Slab Geometry

    Institute of Scientific and Technical Information of China (English)

    沈智军; 袁光伟; 沈隆钧

    2002-01-01

    @@ There are some common numerical methods for solving neutron transport equation, which including the well-known discrete ordinates method, PN approximation and integral transport methods[1]. There exists certain singularities in the solution of transport equation near the boundary and interface[2]. It gives rise to the difficulty in the construction of high order accurate numerical methods. The numerical solution obtained by now can not attain the second order convergent accuracy[3,4].

  7. Impulsive neutral fractional integro-differential equations with state dependent delays and integral condition

    Directory of Open Access Journals (Sweden)

    Jaydev Dabas

    2013-12-01

    Full Text Available In this article, we establish the existence of a solution for an impulsive neutral fractional integro-differential state dependent delay equation subject to an integral boundary condition. The existence results are proved by applying the classical fixed point theorems. An example is presented to demonstrate the application of the results established.

  8. Impulsive neutral fractional integro-differential equations with state dependent delays and integral condition

    OpenAIRE

    Jaydev Dabas; Ganga Ram Gautam

    2013-01-01

    In this article, we establish the existence of a solution for an impulsive neutral fractional integro-differential state dependent delay equation subject to an integral boundary condition. The existence results are proved by applying the classical fixed point theorems. An example is presented to demonstrate the application of the results established.

  9. A CLASS OF SINGULARLY PERTURBED INITIAL BOUNDARY PROBLEM FOR REACTION DIFFUSION EQUATION

    Institute of Scientific and Technical Information of China (English)

    Xie Feng

    2003-01-01

    The singularly perturbed initial boundary value problem for a class of reaction diffusion equation isconsidered. Under appropriate conditions, the existence-uniqueness and the asymptotic behavior of the solu-tion are showed by using the fixed-point theorem.

  10. Existence of the Optimal Control for Stochastic Boundary Control Problems Governed by Semilinear Parabolic Equations

    OpenAIRE

    2014-01-01

    We study an optimal control problem governed by a semilinear parabolic equation, whose control variable is contained only in the boundary condition. An existence theorem for the optimal control is obtained.

  11. Asymptotic behavior of solutions to nonlinear parabolic equation with nonlinear boundary conditions

    Directory of Open Access Journals (Sweden)

    Diabate Nabongo

    2008-01-01

    Full Text Available We show that solutions of a nonlinear parabolic equation of second order with nonlinear boundary conditions approach zero as t approaches infinity. Also, under additional assumptions, the solutions behave as a function determined here.

  12. On a periodic boundary value problem for second-order linear functional differential equations

    Czech Academy of Sciences Publication Activity Database

    Mukhigulashvili, Sulkhan

    2005-01-01

    Roč. 3, - (2005), s. 247-261. ISSN 1687-2762 Institutional research plan: CEZ:AV0Z10190503 Keywords : second order linear functional differential equation * periodic boundary value problem * unique solvability Subject RIV: BA - General Mathematics

  13. On solvability of some boundary value problems for a biharmonic equation with periodic conditions

    Science.gov (United States)

    Karachik, Valery V.; Massanov, Saparbay K.; Turmetov, Batirkhan Kh.

    2016-08-01

    In the paper we study questions about solvability of some boundary value problems with periodic conditions for an inhomogeneous biharmonic equation. The exact conditions for solvability of the problems are found.

  14. Periodic solutions of a non-linear wave equation with homogeneous boundary conditions

    Energy Technology Data Exchange (ETDEWEB)

    Rudakov, I A [M.V. Lomonosov Moscow State University, Moscow (Russian Federation)

    2006-02-28

    We prove the existence of time-periodic solutions of a non-linear wave equation with homogeneous boundary conditions. The non-linear term either has polynomial growth or satisfies a 'non-resonance' condition.

  15. Antiperiodic Boundary Value Problems for Second-Order Impulsive Ordinary Differential Equations

    Directory of Open Access Journals (Sweden)

    2009-02-01

    Full Text Available We consider a second-order ordinary differential equation with antiperiodic boundary conditions and impulses. By using Schaefer's fixed-point theorem, some existence results are obtained.

  16. Existence of the Optimal Control for Stochastic Boundary Control Problems Governed by Semilinear Parabolic Equations

    Directory of Open Access Journals (Sweden)

    Weifeng Wang

    2014-01-01

    Full Text Available We study an optimal control problem governed by a semilinear parabolic equation, whose control variable is contained only in the boundary condition. An existence theorem for the optimal control is obtained.

  17. Existence of Solutions for p-Laplace Equations Subjected to Neumann Boundary Value Problem

    Institute of Scientific and Technical Information of China (English)

    HU Zhi-gang; RUI Wen-juan; LIU Wen-bing

    2006-01-01

    The existence of solutions for one dimensional p-Laplace equation (φp(u'))'=f(t,u,u') with t∈(0,1) and φp(s)=│s│p-2s,s≠0 subjected to Neumann boundary value problem at u'(0)=0,u'(1)=0. By using the degree theory, the sufficient conditions of the existence of solutions for p-Laplace equation subjected to Neumann boundary value condition are established.

  18. Sobolev type equations of time-fractional order with periodical boundary conditions

    Science.gov (United States)

    Plekhanova, Marina

    2016-08-01

    The existence of a unique local solution for a class of time-fractional Sobolev type partial differential equations endowed by the Cauchy initial conditions and periodical with respect to every spatial variable boundary conditions on a parallelepiped is proved. General results are applied to study of the unique solvability for the initial boundary value problem to Benjamin-Bona-Mahony-Burgers and Allair partial differential equations.

  19. Partial differential equations of mathematical physics and integral equations

    CERN Document Server

    Guenther, Ronald B

    1996-01-01

    This book was written to help mathematics students and those in the physical sciences learn modern mathematical techniques for setting up and analyzing problems. The mathematics used is rigorous, but not overwhelming, while the authors carefully model physical situations, emphasizing feedback among a beginning model, physical experiments, mathematical predictions, and the subsequent refinement and reevaluation of the physical model itself. Chapter 1 begins with a discussion of various physical problems and equations that play a central role in applications. The following chapters take up the t

  20. On singular solutions of a magnetohydrodynamic nonlinear boundary layer equation

    Directory of Open Access Journals (Sweden)

    Mohammed Guedda

    2007-05-01

    Full Text Available This paper concerns the singular solutions of the equation $$ f''' +kappa ff''-eta {f'}^2 = 0, $$ where $eta < 0$ and $kappa = 0$ or 1. This equation arises when modelling heat transfer past a vertical flat plate embedded in a saturated porous medium with an applied magnetic field. After suitable normalization, $f'$ represents the velocity parallel to the surface or the non-dimensional fluid temperature. Our interest is in solutions which develop a singularity at some point (the blow-up point. In particular, we shall examine in detail the behavior of $f$ near the blow-up point.

  1. A simple multigrid scheme for solving the Poisson equation with arbitrary domain boundaries

    CERN Document Server

    Guillet, Thomas

    2011-01-01

    We present a new multigrid scheme for solving the Poisson equation with Dirichlet boundary conditions on a Cartesian grid with irregular domain boundaries. This scheme was developed in the context of the Adaptive Mesh Refinement (AMR) schemes based on a graded-octree data structure. The Poisson equation is solved on a level-by-level basis, using a "one-way interface" scheme in which boundary conditions are interpolated from the previous coarser level solution. Such a scheme is particularly well suited for self-gravitating astrophysical flows requiring an adaptive time stepping strategy. By constructing a multigrid hierarchy covering the active cells of each AMR level, we have designed a memory-efficient algorithm that can benefit fully from the multigrid acceleration. We present a simple method for capturing the boundary conditions across the multigrid hierarchy, based on a second-order accurate reconstruction of the boundaries of the multigrid levels. In case of very complex boundaries, small scale features ...

  2. Riemann boundary value problem for triharmonic equation in higher space.

    Science.gov (United States)

    Gu, Longfei

    2014-01-01

    We mainly deal with the boundary value problem for triharmonic function with value in a universal Clifford algebra: Δ(3)[u](x) = 0, x ∈ R (n)\\∂Ω, u (+)(x) = u (-)(x)G(x) + g(x), x ∈ ∂Ω, (D (j) u)(+)(x) = (D (j) u)(-)(x)A j + f j (x), x ∈ ∂Ω, u(∞) = 0, where (j = 1,…, 5)  ∂Ω is a Lyapunov surface in R (n) , D = ∑ k=1 (n) e k (∂/∂x k) is the Dirac operator, and u(x) = ∑ A e A u A (x) are unknown functions with values in a universal Clifford algebra Cl(V n,n). Under some hypotheses, it is proved that the boundary value problem has a unique solution. PMID:25114963

  3. Positivity for equations involving polyharmonic operators with Dirichlet boundary conditions

    NARCIS (Netherlands)

    Grunau, H.-Ch.; Sweers, G.

    2001-01-01

    Cranston, Fabes and Zhao ([26], [5]) established the uniform bound sup x; y 2 x 6= y R G1;n (x; z)G1;n (z; y) dz G1;n (x; y) M < 1; (1) where G1;n (x; y) is the Green function for the Laplacian - with Dirichlet boundary conditions on a Lipschitz domain - Rn with n 3 (see [27] for n = 2). This

  4. Solvability of functional quadratic integral equations with perturbation

    OpenAIRE

    Mohamed M. A. Metwali

    2013-01-01

    We study the existence of solutions of the functional quadratic integral equation with a perturbation term in the space of Lebesgue integrable functions on an unbounded interval by using the Krasnoselskii fixed point theory and the measure of weak noncompactness.

  5. Existence of Solutions of Periodic Boundary Value Problems for Impulsive Functional Duffing Equations at Nonresonance Case

    Directory of Open Access Journals (Sweden)

    Liu Yuji

    2008-01-01

    Full Text Available Abstract This paper deals with the existence of solutions of the periodic boundary value problem of the impulsive Duffing equations: . Sufficient conditions are established for the existence of at least one solution of above-mentioned boundary value problem. Our method is based upon Schaeffer's fixed-point theorem. Examples are presented to illustrate the efficiency of the obtained results.

  6. Recovery of small electromagnetic inhomogeneities from boundary measurements in time-dependent Maxwell's equations

    CERN Document Server

    Daveau, Christian

    2007-01-01

    We consider for the time-dependent Maxwell's equations the inverse problem of identifying locations and certain properties of small electromagnetic inhomogeneities in a homogeneous background medium from dynamic measurements of the tangential component of the magnetic field on the boundary (or a part of the boundary) of a domain.

  7. The immersed boundary method for the (2D) incompressible Navier-Stokes equations

    NARCIS (Netherlands)

    Meûlen, R.J.R. vander

    2006-01-01

    Immersed Boundary Methods (IBMs) are a class of methods in Computational Fluid Dynamics where the grids do not conform to the shape of the body. Instead they employ Cartesian meshes and alternative ways to incorporate the boundary conditions in the (discrete) governing equations. The simple grids an

  8. m-POINT BOUNDARY VALUE PROBLEM FOR SECOND ORDER IMPULSIVE DIFFERENTIAL EQUATION AT RESONANCE

    Institute of Scientific and Technical Information of China (English)

    2012-01-01

    In his paper,we obtain a general theorem concerning the existence of solutions to an m-point boundary value problem for the second-order differential equation with impulses.Moreover,the result can also be applied to study the usual m-point boundary value problem at resonance without impulses.

  9. GENERAL DECAY OF A TRANSMISSION PROBLEM FOR KIRCHHOFF TYPE WAVE EQUATIONS WITH BOUNDARY MEMORY CONDITION

    Institute of Scientific and Technical Information of China (English)

    Sun Hye PARK

    2014-01-01

    In this paper, we investigate the influence of boundary dissipation on the de-cay property of solutions for a transmission problem of Kirchhoff type wave equation with boundary memory condition. By introducing suitable energy and Lyapunov functionals, we establish a general decay estimate for the energy, which depends on the behavior of relaxation function.

  10. Finite-Volume Analysis for the Cahn-Hilliard equation with Dynamic boundary conditions

    OpenAIRE

    Nabet, Flore

    2014-01-01

    This work is devoted to the convergence analysis of a finite-volume approximation of the 2D Cahn-Hilliard equation with dynamic boundary conditions. The method that we propose couples a 2d-finite-volume method in a bounded, smooth domain and a 1d-finite-volume method on its boundary. We prove convergence of the sequence of approximate solutions.

  11. The BV solution of the parabolic equation with degeneracy on the boundary

    Directory of Open Access Journals (Sweden)

    Zhan Huashui

    2016-12-01

    Full Text Available Consider a parabolic equation which is degenerate on the boundary. By the degeneracy, to assure the well-posedness of the solutions, only a partial boundary condition is generally necessary. When 1 ≤ α < p – 1, the existence of the local BV solution is proved. By choosing some kinds of test functions, the stability of the solutions based on a partial boundary condition is established.

  12. Sliding periodic boundary conditions for lattice Boltzmann and lattice kinetic equations

    OpenAIRE

    Adhikari, R.; Desplat, J. -C.; Stratford, K.

    2005-01-01

    We present a method to impose linear shear flow in discrete-velocity kinetic models of hydrodynamics through the use of sliding periodic boundary conditions. Our method is derived by an explicit coarse-graining of the Lees-Edwards boundary conditions for Couette flow in molecular dynamics, followed by a projection of the resulting equations onto the subspace spanned by the discrete velocities of the lattice Boltzmann method. The boundary conditions are obtained without resort to perturbative ...

  13. RADIATION BOUNDARY CONDITIONS FOR MAXWELL'S EQUATIONS: A REVIEW OF ACCURATE TIME-DOMAIN FORMULATIONS

    Institute of Scientific and Technical Information of China (English)

    Thomas Hagstrom; Stephen Lau

    2007-01-01

    We review time-domain formulations of radiation boundary conditions for Maxwell's equations, focusing on methods which can deliver arbitrary accuracy at acceptable computational cost. Examples include fast evaluations of nonlocal conditions on symmetric and general boundaries, methods based on identifying and evaluating equivalent sources, and local approximations such as the perfectly matched layer and sequences of local boundary conditions. Complexity estimates are derived to assess work and storage requirements as a function of wavelength and simulation time.

  14. Run-and-tumble particles, telegrapher’s equation and absorption problems with partially reflecting boundaries

    International Nuclear Information System (INIS)

    Absorption problems of run-and-tumble particles, described by the telegrapher's equation, are analyzed in one space dimension considering partially reflecting boundaries. Exact expressions for the probability distribution function in the Laplace domain and for the mean time to absorption are given, discussing some interesting limits (Brownian and wave limit, large volume limit) and different case studies (semi-infinite segment, equal and symmetric boundaries, totally/partially reflecting boundaries). (paper)

  15. Initial-boundary value problems for quasilinear dispersive equations posed on a bounded interval

    Directory of Open Access Journals (Sweden)

    Andrei V. Faminskii

    2010-01-01

    Full Text Available This paper studies nonhomogeneous initial-boundary value problems for quasilinear one-dimensional odd-order equations posed on a bounded interval. For reasonable initial and boundary conditions we prove existence and uniqueness of global weak and regular solutions. Also we show the exponential decay of the obtained solution with zero boundary conditions and right-hand side, and small initial data.

  16. Boundary value problems for the Helmholtz equation in a half-plane

    OpenAIRE

    Chandler-Wilde, SN

    1994-01-01

    The Dirichlet and impedance boundary value problems for the Helmholtz equation in a half-plane with bounded continuous boundary data are studied. For the Dirichlet problem the solution can be constructed explicitly. We point out that, for wavenumbers k > 0, the solution, although it satisfies a limiting absorption principle, may increase in magnitude with distance from the boundary. Using the explicit solution we propose a novel radiation condition which we utilise in formulating the impedanc...

  17. Implementation of higher-order absorbing boundary conditions for the Einstein equations

    International Nuclear Information System (INIS)

    We present an implementation of absorbing boundary conditions for the Einstein equations based on the recent work of Buchman and Sarbach. In this paper, we assume that spacetime may be linearized about Minkowski space close to the outer boundary, which is taken to be a coordinate sphere. We reformulate the boundary conditions as conditions on the gauge-invariant Regge-Wheeler-Zerilli scalars. Higher-order radial derivatives are eliminated by rewriting the boundary conditions as a system of ODEs for a set of auxiliary variables intrinsic to the boundary. From these we construct boundary data for a set of well-posed constraint-preserving boundary conditions for the Einstein equations in a first-order generalized harmonic formulation. This construction has direct applications to outer boundary conditions in simulations of isolated systems (e.g., binary black holes) as well as to the problem of Cauchy-perturbative matching. As a test problem for our numerical implementation, we consider linearized multipolar gravitational waves in TT gauge, with angular momentum numbers l = 2 (Teukolsky waves), 3 and 4. We demonstrate that the perfectly absorbing boundary condition ΒL of order L = l yields no spurious reflections to linear order in perturbation theory. This is in contrast to the lower-order absorbing boundary conditions ΒL with L 0 boundary condition that imposes the vanishing of the Newman-Penrose scalar Ψ0.

  18. Approximation of 2D Euler Equations by the Second-Grade Fluid Equations with Dirichlet Boundary Conditions

    Science.gov (United States)

    Lopes Filho, Milton C.; Nussenzveig Lopes, Helena J.; Titi, Edriss S.; Zang, Aibin

    2015-06-01

    The second-grade fluid equations are a model for viscoelastic fluids, with two parameters: α > 0, corresponding to the elastic response, and , corresponding to viscosity. Formally setting these parameters to 0 reduces the equations to the incompressible Euler equations of ideal fluid flow. In this article we study the limits of solutions of the second-grade fluid system, in a smooth, bounded, two-dimensional domain with no-slip boundary conditions. This class of problems interpolates between the Euler- α model (), for which the authors recently proved convergence to the solution of the incompressible Euler equations, and the Navier-Stokes case ( α = 0), for which the vanishing viscosity limit is an important open problem. We prove three results. First, we establish convergence of the solutions of the second-grade model to those of the Euler equations provided , as α → 0, extending the main result in (Lopes Filho et al., Physica D 292(293):51-61, 2015). Second, we prove equivalence between convergence (of the second-grade fluid equations to the Euler equations) and vanishing of the energy dissipation in a suitably thin region near the boundary, in the asymptotic regime , as α → 0. This amounts to a convergence criterion similar to the well-known Kato criterion for the vanishing viscosity limit of the Navier-Stokes equations to the Euler equations. Finally, we obtain an extension of Kato's classical criterion to the second-grade fluid model, valid if , as . The proof of all these results relies on energy estimates and boundary correctors, following the original idea by Kato.

  19. Integral Formulation of the Boundary Value Problem in Waveguides.

    Science.gov (United States)

    Sancho, M.

    1980-01-01

    Presents an integral approach to the boundary value problem in waveguides deduced from the Kirchoff's integral formulation of the electromagnetic field. Also, the basis for the numerical solution of more general problems is given, including the example of the isosceles right triangular guide. (Author/SK)

  20. The solution of integral equation of plasma diffusion problems

    International Nuclear Information System (INIS)

    Two topics which concern with the solution of one dimensional diffusion equation are discussed in this paper. The solution of the one dimensional diffusion equation with one initial condition is discussed in the first topic, and the solution of the same equation with two initial conditions is discussed in the second topic. Both solutions are presented in integral forms. (author)

  1. The boundary value problems for the scalar Oseen equation

    Czech Academy of Sciences Publication Activity Database

    Medková, Dagmar; Skopin, E.; Varnhorn, W.

    2012-01-01

    Roč. 285, 17-18 (2012), s. 2208-2221. ISSN 0025-584X R&D Projects: GA ČR(CZ) GAP201/11/1304 Institutional support: RVO:67985840 Keywords : scalar Oseen equation * Dirichlet problem * Neumann problem Subject RIV: BA - General Mathematics Impact factor: 0.576, year: 2012 http://onlinelibrary.wiley.com/doi/10.1002/mana.201100219/abstract

  2. A stable penalty method for the compressible Navier-Stokes equations: I. Open boundary conditions

    DEFF Research Database (Denmark)

    Hesthaven, Jan; Gottlieb, D.

    1996-01-01

    The purpose of this paper is to present asymptotically stable open boundary conditions for the numerical approximation of the compressible Navier-Stokes equations in three spatial dimensions. The treatment uses the conservation form of the Navier-Stokes equations and utilizes linearization and...

  3. On The Stabilization of the Linear Kawahara Equation with Periodic Boundary Conditions

    Directory of Open Access Journals (Sweden)

    Patricia N. da Silva

    2015-03-01

    Full Text Available We study the stabilization of global solutions of the linear Kawahara equation (K with periodic boundary conditions under the effect of a localized damping mechanism. The Kawahara equation is a model for small amplitude long waves. Using separation of variables, the Ingham inequality, multiplier techniques and compactness arguments we prove the exponential decay of the solutions of the (K model.

  4. On The Stabilization of the Linear Kawahara Equation with Periodic Boundary Conditions

    OpenAIRE

    Patricia N. da Silva; Carlos F. Vasconcellos

    2015-01-01

    We study the stabilization of global solutions of the linear Kawahara equation (K) with periodic boundary conditions under the effect of a localized damping mechanism. The Kawahara equation is a model for small amplitude long waves. Using separation of variables, the Ingham inequality, multiplier techniques and compactness arguments we prove the exponential decay of the solutions of the (K) model.

  5. Strang-type preconditioners for solving fractional diffusion equations by boundary value methods

    NARCIS (Netherlands)

    Gu, Xian-Ming; Huang, Ting-Zhu; Zhao, Xi-Le; Li, Hou-Biao; Li, Liang

    2015-01-01

    The finite difference scheme with the shifted Grünwarld formula is employed to semi-discrete the fractional diffusion equations. This spatial discretization can reduce to the large system of ordinary differential equations (ODEs) with initial values. Recently, boundary value method (BVM) was develop

  6. On the integrability of nonlinear partial differential equations

    CERN Document Server

    Dorren, H J S

    1998-01-01

    We investigate the integrability of Nonlinear Partial Differential Equations (NPDEs). The concepts are developed by firstly discussing the integrability of the KdV equation. We proceed by generalizing the ideas introduced for the KdV equation to other NPDEs. The method is based upon a linearization principle which can be applied on nonlinearities which have a polynomial form. We illustrate the potential of the method by finding solutions of the (coupled) nonlinear Schrödinger equation and the Manakov equation which play an important role in optical fiber communication. Finally, it is shown that the method can also be generalized to higher-dimensions.

  7. Shape integral method for magnetospheric shapes. [boundary layer calculations

    Science.gov (United States)

    Michel, F. C.

    1979-01-01

    A method is developed for calculating the shape of any magnetopause to arbitrarily high precision. The method uses an integral equation which is evaluated for a trial shape. The resulting values of the integral equation as a function of auxiliary variables indicate how close one is to the desired solution. A variational method can then be used to improve the trial shape. Some potential applications are briefly mentioned.

  8. THE DIFFERENTIAL INTEGRAL EQUATIONS ON SMOOTH CLOSED ORIENTABLE MANIFOLDS

    Institute of Scientific and Technical Information of China (English)

    2001-01-01

    Using integration by parts and Stokes' formula, the authors give a new definition of Hadamard principal value of higher order singular integrals with Bochner-Martinelli kernel on smooth closed orientable manifolds in Cn. The Plemelj formula and composite formula of higher order singular integral arc obtained. Differential integral equations on smooth closed orientable manifolds are treated by using the composite formula.

  9. Momentum conservation law in explicit symplectic integrators for a nonlinear Schrödinger-type equation

    International Nuclear Information System (INIS)

    This paper discusses the momentum conservation law in a symplectic integrator for a nonlinear Schrödinger-type equation. We show that the total momentum, which is formally expressed by a polynomial of a discrete variable, is conserved under cyclic boundary conditions. We also perform numerical simulations to demonstrate the validity and numerical convergence of our expression for the total momentum. As a result, the symplectic integrator simultaneously satisfies the energy, density and momentum conservation laws. (author)

  10. Existence and uniqueness for one-phase Stefan problems of non-classical heat equations with temperature boundary condition at a fixed face

    Directory of Open Access Journals (Sweden)

    Adriana C. Briozzo

    2006-02-01

    Full Text Available We prove the existence and uniqueness, local in time, of a solution for a one-phase Stefan problem of a non-classical heat equation for a semi-infinite material with temperature boundary condition at the fixed face. We use the Friedman-Rubinstein integral representation method and the Banach contraction theorem in order to solve an equivalent system of two Volterra integral equations.

  11. Well posed constraint-preserving boundary conditions for the linearized Einstein equations

    CERN Document Server

    Calabrese, G; Reula, O; Sarbach, O; Tiglio, M H; Calabrese, Gioel; Pullin, Jorge; Reula, Oscar; Sarbach, Olivier; Tiglio, Manuel

    2003-01-01

    In the Cauchy problem of general relativity one considers initial data that satisfies certain constraints. The evolution equations guarantee that the evolved variables will satisfy the constraints at later instants of time. This is only true within the domain of dependence of the initial data. If one wishes to consider situations where the evolutions are studied for longer intervals than the size of the domain of dependence, as is usually the case in three dimensional numerical relativity, one needs to give boundary data. The boundary data should be specified in such a way that the constraints are satisfied everywhere, at all times. In this paper we address this problem for the case of general relativity linearized around Minkowski space using the generalized Einstein-Christoffel symmetric hyperbolic system of evolution equations. We study the evolution equations for the constraints, specify boundary conditions for them that make them well posed and further choose these boundary conditions in such a way that ...

  12. Exact integral equation for the renormalized Fermi surface

    OpenAIRE

    Ledowski, Sascha; Kopietz, Peter

    2002-01-01

    The true Fermi surface of a fermionic many-body system can be viewed as a fixed point manifold of the renormalization group (RG). Within the framework of the exact functional RG we show that the fixed point condition implies an exact integral equation for the counterterm which is needed for a self-consistent calculation of the Fermi surface. In the simplest approximation, our integral equation reduces to the self-consistent Hartree-Fock equation for the counterterm.

  13. Boundary element method for the solution of the diffusion equation in cylindrical symmetry

    International Nuclear Information System (INIS)

    Equations for the solution of the diffusion equation in plane Cartesian geometry with the Boundary Element method was derived. The equation for the axi-symmetric case were set and included in the computer program. The results were compared to those obtained by the Finite Difference method. Comparing the results some advantages of the proposed method can be observed, with implications on the multidimensional problems. (author)

  14. THE INITIAL BOUNDARY VALUE PROBLEM FOR QUASI-LINEAR SCHR(O)DINGER-POISSON EQUATIONS

    Institute of Scientific and Technical Information of China (English)

    2006-01-01

    In this article, the author studies the initial-(Dirichlet) boundary problem for a high-field version of the Schr(o)dinger-Poisson equations, which include a nonlinear term in the Poisson equation corresponding to a field-dependent dielectric constant and an effective potential in the Schrodinger equations on the unit cube. A global existence and uniqueness is established for a solution to this problem.

  15. Incompressible limit of the compressible magnetohydrodynamic equations with periodic boundary conditions

    OpenAIRE

    Jiang, Song; Ju, Qiangchang; Li, Fucai

    2010-01-01

    This paper is concerned with the incompressible limit of the compressible magnetohydrodynamic equations with periodic boundary conditions. It is rigorously shown that the weak solutions of the compressible magnetohydrodynamic equations converge to the strong solution of the viscous or inviscid incompressible magnetohydrodynamic equations as long as the latter exists both for the well-prepared initial data and general initial data. Furthermore, the convergence rates are also obtained in the ca...

  16. Existence and uniqueness of entropy solution to initial boundary value problem for the inviscid Burgers equation

    CERN Document Server

    Zhu, C

    2003-01-01

    This paper is concerned with the existence and uniqueness of the entropy solution to the initial boundary value problem for the inviscid Burgers equation. To apply the method of vanishing viscosity to study the existence of the entropy solution, we first introduce the initial boundary value problem for the viscous Burgers equation, and as in Evans (1998 Partial Differential Equations (Providence, RI: American Mathematical Society) and Hopf (1950 Commun. Pure Appl. Math. 3 201-30), give the formula of the corresponding viscosity solutions by Hopf-Cole transformation. Secondly, we prove the convergence of the viscosity solution sequences and verify that the limiting function is an entropy solution. Finally, we give an example to show how our main result can be applied to solve the initial boundary value problem for the Burgers equation.

  17. Existence and uniqueness of entropy solution to initial boundary value problem for the inviscid Burgers equation

    International Nuclear Information System (INIS)

    This paper is concerned with the existence and uniqueness of the entropy solution to the initial boundary value problem for the inviscid Burgers equation. To apply the method of vanishing viscosity to study the existence of the entropy solution, we first introduce the initial boundary value problem for the viscous Burgers equation, and as in Evans (1998 Partial Differential Equations (Providence, RI: American Mathematical Society) and Hopf (1950 Commun. Pure Appl. Math. 3 201-30), give the formula of the corresponding viscosity solutions by Hopf-Cole transformation. Secondly, we prove the convergence of the viscosity solution sequences and verify that the limiting function is an entropy solution. Finally, we give an example to show how our main result can be applied to solve the initial boundary value problem for the Burgers equation

  18. Existence and uniqueness of entropy solution to initial boundary value problem for the inviscid Burgers equation

    Energy Technology Data Exchange (ETDEWEB)

    Zhu, Changjiang; Duan, Renjun [Laboratory of Nonlinear Analysis, Department of Mathematics, Central China Normal University, Wuhan 430079, People' s Republic of China (China)

    2003-02-28

    This paper is concerned with the existence and uniqueness of the entropy solution to the initial boundary value problem for the inviscid Burgers equation. To apply the method of vanishing viscosity to study the existence of the entropy solution, we first introduce the initial boundary value problem for the viscous Burgers equation, and as in Evans (1998 Partial Differential Equations (Providence, RI: American Mathematical Society) and Hopf (1950 Commun. Pure Appl. Math. 3 201-30), give the formula of the corresponding viscosity solutions by Hopf-Cole transformation. Secondly, we prove the convergence of the viscosity solution sequences and verify that the limiting function is an entropy solution. Finally, we give an example to show how our main result can be applied to solve the initial boundary value problem for the Burgers equation.

  19. Calculation of unsteady transonic flows using the integral equation method

    Science.gov (United States)

    Nixon, D.

    1978-01-01

    The basic integral equations for a harmonically oscillating airfoil in a transonic flow with shock waves are derived; the reduced frequency is assumed to be small. The problems associated with shock wave motion are treated using a strained coordinate system. The integral equation is linear and consists of both line integrals and surface integrals over the flow field which are evaluated by quadrature. This leads to a set of linear algebraic equations that can be solved directly. The shock motion is obtained explicitly by enforcing the condition that the flow is continuous except at a shock wave. Results obtained for both lifting and nonlifting oscillatory flows agree satisfactorily with other accurate results.

  20. Mixed initial-boundary value problem for equations of motion of Kelvin-Voigt fluids

    Science.gov (United States)

    Baranovskii, E. S.

    2016-07-01

    The initial-boundary value problem for equations of motion of Kelvin-Voigt fluids with mixed boundary conditions is studied. The no-slip condition is used on some portion of the boundary, while the impermeability condition and the tangential component of the surface force field are specified on the rest of the boundary. The global-in-time existence of a weak solution is proved. It is shown that the solution is unique and depends continuously on the field of external forces, the field of surface forces, and initial data.

  1. Mixed Boundary Value Problems for Stationary Magnetohydrodynamic Equations of a Viscous Heat-Conducting Fluid

    Science.gov (United States)

    Alekseev, Gennady

    2016-04-01

    We consider the boundary value problem for stationary magnetohydrodynamic equations of electrically and heat conducting fluid under inhomogeneous mixed boundary conditions for electromagnetic field and temperature and Dirichlet condition for the velocity. The problem describes the thermoelectromagnetic flow of a viscous fluid in 3D bounded domain with the boundary consisting of several parts with different thermo- and electrophysical properties. The global solvability of the boundary value problem is proved and the apriori estimates of the solution are derived. The sufficient conditions on the data are established which provide a local uniqueness of the solution.

  2. Homotopic Mapping Solving Method for Perturbed Mechanism of Western Boundary Undercurrents in Equator Pacific

    Institute of Scientific and Technical Information of China (English)

    MO Jiaqi; WANG Hui; LIN Wantao

    2006-01-01

    The perturbed boundary undercurrent is an exceptional event in the tropical atmosphere and ocean. It is a complicated nonlinear system. Its appearance badly affects not only natural conditions such as climate and environment,but also global economic development and human living, and brings about many calamities. Thus there is very attractive study on its rules in the international academic circles. Many scholars made more studies on its local and whole behaviors using different methods, such as self-anamnestic principle, Fokker-Plank Equation method, higher order singular pedigree and predictable study, rapid change on boundary, indeterminate adaptive control, multi-cogradient method and so on. Nonlinear perturbed theory and approximate method are very attractive studies in the international academic circles. Many scholars considered a class of nonlinear problems for the ordinary differential equation, the reaction diffusion equations, the boundary value of elliptic equation, the initial boundary value of hyperbolic equation, the shock layer solution of nonlinear equation and so on. In this paper, a class of perturbed mechanism for the western boundary undercurrents in the equator Pacific is considered. Under suitable conditions, using a homotopic mapping theory and method, we obtain a simple and rapid arbitrary order approximate solution for the corresponding nonlinear system. For example, a special case shows that using the homotopic mapping method, there is a high accuracy for the computed value.It is also provided from the results that the solution for homotopic mapping solving method can be used for analyzing operator for perturbed mechanism of western boundary undercurrents in the equator Pacific.

  3. Existence of solutions to boundary value problems arising from the fractional advection dispersion equation

    Directory of Open Access Journals (Sweden)

    Lingju Kong

    2013-04-01

    Full Text Available We study the existence of multiple solutions to the boundary value problem $$displaylines{ frac{d}{dt}Big(frac12{}_0D_t^{-eta}(u'(t+frac12{}_tD_T^{-eta}(u'(t Big+lambda abla F(t,u(t=0,quad tin [0,T],cr u(0=u(T=0, }$$ where $T>0$, $lambda>0$ is a parameter, $0leqeta<1$, ${}_0D_t^{-eta}$ and ${}_tD_T^{-eta}$ are, respectively, the left and right Riemann-Liouville fractional integrals of order $eta$, $F: [0,T]imesmathbb{R}^Nomathbb{R}$ is a given function. Our interest in the above system arises from studying the steady fractional advection dispersion equation. By applying variational methods, we obtain sufficient conditions under which the above equation has at least three solutions. Our results are new even for the special case when $eta=0$. Examples are provided to illustrate the applicability of our results.

  4. Difference equations for correlation functions of A{sup (1)}{sub n-1}-face model with boundary reflection

    Energy Technology Data Exchange (ETDEWEB)

    Quano, Yas-Hiro [Department of Medical Electronics, Suzuka University of Medical Science, Kishioka-cho, Suzuka (Japan)]. E-mail: quanoy@suzuka-u.ac.jp

    2001-10-19

    The A{sup (1)}{sub n-1}-face model with boundary reflection is considered on the basis of the boundary corner transfer matrix bootstrap. We construct the fused boundary Boltzmann weights to determine the normalization factor. We derive difference equations of the quantum Knizhnik-Zamolodchikov type for correlation functions of the boundary model. The simplest difference equations are solved in the case of the free boundary condition. (author)

  5. Time-independent integral equation for Maxwell's system. Application of radar cross section computation

    International Nuclear Information System (INIS)

    We prove that the scattering operator for the wave equation in the exterior of an non-homogeneous obstacle exists. Its distribution kernel is represented by a time-dependent boundary integral equation. A space-time integral variational formulation is developed for determining the current induced by the scattering of an electromagnetic wave by an homogeneous object. The discrete approximation of the variational problem using a finite element method in both space and time leads to stable convergent schemes, giving a numerical code for perfectly conducting cylinders. (author)

  6. Integrable discretisations for a class of nonlinear Schrödinger equations on Grassmann algebras

    International Nuclear Information System (INIS)

    Integrable discretisations for a class of coupled (super) nonlinear Schrödinger (NLS) type of equations are presented. The class corresponds to a Lax operator with entries in a Grassmann algebra. Elementary Darboux transformations are constructed. As a result, Grassmann generalisations of the Toda lattice and the NLS dressing chain are obtained. The compatibility (Bianchi commutativity) of these Darboux transformations leads to integrable Grassmann generalisations of the difference Toda and NLS equations. The resulting systems will have discrete Lax representations provided by the set of two consistent elementary Darboux transformations. For the two discrete systems obtained, initial value and initial-boundary problems are formulated.

  7. Integrable discretisations for a class of nonlinear Schrödinger equations on Grassmann algebras

    Energy Technology Data Exchange (ETDEWEB)

    Grahovski, Georgi G., E-mail: G.Grahovski@leeds.ac.uk [Department of Applied Mathematics, University of Leeds, Woodhouse Lane, Leeds, LS2 9JT (United Kingdom); Institute of Nuclear Research and Nuclear Energy, Bulgarian Academy of Sciences, 72 Tsarigradsko chausee, Sofia 1784 (Bulgaria); Mikhailov, Alexander V., E-mail: A.V.Mikhailov@leeds.ac.uk [Department of Applied Mathematics, University of Leeds, Woodhouse Lane, Leeds, LS2 9JT (United Kingdom)

    2013-12-17

    Integrable discretisations for a class of coupled (super) nonlinear Schrödinger (NLS) type of equations are presented. The class corresponds to a Lax operator with entries in a Grassmann algebra. Elementary Darboux transformations are constructed. As a result, Grassmann generalisations of the Toda lattice and the NLS dressing chain are obtained. The compatibility (Bianchi commutativity) of these Darboux transformations leads to integrable Grassmann generalisations of the difference Toda and NLS equations. The resulting systems will have discrete Lax representations provided by the set of two consistent elementary Darboux transformations. For the two discrete systems obtained, initial value and initial-boundary problems are formulated.

  8. Master equations and the theory of stochastic path integrals

    CERN Document Server

    Weber, Markus F

    2016-01-01

    This review provides a pedagogic and self-contained introduction to master equations and to their representation by path integrals. We discuss analytical and numerical methods for the solution of master equations, keeping our focus on methods that are applicable even when stochastic fluctuations are strong. The reviewed methods include the generating function technique and the Poisson representation, as well as novel ways of mapping the forward and backward master equations onto linear partial differential equations (PDEs). Spectral methods, WKB approximations, and a variational approach have been proposed for the analysis of the PDE obeyed by the generating function. After outlining these methods, we solve the derived PDEs in terms of two path integrals. The path integrals provide distinct exact representations of the conditional probability distribution solving the master equations. We exemplify both path integrals in analysing elementary chemical reactions. Furthermore, we review a method for the approxima...

  9. POSITIVE SOLUTIONS OF A NONLINEAR THREE-POINT EIGENVALUE PROBLEM WITH INTEGRAL BOUNDARY CONDITIONS

    Directory of Open Access Journals (Sweden)

    FAOUZI HADDOUCHI

    2015-11-01

    Full Text Available In this paper, we study the existence of positive solutions of a three-point integral boundary value problem (BVP for the following second-order differential equation u''(t + \\lambda a(tf(u(t = 0; 0 0 is a parameter, 0 <\\eta < 1, 0 <\\alpha < 1/{\\eta}. . By using the properties of the Green's function and Krasnoselskii's fixed point theorem on cones, the eigenvalue intervals of the nonlinear boundary value problem are considered, some sufficient conditions for the existence of at least one positive solutions are established.

  10. Radiation boundary condition and anisotropy correction for finite difference solutions of the Helmholtz equation

    Science.gov (United States)

    Tam, Christopher K. W.; Webb, Jay C.

    1994-01-01

    In this paper finite-difference solutions of the Helmholtz equation in an open domain are considered. By using a second-order central difference scheme and the Bayliss-Turkel radiation boundary condition, reasonably accurate solutions can be obtained when the number of grid points per acoustic wavelength used is large. However, when a smaller number of grid points per wavelength is used excessive reflections occur which tend to overwhelm the computed solutions. Excessive reflections are due to the incompability between the governing finite difference equation and the Bayliss-Turkel radiation boundary condition. The Bayliss-Turkel radiation boundary condition was developed from the asymptotic solution of the partial differential equation. To obtain compatibility, the radiation boundary condition should be constructed from the asymptotic solution of the finite difference equation instead. Examples are provided using the improved radiation boundary condition based on the asymptotic solution of the governing finite difference equation. The computed results are free of reflections even when only five grid points per wavelength are used. The improved radiation boundary condition has also been tested for problems with complex acoustic sources and sources embedded in a uniform mean flow. The present method of developing a radiation boundary condition is also applicable to higher order finite difference schemes. In all these cases no reflected waves could be detected. The use of finite difference approximation inevita bly introduces anisotropy into the governing field equation. The effect of anisotropy is to distort the directional distribution of the amplitude and phase of the computed solution. It can be quite large when the number of grid points per wavelength used in the computation is small. A way to correct this effect is proposed. The correction factor developed from the asymptotic solutions is source independent and, hence, can be determined once and for all. The

  11. Multicomponent integrable wave equations: II. Soliton solutions

    International Nuclear Information System (INIS)

    The Darboux-dressing transformations developed in Degasperis and Lombardo (2007 J. Phys. A: Math. Theor. 40 961-77) are here applied to construct soliton solutions for a class of boomeronic-type equations. The vacuum (i.e. vanishing) solution and the generic plane wave solution are both dressed to yield one-soliton solutions. The formulae are specialized to the particularly interesting case of the resonant interaction of three waves, a well-known model which is of boomeronic type. For this equation a novel solution which describes three locked dark pulses (simulton) is introduced.

  12. Multicomponent integrable wave equations: II. Soliton solutions

    Energy Technology Data Exchange (ETDEWEB)

    Degasperis, A [Dipartimento di Fisica, Universita di Roma ' La Sapienza' , and Istituto Nazionale di Fisica Nucleare, Sezione di Roma, Rome (Italy); Lombardo, S [School of Mathematics, University of Manchester, Alan Turing Building, Upper Brook Street, Manchester M13 9EP (United Kingdom)], E-mail: antonio.degasperis@roma1.infn.it, E-mail: sara.lombardo@manchester.ac.uk, E-mail: sara@few.vu.nl

    2009-09-25

    The Darboux-dressing transformations developed in Degasperis and Lombardo (2007 J. Phys. A: Math. Theor. 40 961-77) are here applied to construct soliton solutions for a class of boomeronic-type equations. The vacuum (i.e. vanishing) solution and the generic plane wave solution are both dressed to yield one-soliton solutions. The formulae are specialized to the particularly interesting case of the resonant interaction of three waves, a well-known model which is of boomeronic type. For this equation a novel solution which describes three locked dark pulses (simulton) is introduced.

  13. Entropy Stable Wall Boundary Conditions for the Compressible Navier-Stokes Equations

    Science.gov (United States)

    Parsani, Matteo; Carpenter, Mark H.; Nielsen, Eric J.

    2014-01-01

    Non-linear entropy stability and a summation-by-parts framework are used to derive entropy stable wall boundary conditions for the compressible Navier-Stokes equations. A semi-discrete entropy estimate for the entire domain is achieved when the new boundary conditions are coupled with an entropy stable discrete interior operator. The data at the boundary are weakly imposed using a penalty flux approach and a simultaneous-approximation-term penalty technique. Although discontinuous spectral collocation operators are used herein for the purpose of demonstrating their robustness and efficacy, the new boundary conditions are compatible with any diagonal norm summation-by-parts spatial operator, including finite element, finite volume, finite difference, discontinuous Galerkin, and flux reconstruction schemes. The proposed boundary treatment is tested for three-dimensional subsonic and supersonic flows. The numerical computations corroborate the non-linear stability (entropy stability) and accuracy of the boundary conditions.

  14. Entropy Stable Wall Boundary Conditions for the Three-Dimensional Compressible Navier-Stokes Equations

    Science.gov (United States)

    Parsani, Matteo; Carpenter, Mark H.; Nielsen, Eric J.

    2015-01-01

    Non-linear entropy stability and a summation-by-parts framework are used to derive entropy stable wall boundary conditions for the three-dimensional compressible Navier-Stokes equations. A semi-discrete entropy estimate for the entire domain is achieved when the new boundary conditions are coupled with an entropy stable discrete interior operator. The data at the boundary are weakly imposed using a penalty flux approach and a simultaneous-approximation-term penalty technique. Although discontinuous spectral collocation operators on unstructured grids are used herein for the purpose of demonstrating their robustness and efficacy, the new boundary conditions are compatible with any diagonal norm summation-by-parts spatial operator, including finite element, finite difference, finite volume, discontinuous Galerkin, and flux reconstruction/correction procedure via reconstruction schemes. The proposed boundary treatment is tested for three-dimensional subsonic and supersonic flows. The numerical computations corroborate the non-linear stability (entropy stability) and accuracy of the boundary conditions.

  15. A device adaptive inflow boundary condition for Wigner equations of quantum transport

    Energy Technology Data Exchange (ETDEWEB)

    Jiang, Haiyan [Department of Applied Mathematics, Beijing Institute of Technology, Beijing 100081 (China); Lu, Tiao [HEDPS and CAPT, LMAM and School of Mathematical Sciences, Peking University, Beijing 100871 (China); Cai, Wei, E-mail: wcai@uncc.edu [Department of Mathematics and Statistics, University of North Carolina at Charlotte, Charlotte, NC 28223-0001 (United States)

    2014-02-01

    In this paper, an improved inflow boundary condition is proposed for Wigner equations in simulating a resonant tunneling diode (RTD), which takes into consideration the band structure of the device. The original Frensley inflow boundary condition prescribes the Wigner distribution function at the device boundary to be the semi-classical Fermi–Dirac distribution for free electrons in the device contacts without considering the effect of the quantum interaction inside the quantum device. The proposed device adaptive inflow boundary condition includes this effect by assigning the Wigner distribution to the value obtained from the Wigner transform of wave functions inside the device at zero external bias voltage, thus including the dominant effect on the electron distribution in the contacts due to the device internal band energy profile. Numerical results on computing the electron density inside the RTD under various incident waves and non-zero bias conditions show much improvement by the new boundary condition over the traditional Frensley inflow boundary condition.

  16. Solvability of Nonlinear Langevin Equation Involving Two Fractional Orders with Dirichlet Boundary Conditions

    Directory of Open Access Journals (Sweden)

    Bashir Ahmad

    2010-01-01

    Full Text Available We study a Dirichlet boundary value problem for Langevin equation involving two fractional orders. Langevin equation has been widely used to describe the evolution of physical phenomena in fluctuating environments. However, ordinary Langevin equation does not provide the correct description of the dynamics for systems in complex media. In order to overcome this problem and describe dynamical processes in a fractal medium, numerous generalizations of Langevin equation have been proposed. One such generalization replaces the ordinary derivative by a fractional derivative in the Langevin equation. This gives rise to the fractional Langevin equation with a single index. Recently, a new type of Langevin equation with two different fractional orders has been introduced which provides a more flexible model for fractal processes as compared with the usual one characterized by a single index. The contraction mapping principle and Krasnoselskii's fixed point theorem are applied to prove the existence of solutions of the problem in a Banach space.

  17. Dynamics and local boundary properties of the dawn-side magnetopause under conditions observed by Equator-S

    Directory of Open Access Journals (Sweden)

    M. W. Dunlop

    Full Text Available Magnetic field measurements, taken by the magnetometer experiment (MAM on board the German Equator-S spacecraft, have been used to identify and categorise 131 crossings of the dawn-side magnetopause at low latitude, providing unusual, long duration coverage of the adjacent magnetospheric regions and near magnetosheath. The crossings occurred on 31 orbits, providing unbiased coverage over the full range of local magnetic shear from 06:00 to 10:40 LT. Apogee extent places the spacecraft in conditions associated with intermediate, rather than low, solar wind dynamic pressure, as it processes into the flank region. The apogee of the spacecraft remains close to the magnetopause for mean solar wind pressure. The occurrence of the magnetopause encounters are summarised and are found to compare well with predicted boundary location, where solar wind conditions are known. Most scale with solar wind pressure. Magnetopause shape is also documented and we find that the magnetopause orientation is consistently sunward of a model boundary and is not accounted for by IMF or local magnetic shear conditions. A number of well-established crossings, particularly those at high magnetic shear, or exhibiting unusually high-pressure states, were observed and have been analysed for their boundary characteristics and some details of their boundary and near magnetosheath properties are discussed. Of particular note are the occurrence of mirror-like signatures in the adjacent magnetosheath during a significant fraction of the encounters and a high number of multiple crossings over a long time period. The latter is facilitated by the spacecraft orbit which is designed to remain in the near magnetosheath for average solar wind pressure. For most encounters, a well-ordered, tangential (draped magnetosheath field is observed and there is little evidence of large deviations in local boundary orientations. Two passes corresponding to close conjunctions of the Geotail spacecraft

  18. Solutions of Laplace's equation with simple boundary conditions, and their applications for capacitors with multiple symmetries

    CERN Document Server

    Morales, Mayckol; Herrera, William J

    2015-01-01

    We find solutions of Laplace's equation with specific boundary conditions (in which such solutions take either the value zero or unity in each surface) using a generic curvilinear system of coordinates. Such purely geometrical solutions (that we shall call Basic Harmonic Functions BHF's) are utilized to obtain a more general class of solutions for Laplace's equation, in which the functions take arbitrary constant values on the boundaries. On the other hand, the BHF's are also used to obtain the capacitance of many electrostatic configurations of conductors. This method of finding solutions of Laplace's equation and capacitances with multiple symmetries is particularly simple, owing to the fact that the method of separation of variables becomes much simpler under the boundary conditions that lead to the BHF's. Examples of application in complex symmetries are given. Then, configurations of succesive embedding of conductors are also examined. In addition, expressions for electric fields between two conductors a...

  19. An energy absorbing far-field boundary condition for the elastic wave equation

    Energy Technology Data Exchange (ETDEWEB)

    Petersson, N A; Sjogreen, B

    2008-07-15

    The authors present an energy absorbing non-reflecting boundary condition of Clayton-Engquist type for the elastic wave equation together with a discretization which is stable for any ratio of compressional to shear wave speed. They prove stability for a second order accurate finite-difference discretization of the elastic wave equation in three space dimensions together with a discretization of the proposed non-reflecting boundary condition. The stability proof is based on a discrete energy estimate and is valid for heterogeneous materials. The proof includes all six boundaries of the computational domain where special discretizations are needed at the edges and corners. The stability proof holds also when a free surface boundary condition is imposed on some sides of the computational domain.

  20. Numerical methods for solving some boundary problems associated to the Schroedinger equation

    International Nuclear Information System (INIS)

    Full text: The Schroedinger equation is the basic equation of nonrelativistic quantum mechanics. From mathematical point of view this is a partial differential equation of parabolic type, containing a temporal variable and one or more space variables. By separating the temporal part one arrives to the stationary equation, which is a partial differential equation of elliptic type. In some conditions, this equation can be transformed in a system of ordinary differential equations with boundary conditions specific to the problem. The natural domain of definition is unbound. For the numerical solution one should restrict to a finite domain. The adequate treatment of boundary conditions on this computational domain is decisive for the accuracy of the results. Concerning the stationary equation, of interest is the radial equation which is obtained by separating the angular part from the radial one. The interval of the independent variable r (which from physical point of view is [0, ∞) is truncated to a segment whose first end is in the vicinity of the origin (rmin), and the second in the asymptotic region (ras). The influence of the two excluded subintervals [0, rmin] and [ras, ∞), should be carefully taken into account. Mathematically, the essential peculiarity is that the first interval is a singular zone, while the second is an asymptotic zone and thus the numerical solution in these regions should be constructed by considering these features. The thesis presents a systematic approach to the numerical solution of the Schroedinger equation in the two extreme zones. Our aim was to obtain stable and rapid algorithms with controlled accuracy. The used procedures have been perturbative, exponential fitting, convergence acceleration and summation methods. The case of systems of coupled equations has been also discussed. In the case of the time dependent equation, the limitation to a finite spatial domain leads to reflexions which affect the propagated wavefunction and

  1. Integrability Estimates for Gaussian Rough Differential Equations

    CERN Document Server

    Cass, Thomas; Lyons, Terry

    2011-01-01

    We derive explicit tail-estimates for the Jacobian of the solution flow of stochastic differential equations driven by Gaussian rough paths. In particular, we deduce that the Jacobian has finite moments of all order for a wide class of Gaussian process including fractional Brownian motion with Hurst parameter H>1/4. We remark on the relevance of such estimates to a number of significant open problems.

  2. On Neumann boundary value problems for some quasilinear elliptic equations

    Directory of Open Access Journals (Sweden)

    Paul A. Binding

    1997-01-01

    Full Text Available function $a(x$ on the existence of positive solutions to the problem $$left{ eqalign{ -{ m div},(|abla u|^{p-2}abla u&= lambda a(x|u|^{p-2}u+b(x|u|^{gamma-2}u, quad xinOmega, cr x{partial u overpartial n}&=0, quad xinpartialOmega,,} ight. $$ where $Omega$ is a smooth bounded domain in $R^n$, $b$ changes sign, $1integral condition, then there exists some $lambda^*$ such that $lambda^* int_Omega a(x, dx<0$ and, for $lambda$ strictly between 0 and $lambda^*$, the problem has a positive solution. (ii if $int_Omega a(x, dx=0$, then the problem has a positive solution for small $lambda$ provided that $int_Omega b(x,dx<0$.

  3. 跨共振的周期-积分边值问题%Periodic-Integral Boundary Value Problems across Resonance

    Institute of Scientific and Technical Information of China (English)

    宋新; 杨雪

    2011-01-01

    研究二阶微分方程周期-积分边值问题,应用最优控制理论给出了跨多个共振情形下的二阶微分方程周期-积分边值问题唯一可解的最优条件.%The periodic-integral boundary value problems for second order differential equations were considered. On the basis of optimal control theory method, we gave an optimal condition of the unique solvability to the periodic-integral boundary value problems for second order differential equations across multiple resonance.

  4. Explicit construction of solutions to the Burgers equation with discontinuous initial-boundary conditions

    OpenAIRE

    Désilles, Anya; Frankowska, Hélène

    2013-01-01

    International audience Abstract. A solution of the initial-boundary value problem on the strip (0,1) × [0, 1] for scalar conservation laws with strictly convex flux can be obtained by considering gradients of the unique solution V to an associated Hamilton-Jacobi equation (with appropriately defined initial and boundary conditions). It was shown in Frankowska (2010) that V can be expressed as the minimum of three value functions arising in calculus of variations problems that, in turn, can...

  5. Nonlinear boundary value problems for first order impulsive integro-differential equations

    Directory of Open Access Journals (Sweden)

    Xinzhi Liu

    1989-01-01

    Full Text Available In this paper, we investigate a class of first order impulsive integro-differential equations subject to certain nonlinear boundary conditions and prove, with the help of upper and lower solutions, that the problem has a solution lying between the upper and lower solutions. We also develop monotone iterative technique and show the existence of multiple solutions of a class of periodic boundary value problems.

  6. On Initial-Boundary Value Problem of Stochastic Heat Equation in a Lipschitz Cylinder

    CERN Document Server

    Chang, Tongkeun; Yang, Minsuk

    2011-01-01

    We consider the initial boundary value problem of non-homogeneous stochastic heat equation. The derivative of the solution with respect to time receives heavy random perturbation. The space boundary is Lipschitz and we impose non-zero cylinder condition. We prove a regularity result after finding suitable spaces for the solution and the pre-assigned datum in the problem. The tools from potential theory, harmonic analysis and probability are used. Some Lemmas are as important as the main Theorem.

  7. Exact Synchronization for a Coupled System of Wave Equations with Dirichlet Boundary Controls

    Institute of Scientific and Technical Information of China (English)

    Tatsien LI; Bopeng RAO

    2013-01-01

    In this paper,the exact synchronization for a coupled system of wave equations with Dirichlet boundary controls and some related concepts are introduced.By means of the exact null controllability of a reduced coupled system,under certain conditions of compatibility,the exact synchronization,the exact synchronization by groups,and the exact null controllability and synchronization by groups are all realized by suitable boundary controls.

  8. Quenching rates for heat equations with coupled singular nonlinear boundary flux

    Institute of Scientific and Technical Information of China (English)

    2008-01-01

    We study finite time quenching for heat equations coupled via singular nonlinear bound-ary flux. A criterion is proposed to identify the simultaneous and non-simultaneous quenchings. In particular, three kinds of simultaneous quenching rates are obtained for different nonlinear exponent re-gions and appropriate initial data. This extends an original work by Pablo, Quir′os and Rossi for a heat system with coupled inner absorption terms subject to homogeneous Neumann boundary conditions.

  9. Asymptotic stability and blow up for a semilinear damped wave equation with dynamic boundary conditions

    KAUST Repository

    Gerbi, Stéphane

    2011-12-01

    In this paper we consider a multi-dimensional wave equation with dynamic boundary conditions, related to the KelvinVoigt damping. Global existence and asymptotic stability of solutions starting in a stable set are proved. Blow up for solutions of the problem with linear dynamic boundary conditions with initial data in the unstable set is also obtained. © 2011 Elsevier Ltd. All rights reserved.

  10. Regularized boundary integral representations for dislocations and cracks in smart media

    International Nuclear Information System (INIS)

    This paper presents a complete set of singularity-reduced boundary integral relations for isolated discontinuities embedded in three-dimensional infinite media. The development is carried out within a broad context that allows the treatment of a well-known class of smart media such as linear piezoelectric, linear piezomagnetic and linear piezoelectromagnetic materials. In addition, resulting boundary integral representations are applicable to general discontinuities of arbitrary geometry and possessing a general jump distribution. The latter aspect allows the treatment of two special kinds of discontinuities: dislocations and cracks. The most attractive feature of the current development is that all integral relations for field quantities such as state variables and their gradients, the body flux, and the generalized interaction energy produced by dislocations are expressed only in terms of line integrals over the dislocation loops and, for cracks, the key governing boundary integral equation is established in a symmetric weak form and contains only weakly singular kernels of O(1/r). Results for the former case are fundamental and useful in the context of dislocation mechanics and modeling while the resulting weakly singular, weak form integral equation constitutes a basis for the development of a well-known numerical technique, called a symmetric Galerkin boundary element method (SGBEM), for analysis of cracked bodies. The weakly singular nature of such an integral equation allows low order interpolations to be used in the numerical approximation. The key ingredient for achieving such development of integral representations is the use of certain special decompositions in the derivative-transferring process via Stokes's theorem. Existence of such decompositions is ensured by a careful consideration of the singularity nature of the kernels, and a particular solution of the weakly singular functions involved is obtained by solving a system of partial differential

  11. Integral conditions in the theory of the Beltrami equations

    CERN Document Server

    Ryazanov, V; Yakubov, E

    2010-01-01

    It is shown that many recent and new results on the existence of ACL homeomorphic solutions for the degenerate Beltrami equations with integral constraints follow from our extension of the well--known Lehto existence theorem.

  12. Nonlinear wave equations, covariant exterior derivative, Painleve test, and integrability

    International Nuclear Information System (INIS)

    A class of nonlinear wave equations is derived which can be written as the vanishing of a covariant exterior derivative. The Painleve test is performed and the connection with integrability is discussed

  13. Existence and uniqueness theorem for a Hammerstein nonlinear integral equation

    OpenAIRE

    A. Kh. Khachatryan; Kh. A. Khachatryan

    2011-01-01

    The existence of a solution, as well as some properties of the obtained solution for a Hammerstein type nonlinear integral equation have been investigated. For a certain class of functions the uniqueness theorem has also been proved.

  14. LETTER TO THE EDITOR: Deconstructing an integrable lattice equation

    Science.gov (United States)

    Ramani, A.; Joshi, N.; Grammaticos, B.; Tamizhmani, T.

    2006-02-01

    We show that an integrable lattice equation, obtained by J Hietarinta using the 'consistency around a cube' method without the tetrahedron assumption, is indeed solvable by linearization. We also present its nonautonomous extension.

  15. Generalised Fourier integral operator methods for hyperbolic equations with singularities

    OpenAIRE

    Garetto, Claudia; Oberguggenberger, Michael

    2011-01-01

    This article addresses linear hyperbolic partial differential equations and pseudodifferential equations with strongly singular coefficients and data, modelled as members of algebras of generalised functions. We employ the recently developed theory of generalised Fourier integral operators to construct parametrices for the solutions and to describe propagation of singularities in this setting. As required tools, the construction of generalised solutions to eikonal and transport equations is g...

  16. Exact Solutions and Conservation Laws for a New Integrable Equation

    International Nuclear Information System (INIS)

    In this work we study a generalization of an integrable equation proposed by Qiao and Liu from the point of view of the theory of symmetry reductions in partial differential equations. Among the solutions we obtain a travelling wave with decaying velocity and a smooth soliton solution. We determine the subclass of these equations which are quasi-self-adjoint and we get a nontrivial conservation law.

  17. On the analytical solution of the SN equation in a rectangle assuming an exponential exiting angular flux boundary

    International Nuclear Information System (INIS)

    In this work, we report an analytical solution for the set of SN equations for the angular flux, in a rectangle, using the double Laplace transform technique. Its main idea comprehends the steps: application of the Laplace transform in one space variable, solution of the resulting equation by the LTSN method and reconstruction of the double Laplace transformed angular flux using the inversion theorem of the Laplace transform. We must emphasize that we perform the Laplace inversion by the LTSN method in the x direction, meanwhile we evaluate the inversion in the y direction performing the calculation of the corresponding line integral solution by the Stefest method. We have also to figure out that the application of Laplace transform to this type of boundary value problem introduces additional unknown functions associated to the partial derivatives of the angular flux at boundary. Based on the good results attained by the nodal LTSN method, we assume that the angular flux at boundary is also approximated by an exponential function. By analytical we mean that no approximation is done along the solution derivation except for the exponential hypothesis for the exiting angular flux at boundary. For sake of completeness, we report numerical comparisons of the obtained results against the ones of the literature. (author)

  18. Integral conditions for nonoscillation of second order nonlinear differential equations

    Czech Academy of Sciences Publication Activity Database

    Cecchi, M.; Došlá, Z.; Marini, M.; Vrkoč, Ivo

    2006-01-01

    Roč. 64, č. 6 (2006), s. 1278-1289. ISSN 0362-546X R&D Projects: GA AV ČR(CZ) IAA1163401 Institutional research plan: CEZ:AV0Z10190503 Keywords : change of integration * half-linear differential equation * nonlinear differential equation Subject RIV: BA - General Mathematics Impact factor: 0.677, year: 2006

  19. Radial symmetry and monotonicity for an integral equation

    Science.gov (United States)

    Ma, Li; Chen, Dezhong

    2008-06-01

    In this paper we study radial symmetry and monotonicity of positive solutions of an integral equation arising from some higher-order semilinear elliptic equations in the whole space . Instead of the usual method of moving planes, we use a new Hardy-Littlewood-Sobolev (HLS) type inequality for the Bessel potentials to establish the radial symmetry and monotonicity results.

  20. Singular perturbation for the weakly nonlinear reaction diffusion equation with boundary perturbation

    Institute of Scientific and Technical Information of China (English)

    MO Jia-qi

    2008-01-01

    In this paper, a class of nonlinear singularly perturbed initial boundary value problems for reaction diffusion equations with boundary perturbation are considered under suitable conditions. Firstly, by dint of the regular perturbation method, the outer solution of the original problem is obtained. Secondly, by using the stretched variable and the expansion theory of power series the initial layer of the solution is constructed. And then, by using the theory of differential inequalities, the asymptotic behavior of the solution for the initial boundary value problems is studied. Finally, using some relational inequalities the existence and uniqueness of solution for the original problem and the uniformly valid asymptotic estimation are discussed.

  1. THE NONLINEAR NONLOCAL SINGULARLY PERTURBED PROBLEMS FOR REACTION DIFFUSION EQUATIONS WITH A BOUNDARY PERTURBATION

    Institute of Scientific and Technical Information of China (English)

    Jingsun Yao; Jiaqi Mo

    2005-01-01

    The nonlinear nonlocal singularly perturbed initial boundary value problems for reaction diffusion equations with a boundary perturbation is considered. Under suitable conditions, the outer solution of the original problem is obtained. Using the stretched variable, the composing expansion method and the expanding theory of power series the initial layer is constructed. And then using the theory of differential inequalities the asymptotic behavior of solution for the initial boundary value problems is studied. Finally the existence and uniqueness of solution for the original problem and the uniformly valid asymptotic estimation are discussed.

  2. A CLASS OF NONLINEAR NONLOCAL SINGULARLY PERTURBED PROBLEMS FOR REACTION DIFFUSION EQUATIONS WITH BOUNDARY PERTURBATION

    Institute of Scientific and Technical Information of China (English)

    MO Jia-qi; WANG Hui; LIN Wan-tao

    2005-01-01

    A class of nonlinear nonlocal for singularly perturbed Robin initial boundary value problems for reaction diffusion equations with boundary perturbation is considered. Under suitable conditions, first, the outer solution of the original problem was obtained. Secondly, using the stretched variable, the composing expansion method and the expanding theory of power series the initial layer was constructed. Finally, using the theory of differential inequalities the asymptotic behavior of solution for the initial boundary value problems was studied, and educing some relational inequalities the existence and uniqueness of solution for the original problem and the uniformly valid asymptotic estimation were discussed.

  3. The Focusing NLS Equation on the Half-Line with Periodic Boundary Conditions

    OpenAIRE

    Kamvissis, S.; Fokas, A. S.

    2007-01-01

    We consider the Dirichlet problem for the focusing NLS equation on the half-line, with given Schwartz initial data and boundary data $q(0,t)$ equal to an exponentially decaying perturbation $u(t)$ of the periodic boundary data $ a e^{2i\\omega t + i \\epsilon}$ at $x=0.$ It is known from PDE theory that this is a well-posed problem (for fixed initial data and fixed $u$). On the other hand, the associated inverse scattering transform formalism involves the Neumann boundary value for $x=0$. Thus ...

  4. Second-order boundary estimates for solutions to singular elliptic equations in borderline cases

    Directory of Open Access Journals (Sweden)

    Claudia Anedda

    2011-04-01

    Full Text Available Let $Omegasubset R^N$ be a bounded smooth domain. We investigate the effect of the mean curvature of the boundary $partialOmega$ on the behaviour of the solution to the homogeneous Dirichlet boundary value problem for the equation $Delta u+f(u=0$. Under appropriate growth conditions on $f(t$ as $t$ approaches zero, we find asymptotic expansions up to the second order of the solution in terms of the distance from $x$ to the boundary $partialOmega$.

  5. Integrable achiral D5-brane reflections and asymptotic Bethe equations

    CERN Document Server

    Correa, Diego H; Young, Charles A S

    2011-01-01

    We study the reflection of magnons from a D5-brane in the framework of the AdS/CFT correspondence. We consider two possible orientations of the D5-brane with respect to the reference vacuum state, namely vacuum states aligned along "vertical" and "horizontal" directions. We show that the reflections are of the achiral type. We also show that the reflection matrices satisfy the boundary Yang-Baxter equations for both orientations. In the horizontal case the reflection matrix can be interpreted in terms of a bulk S-matrix, S(p, -p), and factorizability of boundary scattering therefore follows from that of bulk scattering. Finally, we solve the nested coordinate Bethe ansatz for the system in the vertical case to find the Bethe equations. In the horizontal case, the Bethe equations are of the same form as those for the closed string.

  6. Existence and asymptotic behavior of the wave equation with dynamic boundary conditions

    KAUST Repository

    Graber, Philip Jameson

    2012-03-07

    The goal of this work is to study a model of the strongly damped wave equation with dynamic boundary conditions and nonlinear boundary/interior sources and nonlinear boundary/interior damping. First, applying the nonlinear semigroup theory, we show the existence and uniqueness of local in time solutions. In addition, we show that in the strongly damped case solutions gain additional regularity for positive times t>0. Second, we show that under some restrictions on the initial data and if the interior source dominates the interior damping term and if the boundary source dominates the boundary damping, then the solution grows as an exponential function. Moreover, in the absence of the strong damping term, we prove that the solution ceases to exists and blows up in finite time. © 2012 Springer Science+Business Media, LLC.

  7. Imposing Neumann boundary condition on cosmological perturbation equations and trajectories of particles

    Science.gov (United States)

    Shenavar, Hossein

    2016-03-01

    We impose Neumann boundary condition to solve cosmological perturbation equations and we derive a modified Friedmann equation and a new lensing equation. To check the new lensing equation and the value of Neumann constant, a sample that contains ten strong lensing systems is surveyed. Except for one lens, masses of the other lenses are found to be within the constrains of the observational data. Furthermore, we argue that by using the concept of geometrodynamic clocks it is possible to modify the equation of motion of massive particles too. Also, a sample that includes 101 HSB and LSB galaxies is used to re-estimate the value of the Neumann constant and we found that this value is consistent with the prior evaluation from Friedmann and lensing equations. Finally, the growth of structure is studied by a Newtonian approach which resulted in a more rapid rate of the structure formation in matter dominated era.

  8. Navier-Stokes equations in three-dimensional thin domains with various boundary conditions

    OpenAIRE

    Temam, R.; M. Ziane

    1996-01-01

    In this work we develop methods for studying the Navier-Stokes equations in thin domains. We consider various boundary conditions and establish the global existence of strong solutions when the initial data belong to "large sets." Our work was inspired by the recent interesting results of G. Raugel and G. Sell [22, 23, 24] which, in the periodic case, give global existence for smooth solutions of the 3D Navier-Stokes equations in thin domains for large sets of initial condition...

  9. Asymptotic Behavior of the 3D Compressible Euler Equations with Nonlinear Damping and Slip Boundary Condition

    OpenAIRE

    Huimin Yu

    2012-01-01

    The asymptotic behavior (as well as the global existence) of classical solutions to the 3D compressible Euler equations are considered. For polytropic perfect gas $(P(\\rho )={P}_{0}{\\rho }^{\\gamma })$ , time asymptotically, it has been proved by Pan and Zhao (2009) that linear damping and slip boundary effect make the density satisfying the porous medium equation and the momentum obeying the classical Darcy's law. In this paper, we use a more general method and extend this resu...

  10. Reconstruction of the boundary condition for the heat conduction equation of fractional order

    OpenAIRE

    Brociek Rafal; Slota Damian

    2015-01-01

    This paper describes reconstruction of the heat transfer coefficient occurring in the boundary condition of the third kind for the time fractional heat conduction equation. Fractional derivative with respect to time, occurring in considered equation, is defined as the Caputo derivative. Additional information for the considered inverse problem is given by the temperature measurements at selected points of the domain. The direct problem is solved by using th...

  11. Strang-type preconditioners for solving fractional diffusion equations by boundary value methods

    OpenAIRE

    Gu, Xian-Ming; Huang, Ting-Zhu; Zhao, Xi-Le; Li, Hou-Biao; Li, Liang

    2013-01-01

    The finite difference scheme with the shifted Gr\\"{u}nwarld formula is employed to semi-discrete the fractional diffusion equations. This spatial discretization can reduce to the large system of ordinary differential equations (ODEs) with initial values. Recently, boundary value method (BVM) was developed as a popular algorithm for solving large systems of ODEs. This method requires the solutions of one or more nonsymmetric, large and sparse linear systems. In this paper, the GMRES method wit...

  12. A boundary value approach for solving three-dimensional elliptic and hyperbolic partial differential equations

    OpenAIRE

    Biala, T. A.; Jator, S. N.

    2015-01-01

    In this article, the boundary value method is applied to solve three dimensional elliptic and hyperbolic partial differential equations. The partial derivatives with respect to two of the spatial variables (y, z) are discretized using finite difference approximations to obtain a large system of ordinary differential equations (ODEs) in the third spatial variable (x). Using interpolation and collocation techniques, a continuous scheme is developed and used to obtain discrete methods which are ...

  13. An efficient computer based wavelets approximation method to solve Fuzzy boundary value differential equations

    Science.gov (United States)

    Alam Khan, Najeeb; Razzaq, Oyoon Abdul

    2016-03-01

    In the present work a wavelets approximation method is employed to solve fuzzy boundary value differential equations (FBVDEs). Essentially, a truncated Legendre wavelets series together with the Legendre wavelets operational matrix of derivative are utilized to convert FB- VDE into a simple computational problem by reducing it into a system of fuzzy algebraic linear equations. The capability of scheme is investigated on second order FB- VDE considered under generalized H-differentiability. Solutions are represented graphically showing competency and accuracy of this method.

  14. Localization of the eigenvalues of linear integral equations with applications to linear ordinary differential equations.

    Science.gov (United States)

    Sloss, J. M.; Kranzler, S. K.

    1972-01-01

    The equivalence of a considered integral equation form with an infinite system of linear equations is proved, and the localization of the eigenvalues of the infinite system is expressed. Error estimates are derived, and the problems of finding upper bounds and lower bounds for the eigenvalues are solved simultaneously.

  15. Smoothed Particle Hydrodynamics Continuous Boundary Force method for Navier-Stokes equations subject to a Robin boundary condition

    Science.gov (United States)

    Pan, Wenxiao; Bao, Jie; Tartakovsky, Alexandre

    2013-11-01

    A Continuous Boundary Force (CBF) method was developed for implementing Robin (Navier) boundary condition (BC) that can describe no-slip or slip conditions (slip length from zero to infinity) at the fluid-solid interface. In the CBF method the Robin BC is replaced by a homogeneous Neumann BC and an additional volumetric source term in the governing momentum equation. The formulation is derived based on an approximation of the sharp boundary with a diffuse interface of finite thickness, across which the BC is reformulated by means of a smoothed characteristic function. The CBF method is easy to be implemented in Lagrangian particle-based methods. We first implemented it in smoothed particle hydrodynamics (SPH) to solve numerically the Navier-Stokes equations subject to spatial-independent or dependent Robin BC in two and three dimensions. The numerical accuracy and convergence is examined through comparisons with the corresponding finite difference or finite element solutions. The CBF method is further implemented in smoothed dissipative particle dynamics (SDPD), a mesoscale scheme, for modeling slip flows commonly existent in micro/nano channels and microfluidic devices. The authors acknowledge the funding support by the ASCR Program of the Office of Science, U.S. Department of Energy.

  16. Integral equation theory for nematic fluids

    Directory of Open Access Journals (Sweden)

    M.F. Holovko

    2010-01-01

    Full Text Available The traditional formalism in liquid state theory based on the calculation of the pair distribution function is generalized and reviewed for nematic fluids. The considered approach is based on the solution of orientationally inhomogeneous Ornstein-Zernike equation in combination with the Triezenberg-Zwanzig-Lovett-Mou-Buff-Wertheim equation. It is shown that such an approach correctly describes the behavior of correlation functions of anisotropic fluids connected with the presence of Goldstone modes in the ordered phase in the zero-field limit. We focus on the discussions of analytical results obtained in collaboration with T.G. Sokolovska in the framework of the mean spherical approximation for Maier-Saupe nematogenic model. The phase behavior of this model is presented. It is found that in the nematic state the harmonics of the pair distribution function connected with the correlations of the director transverse fluctuations become long-range in the zero-field limit. It is shown that such a behavior of distribution function of nematic fluid leads to dipole-like and quadrupole-like long-range asymptotes for effective interaction between colloids solved in nematic fluids, predicted before by phenomenological theories.

  17. Open two-species exclusion processes with integrable boundaries

    International Nuclear Information System (INIS)

    We give a complete classification of integrable Markovian boundary conditions for the asymmetric simple exclusion process with two species (or classes) of particles. Some of these boundary conditions lead to non-vanishing particle currents for each species. We explain how the stationary state of all these models can be expressed in a matrix product form, starting from two key components, the Zamolodchikov–Faddeev and Ghoshal–Zamolodchikov relations. This statement is illustrated by studying in detail a specific example, for which the matrix ansatz (involving nine generators) is explicitly constructed and physical observables (such as currents, densities) calculated. (paper)

  18. Dirac equation for particles with arbitrary half-integral spin

    Science.gov (United States)

    Guseinov, I. I.

    2011-11-01

    The sets of ? -component irreducible and Clifford algebraic Hermitian and unitary matrices through the two-component Pauli matrices are suggested, where s = 1/2, 3/2, 5/2, … . Using these matrix sets, the eigenvalues of which are ? , the ? -component generalized Dirac equation for a description of arbitrary half-integral spin particles is constructed. In accordance with the correspondence principle, the generalized Dirac equation suggested arises from the condition of relativistic invariance. This equation is reduced to the sets of two-component matrix equations the number of which is equal to ? . The new relativistic invariant equation of motion leads to an equation of continuity with a positive-definite probability density and also to the Klein-Gordon equation. This relativistic equation is causal in the presence of an external electromagnetic field interaction. It is shown that, in the case of nonrelativistic limit, the relativistic equation presented is reduced to the Pauli equation describing the motion of half-integral spin particle in the electromagnetic field.

  19. Analytical solutions to the fractional advection-diffusion equation with time-dependent pulses on the boundary

    Science.gov (United States)

    Rubbab, Qammar; Mirza, Itrat Abbas; Qureshi, M. Zubair Akbar

    2016-07-01

    The time-fractional advection-diffusion equation with Caputo-Fabrizio fractional derivatives (fractional derivatives without singular kernel) is considered under the time-dependent emissions on the boundary and the first order chemical reaction. The non-dimensional problem is formulated by using suitable dimensionless variables and the fundamental solutions to the Dirichlet problem for the fractional advection-diffusion equation are determined using the integral transforms technique. The fundamental solutions for the ordinary advection-diffusion equation, fractional and ordinary diffusion equation are obtained as limiting cases of the previous model. Using Duhamel's principle, the analytical solutions to the Dirichlet problem with time-dependent boundary pulses have been obtained. The influence of the fractional parameter and of the drift parameter on the solute concentration in various spatial positions was analyzed by numerical calculations. It is found that the variation of the fractional parameter has a significant effect on the solute concentration, namely, the memory effects lead to the retardation of the mass transport.

  20. Fringe integral equation method for a truncated grounded dielectric slab

    DEFF Research Database (Denmark)

    Jørgensen, Erik; Maci, S.; Toccafondi, A.

    2001-01-01

    The problem of scattering by a semi-infinite grounded dielectric slab illuminated by an arbitrary incident TMz polarized electric field is studied by solving a new set of “fringe” integral equations (F-IEs), whose functional unknowns are physically associated to the wave diffraction processes...... occurring at the truncation. The F-IEs are obtained by subtracting from the surface/surface integral equations pertinent to the truncated slab, an auxiliary set of equations obtained for the canonical problem of an infinite grounded slab illuminated by the same source. The F-IEs are solved by the method of...

  1. A Hierarchy of Integrable Lattice Soliton Equations and New Integrable Symplectic Map

    International Nuclear Information System (INIS)

    Starting from a discrete spectral problem, a hierarchy of integrable lattice soliton equations is derived. It is shown that the hierarchy is completely integrable in the Liouville sense and possesses discrete bi-Hamiltonian structure. A new integrable symplectic map and finite-dimensional integrable systems are given by nonlinearization method. The binary Bargmann constraint gives rise to a Baecklund transformation for the resulting integrable lattice equations. At last, conservation laws of the hierarchy are presented.

  2. A Hierarchy of Integrable Lattice Soliton Equations and New Integrable Symplectic Map

    Institute of Scientific and Technical Information of China (English)

    2006-01-01

    Starting from a discrete spectral problem, a hierarchy of integrable lattice soliton equations is derived. It is shown that the hierarchy is completely integrable in the Liouville sense and possesses discrete bi-Hamiltonian structure.A new integrable symplectic map and finite-dimensional integrable systems are given by nonlinearization method. The binary Bargmann constraint gives rise to a B(a)cklund transformation for the resulting integrable lattice equations. At last, conservation laws of the hierarchy are presented.

  3. Nonlinear integral equations for shape reconstruction in the inverse interior scattering problem

    International Nuclear Information System (INIS)

    In this paper, we consider the inverse scattering problem of recovering the shape of a perfectly conducting cavity from one source and several measurements placed on a curve inside the cavity. Under restrictive assumptions on the size of the cavity, a uniqueness theorem for finitely many excitations is given. Based on a system of nonlinear and ill-posed integral equations for the unknown boundary, we apply a regularized Newton iterative approach to find the boundary. We present the mathematical foundation of the method and give several numerical examples to show the viability of the method

  4. Vanishing viscosity limits for the degenerate lake equations with Navier boundary conditions

    CERN Document Server

    Jiu, Quansen; Wu, Jiahong

    2011-01-01

    The paper is concerned with the vanishing viscosity limit of the two-dimensional degenerate viscous lake equations when the Navier slip conditions are prescribed on the impermeable boundary of a simply connected bounded regular domain. When the initial vorticity is in the Lebesgue space $L^q$ with $2equations possess a unique global solution and the solution converges to a corresponding weak solution of the inviscid lake equations. In the special case when the vorticity is in $L^\\infty$, an explicit convergence rate is obtained.

  5. Positive Solutions of Singular Boundary Value Problem of Negative Exponent Emden–Fowler Equation

    Indian Academy of Sciences (India)

    Yuxia Wang; Xiyu Liu

    2003-05-01

    This paper investigates the existence of positive solutions of a singular boundary value problem with negative exponent similar to standard Emden–Fowler equation. A necessary and sufficient condition for the existence of [0, 1] positive solutions as well as 1[0, 1] positive solutions is given by means of the method of lower and upper solutions with the Schauder fixed point theorem.

  6. Multiple Positive Solutions of Boundary Value Problems for Systems of Nonlinear Third-Order Differential Equations

    Institute of Scientific and Technical Information of China (English)

    Yaohong LI; Xiaoyan ZHANG

    2013-01-01

    In this paper,we consider boundary value problems for systems of nonlinear thirdorder differential equations.By applying the fixed point theorems of cone expansion and compression of norm type and Leggett-Williams fixed point theorem,the existence of multiple positive solutions is obtained.As application,we give some examples to demonstrate our results.

  7. Quasilinearization for the periodic boundary value problem for systems of impulsive differential equations

    Directory of Open Access Journals (Sweden)

    2006-01-01

    Full Text Available The method of generalized quasilinearization for the system of nonlinear impulsive differential equations with periodic boundary conditions is studied. As a byproduct, the result for the system without impulses can be obtained, which is a new result as well.

  8. Non-perturbative solution of free-convective boundary-layer equation by Adomian decomposition method

    Energy Technology Data Exchange (ETDEWEB)

    Kechil, Seripah Awang [Department of Mathematics, Universiti Teknologi MARA, 40450 Shah Alam Selangor (Malaysia); Hashim, Ishak [School of Mathematical Sciences, Universiti Kebangsaan Malaysia, 43600 UKM Bangi Selangor (Malaysia)]. E-mail: ishak_h@ukm.my

    2007-03-19

    A free-convective boundary layer flow modeled by a system of nonlinear ordinary differential equations is considered. The system is solved using the Adomian decomposition method (ADM) which yields an analytic solution in the form of a rapidly convergent infinite series with easily computable terms. The analytical solutions and the pertinent features of the illustrations show the efficiency of the method.

  9. Initial boundary value problems of nonlinear wave equations in an exterior domain

    International Nuclear Information System (INIS)

    In this paper, we investigate the existence and uniqueness of the global solutions to the initial boundary value problems of nonlinear wave equations in an exterior domain. When the space dimension n >= 3, the unique global solution of the above problem is obtained for small initial data, even if the nonlinear term is fully nonlinear and contains the unknown function itself. (author). 10 refs

  10. Weak versus strong no-slip boundary conditions for the Navier-Stokes equations

    OpenAIRE

    Abbas, Qaisar; Nordström, Jan

    2010-01-01

    Original Publication:Qaisar Abbas and Jan Nordström, Weak versus strong no-slip boundary conditions for the Navier-Stokes equations, 2010, Engineering Applications of Computational Fluid Mechanics, (4), 29-38. Copyright: Engineering Applications of Computational Fluid Mechanics

  11. On two-point boundary value problems for second order singular functional differential equations

    Czech Academy of Sciences Publication Activity Database

    Kiguradze, I.; Půža, Bedřich

    2005-01-01

    Roč. 12, 3-4 (2005), s. 271-294. ISSN 0793-1786 Institutional research plan: CEZ:AV0Z10190503 Keywords : second order singular functional differential equation * two-point boundary value problem * solvability Subject RIV: BA - General Mathematics

  12. Positive solutions of a boundary-value problem for second order ordinary differential equations

    Directory of Open Access Journals (Sweden)

    George L. Karakostas

    2000-06-01

    Full Text Available The existence of positive solutions of a two-point boundary value problem for a second order differential equation is investigated. By using indices of convergence of the nonlinearities at zero and at positive infinity, we providea priori upper and lower bounds for the slope of the solutions are also provided.

  13. A coupled boundary element-finite difference solution of the elliptic modified mild slope equation

    DEFF Research Database (Denmark)

    Naserizadeh, R.; Bingham, Harry B.; Noorzad, A.

    2011-01-01

    The modified mild slope equation of [5] is solved using a combination of the boundary element method (BEM) and the finite difference method (FDM). The exterior domain of constant depth and infinite horizontal extent is solved by a BEM using linear or quadratic elements. The interior domain with...

  14. Inertial Manifolds for the 3D Cahn-Hilliard Equations with Periodic Boundary Conditions

    OpenAIRE

    Kostianko, A; Zelik, S.

    2014-01-01

    The existence of an inertial manifold for the 3D Cahn-Hilliard equation with periodic boundary conditions is verified using the proper extension of the so-called spatial averaging principle introduced by G. Sell and J. Mallet-Paret. Moreover, the extra regularity of this manifold is also obtained.

  15. Quasi-linear dynamics in nonlinear Schr\\" odinger equation with periodic boundary conditions

    OpenAIRE

    Erdogan, M. Burak; Zharnitsky, Vadim

    2007-01-01

    It is shown that a large subset of initial data with finite energy ($L^2$ norm)evolves nearly linearly in nonlinear Schr\\" odinger equation with periodic boundary conditions. These new solutions are not perturbations of the known ones such as solitons, semiclassical or weakly linear solutions.

  16. Periodic solutions of the quasilinear equation of forced beam vibrations with homogeneous boundary conditions

    Science.gov (United States)

    Rudakov, I. A.

    2015-10-01

    We prove the existence of a countable family of time-periodic solutions of the quasilinear equation of beam vibrations with homogeneous boundary conditions and time-periodic right-hand side in the case when the non-linear term has power growth.

  17. Fredholm-type theorem for boundary value problems for systems of nonlinear functional differential equations

    Czech Academy of Sciences Publication Activity Database

    Hakl, Robert; Zamora, M.

    -, May 2014 (2014), s. 113. ISSN 1687-2770 Institutional support: RVO:67985840 Keywords : functional-differential equations * boundary value problems * existence of solutions Subject RIV: BA - General Mathematics Impact factor: 1.014, year: 2014 http://www.boundaryvalueproblems.com/content/2014/1/113

  18. THE ANALYTICAL SOLUTION FOR SEDIMENT REACTION AND DIFFUSION EQUATION WITH GENERALIZED INITIAL-BOUNDARY CONDITIONS

    Institute of Scientific and Technical Information of China (English)

    熊岳山; 韦永康

    2001-01-01

    The sediment reaction and diffusion equation with generalized initial and boundary condition is studied. By using Laplace transform and Jordan lemma , an analytical solution is got, which is an extension of analytical solution provided by Cheng Kwokming James ( only diffusion was considered in analytical solution of Cheng ). Some problems arisen in the computation of analytical solution formula are also analysed.

  19. Nonlocal boundary value problem for strongly singular higher-order linear functional-differential equations

    Czech Academy of Sciences Publication Activity Database

    Mukhigulashvili, Sulkhan

    -, č. 33 (2013), s. 1-38 ISSN 1417-3875 Institutional support: RVO:67985840 Keywords : higher order linear differential equation * nonlocal boundary conditions * deviating argument Subject RIV: BA - General Mathematics Impact factor: 0.638, year: 2013 http://www.emis.de/journals/EJQTDE/p2108. pdf

  20. The nonlocal boundary value problems for strongly singular higher-order nonlinear functional-differential equations

    Czech Academy of Sciences Publication Activity Database

    Mukhigulashvili, Sulkhan

    -, č. 35 (2015), s. 23-50. ISSN 1126-8042 Institutional support: RVO:67985840 Keywords : higher order functional differential equations * Dirichlet boundary value problem * strong singularity Subject RIV: BA - General Mathematics http://ijpam.uniud.it/online_issue/201535/03-Mukhigulashvili. pdf

  1. Burgers equation with no-flux boundary conditions and its application for complete fluid separation

    Science.gov (United States)

    Watanabe, Shinya; Matsumoto, Sohei; Higurashi, Tomohiro; Ono, Naoki

    2016-09-01

    Burgers equation in a one-dimensional bounded domain with no-flux boundary conditions at both ends is proven to be exactly solvable. Cole-Hopf transformation converts not only the governing equation to the heat equation with an extra damping but also the nonlinear mixed boundary conditions to Dirichlet boundary conditions. The average of the solution v bar is conserved. Consequently, from an arbitrary initial condition, solutions converge to the equilibrium solution which is unique for the given v bar. The problem arises naturally as a continuum limit of a network of certain micro-devices. Each micro-device imperfectly separates a target fluid component from a mixture of more than one component, and its input-output concentration relationships are modeled by a pair of quadratic maps. The solvability of the initial boundary value problem is used to demonstrate that such a network acts as an ideal macro-separator, separating out the target component almost completely. Another network is also proposed which leads to a modified Burgers equation with a nonlinear diffusion coefficient.

  2. Solvability of Boundary Value Problem at Resonance for Third-Order Functional Differential Equations

    Indian Academy of Sciences (India)

    Pinghua Yang; Zengji Du; Weigao Ge

    2008-05-01

    This paper is devoted to the study of boundary value problem of third-order functional differential equations. We obtain some existence results for the problem at resonance under the condition that the nonlinear terms is bounded or generally unbounded. In this paper we mainly use the topological degree theory.

  3. The focal boundary value problem for strongly singular higher-order nonlinear functional-differential equations

    Czech Academy of Sciences Publication Activity Database

    Mukhigulashvili, Sulkhan; Půža, B.

    2015-01-01

    Roč. 2015, January (2015), s. 17. ISSN 1687-2770 Institutional support: RVO:67985840 Keywords : higher order nonlinear functional-differential equations * two-point right-focal boundary value problem * strong singularity Subject RIV: BA - General Mathematics Impact factor: 1.014, year: 2014 http://link.springer.com/article/10.1186%2Fs13661-014-0277-1

  4. Analysis of Blasius Equation for Flat-Plate Flow with Infinite Boundary Value

    DEFF Research Database (Denmark)

    Miansari, M. O.; Miansari, M. E.; Barari, Amin; Domairry, G.

    2010-01-01

    and write the nonlinear differential equation in the state space format, and then solve the initial value problem instead of boundary value problem. The significance of linear part is a key factor in convergence. A first seen linear part may lead to an unstable solution, therefore an extra term is...

  5. Geometry and a priori estimates for free boundary problems of the Euler's equation

    OpenAIRE

    Shatah, Jalal; Zeng, Chongchun

    2006-01-01

    In this paper we derive estimates to the free boundary problem for the Euler equation with surface tension, and without surface tension provided the Rayleigh-Taylor sign condition holds. We prove that as the surface tension tends to zero, when the Rayleigh-Taylor condition is satisfied, solutions converge to the Euler flow with zero surface tension.

  6. Third-order operator-differential equations with discontinuous coefficients and operators in the boundary conditions

    Directory of Open Access Journals (Sweden)

    Araz R. Aliev

    2013-10-01

    Full Text Available We study a third-order operator-differential equation on the semi-axis with a discontinuous coefficient and boundary conditions which include an abstract linear operator. Sufficient conditions for the well-posed and unique solvability are found by means of properties of the operator coefficients in a Sobolev-type space.

  7. PERIODIC BOUNDARY VALUE PROBLEM FOR NONLINEAR INTEGRO-DIFFERENTIAL EQUATION OF MIXED TYPE ON TIME SCALES

    Institute of Scientific and Technical Information of China (English)

    Yepeng Xing; Qiong Wang; Valery G. Romanovski

    2009-01-01

    We prove several new comparison results and develop the monotone iterative tech-nique to show the existence of extremal solutions to a kind of periodic boundary value problem (PBVP) for nonlinear integro-differential equation of mixed type on time scales.

  8. Three-Point Boundary Value Problems of Nonlinear Second-Order q-Difference Equations Involving Different Numbers of q

    Directory of Open Access Journals (Sweden)

    Thanin Sitthiwirattham

    2013-01-01

    Full Text Available We study a new class of three-point boundary value problems of nonlinear second-order q-difference equations. Our problems contain different numbers of q in derivatives and integrals. By using a variety of fixed point theorems (such as Banach’s contraction principle, Boyd and Wong fixed point theorem for nonlinear contractions, Krasnoselskii’s fixed point theorem, and Leray-Schauder nonlinear alternative and Leray-Schauder degree theory, some new existence and uniqueness results are obtained. Illustrative examples are also presented.

  9. Surfaces specified by integrable systems of partial differential equations determined by structure equations and Lax pair

    Science.gov (United States)

    Bracken, Paul

    2010-04-01

    A system of evolution equations can be developed from the structure equations for a submanifold embedded in a three-dimensional space. It is seen how these same equations can be obtained from a generalized matrix Lax pair provided a single constraint equation is imposed. This can be done in Euclidean space as well as in Minkowski space. The integrable systems which result from this process can be thought of as generalizing the SO(3) and SO(2,1) Lax pairs which have been studied previously.

  10. An algorithm of computing inhomogeneous differential equations for definite integrals

    OpenAIRE

    Nakayama, Hiromasa; Nishiyama, Kenta

    2010-01-01

    We give an algorithm to compute inhomogeneous differential equations for definite integrals with parameters. The algorithm is based on the integration algorithm for $D$-modules by Oaku. Main tool in the algorithm is the Gr\\"obner basis method in the ring of differential operators.

  11. Long-term behavior of reaction-diffusion equations with nonlocal boundary conditions on rough domains

    Science.gov (United States)

    Gal, Ciprian G.; Warma, Mahamadi

    2016-08-01

    We investigate the long term behavior in terms of finite dimensional global and exponential attractors, as time goes to infinity, of solutions to a semilinear reaction-diffusion equation on non-smooth domains subject to nonlocal Robin boundary conditions, characterized by the presence of fractional diffusion on the boundary. Our results are of general character and apply to a large class of irregular domains, including domains whose boundary is Hölder continuous and domains which have fractal-like geometry. In addition to recovering most of the existing results on existence, regularity, uniqueness, stability, attractor existence, and dimension, for the well-known reaction-diffusion equation in smooth domains, the framework we develop also makes possible a number of new results for all diffusion models in other non-smooth settings.

  12. Attractor of Beam Equation with Structural Damping under Nonlinear Boundary Conditions

    Directory of Open Access Journals (Sweden)

    Danxia Wang

    2015-01-01

    Full Text Available Simultaneously, considering the viscous effect of material, damping of medium, and rotational inertia, we study a kind of more general Kirchhoff-type extensible beam equation utt-uxxtt+uxxxx-σ(∫0l‍(ux2dxuxx-ϕ(∫0l‍(ux2dxuxxt=q(x, in  [0,L]×R+ with the structural damping and the rotational inertia term. Little attention is paid to the longtime behavior of the beam equation under nonlinear boundary conditions. In this paper, under nonlinear boundary conditions, we prove not only the existence and uniqueness of global solutions by prior estimates combined with some inequality skills, but also the existence of a global attractor by the existence of an absorbing set and asymptotic compactness of corresponding solution semigroup. In addition, the same results also can be proved under the other nonlinear boundary conditions.

  13. Fundamental solutions for transversely isotropic magneto-electro-elastic media and boundary integral formulation

    Institute of Scientific and Technical Information of China (English)

    丁皓江; 江爱民

    2003-01-01

    To obtain the fundamental solutions for computation of magneto-electro-elastic media by the boundary element method, the general solutions in the case of distinct eigenvalues are derived and expressed in five harmonic functions from the governing equations and the strict differential operator theorem. On the basis of these general solutions, the fundamental solution of infinite magneto-electro-elastic solid are obtained with the method of trial-and-error. Finally, the boundary integral formulation is derived and the corresponding boundary element method program is implemented to perform two numerical calculations(a column under uni-axial tension, uniform electric displacement or uniform magnetic induction, an annular plate simply-supported on outer and inner surfaces under axial loads). The numerical results agree well with the analytical ones.

  14. Generalized analytical solution for advection-dispersion equation in finite spatial domain with arbitrary time-dependent inlet boundary condition

    Directory of Open Access Journals (Sweden)

    J.-S. Chen

    2011-04-01

    Full Text Available This study presents a generalized analytical solution for one-dimensional solute transport in finite spatial domain subject to arbitrary time-dependent inlet boundary condition. The governing equation includes terms accounting for advection, hydrodynamic dispersion, linear equilibrium sorption and first order decay processes. The generalized analytical solution is derived by using the Laplace transform with respect to time and the generalized integral transform technique with respect to the spatial coordinate. Several special cases are presented and compared to illustrate the robustness of the derived generalized analytical solution. Result shows an excellent agreement. The analytical solutions of the special cases derived in this study have practical applications. Moreover, the derived generalized solution which consists an integral representation is evaluated by the numerical integration to extend its usage. The developed generalized solution offers a convenient tool for further development of analytical solution of specified time-dependent inlet boundary conditions or numerical evaluation of the concentration field for arbitrary time-dependent inlet boundary problem.

  15. Estimation of parameters of boundary value problems for linear ordinary differential equations with uncertain data

    CERN Document Server

    Nakonechnyi, Olexandr; Shestopalov, Yury

    2009-01-01

    In this paper we construct optimal, in certain sense, estimates of values of linear functionals on solutions to two-point boundary value problems (BVPs) for systems of linear first-order ordinary differential equations from observations which are linear transformations of the same solutions perturbed by additive random noises. It is assumed here that right-hand sides of equations and boundary data as well as statistical characteristics of random noises in observations are not known and belong to certain given sets in corresponding functional spaces. This leads to the necessity of introducing minimax statement of an estimation problem when optimal estimates are defined as linear, with respect to observations, estimates for which the maximum of mean square error of estimation taken over the above-mentioned sets attains minimal value. Such estimates are called minimax estimates. We establish that the minimax estimates are expressed via solutions of some systems of differential equations of special type. Similar ...

  16. $A_n^{(1)}$ affine Toda field theories with integrable boundary conditions revisited

    CERN Document Server

    Doikou, Anastasia

    2008-01-01

    Generic classically integrable boundary conditions for the $A_{n}^{(1)}$ affine Toda field theories (ATFT) are investigated. The present analysis relies primarily on the underlying algebra, defined by the classical version of the reflection equation. We use as a prototype example the first non-trivial model of the hierarchy i.e. the $A_2^{(1)}$ ATFT, however our results may be generalized for any $A_{n}^{(1)}$ ($n>1$). We assume here two distinct types of boundary conditions called some times soliton preserving (SP), and soliton non-preserving (SNP) associated to two distinct algebras, i.e. the reflection algebra and the ($q$) twisted Yangian respectively. The boundary local integrals of motion are then systematically extracted from the asymptotic expansion of the associated transfer matrix. In the case of SNP boundary conditions we recover previously known results. The other type of boundary conditions (SP), associated to the reflection algebra, are novel in this context and lead to a different set of conser...

  17. A treecode-accelerated boundary integral Poisson-Boltzmann solver for electrostatics of solvated biomolecules

    Science.gov (United States)

    Geng, Weihua; Krasny, Robert

    2013-08-01

    We present a treecode-accelerated boundary integral (TABI) solver for electrostatics of solvated biomolecules described by the linear Poisson-Boltzmann equation. The method employs a well-conditioned boundary integral formulation for the electrostatic potential and its normal derivative on the molecular surface. The surface is triangulated and the integral equations are discretized by centroid collocation. The linear system is solved by GMRES iteration and the matrix-vector product is carried out by a Cartesian treecode which reduces the cost from O(N2) to O(NlogN), where N is the number of faces in the triangulation. The TABI solver is applied to compute the electrostatic solvation energy in two cases, the Kirkwood sphere and a solvated protein. We present the error, CPU time, and memory usage, and compare results for the Poisson-Boltzmann and Poisson equations. We show that the treecode approximation error can be made smaller than the discretization error, and we compare two versions of the treecode, one with uniform clusters and one with non-uniform clusters adapted to the molecular surface. For the protein test case, we compare TABI results with those obtained using the grid-based APBS code, and we also present parallel TABI simulations using up to eight processors. We find that the TABI solver exhibits good serial and parallel performance combined with relatively simple implementation, efficient memory usage, and geometric adaptability.

  18. Boundary Integral Equation Methods for Conical Diffraction and Short Waves

    OpenAIRE

    Goray, Leonid; Schmidt, Gunther

    2014-01-01

    This work is part of research that has been pursued by the authors over a long period of time for the purpose of developing accurate and fast numerical algorithms, including the commercial packages PCGrate and DiPoG designed to model multilayered gratings having mostly one-dimensional periodicity, including roughness, and working in all, including the shortest, optical wavelength ranges at arbitrary optical mounts. The main purpose of this Chapter is to present a complete description in gener...

  19. Convergence of a Boundary Integral Method for Water Waves

    OpenAIRE

    Beale, J. Thomas; Hou, Thomas Y.; Lowengrub, John

    1996-01-01

    We prove nonlinear stability and convergence of certain boundary integral methods for time-dependent water waves in a two-dimensional, inviscid, irrotational, incompressible fluid, with or without surface tension. The methods are convergent as long as the underlying solution remains fairly regular (and a sign condition holds in the case without surface tension). Thus, numerical instabilities are ruled out even in a fully nonlinear regime. The analysis is based on delicate energy estimates, fo...

  20. Algebraic Integrability of Lotka-Volterra equations in three dimensions

    CERN Document Server

    Constandinides, Kyriacos

    2009-01-01

    We examine the algebraic complete integrability of Lotka-Volterra equations in three dimensions. We restrict our attention to Lotka-Volterra systems defined by a skew symmetric matrix. We obtain a complete classification of such systems. The classification is obtained using Painleve analysis and more specifically by the use of Kowalevski exponents. The imposition of certain integrality conditions on the Kowalevski exponents gives necessary conditions for the algebraic integrability of the corresponding systems. We also show that the conditions are sufficient.

  1. Exact solutions to quantum field theories and integrable equations

    OpenAIRE

    Marshakov, A.

    1996-01-01

    The exact solutions to quantum string and gauge field theories are discussed and their formulation in the framework of integrable systems is presented. In particular I consider in detail several examples of appearence of solutions to the first-order integrable equations of hydrodynamical type and stress that all known examples can be treated as partial solutions to the same problem in the theory of integrable systems.

  2. The second boundary value problem for equations of viscoelastic diffusion in polymers

    CERN Document Server

    Vorotnikov, Dmitry A

    2009-01-01

    The classical approach to diffusion processes is based on Fick's law that the flux is proportional to the concentration gradient. Various phenomena occurring during propagation of penetrating liquids in polymers show that this type of diffusion exhibits anomalous behavior and contradicts the just mentioned law. However, they can be explained in the framework of non-Fickian diffusion theories based on viscoelasticity of polymers. Initial-boundary value problems for viscoelastic diffusion equations have been studied by several authors. Most of the studies are devoted to the Dirichlet BVP (the concentration is given on the boundary of the domain). In this chapter we study the second BVP, i.e. when the normal component of the concentration flux is prescribed on the boundary, which is more realistic in many physical situations. We establish existence of weak solutions to this problem. We suggest some conditions on the coefficients and boundary data under which all the solutions tend to the homogeneous state as tim...

  3. Effect of the Lower Boundary Position of the Fourier Equation on the Soil Energy Balance

    Institute of Scientific and Technical Information of China (English)

    孙菽芬; 张霞

    2004-01-01

    In this study, the effect of the lower boundary position selection for the Fourier equation on heat transfer and energy balance in soil is evaluated. A detailed numerical study shows that the proper position of the lower boundary is critical when solving the Fourier equation by using zero heat flux as the lower boundary condition. Since the position defines the capacity of soil as a heat sink or source, which absorbs and stores radiation energy from the sky in summer and then releases the energy to the atmosphere in winter, and regulates the deep soil temperature distribution, the depth of the position greatly influences the heat balance within the soil as well as the interaction between the soil and the atmosphere. Based on physical reasoning and the results of numerical simulation, the proper depth of the position should be equal to approximately 3 times of the annual heat wave damping depth. For most soils, the proper lower boundary depth for the Fourier equation should be around 8 m to 15 m, depending on soil texture.

  4. Boundary value problem for one-dimensional fractional differential advection-dispersion equation

    Directory of Open Access Journals (Sweden)

    Khasambiev Mokhammad Vakhaevich

    2014-07-01

    Full Text Available An equation commonly used to describe solute transport in aquifers has attracted more attention in recent years. After a formal study of some aspects of the advection-diffusion equation, basically from the mathematical point of view with the solution of a differential equation with fractional derivative, the main interest to this problem shifted onto physical aspects of the dynamical system, such as the total energy and the dynamical response. In this regard it should be pointed out that the interaction with environment is expressed in terms of stochastic arrow of time. This allows one also to reach a progress in one more issue. Formerly the equation of advection-diffusion was not obtained from any physical principles. However, mainly the success concerns linear fractional systems. In fact, there are many cases in which linear treatments are not sufficient. The more general systems described by nonlinear fractional differential equations have not been studied enough. The ordinary calculus brings out clearly that essentially new phenomena occur in nonlinear systems, which generally cannot occur in linear systems. Due to vast range of application of the fractional advection-dispersion equation, a lot of work has been done to find numerical solution and fundamental solution of this equation. The research on the analytical solution of initial-boundary problem for space-fractional advection-dispersion equation is relatively new and is still at an early stage of development. In this paper, we will take use of the method of variable separation to solve space-fractional advection-dispersion equation with initial boundary data.

  5. Velocity boundary conditions for vorticity formulations of the incompressible Navier-Stokes equations

    Energy Technology Data Exchange (ETDEWEB)

    Kempka, S.N.; Strickland, J.H.; Glass, M.W.; Peery, J.S. [Sandia National Labs., Albuquerque, NM (United States); Ingber, M.S. [Univ. of New Mexico, Albuquerque, NM (United States)

    1995-04-01

    formulation to satisfy velocity boundary conditions for the vorticity form of the incompressible, viscous fluid momentum equations is presented. The tangential and normal components of the velocity boundary condition are satisfied simultaneously by creating vorticity adjacent to boundaries. The newly created vorticity is determined using a kinematical formulation which is a generalization of Helmholtz` decomposition of a vector field. Though it has not been generally recognized, these formulations resolve the over-specification issue associated with creating voracity to satisfy velocity boundary conditions. The generalized decomposition has not been widely used, apparently due to a lack of a useful physical interpretation. An analysis is presented which shows that the generalized decomposition has a relatively simple physical interpretation which facilitates its numerical implementation. The implementation of the generalized decomposition is discussed in detail. As an example the flow in a two-dimensional lid-driven cavity is simulated. The solution technique is based on a Lagrangian transport algorithm in the hydrocode ALEGRA. ALEGRA`s Lagrangian transport algorithm has been modified to solve the vorticity transport equation and the generalized decomposition, thus providing a new, accurate method to simulate incompressible flows. This numerical implementation and the new boundary condition formulation allow vorticity-based formulations to be used in a wider range of engineering problems.

  6. A Numerical Method for Solving the 3D Unsteady Incompressible Navier-Stokes Equations in Curvilinear Domains with Complex Immersed Boundaries.

    Science.gov (United States)

    Ge, Liang; Sotiropoulos, Fotis

    2007-08-01

    A novel numerical method is developed that integrates boundary-conforming grids with a sharp interface, immersed boundary methodology. The method is intended for simulating internal flows containing complex, moving immersed boundaries such as those encountered in several cardiovascular applications. The background domain (e.g the empty aorta) is discretized efficiently with a curvilinear boundary-fitted mesh while the complex moving immersed boundary (say a prosthetic heart valve) is treated with the sharp-interface, hybrid Cartesian/immersed-boundary approach of Gilmanov and Sotiropoulos [1]. To facilitate the implementation of this novel modeling paradigm in complex flow simulations, an accurate and efficient numerical method is developed for solving the unsteady, incompressible Navier-Stokes equations in generalized curvilinear coordinates. The method employs a novel, fully-curvilinear staggered grid discretization approach, which does not require either the explicit evaluation of the Christoffel symbols or the discretization of all three momentum equations at cell interfaces as done in previous formulations. The equations are integrated in time using an efficient, second-order accurate fractional step methodology coupled with a Jacobian-free, Newton-Krylov solver for the momentum equations and a GMRES solver enhanced with multigrid as preconditioner for the Poisson equation. Several numerical experiments are carried out on fine computational meshes to demonstrate the accuracy and efficiency of the proposed method for standard benchmark problems as well as for unsteady, pulsatile flow through a curved, pipe bend. To demonstrate the ability of the method to simulate flows with complex, moving immersed boundaries we apply it to calculate pulsatile, physiological flow through a mechanical, bileaflet heart valve mounted in a model straight aorta with an anatomical-like triple sinus. PMID:19194533

  7. An initial-boundary value problem in a strip for a two-dimensional equation of Zakharov-Kuznetsov type

    OpenAIRE

    Faminskii, Andrei V.

    2013-01-01

    An initial-boundary value problem in a strip with homogeneous Diriclet boundary conditions for two-dimensional generalized Zakharov-Kuznetsov equation is considered. In particular, dissipative and absorbing degenerate terms can be supplemented to the original Zakharov-Kuznetsov equation. Results on global existence, uniqueness and long-time decay of weak silutions are established.

  8. On the physical solutions to the heat equation subjected to nonlinear boundary conditions

    International Nuclear Information System (INIS)

    This work consists of a discussion on the physical solutions to the steady-state heat transfer equation, when it is subjected to nonlinear boundary conditions. It will be presented a functional, whose minimum occurs for the (unique) physical solution to the condidered heat transfer problem, suitable for a large class of typical (nonlinear) boundary conditions (representing the radiative/convective loss from the body to the environment). It will be demonstrated that these problems admit-always one, and only one, physical solution (which represents the absolute temperature). (author)

  9. Commentary on local and boundary regularity of weak solutions to Navier-Stokes equations

    Directory of Open Access Journals (Sweden)

    Zdenek Skalak

    2004-01-01

    Full Text Available We present results on local and boundary regularity for weak solutions to the Navier-Stokes equations. Beginning with the regularity criterion proved recently in [14] for the Cauchy problem, we show that this criterion holds also locally. This is also the case for some other results: We present two examples concerning the regularity of weak solutions stemming from the regularity of two components of the vorticity ([2] or from the regularity of the pressure ([3]. We conclude by presenting regularity criteria near the boundary based on the results in [10] and [16].

  10. Global existence and exponential growth for a viscoelastic wave equation with dynamic boundary conditions

    KAUST Repository

    Gerbi, Stéphane

    2013-01-15

    The goal of this work is to study a model of the wave equation with dynamic boundary conditions and a viscoelastic term. First, applying the Faedo-Galerkin method combined with the fixed point theorem, we show the existence and uniqueness of a local in time solution. Second, we show that under some restrictions on the initial data, the solution continues to exist globally in time. On the other hand, if the interior source dominates the boundary damping, then the solution is unbounded and grows as an exponential function. In addition, in the absence of the strong damping, then the solution ceases to exist and blows up in finite time.

  11. Field Method for Integrating the First Order Differential Equation

    Institute of Scientific and Technical Information of China (English)

    JIA Li-qun; ZHENG Shi-wang; ZHANG Yao-yu

    2007-01-01

    An important modern method in analytical mechanics for finding the integral, which is called the field-method, is used to research the solution of a differential equation of the first order. First, by introducing an intermediate variable, a more complicated differential equation of the first order can be expressed by two simple differential equations of the first order, then the field-method in analytical mechanics is introduced for solving the two differential equations of the first order. The conclusion shows that the field-method in analytical mechanics can be fully used to find the solutions of a differential equation of the first order, thus a new method for finding the solutions of the first order is provided.

  12. Transonic integro-differential and integral equations with artificial viscosity

    International Nuclear Information System (INIS)

    From the two-dimensional transonic small disturbance equation, with artificial viscosity, we derive an integro-differential equation and subsequently an integral equation. The computational domain is discretized into rectangular elements and functions of the dependent variable and its derivatives are assumed to be constant in each element. The resulting nonlinear algebraic systems are solved by Jacobi iteration. The method is tested for parabolic-arc and NACA0012 airfoils. Convergence is fast and the solutions compare well with finite-difference results, despite the use of a comparably small number of nodes. (author). 18 refs, 7 figs

  13. Recursive integral equations with positive kernel for lattice calculations

    International Nuclear Information System (INIS)

    A Kirkwood-Salzburg integral equation, with positive defined kernel, for the states of lattice models of statistical mechanics and quantum field theory is derived. The equation is defined in the thermodynamic limit, and its iterative solution is convergent. Moreover, positivity leads to an exact a priori bound on the iteration. The equation's relevance as a reliable algorithm for lattice calculations is therefore suggested, and it is illustrated with a simple application. It should provide a viable alternative to Monte Carlo methods for models of statistical mechanics and lattice gauge theories. 10 refs

  14. An integral equation solution for multistage turbomachinery design calculations

    Science.gov (United States)

    Mcfarland, Eric R.

    1993-01-01

    A method was developed to calculate flows in multistage turbomachinery. The method is an extension of quasi-three-dimensional blade-to-blade solution methods. Governing equations for steady compressible inviscid flow are linearized by introducing approximations. The linearized flow equations are solved using integral equation techniques. The flows through both stationary and rotating blade rows are determined in a single calculation. Multiple bodies can be modelled for each blade row, so that arbitrary blade counts can be analyzed. The method's benefits are its speed and versatility.

  15. Normalized RBF networks: application to a system of integral equations

    International Nuclear Information System (INIS)

    Linear integral and integro-differential equations of Fredholm and Volterra types have been successfully treated using radial basis function (RBF) networks in previous works. This paper deals with the case of a system of integral equations of Fredholm and Volterra types with a normalized radial basis function (NRBF) network. A novel learning algorithm is developed for the training of NRBF networks in which the BFGS backpropagation (BFGS-BP) least-squares optimization method as a recursive model is used. In the approach presented here, a trial solution is given by an NRBF network of incremental architecture with a set of unknown parameters. Detailed learning algorithms and concrete examples are also included

  16. An approximate method to study the one-velocity neutron integral transport equation

    International Nuclear Information System (INIS)

    An approximate method to study the monokinetic linear transport equation is outlined, starting from its integral form, rather than the integro-differential one. The approximate solution may be deduced either analytically, in simple cases, or numerically by means of typical space discretization techniques, through a system of second-order differential equations, associated with proper boundary conditions. Both the system and the boundary conditions may be matched with the standard neutron diffusion multigroup ones, by means of a proper correspondence of the coefficients and of the unknowns. The slab and the radially-symmetric sphere are then analysed in detail. It is shown how, in the plane case, the present approximation is perfectly equivalent to the well-known discrete ordinate one. For curved geometries no such equivalence exists, and it is in these cases that the application of the method at hand looks promising, in order to avoid complications and numerical problems in practical applications. (author)

  17. A diffusive Fisher-KPP equation with free boundaries and time-periodic advections

    OpenAIRE

    Sun, Ningkui; Lou, Bendong; Zhou, Maolin

    2016-01-01

    We consider a reaction-diffusion-advection equation of the form: $u_t=u_{xx}-\\beta(t)u_x+f(t,u)$ for $x\\in (g(t),h(t))$, where $\\beta(t)$ is a $T$-periodic function representing the intensity of the advection, $f(t,u)$ is a Fisher-KPP type of nonlinearity, $T$-periodic in $t$, $g(t)$ and $h(t)$ are two free boundaries satisfying Stefan conditions. This equation can be used to describe the population dynamics in time-periodic environment with advection. Its homogeneous version (that is, both $...

  18. A boundary value approach for solving three-dimensional elliptic and hyperbolic partial differential equations.

    Science.gov (United States)

    Biala, T A; Jator, S N

    2015-01-01

    In this article, the boundary value method is applied to solve three dimensional elliptic and hyperbolic partial differential equations. The partial derivatives with respect to two of the spatial variables (y, z) are discretized using finite difference approximations to obtain a large system of ordinary differential equations (ODEs) in the third spatial variable (x). Using interpolation and collocation techniques, a continuous scheme is developed and used to obtain discrete methods which are applied via the Block unification approach to obtain approximations to the resulting large system of ODEs. Several test problems are investigated to elucidate the solution process. PMID:26543723

  19. Diffusive KPP Equations with Free Boundaries in Time Almost Periodic Environments: II. Spreading Speeds and Semi-Wave

    OpenAIRE

    Li, Fang; Liang, Xing; Shen, Wenxian

    2015-01-01

    In this series of papers, we investigate the spreading and vanishing dynamics of time almost periodic diffusive KPP equations with free boundaries. Such equations are used to characterize the spreading of a new species in time almost periodic environments with free boundaries representing the spreading fronts. In the first part of the series, we showed that a spreading-vanishing dichotomy occurs for such free boundary problems (see [16]). In this second part of the series, we investigate the ...

  20. The Edwards-Wilkinson equation with drift and Neumann boundary conditions

    International Nuclear Information System (INIS)

    The well-known scaling of the Edwards-Wilkinson equation is essentially determined by dimensional analysis. In a range of experimental setups, be it due to the presence of an electrical or a gravitational field, or the indirect effect of other terms or an expansion, an additional drift term has to be considered. Once the drift term is added, more sophisticated reasoning is required to determine the scaling, which initially suggests that the drift term dominates the diffusion. However, the diffusion term is dangerously irrelevant and the resulting scaling in fact non-trivial. In order to assess the universality of the Edwards-Wilkinson equation with drift and to describe a physically more relevant situation, we compare the scaling obtained with Neumann boundary conditions to the published case with Dirichlet boundary conditions.

  1. Mixed direct-iterative methods for boundary integral formulations of continuum dielectric solvation models

    Energy Technology Data Exchange (ETDEWEB)

    Corcelli, S.A.; Kress, J.D.; Pratt, L.R.

    1995-08-07

    This paper develops and characterizes mixed direct-iterative methods for boundary integral formulations of continuum dielectric solvation models. We give an example, the Ca{sup ++}{hor_ellipsis}Cl{sup {minus}} pair potential of mean force in aqueous solution, for which a direct solution at thermal accuracy is difficult and, thus for which mixed direct-iterative methods seem necessary to obtain the required high resolution. For the simplest such formulations, Gauss-Seidel iteration diverges in rare cases. This difficulty is analyzed by obtaining the eigenvalues and the spectral radius of the non-symmetric iteration matrix. This establishes that those divergences are due to inaccuracies of the asymptotic approximations used in evaluation of the matrix elements corresponding to accidental close encounters of boundary elements on different atomic spheres. The spectral radii are then greater than one for those diverging cases. This problem is cured by checking for boundary element pairs closer than the typical spatial extent of the boundary elements and for those cases performing an ``in-line`` Monte Carlo integration to evaluate the required matrix elements. These difficulties are not expected and have not been observed for the thoroughly coarsened equations obtained when only a direct solution is sought. Finally, we give an example application of hybrid quantum-classical methods to deprotonation of orthosilicic acid in water.

  2. Transient Growth Analysis of Compressible Boundary Layers with Parabolized Stability Equations

    Science.gov (United States)

    Paredes, Pedro; Choudhari, Meelan M.; Li, Fei; Chang, Chau-Lyan

    2016-01-01

    The linear form of parabolized linear stability equations (PSE) is used in a variational approach to extend the previous body of results for the optimal, non-modal disturbance growth in boundary layer flows. This methodology includes the non-parallel effects associated with the spatial development of boundary layer flows. As noted in literature, the optimal initial disturbances correspond to steady counter-rotating stream-wise vortices, which subsequently lead to the formation of stream-wise-elongated structures, i.e., streaks, via a lift-up effect. The parameter space for optimal growth is extended to the hypersonic Mach number regime without any high enthalpy effects, and the effect of wall cooling is studied with particular emphasis on the role of the initial disturbance location and the value of the span-wise wavenumber that leads to the maximum energy growth up to a specified location. Unlike previous predictions that used a basic state obtained from a self-similar solution to the boundary layer equations, mean flow solutions based on the full Navier-Stokes (NS) equations are used in select cases to help account for the viscous-inviscid interaction near the leading edge of the plate and also for the weak shock wave emanating from that region. These differences in the base flow lead to an increasing reduction with Mach number in the magnitude of optimal growth relative to the predictions based on self-similar mean-flow approximation. Finally, the maximum optimal energy gain for the favorable pressure gradient boundary layer near a planar stagnation point is found to be substantially weaker than that in a zero pressure gradient Blasius boundary layer.

  3. The pentabox Master Integrals with the Simplified Differential Equations approach

    Science.gov (United States)

    Papadopoulos, Costas G.; Tommasini, Damiano; Wever, Christopher

    2016-04-01

    We present the calculation of massless two-loop Master Integrals relevant to five-point amplitudes with one off-shell external leg and derive the complete set of planar Master Integrals with five on-mass-shell legs, that contribute to many 2 → 3 amplitudes of interest at the LHC, as for instance three jet production, γ , V, H + 2 jets etc., based on the Simplified Differential Equations approach.

  4. The Pentabox Master Integrals with the Simplified Differential Equations approach

    CERN Document Server

    Papadopoulos, Costas G; Wever, Christopher

    2015-01-01

    We present the calculation of massless two-loop Master Integrals relevant to five-point amplitudes with one off-shell external leg and derive the complete set of planar Master Integrals with five on-mass-shell legs, that contribute to many $2\\to 3$ amplitudes of interest at the LHC, as for instance three jet production, $\\gamma, V, H +2$ jets etc., based on the Simplified Differential Equations approach.

  5. Numerical Integration of the Transport Equation For Infinite Homogeneous Media

    International Nuclear Information System (INIS)

    The transport equation for neutrons in infinite homogeneous media is solved by direct numerical integration. Accounts are taken to the anisotropy and the inelastic scattering. The integration has been performed by means of the trapezoidal rule and the length of the energy intervals are constant in lethargy scale. The machine used is a Ferranti Mercury computer. Results are given for water, heavy water, aluminium water mixture and iron-aluminium-water mixture

  6. Extension to the integral transport equation of an iterative method

    Energy Technology Data Exchange (ETDEWEB)

    Jehouani, A. [EPRA, Department of Physics, Faculty of Sciences, Semlalia, PO Box 2390, 40000 Marrakech (Morocco)]. E-mail: jehouani@ucam.ac.ma; Elmorabiti, A. [Centre d' Etudes Nucleaires de Maamoura, CNESTEN (Morocco); Ghassoun, J. [EPRA, Department of Physics, Faculty of Sciences, Semlalia, PO Box 2390, 40000 Marrakech (Morocco)

    2006-09-15

    In this paper we describe an extension to the neutron integral transport equation of an iterative method. Indeed an iterative scheme is used for both energy and space including external iteration for the multiplication factor and internal iteration for flux calculations. The Monte Carlo method is used to evaluate the spatial transfer integrals. The results were compared with those obtained by using the APOLLO2 code for a cylindrical cell.

  7. Extension to the integral transport equation of an iterative method

    International Nuclear Information System (INIS)

    In this paper we describe an extension to the neutron integral transport equation of an iterative method. Indeed an iterative scheme is used for both energy and space including external iteration for the multiplication factor and internal iteration for flux calculations. The Monte Carlo method is used to evaluate the spatial transfer integrals. The results were compared with those obtained by using the APOLLO2 code for a cylindrical cell

  8. Continuum and Discrete Initial-Boundary Value Problems and Einstein's Field Equations

    Directory of Open Access Journals (Sweden)

    Olivier Sarbach

    2012-08-01

    Full Text Available Many evolution problems in physics are described by partial differential equations on an infinite domain; therefore, one is interested in the solutions to such problems for a given initial dataset. A prominent example is the binary black-hole problem within Einstein's theory of gravitation, in which one computes the gravitational radiation emitted from the inspiral of the two black holes, merger and ringdown. Powerful mathematical tools can be used to establish qualitative statements about the solutions, such as their existence, uniqueness, continuous dependence on the initial data, or their asymptotic behavior over large time scales. However, one is often interested in computing the solution itself, and unless the partial differential equation is very simple, or the initial data possesses a high degree of symmetry, this computation requires approximation by numerical discretization. When solving such discrete problems on a machine, one is faced with a finite limit to computational resources, which leads to the replacement of the infinite continuum domain with a finite computer grid. This, in turn, leads to a discrete initial-boundary value problem. The hope is to recover, with high accuracy, the exact solution in the limit where the grid spacing converges to zero with the boundary being pushed to infinity. The goal of this article is to review some of the theory necessary to understand the continuum and discrete initial boundary-value problems arising from hyperbolic partial differential equations and to discuss its applications to numerical relativity; in particular, we present well-posed initial and initial-boundary value formulations of Einstein's equations, and we discuss multi-domain high-order finite difference and spectral methods to solve them.

  9. Robin boundary value problems for elliptic operational differential equations with variable coefficients

    Directory of Open Access Journals (Sweden)

    Rabah Haoua

    2015-04-01

    Full Text Available In this article we give some new results on abstract second-order differential equations of elliptic type with variable operator coefficients and general Robin boundary conditions, in the framework of Holder spaces. We assume that the family of variable coefficients verify the well known Labbas-Terreni assumption used in the sum theory. We use Dunford calculus, interpolation spaces and the semigroup theory to obtain existence, uniqueness and maximal regularity results for the solution of the problem.

  10. On the Spectrum of Neutron Transport Equations with Reflecting Boundary Conditions

    OpenAIRE

    Song, Degong

    2000-01-01

    This dissertation is devoted to investigating the time dependent neutron transport equations with reflecting boundary conditions. Two typical geometries --- slab geometry and spherical geometry --- are considered in the setting of L^p including L^1. Some aspects of the spectral properties of the transport operator A and the strongly continuous semigroup T(t) generated by A are studied. It is shown under fairly general assumptions that the accumulation points of { m Pas...

  11. Maximum Principles and Boundary Value Problems for First-Order Neutral Functional Differential Equations

    Directory of Open Access Journals (Sweden)

    Domoshnitsky Alexander

    2009-01-01

    Full Text Available We obtain the maximum principles for the first-order neutral functional differential equation where , and are linear continuous operators, and are positive operators, is the space of continuous functions, and is the space of essentially bounded functions defined on . New tests on positivity of the Cauchy function and its derivative are proposed. Results on existence and uniqueness of solutions for various boundary value problems are obtained on the basis of the maximum principles.

  12. Adaptive finite difference methods for nonlinear elliptic and parabolic partial differential equations with free boundaries

    OpenAIRE

    Oberman, Adam M.; Zwiers, Ian

    2014-01-01

    Monotone finite difference methods provide stable convergent discretizations of a class of degenerate elliptic and parabolic Partial Differential Equations (PDEs). These methods are best suited to regular rectangular grids, which leads to low accuracy near curved boundaries or singularities of solutions. In this article we combine monotone finite difference methods with an adaptive grid refinement technique to produce a PDE discretization and solver which is applied to a broad class of equati...

  13. Uniform regularity for the Navier-Stokes equation with Navier boundary condition

    OpenAIRE

    Masmoudi, Nader; Rousset, Frederic

    2010-01-01

    We prove that there exists an interval of time which is uniform in the vanishing viscosity limit and for which the Navier-Stokes equation with Navier boundary condition has a strong solution. This solution is uniformly bounded in a conormal Sobolev space and has only one normal derivative bounded in $L^\\infty$. This allows to get the vanishing viscosity limit to the incompressible Euler system from a strong compactness argument.

  14. On a two-point boundary value problem for second order ordinary differential equations with singularities

    Czech Academy of Sciences Publication Activity Database

    Lomtatidze, Alexander; Malaguti, L.

    2003-01-01

    Roč. 52, č. 6 (2003), s. 1553-1567. ISSN 0362-546X R&D Projects: GA ČR GA201/99/0295 Institutional research plan: CEZ:AV0Z1019905; CEZ:AV0Z1019905 Keywords : second order singular differential equation * two-point boundary value problem * lower and upper functions Subject RIV: BA - General Mathematics Impact factor: 0.354, year: 2003

  15. Solutions to fourth-order random differential equations with periodic boundary conditions

    Directory of Open Access Journals (Sweden)

    Xiaoling Han

    2012-12-01

    Full Text Available Existence of solutions and of extremal random solutions are proved for periodic boundary-value problems of fourth-order ordinary random differential equations. Our investigation is done in the space of continuous real-valued functions defined on closed and bounded intervals. Also we study the applications of the random version of a nonlinear alternative of Leray-Schauder type and an algebraic random fixed point theorem by Dhage.

  16. Stable and chaotic solutions of the complex Ginzburg-Landau equation with periodic boundary conditions

    OpenAIRE

    Scheuer, Jacob; Malomed, Boris A.

    2001-01-01

    We study, analytically and numerically, the dynamical behavior of the solutions of the complex Ginzburg-Landau equation with diffraction but without diffusion, which governs the spatial evolution of the field in an active nonlinear laser cavity. Accordingly, the solutions are subject to periodic boundary conditions. The analysis reveals regions of stable stationary solutions in the model’s parameter space, and a wide range of oscillatory and chaotic behaviors. Close to the first bifurca...

  17. PERIODIC BOUNDARY VALUE PROBLEM AND CAUCHY PROBLEM OF THE GENERALIZED CUBIC DOUBLE DISPERSION EQUATION

    Institute of Scientific and Technical Information of China (English)

    Chen Guowang; Xue Hongxia

    2008-01-01

    In this article, the existence, uniqueness and regularities of the global gener-alized solution and global classical solution for the periodic boundary value problem and the Cauchy problem of the general cubic double dispersion equation utt -uxx-auxxtt+bux4 - duxxt= f(u)xx are proved, and the sufficient conditions of blow-up of the solutions for the Cauchy problems in finite time are given.

  18. Existence of solutions for boundary value problem of fractional order impulsive differential equations systems

    Directory of Open Access Journals (Sweden)

    Weihua JIANG

    2015-04-01

    Full Text Available By defining appropriate linear space and norm, giving the appropriate operator, using the contraction mapping principle and krasnoselskii fixed point theorem respectively, the existence and uniqueness of solutions for boundary value problem of fractional order impulsive differential equations systems are investigated under certain condition that nonlinear term and pulse value are satisfied. An example is given to illustrate that the required conditions can be satisfied.

  19. The second boundary value problem for equations of viscoelastic diffusion in polymers

    OpenAIRE

    Vorotnikov, Dmitry A.

    2009-01-01

    The classical approach to diffusion processes is based on Fick's law that the flux is proportional to the concentration gradient. Various phenomena occurring during propagation of penetrating liquids in polymers show that this type of diffusion exhibits anomalous behavior and contradicts the just mentioned law. However, they can be explained in the framework of non-Fickian diffusion theories based on viscoelasticity of polymers. Initial-boundary value problems for viscoelastic diffusion equat...

  20. Wavelet multiresolution analyses adapted for the fast solution of boundary value ordinary differential equations

    Science.gov (United States)

    Jawerth, Bjoern; Sweldens, Wim

    1993-01-01

    We present ideas on how to use wavelets in the solution of boundary value ordinary differential equations. Rather than using classical wavelets, we adapt their construction so that they become (bi)orthogonal with respect to the inner product defined by the operator. The stiffness matrix in a Galerkin method then becomes diagonal and can thus be trivially inverted. We show how one can construct an O(N) algorithm for various constant and variable coefficient operators.

  1. Navier–Stokes equations on a flat cylinder with vorticity production on the boundary

    International Nuclear Information System (INIS)

    We study a singular version of the incompressible two-dimensional Navier–Stokes (NS) system on a flat cylinder C, with Neumann conditions for the vorticity and a vorticity production term on the boundary ∂C to restore the no-slip boundary condition for the velocity u|∂C = 0. The problem is formulated as an infinite system of coupled ordinary differential equations (ODEs) for the Neumann Fourier modes. For a general class of initial data we prove existence and uniqueness of the solution, and equivalence to the usual NS system. The main tool in the proofs is a suitable decay of the modes, obtained by the explicit form of the ODEs. We finally show that the resulting expansions of the velocity u and of its first and second space derivatives converge and define continuous functions up to the boundary

  2. A Cartesian grid embedded boundary method for Poisson`s equation on irregular domains

    Energy Technology Data Exchange (ETDEWEB)

    Johansen, H. [Univ. of California, Berkeley, CA (United States). Dept. of Mechanical Engineering; Colella, P. [Lawrence Berkeley National Lab., CA (United States). Center for Computational Sciences and Engineering

    1997-01-31

    The authors present a numerical method for solving Poisson`s equation, with variable coefficients and Dirichlet boundary conditions, on two-dimensional regions. The approach uses a finite-volume discretization, which embeds the domain in a regular Cartesian grid. They treat the solution as a cell-centered quantity, even when those centers are outside the domain. Cells that contain a portion of the domain boundary use conservation differencing of second-order accurate fluxes, on each cell volume. The calculation of the boundary flux ensures that the conditioning of the matrix is relatively unaffected by small cell volumes. This allows them to use multi-grid iterations with a simple point relaxation strategy. They have combined this with an adaptive mesh refinement (AMR) procedure. They provide evidence that the algorithm is second-order accurate on various exact solutions, and compare the adaptive and non-adaptive calculations.

  3. Accuracy of the Frensley inflow boundary condition for Wigner equations in simulating resonant tunneling diodes

    International Nuclear Information System (INIS)

    In this paper, the accuracy of the Frensley inflow boundary condition of the Wigner equation is analyzed in computing the I-V characteristics of a resonant tunneling diode (RTD). It is found that the Frensley inflow boundary condition for incoming electrons holds only exactly infinite away from the active device region and its accuracy depends on the length of contacts included in the simulation. For this study, the non-equilibrium Green's function (NEGF) with a Dirichlet to Neumann mapping boundary condition is used for comparison. The I-V characteristics of the RTD are found to agree between self-consistent NEGF and Wigner methods at low bias potentials with sufficiently large GaAs contact lengths. Finally, the relation between the negative differential conductance (NDC) of the RTD and the sizes of contact and buffer in the RTD is investigated using both methods.

  4. Boundary conditions for the solution of the 3-dimensional Poisson equation in open metallic enclosures

    CERN Document Server

    Biswas, Debabrata; Kumar, Raghwendra

    2015-01-01

    Numerical solution of the Poisson equation in metallic enclosures, open at one or more ends, is important in many practical situations such as High Power Microwave (HPM) or photo-cathode devices. It requires imposition of a suitable boundary condition at the open end. In this paper, methods for solving the Poisson equation are investigated for various charge densities and aspect ratios of the open ends. It is found that a mixture of second order and third order local asymptotic boundary condition (ABC) is best suited for large aspect ratios while a proposed non-local matching method, based on the solution of the Laplace equation, scores well when the aspect ratio is near unity for all charge density variations, including ones where the centre of charge is close to an open end or the charge density is non-localized. The two methods complement each other and can be used in electrostatic calculations where the computational domain needs to be terminated at the open boundaries of the metallic enclosure.

  5. Boundary conditions for the solution of the three-dimensional Poisson equation in open metallic enclosures

    Energy Technology Data Exchange (ETDEWEB)

    Biswas, Debabrata; Singh, Gaurav; Kumar, Raghwendra [Bhabha Atomic Research Centre, Mumbai 400085 (India)

    2015-09-15

    Numerical solution of the Poisson equation in metallic enclosures, open at one or more ends, is important in many practical situations, such as high power microwave or photo-cathode devices. It requires imposition of a suitable boundary condition at the open end. In this paper, methods for solving the Poisson equation are investigated for various charge densities and aspect ratios of the open ends. It is found that a mixture of second order and third order local asymptotic boundary conditions is best suited for large aspect ratios, while a proposed non-local matching method, based on the solution of the Laplace equation, scores well when the aspect ratio is near unity for all charge density variations, including ones where the centre of charge is close to an open end or the charge density is non-localized. The two methods complement each other and can be used in electrostatic calculations where the computational domain needs to be terminated at the open boundaries of the metallic enclosure.

  6. Boundary conditions for the solution of the three-dimensional Poisson equation in open metallic enclosures

    Science.gov (United States)

    Biswas, Debabrata; Singh, Gaurav; Kumar, Raghwendra

    2015-09-01

    Numerical solution of the Poisson equation in metallic enclosures, open at one or more ends, is important in many practical situations, such as high power microwave or photo-cathode devices. It requires imposition of a suitable boundary condition at the open end. In this paper, methods for solving the Poisson equation are investigated for various charge densities and aspect ratios of the open ends. It is found that a mixture of second order and third order local asymptotic boundary conditions is best suited for large aspect ratios, while a proposed non-local matching method, based on the solution of the Laplace equation, scores well when the aspect ratio is near unity for all charge density variations, including ones where the centre of charge is close to an open end or the charge density is non-localized. The two methods complement each other and can be used in electrostatic calculations where the computational domain needs to be terminated at the open boundaries of the metallic enclosure.

  7. Hermite-DG methods for pdf equations modelling particle transport and deposition in turbulent boundary layers

    International Nuclear Information System (INIS)

    A novel methodology is presented for the numerical treatment of multi-dimensional pdf (probability density function) models used to study particle transport in turbulent boundary layers. A system of coupled Fokker–Planck type equations is constructed to describe the transport of phase-space conditioned moments of particle and fluid velocities, both streamwise and wall-normal. This system, unlike conventional moment-based transport equations, allows for an exact treatment of particle deposition at the flow boundary and provides an efficient way to handle the 5-dimensional phase-space domain. Moreover, the equations in the system are linear and can be solved in a sequential fashion; there is no closure problem to address. A hybrid Hermite-Discontinuous Galerkin scheme is developed to treat the system. The choice of Hermite basis functions in combination with an iterative scaling approach permits the efficient computation of solutions to high accuracy. Results demonstrate the effectiveness of the methodology in resolving the extreme gradients characteristic of distributions near an absorbing boundary.

  8. Numerical solution for Laplace equation with mixed boundary condition for ship problem in the sea

    Science.gov (United States)

    Silalahi, Fitriani Tupa R.; Budhi, Wono Setya; Adytia, Didit; van Groesen, E.

    2015-09-01

    One interesting phenomena is investigating the movement of ships at the sea. To start with the investigation in modelling of this problem, we will assume that the ship is only a one-dimensional object that is floating on the sea surface. Similarly, we assume that the water flow is uniform in parallel directions to the ship. Therefore, we simply use the two-dimensional Laplace equation in this problem. In the section that describes the surface of sea, Neumann boundary condition is imposed in part related to the ship and the Dirichlet boundary condition for others. Then on the other three boundaries, we imposed the Neumann boundary condition by assuming that the water does not flow on the bottom, and both end. The model is solved by numerical solution using the finite element method. Velocity potential solution on the whole domain is demonstrated as a result of the implementation of the finite element method. In this paper, we initiate an investigation with assuming that the ship is on the water so that the domain of the Laplace equation is rectangular. Then we assume the drift ship. Furthermore, we also study the dependence of width and depth of the domain to the velocity potential.

  9. Stability analysis of Boundary Layer in Poiseuille Flow Through A Modified Orr-Sommerfeld Equation

    CERN Document Server

    Monwanou, A V; Orou, J B Chabi; 10.5539/apr.v4n4p138

    2013-01-01

    For applications regarding transition prediction, wing design and control of boundary layers, the fundamental understanding of disturbance growth in the flat-plate boundary layer is an important issue. In the present work we investigate the stability of boundary layer in Poiseuille flow. We normalize pressure and time by inertial and viscous effects. The disturbances are taken to be periodic in the spanwise direction and time. We present a set of linear governing equations for the parabolic evolution of wavelike disturbances. Then, we derive modified Orr-Sommerfeld equations that can be applied in the layer. Contrary to what one might think, we find that Squire's theorem is not applicable for the boundary layer. We find also that normalization by inertial or viscous effects leads to the same order of stability or instability. For the 2D disturbances flow ($\\theta=0$), we found the same critical Reynolds number for our two normalizations. This value coincides with the one we know for neutral stability of the k...

  10. Quadrature of one integral solution of transport slowing down equation

    International Nuclear Information System (INIS)

    Analytical solution of neutron slowing down equation leads to numerical calculation of neutron flux for estimated energy values. This paper shows that Newton interpolation for obtaining tabulated values of effective cross sections and application of Gauss quadrature formula for numerical integration are valid

  11. Quadrature of an integral solution of slowing down transport equation

    International Nuclear Information System (INIS)

    Full text: Numerical solution of the slowing-down equation is derived from analytical solution for neutron flux for given values of energy. It is shown in this paper that the application of Newton interpolation for obtaining effective cross sections as well as Gauss quadrature formula for numerical integration is justified

  12. A Contracted Path Integral Solution of the Discrete Master Equation

    OpenAIRE

    Helbing, Dirk

    1998-01-01

    A new representation of the exact time dependent solution of the discrete master equation is derived. This representation can be considered as contraction of the path integral solution of Haken. It allows the calculation of the probability distribution of the occurence time for each path and is suitable as basis of new computational solution methods.

  13. ADAPTIVE INTERVAL WAVELET PRECISE INTEGRATION METHOD FOR PARTIAL DIFFERENTIAL EQUATIONS

    Institute of Scientific and Technical Information of China (English)

    MEI Shu-li; LU Qi-shao; ZHANG Sen-wen; JIN Li

    2005-01-01

    The quasi-Shannon interval wavelet is constructed based on the interpolation wavelet theory, and an adaptive precise integration method, which is based on extrapolation method is presented for nonlinear ordinary differential equations (ODEs). And then, an adaptive interval wavelet precise integration method (AIWPIM) for nonlinear partial differential equations(PDEs) is proposed. The numerical results show that the computational precision of AIWPIM is higher than that of the method constructed by combining the wavelet and the 4th Runge-Kutta method, and the computational amounts of these two methods are almost equal. For convenience, the Burgers equation is taken as an example in introducing this method, which is also valid for more general cases.

  14. Determination of a source term for a time fractional diffusion equation with an integral type over-determining condition

    Directory of Open Access Journals (Sweden)

    Timurkhan S. Aleroev

    2013-12-01

    Full Text Available We consider a linear heat equation involving a fractional derivative in time, with a nonlocal boundary condition. We determine a source term independent of the space variable, and the temperature distribution for a problem with an over-determining condition of integral type. We prove the existence and uniqueness of the solution, and its continuous dependence on the data.

  15. Numerical Approximation of Boundary Control for the Wave Equation - with Application to an Inverse Problem

    DEFF Research Database (Denmark)

    Mariegaard, Jesper Sandvig

    We consider a control problem for the wave equation: Given the initial state, find a specific boundary condition, called a control, that steers the system to a desired final state. The Hilbert uniqueness method (HUM) is a mathematical method for the solution of such control problems. It builds on...... the duality between the control system and its adjoint system, and these systems are connected via a so-called controllability operator. In this project, we are concerned with the numerical approximation of HUM control for the one-dimensional wave equation. We study two semi-discretizations of the...... wave equation: a linear finite element method (L-FEM) and a discontinuous Galerkin-FEM (DG-FEM). The controllability operator is discretized with both L-FEM and DG-FEM to obtain a HUM matrix. We show that formulating HUM in a sine basis is beneficial for several reasons: (i) separation of low and high...

  16. On energy boundary layer equations in power law non-Newtonian fluids

    Institute of Scientific and Technical Information of China (English)

    郑连存; 张欣欣

    2008-01-01

    The hear transfer mechanism and the constitutive models for energy boundary layer in power law fluids were investigated.Two energy transfer constitutive equations models were proposed based on the assumption of similarity of velocity field momentum diffusion and temperature field heat transfer.The governing systems of partial different equations were transformed into ordinary differential equations respectively by using the similarity transformation group.One model was assumed that Prandtl number is a constant,and the other model was assumed that viscosity diffusion is analogous to thermal diffusion.The solutions were presented analytically and numerically by using the Runge-Kutta formulas and shooting technique and the associated transfer characteristics were discussed.

  17. Brunet-Derrida particle systems, free boundary problems and Wiener-Hopf equations

    CERN Document Server

    Durrett, Rick

    2009-01-01

    We consider a branching-selection system in $\\rr$ with $N$ particles which give birth independently at rate 1 and where after each birth the leftmost particle is erased, keeping the number of particles constant. We show that, as $N\\to\\infty$, the empirical measure process associated to the system converges in distribution to a deterministic measure-valued process whose densities solve a free boundary integro-differential equation. We also show that this equation has a unique traveling wave solution traveling at speed $c$ or no such solution depending on whether $c>a$ or $c\\leq a$, where $a$ is the asymptotic speed of the branching random walk obtained by ignoring the removal of the leftmost particles in our process. The traveling wave solutions correspond to solutions of Wiener-Hopf equations.

  18. Difference Scheme for Hyperbolic Heat Conduction Equation with Pulsed Heating Boundary

    Institute of Scientific and Technical Information of China (English)

    LiJi; ZhangZhengfang; 等

    2000-01-01

    In this paper,the authors have analyzed the second-order hyperbolic differential equation with pulsed heating boundary,and tried to solve this kind of equation with numerical method.After analyzing the error of the difference scheme and comparing the numerical results with the theoretical results.The authors present some examples to show how the thermal wave propagates in materials.By analyzing the calculation results,the conditions to observe the thermal wave phenomena in the laboratory are conferred.Finally the heat transfter in a complex combined structure is calculated with this method.The result is quite different from the calculated result from the aprabolic heat conduction equation.

  19. Nonlinear fractional differential equations and inclusions of arbitrary order and multi-strip boundary conditions

    Directory of Open Access Journals (Sweden)

    Bashir Ahmad

    2012-06-01

    Full Text Available We study boundary value problems of nonlinear fractional differential equations and inclusions of order $q in (m-1, m]$, $m ge 2$ with multi-strip boundary conditions. Multi-strip boundary conditions may be regarded as the generalization of multi-point boundary conditions. Our problem is new in the sense that we consider a nonlocal strip condition of the form: $$ x(1=sum_{i=1}^{n-2}alpha_i int^{eta_i}_{zeta_i} x(sds, $$ which can be viewed as an extension of a multi-point nonlocal boundary condition: $$ x(1=sum_{i=1}^{n-2}alpha_i x(eta_i. $$ In fact, the strip condition corresponds to a continuous distribution of the values of the unknown function on arbitrary finite segments $(zeta_i,eta_i$ of the interval $[0,1]$ and the effect of these strips is accumulated at $x=1$. Such problems occur in the applied fields such as wave propagation and geophysics. Some new existence and uniqueness results are obtained by using a variety of fixed point theorems. Some illustrative examples are also discussed.

  20. A boundary integral method for a dynamic, transient mode I crack problem with viscoelastic cohesive zone

    KAUST Repository

    Leise, Tanya L.

    2009-08-19

    We consider the problem of the dynamic, transient propagation of a semi-infinite, mode I crack in an infinite elastic body with a nonlinear, viscoelastic cohesize zone. Our problem formulation includes boundary conditions that preclude crack face interpenetration, in contrast to the usual mode I boundary conditions that assume all unloaded crack faces are stress-free. The nonlinear viscoelastic cohesive zone behavior is motivated by dynamic fracture in brittle polymers in which crack propagation is preceeded by significant crazing in a thin region surrounding the crack tip. We present a combined analytical/numerical solution method that involves reducing the problem to a Dirichlet-to-Neumann map along the crack face plane, resulting in a differo-integral equation relating the displacement and stress along the crack faces and within the cohesive zone. © 2009 Springer Science+Business Media B.V.

  1. Integrable Rosochatius Deformations for an Integrable Couplings of CKdV Equation Hierarchy

    International Nuclear Information System (INIS)

    We propose a method to construct the integrable Rosochatius deformations for an integrable couplings equations hierarchy. As applications, the integrable Rosochatius deformations of the coupled CKdV hierarchy with self-consistent sources and its Lax representation are presented. (general)

  2. Euler integral symmetries for a deformed Heun equation and symmetries of the Painlevé PVI equation

    Science.gov (United States)

    Kazakov, A. Ya.; Slavyanov, S. Yu.

    2008-05-01

    Euler integral transformations relate solutions of ordinary linear differential equations and generate integral representations of the solutions in a number of cases or relations between solutions of constrained equations (Euler symmetries) in some other cases. These relations lead to the corresponding symmetries of the monodromy matrices. We discuss Euler symmetries in the case of the simplest Fuchsian system that is equivalent to a deformed Heun equation, which is in turn related to the Painlevé PVI equation. The existence of integral symmetries of the deformed Heun equation leads to the corresponding symmetries of the PVI equation.

  3. Integrability of a deterministic cellular automaton driven by stochastic boundaries

    Science.gov (United States)

    Prosen, Tomaž; Mejía-Monasterio, Carlos

    2016-05-01

    We propose an interacting many-body space–time-discrete Markov chain model, which is composed of an integrable deterministic and reversible cellular automaton (rule 54 of Bobenko et al 1993 Commun. Math. Phys. 158 127) on a finite one-dimensional lattice {({{{Z}}}2)}× n, and local stochastic Markov chains at the two lattice boundaries which provide chemical baths for absorbing or emitting the solitons. Ergodicity and mixing of this many-body Markov chain is proven for generic values of bath parameters, implying the existence of a unique nonequilibrium steady state. The latter is constructed exactly and explicitly in terms of a particularly simple form of matrix product ansatz which is termed a patch ansatz. This gives rise to an explicit computation of observables and k-point correlations in the steady state as well as the construction of a nontrivial set of local conservation laws. The feasibility of an exact solution for the full spectrum and eigenvectors (decay modes) of the Markov matrix is suggested as well. We conjecture that our ideas can pave the road towards a theory of integrability of boundary driven classical deterministic lattice systems.

  4. Rational first integrals of geodesic equations and generalised hidden symmetries

    CERN Document Server

    Aoki, Arata; Tomoda, Kentaro

    2016-01-01

    We discuss novel generalisations of Killing tensors, which are introduced by considering rational first integrals of geodesic equations. We introduce the notion of inconstructible generalised Killing tensors, which cannot be constructed from ordinary Killing tensors. Moreover, we introduce inconstructible rational first integrals, which are constructed from inconstructible generalised Killing tensors, and provide a method for checking the inconstructibility of a rational first integral. Using the method, we show that the rational first integral of the Collinson-O'Donnell solution is not inconstructible. We also provide several examples of metrics admitting an inconstructible rational first integral in two and four dimensions, by using the Maciejewski-Przybylska system. Furthermore, we attempt to generalise other hidden symmetries such as Killing-Yano tensors.

  5. Kwong-Wong-type integral equation on time scales

    Directory of Open Access Journals (Sweden)

    Baoguo Jia

    2011-09-01

    Full Text Available Consider the second-order nonlinear dynamic equation $$ [r(tx^Delta(ho(t]^Delta+p(tf(x(t=0, $$ where $p(t$ is the backward jump operator. We obtain a Kwong-Wong-type integral equation, that is: If $x(t$ is a nonoscillatory solution of the above equation on $[T_0,infty$, then the integral equation $$ frac{r^sigma(tx^Delta(t}{f(x^sigma(t} =P^sigma(t+int^infty_{sigma(t}frac{r^sigma(s [int^1_0f'(x_h(sdh][x^Delta(s]^2}{f(x(s f(x^sigma(s}Delta s $$ is satisfied for $tgeq T_0$, where $P^sigma(t=int^infty_{sigma(t}p(sDelta s$, and $x_h(s=x(s+hmu(sx^Delta(s$. As an application, we show that the superlinear dynamic equation $$ [r(tx^{Delta}(ho(t]^Delta+p(tf(x(t=0, $$ is oscillatory, under certain conditions.

  6. Coarse projective kMC integration: forward/reverse initial and boundary value problems

    International Nuclear Information System (INIS)

    In 'equation-free' multiscale computation a dynamic model is given at a fine, microscopic level; yet we believe that its coarse-grained, macroscopic dynamics can be described by closed equations involving only coarse variables. These variables are typically various low-order moments of the distributions evolved through the microscopic model. We consider the problem of integrating these unavailable equations by acting directly on kinetic Monte Carlo microscopic simulators, thus circumventing their derivation in closed form. In particular, we use projective multi-step integration to solve the coarse initial value problem forward in time as well as backward in time (under certain conditions). Macroscopic trajectories are thus traced back to unstable, source-type, and even sometimes saddle-like stationary points, even though the microscopic simulator only evolves forward in time. We also demonstrate the use of such projective integrators in a shooting boundary value problem formulation for the computation of 'coarse limit cycles' of the macroscopic behavior, and the approximation of their stability through estimates of the leading 'coarse Floquet multipliers'

  7. Hamiltonian time integrators for Vlasov-Maxwell equations

    International Nuclear Information System (INIS)

    Hamiltonian time integrators for the Vlasov-Maxwell equations are developed by a Hamiltonian splitting technique. The Hamiltonian functional is split into five parts, which produces five exactly solvable subsystems. Each subsystem is a Hamiltonian system equipped with the Morrison-Marsden-Weinstein Poisson bracket. Compositions of the exact solutions provide Poisson structure preserving/Hamiltonian methods of arbitrary high order for the Vlasov-Maxwell equations. They are then accurate and conservative over a long time because of the Poisson-preserving nature

  8. A FAST NUMERICAL METHOD FOR INTEGRAL EQUATIONS OF THE FIRST KIND WITH LOGARITHMIC KERNEL USING MESH GRADING

    Institute of Scientific and Technical Information of China (English)

    Qi-yuan Chen; Tao Tang; Zhen-huan Teng

    2004-01-01

    The aim of this paper is to develop a fast numerical method for two-dimensional boundary integral equations of the first kind with logarithm kernels when the boundary of the domain is smooth and closed. In this case, the use of the conventional boundary element methods gives linear systems with dense matrix. In this paper, we demonstrate that the dense matrix can be replaced by a sparse one if appropriate graded meshes are used in the quadrature rules. It will be demonstrated that this technique can increase the numerical efficiency significantly.

  9. Lagrangian structures, integrability and chaos for 3D dynamical equations

    International Nuclear Information System (INIS)

    In this paper, we consider the general setting for constructing action principles for three-dimensional first-order autonomous equations. We present the results for some integrable and non-integrable cases of the Lotka-Volterra equation, and show Lagrangian descriptions which are valid for systems satisfying Shil'nikov criteria on the existence of strange attractors, though chaotic behaviour has not been verified up to now. The Euler-Lagrange equations we get for these systems usually present 'time reparametrization' invariance, though other kinds of invariance may be found according to the kernel of the associated symplectic 2-form. The formulation of a Hamiltonian structure (Poisson brackets and Hamiltonians) for these systems from the Lagrangian viewpoint leads to a method of finding new constants of the motion starting from known ones, which is applied to some systems found in the literature known to possess a constant of the motion, to find the other and thus showing their integrability. In particular, we show that the so-called ABC system is completely integrable if it possesses one constant of the motion

  10. Superplastic Constitutive Equation Including Percentage of High-Angle Grain Boundaries as a Microstructural Parameter

    Science.gov (United States)

    Wang, K.; Liu, F. C.; Xue, P.; Wang, D.; Xiao, B. L.; Ma, Z. Y.

    2016-01-01

    Fifteen Al-Mg-Sc samples with subgrain/grain sizes in the range of 1.8 to 4.9 μm were prepared through the processing methods of friction stir processing (FSP), equal-channel-angular pressing (ECAP), rolling, annealing, and combinations of the above. The percentages of high-angle grain boundaries (HAGBs) of these fine-grained alloys were distributed from 39 to 97 pct. The samples processed through FSP had a higher percentage of HAGBs compared to other samples. Superplasticity was achieved in all fifteen samples, but the FSP samples exhibited better superplasticity than other samples because their fine equiaxed grains, which were mostly surrounded by HAGBs, were conducive to the occurrence of grain boundary sliding (GBS) during superplastic deformation. The dominant deformation mechanism was the same for all fifteen samples, i.e., GBS controlled by grain boundary diffusion. However, the subgrains were the GBS units for the rolled or ECAP samples, which contained high percentages of unrecrystallized grains, whereas the fine grains were the GBS units for the FSP samples. Superplastic data analysis revealed that the dimensionless A in the classical constitutive equation for superplasticity of fine-grained Al alloys was not a constant, but increased with an increase in the percentage of HAGBs, demonstrating that the enhanced superplastic deformation kinetics can be ascribed to the high percentage of HAGBs. A modified superplastic constitutive equation with the percentage of HAGBs as a new microstructural parameter was established.

  11. A broadband fast multipole accelerated boundary element method for the three dimensional Helmholtz equation.

    Science.gov (United States)

    Gumerov, Nail A; Duraiswami, Ramani

    2009-01-01

    The development of a fast multipole method (FMM) accelerated iterative solution of the boundary element method (BEM) for the Helmholtz equations in three dimensions is described. The FMM for the Helmholtz equation is significantly different for problems with low and high kD (where k is the wavenumber and D the domain size), and for large problems the method must be switched between levels of the hierarchy. The BEM requires several approximate computations (numerical quadrature, approximations of the boundary shapes using elements), and these errors must be balanced against approximations introduced by the FMM and the convergence criterion for iterative solution. These different errors must all be chosen in a way that, on the one hand, excess work is not done and, on the other, that the error achieved by the overall computation is acceptable. Details of translation operators for low and high kD, choice of representations, and BEM quadrature schemes, all consistent with these approximations, are described. A novel preconditioner using a low accuracy FMM accelerated solver as a right preconditioner is also described. Results of the developed solvers for large boundary value problems with 0.0001 less, similarkD less, similar500 are presented and shown to perform close to theoretical expectations. PMID:19173406

  12. Advances in the study of boundary value problems for nonlinear integrable PDEs

    International Nuclear Information System (INIS)

    In this review I summarize some of the most significant advances of the last decade in the analysis and solution of boundary value problems for integrable partial differential equations (PDEs) in two independent variables. These equations arise widely in mathematical physics, and in order to model realistic applications, it is essential to consider bounded domain and inhomogeneous boundary conditions. I focus specifically on a general and widely applicable approach, usually referred to as the unified transform or Fokas transform, that provides a substantial generalization of the classical inverse scattering transform. This approach preserves the conceptual efficiency and aesthetic appeal of the more classical transform approaches, but presents a distinctive and important difference. While the inverse scattering transform follows the ‘separation of variables’ philosophy, albeit in a nonlinear setting, the unified transform is based on the idea of synthesis, rather than separation, of variables. I will outline the main ideas in the case of linear evolution equations, and then illustrate their generalization to certain nonlinear cases of particular significance. (invited article)

  13. Advances in the study of boundary value problems for nonlinear integrable PDEs

    Science.gov (United States)

    Pelloni, Beatrice

    2015-02-01

    In this review I summarize some of the most significant advances of the last decade in the analysis and solution of boundary value problems for integrable partial differential equations (PDEs) in two independent variables. These equations arise widely in mathematical physics, and in order to model realistic applications, it is essential to consider bounded domain and inhomogeneous boundary conditions. I focus specifically on a general and widely applicable approach, usually referred to as the unified transform or Fokas transform, that provides a substantial generalization of the classical inverse scattering transform. This approach preserves the conceptual efficiency and aesthetic appeal of the more classical transform approaches, but presents a distinctive and important difference. While the inverse scattering transform follows the ‘separation of variables’ philosophy, albeit in a nonlinear setting, the unified transform is based on the idea of synthesis, rather than separation, of variables. I will outline the main ideas in the case of linear evolution equations, and then illustrate their generalization to certain nonlinear cases of particular significance.

  14. Computation of Green function of the Schroedinger-like partial differential equations by the numerical functional integration

    International Nuclear Information System (INIS)

    A new method for numerical solution of the boundary problem for Schroedinger-like partial differential equations in Rn is elaborated. The method is based on representation of multidimensional Green function in the form of multiple functional integral and on the use of approximation formulas which are constructed for such integrals. The convergence of approximations to the exact value is proved, the remainder of the formulas is estimated. Method reduces the initial differential problem to quadratures. 16 refs.; 7 tabs

  15. Introduction to stochastic analysis integrals and differential equations

    CERN Document Server

    Mackevicius, Vigirdas

    2013-01-01

    This is an introduction to stochastic integration and stochastic differential equations written in an understandable way for a wide audience, from students of mathematics to practitioners in biology, chemistry, physics, and finances. The presentation is based on the naïve stochastic integration, rather than on abstract theories of measure and stochastic processes. The proofs are rather simple for practitioners and, at the same time, rather rigorous for mathematicians. Detailed application examples in natural sciences and finance are presented. Much attention is paid to simulation diffusion pro

  16. Lotka-Volterra type equations and their explicit integration

    CERN Document Server

    Gervais, Jean-Loup; Jean-Loup Gervais; Mikhail V Saveliev

    1994-01-01

    In the present note we give an explicit integration of some two--dimensionalised Lotka--Volterra type equations associated with simple Lie algebras, other than the familiar A_n case, possessing a representation without branching. This allows us, in particular, to treat the first fundamental representations of A_r, B_r, C_r, and G_2 on the same footing.

  17. Accuracy of surrogate solutions of integral equations by feedforward networks

    Czech Academy of Sciences Publication Activity Database

    Kůrková, Věra

    Cham : Springer, 2014 - (Kóczy, L.; Pozna, C.; Kacprzyk, J.), s. 91-102 ISBN 978-3-319-03205-4. - (Studies in Computational Intelligence . 530) R&D Projects: GA ČR GAP202/11/1368; GA MŠk(CZ) LD13002 Institutional support: RVO:67985807 Keywords : surrogate modeling by neural networks * approximate solutions of integral equations * feedforward neural networks * model complexity * rates of approximation Subject RIV: IN - Informatics, Computer Science

  18. Existence in optimal control problems of certain Fredholm integral equations

    Czech Academy of Sciences Publication Activity Database

    Roubíček, Tomáš; Schmidt, W. H.

    2001-01-01

    Roč. 30, č. 3 (2001), s. 303-322. ISSN 0324-8569 R&D Projects: GA AV ČR IAA1075005; GA ČR GA201/00/0768; GA MŠk 11320007 Institutional research plan: AV0Z1075907 Keywords : optimal control * integral equation * Young measures Subject RIV: BA - General Mathematics Impact factor: 0.154, year: 2001

  19. Deterministic methods to solve the integral transport equation in neutronic

    International Nuclear Information System (INIS)

    We present a synthesis of the methods used to solve the integral transport equation in neutronic. This formulation is above all used to compute solutions in 2D in heterogeneous assemblies. Three kinds of methods are described: - the collision probability method; - the interface current method; - the current coupling collision probability method. These methods don't seem to be the most effective in 3D. (author). 9 figs

  20. DIRICHLET BOUNDARY VALUE PROBLEMS FOR SECOND-ORDER QUASI-LINEAR DIFFERENTIAL EQUATIONS WITH CHANGING SIGN NONLINEARITIES

    Institute of Scientific and Technical Information of China (English)

    2006-01-01

    This paper is concerned with the existence of positive solutions of two-point Dirichlet singular and nonsingular boundary problems for second-order quasi-linear differential equations with changing sign nonlinearities.