Source author record

Tao Xiong

Tao Xiong appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

25works
12topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

25 published item(s)

preprint2022arXiv

sqSGD: Locally Private and Communication Efficient Federated Learning

Federated learning (FL) is a technique that trains machine learning models from decentralized data sources. We study FL under local notions of privacy constraints, which provides strong protection against sensitive data disclosures via obfuscating the data before leaving the client. We identify two major concerns in designing practical privacy-preserving FL algorithms: communication efficiency and high-dimensional compatibility. We then develop a gradient-based learning algorithm called \emph{sqSGD} (selective quantized stochastic gradient descent) that addresses both concerns. The proposed algorithm is based on a novel privacy-preserving quantization scheme that uses a constant number of bits per dimension per client. Then we improve the base algorithm in three ways: first, we apply a gradient subsampling strategy that simultaneously offers better training performance and smaller communication costs under a fixed privacy budget. Secondly, we utilize randomized rotation as a preprocessing step to reduce quantization error. Thirdly, an adaptive gradient norm upper bound shrinkage strategy is adopted to improve accuracy and stabilize training. Finally, the practicality of the proposed framework is demonstrated on benchmark datasets. Experiment results show that sqSGD successfully learns large models like LeNet and ResNet with local privacy constraints. In addition, with fixed privacy and communication level, the performance of sqSGD significantly dominates that of various baseline algorithms.

preprint2021arXiv

High order well-balanced asymptotic preserving finite difference WENO schemes for the shallow water equations in all Froude numbers

In this paper, high order semi-implicit well-balanced and asymptotic preserving finite difference WENO schemes are proposed for the shallow water equations with a non-flat bottom topography. We consider the Froude number ranging from O(1) to 0, which in the zero Froude limit becomes the "lake equations" for balanced flow without gravity waves. We apply a well-balanced finite difference WENO reconstruction, coupled with a stiffly accurate implicit-explicit (IMEX) Runge-Kutta time discretization. The resulting semi-implicit scheme can be shown to be well-balanced, asymptotic preserving (AP) and asymptotically accurate (AA) at the same time. Both one- and two-dimensional numerical results are provided to demonstrate the high order accuracy, AP property and good performance of the proposed methods in capturing small perturbations of steady state solutions.

preprint2020arXiv

Conservative discontinuous Galerkin/Hermite Spectral Method for the Vlasov-Poisson System

We propose a class of conservative discontinuous Galerkin methods for the Vlasov-Poisson system written as a hyperbolic system using Hermite polynomials in the velocity variable. These schemes are designed to be systematically as accurate as one wants with provable conservation of mass and possibly total energy. Such properties in general are hard to achieve within other numerical method frameworks for simulating the Vlasov-Poisson system. The proposed scheme employs discontinuous Galerkin discretization for both the Vlasov and the Poisson equations, resulting in a consistent description of the distribution function and electric field. Numerical simulations are performed to verify the order of accuracy and conservation properties.

preprint2018arXiv

Sensitivity based Neural Networks Explanations

Although neural networks can achieve very high predictive performance on various different tasks such as image recognition or natural language processing, they are often considered as opaque "black boxes". The difficulty of interpreting the predictions of a neural network often prevents its use in fields where explainability is important, such as the financial industry where regulators and auditors often insist on this aspect. In this paper, we present a way to assess the relative input features importance of a neural network based on the sensitivity of the model output with respect to its input. This method has the advantage of being fast to compute, it can provide both global and local levels of explanations and is applicable for many types of neural network architectures. We illustrate the performance of this method on both synthetic and real data and compare it with other interpretation techniques. This method is implemented into an open-source Python package that allows its users to easily generate and visualize explanations for their neural networks.

preprint2016arXiv

Conservative Multi-Dimensional Semi-Lagrangian Finite Difference Scheme: Stability and Applications to the Kinetic and Fluid Simulations

