Mobility and dynamic erosion process of granular flow:insights from numerical investigation using material point method

In order to understand the dynamics of granular flow on an erodible base soil,in this paper,a series of material point method-based granular column collapse tests were conducted to investigate numerically the mobility and dynamic erosion process of granular flow subjected to the complex settings,i.e.,the aspect ratio,granular mass,friction and dilatancy resistance,gravity and presence of water.A set of power scaling laws were proposed to describe the final deposit characteristics of granular flow by the relations of the normalized run-out distance and the normalized final height of granular flow against the aspect ratio,being greatly affected by the complex geological settings,e.g.,granular mass,the friction and dilatancy resistance of granular soil,and presence of water in granular flow.An index of the coefficient of friction of granular soil was defined as a ratio of the target coefficient of friction over the initial coefficient of friction to quantify the scaling extent of friction change(i.e.,friction strengthening or weakening).There is a characteristic aspect ratio of granular column corresponding to the maximum mobility of granular flow with the minimum index of the apparent coefficient of friction.The index of the repose coefficient of friction of granular flow decreased gradually with the increase in aspect ratio because higher potential energy of granular column at a larger aspect ratio causes a larger kinetic energy of granular soil to weaken the friction of granular soil as a kind of velocity-related friction weakening.An increase in granular mass reduces gradually the indexes of the apparent and repose coefficients of friction of granular soil to enhance the mobility of granular flow.The mobility of granular flow increases gradually with the decrease in friction angle or increase in dilatancy angle of granular soil.However,the increase of gravity accelerates granular flow but showing the same final deposit profile without any dependence on gravity.The mobility of granular flow increases gradually by lowering the indexes of the apparent and repose coefficients of friction of granular flow while changing the surroundings,in turn,the dry soil,submerged soil and saturated soil,implying a gradually increased excessive mobility of granular flow with the friction weakening of granular soil.Presence of water in granular flow may be a potential catalyzer to yield a long run-out granular flow,as revealed in comparison of water-absent and water-present granular flows.In addition,the dynamic erosion and entrainment of based soil induced by granular flow subjected to the complex geological settings,i.e.,the aspect ratio,granular mass,gravity,friction and dilatancy resistance,and presence of water,were comprehensively investigated as well.

Numerical manifold method for thermo-mechanical coupling simulation of fractured rock mass

As a calculation method based on the Galerkin variation,the numerical manifold method(NMM)adopts a double covering system,which can easily deal with discontinuous deformation problems and has a high calculation accuracy.Aiming at the thermo-mechanical(TM)coupling problem of fractured rock masses,this study uses the NMM to simulate the processes of crack initiation and propagation in a rock mass under the influence of temperature field,deduces related system equations,and proposes a penalty function method to deal with boundary conditions.Numerical examples are employed to confirm the effectiveness and high accuracy of this method.By the thermal stress analysis of a thick-walled cylinder(TWC),the simulation of cracking in the TWC under heating and cooling conditions,and the simulation of thermal cracking of the SwedishÄspöPillar Stability Experiment(APSE)rock column,the thermal stress,and TM coupling are obtained.The numerical simulation results are in good agreement with the test data and other numerical results,thus verifying the effectiveness of the NMM in dealing with thermal stress and crack propagation problems of fractured rock masses.

Generalized nth-Order Perturbation Method Based on Loop Subdivision Surface Boundary Element Method for Three-Dimensional Broadband Structural Acoustic Uncertainty Analysis

In this paper,a generalized nth-order perturbation method based on the isogeometric boundary element method is proposed for the uncertainty analysis of broadband structural acoustic scattering problems.The Burton-Miller method is employed to solve the problem of non-unique solutions that may be encountered in the external acoustic field,and the nth-order discretization formulation of the boundary integral equation is derived.In addition,the computation of loop subdivision surfaces and the subdivision rules are introduced.In order to confirm the effectiveness of the algorithm,the computed results are contrasted and analyzed with the results under Monte Carlo simulations(MCs)through several numerical examples.

Impact dynamics of granular flow on rigid barriers:insights from numerical investigation using material point method

