Solution of Singularly Perturbed Differential Difference Equations Using Higher Order Finite Differences

This paper is concerned with the numerical solution of the singularly perturbed differential-difference equations with small shifts called delay and advanced parameters. A fourth order finite difference method with a fitting factor is proposed for the solution of the singularly perturbed differential-difference equations with mixed shifts. The delay and advanced shifts are managed by Taylor series and an asymptotically equivalent singularly perturbed two-point boundary value problem is obtained. A fitting factor is introduced in the fourth order finite difference scheme for the problem which takes care of the small values of the perturbation parameter. This fitting factor is obtained from the asymptotic solution of singular perturbations. Thomas algorithm is used to solve the discrete system of the difference scheme. Convergence of the proposed method is analyzed. Maximum absolute errors in comparison with the several numerical experiments aretabulated to illustrate the proposed method.

The mode perturbation method (MPM) is suitable and efficient for solving the eigenvalue problem of a nonuniform soil deposit whose property varies with depth. However, results of the MPM do not always converge to the exact solution, when the variation of soil deposit property is discontinuous. This discontinuity is typical because soil is usually made up of sedimentary layers of different geologic materials. Based on the energy integral of the variational principle, a new mode perturbation method, the energy-based mode perturbation method (EMPM), is proposed to address the convergence of the perturbation solution on the natural frequencies and the corresponding mode shapes and is able to find solution whether the soil properties are continuous or not. First, the variational principle is used to transform the variable coefficient differential equation into an equivalent energy integral equation. Then, the natural mode shapes of the uniform shear beam with same height and boundary conditions are used as Ritz function. The EMPM transforms the energy integral equation into a set of nonlinear algebraic equations which significantly simplifies the eigenvalue solution of the soil layer with variable properties. Finally, the accuracy and convergence of this new method are illustrated with two case study examples. Numerical results show that the EMPM is more accurate and convergent than the MPM. As for the mode shapes of the uniform shear beam included in the EMPM, the additional 8 modes of vibration are sufficient in engineering applications.

In recent years, there has been a tendency towards renewable energy sources considering the damages caused by non-renewable energy resources to nature and humans. One of the renewable energy sources is wind and energy is obtained with the help of wind turbines. To determine the behavior of wind turbines under earthquake loads, dynamic characteristics are required. In this study, the differential transformation method is proposed to determine the free vibration analysis of wind turbines with a variable cross-section. The wind turbine is modeled as an equivalent variable continuous flexural beam and blade weight is considered as a point mass at the top of the structures. The differential equation representing the free vibration of the wind turbine is transformed into an algebraic equation with the help of differential transformation method and the angular frequencies and the mode shapes of the wind turbine are obtained by the help of the differential transformation method. In the study, a sample taken from the literature was solved with the presented method and the suitability of the method was investigated. The same wind turbine example also modeled by finite element modelling software, ABAQUS. Results of the finite element model and differential transformation method are compared with each other and the results are in good agreement.

In this paper, we use Chebyshev polynomials to seek the numerical solution of a class of multi-variable order fractional differential equation (MVODEs) that the fractional derivative is described in the Caputo-Prabhakar sense. Using operational matrices, the original equations are transferred to a system of algebraic equations. By solving the system of equations, the numerical solutions are acquired that this system may be solved numerically using an iterative algorithm. The effectiveness and convergence analysis of the numerical scheme is illustrated through four numerical examples.

In this paper, natural convection heat transfer over a vertical plate in a Darcy porous medium saturated with a nanofluid subject to heat generation/absorption was theoretically studied. The governing partial differential equations were transformed to a set of ordinary differential equations using similarity transformations and solved using finite difference method. The influence of parametric variation of the Brownian motion parameter, thermophoresis parameter and heat generation/absorption parameter on velocity, temperature and nanoparticles concentration profiles was graphically shown. Impact of non-dimensional parameters on the reduced Nusselt number and reduced Sherwood number was also investigated. The results showed that an increase in the heat generation/absorption parameter would increase temperature and velocity profiles; but, it would decrease concentration profiles. Increase of thermophoresis parameter increased magnitude of concentration profiles while not showing any significant effect on velocity and temperature profiles. The results also indicated that increase of Brownian motion parameter did not demonstrate any significant effect on the magnitude of velocity and temperature profiles. It was found that an increase in the heat generation/absorption parameter decreased the reduced Nusselt number whereas it increased the reduced Sherwood number. For negative values of the Brownian motion parameter, increase of the thermophoresis parameter increased the reduced Nusselt and Sherwood numbers.