In this paper, we propose a mass conservative semi-Lagrangian finite difference scheme for multi-dimensional problems without dimensional splitting. The semi-Lagrangian scheme, based on tracing characteristics backward in time from grid points, does not necessarily conserve the total mass. To ensure mass conservation, we propose a conservative correction procedure based on a flux difference form. Such procedure guarantees local mass conservation, while introducing time step constraints for stability. We theoretically investigate such stability constraints from an ODE point of view by assuming exact evaluation of spatial differential operators and from the Fourier analysis for linear PDEs. The scheme is tested by classical two dimensional linear passive-transport problems, such as linear advection, rotation and swirling deformation. The scheme is applied to solve the nonlinear Vlasov-Poisson system using a a high order tracing mechanism proposed in [Qiu and Russo, 2016]. Such high order characteristics tracing scheme is generalized to the nonlinear guiding center Vlasov model and incompressible Euler system. The effectiveness of the proposed conservative semi-Lagrangian scheme is demonstrated numerically by our extensive numerical tests.

preprint2016arXiv

High Order Hierarchical Asymptotic Preserving Nodal Discontinuous Galerkin IMEX Schemes For The BGK Equation

A class of high order asymptotic preserving (AP) schemes has been developed for the BGK equation in Xiong et. al. (2015) [37], which is based on the micro-macro formulation of the equation. The nodal discontinuous Galerkin (NDG) method with Lagrangian basis functions for spatial discretization and globally stiffly accurate implicit-explicit (IMEX) Runge-Kutta (RK) scheme as time discretization are introduced with asymptotic preserving properties. However, it is only necessary to solve the kinetic equation when the hydrodynamic description breaks down. Motivated by the recent work in Filbet and Rey (2015) [23], it is more naturally to construct a hierarchy scheme under the NDG-IMEX framework without hybridization, as the formal analysis in [37] shows that when $ε$ is small, the NDG-IMEX scheme becomes a local discontinuous Galerkin (LDG) scheme for the compressible Navier-Stokes equations, and when $ε=0$ it is a discontinuous Galerkin (DG) scheme for the compressible Euler equations. Moveover, we propose to combine the kinetic regime with the hydrodynamic regime including both the compressible Euler and Navier-Stokes equations. Numerical experiments demonstrate very decent performance of the new approach. In our numerics, all three regimes are clearly divided, leading to great savings in terms of the computational cost.

preprint2015arXiv

On $L_{n}$-Injective Modules and $L_{n}$-Injective Dimensions

Let $R$ be a ring, and $n$ a fixed nonnegative integer. An $R$-module $W$ is called $L_{n}$-injective if ${\rm Ext}_{R}^{1}(M,W)=0$ for any $R$-module $M$ with flat dimension at most $n$. In this paper, we prove first that ($\mathcal{F}_{n},\mathcal{L}_{n}$) is a complete hereditary cotorsion theory, where $\mathcal{F}_n$ (resp. $\mathcal{L}_n$) denotes the class of all $R$-modules with flat dimension at most $n$ (resp. $L_{n}$-injective $R$-modules). Then we introduce the $L_{n}$-injective dimension of a module and $L_n$-global dimension of a ring. Finally, over rings with weak global dimension $\leq n$, perfect rings, and $L_n$-hereditary rings, more properties and applications of $L_{n}$-injective modules, $L_{n}$-injective dimensions of modules and $\mathcal{L}_{n}$-global dimensions of rings are given.

preprint2014arXiv

A Maximum-Principle-Satisfying High-order Finite Volume Compact WENO Scheme for Scalar Conservation Laws

In this paper, a maximum-principle-satisfying finite volume compact scheme is proposed for solving scalar hyperbolic conservation laws. The scheme combines WENO schemes (Weighted Essentially Non-Oscillatory) with a class of compact schemes under a finite volume framework, in which the nonlinear WENO weights are coupled with lower order compact stencils. The maximum-principle-satisfying polynomial rescaling limiter in [Zhang and Shu, JCP, 2010] is adopted to construct the present schemes at each stage of an explicit Runge-Kutta method, without destroying high order accuracy and conservativity. Numerical examples for one and two dimensional problems including incompressible flows are presented to assess the good performance, maximum principle preserving, essentially non-oscillatory and highly accurate resolution of the proposed method.

