5,046 research outputs found
Symbolic computation of exact solutions expressible in hyperbolic and elliptic functions for nonlinear PDEs
Algorithms are presented for the tanh- and sech-methods, which lead to
closed-form solutions of nonlinear ordinary and partial differential equations
(ODEs and PDEs). New algorithms are given to find exact polynomial solutions of
ODEs and PDEs in terms of Jacobi's elliptic functions.
For systems with parameters, the algorithms determine the conditions on the
parameters so that the differential equations admit polynomial solutions in
tanh, sech, combinations thereof, Jacobi's sn or cn functions. Examples
illustrate key steps of the algorithms.
The new algorithms are implemented in Mathematica. The package
DDESpecialSolutions.m can be used to automatically compute new special
solutions of nonlinear PDEs. Use of the package, implementation issues, scope,
limitations, and future extensions of the software are addressed.
A survey is given of related algorithms and symbolic software to compute
exact solutions of nonlinear differential equations.Comment: 39 pages. Software available from Willy Hereman's home page at
http://www.mines.edu/fs_home/whereman
Statistics of non-linear stochastic dynamical systems under L\'evy noises by a convolution quadrature approach
This paper describes a novel numerical approach to find the statistics of the
non-stationary response of scalar non-linear systems excited by L\'evy white
noises. The proposed numerical procedure relies on the introduction of an
integral transform of Wiener-Hopf type into the equation governing the
characteristic function. Once this equation is rewritten as partial
integro-differential equation, it is then solved by applying the method of
convolution quadrature originally proposed by Lubich, here extended to deal
with this particular integral transform. The proposed approach is relevant for
two reasons: 1) Statistics of systems with several different drift terms can be
handled in an efficient way, independently from the kind of white noise; 2) The
particular form of Wiener-Hopf integral transform and its numerical evaluation,
both introduced in this study, are generalizations of fractional
integro-differential operators of potential type and Gr\"unwald-Letnikov
fractional derivatives, respectively.Comment: 20 pages, 5 figure
A space-time pseudospectral discretization method for solving diffusion optimal control problems with two-sided fractional derivatives
We propose a direct numerical method for the solution of an optimal control
problem governed by a two-side space-fractional diffusion equation. The
presented method contains two main steps. In the first step, the space variable
is discretized by using the Jacobi-Gauss pseudospectral discretization and, in
this way, the original problem is transformed into a classical integer-order
optimal control problem. The main challenge, which we faced in this step, is to
derive the left and right fractional differentiation matrices. In this respect,
novel techniques for derivation of these matrices are presented. In the second
step, the Legendre-Gauss-Radau pseudospectral method is employed. With these
two steps, the original problem is converted into a convex quadratic
optimization problem, which can be solved efficiently by available methods. Our
approach can be easily implemented and extended to cover fractional optimal
control problems with state constraints. Five test examples are provided to
demonstrate the efficiency and validity of the presented method. The results
show that our method reaches the solutions with good accuracy and a low CPU
time.Comment: This is a preprint of a paper whose final and definite form is with
'Journal of Vibration and Control', available from
[http://journals.sagepub.com/home/jvc]. Submitted 02-June-2018; Revised
03-Sept-2018; Accepted 12-Oct-201
Diffusive approximation of a time-fractional Burger's equation in nonlinear acoustics
A fractional time derivative is introduced into the Burger's equation to
model losses of nonlinear waves. This term amounts to a time convolution
product, which greatly penalizes the numerical modeling. A diffusive
representation of the fractional derivative is adopted here, replacing this
nonlocal operator by a continuum of memory variables that satisfy local-in-time
ordinary differential equations. Then a quadrature formula yields a system of
local partial differential equations, well-suited to numerical integration. The
determination of the quadrature coefficients is crucial to ensure both the
well-posedness of the system and the computational efficiency of the diffusive
approximation. For this purpose, optimization with constraint is shown to be a
very efficient strategy. Strang splitting is used to solve successively the
hyperbolic part by a shock-capturing scheme, and the diffusive part exactly.
Numerical experiments are proposed to assess the efficiency of the numerical
modeling, and to illustrate the effect of the fractional attenuation on the
wave propagation.Comment: submitted to Siam SIA
The LifeV library: engineering mathematics beyond the proof of concept
LifeV is a library for the finite element (FE) solution of partial
differential equations in one, two, and three dimensions. It is written in C++
and designed to run on diverse parallel architectures, including cloud and high
performance computing facilities. In spite of its academic research nature,
meaning a library for the development and testing of new methods, one
distinguishing feature of LifeV is its use on real world problems and it is
intended to provide a tool for many engineering applications. It has been
actually used in computational hemodynamics, including cardiac mechanics and
fluid-structure interaction problems, in porous media, ice sheets dynamics for
both forward and inverse problems. In this paper we give a short overview of
the features of LifeV and its coding paradigms on simple problems. The main
focus is on the parallel environment which is mainly driven by domain
decomposition methods and based on external libraries such as MPI, the Trilinos
project, HDF5 and ParMetis.
Dedicated to the memory of Fausto Saleri.Comment: Review of the LifeV Finite Element librar
Two Reliable Methods for Solving the (3 + 1)-Dimensional Space-Time Fractional Jimbo-Miwa Equation
We investigate methods for obtaining exact solutions of the (3 + 1)-dimensional nonlinear space-time fractional Jimbo-Miwa equation in the sense of the modified Riemann-Liouville derivative. The methods employed to analytically solve the equation are the G′/G,1/G-expansion method and the novel G′/G-expansion method. To the best of our knowledge, there are no researchers who have applied these methods to obtain exact solutions of the equation. The application of the methods is simple, elegant, efficient, and trustworthy. In particular, applying the novel G′/G-expansion method to the equation, we obtain more exact solutions than using other existing methods such as the G′/G-expansion method and the exp-Φ(ξ)-expansion method. The exact solutions of the equation, obtained using the two methods, can be categorized in terms of hyperbolic, trigonometric, and rational functions. Some of the results obtained by the two methods are new and reported here for the first time. In addition, the obtained exact explicit solutions of the equation characterize many physical meanings such as soliton solitary wave solutions, periodic wave solutions, and singular multiple-soliton solutions
A Shifted Jacobi-Gauss Collocation Scheme for Solving Fractional Neutral Functional-Differential Equations
The shifted Jacobi-Gauss collocation (SJGC) scheme is proposed and implemented to solve the fractional neutral functional-differential equations with proportional delays. The technique we have proposed is based upon shifted Jacobi polynomials with the Gauss quadrature integration technique. The main advantage of the shifted Jacobi-Gauss scheme is to reduce solving the generalized fractional neutral functional-differential equations to a system of algebraic equations in the unknown expansion. Reasonable numerical results are achieved by choosing few shifted Jacobi-Gauss collocation nodes. Numerical results demonstrate the accuracy, and versatility of the proposed algorithm
- …