Numerical study of a mathematical model of internal erosion of soil

Science.gov (United States)

Sibin, A.

2017-10-01

The process of internal erosion in a three-phase saturated soil is studied. A mathematical model describing the process consists of the equations of mass conservation, Darcy’s law and equation for capillary pressure. The original system of equations is reduced to a system of three equations for porosity, pressure and water saturation. Obtained equation for the water saturation is degenerate. The degenerate problem in an one-dimensional domain is solved numerically using the finite-difference method.

Mathematical and numerical analysis of hyper-elastic systems and introduction of plasticity

International Nuclear Information System (INIS)

Kluth, G.

2008-12-01

The goal is to model mathematically and numerically the dynamic phenomenons for solids in finite plasticity. We suggest a model that we call hyper-elasto-plastic based on hyper-elastic systems of conservation laws and on the use of an equation of state that we have constructed so as to achieve the plastic yield criterion of Von Mises. This model gives exact (analytic) solutions with shock split to flyer-plate experiments. The mathematical analysis of this model is done (hyperbolicity, characteristic fields, involutions and entropy). In the numerical part, we give 1D and 2D Lagrangian schemes which satisfy an entropy criterion. Moreover, thanks to a special discretization of the equations on deformation gradient, we satisfy some discrete involutions. In this work, the degeneracy of the solid model into hydrodynamic models is studied at the continuous level, and achieved at the numerical one. On different problems, we show the validity of our model and our numerical schemes. (author)

Partial Differential Equations Modeling and Numerical Simulation

CERN Document Server

Glowinski, Roland

2008-01-01

This book is dedicated to Olivier Pironneau. For more than 250 years partial differential equations have been clearly the most important tool available to mankind in order to understand a large variety of phenomena, natural at first and then those originating from human activity and technological development. Mechanics, physics and their engineering applications were the first to benefit from the impact of partial differential equations on modeling and design, but a little less than a century ago the Schrödinger equation was the key opening the door to the application of partial differential equations to quantum chemistry, for small atomic and molecular systems at first, but then for systems of fast growing complexity. Mathematical modeling methods based on partial differential equations form an important part of contemporary science and are widely used in engineering and scientific applications. In this book several experts in this field present their latest results and discuss trends in the numerical analy...

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.

Numerical solution of ordinary differential equations. For classical, relativistic and nano systems

International Nuclear Information System (INIS)

Greenspan, D.

2006-01-01

An up-to-date survey on numerical solutions with theory, intuition and applications. Ordinary differential equations (ODE) play a significant role in mathematics, physics and engineering sciences, and thus are part of relevant college and university courses. Many problems, however, both traditional and modern, do not possess exact solutions, and must be treated numerically. Usually this is done with software packages, but for this to be efficient requires a sound understanding of the mathematics involved. This work meets the need for an affordable textbook that helps in understanding numerical solutions of ODE. Carefully structured by an experienced textbook author, it provides a survey of ODE for various applications, both classical and modern, including such special applications as relativistic and nano systems. The examples are carefully explained and compiled into an algorithm, each of which is presented generically, independent of a specific programming language, while each chapter is rounded off with exercises. The text meets the demands of MA200 courses and of the newly created Numerical Solution of Differential Equations courses, making it ideal for both students and lecturers in physics, mathematics, mechanical engineering, electrical engineering, as well as for physicists, mathematicians, engineers, and electrical engineers. From the Contents - Euler's Method - Runge-Kutta Methods - The Method of Taylor Expansions - Large Second Order Systems with Application to Nano Systems - Completely Conservative, Covariant Numerical Methodology - Instability - Numerical Solution of Tridiagonal Linear Algebraic Systems and Related Nonlinear Systems - Approximate Solution of Boundary Value Problems - Special Relativistic Motion - Special Topics - Appendix: Basic Matrix Operations - Bibliography. (orig.) (orig.)

Large quantum systems: a mathematical and numerical perspective

International Nuclear Information System (INIS)

Lewin, M.

2009-06-01

This thesis is devoted to the mathematical study of variational models for large quantum systems. The mathematical methods are that of nonlinear analysis, calculus of variations, partial differential equations, spectral theory, and numerical analysis. The first part contains some results on finite systems. We study several approximations of the N-body Schroedinger equation for electrons in an atom or a molecule, and then the so-called Hartree-Fock- Bogoliubov model for a system of fermions interacting via the gravitational force. In a second part, we propose a new method allowing to prove the existence of the thermodynamic limit of Coulomb quantum systems. Then, we construct two Hartree-Fock-type models for infinite systems. The first is a relativistic theory deduced from Quantum Electrodynamics, allowing to describe the behavior of electrons, coupled to that of Dirac's vacuum which can become polarized. The second model describes a nonrelativistic quantum crystal in the presence of a charged defect. A new numerical method is also proposed. The last part of the thesis is devoted to spectral pollution, a phenomenon which is observed when trying to approximate eigenvalues in a gap of the essential spectrum of a self-adjoint operator, for instance for periodic Schroedinger operators or Dirac operators. (author)

Partial differential equations of mathematical physics

CERN Document Server

Sobolev, S L

1964-01-01

Partial Differential Equations of Mathematical Physics emphasizes the study of second-order partial differential equations of mathematical physics, which is deemed as the foundation of investigations into waves, heat conduction, hydrodynamics, and other physical problems. The book discusses in detail a wide spectrum of topics related to partial differential equations, such as the theories of sets and of Lebesgue integration, integral equations, Green's function, and the proof of the Fourier method. Theoretical physicists, experimental physicists, mathematicians engaged in pure and applied math

Mathematics and Maxwell's equations

International Nuclear Information System (INIS)

Boozer, Allen H

2010-01-01

The universality of mathematics and Maxwell's equations is not shared by specific plasma models. Computations become more reliable, efficient and transparent if specific plasma models are used to obtain only the information that would otherwise be missing. Constraints of high universality, such as those from mathematics and Maxwell's equations, can be obscured or lost by integrated computations. Recognition of subtle constraints of high universality is important for (1) focusing the design of control systems for magnetic field errors in tokamaks from perturbations that have little effect on the plasma to those that do, (2) clarifying the limits of applicability to astrophysics of computations of magnetic reconnection in fields that have a double periodicity or have B-vector =0 on a surface, as in a Harris sheet. Both require a degree of symmetry not expected in natural systems. Mathematics and Maxwell's equations imply that neighboring magnetic field lines characteristically separate exponentially with distance along a line. This remarkably universal phenomenon has been largely ignored, though it defines a trigger for reconnection through a critical magnitude of exponentiation. These and other examples of the importance of making distinctions and understanding constraints of high universality are explained.

Review on mathematical basis for thermal conduction equation

Energy Technology Data Exchange (ETDEWEB)

Park, D. G.; Kim, H. M

2007-10-15

In the view point of thermal conductivity measurement technology, It is very useful to understand mathematical theory of thermal conduction equation in order to evaluation of measurement data and to solve diverse technical problem in measurement. To approach this mathematical theory, thermal conduction equation is derived by Fourier thermal conduction law. Since thermal conduction equation depends on the Lapacian operator basically, mathematical meaning of Lapalacian and various diffusion equation including Laplacian have been studied. Stum-Liouville problem and Bessel function were studied in this report to understand analytical solution of various diffusion equation.

Review on mathematical basis for thermal conduction equation

International Nuclear Information System (INIS)

Park, D. G.; Kim, H. M.

2007-10-01

In the view point of thermal conductivity measurement technology, It is very useful to understand mathematical theory of thermal conduction equation in order to evaluation of measurement data and to solve diverse technical problem in measurement. To approach this mathematical theory, thermal conduction equation is derived by Fourier thermal conduction law. Since thermal conduction equation depends on the Lapacian operator basically, mathematical meaning of Lapalacian and various diffusion equation including Laplacian have been studied. Stum-Liouville problem and Bessel function were studied in this report to understand analytical solution of various diffusion equation

A numerical methodology for the Painlevé equations

KAUST Repository

Fornberg, Bengt

2011-07-01

The six Painlevé transcendents PI-PVI have both applications and analytic properties that make them stand out from most other classes of special functions. Although they have been the subject of extensive theoretical investigations for about a century, they still have a reputation for being numerically challenging. In particular, their extensive pole fields in the complex plane have often been perceived as \\'numerical mine fields\\'. In the present work, we note that the Painlevé property in fact provides the opportunity for very fast and accurate numerical solutions throughout such fields. When combining a Taylor/Padé-based ODE initial value solver for the pole fields with a boundary value solver for smooth regions, numerical solutions become available across the full complex plane. We focus here on the numerical methodology, and illustrate it for the PI equation. In later studies, we will concentrate on mathematical aspects of both the PI and the higher Painlevé transcendents. © 2011 Elsevier Inc.

Attractors for equations of mathematical physics

CERN Document Server

Chepyzhov, Vladimir V

2001-01-01

One of the major problems in the study of evolution equations of mathematical physics is the investigation of the behavior of the solutions to these equations when time is large or tends to infinity. The related important questions concern the stability of solutions or the character of the instability if a solution is unstable. In the last few decades, considerable progress in this area has been achieved in the study of autonomous evolution partial differential equations. For a number of basic evolution equations of mathematical physics, it was shown that the long time behavior of their soluti

DEVELOPMENT OF THE MATHEMATICAL INTUITION OF STUDENTS IN TRAINING THE INVERSE PROBLEMS FOR DIFFERENTIAL EQUATIONS

Directory of Open Access Journals (Sweden)

Виктор Семенович Корнилов

2017-12-01

