Systems-based decomposition schemes for the approximate solution of multi-term fractional differential equations

We give a comparison of the efficiency of three alternative decomposition schemes for the approximate solution of multi-term fractional differential equations using the Caputo form of the fractional derivative. The schemes we compare are based on conversion of the original problem into a system of equations. We review alternative approaches and consider how the most appropriate numerical scheme may be chosen to solve a particular equation

Continuous approximations of a class of piece-wise continuous systems

In this paper we provide a rigorous mathematical foundation for continuous approximations of a class of systems with piece-wise continuous functions. By using techniques from the theory of differential inclusions, the underlying piece-wise functions can be locally or globally approximated. The approximation results can be used to model piece-wise continuous-time dynamical systems of integer or fractional-order. In this way, by overcoming the lack of numerical methods for diffrential equations of fractional-order with discontinuous right-hand side, unattainable procedures for systems modeled by this kind of equations, such as chaos control, synchronization, anticontrol and many others, can be easily implemented. Several examples are presented and three comparative applications are studied.Comment: IJBC, accepted (examples revised

Analysis, simulation and design of nonlinear RF circuits

The PhD project consists of two parts. The first part concerns the development of Computer Aided Design (CAD) algorithms for high-frequency circuits. Novel Padébased algorithms for numerical integration of ODEs as arise in high-frequency circuits are proposed. Both single- and multi-step methods are introduced. A large part of this section of the research is concerned with the application of Filon-type integration techniques to circuits subject to modulated signals. Such methods are tested with analog and digital modulated signals and are seen to be very effective. The results confirm that these methods are more accurate than the traditional trapezoidal rule and Runge-Kutta methods. The second part of the research is concerned with the analysis, simulation and design of RF circuits with emphasis on injection-locked frequency dividers (ILFD) and digital delta-sigma modulators (DDSM). Both of these circuits are employed in fractional-N frequency synthesizers. Several simulation methods are proposed to capture the locking range of an ILFD, such as the Warped Multi-time Partial Differential Equation (WaMPDE) and the Multiple-Phase-Condition Envelope Following (MPCENV) methods. The MPCENV method is the more efficient and accurate simulation technique and it is recommended to obviate the need for expensive experiments. The Multi-stAge noise Shaping (MASH) digital delta-sigma modulator (DDSM) is simulated in MATLAB and analysed mathematically. A novel structure employing multimoduli, termed the MM-MASH, is proposed. The goal in this design work is to reduce the noise level in the useful frequency band of the modulator. The success of the novel structure in achieving this aim is confirmed with simulations

All-propulsion design of the drag-free and attitude control of the European satellite GOCE

This paper concerns the drag-free and attitude control (DFAC) of the European Gravity field and steady-state Ocean Circulation Explorer satellite (GOCE), during the science phase. GOCE aims to determine the Earth's gravity field with high accuracy and spatial resolution, through complementary space techniques such as gravity gradiometry and precise orbit determination. Both techniques rely on accurate attitude and drag-free control, especially in the gradiometer measurement bandwidth (5-100mHz), where non-gravitational forces must be counteracted down to micronewton, and spacecraft attitude must track the local orbital reference frame with micro-radian accuracy. DFAC aims to enable the gravity gradiometer to operate so as to determine the Earth's gravity field especially in the so-called measurement bandwidth (5-100mHz), making use of ion and micro-thruster actuators. The DFAC unit has been designed entirely on a simplified discrete-time model (Embedded Model) derived from the fine dynamics of the spacecraft and its environment; the relevant control algorithms are implemented and tuned around the Embedded Model, which is the core of the control unit. The DFAC has been tested against uncertainties in spacecraft and environment and its code has been the preliminary model for final code development. The DFAC assumes an all-propulsion command authority, partly abandoned by the actual GOCE control system because of electric micro-propulsion not being fully developed. Since all-propulsion authority is expected to be imperative for future scientific and observation missions, design and simulated results are believed to be of interest to the space communit

An efficient numerical approach for solving two-point fractional order nonlinear boundary value problems with Robin boundary conditions

This article proposes new strategies for solving two-point Fractional order Nonlinear Boundary Value Problems (FNBVPs) with Robin Boundary Conditions (RBCs). In the new numerical schemes, a two-point FNBVP is transformed into a system of Fractional order Initial Value Problems (FIVPs) with unknown Initial Conditions (ICs). To approximate ICs in the system of FIVPs, we develop nonlinear shooting methods based on Newton's method and Halley's method using the RBC at the right end point. To deal with FIVPs in a system, we mainly employ High-order Predictor-Corrector Methods (HPCMs) with linear interpolation and quadratic interpolation (Nguyen and Jang in Fract. Calc. Appl. Anal. 20(2):447-476, 2017) into Volterra integral equations which are equivalent to FIVPs. The advantage of the proposed schemes with HPCMs is that even though they are designed for solving two-point FNBVPs, they can handle both linear and nonlinear two-point Fractional order Boundary Value Problems (FBVPs) with RBCs and have uniform convergence rates of HPCMs, O(h(2)) and O(h(3)) for shooting techniques with Newton's method and Halley's method, respectively. A variety of numerical examples are demonstrated to confirm the effectiveness and performance of the proposed schemes. Also we compare the accuracy and performance of our schemes with another method

Solution of 3-dimensional time-dependent viscous flows. Part 2: Development of the computer code

There is considerable interest in developing a numerical scheme for solving the time dependent viscous compressible three dimensional flow equations to aid in the design of helicopter rotors. The development of a computer code to solve a three dimensional unsteady approximate form of the Navier-Stokes equations employing a linearized block emplicit technique in conjunction with a QR operator scheme is described. Results of calculations of several Cartesian test cases are presented. The computer code can be applied to more complex flow fields such as these encountered on rotating airfoils