In order to advance the understanding of the impact dynamics of granular flow in complex geological settings,this paper studied the impact dynamics of granular flow on rigid barriers with a number of Material Point Method(MPM)numerical tests.The impact behavior of granular flow on a rigid barrier was characterized by the initial dynamic impact stage,dynamic surge impact stage,compression impact stage and static stage of granular flow,where the impact force of granular flow was comprised of the dynamic and static forces of granular flow.The impact behavior of granular flow on a rigid barrier was characterized by the states of the fast or slow impact behavior of granular flow.The angle of slope and aspect ratio of granular soil greatly affected the impact behavior of granular flow on a column rigid barrier,where a power model was proposed to quantify the residual(Fnr)-over-maximum(Fnmax)normal impact force ratio of granular flow Fnr⁄Fnmax incorporating the effects of the angle of slope and aspect ratio of granular soil.With the increase of the column rigid barrier up to the semi-infinite column rigid barrier,the impact dynamics of granular flow gradually increased up to a maximum by progressively transforming the overflow into the dynamic surge impact of the incoming flow on the rigid barrier to capture more granular soil of granular flow against the rigid barrier.Presence of water in granular flow,i.e.,a mixture of solid and liquid in granular flow,yielded a dynamic coupling contribution of the solid and liquid,being accompanied by the whole dynamic process of granular flow,on the impact behavior of granular flow on a rigid barrier,where the liquid-phase material of granular flow,i.e.,the water,was predominant to contribute on the normal impact force of granular flow in comparison with the solid-phase material of granular flow.In addition,other factors,e.g.,the shape of rigid barrier(i.e.,the column barrier,arch barrier and circle barrier),and the gravity(i.e.,in the gravitational environment of the Moon,Earth and Mars),greatly affected the impact behavior of granular flow on a rigid barrier as well.

Numerical simulation study on the mold strength of magnetic mold casting based on a coupled electromagnetic-structural method

The properties of the magnetic mold in magnetic mold casting directly determine the quality of the final cast parts.In this study,the magnetic mold properties in magnetic mold casting,were studied utilizing a coupled electromagnetic-structural method through numerical simulation.This study investigated key factors including equivalent stress,the distribution of tensile and compressive stresses,and the area ratio of tensile stress.It compared molds made entirely of magnetic materials with those made partially of magnetic materials.Simulation results indicate that as current increases from 4 A to 8 A,both the initial magnetic mold and the material-replaced magnetic mold initially show an increasing trend in equivalent stress,tensile-compressive stress,and the area ratio of tensile stress,peaking at 6 A before declining.After material replacement,the area ratio of tensile stress at 6 A decreases to 19.84%,representing a reduction of 29.72%.Magnetic molds comprising a combination of magnetic and non-magnetic materials exhibit sufficient strength and a reduced area ratio of tensile stress compared to those made entirely from magnetic materials.This study provides valuable insights for optimizing magnetic mold casting processes and offers practical guidance for advancing the application of magnetic molds.

Numerical Simulation of Oil-Water Two-Phase Flow in Low Permeability Tight Reservoirs Based on Weighted Least Squares Meshless Method

In response to the complex characteristics of actual low-permeability tight reservoirs,this study develops a meshless-based numerical simulation method for oil-water two-phase flow in these reservoirs,considering complex boundary shapes.Utilizing radial basis function point interpolation,the method approximates shape functions for unknown functions within the nodal influence domain.The shape functions constructed by the aforementioned meshless interpolation method haveδ-function properties,which facilitate the handling of essential aspects like the controlled bottom-hole flow pressure in horizontal wells.Moreover,the meshless method offers greater flexibility and freedom compared to grid cell discretization,making it simpler to discretize complex geometries.A variational principle for the flow control equation group is introduced using a weighted least squares meshless method,and the pressure distribution is solved implicitly.Example results demonstrate that the computational outcomes of the meshless point cloud model,which has a relatively small degree of freedom,are in close agreement with those of the Discrete Fracture Model(DFM)employing refined grid partitioning,with pressure calculation accuracy exceeding 98.2%.Compared to high-resolution grid-based computational methods,the meshless method can achieve a better balance between computational efficiency and accuracy.Additionally,the impact of fracture half-length on the productivity of horizontal wells is discussed.The results indicate that increasing the fracture half-length is an effective strategy for enhancing production from the perspective of cumulative oil production.

Comparative Studies between Picard’s and Taylor’s Methods of Numerical Solutions of First Ordinary Order Differential Equations Arising from Real-Life Problems

To solve the first-order differential equation derived from the problem of a free-falling object and the problem arising from Newton’s law of cooling, the study compares the numerical solutions obtained from Picard’s and Taylor’s series methods. We have carried out a descriptive analysis using the MATLAB software. Picard’s and Taylor’s techniques for deriving numerical solutions are both strong mathematical instruments that behave similarly. All first-order differential equations in standard form that have a constant function on the right-hand side share this similarity. As a result, we can conclude that Taylor’s approach is simpler to use, more effective, and more accurate. We will contrast Rung Kutta and Taylor’s methods in more detail in the following section.

Regular perturbation solution of Couette flow(non-Newtonian)between two parallel porous plates: a numerical analysis with irreversibility

