Loading...
Search
Toggle navigation
JCM
Home
About Journal
Information for Authors
Editorial Board
Subscription
Editorial Office
Table of Content
15 May 2023, Volume 41 Issue 3
Previous Issue
A FINITE VOLUME METHOD PRESERVING MAXIMUM PRINCIPLE FOR THE CONJUGATE HEAT TRANSFER PROBLEMS WITH GENERAL INTERFACE CONDITIONS
Huifang Zhou, Zhiqiang Sheng, Guangwei Yuan
2023, 41(3): 345369. DOI:
10.4208/jcm.2107m20200266
Asbtract
(
39
)
PDF
References

Related Articles

Metrics
In this paper, we present a unified finite volume method preserving discrete maximum principle (DMP) for the conjugate heat transfer problems with general interface conditions. We prove the existence of the numerical solution and the DMPpreserving property. Numerical experiments show that the nonlinear iteration numbers of the scheme in [24] increase rapidly when the interfacial coefficients decrease to zero. In contrast, the nonlinear iteration numbers of the unified scheme do not increase when the interfacial coefficients decrease to zero, which reveals that the unified scheme is more robust than the scheme in [24]. The accuracy and DMPpreserving property of the scheme are also verified in the numerical experiments.
EXPONENTIAL TIME DIFFERENCINGPADé FINITE ELEMENT METHOD FOR NONLINEAR CONVECTIONDIFFUSIONREACTION EQUATIONS WITH TIME CONSTANT DELAY
Haishen Dai, Qiumei Huang, Cheng Wang
2023, 41(3): 370394. DOI:
10.4208/jcm.2107m20210051
Asbtract
(
33
)
PDF
References

Related Articles

Metrics
In this paper, ETD3Padé and ETD4Padé Galerkin finite element methods are proposed and analyzed for nonlinear delayed convectiondiffusionreaction equations with Dirichlet boundary conditions. An ETDbased RK is used for time integration of the corresponding equation. To overcome a wellknown difficulty of numerical instability associated with the computation of the exponential operator, the Padé approach is used for such an exponential operator approximation, which in turn leads to the corresponding ETDPadé schemes. An unconditional
L
^{2}
numerical stability is proved for the proposed numerical schemes, under a global Lipshitz continuity assumption. In addition, optimal rate error estimates are provided, which gives the convergence order of
^{O}
(
k
^{3}
+
h
^{r}
) (ETD3 Padé) or
O
(
k
^{4}
+
h
^{r}
) (ETD4Padé) in the
L
^{2}
norm, respectively. Numerical experiments are presented to demonstrate the robustness of the proposed numerical schemes.
TWO NOVEL CLASSES OF ARBITRARY HIGHORDER STRUCTUREPRESERVING ALGORITHMS FOR CANONICAL HAMILTONIAN SYSTEMS
Yonghui Bo, Wenjun Cai, Yushun Wang
2023, 41(3): 395414. DOI:
10.4208/jcm.2108m20210076
Asbtract
(
15
)
PDF
References

Related Articles

Metrics
In this paper, we systematically construct two classes of structurepreserving schemes with arbitrary order of accuracy for canonical Hamiltonian systems. The one class is the symplectic scheme, which contains two new families of parameterized symplectic schemes that are derived by basing on the generating function method and the symmetric composition method, respectively. Each member in these schemes is symplectic for any fixed parameter. A more general form of generating functions is introduced, which generalizes the three classical generating functions that are widely used to construct symplectic algorithms. The other class is a novel family of energy and quadratic invariants preserving schemes, which is devised by adjusting the parameter in parameterized symplectic schemes to guarantee energy conservation at each time step. The existence of the solutions of these schemes is verified. Numerical experiments demonstrate the theoretical analysis and conservation of the proposed schemes.
STOCHASTIC VARIATIONAL INEQUALITY APPROACHES TO THE STOCHASTIC GENERALIZED NASH EQUILIBRIUM WITH SHARED CONSTRAINTS
Yanfang Zhang
2023, 41(3): 415436. DOI:
10.4208/jcm.2109m20200099
Asbtract
(
10
)
PDF
References

Related Articles

Metrics
In this paper, we consider the generalized Nash equilibrium with shared constraints in the stochastic environment, and we call it the stochastic generalized Nash equilibrium. The stochastic variational inequalities are employed to solve this kind of problems, and the expected residual minimization model and the conditional valueatrisk formulations defined by the residual function for the stochastic variational inequalities are discussed. We show the risk for different kinds of solutions for the stochastic generalized Nash equilibrium by the conditional valueatrisk formulations. The properties of the stochastic quadratic generalized Nash equilibrium are shown. The smoothing approximations for the expected residual minimization formulation and the conditional valueatrisk formulation are employed. Moreover, we establish the gradient consistency for the measurable smoothing functions and the integrable functions under some suitable conditions, and we also analyze the properties of the formulations. Numerical results for the applications arising from the electricity market model illustrate that the solutions for the stochastic generalized Nash equilibrium given by the ERM model have good properties, such as robustness, low risk and so on.
THE WASSERSTEINFISHERRAO METRIC FOR WAVEFORM BASED EARTHQUAKE LOCATION
Datong Zhou, Jing Chen, Hao Wu, Dinghui Yang, Lingyun Qiu
2023, 41(3): 437458. DOI:
10.4208/jcm.2109m20210045
Asbtract
(
10
)
PDF
References

