Two-dimensional finite element modeling of seepage-erosion coupled process within unsaturated soil slopes

Fines migration along with rainfall infiltration is a possible cause of failures of slopes composed of loose deposits.To investigate the intrinsic mechanisms,a rigid mathematical model which can fully capture the multi-phasic and multiphysical process is necessary.In this research,the macro and micro physical phenomena of fines migration process within deposited soil slopes under rainfall infiltration were summarized.Based on the mixture theory,a seepage-erosion model for unsaturated erodible soils capable to capture these phenomena mathematically was built based on a rigid theoretical framework.The model was used to simulate a set of rainfall flume tests involving fines migration phenomena with the finite element method.Two distinct slope failure modes observed experimentally,which were induced by the soil erosion-deposition properties,can be well reproduced by our numerical model.The seepage-erosion coupled process during the rainfall infiltration,as well as the intrinsic mechanism responsible for the slope failures,was illustrated in detail based on the numerical results.It was shown that the fines migration process can affect the hydro-mechanical response within unsaturated slopes significantly,and therefore special attention should be paid to those soil slopes susceptible to internal erosion.

Two-dimensional finite element mesh generation algorithm for electromagnetic field calculation

Two-dimensional finite element mesh generation algorithm for electromagnetic field calculation is proposed in this paper to improve the efficiency and accuracy of electromagnetic calculation. An image boundary extraction algorithm is developed to map the image on the geometric domain. Identification algorithm for the location of nodes in polygon area is proposed to determine the state of the node. To promote the average quality of the mesh and the efficiency of mesh generation, a novel force-based mesh smoothing algorithm is proposed. One test case and a typical electromagnetic calculation are used to testify the effectiveness and efficiency of the proposed algorithm. The results demonstrate that the proposed algorithm can produce a high-quality mesh with less iteration.

Two-Dimensional Finite Element Method Analysis Effect of the Recombination Velocity at the Grain Boundaries on the Characteristics of a Polycrystalline Silicon Solar Cell

To take into account the variation of the recombination velocity at the grain boundaries, we present in this paper a new approach of characterization of the solar cells, based on the two dimensional finite element method. The results of this study on a bifacial polycrystalline silicon solar cell, modelled in the rectangular form, highlighting the effects of the boundary recombination velocity (Sgb) on the solar cell electrical parameters. The photogenerated excess carrier’s density, the photocurrent density;the phototovoltage and the current-voltage characteristics are analyzed, namely. A good agreement with the results given in the literature is observed.

A Deep Learning Approach to Shape Optimization Problems for Flexoelectric Materials Using the Isogeometric Finite Element Method

A new approach for flexoelectricmaterial shape optimization is proposed in this study.In this work,a proxymodel based on artificial neural network(ANN)is used to solve the parameter optimization and shape optimization problems.To improve the fitting ability of the neural network,we use the idea of pre-training to determine the structure of the neural network and combine different optimizers for training.The isogeometric analysis-finite element method(IGA-FEM)is used to discretize the flexural theoretical formulas and obtain samples,which helps ANN to build a proxy model from the model shape to the target value.The effectiveness of the proposed method is verified through two numerical examples of parameter optimization and one numerical example of shape optimization.

Finite Element Method Simulation of Wellbore Stability under Different Operating and Geomechanical Conditions

The variation of the principal stress of formations with the working and geo-mechanical conditions can trigger wellbore instabilities and adversely affect the well completion.A finite element model,based on the theory of poro-elasticity and the Mohr-Coulomb rock damage criterion,is used here to analyze such a risk.The changes in wellbore stability before and after reservoir acidification are simulated for different pressure differences.The results indicate that the risk of wellbore instability grows with an increase in the production-pressure difference regardless of whether acidification is completed or not;the same is true for the instability area.After acidizing,the changes in the main geomechanical parameters(i.e.,elastic modulus,Poisson’s ratio,and rock strength)cause the maximum wellbore instability coefficient to increase.

Modularized and Parametric Modeling Technology for Finite Element Simulations of Underground Engineering under Complicated Geological Conditions

The surrounding geological conditions and supporting structures of underground engineering are often updated during construction,and these updates require repeated numerical modeling.To improve the numerical modeling efficiency of underground engineering,a modularized and parametric modeling cloud server is developed by using Python codes.The basic framework of the cloud server is as follows:input the modeling parameters into the web platform,implement Rhino software and FLAC3D software to model and run simulations in the cloud server,and return the simulation results to the web platform.The modeling program can automatically generate instructions that can run the modeling process in Rhino based on the input modeling parameters.The main modules of the modeling program include modeling the 3D geological structures,the underground engineering structures,and the supporting structures as well as meshing the geometric models.In particular,various cross-sections of underground caverns are crafted as parametricmodules in themodeling program.Themodularized and parametric modeling program is used for a finite element simulation of the underground powerhouse of the Shuangjiangkou Hydropower Station.This complicatedmodel is rapidly generated for the simulation,and the simulation results are reasonable.Thus,this modularized and parametric modeling program is applicable for three-dimensional finite element simulations and analyses.