preprint2014arXiv

A Positivity-preserving High Order Finite Volume Compact-WENO Scheme for Compressible Euler Equations

In this paper, a positivity-preserving fifth-order finite volume compact-WENO scheme is proposed for solving compressible Euler equations. As we know conservative compact finite volume schemes have high resolution properties while WENO (Weighted Essentially Non-Oscillatory) schemes are essentially non-oscillatory near flow discontinuities. We extend the main idea of WENO schemes to some classical compact finite volume schemes [32], where lower order compact stencils are combined with WENO nonlinear weights to get a higher order finite volume compact-WENO scheme. The newly developed positivity-preserving limiter [46,44] is used to preserve positive density and internal energy for compressible Euler equations of fluid dynamics. The HLLC (Harten, Lax, and van Leer with Contact) approximate Riemann solver [39,2] is used to get the numerical flux at the cell interfaces. Numerical tests are presented to demonstrate the high-order accuracy, positivity-preserving, high-resolution and robustness of the proposed scheme.

preprint2014arXiv

A PSO and Pattern Search based Memetic Algorithm for SVMs Parameters Optimization

Addressing the issue of SVMs parameters optimization, this study proposes an efficient memetic algorithm based on Particle Swarm Optimization algorithm (PSO) and Pattern Search (PS). In the proposed memetic algorithm, PSO is responsible for exploration of the search space and the detection of the potential regions with optimum solutions, while pattern search (PS) is used to produce an effective exploitation on the potential regions obtained by PSO. Moreover, a novel probabilistic selection strategy is proposed to select the appropriate individuals among the current population to undergo local refinement, keeping a well balance between exploration and exploitation. Experimental results confirm that the local refinement with PS and our proposed selection strategy are effective, and finally demonstrate effectiveness and robustness of the proposed PSO-PS based MA for SVMs parameters optimization.

preprint2014arXiv

Analysis of Asymptotic Preserving DG-IMEX Schemes for Linear Kinetic Transport Equations in a Diffusive Scaling

In this paper, some theoretical aspects will be addressed for the asymptotic preserving DG-IMEX schemes recently proposed in [J. Jang, F. Li, J.-M. Qiu and T. Xiong, submitted, arxiv:1306.0227] for kinetic transport equations under a diffusive scaling. We will focus on the methods that are based on discontinuous Galerkin (DG) spatial discretizations with the $P^k$ polynomial space and a first order IMEX temporal discretization, and apply them to two linear models: the telegraph equation and the one-group transport equation in slab geometry. In particular, we will establish uniform numerical stability with respect to Knudsen number $\varepsilon$ using energy methods, as well as error estimates for any given $\varepsilon$. When $\varepsilon\rightarrow 0$, a rigorous asymptotic analysis of the schemes is also obtained. Though the methods and the analysis are presented for one dimension in space, they can be generalized to higher dimensions directly.

preprint2014arXiv

Beyond One-Step-Ahead Forecasting: Evaluation of Alternative Multi-Step-Ahead Forecasting Models for Crude Oil Prices

An accurate prediction of crude oil prices over long future horizons is challenging and of great interest to governments, enterprises, and investors. This paper proposes a revised hybrid model built upon empirical mode decomposition (EMD) based on the feed-forward neural network (FNN) modeling framework incorporating the slope-based method (SBM), which is capable of capturing the complex dynamic of crude oil prices. Three commonly used multi-step-ahead prediction strategies proposed in the literature, including iterated strategy, direct strategy, and MIMO (multiple-input multiple-output) strategy, are examined and compared, and practical considerations for the selection of a prediction strategy for multi-step-ahead forecasting relating to crude oil prices are identified. The weekly data from the WTI (West Texas Intermediate) crude oil spot price are used to compare the performance of the alternative models under the EMD-SBM-FNN modeling framework with selected counterparts. The quantitative and comprehensive assessments are performed on the basis of prediction accuracy and computational cost. The results obtained in this study indicate that the proposed EMD-SBM-FNN model using the MIMO strategy is the best in terms of prediction accuracy with accredited computational load.