Related Articles

Metrics
In this paper, we apply the WassersteinFisherRao (WFR) metric from the unbalanced optimal transport theory to the earthquake location problem. Compared with the quadratic Wasserstein (
W
_{2}
) metric from the classical optimal transport theory, the advantage of this method is that it retains the important amplitude information as a new constraint, which avoids the problem of the degeneration of the optimization objective function near the real earthquake hypocenter and origin time. As a result, the deviation of the global minimum of the optimization objective function based on the WFR metric from the true solution can be much smaller than the results based on the
W
_{2}
metric when there exists strong data noise. Thus, we develop an accurate earthquake location method under strong data noise. Many numerical experiments verify our conclusions.
ANISOTROPIC
EQ
_{1}
^{ROT}
FINITE ELEMENT APPROXIMATION FOR A MULTITERM TIMEFRACTIONAL MIXED SUBDIFFUSION AND DIFFUSIONWAVE EQUATION
Huijun Fan, Yanmin Zhao, Fenling Wang, Yanhua Shi, Fawang Liu
2023, 41(3): 459482. DOI:
10.4208/jcm.2110m20210180
Asbtract
(
12
)
PDF
References

Related Articles

Metrics
By employing
EQ
_{1}
^{rot}
nonconforming finite element, the numerical approximation is presented for multiterm timefractional mixed subdiffusion and diffusionwave equation on anisotropic meshes. Comparing with the multiterm timefractional subdiffusion equation or diffusionwave equation, the mixed case contains a special timespace coupled derivative, which leads to many difficulties in numerical analysis. Firstly, a fully discrete scheme is established by using nonconforming finite element method (FEM) in spatial direction and L1 approximation coupled with CrankNicolson (L1CN) scheme in temporal direction. Furthermore, the fully discrete scheme is proved to be unconditional stable. Besides, convergence and superclose results are derived by using the properties of
EQ
_{1}
^{rot}
nonconforming finite element. What's more, the global superconvergence is obtained via the interpolation postprocessing technique. Finally, several numerical results are provided to demonstrate the theoretical analysis on anisotropic meshes.
INVERSE CONDUCTIVITY PROBLEM WITH INTERNAL DATA
Faouzi Triki, Tao Yin
2023, 41(3): 483501. DOI:
10.4208/jcm.2111m20210093
Asbtract
(
9
)
PDF
References

Related Articles

Metrics
This paper concerns the reconstruction of a scalar coefficient of a secondorder elliptic equation in divergence form posed on a bounded domain from internal data. This problem finds applications in multiwave imaging, greedy methods to approximate parameterdependent elliptic problems, and image treatment with partial differential equations. We first show that the inverse problem for smooth coefficients can be rewritten as a linear transport equation. Assuming that the coefficient is known near the boundary, we study the wellposedness of associated transport equation as well as its numerical resolution using discontinuous Galerkin method. We propose a regularized transport equation that allow us to derive rigorous convergence rates of the numerical method in terms of the order of the polynomial approximation as well as the regularization parameter. We finally provide numerical examples for the inversion assuming a lower regularity of the coefficient, and using synthetic data.
THEORETICAL ANALYSIS OF THE REPRODUCING KERNEL GRADIENT SMOOTHING INTEGRATION TECHNIQUE IN GALERKIN MESHLESS METHODS
Xiaolin Li
2023, 41(3): 502525. DOI:
10.4208/jcm.2201m20210361
Asbtract
(
9
)
PDF
References

Related Articles

Metrics
Numerical integration poses greater challenges in Galerkin meshless methods than finite element methods owing to the nonpolynomial feature of meshless shape functions. The reproducing kernel gradient smoothing integration (RKGSI) is one of the optimal numerical integration techniques in Galerkin meshless methods with minimum integration points. In this paper, properties, quadrature rules and the effect of the RKGSI on meshless methods are analyzed. The existence, uniqueness and error estimates of the solution of Galerkin meshless methods under numerical integration with the RKGSI are established. A procedure on how to choose quadrature rules to recover the optimal convergence rate is presented.
ANNEALED IMPORTANCE SAMPLING FOR ISING MODELS WITH MIXED BOUNDARY CONDITIONS
Lexing Ying
2023, 41(3): 526534. DOI:
10.4208/jcm.2211m20220172
Asbtract
(
8
)
PDF
References

Related Articles

Metrics
This note introduces a method for sampling Ising models with mixed boundary conditions. As an application of annealed importance sampling and the SwendsenWang algorithm, the method adopts a sequence of intermediate distributions that keeps the temperature fixed but turns on the boundary condition gradually. The numerical results show that the variance of the sample weights is relatively small.
Current Issue
Earlier Issues
Advanced Search
Most Read Articles
Most Downloaded Articles
Visit Hong Kong Site