A stable implicit nodal integration-based particle finite element method(N-PFEM)for modelling saturated soil dynamics

In this study,we present a novel nodal integration-based particle finite element method(N-PFEM)designed for the dynamic analysis of saturated soils.Our approach incorporates the nodal integration technique into a generalised Hellinger-Reissner(HR)variational principle,creating an implicit PFEM formulation.To mitigate the volumetric locking issue in low-order elements,we employ a node-based strain smoothing technique.By discretising field variables at the centre of smoothing cells,we achieve nodal integration over cells,eliminating the need for sophisticated mapping operations after re-meshing in the PFEM.We express the discretised governing equations as a min-max optimisation problem,which is further reformulated as a standard second-order cone programming(SOCP)problem.Stresses,pore water pressure,and displacements are simultaneously determined using the advanced primal-dual interior point method.Consequently,our numerical model offers improved accuracy for stresses and pore water pressure compared to the displacement-based PFEM formulation.Numerical experiments demonstrate that the N-PFEM efficiently captures both transient and long-term hydro-mechanical behaviour of saturated soils with high accuracy,obviating the need for stabilisation or regularisation techniques commonly employed in other nodal integration-based PFEM approaches.This work holds significant implications for the development of robust and accurate numerical tools for studying saturated soil dynamics.

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%.

Collapse Behavior of Pipe-Framed Greenhouses with and without Reinforcement under Snow Loading:A 3-D Finite Element Analysis

The present paper first investigates the collapse behavior of a conventional pipe-framed greenhouse under snow loading based on a 3-D finite element analysis,in which both geometrical and material non-linearities are considered.Three snow load distribution patterns related to the wind-driven snow particle movement are used in the analysis.It is found that snow load distribution affects the deformation and collapse behavior of the pipe-framed greenhouse significantly.The results obtained in this study are consistent with the actual damage observed.Next,discussion is made of the effects of reinforcements by adding members to the basic frame on the strength of the whole structure,in which seven kinds of reinforcement methods are examined.A buckling analysis is also carried out.The results indicate that the most effective reinforcement method depends on the snow load distribution pattern.

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.

Contribution to the Full 3D Finite Element Modelling of a Hybrid Stepping Motor with and without Current in the Coils

The paper presents our contribution to the full 3D finite element modelling of a hybrid stepping motor using COMSOL Multiphysics software. This type of four-phase motor has a permanent magnet interposed between the two identical and coaxial half stators. The calculation of the field with or without current in the windings (respectively with or without permanent magnet) is done using a mixed formulation with strong coupling. In addition, the local high saturation of the ferromagnetic material and the radial and axial components of the magnetic flux are taken into account. The results obtained make it possible to clearly observe, as a function of the intensity of the bus current or the remanent induction, the saturation zones, the lines, the orientations and the magnetic flux densities. 3D finite element modelling provide more accurate numerical data on the magnetic field through multiphysics analysis. This analysis considers the actual operating conditions and leads to the design of an optimized machine structure, with or without current in the windings and/or permanent magnet.

TWO-DIMENSIONAL FINITE ELEMENT METHOD FOR DETERMINING NONLINEAR WAVE FORCES ON LARGE SUBMERGED BODIES

In this paper, we first develop the far field asymptotic solutions of the second-order scattering waves for the vertical plane problem taking the second-order Stokes waves as the incident waves. The asymptotic solutions satisfy the Laplace equation, the sea bed and free surface boundary conditions and are the out-going waves. Then the radiation conditions of the second-order mattering waves are derived by using the asymptotic solutions. By using the two-dimensinal finite clement method with the radiation conditions imposed on the ar- tificial boundaries, the computer program, known as 'NWF2', for determining nonlinear wave forces on large submerged bodies has been written. As a numerical example, nonlinear wave forces on a semi-circu- lar cylinder lying on the sea bed arc presented.

Two-Dimensional Large Deformation Finite Element Analysisfor the Pulling-up of Plate Anchor

Based on mesh regeneration and stress interpolation from an old mesh to a new one, a large deformation finite element model is developed for the study of the behaviour of circular plate anchors subjected to uplift loading. For the deterruination of the distributions of stress components across a clay foundation, the Recovery by Equilibrium in Patches is extended to plastic analyses. ABAQUS, a commercial finite element package, is customized and linked into our program so as to keep automatic and efficient running of large deformation calculation. The quality of stress interpolation is testified by evaluations of Tresca stress and nodal reaction forces. The complete pulling-up processes of plate anchors buried in homogeneous clay arc simulated, and typical pulling force-displacement responses of a deep anchor and a shallow anchor are compared. Different from the results of previous studies, large deformation analysis is of the capability of estimating the breakaway between the anchor bottom and soils. For deep anchors, the variation of mobilized uplift resistance with anchor settlement is composed of three stages, and the initial buried depths of anchors affect the separation embedment slightly. The uplift bearing capacity of deep anchors is usually higher than that of shallow anchors.

