Numerical solution of singularly perturbed convection–diffusion problem using parameter uniform B-spline collocation method

AbstractThis paper is concerned with a numerical scheme to solve a singularly perturbed convection–diffusion problem. The solution of this problem exhibits the boundary layer on the right-hand side of the domain due to the presence of singular perturbation parameter ɛ. The scheme involves B-spline collocation method and appropriate piecewise-uniform Shishkin mesh. Bounds are established for the derivative of the analytical solution. Moreover, the present method is boundary layer resolving as well as second-order uniformly convergent in the maximum norm. A comprehensive analysis has been given to prove the uniform convergence with respect to singular perturbation parameter. Several numerical examples are also given to demonstrate the efficiency of B-spline collocation method and to validate the theoretical aspects

On the design and implementation of a hybrid numerical method for singularly perturbed two-point boundary value problems

>Magister Scientiae - MScWith the development of technology seen in the last few decades, numerous solvers have been developed to provide adequate solutions to the problems that model different aspects of science and engineering. Quite often, these solvers are tailor-made for specific classes of problems. Therefore, more of such must be developed to accompany the growing need for mathematical models that help in the understanding of the contemporary world. This thesis treats two point boundary value singularly perturbed problems. The solution to this type of problem undergoes steep changes in narrow regions (called boundary or internal layer regions) thus rendering the classical numerical procedures inappropriate. To this end, robust numerical methods such as finite difference methods, in particular fitted mesh and fitted operator methods have extensively been used. While the former consists of transforming the continuous problem into a discrete one on a non-uniform mesh, the latter involves a special discretisation of the problem on a uniform mesh and are known to be more accurate. Both classes of methods are suitably designed to accommodate the rapid change(s) in the solution. Quite often, finite difference methods on piece-wise uniform meshes (of Shishkin-type) are adopted. However, methods based on such non-uniform meshes, though layer-resolving, are not easily extendable to higher dimensions. This work aims at investigating the possibility of capitalising on the advantages of both fitted mesh and fitted operator methods. Theoretical results are confirmed by extensive numerical simulations

Higher order numerical methods for singular perturbation problems

Philosophiae Doctor - PhDIn recent years, there has been a great interest towards the higher order numerical methods for singularly perturbed problems. As compared to their lower order counterparts, they provide better accuracy with fewer mesh points. Construction and/or implementation of direct higher order methods is usually very complicated. Thus a natural choice is to use some convergence acceleration techniques, e.g., Richardson extrapolation, defect correction, etc. In this thesis, we will consider various classes of problems described by singularly perturbed ordinary and partial differential equations. For these problems, we design some novel numerical methods and attempt to increase their accuracy as well as the order of convergence. We also do the same for existing numerical methods in some instances. We find that, even though the Richardson extrapolation technique always improves the accuracy, it does not perform equally well when applied to different methods for certain classes of problems. Moreover, while in some cases it improves the order of convergence, in other cases it does not. These issues are discussed in this thesis for linear and nonlinear singularly perturbed ODEs as well as PDEs. Extrapolation techniques are analyzed thoroughly in all the cases, whereas the limitations of the defect correction approach for certain problems is indicated at the end of the thesis.South Afric

Robust computational methods for two-parameter singular perturbation problems

Magister Scientiae - MScThis thesis is concerned with singularly perturbed two-parameter problems. We study a tted nite difference method as applied on two different meshes namely a piecewise mesh (of Shishkin type) and a graded mesh (of Bakhvalov type) as well as a tted operator nite di erence method. We notice that results on Bakhvalov mesh are better than those on Shishkin mesh. However, piecewise uniform meshes provide a simpler platform for analysis and computations. Fitted operator methods are even simpler in these regards due to the ease of operating on uniform meshes. Richardson extrapolation is applied on one of the tted mesh nite di erence method (those based on Shishkin mesh) as well as on the tted operator nite di erence method in order to improve the accuracy and/or the order of convergence. This is our main contribution to this eld and in fact we have achieved very good results after extrapolation on the tted operator finitete difference method. Extensive numerical computations are carried out on to confirm the theoretical results.South Afric

Technical Evaluation Report for Symposium AVT-147: Computational Uncertainty in Military Vehicle Design

