Comparison of Fully 3-D,Simplified 3-D and 2-D Numerical Simulations on a Circulating Fluidized Bed Boiler

Based on the fully three-dimensional(3-D)and two-dimensional(2-D)comprehensive CFD(Computational Fluid Dynamics)combustion models for a circulating fluidized bed boiler,a simplified 3-D computational domain considering the corrections of furnace side wall openings is proposed.It aims to compensate for the deficiencies of the large amount of computation in the fully 3-D model and improve the air and gas flow treatments at the openings in the simplified 2-D model.Three different computational domains,named as the fully 3-D model,simplified 3-D model and 2-D model,were implemented to perform a comparative CFD analysis in an ultra-supercritical circulating fluidized bed boiler including the hydrodynamics,penetration depth of secondary air,temperature and species distribution.The simulation results computed by the simplified3-D model yield better agreement with the fully 3-D simulation results than those of the 2-D model.The simplified 3-D model is recommended as an alternative computational domain for the conventional 2-D model in the numerical simulation of large-scale circulating fluidized bed boiler.

Numerical Simulation of Thermo-Hydro-Mechanical-Chemical Response Caused by CO_(2)Injection into Saline Geological Formations:A Case Study from the Ordos Project,China

Thermo-hydro-mechanical-chemical(THMC)interactions are prevalent during CO_(2)geological sequestration(CGS).In this study,a sequential coupling THMC numerical simulation program was constructed,which can be used to explore the following issues of CGS:fluid and heat flow,solute transport;stresses,displacements and rock failures related to geo-mechanical effects;equilibrium and kinetic chemical reactions;chemical damage to mechanical properties of the rock.Then,the coupling program was applied to the Ordos CGS Project to study the formation response under the multifield interaction caused by CO_(2)injection.The simulation results show that the mechanical process dominates the short CO_(2)injection period.Specifically,the formation’s permeability near the injection well increases by 43%,due to the reduction of effective stress,which significantly promotes the lateral migration of CO_(2).When the injection rate exceeds 0.15 million tons per year,the cohesion of the reservoir rock is not enough to resist the shear force inside the rock and rock failure may occur.During the subsequent long-term sequestration period(200 years),the influence of mineral reactions gradually increases.Due to calcite dissolution,the shear modulus of caprock is predicted to decrease by 7.6%,which will to some extent increase the risk of rock failure.

2-D NUMERICAL SIMULATION OF FLOODING EFFECTS CAUSED BY SOUTH-TO-NORTH WATER TRANSFER PROJECT

Since the General Channel designed for the South-to-North Water Transfer Project in China has to cross many rivers and streams flowing from west to east, there are potentially serious effects additional flooding on the western side of the project alignment. Therefore, a 2-D numerical model for forecasting basin flood disasters was established and verified using historical flood data. The model was applied to researching the interaction between the proposed Project and flooding events for 5 streams in the Anyang River reach as a representative case study. Simulated results indicate that the model could correctly forecast the flood, submerged area and depths, and water surface elevations along the left side of the channel. The discharge capacity and location of hydraulic structures in the transfer canal alignment were analyzed. Then adjustments to the dimensions and positioning of proposed hydraulic structures were recommended at intersections, especially the addition of a channel to transfer flood water from one stream to another, which can effectively limit the sluice and protect the Anyang City from flooding.

Numerical simulation study on shear resistance of anchorage joints considering tensile-shear fracture criterion of 2G-NPR bolt

2G-NPR bolt (the 2nd generation Negative Poisson’s Ratio bolt) is a new type of bolt with high strength, high toughness and no yield platform. It has signifcant efects on improving the shear strength of jointed rock mass and controlling the stability of surrounding rock. To achieve an accurate simulation of bolted joint shear tests, we have studied a numerical simulation method that takes into account the 2G-NPR bolt's tensile–shear fracture criterion. Firstly, the indoor experimental study on the tensile–shear mechanical properties of 2G-NPR bolt is carried out to explore its mechanical properties under diferent tensile–shear angles, and the fracture criterion of 2G-NPR bolt considering the tensile–shear angle is established. Then, a three-dimensional numerical simulation method considering the tensile–shear mechanical constitutive and fracture criterion of 2G-NPR bolt, the elastoplastic mechanical behavior of surrounding rock and the damage and deterioration of grouting body is proposed. The feasibility and accuracy of the method are verifed by comparing with the indoor shear test results of 2G-NPR bolt anchorage joints. Finally, based on the numerical simulation results, the deformation and stress of the bolt, the distribution of the plastic zone of the rock mass, the stress distribution and the damage of the grouting body are analyzed in detail. The research results can provide a good reference value for the practical engineering application and shear mechanical performance analysis of 2G-NPR bolt.

