Probabilistic Analysis of Slope Using Finite Element Approach and Limit Equilibrium Approach around Amalpata Landslide of West Central, Nepal

The stability study of the ongoing and recurring Amalpata landslide in Baglung in Nepal’s Gandaki Province is presented in this research. The impacted slope is around 200 meters high, with two terraces that have different slope inclinations. The lower bench, located above the basement, consistently fails and sets others up for failure. The fluctuating water level of the slope, which travels down the slope masses, exacerbates the slide problem. The majority of these rocks are Amalpata landslide area experiences several structural disruptions. The area’s stability must be evaluated in order to prevent and control more harm from occurring to the nearby agricultural land and people living along the slope. The slopes’ failures increase the damages of house existing in nearby area and the erosion of the slope. Two modeling techniques the finite element approach and the limit equilibrium method were used to simulate the slope. The findings show that, in every case, the terrace above the basement is where the majority of the stress is concentrated, with a safety factor of near unity. Using probabilistic slope stability analysis, the failure probability was predicted to be between 98.90% and 100%.

A Full Predictor-Corrector Finite Element Method for the One-Dimensional Heat Equation with Time-Dependent Singularities

The energy norm convergence rate of the finite element solution of the heat equation is reduced by the time-regularity of the exact solution. This paper presents an adaptive finite element treatment of time-dependent singularities on the one-dimensional heat equation. The method is based on a Fourier decomposition of the solution and an extraction formula of the coefficients of the singularities coupled with a predictor-corrector algorithm. The method recovers the optimal convergence rate of the finite element method on a quasi-uniform mesh refinement. Numerical results are carried out to show the efficiency of the method.

PERTURBATIONAL SOLUTIONS FOR FUZZY-STOCHASTIC FINITE ELEMENT EQUILIBRIUM EQUATIONS (FSFEEE)

In this paper, the random interval equilibrium equations (RIEE) is obtained by lambda-level cutting the fuzzy-stochastic finite element equilibrium equations (FSFEEE). Based on the relations between the variables of equilibrium equations, solving RIEE is transformed into solving two kinds of general random equilibrium equations (GREE). Then the recursive equations of evaluating the random interval displacement is derived from the small-parameter perturbation theory. The computational formulae of statistical characteristic of the fuzzy random displacements, the fuzzy random strains and the fuzzy random stresses are also deduced in detail.

A NEW SOLUTION TO STRUCTURAL FUZZY FINITE ELEMENT EQUILIBRIUM EQUATIONS (SFFEEE)

A speedy accurate solution to structural fuzzy finite element equilibrium equations (SFFEEE), by combining the definition of the solution of interval equations with the mechanical meaning of the structural finite element equilibrium equations (SFEEE), was put forward. The fuzzification of the SFFEEE, which is discussed in this paper, originates from that of material property, structural boundary conditions and external loading. The computing quantity of this solution is almost equal to that of the general finite element method (GFEM).

The Finite Element Solutions to the Semiconductor Equations

In this paper, the approximation of stationary equations of the semiconductor devices with mixed boundary conditions is considered. Two schemes are proposed for the system. One is Glerkin discrete scheme, the other is hybrid variable discrete scheme. A convergence analysis is also given.

Error estimates of H^1-Galerkin mixed finite element method for Schrdinger equation

An H^1-Galerkin mixed finite element method is discussed for a class of second order SchrSdinger equation. Optimal error estimates of semidiscrete schemes are derived for problems in one space dimension. At the same time, optimal error estimates are derived for fully discrete schemes. And it is showed that the H1-Galerkin mixed finite element approximations have the same rate of convergence as in the classical mixed finite element methods without requiring the LBB consistency condition.

A lumped mass nonconforming finite element method for nonlinear parabolic integro-differential equations on anisotropic meshes

A lumped mass approximation scheme of a low order Crouzeix-Raviart type noncon- forming triangular finite element is proposed to a kind of nonlinear parabolic integro-differential equations. The L2 error estimate is derived on anisotropic meshes without referring to the traditional nonclassical elliptic projection.