The unavailability of wasted energy due to the irreversibility in the process is called the entropy generation.An irreversible process is a process in which the entropy of the system is increased.The second law of thermodynamics is used to define whether the given system is reversible or irreversible.Here,our focus is how to reduce the entropy of the system and maximize the capability of the system.There are many methods for maximizing the capacity of heat transport.The constant pressure gradient or motion of the wall can be used to increase the heat transfer rate and minimize the entropy.The objective of this study is to analyze the heat and mass transfer of an Eyring-Powell fluid in a porous channel.For this,we choose two different fluid models,namely,the plane and generalized Couette flows.The flow is generated in the channel due to a pressure gradient or with the moving of the upper lid.The present analysis shows the effects of the fluid parameters on the velocity,the temperature,the entropy generation,and the Bejan number.The nonlinear boundary value problem of the flow problem is solved with the help of the regular perturbation method.To validate the perturbation solution,a numerical solution is also obtained with the help of the built-in command NDSolve of MATHEMATICA 11.0.The velocity profile shows the shear thickening behavior via first-order Eyring-Powell parameters.It is also observed that the profile of the Bejan number has a decreasing trend against the Brinkman number.Whenηi→0(i=1,2,3),the Eyring-Powell fluid is transformed into a Newtonian fluid.

A COUPLING DIFFERENCE SCHEME FOR THE NUMERICAL SOLUTION OF A SINGULAR PERTURBATION PROBLEM

A singularly perturbed problem without turning points was considered. On a special discretization mesh, a coupling difference scheme, resulting from central difference scheme and Abrahamsson-Keller-Kreiss box scheme, was proposed and the second order convergence, uniform in the small parameter, was proved. Finally, numerical results were provided.

A Uniformly Convergent Numerical Method Using Weak Formulation for Singularly Perturbed Differential Equations

A numerical method using weak formulation is proposed to solve singularly perturbed differential equations. The numerical method is applied to both linear and nonlinear perturbation problems. A linear differential equation is solved using its weak formulation with a test space composed of exponential functions matching boundary layers. A nonlinear singular perturbation problem is converted into a system of linear differentiation equations. Then each linear differential equation is solved iteratively. The uniform convergence, which is independent of the singular perturbation parameter, is numerically verified.

Numerical modeling of stress perturbations caused by geometric changes of salt bodies

We simulated the stress changes around a salt basin using a static salt structure model under compressive stress background to investigate the stress perturbation caused by different salt body shapes and amplitudes. We designed a two-layer salt model with three bulges and sags using finite element methods to calculate the stress perturbation around the salt. The results show that salt shape is closely related to the stress perturbation in the sediments around the salt, and the fluctuations of the bulge and sag(smooth or steep) can also affect the stress perturbation magnitude. Extrusion horizontal stress, normal stress, and out-of-plane stress on the plane would occur near the salt uplift in a compressive tectonic stress environment. In contrast, tensile horizontal stress, out-of-plane stress, and vertical stress would occur near the salt sag. In addition, smoother bulges are associated with smaller produced stress perturbations, and steeper sags are associated with a greater reduction of stress perturbation in the sediment. The stress of a salt structure in western Kelasu of the Kuqa depression was simulated and the applicability of previous conclusions regarding this structure was verified. These conclusions provide scientific basis for the prediction of stress perturbations around salt basin systems.

A Numerical Method for Nonlinear Singularly Perturbed Multi-Point Boundary Value Problem

We consider a uniform finite difference method for nonlinear singularly perturbed multi-point boundary value problem on Shishkin mesh. The problem is discretized using integral identities, interpolating quadrature rules, exponential basis functions and remainder terms in integral form. We show that this method is the first order convergent in the discrete maximum norm for original problem (independent of the perturbation parameter ε). To illustrate the theoretical results, we solve test problem and we also give the error distributions in the solution in Table 1 and Figures 1-3.

AN APPLICATION OF THE EXPONENTIAL CUBIC SPLINES TO NUMERICAL SOLUTION OF A SELF-ADJOINT PERTURBATION PROBLEM

We present a class of the second order optimal splines difference schemes derived from ex- ponential cubic splines for self-adjoint singularly perturbed 2-point boundary value problem. We prove an optimal error estimate and give illustrative numerical example.

Calculation of effective temperature for pavement rutting using numerical simulation methods

In order to predict the long-term rutting of asphalt pavement, the effective temperature for pavement rutting is calculated using the numerical simulation method. The transient temperature field of asphalt pavement was simulated based on actual meteorological data of Nanjing. 24-hour rutting development under a transient temperature field was calculated in each month. The rutting depth accumulated under the static temperature field was also estimated and the relationship between constant temperature parameters was analyzed. Then the effective temperature for pavement rutting was determined based on the rutting equivalence principle. The results show that the monthly effective temperature is above 40 t in July and August, while in June and September it ranges from 30 to 40 Rutting development can be ignored when the monthly effective temperature is less than 30 t. The yearly effective temperature for rutting in Nanjing is around 38. 5 t. The long-term rutting prediction model based on the effective temperature can reflect the influences of meteorological factors and traffic time distribution.

