373 research outputs found
An approximate solution of the MHD Falkner-Skan flow by Hermite functions pseudospectral method
Based on a new approximation method, namely pseudospectral method, a solution
for the three order nonlinear ordinary differential laminar boundary layer
Falkner-Skan equation has been obtained on the semi-infinite domain. The
proposed approach is equipped by the orthogonal Hermite functions that have
perfect properties to achieve this goal. This method solves the problem on the
semi-infinite domain without truncating it to a finite domain and transforming
domain of the problem to a finite domain. In addition, this method reduces
solution of the problem to solution of a system of algebraic equations. We also
present the comparison of this work with numerical results and show that the
present method is applicable.Comment: 15 pages, 4 figures; Published online in the journal of
"Communications in Nonlinear Science and Numerical Simulation
An approximation algorithm for the solution of the nonlinear Lane-Emden type equations arising in astrophysics using Hermite functions collocation method
In this paper we propose a collocation method for solving some well-known
classes of Lane-Emden type equations which are nonlinear ordinary differential
equations on the semi-infinite domain. They are categorized as singular initial
value problems. The proposed approach is based on a Hermite function
collocation (HFC) method. To illustrate the reliability of the method, some
special cases of the equations are solved as test examples. The new method
reduces the solution of a problem to the solution of a system of algebraic
equations. Hermite functions have prefect properties that make them useful to
achieve this goal. We compare the present work with some well-known results and
show that the new method is efficient and applicable.Comment: 34 pages, 13 figures, Published in "Computer Physics Communications
Numerical approximations for population growth model by Rational Chebyshev and Hermite Functions collocation approach: A comparison
This paper aims to compare rational Chebyshev (RC) and Hermite functions (HF)
collocation approach to solve the Volterra's model for population growth of a
species within a closed system. This model is a nonlinear integro-differential
equation where the integral term represents the effect of toxin. This approach
is based on orthogonal functions which will be defined. The collocation method
reduces the solution of this problem to the solution of a system of algebraic
equations. We also compare these methods with some other numerical results and
show that the present approach is applicable for solving nonlinear
integro-differential equations.Comment: 18 pages, 5 figures; Published online in the journal of "Mathematical
Methods in the Applied Sciences
Concurrent Learning Adaptive Model Predictive Control with Pseudospectral Implementation
This paper presents a control architecture in which a direct adaptive control
technique is used within the model predictive control framework, using the
concurrent learning based approach, to compensate for model uncertainties. At
each time step, the control sequences and the parameter estimates are both used
as the optimization arguments, thereby undermining the need for switching
between the learning phase and the control phase, as is the case with
hybrid-direct-indirect control architectures. The state derivatives are
approximated using pseudospectral methods, which are vastly used for numerical
optimal control problems. Theoretical results and numerical simulation examples
are used to establish the effectiveness of the architecture.Comment: 21 pages, 13 figure
Discontinuous Galerkin methods for general-relativistic hydrodynamics: formulation and application to spherically symmetric spacetimes
We have developed the formalism necessary to employ the
discontinuous-Galerkin approach in general-relativistic hydrodynamics. The
formalism is firstly presented in a general 4-dimensional setting and then
specialized to the case of spherical symmetry within a 3+1 splitting of
spacetime. As a direct application, we have constructed a one-dimensional code,
EDGES, which has been used to asses the viability of these methods via a series
of tests involving highly relativistic flows in strong gravity. Our results
show that discontinuous Galerkin methods are able not only to handle strong
relativistic shock waves but, at the same time, to attain very high orders of
accuracy and exponential convergence rates in smooth regions of the flow. Given
these promising prospects and their affinity with a pseudospectral solution of
the Einstein equations, discontinuous Galerkin methods could represent a new
paradigm for the accurate numerical modelling in relativistic astrophysics.Comment: 24 pages, 19 figures. Small changes; matches version to appear in PR
Spectral Methods for Numerical Relativity. The Initial Data Problem
Numerical relativity has traditionally been pursued via finite differencing.
Here we explore pseudospectral collocation (PSC) as an alternative to finite
differencing, focusing particularly on the solution of the Hamiltonian
constraint (an elliptic partial differential equation) for a black hole
spacetime with angular momentum and for a black hole spacetime superposed with
gravitational radiation. In PSC, an approximate solution, generally expressed
as a sum over a set of orthogonal basis functions (e.g., Chebyshev
polynomials), is substituted into the exact system of equations and the
residual minimized. For systems with analytic solutions the approximate
solutions converge upon the exact solution exponentially as the number of basis
functions is increased. Consequently, PSC has a high computational efficiency:
for solutions of even modest accuracy we find that PSC is substantially more
efficient, as measured by either execution time or memory required, than finite
differencing; furthermore, these savings increase rapidly with increasing
accuracy. The solution provided by PSC is an analytic function given
everywhere; consequently, no interpolation operators need to be defined to
determine the function values at intermediate points and no special
arrangements need to be made to evaluate the solution or its derivatives on the
boundaries. Since the practice of numerical relativity by finite differencing
has been, and continues to be, hampered by both high computational resource
demands and the difficulty of formulating acceptable finite difference
alternatives to the analytic boundary conditions, PSC should be further pursued
as an alternative way of formulating the computational problem of finding
numerical solutions to the field equations of general relativity.Comment: 15 pages, 5 figures, revtex, submitted to PR
- …