Numerical Simulation of Low Cycle Fatigue Behavior of Ti_(2)AlNb Alloy Subcomponents

Many titanium alloy subcomponents are subjected to fatigue loading in aerospace engineering,resulting in fatigue failure.The fatigue behavior of Ti_(2)AlNb alloy subcomponents was investigated based on the Seeger fatigue life theory and the improved Lemaitre damage evolution theory.Firstly,the finite element models of the standard openhole specimen and Y-section subcomponents have been established by ABAQUS.The damage model parameters were determined by fatigue tests,and the reliability of fatigue life simulation results of the Ti_(2)AlNb alloy standard open-hole specimen was verified.Meanwhile,the fatigue life of Ti_(2)AlNb alloy Y-section subcomponents was predicted.Under the same initial conditions,the average error of fatigue life predicted by two different models was 20.6%.Finally,the effects of loading amplitude,temperature,and Y-interface angle on fatigue properties of Ti_(2)AlNb Y-section subcomponents were investigated.These results provide a new idea for evaluating the fatigue life of various Ti_(2)AlNb alloy subcomponents.

A Numerical Simulation Study of Temperature and Pressure Effects on the Breakthrough Pressure of CO_(2)in Unsaturated Low-permeability Rock Core

Gas breakthrough pressure is a key parameter to evaluate the sealing capacity of caprock,and it also plays important roles in safety and capacity of CO_(2)geological storage.Based on the published experimental results,we present numerical simulations on CO_(2)breakthrough pressure in unsaturated low-permeability rock under 9 multiple P-T conditions(which can keep CO_(2)in gaseous,liquid and supercritical states)and thus,a numerical method which can be used to accurately predict CO_(2)breakthrough pressure on rock-core scale is proposed.The simulation results show that CO_(2)breakthrough pressure and breakthrough time are exponential correlated with P-T conditions.Meanwhile,pressure has stronger effects on experimental results than that of temperature.Moreover,we performed sensitivity studies on the pore distribution indexλ(0.6,0.7,0.8,and 0.9)in van Genuchten-Muale model.Results show that with the increase ofλ,CO_(2)breakthrough pressure and breakthrough time both show decreasing trends.In other words,the larger the value ofλis,the better the permeability of the caprock is,and the worse the CO_(2)sealing capacity is.The numerical method established in this study can provide an important reference for the prediction of gas breakthrough pressure on rock-core scale and for related numerical studies.

Effects of CH_(4)/CO_(2) multi-component gas on components and properties of tight oil during CO_(2) utilization and storage: Physical experiment and composition numerical simulation

An essential technology of carbon capture, utilization and storage-enhanced oil recovery (CCUS-EOR) for tight oil reservoirs is CO_(2) huff-puff followed by associated produced gas reinjection. In this paper, the effects of multi-component gas on the properties and components of tight oil are studied. First, the core displacement experiments using the CH_(4)/CO_(2) multi-component gas are conducted to determine the oil displacement efficiency under different CO_(2) and CH_(4) ratios. Then, a viscometer and a liquid density balance are used to investigate the change characteristics of oil viscosity and density after multi-component gas displacement with different CO_(2) and CH_(4) ratios. In addition, a laboratory scale numerical model is established to validate the experimental results. Finally, a composition model of multi-stage fractured horizontal well in tight oil reservoir considering nano-confinement effects is established to investigate the effects of multi-component gas on the components of produced dead oil and formation crude oil. The experimental results show that the oil displacement efficiency of multi-component gas displacement is greater than that of single-component gas displacement. The CH_(4) decreases the viscosity and density of light oil, while CO_(2) decreases the viscosity but increases the density. And the numerical simulation results show that CO_(2) extracts more heavy components from the liquid phase into the vapor phase, while CH_(4) extracts more light components from the liquid phase into the vapor phase during cyclic gas injection. The multi-component gas can extract both the light components and the heavy components from oil, and the balanced production of each component can be achieved by using multi-component gas huff-puff.