This paper presents a numerical investigation for laminar forced convection flow of a radiating gas in a rectangular duct with a solid element that makes a backward facing step. The fluid is treated as a gray, absorbing, emitting and scattering medium. The governing differential equations consisting the continuity, momentum and energy are solved numerically by the computational fluid dynamics techniques. Since the present problem is a conjugate one and both gas and solid elements are considered in the computational domain, simultaneously, the numerical solution of Laplace equation is obtained in the solid element for temperature calculation in this area. Discretized forms of these equations are obtained using the finite volume method and solved by the SIMPLE algorithm. The radiative transfer equation (RTE) is also solved numerically by the discrete ordinate method (DOM) for computation of the radiative term in the gas energy equation.The streamline and isotherm plots in the gas flow and the distributions of convective, radiative and total Nusselt numbers along the solid-gas interface are presented. Besides, the effects of radiation conduction parameter and also solid to gas conduction ratio as two important parameters on thermo hydrodynamic characteristics of such thermal system are explored. It is revealed that the radiative Nusselt number on the interface surface is much affected by RC parameter but the radiation conduction parameter has not considerable effect on the convective Nusselt number. Comparison between the present numerical results with those obtained by other investigators for the case of non-conjugate problems shows good consistency.

In the present paper we consider initial value problems for Hilfer fractional differential equations and for system of Hilfer fractional differential equations. By using equivalent integral equations and some fixed point theorems, we study the local existence of solutions. We extend these local existence results globally with the help of continuation theorems and generalized Gronwall inequality.

High-order adaptive mesh refinement offers the potential for accurate and efficient resolution of problems in fluid dynamics and other fields where a wide range of length scales is present. A critical requirement for the interface closures used with these methods is stability in the context of hyperbolic systems of partial differential equations. In this study, a class of energy-stable high-order finite-difference interface closures is presented for grids with step resolution changes in one dimension. Asymptotic stability in time for these schemes is achieved by imposing a summation-by-parts condition on the interface closure, which is thus also nondissipative. Interface closures compatible with interior fourth- and sixth-order explicit, and fourth-order implicit centered schemes are presented. Validation tests include linear and nonlinear problems in one and in two dimensions with tensor-product grid refinement. A second class of stable high-order interface closures is presented for two-dimensional cell-centered grids with patch-refinement and step-changes in resolution. For these grids, coarse and fine nodes are not aligned at the mesh interfaces, resulting in hanging nodes. Stability is achieved by again imposing a summation-by-parts condition, resulting in nondissipative closures, at the cost of accuracy at corner interfaces. Interface stencils for an explicit fourth-order finite-difference scheme are presented for each geometry. Validation tests confirm the stability and accuracy of these closures for linear and nonlinear problems. The Richtmyer-Meshkov instability is investigated using a novel first-order perturbation of the two-dimensional Navier-Stokes equations about a shock-resolved base flow. The computational domain is efficiently resolved using the one-dimensional fourth-order interface scheme. Results are compared to analytic models of the instability, showing agreement with predicted asymptotic growth rates in the inviscid range, while significant discrepancies are noted in the transient growth phase. Viscous effects are found to be poorly predicted by existing models.

In this paper, we present a comparative study between the modified variational iteration method (MVIM) and a hybrid of Fourier transform and variational iteration method (FTVIM). The study outlines the efficiency and convergence of the two methods. The analysis is illustrated by investigating four singular partial differential equations with variable coefficients. The solution of singular partial differential equations usually needs a coordinate transformation in order to discard the singularity of the partial differential equation. Most often this transformation is not applicable and even does not exist. Therefore in this case the solution for the singular partial differential equation does not exist. In the present study the results of simulation for the singular partial differential equations with variable coefficients using the Fourier transform variational iteration method are compared with the results of simulation using the modified variational iteration method. The comparison shows that the effectiveness and accuracy of Fourier transform variational iteration method is more than that of the modified variational iteration method for the simulation of singular partial differential equations.

Commentary on: [L 45] Perturbations of Discontinuous Solutions of Non-Linear Systems of Differential Equations [L 48] An Ordinary Differential Equation with an Interval of Stability, a Separation Point, and an Interval of Instability [L 60] (with J. J. Levin), Singular Perturbations of Non-Linear Systems of Differential Equations and an Associated Boundary Layer Equation, [L 62] (with L. Flatto), Periodic Solutions