The complexity of modern military systems, as well as the cost and difficulty associated with experimentally verifying system and subsystem design makes the use of high-fidelity based simulation a future alternative for design and development. The predictive ability of such simulations such as computational fluid dynamics (CFD) and computational structural mechanics (CSM) have matured significantly. However, for numerical simulations to be used with confidence in design and development, quantitative measures of uncertainty must be available. The AVT 147 Symposium has been established to compile state-of-the art methods of assessing computational uncertainty, to identify future research and development needs associated with these methods, and to present examples of how these needs are being addressed and how the methods are being applied. Papers were solicited that address uncertainty estimation associated with high fidelity, physics-based simulations. The solicitation included papers that identify sources of error and uncertainty in numerical simulation from either the industry perspective or from the disciplinary or cross-disciplinary research perspective. Examples of the industry perspective were to include how computational uncertainty methods are used to reduce system risk in various stages of design or development

Modelling two-phase flow and transport effects of multi-component fuels

Three novel multicomponent fuel spray droplet evaporation models are developed by employing the theory of continuous thermodynamics(CT) with the aim of applying them in the design and analysis of various energy conversion devices such as, aircraft jet engines, liquid-fuel rocket engines, diesel engines, and industrial furnaces. The CT methodology seeks to represent complex mixtures - for example,aviation kerosene or JP8 that typically comprise blends of a large number of chemical compounds by using probability distribution functions (PDFs). The components of JP8, which is constituted by the homologous series of paraffin, naphthene, and aromatic hydrocarbons; are each represented by the Pearson-Shultz type three-parameter gamma PDF, where the three (shape, scale, and origin) parameters characterise changes in the mixture composition. The phase transition of the liquid droplet due to evaporation is modelled using both low-pressure (LP) and high-pressure (HP) vapour-liquid equilibrium (VLE) models employing various mixing and combining rules by applying a general cubic equation of state (CEOS). Interestingly enough, the phase transition of the liquid fuel into vapour mixture is characterised by a change in the PDF scale parameter alone. Once the description of the fuel mixture is complete, the traditional species and energy transport equations both for the liquid and vapour phases respectively, are re-written using the composition PDF moments under Lagrangian and Eulerian frameworks. In order to solve the governing equations for the three droplet evaporation models, which characteristically involve phase change and a moving interface, a novel fully Adaptive Method Of Lines using B-Spline Collocation (AMOLBSC) is developed. The models are tested at various pressures, temperatures and convective conditions, including at a lean, premixed, prevaporised (LPP) combustor operating condition. In general, the computational results at an ambient pressure close to atmospheric showed good to excellent agreement against available experimental data in the literature. However, for ambient conditions with elevated-high pressures and temperatures only models that employ the HP formulation gave reliable results. In particular, when the liquid is at or near its critical pressure and temperature it is characterised by faster vaporisation and shorter droplet lifetime, including some evidence of liquid mass diffusion. The liquid model that incorporates the effects of liquid core circulation using semiempirical approximation and adaptive mesh refinement (AMR) technique is the most accurate and computationally efficient, although further work is required to establish its ranges of applicability.EThOS - Electronic Theses Online ServiceGBUnited Kingdo

Research in progress and other activities of the Institute for Computer Applications in Science and Engineering

This report summarizes research conducted at the Institute for Computer Applications in Science and Engineering in applied mathematics and computer science during the period April 1, 1993 through September 30, 1993. The major categories of the current ICASE research program are: (1) applied and numerical mathematics, including numerical analysis and algorithm development; (2) theoretical and computational research in fluid mechanics in selected areas of interest to LaRC, including acoustic and combustion; (3) experimental research in transition and turbulence and aerodynamics involving LaRC facilities and scientists; and (4) computer science

Numerical Computation, Data Analysis and Software in Mathematics and Engineering

The present book contains 14 articles that were accepted for publication in the Special Issue “Numerical Computation, Data Analysis and Software in Mathematics and Engineering” of the MDPI journal Mathematics. The topics of these articles include the aspects of the meshless method, numerical simulation, mathematical models, deep learning and data analysis. Meshless methods, such as the improved element-free Galerkin method, the dimension-splitting, interpolating, moving, least-squares method, the dimension-splitting, generalized, interpolating, element-free Galerkin method and the improved interpolating, complex variable, element-free Galerkin method, are presented. Some complicated problems, such as tge cold roll-forming process, ceramsite compound insulation block, crack propagation and heavy-haul railway tunnel with defects, are numerically analyzed. Mathematical models, such as the lattice hydrodynamic model, extended car-following model and smart helmet-based PLS-BPNN error compensation model, are proposed. The use of the deep learning approach to predict the mechanical properties of single-network hydrogel is presented, and data analysis for land leasing is discussed. This book will be interesting and useful for those working in the meshless method, numerical simulation, mathematical model, deep learning and data analysis fields