2-D NUMERICAL SIMULATION OF CRUSH BEDROCK RIVER

In this paper, the erosion-resisting coefficient was introduced to computebed deformation in a crush bedrock river. In the case of crush bedrock, there has been no propercontrol equation to describe bed stability, which leads to difficulty in calculation of the beddeformation with conventional methods. The data from field survey were used to give thee-rosion-resisting capability with an appropriate coefficient. After the determination oflongitudinal distribution expressed by polynomial regression and transversal distribution expressedby normal distribution function, the plane distribution of erosion-resisting coefficient in a crushbedrock river was obtained. With the computational results from a 2-D horizontal flow mathematicalmodel, the erosion-resisting coefficient and controlling condition of local stability were employedto compute the values of bed deformation when riverbed is stable. The above method was applied in acase study, and the computational results of flow and bed deformations are in good a-greement withphysical model test data.

2-D NUMERICAL SIMULATION OF RADIONUCLIDE TRANSPORT IN THE LOWER YANGTZE RIVER

The assessment of the radiological impact of the liquid discharges from nuclear power plants is a major issue for the envi- ronmental protection. In this study, a numerical model for the radionuclide transport in the aquatic environment is built, based on the hydrodynamic equations, including the complete set of Saint-Venant equations, the sediment transport equations, with consideration of several different particle sizes and the deposition and erosion of the suspended sediments, and the radionuclide transport equations The exchanges of radionuclides between water, suspended matter and bed sediments are described in terms of kinetic transfer coe- fficients. The model is used to simulate the transport of the radionuclides discharged from a planned nuclear power plant project to be sited along the lower Yangtze River. From the model results, one may see the detailed temporal-spatial evolution of the radio- nuclide contamination in the solution, in the suspended matter and in the bed sediments. The model can be used as a basic tool for studying the environmental impacts of the liquid discharges from nuclear facilities on a river system.

Two Improvements on Numerical Simulation of 2-DOF Vortex-Induced Vibration with Low Mass Ratio

Till now, there have been lots of researches on numerical simulation of vortex-induced vibration. Acceptable results have been obtained for fixed cylinders with low Reynolds number. However, for responses of 2-DOF vortex-induced vibration with low mass ratio, the accuracy is not satisfactory, especially for the maximum amplitudes. In Jauvtis and Williamson＇s work, the maximum amplitude of the cylinder with low mass ratio m＊=2.6 can reach as large as 1.5D to be called as the ＂super-upper branch＂, but from current literatures, few simulation results can achieve such value, even fail to capture the upper branch. Besides, it is found that the amplitude decays too fast in the lower branch with the RANS-based turbulence model. The reason is likely to be the defects of the turbulence model itself in the prediction of unsteady separated flows as well as the unreasonable setting of the numerical simulation parameters. Aiming at above issues, a modified turbulence model is proposed in this paper, and the effect of the acceleration of flow field on the response of vortex-induced vibration is studied based on OpenFOAM. By analyzing the responses of amplitude, phase and trajectory, frequency and vortex mode, it is proved that the vortex-induced vibration can be predicted accurately with the modified turbulence model under appropriate flow field acceleration.

Numerical simulation study on the effect of temperature on the restricted diffusion in porous media

