Loading...
Search
Toggle navigation
JCM
Home
About Journal
Information for Authors
Editorial Board
Subscription
Editorial Office
Table of Content
15 July 2018, Volume 36 Issue 4
Previous Issue
Next Issue
A WEAK GALERKIN FINITE ELEMENT METHOD FOR THE LINEAR ELASTICITY PROBLEM IN MIXED FORM
Ruishu Wang, Ran Zhang
2018, 36(4): 469-491. DOI:
10.4208/jcm.1701-m2016-0733
Asbtract
(
175
)
PDF
References
|
Related Articles
|
Metrics
In this paper, we use the weak Galerkin (WG) finite element method to solve the mixed form linear elasticity problem. In the mixed form, we get the discrete of proximation of the stress tensor and the displacement field. For the WG methods, we define the weak function and the weak differential operator in an optimal polynomial approximation spaces. The optimal error estimates are given and numerical results are presented to demonstrate the efficiency and the accuracy of the weak Galerkin finite element method.
BLOCK-CENTERED FINITE DIFFERENCE METHODS FOR NON-FICKIAN FLOW IN POROUS MEDIA
Xiaoli Li, Hongxing Rui
2018, 36(4): 492-516. DOI:
10.4208/jcm.1701-m2016-0628
Asbtract
(
169
)
PDF
References
|
Related Articles
|
Metrics
In this article, two block-centered finite difference schemes are introduced and analyzed to solve the parabolic integro-differential equation arising in modeling non-Fickian flow in porous media. One scheme is Euler backward scheme with first order accuracy in time increment while the other is Crank-Nicolson scheme with second order accuracy in time increment. Stability analysis and second-order error estimates in spatial meshsize for both pressure and velocity in discrete
L
2
norms are established on non-uniform rectangular grid. Numerical experiments using the schemes show that the convergence rates are in agreement with the theoretical analysis.
HIGH ORDER STABLE MULTI-DOMAIN HYBRID RKDG AND WENO-FD METHODS
Fan Zhang, Tiegang Liu, Jian Cheng
2018, 36(4): 517-541. DOI:
10.4208/jcm.1702-m2016-0707
Asbtract
(
105
)
PDF
References
|
Related Articles
|
Metrics
Recently, a kind of high order hybrid methods based on Runge-Kutta discontinuous Galerkin (RKDG) method and weighted essentially non-oscillatory finite difference (WENO-FD) scheme was proposed. Those methods are computationally efficient, however stable problems might sometimes be encountered in practical applications. In this work, we first analyze the linear stabilities of those methods based on the Heuristic theory. We find that the conservative hybrid method is linearly unstable if the numerical flux at the coupling interface is chosen to be ‘downstream’. Then we introduce two ways of healing this defect. One is to choose the numerical flux at the coupling interface to be ‘upstream’. The other is to employ a slope limiter function to enforce the hybrid method satisfying the local total variation diminishing (TVD) condition. In the end, numerical experiments are provided to validate the effectiveness of the proposed methods.
QUASI-NEWTON WAVEFORM RELAXATION BASED ON ENERGY METHOD
Yaolin Jiang, Zhen Miao
2018, 36(4): 542-562. DOI:
10.4208/jcm.1702-m2016-0700
Asbtract
(
121
)
PDF
References
|
Related Articles
|
Metrics
A quasi-Newton waveform relaxation (WR) algorithm for semi-linear reaction-diffusion equations is presented at first in this paper. Using the idea of energy estimate, a general proof method for convergence of the continuous case and the discrete case of quasi-Newton WR is given, which appears to be the superlinear rate. The semi-linear wave equation and semi-linear coupled equations can similarly be solved by quasi-Newton WR algorithm and be proved as convergent with the energy inequalities. Finally several parallel numerical experiments are implemented to confirm the effectiveness of the above theories.
ANOMALOUS DIFFUSION IN FINITE LENGTH FINGERS COMB FRAME WITH THE EFFECTS OF TIME AND SPACE RIESZ FRACTIONAL CATTANEO-CHRISTOV FLUX AND POISEUILLE FLOW
Lin Liu, Liancun Zheng, Fawang Liu, Xinxin Zhang
2018, 36(4): 563-578. DOI:
10.4208/jcm.1702-m2016-0627
Asbtract
(
89
)
PDF
References
|
Related Articles
|
Metrics
This paper presents an investigation on the anomalous diffusion in finite length fingers comb frame, the time and space Riesz fractional Cattaneo-Christov flux is introduced with the Oldroyds' upper convective derivative and the effect of Poiseuille flow is also taken into account. Formulated governing equation possesses the coexisting characteristics of parabolicity and hyperbolicity. Numerical solution is obtained by the
L
1
-scheme and shifted Grünwald formulae, which is verified by introducing a source item to construct an exact solution. The effects, such as time and space fractional parameters, relaxation parameter and the ratio of the pressure gradient and viscosity coefficient, on the spatial and temporal evolution of particles distribution and dynamic characteristics are shown graphically and analyzed in detail.
A MODIFIED PRECONDITIONER FOR PARAMETERIZED INEXACT UZAWA METHOD FOR INDEFINITE SADDLE POINT PROBLEMS
Xinhui Shao, Chen Li, Tie Zhang, Changjun Li
2018, 36(4): 579-590. DOI:
10.4208/jcm.1702-m2016-0665
Asbtract
(
95
)
PDF
References
|
Related Articles
|
Metrics
The preconditioner for parameterized inexact Uzawa methods have been used to solve some indefinite saddle point problems. Firstly, we modify the preconditioner by making it more generalized, then we use theoretical analyses to show that the iteration method converges under certain conditions. Moreover, we discuss the optimal parameter and matrices based on these conditions. Finally, we propose two improved methods. Numerical experiments are provided to show the effectiveness of the modified preconditioner. All methods have fantastic convergence rates by choosing the optimal parameter and matrices.
HIGH ORDER COMPACT MULTISYMPLECTIC SCHEME FOR COUPLED NONLINEAR SCHRÖDINGER-KDV EQUATIONS
Lan Wang, Yushun Wang
2018, 36(4): 591-604. DOI:
10.4208/jcm.1702-m2016-0789
Asbtract
(
106
)
PDF
References
|
Related Articles
|
Metrics
In this paper, a novel multisymplectic scheme is proposed for the coupled nonlinear Schrödinger-KdV (CNLS-KdV) equations. The CNLS-KdV equations are rewritten into the multisymplectic Hamiltonian form by introducing some canonical momenta. To simulate the problem efficiently, the CNLS-KdV equations are approximated by a high order compact method in space which preserves
N
semi-discrete multisymplectic conservation laws. We then discretize the semi-discrete system by using a symplectic midpoint scheme in time. Thus, a full-discrete multisymplectic scheme is obtained for the CNLS-KdV equations. The conservation laws of the full-discrete scheme are analyzed. Some numerical experiments are presented to further verify the convergence and conservation laws of the new scheme.
A NEW BOUNDARY CONDITION FOR RATE-TYPE NON-NEWTONIAN DIFFUSIVE MODELS AND THE STABLE MAC SCHEME
Kun Li, Youngju Lee, Christina Starkey
2018, 36(4): 605-626. DOI:
10.4208/jcm.1703-m2015-0359
Asbtract
(
126
)
PDF
References
|
Related Articles
|
Metrics
We present a new Dirichlet boundary condition for the rate-type non-Newtonian diffusive constitutive models. The newly proposed boundary condition is compared with two such well-known and popularly used boundary conditions as the pure Neumann condition[1] and the Dirichlet condition by Sureshkumar and Beris[2]. Our condition is demonstrated to be more stable and robust in a number of numerical test cases. A new Dirichlet boundary condition is implemented in the framework of the finite difference Marker and Cell (MAC) method. In this paper, we also present an energy-stable finite difference MAC scheme that preserves the positivity for the conformation tensor and show how the addition of the diffusion helps the energy-stability in a finite difference MAC scheme-setting.
Current Issue
Earlier Issues
Advanced Search
Most Read Articles
Most Downloaded Articles
Visit Hong Kong Site