Numerical analysis of matrix swelling and its effect on microstructure of digital coal and its associated permeability during CO_(2)-ECBM process based on X-ray CT data

Matrix swelling effect will cause the change of microstructure of coal reservoir and its permeability,which is the key factor affecting the engineering effect of CO_(2)-ECBM technology.The Sihe and Yuwu collieries are taken as research objects.Firstly,visualization reconstruction of coal reservoir is realized.Secondly,the evolution of the pore/fracture structures under different swelling contents is discussed.Then,the influence of matrix phase with different swelling contents on permeability is discussed.Finally,the mechanism of swelling effect during the CO_(2)-ECBM process is further discussed.The results show that the intra-matrix pores and matrix-edge fractures are the focus of this study,and the contacting area between matrix and pore/fracture is the core area of matrix swelling.The number of matrix particles decreases with the increase of size,and the distribution of which is isolated with small size and interconnected with large size.The swelling effect of matrix particles with larger size has a great influence on the pore/fracture structures.The number of connected pores/fractures is limited and only interconnected in a certain direction.With the increase of matrix swelling content,the number,porosity,width,fractal dimension,surface area and volume of pores/fractures decrease,and their negative contribution to absolute permeability increases from 0.368% to 0.633% and 0.868%-1.404%,respectively.With the increase of swelling content,the number of intra-matrix pores gradually decreases and the pore radius becomes shorter during the CO_(2)-ECBM process.The matrix continuously expands to the connected fractures,and the width of connected fractures gradually shorten.Under the influence of matrix swelling,the bending degree of fluid flow increases gradually,so the resistance of fluid migration increases and the permeability gradually decreases.This study shows that the matrix swelling effect is the key factor affecting CBM recovery,and the application of this effect in CO_(2)-ECBM process can be discussed.

Numerical Analysis of a Spiral-groove Dry-gas Seal Considering Micro-scale Effects

A dry-gas seal system is a non-contact seal technology that is widely used in different industrial applications.Spiral-groove dry-gas seal utilizes fluid dynamic pressure effects to realize the seal and lubrication processes,while forming a high pressure gas film between two sealing faces due to the deceleration of the gas pumped in or out.There is little research into the effects and the influence on seal performance,if the grooves and the gas film are at the micro-scale.This paper investigates the micro-scale effects on spiral-groove dry-gas seal performance in a numerical solution of a corrected Reynolds equation.The Reynolds equation is discretized by means of the finite difference method with the second order scheme and solved by the successive-over-relaxation（SOR） iterative method.The Knudsen number of the flow in the sealing gas film is changed from 0.005 to 0.120 with a variation of film depth and sealing pressure.The numerical results show that the average pressure in the gas film and the sealed gas leakage increase due to micro-scale effects.The open force is enlarged,while the gas film stiffness is significantly decreased due to micro-scale effects.The friction torque and power consumption remain constant,even in low sealing pressure and spin speed conditions.In this paper,the seal performance at different rotor face spin speeds is also described.The proposed research clarifies the micro-scale effects in a spiral-groove dry-gas seal and their influence on seal performance,which is expected to be useful for the improvement of the design of dry-gas seal systems operating in the slip flow regime.

Numerical analysis and geophysical monitoring for stability assessment of the Northwest tailings dam at Westwood Mine

The Westwood Mine aims to reuse the tailings storage facility #1(TSF #1) for solid waste storage, but,downstream of the Northwest dike is considered critical in terms of stability. This paper uses numerical modeling along with geophysical monitoring for assessing the Northwest dike stability during the restoration phase. The impact of waste rock deposition in the upstream TSF #1 is considered. The geophysical monitoring is based on electrical resistivity methods and was used to investigate the internal structure of the dike embankment in different deposition stages. The numerical simulations were performed with SLOPE/W code. The results show a factor of safety well above the minimum recommended value of 1.5. Geophysical monitoring revealed a vertical variation in the electrical resistivity across the dike, which indicates a multilayer structure of the embankment. Without any current in situ data, the geophysical monitoring helped estimating the nature of the materials used and the internal structure of the embankment. These interpretations were validated by geological observation of geotechnical log of the embankment. Based on this study, it is recommended that the water polishing pond be partly filled before waste rock is deposited in TSF #1. In addition, to ensure the stability of the dike, the piezometric head monitoring prior to and during waste rock deposition is recommended.

Shaking table test and numerical analysis of a 1:12 scale model of a special concentrically braced steel frame with pinned connections

