We present a new algorithm to simulate unsteady viscoelastic flows in abrupt contraction channels. In our approach we split the viscoelastic terms of the Oldroyd-B constitutive equation using Duhamel s formula and discretize the resulting PDEs using a semi-implicit finite difference method based on a Lax–Wendroff method for hyperbolic terms. In particular, we leave a small residual elastic term...

High order fast sweeping methods have been developed recently in the literature to solve static Hamilton-Jacobi equations efficiently. Comparing with the first order fast sweeping methods, the high order fast sweeping methods are more accurate, but they often require additional numerical boundary treatment for several grid points near the boundary because of the wider numerical stencil. It is p...

In this paper, a new scheme of arbitrary high order accuracy in both space and time is proposed to solve hyperbolic conservative laws. Based on the idea of flux vector splitting(FVS) scheme, we split all the space and time derivatives in the Taylor expansion of the numerical flux into two parts: one part with positive eigenvalues, another part with negative eigenvalues. According to a Lax-Wendr...

To gain insight into cardio-arterial interactions, a coupled left ventricle-systemic artery (LV-SA) model is developed that incorporates a three-dimensional finite-strain left ventricle (LV), and a physiologically-based one-dimensional model for the systemic arteries (SA). The coupling of the LV model and the SA model is achieved by matching the pressure and the flow rate at the aortic root, i....

The goal of this article is to design robust and simple first order explicit solvers for one-dimensional nonconservative hyperbolic systems. These solvers are intended to be used as the basis for higher order methods for one or multidimensional problems. The starting point for the development of these solvers is the general definition of a Roe linearization introduced by Toumi in 1992 based on ...

We develop a theoretical tool to examine the properties of numerical schemes for advection equations. To magnify the defects of a scheme we consider a convection-reaction equation 0021-9 doi:10. q Th * Co E-m ut þ ðjujq=qÞx 1⁄4 u; u; x 2 R; t 2 Rþ; q > 1: It is shown that, if a numerical scheme for the advection part is performed with a splitting method, the intrinsic properties of the scheme a...

It is shown that for quasi-linear hyperbolic systems of the conservation form Wt = —Fx = —AWX, it is possible to build up relatively simple finite-difference numerical schemes accurate to 3rd and 4th order provided that the matrix A satisfies commutativity relations with its partial-derivative-matrices. These schemes generalize the Lax-Wendroff 2nd order scheme, and are written down explicitly....

A 3-D model for atmospheric pollutant transport is proposed considering a set of coupled convection–diffusion–reaction equations. The convective phenomenon is mainly produced by a wind field obtained from a 3-D mass consistent model. In particular, the modelling of oxidation and hydrolysis of sulphur and nitrogen oxides released to the surface layer is carried out by using a linear module of ch...

ABSTRACT The conservative high-order accurate spectral difference method is presented for simulation of rotating shallowwater equations. The method is formulated using Lagrange interpolations on Gauss-Lobatto points for the desired order of accuracy without suffering numerical dissipation and dispersion errors. The optimal third-order total variation diminishing (TVD) Runge-Kutta algorithm is u...