preprint2014arXiv

Does Restraining End Effect Matter in EMD-Based Modeling Framework for Time Series Prediction? Some Experimental Evidences

Following the "decomposition-and-ensemble" principle, the empirical mode decomposition (EMD)-based modeling framework has been widely used as a promising alternative for nonlinear and nonstationary time series modeling and prediction. The end effect, which occurs during the sifting process of EMD and is apt to distort the decomposed sub-series and hurt the modeling process followed, however, has been ignored in previous studies. Addressing the end effect issue, this study proposes to incorporate end condition methods into EMD-based decomposition and ensemble modeling framework for one- and multi-step ahead time series prediction. Four well-established end condition methods, Mirror method, Coughlin's method, Slope-based method, and Rato's method, are selected, and support vector regression (SVR) is employed as the modeling technique. For the purpose of justification and comparison, well-known NN3 competition data sets are used and four well-established prediction models are selected as benchmarks. The experimental results demonstrated that significant improvement can be achieved by the proposed EMD-based SVR models with end condition methods. The EMD-SBM-SVR model and EMD-Rato-SVR model, in particular, achieved the best prediction performances in terms of goodness of forecast measures and equality of accuracy of competing forecasts test.

preprint2014arXiv

Exploring gender differences on general and specific computer self-efficacy in mobile learning adoption

Reasons for contradictory findings regarding the gender moderate effect on computer self-efficacy in the adoption of e-learning/mobile learning are limited. Recognizing the multilevel nature of the computer self-efficacy (CSE), this study attempts to explore gender differences in the adoption of mobile learning, by extending the Technology Acceptance Model (TAM) with general and specific CSE. Data collected from 137 university students were tested against the research model using the structural equation modeling approach. The results suggest that there are significant gender differences in perceptions of general CSE, perceived ease of use and behavioral intention to use but no significant differences in specific CSE, perceived usefulness. Additionally, the findings reveal that specific CSE is more salient than general CSE in influencing perceived ease of use while general CSE seems to be the salient factor on perceived usefulness for both female and male combined. Moreover, general CSE was salient to determine the behavioral intention to use indirectly for female despite lower perception of general CSE than male's, and specific CSE exhibited stronger indirect effect on behavioral intention to use than general CSE for female despite similar perception of specific CSE as males'. These findings provide important implications for mobile learning adoption and usage.

preprint2014arXiv

High Order Asymptotic Preserving Nodal Discontinuous Galerkin IMEX Schemes for the BGK Equation

In this paper, we develop high-order asymptotic preserving (AP) schemes for the BGK equation in a hyperbolic scaling, which leads to the macroscopic models such as the Euler and compressible Navier-Stokes equations in the asymptotic limit. Our approaches are based on the so-called micro-macro formulation of the kinetic equation which involves a natural decomposition of the problem to the equilibrium and the non-equilibrium parts. The proposed methods are formulated for the BGK equation with constant or spatially variant Knudsen number. The new ingredients for the proposed methods to achieve high order accuracy are the following: we introduce discontinuous Galerkin (DG) discretization of arbitrary order of accuracy with nodal Lagrangian basis functions in space; we employ a high order globally stiffly accurate implicit-explicit (IMEX) Runge-Kutta (RK) scheme as time discretization. Two versions of the schemes are proposed: Scheme I is a direct formulation based on the micro-macro decomposition of the BGK equation, while Scheme II, motivated by the asymptotic analysis for the continuous problem, utilizes certain properties of the projection operator. Compared with Scheme I, Scheme II not only has better computational efficiency (the computational cost is reduced by half roughly), but also allows the establishment of a formal asymptotic analysis. Specifically, it is demonstrated that when $0<ε\ll 1$, Scheme II, up to $\mathcal{O}(ε^2)$, becomes a local DG discretization with an explicit RK method for the macroscopic compressible Navier-Stokes equations, a method in a similar spirit to the ones in [Bassi \& Rabey 1997, Cockburn \& Shu 1998]. Numerical results are presented for a wide range of Knudsen number to illustrate the effectiveness and high order accuracy of the methods.

