This paper considers pricing European options under the well-known of SVJ model of Bates and related computational methods. According to the no-arbitrage principle, we first derive a partial differential equation that...This paper considers pricing European options under the well-known of SVJ model of Bates and related computational methods. According to the no-arbitrage principle, we first derive a partial differential equation that the value of any European contingent claim should satisfy, where the asset price obeys the SVJ model. This equation is numerically solved by using the implicit- explicit backward difference method and time semi-discretization. In order to explain the validity of our method, the stability of time semi-discretization scheme is also proved. Finally, we use a simulation example to illustrate the efficiency of the method.展开更多
Numerical diffusion and oscillatory behavior characteristics are averted applying numerical solutions of advection-diffusion equation are themselves immensely sophisticated. In this paper, two numerical methods have b...Numerical diffusion and oscillatory behavior characteristics are averted applying numerical solutions of advection-diffusion equation are themselves immensely sophisticated. In this paper, two numerical methods have been used to solve the advection diffusion equation. We use an explicit finite difference scheme for the advection diffusion equation and semi-discretization on the spatial variable for advection-diffusion equation yields a system of ordinary differential equations solved by Euler’s method. Numerical assessment has been executed with specified initial and boundary conditions, for which the exact solution is known. We compare the solutions of the advection diffusion equation as well as error analysis for both schemes.展开更多
Using the weight coefficient method, we first discuss semi-discrete Hilbert-type inequalities, and then discuss boundedness of integral and discrete operators and operator norm estimates based on Hilbert-type inequali...Using the weight coefficient method, we first discuss semi-discrete Hilbert-type inequalities, and then discuss boundedness of integral and discrete operators and operator norm estimates based on Hilbert-type inequalities in weighted Lebesgue space and weighted normed sequence space.展开更多
In this paper, the improved Jacobian elliptic function expansion approach is extended and applied to constructing discrete solutions of the semi-discrete coupled modified Korteweg de Vries (mKdV) equations with the ...In this paper, the improved Jacobian elliptic function expansion approach is extended and applied to constructing discrete solutions of the semi-discrete coupled modified Korteweg de Vries (mKdV) equations with the aid of the symbolic computation system Maple. Some new discrete Jacobian doubly periodic solutions are obtained. When the modulus m →1, these doubly periodic solutions degenerate into the corresponding solitary wave solutions, including kink-type, bell-type and other types of excitations.展开更多
Fluidelastic instability is destructive in tube bundles subjected to cross flow.Flow channel model proposed by Leaver and Weaver is well used for modeling this problem.However,as the tube motion is supposed to be harm...Fluidelastic instability is destructive in tube bundles subjected to cross flow.Flow channel model proposed by Leaver and Weaver is well used for modeling this problem.However,as the tube motion is supposed to be harmonic,it may not simulate the general dynamic behaviors of tubes.To improve this,a model with arbitrary tube motion is proposed by Hassan and Hayder.While,due to involving in the time delay term,the stability problem cannot be solved by the eigenvalue scheme,and time domain responses of the tube have to be obtained to assess the instability threshold.To overcome this weakness,a new approach based on semi-discretizing method(SDM)is proposed in this study to make the instability threshold be predicted by eigenvalues directly.The motion equation of tube is built with considering the arbitrary tube motion and the time delay between fluid flow and tube vibration.A time delay integral term is derived and the SDM is employed to construct a transfer matrix,which transforms the infinite dimensional eigenvalue problem into a finite one.Hence the stability problem become solvable accordingly.With the proposed method,the instability threshold of a typical square tube array model is predicted,and the influences of system parameters on stability are also discussed.With comparing with prior works,it shows significant efficiency improvement in prediction of the instability threshold of tube bundles.展开更多
A numerical method for American options pricing on assets under the Heston stochastic volatility model is developed.A preliminary transformation is applied to remove the mixed derivative term avoiding known numerical ...A numerical method for American options pricing on assets under the Heston stochastic volatility model is developed.A preliminary transformation is applied to remove the mixed derivative term avoiding known numerical drawbacks and reducing computational costs.Free boundary is treated by the penalty method.Transformed nonlinear partial differential equation is solved numerically by using the method of lines.For full discretization the exponential time differencing method is used.Numerical analysis establishes the stability and positivity of the proposed method.The numerical convergence behaviour and effectiveness are investigated in extensive numerical experiments.展开更多
In this paper a numerical investigation on the non-Newtonian flow problem is conducted, in order to shed further light on the mathematical and virtual test methods in the auto-crash safety analysis. The accurate mathe...In this paper a numerical investigation on the non-Newtonian flow problem is conducted, in order to shed further light on the mathematical and virtual test methods in the auto-crash safety analysis. The accurate mathematical prediction would supply ultimate research tool for the passive safety analysis in such a scale.展开更多
A nonconforming mixed finite element method for nonlinear hyperbolic equations is discussed. Existence and uniqueness of the solution to the discrete problem are proved. Priori estimates of optimal order are derived f...A nonconforming mixed finite element method for nonlinear hyperbolic equations is discussed. Existence and uniqueness of the solution to the discrete problem are proved. Priori estimates of optimal order are derived for both the displacement and the stress.展开更多
The article gives a semi-discrete method for solving high-dimension wave equationBy the method, high-dimension wave equation is converted by, means of diseretizationinto I-D wave equation system which is well-posed. T...The article gives a semi-discrete method for solving high-dimension wave equationBy the method, high-dimension wave equation is converted by, means of diseretizationinto I-D wave equation system which is well-posed. The convergence of the semidijcrete method is given. The numerical calculating resulis show that the speed of convergence is high.展开更多
This paper considers weak Galerkin finite element approximations on polygonal/polyhedral meshes for a quasistatic Maxwell viscoelastic model.The spatial discretization uses piecewise polynomials of degree k(k≥1)for t...This paper considers weak Galerkin finite element approximations on polygonal/polyhedral meshes for a quasistatic Maxwell viscoelastic model.The spatial discretization uses piecewise polynomials of degree k(k≥1)for the stress approximation,degree k+1 for the velocity approximation,and degree k for the numerical trace of velocity on the inter-element boundaries.The temporal discretization in the fully discrete method adopts a backward Euler difference scheme.We show the existence and uniqueness of the semi-discrete and fully discrete solutions,and derive optimal a priori error estimates.Numerical examples are provided to support the theoretical analysis.展开更多
A numerical two-dimensional shallow water method was based on method for solving the equations was presented. This the third-order genuinely multidimensional semi-discrete central scheme for spatial discretization an...A numerical two-dimensional shallow water method was based on method for solving the equations was presented. This the third-order genuinely multidimensional semi-discrete central scheme for spatial discretization and the optimal third-order Strong Stability Preserving (SSP) Runge-Kutta method for time integration. The third-order compact Central Weighted Essentially Non-Oscillatory (CWENO) reconstruction was adopted to guarantee the non-oscillatory behavior of the presented scheme and improve the resolution. Two kinds of source terms were considered in this work. They were evaluated using different approaches. The resulting scheme does not require Riemann solvers or characteristic decomposition, hence it retains all the attractive features of central schemes such as simplicity and high resolution. To evaluate the performance of the presented scheme, several numerical examples were tested. The results demonstrate that our method is efficient, stable and robust.展开更多
In this work,the Lie point symmetries of the inhomogeneous Toda lattice equation are obtained by semi-discrete exterior calculus,which is a semi-discrete version of Harrison and Estabrook’s geometric approach.A four-...In this work,the Lie point symmetries of the inhomogeneous Toda lattice equation are obtained by semi-discrete exterior calculus,which is a semi-discrete version of Harrison and Estabrook’s geometric approach.A four-dimensional Lie algebra and its one-,two-and three-dimensional subalgebras are given.Two similarity reductions of the inhomogeneous Toda lattice equation are obtained by using the symmetry vectors.展开更多
Semi-discrete and fully discrete mixedfinite element methods are consid-ered for Maxwell-model-based problems of wave propagation in linear viscoelastic solid.This mixedfinite element framework allows the use of a large...Semi-discrete and fully discrete mixedfinite element methods are consid-ered for Maxwell-model-based problems of wave propagation in linear viscoelastic solid.This mixedfinite element framework allows the use of a large class of exist-ing mixed conformingfinite elements for elasticity in the spatial discretization.In the fully discrete scheme,a Crank-Nicolson scheme is adopted for the approximation of the temporal derivatives of stress and velocity variables.Error estimates of the semi-discrete and fully discrete schemes,as well as an unconditional stability result for the fully discrete scheme,are derived.Numerical experiments are provided to verify the theoretical results.展开更多
Whirling has been adopted for the cost-effective machining of blade-shape components in addition to traditional end milling and flank milling processes.To satisfy the requirements of rotary forming in the blade whirli...Whirling has been adopted for the cost-effective machining of blade-shape components in addition to traditional end milling and flank milling processes.To satisfy the requirements of rotary forming in the blade whirling process,the workpiece must be clamped at both ends in suspension and rotated slowly during machining,which complicates the dynamics.This study aims to identify the dynamic characteristics within the blade whirling operation and present strategies for stability prediction.In this study,the dynamic characteristics of a whirling system are modeled by assuming symmetric and asymmetric parameters.Theoretical prediction frequency response function(FRF)results are compared with experimental results.Moreover,semi-discretization stability lobe diagrams(SLDs)obtained using the dynamic parameters of these models are investigated experimentally.The results show that the asymmetric model is more suitable for describing the whirling system,whereas the symmetric model presents limitations associated with the frequency range and location of measuring points.Finally,a set of airfoil propeller blade whirling operations is conducted to verify the prediction accuracy.展开更多
In this paper, nonconforming quasi-Wilson finite element approximation to a class of nonlinear sine-Gordan equations is discussed. Based on the known higher accuracy results of bilinear element and different technique...In this paper, nonconforming quasi-Wilson finite element approximation to a class of nonlinear sine-Gordan equations is discussed. Based on the known higher accuracy results of bilinear element and different techniques from the existing literature, it is proved that the inner product △↓(u - Ih^1u), △↓vh) and the consistency error can be estimated as order O(h^2) in broken H^1 - norm/L^2 - norm when u ∈ H^3(Ω)/H^4(Ω), where Ih^1u is the bilinear interpolation of u, Vh belongs to the quasi-Wilson finite element space. At the same time, the superclose result with order O(h^2) for semi-discrete scheme under generalized rectangular meshes is derived. Furthermore, a fully-discrete scheme is proposed and the corresponding error estimate of order O(h^2 + τ^2) is obtained for the rectangular partition when u ∈ H^4(Ω), which is as same as that of the bilinear element with ADI scheme and one order higher than that of the usual analysis on nonconforming finite elements.展开更多
EQrot nonconforming finite element approximation to a class of nonlinear dual phase lagging heat conduction equations is discussed for semi-discrete and fully-discrete schemes. By use of a special property, that is, t...EQrot nonconforming finite element approximation to a class of nonlinear dual phase lagging heat conduction equations is discussed for semi-discrete and fully-discrete schemes. By use of a special property, that is, the consistency error of this element is of order O(h2) one order higher than its interpolation error O(h), the superclose results of order O(h2) in broken Hi-norm are obtained. At the same time, the global superconvergence in broken Hi-norm is deduced by interpolation postprocessing technique. Moreover, the extrapolation result with order O(h4) is derived by constructing a new interpolation postprocessing operator and extrapolation scheme based on the known asymptotic expansion formulas of EQrot element. Finally, optimal error estimate is gained for a proposed fully-discrete scheme by different approaches from the previous literature.展开更多
The lowest order H1-Galerkin mixed finite element method (for short MFEM) is proposed for a class of nonlinear sine-Gordon equations with the simplest bilinear rectangular element and zero order Raviart- Thomas elem...The lowest order H1-Galerkin mixed finite element method (for short MFEM) is proposed for a class of nonlinear sine-Gordon equations with the simplest bilinear rectangular element and zero order Raviart- Thomas element. Base on the interpolation operator instead of the traditional Ritz projection operator which is an indispensable tool in the traditional FEM analysis, together with mean-value technique and high accuracy analysis, the superclose properties of order O(h2)/O(h2 + r2) in Hi-norm and H(div; Ω)-norm axe deduced for the semi-discrete and the fully-discrete schemes, where h, r- denote the mesh size and the time step, respectively, which improve the results in the previous literature.展开更多
文摘This paper considers pricing European options under the well-known of SVJ model of Bates and related computational methods. According to the no-arbitrage principle, we first derive a partial differential equation that the value of any European contingent claim should satisfy, where the asset price obeys the SVJ model. This equation is numerically solved by using the implicit- explicit backward difference method and time semi-discretization. In order to explain the validity of our method, the stability of time semi-discretization scheme is also proved. Finally, we use a simulation example to illustrate the efficiency of the method.
文摘Numerical diffusion and oscillatory behavior characteristics are averted applying numerical solutions of advection-diffusion equation are themselves immensely sophisticated. In this paper, two numerical methods have been used to solve the advection diffusion equation. We use an explicit finite difference scheme for the advection diffusion equation and semi-discretization on the spatial variable for advection-diffusion equation yields a system of ordinary differential equations solved by Euler’s method. Numerical assessment has been executed with specified initial and boundary conditions, for which the exact solution is known. We compare the solutions of the advection diffusion equation as well as error analysis for both schemes.
基金Supported by Guangdong Basic and Applied Basic Research Foundation(Grant No.2022A1515012429)Guangzhou Huashang College Research Team Project(Grant No.2021HSKT03)。
文摘Using the weight coefficient method, we first discuss semi-discrete Hilbert-type inequalities, and then discuss boundedness of integral and discrete operators and operator norm estimates based on Hilbert-type inequalities in weighted Lebesgue space and weighted normed sequence space.
基金Project supported by the National Natural Science Foundation of China (Grant No 10272071) and the Natural Science Foundation of Zhejiang Lishui University of China (Grant Nos KZ05004 and KY06024).
文摘In this paper, the improved Jacobian elliptic function expansion approach is extended and applied to constructing discrete solutions of the semi-discrete coupled modified Korteweg de Vries (mKdV) equations with the aid of the symbolic computation system Maple. Some new discrete Jacobian doubly periodic solutions are obtained. When the modulus m →1, these doubly periodic solutions degenerate into the corresponding solitary wave solutions, including kink-type, bell-type and other types of excitations.
基金The support from the National Natural Science Foundation of China(No.11672179)is greatly acknowledged.
文摘Fluidelastic instability is destructive in tube bundles subjected to cross flow.Flow channel model proposed by Leaver and Weaver is well used for modeling this problem.However,as the tube motion is supposed to be harmonic,it may not simulate the general dynamic behaviors of tubes.To improve this,a model with arbitrary tube motion is proposed by Hassan and Hayder.While,due to involving in the time delay term,the stability problem cannot be solved by the eigenvalue scheme,and time domain responses of the tube have to be obtained to assess the instability threshold.To overcome this weakness,a new approach based on semi-discretizing method(SDM)is proposed in this study to make the instability threshold be predicted by eigenvalues directly.The motion equation of tube is built with considering the arbitrary tube motion and the time delay between fluid flow and tube vibration.A time delay integral term is derived and the SDM is employed to construct a transfer matrix,which transforms the infinite dimensional eigenvalue problem into a finite one.Hence the stability problem become solvable accordingly.With the proposed method,the instability threshold of a typical square tube array model is predicted,and the influences of system parameters on stability are also discussed.With comparing with prior works,it shows significant efficiency improvement in prediction of the instability threshold of tube bundles.
基金This work has been supported by the Spanish Ministerio de Economía,Industria y Competitividad(MINECO),the Agencia Estatal de Investigación(AEI)and Fondo Europeo de Desarrollo Regional(FEDER UE)grant MTM2017-89664-P.
文摘A numerical method for American options pricing on assets under the Heston stochastic volatility model is developed.A preliminary transformation is applied to remove the mixed derivative term avoiding known numerical drawbacks and reducing computational costs.Free boundary is treated by the penalty method.Transformed nonlinear partial differential equation is solved numerically by using the method of lines.For full discretization the exponential time differencing method is used.Numerical analysis establishes the stability and positivity of the proposed method.The numerical convergence behaviour and effectiveness are investigated in extensive numerical experiments.
基金supported by the National Natural Science Foundation of China (Grant No.10871225) the Pujing Talents’ Project of Shanghai (Grant No.PJ[2006]118)the E-Institutes of Shanghai Municipal Education Commission (Grant No.E03004)
文摘In this paper a numerical investigation on the non-Newtonian flow problem is conducted, in order to shed further light on the mathematical and virtual test methods in the auto-crash safety analysis. The accurate mathematical prediction would supply ultimate research tool for the passive safety analysis in such a scale.
文摘A nonconforming mixed finite element method for nonlinear hyperbolic equations is discussed. Existence and uniqueness of the solution to the discrete problem are proved. Priori estimates of optimal order are derived for both the displacement and the stress.
文摘The article gives a semi-discrete method for solving high-dimension wave equationBy the method, high-dimension wave equation is converted by, means of diseretizationinto I-D wave equation system which is well-posed. The convergence of the semidijcrete method is given. The numerical calculating resulis show that the speed of convergence is high.
基金This work was supported by the National Natural Science Foundation of China(Grant No.12171340).
文摘This paper considers weak Galerkin finite element approximations on polygonal/polyhedral meshes for a quasistatic Maxwell viscoelastic model.The spatial discretization uses piecewise polynomials of degree k(k≥1)for the stress approximation,degree k+1 for the velocity approximation,and degree k for the numerical trace of velocity on the inter-element boundaries.The temporal discretization in the fully discrete method adopts a backward Euler difference scheme.We show the existence and uniqueness of the semi-discrete and fully discrete solutions,and derive optimal a priori error estimates.Numerical examples are provided to support the theoretical analysis.
基金Project supported by the National Natural Science Foundation of China (Grant No: 60134010).
文摘A numerical two-dimensional shallow water method was based on method for solving the equations was presented. This the third-order genuinely multidimensional semi-discrete central scheme for spatial discretization and the optimal third-order Strong Stability Preserving (SSP) Runge-Kutta method for time integration. The third-order compact Central Weighted Essentially Non-Oscillatory (CWENO) reconstruction was adopted to guarantee the non-oscillatory behavior of the presented scheme and improve the resolution. Two kinds of source terms were considered in this work. They were evaluated using different approaches. The resulting scheme does not require Riemann solvers or characteristic decomposition, hence it retains all the attractive features of central schemes such as simplicity and high resolution. To evaluate the performance of the presented scheme, several numerical examples were tested. The results demonstrate that our method is efficient, stable and robust.
基金Supported by National Natural Science Foundation of China under Grant Nos.11375030,11472315Department of Science and Technology of Henan Province under Grant No.162300410223Beijing Finance Funds of Natural Science Program for Excellent Talents under Grant No.2014000026833ZK19
文摘In this work,the Lie point symmetries of the inhomogeneous Toda lattice equation are obtained by semi-discrete exterior calculus,which is a semi-discrete version of Harrison and Estabrook’s geometric approach.A four-dimensional Lie algebra and its one-,two-and three-dimensional subalgebras are given.Two similarity reductions of the inhomogeneous Toda lattice equation are obtained by using the symmetry vectors.
基金supported in part by National Natural Science Foundation of China(No.11771312).
文摘Semi-discrete and fully discrete mixedfinite element methods are consid-ered for Maxwell-model-based problems of wave propagation in linear viscoelastic solid.This mixedfinite element framework allows the use of a large class of exist-ing mixed conformingfinite elements for elasticity in the spatial discretization.In the fully discrete scheme,a Crank-Nicolson scheme is adopted for the approximation of the temporal derivatives of stress and velocity variables.Error estimates of the semi-discrete and fully discrete schemes,as well as an unconditional stability result for the fully discrete scheme,are derived.Numerical experiments are provided to verify the theoretical results.
基金the Shandong Provincial Natural Science Foundation,China(Grant No.ZR2017MEE021)for the financial support toward this research.
文摘Whirling has been adopted for the cost-effective machining of blade-shape components in addition to traditional end milling and flank milling processes.To satisfy the requirements of rotary forming in the blade whirling process,the workpiece must be clamped at both ends in suspension and rotated slowly during machining,which complicates the dynamics.This study aims to identify the dynamic characteristics within the blade whirling operation and present strategies for stability prediction.In this study,the dynamic characteristics of a whirling system are modeled by assuming symmetric and asymmetric parameters.Theoretical prediction frequency response function(FRF)results are compared with experimental results.Moreover,semi-discretization stability lobe diagrams(SLDs)obtained using the dynamic parameters of these models are investigated experimentally.The results show that the asymmetric model is more suitable for describing the whirling system,whereas the symmetric model presents limitations associated with the frequency range and location of measuring points.Finally,a set of airfoil propeller blade whirling operations is conducted to verify the prediction accuracy.
文摘In this paper, nonconforming quasi-Wilson finite element approximation to a class of nonlinear sine-Gordan equations is discussed. Based on the known higher accuracy results of bilinear element and different techniques from the existing literature, it is proved that the inner product △↓(u - Ih^1u), △↓vh) and the consistency error can be estimated as order O(h^2) in broken H^1 - norm/L^2 - norm when u ∈ H^3(Ω)/H^4(Ω), where Ih^1u is the bilinear interpolation of u, Vh belongs to the quasi-Wilson finite element space. At the same time, the superclose result with order O(h^2) for semi-discrete scheme under generalized rectangular meshes is derived. Furthermore, a fully-discrete scheme is proposed and the corresponding error estimate of order O(h^2 + τ^2) is obtained for the rectangular partition when u ∈ H^4(Ω), which is as same as that of the bilinear element with ADI scheme and one order higher than that of the usual analysis on nonconforming finite elements.
基金Supported by the National Natural Science Foundation of China (Nos. 10971203 11101381)+3 种基金Tianyuan Mathe-matics Foundation of National Natural Science Foundation of China (No. 11026154)Natural Science Foundation of Henan Province (No. 112300410026)Natural Science Foundation of the Education Department of Henan Province (Nos. 2011A110020 12A110021)
文摘EQrot nonconforming finite element approximation to a class of nonlinear dual phase lagging heat conduction equations is discussed for semi-discrete and fully-discrete schemes. By use of a special property, that is, the consistency error of this element is of order O(h2) one order higher than its interpolation error O(h), the superclose results of order O(h2) in broken Hi-norm are obtained. At the same time, the global superconvergence in broken Hi-norm is deduced by interpolation postprocessing technique. Moreover, the extrapolation result with order O(h4) is derived by constructing a new interpolation postprocessing operator and extrapolation scheme based on the known asymptotic expansion formulas of EQrot element. Finally, optimal error estimate is gained for a proposed fully-discrete scheme by different approaches from the previous literature.
基金Supported in part by the National Natural Science Foundation of China under Grant Nos.11671369,11271340the Natural Science Foundation of the Education Department of Henan Province under Grant Nos.14A110009,16A110022
文摘The lowest order H1-Galerkin mixed finite element method (for short MFEM) is proposed for a class of nonlinear sine-Gordon equations with the simplest bilinear rectangular element and zero order Raviart- Thomas element. Base on the interpolation operator instead of the traditional Ritz projection operator which is an indispensable tool in the traditional FEM analysis, together with mean-value technique and high accuracy analysis, the superclose properties of order O(h2)/O(h2 + r2) in Hi-norm and H(div; Ω)-norm axe deduced for the semi-discrete and the fully-discrete schemes, where h, r- denote the mesh size and the time step, respectively, which improve the results in the previous literature.