Abstract methods in partial differential equations

CERN Document Server

Carroll, Robert W

2012-01-01

Detailed, self-contained treatment examines modern abstract methods in partial differential equations, especially abstract evolution equations. Suitable for graduate students with some previous exposure to classical partial differential equations. 1969 edition.

Numerical Simulation of a Seaway with Breaking

Science.gov (United States)

Dommermuth, Douglas; O'Shea, Thomas; Brucker, Kyle; Wyatt, Donald

2012-11-01

The focus of this presentation is to describe the recent efforts to simulate a fully non-linear seaway with breaking by using a high-order spectral (HOS) solution of the free-surface boundary value problem to drive a three-dimensional Volume of Fluid (VOF) solution. Historically, the two main types of simulations to simulate free-surface flows are the boundary integral equations method (BIEM) and high-order spectral (HOS) methods. BIEM calculations fail at the point at which the surface impacts upon itself, if not sooner, and HOS methods can only simulate a single valued free-surface. Both also employ a single-phase approximation in which the effects of the air on the water are neglected. Due to these limitations they are unable to simulate breaking waves and air entrainment. The Volume of Fluid (VOF) method on the other hand is suitable for modeling breaking waves and air entrainment. However it is computationally intractable to generate a realistic non-linear sea-state. Here, we use the HOS solution to quickly drive, or nudge, the VOF solution into a non-linear state. The computational strategies, mathematical formulation, and numerical implementation will be discussed. The results of the VOF simulation of a seaway with breaking will also be presented, and compared to the single phase, single valued HOS results.

A generalized fractional sub-equation method for fractional differential equations with variable coefficients

International Nuclear Information System (INIS)

Tang, Bo; He, Yinnian; Wei, Leilei; Zhang, Xindong

2012-01-01

In this Letter, a generalized fractional sub-equation method is proposed for solving fractional differential equations with variable coefficients. Being concise and straightforward, this method is applied to the space–time fractional Gardner equation with variable coefficients. As a result, many exact solutions are obtained including hyperbolic function solutions, trigonometric function solutions and rational solutions. It is shown that the considered method provides a very effective, convenient and powerful mathematical tool for solving many other fractional differential equations in mathematical physics. -- Highlights: ► Study of fractional differential equations with variable coefficients plays a role in applied physical sciences. ► It is shown that the proposed algorithm is effective for solving fractional differential equations with variable coefficients. ► The obtained solutions may give insight into many considerable physical processes.

The parabolic equation method for outdoor sound propagation

DEFF Research Database (Denmark)

Arranz, Marta Galindo

The parabolic equation method is a versatile tool for outdoor sound propagation. The present study has focused on the Cranck-Nicolson type Parabolic Equation method (CNPE). Three different applications of the CNPE method have been investigated. The first two applications study variations of the g......The parabolic equation method is a versatile tool for outdoor sound propagation. The present study has focused on the Cranck-Nicolson type Parabolic Equation method (CNPE). Three different applications of the CNPE method have been investigated. The first two applications study variations...

Numerical method for the eigenvalue problem and the singular equation by using the multi-grid method and application to ordinary differential equation

International Nuclear Information System (INIS)

Kanki, Takashi; Uyama, Tadao; Tokuda, Shinji.

1995-07-01

In the numerical method to compute the matching data which are necessary for resistive MHD stability analyses, it is required to solve the eigenvalue problem and the associated singular equation. An iterative method is developed to solve the eigenvalue problem and the singular equation. In this method, the eigenvalue problem is replaced with an equivalent nonlinear equation and a singular equation is derived from Newton's method for the nonlinear equation. The multi-grid method (MGM), a high speed iterative method, can be applied to this method. The convergence of the eigenvalue and the eigenvector, and the CPU time in this method are investigated for a model equation. It is confirmed from the numerical results that this method is effective for solving the eigenvalue problem and the singular equation with numerical stability and high accuracy. It is shown by improving the MGM that the CPU time for this method is 50 times shorter than that of the direct method. (author)

Partial differential equations methods, applications and theories

CERN Document Server

Hattori, Harumi

2013-01-01