This paper describes shaking table tests of a 1：12 scale model of a special concentrically braced steel frame with pinned connections, which was fabricated according to a one-bay braced frame selected from a typical main factory building of a large thermal power plant. In order to investigate the seismic performance of this type of structure, several ground motion accelerations with different levels for seismic intensity Ⅷ, based on the Chinese Code for Seismic Design of Buildings, were selected to excite the model. The results show that the design methods of the members and the connections are adequate and that the structural system will perform well in regions of high seismicity. In addition to the tests, numerical simulations were also conducted and the results showed good agreement with the test results. Thus, the numerical model is shown to be accurate and the beam element can be used to model this structural system.

Numerical analysis of surface subsidence in asymmetric parallel highway tunnels

Tunnelling related hazards are very common in the Himalayan terrain and a number of such instances have been reported. Several twin tunnels are being planned for transportation purposes which will require good understanding for prediction of tunnel deformation and surface settlement during the engineering life of the structure. The deformational behaviour, design of sequential excavation and support of any jointed rock mass are challenging during underground construction. We have raised several commonly assumed issues while performing stability analysis of underground opening at shallow depth. For this purpose, Kainchi-mod Nerchowck twin tunnels(Himachal Pradesh, India) are taken for in-depth analysis of the stability of two asymmetric tunnels to address the influence of topography, twin tunnel dimension and geometry. The host rock encountered during excavation is composed mainly of moderately to highly jointed grey sandstone, maroon sandstone and siltstones. In contrast to equidimensional tunnels where the maximum subsidence is observed vertically above the centreline of the tunnel, the result from the present study shows shifting of the maximum subsidence away from the tunnel centreline. The maximum subsidence of 0.99 mm is observed at 4.54 m left to the escape tunnel centreline whereas the maximum subsidence of 3.14 mm is observed at 8.89 m right to the main tunnel centreline. This shifting clearly indicates the influence of undulating topography and inequidimensional noncircular tunnel.

Numerical analysis of dynamic variation of weld pool geometry in fully-penetrated TIG welding

A mathematical model is developed for numerical analysis of thermal process in TIG welding with a moving arc, which is considered the double-elliptic distribution for both arc heat flux and arc pressure. An adjusting factor is introduced into the expression of arc pressure. The domain within which the arc heat flux is distributed non-symmetrically due to arc moving is selected appropriately, and three conditions for the domain to meet are described. The latent heat is taken into consideration by liquid fraction method. The dynamic development of weld pool geometry during TIG welding is analyzed numerically, and the effect of arc moving on the weld pool geometry is discussed. The experimental results show that the numerical analysis accuracy is obviously improved through taking the above-mentioned measures.

Numerical Analysis of Interaction Between Pile-Supported Pier and Bank Slope

Two and three-dimensional finite element analysis programs for pile-soil interaction are compiled. Duncan-Chang's Model is used. The construction sequence of the pier is modeled. The pile-soil interface element is used. The influence of the combination type of piles on the deformation of bank slope and pile behaviour is analyzed. Different designs of a pile-supported pier are compared thoroughly. Calculation results show that the stresses and displacements of the pile are directly related to the distance from the bank slope and the direction of inclination. An inclined prop pile set in the rear platform would remarkably reduce the stresses of piles and the displacement of the pier.

Numerical Analysis of Hydrodynamic Pressure Induced by Fluid-Solid Impact

As a further development of the authors' work (Huang and Qian, 1993), in this paper a new numerical method based on the time domain boundary element technique is proposed for solving fluid-solid coupling problems, in which a rigid body impacts normally on the calm surface of a half-space fluid. A fundamental solution to the half-space potential flow problem is first derived with the method of images. Then, an equivalent boundary integral equation in the Laplace transform domain is established by means of Green's second identity. Through the inverse Laplace transform and discretization in both time and boundary of the fluid region, the numerical calculation for the problem under consideration has been carried out. Several examples demonstrate that the present method is more efficient than existing ones, from which it is also seen that the shape of the impacting body has a considerable effect on the total impact force.

Experimental and Numerical Analysis of Restrained Early Age Cracking based on Electrical Resistivity Using Eccentric Sample

Crack potential and hydration processes of the cement pastes were monitored using an upto-date eccentric steel cracking frame(ESCF), associated with the non-contact electrical resistivity apparatus, independently. The objective of employing the ESCF is to give a new method determining cracks of concrete at early age. The findings indicate that the lowest water-cement ratio paste reveals highest resistivity values, compasses an earlier inflection point and obtained higher stress. The eccentric restrained cracking test exhibited that lower water-cement ratio paste cracked at the earliest time, accordingly confirms cracking tendency is the highest. Tensile strength test and stresses utilizing ABAQUS simulation was performed. The crack initiation ages obtained are consistent with the experimental program results, which indicates that ABAQUS numerical analysis can well be utilized to predict the crack tendency of cement.

NUMERICAL ANALYSIS AND EXPERI-MENT OF UNSTEADY THERMAL FIELD OF ROTOR PLATE FOR EDDY CURRENT RETARDER

