Search

We develop in this paper a new regularized flow dynamic approach to construct efficient numerical schemes for Wasserstein gradient flows in Lagrangian coordinates. Instead of approximating the Wasserstein distance which needs to solve constrained minimization problems, we reformulate the problem using the Benamou-Brenier's flow dynamic approach, leading to algorithms which only need to solve unconstrained minimization problem in $L^2$ distance. Our schemes automatically inherit some essential properties of Wasserstein gradient systems such as positivity-preserving, mass conservative and energy dissipation. We present ample numerical simulations of Porous-Medium equations, Keller-Segel equations and Aggregation equations to validate the accuracy and stability of the proposed schemes. Compared to numerical schemes in Eulerian coordinates, our new schemes can capture sharp interfaces for various Wasserstein gradient flows using relatively smaller number of unknowns.

We present an optimal rate convergence analysis for a second order accurate in time, fully discrete finite difference scheme for the Cahn-Hilliard-Navier-Stokes (CHNS) system, combined with logarithmic Flory-Huggins energy potential. The numerical scheme has been recently proposed, and the positivity-preserving property of the logarithmic arguments, as well as the total energy stability, have been theoretically justified. In this paper, we rigorously prove second order convergence of the proposed numerical scheme, in both time and space. Since the CHNS is a coupled system, the standard $\ell^\infty (0, T; \ell^2) \cap \ell^2 (0, T; H_h^2)$ error estimate could not be easily derived, due to the lack of regularity to control the numerical error associated with the coupled terms. Instead, the $\ell^\infty (0, T; H_h^1) \cap \ell^2 (0, T; H_h^3)$ error analysis for the phase variable and the $\ell^\infty (0, T; \ell^2)$ analysis for the velocity vector, which shares the same regularity as the energy estimate, is more suitable to pass through the nonlinear analysis for the error terms associated with the coupled physical process. Furthermore, the highly nonlinear and singular nature of the logarithmic error terms makes the convergence analysis even more challenging, since a uniform distance between the numerical solution and the singular limit values of is needed for the associated error estimate. Many highly non-standard estimates, such as a higher order asymptotic expansion of the numerical solution (up to the third order accuracy in time and fourth order in space), combined with a rough error estimate (to establish the maximum norm bound for the phase variable), as well as a refined error estimate, have to be carried out to conclude the desired convergence result.

This paper is concerned with the extraction of the smallest eigenvalue and the corresponding eigenvector of a symmetric positive definite matrix pencil. We reveal implicit convexity of the eigenvalue problem in Euclidean space. A provable accelerated eigensolver based on preconditioning and implicit convexity (EPIC) is proposed. Theoretical analysis shows the acceleration of EPIC with the rate of convergence resembling the expected rate of convergence of the well-known locally optimal preconditioned conjugate gradient (LOPCG). A complete proof of the expected rate of convergence of LOPCG is elusive so far. Numerical results confirm our theoretical findings of EPIC.

The analysis of the acceleration behavior of gradient-based eigensolvers with preconditioning presents a substantial theoretical challenge. In this work, we present a novel framework for preconditioning on Riemannian manifolds and introduce a metric, the leading angle, to evaluate preconditioners for symmetric eigenvalue problems. We extend the locally optimal Riemannian accelerated gradient method for Riemannian convex optimization to develop the Riemannian Acceleration with Preconditioning (RAP) method for symmetric eigenvalue problems, thereby providing theoretical evidence to support its acceleration. Our analysis of the Schwarz preconditioner for elliptic eigenvalue problems demonstrates that RAP achieves a convergence rate of $1-C\kappa^{-1/2}$, which is an improvement over the preconditioned steepest descent method's rate of $1-C\kappa^{-1}$. The exponent in $\kappa^{-1/2}$ is sharp, and numerical experiments confirm our theoretical findings.