preprint2014arXiv

High order maximum principle preserving discontinuous Galerkin method for convection-diffusion equations

In this paper, we propose to apply the parametrized maximum-principle-preserving (MPP) flux limiter in [Xiong et. al., JCP, 2013] to the discontinuous Galerkin (DG) method for solving the convection-diffusion equations. The feasibility of applying the MPP flux limiters to the DG solution of convection-diffusion problem is based on the fact that the cell averages for the DG solutions are updated in a conservative fashion (by using flux difference) even in the presence of diffusion terms. The main purpose of this paper is to address the difficulty of obtaining higher than second order accuracy while maintaining a discrete maximum principle for the DG method solving convection diffusion equations. We found that the proposed MPP flux limiter can be applied to arbitrarily high order DG method. Numerical evidence is presented to show that the proposed MPP flux limiter method does not adversely affect the desired high order accuracy, nor does it require restrictive time steps. Numerical experiments including incompressible Navier-Stokes equations demonstrate the high order accuracy preserving, the MPP performance, and the robustness of the proposed method.

preprint2014arXiv

High order maximum principle preserving finite volume method for convection dominated problems

In this paper, we investigate the application of the maximum principle preserving (MPP) parametrized flux limiters to the high order finite volume scheme with Runge-Kutta time discretization for solving convection dominated problems. Such flux limiter was originally proposed in [Xu, Math. Comp., 2013] and further developed in [Xiong et. al., J. Comp. Phys., 2013] for finite difference WENO schemes with Runge-Kutta time discretization for convection equations. The main idea is to limit the temporal integrated high order numerical flux toward a first order MPP monotone flux. In this paper, we generalize such flux limiter to high order finite volume methods solving convection-dominated problems, which is easy to implement and introduces little computational overhead. More importantly, for the first time in the finite volume setting, we provide a general proof that the proposed flux limiter maintains high order accuracy of the original WENO scheme for linear advection problems without any additional time step restriction. For general nonlinear convection-dominated problems, we prove that the proposed flux limiter introduces up to $O(Δx^3+Δt^3)$ modification to the high order temporal integrated flux in the original WENO scheme without extra time step constraint. We also numerically investigate the preservation of up to ninth order accuracy of the proposed flux limiter in a general setting. The advantage of the proposed method is demonstrated through various numerical experiments.

preprint2014arXiv

High Order Maximum Principle Preserving Semi-Lagrangian Finite Difference WENO schemes for the Vlasov Equation

In this paper, we propose the parametrized maximum principle preserving (MPP) flux limiter, originally developed in [Z. Xu, Math. Comp., (2013), in press], to the semi- Lagrangian finite difference weighted essentially non-oscillatory scheme for solving the Vlasov equation. The MPP flux limiter is proved to maintain up to fourth order accuracy for the semi-Lagrangian finite difference scheme without any time step restriction. Numerical studies on the Vlasov-Poisson system demonstrate the performance of the proposed method and its ability in preserving the positivity of the probability distribution function while maintaining the high order accuracy.

preprint2014arXiv

Interval Forecasting of Electricity Demand: A Novel Bivariate EMD-based Support Vector Regression Modeling Framework