The physical model based on heat transfer theory and virtual boundary method for analyzing unsteady thermal field of rotor plate for eddy current retarder used in automobile is established and boundary conditions are also defined. The finite element governing equation is derived by Galerkin method. The time differential item is discrete based on Galerkin format that is stable at any condition. And a new style of varying time step method is used in iteration process. The thermal field on the rotor plate at the radial and axle directions is analyzed and varying temperature at appointed points on two side-surfaces is measured. The testing and analytical data are uniform approximately. Finite element method can be used for estimating thermal field of the rotor plate at initial design stage of eddy current retarder.

NUMERICAL ANALYSIS ON THREE-DIMENSIONAL FLOW FIELD OF TURBINE IN TORQUE CONVERTER

Three-dimensional flow field of turbine in torque converter is simulated by numerical calculation in order to improve the performance of torque converter. Calculation model of a torque converter is presented based on the mixing-plane technology. In the calculation of flow field, the 3D N-S equations are separated by finite-volume method and solved by semi-implicit method for pressure-linked equations（SIMPLE）. Based on flow field calculation, the flow field of turbine is simulated. The velocity and pressure in the flow field of turbine are analyzed. The external performance of the torque converter is also calculated. Results of flow simulation show that there are secondary flow, off flow and velocity gradient in turbine passage. The validity of numerical simulation is verified by comparing the results of numerical simulation with experiment data.

Numerical Analysis of Deteriorated Sub-sea Pipelines under Environmental Loads

The significant point is the bidirectional interaction technique in FSI analysis while investigating subsea corrosion effect. By this way, pipe environment is accurately modelled and fluid effects are also considered. The effect of external corrosion defects on structural behaviour of a pipeline is studied by creating a nonlinear numerical model based on the finite element method according to ABAQUS analysis program. Corrosion losses of sections are obtained from experimental results and applied to the model. Numerical model is formed by a span of sub-sea pipeline that is subjected to environmental loads. Seismic and wind-generated irregular wave loads are considered as environmental loads. Irregular wave is represented with equivalent eight regular waves via FFT. The pipe is modelled according to two different types which are non-corroded（intact） and corroded（deteriorated） to demonstrate corrosion effects on it. The visible type of corrosion in marine environment is named ‘pitting＇ corrosion, in which the material loss is locally interpenetrated over the surface. By considering this situation, the corroded and non-corroded pipes are modelled as 3D solid elements. The main point is revealing how the subsea corrosion affects the structural behaviour of pipelines on the basis of implementation of experimental results to a model structure due to changes of stresses and displacement.

NUMERICAL ANALYSIS OF GASEOUS FLOW IN MICRO-CHANNELS

The algorithm of gaseous flow in bi-dimensional micro-channels is set up andthe corresponding program based on micro-flow theory is presented. Gaseous flow in micro-channels isnumerically analyzed and the pressure drop along the duct as well .as the velocity profile in themicro-channels is obtained. The numerical results agreed well with the experimental results in thereferences. Moreover, the effects of Kn, sigma_v and Re on the velocity profiles are analyzed. It isfound that for Kn>0.001, with increasing Kn number, the slip velocity on the wall boundaryincreases; the tangential momentum coefficient sigma_v affects the slip velocity greatly. The slipvelocity increases with decreasing a, In the slip flow regime and for low Re numbers, the slipvelocity is little influenced by the Re number.

Numerical analysis for the effcacy of nasal surgery in obstructive sleep apnea hypopnea syndrome

In the present study, we reconstructed upper airway and soft palate models of 3 obstructive sleep apnea-hypopnea syndrome（OSAHS） patients with nasal obstruction. The airflow distribution and movement of the soft palate before and after surgery were described by a numerical simulation method. The curative effect of nasal surgery was evaluated for the three patients with OSAHS. The degree of nasal obstruction in the 3 patients was improved after surgery. For 2 patients with mild OSAHS, the upper airway resistance and soft palate displacement were reduced after surgery. These changes contributed to the mitigation of respiratory airflow limitation. For the patient with severe OSAHS, the upper airway resistance and soft palate displacement increased after surgery, which aggravated the airway obstruction. The effcacy of nasal surgery for patients with OSAHS is determined by the degree of improvement in nasal obstruction and whether the effects on the pharynx are beneficial. Numerical simulation results are consistent with the polysomnogram（PSG） test results, chief complaints, and clinical findings, and can indirectly reflect the degree of nasal patency and improvement of snoring symptoms, and further,provide a theoretical basis to solve relevant clinical problems.

Numerical Analysis of Magnetic-Shielding Effectiveness for Magnetic Resonant Wireless Power Transfer System