In high-speed rotating channels, significant compressive effects are observed, resulting in distinct flow characteristics compared to incompressible flows. In this study, we employed a finite volume method based on the simple algorithm to solve for low-speed compressible laminar flow within rotating channels using an orthogonal uniform grid. The governing equations include the full Navier-Stokes equations and the energy equation. Contrary to stationary channel, the alterations in flow within rotating channel are primarily influenced by the compressive effects of centrifugal force and the compressibility of fluid within the flow's normal section. The first effect involves a reduction in the velocity due to centrifugal force, leading to an increasing influence of the Coriolis force compared to inertial forces along the flow direction. This trend in axial changes aligns closely with the increase in rotation speed. The second effect arises from the increase in Mach number and the Coriolis compression, resulting in slight density differences within the cross-section. Strong centrifugal forces generate significant centrifugal additional force (buoyancy force). Consequently, under the same local rotation number, the velocity profiles of the mainstream experience considerable changes. Additionally, higher Mach number significantly impact wall shear stress, with the leading side being notably affected. For instance, at a cross-sectional Ro = 0.6 and Ma = 0.035, the dimensionless shear stress on the leading side decreased by 13%. Furthermore, while an increase in Mach number has minimal impact on the cross-sectional secondary flow structure, changes in mainstream velocity profiles influence secondary flow intensity, resulting in an enhanced velocity peak and a shift towards the trailing side.

In the study of rotating channel flow, the key dimensionless parameters typically include the Reynolds number, rotation number, Prandtl number and buoyancy number. Our research focused on comparing the flow characteristics between the enlarged model, analyzed under the rotating similarity theory, and the original channel flow. Significantly different flow behaviors were observed between these two cases. Through theoretical derivation and dimensional analysis, we identified a new significant parameter - the centrifugal work number (CW). This parameter characterizes the ratio of centrifugal work to gas enthalpy in the rotating channel and plays a crucial role in measuring the compressibility of fluids within the rotating channel. Additionally, we utilized large eddy simulation(LES) to validate the impact of the centrifugal work ratio on the flow state of the rotating channel, thus enhancing the similarity theory of rotating channel compressible flow., Comment: 14 pages, 8 figures

The Keller-Segel equations are widely used for describing chemotaxis in biology. Recently, a new fully discrete scheme for this model was proposed in [46], mass conservation, positivity and energy decay were proved for the proposed scheme, which are important properties of the original system. In this paper, we establish the error estimates of this scheme. Then, based on the error estimates, we derive the finite-time blowup of nonradial numerical solutions under some conditions on the mass and the moment of the initial data.

We propose and analyze a first-order finite difference scheme for the functionalized Cahn-Hilliard (FCH) equation with a logarithmic Flory-Huggins potential. The semi-implicit numerical scheme is designed based on a suitable convex-concave decomposition of the FCH free energy. We prove unique solvability of the numerical algorithm and verify its unconditional energy stability without any restriction on the time step size. Thanks to the singular nature of the logarithmic part in the Flory-Huggins potential near the pure states $\pm 1$, we establish the so-called positivity-preserving property for the phase function at a theoretic level. As a consequence, the numerical solutions will never reach the singular values $\pm 1$ in the point-wise sense and the fully discrete scheme is well defined at each time step. Next, we present a detailed optimal rate convergence analysis and derive error estimates in $l^{\infty}(0,T;L_h^2)\cap l^2(0,T;H^3_h)$ under a linear refinement requirement $\Delta t\leq C_1 h$. To achieve the goal, a higher order asymptotic expansion (up to the second order temporal and spatial accuracy) based on the Fourier projection is utilized to control the discrete maximum norm of solutions to the numerical scheme. We show that if the exact solution to the continuous problem is strictly separated from the pure states $\pm 1$, then the numerical solutions can be kept away from $\pm 1$ by a positive distance that is uniform with respect to the size of the time step and the grid. Finally, a few numerical experiments are presented. Convergence test is performed to demonstrate the accuracy and robustness of the proposed numerical scheme. Pearling bifurcation, meandering instability and spinodal decomposition are observed in the numerical simulations., Comment: 46 pages, 11 figures. The linear refinement condition on the time step for the optimal rate convergence analysis was improved. Typos were fixed