Highly accurate interval forecasting of electricity demand is fundamental to the success of reducing the risk when making power system planning and operational decisions by providing a range rather than point estimation. In this study, a novel modeling framework integrating bivariate empirical mode decomposition (BEMD) and support vector regression (SVR), extended from the well-established empirical mode decomposition (EMD) based time series modeling framework in the energy demand forecasting literature, is proposed for interval forecasting of electricity demand. The novelty of this study arises from the employment of BEMD, a new extension of classical empirical model decomposition (EMD) destined to handle bivariate time series treated as complex-valued time series, as decomposition method instead of classical EMD only capable of decomposing one-dimensional single-valued time series. This proposed modeling framework is endowed with BEMD to decompose simultaneously both the lower and upper bounds time series, constructed in forms of complex-valued time series, of electricity demand on a monthly per hour basis, resulting in capturing the potential interrelationship between lower and upper bounds. The proposed modeling framework is justified with monthly interval-valued electricity demand data per hour in Pennsylvania-New Jersey-Maryland Interconnection, indicating it as a promising method for interval-valued electricity demand forecasting.

preprint2014arXiv

Multi-Step-Ahead Time Series Prediction using Multiple-Output Support Vector Regression

Accurate time series prediction over long future horizons is challenging and of great interest to both practitioners and academics. As a well-known intelligent algorithm, the standard formulation of Support Vector Regression (SVR) could be taken for multi-step-ahead time series prediction, only relying either on iterated strategy or direct strategy. This study proposes a novel multiple-step-ahead time series prediction approach which employs multiple-output support vector regression (M-SVR) with multiple-input multiple-output (MIMO) prediction strategy. In addition, the rank of three leading prediction strategies with SVR is comparatively examined, providing practical implications on the selection of the prediction strategy for multi-step-ahead forecasting while taking SVR as modeling technique. The proposed approach is validated with the simulated and real datasets. The quantitative and comprehensive assessments are performed on the basis of the prediction accuracy and computational cost. The results indicate that: 1) the M-SVR using MIMO strategy achieves the best accurate forecasts with accredited computational load, 2) the standard SVR using direct strategy achieves the second best accurate forecasts, but with the most expensive computational cost, and 3) the standard SVR using iterated strategy is the worst in terms of prediction accuracy, but with the least computational cost.

preprint2014arXiv

Multiple-output support vector regression with a firefly algorithm for interval-valued stock price index forecasting

Highly accurate interval forecasting of a stock price index is fundamental to successfully making a profit when making investment decisions, by providing a range of values rather than a point estimate. In this study, we investigate the possibility of forecasting an interval-valued stock price index series over short and long horizons using multi-output support vector regression (MSVR). Furthermore, this study proposes a firefly algorithm (FA)-based approach, built on the established MSVR, for determining the parameters of MSVR (abbreviated as FA-MSVR). Three globally traded broad market indices are used to compare the performance of the proposed FA-MSVR method with selected counterparts. The quantitative and comprehensive assessments are performed on the basis of statistical criteria, economic criteria, and computational cost. In terms of statistical criteria, we compare the out-of-sample forecasting using goodness-of-forecast measures and testing approaches. In terms of economic criteria, we assess the relative forecast performance with a simple trading strategy. The results obtained in this study indicate that the proposed FA-MSVR method is a promising alternative for forecasting interval-valued financial time series.

preprint2014arXiv

Parametrized Positivity Preserving Flux Limiters for the High Order Finite Difference WENO Scheme Solving Compressible Euler Equations

In this paper, we develop parametrized positivity satisfying flux limiters for the high order finite difference Runge-Kutta weighted essentially non-oscillatory (WENO) scheme solving compressible Euler equations to maintain positive density and pressure. Negative density and pressure, which often leads to simulation blow-ups or nonphysical solutions, emerges from many high resolution computations in some extreme cases. The methodology we propose in this paper is a nontrivial generalization of the parametrized maximum principle preserving flux limiters for high order finite difference schemes solving scalar hyperbolic conservation laws [22, 10, 20]. To preserve the maximum principle, the high order flux is limited towards a first order monotone flux, where the limiting procedures are designed by decoupling linear maximum principle constraints. High order schemes with such flux limiters are shown to preserve the high order accuracy via local truncation error analysis and by extensive numerical experiments with mild CFL constraints. The parametrized flux limiting approach is generalized to the Euler system to preserve the positivity of density and pressure of numerical solutions via decoupling some nonlinear constraints. Compared with existing high order positivity preserving approaches [24, 26, 25], our proposed algorithm is positivity preserving by the design; it is computationally efficient and maintains high order spatial and temporal accuracy in our extensive numerical tests. Numerical tests are performed to demonstrate the efficiency and effectiveness of the proposed new algorithm.