Application of ALE Method on the Numerical Simulation of Reinforced Concrete Penetration

LS-DYNA program and the principle of ALE method were introduced, and the target features of the reinforced concrete penetration were analyzed by using the D material model and the ALE method. A numerical simulation has been done to show the penetration visually and veritably. The simulation results are analyzed carefully and explicitly prove their significance to the research of reinforced concrete penetration.

A numerical study of two-and three-dimensional detonation dynamics of pulse detonation engine by the CE/SE method

In this paper, the CE/SE method is developed to simulate the two- and three-dimensional flow-field of Pulse Detonation Engine (PDE). The conservation equations with stiff source terms for chemical reaction are solved in two steps. The detailed analysis of computational results of a PDE with a single detonation tube and a PDE with five detonation tubes are given in this paper. Complex wave systems are observed inside and outside a PDE. For a PDE with 5 detonation tubes, there is a big bow shock produced from a number of little shocks near the open ends of tubes. A lot of vortexes interact with shocks and a large expansion wave propagates forward and backward with respect to the PDE in a semi-oval shape.

Numerical simulation of the welding deformation for the side sill of the bogie frame based on local-global method

Considering the limitation of computational capacity, a new finite element solution is used to simulate the welding deformation of the side sill of railroad car＇ s bogie frame based on the local-global method. Firstly, a volumetric heat source defined by a double ellipsoid is adopted to simulate the thermal distributions of the arc welding process. And then, the local models extracted from the global model are computed with refined meshes. On these bases, the global distortions of the subject studied are ascertained by transferring the inner forces of computed local models to the global model. It indicates that the local-global method is feasible for simulating the large welded structures by comparing the computed results with the corresponding actual measured values. The work provides basis for optimizing the welding sequence and clamping conditions, and has theoretical values and engineering significance in the integral design, manufacturing technique selection of the bogie frame, as well as other kinds of large welded structures.

PERTURBATION FINITE VOLUME METHOD FOR CONVECTIVE-DIFFUSION INTEGRAL EQUATION

A perturbation finite volume(PFV)method for the convective-diffusion integral equa- tion is developed in this paper.The PFV scheme is an upwind and mixed scheme using any higher-order interpolation and second-order integration approximations,with the least nodes similar to the standard three-point schemes,that is,the number of the nodes needed is equal to unity plus the face-number of the control volume.For instance,in the two-dimensional(2-D)case,only four nodes for the triangle grids and five nodes for the Cartesian grids are utilized,respectively.The PFV scheme is applied on a number of 1-D linear and nonlinear problems,2-D and 3-D flow model equations.Comparing with other standard three-point schemes,the PFV scheme has much smaller numerical diffusion than the first-order upwind scheme(UDS).Its numerical accuracies are also higher than the second-order central scheme(CDS),the power-law scheme(PLS)and QUICK scheme.

Surplus Space Method: A New Numerical Model for Prediction of Shallow-seated Magmatic Bodies

Based on the data of field measurement and drilling in the Tongling area, a series of numerical simulations are carried out by using the 'Surplus Space Method' (SSM), which is first put forward in this paper and applied to predict the shallow-seated magmatic bodies. The results of the numerical simulations show the existence and the 3-D shape of a conical magmatic structure at a depth of-1000 m beneath the center of the area: its top offsets southwards and bifurcates to several branches, while its lower part stretches northeastwards and contracts rapidly to a point at about -1000 m depth. This point is reckoned to be a 'sink' of magma system, transferring ore materials and heat energy from the deep magma chamber to the sub-surface apophyses. The preliminary application of the SSM proves that it may be developed as a new detection means for determining the existence of shallow-seated magmatic bodies and analyzing their three-dimensional features.

Stress initialization methods for dynamic numerical simulation of rock mass with high in-situ stress

In the context of deep rock engineering,the in-situ stress state is of major importance as it plays an important role in rock dynamic response behavior.Thus,stress initialization becomes crucial and is the first step for the dynamic response simulation of rock mass in a high in-situ stress field.In this paper,stress initialization methods,including their principles and operating procedures for reproducing steady in-situ stress state in LS-DYNA,are first introduced.Then the most popular four methods,i.e.,explicit dynamic relaxation(DR)method,implicit-explicit sequence method,Dynain file method and quasi-static method,are exemplified through a case analysis by using the RHT and plastic hardening rock material models to simulate rock blasting under in-situ stress condition.Based on the simulations,it is concluded that the stress initialization results obtained by implicit-explicit sequence method and dynain file method are closely related to the rock material model,and the explicit DR method has an obvious advantage in solution time when compared to other methods.Besides that,it is recommended to adopt two separate analyses for the whole numerical simulation of rock mass under the combined action of in-situ stress and dynamic disturbance.