Blood flow is modeled as non-Newtonian micropolar fluid. The non-linear governing equations of continuum and momentum in the cylindrical coordinate are being discretized using a finite difference approach and have been solved iteratively ,through Crank-Nicolson method. The blood velocity distribution, volumetric flow rate and Resistance to blood flow at the stenosis throat are computed for vari...

in this article we discussed the numerical solution of Burgers’ equation using multigrid method. We used implicit method for time discritization and Crank-Nicolson scheme for space discritization for fully discrete scheme. For improvement we used Multigrid method in fully discrete solution. And also Multigrid method accelerates convergence of a basics iterative method by global correction. Nume...

We develop systematically a numerical approximation strategy to discretize a hydrodynamic phase field model for a binary fluid mixture of two immiscible viscous fluids, derived using the generalized Onsager principle that warrants not only the variational structure but also the energy dissipation property. We first discretize the governing equations in space to arrive at a semi-discretized, tim...

Transparent boundary conditions (TBCs) for general Schrödinger– type equations on a bounded domain can be derived explicitly under the assumption that the given potential V is constant on the exterior of that domain. In 1D these boundary conditions are non–local in time (of memory type). Existing discretizations of these TBCs have accuracy problems and render the overall Crank–Nicolson finite d...

We propose and analyze a partitioned numerical method for the fully evolutionary Stokes-Darcy equations that model the coupling between surface and groundwater flows. The proposed method uncouples the surface from the groundwater flow by using the implicit-explicit combination of the Crank-Nicolson and Leapfrog methods for the discretization in time with added stabilization terms. We prove that...

The standard ‘parabolic’ approximation to the Helmholtz equation is used in order to model long-range propagation of sound in the sea in the presence of cylindrical symmetry in a domain with a rigid bottom of variable topography. The rigid bottom is modeled by a homogeneous Neumann condition and a paraxial approximation thereof proposed by Abrahamsson and Kreiss. The resulting initial-boundary-...

In this paper, we derive the exact rate of convergence of some approximation schemes associated to scalar stochastic differential equations driven by a fractional Brownian motion with Hurst index H . We consider two cases. If H > 1/2, the exact rate of convergence of the Euler scheme is determined. We show that the error of the Euler scheme converges almost surely to a random variable, which in...

An unconditionally stable finite difference scheme for systems whose dynamics are described by a fourth-order partial differential equation is developed with the use of a regular hexagonal grid. The scheme is motivated by the well-known Crank-Nicolson discretization that was originally developed for second-order systems and it is used to develop a discrete in time and space model of a deformabl...

We consider the fractional cable equation. For solution of fractional Cable equation involving Caputo fractional derivative, a new difference scheme is constructed based on Crank Nicholson difference scheme. We prove that the proposed method is unconditionally stable by using spectral stability technique.

In this paper convergence properties are discussed for some locally one-dimensional (LOD) splitting methods applied to linear parabolic initialboundary value problems. We shall consider unconditional convergence, where both the stepsize in time and the meshwidth in space tend to zero, independently of each other.