New approach to the Lax‐Wendroff modified differential equation for linear and nonlinear advection

The paper presents an enhanced analysis of the Lax‐Wendroff difference scheme—up to the eighth‐order with respect to time and space derivatives—of the modified‐partial differential equation (MDE) of the constant‐wind‐speed advection equation. The modified equation has been so far derived mainly as a fourth‐order equation. The Π ‐form of the first differential approximation (differential approximation or equivalent equation) derived by expressing the time derivatives in terms of the space derivatives is used for presenting the MDE. The obtained coefficients at higher order derivatives are analyzed for indications of the character of the dissipative and dispersive errors. The authors included a part of the stencil applied for determining the modified differential equation up to the eighth‐order of the analyzed modified differential equation for the second‐order Lax‐Wendroff scheme. Neither the derived coefficients at the space derivatives of order p ∈ (7 – 8) in the modified differential equation for the Lax‐Wendroff difference scheme nor the results of analyses on the basis of these coefficients of the group velocity, phase shift errors, or dispersive and dissipative features of the scheme have been published. The MDEs for 2 two‐step variants of the Lax‐Wendroff type difference schemes and the MacCormack predictor–corrector scheme (see MacCormack's study) constructed for the scalar hyperbolic conservation laws are also presented in this paper. The analysis of the inviscid Burgers equation solution with the initial condition in a form of a shock wave has been discussed on their basis. The inviscid Burgers equation with the source is also presented. The theory of MDE started to develop after the paper of C. W. Hirt was published in 1968.     

A Steffensen's type method in Banach spaces with applications on boundary-value problems

In this paper, a modified Steffensen's type iterative scheme for the numerical solution of a system of nonlinear equations is studied. Two convergence theorems are presented. The numerical solution of boundary-value problems by the multiple shooting method using the proposed iterative scheme is analyzed.     

Existence of solution of a finite volume scheme preserving maximum principle for diffusion equations

In this article, a cell‐centered finite volume scheme preserving maximum principle for diffusion equations with scalar coefficients is developed. The construction of the scheme consists of three steps: at first the discrete normal flux is obtained by a linear combination of two single‐sided fluxes, then the tangential term of the normal flux is modified by using a nonlinear combination of two single‐sided tangential fluxes, finally the auxiliary unknowns in the tangential fluxes are calculated by the convex combinations of the cell‐centered unknowns. It is proved that this nonlinear scheme satisfies the discrete maximum principle (DMP). Moreover, the existence of a solution of the nonlinear scheme is proved by using the Brouwer's fixed point theorem and the bounded estimates. Numerical experiments are presented to show that the scheme not only satisfies DMP, but also obtains the second‐order accuracy and conservation.© 2017 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 34: 80–96, 2018     

A high‐order finite difference method for 1D nonhomogeneous heat equations

In this article a sixth‐order approximation method (in both temporal and spatial variables) for solving nonhomogeneous heat equations is proposed. We first develop a sixth‐order finite difference approximation scheme for a two‐point boundary value problem, and then heat equation is approximated by a system of ODEs defined on spatial grid points. The ODE system is discretized to a Sylvester matrix equation via boundary value method. The obtained algebraic system is solved by a modified Bartels‐Stewart method. The proposed approach is unconditionally stable. Numerical results are provided to illustrate the accuracy and efficiency of our approximation method along with comparisons with those generated by the standard second‐order Crank‐Nicolson scheme as well as Sun‐Zhang's recent fourth‐order method. © 2008 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2009     

The Cayley Transform as a Time Discretization Scheme

We interpret the Cayley transform of linear (finite- or infinite-dimensional) state space systems as a numerical integration scheme of Crank–Nicolson type. The scheme is known as Tustin's method in the engineering literature, and it has the following important Hamiltonian integrator property: if Tustin's method is applied to a conservative (continuous time) linear system, then the resulting (discrete time) linear system is conservative in the discrete time sense. The purpose of this paper is to study the convergence of this integration scheme from the input/output point of view.     

Newton–Milstein scheme for stochastic differential equations and its fast uniform convergence

We show that a modified Milstein scheme combined with explicit Newton’s method enables us to construct fast converging sequences of approximate solutions of stochastic differential equations. The fast uniform convergence of our Newton–Milstein scheme follows from Amano’s probabilistic second-order error estimate, which had been an open problem since 1991. The Newton–Milstein scheme, which is based on a modified Milstein scheme and the symbolic Newton’s method, will be classified as a numerical and computer algebraic hybrid method and it may give a new possibility to the study of computer algebraic method in stochastic analysis.     

Nonasymptotic analysis of robust regression with modified Huber's loss

To achieve robustness against the outliers or heavy-tailed sampling distribution, we consider an Ivanov regularized empirical risk minimization scheme associated with a modified Huber's loss for nonparametric regression in reproducing kernel Hilbert space. By tuning the scaling and regularization parameters in accordance with the sample size, we develop nonasymptotic concentration results for such an adaptive estimator. Specifically, we establish the best convergence rates for prediction error when the conditional distribution satisfies a weak moment condition.     

An explicit fourth-order staggered finite-difference time-domain method for Maxwell's equations

A new explicit fourth-order accurate staggered finite-difference time-domain (FDTD) scheme is proposed and applied to electromagnetic wave problems. It is fourth-order accurate in both space and time, conditionally stable, and highly efficient (with respect to Yee's scheme) and still retains much of the original simplicity of Yee's scheme. Both extension to perfectly matched layers and modification to deal with dielectric interfaces and perfectly conducting boundaries of the scheme have also been presented. Numerical examples are shown to illustrate the efficiency of the method.     

A modified method of characteristics incorporating streamline diffusion

A hybrid finite-element method, combining ideas from a modified method of characteristics and the streamline diffusion method, delivers accurate solutions to the advection–diffusion equation. An error analysis for the case of tensorial diffusion shows that the lowest-order version of the scheme, which allows one to use a symmetric linear solvers at each time step, possesses first-order accuracy in time and space. Numerical experiments demonstrate the scheme's ability to model advection-dominated transport of solute plumes without distorting sharp fronts. © 1995 John Wiley & Sons, Inc.     

A Moreau-type Variational Integrator

In this paper we derive a variational integrator for nonsmooth mechanical systems by discretizing the principle of virtual action with finite elements in time. After the discretization with local finite elements, the constitutive laws for the contact forces are introduced as in Moreau's time stepping scheme. This derivation shows exemplary how variational integrators for systems with frictional unilateral constraints can be derived. The long-time energy behavior of the presented scheme is compared with the behavior of Moreau's stepping scheme on an example system. (© 2016 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Application of MFCAV Riemann Solver to Maire's Cell-Centered Lagrangian Method

In this paper, we apply arbitrary Riemann solvers, which may not satisfy the Maire's requirement, to the Maire's node-based Lagrangian scheme developed in [P. H. Maire et al., SIAM J. Sci. Comput, 29 (2007), 1781-1824]. In particular, we apply the so-called Multi-Fluid Channel on Averaged Volume (MFCAV) Riemann solver and a Riemann solver that adaptively combines the MFCAV solver with other more dissipative Riemann solvers to the Maire's scheme. It is noted that neither of the two solvers satisfies the Maire's requirement. Numerical experiments are presented to demonstrate that the application of the two Riemann solvers is successful.     

Optimal error estimates and modified energy conservation identities of the ADI-FDTD scheme on staggered grids for 3D Maxwell’s equations

This paper is concerned with the optimal error estimates and energy conservation properties of the alternating direction implicit finite-difference time-domain(ADI-FDTD) method which is a popular scheme for solving the 3D Maxwell’s equations.Precisely,for the case with a perfectly electric conducting(PEC) boundary condition we establish the optimal second-order error estimates in both space and time in the discrete H 1-norm for the ADI-FDTD scheme,and prove the approximate divergence preserving property that if the divergence of the initial electric and magnetic fields are zero,then the discrete L 2-norm of the discrete divergence of the ADI-FDTD solution is approximately zero with the second-order accuracy in both space and time.The key ingredient is two new discrete modified energy norms which are second-order in time perturbations of two new energy conservation laws for the Maxwell’s equations introduced in this paper.Furthermore,we prove that,in addition to two known discrete modified energy identities which are second-order in time perturbations of two known energy conservation laws,the ADI-FDTD scheme also satisfies two new discrete modified energy identities which are second-order in time perturbations of the two new energy conservation laws.This means that the ADI-FDTD scheme is unconditionally stable under the four discrete modified energy norms.Experimental results which confirm the theoretical results are presented.     

Analysis of the energy‐conserved S‐FDTD scheme for variable coefficient Maxwell's equations in disk domains

In this paper, we analyze the energy‐conserved splitting finite‐difference time‐domain (FDTD) scheme for variable coefficient Maxwell's equations in two‐dimensional disk domains. The approach is energy‐conserved, unconditionally stable, and effective. We strictly prove that the EC‐S‐FDTD scheme for the variable coefficient Maxwell's equations in disk domains is of second order accuracy both in time and space. It is also strictly proved that the scheme is energy‐conserved, and the discrete divergence‐free is of second order convergence. Numerical experiments confirm the theoretical results, and practical test is simulated as well to demonstrate the efficiency of the proposed EC‐S‐FDTD scheme. Copyright © 2015 John Wiley & Sons, Ltd.     

Numerical analysis for generalized Forchheimer flows of slightly compressible fluids

In this article, we will consider the generalized Forchheimer flows for slightly compressible fluids. Using Muskat's and Ward's general form of Forchheimer equations, we describe the fluid dynamics by a nonlinear degenerate parabolic equation for density. The long‐time numerical approximation of the nonlinear degenerate parabolic equation with time dependent boundary conditions is studied. The stability for all time is established in a continuous time scheme and a discrete backward Euler scheme. A Gronwall's inequality‐type is used to study the asymptotic behavior of the solution. Error estimates for the solution are derived for both continuous and discrete time procedures. Numerical experiments confirm the theoretical analysis regarding convergence rates.     

An Integrated Research on Children's Construction of Meaningful,Symbolic, Partitioning-related Conceptions and the Teacher's Role in Fostering That Learning

A teaching experiment was conducted with two fourth graders to study the co-emergence of teaching and children's construction of fraction knowledge. The children's learning, i.e., modifications in their fraction schemes, was fostered through working on tasks in a computer microworld. The children advanced from thinking about a unit fraction as one of several equal parts in a whole (the equipartitioning scheme) to operating with a unit fraction as a symbolized, iterable part the magnitude of which is based on the numerosity of the partitioned whole (the partitive fraction scheme). The paper interweaves an analysis of children's construction of partitioning-related symbolic conceptions of fractions with an analysis of the teaching—planning and using tasks—that fosters such an advancement by introducing fraction words and numerals in the context of the children's partitioning activities.     

大气环流闭合方程组的总能量守恒格式及全球预报误差的严格估计

本文讨论大气环流闭合方程组,由于同时考虑了热传导效应,内摩擦效应及表达动能向内能转化的耗散项,因此符合总能量守恒律,文中对这一方程组建立了加权平均守恒型差分格式,并证明当选择最优参数时,它满足离散形式的总能量守恒律,通常的二次守恒格式是其次优的情况,文中还综合应用了Jessen不等式,Hardy不等式等等,从而严格证明了在一定条件下,存在t₀>0,当t     

Analysis of DG approximations for Stokes problem based on velocity‐pseudostress formulation

In this article, we first discuss the well posedness of a modified LDG scheme of Stokes problem, considering a velocity‐pseudostress formulation. The difficulty here relies on the fact that the application of classical Babu?ka‐Brezzi theory is not easy, so we proceed in a nonstandard way. For uniqueness, we apply a discrete version of Fredholm's alternative theorem, while the a priori error analysis is done introducing suitable projections of exact solution. As a result, we prove that the method is convergent, and under suitable regularity assumptions on the exact solution, the optimal rate of convergence is guaranteed. Next, we explore two stabilizations to the previous scheme, by adding least squares type terms. For these cases, well posedness and the a priori error estimates are proved by the application of standard theory. We end this work with some numerical experiments considering our third scheme, whose results are in agreement with the theoretical properties we deduce.© 2017 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 33: 1540–1564, 2017     

A Finite Volume Scheme for Savage-Hutter Equations on Unstructured Grids

A Godunov-type finite volume scheme on unstructured grids is proposed to numerically solve the Savage-Hutter equations in curvilinear coordinate. We show the direct observation that the model isn't a Galilean invariant system. At the cell boundary, the modified Harten-Lax-van Leer (HLL) approximate Riemann solver is adopted to calculate the numerical flux. The modified HLL flux is not troubled by the lack of Galilean invariance of the model and it is helpful to handle discontinuities at free interface. Rigidly the system is not always a hyperbolic system due to the dependence of flux on the velocity gradient. Even so, our numerical results still show quite good agreements with reference solutions. The simulations for granular avalanche flows with shock waves indicate that the scheme is applicable.     

一个安全有效的智能卡远程用户认证方案

2000年,Hwang和Li提出了一个新的智能卡远程用户认证方案,随后Chan和Cheng对该方案进行了成功的攻击.最近Shen,Lin和Hwang针对该方案提出了一种不同的攻击方法,并提供了一个改进方案用于抵御这些攻击.2003年,Leung等认为Shen-Lin-Hwang改进方案仍然不能抵御Chan和Cheng的攻击,他们用改进后的Chang-Hwang攻击方法进行了攻击.文中主要在Hwang-Li方案的基础上,提出了一个新的远程用户认证方案,该方案主要在注册阶段和登录阶段加强了安全性,抵御了类似Chan-Cheng和Chang-Hwang的攻击.     

A modified Green's function for the internal gravitational wave equation in a layer of a stratified medium with a constant shear flow

The construction of a modified Green's function for the internal gravitational wave (IGW) equation in a layer of a stratified medium when there are constant mean shear flows is considered and the basic properties of the corresponding eigenvalue problems and the modified eigenfunctions and eigenvalues are investigated. It is shown that each mode of the modified Green's function consists of a sum of three terms describing (1) the IGWs that propagate from the source, (2) the effects of a time varying source, localized in a certain neighbourhood of it, and (3) the effects of the displacement of the fluid (an internal discontinuity) caused by the source. The resulting expressions are analysed out for a constant and oscillating source of the generation of IGWs in which each of the terms of Green's function is represented in the form of simple quadratures.