Full Text Available In article attention to that fact that at students of higher educational institutions of the physical and mathematical and natural-science directions of preparation when training in the reverse tasks for differential equations the mathematical intuition which is an important component of their creative potential develops is paid. The mathematical intuition helps students to comprehend a physical sense of the researched application-oriented task, to select effective methods of mathematical physics for the decision of the reverse task for differential equations.The mathematical intuition of students develops in many respects in case of the decision of different educational jobs. Among such educational jobs: creation of system of integrable equations of the reverse task for differential equations, the proof of the conditional correctness of the decision of the reverse task for differential equations, creation of the difference analog of the reverse task for a differential equation; finding of the numerical decision of the reverse task, the proof of convergence of approximate solution of the reverse task to the exact decision, reasons for the idea of the proof of a correctness (the conditional correctness of the decision of the reverse task for differential equations, a statement of logical outputs of application-oriented or humanitarian character on the basis of the conducted research of the reverse task and other educational jobs.In the course of such training students create system of fundamental knowledge in the field of the reverse and incorrect tasks, acquire new scientific knowledge in the field of applied and calculus mathematics, but, obviously, and develop a mathematical intuition.

Mathematical and numerical methods for Vlasov-Maxwell equations: the contributions of data mining

International Nuclear Information System (INIS)

Assous, F.; Chaskalovic, J.

2014-01-01

There exist a lot of formulations that can model plasma physics or particle accelerators problems as the Vlasov- Maxwell equations. This paper deals with the applications of data mining techniques in the evaluation of numerical solutions of Vlasov-Maxwell models. This is part of the topic of characterizing the model and approximation errors via learning techniques. We give two examples of application. The first one aims at comparing two Vlasov-Maxwell approximate models. In the second one, a scheme based on data mining techniques is proposed to characterize the errors between a P1 and a P2 finite element Particle-In-Cell approach. Beyond these examples, this original approach should operate in all cases where intricate numerical simulations like for the Vlasov-Maxwell equations take a central part. (authors)

assessment of concentration of air pollutants using analytical and numerical solution of the atmospheric diffusion equation

International Nuclear Information System (INIS)

Esmail, S.F.H.

2011-01-01

The mathematical formulation of numerous physical problems a results in differential equations actually partial or ordinary differential equations.In our study we are interested in solutions of partial differential equations.The aim of this work is to calculate the concentrations of the pollution, by solving the atmospheric diffusion equation(ADE) using different mathematical methods of solution. It is difficult to solve the general form of ADE analytically, so we use some assumptions to get its solution.The solutions of it depend on the eddy diffusivity profiles(k) and the wind speed u. We use some physical assumptions to simplify its formula and solve it. In the present work, we solve the ADE analytically in three dimensions using Green's function method, Laplace transform method, normal mode method and these separation of variables method. Also, we use ADM as a numerical method. Finally, comparisons are made with the results predicted by the previous methods and the observed data.

Comparison of Numerical Approaches to a Steady-State Landscape Equation

Science.gov (United States)

Bachman, S.; Peckham, S.

2008-12-01

A mathematical model of an idealized fluvial landscape has been developed, in which a land surface will evolve to preserve dendritic channel networks as the surface is lowered. The physical basis for this model stems from the equations for conservation of mass for water and sediment. These equations relate the divergence of the 2D vector fields showing the unit-width discharge of water and sediment to the excess rainrate and tectonic uplift on the land surface. The 2D flow direction is taken to be opposite to the water- surface gradient vector. These notions are combined with a generalized Manning-type flow resistance formula and a generalized sediment transport law to give a closed mathematical system that can, in principle, be solved for all variables of interest: discharge of water and sediment, land surface height, vertically- averaged flow velocity, water depth, and shear stress. The hydraulic geometry equations (Leopold et. al, 1964, 1995) are used to incorporate width, depth, velocity, and slope of river channels as powers of the mean-annual river discharge. Combined, they give the unit- width discharge of the stream as a power, γ, of the water surface slope. The simplified steady-state model takes into account three components among those listed above: conservation of mass for water, flow opposite the gradient, and a slope-discharge exponent γ = -1 to reflect mature drainage networks. The mathematical representation of this model appears as a second-order hyperbolic partial differential equation (PDE) where the diffusivity is inversely proportional to the square of the local surface slope. The highly nonlinear nature of this PDE has made it very difficult to solve both analytically and numerically. We present simplistic analytic solutions to this equation which are used to test the validity of the numerical algorithms. We also present three such numerical approaches which have been used in solving the differential equation. The first is based on a

Numerical solution of Boltzmann's equation

International Nuclear Information System (INIS)

Sod, G.A.

1976-04-01

The numerical solution of Boltzmann's equation is considered for a gas model consisting of rigid spheres by means of Hilbert's expansion. If only the first two terms of the expansion are retained, Boltzmann's equation reduces to the Boltzmann-Hilbert integral equation. Successive terms in the Hilbert expansion are obtained by solving the same integral equation with a different source term. The Boltzmann-Hilbert integral equation is solved by a new very fast numerical method. The success of the method rests upon the simultaneous use of four judiciously chosen expansions; Hilbert's expansion for the distribution function, another expansion of the distribution function in terms of Hermite polynomials, the expansion of the kernel in terms of the eigenvalues and eigenfunctions of the Hilbert operator, and an expansion involved in solving a system of linear equations through a singular value decomposition. The numerical method is applied to the study of the shock structure in one space dimension. Numerical results are presented for Mach numbers of 1.1 and 1.6. 94 refs, 7 tables, 1 fig

Partial differential equations of mathematical physics and integral equations

CERN Document Server

Guenther, Ronald B

1996-01-01

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

Mathematical modelling and numerical resolution of multi-phase compressible fluid flows problems

International Nuclear Information System (INIS)

Lagoutiere, Frederic

2000-01-01

This work deals with Eulerian compressible multi-species fluid dynamics, the species being either mixed or separated (with interfaces). The document is composed of three parts. The first parts devoted to the numerical resolution of model problems: advection equation, Burgers equation, and Euler equations, in dimensions one and two. The goal is to find a precise method, especially for discontinuous initial conditions, and we develop non dissipative algorithms. They are based on a downwind finite-volume discretization under some stability constraints. The second part treats of the mathematical modelling of fluids mixtures. We construct and analyse a set of multi-temperature and multi-pressure models that are entropy, symmetrizable, hyperbolic, not ever conservative. In the third part, we apply the ideas developed in the first part (downwind discretization) to the numerical resolution of the partial differential problems we have constructed for fluids mixtures in the second part. We present some numerical results in dimensions one and two. (author) [fr

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…

Mathematical and numerical analysis of systems of compressible hydrodynamics and photonics with polar coordinates

International Nuclear Information System (INIS)

Meltz, Bertrand

2015-01-01

This thesis deals with the mathematical and numerical analysis of the systems of compressible hydrodynamics and radiative transfer. More precisely, we study the derivation of numerical methods with 2D polar coordinates (one for the radius, one for the angle) where equations are discretized on regular polar grids. On one hand, these methods are well-suited for the simulation of flows with polar symmetries since they preserve these symmetries by construction. On the other hand, such coordinates systems introduce geometrical singularities as well as geometrical source terms which must be carefully treated. The first part of this document is devoted to the study of hydrodynamics equations, or Euler equations. We propose a new class of arbitrary high-order numerical schemes in both space and time and rely on directional splitting methods for the resolution of 2D equations. Each sub-system is solved using a Lagrange+Remap solver. We study the influence of the r=0 geometrical singularities of the cylindrical and spherical coordinates systems on the precision of the 2D numerical solutions. The second part of this document is devoted to the study of radiative transfer equations. In these equations, the unknowns depend on a large number of variables and a stiff source term is involved. The main difficulty consists in capturing the correct asymptotic behavior on coarse grids. We first construct a class of models where the radiative intensity is projected on a truncated spherical harmonics basis in order to lower the number of mathematical dimensions. Then we propose an Asymptotic Preserving scheme built in polar coordinates and we show that the scheme capture the correct diffusion limit in the radial direction as well as in the polar direction. (author) [fr

A practical course in differential equations and mathematical modeling

CERN Document Server

Ibragimov , Nail H

2009-01-01

A Practical Course in Differential Equations and Mathematical Modelling is a unique blend of the traditional methods of ordinary and partial differential equations with Lie group analysis enriched by the author's own theoretical developments. The book which aims to present new mathematical curricula based on symmetry and invariance principles is tailored to develop analytic skills and working knowledge in both classical and Lie's methods for solving linear and nonlinear equations. This approach helps to make courses in differential equations, mathematical modelling, distributions and fundame

Trends in differential equations and applications

CERN Document Server

Neble, María; Galván, José

2016-01-01

This work collects the most important results presented at the Congress on Differential Equations and Applications/Congress on Applied Mathematics (CEDYA/CMA) in Cádiz (Spain) in 2015. It supports further research in differential equations, numerical analysis, mechanics, control and optimization. In particular, it helps readers gain an overview of specific problems of interest in the current mathematical research related to different branches of applied mathematics. This includes the analysis of nonlinear partial differential equations, exact solutions techniques for ordinary differential equations, numerical analysis and numerical simulation of some models arising in experimental sciences and engineering, control and optimization, and also trending topics on numerical linear Algebra, dynamical systems, and applied mathematics for Industry. This volume is mainly addressed to any researcher interested in the applications of mathematics, especially in any subject mentioned above. It may be also useful to PhD s...

Use of artificial bee colonies algorithm as numerical approximation of differential equations solution

Science.gov (United States)

Fikri, Fariz Fahmi; Nuraini, Nuning

2018-03-01

The differential equation is one of the branches in mathematics which is closely related to human life problems. Some problems that occur in our life can be modeled into differential equations as well as systems of differential equations such as the Lotka-Volterra model and SIR model. Therefore, solving a problem of differential equations is very important. Some differential equations are difficult to solve, so numerical methods are needed to solve that problems. Some numerical methods for solving differential equations that have been widely used are Euler Method, Heun Method, Runge-Kutta and others. However, some of these methods still have some restrictions that cause the method cannot be used to solve more complex problems such as an evaluation interval that we cannot change freely. New methods are needed to improve that problems. One of the method that can be used is the artificial bees colony algorithm. This algorithm is one of metaheuristic algorithm method, which can come out from local search space and do exploration in solution search space so that will get better solution than other method.

Trajectory attractors of equations of mathematical physics

International Nuclear Information System (INIS)

Vishik, Marko I; Chepyzhov, Vladimir V

2011-01-01

In this survey the method of trajectory dynamical systems and trajectory attractors is described, and is applied in the study of the limiting asymptotic behaviour of solutions of non-linear evolution equations. This method is especially useful in the study of dissipative equations of mathematical physics for which the corresponding Cauchy initial-value problem has a global (weak) solution with respect to the time but the uniqueness of this solution either has not been established or does not hold. An important example of such an equation is the 3D Navier-Stokes system in a bounded domain. In such a situation one cannot use directly the classical scheme of construction of a dynamical system in the phase space of initial conditions of the Cauchy problem of a given equation and find a global attractor of this dynamical system. Nevertheless, for such equations it is possible to construct a trajectory dynamical system and investigate a trajectory attractor of the corresponding translation semigroup. This universal method is applied for various types of equations arising in mathematical physics: for general dissipative reaction-diffusion systems, for the 3D Navier-Stokes system, for dissipative wave equations, for non-linear elliptic equations in cylindrical domains, and for other equations and systems. Special attention is given to using the method of trajectory attractors in approximation and perturbation problems arising in complicated models of mathematical physics. Bibliography: 96 titles.

How students process equations in solving quantitative synthesis problems? Role of mathematical complexity in students’ mathematical performance

Directory of Open Access Journals (Sweden)

Bashirah Ibrahim

2017-10-01

Full Text Available We examine students’ mathematical performance on quantitative “synthesis problems” with varying mathematical complexity. Synthesis problems are tasks comprising multiple concepts typically taught in different chapters. Mathematical performance refers to the formulation, combination, and simplification of equations. Generally speaking, formulation and combination of equations require conceptual reasoning; simplification of equations requires manipulation of equations as computational tools. Mathematical complexity is operationally defined by the number and the type of equations to be manipulated concurrently due to the number of unknowns in each equation. We use two types of synthesis problems, namely, sequential and simultaneous tasks. Sequential synthesis tasks require a chronological application of pertinent concepts, and simultaneous synthesis tasks require a concurrent application of the pertinent concepts. A total of 179 physics major students from a second year mechanics course participated in the study. Data were collected from written tasks and individual interviews. Results show that mathematical complexity negatively influences the students’ mathematical performance on both types of synthesis problems. However, for the sequential synthesis tasks, it interferes only with the students’ simplification of equations. For the simultaneous synthesis tasks, mathematical complexity additionally impedes the students’ formulation and combination of equations. Several reasons may explain this difference, including the students’ different approaches to the two types of synthesis problems, cognitive load, and the variation of mathematical complexity within each synthesis type.

Euler's pioneering equation the most beautiful theorem in mathematics

CERN Document Server

Wilson, Robin

2018-01-01

In 1988 The Mathematical Intelligencer, a quarterly mathematics journal, carried out a poll to find the most beautiful theorem in mathematics. Twenty-four theorems were listed and readers were invited to award each a 'score for beauty'. While there were many worthy competitors, the winner was 'Euler's equation'. In 2004 Physics World carried out a similar poll of 'greatest equations', and found that among physicists Euler's mathematical result came second only to Maxwell's equations. The Stanford mathematician Keith Devlin reflected the feelings of many in describing it as "like a Shakespearian sonnet that captures the very essence of love, or a painting which brings out the beauty of the human form that is far more than just skin deep, Euler's equation reaches down into the very depths of existence."

Multi-band effective mass approximations advanced mathematical models and numerical techniques

CERN Document Server

Koprucki, Thomas

2014-01-01

This book addresses several mathematical models from the most relevant class of kp-Schrödinger systems. Both mathematical models and state-of-the-art numerical methods for adequately solving the arising systems of differential equations are presented. The operational principle of modern semiconductor nano structures, such as quantum wells, quantum wires or quantum dots, relies on quantum mechanical effects. The goal of numerical simulations using quantum mechanical models in the development of semiconductor nano structures is threefold: First they are needed for a deeper understanding of experimental data and of the operational principle. Secondly, they allow us to predict and optimize in advance the qualitative and quantitative properties of new devices in order to minimize the number of prototypes needed. Semiconductor nano structures are embedded as an active region in semiconductor devices. Thirdly and finally, the results of quantum mechanical simulations of semiconductor nano structures can be used wit...

Numerical solution of newton´s cooling differential equation by the methods of euler and runge-kutta

Directory of Open Access Journals (Sweden)

Andresa Pescador

2016-04-01

Full Text Available This article presents the first-order differential equations, which are a very important branch of mathematics as they have a wide applicability, in mathematics, as in physics, biology and economy. The objective of this study was to analyze the resolution of the equation that defines the cooling Newton's law. Verify its behavior using some applications that can be used in the classroom as an auxiliary instrument to the teacher in addressing these contents bringing answers to the questions of the students and motivating them to build their knowledge. It attempted to its resolution through two numerical methods, Euler method and Runge -Kutta method. Finally, there was a comparison of the approach of the solution given by the numerical solution with the analytical resolution whose solution is accurate.

Dispersive traveling wave solutions of the Equal-Width and Modified Equal-Width equations via mathematical methods and its applications

Science.gov (United States)

Lu, Dianchen; Seadawy, Aly R.; Ali, Asghar

2018-06-01

The Equal-Width and Modified Equal-Width equations are used as a model in partial differential equations for the simulation of one-dimensional wave transmission in nonlinear media with dispersion processes. In this article we have employed extend simple equation method and the exp(-varphi(ξ)) expansion method to construct the exact traveling wave solutions of equal width and modified equal width equations. The obtained results are novel and have numerous applications in current areas of research in mathematical physics. It is exposed that our method, with the help of symbolic computation, provides a effective and powerful mathematical tool for solving different kind nonlinear wave problems.

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.

Generalized Bondi-Sachs equations for characteristic formalism of numerical relativity

Science.gov (United States)

Cao, Zhoujian; He, Xiaokai

2013-11-01

The Cauchy formalism of numerical relativity has been successfully applied to simulate various dynamical spacetimes without any symmetry assumption. But discovering how to set a mathematically consistent and physically realistic boundary condition is still an open problem for Cauchy formalism. In addition, the numerical truncation error and finite region ambiguity affect the accuracy of gravitational wave form calculation. As to the finite region ambiguity issue, the characteristic extraction method helps much. But it does not solve all of the above issues. Besides the above problems for Cauchy formalism, the computational efficiency is another problem. Although characteristic formalism of numerical relativity suffers the difficulty from caustics in the inner near zone, it has advantages in relation to all of the issues listed above. Cauchy-characteristic matching (CCM) is a possible way to take advantage of characteristic formalism regarding these issues and treat the inner caustics at the same time. CCM has difficulty treating the gauge difference between the Cauchy part and the characteristic part. We propose generalized Bondi-Sachs equations for characteristic formalism for the Cauchy-characteristic matching end. Our proposal gives out a possible same numerical evolution scheme for both the Cauchy part and the characteristic part. And our generalized Bondi-Sachs equations have one adjustable gauge freedom which can be used to relate the gauge used in the Cauchy part. Then these equations can make the Cauchy part and the characteristic part share a consistent gauge condition. So our proposal gives a possible new starting point for Cauchy-characteristic matching.

Numerical solution of the Black-Scholes equation using cubic spline wavelets

Science.gov (United States)

Černá, Dana

2016-12-01

The Black-Scholes equation is used in financial mathematics for computation of market values of options at a given time. We use the θ-scheme for time discretization and an adaptive scheme based on wavelets for discretization on the given time level. Advantages of the proposed method are small number of degrees of freedom, high-order accuracy with respect to variables representing prices and relatively small number of iterations needed to resolve the problem with a desired accuracy. We use several cubic spline wavelet and multi-wavelet bases and discuss their advantages and disadvantages. We also compare an isotropic and anisotropic approach. Numerical experiments are presented for the two-dimensional Black-Scholes equation.

Numerical Methods for Partial Differential Equations

CERN Document Server

Guo, Ben-yu

1987-01-01

These Proceedings of the first Chinese Conference on Numerical Methods for Partial Differential Equations covers topics such as difference methods, finite element methods, spectral methods, splitting methods, parallel algorithm etc., their theoretical foundation and applications to engineering. Numerical methods both for boundary value problems of elliptic equations and for initial-boundary value problems of evolution equations, such as hyperbolic systems and parabolic equations, are involved. The 16 papers of this volume present recent or new unpublished results and provide a good overview of current research being done in this field in China.

Theoretical and numerical study of the equations of Vlasov-Maxwell in the covariant formalism

International Nuclear Information System (INIS)

Back, A.

2011-11-01

A new point of view is proposed for the simulation of plasmas using the kinetic model which links the equations of Vlasov for the distribution of particles and the equations of Maxwell for the electromagnetic contribution of fields. We use the following principle: the equations of Physics are mathematical objects which put in relation geometrical objects. To preserve the geometrical properties of the various objects in an equation, we use, for the theoretical and numerical study, the differential geometry. All the equations of Physics can be written with differential forms and this point of view is not dependent on the choice of coordinates. We propose then a discretization of the differential forms by using B-Splines. To be coherent with the theory, we also propose a discretization of the various operations of the differential geometry. We test our scheme, first on the equations of Maxwell with several boundary conditions and since it does not depend on the system of coordinates, we also test it when we change coordinates. Finally, we apply the same method to the equations of Vlasov-Poisson in one-dimension and we propose several numerical schemes. (author)

Correlation of Numerical Anxiety and Mathematics Performance

Directory of Open Access Journals (Sweden)

Michael Howard D. Morada

2015-12-01

Full Text Available It has been observed that most students had negative view towards mathematics and as a result, they also performed poorly.As such, it is imperative for every math teacher to understand the reasons behind this negative view to improve their student’s performance. This observation led the researcher to conduct a study on Correlation of Mathematics Performance and Anxiety of third and fourth year students for school year 2012-2013 across the different programs.This study determined the numerical anxiety level and mathematics performance of the respondents along age, gender and programs. The study revealed that students, regardless of age had passing performance. However, female and male students had fair and passing mathematics performance, respectively. Students from College of Business Education, Teacher Education and Computer Studies had fair performance while those from Marine Transportation, Criminal Justice Education and Engineering had passing performance. The study also revealed that students across different variables had moderate numerical anxiety level. Furthermore, it was found out that mathematics performance is significantly related to numerical anxiety. However, the relationship was inverse and small.

An Equation-Type Approach for the Numerical Solution of the Partial Differential Equations Governing Transport Phenomena in Porous Media

KAUST Repository

Sun, Shuyu

2012-06-02

A new technique for the numerical solution of the partial differential equations governing transport phenomena in porous media is introduced. In this technique, the governing equations as depicted from the physics of the problem are used without extra manipulations. In other words, there is no need to reduce the number of governing equations by some sort of mathematical manipulations. This technique enables the separation of the physics part of the problem and the solver part, which makes coding more robust and could be used in several other applications with little or no modifications (e.g., multi-phase flow in porous media). In this method, one abandons the need to construct the coefficient matrix for the pressure equation. Alternatively, the coefficients are automatically generated within the solver routine. We show examples of using this technique to solving several flow problems in porous media.

An Equation-Type Approach for the Numerical Solution of the Partial Differential Equations Governing Transport Phenomena in Porous Media

KAUST Repository

Sun, Shuyu; Salama, Amgad; El-Amin, Mohamed

2012-01-01

A new technique for the numerical solution of the partial differential equations governing transport phenomena in porous media is introduced. In this technique, the governing equations as depicted from the physics of the problem are used without extra manipulations. In other words, there is no need to reduce the number of governing equations by some sort of mathematical manipulations. This technique enables the separation of the physics part of the problem and the solver part, which makes coding more robust and could be used in several other applications with little or no modifications (e.g., multi-phase flow in porous media). In this method, one abandons the need to construct the coefficient matrix for the pressure equation. Alternatively, the coefficients are automatically generated within the solver routine. We show examples of using this technique to solving several flow problems in porous media.

Mathematics and Maxwell's equations

Energy Technology Data Exchange (ETDEWEB)

Boozer, Allen H, E-mail: ahb17@columbia.ed [Department of Applied Physics and Applied Mathematics, Columbia University, New York, NY 10027 (United States)

2010-12-15

The universality of mathematics and Maxwell's equations is not shared by specific plasma models. Computations become more reliable, efficient and transparent if specific plasma models are used to obtain only the information that would otherwise be missing. Constraints of high universality, such as those from mathematics and Maxwell's equations, can be obscured or lost by integrated computations. Recognition of subtle constraints of high universality is important for (1) focusing the design of control systems for magnetic field errors in tokamaks from perturbations that have little effect on the plasma to those that do, (2) clarifying the limits of applicability to astrophysics of computations of magnetic reconnection in fields that have a double periodicity or have B-vector =0 on a surface, as in a Harris sheet. Both require a degree of symmetry not expected in natural systems. Mathematics and Maxwell's equations imply that neighboring magnetic field lines characteristically separate exponentially with distance along a line. This remarkably universal phenomenon has been largely ignored, though it defines a trigger for reconnection through a critical magnitude of exponentiation. These and other examples of the importance of making distinctions and understanding constraints of high universality are explained.

Mathematical and numerical analysis of PN models for photons transport problems

International Nuclear Information System (INIS)

Valentin, Xavier

2015-01-01

Computational costs for direct numerical simulations of photon transport problems are very high in terms of CPU time and memory. One way to tackle this issue is to develop reduced models that a cheaper to solve numerically. There exists number of these models: moments models, discrete ordinates models (S N ), diffusion-like models... In this thesis, we focus on P N models in which the transport operator is approached by mean of a truncated development on the spherical harmonics basis. These models are arbitrary accurate in the angular dimension and are rotationally invariants (in multiple space dimensions). The latter point is fundamental when one wants to simulate inertial confinement fusion (ICF) experiments where the spherical symmetry plays an important part in the accuracy of the numerical solutions. We study the mathematical structure of the PN models and construct a new numerical method in the special case of a one dimensional space dimension with spherical symmetry photon transport problems. We first focus on a linear transport problem in the vacuum. Even in this simple case, it appears in the P N equations geometrical source terms that are stiff in the neighborhood of r = 0 and thus hard to discretize. Existing numerical methods are not satisfactory for multiple reasons: (1) inaccuracy in the neighborhood of r = 0 ('flux-dip'), (2) do not capture steady states (well-balanced scheme), (3) no stability proof. Following recent works, we develop a new well-balanced scheme for which we show the L 2 stability. We then extend the scheme for photon transport problems within a no moving media, the linear Boltzmann equation, and interest ourselves on its behavior in the diffusion limit (asymptotic-preserving property). In a second part, we consider radiation hydrodynamics problems. Since modelization of these problems is still under discussion in the literature, we compare a set of existing models by mean of mathematical analysis and establish a hierarchy

NUMERICAL HOPF BIFURCATION OF DELAY-DIFFERENTIAL EQUATIONS

Institute of Scientific and Technical Information of China (English)

无

2006-01-01

In this paper we consider the numerical solution of some delay differential equations undergoing a Hopf bifurcation. We prove that if the delay differential equations have a Hopf bifurcation point atλ=λ*, then the numerical solution of the equation also has a Hopf bifurcation point atλh =λ* + O(h).

European Conference on Numerical Mathematics and Advanced Applications

CERN Document Server

Manguoğlu, Murat; Tezer-Sezgin, Münevver; Göktepe, Serdar; Uğur, Ömür

2016-01-01

The European Conference on Numerical Mathematics and Advanced Applications (ENUMATH), held every 2 years, provides a forum for discussing recent advances in and aspects of numerical mathematics and scientific and industrial applications. The previous ENUMATH meetings took place in Paris (1995), Heidelberg (1997), Jyvaskyla (1999), Ischia (2001), Prague (2003), Santiago de Compostela (2005), Graz (2007), Uppsala (2009), Leicester (2011) and Lausanne (2013). This book presents a selection of invited and contributed lectures from the ENUMATH 2015 conference, which was organised by the Institute of Applied Mathematics (IAM), Middle East Technical University, Ankara, Turkey, from September 14 to 18, 2015. It offers an overview of central recent developments in numerical analysis, computational mathematics, and applications in the form of contributions by leading experts in the field.

Numerical solution of kinetics equation for point defects accumulation in metals under irradiation

International Nuclear Information System (INIS)

Aldzhambekova, G.T.; Iskakov, B.M.

1999-01-01

In the report the mathematical model, describing processes of generation and accumulation of defects in solids under irradiation is considered. The equations of this model take into account the velocity of Frenkel pairs generation, the mutual recombination of vacancies and the interstitials, as well as velocity of defects absorption by discharge channeling of vacancies and interstitials. By Runge-Kutta method the numerical solution of the model was carried out

Convergence order vs. parallelism in the numerical simulation of the bidomain equations

International Nuclear Information System (INIS)

Sharomi, Oluwaseun; Spiteri, Raymond J

2012-01-01

The propagation of electrical activity in the human heart can be modelled mathematically by the bidomain equations. The bidomain equations represent a multi-scale reaction-diffusion model that consists of a set of ordinary differential equations governing the dynamics at the cellular level coupled with a set of partial differential equations governing the dynamics at the tissue level. Significant computation is generally required to generate clinically useful data from the bidomain equations. Contemporary developments in computer architecture, in particular multi- and many-core computers and graphics processing units, have made such computations feasible. However, the zeal to take advantage to parallel architectures has typically caused another important aspect of numerical methods for the solution of differential equations to be overlooked, namely the convergence order. It is well known that higher-order methods are generally more efficient than lower-order ones when solutions are smooth and relatively high accuracy is desired. In these situations, serial implementations of high-order methods may remain surprisingly competitive with parallel implementations of low-order methods. In this paper, we examine the effect of order on the numerical solution of the bidomain equations in parallel. We find that high-order methods, in particular high-order time-integration methods with relatively better stability properties, tend to outperform their low-order counterparts, even when the latter are run in parallel. In other words, increasing integration order often trumps increasing available computational resources, especially when relatively high accuracy is desired.

Numerical Analysis of Partial Differential Equations

CERN Document Server

Lui, S H

2011-01-01

A balanced guide to the essential techniques for solving elliptic partial differential equations Numerical Analysis of Partial Differential Equations provides a comprehensive, self-contained treatment of the quantitative methods used to solve elliptic partial differential equations (PDEs), with a focus on the efficiency as well as the error of the presented methods. The author utilizes coverage of theoretical PDEs, along with the nu merical solution of linear systems and various examples and exercises, to supply readers with an introduction to the essential concepts in the numerical analysis

Various types of numerical schema for the one-dimensional spherical geometry transport equation

International Nuclear Information System (INIS)

Jaber, Abdelouhab.

1981-07-01

Mathematical and numerical studies of new schemas possessing high accuracy spatial variable properties are described and the corresponding studies presented. In order to do this, the [0,R] x [-1,+1] rectangle is decomposad into Ksub(ij) = [rsub(i),rsub(i+1)] x [μsub(j),μsub(j+1) ] rectangles. Continuous finite element methods employing polynominals of degree 1 in μ and degree 2 in r are defined for each elements. In chapter I, different ways of rendering the particular equation (for μ = -1) discrete are studied. In chapter II, numerical schemas are described and their stability investigated. In chapter III, error estimation theories are exposed and numerical results for different second members, S, given [fr

Mathematical and numerical analysis of a few hydrodynamic and kinetic models of plasma physics

International Nuclear Information System (INIS)

Buet, C.

2005-01-01

My research work deals mainly with the mathematical modelling and the numerical simulation of plasma physics. This document is divided into 3 parts. The first one is a summary of the works done for the numerical solving of collision operators. The common thread of this part is obtaining numerical schemes preserving operators' properties namely physical invariants like mass, momentum and energy, equilibrium states and entropy decrease. These properties are generally checked formally for continuous operators, may give rise to some difficulties for discrete operators. In the second part I present a summary of the works regarding moments methods applied to radiative transfer and the numerical issues dealing with their discretization. The common thread of this part is how to get numerical schemes preserving asymptotic scattering and invariant domains for Lorentz models and also for non-linear telegraph-type equations involved in radiative transfer or electronic plasma. In the third part I present 2 themes linked to collision operators: multi-fluid ionization and the non-existence of linear monotone schemes for some linear parabolic equations

A rapid numerical method for solving Serre-Green-Naghdi equations describing long free surface gravity waves

Science.gov (United States)

Favrie, N.; Gavrilyuk, S.

2017-07-01

A new numerical method for solving the Serre-Green-Naghdi (SGN) equations describing dispersive waves on shallow water is proposed. From the mathematical point of view, the SGN equations are the Euler-Lagrange equations for a ‘master’ lagrangian submitted to a differential constraint which is the mass conservation law. One major numerical challenge in solving the SGN equations is the resolution of an elliptic problem at each time instant. This is the most time-consuming part of the numerical method. The idea is to replace the ‘master’ lagrangian by a one-parameter family of ‘augmented’ lagrangians, depending on a greater number of variables, for which the corresponding Euler-Lagrange equations are hyperbolic. In such an approach, the ‘master’ lagrangian is recovered by the augmented lagrangian in some limit (for example, when the corresponding parameter is large). The choice of such a family of augmented lagrangians is proposed and discussed. The corresponding hyperbolic system is numerically solved by a Godunov type method. Numerical solutions are compared with exact solutions to the SGN equations. It appears that the computational time in solving the hyperbolic system is much lower than in the case where the elliptic operator is inverted. The new method is applied, in particular, to the study of ‘Favre waves’ representing non-stationary undular bores produced after reflection of the fluid flow with a free surface at an immobile wall.

A first course in ordinary differential equations analytical and numerical methods

CERN Document Server

Hermann, Martin

2014-01-01

This book presents a modern introduction to analytical and numerical techniques for solving ordinary differential equations (ODEs). Contrary to the traditional format—the theorem-and-proof format—the book is focusing on analytical and numerical methods. The book supplies a variety of problems and examples, ranging from the elementary to the advanced level, to introduce and study the mathematics of ODEs. The analytical part of the book deals with solution techniques for scalar first-order and second-order linear ODEs, and systems of linear ODEs—with a special focus on the Laplace transform, operator techniques and power series solutions. In the numerical part, theoretical and practical aspects of Runge-Kutta methods for solving initial-value problems and shooting methods for linear two-point boundary-value problems are considered. The book is intended as a primary text for courses on the theory of ODEs and numerical treatment of ODEs for advanced undergraduate and early graduate students. It is assumed t...

Algorithm for research of mathematical physics equations symmetries. Symmetries of the free Schroedinger equation

International Nuclear Information System (INIS)

Kotel'nikov, G.A.

1994-01-01

An algorithm id proposed for research the symmetries of mathematical physics equation. The application of this algorithm to the Schroedinger equation permitted to establish, that in addition to the known symmetry the Schroedinger equation possesses also the relativistic symmetry

Numerical solution of distributed order fractional differential equations

Science.gov (United States)

Katsikadelis, John T.

2014-02-01

In this paper a method for the numerical solution of distributed order FDEs (fractional differential equations) of a general form is presented. The method applies to both linear and nonlinear equations. The Caputo type fractional derivative is employed. The distributed order FDE is approximated with a multi-term FDE, which is then solved by adjusting appropriately the numerical method developed for multi-term FDEs by Katsikadelis. Several example equations are solved and the response of mechanical systems described by such equations is studied. The convergence and the accuracy of the method for linear and nonlinear equations are demonstrated through well corroborated numerical results.

A general theory for the construction of best-fit correlation equations for multi-dimensioned numerical data

International Nuclear Information System (INIS)

Moore, S.E.; Moffat, D.G.

2007-01-01

A general theory for the construction of best-fit correlation equations for multi-dimensioned sets of numerical data is presented. This new theory is based on the mathematics of n-dimensional surfaces and goodness-of-fit statistics. It is shown that orthogonal best-fit analytical trend lines for each of the independent parameters of the data can be used to construct an overall best-fit correlation equation that satisfies both physical boundary conditions and best-of-fit statistical measurements. Application of the theory is illustrated by fitting a three-parameter set of numerical finite-element maximum-stress data obtained earlier by Dr. Moffat for pressure vessel nozzles and/or piping system branch connections

Integrating spatial and numerical structure in mathematical patterning

Science.gov (United States)

Ni’mah, K.; Purwanto; Irawan, E. B.; Hidayanto, E.

2018-03-01

This paper reports a study monitoring the integrating spatial and numerical structure in mathematical patterning skills of 30 students grade 7th of junior high school. The purpose of this research is to clarify the processes by which learners construct new knowledge in mathematical patterning. Findings indicate that: (1) students are unable to organize the structure of spatial and numerical, (2) students were only able to organize the spatial structure, but the numerical structure is still incorrect, (3) students were only able to organize numerical structure, but its spatial structure is still incorrect, (4) students were able to organize both of the spatial and numerical structure.

Mathematical and Numerical Modeling in Maritime Geomechanics

Directory of Open Access Journals (Sweden)

Miguel Martín Stickle

2012-04-01

Full Text Available A theoretical and numerical framework to model the foundation of marine offshore structures is presented. The theoretical model is composed by a system of partial differential equations describing coupling between seabed solid skeleton and pore fluids (water, air, oil,... combined with a system of ordinary differential equations describing the specific constitutive relation of the seabed soil skeleton. Once the theoretical model is described, the finite element numerical procedure to achieve an approximate solution of the overning equations is outlined. In order to validate the proposed theoretical and numerical framework the seaward tilt mechanism induced by the action of breaking waves over a vertical breakwater is numerically reproduced. The results numerically attained are in agreement with the main conclusions drawn from the literature associated with this failure mechanism.

Numerical solving of equations in the work of José Mariano Vallejo

Science.gov (United States)

Pacheco Castelao, José-Miguel; Pérez-Fern; ández, F. Javier; Suárez Alemán, Carlos-Oswaldo

2007-09-01

The progress of Mathematics during the nineteenth century was characterised both by an enormous acquisition of new knowledge and by the attempts to introduce rigour in reasoning patterns and mathematical writing. Cauchy's presentation of Mathematical Analysis was not immediately accepted, and many writers, though aware of that new style, did not use it in their own mathematical production. This paper is devoted to an episode of this sort that took place in Spain during the first half of the century: It deals with the presentation of a method for numerically solving algebraic equations by José Mariano Vallejo, a late Spanish follower of the Enlightenment ideas, politician, writer, and mathematician who published it in the fourth (1840) edition of his book Compendio de Mathemáticas Puras y Mistas, claiming to have discovered it on his own. Vallejo's main achievement was to write down the whole procedure in a very careful way taking into account the different types of roots, although he paid little attention to questions such as convergence checks and the fulfilment of the hypotheses of Rolle's Theorem. For sure this lack of mathematical care prevented Vallejo to occupy a place among the forerunners of Computational Algebra.

Stochastic optimal control, forward-backward stochastic differential equations and the Schroedinger equation

Energy Technology Data Exchange (ETDEWEB)

Paul, Wolfgang; Koeppe, Jeanette [Institut fuer Physik, Martin Luther Universitaet, 06099 Halle (Germany); Grecksch, Wilfried [Institut fuer Mathematik, Martin Luther Universitaet, 06099 Halle (Germany)

2016-07-01

The standard approach to solve a non-relativistic quantum problem is through analytical or numerical solution of the Schroedinger equation. We show a way to go around it. This way is based on the derivation of the Schroedinger equation from conservative diffusion processes and the establishment of (several) stochastic variational principles leading to the Schroedinger equation under the assumption of a kinematics described by Nelson's diffusion processes. Mathematically, the variational principle can be considered as a stochastic optimal control problem linked to the forward-backward stochastic differential equations of Nelson's stochastic mechanics. The Hamilton-Jacobi-Bellmann equation of this control problem is the Schroedinger equation. We present the mathematical background and how to turn it into a numerical scheme for analyzing a quantum system without using the Schroedinger equation and exemplify the approach for a simple 1d problem.

Mathematical and numerical aspects of low mach number flows

Energy Technology Data Exchange (ETDEWEB)

Schochet, St.; Bresch, D.; Grenier, E.; Alazard, T.; Gordner, A.; Sankaran, V.; Massot, M.; Sery, R.; Pebay, P.; Lunch, O.; Mazhorova, O.; Turkel, O.E.; Faille, I.; Danchin, R.; Allain, O.; Birken, P.; Lafitte, O.; Kloczko, T.; Frick, W.; Bui, T.; Dellacherie, S.; Klein, R.; Roe, Ph.; Accary, G.; Braack, M.; Picano, F.; Cadiou, A.; Dinescu, C.; Lesage, A.C.; Wesseling, P.; Heuveline, V.; Jobelin, M.; Weisman, C.; Merkle, C.

2004-07-01

Low Mach number flows represent a significant part of the various flows encountered in geophysics, industry or every day life. Paradoxically, the mathematical analysis of the equations governing these flows is difficult and on the practical side, the research of numerical algorithms valid for all flow speeds is continuing to be a challenge. However, in the last decade, both from the theoretical and the numerical sides, significant progresses were made in the understanding and analysis of the equations governing these flows. This conference intends to provide an up-to-date inventory of recent mathematical and numerical results in the analysis of these flows by bringing together both mathematicians and numericists active in this area. In the framework of the conference, a numerical workshop is organized which proposes to compute several challenging low Mach number flows: liquid flow around non-cavitating and cavitating NACA0015 hydrofoil, natural convection with large temperature differences, free convection, free surface flow, vessel pressurization. This document brings together the descriptions of the test cases of the numerical workshop and the abstracts of the conference papers: A 3D high order finite volume method for the prediction of near-critical fluid flows (G. ACCARY, I. RASPO, P. BONTOUX, B. ZAPPOLI); low Mach number limit of the non-isentropic Navier-Stokes equations (T. ALAZARD); simulation of cavitation rolls past a forward step with a bubble model (O. ALLAIN, N. BLASKA, C. LECA); flux preconditioning methods and fire events (P. BIRKEN, A. MEISTER); an adaptive finite element solver for compressible flows: application to heat-driven cavity benchmarks in 2D and 3D (M. BRAACK); comparison of various implicit, explicit, centered and upwind schemes for the simulation of compressed flows on moving mesh (A. CADIOU, M. BUFFAT, L. Le PENVEN, C. Le RIBAULT); low Mach number limit for viscous compressible flows (R. DANCHIN); some Properties of the low Mach number

Dimensional analysis yields the general second-order differential equation underlying many natural phenomena: the mathematical properties of a phenomenon's data plot then specify a unique differential equation for it.

Science.gov (United States)

Kepner, Gordon R

2014-08-27

This study uses dimensional analysis to derive the general second-order differential equation that underlies numerous physical and natural phenomena described by common mathematical functions. It eschews assumptions about empirical constants and mechanisms. It relies only on the data plot's mathematical properties to provide the conditions and constraints needed to specify a second-order differential equation that is free of empirical constants for each phenomenon. A practical example of each function is analyzed using the general form of the underlying differential equation and the observable unique mathematical properties of each data plot, including boundary conditions. This yields a differential equation that describes the relationship among the physical variables governing the phenomenon's behavior. Complex phenomena such as the Standard Normal Distribution, the Logistic Growth Function, and Hill Ligand binding, which are characterized by data plots of distinctly different sigmoidal character, are readily analyzed by this approach. It provides an alternative, simple, unifying basis for analyzing each of these varied phenomena from a common perspective that ties them together and offers new insights into the appropriate empirical constants for describing each phenomenon.

Mathematical and numerical foundations of turbulence models and applications

CERN Document Server

Chacón Rebollo, Tomás

2014-01-01

With applications to climate, technology, and industry, the modeling and numerical simulation of turbulent flows are rich with history and modern relevance. The complexity of the problems that arise in the study of turbulence requires tools from various scientific disciplines, including mathematics, physics, engineering, and computer science. Authored by two experts in the area with a long history of collaboration, this monograph provides a current, detailed look at several turbulence models from both the theoretical and numerical perspectives. The k-epsilon, large-eddy simulation, and other models are rigorously derived and their performance is analyzed using benchmark simulations for real-world turbulent flows. Mathematical and Numerical Foundations of Turbulence Models and Applications is an ideal reference for students in applied mathematics and engineering, as well as researchers in mathematical and numerical fluid dynamics. It is also a valuable resource for advanced graduate students in fluid dynamics,...

Mathematical Methods for Engineers and Scientists 2 Vector Analysis, Ordinary Differential Equations and Laplace Transforms

CERN Document Server

Tang, Kwong-Tin

2007-01-01

Pedagogical insights gained through 30 years of teaching applied mathematics led the author to write this set of student-oriented books. Topics such as complex analysis, matrix theory, vector and tensor analysis, Fourier analysis, integral transforms, ordinary and partial differential equations are presented in a discursive style that is readable and easy to follow. Numerous clearly stated, completely worked out examples together with carefully selected problem sets with answers are used to enhance students' understanding and manipulative skill. The goal is to make students comfortable and confident in using advanced mathematical tools in junior, senior, and beginning graduate courses.

Mathematical Methods for Engineers and Scientists 3 Fourier Analysis, Partial Differential Equations and Variational Methods

CERN Document Server

Tang, Kwong-Tin

2007-01-01

Pedagogical insights gained through 30 years of teaching applied mathematics led the author to write this set of student oriented books. Topics such as complex analysis, matrix theory, vector and tensor analysis, Fourier analysis, integral transforms, ordinary and partial differential equations are presented in a discursive style that is readable and easy to follow. Numerous clearly stated, completely worked out examples together with carefully selected problem sets with answers are used to enhance students' understanding and manipulative skill. The goal is to make students comfortable and confident in using advanced mathematical tools in junior, senior, and beginning graduate courses.

Equations in mathematical physics a practical course

CERN Document Server

Pikulin, Victor P

2001-01-01

This handbook is addressed to students of technology institutf's where a course on mathematical physics of relatively reduced volume is offered, as well as to engineers and scientists. The aim of the handbook is to treat (demonstrate) the basic methods for solving the simplest problems of classical mathematical physics. The most basic among the methods considered hrre i8 the superposition method. It allows one, based on particular linearly indepmdent HolutionH (solution "atoms"), to obtain the solution of a given problem. To that end the "Hupply" of solution atoms must be complete. This method is a development of the well-known method of particular solutions from the theory of ordinar~' differelltial equations. In contrast to the case of ordinary differential equations, where the number of linearly independent 80lutions is always finite, for a linear partial differrntial equation a complete "supply" of solution atoms is always infinite. This infinite set of Holutions may be discrete (for example, for regular ...

Numerical Solution of Heun Equation Via Linear Stochastic Differential Equation

Directory of Open Access Journals (Sweden)

Hamidreza Rezazadeh

2014-05-01

Full Text Available In this paper, we intend to solve special kind of ordinary differential equations which is called Heun equations, by converting to a corresponding stochastic differential equation(S.D.E.. So, we construct a stochastic linear equation system from this equation which its solution is based on computing fundamental matrix of this system and then, this S.D.E. is solved by numerically methods. Moreover, its asymptotic stability and statistical concepts like expectation and variance of solutions are discussed. Finally, the attained solutions of these S.D.E.s compared with exact solution of corresponding differential equations.

An introduction to differential equations using MATLAB

CERN Document Server

Butt, Rizwan

2016-01-01

An Introduction to Differential Equations using MATLAB exploits the symbolic, numerical, and graphical capabilitiesof MATLAB to develop a thorough understanding of differential equations algorithms. This book provides the readerwith numerous applications, m-files, and practical examples to problems. Balancing theoretical concepts withcomputational speed and accuracy, the book includes numerous short programs in MATLAB that can be used to solveproblems involving first-and higher-order differential equations, Laplace transforms, linear systems of differentialequations, numerical solutions of differential equations, computer graphics, and more. The author emphasizes thebasic ideas of analytical and numerical techniques and the uses of modern mathematical software (MATLAB) ratherthan relying only on complex mathematical derivations to engineers, mathematician, computer scientists, andphysicists or for use as a textbook in applied or computational courses.A CD-ROM with all the figures, codes, solutions, appendices...

On the Existence and the Applications of Modified Equations for Stochastic Differential Equations

KAUST Repository

Zygalakis, K. C.

2011-01-01

In this paper we describe a general framework for deriving modified equations for stochastic differential equations (SDEs) with respect to weak convergence. Modified equations are derived for a variety of numerical methods, such as the Euler or the Milstein method. Existence of higher order modified equations is also discussed. In the case of linear SDEs, using the Gaussianity of the underlying solutions, we derive an SDE which the numerical method solves exactly in the weak sense. Applications of modified equations in the numerical study of Langevin equations is also discussed. © 2011 Society for Industrial and Applied Mathematics.

Numerical Treatment of the Boltzmann Equation for Self-Propelled Particle Systems

Directory of Open Access Journals (Sweden)

Florian Thüroff

2014-11-01

Full Text Available Kinetic theories constitute one of the most promising tools to decipher the characteristic spatiotemporal dynamics in systems of actively propelled particles. In this context, the Boltzmann equation plays a pivotal role, since it provides a natural translation between a particle-level description of the system’s dynamics and the corresponding hydrodynamic fields. Yet, the intricate mathematical structure of the Boltzmann equation substantially limits the progress toward a full understanding of this equation by solely analytical means. Here, we propose a general framework to numerically solve the Boltzmann equation for self-propelled particle systems in two spatial dimensions and with arbitrary boundary conditions. We discuss potential applications of this numerical framework to active matter systems and use the algorithm to give a detailed analysis to a model system of self-propelled particles with polar interactions. In accordance with previous studies, we find that spatially homogeneous isotropic and broken-symmetry states populate two distinct regions in parameter space, which are separated by a narrow region of spatially inhomogeneous, density-segregated moving patterns. We find clear evidence that these three regions in parameter space are connected by first-order phase transitions and that the transition between the spatially homogeneous isotropic and polar ordered phases bears striking similarities to liquid-gas phase transitions in equilibrium systems. Within the density-segregated parameter regime, we find a novel stable limit-cycle solution of the Boltzmann equation, which consists of parallel lanes of polar clusters moving in opposite directions, so as to render the overall symmetry of the system’s ordered state nematic, despite purely polar interactions on the level of single particles.

Numerical solutions of the Vlasov equation

International Nuclear Information System (INIS)

Satofuka, Nobuyuki; Morinishi, Koji; Nishida, Hidetoshi

1985-01-01

A numerical procedure is derived for the solutions of the one- and two-dimensional Vlasov-Poisson system equations. This numerical procedure consists of the phase space discretization and the integration of the resulting set of ordinary differential equations. In the phase space discretization, derivatives with respect to the phase space variable are approximated by a weighted sum of the values of the distribution function at properly chosen neighboring points. Then, the resulting set of ordinary differential equations is solved by using an appropriate time integration scheme. The results for linear Landau damping, nonlinear Landau damping and counter-streaming plasmas are investigated and compared with those of the splitting scheme. The proposed method is found to be very accurate and efficient. (author)

Numerical solution of stiff burnup equation with short half lived nuclides by the Krylov subspace method

International Nuclear Information System (INIS)

Yamamoto, Akio; Tatsumi, Masahiro; Sugimura, Naoki

2007-01-01

The Krylov subspace method is applied to solve nuclide burnup equations used for lattice physics calculations. The Krylov method is an efficient approach for solving ordinary differential equations with stiff nature such as the nuclide burnup with short lived nuclides. Some mathematical fundamentals of the Krylov subspace method and its application to burnup equations are discussed. Verification calculations are carried out in a PWR pin-cell geometry with UO 2 fuel. A detailed burnup chain that includes 193 fission products and 28 heavy nuclides is used in the verification calculations. Shortest half life found in the present burnup chain is approximately 30 s ( 106 Rh). Therefore, conventional methods (e.g., the Taylor series expansion with scaling and squaring) tend to require longer computation time due to numerical stiffness. Comparison with other numerical methods (e.g., the 4-th order Runge-Kutta-Gill) reveals that the Krylov subspace method can provide accurate solution for a detailed burnup chain used in the present study with short computation time. (author)

Numerical solution of neutron transport equations in discrete ordinates and slab geometry

International Nuclear Information System (INIS)

Serrano Pedraza, F.

1985-01-01

An unified formalism to solve numerically, between other equation, the neutron transport in discrete ordinates, slab geometry, several energy groups and independents of time, has been developed recently. Such a formalism cover some of the conventional schemes as diamond difference, (WDD) characteristic step (SC) lineal characteristic (LC), quadratic characteristic (QC) and lineal discontinuous. Unified formation gives before hand the convergence order of the previously selected scheme. In fact it allows besides to generate a big amount of numerical schemes, with which is also possible to solve numerical equations as soon as neutron transport. The essential purpose of this work was to solve the neutron transport equations in slab geometry and discrete ordinates considering several energy groups without to take under advisement time dependence based in the above mentioned unified formalism. To reach this purpose it was necesary to design a computer code with the name TNOD1 (Neutron transport in discrete ordinates and 1 dimension) which includes each one of the schemes already pointed out. there exist two numerical schemes, also recently developed, quadratic continuous (QC) and cubic continuous (CN), although covered by unified formalism, it has been possible to include them inside this computer code without make substantial changes in its structure. In chapter I, derivative of neutron transport equation independent of time is taken, for angular flux, including boundary conditions and discontinuity. In chapter II the neutron transport equations are obtained in multigroups, independents of time, for approximation of discrete ordinates. Description of theory related with unified formalism and its relationship with mentioned discretization schemes is presented in chapter III. Chapter IV describes the computer code developed and finally, in chapter V different numerical results obtained with TNOD1 program are shown. In Appendix A theorems and mathematical arguments used

A Numerical Comparison of Soave Redlich Kwong and Peng-Robinson Equations of State for Predicting Hydrocarbons’ Thermodynamic Properties

Directory of Open Access Journals (Sweden)

B. Hussain

2018-02-01

Full Text Available Mixture phase equilibrium and thermodynamic properties have a significant role in industry. Numerical analysis of flash calculation generates an appropriate solution for the problem. In this research, a comparison of Soave Redlich Kwong (SRK and Peng-Robinson (PR equations of state predicting the thermodynamic properties of a mixture of hydrocarbon and related compounds in a critical region at phase equilibrium is performed. By applying mathematical modeling of both equations of states, the behavior of binary gases mixtures is monitored. The numerical analysis of isothermal flash calculations is applied to study the pressure behavior with volume and mole fraction. The approach used in this research shows considerable convergence with experimental results available in the literature.

Numerical solutions of diffusive logistic equation

International Nuclear Information System (INIS)

Afrouzi, G.A.; Khademloo, S.

2007-01-01

In this paper we investigate numerically positive solutions of a superlinear Elliptic equation on bounded domains. The study of Diffusive logistic equation continues to be an active field of research. The subject has important applications to population migration as well as many other branches of science and engineering. In this paper the 'finite difference scheme' will be developed and compared for solving the one- and three-dimensional Diffusive logistic equation. The basis of the analysis of the finite difference equations considered here is the modified equivalent partial differential equation approach, developed from many authors these years

Distribution of the Discretization and Algebraic Error in Numerical Solution of Partial Differential Equations

Czech Academy of Sciences Publication Activity Database

Papež, Jan; Liesen, J.; Strakoš, Z.

2014-01-01

Roč. 449, 15 May (2014), s. 89-114 ISSN 0024-3795 R&D Projects: GA AV ČR IAA100300802; GA ČR GA201/09/0917 Grant - others:GA MŠk(CZ) LL1202; GA UK(CZ) 695612 Institutional support: RVO:67985807 Keywords : numerical solution of partial differential equations * finite element method * adaptivity * a posteriori error analysis * discretization error * algebra ic error * spatial distribution of the error Subject RIV: BA - General Mathematics Impact factor: 0.939, year: 2014

Mathematical and numerical analysis of hyper-elastic systems and introduction of plasticity; Analyse mathematique et numerique de systemes hyperelastiques et introduction de la plasticite

Energy Technology Data Exchange (ETDEWEB)

Kluth, G

2008-12-15

The goal is to model mathematically and numerically the dynamic phenomenons for solids in finite plasticity. We suggest a model that we call hyper-elasto-plastic based on hyper-elastic systems of conservation laws and on the use of an equation of state that we have constructed so as to achieve the plastic yield criterion of Von Mises. This model gives exact (analytic) solutions with shock split to flyer-plate experiments. The mathematical analysis of this model is done (hyperbolicity, characteristic fields, involutions and entropy). In the numerical part, we give 1D and 2D Lagrangian schemes which satisfy an entropy criterion. Moreover, thanks to a special discretization of the equations on deformation gradient, we satisfy some discrete involutions. In this work, the degeneracy of the solid model into hydrodynamic models is studied at the continuous level, and achieved at the numerical one. On different problems, we show the validity of our model and our numerical schemes. (author)

The essential harmony in the classical equations of mathematical physics

CERN Document Server

Nucci, M C

2001-01-01

The possibility to transform any system of linear ordinary differential equations into a system of constant coefficient equations is demonstrated using Lie theory. Some examples relate the classical equations of mathematical physics to the simple harmonic oscillator. The roles of the third order form of the Ermakov-Pinney equation and of Fleischen-von Weltunter systems are explained.

Analysis of numerical solutions for Bateman equations

International Nuclear Information System (INIS)

Loch, Guilherme G.; Bevilacqua, Joyce S.

2013-01-01

The implementation of stable and efficient numerical methods for solving problems involving nuclear transmutation and radioactive decay chains is the main scope of this work. The physical processes associated with irradiations of samples in particle accelerators, or the burning spent nuclear fuel in reactors, or simply the natural decay chains, can be represented by a set of first order ordinary differential equations with constant coefficients, for instance, the decay radioactive constants of each nuclide in the chain. Bateman proposed an analytical solution for a particular case of a linear chain with n nuclides decaying in series and with different decay constants. For more complex and realistic applications, the construction of analytical solutions is not viable and the introduction of numerical techniques is imperative. However, depending on the magnitudes of the decay radioactive constants, the matrix of coefficients could be almost singular, generating unstable and non convergent numerical solutions. In this work, different numerical strategies for solving systems of differential equations were implemented, the Runge-Kutta 4-4, Adams Predictor-Corrector (PC2) and the Rosenbrock algorithm, this last one more specific for stiff equations. Consistency, convergence and stability of the numerical solutions are studied and the performance of the methods is analyzed for the case of the natural decay chain of Uranium-235 comparing numerical with analytical solutions. (author)

New Numerical Treatment for Solving the KDV Equation

Directory of Open Access Journals (Sweden)

khalid ali

2017-01-01

Full Text Available In the present article, a numerical method is proposed for the numerical solution of the KdV equation by using collocation method with the modified exponential cubic B-spline. In this paper we convert the KdV equation to system of two equations. The method is shown to be unconditionally stable using von-Neumann technique. To test accuracy the error norms2L, ?L are computed. Three invariants of motion are predestined to determine the preservation properties of the problem, and the numerical scheme leads to careful and active results. Furthermore, interaction of two and three solitary waves is shown. These results show that the technique introduced here is easy to apply.

Intuitive physics knowledge, physics problem solving and the role of mathematical equations

Directory of Open Access Journals (Sweden)

Laura Buteler

2012-09-01

Full Text Available The present work explores the role that mathematical equations play in modifying students’ physical intuition (diSessa, 1993. The work is carried out assuming that students achieve a great deal of the refinement in their physical intuitions during problem solving (Sherin, 2006. The study is guided by the question of how the use of mathematical equations contributes to this refinement. The authors aim at expanding on Sherin´s (2006 hypothesis, suggesting a more bounding relation between physical intuitions and mathematics. In this scenario, intuitions play a more compelling role in “deciding” which equations are acceptable and which are not. Our hypothesis is constructed on the basis of three cases: the first published by Sherin (2006 and two more from registries of our own. The three cases are compared and analyzed in relation to the role of mathematical equations in refining – or not – the intuitive knowledge students bring to play during problem solving.

Mathematical modeling and the two-phase constitutive equations

International Nuclear Information System (INIS)

Boure, J.A.

1975-01-01

The problems raised by the mathematical modeling of two-phase flows are summarized. The models include several kinds of equations, which cannot be discussed independently, such as the balance equations and the constitutive equations. A review of the various two-phase one-dimensional models proposed to date, and of the constitutive equations they imply, is made. These models are either mixture models or two-fluid models. Due to their potentialities, the two-fluid models are discussed in more detail. To avoid contradictions, the form of the constitutive equations involved in two-fluid models must be sufficiently general. A special form of the two-fluid models, which has particular advantages, is proposed. It involves three mixture balance equations, three balance equations for slip and thermal non-equilibriums, and the necessary constitutive equations [fr

The Influence of Pre-University Students' Mathematics Test Anxiety and Numerical Anxiety on Mathematics Achievement

Science.gov (United States)

Seng, Ernest Lim Kok

2015-01-01

This study examines the relationship between mathematics test anxiety and numerical anxiety on students' mathematics achievement. 140 pre-university students who studied at one of the institutes of higher learning were being investigated. Gender issue pertaining to mathematics anxieties was being addressed besides investigating the magnitude of…

Spurious Numerical Solutions Of Differential Equations

Science.gov (United States)

Lafon, A.; Yee, H. C.

1995-01-01

Paper presents detailed study of spurious steady-state numerical solutions of differential equations that contain nonlinear source terms. Main objectives of this study are (1) to investigate how well numerical steady-state solutions of model nonlinear reaction/convection boundary-value problem mimic true steady-state solutions and (2) to relate findings of this investigation to implications for interpretation of numerical results from computational-fluid-dynamics algorithms and computer codes used to simulate reacting flows.

Numerical resolution of Navier-Stokes equations coupled to the heat equation

International Nuclear Information System (INIS)

Zenouda, Jean-Claude

1970-08-01

The author proves a uniqueness theorem for the time dependent Navier-Stokes equations coupled with heat flow in the two-dimensional case. He studies stability and convergence of several finite - difference schemes to solve these equations. Numerical experiments are done in the case of a square domain. (author) [fr

Intentional and Automatic Numerical Processing as Predictors of Mathematical Abilities in Primary School Children

Directory of Open Access Journals (Sweden)

Violeta ePina

2015-03-01

Full Text Available Previous studies have suggested that numerical processing relates to mathematical performance, but it seems that such relationship is more evident for intentional than for automatic numerical processing. In the present study we assessed the relationship between the two types of numerical processing and specific mathematical abilities in a sample of 109 children in grades 1 to 6. Participants were tested in an ample range of mathematical tests and also performed both a numerical and a size comparison task. The results showed that numerical processing related to mathematical performance only when inhibitory control was involved in the comparison tasks. Concretely, we found that intentional numerical processing, as indexed by the numerical distance effect in the numerical comparison task, was related to mathematical reasoning skills only when the task-irrelevant dimension (the physical size was incongruent; whereas automatic numerical processing, indexed by the congruency effect in the size comparison task, was related to mathematical calculation skills only when digits were separated by small distance. The observed double dissociation highlights the relevance of both intentional and automatic numerical processing in mathematical skills, but when inhibitory control is also involved.

Numerical analysis for multi-group neutron-diffusion equation using Radial Point Interpolation Method (RPIM)

International Nuclear Information System (INIS)

Kim, Kyung-O; Jeong, Hae Sun; Jo, Daeseong

2017-01-01

Highlights: • Employing the Radial Point Interpolation Method (RPIM) in numerical analysis of multi-group neutron-diffusion equation. • Establishing mathematical formation of modified multi-group neutron-diffusion equation by RPIM. • Performing the numerical analysis for 2D critical problem. - Abstract: A mesh-free method is introduced to overcome the drawbacks (e.g., mesh generation and connectivity definition between the meshes) of mesh-based (nodal) methods such as the finite-element method and finite-difference method. In particular, the Point Interpolation Method (PIM) using a radial basis function is employed in the numerical analysis for the multi-group neutron-diffusion equation. The benchmark calculations are performed for the 2D homogeneous and heterogeneous problems, and the Multiquadrics (MQ) and Gaussian (EXP) functions are employed to analyze the effect of the radial basis function on the numerical solution. Additionally, the effect of the dimensionless shape parameter in those functions on the calculation accuracy is evaluated. According to the results, the radial PIM (RPIM) can provide a highly accurate solution for the multiplication eigenvalue and the neutron flux distribution, and the numerical solution with the MQ radial basis function exhibits the stable accuracy with respect to the reference solutions compared with the other solution. The dimensionless shape parameter directly affects the calculation accuracy and computing time. Values between 1.87 and 3.0 for the benchmark problems considered in this study lead to the most accurate solution. The difference between the analytical and numerical results for the neutron flux is significantly increased in the edge of the problem geometry, even though the maximum difference is lower than 4%. This phenomenon seems to arise from the derivative boundary condition at (x,0) and (0,y) positions, and it may be necessary to introduce additional strategy (e.g., the method using fictitious points and

Advanced Topics in Computational Partial Differential Equations: Numerical Methods and Diffpack Programming

International Nuclear Information System (INIS)

Katsaounis, T D

2005-01-01

The scope of this book is to present well known simple and advanced numerical methods for solving partial differential equations (PDEs) and how to implement these methods using the programming environment of the software package Diffpack. A basic background in PDEs and numerical methods is required by the potential reader. Further, a basic knowledge of the finite element method and its implementation in one and two space dimensions is required. The authors claim that no prior knowledge of the package Diffpack is required, which is true, but the reader should be at least familiar with an object oriented programming language like C++ in order to better comprehend the programming environment of Diffpack. Certainly, a prior knowledge or usage of Diffpack would be a great advantage to the reader. The book consists of 15 chapters, each one written by one or more authors. Each chapter is basically divided into two parts: the first part is about mathematical models described by PDEs and numerical methods to solve these models and the second part describes how to implement the numerical methods using the programming environment of Diffpack. Each chapter closes with a list of references on its subject. The first nine chapters cover well known numerical methods for solving the basic types of PDEs. Further, programming techniques on the serial as well as on the parallel implementation of numerical methods are also included in these chapters. The last five chapters are dedicated to applications, modelled by PDEs, in a variety of fields. In summary, the book focuses on the computational and implementational issues involved in solving partial differential equations. The potential reader should have a basic knowledge of PDEs and the finite difference and finite element methods. The examples presented are solved within the programming framework of Diffpack and the reader should have prior experience with the particular software in order to take full advantage of the book. Overall

The universe in zero words the story of mathematics as told through equations

CERN Document Server

Mackenzie, Dana

2012-01-01

Most popular books about science, and even about mathematics, tiptoe around equations as if they were something to be hidden from the reader's tender eyes. Dana Mackenzie starts from the opposite premise: He celebrates equations. No history of art would be complete without pictures. Why, then, should a history of mathematics--the universal language of science--keep the masterpieces of the subject hidden behind a veil? The Universe in Zero Words tells the history of twenty-four great and beautiful equations that have shaped mathematics, science, and society--from the elementary

Numerical integration of asymptotic solutions of ordinary differential equations

Science.gov (United States)

Thurston, Gaylen A.

1989-01-01

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

Numerical analysis

CERN Document Server

Rao, G Shanker

2006-01-01

About the Book: This book provides an introduction to Numerical Analysis for the students of Mathematics and Engineering. The book is designed in accordance with the common core syllabus of Numerical Analysis of Universities of Andhra Pradesh and also the syllabus prescribed in most of the Indian Universities. Salient features: Approximate and Numerical Solutions of Algebraic and Transcendental Equation Interpolation of Functions Numerical Differentiation and Integration and Numerical Solution of Ordinary Differential Equations The last three chapters deal with Curve Fitting, Eigen Values and Eigen Vectors of a Matrix and Regression Analysis. Each chapter is supplemented with a number of worked-out examples as well as number of problems to be solved by the students. This would help in the better understanding of the subject. Contents: Errors Solution of Algebraic and Transcendental Equations Finite Differences Interpolation with Equal Intervals Interpolation with Unequal Int...

Numerical methods of mathematical optimization with Algol and Fortran programs

CERN Document Server

Künzi, Hans P; Zehnder, C A; Rheinboldt, Werner

1971-01-01

Numerical Methods of Mathematical Optimization: With ALGOL and FORTRAN Programs reviews the theory and the practical application of the numerical methods of mathematical optimization. An ALGOL and a FORTRAN program was developed for each one of the algorithms described in the theoretical section. This should result in easy access to the application of the different optimization methods.Comprised of four chapters, this volume begins with a discussion on the theory of linear and nonlinear optimization, with the main stress on an easily understood, mathematically precise presentation. In addition

Numerical methods and analysis of multiscale problems

CERN Document Server

Madureira, Alexandre L

2017-01-01

This book is about numerical modeling of multiscale problems, and introduces several asymptotic analysis and numerical techniques which are necessary for a proper approximation of equations that depend on different physical scales. Aimed at advanced undergraduate and graduate students in mathematics, engineering and physics – or researchers seeking a no-nonsense approach –, it discusses examples in their simplest possible settings, removing mathematical hurdles that might hinder a clear understanding of the methods. The problems considered are given by singular perturbed reaction advection diffusion equations in one and two-dimensional domains, partial differential equations in domains with rough boundaries, and equations with oscillatory coefficients. This work shows how asymptotic analysis can be used to develop and analyze models and numerical methods that are robust and work well for a wide range of parameters.

Numerical methods and analysis of the nonlinear Vlasov equation on unstructured meshes of phase space

International Nuclear Information System (INIS)

Besse, Nicolas

2003-01-01

This work is dedicated to the mathematical and numerical studies of the Vlasov equation on phase-space unstructured meshes. In the first part, new semi-Lagrangian methods are developed to solve the Vlasov equation on unstructured meshes of phase space. As the Vlasov equation describes multi-scale phenomena, we also propose original methods based on a wavelet multi-resolution analysis. The resulting algorithm leads to an adaptive mesh-refinement strategy. The new massively-parallel computers allow to use these methods with several phase-space dimensions. Particularly, these numerical schemes are applied to plasma physics and charged particle beams in the case of two-, three-, and four-dimensional Vlasov-Poisson systems. In the second part we prove the convergence and give error estimates for several numerical schemes applied to the Vlasov-Poisson system when strong and classical solutions are considered. First we show the convergence of a semi-Lagrangian scheme on an unstructured mesh of phase space, when the regularity hypotheses for the initial data are minimal. Then we demonstrate the convergence of classes of high-order semi-Lagrangian schemes in the framework of the regular classical solution. In order to reconstruct the distribution function, we consider symmetrical Lagrange polynomials, B-Splines and wavelets bases. Finally we prove the convergence of a semi-Lagrangian scheme with propagation of gradients yielding a high-order and stable reconstruction of the solution. (author) [fr

Solution of the mathematical adjoint equations for an interface current nodal formulation

International Nuclear Information System (INIS)

Yang, W.S.; Taiwo, T.A.; Khalil, H.

1994-01-01

Two techniques for solving the mathematical adjoint equations of an interface current nodal method are described. These techniques are the ''similarity transformation'' procedure and a direct solution scheme. A theoretical basis is provided for the similarity transformation procedure originally proposed by Lawrence. It is shown that the matrices associated with the mathematical and physical adjoint equations are similar to each other for the flat transverse leakage approximation but not for the quadratic leakage approximation. It is also shown that a good approximate solution of the mathematical adjoint for the quadratic transverse leakage approximation is obtained by applying the similarity transformation for the flat transverse leakage approximation to the physical adjoint solution. The direct solution scheme, which was developed as an alternative to the similarity transformation procedure, yields the correct mathematical adjoint solution for both flat and quadratic transverse leakage approximations. In this scheme, adjoint nodal equations are cast in a form very similar to that of the forward equations by employing a linear transformation of the adjoint partial currents. This enables the use of the forward solution algorithm with only minor modifications for solving the mathematical adjoint equations. By using the direct solution scheme as a reference method, it is shown that while the results computed with the similarity transformation procedure are approximate, they are sufficiently accurate for calculations of global and local reactivity changes resulting from coolant voiding in a liquid-metal reactor

Investigating Students' Mathematical Difficulties with Quadratic Equations

Science.gov (United States)

O'Connor, Bronwyn Reid; Norton, Stephen

2016-01-01

This paper examines the factors that hinder students' success in working with and understanding the mathematics of quadratic equations using a case study analysis of student error patterns. Twenty-five Year 11 students were administered a written test to examine their understanding of concepts and procedures associated with this topic. The…

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...

Numerical instability of time-discretized one-point kinetic equations

International Nuclear Information System (INIS)

Hashimoto, Kengo; Ikeda, Hideaki; Takeda, Toshikazu

2000-01-01

The one-point kinetic equations with numerical errors induced by the explicit, implicit and Crank-Nicolson integration methods are derived. The zero-power transfer functions based on the present equations are demonstrated to investigate the numerical stability of the discretized systems. These demonstrations indicate unconditional stability for the implicit and Crank-Nicolson methods but present the possibility of numerical instability for the explicit method. An upper limit of time mesh spacing for the stability is formulated and several numerical calculations are made to confirm the validity of this formula

Non-linear wave equations:Mathematical techniques

International Nuclear Information System (INIS)

1978-01-01

An account of certain well-established mathematical methods, which prove useful to deal with non-linear partial differential equations is presented. Within the strict framework of Functional Analysis, it describes Semigroup Techniques in Banach Spaces as well as variational approaches towards critical points. Detailed proofs are given of the existence of local and global solutions of the Cauchy problem and of the stability of stationary solutions. The formal approach based upon invariance under Lie transformations deserves attention due to its wide range of applicability, even if the explicit solutions thus obtained do not allow for a deep analysis of the equations. A compre ensive introduction to the inverse scattering approach and to the solution concept for certain non-linear equations of physical interest are also presented. A detailed discussion is made about certain convergence and stability problems which arise in importance need not be emphasized. (author) [es

The MATH--Open Source Application for Easier Learning of Numerical Mathematics

Science.gov (United States)

Glaser-Opitz, Henrich; Budajová, Kristina

2016-01-01

The article introduces a software application (MATH) supporting an education of Applied Mathematics, with focus on Numerical Mathematics. The MATH is an easy to use tool supporting various numerical methods calculations with graphical user interface and integrated plotting tool for graphical representation written in Qt with extensive use of Qwt…

Numerical analysis of systems of ordinary and stochastic differential equations

CERN Document Server

Artemiev, S S

1997-01-01

This text deals with numerical analysis of systems of both ordinary and stochastic differential equations. It covers numerical solution problems of the Cauchy problem for stiff ordinary differential equations (ODE) systems by Rosenbrock-type methods (RTMs).

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

Numerical optimization using flow equations

Science.gov (United States)

Punk, Matthias

2014-12-01

We develop a method for multidimensional optimization using flow equations. This method is based on homotopy continuation in combination with a maximum entropy approach. Extrema of the optimizing functional correspond to fixed points of the flow equation. While ideas based on Bayesian inference such as the maximum entropy method always depend on a prior probability, the additional step in our approach is to perform a continuous update of the prior during the homotopy flow. The prior probability thus enters the flow equation only as an initial condition. We demonstrate the applicability of this optimization method for two paradigmatic problems in theoretical condensed matter physics: numerical analytic continuation from imaginary to real frequencies and finding (variational) ground states of frustrated (quantum) Ising models with random or long-range antiferromagnetic interactions.

Transient three-dimensional thermal-hydraulic analysis of nuclear reactor fuel rod arrays: general equations and numerical scheme

International Nuclear Information System (INIS)

Wnek, W.J.; Ramshaw, J.D.; Trapp, J.A.; Hughes, E.D.; Solbrig, C.W.

1975-11-01

A mathematical model and a numerical solution scheme for thermal-hydraulic analysis of fuel rod arrays are given. The model alleviates the two major deficiencies associated with existing rod array analysis models, that of a correct transverse momentum equation and the capability of handling reversing and circulatory flows. Possible applications of the model include steady state and transient subchannel calculations as well as analysis of flows in heat exchangers, other engineering equipment, and porous media

Numerical Treatment of Degenerate Diffusion Equations via Feller's Boundary Classification, and Applications

Science.gov (United States)

Cacio, Emanuela; Cohn, Stephen E.; Spigler, Renato

2011-01-01

A numerical method is devised to solve a class of linear boundary-value problems for one-dimensional parabolic equations degenerate at the boundaries. Feller theory, which classifies the nature of the boundary points, is used to decide whether boundary conditions are needed to ensure uniqueness, and, if so, which ones they are. The algorithm is based on a suitable preconditioned implicit finite-difference scheme, grid, and treatment of the boundary data. Second-order accuracy, unconditional stability, and unconditional convergence of solutions of the finite-difference scheme to a constant as the time-step index tends to infinity are further properties of the method. Several examples, pertaining to financial mathematics, physics, and genetics, are presented for the purpose of illustration.

Mathematical and numerical models for eddy currents and magnetostatics with selected applications

CERN Document Server

Rappaz, Jacques

2013-01-01

This monograph addresses fundamental aspects of mathematical modeling and numerical solution methods of electromagnetic problems involving low frequencies, i.e. magnetostatic and eddy current problems which are rarely presented in the applied mathematics literature. In the first part, the authors introduce the mathematical models in a realistic context in view of their use for industrial applications. Several geometric configurations of electric conductors leading to different mathematical models are carefully derived and analyzed, and numerical methods for the solution of the obtained problem

Numerical solution of boundary-integral equations for molecular electrostatics.

Science.gov (United States)

Bardhan, Jaydeep P

2009-03-07

Numerous molecular processes, such as ion permeation through channel proteins, are governed by relatively small changes in energetics. As a result, theoretical investigations of these processes require accurate numerical methods. In the present paper, we evaluate the accuracy of two approaches to simulating boundary-integral equations for continuum models of the electrostatics of solvation. The analysis emphasizes boundary-element method simulations of the integral-equation formulation known as the apparent-surface-charge (ASC) method or polarizable-continuum model (PCM). In many numerical implementations of the ASC/PCM model, one forces the integral equation to be satisfied exactly at a set of discrete points on the boundary. We demonstrate in this paper that this approach to discretization, known as point collocation, is significantly less accurate than an alternative approach known as qualocation. Furthermore, the qualocation method offers this improvement in accuracy without increasing simulation time. Numerical examples demonstrate that electrostatic part of the solvation free energy, when calculated using the collocation and qualocation methods, can differ significantly; for a polypeptide, the answers can differ by as much as 10 kcal/mol (approximately 4% of the total electrostatic contribution to solvation). The applicability of the qualocation discretization to other integral-equation formulations is also discussed, and two equivalences between integral-equation methods are derived.

Technological pedagogical content knowledge of junior high school mathematics teachers in teaching linear equation

Science.gov (United States)

Wati, S.; Fitriana, L.; Mardiyana

2018-04-01

Linear equation is one of the topics in mathematics that are considered difficult. Student difficulties of understanding linear equation can be caused by lack of understanding this concept and the way of teachers teach. TPACK is a way to understand the complex relationships between teaching and content taught through the use of specific teaching approaches and supported by the right technology tools. This study aims to identify TPACK of junior high school mathematics teachers in teaching linear equation. The method used in the study was descriptive. In the first phase, a survey using a questionnaire was carried out on 45 junior high school mathematics teachers in teaching linear equation. While in the second phase, the interview involved three teachers. The analysis of data used were quantitative and qualitative technique. The result PCK revealed teachers emphasized developing procedural and conceptual knowledge through reliance on traditional in teaching linear equation. The result of TPK revealed teachers’ lower capacity to deal with the general information and communications technologies goals across the curriculum in teaching linear equation. The result indicated that PowerPoint constitutes TCK modal technological capability in teaching linear equation. The result of TPACK seems to suggest a low standard in teachers’ technological skills across a variety of mathematics education goals in teaching linear equation. This means that the ability of teachers’ TPACK in teaching linear equation still needs to be improved.

Exact Solutions for Nonlinear Differential Difference Equations in Mathematical Physics

Directory of Open Access Journals (Sweden)

Khaled A. Gepreel

2013-01-01

Full Text Available We modified the truncated expansion method to construct the exact solutions for some nonlinear differential difference equations in mathematical physics via the general lattice equation, the discrete nonlinear Schrodinger with a saturable nonlinearity, the quintic discrete nonlinear Schrodinger equation, and the relativistic Toda lattice system. Also, we put a rational solitary wave function method to find the rational solitary wave solutions for some nonlinear differential difference equations. The proposed methods are more effective and powerful to obtain the exact solutions for nonlinear difference differential equations.

Numerical study of fractional nonlinear Schrodinger equations

KAUST Repository

Klein, Christian

2014-10-08

Using a Fourier spectral method, we provide a detailed numerical investigation of dispersive Schrödinger-type equations involving a fractional Laplacian in an one-dimensional case. By an appropriate choice of the dispersive exponent, both mass and energy sub- and supercritical regimes can be identified. This allows us to study the possibility of finite time blow-up versus global existence, the nature of the blow-up, the stability and instability of nonlinear ground states and the long-time dynamics of solutions. The latter is also studied in a semiclassical setting. Moreover, we numerically construct ground state solutions of the fractional nonlinear Schrödinger equation.

Differential equations a dynamical systems approach ordinary differential equations

CERN Document Server

Hubbard, John H

1991-01-01

This is a corrected third printing of the first part of the text Differential Equations: A Dynamical Systems Approach written by John Hubbard and Beverly West. The authors' main emphasis in this book is on ordinary differential equations. The book is most appropriate for upper level undergraduate and graduate students in the fields of mathematics, engineering, and applied mathematics, as well as the life sciences, physics and economics. Traditional courses on differential equations focus on techniques leading to solutions. Yet most differential equations do not admit solutions which can be written in elementary terms. The authors have taken the view that a differential equations defines functions; the object of the theory is to understand the behavior of these functions. The tools the authors use include qualitative and numerical methods besides the traditional analytic methods. The companion software, MacMath, is designed to bring these notions to life.

The numerical simulation of convection delayed dominated diffusion equation

Directory of Open Access Journals (Sweden)

Mohan Kumar P. Murali

2016-01-01

Full Text Available In this paper, we propose a fitted numerical method for solving convection delayed dominated diffusion equation. A fitting factor is introduced and the model equation is discretized by cubic spline method. The error analysis is analyzed for the consider problem. The numerical examples are solved using the present method and compared the result with the exact solution.

Representations of Numerical and Non-Numerical Magnitude Both Contribute to Mathematical Competence in Children

Science.gov (United States)

Lourenco, Stella F.; Bonny, Justin W.

2017-01-01

A growing body of evidence suggests that non-symbolic representations of number, which humans share with nonhuman animals, are functionally related to uniquely human mathematical thought. Other research suggesting that numerical and non-numerical magnitudes not only share analog format but also form part of a general magnitude system raises…

A new numerical approximation of the fractal ordinary differential equation

Science.gov (United States)

Atangana, Abdon; Jain, Sonal

2018-02-01

The concept of fractal medium is present in several real-world problems, for instance, in the geological formation that constitutes the well-known subsurface water called aquifers. However, attention has not been quite devoted to modeling for instance, the flow of a fluid within these media. We deem it important to remind the reader that the concept of fractal derivative is not to represent the fractal sharps but to describe the movement of the fluid within these media. Since this class of ordinary differential equations is highly complex to solve analytically, we present a novel numerical scheme that allows to solve fractal ordinary differential equations. Error analysis of the method is also presented. Application of the method and numerical approximation are presented for fractal order differential equation. The stability and the convergence of the numerical schemes are investigated in detail. Also some exact solutions of fractal order differential equations are presented and finally some numerical simulations are presented.

Numerical capacities as domain-specific predictors beyond early mathematics learning: a longitudinal study.

Science.gov (United States)

Reigosa-Crespo, Vivian; González-Alemañy, Eduardo; León, Teresa; Torres, Rosario; Mosquera, Raysil; Valdés-Sosa, Mitchell

2013-01-01

The first aim of the present study was to investigate whether numerical effects (Numerical Distance Effect, Counting Effect and Subitizing Effect) are domain-specific predictors of mathematics development at the end of elementary school by exploring whether they explain additional variance of later mathematics fluency after controlling for the effects of general cognitive skills, focused on nonnumerical aspects. The second aim was to address the same issues but applied to achievement in mathematics curriculum that requires solutions to fluency in calculation. These analyses assess whether the relationship found for fluency are generalized to mathematics content beyond fluency in calculation. As a third aim, the domain specificity of the numerical effects was examined by analyzing whether they contribute to the development of reading skills, such as decoding fluency and reading comprehension, after controlling for general cognitive skills and phonological processing. Basic numerical capacities were evaluated in children of 3(rd) and 4(th) grades (n=49). Mathematics and reading achievements were assessed in these children one year later. Results showed that the size of the Subitizing Effect was a significant domain-specific predictor of fluency in calculation and also in curricular mathematics achievement, but not in reading skills, assessed at the end of elementary school. Furthermore, the size of the Counting Effect also predicted fluency in calculation, although this association only approached significance. These findings contrast with proposals that the core numerical competencies measured by enumeration will bear little relationship to mathematics achievement. We conclude that basic numerical capacities constitute domain-specific predictors and that they are not exclusively "start-up" tools for the acquisition of Mathematics; but they continue modulating this learning at the end of elementary school.

Spurious solutions in few-body equations. II. Numerical investigations

International Nuclear Information System (INIS)

Adhikari, S.K.

1979-01-01

A recent analytic study of spurious solutions in few-body equations by Adhikari and Gloeckle is here complemented by numerical investigations. As proposed by Adhikari and Gloeckle we study numerically the spurious solutions in the three-body Weinberg type equations and draw some general conclusions about the existence of spurious solutions in three-body equations with the Weinberg kernel and in other few-body formulations. In particular we conclude that for most of the potentials we encounter in problems of nuclear physics the three-body Weinberg type equation will not have a spurious solution which may interfere with the bound state or scattering calculation. Hence, if proven convenient, the three-body Weinberg type equation can be used in practical calculations. The same conclusion is true for the three-body channel coupling array scheme of Kouri, Levin, and Tobocman. In the case of the set of six coupled four-body equations proposed by Rosenberg et al. and the set of the Bencze-Redish-Sloan equations a careful study of the possible spurious solutions is needed before using these equations in practical calculations

Numerical Magnitude Representation in Children With Mathematical Difficulties With or Without Reading Difficulties.

Science.gov (United States)

Tobia, Valentina; Fasola, Anna; Lupieri, Alice; Marzocchi, Gian Marco

2016-01-01

This study aimed to explore the spatial numerical association of response codes (SNARC), the flanker, and the numerical distance effects in children with mathematical difficulties. From a sample of 720 third, fourth, and fifth graders, 60 children were selected and divided into the following three groups: typically developing children (TD; n = 29), children with mathematical difficulties only (MD only; n = 21), and children with mathematical and reading difficulties (MD+RD; n = 10). Children were tested with a numerical Eriksen task that was built to assess SNARC, numerical distance, and flanker (first and second order congruency) effects. Children with MD only showed stronger SNARC and second order congruency effects than did TD children, whereas the numerical distance effects were similar across the three groups. Finally, the first order congruency effect was associated with reading difficulties. These results showed that children with mathematical difficulties with or without reading difficulties were globally more impaired when spatial incompatibilities were presented. © Hammill Institute on Disabilities 2014.

Mathematics Literacy of Secondary Students in Solving Simultanenous Linear Equations

Science.gov (United States)

Sitompul, R. S. I.; Budayasa, I. K.; Masriyah

2018-01-01

This study examines the profile of secondary students’ mathematical literacy in solving simultanenous linear equations problems in terms of cognitive style of visualizer and verbalizer. This research is a descriptive research with qualitative approach. The subjects in this research consist of one student with cognitive style of visualizer and one student with cognitive style of verbalizer. The main instrument in this research is the researcher herself and supporting instruments are cognitive style tests, mathematics skills tests, problem-solving tests and interview guidelines. Research was begun by determining the cognitive style test and mathematics skill test. The subjects chosen were given problem-solving test about simultaneous linear equations and continued with interview. To ensure the validity of the data, the researcher conducted data triangulation; the steps of data reduction, data presentation, data interpretation, and conclusion drawing. The results show that there is a similarity of visualizer and verbalizer-cognitive style in identifying and understanding the mathematical structure in the process of formulating. There are differences in how to represent problems in the process of implementing, there are differences in designing strategies and in the process of interpreting, and there are differences in explaining the logical reasons.

Numerical algebra, matrix theory, differential-algebraic equations and control theory festschrift in honor of Volker Mehrmann

CERN Document Server

Bollhöfer, Matthias; Kressner, Daniel; Mehl, Christian; Stykel, Tatjana

2015-01-01

This edited volume highlights the scientific contributions of Volker Mehrmann, a leading expert in the area of numerical (linear) algebra, matrix theory, differential-algebraic equations and control theory. These mathematical research areas are strongly related and often occur in the same real-world applications. The main areas where such applications emerge are computational engineering and sciences, but increasingly also social sciences and economics. This book also reflects some of Volker Mehrmann's major career stages. Starting out working in the areas of numerical linear algebra (his first full professorship at TU Chemnitz was in "Numerical Algebra," hence the title of the book) and matrix theory, Volker Mehrmann has made significant contributions to these areas ever since. The highlights of these are discussed in Parts I and II of the present book. Often the development of new algorithms in numerical linear algebra is motivated by problems in system and control theory. These and his later major work on ...

Numerical Capacities as Domain-Specific Predictors beyond Early Mathematics Learning: A Longitudinal Study

Science.gov (United States)

Reigosa-Crespo, Vivian; González-Alemañy, Eduardo; León, Teresa; Torres, Rosario; Mosquera, Raysil; Valdés-Sosa, Mitchell

2013-01-01

The first aim of the present study was to investigate whether numerical effects (Numerical Distance Effect, Counting Effect and Subitizing Effect) are domain-specific predictors of mathematics development at the end of elementary school by exploring whether they explain additional variance of later mathematics fluency after controlling for the effects of general cognitive skills, focused on nonnumerical aspects. The second aim was to address the same issues but applied to achievement in mathematics curriculum that requires solutions to fluency in calculation. These analyses assess whether the relationship found for fluency are generalized to mathematics content beyond fluency in calculation. As a third aim, the domain specificity of the numerical effects was examined by analyzing whether they contribute to the development of reading skills, such as decoding fluency and reading comprehension, after controlling for general cognitive skills and phonological processing. Basic numerical capacities were evaluated in children of 3rd and 4th grades (n=49). Mathematics and reading achievements were assessed in these children one year later. Results showed that the size of the Subitizing Effect was a significant domain-specific predictor of fluency in calculation and also in curricular mathematics achievement, but not in reading skills, assessed at the end of elementary school. Furthermore, the size of the Counting Effect also predicted fluency in calculation, although this association only approached significance. These findings contrast with proposals that the core numerical competencies measured by enumeration will bear little relationship to mathematics achievement. We conclude that basic numerical capacities constitute domain-specific predictors and that they are not exclusively “start-up” tools for the acquisition of Mathematics; but they continue modulating this learning at the end of elementary school. PMID:24255710

A numerical scheme for the generalized Burgers–Huxley equation

Directory of Open Access Journals (Sweden)

Brajesh K. Singh

2016-10-01

Full Text Available In this article, a numerical solution of generalized Burgers–Huxley (gBH equation is approximated by using a new scheme: modified cubic B-spline differential quadrature method (MCB-DQM. The scheme is based on differential quadrature method in which the weighting coefficients are obtained by using modified cubic B-splines as a set of basis functions. This scheme reduces the equation into a system of first-order ordinary differential equation (ODE which is solved by adopting SSP-RK43 scheme. Further, it is shown that the proposed scheme is stable. The efficiency of the proposed method is illustrated by four numerical experiments, which confirm that obtained results are in good agreement with earlier studies. This scheme is an easy, economical and efficient technique for finding numerical solutions for various kinds of (nonlinear physical models as compared to the earlier schemes.

Symbolic-Numeric Integration of the Dynamical Cosserat Equations

KAUST Repository

Lyakhov, Dmitry A.

2017-08-29

We devise a symbolic-numeric approach to the integration of the dynamical part of the Cosserat equations, a system of nonlinear partial differential equations describing the mechanical behavior of slender structures, like fibers and rods. This is based on our previous results on the construction of a closed form general solution to the kinematic part of the Cosserat system. Our approach combines methods of numerical exponential integration and symbolic integration of the intermediate system of nonlinear ordinary differential equations describing the dynamics of one of the arbitrary vector-functions in the general solution of the kinematic part in terms of the module of the twist vector-function. We present an experimental comparison with the well-established generalized \\\\alpha -method illustrating the computational efficiency of our approach for problems in structural mechanics.

Symbolic-Numeric Integration of the Dynamical Cosserat Equations

KAUST Repository

Lyakhov, Dmitry A.; Gerdt, Vladimir P.; Weber, Andreas G.; Michels, Dominik L.

2017-01-01

We devise a symbolic-numeric approach to the integration of the dynamical part of the Cosserat equations, a system of nonlinear partial differential equations describing the mechanical behavior of slender structures, like fibers and rods. This is based on our previous results on the construction of a closed form general solution to the kinematic part of the Cosserat system. Our approach combines methods of numerical exponential integration and symbolic integration of the intermediate system of nonlinear ordinary differential equations describing the dynamics of one of the arbitrary vector-functions in the general solution of the kinematic part in terms of the module of the twist vector-function. We present an experimental comparison with the well-established generalized \\alpha -method illustrating the computational efficiency of our approach for problems in structural mechanics.

Jacobi Elliptic Solutions for Nonlinear Differential Difference Equations in Mathematical Physics

Directory of Open Access Journals (Sweden)

Khaled A. Gepreel

2012-01-01

Full Text Available We put a direct new method to construct the rational Jacobi elliptic solutions for nonlinear differential difference equations which may be called the rational Jacobi elliptic functions method. We use the rational Jacobi elliptic function method to construct many new exact solutions for some nonlinear differential difference equations in mathematical physics via the lattice equation and the discrete nonlinear Schrodinger equation with a saturable nonlinearity. The proposed method is more effective and powerful to obtain the exact solutions for nonlinear differential difference equations.

The numerical solution of linear multi-term fractional differential equations: systems of equations

Science.gov (United States)

Edwards, John T.; Ford, Neville J.; Simpson, A. Charles

2002-11-01

In this paper, we show how the numerical approximation of the solution of a linear multi-term fractional differential equation can be calculated by reduction of the problem to a system of ordinary and fractional differential equations each of order at most unity. We begin by showing how our method applies to a simple class of problems and we give a convergence result. We solve the Bagley Torvik equation as an example. We show how the method can be applied to a general linear multi-term equation and give two further examples.

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

KAUST Repository

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

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

Numerical investigation of sixth order Boussinesq equation

Science.gov (United States)

Kolkovska, N.; Vucheva, V.

2017-10-01

We propose a family of conservative finite difference schemes for the Boussinesq equation with sixth order dispersion terms. The schemes are of second order of approximation. The method is conditionally stable with a mild restriction τ = O(h) on the step sizes. Numerical tests are performed for quadratic and cubic nonlinearities. The numerical experiments show second order of convergence of the discrete solution to the exact one.

Numerical Solution of Stochastic Nonlinear Fractional Differential Equations

KAUST Repository

El-Beltagy, Mohamed A.

2015-01-07

Using Wiener-Hermite expansion (WHE) technique in the solution of the stochastic partial differential equations (SPDEs) has the advantage of converting the problem to a system of deterministic equations that can be solved efficiently using the standard deterministic numerical methods [1]. WHE is the only known expansion that handles the white/colored noise exactly. This work introduces a numerical estimation of the stochastic response of the Duffing oscillator with fractional or variable order damping and driven by white noise. The WHE technique is integrated with the Grunwald-Letnikov approximation in case of fractional order and with Coimbra approximation in case of variable-order damping. The numerical solver was tested with the analytic solution and with Monte-Carlo simulations. The developed mixed technique was shown to be efficient in simulating SPDEs.

Numerical Solution of Stochastic Nonlinear Fractional Differential Equations

KAUST Repository

El-Beltagy, Mohamed A.; Al-Juhani, Amnah

2015-01-01

Using Wiener-Hermite expansion (WHE) technique in the solution of the stochastic partial differential equations (SPDEs) has the advantage of converting the problem to a system of deterministic equations that can be solved efficiently using the standard deterministic numerical methods [1]. WHE is the only known expansion that handles the white/colored noise exactly. This work introduces a numerical estimation of the stochastic response of the Duffing oscillator with fractional or variable order damping and driven by white noise. The WHE technique is integrated with the Grunwald-Letnikov approximation in case of fractional order and with Coimbra approximation in case of variable-order damping. The numerical solver was tested with the analytic solution and with Monte-Carlo simulations. The developed mixed technique was shown to be efficient in simulating SPDEs.

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.

Numerical Analysis of Partial Differential Equations

CERN Document Server

Lions, Jacques-Louis

2011-01-01

S. Albertoni: Alcuni metodi di calcolo nella teoria della diffusione dei neutroni.- I. Babuska: Optimization and numerical stability in computations.- J.H. Bramble: Error estimates in elliptic boundary value problems.- G. Capriz: The numerical approach to hydrodynamic problems.- A. Dou: Energy inequalities in an elastic cylinder.- T. Doupont: On the existence of an iterative method for the solution of elliptic difference equation with an improved work estimate.- J. Douglas, J.R. Cannon: The approximation of harmonic and parabolic functions of half-spaces from interior data.- B.E. Hubbard: Erro

Introduction to computation and modeling for differential equations

CERN Document Server

Edsberg, Lennart

2008-01-01

An introduction to scientific computing for differential equationsIntroduction to Computation and Modeling for Differential Equations provides a unified and integrated view of numerical analysis, mathematical modeling in applications, and programming to solve differential equations, which is essential in problem-solving across many disciplines, such as engineering, physics, and economics. This book successfully introduces readers to the subject through a unique ""Five-M"" approach: Modeling, Mathematics, Methods, MATLAB, and Multiphysics. This approach facilitates a thorough understanding of h

New Numerical Solution of von Karman Equation of Lengthwise Rolling

Directory of Open Access Journals (Sweden)

Rudolf Pernis

2015-01-01

Full Text Available The calculation of average material contact pressure to rolls base on mathematical theory of rolling process given by Karman equation was solved by many authors. The solutions reported by authors are used simplifications for solution of Karman equation. The simplifications are based on two cases for approximation of the circular arch: (a by polygonal curve and (b by parabola. The contribution of the present paper for solution of two-dimensional differential equation of rolling is based on description of the circular arch by equation of a circle. The new term relative stress as nondimensional variable was defined. The result from derived mathematical models can be calculated following variables: normal contact stress distribution, front and back tensions, angle of neutral point, coefficient of the arm of rolling force, rolling force, and rolling torque during rolling process. Laboratory cold rolled experiment of CuZn30 brass material was performed. Work hardening during brass processing was calculated. Comparison of theoretical values of normal contact stress with values of normal contact stress obtained from cold rolling experiment was performed. The calculations were not concluded with roll flattening.

Numerical models for differential problems

CERN Document Server

Quarteroni, Alfio

2017-01-01

In this text, we introduce the basic concepts for the numerical modelling of partial differential equations. We consider the classical elliptic, parabolic and hyperbolic linear equations, but also the diffusion, transport, and Navier-Stokes equations, as well as equations representing conservation laws, saddle-point problems and optimal control problems. Furthermore, we provide numerous physical examples which underline such equations. We then analyze numerical solution methods based on finite elements, finite differences, finite volumes, spectral methods and domain decomposition methods, and reduced basis methods. In particular, we discuss the algorithmic and computer implementation aspects and provide a number of easy-to-use programs. The text does not require any previous advanced mathematical knowledge of partial differential equations: the absolutely essential concepts are reported in a preliminary chapter. It is therefore suitable for students of bachelor and master courses in scientific disciplines, an...

Numerical solution of second-order stochastic differential equations with Gaussian random parameters

Directory of Open Access Journals (Sweden)

Rahman Farnoosh

2014-07-01

Full Text Available In this paper, we present the numerical solution of ordinary differential equations (or SDEs, from each orderespecially second-order with time-varying and Gaussian random coefficients. We indicate a complete analysisfor second-order equations in specially case of scalar linear second-order equations (damped harmonicoscillators with additive or multiplicative noises. Making stochastic differential equations system from thisequation, it could be approximated or solved numerically by different numerical methods. In the case oflinear stochastic differential equations system by Computing fundamental matrix of this system, it could becalculated based on the exact solution of this system. Finally, this stochastic equation is solved by numericallymethod like E.M. and Milstein. Also its Asymptotic stability and statistical concepts like expectationand variance of solutions are discussed.

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.

Elliptic differential equations theory and numerical treatment

CERN Document Server

Hackbusch, Wolfgang

2017-01-01

This book simultaneously presents the theory and the numerical treatment of elliptic boundary value problems, since an understanding of the theory is necessary for the numerical analysis of the discretisation. It first discusses the Laplace equation and its finite difference discretisation before addressing the general linear differential equation of second order. The variational formulation together with the necessary background from functional analysis provides the basis for the Galerkin and finite-element methods, which are explored in detail. A more advanced chapter leads the reader to the theory of regularity. Individual chapters are devoted to singularly perturbed as well as to elliptic eigenvalue problems. The book also presents the Stokes problem and its discretisation as an example of a saddle-point problem taking into account its relevance to applications in fluid dynamics.

Direct numerical methods of mathematical modeling in mechanical structural design

International Nuclear Information System (INIS)

Sahili, Jihad; Verchery, Georges; Ghaddar, Ahmad; Zoaeter, Mohamed

2002-01-01

Full text.Structural design and numerical methods are generally interactive; requiring optimization procedures as the structure is analyzed. This analysis leads to define some mathematical terms, as the stiffness matrix, which are resulting from the modeling and then used in numerical techniques during the dimensioning procedure. These techniques and many others involve the calculation of the generalized inverse of the stiffness matrix, called also the 'compliance matrix'. The aim of this paper is to introduce first, some different existing mathematical procedures, used to calculate the compliance matrix from the stiffness matrix, then apply direct numerical methods to solve the obtained system with the lowest computational time, and to compare the obtained results. The results show a big difference of the computational time between the different procedures

International Conference on Differential Equations and Mathematical Physics

CERN Document Server

Saitō, Yoshimi

1987-01-01

The meeting in Birmingham, Alabama, provided a forum for the discussion of recent developments in the theory of ordinary and partial differential equations, both linear and non-linear, with particular reference to work relating to the equations of mathematical physics. The meeting was attended by about 250 mathematicians from 22 countries. The papers in this volume all involve new research material, with at least outline proofs; some papers also contain survey material. Topics covered include: Schrödinger theory, scattering and inverse scattering, fluid mechanics (including conservative systems and inertial manifold theory attractors), elasticity, non-linear waves, and feedback control theory.

Numerical Solutions for Convection-Diffusion Equation through Non-Polynomial Spline

Directory of Open Access Journals (Sweden)

Ravi Kanth A.S.V.

2016-01-01

Full Text Available In this paper, numerical solutions for convection-diffusion equation via non-polynomial splines are studied. We purpose an implicit method based on non-polynomial spline functions for solving the convection-diffusion equation. The method is proven to be unconditionally stable by using Von Neumann technique. Numerical results are illustrated to demonstrate the efficiency and stability of the purposed method.

Mathematical methods for physical and analytical chemistry

CERN Document Server

Goodson, David Z

2011-01-01

Mathematical Methods for Physical and Analytical Chemistry presents mathematical and statistical methods to students of chemistry at the intermediate, post-calculus level. The content includes a review of general calculus; a review of numerical techniques often omitted from calculus courses, such as cubic splines and Newton's method; a detailed treatment of statistical methods for experimental data analysis; complex numbers; extrapolation; linear algebra; and differential equations. With numerous example problems and helpful anecdotes, this text gives chemistry students the mathematical

Mathematics in physics and engineering

CERN Document Server

Irving, J; Massey, H S W; Brueckner, Keith A

1959-01-01

Mathematics in Physics and Engineering describes the analytical and numerical (desk-machine) methods that arise in pure and applied science, including wave equations, Bessel and Legendre functions, and matrices. The manuscript first discusses partial differential equations, as well as the method of separation of variables, three-dimensional wave equation, diffusion or heat flow equation, and wave equation in plane and cylindrical polar coordinates. The text also ponders on Frobenius' and other methods of solution. Discussions focus on hypergeometric equation, Bessel's equation, confluent hyper

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.

A systematic literature review of Burgers' equation with recent advances

Science.gov (United States)

Bonkile, Mayur P.; Awasthi, Ashish; Lakshmi, C.; Mukundan, Vijitha; Aswin, V. S.

2018-06-01

Even if numerical simulation of the Burgers' equation is well documented in the literature, a detailed literature survey indicates that gaps still exist for comparative discussion regarding the physical and mathematical significance of the Burgers' equation. Recently, an increasing interest has been developed within the scientific community, for studying non-linear convective-diffusive partial differential equations partly due to the tremendous improvement in computational capacity. Burgers' equation whose exact solution is well known, is one of the famous non-linear partial differential equations which is suitable for the analysis of various important areas. A brief historical review of not only the mathematical, but also the physical significance of the solution of Burgers' equation is presented, emphasising current research strategies, and the challenges that remain regarding the accuracy, stability and convergence of various schemes are discussed. One of the objectives of this paper is to discuss the recent developments in mathematical modelling of Burgers' equation and thus open doors for improvement. No claim is made that the content of the paper is new. However, it is a sincere effort to outline the physical and mathematical importance of Burgers' equation in the most simplified ways. We throw some light on the plethora of challenges which need to be overcome in the research areas and give motivation for the next breakthrough to take place in a numerical simulation of ordinary / partial differential equations.

Numerical solution of special ultra-relativistic Euler equations using central upwind scheme

Science.gov (United States)

Ghaffar, Tayabia; Yousaf, Muhammad; Qamar, Shamsul

2018-06-01

This article is concerned with the numerical approximation of one and two-dimensional special ultra-relativistic Euler equations. The governing equations are coupled first-order nonlinear hyperbolic partial differential equations. These equations describe perfect fluid flow in terms of the particle density, the four-velocity and the pressure. A high-resolution shock-capturing central upwind scheme is employed to solve the model equations. To avoid excessive numerical diffusion, the considered scheme avails the specific information of local propagation speeds. By using Runge-Kutta time stepping method and MUSCL-type initial reconstruction, we have obtained 2nd order accuracy of the proposed scheme. After discussing the model equations and the numerical technique, several 1D and 2D test problems are investigated. For all the numerical test cases, our proposed scheme demonstrates very good agreement with the results obtained by well-established algorithms, even in the case of highly relativistic 2D test problems. For validation and comparison, the staggered central scheme and the kinetic flux-vector splitting (KFVS) method are also implemented to the same model. The robustness and efficiency of central upwind scheme is demonstrated by the numerical results.

Inferring Mathematical Equations Using Crowdsourcing.

Directory of Open Access Journals (Sweden)

Szymon Wasik

Full Text Available Crowdsourcing, understood as outsourcing work to a large network of people in the form of an open call, has been utilized successfully many times, including a very interesting concept involving the implementation of computer games with the objective of solving a scientific problem by employing users to play a game-so-called crowdsourced serious games. Our main objective was to verify whether such an approach could be successfully applied to the discovery of mathematical equations that explain experimental data gathered during the observation of a given dynamic system. Moreover, we wanted to compare it with an approach based on artificial intelligence that uses symbolic regression to find such formulae automatically. To achieve this, we designed and implemented an Internet game in which players attempt to design a spaceship representing an equation that models the observed system. The game was designed while considering that it should be easy to use for people without strong mathematical backgrounds. Moreover, we tried to make use of the collective intelligence observed in crowdsourced systems by enabling many players to collaborate on a single solution. The idea was tested on several hundred players playing almost 10,000 games and conducting a user opinion survey. The results prove that the proposed solution has very high potential. The function generated during weeklong tests was almost as precise as the analytical solution of the model of the system and, up to a certain complexity level of the formulae, it explained data better than the solution generated automatically by Eureqa, the leading software application for the implementation of symbolic regression. Moreover, we observed benefits of using crowdsourcing; the chain of consecutive solutions that led to the best solution was obtained by the continuous collaboration of several players.

Inferring Mathematical Equations Using Crowdsourcing.

Science.gov (United States)

Wasik, Szymon; Fratczak, Filip; Krzyskow, Jakub; Wulnikowski, Jaroslaw

2015-01-01

Crowdsourcing, understood as outsourcing work to a large network of people in the form of an open call, has been utilized successfully many times, including a very interesting concept involving the implementation of computer games with the objective of solving a scientific problem by employing users to play a game-so-called crowdsourced serious games. Our main objective was to verify whether such an approach could be successfully applied to the discovery of mathematical equations that explain experimental data gathered during the observation of a given dynamic system. Moreover, we wanted to compare it with an approach based on artificial intelligence that uses symbolic regression to find such formulae automatically. To achieve this, we designed and implemented an Internet game in which players attempt to design a spaceship representing an equation that models the observed system. The game was designed while considering that it should be easy to use for people without strong mathematical backgrounds. Moreover, we tried to make use of the collective intelligence observed in crowdsourced systems by enabling many players to collaborate on a single solution. The idea was tested on several hundred players playing almost 10,000 games and conducting a user opinion survey. The results prove that the proposed solution has very high potential. The function generated during weeklong tests was almost as precise as the analytical solution of the model of the system and, up to a certain complexity level of the formulae, it explained data better than the solution generated automatically by Eureqa, the leading software application for the implementation of symbolic regression. Moreover, we observed benefits of using crowdsourcing; the chain of consecutive solutions that led to the best solution was obtained by the continuous collaboration of several players.

BOOK REVIEW: Advanced Topics in Computational Partial Differential Equations: Numerical Methods and Diffpack Programming

Science.gov (United States)

Katsaounis, T. D.

2005-02-01

The scope of this book is to present well known simple and advanced numerical methods for solving partial differential equations (PDEs) and how to implement these methods using the programming environment of the software package Diffpack. A basic background in PDEs and numerical methods is required by the potential reader. Further, a basic knowledge of the finite element method and its implementation in one and two space dimensions is required. The authors claim that no prior knowledge of the package Diffpack is required, which is true, but the reader should be at least familiar with an object oriented programming language like C++ in order to better comprehend the programming environment of Diffpack. Certainly, a prior knowledge or usage of Diffpack would be a great advantage to the reader. The book consists of 15 chapters, each one written by one or more authors. Each chapter is basically divided into two parts: the first part is about mathematical models described by PDEs and numerical methods to solve these models and the second part describes how to implement the numerical methods using the programming environment of Diffpack. Each chapter closes with a list of references on its subject. The first nine chapters cover well known numerical methods for solving the basic types of PDEs. Further, programming techniques on the serial as well as on the parallel implementation of numerical methods are also included in these chapters. The last five chapters are dedicated to applications, modelled by PDEs, in a variety of fields. The first chapter is an introduction to parallel processing. It covers fundamentals of parallel processing in a simple and concrete way and no prior knowledge of the subject is required. Examples of parallel implementation of basic linear algebra operations are presented using the Message Passing Interface (MPI) programming environment. Here, some knowledge of MPI routines is required by the reader. Examples solving in parallel simple PDEs using

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

Numerical study of fractional nonlinear Schrodinger equations

KAUST Repository

Klein, Christian; Sparber, Christof; Markowich, Peter A.

2014-01-01

Using a Fourier spectral method, we provide a detailed numerical investigation of dispersive Schrödinger-type equations involving a fractional Laplacian in an one-dimensional case. By an appropriate choice of the dispersive exponent, both mass

Handbook of integral equations

CERN Document Server

Polyanin, Andrei D

2008-01-01

This handbook contains over 2,500 integral equations with solutions as well as analytical and numerical methods for solving linear and nonlinear equations. It explores Volterra, Fredholm, WienerHopf, Hammerstein, Uryson, and other equations that arise in mathematics, physics, engineering, the sciences, and economics. This second edition includes new chapters on mixed multidimensional equations and methods of integral equations for ODEs and PDEs, along with over 400 new equations with exact solutions. With many examples added for illustrative purposes, it presents new material on Volterra, Fredholm, singular, hypersingular, dual, and nonlinear integral equations, integral transforms, and special functions.

Differential equations and applications recent advances

CERN Document Server

2014-01-01

Differential Equations and Applications : Recent Advances focus on the latest developments in Nonlinear Dynamical Systems, Neural Networks, Fluid Dynamics, Fractional Differential Systems, Mathematical Modelling and Qualitative Theory. Different aspects such as Existence, Stability, Controllability, Viscosity and Numerical Analysis for different systems have been discussed in this book. This book will be of great interest and use to researchers in Applied Mathematics, Engineering and Mathematical Physics.

Quasi-gas dynamic equations

CERN Document Server

Elizarova, Tatiana G

2009-01-01

This book presents two interconnected mathematical models generalizing the Navier-Stokes system. The models, called the quasi-gas-dynamic and quasi-hydrodynamic equations, are then used as the basis of numerical methods solving gas- and fluid-dynamic problems.

Theoretical numerical analysis a functional analysis framework

CERN Document Server

Atkinson, Kendall

2005-01-01

This textbook prepares graduate students for research in numerical analysis/computational mathematics by giving to them a mathematical framework embedded in functional analysis and focused on numerical analysis. This helps the student to move rapidly into a research program. The text covers basic results of functional analysis, approximation theory, Fourier analysis and wavelets, iteration methods for nonlinear equations, finite difference methods, Sobolev spaces and weak formulations of boundary value problems, finite element methods, elliptic variational inequalities and their numerical solu

Brain correlates of mathematical competence in processing mathematical representations

Directory of Open Access Journals (Sweden)

Roland H. Grabner

2011-11-01

Full Text Available The ability to extract numerical information from different representation formats (e.g., equations, tables, or diagrams is a key component of mathematical competence but little is known about its neural correlate. Previous studies comparing mathematically less and more competent adults have focused on mental arithmetic and reported differences in left angular gyrus activity which were interpreted to reflect differential reliance on arithmetic fact retrieval during problem solving. The aim of the present functional magnetic resonance imaging (fMRI study was to investigate the brain correlates of mathematical competence in a task requiring the processing of typical mathematical representations. Twenty-eight adults of lower and higher mathematical competence worked on a representation matching task in which they had to evaluate whether the numerical information of a symbolic equation matches that of a bar chart. Two task conditions without and one condition with arithmetic demands were administered. Both competence groups performed equally well in the non-arithmetic conditions and only differed in accuracy in the condition requiring calculation. Activation contrasts between the groups revealed consistently stronger left angular gyrus activation in the more competent individuals across all three task conditions. The finding of competence-related activation differences independently of arithmetic demands suggests that more and less competent individuals differ in a cognitive process other than arithmetic fact retrieval. Specifically, it is argued that the stronger left angular gyrus activity in the more competent adults may reflect their higher proficiency in processing mathematical symbols. Moreover, the study demonstrates competence-related parietal activation differences that were not accompanied by differential experimental performance.

Mathematical analysis of partial differential equations modeling electrostatic MEMS

CERN Document Server

Esposito, Pierpaolo; Guo, Yujin

2010-01-01

Micro- and nanoelectromechanical systems (MEMS and NEMS), which combine electronics with miniature-size mechanical devices, are essential components of modern technology. It is the mathematical model describing "electrostatically actuated" MEMS that is addressed in this monograph. Even the simplified models that the authors deal with still lead to very interesting second- and fourth-order nonlinear elliptic equations (in the stationary case) and to nonlinear parabolic equations (in the dynamic case). While nonlinear eigenvalue problems-where the stationary MEMS models fit-are a well-developed

Approximate numerical abilities and mathematics: Insight from correlational and experimental training studies.

Science.gov (United States)

Hyde, D C; Berteletti, I; Mou, Y

2016-01-01

Humans have the ability to nonverbally represent the approximate numerosity of sets of objects. The cognitive system that supports this ability, often referred to as the approximate number system (ANS), is present in early infancy and continues to develop in precision over the life span. It has been proposed that the ANS forms a foundation for uniquely human symbolic number and mathematics learning. Recent work has brought two types of evidence to bear on the relationship between the ANS and human mathematics: correlational studies showing individual differences in approximate numerical abilities correlate with individual differences in mathematics achievement and experimental studies showing enhancing effects of nonsymbolic approximate numerical training on exact, symbolic mathematical abilities. From this work, at least two accounts can be derived from these empirical data. It may be the case that the ANS and mathematics are related because the cognitive and brain processes responsible for representing numerical quantity in each format overlap, the Representational Overlap Hypothesis, or because of commonalities in the cognitive operations involved in mentally manipulating the representations of each format, the Operational Overlap hypothesis. The two hypotheses make distinct predictions for future work to test. © 2016 Elsevier B.V. All rights reserved.

Applied mathematics

International Nuclear Information System (INIS)

Nedelec, J.C.

1988-01-01

The 1988 progress report of the Applied Mathematics center (Polytechnic School, France), is presented. The research fields of the Center are the scientific calculus, the probabilities and statistics and the video image synthesis. The research topics developed are: the analysis of numerical methods, the mathematical analysis of the physics and mechanics fundamental models, the numerical solution of complex models related to the industrial problems, the stochastic calculus and the brownian movement, the stochastic partial differential equations, the identification of the adaptive filtering parameters, the discrete element systems, statistics, the stochastic control and the development, the image synthesis techniques for education and research programs. The published papers, the congress communications and the thesis are listed [fr

Large quantum systems: a mathematical and numerical perspective; Systemes quantiques a grand nombre de particules: une perspective mathematique et numerique

Energy Technology Data Exchange (ETDEWEB)

Lewin, M.

2009-06-15

This thesis is devoted to the mathematical study of variational models for large quantum systems. The mathematical methods are that of nonlinear analysis, calculus of variations, partial differential equations, spectral theory, and numerical analysis. The first part contains some results on finite systems. We study several approximations of the N-body Schroedinger equation for electrons in an atom or a molecule, and then the so-called Hartree-Fock- Bogoliubov model for a system of fermions interacting via the gravitational force. In a second part, we propose a new method allowing to prove the existence of the thermodynamic limit of Coulomb quantum systems. Then, we construct two Hartree-Fock-type models for infinite systems. The first is a relativistic theory deduced from Quantum Electrodynamics, allowing to describe the behavior of electrons, coupled to that of Dirac's vacuum which can become polarized. The second model describes a nonrelativistic quantum crystal in the presence of a charged defect. A new numerical method is also proposed. The last part of the thesis is devoted to spectral pollution, a phenomenon which is observed when trying to approximate eigenvalues in a gap of the essential spectrum of a self-adjoint operator, for instance for periodic Schroedinger operators or Dirac operators. (author)

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

Analytical-numerical solution of a nonlinear integrodifferential equation in econometrics

Science.gov (United States)

Kakhktsyan, V. M.; Khachatryan, A. Kh.

2013-07-01

A mixed problem for a nonlinear integrodifferential equation arising in econometrics is considered. An analytical-numerical method is proposed for solving the problem. Some numerical results are presented.

Numerical simulation of fractional Cable equation of spiny neuronal dendrites

Directory of Open Access Journals (Sweden)

N.H. Sweilam

2014-03-01

Full Text Available In this article, numerical study for the fractional Cable equation which is fundamental equations for modeling neuronal dynamics is introduced by using weighted average of finite difference methods. The stability analysis of the proposed methods is given by a recently proposed procedure similar to the standard John von Neumann stability analysis. A simple and an accurate stability criterion valid for different discretization schemes of the fractional derivative and arbitrary weight factor is introduced and checked numerically. Numerical results, figures, and comparisons have been presented to confirm the theoretical results and efficiency of the proposed method.

Numerical magnitude processing in abacus-trained children with superior mathematical ability: an EEG study.

Science.gov (United States)

Huang, Jian; Du, Feng-lei; Yao, Yuan; Wan, Qun; Wang, Xiao-Song; Chen, Fei-Yan

2015-08-01

Distance effect has been regarded as the best established marker of basic numerical magnitude processes and is related to individual mathematical abilities. A larger behavioral distance effect is suggested to be concomitant with lower mathematical achievement in children. However, the relationship between distance effect and superior mathematical abilities is unclear. One could get superior mathematical abilities by acquiring the skill of abacus-based mental calculation (AMC), which can be used to solve calculation problems with exceptional speed and high accuracy. In the current study, we explore the relationship between distance effect and superior mathematical abilities by examining whether and how the AMC training modifies numerical magnitude processing. Thus, mathematical competencies were tested in 18 abacus-trained children (who accepted the AMC training) and 18 non-trained children. Electroencephalography (EEG) waveforms were recorded when these children executed numerical comparison tasks in both Arabic digit and dot array forms. We found that: (a) the abacus-trained group had superior mathematical abilities than their peers; (b) distance effects were found both in behavioral results and on EEG waveforms; (c) the distance effect size of the average amplitude on the late negative-going component was different between groups in the digit task, with a larger effect size for abacus-trained children; (d) both the behavioral and EEG distance effects were modulated by the notation. These results revealed that the neural substrates of magnitude processing were modified by AMC training, and suggested that the mechanism of the representation of numerical magnitude for children with superior mathematical abilities was different from their peers. In addition, the results provide evidence for a view of non-abstract numerical representation.

Research in applied mathematics, numerical analysis, and computer science

Science.gov (United States)

1984-01-01

Research conducted at the Institute for Computer Applications in Science and Engineering (ICASE) in applied mathematics, numerical analysis, and computer science is summarized and abstracts of published reports are presented. The major categories of the ICASE research program are: (1) numerical methods, with particular emphasis on the development and analysis of basic numerical algorithms; (2) control and parameter identification; (3) computational problems in engineering and the physical sciences, particularly fluid dynamics, acoustics, and structural analysis; and (4) computer systems and software, especially vector and parallel computers.

Research in progress in applied mathematics, numerical analysis, and computer science

Science.gov (United States)

1990-01-01

Research conducted at the Institute in Science and Engineering in applied mathematics, numerical analysis, and computer science is summarized. The Institute conducts unclassified basic research in applied mathematics in order to extend and improve problem solving capabilities in science and engineering, particularly in aeronautics and space.

Geophysical interpretation using integral equations

CERN Document Server

Eskola, L

1992-01-01

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

Uncertain differential equations

CERN Document Server

Yao, Kai

2016-01-01

This book introduces readers to the basic concepts of and latest findings in the area of differential equations with uncertain factors. It covers the analytic method and numerical method for solving uncertain differential equations, as well as their applications in the field of finance. Furthermore, the book provides a number of new potential research directions for uncertain differential equation. It will be of interest to researchers, engineers and students in the fields of mathematics, information science, operations research, industrial engineering, computer science, artificial intelligence, automation, economics, and management science.

Stability of Functional Differential Equations

CERN Document Server

Lemm, Jeffrey M

1986-01-01

This book provides an introduction to the structure and stability properties of solutions of functional differential equations. Numerous examples of applications (such as feedback systrems with aftereffect, two-reflector antennae, nuclear reactors, mathematical models in immunology, viscoelastic bodies, aeroautoelastic phenomena and so on) are considered in detail. The development is illustrated by numerous figures and tables.

Numerical solution of plasma fluid equations using locally refined grids

International Nuclear Information System (INIS)

Colella, P.

1997-01-01

This paper describes a numerical method for the solution of plasma fluid equations on block-structured, locally refined grids. The plasma under consideration is typical of those used for the processing of semiconductors. The governing equations consist of a drift-diffusion model of the electrons and an isothermal model of the ions coupled by Poisson's equation. A discretization of the equations is given for a uniform spatial grid, and a time-split integration scheme is developed. The algorithm is then extended to accommodate locally refined grids. This extension involves the advancement of the discrete system on a hierarchy of levels, each of which represents a degree of refinement, together with synchronization steps to ensure consistency across levels. A brief discussion of a software implementation is followed by a presentation of numerical results

Stable Numerical Approach for Fractional Delay Differential Equations

Science.gov (United States)

Singh, Harendra; Pandey, Rajesh K.; Baleanu, D.

2017-12-01

In this paper, we present a new stable numerical approach based on the operational matrix of integration of Jacobi polynomials for solving fractional delay differential equations (FDDEs). The operational matrix approach converts the FDDE into a system of linear equations, and hence the numerical solution is obtained by solving the linear system. The error analysis of the proposed method is also established. Further, a comparative study of the approximate solutions is provided for the test examples of the FDDE by varying the values of the parameters in the Jacobi polynomials. As in special case, the Jacobi polynomials reduce to the well-known polynomials such as (1) Legendre polynomial, (2) Chebyshev polynomial of second kind, (3) Chebyshev polynomial of third and (4) Chebyshev polynomial of fourth kind respectively. Maximum absolute error and root mean square error are calculated for the illustrated examples and presented in form of tables for the comparison purpose. Numerical stability of the presented method with respect to all four kind of polynomials are discussed. Further, the obtained numerical results are compared with some known methods from the literature and it is observed that obtained results from the proposed method is better than these methods.

The Numerical Solution of an Abelian Ordinary Differential Equation ...

African Journals Online (AJOL)

In this paper we present a relatively new technique call theNew Hybrid of Adomian decomposition method (ADM) for solution of an Abelian Differential equation. The numerical results of the equation have been obtained in terms of convergent series with easily computable component. These methods are applied to solve ...

Brain stimulation, mathematical, and numerical training: Contribution of core and noncore skills.

Science.gov (United States)

Looi, C Y; Cohen Kadosh, R

2016-01-01

Mathematical abilities that are correlated with various life outcomes vary across individuals. One approach to improve mathematical abilities is by understanding the underlying cognitive functions. Theoretical and experimental evidence suggest that mathematical abilities are subserved by "core" and "noncore" skills. Core skills are commonly regarded as the "innate" capacity to attend to and process numerical information, while noncore skills are those that are important for mathematical cognition, but are not exclusive to the mathematical domain such as executive functions, spatial skills, and attention. In recent years, mathematical training has been combined with the application of noninvasive brain stimulation to further enhance training outcomes. However, the development of more strategic training paradigms is hindered by the lack of understanding on the contributory nature of core and noncore skills and their neural underpinnings. In the current review, we will examine the effects of brain stimulation with focus on transcranial electrical stimulation on core and noncore skills, and its impact on mathematical and numerical training. We will conclude with a discussion on the theoretical and experimental implications of these studies and directions for further research. © 2016 Elsevier B.V. All rights reserved.

Computational mathematics in China

CERN Document Server

Shi, Zhong-Ci

1994-01-01

This volume describes the most significant contributions made by Chinese mathematicians over the past decades in various areas of computational mathematics. Some of the results are quite important and complement Western developments in the field. The contributors to the volume range from noted senior mathematicians to promising young researchers. The topics include finite element methods, computational fluid mechanics, numerical solutions of differential equations, computational methods in dynamical systems, numerical algebra, approximation, and optimization. Containing a number of survey articles, the book provides an excellent way for Western readers to gain an understanding of the status and trends of computational mathematics in China.

Numerical approximation abilities correlate with and predict informal but not formal mathematics abilities.

Science.gov (United States)

Libertus, Melissa E; Feigenson, Lisa; Halberda, Justin

2013-12-01

Previous research has found a relationship between individual differences in children's precision when nonverbally approximating quantities and their school mathematics performance. School mathematics performance emerges from both informal (e.g., counting) and formal (e.g., knowledge of mathematics facts) abilities. It remains unknown whether approximation precision relates to both of these types of mathematics abilities. In the current study, we assessed the precision of numerical approximation in 85 3- to 7-year-old children four times over a span of 2years. In addition, at the final time point, we tested children's informal and formal mathematics abilities using the Test of Early Mathematics Ability (TEMA-3). We found that children's numerical approximation precision correlated with and predicted their informal, but not formal, mathematics abilities when controlling for age and IQ. These results add to our growing understanding of the relationship between an unlearned nonsymbolic system of quantity representation and the system of mathematics reasoning that children come to master through instruction. Copyright © 2013 Elsevier Inc. All rights reserved.

Mathematical and numerical modelling of fluids at Nano-metric scales

International Nuclear Information System (INIS)

Joubaud, R.

2012-01-01

This work presents some contributions to the mathematical and numerical modelling of fluids at Nano-metric scales. We are interested in two levels of modelling. The first level consists in an atomic description. We consider the problem of computing the shear viscosity of a fluid from a microscopic description. More precisely, we study the mathematical properties of the nonequilibrium Langevin dynamics allowing to compute the shear viscosity. The second level of description is a continuous description, and we consider a class of continuous models for equilibrium electrolytes, which incorporate on the one hand a confinement by charged solid objects and on the other hand non-ideality effects stemming from electrostatic correlations and steric exclusion phenomena due to the excluded volume effects. First, we perform the mathematical analysis of the case where the free energy is a convex function (mild non-ideality). Second, we consider numerically the case where the free energy is a non convex function (strong non-ideality) leading in particular to phase separation. (author)

Sparse grid spectral methods for the numerical solution of partial differential equations with periodic boundary conditions

International Nuclear Information System (INIS)

Kupka, F.

1997-11-01

This thesis deals with the extension of sparse grid techniques to spectral methods for the solution of partial differential equations with periodic boundary conditions. A review on boundary and initial-boundary value problems and a discussion on numerical resolution is used to motivate this research. Spectral methods are introduced by projection techniques, and by three model problems: the stationary and the transient Helmholtz equations, and the linear advection equation. The approximation theory on the hyperbolic cross is reviewed and its close relation to sparse grids is demonstrated. This approach extends to non-periodic problems. Various Sobolev spaces with dominant mixed derivative are introduced to provide error estimates for Fourier approximation and interpolation on the hyperbolic cross and on sparse grids by means of Sobolev norms. The theorems are immediately applicable to the stability and convergence analysis of sparse grid spectral methods. This is explicitly demonstrated for the three model problems. A variant of the von Neumann condition is introduced to simplify the stability analysis of the time-dependent model problems. The discrete Fourier transformation on sparse grids is discussed together with its software implementation. Results on numerical experiments are used to illustrate the performance of the new method with respect to the smoothness properties of each example. The potential of the method in mathematical modelling is estimated and generalizations to other sparse grid methods are suggested. The appendix includes a complete Fortran90 program to solve the linear advection equation by the sparse grid Fourier collocation method and a third-order Runge-Kutta routine for integration in time. (author)

Numerical solutions of ordinary and partial differential equations in the frequency domain

International Nuclear Information System (INIS)

Hazi, G.; Por, G.

1997-01-01

Numerical problems during the noise simulation in a nuclear power plant are discussed. The solutions of ordinary and partial differential equations are studied in the frequency domain. Numerical methods by the transfer function method are applied. It is shown that the correctness of the numerical methods is limited for ordinary differential equations in the frequency domain. To overcome the difficulties, step-size selection is suggested. (author)

An Adjoint-based Numerical Method for a class of nonlinear Fokker-Planck Equations

KAUST Repository

Festa, Adriano; Gomes, Diogo A.; Machado Velho, Roberto

2017-01-01

Here, we introduce a numerical approach for a class of Fokker-Planck (FP) equations. These equations are the adjoint of the linearization of Hamilton-Jacobi (HJ) equations. Using this structure, we show how to transfer the properties of schemes for HJ equations to the FP equations. Hence, we get numerical schemes with desirable features such as positivity and mass-preservation. We illustrate this approach in examples that include mean-field games and a crowd motion model.

An Adjoint-based Numerical Method for a class of nonlinear Fokker-Planck Equations

KAUST Repository

Festa, Adriano

2017-03-22

Here, we introduce a numerical approach for a class of Fokker-Planck (FP) equations. These equations are the adjoint of the linearization of Hamilton-Jacobi (HJ) equations. Using this structure, we show how to transfer the properties of schemes for HJ equations to the FP equations. Hence, we get numerical schemes with desirable features such as positivity and mass-preservation. We illustrate this approach in examples that include mean-field games and a crowd motion model.

A course in mathematical methods for physicists

CERN Document Server

Herman, Russell L

2014-01-01

Based on the author’s junior-level undergraduate course, this introductory textbook is designed for a course in mathematical physics. Focusing on the physics of oscillations and waves, A Course in Mathematical Methods for Physicists helps students understand the mathematical techniques needed for their future studies in physics. It takes a bottom-up approach that emphasizes physical applications of the mathematics. The book offers: •A quick review of mathematical prerequisites, proceeding to applications of differential equations and linear algebra •Classroom-tested explanations of complex and Fourier analysis for trigonometric and special functions •Coverage of vector analysis and curvilinear coordinates for solving higher dimensional problems •Sections on nonlinear dynamics, variational calculus, numerical solutions of differential equations, and Green's functions

on the properties of solutions and some applications on the TOV differential equation with a model of nuclear equation of state

International Nuclear Information System (INIS)

Esmail, S.F.H.

2006-01-01

the mathematical formulation of numerous physical problems results in differential equations actually non-linear differential equations . in our study we are interested in solutions of differential equations which describe the structure of neutron star in non-relativistic and relativistic cases. the aim of this work is to determine the mass and the radius of a neutron star, by solving the tolmann-oppenheimer-volkoff (TOV) differential equation using different models of the nuclear equation of state (EOS). analytically solutions are obtained for a simple form of the nuclear equation of state of Clayton model and poly trope model. for a more realistic equation of state the TOV differential equation is solved numerically using rung -Kutta method

Mathematical well-posedness of a two-fluid equations for bubbly two-phase flows

International Nuclear Information System (INIS)

Okawa, Tomio; Kataoka, Isao

2000-01-01

It is widely known that two-fluid equations used in most engineering applications do not satisfy the necessary condition for being mathematical well-posed as initial-value problems. In the case of stratified two-phase flows, several researchers have revealed that differential models satisfying the necessary condition are to be derived if the pressure difference between the phases is related to the spatial gradient of the void fraction through the effects of gravity or surface tension. While, in the case of dispersed two-phase flows, no physically reasonable method to derive mathematically well-posed two-fluid model has been proposed. In the present study, particularly focusing on the effect of interfacial pressure terms, we derived the mathematically closed form of the volume-averaged two-fluid model for bubbly two-phase flows. As a result of characteristic analyses, it was shown that the proposed two-fluid equations satisfy the necessary condition of mathematical well-posedness if the void fraction is sufficiently small. (author)

Computational mathematics and mathematical computer software. Vychislitel'naia matematika i matematicheskoe obespechenie EVM

Energy Technology Data Exchange (ETDEWEB)

Tikhonov, A.N.; Samarskii, A.A.

1985-01-01

Various aspects of mathematical modeling and problem-oriented computer software are examined with reference to numerical methods in mathematical physics, methods for solving inverse problems, development of automatic systems for experimental data processing, and mathematical modeling in plasma physics. Papers are presented on some properties of difference schemes in one-dimensional gas dynamics, an algorithm for processing signals reflected from multipoint targets, and the application of simplified Navier-Stokes equations for calculating flow of a viscous gas past long bodies.

CSR Fields: Direct Numerical Solution of the Maxwell's Equation

International Nuclear Information System (INIS)

Novokhatski, Alexander

2011-01-01

We discuss the properties of the coherent electromagnetic fields of a very short, ultra-relativistic bunch in a rectangular vacuum chamber inside a bending magnet. The analysis is based on the results of a direct numerical solution of Maxwell's equations together with Newton's equations. We use a new dispersion-free time-domain algorithm which employs a more efficient use of finite element mesh techniques and hence produces self-consistent and stable solutions for very short bunches. We investigate the fine structure of the CSR fields including coherent edge radiation. This approach should be useful in the study of existing and future concepts of particle accelerators and ultrafast coherent light sources. The coherent synchrotron radiation (CSR) fields have a strong action on the beam dynamics of very short bunches, which are moving in the bends of all kinds of magnetic elements. They are responsible for additional energy loss and energy spread; micro bunching and beam emittance growth. These fields may bound the efficiency of damping rings, electron-positron colliders and ultrafast coherent light sources, where high peak currents and very short bunches are envisioned. This is relevant to most high-brightness beam applications. On the other hand these fields together with transition radiation fields can be used for beam diagnostics or even as a powerful resource of THz radiation. A history of the study of CSR and a good collection of references can be found in (1). Electromagnetic theory suggests several methods on how to calculate CSR fields. The most popular method is to use Lienard-Wiechert potentials. Other approach is to solve numerically the approximate equations, which are a Schrodinger type equation. These numerical methods are described in (2). We suggest that a direct solution of Maxwell's equations together with Newton's equations can describe the detailed structure of the CSR fields (3).

Numerical analysis

CERN Document Server

Brezinski, C

2012-01-01

Numerical analysis has witnessed many significant developments in the 20th century. This book brings together 16 papers dealing with historical developments, survey papers and papers on recent trends in selected areas of numerical analysis, such as: approximation and interpolation, solution of linear systems and eigenvalue problems, iterative methods, quadrature rules, solution of ordinary-, partial- and integral equations. The papers are reprinted from the 7-volume project of the Journal of Computational and Applied Mathematics on '/homepage/sac/cam/na2000/index.html<

A mathematical solution for the parameters of three interfering resonances

Science.gov (United States)

Han, X.; Shen, C. P.

2018-04-01

The multiple-solution problem in determining the parameters of three interfering resonances from a fit to an experimentally measured distribution is considered from a mathematical viewpoint. It is shown that there are four numerical solutions for a fit with three coherent Breit-Wigner functions. Although explicit analytical formulae cannot be derived in this case, we provide some constraint equations between the four solutions. For the cases of nonrelativistic and relativistic Breit-Wigner forms of amplitude functions, a numerical method is provided to derive the other solutions from that already obtained, based on the obtained constraint equations. In real experimental measurements with more complicated amplitude forms similar to Breit-Wigner functions, the same method can be deduced and performed to get numerical solutions. The good agreement between the solutions found using this mathematical method and those directly from the fit verifies the correctness of the constraint equations and mathematical methodology used. Supported by National Natural Science Foundation of China (NSFC) (11575017, 11761141009), the Ministry of Science and Technology of China (2015CB856701) and the CAS Center for Excellence in Particle Physics (CCEPP)

Mathematical models and numerical simulation in electromagnetism

CERN Document Server

Bermúdez, Alfredo; Salgado, Pilar

2014-01-01

The book represents a basic support for a master course in electromagnetism oriented to numerical simulation. The main goal of the book is that the reader knows the boundary-value problems of partial differential equations that should be solved in order to perform computer simulation of electromagnetic processes. Moreover it includes a part devoted to electric circuit theory  based on ordinary differential equations. The book is mainly oriented to electric engineering applications, going from the general to the specific, namely, from the full Maxwell’s equations to the particular cases of electrostatics, direct current, magnetostatics and eddy currents models. Apart from standard exercises related to analytical calculus, the book includes some others oriented to real-life applications solved with MaxFEM free simulation software.

Mathematical Modelling of Intraretinal Oxygen Partial Pressure

African Journals Online (AJOL)

Erah

The system of non-linear differential equations was solved numerically using Runge-kutta. Nystroms method. ... artery occlusion. Keywords: Mathematical modeling, Intraretinal oxygen pressure, Retinal capillaries, Oxygen ..... Mass transfer,.

Equations in mathematical physics a practical course

CERN Document Server

Pikulin, Victor P

2001-01-01

Many physical processes in fields such as mechanics, thermodynamics, electricity, magnetism or optics are described by means of partial differential equations. The aim of the present book is to demonstrate the basic methods for solving the classical linear problems in mathematical physics of elliptic, parabolic and hyperbolic type. In particular, the methods of conformal mappings, Fourier analysis and Green`s functions are considered, as well as the perturbation method and integral transformation method, among others. Every chapter contains concrete examples with a detailed analysis of their solution. The book is intended as a textbook for students in mathematical physics, but will also serve as a handbook for scientists and engineers.   ------------ [A] manual for future engineers must strongly differ from the textbook for pure mathematicians, and the book by Pikulin and Pohozaev is the good example. (…) The purpose (…)  is to offer quick access to the principal facts (…) This well written book is a...

Numerical bifurcation analysis of a class of nonlinear renewal equations

NARCIS (Netherlands)

Breda, Dimitri; Diekmann, Odo; Liessi, Davide; Scarabel, Francesca

2016-01-01

We show, by way of an example, that numerical bifurcation tools for ODE yield reliable bifurcation diagrams when applied to the pseudospectral approximation of a one-parameter family of nonlinear renewal equations. The example resembles logistic-and Ricker-type population equations and exhibits

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.

Numerical solutions of multi-order fractional differential equations by Boubaker polynomials

Directory of Open Access Journals (Sweden)

Bolandtalat A.

2016-01-01

Full Text Available In this paper, we have applied a numerical method based on Boubaker polynomials to obtain approximate numerical solutions of multi-order fractional differential equations. We obtain an operational matrix of fractional integration based on Boubaker polynomials. Using this operational matrix, the given problem is converted into a set of algebraic equations. Illustrative examples are are given to demonstrate the efficiency and simplicity of this technique.

A Numerical Simulation for a Deterministic Compartmental ...

African Journals Online (AJOL)

In this work, an earlier deterministic mathematical model of HIV/AIDS is revisited and numerical solutions obtained using Eulers numerical method. Using hypothetical values for the parameters, a program was written in VISUAL BASIC programming language to generate series for the system of difference equations from the ...

International Winter Workshop on Differential Equations and Numerical Analysis

CERN Document Server

Miller, John; Narasimhan, Ramanujam; Mathiazhagan, Paramasivam; Victor, Franklin

2016-01-01

This book offers an ideal introduction to singular perturbation problems, and a valuable guide for researchers in the field of differential equations. It also includes chapters on new contributions to both fields: differential equations and singular perturbation problems. Written by experts who are active researchers in the related fields, the book serves as a comprehensive source of information on the underlying ideas in the construction of numerical methods to address different classes of problems with solutions of different behaviors, which will ultimately help researchers to design and assess numerical methods for solving new problems. All the chapters presented in the volume are complemented by illustrations in the form of tables and graphs.

Physical mathematics

CERN Document Server

Cahill, Kevin

2013-01-01

Unique in its clarity, examples and range, Physical Mathematics explains as simply as possible the mathematics that graduate students and professional physicists need in their courses and research. The author illustrates the mathematics with numerous physical examples drawn from contemporary research. In addition to basic subjects such as linear algebra, Fourier analysis, complex variables, differential equations and Bessel functions, this textbook covers topics such as the singular-value decomposition, Lie algebras, the tensors and forms of general relativity, the central limit theorem and Kolmogorov test of statistics, the Monte Carlo methods of experimental and theoretical physics, the renormalization group of condensed-matter physics and the functional derivatives and Feynman path integrals of quantum field theory.

Differential equations for dummies

CERN Document Server

Holzner, Steven

2008-01-01

The fun and easy way to understand and solve complex equations Many of the fundamental laws of physics, chemistry, biology, and economics can be formulated as differential equations. This plain-English guide explores the many applications of this mathematical tool and shows how differential equations can help us understand the world around us. Differential Equations For Dummies is the perfect companion for a college differential equations course and is an ideal supplemental resource for other calculus classes as well as science and engineering courses. It offers step-by-step techniques, practical tips, numerous exercises, and clear, concise examples to help readers improve their differential equation-solving skills and boost their test scores.

Conservative numerical schemes for Euler-Lagrange equations

Energy Technology Data Exchange (ETDEWEB)

Vazquez, L. [Universidad Complutense, Madrid (Spain). Dept. de Matematica Aplicada; Jimenez, S. [Universidad Alfonso X El Sabio, Madrid (Spain). Dept. de Matematica Aplicada

1999-05-01

As a preliminary step to study magnetic field lines, the authors seek numerical schemes that reproduce at discrete level the significant feature of the continuous model, based on an underling Lagrangian structure. The resulting scheme give discrete counterparts of the variation law for the energy as well of as the Euler-Lagrange equations and their symmetries.

Numerical stability in problems of linear algebra.

Science.gov (United States)

Babuska, I.

1972-01-01

Mathematical problems are introduced as mappings from the space of input data to that of the desired output information. Then a numerical process is defined as a prescribed recurrence of elementary operations creating the mapping of the underlying mathematical problem. The ratio of the error committed by executing the operations of the numerical process (the roundoff errors) to the error introduced by perturbations of the input data (initial error) gives rise to the concept of lambda-stability. As examples, several processes are analyzed from this point of view, including, especially, old and new processes for solving systems of linear algebraic equations with tridiagonal matrices. In particular, it is shown how such a priori information can be utilized as, for instance, a knowledge of the row sums of the matrix. Information of this type is frequently available where the system arises in connection with the numerical solution of differential equations.

Mathematic study and numerical implantation of the Vlasov-Darwin model

International Nuclear Information System (INIS)

Sonnendrucker, E.

1994-12-01

Numerical simulation of some phenomena in plasma physics, or more generally in electromagnetism, can be more easily done using approximate models of Maxwell equations such as the Darwin model in which the transverse part of the displacement current in the Ampere equation is neglected, or such as the static model in which the time derivatives are neglected. In this note, the Darwin model is presented first, and then an asymptotic analysis of Maxwell equations is given with limit conditions of perfect conductor on one part of the side, and Silver-Muller absorbing conditions on the other part. This allows to obtain a variational formulation for the Darwin model which is a good approximation of Maxwell equations. A variational formulation for the quasi-static model is also obtained. In a second part this implantation is described using a 2-D finite element method coupled with a particulate method for the Vlasov equations which leads to numerical results allowing a determination of the different models application. (J.S.). 2 refs

Numerical solution of modified differential equations based on symmetry preservation.

Science.gov (United States)

Ozbenli, Ersin; Vedula, Prakash

2017-12-01

In this paper, we propose a method to construct invariant finite-difference schemes for solution of partial differential equations (PDEs) via consideration of modified forms of the underlying PDEs. The invariant schemes, which preserve Lie symmetries, are obtained based on the method of equivariant moving frames. While it is often difficult to construct invariant numerical schemes for PDEs due to complicated symmetry groups associated with cumbersome discrete variable transformations, we note that symmetries associated with more convenient transformations can often be obtained by appropriately modifying the original PDEs. In some cases, modifications to the original PDEs are also found to be useful in order to avoid trivial solutions that might arise from particular selections of moving frames. In our proposed method, modified forms of PDEs can be obtained either by addition of perturbation terms to the original PDEs or through defect correction procedures. These additional terms, whose primary purpose is to enable symmetries with more convenient transformations, are then removed from the system by considering moving frames for which these specific terms go to zero. Further, we explore selection of appropriate moving frames that result in improvement in accuracy of invariant numerical schemes based on modified PDEs. The proposed method is tested using the linear advection equation (in one- and two-dimensions) and the inviscid Burgers' equation. Results obtained for these tests cases indicate that numerical schemes derived from the proposed method perform significantly better than existing schemes not only by virtue of improvement in numerical accuracy but also due to preservation of qualitative properties or symmetries of the underlying differential equations.

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.

Consistent three-equation model for thin films

Science.gov (United States)

Richard, Gael; Gisclon, Marguerite; Ruyer-Quil, Christian; Vila, Jean-Paul

2017-11-01

Numerical simulations of thin films of newtonian fluids down an inclined plane use reduced models for computational cost reasons. These models are usually derived by averaging over the fluid depth the physical equations of fluid mechanics with an asymptotic method in the long-wave limit. Two-equation models are based on the mass conservation equation and either on the momentum balance equation or on the work-energy theorem. We show that there is no two-equation model that is both consistent and theoretically coherent and that a third variable and a three-equation model are required to solve all theoretical contradictions. The linear and nonlinear properties of two and three-equation models are tested on various practical problems. We present a new consistent three-equation model with a simple mathematical structure which allows an easy and reliable numerical resolution. The numerical calculations agree fairly well with experimental measurements or with direct numerical resolutions for neutral stability curves, speed of kinematic waves and of solitary waves and depth profiles of wavy films. The model can also predict the flow reversal at the first capillary trough ahead of the main wave hump.

Validation of a numerical algorithm based on transformed equations

International Nuclear Information System (INIS)

Xu, H.; Barron, R.M.; Zhang, C.

2003-01-01

Generally, a typical equation governing a physical process, such as fluid flow or heat transfer, has three types of terms that involve partial derivatives, namely, the transient term, the convective terms and the diffusion terms. The major difficulty in obtaining numerical solutions of these partial differential equations is the discretization of the convective terms. The transient term is usually discretized using the first-order forward or backward differencing scheme. The diffusion terms are usually discretized using the central differencing scheme and no difficulty arises since these terms involve second-order spatial derivatives of the flow variables. The convective terms are non-linear and contain first-order spatial derivatives. The main difference between various numerical algorithms is the discretization of the convective terms. In the present study, an alternative approach to discretizing the governing equations is presented. In this algorithm, the governing equations are first transformed by introducing an exponential function to eliminate the convective terms in the equations. The proposed algorithm is applied to simulate some fluid flows with exact solutions to validate the proposed algorithm. The fluid flows used in this study are a self-designed quasi-fluid flow problem, stagnation in plane flow (Hiemenz flow), and flow between two concentric cylinders. The comparisons with the power-law scheme indicate that the proposed scheme exhibits better performance. (author)

Effects of finger counting on numerical development – the opposing views of neurocognition and mathematics education

Directory of Open Access Journals (Sweden)

Korbinian eMoeller

2011-11-01

Full Text Available Usually children learn the basic principles of number and arithmetic by the help of finger-based representations. However, whether the reliance on finger-based representations is only beneficial or whether it may even become detrimental is the subject of an ongoing debate between neuro-cognitive and mathematics education researchers. From the neuro-cognitive perspective finger counting provides multi-sensory input conveying both cardinal and ordinal aspects of numbers. Recent data indicate that children with good finger-based numerical representations show better arithmetic skills and that training finger gnosis enhances mathematical skills. From this neuro-cognitive researchers conclude that elaborate finger-based numerical representations are beneficial for later numerical development.However, mathematics education research recommends fostering mental numerical representations so as to induce children to abandon finger-counting. More precisely mathematics education recommends moving from finger counting to concrete structured representations and then, finally, to mental representations of numbers.Taken together, there is obviously an important debate between the neuro-cognitve and mathematics education research concerning the benefits or detriments of finger-based strategies for numerical development. In the present review, the rationale of both lines of evidence will be presented and discussed.

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.

Theoretical and numerical method in aeroacoustics

Directory of Open Access Journals (Sweden)

Nicuşor ALEXANDRESCU

2010-06-01

Full Text Available The paper deals with the mathematical and numerical modeling of the aerodynamic noisegenerated by the fluid flow interaction with the solid structure of a rotor blade.Our analysis use Lighthill’s acoustic analogy. Lighthill idea was to express the fundamental equationsof motion into a wave equation for acoustic fluctuation with a source term on the right-hand side. Theobtained wave equation is solved numerically by the spatial discretization. The method is applied inthe case of monopole source placed in different points of blade surfaces to find this effect of noisepropagation.

Simple Numerical Schemes for the Korteweg-deVries Equation

International Nuclear Information System (INIS)

McKinstrie, C. J.; Kozlov, M.V.

2000-01-01

Two numerical schemes, which simulate the propagation of dispersive non-linear waves, are described. The first is a split-step Fourier scheme for the Korteweg-de Vries (KdV) equation. The second is a finite-difference scheme for the modified KdV equation. The stability and accuracy of both schemes are discussed. These simple schemes can be used to study a wide variety of physical processes that involve dispersive nonlinear waves

Simple Numerical Schemes for the Korteweg-deVries Equation

Energy Technology Data Exchange (ETDEWEB)

C. J. McKinstrie; M. V. Kozlov

2000-12-01

Two numerical schemes, which simulate the propagation of dispersive non-linear waves, are described. The first is a split-step Fourier scheme for the Korteweg-de Vries (KdV) equation. The second is a finite-difference scheme for the modified KdV equation. The stability and accuracy of both schemes are discussed. These simple schemes can be used to study a wide variety of physical processes that involve dispersive nonlinear waves.

Teaching Mathematical Modelling for Earth Sciences via Case Studies

Science.gov (United States)

Yang, Xin-She

2010-05-01

Mathematical modelling is becoming crucially important for earth sciences because the modelling of complex systems such as geological, geophysical and environmental processes requires mathematical analysis, numerical methods and computer programming. However, a substantial fraction of earth science undergraduates and graduates may not have sufficient skills in mathematical modelling, which is due to either limited mathematical training or lack of appropriate mathematical textbooks for self-study. In this paper, we described a detailed case-study-based approach for teaching mathematical modelling. We illustrate how essential mathematical skills can be developed for students with limited training in secondary mathematics so that they are confident in dealing with real-world mathematical modelling at university level. We have chosen various topics such as Airy isostasy, greenhouse effect, sedimentation and Stokes' flow,free-air and Bouguer gravity, Brownian motion, rain-drop dynamics, impact cratering, heat conduction and cooling of the lithosphere as case studies; and we use these step-by-step case studies to teach exponentials, logarithms, spherical geometry, basic calculus, complex numbers, Fourier transforms, ordinary differential equations, vectors and matrix algebra, partial differential equations, geostatistics and basic numeric methods. Implications for teaching university mathematics for earth scientists for tomorrow's classroom will also be discussed. Refereces 1) D. L. Turcotte and G. Schubert, Geodynamics, 2nd Edition, Cambridge University Press, (2002). 2) X. S. Yang, Introductory Mathematics for Earth Scientists, Dunedin Academic Press, (2009).

Lectures on the Numerical Solution of Partial Differential Equations.

Science.gov (United States)

1981-12-01

Mathematics Rockefeller University Professor F. Brezzi New York, New York 10021 Laboratorio di Analisi Numerica Universita di Pavia Professor Amiram Harten...equations and to control the spacing of the points sj. In the MFE process the grid points move with the solution and cluster atound areas of roughness...149-159. [28] Fichera, G.: Analisi essistenziale per le soluzioni die problemi al contorno misti relativi alle equazione ed ai sistemi di equazioni

Random ordinary differential equations and their numerical solution

CERN Document Server

Han, Xiaoying

2017-01-01

This book is intended to make recent results on the derivation of higher order numerical schemes for random ordinary differential equations (RODEs) available to a broader readership, and to familiarize readers with RODEs themselves as well as the closely associated theory of random dynamical systems. In addition, it demonstrates how RODEs are being used in the biological sciences, where non-Gaussian and bounded noise are often more realistic than the Gaussian white noise in stochastic differential equations (SODEs).   RODEs are used in many important applications and play a fundamental role in the theory of random dynamical systems.  They can be analyzed pathwise with deterministic calculus, but require further treatment beyond that of classical ODE theory due to the lack of smoothness in their time variable. Although classical numerical schemes for ODEs can be used pathwise for RODEs, they rarely attain their traditional order since the solutions of RODEs do not have sufficient smoothness to have Taylor ...

Theoretical stability in coefficient inverse problems for general hyperbolic equations with numerical reconstruction

Science.gov (United States)

Yu, Jie; Liu, Yikan; Yamamoto, Masahiro

2018-04-01

In this article, we investigate the determination of the spatial component in the time-dependent second order coefficient of a hyperbolic equation from both theoretical and numerical aspects. By the Carleman estimates for general hyperbolic operators and an auxiliary Carleman estimate, we establish local Hölder stability with either partial boundary or interior measurements under certain geometrical conditions. For numerical reconstruction, we minimize a Tikhonov functional which penalizes the gradient of the unknown function. Based on the resulting variational equation, we design an iteration method which is updated by solving a Poisson equation at each step. One-dimensional prototype examples illustrate the numerical performance of the proposed iteration.

A asymptotic numerical method for the steady-state convection diffusion equation

International Nuclear Information System (INIS)

Wu Qiguang

1988-01-01

In this paper, A asymptotic numerical method for the steady-state Convection diffusion equation is proposed, which need not take very fine mesh size in the neighbourhood of the boundary layer. Numerical computation for model problem show that we can obtain the numerical solution in the boundary layer with moderate step size

Mathematics for electronic technology

CERN Document Server

Howson, D P

1975-01-01

Mathematics for Electronic Technology is a nine-chapter book that begins with the elucidation of the introductory concepts related to use of mathematics in electronic engineering, including differentiation, integration, partial differentiation, infinite series, vectors, vector algebra, and surface, volume and line integrals. Subsequent chapters explore the determinants, differential equations, matrix analysis, complex variable, topography, graph theory, and numerical analysis used in this field. The use of Fourier method for harmonic analysis and the Laplace transform is also described. The ma

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)

Attempts at a numerical realisation of stochastic differential equations containing Preisach operator

International Nuclear Information System (INIS)

McCarthy, S; Rachinskii, D

2011-01-01

We describe two Euler type numerical schemes obtained by discretisation of a stochastic differential equation which contains the Preisach memory operator. Equations of this type are of interest in areas such as macroeconomics and terrestrial hydrology where deterministic models containing the Preisach operator have been developed but do not fully encapsulate stochastic aspects of the area. A simple price dynamics model is presented as one motivating example for our studies. Some numerical evidence is given that the two numerical schemes converge to the same limit as the time step decreases. We show that the Preisach term introduces a damping effect which increases on the parts of the trajectory demonstrating a stronger upwards or downwards trend. The results are preliminary to a broader programme of research of stochastic differential equations with the Preisach hysteresis operator.

The Effect of Tutoring With Nonstandard Equations for Students With Mathematics Difficulty.

Science.gov (United States)

Powell, Sarah R; Driver, Melissa K; Julian, Tyler E

2015-01-01

Students often misinterpret the equal sign (=) as operational instead of relational. Research indicates misinterpretation of the equal sign occurs because students receive relatively little exposure to equations that promote relational understanding of the equal sign. No study, however, has examined effects of nonstandard equations on the equation solving and equal-sign understanding of students with mathematics difficulty (MD). In the present study, second-grade students with MD (n = 51) were randomly assigned to standard equations tutoring, combined tutoring (standard and nonstandard equations), and no-tutoring control. Combined tutoring students demonstrated greater gains on equation-solving assessments and equal-sign tasks compared to the other two conditions. Standard tutoring students demonstrated improved skill on equation solving over control students, but combined tutoring students' performance gains were significantly larger. Results indicate that exposure to and practice with nonstandard equations positively influence student understanding of the equal sign. © Hammill Institute on Disabilities 2013.

Numerical treatment of linearized equations describing inhomogeneous collisionless plasmas

International Nuclear Information System (INIS)

Lewis, H.R.

1979-01-01

The equations governing the small-signal response of spatially inhomogeneous collisionless plasmas have practical significance in physics, for example in controlled thermonuclear fusion research. Although the solutions are very complicated and the equations are different to solve numerically, effective methods for them are being developed which are applicable when the equilibrium involves only one nonignorable coordinate. The general theoretical framework probably will provide a basis for progress when there are two or three nonignorable coordinates

Numerical solution of a reaction-diffusion equation

International Nuclear Information System (INIS)

Moyano, Edgardo A.; Scarpettini, Alberto F.

2000-01-01

The purpose of the present work to continue the observations and the numerical experiences on a reaction-diffusion model, that is a simplified form of the neutronic flux equation. The model is parabolic, nonlinear, with Dirichlet boundary conditions. The purpose is to approximate non trivial solutions, asymptotically stables for t → ∞, that is solutions that tend to the elliptic problem, in the Lyapunov sense. It belongs to the so-called reaction-diffusion equations of semi linear kind, that is, linear equations in the heat operator and they have a nonlinear reaction function, in this case f (u, a, b) = u (a - b u), being u concentration, a and b parameters. The study of the incidence of these parameters take an interest to the neutronic flux physics. So that we search non trivial, positive and bounded solutions. The used algorithm is based on the concept of monotone and ordered sequences, and on the existence theorem of Amann and Sattinger. (author)

Iterated Crank-Nicolson method for hyperbolic and parabolic equations in numerical relativity

International Nuclear Information System (INIS)

Leiler, Gregor; Rezzolla, Luciano

2006-01-01

The iterated Crank-Nicolson is a predictor-corrector algorithm commonly used in numerical relativity for the solution of both hyperbolic and parabolic partial differential equations. We here extend the recent work on the stability of this scheme for hyperbolic equations by investigating the properties when the average between the predicted and corrected values is made with unequal weights and when the scheme is applied to a parabolic equation. We also propose a variant of the scheme in which the coefficients in the averages are swapped between two corrections leading to systematically larger amplification factors and to a smaller numerical dispersion

Introduction to precise numerical methods

CERN Document Server

Aberth, Oliver

2007-01-01

Precise numerical analysis may be defined as the study of computer methods for solving mathematical problems either exactly or to prescribed accuracy. This book explains how precise numerical analysis is constructed. The book also provides exercises which illustrate points from the text and references for the methods presented. All disc-based content for this title is now available on the Web. · Clearer, simpler descriptions and explanations ofthe various numerical methods· Two new types of numerical problems; accurately solving partial differential equations with the included software and computing line integrals in the complex plane.

Approximate solutions to Mathieu's equation

Science.gov (United States)

Wilkinson, Samuel A.; Vogt, Nicolas; Golubev, Dmitry S.; Cole, Jared H.

2018-06-01

Mathieu's equation has many applications throughout theoretical physics. It is especially important to the theory of Josephson junctions, where it is equivalent to Schrödinger's equation. Mathieu's equation can be easily solved numerically, however there exists no closed-form analytic solution. Here we collect various approximations which appear throughout the physics and mathematics literature and examine their accuracy and regimes of applicability. Particular attention is paid to quantities relevant to the physics of Josephson junctions, but the arguments and notation are kept general so as to be of use to the broader physics community.

Numerical Calculation of Transport Based on the Drift-Kinetic Equation for Plasmas in General Toroidal Magnetic Geometry: Numerical Methods

International Nuclear Information System (INIS)

Reynolds, J. M.; Lopez-Bruna, D.

2009-01-01

In this report we continue with the description of a newly developed numerical method to solve the drift kinetic equation for ions and electrons in toroidal plasmas. Several numerical aspects, already outlined in a previous report [Informes Tecnicos Ciemat 1165, mayo 2009], will be treated now in more detail. Aside from discussing the method in the context of other existing codes, various aspects will be now explained from the viewpoint of numerical methods: the way to solve convection equations, the adopted boundary conditions, the real-space meshing procedures along with a new software developed to build them, and some additional questions related with the parallelization and the numerical integration. (Author) 16 refs

A numerical study of adaptive space and time discretisations for Gross–Pitaevskii equations

Science.gov (United States)

Thalhammer, Mechthild; Abhau, Jochen

2012-01-01

As a basic principle, benefits of adaptive discretisations are an improved balance between required accuracy and efficiency as well as an enhancement of the reliability of numerical computations. In this work, the capacity of locally adaptive space and time discretisations for the numerical solution of low-dimensional nonlinear Schrödinger equations is investigated. The considered model equation is related to the time-dependent Gross–Pitaevskii equation arising in the description of Bose–Einstein condensates in dilute gases. The performance of the Fourier-pseudo spectral method constrained to uniform meshes versus the locally adaptive finite element method and of higher-order exponential operator splitting methods with variable time stepsizes is studied. Numerical experiments confirm that a local time stepsize control based on a posteriori local error estimators or embedded splitting pairs, respectively, is effective in different situations with an enhancement either in efficiency or reliability. As expected, adaptive time-splitting schemes combined with fast Fourier transform techniques are favourable regarding accuracy and efficiency when applied to Gross–Pitaevskii equations with a defocusing nonlinearity and a mildly varying regular solution. However, the numerical solution of nonlinear Schrödinger equations in the semi-classical regime becomes a demanding task. Due to the highly oscillatory and nonlinear nature of the problem, the spatial mesh size and the time increments need to be of the size of the decisive parameter 0Fourier pseudo-spectral and the finite element method. Nevertheless, for smaller parameter values locally adaptive time discretisations facilitate to determine the time stepsizes sufficiently small in order that the numerical approximation captures correctly the behaviour of the analytical solution. Further illustrations for Gross–Pitaevskii equations with a focusing nonlinearity or a sharp Gaussian as initial condition, respectively

Numerical solution of integral equations, describing mass spectrum of vector mesons

International Nuclear Information System (INIS)

Zhidkov, E.P.; Nikonov, E.G.; Sidorov, A.V.; Skachkov, N.B.; Khoromskij, B.N.

1988-01-01

The description of the numerical algorithm for solving quasipotential integral equation in impulse space is presented. The results of numerical computations of the vector meson mass spectrum and the leptonic decay width are given in comparison with the experimental data

Numerical Solution of Differential Algebraic Equations and Applications

DEFF Research Database (Denmark)

Thomsen, Per Grove

2005-01-01

These lecture notes have been written as part of a special course on the numerical solution of Differential Algebraic Equations and applications . The course was held at IMM in the spring of 2005. The authors of the different chapters have all taken part in the course and the chapters are written...

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.

Field data-based mathematical modeling by Bode equations and vector fitting algorithm for renewable energy applications

Science.gov (United States)

W. Hasan, W. Z.

2018-01-01

The power system always has several variations in its profile due to random load changes or environmental effects such as device switching effects when generating further transients. Thus, an accurate mathematical model is important because most system parameters vary with time. Curve modeling of power generation is a significant tool for evaluating system performance, monitoring and forecasting. Several numerical techniques compete to fit the curves of empirical data such as wind, solar, and demand power rates. This paper proposes a new modified methodology presented as a parametric technique to determine the system’s modeling equations based on the Bode plot equations and the vector fitting (VF) algorithm by fitting the experimental data points. The modification is derived from the familiar VF algorithm as a robust numerical method. This development increases the application range of the VF algorithm for modeling not only in the frequency domain but also for all power curves. Four case studies are addressed and compared with several common methods. From the minimal RMSE, the results show clear improvements in data fitting over other methods. The most powerful features of this method is the ability to model irregular or randomly shaped data and to be applied to any algorithms that estimating models using frequency-domain data to provide state-space or transfer function for the model. PMID:29351554

Field data-based mathematical modeling by Bode equations and vector fitting algorithm for renewable energy applications.

Science.gov (United States)

Sabry, A H; W Hasan, W Z; Ab Kadir, M Z A; Radzi, M A M; Shafie, S

2018-01-01

The power system always has several variations in its profile due to random load changes or environmental effects such as device switching effects when generating further transients. Thus, an accurate mathematical model is important because most system parameters vary with time. Curve modeling of power generation is a significant tool for evaluating system performance, monitoring and forecasting. Several numerical techniques compete to fit the curves of empirical data such as wind, solar, and demand power rates. This paper proposes a new modified methodology presented as a parametric technique to determine the system's modeling equations based on the Bode plot equations and the vector fitting (VF) algorithm by fitting the experimental data points. The modification is derived from the familiar VF algorithm as a robust numerical method. This development increases the application range of the VF algorithm for modeling not only in the frequency domain but also for all power curves. Four case studies are addressed and compared with several common methods. From the minimal RMSE, the results show clear improvements in data fitting over other methods. The most powerful features of this method is the ability to model irregular or randomly shaped data and to be applied to any algorithms that estimating models using frequency-domain data to provide state-space or transfer function for the model.

Field data-based mathematical modeling by Bode equations and vector fitting algorithm for renewable energy applications.

Directory of Open Access Journals (Sweden)

A H Sabry

Full Text Available The power system always has several variations in its profile due to random load changes or environmental effects such as device switching effects when generating further transients. Thus, an accurate mathematical model is important because most system parameters vary with time. Curve modeling of power generation is a significant tool for evaluating system performance, monitoring and forecasting. Several numerical techniques compete to fit the curves of empirical data such as wind, solar, and demand power rates. This paper proposes a new modified methodology presented as a parametric technique to determine the system's modeling equations based on the Bode plot equations and the vector fitting (VF algorithm by fitting the experimental data points. The modification is derived from the familiar VF algorithm as a robust numerical method. This development increases the application range of the VF algorithm for modeling not only in the frequency domain but also for all power curves. Four case studies are addressed and compared with several common methods. From the minimal RMSE, the results show clear improvements in data fitting over other methods. The most powerful features of this method is the ability to model irregular or randomly shaped data and to be applied to any algorithms that estimating models using frequency-domain data to provide state-space or transfer function for the model.

Cross Coursing in Mathematics: Physical Modelling in Differential Equations Crossing to Discrete Dynamical Systems

Science.gov (United States)

Winkel, Brian

2012-01-01

We give an example of cross coursing in which a subject or approach in one course in undergraduate mathematics is used in a completely different course. This situation crosses falling body modelling in an upper level differential equations course into a modest discrete dynamical systems unit of a first-year mathematics course. (Contains 1 figure.)

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.

Research in progress in applied mathematics, numerical analysis, fluid mechanics, and computer science

Science.gov (United States)

1994-01-01

This report summarizes research conducted at the Institute for Computer Applications in Science and Engineering in applied mathematics, fluid mechanics, and computer science during the period October 1, 1993 through March 31, 1994. The major categories of the current ICASE research program are: (1) applied and numerical mathematics, including numerical analysis and algorithm development; (2) theoretical and computational research in fluid mechanics in selected areas of interest to LaRC, including acoustics and combustion; (3) experimental research in transition and turbulence and aerodynamics involving LaRC facilities and scientists; and (4) computer science.

The Liquid Film Flow with Evaporation: Numerical Modelling

Directory of Open Access Journals (Sweden)

Rezanova Ekaterina

2016-01-01

Full Text Available The flow of thin liquid layer on an inclined substrate is investigated numerically. The mathematical modelling is based on the Oberbeck-Boussinesq equations and the generalized conditions on the thermocapillary boundary simplified during the parametrical analysis. In the framework of the long-wave approximation the evolution equation which determines the thickness of the liquid layer in the case of moderate Reynolds numbers is derived. The results of numerical modelling of the liquid flow with evaporation at the interface are obtained.

A numerical study of adaptive space and time discretisations for Gross-Pitaevskii equations.

Science.gov (United States)

Thalhammer, Mechthild; Abhau, Jochen

2012-08-15

As a basic principle, benefits of adaptive discretisations are an improved balance between required accuracy and efficiency as well as an enhancement of the reliability of numerical computations. In this work, the capacity of locally adaptive space and time discretisations for the numerical solution of low-dimensional nonlinear Schrödinger equations is investigated. The considered model equation is related to the time-dependent Gross-Pitaevskii equation arising in the description of Bose-Einstein condensates in dilute gases. The performance of the Fourier-pseudo spectral method constrained to uniform meshes versus the locally adaptive finite element method and of higher-order exponential operator splitting methods with variable time stepsizes is studied. Numerical experiments confirm that a local time stepsize control based on a posteriori local error estimators or embedded splitting pairs, respectively, is effective in different situations with an enhancement either in efficiency or reliability. As expected, adaptive time-splitting schemes combined with fast Fourier transform techniques are favourable regarding accuracy and efficiency when applied to Gross-Pitaevskii equations with a defocusing nonlinearity and a mildly varying regular solution. However, the numerical solution of nonlinear Schrödinger equations in the semi-classical regime becomes a demanding task. Due to the highly oscillatory and nonlinear nature of the problem, the spatial mesh size and the time increments need to be of the size of the decisive parameter [Formula: see text], especially when it is desired to capture correctly the quantitative behaviour of the wave function itself. The required high resolution in space constricts the feasibility of numerical computations for both, the Fourier pseudo-spectral and the finite element method. Nevertheless, for smaller parameter values locally adaptive time discretisations facilitate to determine the time stepsizes sufficiently small in order that

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.

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...

General existence principles for Stieltjes differential equations with applications to mathematical biology

Science.gov (United States)

López Pouso, Rodrigo; Márquez Albés, Ignacio

2018-04-01

Stieltjes differential equations, which contain equations with impulses and equations on time scales as particular cases, simply consist on replacing usual derivatives by derivatives with respect to a nondecreasing function. In this paper we prove new existence results for functional and discontinuous Stieltjes differential equations and we show that such general results have real world applications. Specifically, we show that Stieltjes differential equations are specially suitable to study populations which exhibit dormant states and/or very short (impulsive) periods of reproduction. In particular, we construct two mathematical models for the evolution of a silkworm population. Our first model can be explicitly solved, as it consists on a linear Stieltjes equation. Our second model, more realistic, is nonlinear, discontinuous and functional, and we deduce the existence of solutions by means of a result proven in this paper.

Recent developments of mathematical fluid mechanics

CERN Document Server

Giga, Yoshikazu; Kozono, Hideo; Okamoto, Hisashi; Yamazaki, Masao

2016-01-01

The book addresses recent developments of the mathematical research on the Navier-Stokes and Euler equations as well as on related problems. In particular, there are covered:   1) existence, uniqueness, and the regularity of weak solutions; 2) stability of the motion in rest and the asymptotic behavior of solutions; 3) singularity and blow-up of weak and strong solutions; 4) vorticity and energy conservation; 5) motions of rotating fluids, or of fluids surrounding a rotating body; 6) free boundary problems; 7) maximal regularity theory and other abstract results for mathematical fluid mechanics.   For this quarter century, these topics have been playing a central role in both pure and applied mathematics and having a great influence to the developm ent of the functional analysis, harmonic analysis and numerical analysis whose tools make a a substantial contribution to the investigation of nonlinear partial differential equations, particularly the Navier-Stokes and the Euler equations.      There are 24...

Mathematical Model of Age Aggression

OpenAIRE

Golovinski, P. A.

2013-01-01

We formulate a mathematical model of competition for resources between representatives of different age groups. A nonlinear kinetic integral-differential equation of the age aggression describes the process of redistribution of resources. It is shown that the equation of the age aggression has a stationary solution, in the absence of age-dependency in the interaction of different age groups. A numerical simulation of the evolution of resources for different initial distributions has done. It ...

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.

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

CERN Document Server

Constanda, Christian; Hamill, William

2016-01-01

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

Mathematical methods in the solution of the the Hamilton-Darwin and the Takagi-Taupin equations

International Nuclear Information System (INIS)

Werner, S.A.; Berliner, R.R.; Arif, M.; Missouri Univ., Columbia

1986-01-01

The diffraction of neutrons by a single crystal is intrinsically a multiple scattering problem. For an ideally imperfect mosaic crystal the Hamilton-Darwin transfer equations describe the coupling of the incident and diffracted beams; whereas, for a perfect crystal one must use the dynamical theory of diffraction, which can be recast in the form of two coupled partial differential equations commonly referred to as the Takagi-Taupin equations. From a mathematical point of view these two problems are equivalent, although the physical manifestations of the solutions are quite different. For the occasion of Professor Shull's seventieth birthday celebration, we bring together in this paper some of the mathematical techniques which we have found useful in elucidating the subtleties of the Bragg diffraction of neutron by crystals. (orig.)

Numerical study of the Kadomtsev-Petviashvili equation and dispersive shock waves

Science.gov (United States)

Grava, T.; Klein, C.; Pitton, G.

2018-02-01

A detailed numerical study of the long time behaviour of dispersive shock waves in solutions to the Kadomtsev-Petviashvili (KP) I equation is presented. It is shown that modulated lump solutions emerge from the dispersive shock waves. For the description of dispersive shock waves, Whitham modulation equations for KP are obtained. It is shown that the modulation equations near the soliton line are hyperbolic for the KPII equation while they are elliptic for the KPI equation leading to a focusing effect and the formation of lumps. Such a behaviour is similar to the appearance of breathers for the focusing nonlinear Schrödinger equation in the semiclassical limit.

Dynamic transition between fixed- and mobile-bed: mathematical and numerical aspects

Science.gov (United States)

Zugliani, Daniel; Pasqualini, Matteo; Rosatti, Giorgio

2017-04-01

Free-surface flows with high sediment transport (as debris flow or hyper-concentrated flow) are composed by a mixture of fluid and solid phase, usually water and sediment. When these flows propagate over loose beds, particles constituting the mixture of water and sediments strongly interact with the ones forming the bed, leading to erosion or deposition. However, there are lots of other situations when the mixture flows over rigid bedrocks or over artificially paved transects, so there is no mass exchange between bed and mixture. The two situations are usually referred to as, respectively, mobile- and fixed-bed conditions. From a mathematical point of view, the systems of Partial Differential Equations (PDEs) that describe these flows derive from mass and momentum balance of both phases, but, the two resulting PDEs systems are different. The main difference concerns the concentration: in the mobile-bed condition, the concentration is linked to the local flow conditions by means of a suitable rheological relation, while in the fixed-bed case, the concentration is an unknown of the problem. It is quite common that a free surface flow with high sediment transport, in its path, encounters both conditions. In the recent work of Rosatti & Zugliani 2015, the mathematical and numerical description of the transition between fixed- and mobile-bed was successfully resolved, for the case of low sediment transport phenomena, by the introduction of a suitable erodibility variable and satisfactory results were obtained. The main disadvantage of the approach is related to the erodibility variable, that changes in space, based on bed characteristics, but remains constant in time. However, the nature of the bed can change dynamically as result of deposition over fixed bed or high erosion over mobile bed. With this work, we extend the applicability of the mentioned approach to the more complex PDEs describing the hyper-concentrated flow. Moreover, we introduce a strategy that allows

Numerical solution of conservation equations in the transient model for the system thermal - hydraulics in the Korsar computer code

International Nuclear Information System (INIS)

Yudov, Y.V.

2001-01-01

The functional part of the KORSAR computer code is based on the computational unit for the reactor system thermal-hydraulics and other thermal power systems with water cooling. The two-phase flow dynamics of the thermal-hydraulic network is modelled by KORSAR in one-dimensional two-fluid (non-equilibrium and nonhomogeneous) approximation with the same pressure of both phases. Each phase is characterized by parameters averaged over the channel sections, and described by the conservation equations for mass, energy and momentum. The KORSAR computer code relies upon a novel approach to mathematical modelling of two-phase dispersed-annular flows. This approach allows a two-fluid model to differentiate the effects of the liquid film and droplets in the gas core on the flow characteristics. A semi-implicit numerical scheme has been chosen for deriving discrete analogs the conservation equations in KORSAR. In the semi-implicit numerical scheme, solution of finite-difference equations is reduced to the problem of determining the pressure field at a new time level. For the one-channel case, the pressure field is found from the solution of a system of linear algebraic equations by using the tri-diagonal matrix method. In the branched network calculation, the matrix of coefficients in the equations describing the pressure field is no longer tri-diagonal but has a sparseness structure. In this case, the system of linear equations for the pressure field can be solved with any of the known classical methods. Such an approach is implemented in the existing best-estimate thermal-hydraulic computer codes (TRAC, RELAP5, etc.) For the KORSAR computer code, we have developed a new non-iterative method for calculating the pressure field in the network of any topology. This method is based on the tri-diagonal matrix method and performs well when solving the thermal-hydraulic network problems. (author)

Numerical simulation of the regularized long wave equation by He's homotopy perturbation method

International Nuclear Information System (INIS)

Inc, Mustafa; Ugurlu, Yavuz

2007-01-01

In this Letter, we present the homotopy perturbation method (shortly HPM) for obtaining the numerical solution of the RLW equation. We obtain the exact and numerical solutions of the Regularized Long Wave (RLW) equation for certain initial condition. The initial approximation can be freely chosen with possible unknown constants which can be determined by imposing the boundary and initial conditions. Comparison of the results with those of other methods have led us to significant consequences. The numerical solutions are compared with the known analytical solutions

Partial differential equations II elements of the modern theory equations with constant coefficients

CERN Document Server

Shubin, M

1994-01-01

This book, the first printing of which was published as Volume 31 of the Encyclopaedia of Mathematical Sciences, contains a survey of the modern theory of general linear partial differential equations and a detailed review of equations with constant coefficients. Readers will be interested in an introduction to microlocal analysis and its applications including singular integral operators, pseudodifferential operators, Fourier integral operators and wavefronts, a survey of the most important results about the mixed problem for hyperbolic equations, a review of asymptotic methods including short wave asymptotics, the Maslov canonical operator and spectral asymptotics, a detailed description of the applications of distribution theory to partial differential equations with constant coefficients including numerous interesting special topics.

Numerical computations with GPUs

CERN Document Server

Kindratenko, Volodymyr

2014-01-01

This book brings together research on numerical methods adapted for Graphics Processing Units (GPUs). It explains recent efforts to adapt classic numerical methods, including solution of linear equations and FFT, for massively parallel GPU architectures. This volume consolidates recent research and adaptations, covering widely used methods that are at the core of many scientific and engineering computations. Each chapter is written by authors working on a specific group of methods; these leading experts provide mathematical background, parallel algorithms and implementation details leading to

Conservation form of the equations of fluid dynamics in general nonsteady coordinates

Science.gov (United States)

Zhang, H.; Camarero, R.; Kahawita, R.

1985-11-01

Many of the differential equations arising in fluid dynamics may be stated in conservation-law form. A number of investigations have been conducted with the aim to derive the conservation-law form of the Navier-Stokes equations in general nonsteady coordinate systems. The present note has the objective to illustrate a mathematical methodology with which such forms of the equations may be derived in an easier and more general fashion. For numerical applications, the scalar form of the equations is eventually provided. Attention is given to the conservation form of equations in curvilinear coordinates and numerical considerations.

Conservation form of the equations of fluid dynamics in general nonsteady coordinates

International Nuclear Information System (INIS)

Zhang, H.; Camarero, R.; Kahawita, R.

1985-01-01

Many of the differential equations arising in fluid dynamics may be stated in conservation-law form. A number of investigations have been conducted with the aim to derive the conservation-law form of the Navier-Stokes equations in general nonsteady coordinate systems. The present note has the objective to illustrate a mathematical methodology with which such forms of the equations may be derived in an easier and more general fashion. For numerical applications, the scalar form of the equations is eventually provided. Attention is given to the conservation form of equations in curvilinear coordinates and numerical considerations. 6 references

Numerical integration of the Teukolsky equation in the time domain

International Nuclear Information System (INIS)

Pazos-Avalos, Enrique; Lousto, Carlos O.

2005-01-01

We present a fourth-order convergent (2+1)-dimensional, numerical formalism to solve the Teukolsky equation in the time domain. Our approach is first to rewrite the Teukolsky equation as a system of first-order differential equations. In this way we get a system that has the form of an advection equation. This is then used in combination with a series expansion of the solution in powers of time. To obtain a fourth-order scheme we kept terms up to fourth derivative in time and use the advectionlike system of differential equations to substitute the temporal derivatives by spatial derivatives. This scheme is applied to evolve gravitational perturbations in the Schwarzschild and Kerr backgrounds. Our numerical method proved to be stable and fourth-order convergent in r* and θ directions. The correct power-law tail, ∼1/t 2l+3 , for general initial data, and ∼1/t 2l+4 , for time-symmetric data, was found in our runs. We noted that it is crucial to resolve accurately the angular dependence of the mode at late times in order to obtain these values of the exponents in the power-law decay. In other cases, when the decay was too fast and round-off error was reached before a tail was developed, then the quasinormal modes frequencies provided a test to determine the validity of our code

Undergraduate Students' Perceptions of Collaborative Learning in a Differential Equations Mathematics Course

Science.gov (United States)

Hajra, Sayonita Ghosh; Das, Ujjaini

2015-01-01

This paper uses collaborative learning strategies to examine students' perceptions in a differential equations mathematics course. Students' perceptions were analyzed using three collaborative learning strategies including collaborative activity, group-quiz and online discussion. The study results show that students identified both strengths and…

Minimal features of a computer and its basic software to executs NEPTUNIX 2 numerical step

International Nuclear Information System (INIS)

Roux, Pierre.

1982-12-01

NEPTUNIX 2 is a package which carries out the simulation of complex processes described by numerous non linear algebro-differential equations. Main features are: non linear or time dependent parameters, implicit form, stiff systems, dynamic change of equations leading to discontinuities on some variables. Thus the mathematical model is built with an equation set F(x,x',t,l) = 0, where t is the independent variable, x' the derivative of x and l an ''algebrized'' logical variable. The NEPTUNIX 2 package is divided into two successive major steps: a non numerical step and a numerical step. The non numerical step must be executed on a series 370 IBM computer or a compatible computer. This step generates a FORTRAN language model picture fitted for the computer carrying out the numerical step. The numerical step consists in building and running a mathematical model simulator. This execution step of NEPTUNIX 2 has been designed in order to be transportable on many computers. The present manual describes minimal features of such host computer used for executing the NEPTUNIX 2 numerical step [fr

The Navier-Stokes Equations Theory and Numerical Methods

CERN Document Server

Masuda, Kyûya; Rautmann, Reimund; Solonnikov, Vsevolod

1990-01-01

These proceedings contain original (refereed) research articles by specialists from many countries, on a wide variety of aspects of Navier-Stokes equations. Additionally, 2 survey articles intended for a general readership are included: one surveys the present state of the subject via open problems, and the other deals with the interplay between theory and numerical analysis.

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

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.

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.

Group Classification of a General Bond-Option Pricing Equation of Mathematical Finance

OpenAIRE

Motsepa, Tanki; Khalique, Chaudry Masood; Molati, Motlatsi

2014-01-01

We carry out group classification of a general bond-option pricing equation. We show that the equation admits a three-dimensional equivalence Lie algebra. We also show that some of the values of the constants which result from group classification give us well-known models in mathematics of finance such as Black-Scholes, Vasicek, and Cox-Ingersoll-Ross. For all such values of these arbitrary constants we obtain Lie point symmetries. Symmetry reductions are then obtained and group invariant so...

Soliton equations and Hamiltonian systems

CERN Document Server

Dickey, L A

2002-01-01

The theory of soliton equations and integrable systems has developed rapidly during the last 30 years with numerous applications in mechanics and physics. For a long time, books in this field have not been written but the flood of papers was overwhelming: many hundreds, maybe thousands of them. All this output followed one single work by Gardner, Green, Kruskal, and Mizura on the Korteweg-de Vries equation (KdV), which had seemed to be merely an unassuming equation of mathematical physics describing waves in shallow water. Besides its obvious practical use, this theory is attractive also becau

Numerical Schemes for Rough Parabolic Equations

Energy Technology Data Exchange (ETDEWEB)

Deya, Aurelien, E-mail: deya@iecn.u-nancy.fr [Universite de Nancy 1, Institut Elie Cartan Nancy (France)

2012-04-15

This paper is devoted to the study of numerical approximation schemes for a class of parabolic equations on (0,1) perturbed by a non-linear rough signal. It is the continuation of Deya (Electron. J. Probab. 16:1489-1518, 2011) and Deya et al. (Probab. Theory Relat. Fields, to appear), where the existence and uniqueness of a solution has been established. The approach combines rough paths methods with standard considerations on discretizing stochastic PDEs. The results apply to a geometric 2-rough path, which covers the case of the multidimensional fractional Brownian motion with Hurst index H>1/3.

Rotationally symmetric numerical solutions to the sine-Gordon equation

DEFF Research Database (Denmark)

Olsen, O. H.; Samuelsen, Mogens Rugholm

1981-01-01

We examine numerically the properties of solutions to the spherically symmetric sine-Gordon equation given an initial profile which coincides with the one-dimensional breather solution and refer to such solutions as ring waves. Expanding ring waves either exhibit a return effect or expand towards...

Numerical solution for multi-term fractional (arbitrary) orders differential equations

OpenAIRE

El-Sayed, A. M. A.; El-Mesiry, A. E. M.; El-Saka, H. A. A.

2004-01-01

Our main concern here is to give a numerical scheme to solve a nonlinear multi-term fractional (arbitrary) orders differential equation. Some results concerning the existence and uniqueness have been also obtained.

Applications of the Soave-Redlich-Kwong Equations of State Using Mathematic

Science.gov (United States)

Sun, Lanyi; Zhai, Cheng; Zhang, Hui

The application of the Peng-Robinson equations of state (PR EOS) using Matlab and Mathematic has already been demonstrated. In this paper, using Mathematic to solve Soave-Redlich-Kwong (SRK) EOS, as well as the estimation of pure component properties, plotting of vapor-liquid equilibrium (VLE) diagram and calculation of chemical equilibrium, is presented. First the SRK EOS is used to predict several pure-component properties, such as liquid and gas molar volumes for isobutane. The vapor-liquid isobaric diagram is then plotted for a binary mixture composed of n-pentane and n-hexane under the pressures of 2*10^5 and 8*10^5 Pa respectively.

Numerical computation of soliton dynamics for NLS equations in a driving potential

Directory of Open Access Journals (Sweden)

Marco Caliari

2010-06-01

Full Text Available We provide numerical computations for the soliton dynamics of the nonlinear Schrodinger equation with an external potential. After computing the ground state solution r of a related elliptic equation we show that, in the semi-classical regime, the center of mass of the solution with initial datum built upon r is driven by the solution to $ddot x=- abla V(x$. Finally, we provide examples and analyze the numerical errors in the two dimensional case when V is a harmonic potential.

On the numerical treatment of the Griffin-Hill-Eheeler equation

International Nuclear Information System (INIS)

Galleti, D.; Toledo Piza, A.F.R.

The precision attainable in the numerical treatment of the Griffin-Hill-Wheeler equation is studied in a solvable model. Trucation errors related to the generator coordinate kinematics are exhibited and briefly discussed [pt

Numerical soliton-like solutions of the potential Kadomtsev-Petviashvili equation by the decomposition method

International Nuclear Information System (INIS)

Kaya, Dogan; El-Sayed, Salah M.

2003-01-01

In this Letter we present an Adomian's decomposition method (shortly ADM) for obtaining the numerical soliton-like solutions of the potential Kadomtsev-Petviashvili (shortly PKP) equation. We will prove the convergence of the ADM. We obtain the exact and numerical solitary-wave solutions of the PKP equation for certain initial conditions. Then ADM yields the analytic approximate solution with fast convergence rate and high accuracy through previous works. The numerical solutions are compared with the known analytical solutions

Mathematical geophysics an introduction to rotating fluids and the Navier-Stokes equations

CERN Document Server

Chemin, Jean-Yves; Gallagher, Isabelle; Grenier, Emmanuel

2006-01-01

Aimed at graduate students and researchers in mathematics, engineering, oceanography, meteorology and mechanics, this text provides a detailed introduction to the physical theory of rotating fluids, a significant part of geophysical fluid dynamics. The Navier-Stokes equations are examined in both incompressible and rapidly rotating forms.

Mathematical and numerical modeling of early atherosclerotic lesions***

Directory of Open Access Journals (Sweden)

Raoult Annie

2010-12-01

Full Text Available This article is devoted to the construction of a mathematical model describing the early formation of atherosclerotic lesions. The early stage of atherosclerosis is an inflammatory process that starts with the penetration of low density lipoproteins in the intima and with their oxidation. This phenomenon is closely linked to the local blood flow dynamics. Extending a previous work [5] that was mainly restricted to a one-dimensional setting, we couple a simple lesion growth model relying on the biomolecular process that takes place in the intima with blood flow dynamics and mass transfer. We perform numerical simulations on a two-dimensional geometry taken from [6,7] that mimicks a carotid artery deformed by a perivascular cast and we compare the numerical results with experimental data.

An uncoupling strategy for numerically solving the dynamic thermoelasticity equations

International Nuclear Information System (INIS)

Moura, C.A. de; Feijoo, R.A.

1981-01-01

The dynamic equations of coupled linear thermoelasticity are presented. A numerical algorithm which combines finite-element space approximation with a two-step time discretization in such a way as to reach significant computational savings is presented: It features a strategy for independently calculating the displacement and temperature fields through equations that nevertheless remain coupled. The scheme convergence was shown to be optimal and its machine performance, as ilustrated by some examples, fairly satisfactory. (Author) [pt

Numerical treatments for solving nonlinear mixed integral equation

Directory of Open Access Journals (Sweden)

M.A. Abdou

2016-12-01

Full Text Available We consider a mixed type of nonlinear integral equation (MNLIE of the second kind in the space C[0,T]×L2(Ω,T<1. The Volterra integral terms (VITs are considered in time with continuous kernels, while the Fredholm integral term (FIT is considered in position with singular general kernel. Using the quadratic method and separation of variables method, we obtain a nonlinear system of Fredholm integral equations (NLSFIEs with singular kernel. A Toeplitz matrix method, in each case, is then used to obtain a nonlinear algebraic system. Numerical results are calculated when the kernels take a logarithmic form or Carleman function. Moreover, the error estimates, in each case, are then computed.

Mathematical modeling of efficient protocols to control glioma growth.

Science.gov (United States)

Branco, J R; Ferreira, J A; de Oliveira, Paula

2014-09-01

In this paper we propose a mathematical model to describe the evolution of glioma cells taking into account the viscoelastic properties of brain tissue. The mathematical model is established considering that the glioma cells are of two phenotypes: migratory and proliferative. The evolution of the migratory cells is described by a diffusion-reaction equation of non Fickian type deduced considering a mass conservation law with a non Fickian migratory mass flux. The evolution of the proliferative cells is described by a reaction equation. A stability analysis that leads to the design of efficient protocols is presented. Numerical simulations that illustrate the behavior of the mathematical model are included. Copyright © 2014 Elsevier Inc. All rights reserved.

Wind laws for shockless initialization. [numerical forecasting model

Science.gov (United States)

Ghil, M.; Shkoller, B.

1976-01-01

A system of diagnostic equations for the velocity field, or wind laws, was derived for each of a number of models of large-scale atmospheric flow. The derivation in each case is mathematically exact and does not involve any physical assumptions not already present in the prognostic equations, such as nondivergence or vanishing of derivatives of the divergence. Therefore, initial states computed by solving these diagnostic equations should be compatible with the type of motion described by the prognostic equations of the model and should not generate initialization shocks when inserted into the model. Numerical solutions of the diagnostic system corresponding to a barotropic model are exhibited. Some problems concerning the possibility of implementing such a system in operational numerical weather prediction are discussed.

Fractional Calculus in Hydrologic Modeling: A Numerical Perspective

Energy Technology Data Exchange (ETDEWEB)

David A. Benson; Mark M. Meerschaert; Jordan Revielle

2012-01-01

Fractional derivatives can be viewed either as a handy extension of classical calculus or, more deeply, as mathematical operators defined by natural phenomena. This follows the view that the diffusion equation is defined as the governing equation of a Brownian motion. In this paper, we emphasize that fractional derivatives come from the governing equations of stable Levy motion, and that fractional integration is the corresponding inverse operator. Fractional integration, and its multi-dimensional extensions derived in this way, are intimately tied to fractional Brownian (and Levy) motions and noises. By following these general principles, we discuss the Eulerian and Lagrangian numerical solutions to fractional partial differential equations, and Eulerian methods for stochastic integrals. These numerical approximations illuminate the essential nature of the fractional calculus.

Mathematical oncology 2013

CERN Document Server

Gandolfi, Alberto

2014-01-01

With chapters on free boundaries, constitutive equations, stochastic dynamics, nonlinear diffusion–consumption, structured populations, and applications of optimal control theory, this volume presents the most significant recent results in the field of mathematical oncology. It highlights the work of world-class research teams, and explores how different researchers approach the same problem in various ways. Tumors are complex entities that present numerous challenges to the mathematical modeler. First and foremost, they grow. Thus their spatial mean field description involves a free boundary problem. Second, their interiors should be modeled as nontrivial porous media using constitutive equations. Third, at the end of anti-cancer therapy, a small number of malignant cells remain, making the post-treatment dynamics inherently stochastic. Fourth, the growth parameters of macroscopic tumors are non-constant, as are the parameters of anti-tumor therapies. Changes in these parameters may induce phenomena that a...

Zdeněk Kopal: Numerical Analyst

Science.gov (United States)

Křížek, M.

2015-07-01

We give a brief overview of Zdeněk Kopal's life, his activities in the Czech Astronomical Society, his collaboration with Vladimír Vand, and his studies at Charles University, Cambridge, Harvard, and MIT. Then we survey Kopal's professional life. He published 26 monographs and 20 conference proceedings. We will concentrate on Kopal's extensive monograph Numerical Analysis (1955, 1961) that is widely accepted to be the first comprehensive textbook on numerical methods. It describes, for instance, methods for polynomial interpolation, numerical differentiation and integration, numerical solution of ordinary differential equations with initial or boundary conditions, and numerical solution of integral and integro-differential equations. Special emphasis will be laid on error analysis. Kopal himself applied numerical methods to celestial mechanics, in particular to the N-body problem. He also used Fourier analysis to investigate light curves of close binaries to discover their properties. This is, in fact, a problem from mathematical analysis.

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

Science.gov (United States)

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

2013-03-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 techniques can also be extended to other kinds of the multi-term fractional time-space models with fractional Laplacian.

A Predictor-Corrector Approach for the Numerical Solution of Fractional Differential Equations

Science.gov (United States)

Diethelm, Kai; Ford, Neville J.; Freed, Alan D.; Gray, Hugh R. (Technical Monitor)

2002-01-01

We discuss an Adams-type predictor-corrector method for the numerical solution of fractional differential equations. The method may be used both for linear and for nonlinear problems, and it may be extended to multi-term equations (involving more than one differential operator) too.

Mathematical analysis, approximation theory and their applications

CERN Document Server

Gupta, Vijay

2016-01-01

Designed for graduate students, researchers, and engineers in mathematics, optimization, and economics, this self-contained volume presents theory, methods, and applications in mathematical analysis and approximation theory. Specific topics include: approximation of functions by linear positive operators with applications to computer aided geometric design, numerical analysis, optimization theory, and solutions of differential equations. Recent and significant developments in approximation theory, special functions and q-calculus along with their applications to mathematics, engineering, and social sciences are discussed and analyzed. Each chapter enriches the understanding of current research problems and theories in pure and applied research.

Nonlinear reaction-diffusion equations with delay: some theorems, test problems, exact and numerical solutions

Science.gov (United States)

Polyanin, A. D.; Sorokin, V. G.

2017-12-01

The paper deals with nonlinear reaction-diffusion equations with one or several delays. We formulate theorems that allow constructing exact solutions for some classes of these equations, which depend on several arbitrary functions. Examples of application of these theorems for obtaining new exact solutions in elementary functions are provided. We state basic principles of construction, selection, and use of test problems for nonlinear partial differential equations with delay. Some test problems which can be suitable for estimating accuracy of approximate analytical and numerical methods of solving reaction-diffusion equations with delay are presented. Some examples of numerical solutions of nonlinear test problems with delay are considered.

Matrix-oriented implementation for the numerical solution of the partial differential equations governing flows and transport in porous media

KAUST Repository

Sun, Shuyu

2012-09-01

In this paper we introduce a new technique for the numerical solution of the various partial differential equations governing flow and transport phenomena in porous media. This method is proposed to be used in high level programming languages like MATLAB, Python, etc., which show to be more efficient for certain mathematical operations than for others. The proposed technique utilizes those operations in which these programming languages are efficient the most and keeps away as much as possible from those inefficient, time-consuming operations. In particular, this technique is based on the minimization of using multiple indices looping operations by reshaping the unknown variables into one-dimensional column vectors and performing the numerical operations using shifting matrices. The cell-centered information as well as the face-centered information are shifted to the adjacent face-center and cell-center, respectively. This enables the difference equations to be done for all the cells at once using matrix operations rather than within loops. Furthermore, for results post-processing, the face-center information can further be mapped to the physical grid nodes for contour plotting and stream lines constructions. In this work we apply this technique to flow and transport phenomena in porous media. © 2012 Elsevier Ltd.

Mathematical and numerical study of nonlinear boundary problems related to plasma physics

International Nuclear Information System (INIS)

Sermange, M.

1982-06-01

After the study of some equations based on the Hodgkin-Huxley model, the work presented here is concerned with nonlinear boundary problems in MHD. They are gathered in two subjects: equilibrium equations and stability equations. The axisymmetric MHD equilibrium equations with free boundary have been studied by different authors, particularly the existence, regularity, unicity and non-unicity. Here, bifurcation, convergence of calculation methods existence of solutions in a discontinuous frame are studied. MHD stability can be determined by the principle of Bernstein et al; the mathematical work concerned here bears on the equivalence, in the case of two-dimensional or axisymmetric stability, between this model and a scalar eigenvalue problem which is introduced. At last, modules for computing MHD equilibrium for the simulation of plasma confinement in a tokamak are described [fr

Group Classification of a General Bond-Option Pricing Equation of Mathematical Finance

Directory of Open Access Journals (Sweden)

Tanki Motsepa

2014-01-01

Full Text Available We carry out group classification of a general bond-option pricing equation. We show that the equation admits a three-dimensional equivalence Lie algebra. We also show that some of the values of the constants which result from group classification give us well-known models in mathematics of finance such as Black-Scholes, Vasicek, and Cox-Ingersoll-Ross. For all such values of these arbitrary constants we obtain Lie point symmetries. Symmetry reductions are then obtained and group invariant solutions are constructed for some cases.

Hamiltonian partial differential equations and applications

CERN Document Server

Nicholls, David; Sulem, Catherine

2015-01-01

This book is a unique selection of work by world-class experts exploring the latest developments in Hamiltonian partial differential equations and their applications. Topics covered within are representative of the field’s wide scope, including KAM and normal form theories, perturbation and variational methods, integrable systems, stability of nonlinear solutions as well as applications to cosmology, fluid mechanics and water waves. The volume contains both surveys and original research papers and gives a concise overview of the above topics, with results ranging from mathematical modeling to rigorous analysis and numerical simulation. It will be of particular interest to graduate students as well as researchers in mathematics and physics, who wish to learn more about the powerful and elegant analytical techniques for Hamiltonian partial differential equations.

A Study of Enhanced, Higher Order Boussinesq-Type Equations and Their Numerical Modelling

DEFF Research Database (Denmark)

Banijamali, Babak

model is designated for the solution of higher-order Boussinesq-type equations, formulated in terms of the horizontal velocity at an arbitrary depth vector. Various discretisation techniques and grid definitions have been considered in this endeavour, undertaking a detailed analysis of the selected......This project has encompassed efforts in two separate veins: on the one hand, the acquiring of highly accurate model equations of the Boussinesq-type, and on the other hand, the theoretical and practical work in implementing such equations in the form of conventional numerical models, with obvious...... potential for applications to the realm of numerical modelling in coastal engineering. The derivation and analysis of several forms of higher-order in dispersion and non-linearity Boussinesq-type equations have been undertaken, obtaining and investigating the properties of a new and generalised class...

Numerical solution of neutral functional-differential equations with proportional delays

Directory of Open Access Journals (Sweden)

Mehmet Giyas Sakar

2017-07-01

Full Text Available In this paper, homotopy analysis method is improved with optimal determination of auxiliary parameter by use of residual error function for solving neutral functional-differential equations (NFDEs with proportional delays. Convergence analysis and error estimate of method are given. Some numerical examples are solved and comparisons are made with the existing results. The numerical results show that the homotopy analysis method with residual error function is very effective and simple.

Incompressible Navier-Stokes equations. Theory and practice

Energy Technology Data Exchange (ETDEWEB)

Gjesdal, T.

1996-12-31

This paper contains notes from a seminar presented at the Dept. of Mathematics in the University of Bergen, Norway, Oct. 1996. It first introduces the theory of existence and uniqueness of solutions to the incompressible Navier-Stokes equation and defines a well-posed initial-boundary value problem. It then discusses different methods for solving numerically the Navier-Stokes equations in velocity-pressure formulation. The emphasis is on pressure correction methods. 19 refs.

Scattering integral equations and four nucleon problem

International Nuclear Information System (INIS)

Narodetskii, I.M.

1980-01-01

Existing results from the application of integral equation technique to the four-nucleon bound states and scattering are reviewed. The first numerical calculations of the four-body integral equations have been done ten years ago. Yet, it is still widely believed that these equations are too complicated to solve numerically. The purpose of this review is to provide a clear and elementary introduction in the integral equation method and to demonstrate its usefulness in physical applications. The presentation is based on the quasiparticle approach. This permits a simple interpretation of the equations in terms of quasiparticle scattering. The mathematical basis for the quasiparticle approach is the Hilbert-Schmidt method of the Fredholm integral equation theory. The first part of this review contains a detailed discussion of the Hilbert-Schmidt expansion as applied to the 2-particle amplitudes and to the kernel of the four-body equations. The second part contains the discussion of the four-body quasiparticle equations and of the resed forullts obtain bound states and scattering

Higher-order automatic differentiation of mathematical functions

Science.gov (United States)

Charpentier, Isabelle; Dal Cappello, Claude

2015-04-01

Functions of mathematical physics such as the Bessel functions, the Chebyshev polynomials, the Gauss hypergeometric function and so forth, have practical applications in many scientific domains. On the one hand, differentiation formulas provided in reference books apply to real or complex variables. These do not account for the chain rule. On the other hand, based on the chain rule, the automatic differentiation has become a natural tool in numerical modeling. Nevertheless automatic differentiation tools do not deal with the numerous mathematical functions. This paper describes formulas and provides codes for the higher-order automatic differentiation of mathematical functions. The first method is based on Faà di Bruno's formula that generalizes the chain rule. The second one makes use of the second order differential equation they satisfy. Both methods are exemplified with the aforementioned functions.

A numerical study of the integral equations for the laser fields in free-electron lasers

International Nuclear Information System (INIS)

Yoo, J. G.; Park, S. H.; Jeong, Y. U.; Lee, B. C.; Rhee, Y. J.; Cho, S. O.

2004-01-01

The dynamics of the radiation fields in free-electron lasers is investigated on the basis of the integro-differential equations in the one-dimensional formulation. For simple cases we solved the integro-differential equations analytically and numerically to test our numerical procedures developed on the basis of the Filon method. The numerical results showed good agreement with the analytical solutions. To confirm the legitimacy of the numerical package, we carried out numerical studies on the inhomogeneous broadening effects, where no analytic solutions are available, due to the energy spread and the emittance of the electron beam.

Numerical solution of the Navier--Stokes equations at high Reynolds numbers

International Nuclear Information System (INIS)

Shestakov, A.I.

1974-01-01

A numerical method is presented which is designed to solve the Navier-Stokes equations for two-dimensional, incompressible flow. The method is intended for use on problems with high Reynolds numbers for which calculations via finite difference methods have been unattainable or unreliable. The proposed scheme is a hybrid utilizing a time-splitting finite difference method in areas away from the boundaries. In areas neighboring the boundaries, the equations of motion are solved by the newly proposed vortex method by Chorin. The major accomplishment of the new scheme is that it contains a simple way for merging the two methods at the interface of the two subdomains. The proposed algorithm is designed for use on the time-dependent equations but can be used on steady state problems as well. The method is tested on the popular, time-independent, square cavity problem, an example of a separated flow with closed streamlines. Numerical results are presented for a Reynolds number of 10 3 . (auth)

On the numerical solution of the neutron fractional diffusion equation

International Nuclear Information System (INIS)

Maleki Moghaddam, Nader; Afarideh, Hossein; Espinosa-Paredes, Gilberto

2014-01-01

Highlights: • The new version of neutron diffusion equation which established on the fractional derivatives is presented. • The Neutron Fractional Diffusion Equation (NFDE) is solved in the finite differences frame. • NFDE is solved using shifted Grünwald-Letnikov definition of fractional operators. • The results show that “K eff ” strongly depends on the order of fractional derivative. - Abstract: In order to core calculation in the nuclear reactors there is a new version of neutron diffusion equation which is established on the fractional partial derivatives, named Neutron Fractional Diffusion Equation (NFDE). In the NFDE model, neutron flux in each zone depends directly on the all previous zones (not only on the nearest neighbors). Under this circumstance, it can be said that the NFDE has the space history. We have developed a one-dimension code, NFDE-1D, which can simulate the reactor core using arbitrary exponent of differential operators. In this work a numerical solution of the NFDE is presented using shifted Grünwald-Letnikov definition of fractional derivative in finite differences frame. The model is validated with some numerical experiments where different orders of fractional derivative are considered (e.g. 0.999, 0.98, 0.96, and 0.94). The results show that the effective multiplication factor (K eff ) depends strongly on the order of fractional derivative

Numerical Validation of Chemical Compositional Model for Wettability Alteration Processes

Science.gov (United States)

Bekbauov, Bakhbergen; Berdyshev, Abdumauvlen; Baishemirov, Zharasbek; Bau, Domenico

2017-12-01

Chemical compositional simulation of enhanced oil recovery and surfactant enhanced aquifer remediation processes is a complex task that involves solving dozens of equations for all grid blocks representing a reservoir. In the present work, we perform a numerical validation of the newly developed mathematical formulation which satisfies the conservation laws of mass and energy and allows applying a sequential solution approach to solve the governing equations separately and implicitly. Through its application to the numerical experiment using a wettability alteration model and comparisons with existing chemical compositional model's numerical results, the new model has proven to be practical, reliable and stable.

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.

A Finite Difference, Semi-implicit, Equation-of-State Efficient Algorithm for the Compositional Flow Modeling in the Subsurface: Numerical Examples

KAUST Repository

Saavedra, Sebastian

2012-07-01

The mathematical model that has been recognized to have the more accurate approximation to the physical laws govern subsurface hydrocarbon flow in reservoirs is the Compositional Model. The features of this model are adequate to describe not only the performance of a multiphase system but also to represent the transport of chemical species in a porous medium. Its importance relies not only on its current relevance to simulate petroleum extraction processes, such as, Primary, Secondary, and Enhanced Oil Recovery Process (EOR) processes but also, in the recent years, carbon dioxide (CO2) sequestration. The purpose of this study is to investigate the subsurface compositional flow under isothermal conditions for several oil well cases. While simultaneously addressing computational implementation finesses to contribute to the efficiency of the algorithm. This study provides the theoretical framework and computational implementation subtleties of an IMplicit Pressure Explicit Composition (IMPEC)-Volume-balance (VB), two-phase, equation-of-state, approach to model isothermal compositional flow based on the finite difference scheme. The developed model neglects capillary effects and diffusion. From the phase equilibrium premise, the model accounts for volumetric performances of the phases, compressibility of the phases, and composition-dependent viscosities. The Equation of State (EoS) employed to approximate the hydrocarbons behaviour is the Peng Robinson Equation of State (PR-EOS). Various numerical examples were simulated. The numerical results captured the complex physics involved, i.e., compositional, gravitational, phase-splitting, viscosity and relative permeability effects. Regarding the numerical scheme, a phase-volumetric-flux estimation eases the calculation of phase velocities by naturally fitting to phase-upstream-upwinding. And contributes to a faster computation and an efficient programming development.

Ernst Equation and Riemann Surfaces: Analytical and Numerical Methods

International Nuclear Information System (INIS)

Ernst, Frederick J

2007-01-01

metric tensor components. The first two chapters of this book are devoted to some basic ideas: in the introductory chapter 1 the authors discuss the concept of integrability, comparing the integrability of the vacuum Ernst equation with the integrability of nonlinear equations of Korteweg-de Vries (KdV) type, while in chapter 2 they describe various circumstances in which the vacuum Ernst equation has been determined to be relevant, not only in connection with gravitation but also, for example, in the construction of solutions of the self-dual Yang-Mills equations. It is also in this chapter that one of several equivalent linear systems for the Ernst equation is described. The next two chapters are devoted to Dmitry Korotkin's concept of algebro-geometric solutions of a linear system: in chapter 3 the structure of such solutions of the vacuum Ernst equation, which involve Riemann theta functions of hyperelliptic algebraic curves of any genus, is contrasted with the periodic structure of such solutions of the KdV equation. How such solutions can be obtained, for example, by solving a matrix Riemann-Hilbert problem and how the metric tensor of the associated spacetime can be evaluated is described in detail. In chapter 4 the asymptotic behaviour and the similarity structure of the general algebro-geometric solutions of the Ernst equation are described, and the relationship of such solutions to the perhaps more familiar multi-soliton solutions is discussed. The next three chapters are based upon the authors' own published research: in chapter 5 it is shown that a problem involving counter-rotating infinitely thin disks of matter can be solved in terms of genus two Riemann theta functions, while in chapter 6 the authors describe numerical methods that facilitate the construction of such solutions, and in chapter 7 three-dimensional graphs are displayed that depict all metrical fields of the associated spacetime. Finally, in chapter 8, the difficulties associated with

Applied mathematical methods in nuclear thermal hydraulics

International Nuclear Information System (INIS)

Ransom, V.H.; Trapp, J.A.

1983-01-01

Applied mathematical methods are used extensively in modeling of nuclear reactor thermal-hydraulic behavior. This application has required significant extension to the state-of-the-art. The problems encountered in modeling of two-phase fluid transients and the development of associated numerical solution methods are reviewed and quantified using results from a numerical study of an analogous linear system of differential equations. In particular, some possible approaches for formulating a well-posed numerical problem for an ill-posed differential model are investigated and discussed. The need for closer attention to numerical fidelity is indicated

From stochastic processes to numerical methods: A new scheme for solving reaction subdiffusion fractional partial differential equations

Energy Technology Data Exchange (ETDEWEB)

Angstmann, C.N.; Donnelly, I.C. [School of Mathematics and Statistics, UNSW Australia, Sydney NSW 2052 (Australia); Henry, B.I., E-mail: B.Henry@unsw.edu.au [School of Mathematics and Statistics, UNSW Australia, Sydney NSW 2052 (Australia); Jacobs, B.A. [School of Computer Science and Applied Mathematics, University of the Witwatersrand, Johannesburg, Private Bag 3, Wits 2050 (South Africa); DST–NRF Centre of Excellence in Mathematical and Statistical Sciences (CoE-MaSS) (South Africa); Langlands, T.A.M. [Department of Mathematics and Computing, University of Southern Queensland, Toowoomba QLD 4350 (Australia); Nichols, J.A. [School of Mathematics and Statistics, UNSW Australia, Sydney NSW 2052 (Australia)

2016-02-15

We have introduced a new explicit numerical method, based on a discrete stochastic process, for solving a class of fractional partial differential equations that model reaction subdiffusion. The scheme is derived from the master equations for the evolution of the probability density of a sum of discrete time random walks. We show that the diffusion limit of the master equations recovers the fractional partial differential equation of interest. This limiting procedure guarantees the consistency of the numerical scheme. The positivity of the solution and stability results are simply obtained, provided that the underlying process is well posed. We also show that the method can be applied to standard reaction–diffusion equations. This work highlights the broader applicability of using discrete stochastic processes to provide numerical schemes for partial differential equations, including fractional partial differential equations.

Exact and numerical solutions of generalized Drinfeld-Sokolov equations

International Nuclear Information System (INIS)

Ugurlu, Yavuz; Kaya, Dogan

2008-01-01

In this Letter, we consider a system of generalized Drinfeld-Sokolov (gDS) equations which models one-dimensional nonlinear wave processes in two-component media. We find some exact solutions of gDS by using tanh function method and we also obtain a numerical solution by using the Adomian's Decomposition Method (ADM)

Fast numerical solution of KKR-CPA equations: Testing new algorithms

Energy Technology Data Exchange (ETDEWEB)

Bruno, E.; Florio, G.M.; Ginatempo, B.; Giuliano, E.S. (Universita di Messina (Italy))

1994-04-01

Some numerical methods for the solution of KKR-CPA equations are discussed and tested. New, efficient, computational algorithms are proposed, allowing a remarkable reduction of computing time and a good reliability in evaluating spectral quantities. 16 refs., 7 figs.

Mathematical and numerical study of non-linear models used in plasma physics

International Nuclear Information System (INIS)

Ebrard, G.

2005-12-01

We study the interaction of several crossing beams with a plasma in the Laser-Megajoule context. We start from Euler-Maxwell. The formal asymptotic is the Zakharov system. For simplified systems of Klein-Gordon-wave type, we justify an approximation by a Zakharov equation for solutions of large amplitude. We construct a new system that simulates the interaction of 2 beams and present a whole hierarchy of models. We introduce a numerical scheme using the known results on Zakharov-wave equations which are valid for short pulses. We give a scheme which eliminate the backscattering wave. We give some numerical results. Finally, we do several numerical simulations of laser-plasma interaction for the initial value problem and the boundary value problem. (author)

Applied analysis and differential equations

CERN Document Server

Cârj, Ovidiu

2007-01-01

This volume contains refereed research articles written by experts in the field of applied analysis, differential equations and related topics. Well-known leading mathematicians worldwide and prominent young scientists cover a diverse range of topics, including the most exciting recent developments. A broad range of topics of recent interest are treated: existence, uniqueness, viability, asymptotic stability, viscosity solutions, controllability and numerical analysis for ODE, PDE and stochastic equations. The scope of the book is wide, ranging from pure mathematics to various applied fields such as classical mechanics, biomedicine, and population dynamics.

Mathematical, physical and numerical principles essential for models of turbulent mixing

Energy Technology Data Exchange (ETDEWEB)

Sharp, David Howland [Los Alamos National Laboratory; Lim, Hyunkyung [STONY BROOK UNIV; Yu, Yan [STONY BROOK UNIV; Glimm, James G [STONY BROOK UNIV

2009-01-01

We propose mathematical, physical and numerical principles which are important for the modeling of turbulent mixing, especially the classical and well studied Rayleigh-Taylor and Richtmyer-Meshkov instabilities which involve acceleration driven mixing of a fluid discontinuity layer, by a steady accerleration or an impulsive force.

Numerical Solutions of Mechanical Turbulent Filtration Equation Used in Mechatronics and Micro Mechanic

OpenAIRE

Hassan Fathabadi

2013-01-01

In this study, several novel numerical solutions are presented to solve the turbulent filtration equation and its special case called “Non-Newtonian mechanical filtration equation”. The turbulent filtration equation in porous media is a very important equation which has many applications to solve the problems appearing especially in mechatronics, micro mechanic and fluid mechanic. Many applied mechanical problems can be solved using this equation. For example, non-Newtonian mechanical filtrat...

RDTM solution of Caputo time fractional-order hyperbolic telegraph equation

Directory of Open Access Journals (Sweden)

Vineet K. Srivastava

2013-03-01

Full Text Available In this study, a mathematical model has been developed for the second order hyperbolic one-dimensional time fractional Telegraph equation (TFTE. The fractional derivative has been described in the Caputo sense. The governing equations have been solved by a recent reliable semi-analytic method known as the reduced differential transformation method (RDTM. The method is a powerful mathematical technique for solving wide range of problems. Using RDTM method, it is possible to find exact solution as well as closed approximate solution of any ordinary or partial differential equation. Three numerical examples of TFTE have been provided in order to check the effectiveness, accuracy and convergence of the method. The computed results are also depicted graphically.

Generalized Roe's numerical scheme for a two-fluid model

International Nuclear Information System (INIS)

Toumi, I.; Raymond, P.

1993-01-01

This paper is devoted to a mathematical and numerical study of a six equation two-fluid model. We will prove that the model is strictly hyperbolic due to the inclusion of the virtual mass force term in the phasic momentum equations. The two-fluid model is naturally written under a nonconservative form. To solve the nonlinear Riemann problem for this nonconservative hyperbolic system, a generalized Roe's approximate Riemann solver, is used, based on a linearization of the nonconservative terms. A Godunov type numerical scheme is built, using this approximate Riemann solver. 10 refs., 5 figs,

A numerical spectral approach to solve the dislocation density transport equation

International Nuclear Information System (INIS)

Djaka, K S; Taupin, V; Berbenni, S; Fressengeas, C

2015-01-01

A numerical spectral approach is developed to solve in a fast, stable and accurate fashion, the quasi-linear hyperbolic transport equation governing the spatio-temporal evolution of the dislocation density tensor in the mechanics of dislocation fields. The approach relies on using the Fast Fourier Transform algorithm. Low-pass spectral filters are employed to control both the high frequency Gibbs oscillations inherent to the Fourier method and the fast-growing numerical instabilities resulting from the hyperbolic nature of the transport equation. The numerical scheme is validated by comparison with an exact solution in the 1D case corresponding to dislocation dipole annihilation. The expansion and annihilation of dislocation loops in 2D and 3D settings are also produced and compared with finite element approximations. The spectral solutions are shown to be stable, more accurate for low Courant numbers and much less computation time-consuming than the finite element technique based on an explicit Galerkin-least squares scheme. (paper)

Mathematics for plasma physics; Mathematiques pour la physique des plasmas

Energy Technology Data Exchange (ETDEWEB)

Sentis, R. [CEA Bruyeres-le-Chatel, 91 (France)

2011-01-15

The plasma physics is in the heart of the research of the CEA-DAM. Using mathematics in this domain is necessary, particularly for a precise statement of the partial differential equations systems which are on the basis of the numerical simulations. Examples are given concerning hydrodynamics, models for the thermal conduction and laser-plasma interaction. For the bi-temperature compressible Euler model, the mathematical study of the problem has allowed us to understand why the role of the energy equations dealing with ions on one hand and electrons on the other hand are not identical despite the symmetrical appearance of the system. The mathematical study is also necessary to be sure of the existence and uniqueness of the solution

Exact and numerical solutions of generalized Drinfeld-Sokolov 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: dkaya36@yahoo.com

2008-04-14

In this Letter, we consider a system of generalized Drinfeld-Sokolov (gDS) equations which models one-dimensional nonlinear wave processes in two-component media. We find some exact solutions of gDS by using tanh function method and we also obtain a numerical solution by using the Adomian's Decomposition Method (ADM)

Differences between quadratic equations and functions: Indonesian pre-service secondary mathematics teachers’ views

Science.gov (United States)

Aziz, T. A.; Pramudiani, P.; Purnomo, Y. W.

2018-01-01

Difference between quadratic equation and quadratic function as perceived by Indonesian pre-service secondary mathematics teachers (N = 55) who enrolled at one private university in Jakarta City was investigated. Analysis of participants’ written responses and interviews were conducted consecutively. Participants’ written responses highlighted differences between quadratic equation and function by referring to their general terms, main characteristics, processes, and geometrical aspects. However, they showed several obstacles in describing the differences such as inappropriate constraints and improper interpretations. Implications of the study are discussed.

Numerical artifacts in the Generalized Porous Medium Equation: Why harmonic averaging itself is not to blame

Science.gov (United States)

Maddix, Danielle C.; Sampaio, Luiz; Gerritsen, Margot

2018-05-01

The degenerate parabolic Generalized Porous Medium Equation (GPME) poses numerical challenges due to self-sharpening and its sharp corner solutions. For these problems, we show results for two subclasses of the GPME with differentiable k (p) with respect to p, namely the Porous Medium Equation (PME) and the superslow diffusion equation. Spurious temporal oscillations, and nonphysical locking and lagging have been reported in the literature. These issues have been attributed to harmonic averaging of the coefficient k (p) for small p, and arithmetic averaging has been suggested as an alternative. We show that harmonic averaging is not solely responsible and that an improved discretization can mitigate these issues. Here, we investigate the causes of these numerical artifacts using modified equation analysis. The modified equation framework can be used for any type of discretization. We show results for the second order finite volume method. The observed problems with harmonic averaging can be traced to two leading error terms in its modified equation. This is also illustrated numerically through a Modified Harmonic Method (MHM) that can locally modify the critical terms to remove the aforementioned numerical artifacts.

The unique and shared contributions of arithmetic operation understanding and numerical magnitude representation to children's mathematics achievement.

Science.gov (United States)

Wong, Terry Tin-Yau

2017-12-01

The current study examined the unique and shared contributions of arithmetic operation understanding and numerical magnitude representation to children's mathematics achievement. A sample of 124 fourth graders was tested on their arithmetic operation understanding (as reflected by their understanding of arithmetic principles and the knowledge about the application of arithmetic operations) and their precision of rational number magnitude representation. They were also tested on their mathematics achievement and arithmetic computation performance as well as the potential confounding factors. The findings suggested that both arithmetic operation understanding and numerical magnitude representation uniquely predicted children's mathematics achievement. The findings highlight the significance of arithmetic operation understanding in mathematics learning. Copyright © 2017 Elsevier Inc. All rights reserved.

An integral equation-based numerical solver for Taylor states in toroidal geometries

Science.gov (United States)

O'Neil, Michael; Cerfon, Antoine J.

2018-04-01

We present an algorithm for the numerical calculation of Taylor states in toroidal and toroidal-shell geometries using an analytical framework developed for the solution to the time-harmonic Maxwell equations. Taylor states are a special case of what are known as Beltrami fields, or linear force-free fields. The scheme of this work relies on the generalized Debye source representation of Maxwell fields and an integral representation of Beltrami fields which immediately yields a well-conditioned second-kind integral equation. This integral equation has a unique solution whenever the Beltrami parameter λ is not a member of a discrete, countable set of resonances which physically correspond to spontaneous symmetry breaking. Several numerical examples relevant to magnetohydrodynamic equilibria calculations are provided. Lastly, our approach easily generalizes to arbitrary geometries, both bounded and unbounded, and of varying genus.

Numerical simulation of two phase flows in heat exchangers

International Nuclear Information System (INIS)

Grandotto Biettoli, M.

2006-04-01

The author gives an overview of his research activity since 1981. He first gives a detailed presentation of properties and equations of two-phase flows in heat exchangers, and of their mathematical and numerical investigation: semi-local equations (mass conservation, momentum conservation and energy conservation), homogenized conservation equations (mass, momentum and enthalpy conservation, boundary conditions), equation closures, discretization, resolution algorithm, computational aspects and applications. Then, he reports the works performed in the field of turbulent flows, hyperbolic methods, low Mach methods, the Neptune project, and parallel computing

Methods for constructing exact solutions of partial differential equations mathematical and analytical techniques with applications to engineering

CERN Document Server

Meleshko, Sergey V

2005-01-01

Differential equations, especially nonlinear, present the most effective way for describing complex physical processes. Methods for constructing exact solutions of differential equations play an important role in applied mathematics and mechanics. This book aims to provide scientists, engineers and students with an easy-to-follow, but comprehensive, description of the methods for constructing exact solutions of differential equations.

A global numerical solution of the radial Schroedinger equation by second-order perturbation theory

International Nuclear Information System (INIS)

Adam, G.

1979-01-01

A global numerical method, which uses second-order perturbation theory, is described for the solution of the radial Schroedinger equation. The perturbative numerical (PN) solution is derived in two stages: first, the original potential is approximated by a piecewise continuous parabolic function, and second, the resulting Schroedinger equation is solved on each integration step by second-order perturbation theory, starting with a step function reference approximation for the parabolic potential. We get a manageable PN algorithm, which shows an order of accuracy equal to six in the solution of the original Schroedinger equation, and is very stable against round off errors. (author)

Mathematical and computational methods for semiclassical Schrödinger equations

KAUST Repository

Jin, Shi; Markowich, Peter; Sparber, Christof

2011-01-01

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

Numerical simulation of nonlinear continuity equations by evolving diffeomorphisms

KAUST Repository

Carrillo, José A.

2016-09-22

In this paper we present a numerical scheme for nonlinear continuity equations, which is based on the gradient flow formulation of an energy functional with respect to the quadratic transportation distance. It can be applied to a large class of nonlinear continuity equations, whose dynamics are driven by internal energies, given external potentials and/or interaction energies. The solver is based on its variational formulation as a gradient flow with respect to the Wasserstein distance. Positivity of solutions as well as energy decrease of the semi-discrete scheme are guaranteed by its construction. We illustrate this property with various examples in spatial dimension one and two.

Exponential Convergence for Numerical Solution of Integral Equations Using Radial Basis Functions

Directory of Open Access Journals (Sweden)

Zakieh Avazzadeh

2014-01-01

Full Text Available We solve some different type of Urysohn integral equations by using the radial basis functions. These types include the linear and nonlinear Fredholm, Volterra, and mixed Volterra-Fredholm integral equations. Our main aim is to investigate the rate of convergence to solve these equations using the radial basis functions which have normic structure that utilize approximation in higher dimensions. Of course, the use of this method often leads to ill-posed systems. Thus we propose an algorithm to improve the results. Numerical results show that this method leads to the exponential convergence for solving integral equations as it was already confirmed for partial and ordinary differential equations.

Elliptic Diophantine equations a concrete approach via the elliptic logarithm

CERN Document Server

Tzanakis, Nikos

2013-01-01

This book presents in a unified way the beautiful and deep mathematics, both theoretical and computational, on which the explicit solution of an elliptic Diophantine equation is based. It collects numerous results and methods that are scattered in literature. Some results are even hidden behind a number of routines in software packages, like Magma. This book is suitable for students in mathematics, as well as professional mathematicians.

Unsteady analytical solutions to the Poisson–Nernst–Planck equations

International Nuclear Information System (INIS)

Schönke, Johannes

2012-01-01

It is shown that the Poisson–Nernst–Planck equations for a single ion species can be formulated as one equation in terms of the electric field. This previously not analyzed equation shows similarities to the vector Burgers equation and is identical with it in the one dimensional case. Several unsteady exact solutions for one and multidimensional cases are presented. Besides new mathematical insights which these first known unsteady solutions give, they can serve as test cases in computer simulations to analyze numerical algorithms and to verify code. (paper)

A numerical solution of the coupled proton-H atom transport equations for the proton aurora

International Nuclear Information System (INIS)

Basu, B.; Jasperse, J.R.; Grossbard, N.J.

1990-01-01

A numerical code has been developed to solve the coupled proton-H atom linear transport equations for the proton aurora. The transport equations have been simplified by using plane-parallel geometry and the forward-scattering approximations only. Otherwise, the equations and their numerical solutions are exact. Results are presented for the particle fluxes and the energy deposition rates, and they are compared with the previous analytical results that were obtained by using additional simplifying approximations. It is found that although the analytical solutions for the particle fluxes differ somewhat from the numerical solutions, the energy deposition rates calculated by the two methods agree to within a few percent. The accurate particle fluxes given by the numerical code are useful for accurate calculation of the characteristic quantities of the proton aurora, such as the ionization rates and the emission rates

International Conference on Advances in Applied Mathematics

CERN Document Server

Hammami, Mohamed; Masmoudi, Afif

2015-01-01

This contributed volume presents some recent theoretical advances in mathematics and its applications in various areas of science and technology.   Written by internationally recognized scientists and researchers, the chapters in this book are based on talks given at the International Conference on Advances in Applied Mathematics (ICAAM), which took place December 16-19, 2013, in Hammamet, Tunisia.  Topics discussed at the conference included spectral theory, operator theory, optimization, numerical analysis, ordinary and partial differential equations, dynamical systems, control theory, probability, and statistics.  These proceedings aim to foster and develop further growth in all areas of applied mathematics.

Intentional and automatic processing of numerical information in mathematical anxiety: testing the influence of emotional priming.

Science.gov (United States)

Ashkenazi, Sarit

2018-02-05

Current theoretical approaches suggest that mathematical anxiety (MA) manifests itself as a weakness in quantity manipulations. This study is the first to examine automatic versus intentional processing of numerical information using the numerical Stroop paradigm in participants with high MA. To manipulate anxiety levels, we combined the numerical Stroop task with an affective priming paradigm. We took a group of college students with high MA and compared their performance to a group of participants with low MA. Under low anxiety conditions (neutral priming), participants with high MA showed relatively intact number processing abilities. However, under high anxiety conditions (mathematical priming), participants with high MA showed (1) higher processing of the non-numerical irrelevant information, which aligns with the theoretical view regarding deficits in selective attention in anxiety and (2) an abnormal numerical distance effect. These results demonstrate that abnormal, basic numerical processing in MA is context related.

Numerical simulation of trapped dipolar quantum gases: Collapse studies and vortex dynamics

KAUST Repository

Sparber, Christof; Markowich, Peter; Huang, Zhongyi

2010-01-01

We numerically study the three dimensional Gross-Pitaevskii equation for dipolar quantum gases using a time-splitting algorithm. We are mainly concerned with numerical investigations of the possible blow-up of solutions, i.e. collapse of the condensate, and the dynamics of vortices. © American Institute of Mathematical Sciences.

Identifying generalized Fitzhugh-Nagumo equation from a numerical solution of Hodgkin-Huxley model

Directory of Open Access Journals (Sweden)

Nikola V. Georgiev

2003-01-01

Full Text Available An analytic time series in the form of numerical solution (in an appropriate finite time interval of the Hodgkin-Huxley current clamped (HHCC system of four differential equations, well known in the neurophysiology as an exact empirical model of excitation of a giant axon of Loligo, is presented. Then we search for a second-order differential equation of generalized Fitzhugh-Nagumo (GFN type, having as a solution the given single component (action potential of the numerical solution. The given time series is used as a basis for reconstructing orders, powers, and coefficients of the polynomial right-hand sides of GFN equation approximately governing the process of action potential. For this purpose, a new geometrical method for determining phase space dimension of the unknown dynamical system (GFN equation and a specific modification of least squares method for identifying unknown coefficients are developed and applied.

Parameter Estimation for Partial Differential Equations by Collage-Based Numerical Approximation

Directory of Open Access Journals (Sweden)

Xiaoyan Deng

2009-01-01

into a minimization problem of a function of several variables after the partial differential equation is approximated by a differential dynamical system. Then numerical schemes for solving this minimization problem are proposed, including grid approximation and ant colony optimization. The proposed schemes are applied to a parameter estimation problem for the Belousov-Zhabotinskii equation, and the results show that the proposed approximation method is efficient for both linear and nonlinear partial differential equations with respect to unknown parameters. At worst, the presented method provides an excellent starting point for traditional inversion methods that must first select a good starting point.

Analysis and Application of High Resolution Numerical Perturbation Algorithm for Convective-Diffusion Equation

International Nuclear Information System (INIS)

Gao Zhi; Shen Yi-Qing

2012-01-01

The high resolution numerical perturbation (NP) algorithm is analyzed and tested using various convective-diffusion equations. The NP algorithm is constructed by splitting the second order central difference schemes of both convective and diffusion terms of the convective-diffusion equation into upstream and downstream parts, then the perturbation reconstruction functions of the convective coefficient are determined using the power-series of grid interval and eliminating the truncated errors of the modified differential equation. The important nature, i.e. the upwind dominance nature, which is the basis to ensuring that the NP schemes are stable and essentially oscillation free, is firstly presented and verified. Various numerical cases show that the NP schemes are efficient, robust, and more accurate than the original second order central scheme

Mathematical and numerical analysis of plasma stability

International Nuclear Information System (INIS)

Saramito, B.

1987-11-01

Equilibrium of a tokamak plasma is analyzed using two two-dimensional numerical models. Plasma configuration; convection in a cylindrical plasma layer; and tearing instabilities in a flat layer are considered. The finite element code used is explained. The existence of analogous stationary solutions for a problem concerning compressible fluids is shown. Stationary convection created by the equilibrium density gradient is treated. Approximation using fluid equations is employed in the case of convection resulting from the equilibrium temperature gradient. Evolution towards turbulence of incompressible fluid models is followed [fr

Elements of the history of mathematics

CERN Document Server

Bourbaki, Nicolas

1994-01-01

This work gathers together, without substantial modification, the major­ ity of the historical Notes which have appeared to date in my Elements de M atMmatique. Only the flow has been made independent of the Elements to which these Notes were attached; they are therefore, in principle, accessible to every reader who possesses a sound classical mathematical background, of undergraduate standard. Of course, the separate studies which make up this volume could not in any way pretend to sketch, even in a summary manner, a complete and con­ nected history of the development of Mathematics up to our day. Entire parts of classical mathematics such as differential Geometry, algebraic Geometry, the Calculus of variations, are only mentioned in passing; others, such as the theory of analytic functions, that of differential equations or partial differ­ ential equations, are hardly touched on; all the more do these gaps become more numerous and more important as the modern era is reached. It goes without saying that t...

Associations of non-symbolic and symbolic numerical magnitude processing with mathematical competence: a meta-analysis.

Science.gov (United States)

Schneider, Michael; Beeres, Kassandra; Coban, Leyla; Merz, Simon; Susan Schmidt, S; Stricker, Johannes; De Smedt, Bert

2017-05-01

Many studies have investigated the association between numerical magnitude processing skills, as assessed by the numerical magnitude comparison task, and broader mathematical competence, e.g. counting, arithmetic, or algebra. Most correlations were positive but varied considerably in their strengths. It remains unclear whether and to what extent the strength of these associations differs systematically between non-symbolic and symbolic magnitude comparison tasks and whether age, magnitude comparison measures or mathematical competence measures are additional moderators. We investigated these questions by means of a meta-analysis. The literature search yielded 45 articles reporting 284 effect sizes found with 17,201 participants. Effect sizes were combined by means of a two-level random-effects regression model. The effect size was significantly higher for the symbolic (r = .302, 95% CI [.243, .361]) than for the non-symbolic (r = .241, 95% CI [.198, .284]) magnitude comparison task and decreased very slightly with age. The correlation was higher for solution rates and Weber fractions than for alternative measures of comparison proficiency. It was higher for mathematical competencies that rely more heavily on the processing of magnitudes (i.e. mental arithmetic and early mathematical abilities) than for others. The results support the view that magnitude processing is reliably associated with mathematical competence over the lifespan in a wide range of tasks, measures and mathematical subdomains. The association is stronger for symbolic than for non-symbolic numerical magnitude processing. So symbolic magnitude processing might be a more eligible candidate to be targeted by diagnostic screening instruments and interventions for school-aged children and for adults. © 2016 John Wiley & Sons Ltd.

A numerical guide to the solution of the bidomain equations of cardiac electrophysiology

KAUST Repository

Pathmanathan, Pras

2010-06-01

Simulation of cardiac electrical activity using the bidomain equations can be a massively computationally demanding problem. This study provides a comprehensive guide to numerical bidomain modelling. Each component of bidomain simulations-discretisation, ODE-solution, linear system solution, and parallelisation-is discussed, and previously-used methods are reviewed, new methods are proposed, and issues which cause particular difficulty are highlighted. Particular attention is paid to the choice of stimulus currents, compatibility conditions for the equations, the solution of singular linear systems, and convergence of the numerical scheme. © 2010 Elsevier Ltd.

A numerical guide to the solution of the bidomain equations of cardiac electrophysiology

KAUST Repository

Pathmanathan, Pras; Bernabeu, Miguel O.; Bordas, Rafel; Cooper, Jonathan; Garny, Alan; Pitt-Francis, Joe M.; Whiteley, Jonathan P.; Gavaghan, David J.

2010-01-01

Simulation of cardiac electrical activity using the bidomain equations can be a massively computationally demanding problem. This study provides a comprehensive guide to numerical bidomain modelling. Each component of bidomain simulations-discretisation, ODE-solution, linear system solution, and parallelisation-is discussed, and previously-used methods are reviewed, new methods are proposed, and issues which cause particular difficulty are highlighted. Particular attention is paid to the choice of stimulus currents, compatibility conditions for the equations, the solution of singular linear systems, and convergence of the numerical scheme. © 2010 Elsevier Ltd.

Exploring New Physics Frontiers Through Numerical Relativity.

Science.gov (United States)

Cardoso, Vitor; Gualtieri, Leonardo; Herdeiro, Carlos; Sperhake, Ulrich

2015-01-01

The demand to obtain answers to highly complex problems within strong-field gravity has been met with significant progress in the numerical solution of Einstein's equations - along with some spectacular results - in various setups. We review techniques for solving Einstein's equations in generic spacetimes, focusing on fully nonlinear evolutions but also on how to benchmark those results with perturbative approaches. The results address problems in high-energy physics, holography, mathematical physics, fundamental physics, astrophysics and cosmology.

Comparison of numerical approaches to solve a Poincare-covariant Faddeev equation

International Nuclear Information System (INIS)

Alkofer, R.; Eichmann, G.; Krassnigg, A.; Schwinzerl, M.

2006-01-01

Full text: The quark core of Baryons can be described with the help of the numerical solution of the Poincare-Faddeev equation. Hereby the used elements, as e.g. the quark propagator are taken from non-perturbative studies of Landau gauge QCD. Different numerical approaches to solve in this way the relativistic three quark problem are compared and benchmarked results for the efficiency of different algorithms are presented. (author)

Numerical approximation of null controls for the heat equation: Ill-posedness and remedies

International Nuclear Information System (INIS)

Münch, Arnaud; Zuazua, Enrique

2010-01-01

The numerical approximation of exact or trajectory controls for the wave equation is known to be a delicate issue, since the pioneering work of Glowinski–Lions in the nineties, because of the anomalous behavior of the high-frequency spurious numerical waves. Various efficient remedies have been developed and analyzed in the last decade to filter out these high-frequency components: Fourier filtering, Tychonoff's regularization, mixed finite-element methods, multi-grid strategies, etc. Recently convergence rate results have also been obtained. This work is devoted to analyzing this issue for the heat equation, which is the opposite paradigm because of its strong dissipativity and smoothing properties. The existing analytical results guarantee that, at least in some simple situations, as in the finite-difference scheme in 1 − d, the null or trajectory controls for numerical approximation schemes converge. This is due to the intrinsic high-frequency damping of the heat equation that is inherited by its numerical approximation schemes. But when developing numerical simulations the topic appears to be much more subtle and difficult. In fact, efficiently computing the null control for a numerical approximation scheme of the heat equation is a difficult problem in itself. The difficulty is strongly related to the regularizing effect of the heat kernel. The controls of minimal L 2 -norm are characterized as minima of quadratic functionals on the solutions of the adjoint heat equation, or its numerical versions. These functionals are shown to be coercive in very large spaces of solutions, sufficient to guarantee the L 2 character of controls, but very far from being identifiable as energy spaces for the adjoint system. The very weak coercivity of the functionals under consideration makes the approximation problem exponentially ill-posed and the functional framework far from being well adapted to standard techniques in numerical analysis. In practice, the controls of the

Numerical simulation of the regularized long wave equation by He's homotopy perturbation method

Energy Technology Data Exchange (ETDEWEB)

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

2007-09-17

In this Letter, we present the homotopy perturbation method (shortly HPM) for obtaining the numerical solution of the RLW equation. We obtain the exact and numerical solutions of the Regularized Long Wave (RLW) equation for certain initial condition. The initial approximation can be freely chosen with possible unknown constants which can be determined by imposing the boundary and initial conditions. Comparison of the results with those of other methods have led us to significant consequences. The numerical solutions are compared with the known analytical solutions.

Numerical Analysis Of Hooke Jeeves-Runge Kutta To Determine Reaction Rate Equation In Pyrrole Polymerization

International Nuclear Information System (INIS)

Gunawan, Indra; Sulistyo, Harry; Rochmad

2001-01-01

The numerical analysis of Hooke Jeeves Methods combined with Runge Kutta Methods is used to determine the exact model of reaction rate equation of pyrrole polymerization. Chemical polymerization of pyrrole was conducted with FeCI 3 / pyrrole solution at concentration ratio of 1.62 mole / mole and 2.18 mole / mole with varrying temperature of 28, 40, 50, and 60 o C. FeCl 3 acts as an oxidation agent to form pyrrole cation that will polymerize. The numerical analysis was done to examine the exact model of reaction rate equation which is derived from reaction equation of initiation, propagation, and termination. From its numerical analysis, it is found that the pyrrole polymerization follows third order of pyrrole cation concentration

APPROACHED DECISION OF THE DIFFERENTIAL EQUATIONS

Directory of Open Access Journals (Sweden)

Oleksii B. Krasnozhon

2011-02-01

Full Text Available The urgency of the material stated in the article is caused by necessity of development, updating and improvements of methodical operating time on subject matters of issue "Calculus mathematics" which teaching is carried out in conditions of use of information-communication technologies. In the article the program realizations in Mathcad environment of Adams and Runge-Kutt methods of the approached decision of the differential equations are offered; examples on application of the specified methods are brought; the expediency of application of Mathcad environment during mathematical preparation of experts is proved. Perspective directions of the further scientific researches are methodical, mathematical and algorithmic aspects of creation of effective program realizations of numerical methods in Mathcad environment.

Numerical solutions of the monoenergetic neutron transport equation with anisotropic scattering

International Nuclear Information System (INIS)

Dahl, B.

1985-01-01

The Boltzmann equation for monoenergetic neutrons has been solved numerically with high accuracy for homogeneous slabs and spheres with various degree of linear anisotropy. Vacuum boundary conditions are used. The numerical method is based on previous work by Carlvik. Benchmark values of the criticality factor and higher order eigenvalues are given for multiplying systems of thickness or diameter from 10 -5 to 20 mean free paths and with anisotropy coefficients from 0.0 to 0.3. For slab geometry, both even and odd mode eigenvalues are treated. With increasing anisotropy, an increasing number of complex eigenvalues is observer. The total flux is calculated from the eigenvector and tables of the fundamental mode flux are given. Accurate extrapolation distances are derived for various dimensions and anisotropy coefficients from our eigenvalue results on slabs and spheres and from the work by Sanchez on infinite cylinders.The time eigenvalue spectrum in subcritical systems has also been studied. First, the connection between the eigenvalues arising from the time dependent and stationary transport equation is established. Based on this, the spectrum of real time eigenvalues in slabs and spheres is calculated. For spheres, the existence of complex time eigenvalues in the region beyond the value corresponding to the Corngold limit is numerically established. The presence of such eigenvalues has earlier not been proved. It is further shown that the Boltzmann equation for a sphere is significantly simplified when the decay constant is at the Corngold limit. The spectrum of sphere diameters corresponding to this decay constant is calculated for various linear anisotropies, and detailed numerical results are given. (Author)

Relations of Different Types of Numerical Magnitude Representations to Each Other and to Mathematics Achievement

Science.gov (United States)

Fazio, Lisa K.; Bailey, Drew H.; Thompson, Clarissa A.; Siegler, Robert S.

2014-01-01

We examined relations between symbolic and non-symbolic numerical magnitude representations, between whole number and fraction representations, and between these representations and overall mathematics achievement in fifth graders. Fraction and whole number symbolic and non-symbolic numerical magnitude understandings were measured using both…

Nonlinear ordinary differential equations analytical approximation and numerical methods

CERN Document Server

Hermann, Martin

2016-01-01

The book discusses the solutions to nonlinear ordinary differential equations (ODEs) using analytical and numerical approximation methods. Recently, analytical approximation methods have been largely used in solving linear and nonlinear lower-order ODEs. It also discusses using these methods to solve some strong nonlinear ODEs. There are two chapters devoted to solving nonlinear ODEs using numerical methods, as in practice high-dimensional systems of nonlinear ODEs that cannot be solved by analytical approximate methods are common. Moreover, it studies analytical and numerical techniques for the treatment of parameter-depending ODEs. The book explains various methods for solving nonlinear-oscillator and structural-system problems, including the energy balance method, harmonic balance method, amplitude frequency formulation, variational iteration method, homotopy perturbation method, iteration perturbation method, homotopy analysis method, simple and multiple shooting method, and the nonlinear stabilized march...

Formulation and numerical analysis of nonisothermal multiphase flow in porous media

International Nuclear Information System (INIS)

Martinez, M.J.

1995-06-01

A mathematical formulation is presented for describing the transport of air, water and energy through porous media. The development follows a continuum mechanics approach. The theory assumes the existence of various average macroscopic variables which describe the state of the system. Balance equations for mass and energy are formulated in terms of these macroscopic variables. The system is supplemented with constitutive equations relating fluxes to the state variables, and with transport property specifications. Specification of various mixing rules and thermodynamic relations completes the system of equations. A numerical simulation scheme, employing the method of lines, is described for one-dimensional flow. The numerical method is demonstrated on sample problems involving nonisothermal flow of air and water. The implementation is verified by comparison with existing numerical solutions

Numerical time-dependent partial differential equations for scientists and engineers

CERN Document Server

Brio, Moysey; Zakharian, Aramais R

2010-01-01

It is the first text that in addition to standard convergence theory treats other necessary ingredients for successful numerical simulations of physical systems encountered by every practitioner. The book is aimed at users with interests ranging from application modeling to numerical analysis and scientific software development. It is strongly influenced by the authors research in in space physics, electrical and optical engineering, applied mathematics, numerical analysis and professional software development. The material is based on a year-long graduate course taught at the University of Arizona since 1989. The book covers the first two-semesters of a three semester series. The second semester is based on a semester-long project, while the third semester requirement consists of a particular methods course in specific disciplines like computational fluid dynamics, finite element method in mechanical engineering, computational physics, biology, chemistry, photonics, etc. The first three chapters focus on bas...

Visualising magnetic fields numerical equation solvers in action

CERN Document Server

Beeteson, John Stuart

2001-01-01

Visualizing Magnetic Fields: Numerical Equation Solvers in Action provides a complete description of the theory behind a new technique, a detailed discussion of the ways of solving the equations (including a software visualization of the solution algorithms), the application software itself, and the full source code. Most importantly, there is a succinct, easy-to-follow description of each procedure in the code.The physicist Michael Faraday said that the study of magnetic lines of force was greatly influential in leading him to formulate many of those concepts that are now so fundamental to our modern world, proving to him their "great utility as well as fertility." Michael Faraday could only visualize these lines in his mind's eye and, even with modern computers to help us, it has been very expensive and time consuming to plot lines of force in magnetic fields

An introduction to the mathematical theory of the Navier-Stokes equations

CERN Document Server

Galdi, Giovanni P

1994-01-01

Undoubtedly, the Navier-Stokes equations are of basic importance within the context of modern theory of partial differential equations. Although the range of their applicability to concrete problems has now been clearly recognised to be limited, as my dear friend and bright colleague K.R. Ra­ jagopal has showed me by several examples during the past six years, the mathematical questions that remain open are of such a fascinating and challenging nature that analysts and applied mathematicians cannot help being attracted by them and trying to contribute to their resolution. Thus, it is not a coincidence that over the past ten years more than seventy sig­ nificant research papers have appeared concerning the well-posedness of boundary and initial-boundary value problems. In this monograph I shall perform a systematic and up-to-date investiga­ tion of the fundamental properties of the Navier-Stokes equations, including existence, uniqueness, and regularity of solutions and, whenever the region of flow is unbou...

Scilab software as an alternative low-cost computing in solving the linear equations problem

Science.gov (United States)

Agus, Fahrul; Haviluddin

2017-02-01

Numerical computation packages are widely used both in teaching and research. These packages consist of license (proprietary) and open source software (non-proprietary). One of the reasons to use the package is a complexity of mathematics function (i.e., linear problems). Also, number of variables in a linear or non-linear function has been increased. The aim of this paper was to reflect on key aspects related to the method, didactics and creative praxis in the teaching of linear equations in higher education. If implemented, it could be contribute to a better learning in mathematics area (i.e., solving simultaneous linear equations) that essential for future engineers. The focus of this study was to introduce an additional numerical computation package of Scilab as an alternative low-cost computing programming. In this paper, Scilab software was proposed some activities that related to the mathematical models. In this experiment, four numerical methods such as Gaussian Elimination, Gauss-Jordan, Inverse Matrix, and Lower-Upper Decomposition (LU) have been implemented. The results of this study showed that a routine or procedure in numerical methods have been created and explored by using Scilab procedures. Then, the routine of numerical method that could be as a teaching material course has exploited.

Handbook of Nonlinear Partial Differential Equations

CERN Document Server

Polyanin, Andrei D

2011-01-01

New to the Second Edition More than 1,000 pages with over 1,500 new first-, second-, third-, fourth-, and higher-order nonlinear equations with solutions Parabolic, hyperbolic, elliptic, and other systems of equations with solutions Some exact methods and transformations Symbolic and numerical methods for solving nonlinear PDEs with Maple(t), Mathematica(R), and MATLAB(R) Many new illustrative examples and tables A large list of references consisting of over 1,300 sources To accommodate different mathematical backgrounds, the authors avoid wherever possible the use of special terminology. They

Hermite interpolant multiscaling functions for numerical solution of the convection diffusion equations

Directory of Open Access Journals (Sweden)

Elmira Ashpazzadeh

2018-04-01

Full Text Available A numerical technique based on the Hermite interpolant multiscaling functions is presented for the solution of Convection-diusion equations. The operational matrices of derivative, integration and product are presented for multiscaling functions and are utilized to reduce the solution of linear Convection-diusion equation to the solution of algebraic equations. Because of sparsity of these matrices, this method is computationally very attractive and reduces the CPU time and computer memory. Illustrative examples are included to demonstrate the validity and applicability of the new technique.

Solution of wave-like equation based on Haar wavelet

Directory of Open Access Journals (Sweden)

Naresh Berwal

2012-11-01

Full Text Available Wavelet transform and wavelet analysis are powerful mathematical tools for many problems. Wavelet also can be applied in numerical analysis. In this paper, we apply Haar wavelet method to solve wave-like equation with initial and boundary conditions known. The fundamental idea of Haar wavelet method is to convert the differential equations into a group of algebraic equations, which involves a finite number or variables. The results and graph show that the proposed way is quite reasonable when compared to exact solution.

Integral equations and their applications

CERN Document Server

Rahman, M

2007-01-01

For many years, the subject of functional equations has held a prominent place in the attention of mathematicians. In more recent years this attention has been directed to a particular kind of functional equation, an integral equation, wherein the unknown function occurs under the integral sign. The study of this kind of equation is sometimes referred to as the inversion of a definite integral. While scientists and engineers can already choose from a number of books on integral equations, this new book encompasses recent developments including some preliminary backgrounds of formulations of integral equations governing the physical situation of the problems. It also contains elegant analytical and numerical methods, and an important topic of the variational principles. Primarily intended for senior undergraduate students and first year postgraduate students of engineering and science courses, students of mathematical and physical sciences will also find many sections of direct relevance. The book contains eig...

A Note on Iterative Refinement for Seminormal Equations

Czech Academy of Sciences Publication Activity Database

Rozložník, Miroslav; Smoktunowicz, A.; Kopal, J.

2014-01-01

Roč. 75, January (2014), s. 167-174 ISSN 0168-9274 R&D Projects: GA ČR(CZ) GAP108/11/0853 Grant - others:TUL(CZ) SGS 7822/115 Institutional support: RVO:67985807 Keywords : condition number * numerical stability * normal equations Subject RIV: BA - General Mathematics Impact factor: 1.221, year: 2014

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

Mathematics for Nonlinear Phenomena : Analysis and Computation : International Conference in honor of Professor Yoshikazu Giga on his 60th birthday

CERN Document Server

Jimbo, Shuichi

2017-01-01

This volume covers some of the most seminal research in the areas of mathematical analysis and numerical computation for nonlinear phenomena. Collected from the international conference held in honor of Professor Yoshikazu Giga’s 60th birthday, the featured research papers and survey articles discuss partial differential equations related to fluid mechanics,　electromagnetism, surface diffusion, and evolving interfaces. Specific focus is placed on topics such as the solvability of the Navier-Stokes equations and the regularity, stability, and symmetry of their solutions, analysis of a living fluid, stochastic effects and numerics for Maxwell’s equations, nonlinear heat equations in critical spaces, viscosity solutions describing various kinds of interfaces, numerics for evolving interfaces, and a hyperbolic obstacle problem. Also included in this volume are an introduction of Yoshikazu Giga’s extensive academic career and a long list of his published work. Students and researchers in mathematical analy...

Numerical study of traveling-wave solutions for the Camassa-Holm equation

International Nuclear Information System (INIS)

Kalisch, Henrik; Lenells, Jonatan

2005-01-01

We explore numerically different aspects of periodic traveling-wave solutions of the Camassa-Holm equation. In particular, the time evolution of some recently found new traveling-wave solutions and the interaction of peaked and cusped waves is studied

On the numerical evaluation of algebro-geometric solutions to integrable equations

International Nuclear Information System (INIS)

Kalla, C; Klein, C

2012-01-01

Physically meaningful periodic solutions to certain integrable partial differential equations are given in terms of multi-dimensional theta functions associated with real Riemann surfaces. Typical analytical problems in the numerical evaluation of these solutions are studied. In the case of hyperelliptic surfaces efficient algorithms exist even for almost degenerate surfaces. This allows the numerical study of solitonic limits. For general real Riemann surfaces, the choice of a homology basis adapted to the anti-holomorphic involution is important for a convenient formulation of the solutions and smoothness conditions. Since existing algorithms for algebraic curves produce a homology basis not related to automorphisms of the curve, we study symplectic transformations to an adapted basis and give explicit formulae for M-curves. As examples we discuss solutions of the Davey–Stewartson and the multi-component nonlinear Schrödinger equations

Numerical solution of the Schroedinger equation with a polynomial potential

International Nuclear Information System (INIS)

Campoy, G.; Palma, A.

1986-01-01

A numerical method for solving the Schroedinger equation for a potential expressed as a polynomial is proposed. The basic assumption relies on the asymptotic properties of the solution of this equation. It is possible to obtain the energies and the stationary state functions simultaneously. They analyze, in particular, the cases of the quartic anharmonic oscillator and a hydrogen atom perturbed by a quadratic term, obtaining its energy eigenvalues for some values of the perturbation parameter. Together with the Hellmann-Feynman theorem, they use their algorithm to calculate expectation values of x'' for arbitrary positive values of n. 4 tables

Numerical Solution of Multiterm Fractional Differential Equations Using the Matrix Mittag–Leffler Functions

Directory of Open Access Journals (Sweden)

Marina Popolizio

2018-01-01

Full Text Available Multiterm fractional differential equations (MTFDEs nowadays represent a widely used tool to model many important processes, particularly for multirate systems. Their numerical solution is then a compelling subject that deserves great attention, not least because of the difficulties to apply general purpose methods for fractional differential equations (FDEs to this case. In this paper, we first transform the MTFDEs into equivalent systems of FDEs, as done by Diethelm and Ford; in this way, the solution can be expressed in terms of Mittag–Leffler (ML functions evaluated at matrix arguments. We then propose to compute it by resorting to the matrix approach proposed by Garrappa and Popolizio. Several numerical tests are presented that clearly show that this matrix approach is very accurate and fast, also in comparison with other numerical methods.

Foundations of children's numerical and mathematical skills: the roles of symbolic and nonsymbolic representations of numerical magnitude.

Science.gov (United States)

Lyons, Ian M; Ansari, Daniel

2015-01-01

Numerical and mathematical skills are critical predictors of academic success. The last three decades have seen a substantial growth in our understanding of how the human mind and brain represent and process numbers. In particular, research has shown that we share with animals the ability to represent numerical magnitude (the total number of items in a set) and that preverbal infants can process numerical magnitude. Further research has shown that similar processing signatures characterize numerical magnitude processing across species and developmental time. These findings suggest that an approximate system for nonsymbolic (e.g., dot arrays) numerical magnitude representation serves as the basis for the acquisition of cultural, symbolic (e.g., Arabic numerals) representations of numerical magnitude. This chapter explores this hypothesis by reviewing studies that have examined the relation between individual differences in nonsymbolic numerical magnitude processing and symbolic math abilities (e.g., arithmetic). Furthermore, we examine the extent to which the available literature provides strong evidence for a link between symbolic and nonsymbolic representations of numerical magnitude at the behavioral and neural levels of analysis. We conclude that claims that symbolic number abilities are grounded in the approximate system for the nonsymbolic representation of numerical magnitude are not strongly supported by the available evidence. Alternative models and future research directions are discussed. © 2015 Elsevier Inc. All rights reserved.

Numerical implementation of the Dirac equation on hypercube multicomputers

International Nuclear Information System (INIS)

Wells, J.C.

1991-01-01

Motivated by an interest in nonperturbative electromagnetic lepton-pair production in relativistic heavy-ion collisions, we discuss the numerical methods used in implementing a lattice solution of the time-dependent Dirac equation in three-dimensional Cartesian coordinates. Discretization is obtained using the lattice basis-spline collocation method, in which quantum-state vectors and coordinate-space operators are expressed in terms of basis-spline functions, and represented on a spatial lattice. All numerical procedures reduce to a series of matrix-vector operations which we perform on the Intel iPSC/860 hypercube multicomputer. We discuss solutions to the problems of limited node memory and node-to-node communication overhead inherent in using distributed-memory, multiple-instruction, multiple-data parallel computers

The numerical dynamic for highly nonlinear partial differential equations

Science.gov (United States)

Lafon, A.; Yee, H. C.

1992-01-01

Problems associated with the numerical computation of highly nonlinear equations in computational fluid dynamics are set forth and analyzed in terms of the potential ranges of spurious behaviors. A reaction-convection equation with a nonlinear source term is employed to evaluate the effects related to spatial and temporal discretizations. The discretization of the source term is described according to several methods, and the various techniques are shown to have a significant effect on the stability of the spurious solutions. Traditional linearized stability analyses cannot provide the level of confidence required for accurate fluid dynamics computations, and the incorporation of nonlinear analysis is proposed. Nonlinear analysis based on nonlinear dynamical systems complements the conventional linear approach and is valuable in the analysis of hypersonic aerodynamics and combustion phenomena.

Problem-Based Instructional Strategy and Numerical Ability as Determinants of Senior Secondary Achievement in Mathematics

Science.gov (United States)

Badru, Ademola K.

2016-01-01

The study investigated Problem-based Instructional Strategy and Numerical ability as determinants of Senior Secondary Achievement in Mathematics. This study used 4 x 2 x 2 non-randomised control group Pretest-Posttest Quasi-experimental Factorial design. It consisted of two independent variables (treatment and Numerical ability) and one moderating…

Predicting Relationships between Mathematics Anxiety, Mathematics Teaching Anxiety, Self-Efficacy Beliefs towards Mathematics and Mathematics Teaching

Science.gov (United States)

Unlu, Melihan; Ertekin, Erhan; Dilmac, Bulent

2017-01-01

The purpose of the research is to investigate the relationships between self-efficacy beliefs toward mathematics, mathematics anxiety and self-efficacy beliefs toward mathematics teaching, mathematics teaching anxiety variables and testing the relationships between these variables with structural equation model. The sample of the research, which…

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

Numerical solution of the resistive magnetohydrodynamic boundary-layer equations

International Nuclear Information System (INIS)

Glasser, A.H.; Jardin, S.C.; Tesauro, G.

1983-10-01

Three different techniques are presented for numerical solution of the equations governing the boundary layer of resistive magnetohydrodynamic tearing and interchange instabilities in toroidal geometry. Excellent agreement among these methods and with analytical results provides confidence in the correctness of the results. Solutions obtained in regimes where analytical medthods fail indicate a new scaling for the tearing mode as well as the existence of a new regime of stability

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

Finite-element discretization of 3D energy-transport equations for semiconductors

Energy Technology Data Exchange (ETDEWEB)

Gadau, Stephan

2007-07-01

In this thesis a mathematical model was derived that describes the charge and energy transport in semiconductor devices like transistors. Moreover, numerical simulations of these physical processes are performed. In order to accomplish this, methods of theoretical physics, functional analysis, numerical mathematics and computer programming are applied. After an introduction to the status quo of semiconductor device simulation methods and a brief review of historical facts up to now, the attention is shifted to the construction of a model, which serves as the basis of the subsequent derivations in the thesis. Thereby the starting point is an important equation of the theory of dilute gases. From this equation the model equations are derived and specified by means of a series expansion method. This is done in a multi-stage derivation process, which is mainly taken from a scientific paper and which does not constitute the focus of this thesis. In the following phase we specify the mathematical setting and make precise the model assumptions. Thereby we make use of methods of functional analysis. Since the equations we deal with are coupled, we are concerned with a nonstandard problem. In contrary, the theory of scalar elliptic equations is established meanwhile. Subsequently, we are preoccupied with the numerical discretization of the equations. A special finite-element method is used for the discretization. This special approach has to be done in order to make the numerical results appropriate for practical application. By a series of transformations from the discrete model we derive a system of algebraic equations that are eligible for numerical evaluation. Using self-made computer programs we solve the equations to get approximate solutions. These programs are based on new and specialized iteration procedures that are developed and thoroughly tested within the frame of this research work. Due to their importance and their novel status, they are explained and

Numerical solution of modified fokker-planck equation with poissonian input

Czech Academy of Sciences Publication Activity Database

Náprstek, Jiří; Král, Radomil

2010-01-01

Roč. 17, 3/4 (2010), s. 251-268 ISSN 1802-1484 R&D Projects: GA AV ČR(CZ) IAA200710805; GA ČR(CZ) GA103/09/0094 Institutional research plan: CEZ:AV0Z20710524 Keywords : Fokker-Planck equation * poisson ian exciation * numerical solution * transition effects Subject RIV: JN - Civil Engineering

LIE GROUPS AND NUMERICAL SOLUTIONS OF DIFFERENTIAL EQUATIONS: INVARIANT DISCRETIZATION VERSUS DIFFERENTIAL APPROXIMATION

Directory of Open Access Journals (Sweden)

Decio Levi

2013-10-01

Full Text Available We briefly review two different methods of applying Lie group theory in the numerical solution of ordinary differential equations. On specific examples we show how the symmetry preserving discretization provides difference schemes for which the “first differential approximation” is invariant under the same Lie group as the original ordinary differential equation.

Numerical Solution of the Kzk Equation for Pulsed Finite Amplitude Sound Beams in Thermoviscous Fluids

Science.gov (United States)

Lee, Yang-Sub

A time-domain numerical algorithm for solving the KZK (Khokhlov-Zabolotskaya-Kuznetsov) nonlinear parabolic wave equation is developed for pulsed, axisymmetric, finite amplitude sound beams in thermoviscous fluids. The KZK equation accounts for the combined effects of diffraction, absorption, and nonlinearity at the same order of approximation. The accuracy of the algorithm is established via comparison with analytical solutions for several limiting cases, and with numerical results obtained from a widely used algorithm for solving the KZK equation in the frequency domain. The time domain algorithm is used to investigate waveform distortion and shock formation in directive sound beams radiated by pulsed circular piston sources. New results include predictions for the entire process of self-demodulation, and for the effect of frequency modulation on pulse envelope distortion. Numerical results are compared with measurements, and focused sources are investigated briefly.

Direct numerical solution of Poisson's equation in cylindrical (r, z) coordinates

International Nuclear Information System (INIS)

Chao, E.H.; Paul, S.F.; Davidson, R.C.; Fine, K.S.

1997-01-01

A direct solver method is developed for solving Poisson's equation numerically for the electrostatic potential φ(r,z) in a cylindrical region (r wall , 0 wall , z) are specified, and ∂φ/∂z = 0 at the axial boundaries (z = 0, L)

Mathematical and numerical methods for the nonlinear hyperbolic propagation problem: d2GAMMA/dt2 = d/dz [dGAMMA/dt dGAMMA/dz

International Nuclear Information System (INIS)

Stanco, L.; Vaccaro, V.G.; Funk, U.; Krueger, U.; Mika, K.; Wuestefeld, G.

1982-03-01

In the first part of this report a physical model is presented, which describes the deforming of a bunch in a storage ring influenced only by its own space charge field. A system of two differential equations for the density and the momentum of the particles is set up, which is independent of any special machine parameter. Due to the sign of the inductance of the chamber walls and the sign of the dispersion of the revolution frequency, we distinguish between a de-bunching and a self-bunching situation. The de-bunching corresponds to a nonlinear hyperbolic propagation problem well-known in gas dynamics, and the self-bunching to a nonlinear elliptic initial value problem. The second part deals with a mathematical and numerical treatment of an approximate equation for the hyperbolic case. For this nonlinear second order partial differential equation we first present three particular integrals: the solution by separating the variables, the similarity solution, and the solution for a parabolic initial distribution of the density. For a more realistic initial condition, we must resort to other methods: Results are obtained in three different ways, first from a highly accurate Taylor series expansion, second from a common finite difference method, and thirdly from the numerical method of characteristics. The appearance of a shock discontinuity is furthermore established in each of these cases. (orig.)

Simple equation method for nonlinear partial differential equations and its applications

Directory of Open Access Journals (Sweden)

Taher A. Nofal

2016-04-01

Full Text Available In this article, we focus on the exact solution of the some nonlinear partial differential equations (NLPDEs such as, Kodomtsev–Petviashvili (KP equation, the (2 + 1-dimensional breaking soliton equation and the modified generalized Vakhnenko equation by using the simple equation method. In the simple equation method the trial condition is the Bernoulli equation or the Riccati equation. It has been shown that the method provides a powerful mathematical tool for solving nonlinear wave equations in mathematical physics and engineering problems.

Asymptotic preserving error estimates for numerical solutions of compressible Navier-Stokes equations in the low Mach number regime

Czech Academy of Sciences Publication Activity Database

Feireisl, Eduard; Medviďová-Lukáčová, M.; Nečasová, Šárka; Novotný, A.; She, Bangwei

2018-01-01

Roč. 16, č. 1 (2018), s. 150-183 ISSN 1540-3459 R&D Projects: GA ČR GA16-03230S EU Projects: European Commission(XE) 320078 - MATHEF Institutional support: RVO:67985840 Keywords : Navier-Stokes system * finite element numerical method * finite volume numerical method * asymptotic preserving schemes Subject RIV: BA - General Mathematics OBOR OECD: Pure mathematics Impact factor: 1.865, year: 2016 http://epubs.siam.org/doi/10.1137/16M1094233

Numerical simulation of stochastic point kinetic equation in the dynamical system of nuclear reactor

International Nuclear Information System (INIS)

Saha Ray, S.

2012-01-01

Highlights: ► In this paper stochastic neutron point kinetic equations have been analyzed. ► Euler–Maruyama method and Strong Taylor 1.5 order method have been discussed. ► These methods are applied for the solution of stochastic point kinetic equations. ► Comparison between the results of these methods and others are presented in tables. ► Graphs for neutron and precursor sample paths are also presented. -- Abstract: In the present paper, the numerical approximation methods, applied to efficiently calculate the solution for stochastic point kinetic equations () in nuclear reactor dynamics, are investigated. A system of Itô stochastic differential equations has been analyzed to model the neutron density and the delayed neutron precursors in a point nuclear reactor. The resulting system of Itô stochastic differential equations are solved over each time-step size. The methods are verified by considering different initial conditions, experimental data and over constant reactivities. The computational results indicate that the methods are simple and suitable for solving stochastic point kinetic equations. In this article, a numerical investigation is made in order to observe the random oscillations in neutron and precursor population dynamics in subcritical and critical reactors.

Exponential and Bessel fitting methods for the numerical solution of the Schroedinger equation

International Nuclear Information System (INIS)

Raptis, A.D.; Cash, J.R.

1987-01-01

A new method is developed for the numerical integration of the one dimensional radial Schroedinger equation. This method involves using different integration formulae in different parts of the range of integration rather than using the same integration formula throughout. Two new integration formulae are derived, one which integrates Bessel and Neumann functions exactly and another which exactly integrates certain exponential functions. It is shown that, for large r, these new formulae are much more accurate than standard integration methods for the Schroedinger equation. The benefit of using this new approach is demonstrated by considering some numerical examples based on the Lennard-Jones potential. (orig.)

Trends in contemporary mathematics

CERN Document Server

Strickland, Elisabetta

2014-01-01

This book covers a wide spectrum of hot topics and current trends in mathematics, including noncommutative algebra via deformation theory,  optimal transportation, nonlinear potential theory, kinetic theory and gas dynamics, geometric numerical integration, finite simple groups of small essential dimension, optimal control problems, extended Dynkin diagrams, spin glasses, aspherical closed manifolds, Boltzmann systems, birational geometry of projective varieties and directed graphs, nonlinear diffusion, geometric constructions of extremal metrics on complex manifolds, and Pell’s equation in polynomials. The book comprises a selection of contributions by leading international mathematicians who were speakers at the "INdAM Day", an initiative dating back to 2004 at which the most recent developments in contemporary mathematics are presented.

Implementing a Flipped Classroom Approach in a University Numerical Methods Mathematics Course

Science.gov (United States)

Johnston, Barbara M.

2017-01-01

This paper describes and analyses the implementation of a "flipped classroom" approach, in an undergraduate mathematics course on numerical methods. The approach replaced all the lecture contents by instructor-made videos and was implemented in the consecutive years 2014 and 2015. The sequential case study presented here begins with an…

HYDRA-II: A hydrothermal analysis computer code: Volume 1, Equations and numerics

International Nuclear Information System (INIS)

McCann, R.A.

1987-04-01

HYDRA-II is a hydrothermal computer code capable of three-dimensional analysis of coupled conduction, convection, and thermal radiation problems. This code is especially appropriate for simulating the steady-state performance of spent fuel storage systems. The code has been evaluated for this application for the US Department of Energy's Commercial Spent Fuel Management Program. HYDRA-II provides a finite difference solution in Cartesian coordinates to the equations governing the conservation of mass, momentum, and energy. A cylindrical coordinate system may also be used to enclose the Cartesian coordinate system. This exterior coordinate system is useful for modeling cylindrical cask bodies. The difference equations for conservation of momentum are enhanced by the incorporation of directional porosities and permeabilities that aid in modeling solid structures whose dimensions may be smaller than the computational mesh. The equation for conservation of energy permits of modeling of orthotropic physical properties and film resistances. Several automated procedures are available to model radiation transfer within enclosures and from fuel rod to fuel rod. The documentation of HYDRA-II is presented in three separate volumes. This volume, Volume I - Equations and Numerics, describes the basic differential equations, illustrates how the difference equations are formulated, and gives the solution procedures employed. Volume II - User's Manual contains code flow charts, discusses the code structure, provides detailed instructions for preparing an input file, and illustrates the operation of the code by means of a model problem. The final volume, Volume III - Verification/Validation Assessments, presents results of numerical simulations of single- and multiassembly storage systems and comparisons with experimental data. 4 refs

Tree-grass interaction dynamics and pulsed fires : mathematical and numerical studies

OpenAIRE

Tamen, A. T.; Dumont, Y.; Tewa, J. J.; Bowong, S.; Couteron, Pierre

2016-01-01

Savannas are dynamical systems where grasses and trees can either dominate or coexist. Fires are known to be central in the functioning of the savanna biome although their characteristics are expected to vary along the rainfall gradients as observed in Sub-Saharan Africa. In this paper, we model the tree-grass dynamics using impulsive differential equations that consider fires as discrete events. This framework allows us to carry out a comprehensive qualitative mathematical analysis that reve...

Tree-Grass interactions dynamics and Pulse Fires: mathematical and numerical studies

OpenAIRE

Tamen, A. Tchuinté; Dumont, Y.; Bowong, S.; Tewa, J. J.; Couteron, P.

2015-01-01

Savannas are dynamical systems where grasses and trees can either dominate or coexist. Fires are known to be central in the functioning of the savanna biome though their characteristics are expected to vary along the rainfall gradients as observed in Sub-Saharan Africa. In this paper, we model the tree-grass dynamics using impulsive differential equations that consider fires as discrete events. This framework allows us to carry out a comprehensive qualitative mathematical analysis that reveal...

Mathematical modelling of steam generator and design of temperature regulator

Energy Technology Data Exchange (ETDEWEB)

Bogdanovic, S.S. [EE Institute Nikola Tesla, Belgrade (Yugoslavia)

1999-07-01

The paper considers mathematical modelling of once-through power station boiler and numerical algorithm for simulation of the model. Fast and numerically stable algorithm based on the linearisation of model equations and on the simultaneous solving of differential and algebraic equations is proposed. The paper also presents the design of steam temperature regulator by using the method of projective controls. Dynamic behaviour of the system closed with optimal linear quadratic regulator is taken as the reference system. The desired proprieties of the reference system are retained and solutions for superheated steam temperature regulator are determined. (author)

9th and 10th Asian Symposium on Computer Mathematics

CERN Document Server

Lee, Wen-shin; Sato, Yosuke

2014-01-01

This book covers original research and the latest advances in symbolic, algebraic and geometric computation; computational methods for differential and difference equations, symbolic-numerical computation; mathematics software design and implementation; and scientific and engineering applications based on features, invited talks, special sessions and contributed papers presented at the 9th (in Fukuoka, Japan in 2009) and 10th (in Beijing China in 2012) Asian Symposium on Computer Mathematics (ASCM). Thirty selected and refereed articles in the book present the conference participants’ ideas and views on researching mathematics using computers.

Finite difference method and algebraic polynomial interpolation for numerically solving Poisson's equation over arbitrary domains

Directory of Open Access Journals (Sweden)

Tsugio Fukuchi

2014-06-01

Full Text Available The finite difference method (FDM based on Cartesian coordinate systems can be applied to numerical analyses over any complex domain. A complex domain is usually taken to mean that the geometry of an immersed body in a fluid is complex; here, it means simply an analytical domain of arbitrary configuration. In such an approach, we do not need to treat the outer and inner boundaries differently in numerical calculations; both are treated in the same way. Using a method that adopts algebraic polynomial interpolations in the calculation around near-wall elements, all the calculations over irregular domains reduce to those over regular domains. Discretization of the space differential in the FDM is usually derived using the Taylor series expansion; however, if we use the polynomial interpolation systematically, exceptional advantages are gained in deriving high-order differences. In using the polynomial interpolations, we can numerically solve the Poisson equation freely over any complex domain. Only a particular type of partial differential equation, Poisson's equations, is treated; however, the arguments put forward have wider generality in numerical calculations using the FDM.

Numerical Oscillations Analysis for Nonlinear Delay Differential Equations in Physiological Control Systems

Directory of Open Access Journals (Sweden)

Qi Wang

2012-01-01

Full Text Available This paper deals with the oscillations of numerical solutions for the nonlinear delay differential equations in physiological control systems. The exponential θ-method is applied to p′(t=β0ωμp(t−τ/(ωμ+pμ(t−τ−γp(t and it is shown that the exponential θ-method has the same order of convergence as that of the classical θ-method. Several conditions under which the numerical solutions oscillate are derived. Moreover, it is proven that every nonoscillatory numerical solution tends to positive equilibrium of the continuous system. Finally, the main results are illustrated with numerical examples.

Modeling and numerical simulations of the influenced Sznajd model

Science.gov (United States)

Karan, Farshad Salimi Naneh; Srinivasan, Aravinda Ramakrishnan; Chakraborty, Subhadeep

2017-08-01

This paper investigates the effects of independent nonconformists or influencers on the behavioral dynamic of a population of agents interacting with each other based on the Sznajd model. The system is modeled on a complete graph using the master equation. The acquired equation has been numerically solved. Accuracy of the mathematical model and its corresponding assumptions have been validated by numerical simulations. Regions of initial magnetization have been found from where the system converges to one of two unique steady-state PDFs, depending on the distribution of influencers. The scaling property and entropy of the stationary system in presence of varying level of influence have been presented and discussed.

The Contributions of Domain-General and Numerical Factors to Third-Grade Arithmetic Skills and Mathematical Learning Disability

Science.gov (United States)

Cowan, Richard; Powell, Daisy

2014-01-01

Explanations of the marked individual differences in elementary school mathematical achievement and mathematical learning disability (MLD or dyscalculia) have involved domain-general factors (working memory, reasoning, processing speed, and oral language) and numerical factors that include single-digit processing efficiency and multidigit skills…

Numerical method of identification of an unknown source term in a heat equation

Directory of Open Access Journals (Sweden)

Fatullayev Afet Golayo?lu

2002-01-01

Full Text Available A numerical procedure for an inverse problem of identification of an unknown source in a heat equation is presented. Approach of proposed method is to approximate unknown function by polygons linear pieces which are determined consecutively from the solution of minimization problem based on the overspecified data. Numerical examples are presented.

Partial differential equations and their applications

International Nuclear Information System (INIS)

Gauthier-Villars

1998-01-01

This book is dedicated to the French mathematician J.L.Lions. It represents a compilation of articles from about 80 authors. The topics treated are diverse but the more or less commune matter is the study of the characteristics of some partial differential equations. Stability, optimal approximation, numerical resolution, particular applications are among the subjects reviewed. An article deals with the MHD stability of fusion plasmas in tokamaks, another presents the scientific and technical challenges of nuclear energy in France. The latter that contains no equations can be considered as an enjoyable break in a sea of about 40 mathematical articles. (A.C.)

Assessing numerical methods used in nuclear aerosol transport models

International Nuclear Information System (INIS)

McDonald, B.H.

1987-01-01

Several computer codes are in use for predicting the behaviour of nuclear aerosols released into containment during postulated accidents in water-cooled reactors. Each of these codes uses numerical methods to discretize and integrate the equations that govern the aerosol transport process. Computers perform only algebraic operations and generate only numbers. It is in the numerical methods that sense can be made of these numbers and where they can be related to the actual solution of the equations. In this report, the numerical methods most commonly used in the aerosol transport codes are examined as special cases of a general solution procedure, the Method of Weighted Residuals. It would appear that the numerical methods used in the codes are all capable of producing reasonable answers to the mathematical problem when used with skill and care. 27 refs

Mathematics anxiety reduces default mode network deactivation in response to numerical tasks

Directory of Open Access Journals (Sweden)

Belinda ePletzer

2015-04-01

Full Text Available Mathematics anxiety is negatively related to mathematics performance, thereby threatening the professional success. Preoccupation with the emotional content of the stimuli may consume working memory resources, which may be reflected in decreased deactivation of areas associated with the default mode network (DMN activated during self-referential and emotional processing. The common problem is that math anxiety is usually associated with poor math performance, so that any group differences are difficult to interpret.Here we compared the BOLD-response of 18 participants with high (HMAs and 18 participants with low mathematics anxiety (LMAs matched for their mathematical performance to two numerical tasks (number comparison, number bisection. During both tasks, we found stronger deactivation within the DMN in LMAs compared to HMAs, while BOLD-response in task-related activation areas did not differ between HMAs and LMAs. The difference in DMN deactivation between the HMA and LMA group was more pronounced in stimuli with additional requirement on inhibitory functions, but did not differ between number magnitude processing and arithmetic fact retrieval.

Mathematics anxiety reduces default mode network deactivation in response to numerical tasks.

Science.gov (United States)

Pletzer, Belinda; Kronbichler, Martin; Nuerk, Hans-Christoph; Kerschbaum, Hubert H

2015-01-01

Mathematics anxiety is negatively related to mathematics performance, thereby threatening the professional success. Preoccupation with the emotional content of the stimuli may consume working memory resources, which may be reflected in decreased deactivation of areas associated with the default mode network (DMN) activated during self-referential and emotional processing. The common problem is that math anxiety is usually associated with poor math performance, so that any group differences are difficult to interpret. Here we compared the BOLD-response of 18 participants with high (HMAs) and 18 participants with low mathematics anxiety (LMAs) matched for their mathematical performance to two numerical tasks (number comparison, number bisection). During both tasks, we found stronger deactivation within the DMN in LMAs compared to HMAs, while BOLD-response in task-related activation areas did not differ between HMAs and LMAs. The difference in DMN deactivation between the HMA and LMA group was more pronounced in stimuli with additional requirement on inhibitory functions, but did not differ between number magnitude processing and arithmetic fact retrieval.

Numerical analysis of critical two-phase flow in a convergent-divergent nozzle

International Nuclear Information System (INIS)

Romstedt, P.; Werner, W.

1985-01-01

The numerical calculation of critical two-phase flow in a convergent-divergent nozzle is complicated by a singularity of the fluid flow equations at the unknown critical point. This paper describes a method which is able to calculate critical state and its location without any additional assumptions. The critical state is identified by its mathematical properties: characteristics and solvability of linear systems with singular matrix. Because the numerically evaluable mathematical properties are only necessary conditions for the existence of critical flow, some physical ''compatibility-criteria'' (flow velocity equals two-phase sonic velocity, critical flow is independent of downstream flow state variations) are used as a substitute for mathematically sufficient conditions. Numerical results are shown for the critical flow in a LOBI nozzle; the two-phase flow is described by a model with equal phase velocities and thermodynamic non-equilibrium

Applied mathematics and condensed matter; Mathematiques appliquees et matiere condensee

Energy Technology Data Exchange (ETDEWEB)

Bouche, D.; Jollet, F. [CEA Bruyeres-le-Chatel, 91 (France)

2011-01-15

Applied mathematics have always been a key tool in computing the structure of condensed matter. In this paper, we present the most widely used methods, and show the importance of mathematics in their genesis and evolution. After a brief survey of quantum Monte Carlo methods, which try to compute the N electrons wave function, the paper describes the theoretical foundations of N independent particle approximations. We mainly focus on density functional theory (DFT). This theory associated with advanced numerical methods, and high performance computing, has produced significant achievements in the field. This paper presents the foundations of the theory, as well as different numerical methods used to solve DFT equations. (authors)

Solving ordinary differential equations by electrical analogy: a multidisciplinary teaching tool

Science.gov (United States)

Sanchez Perez, J. F.; Conesa, M.; Alhama, I.

2016-11-01

Ordinary differential equations are the mathematical formulation for a great variety of problems in science and engineering, and frequently, two different problems are equivalent from a mathematical point of view when they are formulated by the same equations. Students acquire the knowledge of how to solve these equations (at least some types of them) using protocols and strict algorithms of mathematical calculation without thinking about the meaning of the equation. The aim of this work is that students learn to design network models or circuits in this way; with simple knowledge of them, students can establish the association of electric circuits and differential equations and their equivalences, from a formal point of view, that allows them to associate knowledge of two disciplines and promote the use of this interdisciplinary approach to address complex problems. Therefore, they learn to use a multidisciplinary tool that allows them to solve these kinds of equations, even students of first course of engineering, whatever the order, grade or type of non-linearity. This methodology has been implemented in numerous final degree projects in engineering and science, e.g., chemical engineering, building engineering, industrial engineering, mechanical engineering, architecture, etc. Applications are presented to illustrate the subject of this manuscript.

A numerical scheme using multi-shockpeakons to compute solutions of the Degasperis-Procesi equation

Directory of Open Access Journals (Sweden)

Hakon A. Hoel

2007-07-01

Full Text Available We consider a numerical scheme for entropy weak solutions of the DP (Degasperis-Procesi equation $u_t - u_{xxt} + 4uu_x = 3u_{x}u_{xx}+ uu_{xxx}$. Multi-shockpeakons, functions of the form $$ u(x,t =sum_{i=1}^n(m_i(t -hbox{sign}(x-x_i(ts_i(te^{-|x-x_i(t|}, $$ are solutions of the DP equation with a special property; their evolution in time is described by a dynamical system of ODEs. This property makes multi-shockpeakons relatively easy to simulate numerically. We prove that if we are given a non-negative initial function $u_0 in L^1(mathbb{R}cap BV(mathbb{R}$ such that $u_{0} - u_{0,x}$ is a positive Radon measure, then one can construct a sequence of multi-shockpeakons which converges to the unique entropy weak solution in $mathbb{R}imes[0,T$ for any $T>0$. From this convergence result, we construct a multi-shockpeakon based numerical scheme for solving the DP equation.

Numerical solution of the radionuclide transport equation

International Nuclear Information System (INIS)

Hadermann, J.; Roesel, F.

1983-11-01

A numerical solution of the one-dimensional geospheric radionuclide chain transport equation based on the pseudospectral method is developed. The advantages of this approach are flexibility in incorporating space and time dependent migration parameters, arbitrary boundary conditions and solute rock interactions as well as efficiency and reliability. As an application the authors investigate the impact of non-linear sorption isotherms on migration in crystalline rock. It is shown that non-linear sorption, in the present case a Freundlich isotherm, may reduce concentration at the geosphere outlet by orders of magnitude provided the migration time is comparable or larger than the half-life of the nuclide in question. The importance of fixing dispersivity within the continuum approach is stressed. (Auth.)

Transition behaviours in two coupled Josephson junction equations

International Nuclear Information System (INIS)

Wang Jiazeng; Zhang Xuejuan; You Gongqiang; Zhou Fengyan

2007-01-01

The dynamics of two coupled Josephson junction equations are investigated via mathematical reasoning and numerical simulations. We show that for a fixed coupling K, the whole parameter space can be comparted into three regions: a quenching region, a synchronized running periodic region and a region where these two states coexist. It is further shown that with the increase of the coupling K, the system may transit from a synchronizing state to a quenching state. The characteristic of the critical line K*(b) which separates these two states is mathematically analysed

Numerical bifurcation analysis of delay differential equations arising from physiological modeling.

Science.gov (United States)

Engelborghs, K; Lemaire, V; Bélair, J; Roose, D

2001-04-01

This paper has a dual purpose. First, we describe numerical methods for continuation and bifurcation analysis of steady state solutions and periodic solutions of systems of delay differential equations with an arbitrary number of fixed, discrete delays. Second, we demonstrate how these methods can be used to obtain insight into complex biological regulatory systems in which interactions occur with time delays: for this, we consider a system of two equations for the plasma glucose and insulin concentrations in a diabetic patient subject to a system of external assistance. The model has two delays: the technological delay of the external system, and the physiological delay of the patient's liver. We compute stability of the steady state solution as a function of two parameters, compare with analytical results and compute several branches of periodic solutions and their stability. These numerical results allow to infer two categories of diabetic patients for which the external system has different efficiency.

Mathematics without boundaries surveys in pure mathematics

CERN Document Server

Pardalos, Panos

2014-01-01

The contributions in this volume have been written by eminent scientists from the international mathematical community and present significant advances in several theories, methods and problems of Mathematical Analysis, Discrete Mathematics, Geometry and their Applications. The chapters focus on both old and recent developments in Functional Analysis, Harmonic Analysis, Complex Analysis, Operator Theory, Combinatorics, Functional Equations, Differential Equations as well as a variety of Applications. The book also contains some review works, which could prove particularly useful for a broader audience of readers in Mathematical Sciences, and especially to graduate students looking for the  latest information.

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

The Role of Numerical Methods in the Sensitivity Analysis of a ...

African Journals Online (AJOL)

The mathematical modelling of physiochemical interaction in the framework of industrial and environmental physics which relies on an initial value problem is defined by a first order ordinary differential equation. Two numerical methods of studying sensitivity analysis of physiochemical interaction data are developed.

A numerical solution for a class of time fractional diffusion equations with delay

Directory of Open Access Journals (Sweden)

Pimenov Vladimir G.

2017-09-01

Full Text Available This paper describes a numerical scheme for a class of fractional diffusion equations with fixed time delay. The study focuses on the uniqueness, convergence and stability of the resulting numerical solution by means of the discrete energy method. The derivation of a linearized difference scheme with convergence order O(τ2−α+ h4 in L∞-norm is the main purpose of this study. Numerical experiments are carried out to support the obtained theoretical results.

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

Efficient Numerical Methods for Stochastic Differential Equations in Computational Finance

KAUST Repository

Happola, Juho

2017-09-19

Stochastic Differential Equations (SDE) offer a rich framework to model the probabilistic evolution of the state of a system. Numerical approximation methods are typically needed in evaluating relevant Quantities of Interest arising from such models. In this dissertation, we present novel effective methods for evaluating Quantities of Interest relevant to computational finance when the state of the system is described by an SDE.

Efficient Numerical Methods for Stochastic Differential Equations in Computational Finance

KAUST Repository

Happola, Juho

2017-01-01

Stochastic Differential Equations (SDE) offer a rich framework to model the probabilistic evolution of the state of a system. Numerical approximation methods are typically needed in evaluating relevant Quantities of Interest arising from such models. In this dissertation, we present novel effective methods for evaluating Quantities of Interest relevant to computational finance when the state of the system is described by an SDE.

Applied partial differential equations

CERN Document Server

Logan, J David

2015-01-01

This text presents the standard material usually covered in a one-semester, undergraduate course on boundary value problems and PDEs.  Emphasis is placed on motivation, concepts, methods, and interpretation, rather than on formal theory. The concise treatment of the subject is maintained in this third edition covering all the major ideas: the wave equation, the diffusion equation, the Laplace equation, and the advection equation on bounded and unbounded domains. Methods include eigenfunction expansions, integral transforms, and characteristics. In this third edition, text remains intimately tied to applications in heat transfer, wave motion, biological systems, and a variety other topics in pure and applied science. The text offers flexibility to instructors who, for example, may wish to insert topics from biology or numerical methods at any time in the course. The exposition is presented in a friendly, easy-to-read, style, with mathematical ideas motivated from physical problems. Many exercises and worked e...

Stratified flows with variable density: mathematical modelling and numerical challenges.

Science.gov (United States)

Murillo, Javier; Navas-Montilla, Adrian

2017-04-01

Stratified flows appear in a wide variety of fundamental problems in hydrological and geophysical sciences. They may involve from hyperconcentrated floods carrying sediment causing collapse, landslides and debris flows, to suspended material in turbidity currents where turbulence is a key process. Also, in stratified flows variable horizontal density is present. Depending on the case, density varies according to the volumetric concentration of different components or species that can represent transported or suspended materials or soluble substances. Multilayer approaches based on the shallow water equations provide suitable models but are not free from difficulties when moving to the numerical resolution of the governing equations. Considering the variety of temporal and spatial scales, transfer of mass and energy among layers may strongly differ from one case to another. As a consequence, in order to provide accurate solutions, very high order methods of proved quality are demanded. Under these complex scenarios it is necessary to observe that the numerical solution provides the expected order of accuracy but also converges to the physically based solution, which is not an easy task. To this purpose, this work will focus in the use of Energy balanced augmented solvers, in particular, the Augmented Roe Flux ADER scheme. References: J. Murillo , P. García-Navarro, Wave Riemann description of friction terms in unsteady shallow flows: Application to water and mud/debris floods. J. Comput. Phys. 231 (2012) 1963-2001. J. Murillo B. Latorre, P. García-Navarro. A Riemann solver for unsteady computation of 2D shallow flows with variable density. J. Comput. Phys.231 (2012) 4775-4807. A. Navas-Montilla, J. Murillo, Energy balanced numerical schemes with very high order. The Augmented Roe Flux ADER scheme. Application to the shallow water equations, J. Comput. Phys. 290 (2015) 188-218. A. Navas-Montilla, J. Murillo, Asymptotically and exactly energy balanced augmented flux

Numerical study of a Vlasov equation for systems with interacting particles

Energy Technology Data Exchange (ETDEWEB)

Herrera, Dianela; Curilef, Sergio [Departamento de Física, Universidad Católica del Norte, Avenida Angamos 0610, Antofagasta (Chile)

2015-03-10

We solve numerically the Vlasov equation for the self-gravitating sheet model. We used the method introduced by Cheng and Knorr [Comput Phys 22, 330-351 (1976)]. We discuss the quasi-stationary state for some thermodynamical observables, specifically the kinetic energy, whose trend is depicted for early evolution.

Numerical solution of the Navier-Stokes equations by discontinuous Galerkin method