Finite Element Analysis to Two-Dimensional Nonlinear Sloshing Problems

A two-dimensional nonlinear sloshing problem is analyzed by means of the fully nonlinear theory and time domain second order theory of water waves. Liquid sloshing in a rectangular container Subjected to a horizontal excitation is simulated by the finite element method. Comparisons between the two theories are made based on their numerical results. It is found that good agreement is obtained for the case of small amplitude oscillation and obvious differences occur for large amplitude excitation. Even though, the second order solution can still exhibit typical nonlinear features of nonlinear wave and can be used instead of the fully nonlinear theory.

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.

A Numerical Algorithm Based on Quadratic Finite Element for Two-Dimensional Nonlinear Time Fractional Thermal Diffusion Model

In this article,a high-order scheme,which is formulated by combining the quadratic finite element method in space with a second-order time discrete scheme,is developed for looking for the numerical solution of a two-dimensional nonlinear time fractional thermal diffusion model.The time Caputo fractional derivative is approximated by using the L2-1formula,the first-order derivative and nonlinear term are discretized by some second-order approximation formulas,and the quadratic finite element is used to approximate the spatial direction.The error accuracy O(h3+
t2)is obtained,which is verified by the numerical results.

Adaptive Finite Element Modeling of Marine Controlled-Source Electromagnetic Fields in Two-Dimensional General Anisotropic Media

In this paper, we extend the scope of numerical simulations of marine controlled-source electromagnetic (CSEM) fields in a particular case of anisotropy (dipping anisotropy) to the general case of anisotropy by using an adaptive finite element approach. In comparison to a dipping anisotropy case, the first order spatial derivatives of the strike-parallel components arise in the partial differential equations for generally anisotropic media, which cause a non-symmetric linear system of equations for finite element modeling. The adaptive finite element method is employed to obtain numerical solutions on a sequence of refined unstructured triangular meshes, which allows for arbitrary model geometries including bathymetry and dipping layers. Numerical results of a 2D anisotropic model show both anisotropy strike and dipping angles have great influence on the marine CSEM responses.

Finite Element Method for a Kind of Two-Dimensional Space-Fractional Diffusion Equation with Its Implementation

In this article, we consider a two-dimensional symmetric space-fractional diffusion equation in which the space fractional derivatives are defined in Riesz potential sense. The well-posed feature is guaranteed by energy inequality. To solve the diffusion equation, a fully discrete form is established by employing Crank-Nicolson technique in time and Galerkin finite element method in space. The stability and convergence are proved and the stiffness matrix is given analytically. Three numerical examples are given to confirm our theoretical analysis in which we find that even with the same initial condition, the classical and fractional diffusion equations perform differently but tend to be uniform diffusion at last.

FINITE ELEMENT METHOD FOR SOLVING TWO-DIMENSIONAL DIFFUSION-REACTION EQUATIONS OF BOUNDARY LAYER TYPE IN POROUS CATALYST PELLET

In this paper,finite element method(FEM)is used to solve two-dimensional diffu-sion-reaction equations of boundary layer type.This kind of equations are usually too complicatedand diffcult to be solved by applying the traditional methods used in chemical engineering becauseof the steep gradients of concentration and temperature.But,these difficulties are easy to be over-comed when the FEM is used.The integraded steps of solving this kind of problems by the FEMare presented in this paper.By applying the FEM to the two actual examples,the conclusion can bereached that the FEM has the advantages of simplicity and good accuracy.

Reduced-order proper orthogonal decomposition extrapolating finite volume element format for two-dimensional hyperbolic equations

This paper is concerned with establishing a reduced-order extrapolating fi- nite volume element （FVE） format based on proper orthogonal decomposition （POD） for two-dimensional （2D） hyperbolic equations. For this purpose, a semi discrete variational format relative time and a fully discrete FVE format for the 2D hyperbolic equations are built, and a set of snapshots from the very few FVE solutions are extracted on the first very short time interval. Then, the POD basis from the snapshots is formulated, and the reduced-order POD extrapolating FVE format containing very few degrees of freedom but holding sufficiently high accuracy is built. Next, the error estimates of the reduced-order solutions and the algorithm procedure for solving the reduced-order for- mat are furnished. Finally, a numerical example is shown to confirm the correctness of theoretical conclusions. This means that the format is efficient and feasible to solve the 2D hyperbolic equations.