Search

This paper presents an asymptotic preserving (AP) implicit-explicit (IMEX) scheme for solving the quantum BGK equation using the Hermite spectral method. The distribution function is expanded in a series of Hermite polynomials, with the Gaussian function serving as the weight function. The main challenge in this numerical scheme lies in efficiently expanding the quantum Maxwellian with the Hermite basis functions. To overcome this, we simplify the problem to the calculation of polylogarithms and propose an efficient algorithm to handle it, utilizing the Gauss-Hermite quadrature. Several numerical simulations, including a spatially 2D lid-driven cavity flow, demonstrate the AP property and remarkable efficiency of this method.

We propose a Hermite spectral method for the inelastic Boltzmann equation, which makes two-dimensional periodic problem computation affordable by the hardware nowadays. The new algorithm is based on a Hermite expansion, where the expansion coefficients for the VHS model are reduced into several summations and can be derived exactly. Moreover, a new collision model is built with a combination of the quadratic collision operator and a linearized collision operator, which helps us to balance the computational cost and the accuracy. Various numerical experiments, including spatially two-dimensional simulations, demonstrate the accuracy and efficiency of this numerical scheme.

We introduce a numerical scheme for the full multi-species Boltzmann equation based on Hermite spectral method. With the proper choice of expansion centers for different species, a practical algorithm is derived to evaluate the complicated multi-species binary collision operator. New collision models are built by combining the quadratic collision model and the simple BGK collision model under the framework of the Hermite spectral method, which enables us to balance the computational cost and accuracy. Several numerical experiments are implemented to validate the dramatic efficiency of this new Hermite spectral method. Moreover, we can handle the problems with as many as 100 species, which is far beyond the capability of the state-of-art algorithms.

In this paper, a new numerical method based on adaptive gradient descent optimizers is provided for computing the implied volatility from the Black-Scholes (B-S) option pricing model. It is shown that the new method is more accurate than the close form approximation. Compared with the Newton-Raphson method, the new method obtains a reliable rate of convergence and tends to be less sensitive to the beginning point., Comment: Our implement of Newton-Raphson iteration has defects. After correcting the code implement, we find Newton-Raphson won't be non-convergent. See https://github.com/cloudy-sfu/Newton-Raphson-Implied-Volatility for details

The problem of the new energy economy has become a global hot issue. This study examines the causal relationship between the ratio of thermal power in total power generation (RTPG) and economic growth (GDP) in the western and central China by using the heterogeneous mixed panel Granger causality approach that accounts for both slope heterogeneity and cross-sectional dependence. For the overall panel, the empirical findings support the presence of unidirectional causality running from GDP to RTPG (in northwest China), and from RTPG to GDP (in central). At the provincial level, there is causality from GDP to RTPG in NeiMongol and Ningxia, and causality from RTPG to GDP in Shanxi, Anhui, and Jiangxi. As for the cross regions relationships, we find that GDP (in western) Granger-cause RTPG (in central), and RTPG (in southwest) Granger-cause GDP (in central and northwest). Moreover, panel regressions show the negative impact from GDP to RTPG in the northwest, and RTPG to GDP in the central. However, RTPG has a positive influence on GDP in the northwest. Therefore, to improve economic development without compromising the regions' competitiveness in central and western China, we can adjust the power generation structure, and increase investments in the renewable energy supply and energy efficiency., Comment: The results of this paper cannot be reproduced on the 2018 data. Some updated data obtained from government statistical reports may overturn the conclusions of this article. So we consider changing the method of analysis, conduct further research, and rewrite this article

A two-fluid model is derived from the plasma kinetic equations using the moment model reduction method. The moment method we adopt was recently developed with a globally hyperbolic regularization where the moment model attained is locally well-posed in time. Based on the hyperbolic moment model with well-posedness, the Maxwellian iteration method is utilized to get the closure relations for the resulted two-fluid model. By taking the Shakhov collision operator in the Maxwellian iteration, the two-fluid model inherits the correct Prandtl number from the plasma kinetic equations. The new model is formally the same as the five-moment two-fluid model except for the closure relations, where the pressure tensor is anisotropic and the heat flux is presented. This provides the model the capacity to depict problems with anisotropic pressure tensor and large heat flux.