Magnetic radiation phenomena appear inevitably in the magnetic-resonance wireless power transfer （MR-WPT） system, and regarding this problem the magnetic-shielding scheme is applied to improve the electromagnetic performance in engineering. In this study, the shielding effectiveness of a two-coil MR-WPT system for different material shields is analyzed in theory using Moser＇s formula and Schelkunoff＇s formula. On this basis a candidate magnetic-shielding scheme with a double-layer structure is determined, which has better shielding effectiveness and coils coupling coefficient. Finally, some finite element simulation results validate the correctness of the theoretical analysis, and the shielding effectiveness with the double-layer shield in maximum is 30？dB larger than the one with the single-layer case.

Numerical analysis and experimental studies of membrane wrinkles

Wrinkling patterns and the rules of wrinkling formation and expansion were analyzed in this paper. Direct disturbing method,a numerical method for membrane wrinkling analysis,was established based on finite element programme. The experimental results using photogrammetric procedure were compared with numerical analysis results. It is shown that the reasonable wrinkling deformation can be obtained by rational mesh density and timely removes the initial flaw during analysis. Wrinkles firstly appear on the left and right sides then extend to the membrane center,and the wrinkling amplitude increases with shear displacement. Results from both of the experiments and numerical analysis have a good agreement,which verifies the rationality of the numerical analysis method.

On the Characteristics of Coaxial-Type Microwave Excited Linear Plasma:a Simple Numerical Analysis

To unveil the characteristics and available propagation mechanism of coaxial-type microwave excited line-shape plasma, the effects of parameters including microwave power, working pressure, dielectric constant, and external magnetic field on the plasma distribution were numerically investigated by solving a coupled system of Maxwell＇s equations and continuity equations. Numerical results indicate that high microwave power, relatively high working pressure, low dielectric constant, and shaped magnetic field profiles will help produce a high-density and uniform plasma source. Exciting both ends by microwave contributed to the high-density and uni- form plasma source as well. Possible mechanisms were analyzed by using the polarization model of low temperature plasma. The generation and propagation processes of the line-shape plasma mainly depend on the interaction of three aspects, i.e. the transmitted part, penetration part and absorptive part of the electromagnetic field. The numerical results were qualitatively consistent with available experimental results from literature. More elaborate descriptions of the three aspects and corresponding interactions among them need to be investigated further to improve the properties of the line-shape plasma.

Numerical analysis of external magnetic field for mitigating backward flow momentum in weld pool

The external magnetic field is applied to mitigating backward flow jet of molten metal in weld pool so that humping bead may be suppressed during high speed gas metal arc welding（GMAW）. Therefore, the external magnetic field distribution in workpiece is critical to understand the interaction mechanisms of the external magnetic field with molten metal flow. In this study, the steady state external magnetic field induced by excitation device is numerically analyzed by using the the finite element software ANSYS and the three dimensional static magnetic scalar method. The distribution of external transverse magnetic field By in workpiece and arc area is calculated, and the influence of excitation current and air-gap distance on the distribution of transverse magnetic field By has been discussed. The magnetic field distribution in workpiece is measured by using a Tesla-Meter and compared with the simulated result. It is found that both are in good agreement.

The numerical analysis of thermal effect and weld shaping of DE-GMAW under composite heat source

The action of the composite heat source model in DE-GMAW has been carried out according to the characteristics of the DE-GMA W process, and the distribution of surface and body heat source was analyzed. The weld temperature field distribution has been derived from the appropriate boundary conditions and the thermal physical property parameters by COMSOL software. The effects of the positions of the surface and body heat source on the cross- sectional shape were investigated by studying the experimental welding thermal cycle curve. The simulated results are fully compliant with the experimental data under the same conditions. This illustrates that the composite heat source model is correct. It reflects the thermal mechanism of DE-GMA W process, and reveals the internal influence of weld shaping.

NUMERICAL ANALYSIS OF THREE-DIMENSIONAL FLUID FLOW AND HEAT TRANSFER IN TIG WELD POOL WITH FULL PENETRATION WU Chuansong,CAO Zhenning,WU Lin Harbin Institute of Technology,Harbin,China

A model is established to analyze three-dimensional fluid flow and heat transfer in TICweld pools with full penetration.It considers the deformation of the molten pool surfaceat the condition of full penetrated workpieees,takes the are pressure as the drivingforce of the pool surface deformation,and determines the surface configuration of weldpool based on the dynamic balance of arc pressure,pool gravity and surface tension atdeformed weld pool surface. The SIMPLER algorithm is used to calculate the fluid flowfield and temperature distribution in TIG weld pools of stainless steel workpieces.TIGwelding experiments are made to verify the validity of the model.It shows the calculatedresults by the model are in good agreement with experimental measurements. professor,Dept of Welding Engineering,Harbin Institute of Technology,Harbin 150006,China