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.

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

A SEMI-ANALYTICAL AND SEMI-NUMERICAL METHOD FOR SOLVING 2-D SOUND-STRUCTURE INTERACTION PROBLEMS

Based on the transfer matrix method and the virtual source simulation technique, this paper proposes a novel semi-analytical and semi-numerical method for solving 2-D sound- structure interaction problems under a harmonic excitation.Within any integration segment, as long as its length is small enough,along the circumferential curvilinear coordinate,the non- homogeneous matrix differential equation of an elastic ring of complex geometrical shape can be rewritten in terms of the homogeneous one by the method of extended homogeneous capacity proposed in this paper.For the exterior fluid domain,the multi-circular virtual source simulation technique is adopted.The source density distributed on each virtual circular curve may be ex- panded as the Fourier's series.Combining with the inverse fast Fourier transformation,a higher accuracy and efficiency method for solving 2-D exterior Helmholtz's problems is presented in this paper.In the aspect of solution to the coupling equations,the state vectors of elastic ring induced by the given harmonic excitation and generalized forces of coefficients of the Fourier series can be obtained respectively by using a high precision integration scheme combined with the method of extended homogeneous capacity put forward in this paper.According to the superposition princi- ple and compatibility conditions at the interface between the elastic ring and fluid,the algebraic equation of system can be directly constructed by using the least square approximation.Examples of acoustic radiation from two typical fluid-loaded elastic rings under a harmonic concentrated force are presented.Numerical results show that the method proposed is more efficient than the mixed FE-BE method in common use.

SEMI-ANALYTICAL AND SEMI-NUMERICAL METHOD FOR DYNAMIC ANALYSIS OF FOUNDATION

A semi-analytical and semi-numerical method is proposed for the dynamic analysis of foundations. The Lamb＇s solution and the approximate formulae were used to establish the relation of the contact force and deflection between the foundation and soil. Therefore, the foundation can be separated from soil and analyzed by FEM as for the static cases. The plate can be treated as that the known forces are acting on the upper surface, and the contact pressure from soil can be represented as the deflection. So that only the plate needs to be divided into elements in the analysis. By this method, a series of vibration problems, including various shapes and rigidities of foundations, different excitation frequencies, were analyzed. Furthermore, it can be used for the embedded foundation. The numerical examples show that this method has simplicity, highly accurate and versatile. It is an effective method for the dynamic analysis of foundations.

Numerical Study on Motion Responses of Two Parallel Ships During Underway Replenishment by the Semi-Analytical HOTP Method

Motion responses of two ships advancing parallel in waves with hydrodynamic interactions are investigated in this paper. Within the framework of the frequency-domain potential flow theory, a semi-analytical higher-order translating-pulsating source(HOTP) method is presented to solve the problems of coupled radiation and diffraction potential. The method employs nine-node bi-quadratic curvilinear elements to discretize the boundary integral equations(BIEs) constructed over the mean wetted surface of the two ship hulls. In order to eliminate the numerical oscillation, analytical quadrature formulas are derived and adopted to evaluate the integrals related to the Froudedependent part of the Green’s function along the horizontal direction in the BIEs. Based on the method, a numerical program is originally coded. Through the calculations of hydrodynamic responses of single ships, the numerical implementation is proved successful. Then the validated program is applied in the investigations on the hydrodynamic interactions of two identical Wigley Ⅲ hulls and the underway replenishment of a frigate and a supply ship in waves with and without stagger, respectively. The comparison between the present computed results with experimental data and numerical solutions of other methods shows that the semi-analytical HOTP method is of higher accuracy than the pulsating source Green’s function method with speed correction and better stability than the traditional HOTP method based on Gauss quadrature. In addition, for two ships with obviously different dimensions,the influence of hydrodynamic interactions on the smaller ship is found to be more noticeable than that on the larger ship, which leads to the differences between the motions of frigate with and without the presence of supply ship.

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.

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.

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.

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.

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.

2D numerical manifold method based on quartic uniform B-spline interpolation and its application in thin plate bending

A new numerical manifold （NMM） method is derived on the basis of quartic uniform B-spline interpolation. The analysis shows that the new interpolation function possesses higher-order continuity and polynomial consistency compared with the conven- tional NMM. The stiffness matrix of the new element is well-conditioned. The proposed method is applied for the numerical example of thin plate bending. Based on the prin- ciple of minimum potential energy, the manifold matrices and equilibrium equation are deduced. Numerical results reveal that the NMM has high interpolation accuracy and rapid convergence for the global cover function and its higher-order partial derivatives.

Analogue correction method of errors and its application to numerical weather prediction

In this paper, an analogue correction method of errors （ACE） based on a complicated atmospheric model is further developed and applied to numerical weather prediction （NWP）. The analysis shows that the ACE can effectively reduce model errors by combining the statistical analogue method with the dynamical model together in order that the information of plenty of historical data is utilized in the current complicated NWP model, Furthermore, in the ACE, the differences of the similarities between different historical analogues and the current initial state are considered as the weights for estimating model errors. The results of daily, decad and monthly prediction experiments on a complicated T63 atmospheric model show that the performance of the ACE by correcting model errors based on the estimation of the errors of 4 historical analogue predictions is not only better than that of the scheme of only introducing the correction of the errors of every single analogue prediction, but is also better than that of the T63 model.

Response Sensitivity Analysis of the Dynamic Milling Process Based on the Numerical Integration Method

As one of the bases of gradient-based optimization algorithms, sensitivity analysis is usually required to calculate the derivatives of the system response with respect to the machining parameters. The most widely used approaches for sensitivity analysis are based on time-consuming numerical methods, such as finite difference methods. This paper presents a semi-analytical method for calculation of the sensitivity of the stability boundary in milling. After transforming the delay-differential equation with time-periodic coefficients governing the dynamic milling process into the integral form, the Floquet transition matrix is constructed by using the numerical integration method. Then, the analytical expressions of derivatives of the Floquet transition matrix with respect to the machining parameters are obtained. Thereafter, the classical analytical expression of the sensitivity of matrix eigenvalues is employed to calculate the sensitivity of the stability lobe diagram. The two-degree-of-freedom milling example illustrates the accuracy and efficiency of the proposed method. Compared with the existing methods, the unique merit of the proposed method is that it can be used for analytically computing the sensitivity of the stability boundary in milling, without employing any finite difference methods. Therefore, the high accuracy and high efficiency are both achieved. The proposed method can serve as an effective tool for machining parameter optimization and uncertainty analysis in high-speed milling.