It is of great significance to study how temperature affects the restricted diffusion in pores for an accurate evaluation of reservoir physical properties by using nuclear magnetic resonance(NMR)diffusion-transverse relaxation(D-T2)spectrum under reservoir temperature conditions.In this paper,we simulate the restricted diffusion and twodimensional(2D)NMR D-T2 spectra of water molecules at different temperatures using random-walk method.The one-dimensional(1D)restricted diffusion simulation results show that the diffusion coefficient in the pore at room temperature decays with the diffusion time and eventually reaches a plateau.Under the condition of long-time diffusion,the ratio of restricted diffusion coefficient to bulk diffusion coefficient at different temperatures tends to be the same constant.With the increase in temperature,the simulated D-T2 spectra also gradually move upward.The simulated D-T2 spectra at all temperatures are consistent with the Pade interpolation equation.In addition,the results calculated by Pade interpolation equation demonstrate that the degree of temperature influence on the D-T2 spectrum of rock is quantitatively related to the pore radius,porosity and cementation index.

Numerical simulation of enhancing coalbed methane recovery by injecting CO_2 with heat injection

The technology used to enhance coalbed methane(CBM) recovery by injecting CO_2(CO_2-ECBM) with heat, combining heat injection with CO_2 injection, is still in its infancy; therefore, theoretical studies of this CO_2-ECBM technology should be perused. First, the coupling equations of the di usion–adsorption–seepage–heat transfer fields of gas are established. The displacement processes under di erent pressures and temperatures are simulated by COMSOL. Finally, the displacement effects, a comparison of the CO_2 storage capacity with the CH_4 output and the e ective influencing radius of CO_2 injection are analyzed and discussed. The results show that(1) the displacement pressure and temperature are two key factors influencing the CH_4 output and the CO_2 storage capacity, and the increase in the CO_2 storage capacity is more sensitive to temperature and pressure than the CH_4 output.(2) The gas flow direction is from the injection hole to the discharge hole during the displacement process, and the regions with high velocity are concentrated at the injection hole and the discharge hole.(3) A reduction in the CH_4 concentration and an increase in the CO_2 concentration are obvious during the displacement process.(4) The e ective influencing radius of injecting CO_2 with heat increases with the increase in time and pressure. The relationship between the e ective influencing radius and the injection time of CO_2 has a power exponential function, and there is a linear relationship between the functional coe cient and the injection pressure of CO_2. This numerical simulation study on enhancing CBM recovery by injecting CO_2 with heat can further promote the implementation of CO_2-ECBM project in deep coal seams.

Geomechanics model test and numerical simulation of 2G-NPR bolt support effect in an active fault tunnel

Active faults are a common adverse geological phenomenon that can occur during tunnel excavation and has a very negative impact on the construction and operation of the tunnel.In this paper,the grade IV rock surrounding the cross-fault tunnel with poor geological conditions has been chosen for the study.The support capacity of 2^(nd) Generation-Negative Poisson’s Ratio(2G-NPR)bolt in an active fault tunnel has been carried out on the basis of relevant results obtained from the geomechanical model test and numerical investigations of failure model for existing unsupported fault tunnel.The investigation shows that surrounding rock of the tunnel is prone to shear deformation and crack formation along the fault,as a result,the rock mass on the upper part of the fault slips as a whole.Furthermore,small-scale deformation and loss of blocks are observed around the tunnel;however,the 2G-NPR bolt support is found to be helpful in keeping the overall tunnel intact without any damage and instability.Due to the blocking effect of fault,the stress of the surrounding rock on the upper and lower parts of the fault is significantly different,and the stress at the left shoulder of the tunnel is greater than that at the right shoulder.The asymmetrical arrangement of 2G-NPR bolts can effectively control the asymmetric deformation and instability of the surrounding rock.The present numerical scheme is in good agreement with the model test results,and can reasonably reflect the stress and displacement characteristics of the surrounding rock of the tunnel.In comparison to unsupported and ordinary PR(Poisson’s Ratio)bolt support,2G-NPR bolt can effectively limit the fault slip and control the stability of the surrounding rock of the fault tunnel.The research findings may serve as a guideline for the use of 2G-NPR bolts in fault tunnel support engineering.

Using TOUGH2 numerical simulation to analyse the geothermal formation in Guide basin,China