Failure management plays a significant role in optical networks. It ensures secure operation, mitigates potential risks, and executes proactive protection. Machine learning (ML) is considered to be an extremely powerful technique for performing comprehensive data analysis and complex network management and is widely utilized for failure management in optical networks to revolutionize the conventional manual methods. In this study, the background of failure management is introduced, where typical failure tasks, physical objects, ML algorithms, data source, and extracted information are illustrated in detail. An overview of the applications of ML in failure management is provided in terms of alarm analysis, failure prediction, failure detection, failure localization, and failure identification. Finally, the future directions on ML for failure management are discussed from the perspective of data, model, task, and emerging techniques.

In this paper, we consider the Newton-Schur method in Hilbert space and obtain quadratic convergence. For the symmetric elliptic eigenvalue problem discretized by the standard finite element method and non-overlapping domain decomposition method, we use the Steklov-Poincar\'e operator to reduce the eigenvalue problem on the domain $\Omega$ into the nonlinear eigenvalue subproblem on $\Gamma$, which is the union of subdomain boundaries. We prove that the convergence rate for the Newton-Schur method is $\epsilon_{N}\leq CH^{2}(1+\ln(H/h))^{2}\epsilon^{2}$, where the constant $C$ is independent of the fine mesh size $h$ and coarse mesh size $H$, and $\epsilon_{N}$ and $\epsilon$ are errors after and before one iteration step respectively. Numerical experiments confirm our theoretical analysis.

This article concerns the weak Galerkin mixed finite element method (WG-MFEM) for second order elliptic equations on 2D domains with curved boundary. The Neumann boundary condition is considered since it becomes the essential boundary condition in this case. It is well-known that the discrepancy between the curved physical domain and the polygonal approximation domain leads to a loss of accuracy for discretization with polynomial order $\alpha>1$. The purpose of this paper is two-fold. First, we present a detailed error analysis of the original WG-MFEM for solving problems on curved domains, which exhibits an $O(h^{1/2})$ convergence for all $\alpha\ge 1$. It is a little surprising to see that even the lowest-order WG-MFEM ($\alpha=1$) experiences a loss of accuracy. This is different from known results for the finite element method (FEM) or the mixed FEM, and appears to be a combined effect of the WG-MFEM design and the fact that the outward normal vector on the polygonal approximation domain is different from the one on the curved domain. Second, we propose a remedy to bring the approximation rate back to optimal by employing two techniques. One is a specially designed boundary correction technique. The other is to take full advantage of the nice feature that weak Galerkin discretization can be defined on polygonal meshes, which allows the curved boundary to be better approximated by multiple short edges without increasing the total number of mesh elements. Rigorous analysis shows that a combination of the above two techniques renders optimal convergence for all $\alpha$. Numerical results further confirm this conclusion.

We analyze a fully discrete finite element numerical scheme for the Cahn-Hilliard-Stokes-Darcy system that models two-phase flows in coupled free flow and porous media. To avoid a well-known difficulty associated with the coupling between the Cahn-Hilliard equation and the fluid motion, we make use of the operator-splitting in the numerical scheme, so that these two solvers are decoupled, which in turn would greatly improve the computational efficiency. The unique solvability and the energy stability have been proved in~\cite{CHW2017}. In this work, we carry out a detailed convergence analysis and error estimate for the fully discrete finite element scheme, so that the optimal rate convergence order is established in the energy norm, i.e.,, in the $\ell^\infty (0, T; H^1) \cap \ell^2 (0, T; H^2)$ norm for the phase variables, as well as in the $\ell^\infty (0, T; H^1) \cap \ell^2 (0, T; H^2)$ norm for the velocity variable. Such an energy norm error estimate leads to a cancellation of a nonlinear error term associated with the convection part, which turns out to be a key step to pass through the analysis. In addition, a discrete $\ell^2 (0;T; H^3)$ bound of the numerical solution for the phase variables plays an important role in the error estimate, which is accomplished via a discrete version of Gagliardo-Nirenberg inequality in the finite element setting.