Search

Due to quasicrystals having long-range orientational order but without translational symmetry, traditional numerical methods usually suffer when applied as is. In the past decade, the projection method has emerged as a prominent solver for quasiperiodic problems. Transforming them into a higher dimensional but periodic ones, the projection method facilitates the application of the fast Fourier transform. However, the computational complexity inevitably becomes high which significantly impedes e.g. the generation of the phase diagram since a high-fidelity simulation of a problem whose dimension is doubled must be performed for numerous times. To address the computational challenge of quasiperiodic problems based on the projection method, this paper proposes a multi-component multi-state reduced basis method (MCMS-RBM). Featuring multiple components with each providing reduction functionality for one branch of the problem induced by one part of the parameter domain, the MCMS-RBM does not resort to the parameter domain configurations (e.g. phase diagrams) a priori. It enriches each component in a greedy fashion via a phase-transition guided exploration of the multiple states inherent to the problem. Adopting the empirical interpolation method, the resulting online-efficient method vastly accelerates the generation of a delicate phase diagram to a matter of minutes for a parametrized two-turn-four dimensional Lifshitz-Petrich model with two length scales. Moreover, it furnishes surrogate and equally accurate field variables anywhere in the parameter domain., Comment: 18 pages, 7 figures

Supercapacitors are promising electrochemical energy storage devices due to their prominent performance in rapid charging/discharging rates, long cycle life, stability, etc. Experimental measurement and theoretical prediction on charging timescale for supercapacitors often have large difference. This work develops a matched asymptotic expansion method to derive the charging dynamics of supercapacitors with porous electrodes, in which the supercapacitors are described by the stack-electrode model. Coupling leading-order solutions between every two stacks by continuity of ionic concentration and fluxes leads to an ODE system, which is a generalized equivalent circuit model for zeta potentials, with the potential-dependent nonlinear capacitance and resistance determined by physical parameters of electrolytes, e.g., specific counterion valences for asymmetric electrolytes. Linearized stability analysis on the ODE system after projection is developed to theoretically characterize the charging timescale. The derived asymptotic solutions are numerically verified. Further numerical investigations on the biexponential charging timescales demonstrate that the proposed generalized equivalent circuit model, as well as companion linearized stability analysis, can faithfully capture the charging dynamics of symmetric/asymmetric electrolytes in supercapacitors with porous electrodes., Comment: 26 pages, 6 figures

In this paper, we develop an asymptotic-preserving and energy-conserving (APEC) Particle-In-Cell (PIC) algorithm for the Vlasov-Maxwell system. This algorithm not only guarantees that the asymptotic limiting of the discrete scheme is a consistent and stable discretization of the quasi-neutral limit of the continuous model, but also preserves Gauss's law and energy conservation at the same time, thus it is promising to provide stable simulations of complex plasma systems even in the quasi-neutral regime. The key ingredients for achieving these properties include the generalized Ohm's law for electric field such that the asymptotic-preserving discretization can be achieved, and a proper decomposition of the effects of the electromagnetic fields such that a Lagrange multiplier method can be appropriately employed for correcting the kinetic energy. We investigate the performance of the APEC method with three benchmark tests in one dimension, including the linear Landau damping, the bump-on-tail problem and the two-stream instability. Detailed comparisons are conducted by including the results from the classical explicit leapfrog and the previously developed asymptotic-preserving PIC schemes. Our numerical experiments show that the proposed APEC scheme can give accurate and stable simulations both kinetic and quasi-neutral regimes, demonstrating the attractive properties of the method crossing scales., Comment: 21 pages, 30 figures

The need for accelerating the repeated solving of certain parametrized systems motivates the development of more efficient reduced order methods. The classical reduced basis method is popular due to an offline-online decomposition and a mathematically rigorous {\em a posterior} error estimator which guides a greedy algorithm offline. For nonlinear and nonaffine problems, hyper reduction techniques have been introduced to make this decomposition efficient. However, they may be tricky to implement and often degrade the online computation efficiency. To avoid this degradation, reduced residual reduced over-collocation (R2-ROC) was invented integrating empirical interpolation techniques on the solution snapshots and well-chosen residuals, the collocation philosophy, and the simplicity of evaluating the hyper-reduced well-chosen residuals. In this paper, we introduce an adaptive enrichment strategy for R2-ROC rendering it capable of handling parametric fluid flow problems. Built on top of an underlying Marker and Cell (MAC) scheme, a novel hyper-reduced MAC scheme is therefore presented and tested on Stokes and Navier-Stokes equations demonstrating its high efficiency, stability and accuracy.