The Guide sedimentary basin is located in the northeastern part of Qinghai-Xizang Plateau,which is rich in geothermal resources.However,exploitation of the geothermal resources has so far been limited,because of limited understanding of the resources quantity and storage gained from scientific researches.In this study,using a typical cross section across the basin and taking into account its geothermal and geological conditions,a new waterheat coupled model was built and associated modelling was done by the software TOUGH2.During modelling process,the accuracy and applicability of the model was confirmed through the calibration of relevant parameters for modelling the heat and water transport and the formation of geothermal reservoir across the basin,with particular focus on the Neogene geothermal field.Results show that the groundwater that flows from the basin margins to the center is heated by the Neogene and Paleogene sedimentary rocks with high geothermal gradients.Since the east-west extending fault F1 is conductive,it acts as preferential flow paths which on one hand provide additional and rapid flows to the thermal reservoir;and on the other hand,cool down the thermal water to a certain extent due to the infiltration of shallower water sources in the vicinity of the fault.Furthermore,the estimated geothermal resources quantity is close to that of previous studies.In comparison with the Paleogene rock formations,the Neogene geothermal reservoir shows a better nature in terms of water content,aquifer permeability and resources exploitability,although the resource quantity of the Paleogene reservoir is considerable.

Debris flow simulation 2D(DFS 2D):Numerical modelling of debris flows and calibration of friction parameters

Debris flows are rapid mass movements with a mixture of rock,soil and water.High-intensity rainfall events have triggered multiple debris flows around the globe,making it an important concern from the disaster management perspective.This study presents a numerical model called debris flow simulation 2D(DFS 2D)and applicability of the proposed model is investigated through the values of the model parameters used for the reproduction of an occurred debris flow at Yindongzi gully in China on 13 August 2010.The model can be used to simulate debris flows using three different rheologies and has a userfriendly interface for providing the inputs.Using DFS 2D,flow parameters can be estimated with respect to space and time.The values of the flow resistance parameters of model,dry-Coulomb and turbulent friction,were calibrated through the back analysis and the values obtained are 0.1 and 1000 m/s^(2),respectively.Two new methods of calibration are proposed in this study,considering the crosssectional area of flow and topographical changes induced by the debris flow.The proposed methods of calibration provide an effective solution to the cumulative errors induced by coarse-resolution digital elevation models(DEMs)in numerical modelling of debris flows.The statistical indices such as Willmott's index of agreement,mean-absolute-error,and normalized-root-mean-square-error of the calibrated model are 0.5,1.02 and 1.44,respectively.The comparison between simulated and observed values of topographic changes indicates that DFS 2D provides satisfactory results and can be used for dynamic modelling of debris flows.

AN INTEGRATION METHOD WITH FITTING CUBIC SPLINE FUNCTIONS TO A NUMERICAL MODEL OF 2ND-ORDER SPACE-TIME DIFFERENTIAL REMAINDER——FOR AN IDEAL GLOBAL SIMULATION CASE WITH PRIMITIVE ATMOSPHERIC EQUATIONS

In this paper,the forecasting equations of a 2nd-order space-time differential remainder are deduced from the Navier-Stokes primitive equations and Eulerian operator by Taylor-series expansion.Here we introduce a cubic spline numerical model(Spline Model for short),which is with a quasi-Lagrangian time-split integration scheme of fitting cubic spline/bicubic surface to all physical variable fields in the atmospheric equations on spherical discrete latitude-longitude mesh.A new algorithm of"fitting cubic spline—time step integration—fitting cubic spline—……"is developed to determine their first-and2nd-order derivatives and their upstream points for time discrete integral to the governing equations in Spline Model.And the cubic spline function and its mathematical polarities are also discussed to understand the Spline Model’s mathematical foundation of numerical analysis.It is pointed out that the Spline Model has mathematical laws of"convergence"of the cubic spline functions contracting to the original functions as well as its 1st-order and 2nd-order derivatives.The"optimality"of the 2nd-order derivative of the cubic spline functions is optimal approximation to that of the original functions.In addition,a Hermite bicubic patch is equivalent to operate on a grid for a 2nd-order derivative variable field.Besides,the slopes and curvatures of a central difference are identified respectively,with a smoothing coefficient of 1/3,three-point smoothing of that of a cubic spline.Then the slopes and curvatures of a central difference are calculated from the smoothing coefficient 1/3 and three-point smoothing of that of a cubic spline,respectively.Furthermore,a global simulation case of adiabatic,non-frictional and"incompressible"model atmosphere is shown with the quasi-Lagrangian time integration by using a global Spline Model,whose initial condition comes from the NCEP reanalysis data,along with quasi-uniform latitude-longitude grids and the so-called"shallow atmosphere"Navier-Stokes primitive equations in the spherical coordinates.The Spline Model,which adopted the Navier-Stokes primitive equations and quasi-Lagrangian time-split integration scheme,provides an initial ideal case of global atmospheric circulation.In addition,considering the essentially non-linear atmospheric motions,the Spline Model could judge reasonably well simple points of any smoothed variable field according to its fitting spline curvatures that must conform to its physical interpretation.

