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.

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.

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.

Wavelet Multi-Resolution Interpolation Galerkin Method for Linear Singularly Perturbed Boundary Value Problems

In this study,a wavelet multi-resolution interpolation Galerkin method(WMIGM)is proposed to solve linear singularly perturbed boundary value problems.Unlike conventional wavelet schemes,the proposed algorithm can be readily extended to special node generation techniques,such as the Shishkin node.Such a wavelet method allows a high degree of local refinement of the nodal distribution to efficiently capture localized steep gradients.All the shape functions possess the Kronecker delta property,making the imposition of boundary conditions as easy as that in the finite element method.Four numerical examples are studied to demonstrate the validity and accuracy of the proposedwavelet method.The results showthat the use ofmodified Shishkin nodes can significantly reduce numerical oscillation near the boundary layer.Compared with many other methods,the proposed method possesses satisfactory accuracy and efficiency.The theoretical and numerical results demonstrate that the order of theε-uniform convergence of this wavelet method can reach 5.

Multiparameter Numerical Investigation of Two Types of Moving Interactions Between the Deep-Sea Mining Vehicle Track Plate and Seabed Soil:Digging and Rotating Motions

To ensure the safe performance of deep-sea mining vehicles(DSMVs),it is necessary to study the mechanical characteristics of the interaction between the seabed soil and the track plate.The rotation and digging motions of the track plate are important links in the contact between the driving mechanism of the DSMV and seabed soil.In this study,a numerical simulation is conducted using the coupled Eulerian–Lagrangian(CEL)large deformation numerical method to investigate the interaction between the track plate of the DSMV and the seabed soil under two working conditions:rotating condition and digging condition.First,a soil numerical model is established based on the elastoplastic mechanical characterization using the basic physical and mechanical properties of the seabed soil obtained by in situ sampling.Subsequently,the soil disturbance mechanism and the dynamic mechanical response of the track plate under rotating and digging conditions are obtained through the analysis of the sensitivity of the motion parameters,the grouser structure,the layered soil features and the soil heterogeneity.The results indicate that the above parameters remarkably influence the interaction between the DSMV and the seabed soil.Therefore,it is important to consider the rotating and digging motion of the DSMV in practical engineering to develop a detailed optimization design of the track plate.

Numerical investigation of geostress influence on the grouting reinforcement effectiveness of tunnel surrounding rock mass in fault fracture zones

Grouting is a widely used approach to reinforce broken surrounding rock mass during the construction of underground tunnels in fault fracture zones,and its reinforcement effectiveness is highly affected by geostress.In this study,a numerical manifold method(NMM)based simulator has been developed to examine the impact of geostress conditions on grouting reinforcement during tunnel excavation.To develop this simulator,a detection technique for identifying slurry migration channels and an improved fluid-solid coupling(FeS)framework,which considers the influence of fracture properties and geostress states,is developed and incorporated into a zero-thickness cohesive element(ZE)based NMM(Co-NMM)for simulating tunnel excavation.Additionally,to simulate coagulation of injected slurry,a bonding repair algorithm is further proposed based on the ZE model.To verify the accuracy of the proposed simulator,a series of simulations about slurry migration in single fractures and fracture networks are numerically reproduced,and the results align well with analytical and laboratory test results.Furthermore,these numerical results show that neglecting the influence of geostress condition can lead to a serious over-estimation of slurry migration range and reinforcement effectiveness.After validations,a series of simulations about tunnel grouting reinforcement and tunnel excavation in fault fracture zones with varying fracture densities under different geostress conditions are conducted.Based on these simula-tions,the influence of geostress conditions and the optimization of grouting schemes are discussed.

FFT-Based Numerical Method for Nonlinear Elastic

In theoretical research pertaining to sealing, a contact model must be used to obtain the leakage channel. However, for elastoplastic contact, current numerical methods require a long calculation time. Hyperelastic contact is typically simplifed to a linear elastic contact problem, which must be improved in terms of calculation accuracy. Based on the fast Fourier transform, a numerical method suitable for elastoplastic and hyperelastic frictionless contact that can be used for solving two-dimensional and three-dimensional (3D) contact problems is proposed herein. The nonlinear elastic contact problem is converted into a linear elastic contact problem considering residual deformation (or the equivalent residual deformation). Results from numerical simulations for elastic, elastoplastic, and hyperelastic contact between a hemisphere and a rigid plane are compared with those obtained using the fnite element method to verify the accuracy of the numerical method. Compared with the existing elastoplastic contact numerical methods, the proposed method achieves a higher calculation efciency while ensuring a certain calculation accuracy (i.e., the pressure error does not exceed 15%, whereas the calculation time does not exceed 10 min in a 64 × 64 grid). For hyperelastic contact, the proposed method reduces the dependence of the approximation result on the load, as in a linear elastic approximation. Finally, using the sealing application as an example, the contact and leakage rates between complicated 3D rough surfaces are calculated. Despite a certain error, the simplifed numerical method yields a better approximation result than the linear elastic contact approximation. Additionally, the result can be used as fast solutions in engineering applications.

Numerical Simulation of a Viscoelastic Thinning Process for Preparing Flexible Glasses by Redrawing Method