Highly efficient H^1-Galerkin mixed finite element method (MFEM) for parabolic integro-differential equation

A highly efficient H1-Galerkin mixed finite element method （MFEM） is presented with linear triangular element for the parabolic integro-differential equation. Firstly, some new results about the integral estimation and asymptotic expansions are studied. Then, the superconvergence of order O（h^2） for both the original variable u in H1 （Ω） norm and the flux p = u in H（div, Ω） norm is derived through the interpolation post processing technique. Furthermore, with the help of the asymptotic expansions and a suitable auxiliary problem, the extrapolation solutions with accuracy O（h^3） are obtained for the above two variables. Finally, some numerical results are provided to confirm validity of the theoretical analysis and excellent performance of the proposed method.

HIGH ACCURACY FINITE VOLUME ELEMENT METHOD FOR TWO-POINT BOUNDARY VALUE PROBLEM OF SECOND ORDER ORDINARY DIFFERENTIAL EQUATIONS

In this paper, a high accuracy finite volume element method is presented for two-point boundary value problem of second order ordinary differential equation, which differs from the high order generalized difference methods. It is proved that the method has optimal order error estimate O(h3) in H1 norm. Finally, two examples show that the method is effective.

Error Estimates for Mixed Finite Element Methods for Sobolev Equation

