From Time-Collocated to Leapfrog Fundamental Schemes for ADI and CDI FDTD Methods

The leapfrog schemes have been developed for unconditionally stable alternating-direction implicit (ADI) finite-difference time-domain (FDTD) method, and recently the complying-divergence implicit (CDI) FDTD method. In this paper, the formulations from time-collocated to leapfrog fundamental schemes are presented for ADI and CDI FDTD methods. For the ADI FDTD method, the time-collocated fundamental schemes are implemented using implicit E-E and E-H update procedures, which comprise simple and concise right-hand sides (RHS) in their update equations. From the fundamental implicit E-H scheme, the leapfrog ADI FDTD method is formulated in conventional form, whose RHS are simplified into the leapfrog fundamental scheme with reduced operations and improved efficiency. For the CDI FDTD method, the time-collocated fundamental scheme is presented based on locally one-dimensional (LOD) FDTD method with complying divergence. The formulations from time-collocated to leapfrog schemes are provided, which result in the leapfrog fundamental scheme for CDI FDTD method. Based on their fundamental forms, further insights are given into the relations of leapfrog fundamental schemes for ADI and CDI FDTD methods. The time-collocated fundamental schemes require considerably fewer operations than all conventional ADI, LOD and leapfrog ADI FDTD methods, while the leapfrog fundamental schemes for ADI and CDI FDTD methods constitute the most efficient implicit FDTD schemes to date.

Stability properties of a projector-splitting scheme for dynamical low rank approximation of random parabolic equations

We consider the Dynamical Low Rank (DLR) approximation of random parabolic equations and propose a class of fully discrete numerical schemes. Similarly to the continuous DLR approximation, our schemes are shown to satisfy a discrete variational formulation. By exploiting this property, we establish stability of our schemes: we show that our explicit and semi-implicit versions are conditionally stable under a “parabolic” type CFL condition which does not depend on the smallest singular value of the DLR solution; whereas our implicit scheme is unconditionally stable. Moreover, we show that, in certain cases, the semi-implicit scheme can be unconditionally stable if the randomness in the system is sufficiently small. Furthermore, we show that these schemes can be interpreted as projector-splitting integrators and are strongly related to the scheme proposed in [29, 30], to which our stability analysis applies as well. The analysis is supported by numerical results showing the sharpness of the obtained stability conditions.

Theoretical Investigation of Unconditionally Stable Sequential Algorithms of Coupled Flow and Geomechanics for Hydraulically Fractured Systems

We investigate unconditionally stable sequential algorithms for coupled hydraulically fractured geomechanics and flow systems, which can account for poromechanics behavior within the fractures. We focus on modifying the concepts of the fixed stress and undrained sequential methods properly for the coupled systems by taking appropriate stabilization terms for stability and convergence with energy analyses. Specifically, an apparent fracture stiffness is used for for numerical stabilization. Because this fracture stiffness depends on the fracture length, the stabilization term needs to be updated dynamically, different from the drained bulk modulus used for typical poromechanics problems. For numerical tests, we take the extended finite element method for geomechanics while the piecewise constant finite element method is used for flow within an existing hydraulic fracture. The numerical results support a priori stability analyses.

Flume Tests to Investigate the Initiation Mechanism of Loess Mudflows on the Chinese Loess Plateau

Loess has a strong water sensitivity, so loess landslides often transform into loess mudflows when water is added on the Chinese Loess Plateau, which results in high casualties and property loss of the Chinese government. In this study, a series of flume tests were designed to study the initiation of loess mudflows. The results reveal that the initiation modes of loess mudflows include large-scale mudflow and retrogressive toe sliding (Type A), and small-scale mudflow and retrogressive toe sliding (Type B). A model was used to analyze the test results that describe the effects of water flow on the potential for hillslope failure and liquefaction. It was found that the soil accumulation was unconditionally stable before a loess mudflow was formed, but as the rainfall continued, the water gradually infiltrated the soil, and the soil accumulation changed from unconditionally stable to unconditionally unstable. Thus, this led to different initiation modes during the tests. For Type A, the water preferentially infiltrated into the area with an uneven density and a large amount of water accumulated. The pore water pressure increased quickly and could not dissipate in time, so the loess liquefied. As the liquefaction area continued to expand and became larger, Type A occurred. Relatively speaking, Type B occurs in soil accumulations with relatively uniform densities. These results provide a certain scientific reference for the study of loess mudflows.

Numerical Solution of the coupled Burgers' equation by Trigonometric B-spline Collocation Method

In the present study, the coupled Burgers’ equation is going to be solved numerically by presenting a new technique based on collocation finite element method in which trigonometric cubic and quintic B-splines are used as approximate functions. In order to support the present study, three test problems given with appropriate initial and boundary conditions are studied. The newly obtained results are compared with some of the other published numerical solutions available in the literature. The accuracy of the proposed method is discussed by computing the error norms L₂ and L_{∞}. A linear stability analysis of the approximation obtained by the scheme shows that the method is unconditionally stable.