preprint2014arXiv

Runge-Kutta Central Discontinuous Galerkin BGK Method for the Navier-Stokes Equations

In this paper, we propose a Runge-Kutta (RK) central discontinuous Galerkin (CDG) gas-kinetic BGK method for the Navier-Stokes equations. The proposed method is based on the CDG method defined on two sets of overlapping meshes to avoid discontinuous solutions at cell interfaces, as well as the gas-kinetic BGK model to evaluate fluxes for both convection and diffusion terms. Redundant representation of the numerical solution in the CDG method offers great convenience in the design of gas-kinetic BGK fluxes. Specifically, the evaluation of fluxes at cell interfaces of one set of computational mesh is right inside the cells of the staggered mesh, hence the corresponding particle distribution function for flux evaluation is much simpler than that in existing gas-kinetic BGK methods. As a central scheme, the proposed CDG-BGK has doubled the memory requirement as the corresponding DG scheme; on the other hand, {for the convection part,} the CFL time step constraint of the CDG method for numerical stability is relatively large compared with that for the DG method. Numerical boundary conditions have to be treated with special care. Numerical examples for 1D and 2D viscous flow simulations are presented to validate the accuracy and robustness of the proposed RK CDG-BGK method.

preprint2013arXiv

High Order Asymptotic Preserving DG-IMEX Schemes for Discrete-Velocity Kinetic Equations in a Diffusive Scaling

In this paper, we develop a family of high order asymptotic preserving schemes for some discrete-velocity kinetic equations under a diffusive scaling, that in the asymptotic limit lead to macroscopic models such as the heat equation, the porous media equation, the advection-diffusion equation, and the viscous Burgers equation. Our approach is based on the micro-macro reformulation of the kinetic equation which involves a natural decomposition of the equation to the equilibrium and non-equilibrium parts. To achieve high order accuracy and uniform stability as well as to capture the correct asymptotic limit, two new ingredients are employed in the proposed methods: discontinuous Galerkin spatial discretization of arbitrary order of accuracy with suitable numerical fluxes; high order globally stiffly accurate implicit-explicit Runge-Kutta scheme in time equipped with a properly chosen implicit-explicit strategy. Formal asymptotic analysis shows that the proposed scheme in the limit of epsilon -> 0 is an explicit, consistent and high order discretization for the limiting equation. Numerical results are presented to demonstrate the stability and high order accuracy of the proposed schemes together with their performance in the limit.

preprint2013arXiv

PSO-MISMO Modeling Strategy for Multi-Step-Ahead Time Series Prediction

Multi-step-ahead time series prediction is one of the most challenging research topics in the field of time series modeling and prediction, and is continually under research. Recently, the multiple-input several multiple-outputs (MISMO) modeling strategy has been proposed as a promising alternative for multi-step-ahead time series prediction, exhibiting advantages compared with the two currently dominating strategies, the iterated and the direct strategies. Built on the established MISMO strategy, this study proposes a particle swarm optimization (PSO)-based MISMO modeling strategy, which is capable of determining the number of sub-models in a self-adaptive mode, with varying prediction horizons. Rather than deriving crisp divides with equal-size s prediction horizons from the established MISMO, the proposed PSO-MISMO strategy, implemented with neural networks, employs a heuristic to create flexible divides with varying sizes of prediction horizons and to generate corresponding sub-models, providing considerable flexibility in model construction, which has been validated with simulated and real datasets.