126,937 research outputs found
Status of the differential transformation method
Further to a recent controversy on whether the differential transformation
method (DTM) for solving a differential equation is purely and solely the
traditional Taylor series method, it is emphasized that the DTM is currently
used, often only, as a technique for (analytically) calculating the power
series of the solution (in terms of the initial value parameters). Sometimes, a
piecewise analytic continuation process is implemented either in a numerical
routine (e.g., within a shooting method) or in a semi-analytical procedure
(e.g., to solve a boundary value problem). Emphasized also is the fact that, at
the time of its invention, the currently-used basic ingredients of the DTM
(that transform a differential equation into a difference equation of same
order that is iteratively solvable) were already known for a long time by the
"traditional"-Taylor-method users (notably in the elaboration of software
packages --numerical routines-- for automatically solving ordinary differential
equations). At now, the defenders of the DTM still ignore the, though much
better developed, studies of the "traditional"-Taylor-method users who, in
turn, seem to ignore similarly the existence of the DTM. The DTM has been given
an apparent strong formalization (set on the same footing as the Fourier,
Laplace or Mellin transformations). Though often used trivially, it is easily
attainable and easily adaptable to different kinds of differentiation
procedures. That has made it very attractive. Hence applications to various
problems of the Taylor method, and more generally of the power series method
(including noninteger powers) has been sketched. It seems that its potential
has not been exploited as it could be. After a discussion on the reasons of the
"misunderstandings" which have caused the controversy, the preceding topics are
concretely illustrated.Comment: To appear in Applied Mathematics and Computation, 29 pages,
references and further considerations adde
Multi-Dimensional Astrophysical Structural and Dynamical Analysis I. Development of a Nonlinear Finite Element Approach
A new field of numerical astrophysics is introduced which addresses the
solution of large, multidimensional structural or slowly-evolving problems
(rotating stars, interacting binaries, thick advective accretion disks, four
dimensional spacetimes, etc.). The technique employed is the Finite Element
Method (FEM), commonly used to solve engineering structural problems. The
approach developed herein has the following key features:
1. The computational mesh can extend into the time dimension, as well as
space, perhaps only a few cells, or throughout spacetime.
2. Virtually all equations describing the astrophysics of continuous media,
including the field equations, can be written in a compact form similar to that
routinely solved by most engineering finite element codes.
3. The transformations that occur naturally in the four-dimensional FEM
possess both coordinate and boost features, such that
(a) although the computational mesh may have a complex, non-analytic,
curvilinear structure, the physical equations still can be written in a simple
coordinate system independent of the mesh geometry.
(b) if the mesh has a complex flow velocity with respect to coordinate space,
the transformations will form the proper arbitrary Lagrangian- Eulerian
advective derivatives automatically.
4. The complex difference equations on the arbitrary curvilinear grid are
generated automatically from encoded differential equations.
This first paper concentrates on developing a robust and widely-applicable
set of techniques using the nonlinear FEM and presents some examples.Comment: 28 pages, 9 figures; added integral boundary conditions, allowing
very rapidly-rotating stars; accepted for publication in Ap.
Can Computer Algebra be Liberated from its Algebraic Yoke ?
So far, the scope of computer algebra has been needlessly restricted to exact
algebraic methods. Its possible extension to approximate analytical methods is
discussed. The entangled roles of functional analysis and symbolic programming,
especially the functional and transformational paradigms, are put forward. In
the future, algebraic algorithms could constitute the core of extended symbolic
manipulation systems including primitives for symbolic approximations.Comment: 8 pages, 2-column presentation, 2 figure
Inverse heat conduction problems by using particular solutions
Based on the method of fundamental solutions, we develop in this paper a new computational method to solve two-dimensional transient heat conduction inverse problems. The main idea is to use particular solutions as radial basis functions (PSRBF) for approximation of the solutions to the inverse heat conduction problems. The heat conduction equations are first analyzed in the Laplace transformed domain and the Durbin inversion method is then used to determine the solutions in the time domain. Least-square and singular value decomposition (SVD) techniques are adopted to solve the ill-conditioned linear system of algebraic equations obtained from the proposed PSRBF method. To demonstrate the effectiveness and simplicity of this approach, several numerical examples are given with satisfactory accuracy and stability.Peer reviewe
Spectral methods in general relativistic astrophysics
We present spectral methods developed in our group to solve three-dimensional
partial differential equations. The emphasis is put on equations arising from
astrophysical problems in the framework of general relativity.Comment: 51 pages, elsart (Elsevier Preprint), 19 PostScript figures,
submitted to Journal of Computational & Applied Mathematic
Developing Clean Technology through Approximate Solutions of Mathematical Models
In this paper, the role of mathematical modeling in the development of clean technology has been considered.
One method each for obtaining approximate solutions of mathematical models by ordinary differential equations
and partial differential equations respectively arising from the modeling of systems and physical phenomena has
been considered. The construction of continuous hybrid methods for the numerical approximation of the solutions
of initial value problems of ordinary differential equations as well as homotopy analysis method, an approximate
analytical method, for the solution of nonlinear partial differential equations are discussed
- …