This volume is an introductory level textbook for partial differential equations (PDE's) and suitable for a one-semester undergraduate level or two-semester graduate level course in PDE's or applied mathematics. Chapters One to Five are organized according to the equations and the basic PDE's are introduced in an easy to understand manner. They include the first-order equations and the three fundamental second-order equations, i.e. the heat, wave and Laplace equations. Through these equations we learn the types of problems, how we pose the problems, and the methods of solutions such as the separation of variables and the method of characteristics. The modeling aspects are explained as well. The methods introduced in earlier chapters are developed further in Chapters Six to Twelve. They include the Fourier series, the Fourier and the Laplace transforms, and the Green's functions. The equations in higher dimensions are also discussed in detail. This volume is application-oriented and rich in examples. Going thr...

Differential equations methods and applications

CERN Document Server

Said-Houari, Belkacem

2015-01-01

This book presents a variety of techniques for solving ordinary differential equations analytically and features a wealth of examples. Focusing on the modeling of real-world phenomena, it begins with a basic introduction to differential equations, followed by linear and nonlinear first order equations and a detailed treatment of the second order linear equations. After presenting solution methods for the Laplace transform and power series, it lastly presents systems of equations and offers an introduction to the stability theory. To help readers practice the theory covered, two types of exercises are provided: those that illustrate the general theory, and others designed to expand on the text material. Detailed solutions to all the exercises are included. The book is excellently suited for use as a textbook for an undergraduate class (of all disciplines) in ordinary differential equations. .

Integral Equation Methods for Electromagnetic and Elastic Waves

CERN Document Server

Chew, Weng; Hu, Bin

2008-01-01

Integral Equation Methods for Electromagnetic and Elastic Waves is an outgrowth of several years of work. There have been no recent books on integral equation methods. There are books written on integral equations, but either they have been around for a while, or they were written by mathematicians. Much of the knowledge in integral equation methods still resides in journal papers. With this book, important relevant knowledge for integral equations are consolidated in one place and researchers need only read the pertinent chapters in this book to gain important knowledge needed for integral eq

Nodal spectrum method for solving neutron diffusion equation

International Nuclear Information System (INIS)

Sanchez, D.; Garcia, C. R.; Barros, R. C. de; Milian, D.E.

1999-01-01

Presented here is a new numerical nodal method for solving static multidimensional neutron diffusion equation in rectangular geometry. Our method is based on a spectral analysis of the nodal diffusion equations. These equations are obtained by integrating the diffusion equation in X, Y directions and then considering flat approximations for the current. These flat approximations are the only approximations that are considered in this method, as a result the numerical solutions are completely free from truncation errors. We show numerical results to illustrate the methods accuracy for coarse mesh calculations

Generalized differential transform method to differential-difference equation

International Nuclear Information System (INIS)

Zou Li; Wang Zhen; Zong Zhi

2009-01-01

In this Letter, we generalize the differential transform method to solve differential-difference equation for the first time. Two simple but typical examples are applied to illustrate the validity and the great potential of the generalized differential transform method in solving differential-difference equation. A Pade technique is also introduced and combined with GDTM in aim of extending the convergence area of presented series solutions. Comparisons are made between the results of the proposed method and exact solutions. Then we apply the differential transform method to the discrete KdV equation and the discrete mKdV equation, and successfully obtain solitary wave solutions. The results reveal that the proposed method is very effective and simple. We should point out that generalized differential transform method is also easy to be applied to other nonlinear differential-difference equation.

Advances in iterative methods for nonlinear equations

CERN Document Server

Busquier, Sonia

2016-01-01

This book focuses on the approximation of nonlinear equations using iterative methods. Nine contributions are presented on the construction and analysis of these methods, the coverage encompassing convergence, efficiency, robustness, dynamics, and applications. Many problems are stated in the form of nonlinear equations, using mathematical modeling. In particular, a wide range of problems in Applied Mathematics and in Engineering can be solved by finding the solutions to these equations. The book reveals the importance of studying convergence aspects in iterative methods and shows that selection of the most efficient and robust iterative method for a given problem is crucial to guaranteeing a good approximation. A number of sample criteria for selecting the optimal method are presented, including those regarding the order of convergence, the computational cost, and the stability, including the dynamics. This book will appeal to researchers whose field of interest is related to nonlinear problems and equations...

An irrational trial equation method and its applications

Indian Academy of Sciences (India)

equation method which is different from those direct methods. Liu's key idea is that exact solution to a differential equation can be given by solving an integration. For example, consider a differential equation of u. We always assume that its exact solution satisfies a solvable equation u = F(u). Therefore, our task is just to find.

The modified simplest equation method to look for exact solutions of nonlinear partial differential equations

OpenAIRE

Efimova, Olga Yu.

2010-01-01

The modification of simplest equation method to look for exact solutions of nonlinear partial differential equations is presented. Using this method we obtain exact solutions of generalized Korteweg-de Vries equation with cubic source and exact solutions of third-order Kudryashov-Sinelshchikov equation describing nonlinear waves in liquids with gas bubbles.

Deriving average soliton equations with a perturbative method

International Nuclear Information System (INIS)

Ballantyne, G.J.; Gough, P.T.; Taylor, D.P.

1995-01-01

The method of multiple scales is applied to periodically amplified, lossy media described by either the nonlinear Schroedinger (NLS) equation or the Korteweg--de Vries (KdV) equation. An existing result for the NLS equation, derived in the context of nonlinear optical communications, is confirmed. The method is then applied to the KdV equation and the result is confirmed numerically

Sinc-collocation method for solving the Blasius equation

International Nuclear Information System (INIS)

Parand, K.; Dehghan, Mehdi; Pirkhedri, A.

2009-01-01

Sinc-collocation method is applied for solving Blasius equation which comes from boundary layer equations. It is well known that sinc procedure converges to the solution at an exponential rate. Comparison with Howarth and Asaithambi's numerical solutions reveals that the proposed method is of high accuracy and reduces the solution of Blasius' equation to the solution of a system of algebraic equations.

A Photon Free Method to Solve Radiation Transport Equations

International Nuclear Information System (INIS)

Chang, B

2006-01-01

The multi-group discrete-ordinate equations of radiation transfer is solved for the first time by Newton's method. It is a photon free method because the photon variables are eliminated from the radiation equations to yield a N group XN direction smaller but equivalent system of equations. The smaller set of equations can be solved more efficiently than the original set of equations. Newton's method is more stable than the Semi-implicit Linear method currently used by conventional radiation codes

Modified Method of Simplest Equation Applied to the Nonlinear Schrödinger Equation

Science.gov (United States)

Vitanov, Nikolay K.; Dimitrova, Zlatinka I.

2018-03-01

We consider an extension of the methodology of the modified method of simplest equation to the case of use of two simplest equations. The extended methodology is applied for obtaining exact solutions of model nonlinear partial differential equations for deep water waves: the nonlinear Schrödinger equation. It is shown that the methodology works also for other equations of the nonlinear Schrödinger kind.

Numerov iteration method for second order integral-differential equation

International Nuclear Information System (INIS)

Zeng Fanan; Zhang Jiaju; Zhao Xuan

1987-01-01

In this paper, Numerov iterative method for second order integral-differential equation and system of equations are constructed. Numerical examples show that this method is better than direct method (Gauss elimination method) in CPU time and memoy requireing. Therefore, this method is an efficient method for solving integral-differential equation in nuclear physics

Modified Method of Simplest Equation Applied to the Nonlinear Schrödinger Equation

Directory of Open Access Journals (Sweden)

Vitanov Nikolay K.

2018-03-01

Full Text Available We consider an extension of the methodology of the modified method of simplest equation to the case of use of two simplest equations. The extended methodology is applied for obtaining exact solutions of model nonlinear partial differential equations for deep water waves: the nonlinear Schrödinger equation. It is shown that the methodology works also for other equations of the nonlinear Schrödinger kind.

Five-equation and robust three-equation methods for solution verification of large eddy simulation

Science.gov (United States)

Dutta, Rabijit; Xing, Tao

2018-02-01

This study evaluates the recently developed general framework for solution verification methods for large eddy simulation (LES) using implicitly filtered LES of periodic channel flows at friction Reynolds number of 395 on eight systematically refined grids. The seven-equation method shows that the coupling error based on Hypothesis I is much smaller as compared with the numerical and modeling errors and therefore can be neglected. The authors recommend five-equation method based on Hypothesis II, which shows a monotonic convergence behavior of the predicted numerical benchmark ( S C ), and provides realistic error estimates without the need of fixing the orders of accuracy for either numerical or modeling errors. Based on the results from seven-equation and five-equation methods, less expensive three and four-equation methods for practical LES applications were derived. It was found that the new three-equation method is robust as it can be applied to any convergence types and reasonably predict the error trends. It was also observed that the numerical and modeling errors usually have opposite signs, which suggests error cancellation play an essential role in LES. When Reynolds averaged Navier-Stokes (RANS) based error estimation method is applied, it shows significant error in the prediction of S C on coarse meshes. However, it predicts reasonable S C when the grids resolve at least 80% of the total turbulent kinetic energy.

New exact solutions to MKDV-Burgers equation and (2 + 1)-dimensional dispersive long wave equation via extended Riccati equation method

International Nuclear Information System (INIS)

Kong Cuicui; Wang Dan; Song Lina; Zhang Hongqing

2009-01-01

In this paper, with the aid of symbolic computation and a general ansaetz, we presented a new extended rational expansion method to construct new rational formal exact solutions to nonlinear partial differential equations. In order to illustrate the effectiveness of this method, we apply it to the MKDV-Burgers equation and the (2 + 1)-dimensional dispersive long wave equation, then several new kinds of exact solutions are successfully obtained by using the new ansaetz. The method can also be applied to other nonlinear partial differential equations.

Solutions of hyperbolic equations with the CIP-BS method

International Nuclear Information System (INIS)

Utsumi, Takayuki; Koga, James; Yamagiwa, Mitsuru; Yabe, Takashi; Aoki, Takayuki

2004-01-01

In this paper, we show that a new numerical method, the Constrained Interpolation Profile - Basis Set (CIP-BS) method, can solve general hyperbolic equations efficiently. This method uses a simple polynomial basis set that is easily extendable to any desired higher-order accuracy. The interpolating profile is chosen so that the subgrid scale solution approaches the local real solution owing to the constraints from the spatial derivatives of the master equations. Then, introducing scalar products, the linear and nonlinear partial differential equations are uniquely reduced to the ordinary differential equations for values and spatial derivatives at the grid points. The method gives stable, less diffusive, and accurate results. It is successfully applied to the continuity equation, the Burgers equation, the Korteweg-de Vries equation, and one-dimensional shock tube problems. (author)

Influence of Poroelasticity on the 3D Seismic Response of Complex Geological Media

Directory of Open Access Journals (Sweden)

Wuttke Frank

2017-06-01

Full Text Available Elastic wave propagation in 3D poroelastic geological media with localized heterogeneities, such as an elastic inclusion and a canyon is investigated to visualize the modification of local site responses under consideration of water saturated geomaterial. The extended computational environment herein developed is a direct Boundary Integral Equation Method (BIEM, based on the frequency-dependent fundamental solution of the governing equation in poro-visco elastodynamics. Bardet’s model is introduced in the analysis as the computationally efficient viscoelastic isomorphism to Biot’s equations of dynamic poroelasticity, thus replacing the two-phase material by a complex valued single-phase one. The potential of Bardet’s analogue is illustrated for low frequency vibrations and all simulation results demonstrate the dependency of wave field developed along the free surface on the properties of the soil material.

Fast sweeping method for the factored eikonal equation

Science.gov (United States)

Fomel, Sergey; Luo, Songting; Zhao, Hongkai

2009-09-01

We develop a fast sweeping method for the factored eikonal equation. By decomposing the solution of a general eikonal equation as the product of two factors: the first factor is the solution to a simple eikonal equation (such as distance) or a previously computed solution to an approximate eikonal equation. The second factor is a necessary modification/correction. Appropriate discretization and a fast sweeping strategy are designed for the equation of the correction part. The key idea is to enforce the causality of the original eikonal equation during the Gauss-Seidel iterations. Using extensive numerical examples we demonstrate that (1) the convergence behavior of the fast sweeping method for the factored eikonal equation is the same as for the original eikonal equation, i.e., the number of iterations for the Gauss-Seidel iterations is independent of the mesh size, (2) the numerical solution from the factored eikonal equation is more accurate than the numerical solution directly computed from the original eikonal equation, especially for point sources.

A new sub-equation method applied to obtain exact travelling wave solutions of some complex nonlinear equations

International Nuclear Information System (INIS)

Zhang Huiqun

2009-01-01

By using a new coupled Riccati equations, a direct algebraic method, which was applied to obtain exact travelling wave solutions of some complex nonlinear equations, is improved. And the exact travelling wave solutions of the complex KdV equation, Boussinesq equation and Klein-Gordon equation are investigated using the improved method. The method presented in this paper can also be applied to construct exact travelling wave solutions for other nonlinear complex equations.

Variable-mesh method of solving differential equations

Science.gov (United States)

Van Wyk, R.

1969-01-01

Multistep predictor-corrector method for numerical solution of ordinary differential equations retains high local accuracy and convergence properties. In addition, the method was developed in a form conducive to the generation of effective criteria for the selection of subsequent step sizes in step-by-step solution of differential equations.

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.

Modified Chebyshev Collocation Method for Solving Differential Equations

Directory of Open Access Journals (Sweden)

M Ziaul Arif

2015-05-01

Full Text Available This paper presents derivation of alternative numerical scheme for solving differential equations, which is modified Chebyshev (Vieta-Lucas Polynomial collocation differentiation matrices. The Scheme of modified Chebyshev (Vieta-Lucas Polynomial collocation method is applied to both Ordinary Differential Equations (ODEs and Partial Differential Equations (PDEs cases. Finally, the performance of the proposed method is compared with finite difference method and the exact solution of the example. It is shown that modified Chebyshev collocation method more effective and accurate than FDM for some example given.

Method of controlling chaos in laser equations

International Nuclear Information System (INIS)

Duong-van, M.

1993-01-01

A method of controlling chaotic to laminar flows in the Lorenz equations using fixed points dictated by minimizing the Lyapunov functional was proposed by Singer, Wang, and Bau [Phys. Rev. Lett. 66, 1123 (1991)]. Using different fixed points, we find that the solutions in a chaotic regime can also be periodic. Since the laser equations are isomorphic to the Lorenz equations we use this method to control chaos when the laser is operated over the pump threshold. Furthermore, by solving the laser equations with an occasional proportional feedback mechanism, we recover the essential laser controlling features experimentally discovered by Roy, Murphy, Jr., Maier, Gills, and Hunt [Phys. Rev. Lett. 68, 1259 (1992)

Method of controlling chaos in laser equations

Science.gov (United States)

Duong-van, Minh

1993-01-01

A method of controlling chaotic to laminar flows in the Lorenz equations using fixed points dictated by minimizing the Lyapunov functional was proposed by Singer, Wang, and Bau [Phys. Rev. Lett. 66, 1123 (1991)]. Using different fixed points, we find that the solutions in a chaotic regime can also be periodic. Since the laser equations are isomorphic to the Lorenz equations we use this method to control chaos when the laser is operated over the pump threshold. Furthermore, by solving the laser equations with an occasional proportional feedback mechanism, we recover the essential laser controlling features experimentally discovered by Roy, Murphy, Jr., Maier, Gills, and Hunt [Phys. Rev. Lett. 68, 1259 (1992)].

Lagrange-Noether method for solving second-order differential equations

Institute of Scientific and Technical Information of China (English)

Wu Hui-Bin; Wu Run-Heng

2009-01-01

The purpose of this paper is to provide a new method called the Lagrange-Noether method for solving second-order differential equations. The method is,firstly,to write the second-order differential equations completely or partially in the form of Lagrange equations,and secondly,to obtain the integrals of the equations by using the Noether theory of the Lagrange system. An example is given to illustrate the application of the result.

Equating TIMSS Mathematics Subtests with Nonlinear Equating Methods Using NEAT Design: Circle-Arc Equating Approaches

Science.gov (United States)

Ozdemir, Burhanettin

2017-01-01

The purpose of this study is to equate Trends in International Mathematics and Science Study (TIMSS) mathematics subtest scores obtained from TIMSS 2011 to scores obtained from TIMSS 2007 form with different nonlinear observed score equating methods under Non-Equivalent Anchor Test (NEAT) design where common items are used to link two or more test…

Comparison between results of solution of Burgers' equation and Laplace's equation by Galerkin and least-square finite element methods

Science.gov (United States)

Adib, Arash; Poorveis, Davood; Mehraban, Farid

2018-03-01

In this research, two equations are considered as examples of hyperbolic and elliptic equations. In addition, two finite element methods are applied for solving of these equations. The purpose of this research is the selection of suitable method for solving each of two equations. Burgers' equation is a hyperbolic equation. This equation is a pure advection (without diffusion) equation. This equation is one-dimensional and unsteady. A sudden shock wave is introduced to the model. This wave moves without deformation. In addition, Laplace's equation is an elliptical equation. This equation is steady and two-dimensional. The solution of Laplace's equation in an earth dam is considered. By solution of Laplace's equation, head pressure and the value of seepage in the directions X and Y are calculated in different points of earth dam. At the end, water table is shown in the earth dam. For Burgers' equation, least-square method can show movement of wave with oscillation but Galerkin method can not show it correctly (the best method for solving of the Burgers' equation is discrete space by least-square finite element method and discrete time by forward difference.). For Laplace's equation, Galerkin and least square methods can show water table correctly in earth dam.

Variational linear algebraic equations method

International Nuclear Information System (INIS)

Moiseiwitsch, B.L.

1982-01-01

A modification of the linear algebraic equations method is described which ensures a variational bound on the phaseshifts for potentials having a definite sign at all points. The method is illustrated by the elastic scattering of s-wave electrons by the static field of atomic hydrogen. (author)

A Line-Tau Collocation Method for Partial Differential Equations ...

African Journals Online (AJOL)

This paper deals with the numerical solution of second order linear partial differential equations with the use of the method of lines coupled with the tau collocation method. The method of lines is used to convert the partial differential equation (PDE) to a sequence of ordinary differential equations (ODEs) which is then ...

Applying homotopy analysis method for solving differential-difference equation

International Nuclear Information System (INIS)

Wang Zhen; Zou Li; Zhang Hongqing

2007-01-01

In this Letter, we apply the homotopy analysis method to solving the differential-difference equations. A simple but typical example is applied to illustrate the validity and the great potential of the generalized homotopy analysis method in solving differential-difference equation. Comparisons are made between the results of the proposed method and exact solutions. The results show that the homotopy analysis method is an attractive method in solving the differential-difference equations

On the Inclusion of Difference Equation Problems and Z Transform Methods in Sophomore Differential Equation Classes

Science.gov (United States)

Savoye, Philippe

2009-01-01

In recent years, I started covering difference equations and z transform methods in my introductory differential equations course. This allowed my students to extend the "classical" methods for (ordinary differential equation) ODE's to discrete time problems arising in many applications.

Multi-symplectic Runge-Kutta methods for nonlinear Dirac equations

International Nuclear Information System (INIS)

Hong Jialin; Li Chun

2006-01-01

In this paper, we consider the multi-symplectic Runge-Kutta (MSRK) methods applied to the nonlinear Dirac equation in relativistic quantum physics, based on a discovery of the multi-symplecticity of the equation. In particular, the conservation of energy, momentum and charge under MSRK discretizations is investigated by means of numerical experiments and numerical comparisons with non-MSRK methods. Numerical experiments presented reveal that MSRK methods applied to the nonlinear Dirac equation preserve exactly conservation laws of charge and momentum, and conserve the energy conservation in the corresponding numerical accuracy to the method utilized. It is verified numerically that MSRK methods are stable and convergent with respect to the conservation laws of energy, momentum and charge, and MSRK methods preserve not only the inner geometric structure of the equation, but also some crucial conservative properties in quantum physics. A remarkable advantage of MSRK methods applied to the nonlinear Dirac equation is the precise preservation of charge conservation law

Elliptic equation rational expansion method and new exact travelling solutions for Whitham-Broer-Kaup equations

International Nuclear Information System (INIS)

Chen Yong; Wang Qi; Li Biao

2005-01-01

Based on a new general ansatz and a general subepuation, a new general algebraic method named elliptic equation rational expansion method is devised for constructing multiple travelling wave solutions in terms of rational special function for nonlinear evolution equations (NEEs). We apply the proposed method to solve Whitham-Broer-Kaup equation and explicitly construct a series of exact solutions which include rational form solitary wave solution, rational form triangular periodic wave solutions and rational wave solutions as special cases. In addition, the links among our proposed method with the method by Fan [Chaos, Solitons and Fractals 2004;20:609], are also clarified generally

Incompressible spectral-element method: Derivation of equations

Science.gov (United States)

Deanna, Russell G.

1993-01-01

A fractional-step splitting scheme breaks the full Navier-Stokes equations into explicit and implicit portions amenable to the calculus of variations. Beginning with the functional forms of the Poisson and Helmholtz equations, we substitute finite expansion series for the dependent variables and derive the matrix equations for the unknown expansion coefficients. This method employs a new splitting scheme which differs from conventional three-step (nonlinear, pressure, viscous) schemes. The nonlinear step appears in the conventional, explicit manner, the difference occurs in the pressure step. Instead of solving for the pressure gradient using the nonlinear velocity, we add the viscous portion of the Navier-Stokes equation from the previous time step to the velocity before solving for the pressure gradient. By combining this 'predicted' pressure gradient with the nonlinear velocity in an explicit term, and the Crank-Nicholson method for the viscous terms, we develop a Helmholtz equation for the final velocity.

Weighted particle method for solving the Boltzmann equation

International Nuclear Information System (INIS)

Tohyama, M.; Suraud, E.

1990-01-01

We propose a new, deterministic, method of solution of the nuclear Boltzmann equation. In this Weighted Particle Method two-body collisions are treated by a Master equation for an occupation probability of each numerical particle. We apply the method to the quadrupole motion of 12 C. A comparison with usual stochastic methods is made. Advantages and disadvantages of the Weighted Particle Method are discussed

The modified extended Fan's sub-equation method and its application to (2 + 1)-dimensional dispersive long wave equation

International Nuclear Information System (INIS)

Yomba, Emmanuel

2005-01-01

By using a modified extended Fan's sub-equation method, we have obtained new and more general solutions including a series of non-travelling wave and coefficient function solutions namely: soliton-like solutions, triangular-like solutions, single and combined non-degenerative Jacobi elliptic wave function-like solutions for the (2 + 1)-dimensional dispersive long wave equation. The most important achievement of this method lies on the fact that, we have succeeded in one move to give all the solutions which can be previously obtained by application of at least four methods (method using Riccati equation, or first kind elliptic equation, or auxiliary ordinary equation, or generalized Riccati equation as mapping equation)

A new RBF-Trefftz meshless method for partial differential equations

International Nuclear Information System (INIS)

Cao Leilei; Zhao Ning; Qin Qinghua

2010-01-01

Based on the radial basis functions (RBF) and T-Trefftz solution, this paper presents a new meshless method for numerically solving various partial differential equation systems. First, the analog equation method (AEM) is used to convert the original patial differential equation to an equivalent Poisson's equation. Then, the radial basis functions (RBF) are employed to approxiamate the inhomogeneous term, while the homogeneous solution is obtained by linear combination of a set of T-Trefftz solutions. The present scheme, named RBF-Trefftz has the advantage over the fundamental solution (MFS) method due to the use of nonsingular T-Trefftz solution rather than singular fundamental solutions, so it does not require the artificial boundary. The application and efficiency of the proposed method are validated through several examples which include different type of differential equations, such as Laplace equation, Hellmholtz equation, convectin-diffusion equation and time-dependent equation.

A Numerical Method for Lane-Emden Equations Using Hybrid Functions and the Collocation Method

Directory of Open Access Journals (Sweden)

Changqing Yang

2012-01-01

Full Text Available A numerical method to solve Lane-Emden equations as singular initial value problems is presented in this work. This method is based on the replacement of unknown functions through a truncated series of hybrid of block-pulse functions and Chebyshev polynomials. The collocation method transforms the differential equation into a system of algebraic equations. It also has application in a wide area of differential equations. Corresponding numerical examples are presented to demonstrate the accuracy of the proposed method.

Exp-function method for solving Fisher's equation

Energy Technology Data Exchange (ETDEWEB)

Zhou, X-W [Department of Mathematics, Kunming Teacher' s College, Kunming, Yunnan 650031 (China)], E-mail: km_xwzhou@163.com

2008-02-15

There are many methods to solve Fisher's equation, but each method can only lead to a special solution. In this paper, a new method, namely the exp-function method, is employed to solve the Fisher's equation. The obtained result includes all solutions in open literature as special cases, and the generalized solution with some free parameters might imply some fascinating meanings hidden in the Fisher's equation.

A new Riccati equation rational expansion method and its application to (2 + 1)-dimensional Burgers equation

International Nuclear Information System (INIS)

Wang Qi; Chen Yong; Zhang Hongqing

2005-01-01

In this paper, we present a new Riccati equation rational expansion method to uniformly construct a series of exact solutions for nonlinear evolution equations. Compared with most existing tanh methods and other sophisticated methods, the proposed method not only recover some known solutions, but also find some new and general solutions. The solutions obtained in this paper include rational triangular periodic wave solutions, rational solitary wave solutions and rational wave solutions. The efficiency of the method can be demonstrated on (2 + 1)-dimensional Burgers equation

A new extended elliptic equation rational expansion method and its application to (2 + 1)-dimensional Burgers equation

International Nuclear Information System (INIS)

Wang Baodong; Song Lina; Zhang Hongqing

2007-01-01

In this paper, we present a new elliptic equation rational expansion method to uniformly construct a series of exact solutions for nonlinear partial differential equations. As an application of the method, we choose the (2 + 1)-dimensional Burgers equation to illustrate the method and successfully obtain some new and more general solutions

Low-mode truncation methods in the sine-Gordon equation

International Nuclear Information System (INIS)

Xiong Chuyu.

1991-01-01

In this dissertation, the author studies the chaotic and coherent motions (i.e., low-dimensional chaotic attractor) in some near integrable partial differential equations, particularly the sine-Gordon equation and the nonlinear Schroedinger equation. In order to study the motions, he uses low mode truncation methods to reduce these partial differential equations to some truncated models (low-dimensional ordinary differential equations). By applying many methods available to low-dimensional ordinary differential equations, he can understand the low-dimensional chaotic attractor of PDE's much better. However, there are two important questions one needs to answer: (1) How many modes is good enough for the low mode truncated models to capture the dynamics uniformly? (2) Is the chaotic attractor in a low mode truncated model close to the chaotic attractor in the original PDE? And how close is? He has developed two groups of powerful methods to help to answer these two questions. They are the computation methods of continuation and local bifurcation, and local Lyapunov exponents and Lyapunov exponents. Using these methods, he concludes that the 2N-nls ODE is a good model for the sine-Gordon equation and the nonlinear Schroedinger equation provided one chooses a 'good' basis and uses 'enough' modes (where 'enough' depends on the parameters of the system but is small for the parameter studied here). Therefore, one can use 2N-nls ODE to study the chaos of PDE's in more depth

Insights: A New Method to Balance Chemical Equations.

Science.gov (United States)

Garcia, Arcesio

1987-01-01

Describes a method designed to balance oxidation-reduction chemical equations. Outlines a method which is based on changes in the oxidation number that can be applied to both molecular reactions and ionic reactions. Provides examples and delineates the steps to follow for each type of equation balancing. (TW)

Integration of equations of parabolic type by the method of nets

CERN Document Server

Saul'Yev, V K; Stark, M; Ulam, S

1964-01-01

International Series of Monographs in Pure and Applied Mathematics, Volume 54: Integration of Equations of Parabolic Type by the Method of Nets deals with solving parabolic partial differential equations using the method of nets. The first part of this volume focuses on the construction of net equations, with emphasis on the stability and accuracy of the approximating net equations. The method of nets or method of finite differences (used to define the corresponding numerical method in ordinary differential equations) is one of many different approximate methods of integration of partial diff

Instability of the filtering method for Vlasov's equation

International Nuclear Information System (INIS)

Figua, H.; Bouchut, F.; Fijalkow, E.

1999-01-01

Klimas has introduced a smoothed Fourier-Fourier method. This method consists in convolving the original distribution function with a Gaussian distribution function, and, next, in solving the new system with a transformed splitting algorithm. Unfortunately, a second-order term appears in the new equation. In this work, it is studied how this term affects the numerical equation. In particular it is proven that instability occurs in the linear version of the Vlasov equation obtained by considering only free non-interacting particles. It is proved that the use of Fourier-Fourier transform is a fundamental requirement to solve this new equation. An important property is pointed out concerning the filtered distribution function in the transformed space. (K.A.)

The Extended Fractional Subequation Method for Nonlinear Fractional Differential Equations

OpenAIRE

Zhao, Jianping; Tang, Bo; Kumar, Sunil; Hou, Yanren

2012-01-01

An extended fractional subequation method is proposed for solving fractional differential equations by introducing a new general ansätz and Bäcklund transformation of the fractional Riccati equation with known solutions. Being concise and straightforward, this method is applied to the space-time fractional coupled Burgers’ equations and coupled MKdV equations. As a result, many exact solutions are obtained. It is shown that the considered method provides a very effective, convenient, and powe...

Method of ATMS operators in the formalism of Faddeev equations

International Nuclear Information System (INIS)

Zubarev, D.A.

1991-01-01

The method of ATMS operators is generalized for the case of Faddeev equations. The method to construct effective equations for both elastic scattering and scattering with rearrangement is presented. Properties to obtained equations are considered

Exact solutions for nonlinear evolution equations using Exp-function method

International Nuclear Information System (INIS)

Bekir, Ahmet; Boz, Ahmet

2008-01-01

In this Letter, the Exp-function method is used to construct solitary and soliton solutions of nonlinear evolution equations. The Klein-Gordon, Burger-Fisher and Sharma-Tasso-Olver equations are chosen to illustrate the effectiveness of the method. The method is straightforward and concise, and its applications are promising. The Exp-function method presents a wider applicability for handling nonlinear wave equations

Discrete variational derivative method a structure-preserving numerical method for partial differential equations

CERN Document Server

Furihata, Daisuke

2010-01-01

Nonlinear Partial Differential Equations (PDEs) have become increasingly important in the description of physical phenomena. Unlike Ordinary Differential Equations, PDEs can be used to effectively model multidimensional systems. The methods put forward in Discrete Variational Derivative Method concentrate on a new class of ""structure-preserving numerical equations"" which improves the qualitative behaviour of the PDE solutions and allows for stable computing. The authors have also taken care to present their methods in an accessible manner, which means that the book will be useful to engineer

Convergence of a random walk method for the Burgers equation

International Nuclear Information System (INIS)

Roberts, S.

1985-10-01

In this paper we consider a random walk algorithm for the solution of Burgers' equation. The algorithm uses the method of fractional steps. The non-linear advection term of the equation is solved by advecting ''fluid'' particles in a velocity field induced by the particles. The diffusion term of the equation is approximated by adding an appropriate random perturbation to the positions of the particles. Though the algorithm is inefficient as a method for solving Burgers' equation, it does model a similar method, the random vortex method, which has been used extensively to solve the incompressible Navier-Stokes equations. The purpose of this paper is to demonstrate the strong convergence of our random walk method and so provide a model for the proof of convergence for more complex random walk algorithms; for instance, the random vortex method without boundaries

Exp-function method for solving fractional partial differential equations.

Science.gov (United States)

Zheng, Bin

2013-01-01

We extend the Exp-function method to fractional partial differential equations in the sense of modified Riemann-Liouville derivative based on nonlinear fractional complex transformation. For illustrating the validity of this method, we apply it to the space-time fractional Fokas equation and the nonlinear fractional Sharma-Tasso-Olver (STO) equation. As a result, some new exact solutions for them are successfully established.

An Algebraic Method for Constructing Exact Solutions to Difference-Differential Equations

International Nuclear Information System (INIS)

Wang Zhen; Zhang Hongqing

2006-01-01

In this paper, we present a method to solve difference differential equation(s). As an example, we apply this method to discrete KdV equation and Ablowitz-Ladik lattice equation. As a result, many exact solutions are obtained with the help of Maple including soliton solutions presented by hyperbolic functions sinh and cosh, periodic solutions presented by sin and cos and rational solutions. This method can also be used to other nonlinear difference-differential equation(s).

New application of Exp-function method for improved Boussinesq equation

Energy Technology Data Exchange (ETDEWEB)

Abdou, M.A. [Theoretical Research Group, Physics Department, Faculty of Science, Mansoura University, 35516 Mansoura (Egypt); Department of Physics, Faculty of Education for Girls, Science Departments, King Khalid University, Bisha (Saudi Arabia)], E-mail: m_abdou_eg@yahoo.com; Soliman, A.A. [Department of Mathematics, Faculty of Education (AL-Arish) Suez Canal University, AL-Arish 45111 (Egypt); Department of Mathematics, Teacher' s College (Bisha), King Khalid University, Bisha, PO Box 551 (Saudi Arabia)], E-mail: asoliman_99@yahoo.com; El-Basyony, S.T. [Theoretical Research Group, Physics Department, Faculty of Science, Mansoura University, 35516 Mansoura (Egypt)

2007-10-01

The Exp-function method is used to obtain generalized solitary solutions and periodic solutions for nonlinear evolution equations arising in mathematical physics with the aid of symbolic computation method, namely, the improved Boussinesq equation. The method is straightforward and concise, and its applications is promising for other nonlinear evolution equations in mathematical physics.

A novel method to solve functional differential equations

International Nuclear Information System (INIS)

Tapia, V.

1990-01-01

A method to solve differential equations containing the variational operator as the derivation operation is presented. They are called variational differential equations (VDE). The solution to a VDE should be a function containing the derivatives, with respect to the base space coordinates, of the fields up to a generic order s: a s-th-order function. The variational operator doubles the order of the function on which it acts. Therefore, in order to make compatible the orders of the different terms appearing in a VDE, the solution should be a function containing the derivatives of the fields at all orders. But this takes us again back to the functional methods. In order to avoid this, one must restrict the considerations, in the case of second-order VDEs, to the space of s-th-order functions on which the variational operator acts transitively. These functions have been characterized for a one-dimensional base space for the first- and second-order cases. These functions turn out to be polynomial in the highest-order derivatives of the fields with functions of the lower-order derivatives as coefficients. Then VDEs reduce to a system of coupled partial differential equations for the coefficients above mentioned. The importance of the method lies on the fact that the solutions to VDEs are in a one-to-one correspondence with the solutions of functional differential equations. The previous method finds direct applications in quantum field theory, where the Schroedinger equation plays a central role. Since the Schroedinger equation is reduced to a system of coupled partial differential equations, this provides a nonperturbative scheme for quantum field theory. As an example, the massless scalar field is considered

The improved fractional sub-equation method and its applications to the space–time fractional differential equations in fluid mechanics

International Nuclear Information System (INIS)

Guo, Shimin; Mei, Liquan; Li, Ying; Sun, Youfa

2012-01-01

By introducing a new general ansätz, the improved fractional sub-equation method is proposed to construct analytical solutions of nonlinear evolution equations involving Jumarie's modified Riemann–Liouville derivative. By means of this method, the space–time fractional Whitham–Broer–Kaup and generalized Hirota–Satsuma coupled KdV equations are successfully solved. The obtained results show that the proposed method is quite effective, promising and convenient for solving nonlinear fractional differential equations. -- Highlights: ► We propose a novel method for nonlinear fractional differential equations. ► Two important fractional differential equations in fluid mechanics are solved successfully. ► Some new exact solutions of the fractional differential equations are obtained. ► These solutions will advance the understanding of nonlinear physical phenomena.

New exact solutions of (2 + 1)-dimensional Gardner equation via the new sine-Gordon equation expansion method

International Nuclear Information System (INIS)

Chen Yong; Yan Zhenya

2005-01-01

In this paper (2 + 1)-dimensional Gardner equation is investigated using a sine-Gordon equation expansion method, which was presented via a generalized sine-Gordon reduction equation and a new transformation. As a consequence, it is shown that the method is more powerful to obtain many types of new doubly periodic solutions of (2 + 1)-dimensional Gardner equation. In particular, solitary wave solutions are also given as simple limits of doubly periodic solutions

Variational iteration method for solving coupled-KdV equations

International Nuclear Information System (INIS)

Assas, Laila M.B.

2008-01-01

In this paper, the He's variational iteration method is applied to solve the non-linear coupled-KdV equations. This method is based on the use of Lagrange multipliers for identification of optimal value of a parameter in a functional. This technique provides a sequence of functions which converge to the exact solution of the coupled-KdV equations. This procedure is a powerful tool for solving coupled-KdV equations

Kinetic equation solution by inverse kinetic method

International Nuclear Information System (INIS)

Salas, G.

1983-01-01

We propose a computer program (CAMU) which permits to solve the inverse kinetic equation. The CAMU code is written in HPL language for a HP 982 A microcomputer with a peripheral interface HP 9876 A ''thermal graphic printer''. The CAMU code solves the inverse kinetic equation by taking as data entry the output of the ionization chambers and integrating the equation with the help of the Simpson method. With this program we calculate the evolution of the reactivity in time for a given disturbance

New method for solving three-dimensional Schroedinger equation

International Nuclear Information System (INIS)

Melezhik, V.S.

1992-01-01

A new method is developed for solving the multidimensional Schroedinger equation without the variable separation. To solve the Schroedinger equation in a multidimensional coordinate space X, a difference grid Ω i (i=1,2,...,N) for some of variables, Ω, from X={R,Ω} is introduced and the initial partial-differential equation is reduced to a system of N differential-difference equations in terms of one of the variables R. The arising multi-channel scattering (or eigenvalue) problem is solved by the algorithm based on a continuous analog of the Newton method. The approach has been successfully tested for several two-dimensional problems (scattering on a nonspherical potential well and 'dipole' scatterer, a hydrogen atom in a homogenous magnetic field) and for a three-dimensional problem of the helium-atom bound states. (author)

Solving (2 + 1)-dimensional sine-Poisson equation by a modified variable separated ordinary differential equation method

International Nuclear Information System (INIS)

Ka-Lin, Su; Yuan-Xi, Xie

2010-01-01

By introducing a more general auxiliary ordinary differential equation (ODE), a modified variable separated ordinary differential equation method is presented for solving the (2 + 1)-dimensional sine-Poisson equation. As a result, many explicit and exact solutions of the (2 + 1)-dimensional sine-Poisson equation are derived in a simple manner by this technique. (general)

A method for solving the KDV equation and some numerical experiments

International Nuclear Information System (INIS)

Chang Jinjiang.

1993-01-01

In this paper, by means of difference method for discretization of space partial derivatives of KDV equation, an initial value problem in ordinary differential equations of large dimensions is produced. By using this ordinary differential equations the existence and the uniqueness of the solution of the KDV equation and the conservation of scheme are proved. This ordinary differential equation can be solved by using implicit Runge-Kutta methods, so a new method for finding the numerical solution of the KDV equation is presented. Numerical experiments not only describe in detail the procedure of two solitons collision, soliton reflex and soliton produce, but also show that this method is very effective. (author). 7 refs, 3 figs

A functional-analytic method for the study of difference equations

Directory of Open Access Journals (Sweden)

Siafarikas Panayiotis D

2004-01-01

Full Text Available We will give the generalization of a recently developed functional-analytic method for studying linear and nonlinear, ordinary and partial, difference equations in the and spaces, p∈ℕ, . The method will be illustrated by use of two examples concerning a nonlinear ordinary difference equation known as the Putnam equation, and a linear partial difference equation of three variables describing the discrete Newton law of cooling in three dimensions.

Modified harmonic balance method for the solution of nonlinear jerk equations

Science.gov (United States)

Rahman, M. Saifur; Hasan, A. S. M. Z.

2018-03-01

In this paper, a second approximate solution of nonlinear jerk equations (third order differential equation) can be obtained by using modified harmonic balance method. The method is simpler and easier to carry out the solution of nonlinear differential equations due to less number of nonlinear equations are required to solve than the classical harmonic balance method. The results obtained from this method are compared with those obtained from the other existing analytical methods that are available in the literature and the numerical method. The solution shows a good agreement with the numerical solution as well as the analytical methods of the available literature.

A generalized auxiliary equation method and its application to nonlinear Klein-Gordon and generalized nonlinear Camassa-Holm equations

International Nuclear Information System (INIS)

Yomba, Emmanuel

2008-01-01

With the aid of symbolic computation, a generalized auxiliary equation method is proposed to construct more general exact solutions to two types of NLPDEs. First, we present new family of solutions to a nonlinear Klein-Gordon equation, by using this auxiliary equation method including a new first-order nonlinear ODE with six-degree nonlinear term proposed by Sirendaoreji. Then, we apply an indirect F-function method very close to the F-expansion method to solve the generalized Camassa-Holm equation with fully nonlinear dispersion and fully nonlinear convection C(l,n,p). Taking advantage of the new first-order nonlinear ODE with six degree nonlinear term, this indirect F-function method is used to map the solutions of C(l,n,p) equations to those of that nonlinear ODE. As a result, we can successfully obtain in a unified way, many exact solutions

Solving the discrete KdV equation with homotopy analysis method

International Nuclear Information System (INIS)

Zou, L.; Zong, Z.; Wang, Z.; He, L.

2007-01-01

In this Letter, we apply the homotopy analysis method to differential-difference equations. We take the discrete KdV equation as an example, and successfully obtain double periodic wave solutions and solitary wave solutions. It illustrates the validity and the great potential of the homotopy analysis method in solving discrete KdV equation. Comparisons are made between the results of the proposed method and exact solutions. The results reveal that the proposed method is very effective and convenient

A bivariate Chebyshev spectral collocation quasilinearization method for nonlinear evolution parabolic equations.

Science.gov (United States)

Motsa, S S; Magagula, V M; Sibanda, P

2014-01-01

This paper presents a new method for solving higher order nonlinear evolution partial differential equations (NPDEs). The method combines quasilinearisation, the Chebyshev spectral collocation method, and bivariate Lagrange interpolation. In this paper, we use the method to solve several nonlinear evolution equations, such as the modified KdV-Burgers equation, highly nonlinear modified KdV equation, Fisher's equation, Burgers-Fisher equation, Burgers-Huxley equation, and the Fitzhugh-Nagumo equation. The results are compared with known exact analytical solutions from literature to confirm accuracy, convergence, and effectiveness of the method. There is congruence between the numerical results and the exact solutions to a high order of accuracy. Tables were generated to present the order of accuracy of the method; convergence graphs to verify convergence of the method and error graphs are presented to show the excellent agreement between the results from this study and the known results from literature.

A Bivariate Chebyshev Spectral Collocation Quasilinearization Method for Nonlinear Evolution Parabolic Equations

Directory of Open Access Journals (Sweden)

S. S. Motsa

2014-01-01

Full Text Available This paper presents a new method for solving higher order nonlinear evolution partial differential equations (NPDEs. The method combines quasilinearisation, the Chebyshev spectral collocation method, and bivariate Lagrange interpolation. In this paper, we use the method to solve several nonlinear evolution equations, such as the modified KdV-Burgers equation, highly nonlinear modified KdV equation, Fisher's equation, Burgers-Fisher equation, Burgers-Huxley equation, and the Fitzhugh-Nagumo equation. The results are compared with known exact analytical solutions from literature to confirm accuracy, convergence, and effectiveness of the method. There is congruence between the numerical results and the exact solutions to a high order of accuracy. Tables were generated to present the order of accuracy of the method; convergence graphs to verify convergence of the method and error graphs are presented to show the excellent agreement between the results from this study and the known results from literature.

Symmetrized neutron transport equation and the fast Fourier transform method

International Nuclear Information System (INIS)

Sinh, N.Q.; Kisynski, J.; Mika, J.

1978-01-01

The differential equation obtained from the neutron transport equation by the application of the source iteration method in two-dimensional rectangular geometry is transformed into a symmetrized form with respect to one of the angular variables. The discretization of the symmetrized equation leads to finite difference equations based on the five-point scheme and solved by use of the fast Fourier transform method. Possible advantages of the approach are shown on test calculations

Dirac equation in low dimensions: The factorization method

Energy Technology Data Exchange (ETDEWEB)

Sánchez-Monroy, J.A., E-mail: antosan@if.usp.br [Instituto de Física, Universidade de São Paulo, 05508-090, São Paulo, SP (Brazil); Quimbay, C.J., E-mail: cjquimbayh@unal.edu.co [Departamento de Física, Universidad Nacional de Colombia, Bogotá, D. C. (Colombia); CIF, Bogotá (Colombia)

2014-11-15

We present a general approach to solve the (1+1) and (2+1)-dimensional Dirac equations in the presence of static scalar, pseudoscalar and gauge potentials, for the case in which the potentials have the same functional form and thus the factorization method can be applied. We show that the presence of electric potentials in the Dirac equation leads to two Klein–Gordon equations including an energy-dependent potential. We then generalize the factorization method for the case of energy-dependent Hamiltonians. Additionally, the shape invariance is generalized for a specific class of energy-dependent Hamiltonians. We also present a condition for the absence of the Klein paradox (stability of the Dirac sea), showing how Dirac particles in low dimensions can be confined for a wide family of potentials. - Highlights: • The low-dimensional Dirac equation in the presence of static potentials is solved. • The factorization method is generalized for energy-dependent Hamiltonians. • The shape invariance is generalized for energy-dependent Hamiltonians. • The stability of the Dirac sea is related to the existence of supersymmetric partner Hamiltonians.

The extended Fan's sub-equation method and its application to KdV-MKdV, BKK and variant Boussinesq equations

International Nuclear Information System (INIS)

Yomba, Emmanuel

2005-01-01

An extended Fan's sub-equation method is used for constructing exact travelling wave solutions of nonlinear partial differential equations (NLPDEs). The key idea of this method is to take full advantage of the general elliptic equation involving five parameters which has more new solutions and whose degeneracies can lead to special sub-equations involving three parameters. More new solutions are obtained for KdV-MKdV, Broer-Kaup-Kupershmidt (BKK) and variant Boussinesq equations. Then we present a technique which not only gives us a clear relation among this general elliptic equation and other sub-equations involving three parameters (Riccati equation, first kind elliptic equation, auxiliary ordinary equation, generalized Riccati equation and so on), but also provides an approach to construct new exact solutions to NLPDEs

A method of solving simple harmonic oscillator Schroedinger equation

Science.gov (United States)

Maury, Juan Carlos F.

1995-01-01

A usual step in solving totally Schrodinger equation is to try first the case when dimensionless position independent variable w is large. In this case the Harmonic Oscillator equation takes the form (d(exp 2)/dw(exp 2) - w(exp 2))F = 0, and following W.K.B. method, it gives the intermediate corresponding solution F = exp(-w(exp 2)/2), which actually satisfies exactly another equation, (d(exp 2)/dw(exp 2) + 1 - w(exp 2))F = 0. We apply a different method, useful in anharmonic oscillator equations, similar to that of Rampal and Datta, and although it is slightly more complicated however it is also more general and systematic.

Discontinuous Galerkin finite element methods for hyperbolic differential equations

NARCIS (Netherlands)

van der Vegt, Jacobus J.W.; van der Ven, H.; Boelens, O.J.; Boelens, O.J.; Toro, E.F.

2002-01-01

In this paper a suryey is given of the important steps in the development of discontinuous Galerkin finite element methods for hyperbolic partial differential equations. Special attention is paid to the application of the discontinuous Galerkin method to the solution of the Euler equations of gas

Differential and difference equations a comparison of methods of solution

CERN Document Server

Maximon, Leonard C

2016-01-01

This book, intended for researchers and graduate students in physics, applied mathematics and engineering, presents a detailed comparison of the important methods of solution for linear differential and difference equations - variation of constants, reduction of order, Laplace transforms and generating functions - bringing out the similarities as well as the significant differences in the respective analyses. Equations of arbitrary order are studied, followed by a detailed analysis for equations of first and second order. Equations with polynomial coefficients are considered and explicit solutions for equations with linear coefficients are given, showing significant differences in the functional form of solutions of differential equations from those of difference equations. An alternative method of solution involving transformation of both the dependent and independent variables is given for both differential and difference equations. A comprehensive, detailed treatment of Green’s functions and the associat...

A second order discontinuous Galerkin fast sweeping method for Eikonal equations

Science.gov (United States)

Li, Fengyan; Shu, Chi-Wang; Zhang, Yong-Tao; Zhao, Hongkai

2008-09-01

In this paper, we construct a second order fast sweeping method with a discontinuous Galerkin (DG) local solver for computing viscosity solutions of a class of static Hamilton-Jacobi equations, namely the Eikonal equations. Our piecewise linear DG local solver is built on a DG method developed recently [Y. Cheng, C.-W. Shu, A discontinuous Galerkin finite element method for directly solving the Hamilton-Jacobi equations, Journal of Computational Physics 223 (2007) 398-415] for the time-dependent Hamilton-Jacobi equations. The causality property of Eikonal equations is incorporated into the design of this solver. The resulting local nonlinear system in the Gauss-Seidel iterations is a simple quadratic system and can be solved explicitly. The compactness of the DG method and the fast sweeping strategy lead to fast convergence of the new scheme for Eikonal equations. Extensive numerical examples verify efficiency, convergence and second order accuracy of the proposed method.

Exact solution of some linear matrix equations using algebraic methods

Science.gov (United States)

Djaferis, T. E.; Mitter, S. K.

1977-01-01

A study is done of solution methods for Linear Matrix Equations including Lyapunov's equation, using methods of modern algebra. The emphasis is on the use of finite algebraic procedures which are easily implemented on a digital computer and which lead to an explicit solution to the problem. The action f sub BA is introduced a Basic Lemma is proven. The equation PA + BP = -C as well as the Lyapunov equation are analyzed. Algorithms are given for the solution of the Lyapunov and comment is given on its arithmetic complexity. The equation P - A'PA = Q is studied and numerical examples are given.

The modified simple equation method for solving some fractional ...

Indian Academy of Sciences (India)

... and processes in various areas of natural science. Thus, many effective and powerful methods have been established and improved. In this study, we establish exact solutions of the time fractional biological population model equation and nonlinearfractional Klein–Gordon equation by using the modified simple equation ...

Taylor's series method for solving the nonlinear point kinetics equations

International Nuclear Information System (INIS)

Nahla, Abdallah A.

2011-01-01

Highlights: → Taylor's series method for nonlinear point kinetics equations is applied. → The general order of derivatives are derived for this system. → Stability of Taylor's series method is studied. → Taylor's series method is A-stable for negative reactivity. → Taylor's series method is an accurate computational technique. - Abstract: Taylor's series method for solving the point reactor kinetics equations with multi-group of delayed neutrons in the presence of Newtonian temperature feedback reactivity is applied and programmed by FORTRAN. This system is the couples of the stiff nonlinear ordinary differential equations. This numerical method is based on the different order derivatives of the neutron density, the precursor concentrations of i-group of delayed neutrons and the reactivity. The r th order of derivatives are derived. The stability of Taylor's series method is discussed. Three sets of applications: step, ramp and temperature feedback reactivities are computed. Taylor's series method is an accurate computational technique and stable for negative step, negative ramp and temperature feedback reactivities. This method is useful than the traditional methods for solving the nonlinear point kinetics equations.

Structural equation modeling methods and applications

CERN Document Server

Wang, Jichuan

2012-01-01

A reference guide for applications of SEM using Mplus Structural Equation Modeling: Applications Using Mplus is intended as both a teaching resource and a reference guide. Written in non-mathematical terms, this book focuses on the conceptual and practical aspects of Structural Equation Modeling (SEM). Basic concepts and examples of various SEM models are demonstrated along with recently developed advanced methods, such as mixture modeling and model-based power analysis and sample size estimate for SEM. The statistical modeling program, Mplus, is also featured and provides researchers with a

Solving the Schroedinger equation using the finite difference time domain method

International Nuclear Information System (INIS)

Sudiarta, I Wayan; Geldart, D J Wallace

2007-01-01

In this paper, we solve the Schroedinger equation using the finite difference time domain (FDTD) method to determine energies and eigenfunctions. In order to apply the FDTD method, the Schroedinger equation is first transformed into a diffusion equation by the imaginary time transformation. The resulting time-domain diffusion equation is then solved numerically by the FDTD method. The theory and an algorithm are provided for the procedure. Numerical results are given for illustrative examples in one, two and three dimensions. It is shown that the FDTD method accurately determines eigenfunctions and energies of these systems

Q-Step methods for Newton-Jacobi operator equation | Uwasmusi ...

African Journals Online (AJOL)

The paper considers the Newton-Jacobi operator equation for the solution of nonlinear systems of equations. Special attention is paid to the computational part of this method with particular reference to the q-step methods. Journal of the Nigerian Association of Mathematical Physics Vol. 8 2004: pp. 237-241 ...

Introduction to numerical methods for time dependent differential equations

CERN Document Server

Kreiss, Heinz-Otto

2014-01-01

Introduces both the fundamentals of time dependent differential equations and their numerical solutions Introduction to Numerical Methods for Time Dependent Differential Equations delves into the underlying mathematical theory needed to solve time dependent differential equations numerically. Written as a self-contained introduction, the book is divided into two parts to emphasize both ordinary differential equations (ODEs) and partial differential equations (PDEs). Beginning with ODEs and their approximations, the authors provide a crucial presentation of fundamental notions, such as the t

Methods of mathematical modelling continuous systems and differential equations

CERN Document Server

Witelski, Thomas

2015-01-01

This book presents mathematical modelling and the integrated process of formulating sets of equations to describe real-world problems. It describes methods for obtaining solutions of challenging differential equations stemming from problems in areas such as chemical reactions, population dynamics, mechanical systems, and fluid mechanics. Chapters 1 to 4 cover essential topics in ordinary differential equations, transport equations and the calculus of variations that are important for formulating models. Chapters 5 to 11 then develop more advanced techniques including similarity solutions, matched asymptotic expansions, multiple scale analysis, long-wave models, and fast/slow dynamical systems. Methods of Mathematical Modelling will be useful for advanced undergraduate or beginning graduate students in applied mathematics, engineering and other applied sciences.

Green's function method for perturbed Korteweg-de Vries equation

International Nuclear Information System (INIS)

Cai Hao; Huang Nianning

2003-01-01

The x-derivatives of squared Jost solution are the eigenfunctions with the zero eigenvalue of the linearized equation derived from the perturbed Korteweg-de Vries equation. A method similar to Green's function formalism is introduced to show the completeness of the squared Jost solutions in multi-soliton cases. It is not related to Lax equations directly, and thus it is beneficial to deal with the nonlinear equations with complicated Lax pair

A General Linear Method for Equating with Small Samples

Science.gov (United States)

Albano, Anthony D.

2015-01-01

Research on equating with small samples has shown that methods with stronger assumptions and fewer statistical estimates can lead to decreased error in the estimated equating function. This article introduces a new approach to linear observed-score equating, one which provides flexible control over how form difficulty is assumed versus estimated…

Analytic method for solitary solutions of some partial differential equations

International Nuclear Information System (INIS)

Ugurlu, Yavuz; Kaya, Dogan

2007-01-01

In this Letter by considering an improved tanh function method, we found some exact solutions of the clannish random walker's parabolic equation, the modified Korteweg-de Vries (KdV) equation, and the Sharma-Tasso-Olver (STO) equation with its fission and fusion, the Jaulent-Miodek equation

The Integral Equation Method and the Neumann Problem for the Poisson Equation on NTA Domains

Czech Academy of Sciences Publication Activity Database

Medková, Dagmar

2009-01-01

Roč. 63, č. 21 (2009), s. 227-247 ISSN 0378-620X Institutional research plan: CEZ:AV0Z10190503 Keywords : Poisson equation * Neumann problem * integral equation method Subject RIV: BA - General Mathematics Impact factor: 0.477, year: 2009

Workshop on Numerical Methods for Ordinary Differential Equations

CERN Document Server

Gear, Charles; Russo, Elvira

1989-01-01

Developments in numerical initial value ode methods were the focal topic of the meeting at L'Aquila which explord the connections between the classical background and new research areas such as differental-algebraic equations, delay integral and integro-differential equations, stability properties, continuous extensions (interpolants for Runge-Kutta methods and their applications, effective stepsize control, parallel algorithms for small- and large-scale parallel architectures). The resulting proceedings address many of these topics in both research and survey papers.

Solving the interval type-2 fuzzy polynomial equation using the ranking method

Science.gov (United States)

Rahman, Nurhakimah Ab.; Abdullah, Lazim

2014-07-01

Polynomial equations with trapezoidal and triangular fuzzy numbers have attracted some interest among researchers in mathematics, engineering and social sciences. There are some methods that have been developed in order to solve these equations. In this study we are interested in introducing the interval type-2 fuzzy polynomial equation and solving it using the ranking method of fuzzy numbers. The ranking method concept was firstly proposed to find real roots of fuzzy polynomial equation. Therefore, the ranking method is applied to find real roots of the interval type-2 fuzzy polynomial equation. We transform the interval type-2 fuzzy polynomial equation to a system of crisp interval type-2 fuzzy polynomial equation. This transformation is performed using the ranking method of fuzzy numbers based on three parameters, namely value, ambiguity and fuzziness. Finally, we illustrate our approach by numerical example.

A New Pseudoinverse Matrix Method For Balancing Chemical Equations And Their Stability

International Nuclear Information System (INIS)

Risteski, Ice B.

2008-01-01

In this work is given a new pseudoniverse matrix method for balancing chemical equations. Here offered method is founded on virtue of the solution of a Diophantine matrix equation by using of a Moore-Penrose pseudoinverse matrix. The method has been tested on several typical chemical equations and found to be very successful for the all equations in our extensive balancing research. This method, which works successfully without any limitations, also has the capability to determine the feasibility of a new chemical reaction, and if it is feasible, then it will balance the equation. Chemical equations treated here possess atoms with fractional oxidation numbers. Also, in the present work are introduced necessary and sufficient criteria for stability of chemical equations over stability of their extended matrices

Solution methods for large systems of linear equations in BACCHUS

International Nuclear Information System (INIS)

Homann, C.; Dorr, B.

1993-05-01

The computer programme BACCHUS is used to describe steady state and transient thermal-hydraulic behaviour of a coolant in a fuel element with intact geometry in a fast breeder reactor. In such computer programmes generally large systems of linear equations with sparse matrices of coefficients, resulting from discretization of coolant conservation equations, must be solved thousands of times giving rise to large demands of main storage and CPU time. Direct and iterative solution methods of the systems of linear equations, available in BACCHUS, are described, giving theoretical details and experience with their use in the programme. Besides use of a method of lines, a Runge-Kutta-method, for solution of the partial differential equation is outlined. (orig.) [de

Spline Collocation Method for Nonlinear Multi-Term Fractional Differential Equation

OpenAIRE

Choe, Hui-Chol; Kang, Yong-Suk

2013-01-01

We study an approximation method to solve nonlinear multi-term fractional differential equations with initial conditions or boundary conditions. First, we transform the nonlinear multi-term fractional differential equations with initial conditions and boundary conditions to nonlinear fractional integral equations and consider the relations between them. We present a Spline Collocation Method and prove the existence, uniqueness and convergence of approximate solution as well as error estimatio...

Numerical methods for stochastic partial differential equations with white noise

CERN Document Server

Zhang, Zhongqiang

2017-01-01

This book covers numerical methods for stochastic partial differential equations with white noise using the framework of Wong-Zakai approximation. The book begins with some motivational and background material in the introductory chapters and is divided into three parts. Part I covers numerical stochastic ordinary differential equations. Here the authors start with numerical methods for SDEs with delay using the Wong-Zakai approximation and finite difference in time. Part II covers temporal white noise. Here the authors consider SPDEs as PDEs driven by white noise, where discretization of white noise (Brownian motion) leads to PDEs with smooth noise, which can then be treated by numerical methods for PDEs. In this part, recursive algorithms based on Wiener chaos expansion and stochastic collocation methods are presented for linear stochastic advection-diffusion-reaction equations. In addition, stochastic Euler equations are exploited as an application of stochastic collocation methods, where a numerical compa...

KRYSI, Ordinary Differential Equations Solver with Sdirk Krylov Method

International Nuclear Information System (INIS)

Hindmarsh, A.C.; Norsett, S.P.

2001-01-01

1 - Description of program or function: KRYSI is a set of FORTRAN subroutines for solving ordinary differential equations initial value problems. It is suitable for both stiff and non-stiff systems. When solving the implicit stage equations in the stiff case, KRYSI uses a Krylov subspace iteration method called the SPIGMR (Scaled Preconditioned Incomplete Generalized Minimum Residual) method. No explicit Jacobian storage is required, except where used in pre- conditioning. A demonstration problem is included with a description of two pre-conditioners that are natural for its solution by KRYSI. 2 - Method of solution: KRYSI uses a three-stage, third-order singly diagonally implicit Runge-Kutta (SDIRK) method. In the stiff case, a preconditioned Krylov subspace iteration within a (so-called) inexact Newton iteration is used to solve the system of nonlinear algebraic equations

Analytic method for solitary solutions of some partial differential equations

Energy Technology Data Exchange (ETDEWEB)

Ugurlu, Yavuz [Firat University, Department of Mathematics, 23119 Elazig (Turkey); Kaya, Dogan [Firat University, Department of Mathematics, 23119 Elazig (Turkey)], E-mail: dkaya@firat.edu.tr

2007-10-22

In this Letter by considering an improved tanh function method, we found some exact solutions of the clannish random walker's parabolic equation, the modified Korteweg-de Vries (KdV) equation, and the Sharma-Tasso-Olver (STO) equation with its fission and fusion, the Jaulent-Miodek equation.

Application of Exp-function method for (2 + 1)-dimensional nonlinear evolution equations

International Nuclear Information System (INIS)

Bekir, Ahmet; Boz, Ahmet

2009-01-01

In this paper, the Exp-function method is used to construct solitary and soliton solutions of (2 + 1)-dimensional nonlinear evolution equations. (2 + 1)-dimensional breaking soliton (Calogero) equation, modified Zakharov-Kuznetsov and Konopelchenko-Dubrovsky equations are chosen to illustrate the effectiveness of the method. The method is straightforward and concise, and its applications are promising. The Exp-function method presents a wider applicability for handling nonlinear wave equations.

Partial differential equations with numerical methods

CERN Document Server

Larsson, Stig

2003-01-01

The book is suitable for advanced undergraduate and beginning graduate students of applied mathematics and engineering. The main theme is the integration of the theory of linear PDEs and the numerical solution of such equations. For each type of PDE, elliptic, parabolic, and hyperbolic, the text contains one chapter on the mathematical theory of the differential equation, followed by one chapter on finite difference methods and one on finite element methods. As preparation, the two-point boundary value problem and the initial-value problem for ODEs are discussed in separate chapters. There is also one chapter on the elliptic eigenvalue problem and eigenfunction expansion. The presentation does not presume a deep knowledge of mathematical and functional analysis. Some background on linear functional analysis and Sobolev spaces, and also on numerical linear algebra, is reviewed in two appendices.

Navier-Stokes equations by the finite element method

International Nuclear Information System (INIS)

Portella, P.E.

1984-01-01

A computer program to solve the Navier-Stokes equations by using the Finite Element Method is implemented. The solutions variables investigated are stream-function/vorticity in the steady case and velocity/pressure in the steady state and transient cases. For steady state flow the equations are solved simultaneously by the Newton-Raphson method. For the time dependent formulation, a fractional step method is employed to discretize in time and artificial viscosity is used to preclude spurious oscilations in the solution. The element used is the three node triangle. Some numerical examples are presented and comparisons are made with applications already existent. (Author) [pt

Optimal analytic method for the nonlinear Hasegawa-Mima equation

Science.gov (United States)

Baxter, Mathew; Van Gorder, Robert A.; Vajravelu, Kuppalapalle

2014-05-01

The Hasegawa-Mima equation is a nonlinear partial differential equation that describes the electric potential due to a drift wave in a plasma. In the present paper, we apply the method of homotopy analysis to a slightly more general Hasegawa-Mima equation, which accounts for hyper-viscous damping or viscous dissipation. First, we outline the method for the general initial/boundary value problem over a compact rectangular spatial domain. We use a two-stage method, where both the convergence control parameter and the auxiliary linear operator are optimally selected to minimize the residual error due to the approximation. To do the latter, we consider a family of operators parameterized by a constant which gives the decay rate of the solutions. After outlining the general method, we consider a number of concrete examples in order to demonstrate the utility of this approach. The results enable us to study properties of the initial/boundary value problem for the generalized Hasegawa-Mima equation. In several cases considered, we are able to obtain solutions with extremely small residual errors after relatively few iterations are computed (residual errors on the order of 10-15 are found in multiple cases after only three iterations). The results demonstrate that selecting a parameterized auxiliary linear operator can be extremely useful for minimizing residual errors when used concurrently with the optimal homotopy analysis method, suggesting that this approach can prove useful for a number of nonlinear partial differential equations arising in physics and nonlinear mechanics.

Construct solitary solutions of discrete hybrid equation by Adomian Decomposition Method

International Nuclear Information System (INIS)

Wang Zhen; Zhang Hongqing

2009-01-01

In this paper, we apply the Adomian Decomposition Method to solving the differential-difference equations. A typical example is applied to illustrate the validity and the great potential of the Adomian Decomposition Method in solving differential-difference equation. Kink shaped solitary solution and Bell shaped solitary solution are presented. Comparisons are made between the results of the proposed method and exact solutions. The results show that the Adomian Decomposition Method is an attractive method in solving the differential-difference equations.

Solution of the Schroedinger equation by a spectral method

International Nuclear Information System (INIS)

Feit, M.D.; Fleck, J.A. Jr.; Steiger, A.

1982-01-01

A new computational method for determining the eigenvalues and eigenfunctions of the Schroedinger equation is described. Conventional methods for solving this problem rely on diagonalization of a Hamiltonian matrix or iterative numerical solutions of a time independent wave equation. The new method, in contrast, is based on the spectral properties of solutions to the time-dependent Schroedinger equation. The method requires the computation of a correlation function from a numerical solution psi(r, t). Fourier analysis of this correlation function reveals a set of resonant peaks that correspond to the stationary states of the system. Analysis of the location of these peaks reveals the eigenvalues with high accuracy. Additional Fourier transforms of psi(r, t) with respect to time generate the eigenfunctions. The effectiveness of the method is demonstrated for a one-dimensional asymmetric double well potential and for the two-dimensional Henon--Heiles potential

Various Newton-type iterative methods for solving nonlinear equations

Directory of Open Access Journals (Sweden)

Manoj Kumar

2013-10-01

Full Text Available The aim of the present paper is to introduce and investigate new ninth and seventh order convergent Newton-type iterative methods for solving nonlinear equations. The ninth order convergent Newton-type iterative method is made derivative free to obtain seventh-order convergent Newton-type iterative method. These new with and without derivative methods have efficiency indices 1.5518 and 1.6266, respectively. The error equations are used to establish the order of convergence of these proposed iterative methods. Finally, various numerical comparisons are implemented by MATLAB to demonstrate the performance of the developed methods.

Approximate solution fuzzy pantograph equation by using homotopy perturbation method

Science.gov (United States)

Jameel, A. F.; Saaban, A.; Ahadkulov, H.; Alipiah, F. M.

2017-09-01

In this paper, Homotopy Perturbation Method (HPM) is modified and formulated to find the approximate solution for its employment to solve (FDDEs) involving a fuzzy pantograph equation. The solution that can be obtained by using HPM is in the form of infinite series that converge to the actual solution of the FDDE and this is one of the benefits of this method In addition, it can be used for solving high order fuzzy delay differential equations directly without reduction to a first order system. Moreover, the accuracy of HPM can be detected without needing the exact solution. The HPM is studied for fuzzy initial value problems involving pantograph equation. Using the properties of fuzzy set theory, we reformulate the standard approximate method of HPM and obtain the approximate solutions. The effectiveness of the proposed method is demonstrated for third order fuzzy pantograph equation.

Mathematical and computational methods for semiclassical Schrödinger equations

KAUST Repository

Jin, Shi

2011-04-28

We consider time-dependent (linear and nonlinear) Schrödinger equations in a semiclassical scaling. These equations form a canonical class of (nonlinear) dispersive models whose solutions exhibit high-frequency oscillations. The design of efficient numerical methods which produce an accurate approximation of the solutions, or at least of the associated physical observables, is a formidable mathematical challenge. In this article we shall review the basic analytical methods for dealing with such equations, including WKB asymptotics, Wigner measure techniques and Gaussian beams. Moreover, we shall give an overview of the current state of the art of numerical methods (most of which are based on the described analytical techniques) for the Schrödinger equation in the semiclassical regime. © 2011 Cambridge University Press.

High Weak Order Methods for Stochastic Differential Equations Based on Modified Equations

KAUST Repository

Abdulle, Assyr

2012-01-01

© 2012 Society for Industrial and Applied Mathematics. Inspired by recent advances in the theory of modified differential equations, we propose a new methodology for constructing numerical integrators with high weak order for the time integration of stochastic differential equations. This approach is illustrated with the constructions of new methods of weak order two, in particular, semi-implicit integrators well suited for stiff (meansquare stable) stochastic problems, and implicit integrators that exactly conserve all quadratic first integrals of a stochastic dynamical system. Numerical examples confirm the theoretical results and show the versatility of our methodology.

Adjustment technique without explicit formation of normal equations /conjugate gradient method/

Science.gov (United States)

Saxena, N. K.

1974-01-01

For a simultaneous adjustment of a large geodetic triangulation system, a semiiterative technique is modified and used successfully. In this semiiterative technique, known as the conjugate gradient (CG) method, original observation equations are used, and thus the explicit formation of normal equations is avoided, 'huge' computer storage space being saved in the case of triangulation systems. This method is suitable even for very poorly conditioned systems where solution is obtained only after more iterations. A detailed study of the CG method for its application to large geodetic triangulation systems was done that also considered constraint equations with observation equations. It was programmed and tested on systems as small as two unknowns and three equations up to those as large as 804 unknowns and 1397 equations. When real data (573 unknowns, 965 equations) from a 1858-km-long triangulation system were used, a solution vector accurate to four decimal places was obtained in 2.96 min after 1171 iterations (i.e., 2.0 times the number of unknowns).

Novel Exponentially Fitted Two-Derivative Runge-Kutta Methods with Equation-Dependent Coefficients for First-Order Differential Equations

Directory of Open Access Journals (Sweden)

Yanping Yang

2016-01-01

Full Text Available The construction of exponentially fitted two-derivative Runge-Kutta (EFTDRK methods for the numerical solution of first-order differential equations is investigated. The revised EFTDRK methods proposed, with equation-dependent coefficients, take into consideration the errors produced in the internal stages to the update. The local truncation errors and stability of the new methods are analyzed. The numerical results are reported to show the accuracy of the new methods.

Multi-symplectic Preissmann methods for generalized Zakharov-Kuznetsov equation

International Nuclear Information System (INIS)

Wang Junjie; Yang Kuande; Wang Liantang

2012-01-01

Generalized Zakharov-Kuznetsov equation, a typical nonlinear wave equation, was studied based on the multi-symplectic theory in Hamilton space. The multi-symplectic formulations of generalized Zakharov-Kuznetsov equation with several conservation laws are presented. The multi-symplectic Preissmann method is used to discretize the formulations. The numerical experiment is given, and the results verify the efficiency of the multi-symplectic scheme. (authors)

Multiparameter extrapolation and deflation methods for solving equation systems

Directory of Open Access Journals (Sweden)

A. J. Hughes Hallett

1984-01-01

Full Text Available Most models in economics and the applied sciences are solved by first order iterative techniques, usually those based on the Gauss-Seidel algorithm. This paper examines the convergence of multiparameter extrapolations (accelerations of first order iterations, as an improved approximation to the Newton method for solving arbitrary nonlinear equation systems. It generalises my earlier results on single parameter extrapolations. Richardson's generalised method and the deflation method for detecting successive solutions in nonlinear equation systems are also presented as multiparameter extrapolations of first order iterations. New convergence results are obtained for those methods.

Systems of evolution equations and the singular perturbation method

International Nuclear Information System (INIS)

Mika, J.

Several fundamental theorems are presented important for the solution of linear evolution equations in the Banach space. The algorithm is deduced extending the solution of the system of singularly perturbed evolution equations into an asymptotic series with respect to a small positive parameter. The asymptotic convergence is shown of an approximate solution to the accurate solution. Singularly perturbed evolution equations of the resonance type were analysed. The special role is considered of the asymptotic equivalence of P1 equations obtained as the first order approximation if the spherical harmonics method is applied to the linear Boltzmann equation, and the diffusion equations of the linear transport theory where the small parameter approaches zero. (J.B.)

Analysis of spectral methods for the homogeneous Boltzmann equation

KAUST Repository

Filbet, Francis; Mouhot, Clé ment

2011-01-01

The development of accurate and fast algorithms for the Boltzmann collision integral and their analysis represent a challenging problem in scientific computing and numerical analysis. Recently, several works were devoted to the derivation of spectrally accurate schemes for the Boltzmann equation, but very few of them were concerned with the stability analysis of the method. In particular there was no result of stability except when the method was modified in order to enforce the positivity preservation, which destroys the spectral accuracy. In this paper we propose a new method to study the stability of homogeneous Boltzmann equations perturbed by smoothed balanced operators which do not preserve positivity of the distribution. This method takes advantage of the "spreading" property of the collision, together with estimates on regularity and entropy production. As an application we prove stability and convergence of spectral methods for the Boltzmann equation, when the discretization parameter is large enough (with explicit bound). © 2010 American Mathematical Society.

Analysis of spectral methods for the homogeneous Boltzmann equation

KAUST Repository

Filbet, Francis

2011-04-01

The development of accurate and fast algorithms for the Boltzmann collision integral and their analysis represent a challenging problem in scientific computing and numerical analysis. Recently, several works were devoted to the derivation of spectrally accurate schemes for the Boltzmann equation, but very few of them were concerned with the stability analysis of the method. In particular there was no result of stability except when the method was modified in order to enforce the positivity preservation, which destroys the spectral accuracy. In this paper we propose a new method to study the stability of homogeneous Boltzmann equations perturbed by smoothed balanced operators which do not preserve positivity of the distribution. This method takes advantage of the "spreading" property of the collision, together with estimates on regularity and entropy production. As an application we prove stability and convergence of spectral methods for the Boltzmann equation, when the discretization parameter is large enough (with explicit bound). © 2010 American Mathematical Society.

Waveform relaxation methods for implicit differential equations

NARCIS (Netherlands)

P.J. van der Houwen; W.A. van der Veen

1996-01-01

textabstractWe apply a Runge-Kutta-based waveform relaxation method to initial-value problems for implicit differential equations. In the implementation of such methods, a sequence of nonlinear systems has to be solved iteratively in each step of the integration process. The size of these systems

Fractional Complex Transform and exp-Function Methods for Fractional Differential Equations

Directory of Open Access Journals (Sweden)

Ahmet Bekir

2013-01-01

Full Text Available The exp-function method is presented for finding the exact solutions of nonlinear fractional equations. New solutions are constructed in fractional complex transform to convert fractional differential equations into ordinary differential equations. The fractional derivatives are described in Jumarie's modified Riemann-Liouville sense. We apply the exp-function method to both the nonlinear time and space fractional differential equations. As a result, some new exact solutions for them are successfully established.

Nonlinear Fredholm Integral Equation of the Second Kind with Quadrature Methods

Directory of Open Access Journals (Sweden)

M. Jafari Emamzadeh

2010-06-01

Full Text Available In this paper, a numerical method for solving the nonlinear Fredholm integral equation is presented. We intend to approximate the solution of this equation by quadrature methods and by doing so, we solve the nonlinear Fredholm integral equation more accurately. Several examples are given at the end of this paper

Runge-Kutta Methods for Linear Ordinary Differential Equations

Science.gov (United States)

Zingg, David W.; Chisholm, Todd T.

1997-01-01

Three new Runge-Kutta methods are presented for numerical integration of systems of linear inhomogeneous ordinary differential equations (ODES) with constant coefficients. Such ODEs arise in the numerical solution of the partial differential equations governing linear wave phenomena. The restriction to linear ODEs with constant coefficients reduces the number of conditions which the coefficients of the Runge-Kutta method must satisfy. This freedom is used to develop methods which are more efficient than conventional Runge-Kutta methods. A fourth-order method is presented which uses only two memory locations per dependent variable, while the classical fourth-order Runge-Kutta method uses three. This method is an excellent choice for simulations of linear wave phenomena if memory is a primary concern. In addition, fifth- and sixth-order methods are presented which require five and six stages, respectively, one fewer than their conventional counterparts, and are therefore more efficient. These methods are an excellent option for use with high-order spatial discretizations.

Method of mechanical quadratures for solving singular integral equations of various types

Science.gov (United States)

Sahakyan, A. V.; Amirjanyan, H. A.

2018-04-01

The method of mechanical quadratures is proposed as a common approach intended for solving the integral equations defined on finite intervals and containing Cauchy-type singular integrals. This method can be used to solve singular integral equations of the first and second kind, equations with generalized kernel, weakly singular equations, and integro-differential equations. The quadrature rules for several different integrals represented through the same coefficients are presented. This allows one to reduce the integral equations containing integrals of different types to a system of linear algebraic equations.

Error Analysis of Galerkin's Method for Semilinear Equations

Directory of Open Access Journals (Sweden)

Tadashi Kawanago

2012-01-01

Full Text Available We establish a general existence result for Galerkin's approximate solutions of abstract semilinear equations and conduct an error analysis. Our results may be regarded as some extension of a precedent work (Schultz 1969. The derivation of our results is, however, different from the discussion in his paper and is essentially based on the convergence theorem of Newton’s method and some techniques for deriving it. Some of our results may be applicable for investigating the quality of numerical verification methods for solutions of ordinary and partial differential equations.

Approximate Method for Solving the Linear Fuzzy Delay Differential Equations

Directory of Open Access Journals (Sweden)

S. Narayanamoorthy

2015-01-01

Full Text Available We propose an algorithm of the approximate method to solve linear fuzzy delay differential equations using Adomian decomposition method. The detailed algorithm of the approach is provided. The approximate solution is compared with the exact solution to confirm the validity and efficiency of the method to handle linear fuzzy delay differential equation. To show this proper features of this proposed method, numerical example is illustrated.

Numerical method for the nonlinear Fokker-Planck equation

International Nuclear Information System (INIS)

Zhang, D.S.; Wei, G.W.; Kouri, D.J.; Hoffman, D.K.

1997-01-01

A practical method based on distributed approximating functionals (DAFs) is proposed for numerically solving a general class of nonlinear time-dependent Fokker-Planck equations. The method relies on a numerical scheme that couples the usual path-integral concept to the DAF idea. The high accuracy and reliability of the method are illustrated by applying it to an exactly solvable nonlinear Fokker-Planck equation, and the method is compared with the accurate K-point Stirling interpolation formula finite-difference method. The approach is also used successfully to solve a nonlinear self-consistent dynamic mean-field problem for which both the cumulant expansion and scaling theory have been found by Drozdov and Morillo [Phys. Rev. E 54, 931 (1996)] to be inadequate to describe the occurrence of a long-lived transient bimodality. The standard interpretation of the transient bimodality in terms of the flat region in the kinetic potential fails for the present case. An alternative analysis based on the effective potential of the Schroedinger-like Fokker-Planck equation is suggested. Our analysis of the transient bimodality is strongly supported by two examples that are numerically much more challenging than other examples that have been previously reported for this problem. copyright 1997 The American Physical Society

A functional-analytic method for the study of difference equations

Directory of Open Access Journals (Sweden)

Panayiotis D. Siafarikas

2004-07-01

Full Text Available We will give the generalization of a recently developed functional-analytic method for studying linear and nonlinear, ordinary and partial, difference equations in the ℓp1 and ℓp2 spaces, p∈ℕ, p≥1. The method will be illustrated by use of two examples concerning a nonlinear ordinary difference equation known as the Putnam equation, and a linear partial difference equation of three variables describing the discrete Newton law of cooling in three dimensions.

New finite volume methods for approximating partial differential equations on arbitrary meshes

International Nuclear Information System (INIS)

Hermeline, F.

2008-12-01

This dissertation presents some new methods of finite volume type for approximating partial differential equations on arbitrary meshes. The main idea lies in solving twice the problem to be dealt with. One addresses the elliptic equations with variable (anisotropic, antisymmetric, discontinuous) coefficients, the parabolic linear or non linear equations (heat equation, radiative diffusion, magnetic diffusion with Hall effect), the wave type equations (Maxwell, acoustics), the elasticity and Stokes'equations. Numerous numerical experiments show the good behaviour of this type of method. (author)

Tensor-GMRES method for large sparse systems of nonlinear equations

Science.gov (United States)

Feng, Dan; Pulliam, Thomas H.

1994-01-01

This paper introduces a tensor-Krylov method, the tensor-GMRES method, for large sparse systems of nonlinear equations. This method is a coupling of tensor model formation and solution techniques for nonlinear equations with Krylov subspace projection techniques for unsymmetric systems of linear equations. Traditional tensor methods for nonlinear equations are based on a quadratic model of the nonlinear function, a standard linear model augmented by a simple second order term. These methods are shown to be significantly more efficient than standard methods both on nonsingular problems and on problems where the Jacobian matrix at the solution is singular. A major disadvantage of the traditional tensor methods is that the solution of the tensor model requires the factorization of the Jacobian matrix, which may not be suitable for problems where the Jacobian matrix is large and has a 'bad' sparsity structure for an efficient factorization. We overcome this difficulty by forming and solving the tensor model using an extension of a Newton-GMRES scheme. Like traditional tensor methods, we show that the new tensor method has significant computational advantages over the analogous Newton counterpart. Consistent with Krylov subspace based methods, the new tensor method does not depend on the factorization of the Jacobian matrix. As a matter of fact, the Jacobian matrix is never needed explicitly.

Extended sine-Gordon Equation Method and Its Application to Maccari's System

International Nuclear Information System (INIS)

Song Lina; Zhang Hongqing

2005-01-01

An extended sine-Gordon equation method is proposed to construct exact travelling wave solutions to Maccari's equation based upon a generalized sine-Gordon equation. It is shown that more new travelling wave solutions can be found by this new method, which include bell-shaped soliton solutions, kink-shaped soliton solutions, periodic wave solution, and new travelling waves.

GHM method for obtaining rationalsolutions of nonlinear differential equations.

Science.gov (United States)

Vazquez-Leal, Hector; Sarmiento-Reyes, Arturo

2015-01-01

In this paper, we propose the application of the general homotopy method (GHM) to obtain rational solutions of nonlinear differential equations. It delivers a high precision representation of the nonlinear differential equation using a few linear algebraic terms. In order to assess the benefits of this proposal, three nonlinear problems are solved and compared against other semi-analytic methods or numerical methods. The obtained results show that GHM is a powerful tool, capable to generate highly accurate rational solutions. AMS subject classification 34L30.

Soliton solution for nonlinear partial differential equations by cosine-function method

International Nuclear Information System (INIS)

Ali, A.H.A.; Soliman, A.A.; Raslan, K.R.

2007-01-01

In this Letter, we established a traveling wave solution by using Cosine-function algorithm for nonlinear partial differential equations. The method is used to obtain the exact solutions for five different types of nonlinear partial differential equations such as, general equal width wave equation (GEWE), general regularized long wave equation (GRLW), general Korteweg-de Vries equation (GKdV), general improved Korteweg-de Vries equation (GIKdV), and Coupled equal width wave equations (CEWE), which are the important soliton equations

New extended (G'/G)-expansion method to solve nonlinear evolution equation: the (3 + 1)-dimensional potential-YTSF equation.

Science.gov (United States)

Roshid, Harun-Or-; Akbar, M Ali; Alam, Md Nur; Hoque, Md Fazlul; Rahman, Nizhum

2014-01-01

In this article, a new extended (G'/G) -expansion method has been proposed for constructing more general exact traveling wave solutions of nonlinear evolution equations with the aid of symbolic computation. In order to illustrate the validity and effectiveness of the method, we pick the (3 + 1)-dimensional potential-YTSF equation. As a result, abundant new and more general exact solutions have been achieved of this equation. It has been shown that the proposed method provides a powerful mathematical tool for solving nonlinear wave equations in applied mathematics, engineering and mathematical physics.

On Solution of a Fractional Diffusion Equation by Homotopy Transform Method

International Nuclear Information System (INIS)

Salah, A.; Hassan, S.S.A.

2012-01-01

The homotopy analysis transform method (HATM) is applied in this work in order to find the analytical solution of fractional diffusion equations (FDE). These equations are obtained from standard diffusion equations by replacing a second-order space derivative by a fractional derivative of order α and a first order time derivative by a fractional derivative. Furthermore, some examples are given. Numerical results show that the homotopy analysis transform method is easy to implement and accurate when applied to a fractional diffusion equations.

A Marker Method for the Solution of the Damped Burgers' Equation

International Nuclear Information System (INIS)

Jerome L.V. Lewandowski

2005-01-01

A new method for the solution of the damped Burgers equation is described. The marker method relies on the definition of a convective field associated with the underlying partial differential equation; the information about the approximate solution is associated with the response of an ensemble of markers to this convective field. Some key aspects of the method, such as the selection of the shape function and the initial loading, are discussed in some details. The marker method is applicable to a general class of nonlinear dispersive partial differential equations

students' preference of method of solving simultaneous equations

African Journals Online (AJOL)

Ugboduma,Samuel.O.

substitution method irrespective of their gender for solving simultaneous equations. A recommendation ... advantage given to one method over others. Students' interest .... from two (2) single girls' schools, two (2) single boys schools and ten.

Application of He’s Variational Iteration Method to Nonlinear Helmholtz Equation and Fifth-Order KDV Equation

DEFF Research Database (Denmark)

Miansari, Mo; Miansari, Me; Barari, Amin

2009-01-01

In this article, He’s variational iteration method (VIM), is implemented to solve the linear Helmholtz partial differential equation and some nonlinear fifth-order Korteweg-de Vries (FKdV) partial differential equations with specified initial conditions. The initial approximations can be freely c...

Iterative method of the parameter variation for solution of nonlinear functional equations

International Nuclear Information System (INIS)

Davidenko, D.F.

1975-01-01

The iteration method of parameter variation is used for solving nonlinear functional equations in Banach spaces. The authors consider some methods for numerical integration of ordinary first-order differential equations and construct the relevant iteration methods of parameter variation, both one- and multifactor. They also discuss problems of mathematical substantiation of the method, study the conditions and rate of convergence, estimate the error. The paper considers the application of the method to specific functional equations

New Quasi-Newton Method for Solving Systems of Nonlinear Equations

Czech Academy of Sciences Publication Activity Database

Lukšan, Ladislav; Vlček, Jan

2017-01-01

Roč. 62, č. 2 (2017), s. 121-134 ISSN 0862-7940 R&D Projects: GA ČR GA13-06684S Institutional support: RVO:67985807 Keywords : nonlinear equations * systems of equations * trust-region methods * quasi-Newton methods * adjoint Broyden methods * numerical algorithms * numerical experiments Subject RIV: BA - General Mathematics OBOR OECD: Applied mathematics Impact factor: 0.618, year: 2016 http://hdl.handle.net/10338.dmlcz/146699

Receptor binding kinetics equations: Derivation using the Laplace transform method.

Science.gov (United States)

Hoare, Sam R J

Measuring unlabeled ligand receptor binding kinetics is valuable in optimizing and understanding drug action. Unfortunately, deriving equations for estimating kinetic parameters is challenging because it involves calculus; integration can be a frustrating barrier to the pharmacologist seeking to measure simple rate parameters. Here, a well-known tool for simplifying the derivation, the Laplace transform, is applied to models of receptor-ligand interaction. The method transforms differential equations to a form in which simple algebra can be applied to solve for the variable of interest, for example the concentration of ligand-bound receptor. The goal is to provide instruction using familiar examples, to enable investigators familiar with handling equilibrium binding equations to derive kinetic equations for receptor-ligand interaction. First, the Laplace transform is used to derive the equations for association and dissociation of labeled ligand binding. Next, its use for unlabeled ligand kinetic equations is exemplified by a full derivation of the kinetics of competitive binding equation. Finally, new unlabeled ligand equations are derived using the Laplace transform. These equations incorporate a pre-incubation step with unlabeled or labeled ligand. Four equations for measuring unlabeled ligand kinetics were compared and the two new equations verified by comparison with numerical solution. Importantly, the equations have not been verified with experimental data because no such experiments are evident in the literature. Equations were formatted for use in the curve-fitting program GraphPad Prism 6.0 and fitted to simulated data. This description of the Laplace transform method will enable pharmacologists to derive kinetic equations for their model or experimental paradigm under study. Application of the transform will expand the set of equations available for the pharmacologist to measure unlabeled ligand binding kinetics, and for other time

An effective method for finding special solutions of nonlinear differential equations with variable coefficients

International Nuclear Information System (INIS)

Qin Maochang; Fan Guihong

2008-01-01

There are many interesting methods can be utilized to construct special solutions of nonlinear differential equations with constant coefficients. However, most of these methods are not applicable to nonlinear differential equations with variable coefficients. A new method is presented in this Letter, which can be used to find special solutions of nonlinear differential equations with variable coefficients. This method is based on seeking appropriate Bernoulli equation corresponding to the equation studied. Many well-known equations are chosen to illustrate the application of this method

Exact solitary wave solutions for some nonlinear evolution equations via Exp-function method

International Nuclear Information System (INIS)

Ebaid, A.

2007-01-01

Based on the Exp-function method, exact solutions for some nonlinear evolution equations are obtained. The KdV equation, Burgers' equation and the combined KdV-mKdV equation are chosen to illustrate the effectiveness of the method

[Series: Utilization of Differential Equations and Methods for Solving Them in Medical Physics (2)].

Science.gov (United States)

Murase, Kenya

2015-01-01

In this issue, symbolic methods for solving differential equations were firstly introduced. Of the symbolic methods, Laplace transform method was also introduced together with some examples, in which this method was applied to solving the differential equations derived from a two-compartment kinetic model and an equivalent circuit model for membrane potential. Second, series expansion methods for solving differential equations were introduced together with some examples, in which these methods were used to solve Bessel's and Legendre's differential equations. In the next issue, simultaneous differential equations and various methods for solving these differential equations will be introduced together with some examples in medical physics.

Classical and modular methods applied to Diophantine equations

NARCIS (Netherlands)

Dahmen, S.R.

2008-01-01

Deep methods from the theory of elliptic curves and modular forms have been used to prove Fermat's last theorem and solve other Diophantine equations. These so-called modular methods can often benefit from information obtained by other, classical, methods from number theory; and vice versa. In our

Generalized multiscale finite element methods. nonlinear elliptic equations

KAUST Repository

Efendiev, Yalchin R.; Galvis, Juan; Li, Guanglian; Presho, Michael

2013-01-01

In this paper we use the Generalized Multiscale Finite Element Method (GMsFEM) framework, introduced in [26], in order to solve nonlinear elliptic equations with high-contrast coefficients. The proposed solution method involves linearizing the equation so that coarse-grid quantities of previous solution iterates can be regarded as auxiliary parameters within the problem formulation. With this convention, we systematically construct respective coarse solution spaces that lend themselves to either continuous Galerkin (CG) or discontinuous Galerkin (DG) global formulations. Here, we use Symmetric Interior Penalty Discontinuous Galerkin approach. Both methods yield a predictable error decline that depends on the respective coarse space dimension, and we illustrate the effectiveness of the CG and DG formulations by offering a variety of numerical examples. © 2014 Global-Science Press.

Improved stochastic approximation methods for discretized parabolic partial differential equations

Science.gov (United States)

Guiaş, Flavius

2016-12-01

We present improvements of the stochastic direct simulation method, a known numerical scheme based on Markov jump processes which is used for approximating solutions of ordinary differential equations. This scheme is suited especially for spatial discretizations of evolution partial differential equations (PDEs). By exploiting the full path simulation of the stochastic method, we use this first approximation as a predictor and construct improved approximations by Picard iterations, Runge-Kutta steps, or a combination. This has as consequence an increased order of convergence. We illustrate the features of the improved method at a standard benchmark problem, a reaction-diffusion equation modeling a combustion process in one space dimension (1D) and two space dimensions (2D).

An Evaluation of Kernel Equating: Parallel Equating with Classical Methods in the SAT Subject Tests[TM] Program. Research Report. ETS RR-09-06

Science.gov (United States)

Grant, Mary C.; Zhang, Lilly; Damiano, Michele

2009-01-01

This study investigated kernel equating methods by comparing these methods to operational equatings for two tests in the SAT Subject Tests[TM] program. GENASYS (ETS, 2007) was used for all equating methods and scaled score kernel equating results were compared to Tucker, Levine observed score, chained linear, and chained equipercentile equating…

An introduction to neural network methods for differential equations

CERN Document Server

Yadav, Neha; Kumar, Manoj

2015-01-01

This book introduces a variety of neural network methods for solving differential equations arising in science and engineering. The emphasis is placed on a deep understanding of the neural network techniques, which has been presented in a mostly heuristic and intuitive manner. This approach will enable the reader to understand the working, efficiency and shortcomings of each neural network technique for solving differential equations. The objective of this book is to provide the reader with a sound understanding of the foundations of neural networks, and a comprehensive introduction to neural network methods for solving differential equations together with recent developments in the techniques and their applications. The book comprises four major sections. Section I consists of a brief overview of differential equations and the relevant physical problems arising in science and engineering. Section II illustrates the history of neural networks starting from their beginnings in the 1940s through to the renewed...

Evaluation of Maryland abutment scour equation through selected threshold velocity methods

Science.gov (United States)

Benedict, S.T.

2010-01-01

The U.S. Geological Survey, in cooperation with the Maryland State Highway Administration, used field measurements of scour to evaluate the sensitivity of the Maryland abutment scour equation to the critical (or threshold) velocity variable. Four selected methods for estimating threshold velocity were applied to the Maryland abutment scour equation, and the predicted scour to the field measurements were compared. Results indicated that performance of the Maryland abutment scour equation was sensitive to the threshold velocity with some threshold velocity methods producing better estimates of predicted scour than did others. In addition, results indicated that regional stream characteristics can affect the performance of the Maryland abutment scour equation with moderate-gradient streams performing differently from low-gradient streams. On the basis of the findings of the investigation, guidance for selecting threshold velocity methods for application to the Maryland abutment scour equation are provided, and limitations are noted.

Application of Monte Carlo method to solving boundary value problem of differential equations

International Nuclear Information System (INIS)

Zuo Yinghong; Wang Jianguo

2012-01-01

This paper introduces the foundation of the Monte Carlo method and the way how to generate the random numbers. Based on the basic thought of the Monte Carlo method and finite differential method, the stochastic model for solving the boundary value problem of differential equations is built. To investigate the application of the Monte Carlo method to solving the boundary value problem of differential equations, the model is used to solve Laplace's equations with the first boundary condition and the unsteady heat transfer equation with initial values and boundary conditions. The results show that the boundary value problem of differential equations can be effectively solved with the Monte Carlo method, and the differential equations with initial condition can also be calculated by using a stochastic probability model which is based on the time-domain finite differential equations. Both the simulation results and theoretical analyses show that the errors of numerical results are lowered as the number of simulation particles is increased. (authors)

HAM-Based Adaptive Multiscale Meshless Method for Burgers Equation

Directory of Open Access Journals (Sweden)

Shu-Li Mei

2013-01-01

Full Text Available Based on the multilevel interpolation theory, we constructed a meshless adaptive multiscale interpolation operator (MAMIO with the radial basis function. Using this operator, any nonlinear partial differential equations such as Burgers equation can be discretized adaptively in physical spaces as a nonlinear matrix ordinary differential equation. In order to obtain the analytical solution of the system of ODEs, the homotopy analysis method (HAM proposed by Shijun Liao was developed to solve the system of ODEs by combining the precise integration method (PIM which can be employed to get the analytical solution of linear system of ODEs. The numerical experiences show that HAM is not sensitive to the time step, and so the arithmetic error is mainly derived from the discrete in physical space.

A nondissipative simulation method for the drift kinetic equation

International Nuclear Information System (INIS)

Watanabe, Tomo-Hiko; Sugama, Hideo; Sato, Tetsuya

2001-07-01

With the aim to study the ion temperature gradient (ITG) driven turbulence, a nondissipative kinetic simulation scheme is developed and comprehensively benchmarked. The new simulation method preserving the time-reversibility of basic kinetic equations can successfully reproduce the analytical solutions of asymmetric three-mode ITG equations which are extended to provide a more general reference for benchmarking than the previous work [T.-H. Watanabe, H. Sugama, and T. Sato: Phys. Plasmas 7 (2000) 984]. It is also applied to a dissipative three-mode system, and shows a good agreement with the analytical solution. The nondissipative simulation result of the ITG turbulence accurately satisfies the entropy balance equation. Usefulness of the nondissipative method for the drift kinetic simulations is confirmed in comparisons with other dissipative schemes. (author)

A moving mesh finite difference method for equilibrium radiation diffusion equations

Energy Technology Data Exchange (ETDEWEB)

Yang, Xiaobo, E-mail: xwindyb@126.com [Department of Mathematics, College of Science, China University of Mining and Technology, Xuzhou, Jiangsu 221116 (China); Huang, Weizhang, E-mail: whuang@ku.edu [Department of Mathematics, University of Kansas, Lawrence, KS 66045 (United States); Qiu, Jianxian, E-mail: jxqiu@xmu.edu.cn [School of Mathematical Sciences and Fujian Provincial Key Laboratory of Mathematical Modeling and High-Performance Scientific Computing, Xiamen University, Xiamen, Fujian 361005 (China)

2015-10-01

An efficient moving mesh finite difference method is developed for the numerical solution of equilibrium radiation diffusion equations in two dimensions. The method is based on the moving mesh partial differential equation approach and moves the mesh continuously in time using a system of meshing partial differential equations. The mesh adaptation is controlled through a Hessian-based monitor function and the so-called equidistribution and alignment principles. Several challenging issues in the numerical solution are addressed. Particularly, the radiation diffusion coefficient depends on the energy density highly nonlinearly. This nonlinearity is treated using a predictor–corrector and lagged diffusion strategy. Moreover, the nonnegativity of the energy density is maintained using a cutoff method which has been known in literature to retain the accuracy and convergence order of finite difference approximation for parabolic equations. Numerical examples with multi-material, multiple spot concentration situations are presented. Numerical results show that the method works well for radiation diffusion equations and can produce numerical solutions of good accuracy. It is also shown that a two-level mesh movement strategy can significantly improve the efficiency of the computation.

A moving mesh finite difference method for equilibrium radiation diffusion equations

International Nuclear Information System (INIS)

Yang, Xiaobo; Huang, Weizhang; Qiu, Jianxian

2015-01-01

An efficient moving mesh finite difference method is developed for the numerical solution of equilibrium radiation diffusion equations in two dimensions. The method is based on the moving mesh partial differential equation approach and moves the mesh continuously in time using a system of meshing partial differential equations. The mesh adaptation is controlled through a Hessian-based monitor function and the so-called equidistribution and alignment principles. Several challenging issues in the numerical solution are addressed. Particularly, the radiation diffusion coefficient depends on the energy density highly nonlinearly. This nonlinearity is treated using a predictor–corrector and lagged diffusion strategy. Moreover, the nonnegativity of the energy density is maintained using a cutoff method which has been known in literature to retain the accuracy and convergence order of finite difference approximation for parabolic equations. Numerical examples with multi-material, multiple spot concentration situations are presented. Numerical results show that the method works well for radiation diffusion equations and can produce numerical solutions of good accuracy. It is also shown that a two-level mesh movement strategy can significantly improve the efficiency of the computation

Numerical Solution of Fuzzy Differential Equations by Runge-Kutta Verner Method

Directory of Open Access Journals (Sweden)

T. Jayakumar

2015-01-01

Full Text Available In this paper we study the numerical methods for Fuzzy Differential equations by an application of the Runge-Kutta Verner method for fuzzy differential equations. We prove a convergence result and give numerical examples to illustrate the theory.

A multilevel correction adaptive finite element method for Kohn-Sham equation

Science.gov (United States)

Hu, Guanghui; Xie, Hehu; Xu, Fei

2018-02-01

In this paper, an adaptive finite element method is proposed for solving Kohn-Sham equation with the multilevel correction technique. In the method, the Kohn-Sham equation is solved on a fixed and appropriately coarse mesh with the finite element method in which the finite element space is kept improving by solving the derived boundary value problems on a series of adaptively and successively refined meshes. A main feature of the method is that solving large scale Kohn-Sham system is avoided effectively, and solving the derived boundary value problems can be handled efficiently by classical methods such as the multigrid method. Hence, the significant acceleration can be obtained on solving Kohn-Sham equation with the proposed multilevel correction technique. The performance of the method is examined by a variety of numerical experiments.

Solving the Helmholtz equation in conformal mapped ARROWstructures using homotopy perturbation method

DEFF Research Database (Denmark)

Reck, Kasper; Thomsen, Erik Vilain; Hansen, Ole

2011-01-01

. The solution of the mapped Helmholtz equation is found by solving an infinite series of Poisson equations using two dimensional Fourier series. The solution is entirely based on analytical expressions and is not mesh dependent. The analytical results are compared to a numerical (finite element method) solution......The scalar wave equation, or Helmholtz equation, describes within a certain approximation the electromagnetic field distribution in a given system. In this paper we show how to solve the Helmholtz equation in complex geometries using conformal mapping and the homotopy perturbation method...

Application of the N/D-method to the Roy equations

International Nuclear Information System (INIS)

Heemskerk, A.C.

1978-01-01

The subject of this thesis is the study and numerical application of the Roy equations for elastic pion-pion scattering. These equations were introduced by S.M. Roy in 1971 and express the minimal requirements any ππ partial wave amplitude should satisfy, such as analyticity and crossing symmetry. When combined with the unitarity relations for the partial waves a complete system of equations (the Roy system) results valid in the low-energy region. The author has studied the existance and uniqueness question for the Roy system and developed successful strategies to solve the system by means of iteration. Several solutions are presented and discussed. The main ingredients in the present approach are a) the use of the N/D-method to regularize the singular Roy equations and implement unitarity and b) the introduction of various iteration methods which do not need any derivatives to solve the resulting nonlinear system of equations numerically. (Auth.)

Stability of generalized Runge-Kutta methods for stiff kinetics coupled differential equations

International Nuclear Information System (INIS)

Aboanber, A E

2006-01-01

A stability and efficiency improved class of generalized Runge-Kutta methods of order 4 are developed for the numerical solution of stiff system kinetics equations for linear and/or nonlinear coupled differential equations. The determination of the coefficients required by the method is precisely obtained from the so-called equations of condition which in turn are derived by an approach based on Butcher series. Since the equations of condition are fewer in number, free parameters can be chosen for optimizing any desired feature of the process. A further related coefficient set with different values of these parameters and the region of absolute stability of the method have been introduced. In addition, the A(α) stability properties of the method are investigated. Implementing the method in a personal computer estimated the accuracy and speed of calculations and verified the good performances of the proposed new schemes for several sample problems of the stiff system point kinetics equations with reactivity feedback

An auxiliary differential equation FDTD method for anisotropic magnetized plasmas

International Nuclear Information System (INIS)

Liu Shaobin; Mo Jinjun; Yuan Naichang

2004-01-01

An auxiliary differential equation finite-difference time-domain (ADE-FDTD) methodology for anisotropic magnetized plasmas is derived. The method is based on a difference approximation of the auxiliary differential equation. A comparison with the JEC method is included. The CPU time saving by several times and accuracy of the method are confirmed by computing the reflection and transmission through a magnetized plasma layer with the direction of propagation parallel to the direction of the biasing field

GoFFish: Graph-Oriented Framework for Foresight and Insight Using Scalable Heuristics

Science.gov (United States)

2015-09-01

A. Biem, E. Bouillet, H. Feng, A. Ranganathan , A. Riabov, O. Verscheure, H. Koutsopoulos, and C. Moran, “Ibm infos- phere streams for scalable, real...Systems and Software. Elsevier, 2013, vol. 86, no. 1, pp. 2–11. [5] A. Biem, E. Bouillet, H. Feng, A. Ranganathan , A. Riabov, O. Verscheure, H...Feng, A. Ranganathan , A. Riabov, O. Verscheure, H. Koutsopoulos, and C. Moran. Ibm infosphere streams for scalable, real-time, intelligent

Solving nonlinear evolution equation system using two different methods

Science.gov (United States)

Kaplan, Melike; Bekir, Ahmet; Ozer, Mehmet N.

2015-12-01

This paper deals with constructing more general exact solutions of the coupled Higgs equation by using the (G0/G, 1/G)-expansion and (1/G0)-expansion methods. The obtained solutions are expressed by three types of functions: hyperbolic, trigonometric and rational functions with free parameters. It has been shown that the suggested methods are productive and will be used to solve nonlinear partial differential equations in applied mathematics and engineering. Throughout the paper, all the calculations are made with the aid of the Maple software.

A highly accurate method to solve Fisher's equation

Indian Academy of Sciences (India)

The solution of the Helmholtz equation was approximated by a sixth-order compact finite difference. (CFD6) method in [29]. In [30], a CFD6 scheme has been presented to ... efficiency of the proposed method are reported in §3. Finally .... our discussion, one can apply the proposed method to solve the more general problem.

Integral equation methods for vesicle electrohydrodynamics in three dimensions

Science.gov (United States)

Veerapaneni, Shravan

2016-12-01

In this paper, we develop a new boundary integral equation formulation that describes the coupled electro- and hydro-dynamics of a vesicle suspended in a viscous fluid and subjected to external flow and electric fields. The dynamics of the vesicle are characterized by a competition between the elastic, electric and viscous forces on its membrane. The classical Taylor-Melcher leaky-dielectric model is employed for the electric response of the vesicle and the Helfrich energy model combined with local inextensibility is employed for its elastic response. The coupled governing equations for the vesicle position and its transmembrane electric potential are solved using a numerical method that is spectrally accurate in space and first-order in time. The method uses a semi-implicit time-stepping scheme to overcome the numerical stiffness associated with the governing equations.

Higher-Order Integral Equation Methods in Computational Electromagnetics

DEFF Research Database (Denmark)

Jørgensen, Erik; Meincke, Peter

Higher-order integral equation methods have been investigated. The study has focused on improving the accuracy and efficiency of the Method of Moments (MoM) applied to electromagnetic problems. A new set of hierarchical Legendre basis functions of arbitrary order is developed. The new basis...

Exact solutions of nonlinear fractional differential equations by (G′/G)-expansion method

International Nuclear Information System (INIS)

Bekir Ahmet; Güner Özkan

2013-01-01

In this paper, we use the fractional complex transform and the (G′/G)-expansion method to study the nonlinear fractional differential equations and find the exact solutions. The fractional complex transform is proposed to convert a partial fractional differential equation with Jumarie's modified Riemann—Liouville derivative into its ordinary differential equation. It is shown that the considered transform and method are very efficient and powerful in solving wide classes of nonlinear fractional order equations

Comparative analysis of solution methods of the punctual kinetic equations

International Nuclear Information System (INIS)

Hernandez S, A.

2003-01-01

The following one written it presents a comparative analysis among different analytical solutions for the punctual kinetics equation, which present two variables of interest: a) the temporary behavior of the neutronic population, and b) The temporary behavior of the different groups of precursors of delayed neutrons. The first solution is based on a method that solves the transfer function of the differential equation for the neutronic population, in which intends to obtain the different poles that give the stability of this transfer function. In this section it is demonstrated that the temporary variation of the reactivity of the system can be managed as it is required, since the integration time for this method doesn't affect the result. However, the second solution is based on an iterative method like that of Runge-Kutta or the Euler method where the algorithm was only used to solve first order differential equations giving this way solution to each differential equation that conforms the equations of punctual kinetics. In this section it is demonstrated that only it can obtain a correct temporary behavior of the neutronic population when it is integrated on an interval of very short time, forcing to the temporary variation of the reactivity to change very quick way without one has some control about the time. In both methods the same change is used so much in the reactivity of the system like in the integration times, giving validity to the results graph the one the temporary behavior of the neutronic population vs. time. (Author)

Inverse Free Iterative Methods for Nonlinear Ill-Posed Operator Equations

Directory of Open Access Journals (Sweden)

Ioannis K. Argyros

2014-01-01

ill-posed operator equation F(x=y. The proposed method is a modified form of Tikhonov gradient (TIGRA method considered by Ramlau (2003. The regularization parameter is chosen according to the balancing principle considered by Pereverzev and Schock (2005. The error estimate is derived under a general source condition and is of optimal order. Some numerical examples involving integral equations are also given in this paper.

Operator Splitting Methods for Degenerate Convection-Diffusion Equations I: Convergence and Entropy Estimates

Energy Technology Data Exchange (ETDEWEB)

Holden, Helge; Karlsen, Kenneth H.; Lie, Knut-Andreas

1999-10-01

We present and analyze a numerical method for the solution of a class of scalar, multi-dimensional, nonlinear degenerate convection-diffusion equations. The method is based on operator splitting to separate the convective and the diffusive terms in the governing equation. The nonlinear, convective part is solved using front tracking and dimensional splitting, while the nonlinear diffusion equation is solved by a suitable difference scheme. We verify L{sup 1} compactness of the corresponding set of approximate solutions and derive precise entropy estimates. In particular, these results allow us to pass to the limit in our approximations and recover an entropy solution of the problem in question. The theory presented covers a large class of equations. Important subclasses are hyperbolic conservation laws, porous medium type equations, two-phase reservoir flow equations, and strongly degenerate equations coming from the recent theory of sedimentation-consolidation processes. A thorough numerical investigation of the method analyzed in this paper (and similar methods) is presented in a companion paper. (author)

Application of the Generalized Differential Quadrature Method in Solving Burgers' Equations

International Nuclear Information System (INIS)

Mokhtari, R.; Toodar, A. Samadi; Chegini, N.G.

2011-01-01

The aim of this paper is to obtain numerical solutions of the one-dimensional, two-dimensional and coupled Burgers' equations through the generalized differential quadrature method (GDQM). The polynomial-based differential quadrature (PDQ) method is employed and the obtained system of ordinary differential equations is solved via the total variation diminishing Runge-Kutta (TVD-RK) method. The numerical solutions are satisfactorily coincident with the exact solutions. The method can compete against the methods applied in the literature. (general)

Integral transform method for solving time fractional systems and fractional heat equation

Directory of Open Access Journals (Sweden)

Arman Aghili

2014-01-01

Full Text Available In the present paper, time fractional partial differential equation is considered, where the fractional derivative is defined in the Caputo sense. Laplace transform method has been applied to obtain an exact solution. The authors solved certain homogeneous and nonhomogeneous time fractional heat equations using integral transform. Transform method is a powerful tool for solving fractional singular Integro - differential equations and PDEs. The result reveals that the transform method is very convenient and effective.

Novel method for solution of coupled radial Schrödinger equations

International Nuclear Information System (INIS)

Ershov, S. N.; Vaagen, J. S.; Zhukov, M. V.

2011-01-01

One of the major problems in numerical solution of coupled differential equations is the maintenance of linear independence for different sets of solution vectors. A novel method for solution of radial Schrödinger equations is suggested. It consists of rearrangement of coupled equations in a way that is appropriate to avoid usual numerical instabilities associated with components of the wave function in their classically forbidden regions. Applications of the new method for nuclear structure calculations within the hyperspherical harmonics approach are given.

Predictive equation of state method for heavy materials based on the Dirac equation and density functional theory

Science.gov (United States)

Wills, John M.; Mattsson, Ann E.

2012-02-01

Density functional theory (DFT) provides a formally predictive base for equation of state properties. Available approximations to the exchange/correlation functional provide accurate predictions for many materials in the periodic table. For heavy materials however, DFT calculations, using available functionals, fail to provide quantitative predictions, and often fail to be even qualitative. This deficiency is due both to the lack of the appropriate confinement physics in the exchange/correlation functional and to approximations used to evaluate the underlying equations. In order to assess and develop accurate functionals, it is essential to eliminate all other sources of error. In this talk we describe an efficient first-principles electronic structure method based on the Dirac equation and compare the results obtained with this method with other methods generally used. Implications for high-pressure equation of state of relativistic materials are demonstrated in application to Ce and the light actinides. Sandia National Laboratories is a multi-program laboratory managed andoperated by Sandia Corporation, a wholly owned subsidiary of Lockheed Martin Corporation, for the U.S. Department of Energy's National Nuclear Security Administration under contract DE-AC04-94AL85000.

[Series: Utilization of Differential Equations and Methods for Solving Them in Medical Physics (1)].

Science.gov (United States)

Murase, Kenya

2014-01-01

Utilization of differential equations and methods for solving them in medical physics are presented. First, the basic concept and the kinds of differential equations were overviewed. Second, separable differential equations and well-known first-order and second-order differential equations were introduced, and the methods for solving them were described together with several examples. In the next issue, the symbolic and series expansion methods for solving differential equations will be mainly introduced.

Implicit Boundary Integral Methods for the Helmholtz Equation in Exterior Domains

Science.gov (United States)

2016-06-01

solve the Helmholtz equation as ∂Ω goes through significant change in its shape and topology — applications for which implicit representation of the...boundary-value problems for the wave equation and maxwell’s equations. Russian Math . Surv., 1965. [16] S. Reutskiy. The method of fundamental

The C{sub n} method for approximation of the Boltzmann equation; La methode C{sub n} d'approximation de l'equation de Boltzmann

Energy Technology Data Exchange (ETDEWEB)

Benoist, P; Kavenoky, A [Commissariat a l' Energie Atomique, Saclay (France). Centre d' Etudes Nucleaires

1968-01-15

In a new method of approximation of the Boltzmann equation, one starts from a particular form of the equation which involves only the angular flux at the boundary of the considered medium and where the space variable does not appear explicitly. Expanding in orthogonal polynomials the angular flux of neutrons leaking from the medium and making no assumption about the angular flux within the medium, very good approximations to several classical plane geometry problems, i.e. the albedo of slabs and the transmission by slabs, the extrapolation length of the Milne problem, the spectrum of neutrons reflected by a semi-infinite slowing down medium. The method can be extended to other geometries. (authors) [French] On etablit une nouvelle methode d'approximation pour l'equation de Boltzmann en partant d'une forme particuliere de cette equation qui n'implique que le flux angulaire a la frontiere du milieu et ou les variables d'espace n'apparaissent pas explicitement. Par un developpement en polynomes orthogonaux du flux angulaire sortant du milieu et sans faire d'hypothese sur le flux angulaire a l'interieur du milieu, on obtient de tres bonnes approximations pour plusieurs problemes classiques en geometrie plane: l'albedo et le facteur de transmission des plaques, la longueur d'extrapolation du probleme de Milne, le spectre des neutrons reflechis par un milieu semi-infini ralentisseur. La methode se generalise a d'autres geometries. (auteurs)

On method of solving third-order ordinary differential equations directly using Bernstein polynomials

Science.gov (United States)

Khataybeh, S. N.; Hashim, I.

2018-04-01

In this paper, we propose for the first time a method based on Bernstein polynomials for solving directly a class of third-order ordinary differential equations (ODEs). This method gives a numerical solution by converting the equation into a system of algebraic equations which is solved directly. Some numerical examples are given to show the applicability of the method.

Quadratic algebras in the noncommutative integration method of wave equation

International Nuclear Information System (INIS)

Varaksin, O.L.

1995-01-01

The paper deals with the investigation of applications of the method of noncommutative integration of linear differential equations by partial derivatives. Nontrivial example was taken for integration of three-dimensions wave equation with the use of non-Abelian quadratic algebras

Wielandt method applied to the diffusion equations discretized by finite element nodal methods

International Nuclear Information System (INIS)

Mugica R, A.; Valle G, E. del

2003-01-01

Nowadays the numerical methods of solution to the diffusion equation by means of algorithms and computer programs result so extensive due to the great number of routines and calculations that should carry out, this rebounds directly in the execution times of this programs, being obtained results in relatively long times. This work shows the application of an acceleration method of the convergence of the classic method of those powers that it reduces notably the number of necessary iterations for to obtain reliable results, what means that the compute times they see reduced in great measure. This method is known in the literature like Wielandt method and it has incorporated to a computer program that is based on the discretization of the neutron diffusion equations in plate geometry and stationary state by polynomial nodal methods. In this work the neutron diffusion equations are described for several energy groups and their discretization by means of those called physical nodal methods, being illustrated in particular the quadratic case. It is described a model problem widely described in the literature which is solved for the physical nodal grade schemes 1, 2, 3 and 4 in three different ways: to) with the classic method of the powers, b) method of the powers with the Wielandt acceleration and c) method of the powers with the Wielandt modified acceleration. The results for the model problem as well as for two additional problems known as benchmark problems are reported. Such acceleration method can also be implemented to problems of different geometry to the proposal in this work, besides being possible to extend their application to problems in 2 or 3 dimensions. (Author)

A Laplace transform certified reduced basis method; application to the heat equation and wave equation

OpenAIRE

Knezevic, David; Patera, Anthony T.; Huynh, Dinh Bao Phuong

2010-01-01

We present a certified reduced basis (RB) method for the heat equation and wave equation. The critical ingredients are certified RB approximation of the Laplace transform; the inverse Laplace transform to develop the time-domain RB output approximation and rigorous error bound; a (Butterworth) filter in time to effect the necessary “modal” truncation; RB eigenfunction decomposition and contour integration for Offline–Online decomposition. We present numerical results to demonstrate the accura...

Finite element and discontinuous Galerkin methods for transient wave equations

CERN Document Server

Cohen, Gary

2017-01-01

This monograph presents numerical methods for solving transient wave equations (i.e. in time domain). More precisely, it provides an overview of continuous and discontinuous finite element methods for these equations, including their implementation in physical models, an extensive description of 2D and 3D elements with different shapes, such as prisms or pyramids, an analysis of the accuracy of the methods and the study of the Maxwell’s system and the important problem of its spurious free approximations. After recalling the classical models, i.e. acoustics, linear elastodynamics and electromagnetism and their variational formulations, the authors present a wide variety of finite elements of different shapes useful for the numerical resolution of wave equations. Then, they focus on the construction of efficient continuous and discontinuous Galerkin methods and study their accuracy by plane wave techniques and a priori error estimates. A chapter is devoted to the Maxwell’s system and the important problem ...

An improved front tracking method for the Euler equations

NARCIS (Netherlands)

Witteveen, J.A.S.; Koren, B.; Bakker, P.G.

2007-01-01

An improved front tracking method for hyperbolic conservation laws is presented. The improved method accurately resolves discontinuities as well as continuous phenomena. The method is based on an improved front interaction model for a physically more accurate modeling of the Euler equations, as

Approximate Solutions of Nonlinear Partial Differential Equations by Modified q-Homotopy Analysis Method

Directory of Open Access Journals (Sweden)

Shaheed N. Huseen

2013-01-01

Full Text Available A modified q-homotopy analysis method (mq-HAM was proposed for solving nth-order nonlinear differential equations. This method improves the convergence of the series solution in the nHAM which was proposed in (see Hassan and El-Tawil 2011, 2012. The proposed method provides an approximate solution by rewriting the nth-order nonlinear differential equation in the form of n first-order differential equations. The solution of these n differential equations is obtained as a power series solution. This scheme is tested on two nonlinear exactly solvable differential equations. The results demonstrate the reliability and efficiency of the algorithm developed.

Equation-Method for correcting clipping errors in OFDM signals.

Science.gov (United States)

Bibi, Nargis; Kleerekoper, Anthony; Muhammad, Nazeer; Cheetham, Barry

2016-01-01

Orthogonal frequency division multiplexing (OFDM) is the digital modulation technique used by 4G and many other wireless communication systems. OFDM signals have significant amplitude fluctuations resulting in high peak to average power ratios which can make an OFDM transmitter susceptible to non-linear distortion produced by its high power amplifiers (HPA). A simple and popular solution to this problem is to clip the peaks before an OFDM signal is applied to the HPA but this causes in-band distortion and introduces bit-errors at the receiver. In this paper we discuss a novel technique, which we call the Equation-Method, for correcting these errors. The Equation-Method uses the Fast Fourier Transform to create a set of simultaneous equations which, when solved, return the amplitudes of the peaks before they were clipped. We show analytically and through simulations that this method can, correct all clipping errors over a wide range of clipping thresholds. We show that numerical instability can be avoided and new techniques are needed to enable the receiver to differentiate between correctly and incorrectly received frequency-domain constellation symbols.

Nodal integral method for the neutron diffusion equation in cylindrical geometry

International Nuclear Information System (INIS)

Azmy, Y.Y.

1987-01-01

The nodal methodology is based on retaining a higher a higher degree of analyticity in the process of deriving the discrete-variable equations compared to conventional numerical methods. As a result, extensive numerical testing of nodal methods developed for a wide variety of partial differential equations and comparison of the results to conventional methods have established the superior accuracy of nodal methods on coarse meshes. Moreover, these tests have shown that nodal methods are more computationally efficient than finite difference and finite-element methods in the sense that they require shorter CPU times to achieve comparable accuracy in the solutions. However, nodal formalisms and the final discrete-variable equations they produce are, in general, more complicated than their conventional counterparts. This, together with anticipated difficulties in applying the transverse-averaging procedure in curvilinear coordinates, has limited the applications of nodal methods, so far, to Cartesian geometry, and with additional approximations to hexagonal geometry. In this paper the authors report recent progress in deriving and numerically implementing a nodal integral method (NIM) for solving the neutron diffusion equation in cylindrical r-z geometry. Also, presented are comparisons of numerical solutions to two test problems with those obtained by the Exterminator-2 code, which indicate the superior accuracy of the nodal integral method solutions on much coarser meshes

A Novel Method for Analytical Solutions of Fractional Partial Differential Equations

Directory of Open Access Journals (Sweden)

Mehmet Ali Akinlar

2013-01-01

Full Text Available A new solution technique for analytical solutions of fractional partial differential equations (FPDEs is presented. The solutions are expressed as a finite sum of a vector type functional. By employing MAPLE software, it is shown that the solutions might be extended to an arbitrary degree which makes the present method not only different from the others in the literature but also quite efficient. The method is applied to special Bagley-Torvik and Diethelm fractional differential equations as well as a more general fractional differential equation.

Approximate analytical methods for solving ordinary differential equations

CERN Document Server

Radhika, TSL; Rani, T Raja

2015-01-01

Approximate Analytical Methods for Solving Ordinary Differential Equations (ODEs) is the first book to present all of the available approximate methods for solving ODEs, eliminating the need to wade through multiple books and articles. It covers both well-established techniques and recently developed procedures, including the classical series solution method, diverse perturbation methods, pioneering asymptotic methods, and the latest homotopy methods.The book is suitable not only for mathematicians and engineers but also for biologists, physicists, and economists. It gives a complete descripti

Analysis of a time fractional wave-like equation with the homotopy analysis method

International Nuclear Information System (INIS)

Xu Hang; Cang Jie

2008-01-01

The time fractional wave-like differential equation with a variable coefficient is studied analytically. By using a simple transformation, the governing equation is reduced to two fractional ordinary differential equations. Then the homotopy analysis method is employed to derive the solutions of these equations. The accurate series solutions are obtained. Especially, when h f =h g =-1, these solutions are exactly the same as those results given by the Adomian decomposition method. The present work shows the validity and great potential of the homotopy analysis method for solving nonlinear fractional differential equations. The basic idea described in this Letter is expected to be further employed to solve other similar nonlinear problems in fractional calculus

Inverse operator method for solutions of nonlinear dynamical equations and some typical applications

International Nuclear Information System (INIS)

Fang Jinqing; Yao Weiguang

1993-01-01

The inverse operator method (IOM) is described briefly. We have realized the IOM for the solutions of nonlinear dynamical equations by the mathematics-mechanization (MM) with computers. They can then offer a new and powerful method applicable to many areas of physics. We have applied them successfully to study the chaotic behaviors of some nonlinear dynamical equations. As typical examples, the well-known Lorentz equation, generalized Duffing equation and two coupled generalized Duffing equations are investigated by using the IOM and the MM. The results are in good agreement with those given by Runge-Kutta method. So the IOM realized by the MM is of potential application valuable in nonlinear physics and many other fields

Fast pressure-correction method for incompressible Navier-Stokes equations in curvilinear coordinates

Science.gov (United States)

Aithal, Abhiram; Ferrante, Antonino

2017-11-01

In order to perform direct numerical simulations (DNS) of turbulent flows over curved surfaces and axisymmetric bodies, we have developed the numerical methodology to solve the incompressible Navier-Stokes (NS) equations in curvilinear coordinates for orthogonal meshes. The orthogonal meshes are generated by solving a coupled system of non-linear Poisson equations. The NS equations in orthogonal curvilinear coordinates are discretized in space on a staggered mesh using second-order central-difference scheme and are solved with an FFT-based pressure-correction method. The momentum equation is integrated in time using the second-order Adams-Bashforth scheme. The velocity field is advanced in time by applying the pressure correction to the approximate velocity such that it satisfies the divergence free condition. The novelty of the method stands in solving the variable coefficient Poisson equation for pressure using an FFT-based Poisson solver rather than the slower multigrid methods. We present the verification and validation results of the new numerical method and the DNS results of transitional flow over a curved axisymmetric body.

Stability of numerical method for semi-linear stochastic pantograph differential equations

Directory of Open Access Journals (Sweden)

Yu Zhang

2016-01-01

Full Text Available Abstract As a particular expression of stochastic delay differential equations, stochastic pantograph differential equations have been widely used in nonlinear dynamics, quantum mechanics, and electrodynamics. In this paper, we mainly study the stability of analytical solutions and numerical solutions of semi-linear stochastic pantograph differential equations. Some suitable conditions for the mean-square stability of an analytical solution are obtained. Then we proved the general mean-square stability of the exponential Euler method for a numerical solution of semi-linear stochastic pantograph differential equations, that is, if an analytical solution is stable, then the exponential Euler method applied to the system is mean-square stable for arbitrary step-size h > 0 $h>0$ . Numerical examples further illustrate the obtained theoretical results.

Discontinuous Galerkin Subgrid Finite Element Method for Heterogeneous Brinkman’s Equations

KAUST Repository

Iliev, Oleg P.

2010-01-01

We present a two-scale finite element method for solving Brinkman\\'s equations with piece-wise constant coefficients. This system of equations model fluid flows in highly porous, heterogeneous media with complex topology of the heterogeneities. We make use of the recently proposed discontinuous Galerkin FEM for Stokes equations by Wang and Ye in [12] and the concept of subgrid approximation developed for Darcy\\'s equations by Arbogast in [4]. In order to reduce the error along the coarse-grid interfaces we have added a alternating Schwarz iteration using patches around the coarse-grid boundaries. We have implemented the subgrid method using Deal.II FEM library, [7], and we present the computational results for a number of model problems. © 2010 Springer-Verlag Berlin Heidelberg.

Introduction to partial differential equations and Hilbert space methods

CERN Document Server

Gustafson, Karl E

1997-01-01

Easy-to-use text examines principal method of solving partial differential equations, 1st-order systems, computation methods, and much more. Over 600 exercises, with answers for many. Ideal for a 1-semester or full-year course.

A multi-domain spectral method for time-fractional differential equations

Science.gov (United States)

Chen, Feng; Xu, Qinwu; Hesthaven, Jan S.

2015-07-01

This paper proposes an approach for high-order time integration within a multi-domain setting for time-fractional differential equations. Since the kernel is singular or nearly singular, two main difficulties arise after the domain decomposition: how to properly account for the history/memory part and how to perform the integration accurately. To address these issues, we propose a novel hybrid approach for the numerical integration based on the combination of three-term-recurrence relations of Jacobi polynomials and high-order Gauss quadrature. The different approximations used in the hybrid approach are justified theoretically and through numerical examples. Based on this, we propose a new multi-domain spectral method for high-order accurate time integrations and study its stability properties by identifying the method as a generalized linear method. Numerical experiments confirm hp-convergence for both time-fractional differential equations and time-fractional partial differential equations.

New analytical exact solutions of time fractional KdV-KZK equation by Kudryashov methods

Science.gov (United States)

S Saha, Ray

2016-04-01

In this paper, new exact solutions of the time fractional KdV-Khokhlov-Zabolotskaya-Kuznetsov (KdV-KZK) equation are obtained by the classical Kudryashov method and modified Kudryashov method respectively. For this purpose, the modified Riemann-Liouville derivative is used to convert the nonlinear time fractional KdV-KZK equation into the nonlinear ordinary differential equation. In the present analysis, the classical Kudryashov method and modified Kudryashov method are both used successively to compute the analytical solutions of the time fractional KdV-KZK equation. As a result, new exact solutions involving the symmetrical Fibonacci function, hyperbolic function and exponential function are obtained for the first time. The methods under consideration are reliable and efficient, and can be used as an alternative to establish new exact solutions of different types of fractional differential equations arising from mathematical physics. The obtained results are exhibited graphically in order to demonstrate the efficiencies and applicabilities of these proposed methods of solving the nonlinear time fractional KdV-KZK equation.

Series: Utilization of Differential Equations and Methods for Solving Them in Medical Physics (3).

Science.gov (United States)

Murase, Kenya

2016-01-01

In this issue, simultaneous differential equations were introduced. These differential equations are often used in the field of medical physics. The methods for solving them were also introduced, which include Laplace transform and matrix methods. Some examples were also introduced, in which Laplace transform and matrix methods were applied to solving simultaneous differential equations derived from a three-compartment kinetic model for analyzing the glucose metabolism in tissues and Bloch equations for describing the behavior of the macroscopic magnetization in magnetic resonance imaging.In the next (final) issue, partial differential equations and various methods for solving them will be introduced together with some examples in medical physics.

A Hennart nodal method for the diffusion equation

International Nuclear Information System (INIS)

Lesaint, P.; Noceir, S.; Verwaerde, D.

1995-01-01

A modification of the Hennart nodal method for neutron diffusion problems is presented. The final system of equations obtained by this method is not positive definite. However, a flux elimination technique leads to a simple positive definite system, which can be solved by the traditional iterative methods. Calculations of a two-dimensional International Atomic Energy Agency benchmark problem are performed and compared with results of the original Hennart nodal method and some finite element methods. The high computational efficiency of this modified nodal method is clearly demonstrated

Non-standard finite difference and Chebyshev collocation methods for solving fractional diffusion equation

Science.gov (United States)

Agarwal, P.; El-Sayed, A. A.

2018-06-01

In this paper, a new numerical technique for solving the fractional order diffusion equation is introduced. This technique basically depends on the Non-Standard finite difference method (NSFD) and Chebyshev collocation method, where the fractional derivatives are described in terms of the Caputo sense. The Chebyshev collocation method with the (NSFD) method is used to convert the problem into a system of algebraic equations. These equations solved numerically using Newton's iteration method. The applicability, reliability, and efficiency of the presented technique are demonstrated through some given numerical examples.

Matrix Methods for Solving Hartree-Fock Equations in Atomic Structure Calculations and Line Broadening

Directory of Open Access Journals (Sweden)

Thomas Gomez

2018-04-01

Full Text Available Atomic structure of N-electron atoms is often determined by solving the Hartree-Fock equations, which are a set of integro-differential equations. The integral part of the Hartree-Fock equations treats electron exchange, but the Hartree-Fock equations are not often treated as an integro-differential equation. The exchange term is often approximated as an inhomogeneous or an effective potential so that the Hartree-Fock equations become a set of ordinary differential equations (which can be solved using the usual shooting methods. Because the Hartree-Fock equations are an iterative-refinement method, the inhomogeneous term relies on the previous guess of the wavefunction. In addition, there are numerical complications associated with solving inhomogeneous differential equations. This work uses matrix methods to solve the Hartree-Fock equations as an integro-differential equation. It is well known that a derivative operator can be expressed as a matrix made of finite-difference coefficients; energy eigenvalues and eigenvectors can be obtained by using linear-algebra packages. The integral (exchange part of the Hartree-Fock equation can be approximated as a sum and written as a matrix. The Hartree-Fock equations can be solved as a matrix that is the sum of the differential and integral matrices. We compare calculations using this method against experiment and standard atomic structure calculations. This matrix method can also be used to solve for free-electron wavefunctions, thus improving how the atoms and free electrons interact. This technique is important for spectral line broadening in two ways: it improves the atomic structure calculations, and it improves the motion of the plasma electrons that collide with the atom.

A direct algebraic method applied to obtain complex solutions of some nonlinear partial differential equations

International Nuclear Information System (INIS)

Zhang Huiqun

2009-01-01

By using some exact solutions of an auxiliary ordinary differential equation, a direct algebraic method is described to construct the exact complex solutions for nonlinear partial differential equations. The method is implemented for the NLS equation, a new Hamiltonian amplitude equation, the coupled Schrodinger-KdV equations and the Hirota-Maccari equations. New exact complex solutions are obtained.

A new integral method for solving the point reactor neutron kinetics equations

International Nuclear Information System (INIS)

Li Haofeng; Chen Wenzhen; Luo Lei; Zhu Qian

2009-01-01

A numerical integral method that efficiently provides the solution of the point kinetics equations by using the better basis function (BBF) for the approximation of the neutron density in one time step integrations is described and investigated. The approach is based on an exact analytic integration of the neutron density equation, where the stiffness of the equations is overcome by the fully implicit formulation. The procedure is tested by using a variety of reactivity functions, including step reactivity insertion, ramp input and oscillatory reactivity changes. The solution of the better basis function method is compared to other analytical and numerical solutions of the point reactor kinetics equations. The results show that selecting a better basis function can improve the efficiency and accuracy of this integral method. The better basis function method can be used in real time forecasting for power reactors in order to prevent reactivity accidents.

Flux weighted method for solution of stiff neutron dynamic equations and its application

International Nuclear Information System (INIS)

Li Huiyun; Jiao Huixian

1987-12-01

To analyze reactivity event for nuclear power plants, it is necessary to solve the neutron dynamic equations, which is a group of typical stiff constant differential equations. Very small time steps could only be adopted when the group of equations is solved by common methods. However, a large time steps might be selected if the Flux Weighted Medthod introduced in this paper is used. Generally, weighted factor θ i1 is set as a constant. Naturally, this treatment method can decrease the accuracy of calculation for the increase of the steadiness of solving the equations. An accurate theoretical formula of 4 x 4 matrix of θ i1 is rigorously derived so that the accuracy of calculation is ensured, as well as the steadiness of solved equations is increased. This method have the advantage over classical Runge-kutta Method and other methods. The time steps could be increased by a factor of 1 ∼ 3 orders of magnitude so as to save a lot of computating time. The programe solving neutron dynamic equation, which is prepared by using Flux Weighted Method, could be sued for real time analog of training simulator, as well as for analysis and computation of reactivity event (including rod jumping out event)

Derivation of equation of quasipotential type using the method of Fock-- Podolsky

Energy Technology Data Exchange (ETDEWEB)

Blokhintsev, D I; Rizov, V A; Todorov, I T

1975-12-31

A quasipotential equation is derived for the relativistic Coulomb problem from the equations of motion of quantum electrodynamics using the method of Fock-- Podolsky (Tamm-Dancoff). Relation with an inhomogeneous equation for the 4-point retarded function is exhibited. (auth)

Application of the Exp-function method to the equal-width wave equation

International Nuclear Information System (INIS)

Biazar, J; Ayati, Z

2008-01-01

In this paper, the Exp-function method is used to find an exact solution of the equal-width wave (EW) equation. The method is straightforward and concise, and its applications are promising. It is shown that the Exp-function method, with the help of symbolic computation, provides a very effective and powerful mathematical tool for solving the EW equation.

Transport methods: general. 3. An Additive Angular-Dependent Re-balance Acceleration Method for Neutron Transport Equations

International Nuclear Information System (INIS)

Cho, Nam Zin; Park, Chang Je

2001-01-01

An additive angular-dependent re-balance (AADR) factor acceleration method is described to accelerate the source iteration of discrete ordinates transport calculation. The formulation of the AADR method follows that of the angular-dependent re-balance (ADR) method in that the re-balance factor is defined only on the cell interface and in that the low-order equation is derived by integrating the transport equation (high-order equation) over angular subspaces. But, the re-balance factor is applied additively. While the AADR method is similar to the boundary projection acceleration and the alpha-weighted linear acceleration, it is more general and does have distinct features. The method is easily extendible to DP N and low-order S N re-balancing, and it does not require consistent discretizations between the high- and low-order equations as in diffusion synthetic acceleration. We find by Fourier analysis and numerical results that the AADR method with a chosen form of weighting functions is unconditionally stable and very effective. There also exists an optimal weighting parameter that leads to the smallest spectral radius. The AADR acceleration method described in this paper is simple to implement, unconditionally stable, and very effective. It uses a physically based weighting function with an optimal parameter, leading to the best spectral radius of ρ<0.1865, compared to ρ<0.2247 of DSA. The application of the AADR acceleration method with the LMB scheme on a test problem shows encouraging results

Application of the finite element method to the neutron transport equation

International Nuclear Information System (INIS)

Martin, W.R.

1976-01-01

This paper examines the theoretical and practical application of the finite element method to the neutron transport equation. It is shown that in principle the system of equations obtained by application of the finite element method can be solved with certain physical restrictions concerning the criticality of the medium. The convergence of this approximate solution to the exact solution with mesh refinement is examined, and a non-optical estimate of the convergence rate is obtained analytically. It is noted that the numerical results indicate a faster convergence rate and several approaches to obtain this result analytically are outlined. The practical application of the finite element method involved the development of a computer code capable of solving the neutron transport equation in 1-D plane geometry. Vacuum, reflecting, or specified incoming boundary conditions may be analyzed, and all are treated as natural boundary conditions. The time-dependent transport equation is also examined and it is shown that the application of the finite element method in conjunction with the Crank-Nicholson time discretization method results in a system of algebraic equations which is readily solved. Numerical results are given for several critical slab eigenvalue problems, including anisotropic scattering, and the results compare extremely well with benchmark results. It is seen that the finite element code is more efficient than a standard discrete ordinates code for certain problems. A problem with severe heterogeneities is considered and it is shown that the use of discontinuous spatial and angular elements results in a marked improvement in the results. Finally, time-dependent problems are examined and it is seen that the phenomenon of angular mode separation makes the numerical treatment of the transport equation in slab geometry a considerable challenge, with the result that the angular mesh has a dominant effect on obtaining acceptable solutions

On the evolution equations, solvable through the inverse scattering method

International Nuclear Information System (INIS)

Gerdjikov, V.S.; Khristov, E.Kh.

1979-01-01

The nonlinear evolution equations (NLEE), related to the one-parameter family of Dirac operators are considered in a uniform manner. The class of NLEE solvable through the inverse scatterina method and their conservation laws are described. The description of the hierarchy of Hamiltonian structures and the proof of complete integrability of the NLEE is presented. The class of Baecklund transformations for these NLEE is derived. The general formulae are illustrated by two important examples: the nonlinear Schroedinger equation and the sine-Gordon equation

Application of an analytical method for solution of thermal hydraulic conservation equations

Energy Technology Data Exchange (ETDEWEB)

Fakory, M.R. [Simulation, Systems & Services Technologies Company (S3 Technologies), Columbia, MD (United States)

1995-09-01

An analytical method has been developed and applied for solution of two-phase flow conservation equations. The test results for application of the model for simulation of BWR transients are presented and compared with the results obtained from application of the explicit method for integration of conservation equations. The test results show that with application of the analytical method for integration of conservation equations, the Courant limitation associated with explicit Euler method of integration was eliminated. The results obtained from application of the analytical method (with large time steps) agreed well with the results obtained from application of explicit method of integration (with time steps smaller than the size imposed by Courant limitation). The results demonstrate that application of the analytical approach significantly improves the numerical stability and computational efficiency.

Fibonacci-regularization method for solving Cauchy integral equations of the first kind

Directory of Open Access Journals (Sweden)

Mohammad Ali Fariborzi Araghi

2017-09-01

Full Text Available In this paper, a novel scheme is proposed to solve the first kind Cauchy integral equation over a finite interval. For this purpose, the regularization method is considered. Then, the collocation method with Fibonacci base function is applied to solve the obtained second kind singular integral equation. Also, the error estimate of the proposed scheme is discussed. Finally, some sample Cauchy integral equations stem from the theory of airfoils in fluid mechanics are presented and solved to illustrate the importance and applicability of the given algorithm. The tables in the examples show the efficiency of the method.

The Telegraph Equation and Its Solution by Reduced Differential Transform Method

Directory of Open Access Journals (Sweden)

Vineet K. Srivastava

2013-01-01

Full Text Available One-dimensional second-order hyperbolic telegraph equation was formulated using Ohm’s law and solved by a recent and reliable semianalytic method, namely, the reduced differential transform method (RDTM. Using this method, it is possible to find the exact solution or a closed approximate solution of a differential equation. Three numerical examples have been carried out in order to check the effectiveness, the accuracy, and convergence of the method. The RDTM is a powerful mathematical technique for solving wide range of problems arising in science and engineering fields.

Study on state equation for hydrogen storage measurement by volumetric method

International Nuclear Information System (INIS)

Dai Wei; Xu Jiajing; Wang Chaoyang; Tang Yongjian

2014-01-01

Volumetric measurement technique is one of the most popular methods for determining the amount of hydrogen storage. A new state equation was established which extended the limitations from the ideal gas state equation, the van der Waals equation and the Gou equation. The new state equation was then employed to describe the p-V-T character of hydrogen and investigate the adsorption quantity of hydrogen storage in resorcin-formaldehyde aerogel under different temperatures and pressures. The new equation was used to describe the density of hydrogen under different temperatures and pressures. The results are in good agreement with the experimental data. The differences arising from various underlying physics were carefully analyzed. (authors)

Modified Splitting FDTD Methods for Two-Dimensional Maxwell’s Equations

Directory of Open Access Journals (Sweden)

Liping Gao

2017-01-01

Full Text Available In this paper, we develop a new method to reduce the error in the splitting finite-difference method of Maxwell’s equations. By this method two modified splitting FDTD methods (MS-FDTDI, MS-FDTDII for the two-dimensional Maxwell equations are proposed. It is shown that the two methods are second-order accurate in time and space and unconditionally stable by Fourier methods. By energy method, it is proved that MS-FDTDI is second-order convergent. By deriving the numerical dispersion (ND relations, we prove rigorously that MS-FDTDI has less ND errors than the ADI-FDTD method and the ND errors of ADI-FDTD are less than those of MS-FDTDII. Numerical experiments for computing ND errors and simulating a wave guide problem and a scattering problem are carried out and the efficiency of the MS-FDTDI and MS-FDTDII methods is confirmed.

Iterative methods for compressible Navier-Stokes and Euler equations

Energy Technology Data Exchange (ETDEWEB)

Tang, W.P.; Forsyth, P.A.

1996-12-31

This workshop will focus on methods for solution of compressible Navier-Stokes and Euler equations. In particular, attention will be focused on the interaction between the methods used to solve the non-linear algebraic equations (e.g. full Newton or first order Jacobian) and the resulting large sparse systems. Various types of block and incomplete LU factorization will be discussed, as well as stability issues, and the use of Newton-Krylov methods. These techniques will be demonstrated on a variety of model transonic and supersonic airfoil problems. Applications to industrial CFD problems will also be presented. Experience with the use of C++ for solution of large scale problems will also be discussed. The format for this workshop will be four fifteen minute talks, followed by a roundtable discussion.

Discontinuous Galerkin Subgrid Finite Element Method for Heterogeneous Brinkman’s Equations

KAUST Repository

Iliev, Oleg P.; Lazarov, Raytcho D.; Willems, Joerg

2010-01-01

We present a two-scale finite element method for solving Brinkman's equations with piece-wise constant coefficients. This system of equations model fluid flows in highly porous, heterogeneous media with complex topology of the heterogeneities. We

Simulation electromagnetic scattering on bodies through integral equation and neural networks methods

Science.gov (United States)

Lvovich, I. Ya; Preobrazhenskiy, A. P.; Choporov, O. N.

2018-05-01

The paper deals with the issue of electromagnetic scattering on a perfectly conducting diffractive body of a complex shape. Performance calculation of the body scattering is carried out through the integral equation method. Fredholm equation of the second time was used for calculating electric current density. While solving the integral equation through the moments method, the authors have properly described the core singularity. The authors determined piecewise constant functions as basic functions. The chosen equation was solved through the moments method. Within the Kirchhoff integral approach it is possible to define the scattered electromagnetic field, in some way related to obtained electrical currents. The observation angles sector belongs to the area of the front hemisphere of the diffractive body. To improve characteristics of the diffractive body, the authors used a neural network. All the neurons contained a logsigmoid activation function and weighted sums as discriminant functions. The paper presents the matrix of weighting factors of the connectionist model, as well as the results of the optimized dimensions of the diffractive body. The paper also presents some basic steps in calculation technique of the diffractive bodies, based on the combination of integral equation and neural networks methods.

New analytical exact solutions of time fractional KdV–KZK equation by Kudryashov methods

International Nuclear Information System (INIS)

Saha Ray, S

2016-01-01

In this paper, new exact solutions of the time fractional KdV–Khokhlov–Zabolotskaya–Kuznetsov (KdV–KZK) equation are obtained by the classical Kudryashov method and modified Kudryashov method respectively. For this purpose, the modified Riemann–Liouville derivative is used to convert the nonlinear time fractional KdV–KZK equation into the nonlinear ordinary differential equation. In the present analysis, the classical Kudryashov method and modified Kudryashov method are both used successively to compute the analytical solutions of the time fractional KdV–KZK equation. As a result, new exact solutions involving the symmetrical Fibonacci function, hyperbolic function and exponential function are obtained for the first time. The methods under consideration are reliable and efficient, and can be used as an alternative to establish new exact solutions of different types of fractional differential equations arising from mathematical physics. The obtained results are exhibited graphically in order to demonstrate the efficiencies and applicabilities of these proposed methods of solving the nonlinear time fractional KdV–KZK equation. (paper)

Solving Fokker-Planck Equations on Cantor Sets Using Local Fractional Decomposition Method

Directory of Open Access Journals (Sweden)

Shao-Hong Yan

2014-01-01

Full Text Available The local fractional decomposition method is applied to approximate the solutions for Fokker-Planck equations on Cantor sets with local fractional derivative. The obtained results give the present method that is very effective and simple for solving the differential equations on Cantor set.

Spectral methods for time dependent partial differential equations

Science.gov (United States)

Gottlieb, D.; Turkel, E.

1983-01-01

The theory of spectral methods for time dependent partial differential equations is reviewed. When the domain is periodic Fourier methods are presented while for nonperiodic problems both Chebyshev and Legendre methods are discussed. The theory is presented for both hyperbolic and parabolic systems using both Galerkin and collocation procedures. While most of the review considers problems with constant coefficients the extension to nonlinear problems is also discussed. Some results for problems with shocks are presented.

An Approximate Method for the Acoustic Attenuating VTI Eikonal Equation

KAUST Repository

Hao, Q.

2017-05-26

We present an approximate method to solve the acoustic eikonal equation for attenuating transversely isotropic media with a vertical symmetry axis (VTI). A perturbation method is used to derive the perturbation formula for complex-valued traveltimes. The application of Shanks transform further enhances the accuracy of approximation. We derive both analytical and numerical solutions to the acoustic eikonal equation. The analytic solution is valid for homogeneous VTI media with moderate anellipticity and strong attenuation and attenuation-anisotropy. The numerical solution is applicable for inhomogeneous attenuating VTI media.

An Approximate Method for the Acoustic Attenuating VTI Eikonal Equation

KAUST Repository

Hao, Q.; Alkhalifah, Tariq Ali

2017-01-01

We present an approximate method to solve the acoustic eikonal equation for attenuating transversely isotropic media with a vertical symmetry axis (VTI). A perturbation method is used to derive the perturbation formula for complex-valued traveltimes. The application of Shanks transform further enhances the accuracy of approximation. We derive both analytical and numerical solutions to the acoustic eikonal equation. The analytic solution is valid for homogeneous VTI media with moderate anellipticity and strong attenuation and attenuation-anisotropy. The numerical solution is applicable for inhomogeneous attenuating VTI media.

An Unconditionally Stable Method for Solving the Acoustic Wave Equation

Directory of Open Access Journals (Sweden)

Zhi-Kai Fu

2015-01-01

Full Text Available An unconditionally stable method for solving the time-domain acoustic wave equation using Associated Hermit orthogonal functions is proposed. The second-order time derivatives in acoustic wave equation are expanded by these orthogonal basis functions. By applying Galerkin temporal testing procedure, the time variable can be eliminated from the calculations. The restriction of Courant-Friedrichs-Levy (CFL condition in selecting time step for analyzing thin layer can be avoided. Numerical results show the accuracy and the efficiency of the proposed method.

Conservation properties of numerical integration methods for systems of ordinary differential equations

Science.gov (United States)

Rosenbaum, J. S.

1976-01-01

If a system of ordinary differential equations represents a property conserving system that can be expressed linearly (e.g., conservation of mass), it is then desirable that the numerical integration method used conserve the same quantity. It is shown that both linear multistep methods and Runge-Kutta methods are 'conservative' and that Newton-type methods used to solve the implicit equations preserve the inherent conservation of the numerical method. It is further shown that a method used by several authors is not conservative.

An alternative solver for the nodal expansion method equations - 106

International Nuclear Information System (INIS)

Carvalho da Silva, F.; Carlos Marques Alvim, A.; Senra Martinez, A.

2010-01-01

An automated procedure for nuclear reactor core design is accomplished by using a quick and accurate 3D nodal code, aiming at solving the diffusion equation, which describes the spatial neutron distribution in the reactor. This paper deals with an alternative solver for nodal expansion method (NEM), with only two inner iterations (mesh sweeps) per outer iteration, thus having the potential to reduce the time required to calculate the power distribution in nuclear reactors, but with accuracy similar to the ones found in conventional NEM. The proposed solver was implemented into a computational system which, besides solving the diffusion equation, also solves the burnup equations governing the gradual changes in material compositions of the core due to fuel depletion. Results confirm the effectiveness of the method for practical purposes. (authors)

Sufficient Descent Conjugate Gradient Methods for Solving Convex Constrained Nonlinear Monotone Equations

Directory of Open Access Journals (Sweden)

San-Yang Liu

2014-01-01

Full Text Available Two unified frameworks of some sufficient descent conjugate gradient methods are considered. Combined with the hyperplane projection method of Solodov and Svaiter, they are extended to solve convex constrained nonlinear monotone equations. Their global convergence is proven under some mild conditions. Numerical results illustrate that these methods are efficient and can be applied to solve large-scale nonsmooth equations.

International Conference on Multiscale Methods and Partial Differential Equations.

Energy Technology Data Exchange (ETDEWEB)

Thomas Hou

2006-12-12

The International Conference on Multiscale Methods and Partial Differential Equations (ICMMPDE for short) was held at IPAM, UCLA on August 26-27, 2005. The conference brought together researchers, students and practitioners with interest in the theoretical, computational and practical aspects of multiscale problems and related partial differential equations. The conference provided a forum to exchange and stimulate new ideas from different disciplines, and to formulate new challenging multiscale problems that will have impact in applications.

Application of Legendre spectral-collocation method to delay differential and stochastic delay differential equation

Science.gov (United States)

Khan, Sami Ullah; Ali, Ishtiaq

2018-03-01

Explicit solutions to delay differential equation (DDE) and stochastic delay differential equation (SDDE) can rarely be obtained, therefore numerical methods are adopted to solve these DDE and SDDE. While on the other hand due to unstable nature of both DDE and SDDE numerical solutions are also not straight forward and required more attention. In this study, we derive an efficient numerical scheme for DDE and SDDE based on Legendre spectral-collocation method, which proved to be numerical methods that can significantly speed up the computation. The method transforms the given differential equation into a matrix equation by means of Legendre collocation points which correspond to a system of algebraic equations with unknown Legendre coefficients. The efficiency of the proposed method is confirmed by some numerical examples. We found that our numerical technique has a very good agreement with other methods with less computational effort.

Splitting Method for Solving the Coarse-Mesh Discretized Low-Order Quasi-Diffusion Equations

International Nuclear Information System (INIS)

Hiruta, Hikaru; Anistratov, Dmitriy Y.; Adams, Marvin L.

2005-01-01

In this paper, the development is presented of a splitting method that can efficiently solve coarse-mesh discretized low-order quasi-diffusion (LOQD) equations. The LOQD problem can reproduce exactly the transport scalar flux and current. To solve the LOQD equations efficiently, a splitting method is proposed. The presented method splits the LOQD problem into two parts: (a) the D problem that captures a significant part of the transport solution in the central parts of assemblies and can be reduced to a diffusion-type equation and (b) the Q problem that accounts for the complicated behavior of the transport solution near assembly boundaries. Independent coarse-mesh discretizations are applied: the D problem equations are approximated by means of a finite element method, whereas the Q problem equations are discretized using a finite volume method. Numerical results demonstrate the efficiency of the methodology presented. This methodology can be used to modify existing diffusion codes for full-core calculations (which already solve a version of the D problem) to account for transport effects

Residual power series method for fractional Sharma-Tasso-Olever equation

Directory of Open Access Journals (Sweden)

Amit Kumar

2016-02-01

Full Text Available In this paper, we introduce a modified analytical approximate technique to obtain solution of time fractional Sharma-Tasso-Olever equation. First, we present an alternative framework of the Residual power series method (RPSM which can be used simply and effectively to handle nonlinear fractional differential equations arising in several physical phenomena. This method is basically based on the generalized Taylor series formula and residual error function. A good result is found between our solution and the given solution. It is shown that the proposed method is reliable, efficient and easy to implement on all kinds of fractional nonlinear problems arising in science and technology.

Solution of the generalized Emden-Fowler equations by the hybrid functions method

International Nuclear Information System (INIS)

Tabrizidooz, H R; Marzban, H R; Razzaghi, M

2009-01-01

In this paper, we present a numerical algorithm for solving the generalized Emden-Fowler equations, which have many applications in mathematical physics and astrophysics. The method is based on hybrid functions approximations. The properties of hybrid functions, which consist of block-pulse functions and Lagrange interpolating polynomials, are presented. These properties are then utilized to reduce the computation of the generalized Emden-Fowler equations to a system of nonlinear equations. The method is easy to implement and yields very accurate results.

Cantor-type cylindrical-coordinate method for differential equations with local fractional derivatives

International Nuclear Information System (INIS)

Yang, Xiao-Jun; Srivastava, H.M.; He, Ji-Huan; Baleanu, Dumitru

2013-01-01

In this Letter, we propose to use the Cantor-type cylindrical-coordinate method in order to investigate a family of local fractional differential operators on Cantor sets. Some testing examples are given to illustrate the capability of the proposed method for the heat-conduction equation on a Cantor set and the damped wave equation in fractal strings. It is seen to be a powerful tool to convert differential equations on Cantor sets from Cantorian-coordinate systems to Cantor-type cylindrical-coordinate systems.

The Closest Point Method and Multigrid Solvers for Elliptic Equations on Surfaces

KAUST Repository

Chen, Yujia

2015-01-01

© 2015 Society for Industrial and Applied Mathematics. Elliptic partial differential equations are important from both application and analysis points of view. In this paper we apply the closest point method to solve elliptic equations on general curved surfaces. Based on the closest point representation of the underlying surface, we formulate an embedding equation for the surface elliptic problem, then discretize it using standard finite differences and interpolation schemes on banded but uniform Cartesian grids. We prove the convergence of the difference scheme for the Poisson\\'s equation on a smooth closed curve. In order to solve the resulting large sparse linear systems, we propose a specific geometric multigrid method in the setting of the closest point method. Convergence studies in both the accuracy of the difference scheme and the speed of the multigrid algorithm show that our approaches are effective.

Solution of two group neutron diffusion equation by using homotopy analysis method

International Nuclear Information System (INIS)

Cavdar, S.

2010-01-01

The Homotopy Analysis Method (HAM), proposed in 1992 by Shi Jun Liao and has been developed since then, is based on differential geometry as well as homotopy which is a fundamental concept in topology. It has proved to be useful for obtaining series solutions of many such problems involving algebraic, linear/non-linear, ordinary/partial differential equations, differential-integral equations, differential-difference equations, and coupled equations of them. Briefly, through HAM, it is possible to construct a continuous mapping of an initial guess approximation to the exact solution of the equation of concern. An auxiliary linear operator is chosen to construct such kind of a continuous mapping and an auxiliary parameter is used to ensure the convergence of series solution. We present the solutions of two-group neutron diffusion equation through HAM in this work. We also compare the results with that obtained by other well-known solution analytical and numeric methods.

Methods for the solution of the two-dimensional radiation-transfer equation

International Nuclear Information System (INIS)

Weaver, R.; Mihalas, D.; Olson, G.

1982-01-01

We use the variable Eddington factor (VEF) approximation to solve the time-dependent two-dimensional radiation transfer equation. The transfer equation and its moments are derived for an inertial frame of reference in cylindrical geometry. Using the VEF tensor to close the moment equations, we manipulate them into a combined moment equation that results in an energy equation, which is automatically flux limited. There are two separable facets in this method of solution. First, given the variable Eddington tensor, we discuss the efficient solution of the combined moment matrix equation. The second facet of the problem is the calculation of the variable Eddington tensor. Several options for this calculation, as well as physical limitations on the use of locally-calculated Eddington factors, are discussed

Dispersive optical soliton solutions for the hyperbolic and cubic-quintic nonlinear Schrödinger equations via the extended sinh-Gordon equation expansion method

Science.gov (United States)

Seadawy, Aly R.; Kumar, Dipankar; Chakrabarty, Anuz Kumar

2018-05-01

The (2+1)-dimensional hyperbolic and cubic-quintic nonlinear Schrödinger equations describe the propagation of ultra-short pulses in optical fibers of nonlinear media. By using an extended sinh-Gordon equation expansion method, some new complex hyperbolic and trigonometric functions prototype solutions for two nonlinear Schrödinger equations were derived. The acquired new complex hyperbolic and trigonometric solutions are expressed by dark, bright, combined dark-bright, singular and combined singular solitons. The obtained results are more compatible than those of other applied methods. The extended sinh-Gordon equation expansion method is a more powerful and robust mathematical tool for generating new optical solitary wave solutions for many other nonlinear evolution equations arising in the propagation of optical pulses.

Projection-iteration methods for solving nonlinear operator equations

International Nuclear Information System (INIS)

Nguyen Minh Chuong; Tran thi Lan Anh; Tran Quoc Binh

1989-09-01

In this paper, the authors investigate a nonlinear operator equation in uniformly convex Banach spaces as in metric spaces by using stationary and nonstationary generalized projection-iteration methods. Convergence theorems in the strong and weak sense were established. (author). 7 refs

Solution of fractional differential equations by using differential transform method

International Nuclear Information System (INIS)

Arikoglu, Aytac; Ozkol, Ibrahim

2007-01-01

In this study, we implement a well known transformation technique, Differential Transform Method (DTM), to the area of fractional differential equations. Theorems that never existed before are introduced with their proofs. Also numerical examples are carried out for various types of problems, including the Bagley-Torvik, Ricatti and composite fractional oscillation equations for the application of the method. The results obtained are in good agreement with the existing ones in open literature and it is shown that the technique introduced here is robust, accurate and easy to apply

Solution of fractional differential equations by using differential transform method

Energy Technology Data Exchange (ETDEWEB)

Arikoglu, Aytac [Istanbul Technical University, Faculty of Aeronautics and Astronautics, Department of Aeronautical Engineering, Maslak, TR-34469 Istanbul (Turkey); Ozkol, Ibrahim [Istanbul Technical University, Faculty of Aeronautics and Astronautics, Department of Aeronautical Engineering, Maslak, TR-34469 Istanbul (Turkey)]. E-mail: ozkol@itu.edu.tr

2007-12-15

In this study, we implement a well known transformation technique, Differential Transform Method (DTM), to the area of fractional differential equations. Theorems that never existed before are introduced with their proofs. Also numerical examples are carried out for various types of problems, including the Bagley-Torvik, Ricatti and composite fractional oscillation equations for the application of the method. The results obtained are in good agreement with the existing ones in open literature and it is shown that the technique introduced here is robust, accurate and easy to apply.

Application of the trial equation method for solving some nonlinear ...

Indian Academy of Sciences (India)

Therefore, our aim is just to find the function F. Liu has obtained a number of exact solutions to many nonlinear differential equations when F(u) is a polynomial or a rational function. ... In this study, we apply the trial equation method to seek exact solutions of the ... twice and setting the integration constant to zero, we have.

Perturbation method for periodic solutions of nonlinear jerk equations

International Nuclear Information System (INIS)

Hu, H.

2008-01-01

A Lindstedt-Poincare type perturbation method with bookkeeping parameters is presented for determining accurate analytical approximate periodic solutions of some third-order (jerk) differential equations with cubic nonlinearities. In the process of the solution, higher-order approximate angular frequencies are obtained by Newton's method. A typical example is given to illustrate the effectiveness and simplicity of the proposed method

Convergence of method of lines approximations to partial differential equations

International Nuclear Information System (INIS)

Verwer, J.G.; Sanz-Serna, J.M.

1984-01-01

Many existing numerical schemes for evolutionary problems in partial differential equations (PDEs) can be viewed as method of lines (MOL) schemes. This paper treats the convergence of one-step MOL schemes. The main purpose is to set up a general framework for a convergence analysis applicable to nonlinear problems. The stability materials for this framework are taken from the field of nonlinear stiff ODEs. In this connection, important concepts are the logarithmic matrix norm and C-stability. A nonlinear parabolic equation and the cubic Schroedinger equation are used for illustrating the ideas. (Auth.)

Approximate solution of the transport equation by methods of Galerkin type

International Nuclear Information System (INIS)

Pitkaranta, J.

1977-01-01

Questions of the existence, uniqueness, and convergence of approximate solutions of transport equations by methods of the Galerkin type (where trial and weighting functions are the same) are discussed. The results presented do not exclude the infinite-dimensional case. Two strategies can be followed in the variational approximation of the transport operator: one proceeds from the original form of the transport equation, while the other is based on the partially symmetrized equation. Both principles are discussed in this paper. The transport equation is assumed in a discretized multigroup form

Weak Second Order Explicit Stabilized Methods for Stiff Stochastic Differential Equations

KAUST Repository

Abdulle, Assyr

2013-01-01

We introduce a new family of explicit integrators for stiff Itô stochastic differential equations (SDEs) of weak order two. These numerical methods belong to the class of one-step stabilized methods with extended stability domains and do not suffer from the step size reduction faced by standard explicit methods. The family is based on the standard second order orthogonal Runge-Kutta-Chebyshev (ROCK2) methods for deterministic problems. The convergence, meansquare, and asymptotic stability properties of the methods are analyzed. Numerical experiments, including applications to nonlinear SDEs and parabolic stochastic partial differential equations are presented and confirm the theoretical results. © 2013 Society for Industrial and Applied Mathematics.

Semiconductor device simulation by a new method of solving poisson, Laplace and Schrodinger equations

International Nuclear Information System (INIS)

Sharifi, M. J.; Adibi, A.

2000-01-01

In this paper, we have extended and completed our previous work, that was introducing a new method for finite differentiation. We show the applicability of the method for solving a wide variety of equations such as poisson, Laplace and Schrodinger. These equations are fundamental to the most semiconductor device simulators. In a section, we solve the Shordinger equation by this method in several cases including the problem of finding electron concentration profile in the channel of a HEMT. In another section, we solve the Poisson equation by this method, choosing the problem of SBD as an example. Finally we solve the Laplace equation in two dimensions and as an example, we focus on the VED. In this paper, we have shown that, the method can get stable and precise results in solving all of these problems. Also the programs which have been written based on this method become considerably faster, more clear, and more abstract

New exact solutions of coupled Boussinesq–Burgers equations by Exp-function method

Directory of Open Access Journals (Sweden)

L.K. Ravi

2017-03-01

Full Text Available In the present paper, we build the new analytical exact solutions of a nonlinear differential equation, specifically, coupled Boussinesq–Burgers equations by means of Exp-function method. Then, we analyze the results by plotting the three dimensional soliton graphs for each case, which exhibit the simplicity and effectiveness of the proposed method. The primary purpose of this paper is to employ a new approach, which allows us victorious and efficient derivation of the new analytical exact solutions for the coupled Boussinesq–Burgers equations.

Application of Rational Expansion Method for Differential-Difference Equation

International Nuclear Information System (INIS)

Wang Qi

2011-01-01

In this paper, we applied the rational formal expansion method to construct a series of soliton-like and period-form solutions for nonlinear differential-difference equations. Compared with most existing methods, the proposed method not only recovers some known solutions, but also finds some new and more general solutions. The efficiency of the method can be demonstrated on Toda Lattice and Ablowitz-Ladik Lattice. (general)

Numerical solution of multi group-Two dimensional- Adjoint equation with finite element method

International Nuclear Information System (INIS)

Poursalehi, N.; Khalafi, H.; Shahriari, M.; Minoochehr

2008-01-01

Adjoint equation is used for perturbation theory in nuclear reactor design. For numerical solution of adjoint equation, usually two methods are applied. These are Finite Element and Finite Difference procedures. Usually Finite Element Procedure is chosen for solving of adjoint equation, because it is more use able in variety of geometries. In this article, Galerkin Finite Element method is discussed. This method is applied for numerical solving multi group, multi region and two dimensional (X, Y) adjoint equation. Typical reactor geometry is partitioned with triangular meshes and boundary condition for adjoint flux is considered zero. Finally, for a case of defined parameters, Finite Element Code was applied and results were compared with Citation Code

Finite element method with quadratic quadrilateral unit for solving two dimensional incompressible N-S equation

International Nuclear Information System (INIS)

Tao Ganqiang; Yu Qing; Xiao Xiao

2011-01-01

Viscous and incompressible fluid flow is important for numerous engineering mechanics problems. Because of high non linear and incompressibility for Navier-Stokes equation, it is very difficult to solve Navier-Stokes equation by numerical method. According to its characters of Navier-Stokes equation, quartic derivation controlling equation of the two dimensional incompressible Navier-Stokes equation is set up firstly. The method solves the problem for dealing with vorticity boundary and automatically meets incompressibility condition. Then Finite Element equation for Navier-Stokes equation is proposed by using quadratic quadrilateral unit with 8 nodes in which the unit function is quadratic and non linear.-Based on it, the Finite Element program of quadratic quadrilateral unit with 8 nodes is developed. Lastly, numerical experiment proves the accuracy and dependability of the method and also shows the method has good application prospect in computational fluid mechanics. (authors)

Solution of Fractional Partial Differential Equations in Fluid Mechanics by Extension of Some Iterative Method

Directory of Open Access Journals (Sweden)

A. A. Hemeda

2013-01-01

Full Text Available An extension of the so-called new iterative method (NIM has been used to handle linear and nonlinear fractional partial differential equations. The main property of the method lies in its flexibility and ability to solve nonlinear equations accurately and conveniently. Therefore, a general framework of the NIM is presented for analytical treatment of fractional partial differential equations in fluid mechanics. The fractional derivatives are described in the Caputo sense. Numerical illustrations that include the fractional wave equation, fractional Burgers equation, fractional KdV equation, fractional Klein-Gordon equation, and fractional Boussinesq-like equation are investigated to show the pertinent features of the technique. Comparison of the results obtained by the NIM with those obtained by both Adomian decomposition method (ADM and the variational iteration method (VIM reveals that the NIM is very effective and convenient. The basic idea described in this paper is expected to be further employed to solve other similar linear and nonlinear problems in fractional calculus.

On the economical solution method for a system of linear algebraic equations

Directory of Open Access Journals (Sweden)

Jan Awrejcewicz

2004-01-01

Full Text Available The present work proposes a novel optimal and exact method of solving large systems of linear algebraic equations. In the approach under consideration, the solution of a system of algebraic linear equations is found as a point of intersection of hyperplanes, which needs a minimal amount of computer operating storage. Two examples are given. In the first example, the boundary value problem for a three-dimensional stationary heat transfer equation in a parallelepiped in ℝ3 is considered, where boundary value problems of first, second, or third order, or their combinations, are taken into account. The governing differential equations are reduced to algebraic ones with the help of the finite element and boundary element methods for different meshes applied. The obtained results are compared with known analytical solutions. The second example concerns computation of a nonhomogeneous shallow physically and geometrically nonlinear shell subject to transversal uniformly distributed load. The partial differential equations are reduced to a system of nonlinear algebraic equations with the error of O(hx12+hx22. The linearization process is realized through either Newton method or differentiation with respect to a parameter. In consequence, the relations of the boundary condition variations along the shell side and the conditions for the solution matching are reported.

Series: Utilization of Differential Equations and Methods for Solving Them in Medical Physics (4).

Science.gov (United States)

Murase, Kenya

2016-01-01

Partial differential equations are often used in the field of medical physics. In this (final) issue, the methods for solving the partial differential equations were introduced, which include separation of variables, integral transform (Fourier and Fourier-sine transforms), Green's function, and series expansion methods. Some examples were also introduced, in which the integral transform and Green's function methods were applied to solving Pennes' bioheat transfer equation and the Fourier series expansion method was applied to Navier-Stokes equation for analyzing the wall shear stress in blood vessels.Finally, the author hopes that this series will be helpful for people who engage in medical physics.

Khater method for nonlinear Sharma Tasso-Olever (STO) equation of fractional order

Science.gov (United States)

Bibi, Sadaf; Mohyud-Din, Syed Tauseef; Khan, Umar; Ahmed, Naveed

In this work, we have implemented a direct method, known as Khater method to establish exact solutions of nonlinear partial differential equations of fractional order. Number of solutions provided by this method is greater than other traditional methods. Exact solutions of nonlinear fractional order Sharma Tasso-Olever (STO) equation are expressed in terms of kink, travelling wave, periodic and solitary wave solutions. Modified Riemann-Liouville derivative and Fractional complex transform have been used for compatibility with fractional order sense. Solutions have been graphically simulated for understanding the physical aspects and importance of the method. A comparative discussion between our established results and the results obtained by existing ones is also presented. Our results clearly reveal that the proposed method is an effective, powerful and straightforward technique to work out new solutions of various types of differential equations of non-integer order in the fields of applied sciences and engineering.

Variational formulation and projectional methods for the second order transport equation

International Nuclear Information System (INIS)

Borysiewicz, M.; Stankiewicz, R.

1979-01-01

Herein the variational problem for a second-order boundary value problem for the neutron transport equation is formulated. The projectional methods solving the problem are examined. The approach is compared with that based on the original untransformed form of the neutron transport equation

The modified alternative (G'/G)-expansion method to nonlinear evolution equation: application to the (1+1)-dimensional Drinfel'd-Sokolov-Wilson equation.

Science.gov (United States)

Akbar, M Ali; Mohd Ali, Norhashidah Hj; Mohyud-Din, Syed Tauseef

2013-01-01

Over the years, (G'/G)-expansion method is employed to generate traveling wave solutions to various wave equations in mathematical physics. In the present paper, the alternative (G'/G)-expansion method has been further modified by introducing the generalized Riccati equation to construct new exact solutions. In order to illustrate the novelty and advantages of this approach, the (1+1)-dimensional Drinfel'd-Sokolov-Wilson (DSW) equation is considered and abundant new exact traveling wave solutions are obtained in a uniform way. These solutions may be imperative and significant for the explanation of some practical physical phenomena. It is shown that the modified alternative (G'/G)-expansion method an efficient and advance mathematical tool for solving nonlinear partial differential equations in mathematical physics.

Modulating functions method for parameters estimation in the fifth order KdV equation

KAUST Repository

Asiri, Sharefa M.; Liu, Da-Yan; Laleg-Kirati, Taous-Meriem

2017-01-01

In this work, the modulating functions method is proposed for estimating coefficients in higher-order nonlinear partial differential equation which is the fifth order Kortewegde Vries (KdV) equation. The proposed method transforms the problem into a

Solving Eigenvalue response matrix equations with Jacobian-Free Newton-Krylov methods

International Nuclear Information System (INIS)

Roberts, Jeremy A.; Forget, Benoit

2011-01-01

The response matrix method for reactor eigenvalue problems is motivated as a technique for solving coarse mesh transport equations, and the classical approach of power iteration (PI) for solution is described. The method is then reformulated as a nonlinear system of equations, and the associated Jacobian is derived. A Jacobian-Free Newton-Krylov (JFNK) method is employed to solve the system, using an approximate Jacobian coupled with incomplete factorization as a preconditioner. The unpreconditioned JFNK slightly outperforms PI, and preconditioned JFNK outperforms both PI and Steffensen-accelerated PI significantly. (author)

A constrained Hartree-Fock-Bogoliubov equation derived from the double variational method

International Nuclear Information System (INIS)

Onishi, Naoki; Horibata, Takatoshi.

1980-01-01

The double variational method is applied to the intrinsic state of the generalized BCS wave function. A constrained Hartree-Fock-Bogoliubov equation is derived explicitly in the form of an eigenvalue equation. A method of obtaining approximate overlap and energy overlap integrals is proposed. This will help development of numerical calculations of the angular momentum projection method, especially for general intrinsic wave functions without any symmetry restrictions. (author)

A transformed rational function method and exact solutions to the 3+1 dimensional Jimbo-Miwa equation

International Nuclear Information System (INIS)

Ma Wenxiu; Lee, J.-H.

2009-01-01

A direct approach to exact solutions of nonlinear partial differential equations is proposed, by using rational function transformations. The new method provides a more systematical and convenient handling of the solution process of nonlinear equations, unifying the tanh-function type methods, the homogeneous balance method, the exp-function method, the mapping method, and the F-expansion type methods. Its key point is to search for rational solutions to variable-coefficient ordinary differential equations transformed from given partial differential equations. As an application, the construction problem of exact solutions to the 3+1 dimensional Jimbo-Miwa equation is treated, together with a Baecklund transformation.

Imaginary Time Step Method to Solve the Dirac Equation with Nonlocal Potential

International Nuclear Information System (INIS)

Zhang Ying; Liang Haozhao; Meng Jie

2009-01-01

The imaginary time step (ITS) method is applied to solve the Dirac equation with nonlocal potentials in coordinate space. Taking the nucleus 12 C as an example, even with nonlocal potentials, the direct ITS evolution for the Dirac equation still meets the disaster of the Dirac sea. However, following the recipe in our former investigation, the disaster can be avoided by the ITS evolution for the corresponding Schroedinger-like equation without localization, which gives the convergent results exactly the same with those obtained iteratively by the shooting method with localized effective potentials.

New Efficient Fourth Order Method for Solving Nonlinear Equations

Directory of Open Access Journals (Sweden)

Farooq Ahmad

2013-12-01

Full Text Available In a paper [Appl. Math. Comput., 188 (2 (2007 1587--1591], authors have suggested and analyzed a method for solving nonlinear equations. In the present work, we modified this method by using the finite difference scheme, which has a quintic convergence. We have compared this modified Halley method with some other iterative of fifth-orders convergence methods, which shows that this new method having convergence of fourth order, is efficient.

Critical comparison between equation of motion-Green's function methods and configuration interaction methods: analysis of methods and applications

International Nuclear Information System (INIS)

Freed, K.F.; Herman, M.F.; Yeager, D.L.

1980-01-01

A description is provided of the common conceptual origins of many-body equations of motion and Green's function methods in Liouville operator formulations of the quantum mechanics of atomic and molecular electronic structure. Numerical evidence is provided to show the inadequacies of the traditional strictly perturbative approaches to these methods. Nonperturbative methods are introduced by analogy with techniques developed for handling large configuration interaction calculations and by evaluating individual matrix elements to higher accuracy. The important role of higher excitations is exhibited by the numerical calculations, and explicit comparisons are made between converged equations of motion and configuration interaction calculations for systems where a fundamental theorem requires the equality of the energy differences produced by these different approaches. (Auth.)

On a computer implementation of the block Gauss–Seidel method for normal systems of equations

Directory of Open Access Journals (Sweden)

Alexander I. Zhdanov

2016-12-01

Full Text Available This article focuses on the modification of the block option Gauss-Seidel method for normal systems of equations, which is a sufficiently effective method of solving generally overdetermined, systems of linear algebraic equations of high dimensionality. The main disadvantage of methods based on normal equations systems is the fact that the condition number of the normal system is equal to the square of the condition number of the original problem. This fact has a negative impact on the rate of convergence of iterative methods based on normal equations systems. To increase the speed of convergence of iterative methods based on normal equations systems, for solving ill-conditioned problems currently different preconditioners options are used that reduce the condition number of the original system of equations. However, universal preconditioner for all applications does not exist. One of the effective approaches that improve the speed of convergence of the iterative Gauss–Seidel method for normal systems of equations, is to use its version of the block. The disadvantage of the block Gauss–Seidel method for production systems is the fact that it is necessary to calculate the pseudoinverse matrix for each iteration. We know that finding the pseudoinverse is a difficult computational procedure. In this paper, we propose a procedure to replace the matrix pseudo-solutions to the problem of normal systems of equations by Cholesky. Normal equations arising at each iteration of Gauss–Seidel method, have a relatively low dimension compared to the original system. The results of numerical experimentation demonstrating the effectiveness of the proposed approach are given.

Coupled double-distribution-function lattice Boltzmann method for the compressible Navier-Stokes equations.

Science.gov (United States)

Li, Q; He, Y L; Wang, Y; Tao, W Q

2007-11-01

A coupled double-distribution-function lattice Boltzmann method is developed for the compressible Navier-Stokes equations. Different from existing thermal lattice Boltzmann methods, this method can recover the compressible Navier-Stokes equations with a flexible specific-heat ratio and Prandtl number. In the method, a density distribution function based on a multispeed lattice is used to recover the compressible continuity and momentum equations, while the compressible energy equation is recovered by an energy distribution function. The energy distribution function is then coupled to the density distribution function via the thermal equation of state. In order to obtain an adjustable specific-heat ratio, a constant related to the specific-heat ratio is introduced into the equilibrium energy distribution function. Two different coupled double-distribution-function lattice Boltzmann models are also proposed in the paper. Numerical simulations are performed for the Riemann problem, the double-Mach-reflection problem, and the Couette flow with a range of specific-heat ratios and Prandtl numbers. The numerical results are found to be in excellent agreement with analytical and/or other solutions.

The integral equation method applied to eddy currents

International Nuclear Information System (INIS)

Biddlecombe, C.S.; Collie, C.J.; Simkin, J.; Trowbridge, C.W.

1976-04-01

An algorithm for the numerical solution of eddy current problems is described, based on the direct solution of the integral equation for the potentials. In this method only the conducting and iron regions need to be divided into elements, and there are no boundary conditions. Results from two computer programs using this method for iron free problems for various two-dimensional geometries are presented and compared with analytic solutions. (author)

Test equating, scaling, and linking methods and practices

CERN Document Server

Kolen, Michael J

2014-01-01

This book provides an introduction to test equating, scaling, and linking, including those concepts and practical issues that are critical for developers and all other testing professionals.  In addition to statistical procedures, successful equating, scaling, and linking involves many aspects of testing, including procedures to develop tests, to administer and score tests, and to interpret scores earned on tests. Test equating methods are used with many standardized tests in education and psychology to ensure that scores from multiple test forms can be used interchangeably.  Test scaling is the process of developing score scales that are used when scores on standardized tests are reported. In test linking, scores from two or more tests are related to one another. Linking has received much recent attention, due largely to investigations of linking similarly named tests from different test publishers or tests constructed for different purposes. In recent years, researchers from the education, psychology, and...

Free-complement local-Schrödinger-equation method for solving the Schrödinger equation of atoms and molecules: Basic theories and features

Science.gov (United States)

Nakatsuji, Hiroshi; Nakashima, Hiroyuki

2015-02-01

The free-complement (FC) method is a general method for solving the Schrödinger equation (SE): The produced wave function has the potentially exact structure as the solution of the Schrödinger equation. The variables included are determined either by using the variational principle (FC-VP) or by imposing the local Schrödinger equations (FC-LSE) at the chosen set of the sampling points. The latter method, referred to as the local Schrödinger equation (LSE) method, is integral-free and therefore applicable to any atom and molecule. The purpose of this paper is to formulate the basic theories of the LSE method and explain their basic features. First, we formulate three variants of the LSE method, the AB, HS, and HTQ methods, and explain their properties. Then, the natures of the LSE methods are clarified in some detail using the simple examples of the hydrogen atom and the Hooke's atom. Finally, the ideas obtained in this study are applied to solving the SE of the helium atom highly accurately with the FC-LSE method. The results are very encouraging: we could get the world's most accurate energy of the helium atom within the sampling-type methodologies, which is comparable to those obtained with the FC-VP method. Thus, the FC-LSE method is an easy and yet a powerful integral-free method for solving the Schrödinger equation of general atoms and molecules.

Modeling Solution of Nonlinear Dispersive Partial Differential Equations using the Marker Method

International Nuclear Information System (INIS)

Lewandowski, Jerome L.V.

2005-01-01

A new method for the solution of nonlinear dispersive partial differential equations is described. The marker method relies on the definition of a convective field associated with the underlying partial differential equation; the information about the approximate solution is associated with the response of an ensemble of markers to this convective field. Some key aspects of the method, such as the selection of the shape function and the initial loading, are discussed in some details

A Table Lookup Method for Exact Analytical Solutions of Nonlinear Fractional Partial Differential Equations

Directory of Open Access Journals (Sweden)

Ji Juan-Juan

2017-01-01

Full Text Available A table lookup method for solving nonlinear fractional partial differential equations (fPDEs is proposed in this paper. Looking up the corresponding tables, we can quickly obtain the exact analytical solutions of fPDEs by using this method. To illustrate the validity of the method, we apply it to construct the exact analytical solutions of four nonlinear fPDEs, namely, the time fractional simplified MCH equation, the space-time fractional combined KdV-mKdV equation, the (2+1-dimensional time fractional Zoomeron equation, and the space-time fractional ZKBBM equation. As a result, many new types of exact analytical solutions are obtained including triangular periodic solution, hyperbolic function solution, singular solution, multiple solitary wave solution, and Jacobi elliptic function solution.

Numerical methods for solving the governing equations for a seriated continuum

International Nuclear Information System (INIS)

Narum, R.E.; Noble, C.; Mortensen, G.A.; McFadden, J.H.

1976-09-01

A desire to more accurately predict the behavior of transient two-phase flows has resulted in an investigation of the feasibility of computing unequal phase velocities and unequal phase temperatures. The finite difference forms of a set of equations governing a seriated continuum are presented along with two methods developed for solving the resulting systems of simultaneous nonlinear equations. Results from a one-dimensional computer code are presented to illustrate the capabilities of one of the solution methods

A Novel Method for Analytical Solutions of Fractional Partial Differential Equations

OpenAIRE

Mehmet Ali Akinlar; Muhammet Kurulay

2013-01-01

A new solution technique for analytical solutions of fractional partial differential equations (FPDEs) is presented. The solutions are expressed as a finite sum of a vector type functional. By employing MAPLE software, it is shown that the solutions might be extended to an arbitrary degree which makes the present method not only different from the others in the literature but also quite efficient. The method is applied to special Bagley-Torvik and Diethelm fractional differential equations as...

Perfectly Matched Layer for the Wave Equation Finite Difference Time Domain Method

Science.gov (United States)

Miyazaki, Yutaka; Tsuchiya, Takao

2012-07-01

The perfectly matched layer (PML) is introduced into the wave equation finite difference time domain (WE-FDTD) method. The WE-FDTD method is a finite difference method in which the wave equation is directly discretized on the basis of the central differences. The required memory of the WE-FDTD method is less than that of the standard FDTD method because no particle velocity is stored in the memory. In this study, the WE-FDTD method is first combined with the standard FDTD method. Then, Berenger's PML is combined with the WE-FDTD method. Some numerical demonstrations are given for the two- and three-dimensional sound fields.

Multigrid methods for partial differential equations - a short introduction

International Nuclear Information System (INIS)

Linden, J.; Stueben, K.

1993-01-01

These notes summarize the multigrid methods and emphasis is laid on the algorithmic concepts of multigrid for solving linear and non-linear partial differential equations. In this paper there is brief description of the basic structure of multigrid methods. Detailed introduction is also contained with applications to VLSI process simulation. (A.B.)

A study on linear and nonlinear Schrodinger equations by the variational iteration method

International Nuclear Information System (INIS)

Wazwaz, Abdul-Majid

2008-01-01

In this work, we introduce a framework to obtain exact solutions to linear and nonlinear Schrodinger equations. The He's variational iteration method (VIM) is used for analytic treatment of these equations. Numerical examples are tested to show the pertinent features of this method

Application of the bifurcation method to the modified Boussinesq equation

Directory of Open Access Journals (Sweden)

Shaoyong Li

2014-08-01

Firstly, we give a property of the solutions of the equation, that is, if $1+u(x, t$ is a solution, so is $1-u(x, t$. Secondly, by using the bifurcation method of dynamical systems we obtain some explicit expressions of solutions for the equation, which include kink-shaped solutions, blow-up solutions, periodic blow-up solutions and solitary wave solutions. Some previous results are extended.

The discontinuous finite element method for solving Eigenvalue problems of transport equations

International Nuclear Information System (INIS)

Yang, Shulin; Wang, Ruihong

2011-01-01

In this paper, the multigroup transport equations for solving the eigenvalues λ and K_e_f_f under two dimensional cylindrical coordinate are discussed. Aimed at the equations, the discretizing way combining discontinuous finite element method (DFE) with discrete ordinate method (SN) is developed, and the iterative algorithms and steps are studied. The numerical results show that the algorithms are efficient. (author)

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.

A polynomial expansion method and its application in the coupled Zakharov-Kuznetsov equations

International Nuclear Information System (INIS)

Huang Wenhua

2006-01-01

A polynomial expansion method is presented to solve nonlinear evolution equations. Applying this method, the coupled Zakharov-Kuznetsov equations in fluid system are studied and many exact travelling wave solutions are obtained. These solutions include solitary wave solutions, periodic solutions and rational type solutions

Modulating functions method for parameters estimation in the fifth order KdV equation

KAUST Repository

Asiri, Sharefa M.

2017-07-25

In this work, the modulating functions method is proposed for estimating coefficients in higher-order nonlinear partial differential equation which is the fifth order Kortewegde Vries (KdV) equation. The proposed method transforms the problem into a system of linear algebraic equations of the unknowns. The statistical properties of the modulating functions solution are described in this paper. In addition, guidelines for choosing the number of modulating functions, which is an important design parameter, are provided. The effectiveness and robustness of the proposed method are shown through numerical simulations in both noise-free and noisy cases.

A mixed finite element method for nonlinear diffusion equations

KAUST Repository

Burger, Martin; Carrillo, José ; Wolfram, Marie-Therese

2010-01-01

We propose a mixed finite element method for a class of nonlinear diffusion equations, which is based on their interpretation as gradient flows in optimal transportation metrics. We introduce an appropriate linearization of the optimal transport problem, which leads to a mixed symmetric formulation. This formulation preserves the maximum principle in case of the semi-discrete scheme as well as the fully discrete scheme for a certain class of problems. In addition solutions of the mixed formulation maintain exponential convergence in the relative entropy towards the steady state in case of a nonlinear Fokker-Planck equation with uniformly convex potential. We demonstrate the behavior of the proposed scheme with 2D simulations of the porous medium equations and blow-up questions in the Patlak-Keller-Segel model. © American Institute of Mathematical Sciences.

Comparative analysis among several methods used to solve the point kinetic equations

International Nuclear Information System (INIS)

Nunes, Anderson L.; Goncalves, Alessandro da C.; Martinez, Aquilino S.; Silva, Fernando Carvalho da

2007-01-01

The main objective of this work consists on the methodology development for comparison of several methods for the kinetics equations points solution. The evaluated methods are: the finite differences method, the stiffness confinement method, improved stiffness confinement method and the piecewise constant approximations method. These methods were implemented and compared through a systematic analysis that consists basically of confronting which one of the methods consume smaller computational time with higher precision. It was calculated the relative which function is to combine both criteria in order to reach the goal. Through the analyses of the performance factor it is possible to choose the best method for the solution of point kinetics equations. (author)

Comparative analysis among several methods used to solve the point kinetic equations

Energy Technology Data Exchange (ETDEWEB)

Nunes, Anderson L.; Goncalves, Alessandro da C.; Martinez, Aquilino S.; Silva, Fernando Carvalho da [Universidade Federal, Rio de Janeiro, RJ (Brazil). Coordenacao dos Programas de Pos-graduacao de Engenharia. Programa de Engenharia Nuclear; E-mails: alupo@if.ufrj.br; agoncalves@con.ufrj.br; aquilino@lmp.ufrj.br; fernando@con.ufrj.br

2007-07-01

The main objective of this work consists on the methodology development for comparison of several methods for the kinetics equations points solution. The evaluated methods are: the finite differences method, the stiffness confinement method, improved stiffness confinement method and the piecewise constant approximations method. These methods were implemented and compared through a systematic analysis that consists basically of confronting which one of the methods consume smaller computational time with higher precision. It was calculated the relative which function is to combine both criteria in order to reach the goal. Through the analyses of the performance factor it is possible to choose the best method for the solution of point kinetics equations. (author)

Partial differential equations with variable exponents variational methods and qualitative analysis

CERN Document Server

Radulescu, Vicentiu D

2015-01-01

Partial Differential Equations with Variable Exponents: Variational Methods and Qualitative Analysis provides researchers and graduate students with a thorough introduction to the theory of nonlinear partial differential equations (PDEs) with a variable exponent, particularly those of elliptic type. The book presents the most important variational methods for elliptic PDEs described by nonhomogeneous differential operators and containing one or more power-type nonlinearities with a variable exponent. The authors give a systematic treatment of the basic mathematical theory and constructive meth

Solution of the isotopic depletion equation using decomposition method and analytical solution

Energy Technology Data Exchange (ETDEWEB)

Prata, Fabiano S.; Silva, Fernando C.; Martinez, Aquilino S., E-mail: fprata@con.ufrj.br, E-mail: fernando@con.ufrj.br, E-mail: aquilino@lmp.ufrj.br [Coordenacao dos Programas de Pos-Graduacao de Engenharia (PEN/COPPE/UFRJ), RJ (Brazil). Programa de Engenharia Nuclear

2011-07-01

In this paper an analytical calculation of the isotopic depletion equations is proposed, featuring a chain of major isotopes found in a typical PWR reactor. Part of this chain allows feedback reactions of (n,2n) type. The method is based on decoupling the equations describing feedback from the rest of the chain by using the decomposition method, with analytical solutions for the other isotopes present in the chain. The method was implemented in a PWR reactor simulation code, that makes use of the nodal expansion method (NEM) to solve the neutron diffusion equation, describing the spatial distribution of neutron flux inside the reactor core. Because isotopic depletion calculation module is the most computationally intensive process within simulation systems of nuclear reactor core, it is justified to look for a method that is both efficient and fast, with the objective of evaluating a larger number of core configurations in a short amount of time. (author)

Solution of the isotopic depletion equation using decomposition method and analytical solution

International Nuclear Information System (INIS)

Prata, Fabiano S.; Silva, Fernando C.; Martinez, Aquilino S.

2011-01-01

In this paper an analytical calculation of the isotopic depletion equations is proposed, featuring a chain of major isotopes found in a typical PWR reactor. Part of this chain allows feedback reactions of (n,2n) type. The method is based on decoupling the equations describing feedback from the rest of the chain by using the decomposition method, with analytical solutions for the other isotopes present in the chain. The method was implemented in a PWR reactor simulation code, that makes use of the nodal expansion method (NEM) to solve the neutron diffusion equation, describing the spatial distribution of neutron flux inside the reactor core. Because isotopic depletion calculation module is the most computationally intensive process within simulation systems of nuclear reactor core, it is justified to look for a method that is both efficient and fast, with the objective of evaluating a larger number of core configurations in a short amount of time. (author)

Deterministic methods to solve the integral transport equation in neutronic

International Nuclear Information System (INIS)

Warin, X.

1993-11-01

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

Workshop on Recent Trends in Complex Methods for Partial Differential Equations

CERN Document Server

Celebi, A; Tutschke, Wolfgang

1999-01-01

This volume is a collection of manscripts mainly originating from talks and lectures given at the Workshop on Recent Trends in Complex Methods for Par­ tial Differential Equations held from July 6 to 10, 1998 at the Middle East Technical University in Ankara, Turkey, sponsored by The Scientific and Tech­ nical Research Council of Turkey and the Middle East Technical University. This workshop is a continuation oftwo workshops from 1988 and 1993 at the In­ ternational Centre for Theoretical Physics in Trieste, Italy entitled Functional analytic Methods in Complex Analysis and Applications to Partial Differential Equations. Since classical complex analysis of one and several variables has a long tra­ dition it is of high level. But most of its basic problems are solved nowadays so that within the last few decades it has lost more and more attention. The area of complex and functional analytic methods in partial differential equations, however, is still a growing and flourishing field, in particular as these ...

Solving equations by topological methods

Directory of Open Access Journals (Sweden)

Lech Górniewicz

2005-01-01

Full Text Available In this paper we survey most important results from topological fixed point theory which can be directly applied to differential equations. Some new formulations are presented. We believe that our article will be useful for analysts applying topological fixed point theory in nonlinear analysis and in differential equations.

Solution and study of nodal neutron transport equation applying the LTSN-DiagExp method

International Nuclear Information System (INIS)

Hauser, Eliete Biasotto; Pazos, Ruben Panta; Vilhena, Marco Tullio de; Barros, Ricardo Carvalho de

2003-01-01

In this paper we report advances about the three-dimensional nodal discrete-ordinates approximations of neutron transport equation for Cartesian geometry. We use the combined collocation method of the angular variables and nodal approach for the spatial variables. By nodal approach we mean the iterated transverse integration of the S N equations. This procedure leads to the set of one-dimensional averages angular fluxes in each spatial variable. The resulting system of equations is solved with the LTS N method, first applying the Laplace transform to the set of the nodal S N equations and then obtained the solution by symbolic computation. We include the LTS N method by diagonalization to solve the nodal neutron transport equation and then we outline the convergence of these nodal-LTS N approximations with the help of a norm associated to the quadrature formula used to approximate the integral term of the neutron transport equation. (author)

The application of He's exp-function method to a nonlinear differential-difference equation

International Nuclear Information System (INIS)

Dai Chaoqing; Cen Xu; Wu Shengsheng

2009-01-01

This paper applies He's exp-function method, which was originally proposed to find new exact travelling wave solutions of nonlinear partial differential equations (NPDEs) or coupled nonlinear partial differential equations (CNPDEs), to a nonlinear differential-difference equation, and some new travelling wave solutions are obtained.

Rothe's method for parabolic equations on non-cylindrical domains

Czech Academy of Sciences Publication Activity Database

Dasht, J.; Engström, J.; Kufner, Alois; Persson, L.E.

2006-01-01

Roč. 1, č. 1 (2006), s. 59-80 ISSN 0973-2306 Institutional research plan: CEZ:AV0Z10190503 Keywords : parabolic equations * non-cylindrical domains * Rothe's method * time-discretization Subject RIV: BA - General Mathematics

Exp-function method for constructing exact solutions of Sharma-Tasso-Olver equation

International Nuclear Information System (INIS)

Erbas, Baris; Yusufoglu, Elcin

2009-01-01

In this paper we use the Exp-function method for the analytic treatment of Sharma-Tasso-Olver equation. New solitonary solutions are formally derived. Change of parameters, which drastically changes the characteristics of the equations, is examined. It is shown that the Exp-function method provides a powerful mathematical tool for solving high-dimensional nonlinear evolutions in mathematical physics. The proposed schemes are reliable and manageable.

An Evaluation of Five Linear Equating Methods for the NEAT Design

Science.gov (United States)

Mroch, Andrew A.; Suh, Youngsuk; Kane, Michael T.; Ripkey, Douglas R.

2009-01-01

This study uses the results of two previous papers (Kane, Mroch, Suh, & Ripkey, this issue; Suh, Mroch, Kane, & Ripkey, this issue) and the literature on linear equating to evaluate five linear equating methods along several dimensions, including the plausibility of their assumptions and their levels of bias and root mean squared difference…

Adaptive Finite Volume Method for the Shallow Water Equations on Triangular Grids

Directory of Open Access Journals (Sweden)

Sudi Mungkasi

2016-01-01

Full Text Available This paper presents a numerical entropy production (NEP scheme for two-dimensional shallow water equations on unstructured triangular grids. We implement NEP as the error indicator for adaptive mesh refinement or coarsening in solving the shallow water equations using a finite volume method. Numerical simulations show that NEP is successful to be a refinement/coarsening indicator in the adaptive mesh finite volume method, as the method refines the mesh or grids around nonsmooth regions and coarsens them around smooth regions.

A three operator split-step method covering a larger set of non-linear partial differential equations

Science.gov (United States)

Zia, Haider

2017-06-01

This paper describes an updated exponential Fourier based split-step method that can be applied to a greater class of partial differential equations than previous methods would allow. These equations arise in physics and engineering, a notable example being the generalized derivative non-linear Schrödinger equation that arises in non-linear optics with self-steepening terms. These differential equations feature terms that were previously inaccessible to model accurately with low computational resources. The new method maintains a 3rd order error even with these additional terms and models the equation in all three spatial dimensions and time. The class of non-linear differential equations that this method applies to is shown. The method is fully derived and implementation of the method in the split-step architecture is shown. This paper lays the mathematical ground work for an upcoming paper employing this method in white-light generation simulations in bulk material.

Numerical method for solving linear Fredholm fuzzy integral equations of the second kind

Energy Technology Data Exchange (ETDEWEB)

Abbasbandy, S. [Department of Mathematics, Imam Khomeini International University, P.O. Box 288, Ghazvin 34194 (Iran, Islamic Republic of)]. E-mail: saeid@abbasbandy.com; Babolian, E. [Faculty of Mathematical Sciences and Computer Engineering, Teacher Training University, Tehran 15618 (Iran, Islamic Republic of); Alavi, M. [Department of Mathematics, Arak Branch, Islamic Azad University, Arak 38135 (Iran, Islamic Republic of)

2007-01-15

In this paper we use parametric form of fuzzy number and convert a linear fuzzy Fredholm integral equation to two linear system of integral equation of the second kind in crisp case. We can use one of the numerical method such as Nystrom and find the approximation solution of the system and hence obtain an approximation for fuzzy solution of the linear fuzzy Fredholm integral equations of the second kind. The proposed method is illustrated by solving some numerical examples.

The Closest Point Method and Multigrid Solvers for Elliptic Equations on Surfaces

KAUST Repository

Chen, Yujia; Macdonald, Colin B.

2015-01-01

© 2015 Society for Industrial and Applied Mathematics. Elliptic partial differential equations are important from both application and analysis points of view. In this paper we apply the closest point method to solve elliptic equations on general

Stable multi-domain spectral penalty methods for fractional partial differential equations

Science.gov (United States)

Xu, Qinwu; Hesthaven, Jan S.

2014-01-01

We propose stable multi-domain spectral penalty methods suitable for solving fractional partial differential equations with fractional derivatives of any order. First, a high order discretization is proposed to approximate fractional derivatives of any order on any given grids based on orthogonal polynomials. The approximation order is analyzed and verified through numerical examples. Based on the discrete fractional derivative, we introduce stable multi-domain spectral penalty methods for solving fractional advection and diffusion equations. The equations are discretized in each sub-domain separately and the global schemes are obtained by weakly imposed boundary and interface conditions through a penalty term. Stability of the schemes are analyzed and numerical examples based on both uniform and nonuniform grids are considered to highlight the flexibility and high accuracy of the proposed schemes.

Solution method for the unsteady incompressible Navier-Stokes equations in generalized coordinate systems

International Nuclear Information System (INIS)

Rosenfeld, M.; Kwak, D.; Vinokur, M.

1988-01-01

A solution method based on a fractional step approach is developed for obtaining time-dependent solutions of the three-dimensional, incompressible Navier-Stokes equations in generalized coordinate systems. The governing equations are discretized conservatively by finite volumes using a staggered mesh system. The primitive variable formulation uses the volume fluxes across the faces of each computational cell as dependent variables. This procedure, combined with accurate and consistent approximations of geometric parameters, is done to satisfy the discretized mass conservation equation to machine accuracy as well as to gain favorable convergence properties of the Poisson solver. The discretized equations are second-order-accurate in time and space and no smoothing terms are added. An approximate-factorization scheme is implemented in solving the momentum equations. A novel ZEBRA scheme with four-color ordering is devised for the efficient solution of the Poisson equation. Several two and three-dimensional solutions are compared with other numerical and experimental results to validate the present method. 23 references

New Iterative Method for Fractional Gas Dynamics and Coupled Burger’s Equations

Directory of Open Access Journals (Sweden)

Mohamed S. Al-luhaibi

2015-01-01

Full Text Available This paper presents the approximate analytical solutions to solve the nonlinear gas dynamics and coupled Burger’s equations with fractional time derivative. By using initial values, the explicit solutions of the equations are solved by using a reliable algorithm. Numerical results show that the new iterative method is easy to implement and accurate when applied to time-fractional partial differential equations.

Finding all solutions of nonlinear equations using the dual simplex method

Science.gov (United States)

Yamamura, Kiyotaka; Fujioka, Tsuyoshi

2003-03-01

Recently, an efficient algorithm has been proposed for finding all solutions of systems of nonlinear equations using linear programming. This algorithm is based on a simple test (termed the LP test) for nonexistence of a solution to a system of nonlinear equations using the dual simplex method. In this letter, an improved version of the LP test algorithm is proposed. By numerical examples, it is shown that the proposed algorithm could find all solutions of a system of 300 nonlinear equations in practical computation time.

Numerical simulation for fractional order stationary neutron transport equation using Haar wavelet collocation method

Energy Technology Data Exchange (ETDEWEB)

Saha Ray, S., E-mail: santanusaharay@yahoo.com; Patra, A.

2014-10-15

Highlights: • A stationary transport equation has been solved using the technique of Haar wavelet collocation method. • This paper intends to provide the great utility of Haar wavelets to nuclear science problem. • In the present paper, two-dimensional Haar wavelets are applied. • The proposed method is mathematically very simple, easy and fast. - Abstract: In this paper the numerical solution for the fractional order stationary neutron transport equation is presented using Haar wavelet Collocation Method (HWCM). Haar wavelet collocation method is efficient and powerful in solving wide class of linear and nonlinear differential equations. This paper intends to provide an application of Haar wavelets to nuclear science problems. This paper describes the application of Haar wavelets for the numerical solution of fractional order stationary neutron transport equation in homogeneous medium with isotropic scattering. The proposed method is mathematically very simple, easy and fast. To demonstrate about the efficiency and applicability of the method, two test problems are discussed.

Implementation of compact finite-difference method to parabolized Navier-Stokes equations

International Nuclear Information System (INIS)

Esfahanian, V.; Hejranfar, K.; Darian, H.M.

2005-01-01

The numerical simulation of the Parabolized Navier-Stokes (PNS) equations for supersonic/hypersonic flow field is obtained by using the fourth-order compact finite-difference method. The PNS equations in the general curvilinear coordinates are solved by using the implicit finite-difference algorithm of Beam and Warming. A shock fitting procedure is utilized to obtain the accurate solution in the vicinity of the shock. The computations are performed for hypersonic axisymmetric flow over a blunt cone. The present results for the flow field along with those of the second-order method are presented and accuracy analysis is performed to insure the fourth-order accuracy of the method. (author)

Development of a set of benchmark problems to verify numerical methods for solving burnup equations

International Nuclear Information System (INIS)

Lago, Daniel; Rahnema, Farzad

2017-01-01

Highlights: • Description transmutation chain benchmark problems. • Problems for validating numerical methods for solving burnup equations. • Analytical solutions for the burnup equations. • Numerical solutions for the burnup equations. - Abstract: A comprehensive set of transmutation chain benchmark problems for numerically validating methods for solving burnup equations was created. These benchmark problems were designed to challenge both traditional and modern numerical methods used to solve the complex set of ordinary differential equations used for tracking the change in nuclide concentrations over time due to nuclear phenomena. Given the development of most burnup solvers is done for the purpose of coupling with an established transport solution method, these problems provide a useful resource in testing and validating the burnup equation solver before coupling for use in a lattice or core depletion code. All the relevant parameters for each benchmark problem are described. Results are also provided in the form of reference solutions generated by the Mathematica tool, as well as additional numerical results from MATLAB.

Numerical simulation of GEW equation using RBF collocation method

Directory of Open Access Journals (Sweden)

Hamid Panahipour

2012-08-01

Full Text Available The generalized equal width (GEW equation is solved numerically by a meshless method based on a global collocation with standard types of radial basis functions (RBFs. Test problems including propagation of single solitons, interaction of two and three solitons, development of the Maxwellian initial condition pulses, wave undulation and wave generation are used to indicate the efficiency and accuracy of the method. Comparisons are made between the results of the proposed method and some other published numerical methods.

Isotope decay equations solved by means of a recursive method

International Nuclear Information System (INIS)

Grant, Carlos

2009-01-01

The isotope decay equations have been solved using forward finite differences taking small time steps, among other methods. This is the case of the cell code WIMS, where it is assumed that concentrations of all fissionable isotopes remain constant during the integration interval among other simplifications. Even when the problem could be solved running through a logical tree, all algorithms used for resolution of these equations used an iterative programming formulation. That happened because nearly all computer languages used up to a recent past by the scientific programmers did not support recursion, such as the case of the old versions of FORTRAN or BASIC. Nowadays also an integral form of the depletion equations is used in Monte Carlo simulation. In this paper we propose another programming solution using a recursive algorithm, running through all descendants of each isotope and adding their contributions to all isotopes in each generation. The only assumption made for this solution is that fluxes remain constant during the whole time step. Recursive process is interrupted when a stable isotope was attained or the calculated contributions are smaller than a given precision. These algorithms can be solved by means an exact analytic method that can have some problems when circular loops appear for isotopes with alpha decay, and a more general polynomial method. Both methods are shown. (author)

Exact Solutions for Fractional Differential-Difference Equations by an Extended Riccati Sub-ODE Method

International Nuclear Information System (INIS)

Feng Qinghua

2013-01-01

In this paper, an extended Riccati sub-ODE method is proposed to establish new exact solutions for fractional differential-difference equations in the sense of modified Riemann—Liouville derivative. By a fractional complex transformation, a given fractional differential-difference equation can be turned into another differential-difference equation of integer order. The validity of the method is illustrated by applying it to solve the fractional Hybrid lattice equation and the fractional relativistic Toda lattice system. As a result, some new exact solutions including hyperbolic function solutions, trigonometric function solutions and rational solutions are established. (general)

A Globally Convergent Matrix-Free Method for Constrained Equations and Its Linear Convergence Rate

Directory of Open Access Journals (Sweden)

Min Sun

2014-01-01

Full Text Available A matrix-free method for constrained equations is proposed, which is a combination of the well-known PRP (Polak-Ribière-Polyak conjugate gradient method and the famous hyperplane projection method. The new method is not only derivative-free, but also completely matrix-free, and consequently, it can be applied to solve large-scale constrained equations. We obtain global convergence of the new method without any differentiability requirement on the constrained equations. Compared with the existing gradient methods for solving such problem, the new method possesses linear convergence rate under standard conditions, and a relax factor γ is attached in the update step to accelerate convergence. Preliminary numerical results show that it is promising in practice.

Exact Solutions of Fractional Burgers and Cahn-Hilliard Equations Using Extended Fractional Riccati Expansion Method

Directory of Open Access Journals (Sweden)

Wei Li

2014-01-01

Full Text Available Based on a general fractional Riccati equation and with Jumarie’s modified Riemann-Liouville derivative to an extended fractional Riccati expansion method for solving the time fractional Burgers equation and the space-time fractional Cahn-Hilliard equation, the exact solutions expressed by the hyperbolic functions and trigonometric functions are obtained. The obtained results show that the presented method is effective and appropriate for solving nonlinear fractional differential equations.

Exact Solutions of the Space Time Fractional Symmetric Regularized Long Wave Equation Using Different Methods

Directory of Open Access Journals (Sweden)

Özkan Güner

2014-01-01

Full Text Available We apply the functional variable method, exp-function method, and (G′/G-expansion method to establish the exact solutions of the nonlinear fractional partial differential equation (NLFPDE in the sense of the modified Riemann-Liouville derivative. As a result, some new exact solutions for them are obtained. The results show that these methods are very effective and powerful mathematical tools for solving nonlinear fractional equations arising in mathematical physics. As a result, these methods can also be applied to other nonlinear fractional differential equations.

A negative-norm least-squares method for time-harmonic Maxwell equations

KAUST Repository

Copeland, Dylan M.

2012-04-01

This paper presents and analyzes a negative-norm least-squares finite element discretization method for the dimension-reduced time-harmonic Maxwell equations in the case of axial symmetry. The reduced equations are expressed in cylindrical coordinates, and the analysis consequently involves weighted Sobolev spaces based on the degenerate radial weighting. The main theoretical results established in this work include existence and uniqueness of the continuous and discrete formulations and error estimates for simple finite element functions. Numerical experiments confirm the error estimates and efficiency of the method for piecewise constant coefficients. © 2011 Elsevier Inc.

Lagrange–Galerkin methods for the incompressible Navier-Stokes equations: a review

Directory of Open Access Journals (Sweden)

Bermejo Rodolfo

2016-09-01

Full Text Available We review in this paper the development of Lagrange-Galerkin (LG methods to integrate the incompressible Navier-Stokes equations (NSEs for engineering applications. These methods were introduced in the computational fluid dynamics community in the early eighties of the past century, and at that time they were considered good methods for both their theoretical stability properties and the way of dealing with the nonlinear terms of the equations; however, the numerical experience gained with the application of LG methods to different problems has identified drawbacks of them, such as the calculation of specific integrals that arise in their formulation and the calculation of the ow trajectories, which somehow have hampered the applicability of LG methods. In this paper, we focus on these issues and summarize the convergence results of LG methods; furthermore, we shall briefly introduce a new stabilized LG method suitable for high Reynolds numbers.

Pentadiagonal alternating-direction-implicit finite-difference time-domain method for two-dimensional Schrödinger equation

Science.gov (United States)

Tay, Wei Choon; Tan, Eng Leong

2014-07-01

In this paper, we have proposed a pentadiagonal alternating-direction-implicit (Penta-ADI) finite-difference time-domain (FDTD) method for the two-dimensional Schrödinger equation. Through the separation of complex wave function into real and imaginary parts, a pentadiagonal system of equations for the ADI method is obtained, which results in our Penta-ADI method. The Penta-ADI method is further simplified into pentadiagonal fundamental ADI (Penta-FADI) method, which has matrix-operator-free right-hand-sides (RHS), leading to the simplest and most concise update equations. As the Penta-FADI method involves five stencils in the left-hand-sides (LHS) of the pentadiagonal update equations, special treatments that are required for the implementation of the Dirichlet's boundary conditions will be discussed. Using the Penta-FADI method, a significantly higher efficiency gain can be achieved over the conventional Tri-ADI method, which involves a tridiagonal system of equations.

The Method of Lines Solution of the Regularized Long-Wave Equation Using Runge-Kutta Time Discretization Method

Directory of Open Access Journals (Sweden)

H. O. Bakodah

2013-01-01

Full Text Available A method of lines approach to the numerical solution of nonlinear wave equations typified by the regularized long wave (RLW is presented. The method developed uses a finite differences discretization to the space. Solution of the resulting system was obtained by applying fourth Runge-Kutta time discretization method. Using Von Neumann stability analysis, it is shown that the proposed method is marginally stable. To test the accuracy of the method some numerical experiments on test problems are presented. Test problems including solitary wave motion, two-solitary wave interaction, and the temporal evaluation of a Maxwellian initial pulse are studied. The accuracy of the present method is tested with and error norms and the conservation properties of mass, energy, and momentum under the RLW equation.

Semigroup methods for evolution equations on networks

CERN Document Server

Mugnolo, Delio

2014-01-01

This concise text is based on a series of lectures held only a few years ago and originally intended as an introduction to known results on linear hyperbolic and parabolic equations.  Yet the topic of differential equations on graphs, ramified spaces, and more general network-like objects has recently gained significant momentum and, well beyond the confines of mathematics, there is a lively interdisciplinary discourse on all aspects of so-called complex networks. Such network-like structures can be found in virtually all branches of science, engineering and the humanities, and future research thus calls for solid theoretical foundations.      This book is specifically devoted to the study of evolution equations – i.e., of time-dependent differential equations such as the heat equation, the wave equation, or the Schrödinger equation (quantum graphs) – bearing in mind that the majority of the literature in the last ten years on the subject of differential equations of graphs has been devoted to ellip...

Integral equation models for image restoration: high accuracy methods and fast algorithms

International Nuclear Information System (INIS)

Lu, Yao; Shen, Lixin; Xu, Yuesheng

2010-01-01

Discrete models are consistently used as practical models for image restoration. They are piecewise constant approximations of true physical (continuous) models, and hence, inevitably impose bottleneck model errors. We propose to work directly with continuous models for image restoration aiming at suppressing the model errors caused by the discrete models. A systematic study is conducted in this paper for the continuous out-of-focus image models which can be formulated as an integral equation of the first kind. The resulting integral equation is regularized by the Lavrentiev method and the Tikhonov method. We develop fast multiscale algorithms having high accuracy to solve the regularized integral equations of the second kind. Numerical experiments show that the methods based on the continuous model perform much better than those based on discrete models, in terms of PSNR values and visual quality of the reconstructed images

Analytical Investigation of Beam Deformation Equation using Perturbation, Homotopy Perturbation, Variational Iteration and Optimal Homotopy Asymptotic Methods

DEFF Research Database (Denmark)

Farrokhzad, F.; Mowlaee, P.; Barari, Amin

2011-01-01

The beam deformation equation has very wide applications in structural engineering. As a differential equation, it has its own problem concerning existence, uniqueness and methods of solutions. Often, original forms of governing differential equations used in engineering problems are simplified...... Method (OHAM). The comparisons of the results reveal that these methods are very effective, convenient and quite accurate to systems of non-linear differential equation......., and this process produces noise in the obtained answers. This paper deals with solution of second order of differential equation governing beam deformation using four analytical approximate methods, namely the Homotopy Perturbation Method (HPM), Variational Iteration Method (VIM) and Optimal Homotopy Asymptotic...

Shifted Legendre method with residual error estimation for delay linear Fredholm integro-differential equations

Directory of Open Access Journals (Sweden)

Şuayip Yüzbaşı

2017-03-01

Full Text Available In this paper, we suggest a matrix method for obtaining the approximate solutions of the delay linear Fredholm integro-differential equations with constant coefficients using the shifted Legendre polynomials. The problem is considered with mixed conditions. Using the required matrix operations, the delay linear Fredholm integro-differential equation is transformed into a matrix equation. Additionally, error analysis for the method is presented using the residual function. Illustrative examples are given to demonstrate the efficiency of the method. The results obtained in this study are compared with the known results.

Three-Step Predictor-Corrector of Exponential Fitting Method for Nonlinear Schroedinger Equations

International Nuclear Information System (INIS)

Tang Chen; Zhang Fang; Yan Haiqing; Luo Tao; Chen Zhanqing

2005-01-01

We develop the three-step explicit and implicit schemes of exponential fitting methods. We use the three-step explicit exponential fitting scheme to predict an approximation, then use the three-step implicit exponential fitting scheme to correct this prediction. This combination is called the three-step predictor-corrector of exponential fitting method. The three-step predictor-corrector of exponential fitting method is applied to numerically compute the coupled nonlinear Schroedinger equation and the nonlinear Schroedinger equation with varying coefficients. The numerical results show that the scheme is highly accurate.

Using a Linear Regression Method to Detect Outliers in IRT Common Item Equating

Science.gov (United States)

He, Yong; Cui, Zhongmin; Fang, Yu; Chen, Hanwei

2013-01-01

Common test items play an important role in equating alternate test forms under the common item nonequivalent groups design. When the item response theory (IRT) method is applied in equating, inconsistent item parameter estimates among common items can lead to large bias in equated scores. It is prudent to evaluate inconsistency in parameter…

Multiple travelling wave solutions of nonlinear evolution equations using a unified algebraic method

International Nuclear Information System (INIS)

Fan Engui

2002-01-01

A new direct and unified algebraic method for constructing multiple travelling wave solutions of general nonlinear evolution equations is presented and implemented in a computer algebraic system. Compared with most of the existing tanh methods, the Jacobi elliptic function method or other sophisticated methods, the proposed method not only gives new and more general solutions, but also provides a guideline to classify the various types of the travelling wave solutions according to the values of some parameters. The solutions obtained in this paper include (a) kink-shaped and bell-shaped soliton solutions, (b) rational solutions, (c) triangular periodic solutions and (d) Jacobi and Weierstrass doubly periodic wave solutions. Among them, the Jacobi elliptic periodic wave solutions exactly degenerate to the soliton solutions at a certain limit condition. The efficiency of the method can be demonstrated on a large variety of nonlinear evolution equations such as those considered in this paper, KdV-MKdV, Ito's fifth MKdV, Hirota, Nizhnik-Novikov-Veselov, Broer-Kaup, generalized coupled Hirota-Satsuma, coupled Schroedinger-KdV, (2+1)-dimensional dispersive long wave, (2+1)-dimensional Davey-Stewartson equations. In addition, as an illustrative sample, the properties of the soliton solutions and Jacobi doubly periodic solutions for the Hirota equation are shown by some figures. The links among our proposed method, the tanh method, extended tanh method and the Jacobi elliptic function method are clarified generally. (author)

NUMERICAL METHODS FOR SOLVING THE MULTI-TERM TIME-FRACTIONAL WAVE-DIFFUSION EQUATION

OpenAIRE

Liu, F.; Meerschaert, M.M.; McGough, R.J.; Zhuang, P.; Liu, Q.

2013-01-01

In this paper, the multi-term time-fractional wave-diffusion equations are considered. The multi-term time fractional derivatives are defined in the Caputo sense, whose orders belong to the intervals [0,1], [1,2), [0,2), [0,3), [2,3) and [2,4), respectively. Some computationally effective numerical methods are proposed for simulating the multi-term time-fractional wave-diffusion equations. The numerical results demonstrate the effectiveness of theoretical analysis. These methods and technique...

New method for solving three-dimensional Schroedinger equation

International Nuclear Information System (INIS)

Melezhik, V.S.

1990-01-01

The method derived recently for solving a multidimensional scattering problem is applied to a three-dimensional Schroedinger equation. As compared with direct three-dimensional calculations of finite elements and finite differences, this approach gives sufficiently accurate upper and lower approximations to the helium-atom binding energy, which demonstrates its efficiency. 15 refs.; 1 fig.; 2 tabs

A relaxation-projection method for compressible flows. Part I: The numerical equation of state for the Euler equations

International Nuclear Information System (INIS)

Saurel, Richard; Franquet, Erwin; Daniel, Eric; Le Metayer, Olivier

2007-01-01

A new projection method is developed for the Euler equations to determine the thermodynamic state in computational cells. It consists in the resolution of a mechanical relaxation problem between the various sub-volumes present in a computational cell. These sub-volumes correspond to the ones traveled by the various waves that produce states with different pressures, velocities, densities and temperatures. Contrarily to Godunov type schemes the relaxed state corresponds to mechanical equilibrium only and remains out of thermal equilibrium. The pressure computation with this relaxation process replaces the use of the conventional equation of state (EOS). A simplified relaxation method is also derived and provides a specific EOS (named the Numerical EOS). The use of the Numerical EOS gives a cure to spurious pressure oscillations that appear at contact discontinuities for fluids governed by real gas EOS. It is then extended to the computation of interface problems separating fluids with different EOS (liquid-gas interface for example) with the Euler equations. The resulting method is very robust, accurate, oscillation free and conservative. For the sake of simplicity and efficiency the method is developed in a Lagrange-projection context and is validated over exact solutions. In a companion paper [F. Petitpas, E. Franquet, R. Saurel, A relaxation-projection method for compressible flows. Part II: computation of interfaces and multiphase mixtures with stiff mechanical relaxation. J. Comput. Phys. (submitted for publication)], the method is extended to the numerical approximation of a non-conservative hyperbolic multiphase flow model for interface computation and shock propagation into mixtures

Numerical Simulation of Antennas with Improved Integral Equation Method

International Nuclear Information System (INIS)

Ma Ji; Fang Guang-You; Lu Wei

2015-01-01

Simulating antennas around a conducting object is a challenge task in computational electromagnetism, which is concerned with the behaviour of electromagnetic fields. To analyze this model efficiently, an improved integral equation-fast Fourier transform (IE-FFT) algorithm is presented in this paper. The proposed scheme employs two Cartesian grids with different size and location to enclose the antenna and the other object, respectively. On the one hand, IE-FFT technique is used to store matrix in a sparse form and accelerate the matrix-vector multiplication for each sub-domain independently. On the other hand, the mutual interaction between sub-domains is taken as the additional exciting voltage in each matrix equation. By updating integral equations several times, the whole electromagnetic system can achieve a stable status. Finally, the validity of the presented method is verified through the analysis of typical antennas in the presence of a conducting object. (paper)

Fully Digital Chaotic Differential Equation-based Systems And Methods

KAUST Repository

Radwan, Ahmed Gomaa Ahmed; Zidan, Mohammed A.; Salama, Khaled N.

2012-01-01

Various embodiments are provided for fully digital chaotic differential equation-based systems and methods. In one embodiment, among others, a digital circuit includes digital state registers and one or more digital logic modules configured to obtain a first value from two or more of the digital state registers; determine a second value based upon the obtained first values and a chaotic differential equation; and provide the second value to set a state of one of the plurality of digital state registers. In another embodiment, a digital circuit includes digital state registers, digital logic modules configured to obtain outputs from a subset of the digital shift registers and to provide the input based upon a chaotic differential equation for setting a state of at least one of the subset of digital shift registers, and a digital clock configured to provide a clock signal for operating the digital shift registers.

Fully Digital Chaotic Differential Equation-based Systems And Methods

KAUST Repository

Radwan, Ahmed Gomaa Ahmed

2012-09-06

Various embodiments are provided for fully digital chaotic differential equation-based systems and methods. In one embodiment, among others, a digital circuit includes digital state registers and one or more digital logic modules configured to obtain a first value from two or more of the digital state registers; determine a second value based upon the obtained first values and a chaotic differential equation; and provide the second value to set a state of one of the plurality of digital state registers. In another embodiment, a digital circuit includes digital state registers, digital logic modules configured to obtain outputs from a subset of the digital shift registers and to provide the input based upon a chaotic differential equation for setting a state of at least one of the subset of digital shift registers, and a digital clock configured to provide a clock signal for operating the digital shift registers.

Exact Solution of Space-Time Fractional Coupled EW and Coupled MEW Equations Using Modified Kudryashov Method

International Nuclear Information System (INIS)

Raslan, K. R.; Ali, Khalid K.; EL-Danaf, Talaat S.

2017-01-01

In the present paper, we established a traveling wave solution by using modified Kudryashov method for the space-time fractional nonlinear partial differential equations. The method is used to obtain the exact solutions for different types of the space-time fractional nonlinear partial differential equations such as, the space-time fractional coupled equal width wave equation (CEWE) and the space-time fractional coupled modified equal width wave equation (CMEW), which are the important soliton equations. Both equations are reduced to ordinary differential equations by the use of fractional complex transform and properties of modified Riemann–Liouville derivative. We plot the exact solutions for these equations at different time levels. (paper)

Approximate Analytic and Numerical Solutions to Lane-Emden Equation via Fuzzy Modeling Method

Directory of Open Access Journals (Sweden)

De-Gang Wang

2012-01-01

Full Text Available A novel algorithm, called variable weight fuzzy marginal linearization (VWFML method, is proposed. This method can supply approximate analytic and numerical solutions to Lane-Emden equations. And it is easy to be implemented and extended for solving other nonlinear differential equations. Numerical examples are included to demonstrate the validity and applicability of the developed technique.

Fourier method for three-dimensional partial differential equations in periodic geometry. Application: HELIAC

International Nuclear Information System (INIS)

Shestakov, A.I.; Mirin, A.A.

1984-01-01

A numerical method based on Fourier expansions and finite differences is presented. The method is demonstrated by solving a scalar, three-dimensional elliptic equation arising in MFE research, but has applicability to a wider class of problems. The scheme solves equations whose solutions are expected to be periodic in one or more of the independent variables

On the method of solution of the differential-delay Toda equation

Science.gov (United States)

Villarroel, Javier; Ablowitz, Mark J.

1993-09-01

The method of solution of the Toda differential-delay equation, which is a reduction of the Toda equation in 2+1 dimensions, is described. An important feature of the solution process is to obtain and study a novel Riemann-Hilbert problem. The latter problem requires factorization across an infinite number of strips with a suitable branching structure. Explicit soliton solutions are given.

A Special Variant of the Moment Method for Fredholm Integral Equations of the Second Kind

Directory of Open Access Journals (Sweden)

S. A. Solov’eva

2015-01-01

Full Text Available We consider the linear Fredholm integral equation of the second kind, where the kernel and the free term are smooth functions. We find the unknown function in this class as well.Exact and approximate methods for the solution of linear Fredholm integral equations of the second kind are well developed. However, classical methods do not take into account the structural properties of the kernel and the free term of equation.In this paper we develop and justify a special variant of the moment method to solve this equation, which takes into account the differential properties of initial data. The proposed paper furthers studies of N.S Gabbasov, I.P. Kasakina, and S.A Solov’eva. We use approximation theory, version of the general theory of approximate methods of analysis that Gabdulkhayev B.G suggested, and methods of functional analysis to prove theorems. In addition, we use N.S. Gabbasov’s ideas and methods in papers that are devoted to the Fredholm equations of the first kind, as well as N.S. Gabbasov and S.A Solov’eva’s investigations on the Fredholm equations of the third kind in the space of distributions.The first part of the paper provides a description of the basic function space and elements of the theory of approximation in it.In the second part we propose and theoretically justify a generalized moment method. We have demonstrated that the improvement of differential properties of the initial data improves the approximation accuracy. Since, in practice, the approximate equations are solved, as a rule, only approximately, we prove the stability and causality of the proposed method. The resulting estimate of the paper is in good agreement with the estimate for the ordinary moment method for equations of the second kind in the space of continuous functions.In the final section we have shown that a developed method is optimal in order of accuracy among all polynomial projection methods to solve Fredholm integral equations of the second

An element-free Galerkin (EFG) method for generalized Fisher equations (GFE)

International Nuclear Information System (INIS)

Shi Ting-Yu; Ge Hong-Xia; Cheng Rong-Jun

2013-01-01

A generalized Fisher equation (GFE) relates the time derivative of the average of the intrinsic rate of growth to its variance. The exact mathematical result of the GFE has been widely used in population dynamics and genetics, where it originated. Many researchers have studied the numerical solutions of the GFE, up to now. In this paper, we introduce an element-free Galerkin (EFG) method based on the moving least-square approximation to approximate positive solutions of the GFE from population dynamics. Compared with other numerical methods, the EFG method for the GFE needs only scattered nodes instead of meshing the domain of the problem. The Galerkin weak form is used to obtain the discrete equations, and the essential boundary conditions are enforced by the penalty method. In comparison with the traditional method, numerical solutions show that the new method has higher accuracy and better convergence. Several numerical examples are presented to demonstrate the effectiveness of the method

Implicit methods for the Navier-Stokes equations

Science.gov (United States)

Yoon, S.; Kwak, D.

1990-01-01

Numerical solutions of the Navier-Stokes equations using explicit schemes can be obtained at the expense of efficiency. Conventional implicit methods which often achieve fast convergence rates suffer high cost per iteration. A new implicit scheme based on lower-upper factorization and symmetric Gauss-Seidel relaxation offers very low cost per iteration as well as fast convergence. High efficiency is achieved by accomplishing the complete vectorizability of the algorithm on oblique planes of sweep in three dimensions.

Some New Lie Symmetry Groups of Differential-Difference Equations Obtained from a Simple Direct Method

International Nuclear Information System (INIS)

Zhi Hongyan

2009-01-01

In this paper, based on the symbolic computing system Maple, the direct method for Lie symmetry groups presented by Sen-Yue Lou [J. Phys. A: Math. Gen. 38 (2005) L129] is extended from the continuous differential equations to the differential-difference equations. With the extended method, we study the well-known differential-difference KP equation, KZ equation and (2+1)-dimensional ANNV system, and both the Lie point symmetry groups and the non-Lie symmetry groups are obtained.

A highly accurate finite-difference method with minimum dispersion error for solving the Helmholtz equation

KAUST Repository

Wu, Zedong

2018-04-05

Numerical simulation of the acoustic wave equation in either isotropic or anisotropic media is crucial to seismic modeling, imaging and inversion. Actually, it represents the core computation cost of these highly advanced seismic processing methods. However, the conventional finite-difference method suffers from severe numerical dispersion errors and S-wave artifacts when solving the acoustic wave equation for anisotropic media. We propose a method to obtain the finite-difference coefficients by comparing its numerical dispersion with the exact form. We find the optimal finite difference coefficients that share the dispersion characteristics of the exact equation with minimal dispersion error. The method is extended to solve the acoustic wave equation in transversely isotropic (TI) media without S-wave artifacts. Numerical examples show that the method is is highly accurate and efficient.

A New Method for Constructing Travelling Wave Solutions to the modified Benjamin–Bona–Mahoney Equation

International Nuclear Information System (INIS)

Jun-Mao, Wang; Miao, Zhang; Wen-Liang, Zhang; Rui, Zhang; Jia-Hua, Han

2008-01-01

We present a new method to find the exact travelling wave solutions of nonlinear evolution equations, with the aid of the symbolic computation. Based on this method, we successfully solve the modified Benjamin–Bona–Mahoney equation, and obtain some new solutions which can be expressed by trigonometric functions and hyperbolic functions. It is shown that the proposed method is direct, effective and can be used for many other nonlinear evolution equations in mathematical physics. (general)

A higher order numerical method for time fractional partial differential equations with nonsmooth data

Science.gov (United States)

Xing, Yanyuan; Yan, Yubin

2018-03-01

Gao et al. [11] (2014) introduced a numerical scheme to approximate the Caputo fractional derivative with the convergence rate O (k 3 - α), 0 equation is sufficiently smooth, Lv and Xu [20] (2016) proved by using energy method that the corresponding numerical method for solving time fractional partial differential equation has the convergence rate O (k 3 - α), 0 equation has low regularity and in this case the numerical method fails to have the convergence rate O (k 3 - α), 0 quadratic interpolation polynomials. Based on this scheme, we introduce a time discretization scheme to approximate the time fractional partial differential equation and show by using Laplace transform methods that the time discretization scheme has the convergence rate O (k 3 - α), 0 0 for smooth and nonsmooth data in both homogeneous and inhomogeneous cases. Numerical examples are given to show that the theoretical results are consistent with the numerical results.

A Method for Solving the Voltage and Torque Equations of the Split ...

African Journals Online (AJOL)

Akorede

v′ Voltage applied across the d – axis rotor winding referred ... The embedded MATLAB function and other useful blocks from the ... III. EQUATIONS OF THE SPLIT PHASE INDUCTION MOTOR. The voltage, flux and electromagnetic torque equations are ..... of single phase induction motor using frequency control method ...

An efficient numerical method for solving the Boltzmann equation in multidimensions

Science.gov (United States)

Dimarco, Giacomo; Loubère, Raphaël; Narski, Jacek; Rey, Thomas

2018-01-01

In this paper we deal with the extension of the Fast Kinetic Scheme (FKS) (Dimarco and Loubère, 2013 [26]) originally constructed for solving the BGK equation, to the more challenging case of the Boltzmann equation. The scheme combines a robust and fast method for treating the transport part based on an innovative Lagrangian technique supplemented with conservative fast spectral schemes to treat the collisional operator by means of an operator splitting approach. This approach along with several implementation features related to the parallelization of the algorithm permits to construct an efficient simulation tool which is numerically tested against exact and reference solutions on classical problems arising in rarefied gas dynamic. We present results up to the 3 D × 3 D case for unsteady flows for the Variable Hard Sphere model which may serve as benchmark for future comparisons between different numerical methods for solving the multidimensional Boltzmann equation. For this reason, we also provide for each problem studied details on the computational cost and memory consumption as well as comparisons with the BGK model or the limit model of compressible Euler equations.

Loop equations and bootstrap methods in the lattice

Directory of Open Access Journals (Sweden)

Peter D. Anderson

2017-08-01

Full Text Available Pure gauge theories can be formulated in terms of Wilson Loops by means of the loop equation. In the large-N limit this equation closes in the expectation value of single loops. In particular, using the lattice as a regulator, it becomes a well defined equation for a discrete set of loops. In this paper we study different numerical approaches to solving this equation. Previous ideas gave good results in the strong coupling region. Here we propose an alternative method based on the observation that certain matrices ρˆ of Wilson loop expectation values are positive definite. They also have unit trace (ρˆ⪰0,Trρˆ=1, in fact they can be defined as reduced density matrices in the space of open loops after tracing over color indices and can be used to define an entropy associated with the loss of information due to such trace SWL=−Tr[ρˆln⁡ρˆ]. The condition that such matrices are positive definite allows us to study the weak coupling region which is relevant for the continuum limit. In the exactly solvable case of two dimensions this approach gives very good results by considering just a few loops. In four dimensions it gives good results in the weak coupling region and therefore is complementary to the strong coupling expansion. We compare the results with standard Monte Carlo simulations.

Application of Exp-function method to potential Kadomtsev-Petviashvili equation

International Nuclear Information System (INIS)

Xian Daquan; Dai Zhengde

2009-01-01

Exact periodic kink-wave solution, periodic soliton and doubly periodic solutions for the potential Kadomtsev-Petviashvii (PKP) equation are obtained using Exp-function method with the help of Maple computation.

Domain decomposition method for solving the neutron diffusion equation

International Nuclear Information System (INIS)

Coulomb, F.

1989-03-01

The aim of this work is to study methods for solving the neutron diffusion equation; we are interested in methods based on a classical finite element discretization and well suited for use on parallel computers. Domain decomposition methods seem to answer this preoccupation. This study deals with a decomposition of the domain. A theoretical study is carried out for Lagrange finite elements and some examples are given; in the case of mixed dual finite elements, the study is based on examples [fr

The generalized tanh method to obtain exact solutions of nonlinear partial differential equation

OpenAIRE

Gómez, César

2007-01-01

In this paper, we present the generalized tanh method to obtain exact solutions of nonlinear partial differential equations, and we obtain solitons and exact solutions of some important equations of the mathematical physics.

Green`s function of Maxwell`s equations and corresponding implications for iterative methods

Energy Technology Data Exchange (ETDEWEB)

Singer, B.S. [Macquarie Univ., Sydney (Australia); Fainberg, E.B. [Inst. of Physics of the Earth, Moscow (Russian Federation)

1996-12-31

Energy conservation law imposes constraints on the norm and direction of the Hilbert space vector representing a solution of Maxwell`s equations. In this paper, we derive these constrains and discuss the corresponding implications for the Green`s function of Maxwell`s equations in a dissipative medium. It is shown that Maxwell`s equations can be reduced to an integral equation with a contracting kernel. The equation can be solved using simple iterations. Software based on this algorithm have successfully been applied to a wide range of problems dealing with high contrast models. The matrix corresponding to the integral equation has a well defined spectrum. The equation can be symmetrized and solved using different approaches, for instance one of the conjugate gradient methods.

T-Stability of the Heun Method and Balanced Method for Solving Stochastic Differential Delay Equations

Directory of Open Access Journals (Sweden)

Xiaolin Zhu

2014-01-01

Full Text Available This paper studies the T-stability of the Heun method and balanced method for solving stochastic differential delay equations (SDDEs. Two T-stable conditions of the Heun method are obtained for two kinds of linear SDDEs. Moreover, two conditions under which the balanced method is T-stable are obtained for two kinds of linear SDDEs. Some numerical examples verify the theoretical results proposed.

Solution of second order linear fuzzy difference equation by Lagrange's multiplier method

Directory of Open Access Journals (Sweden)

Sankar Prasad Mondal

2016-06-01

Full Text Available In this paper we execute the solution procedure for second order linear fuzzy difference equation by Lagrange's multiplier method. In crisp sense the difference equation are easy to solve, but when we take in fuzzy sense it forms a system of difference equation which is not so easy to solve. By the help of Lagrange's multiplier we can solved it easily. The results are illustrated by two different numerical examples and followed by two applications.

New generalized and improved (G′/G-expansion method for nonlinear evolution equations in mathematical physics

Directory of Open Access Journals (Sweden)

Hasibun Naher

2014-10-01

Full Text Available In this article, new extension of the generalized and improved (G′/G-expansion method is proposed for constructing more general and a rich class of new exact traveling wave solutions of nonlinear evolution equations. To demonstrate the novelty and motivation of the proposed method, we implement it to the Korteweg-de Vries (KdV equation. The new method is oriented toward the ease of utilize and capability of computer algebraic system and provides a more systematic, convenient handling of the solution process of nonlinear equations. Further, obtained solutions disclose a wider range of applicability for handling a large variety of nonlinear partial differential equations.

Creep-fatigue life prediction method using Diercks equation for Cr-Mo steel

International Nuclear Information System (INIS)

Sonoya, Keiji; Nonaka, Isamu; Kitagawa, Masaki

1990-01-01

For dealing with the situation that creep-fatigue life properties of materials do not exist, a development of the simple method to predict creep-fatigue life properties is necessary. A method to predict the creep-fatigue life properties of Cr-Mo steels is proposed on the basis of D. Diercks equation which correlates the creep-fatigue lifes of SUS 304 steels under various temperatures, strain ranges, strain rates and hold times. The accuracy of the proposed method was compared with that of the existing methods. The following results were obtained. (1) Fatigue strength and creep rupture strength of Cr-Mo steel are different from those of SUS 304 steel. Therefore in order to apply Diercks equation to creep-fatigue prediction for Cr-Mo steel, the difference of fatigue strength was found to be corrected by fatigue life ratio of both steels and the difference of creep rupture strength was found to be corrected by the equivalent temperature corresponding to equal strength of both steels. (2) Creep-fatigue life can be predicted by the modified Diercks equation within a factor of 2 which is nearly as precise as the accuracy of strain range partitioning method. Required test and analysis procedure of this method are not so complicated as strain range partitioning method. (author)

Modulating functions-based method for parameters and source estimation in one-dimensional partial differential equations

KAUST Repository

Asiri, Sharefa M.

2016-10-20

In this paper, modulating functions-based method is proposed for estimating space–time-dependent unknowns in one-dimensional partial differential equations. The proposed method simplifies the problem into a system of algebraic equations linear in unknown parameters. The well-posedness of the modulating functions-based solution is proved. The wave and the fifth-order KdV equations are used as examples to show the effectiveness of the proposed method in both noise-free and noisy cases.

Resolution of the neutron transport equation by a three-dimensional least square method

International Nuclear Information System (INIS)

Varin, Elisabeth

2001-01-01

The knowledge of space and time distribution of neutrons with a certain energy or speed allows the exploitation and control of a nuclear reactor and the assessment of the irradiation dose about an irradiated nuclear fuel storage site. The neutron density is described by a transport equation. The objective of this research thesis is to develop a software for the resolution of this stationary equation in a three-dimensional Cartesian domain by means of a deterministic method. After a presentation of the transport equation, the author gives an overview of the different deterministic resolution approaches, identifies their benefits and drawbacks, and discusses the choice of the Ressel method. The least square method is precisely described and then applied. Numerical benchmarks are reported for validation purposes

Dhage Iteration Method for Generalized Quadratic Functional Integral Equations

Directory of Open Access Journals (Sweden)

Bapurao C. Dhage

2015-01-01

Full Text Available In this paper we prove the existence as well as approximations of the solutions for a certain nonlinear generalized quadratic functional integral equation. An algorithm for the solutions is developed and it is shown that the sequence of successive approximations starting at a lower or upper solution converges monotonically to the solutions of related quadratic functional integral equation under some suitable mixed hybrid conditions. We rely our main result on Dhage iteration method embodied in a recent hybrid fixed point theorem of Dhage (2014 in partially ordered normed linear spaces. An example is also provided to illustrate the abstract theory developed in the paper.

Soliton-like solutions to the GKdV equation by extended mapping method

International Nuclear Information System (INIS)

Wu Ranchao; Sun Jianhua

2007-01-01

In this note, many new exact solutions of the generalized KdV equation, such as rational solutions, periodic solutions like Jacobian elliptic and triangular functions, soliton-like solutions, are constructed by symbolic computation and the extended mapping method, with the auxiliary ordinary equation replaced by a more general one

Formulation and application of optimal homotopty asymptotic method to coupled differential-difference equations.

Science.gov (United States)

Ullah, Hakeem; Islam, Saeed; Khan, Ilyas; Shafie, Sharidan; Fiza, Mehreen

2015-01-01

In this paper we applied a new analytic approximate technique Optimal Homotopy Asymptotic Method (OHAM) for treatment of coupled differential-difference equations (DDEs). To see the efficiency and reliability of the method, we consider Relativistic Toda coupled nonlinear differential-difference equation. It provides us a convenient way to control the convergence of approximate solutions when it is compared with other methods of solution found in the literature. The obtained solutions show that OHAM is effective, simpler, easier and explicit.

New numerical method for solving the solute transport equation

International Nuclear Information System (INIS)

Ross, B.; Koplik, C.M.

1978-01-01

The solute transport equation can be solved numerically by approximating the water flow field by a network of stream tubes and using a Green's function solution within each stream tube. Compared to previous methods, this approach permits greater computational efficiency and easier representation of small discontinuities, and the results are easier to interpret physically. The method has been used to study hypothetical sites for disposal of high-level radioactive waste

Solution of the Multigroup-Diffusion equation by the response matrix method

International Nuclear Information System (INIS)

Oliveira, C.R.E.

1980-10-01

A preliminary analysis of the response matrix method is made, considering its application to the solution of the multigroup diffusion equations. The one-dimensional formulation is presented and used to test some flux expansions, seeking the application of the method to the two-dimensional problem. This formulation also solves the equations that arise from the integro-differential synthesis algorithm. The slow convergence of the power method, used to solve the eigenvalue problem, and its acceleration by means of the Chebyshev polynomial method, are also studied. An algorithm for the estimation of the dominance ratio is presented, based on the residues of two successive iteration vectors. This ratio, which is not known a priori, is fundamental for the efficiency of the method. Some numerical problems are solved, testing the 1D formulation of the response matrix method, its application to the synthesis algorithm and also, at the same time, the algorithm to accelerate the source problem. (Author) [pt

On a method of construction of exact solutions for equations of two-dimensional hydrodynamics of incompressible liquids

International Nuclear Information System (INIS)

Yurov, A.V.; Yurova, A.A.

2006-01-01

The simple algebraic method for construction of exact solutions of two-dimensional hydrodynamic equations of incompressible flow is proposed. This method can be applied both to nonviscous flow (Euler equations) and to viscous flow (Navier-Stokes equations). In the case of nonviscous flow, the problem is reduced to sequential solving of three linear partial differential equations. In the case of viscous flow, the Navier-Stokes equations are reduced to three linear partial differential equations and one differential equation of the first order [ru

Fourth-order numerical solutions of diffusion equation by using SOR method with Crank-Nicolson approach

Science.gov (United States)

Muhiddin, F. A.; Sulaiman, J.

2017-09-01

The aim of this paper is to investigate the effectiveness of the Successive Over-Relaxation (SOR) iterative method by using the fourth-order Crank-Nicolson (CN) discretization scheme to derive a five-point Crank-Nicolson approximation equation in order to solve diffusion equation. From this approximation equation, clearly, it can be shown that corresponding system of five-point approximation equations can be generated and then solved iteratively. In order to access the performance results of the proposed iterative method with the fourth-order CN scheme, another point iterative method which is Gauss-Seidel (GS), also presented as a reference method. Finally the numerical results obtained from the use of the fourth-order CN discretization scheme, it can be pointed out that the SOR iterative method is superior in terms of number of iterations, execution time, and maximum absolute error.

Finite difference method for inner-layer equations in the resistive MagnetoHydroDynamic stability analysis

International Nuclear Information System (INIS)

Tokuda, Shinji; Watanabe, Tomoko.

1996-08-01

The matching problem in resistive MagnetoHydroDynamic stability analysis by the asymptotic matching method has been reformulated as an initial-boundary value problem for the inner-layer equations describing the plasma dynamics in the thin layer around a rational surface. The third boundary conditions at boundaries of a finite interval are imposed on the inner layer equations in the formulation instead of asymptotic conditions at infinities. The finite difference method for this problem has been applied to model equations whose solutions are known in a closed form. It has been shown that the initial value problem and the associated eigenvalue problem for the model equations can be solved by the finite difference method with numerical stability. The formulation presented here enables the asymptotic matching method to be a practical method for the resistive MHD stability analysis. (author)

Local Fractional Series Expansion Method for Solving Wave and Diffusion Equations on Cantor Sets

Directory of Open Access Journals (Sweden)

Ai-Min Yang

2013-01-01

Full Text Available We proposed a local fractional series expansion method to solve the wave and diffusion equations on Cantor sets. Some examples are given to illustrate the efficiency and accuracy of the proposed method to obtain analytical solutions to differential equations within the local fractional derivatives.

Application of coarse-mesh methods to fluid dynamics equations

International Nuclear Information System (INIS)

Romstedt, P.; Werner, W.

1977-01-01

An Asymmetric Weighted Residual (ASWR) method for fluid dynamics equations is described. It leads to local operators with a 7-point Finite Difference (FD) structure, which is independent of the degree of the approximating polynomials. An 1-dimensional problem was solved by both this ASWR-method and a commonly used FD-method. The numerical results demonstrate that the ASWR-method combines high accuracy on a coarse computational mesh with short computing time per space point. The posibility of using fewer space points consequently brings about a considerable reduction in total running time for the ASWR-method as compared with conventional FD-methods. (orig.) [de

Two modified symplectic partitioned Runge-Kutta methods for solving the elastic wave equation

Science.gov (United States)

Su, Bo; Tuo, Xianguo; Xu, Ling

2017-08-01

Based on a modified strategy, two modified symplectic partitioned Runge-Kutta (PRK) methods are proposed for the temporal discretization of the elastic wave equation. The two symplectic schemes are similar in form but are different in nature. After the spatial discretization of the elastic wave equation, the ordinary Hamiltonian formulation for the elastic wave equation is presented. The PRK scheme is then applied for time integration. An additional term associated with spatial discretization is inserted into the different stages of the PRK scheme. Theoretical analyses are conducted to evaluate the numerical dispersion and stability of the two novel PRK methods. A finite difference method is used to approximate the spatial derivatives since the two schemes are independent of the spatial discretization technique used. The numerical solutions computed by the two new schemes are compared with those computed by a conventional symplectic PRK. The numerical results, which verify the new method, are superior to those generated by traditional conventional methods in seismic wave modeling.

On a computer implementation of the block Gauss–Seidel method for normal systems of equations

OpenAIRE

Alexander I. Zhdanov; Ekaterina Yu. Bogdanova

2016-01-01

This article focuses on the modification of the block option Gauss-Seidel method for normal systems of equations, which is a sufficiently effective method of solving generally overdetermined, systems of linear algebraic equations of high dimensionality. The main disadvantage of methods based on normal equations systems is the fact that the condition number of the normal system is equal to the square of the condition number of the original problem. This fact has a negative impact on the rate o...

Generalized conjugate-gradient methods for the Navier-Stokes equations

Science.gov (United States)

Ajmani, Kumud; Ng, Wing-Fai; Liou, Meng-Sing

1991-01-01

A generalized conjugate-gradient method is used to solve the two-dimensional, compressible Navier-Stokes equations of fluid flow. The equations are discretized with an implicit, upwind finite-volume formulation. Preconditioning techniques are incorporated into the new solver to accelerate convergence of the overall iterative method. The superiority of the new solver is demonstrated by comparisons with a conventional line Gauss-Siedel Relaxation solver. Computational test results for transonic flow (trailing edge flow in a transonic turbine cascade) and hypersonic flow (M = 6.0 shock-on-shock phenoena on a cylindrical leading edge) are presented. When applied to the transonic cascade case, the new solver is 4.4 times faster in terms of number of iterations and 3.1 times faster in terms of CPU time than the Relaxation solver. For the hypersonic shock case, the new solver is 3.0 times faster in terms of number of iterations and 2.2 times faster in terms of CPU time than the Relaxation solver.

New Equating Methods and Their Relationships with Levine Observed Score Linear Equating under the Kernel Equating Framework

Science.gov (United States)

Chen, Haiwen; Holland, Paul

2010-01-01

In this paper, we develop a new curvilinear equating for the nonequivalent groups with anchor test (NEAT) design under the assumption of the classical test theory model, that we name curvilinear Levine observed score equating. In fact, by applying both the kernel equating framework and the mean preserving linear transformation of…

Modulating functions-based method for parameters and source estimation in one-dimensional partial differential equations

KAUST Repository

Asiri, Sharefa M.; Laleg-Kirati, Taous-Meriem

2016-01-01

In this paper, modulating functions-based method is proposed for estimating space–time-dependent unknowns in one-dimensional partial differential equations. The proposed method simplifies the problem into a system of algebraic equations linear

Time discretization of the point kinetic equations using matrix exponential method and First-Order Hold

International Nuclear Information System (INIS)

Park, Yujin; Kazantzis, Nikolaos; Parlos, Alexander G.; Chong, Kil To

2013-01-01

Highlights: • Numerical solution for stiff differential equations using matrix exponential method. • The approximation is based on First Order Hold assumption. • Various input examples applied to the point kinetics equations. • The method shows superior useful and effective activity. - Abstract: A system of nonlinear differential equations is derived to model the dynamics of neutron density and the delayed neutron precursors within a point kinetics equation modeling framework for a nuclear reactor. The point kinetic equations are mathematically characterized as stiff, occasionally nonlinear, ordinary differential equations, posing significant challenges when numerical solutions are sought and traditionally resulting in the need for smaller time step intervals within various computational schemes. In light of the above realization, the present paper proposes a new discretization method inspired by system-theoretic notions and technically based on a combination of the matrix exponential method (MEM) and the First-Order Hold (FOH) assumption. Under the proposed time discretization structure, the sampled-data representation of the nonlinear point kinetic system of equations is derived. The performance of the proposed time discretization procedure is evaluated using several case studies with sinusoidal reactivity profiles and multiple input examples (reactivity and neutron source function). It is shown, that by applying the proposed method under a First-Order Hold for the neutron density and the precursor concentrations at each time step interval, the stiffness problem associated with the point kinetic equations can be adequately addressed and resolved. Finally, as evidenced by the aforementioned detailed simulation studies, the proposed method retains its validity and accuracy for a wide range of reactor operating conditions, including large sampling periods dictated by physical and/or technical limitations associated with the current state of sensor and

Telescopic projective methods for parabolic differential equations

CERN Document Server

Gear, C W

2003-01-01

Projective methods were introduced in an earlier paper [C.W. Gear, I.G. Kevrekidis, Projective Methods for Stiff Differential Equations: problems with gaps in their eigenvalue spectrum, NEC Research Institute Report 2001-029, available from http://www.neci.nj.nec.com/homepages/cwg/projective.pdf Abbreviated version to appear in SISC] as having potential for the efficient integration of problems with a large gap between two clusters in their eigenvalue spectrum, one cluster containing eigenvalues corresponding to components that have already been damped in the numerical solution and one corresponding to components that are still active. In this paper we introduce iterated projective methods that allow for explicit integration of stiff problems that have a large spread of eigenvalues with no gaps in their spectrum as arise in the semi-discretization of PDEs with parabolic components.

Telescopic projective methods for parabolic differential equations

International Nuclear Information System (INIS)

Gear, C.W.; Kevrekidis, Ioannis G.

2003-01-01

Projective methods were introduced in an earlier paper [C.W. Gear, I.G. Kevrekidis, Projective Methods for Stiff Differential Equations: problems with gaps in their eigenvalue spectrum, NEC Research Institute Report 2001-029, available from http://www.neci.nj.nec.com/homepages/cwg/projective.pdf Abbreviated version to appear in SISC] as having potential for the efficient integration of problems with a large gap between two clusters in their eigenvalue spectrum, one cluster containing eigenvalues corresponding to components that have already been damped in the numerical solution and one corresponding to components that are still active. In this paper we introduce iterated projective methods that allow for explicit integration of stiff problems that have a large spread of eigenvalues with no gaps in their spectrum as arise in the semi-discretization of PDEs with parabolic components

Analysis of global multiscale finite element methods for wave equations with continuum spatial scales

KAUST Repository

Jiang, Lijian; Efendiev, Yalchin; Ginting, Victor

2010-01-01

In this paper, we discuss a numerical multiscale approach for solving wave equations with heterogeneous coefficients. Our interest comes from geophysics applications and we assume that there is no scale separation with respect to spatial variables. To obtain the solution of these multiscale problems on a coarse grid, we compute global fields such that the solution smoothly depends on these fields. We present a Galerkin multiscale finite element method using the global information and provide a convergence analysis when applied to solve the wave equations. We investigate the relation between the smoothness of the global fields and convergence rates of the global Galerkin multiscale finite element method for the wave equations. Numerical examples demonstrate that the use of global information renders better accuracy for wave equations with heterogeneous coefficients than the local multiscale finite element method. © 2010 IMACS.

Analysis of global multiscale finite element methods for wave equations with continuum spatial scales

KAUST Repository

Jiang, Lijian

2010-08-01

In this paper, we discuss a numerical multiscale approach for solving wave equations with heterogeneous coefficients. Our interest comes from geophysics applications and we assume that there is no scale separation with respect to spatial variables. To obtain the solution of these multiscale problems on a coarse grid, we compute global fields such that the solution smoothly depends on these fields. We present a Galerkin multiscale finite element method using the global information and provide a convergence analysis when applied to solve the wave equations. We investigate the relation between the smoothness of the global fields and convergence rates of the global Galerkin multiscale finite element method for the wave equations. Numerical examples demonstrate that the use of global information renders better accuracy for wave equations with heterogeneous coefficients than the local multiscale finite element method. © 2010 IMACS.

Optimal Homotopy Asymptotic Method for Solving System of Fredholm Integral Equations

Directory of Open Access Journals (Sweden)

Bahman Ghazanfari

2013-08-01

Full Text Available In this paper, optimal homotopy asymptotic method (OHAM is applied to solve system of Fredholm integral equations. The effectiveness of optimal homotopy asymptotic method is presented. This method provides easy tools to control the convergence region of approximating solution series wherever necessary. The results of OHAM are compared with homotopy perturbation method (HPM and Taylor series expansion method (TSEM.

Formulation and Application of Optimal Homotopty Asymptotic Method to Coupled Differential - Difference Equations

Science.gov (United States)

Ullah, Hakeem; Islam, Saeed; Khan, Ilyas; Shafie, Sharidan; Fiza, Mehreen

2015-01-01

In this paper we applied a new analytic approximate technique Optimal Homotopy Asymptotic Method (OHAM) for treatment of coupled differential- difference equations (DDEs). To see the efficiency and reliability of the method, we consider Relativistic Toda coupled nonlinear differential-difference equation. It provides us a convenient way to control the convergence of approximate solutions when it is compared with other methods of solution found in the literature. The obtained solutions show that OHAM is effective, simpler, easier and explicit. PMID:25874457

Numerical solution to generalized Burgers'-Fisher equation using Exp-function method hybridized with heuristic computation.

Science.gov (United States)

Malik, Suheel Abdullah; Qureshi, Ijaz Mansoor; Amir, Muhammad; Malik, Aqdas Naveed; Haq, Ihsanul

2015-01-01

In this paper, a new heuristic scheme for the approximate solution of the generalized Burgers'-Fisher equation is proposed. The scheme is based on the hybridization of Exp-function method with nature inspired algorithm. The given nonlinear partial differential equation (NPDE) through substitution is converted into a nonlinear ordinary differential equation (NODE). The travelling wave solution is approximated by the Exp-function method with unknown parameters. The unknown parameters are estimated by transforming the NODE into an equivalent global error minimization problem by using a fitness function. The popular genetic algorithm (GA) is used to solve the minimization problem, and to achieve the unknown parameters. The proposed scheme is successfully implemented to solve the generalized Burgers'-Fisher equation. The comparison of numerical results with the exact solutions, and the solutions obtained using some traditional methods, including adomian decomposition method (ADM), homotopy perturbation method (HPM), and optimal homotopy asymptotic method (OHAM), show that the suggested scheme is fairly accurate and viable for solving such problems.

Numerical solution to generalized Burgers'-Fisher equation using Exp-function method hybridized with heuristic computation.

Directory of Open Access Journals (Sweden)

Suheel Abdullah Malik

Full Text Available In this paper, a new heuristic scheme for the approximate solution of the generalized Burgers'-Fisher equation is proposed. The scheme is based on the hybridization of Exp-function method with nature inspired algorithm. The given nonlinear partial differential equation (NPDE through substitution is converted into a nonlinear ordinary differential equation (NODE. The travelling wave solution is approximated by the Exp-function method with unknown parameters. The unknown parameters are estimated by transforming the NODE into an equivalent global error minimization problem by using a fitness function. The popular genetic algorithm (GA is used to solve the minimization problem, and to achieve the unknown parameters. The proposed scheme is successfully implemented to solve the generalized Burgers'-Fisher equation. The comparison of numerical results with the exact solutions, and the solutions obtained using some traditional methods, including adomian decomposition method (ADM, homotopy perturbation method (HPM, and optimal homotopy asymptotic method (OHAM, show that the suggested scheme is fairly accurate and viable for solving such problems.

A new generalized algebra method and its application in the (2 + 1) dimensional Boiti-Leon-Pempinelli equation

International Nuclear Information System (INIS)

Ren Yujie; Liu Shutian; Zhang Hongqing

2007-01-01

In the present paper, some types of general solutions of a first-order nonlinear ordinary differential equation with six degree are given and a new generalized algebra method is presented to find more exact solutions of nonlinear differential equations. As an application of the method and the solutions of this equation, we choose the (2 + 1) dimensional Boiti Leon Pempinelli equation to illustrate the validity and advantages of the method. As a consequence, more new types and general solutions are found which include rational solutions and irrational solutions and so on. The new method can also be applied to other nonlinear differential equations in mathematical physics

Weak Second Order Explicit Stabilized Methods for Stiff Stochastic Differential Equations

KAUST Repository

Abdulle, Assyr; Vilmart, Gilles; Zygalakis, Konstantinos C.

2013-01-01

We introduce a new family of explicit integrators for stiff Itô stochastic differential equations (SDEs) of weak order two. These numerical methods belong to the class of one-step stabilized methods with extended stability domains and do not suffer

Steady-state transport equation resolution by particle methods, and numerical results

International Nuclear Information System (INIS)

Mercier, B.

1985-10-01

A method to solve steady-state transport equation has been given. Principles of the method are given. The method is studied in two different cases; estimations given by the theory are compared to numerical results. Results got in 1-D (spherical geometry) and in 2-D (axisymmetric geometry) are given [fr

A Parameter Robust Method for Singularly Perturbed Delay Differential Equations

Directory of Open Access Journals (Sweden)

Erdogan Fevzi

2010-01-01

Full Text Available Uniform finite difference methods are constructed via nonstandard finite difference methods for the numerical solution of singularly perturbed quasilinear initial value problem for delay differential equations. A numerical method is constructed for this problem which involves the appropriate Bakhvalov meshes on each time subinterval. The method is shown to be uniformly convergent with respect to the perturbation parameter. A numerical example is solved using the presented method, and the computed result is compared with exact solution of the problem.

Legendre Wavelet Operational Matrix Method for Solution of Riccati Differential Equation

Directory of Open Access Journals (Sweden)

S. Balaji

2014-01-01

Full Text Available A Legendre wavelet operational matrix method (LWM is presented for the solution of nonlinear fractional-order Riccati differential equations, having variety of applications in quantum chemistry and quantum mechanics. The fractional-order Riccati differential equations converted into a system of algebraic equations using Legendre wavelet operational matrix. Solutions given by the proposed scheme are more accurate and reliable and they are compared with recently developed numerical, analytical, and stochastic approaches. Comparison shows that the proposed LWM approach has a greater performance and less computational effort for getting accurate solutions. Further existence and uniqueness of the proposed problem are given and moreover the condition of convergence is verified.

Numerical Methods for a Class of Differential Algebraic Equations

Directory of Open Access Journals (Sweden)

Lei Ren

2017-01-01

Full Text Available This paper is devoted to the study of some efficient numerical methods for the differential algebraic equations (DAEs. At first, we propose a finite algorithm to compute the Drazin inverse of the time varying DAEs. Numerical experiments are presented by Drazin inverse and Radau IIA method, which illustrate that the precision of the Drazin inverse method is higher than the Radau IIA method. Then, Drazin inverse, Radau IIA, and Padé approximation are applied to the constant coefficient DAEs, respectively. Numerical results demonstrate that the Padé approximation is powerful for solving constant coefficient DAEs.

Approximate solutions for the two-dimensional integral transport equation. The critically mixed methods of resolution

International Nuclear Information System (INIS)

Sanchez, Richard.

1980-11-01

This work is divided into two part the first part (note CEA-N-2165) deals with the solution of complex two-dimensional transport problems, the second one treats the critically mixed methods of resolution. These methods are applied for one-dimensional geometries with highly anisotropic scattering. In order to simplify the set of integral equation provided by the integral transport equation, the integro-differential equation is used to obtain relations that allow to lower the number of integral equation to solve; a general mathematical and numerical study is presented [fr

Convergence acceleration of Navier-Stokes equation using adaptive wavelet method