Numerical Simulation of the Temperature Field in Sintering of TiB_2-BN by SPS

The temperature field during the whole process of SPS sintering TiB2 -BN sample was numerically simulated based on the simplification of the temperature rising process. The result shows that the highest temperature is found in the punch and the heat flows from punch to sample and die in SPS sintering system. In the radial direction, the center temperature of the sample is much higher than the circumference temperature of the sampie. In the axial direction, the center temperature of the sample is lower than the border temperature of the sampie. The temperature difference in the sample is growing bigger in the sintering process and reaches the maximum at the end of the heating-up process.

3-D Finite-Element Numerical Simulation of Centrifugal Induction Electrosalg Casting Solidification Process

A 3-D finite-element numerical simulation model of temperature field for CIESC casting solidification process was developed with the aid of ANSYS software and a series of corresponding experiments were made. The results showed that the good agreement was obtained between the numerical simulation and the experiments. Based on the numerical simulation results, the characteristics of temperature distribution in the castings during CIESC solidification process were analyzed and summarized. According to the G/R-1/2 method and numerical simulation results, there is no any shrinkage defect in the CIESC casting and structure or casting is fine and compact.

Aseismic ridge subduction and flat subduction:Insights from three-dimensional numerical models

Flat subduction can significantly influence the distribution of volcanism,stress state,and surface topography of the overriding plate.However,the mechanisms for inducing flat subduction remain controversial.Previous two-dimensional(2-D)numerical models and laboratory analogue models suggested that a buoyant impactor(aseismic ridge,oceanic plateau,or the like)may induce flat subduction.However,three-dimensional(3-D)systematic studies on the relationship between flat subduction and buoyant blocks are still lacking.Here,we use a 3-D numerical model to investigate the influence of the aseismic ridge,especially its width(which is difficult to consider in 2-D numerical models),on the formation of flat subduction.Our model results suggest that the aseismic ridge needs to be wide and thick enough to induce flat subduction,a condition that is difficult to satisfy on the Earth.We also find that the subduction of an aseismic ridge parallel to the trench or a double aseismic ridge normal to the trench has a similar effect on super-wide aseismic ridge subduction in terms of causing flat subduction,which can explain the flat subduction observed beneath regions such as Chile and Peru.

Experimental and numerical studies of Ca(OH)_(2)/CaO dehydration process in a fixed-bed reactor for thermochemical energy storage

The Ca(OH)_(2)/CaO thermochemical energy storage(TCES)system based on calcium looping has received extensive attention owing to its high energy storage density,prolonged energy storage time,and environmental friendliness.The heat storage process of the Ca(OH)_(2)/CaO TCES system in a mixed heating reactor was evaluated in this study,by employing a combination of direct and indirect heating modes.The dehydration process was studied experimentally,and a numerical model was established and verified based on the experimental results.The dehydration behavior of 500 g of Ca(OH)_(2) powder was investigated in a fixed-bed reactor with mixed heating.The experimental and simulation results indicated that mixed heating causes combined centripetal and horizontal propulsion.Heat input is the main limiting factor in the heat storage process,because the radial advance of the reaction is hindered by the low thermal conductivity of the solid reactant particles.Heat transmission partitions were added to enhance the performance of the reactor.The performance of the modified reactor was compared with that of a conventional reactor.The radial heat transmission partitions in the modified reactor effectively enhance the energy storage rate and reduce the reaction time by 59.5%compared with the reactor without partitions.