The purpose of this paper is to investigate the convergence of the mixed finite element method for the initial-boundary value problem for the Sobolev equation Ut-div{aut + b1 u} = f based on the Raviart-Thomas space Vh × Wh H(div; × L2(). Optimal order estimates are obtained for the approximation of u, ut, the associated velocity p and divp respectively in L(0,T;L2()), L(0,T;L2()), L(0,T;L2()2), and L2(0, T; L2()). Quasi-optimal order estimates are obtained for the approximations of u, ut in L(0, T; L()) and p in L(0,T; L()2).

Local and parallel finite element algorithms based on two-grid discretization for steady Navier-Stokes equations

Local and parallel finite element algorithms based on two-grid discretization for Navier-Stokes equations in two dimension are presented. Its basis is a coarse finite element space on the global domain and a fine finite element space on the subdomain. The local algorithm consists of finding a solution for a given nonlinear problem in the coarse finite element space and a solution for a linear problem in the fine finite element space, then droping the coarse solution of the region near the boundary. By overlapping domain decomposition, the parallel algorithms are obtained. This paper analyzes the error of these algorithms and gets some error estimates which are better than those of the standard finite element method. The numerical experiments are given too. By analyzing and comparing these results, it is shown that these algorithms are correct and high efficient.

Mixed finite element for two-dimensional incompressible convective Brinkman-Forchheimer equations

In this work, the two-dimensional convective Brinkman-Forchheimer equa- tions are considered. The well-posedness for the variational problem and its mixed finite element approximation is established, and the error estimates based on the conforming approximation are obtained. For the computation, a one-step Newton (or semi-Newton) iteration algorithm initialized using a fixed-point iteration is proposed. Finally, numerical experiments using a Taylor-Hood mixed element built on a structured or unstructured triangular mesh are implemented. The numerical results obtained using the algorithm are compared with the analytic data, and are shown to be in very good agreement. Moreover, the lid-driven problem at Reynolds numbers of 100 and 400 is considered and analyzed.

An RKDG finite element method for the one-dimensional inviscid compressible gas dynamics equations in a Lagrangian coordinate

In this paper,Runge-Kutta Discontinuous Galerkin（RKDG） finite element method is presented to solve the onedimensional inviscid compressible gas dynamic equations in a Lagrangian coordinate.The equations are discretized by the DG method in space and the temporal discretization is accomplished by the total variation diminishing Runge-Kutta method.A limiter based on the characteristic field decomposition is applied to maintain stability and non-oscillatory property of the RKDG method.For multi-medium fluid simulation,the two cells adjacent to the interface are treated differently from other cells.At first,a linear Riemann solver is applied to calculate the numerical ？ux at the interface.Numerical examples show that there is some oscillation in the vicinity of the interface.Then a nonlinear Riemann solver based on the characteristic formulation of the equation and the discontinuity relations is adopted to calculate the numerical ？ux at the interface,which suppresses the oscillation successfully.Several single-medium and multi-medium fluid examples are given to demonstrate the reliability and efficiency of the algorithm.

Constitutive equations and finite element implementation of strain localization in sand deformation

A recently proposed model coupling with the solid-fluid of the saturated sand was utilized to study the deformation band. Based on the critical state plasticity model by Borja and Andrade, the hydraulic conductivity tensor was naturally treated as a function of the spatial discretization matrix about the displacement and the stress field, allowing a more realistic representation of the physical phenomenon. The fully Lagrangian form of the Darcy law was resolved by Piola algorithm, and then the flow law was gained, leading to the implementation of a modified model of the saturated sand. Then the criterion for the onset of localization was derived and utilized to detect instability. The constitutive model was implemented in a finite element program coded by FORTRAN, which was used to predict the formation and development of shear bands in plane strain compression of saturated sand. At last, the formation mechanism of the shear band was discussed. It is shown that the model works well, and the simulation sample bifurcates at 1.18% axial strain, which is in a good qualitative agreement with the experiment. The pore pressure greatly affects the onset and development of the deformation band, and it obviously increases around the localization-prone regions with the direction toward the outer side of the normal of the shear band, while the pore stress flows nearly horizontally and is distributed equally far away the shear band region.

A Priori and A Posteriori Error Estimates of Streamline Diffusion Finite Element Method for Optimal Control Problem Governed by Convection Dominated Diffusion Equation

In this paper,we investigate a streamline diffusion finite element approxi- mation scheme for the constrained optimal control problem governed by linear con- vection dominated diffusion equations.We prove the existence and uniqueness of the discretized scheme.Then a priori and a posteriori error estimates are derived for the state,the co-state and the control.Three numerical examples are presented to illustrate our theoretical results.

Unified analysis for stabilized methods of low-order mixed finite elements for stationary Navier-Stokes equations

A unified analysis is presented for the stabilized methods including the pres- sure projection method and the pressure gradient local projection method of conforming and nonconforming low-order mixed finite elements for the stationary Navier-Stokes equa- tions. The existence and uniqueness of the solution and the optimal error estimates are proved.

Low order nonconforming mixed finite element method for nonstationary incompressible Navier-Stokes equations

This paper studies a low order mixed finite element method （FEM） for nonstationary incompressible Navier-Stokes equations. The velocity and pressure are approximated by the nonconforming constrained Q1^4ot element and the piecewise constant, respectively. The superconvergent error estimates of the velocity in the broken H^1-norm and the pressure in the L^2-norm are obtained respectively when the exact solutions are reasonably smooth. A numerical experiment is carried out to confirm the theoretical results.

MIXED FINITE ELEMENT METHOD FOR SOBOLEV EQUATIONS AND ITS ALTERNATING-DIRECTION ITERATIVE SCHEME

In this paper, we study the mixed element method for Sobolev equations. A time-discretization procedure is presented and analysed and the optimal order error estimates are derived.For convenience in practical computation, an alternating-direction iterative scheme of the mixed fi-nite element method is formulated and its stability and converbence are proved for the linear prob-lem. A numerical example is provided at the end of this paper.

Mixed time discontinuous space-time finite element method for convection diffusion equations

A mixed time discontinuous space-time finite element scheme for secondorder convection diffusion problems is constructed and analyzed. Order of the equation is lowered by the mixed finite element method. The low order equation is discretized with a space-time finite element method, continuous in space but discontinuous in time. Stability, existence, uniqueness and convergence of the approximate solutions are proved. Numerical results are presented to illustrate efficiency of the proposed method.

Weak Galerkin Finite Element Method for the Unsteady Stokes Equation

The Weak Galerkin (WG) finite element method for the unsteady Stokes equations in the primary velocity-pressure formulation is introduced in this paper. Optimal-order error estimates are established for the corresponding numerical approximation in an H1 norm for the velocity, and L2 norm for both the velocity and the pressure by use of the Stokes projection.