623 research outputs found
Computational and numerical analysis of differential equations using spectral based collocation method.
Doctoral Degree. University of KwaZulu-Natal, Pietermaritzburg.In this thesis, we develop accurate and computationally eļ¬cient spectral collocation-based methods,
both modiļ¬ed and new, and apply them to solve diļ¬erential equations. Spectral collocation-based
methods are the most commonly used methods for approximating smooth solutions of diļ¬erential
equations deļ¬ned over simple geometries. Procedurally, these methods entail transforming the gov
erning diļ¬erential equation(s) into a system of linear algebraic equations that can be solved directly.
Owing to the complexity of expanding the numerical algorithms to higher dimensions, as reported
in the literature, researchers often transform their models to reduce the number of variables or
narrow them down to problems with fewer dimensions. Such a process is accomplished by making
a series of assumptions that limit the scope of the study. To address this deļ¬ciency, the present
study explores the development of numerical algorithms for solving ordinary and partial diļ¬erential
equations deļ¬ned over simple geometries. The solutions of the diļ¬erential equations considered are
approximated using interpolating polynomials that satisfy the given diļ¬erential equation at se
lected distinct collocation points preferably the Chebyshev-Gauss-Lobatto points. The size of the
computational domain is particularly emphasized as it plays a key role in determining the number
of grid points that are used; a feature that dictates the accuracy and the computational expense of
the spectral method. To solve diļ¬erential equations deļ¬ned on large computational domains much
eļ¬ort is devoted to the development and application of new multidomain approaches, based on
decomposing large spatial domain(s) into a sequence of overlapping subintervals and a large time
interval into equal non-overlapping subintervals. The rigorous analysis of the numerical results con
ļ¬rms the superiority of these multiple domain techniques in terms of accuracy and computational
eļ¬ciency over the single domain approach when applied to problems deļ¬ned over large domains.
The structure of the thesis indicates a smooth sequence of constructing spectral collocation method
algorithms for problems across diļ¬erent dimensions. The process of switching between dimensions
is explained by presenting the work in chronological order from a simple one-dimensional problem
to more complex higher-dimensional problems. The preliminary chapter explores solutions of or
dinary diļ¬erential equations. Subsequent chapters then build on solutions to partial diļ¬erential
i
equations in order of increasing computational complexity. The transition between intermediate
dimensions is demonstrated and reinforced while highlighting the computational complexities in
volved. Discussions of the numerical methods terminate with development and application of a
new method namely; the trivariate spectral collocation method for solving two-dimensional initial
boundary value problems. Finally, the new error bound theorems on polynomial interpolation are
presented with rigorous proofs in each chapter to benchmark the adoption of the diļ¬erent numerical
algorithms. The numerical results of the study conļ¬rm that incorporating domain decomposition
techniques in spectral collocation methods work eļ¬ectively for all dimensions, as we report highly
accurate results obtained in a computationally eļ¬cient manner for problems deļ¬ned on large do
mains. The ļ¬ndings of this study thus lay a solid foundation to overcome major challenges that
numerical analysts might encounter
A numerical study on the viscous fingering instability of immiscible displacement in Hele-Shaw cells
In this thesis, the viscous fingering instability of radial immiscible displacement is analysed numerically using novel mesh-reduction and interface tracking techniques. Using a reduced Hele-Shaw model for the depth averaged lateral flow, viscous fingering instabilities are explored in flow regimes typical of subsurface carbon sequestration involving supercritical CO2 - brine displacements, i.e. with high capillary numbers, low mobility ratios and inhomogeneous permeability/temperature fields.
A high accuracy boundary element method (BEM) is implemented for the solution of homogeneous, finite mobility ratio immiscible displacements. Through efficient, explicit tracking of the sharp fluid-fluid interface, classical fingering processes such as spreading, shielding and splitting are analysed in the late stages of finger growth at low mobility ratios and high capillary numbers. Under these conditions, large differences are found compared with previous high or infinite mobility ratio models and critical events such as plume break-off and coalescence are analysed in much greater detail than has previously been attempted.
For the solution of inhomogeneous mobility problems, a novel meshless radial basis function-finite collocation method is developed that utilises a dynamic quadtree dataset and local enforcement of interface matching conditions. When coupled with the BEM, the numerical scheme allows the analysis of variable permeability effects and the transition in (de)stabilising mechanisms that occurs when the capillary number is increased with a fixed, spatially varying permeability. Finally, thermo-viscous fingering is explored in the context of immiscible flows, with a detailed mechanistic study presented to explain, for the first time, the immiscible thermo-viscous fingering process
Numerical Simulation
Nowadays mathematical modeling and numerical simulations play an important role in life and natural science. Numerous researchers are working in developing different methods and techniques to help understand the behavior of very complex systems, from the brain activity with real importance in medicine to the turbulent flows with important applications in physics and engineering. This book presents an overview of some models, methods, and numerical computations that are useful for the applied research scientists and mathematicians, fluid tech engineers, and postgraduate students
Semi-Infinite Structure Analysis with Bimodular Materials with Infinite Element
The modulus of elasticity of some materials changes under tensile and compressive states is simulated by constructing a typical material nonlinearity in a numerical analysis in this paper. The meshless Finite Block Method (FBM) has been developed to deal with 3D semi-infinite structures in the bimodular materials in this paper. The Lagrange polynomial interpolation is utilized to construct the meshless shape function with the mapping technique to transform the irregular finite domain or semi-infinite physical solids into a normalized domain. A shear modulus strategy is developed to present the nonlinear characteristics of bimodular material. In order to verify the efficiency and accuracy of FBM, the numerical results are compared with both analytical and numerical solutions provided by Finite Element Method (FEM) in four examples
A numerical study on the viscous fingering instability of immiscible displacement in Hele-Shaw cells
In this thesis, the viscous fingering instability of radial immiscible displacement is analysed numerically using novel mesh-reduction and interface tracking techniques. Using a reduced Hele-Shaw model for the depth averaged lateral flow, viscous fingering instabilities are explored in flow regimes typical of subsurface carbon sequestration involving supercritical CO2 - brine displacements, i.e. with high capillary numbers, low mobility ratios and inhomogeneous permeability/temperature fields.
A high accuracy boundary element method (BEM) is implemented for the solution of homogeneous, finite mobility ratio immiscible displacements. Through efficient, explicit tracking of the sharp fluid-fluid interface, classical fingering processes such as spreading, shielding and splitting are analysed in the late stages of finger growth at low mobility ratios and high capillary numbers. Under these conditions, large differences are found compared with previous high or infinite mobility ratio models and critical events such as plume break-off and coalescence are analysed in much greater detail than has previously been attempted.
For the solution of inhomogeneous mobility problems, a novel meshless radial basis function-finite collocation method is developed that utilises a dynamic quadtree dataset and local enforcement of interface matching conditions. When coupled with the BEM, the numerical scheme allows the analysis of variable permeability effects and the transition in (de)stabilising mechanisms that occurs when the capillary number is increased with a fixed, spatially varying permeability. Finally, thermo-viscous fingering is explored in the context of immiscible flows, with a detailed mechanistic study presented to explain, for the first time, the immiscible thermo-viscous fingering process
- ā¦