The need for multiple interactive, real-time simulations using different parameter values has driven the design of fast numerical algorithms with certifiable accuracies. The reduced basis method (RBM) presents itself as such an option. RBM features a mathematically rigorous error estimator which drives the construction of a low-dimensional subspace. A surrogate solution is then sought in this low-dimensional space approximating the parameter-induced high fidelity solution manifold. However when the system is nonlinear or its parameter dependence nonaffine, this efficiency gain degrades tremendously, an inherent drawback of the application of the empirical interpolation method (EIM). In this paper, we augment and extend the EIM approach as a direct solver, as opposed to an assistant, for solving nonlinear partial differential equations on the reduced level. The resulting method, called Reduced Over-Collocation method (ROC), is stable and capable of avoiding the efficiency degradation. Two critical ingredients of the scheme are collocation at about twice as many locations as the number of basis elements for the reduced approximation space, and an efficient error indicator for the strategic building of the reduced solution space. The latter, the main contribution of this paper, results from an adaptive hyper reduction of the residuals for the reduced solution. Together, these two ingredients render the proposed R2-ROC scheme both offline- and online-efficient. A distinctive feature is that the efficiency degradation appearing in traditional RBM approaches that utilize EIM for nonlinear and nonaffine problems is circumvented, both in the offline and online stages. Numerical tests on different families of time-dependent and steady-state nonlinear problems demonstrate the high efficiency and accuracy of our R2-ROC and its superior stability performance., Comment: arXiv admin note: substantial text overlap with arXiv:2009.04812, arXiv:1906.07349

The task of repeatedly solving parametrized partial differential equations (pPDEs) in, e.g. optimization or interactive applications, makes it imperative to design highly efficient and equally accurate surrogate models. The reduced basis method (RBM) presents as such an option. Enabled by a mathematically rigorous error estimator, RBM constructs a low-dimensional subspace of the parameter-induced high fidelity solution manifold from which an approximate solution is computed. It can improve efficiency by several orders of magnitudes leveraging an offline-online decomposition procedure. However, this decomposition, usually through the empirical interpolation method (EIM) when the PDE is nonlinear or its parameter dependence nonaffine, is either challenging to implement, or severely degrades online efficiency. In this paper, we augment and extend the EIM approach as a direct solver, as opposed to an assistant, for solving nonlinear pPDEs on the reduced level. The resulting method, called Reduced Over-Collocation method (ROC), is stable and capable of avoiding the efficiency degradation inherent to a traditional application of EIM. Two critical ingredients of the scheme are collocation at about twice as many locations as the dimension of the reduced solution space, and an efficient L1-norm-based error indicator for the strategic selection of the parameter values to build the reduced solution space. Together, these two ingredients render the proposed L1-ROC scheme both offline- and online-efficient. A distinctive feature is that the efficiency degradation appearing in alternative RBM approaches that utilize EIM for nonlinear and nonaffine problems is circumvented, both in the offline and online stages. Numerical tests on different families of time-dependent and steady-state nonlinear problems demonstrate the high efficiency and accuracy of L1-ROC and its superior stability performance., Comment: 24 pages, 5 figures. arXiv admin note: substantial text overlap with arXiv:1906.07349

The onerous task of repeatedly resolving certain parametrized partial differential equations (pPDEs) in, e.g. the optimization context, makes it imperative to design vastly more efficient numerical solvers without sacrificing any accuracy. The reduced basis method (RBM) presents itself as such an option. With a mathematically rigorous error estimator, RBM seeks a surrogate solution in a carefully-built subspace of the parameter-induced high fidelity solution manifold. It can improve efficiency by several orders of magnitudes leveraging an offline-online decomposition procedure. However, this decomposition, usually through the empirical interpolation method (EIM) when the PDE is nonlinear or its parameter dependence nonaffine, is either challenging to implement, or severely degrading to the online efficiency. In this paper, we augment and extend the EIM approach in the context of solving pPDEs in two different ways, resulting in the Reduced Over-Collocation methods (ROC). These are stable and capable of avoiding the efficiency degradation inherent to a direct application of EIM. There are two ingredients of these methods. First is a strategy to collocate at about twice as many locations as the number of bases for the surrogate space. The second is an efficient approach for the strategic selection of the parameter values to build the reduced solution space for which we study two choices, a recent empirical L1 approach and a new indicator based on the reduced residual. Together, these two ingredients render the schemes, L1-ROC and R2-ROC, online efficient and immune from the efficiency degradation of EIM for nonlinear and nonaffine problems offline and online. Numerical tests on three different families of nonlinear problems demonstrate the high efficiency and accuracy of these new algorithms and their superior stability performance.

In numerical simulations of many charged systems at the micro/nano scale, a common theme is the repeated solution of the Poisson-Boltzmann equation. This task proves challenging, if not entirely infeasible, largely due to the nonlinearity of the equation and the high dimensionality of the physical and parametric domains with the latter emulating the system configuration. In this paper, we for the first time adapt a mathematically rigorous and computationally efficient model order reduction paradigm, the so-called reduced basis method (RBM), to mitigate this challenge. We adopt a finite difference method as the mandatory underlying scheme to produce the {\em truth approximations} of the RBM upon which the fast algorithm is built and its performance is measured against. Numerical tests presented in this paper demonstrate the high efficiency and accuracy of the fast algorithm, the reliability of its error estimation, as well as its capability in effectively capturing the boundary layer.