The objective of this article is investigation of dynamic response of thick multilayer functionally graded (FG) beam under generalized dynamic forces. The plane stress problem is exploited to describe the constitutive equation of thick FG beam to get realistic and accurate response. Applied dynamic forces are assumed to be sinusoidal harmonic, sinusoidal pulse or triangle in time domain and point load. Equations of motion of deep FG beam are derived based on the Hamilton principle from kinematic relations and constitutive equations of plane stress problem. The numerical finite element procedure is adopted to discretize the space domain of structure and transform partial differential equations of motion to ordinary differential equations in time domain. Numerical time integration method is used to solve the system of equations in time domain and find the time responses. Numerical parametric studies are performed to illustrate effects of force type, graduation parameter, geometrical and stacking sequence of layers on the time response of deep multilayer FG beams.

To solve fuzzy differential equations driven by Liu process, three Milstein schemes are proposed in this work, which are explicit Milstein scheme, semi-implicit Milstein scheme and improved Milstein scheme. Improved Milstein scheme is constructed by correcting the error with semi-implicit method, the error is the difference between the exact solution of fuzzy differential equations and the solution derived from Milstein scheme. These numerical methods are proved to have strong convergence with order two. Accompanying the results above, the concept of mean-stability of numerical schemes for fuzzy differential equations is put forward and analysed. For a linear test equation, it is showed the mean-stability region of improved Milstein scheme is bigger than explicit Milstein scheme. Finally, the accuracy and effectiveness of these schemes are confirmed through numerical examples.

Recently, the numerical modelling and simulation for anomalous subdiffusion equation (ASDE), which is a type of fractional partial differential equation( FPDE) and has been found with widely applications in modern engineering and sciences, are attracting more and more attentions. The current dominant numerical method for modelling ASDE is Finite Difference Method (FDM), which is based on a pre-defined grid leading to inherited issues or shortcomings. This paper aims to develop an implicit meshless approach based on the radial basis functions (RBF) for numerical simulation of the non-linear ASDE. The discrete system of equations is obtained by using the meshless shape functions and the strong-forms. The stability and convergence of this meshless approach are then discussed and theoretically proven. Several numerical examples with different problem domains are used to validate and investigate accuracy and efficiency of the newly developed meshless formulation. The results obtained by the meshless formulations are also compared with those obtained by FDM in terms of their accuracy and efficiency. It is concluded that the present meshless formulation is very effective for the modeling and simulation of the ASDE. Therefore, the meshless technique should have good potential in development of a robust simulation tool for problems in engineering and science which are governed by the various types of fractional differential equations.

MHD boundary layer flow in permeable medium along vertical plate with power-law deviation wall temperature in a gyratory system for the effects of Hall and ion-slip currents have investigated numerically. The ordinary differential equations are obtained from considered system of partial differential equations using usual transform. The Nachtsheim- Swigert shooting iteration method along with sixth order Runge-Kutta integration technique has employed to solve the non-dimensional system of ordinary differential equations. The velocities and temperature profiles are displayed graphically for the influence of physical parameters. The numerical values of the rate of change of velocities and temperature have been estimated for concerned physical parameters and represented in table.

Coupling Mortar Finite Element and Boundary Element Methods for 2D Navier-Stokes Equations Courant Element: Before and After Iterative Methods for Solving Stiff Elliptic Problems Straight and Curved Finite Elements of Class C1 and Some Applications to Thin Shell Problems Exact Controllability to Solve the Helmholtz Equation with Absorbing Boundary Conditions Cubic Version of FEM in Elliptic Problems with Interfaces and Singularities Least Squares Mixed Finite Elements Turbulence Modelling in Finite Element Industrial Applications Necessary and Sufficient Conditions for the Numerical Approximation of a Partial Differential Equation Depending on a Small Parameter Efficient Solution Methods for Compressible Flow Computations Parallel Finite Volume Algorithms for Solving the Time-Domain Maxwell Equations on Nonstructured Meshes Coupling Between Nonlinear Maxwell and Heat Equations for an Induction Heating Problem: Modelling and Numerical Methods Solving the 3D Harmonic Maxwell Equations with Finite Elements, Lagrange Multipliers, and Iterative Methods Some Applications of the Hierarchic High Order MITC Finite Elements for Reissner-Mindlin Plates Domain Decomposition for Immiscible Displacement in Single Porosity Systems An Error Estimator for Nonconforming Approximations of a Nonlinear Problem Some Observations on Raviart-Thomas Mixed Finite Elements in p Extension for Parabolic Problems Mixed Finite Element Methods in Fluid Structure Systems A Black-Box Solver for the Solution of General Nonlinear Functional Equations by Mixed FEM A Remark on the Asymptotic Behaviour of Parabolic Variational Inequalities and Their Finite Element Approximation by the Courant Element Domain Decomposition vs. Adaptivity Material Optimization of Composites. Part Contents.