The forming process of the flexible ultrathin glasses(UTG)prepared by the redrawing method was numerically simulated using ANSYS Polyflow software.In the forming process by the redrawing method,temperature,viscosity,transverse and longitudinal velocity distribution of the glasses with different compositions were studied.Furthermore,the influence of these factors on the width and thickness of the flexible glass plate was investigated.It is found that the internal and external heat exchange of glass has a dominant influence on the viscosity variation during the UTG forming process,which is inconsistent with the general viscosity-temperature dependence.The glass that first reaches the lower limit of forming viscosity can significantly resist the shrinking effect caused by surface tension,making the glass wider during the forming.If the original glass width remains unchanged,the glass thickness or feeding speed is reduced,wider and thinner flexible glasses can be produced.

Comparative Analysis for Evaluating Wind Energy Resources Using Intelligent Optimization Algorithms and Numerical Methods

Statistical distributions are used to model wind speed,and the twoparameters Weibull distribution has proven its effectiveness at characterizing wind speed.Accurate estimation of Weibull parameters,the scale(c)and shape(k),is crucial in describing the actual wind speed data and evaluating the wind energy potential.Therefore,this study compares the most common conventional numerical(CN)estimation methods and the recent intelligent optimization algorithms(IOA)to show how precise estimation of c and k affects the wind energy resource assessments.In addition,this study conducts technical and economic feasibility studies for five sites in the northern part of Saudi Arabia,namely Aljouf,Rafha,Tabuk,Turaif,and Yanbo.Results exhibit that IOAs have better performance in attaining optimal Weibull parameters and provided an adequate description of the observed wind speed data.Also,with six wind turbine technologies rating between 1 and 3MW,the technical and economic assessment results reveal that the CN methods tend to overestimate the energy output and underestimate the cost of energy($/kWh)compared to the assessments by IOAs.The energy cost analyses show that Turaif is the windiest site,with an electricity cost of$0.016906/kWh.The highest wind energy output is obtained with the wind turbine having a rated power of 2.5 MW at all considered sites with electricity costs not exceeding$0.02739/kWh.Finally,the outcomes of this study exhibit the potential of wind energy in Saudi Arabia,and its environmental goals can be acquired by harvesting wind energy.

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.

Theoretical and numerical studies of rock breaking mechanism by double disc cutters

A plane mechanical model of rock breaking process by double disc cutter at the center of the cutterhead is established based on contact mechanics to analyze the stress evolution in the rock broken by cutters with different spacings. A continuous-discontinuous coupling numerical method based on zero-thickness cohesive elements is developed to simulate rock breaking using double cutters. The process, mechanism,and characteristics of rock breaking are comprehensively analyzed from five aspects: peak force, breaking form, breaking efficiency, crack mode, and breaking degree. The results show that under the penetrating action of cutters, dense cores are formed due to shear failure under respective cutters. The tensile cracks propagate in the rock, and then rock chips form with increasing penetration depth. When the cutter spacing is increased from 10 to 80 mm, the peak force gradually increases, the rock breaking range increases first and then decreases, the specific energy decreases first and then rises, and the breaking coefficient of intermediate rock decreases from 0.955 to 0.788. The area of rock breaking is positively correlated with the length of the tensile crack. Furthermore, the length of the tensile crack accounts for 14.4%–33.6% of the total crack length.

A Dimension-Splitting Variational Multiscale Element-Free Galerkin Method for Three-Dimensional Singularly Perturbed Convection-Diffusion Problems

By introducing the dimensional splitting(DS)method into the multiscale interpolating element-free Galerkin(VMIEFG)method,a dimension-splitting multiscale interpolating element-free Galerkin(DS-VMIEFG)method is proposed for three-dimensional(3D)singular perturbed convection-diffusion(SPCD)problems.In the DSVMIEFG method,the 3D problem is decomposed into a series of 2D problems by the DS method,and the discrete equations on the 2D splitting surface are obtained by the VMIEFG method.The improved interpolation-type moving least squares(IIMLS)method is used to construct shape functions in the weak form and to combine 2D discrete equations into a global system of discrete equations for the three-dimensional SPCD problems.The solved numerical example verifies the effectiveness of the method in this paper for the 3D SPCD problems.The numerical solution will gradually converge to the analytical solution with the increase in the number of nodes.For extremely small singular diffusion coefficients,the numerical solution will avoid numerical oscillation and has high computational stability.

Numerical Simulation of Aerodynamic Interaction Effects in Coaxial Compound Helicopters

The so-called coaxial compound helicopter features two rigid coaxial rotors,and possesses high-speed capabilities.Nevertheless,the small separation of the coaxial rotors causes severe aerodynamic interactions,which require careful analysis.In the present work,the aerodynamic interaction between the various helicopter components is investigated by means of a numerical method considering both hover and forward flight conditions.While a sliding mesh method is used to deal with the rotating coaxial rotors,the Reynolds-Averaged Navier-Stokes(RANS)equations are solved for the flow field.The Caradonna&Tung(CT)rotor and Harrington-2 coaxial rotor are considered to validate the numerical method.The results show that the aerodynamic interaction of the two rigid coaxial rotors significantly influences hover’s induced velocity and pressure distribution.In addition,the average thrust of an isolated coaxial rotor is smaller than that of the corresponding isolated single rotor.Compared with the isolated coaxial rotor,the existence of the fuselage results in an increment in the thrust of the rotors.Furthermore,these interactions between the components of the considered coaxial compound helicopter decay with an increase in the advance ratio.

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.