Temperature Distributions in Single Cell of Polymer Electrolyte Fuel Cell Simulated by an 1D Multi-plate Heat-Transfer Model and a 3D Numerical Simulation Model

The purpose of this study is to verify an 1D multi-plate heat-transfer model estimating the temperature distribution on the interface between polymer electrolyte membrane and catalyst layer at cathode in single cell of polymer electrolyte fuel cell, which is named as reaction surface in this study, with a 3D numerical simulation model solving many governing equations on the coupling phenomena in the cell. The results from both models/simulations agreed well. The effects of initial operation temperature, flow rate, and relative humidity of supply gas on temperature distribution on the reaction surface were also investigated. It was found in both 1D and 3D simulations that, the temperature rise （i.e., Treact-Tini） of the reaction surface from initial operation temperature at 70℃ was higher than that at 80℃ irrespective of flow rate of supply gas. The effect of relative humidity of supply gas on Treact- Tini near the inlet of the cell was small. Compared to the previous studies conducted under the similar operation conditions, the Treact - Tini calculated by 1D multi-plate heat-transfer model in this study as well as numerical simulation using 3D model was reasonable.

Numerical simulation of protection range in exploiting the upper protective layer with a bow pseudo-incline technique

The developing processes of stress and deformation fields of a protected layer after mining an upper-protective layer with a bow pseudo-incline technique were simulated to locate the protection region. The pressure relief of the protected layer was analyzed after mining the upper-protective layer. The pressure relief angle along the strike and incline were located according to the roles of protection of the deformation and stress pressure-relief of the protective layer after mining. This results show that the upper-protective layer with the bow pseudo-incline technique have an upper and downside pressure relief angle of 85 and 68 degrees respectively; the distribution of strike pressure relief angles along the pseudo-incline working face is uneven and their values range from 38.3 to 51 degrees. The pressure relief angle of the inclined middle location was the largest. The distribution of the protection region of the upper-protective layer with the bow pseudo-incline teelmique located by practical tests and numerical simulation is essentially consistent, compared with the results obtained by these methods.

Study of 3-D Numerical Simulation for Gas Transfer in the Goaf of the Coal Mining

In order to simulate field distribution rules,mathematical models for 3-D air flows and gas transfer in the goaf of the coal mining are established,based on theories of permeability and dynamic dispersion through porous media. A gas dispersion equation in a 3-D field is calculated by use of numerical method on a weighted upstream multi-element balance. Based on data of an example with a U type ventilation mode,surface charts of air pressure distribution and gas concentration are drawn by Graphtool software. Finally,a comparison between actually measured results in the model test and the numerical simulation results is made to proves the numerical implementation feasible.

Numerical and Experimental Study of the 3D Effect on Connecting Arm of Vertical Axis Tidal Current Turbine

Vertical axis tidal current turbine is a promising device to extract energy from ocean current. One of the important components of the turbine is the connecting arm, which can bring about a significant effect on the pressure distribution along the span of the turbine blade, herein we call it 3D effect. However, so far the effect is rarely reported in the research, moreover, in numerical simulation. In the present study, a 3D numerical model of the turbine with the connecting arm was developed by using FLUENT software compiling the UDF（User Defined Function） command. The simulation results show that the pressure distribution along the span of blade with the connecting arm model is significantly different from those without the connecting arm. To facilitate the validation of numerical model, the laboratory experiment has been carried out by using three different types of NACA aerofoil connecting arm and circle section connecting arm. And results show that the turbine with NACA0012 connecting arm has the best start-up performance which is 0.346 m/s and the peak point of power conversion coefficient is around 0.33. A further study has been performed and a conclusion is drawn that the aerofoil and thickness of connecting arm are the most important factors on the power conversion coefficient of the vertical axis tidal current turbine.

Parallel numerical simulations for quantized vortices in Bose-Einstein condensates

We employ the parallel computing technology to study numerically the three-dimensional structure of quantized vortices of Bose-Einstein condensates, For anisotropic cases, the bending process of vortices is described in detail by the decrease of Gross-Pitaevskii energy. A completely straight vortex and the steady and symmetrical multiple-vortex configurations are obtained. We analyse the effect of initial conditions and angular velocity on the number and shape of vortices.

Numerical simulation of characteristics of semidiurnal tidal waves in sea region around Taiwan

POM97, an oceanic model, has been used for the first time to the numerical study on the tidal waves of the as Regions around Taiwan. In this paper, we have got the result that the semidiurnal tidal waves of these area mainly are the co-operating tides which come from the south of 23'N of the western Pacific. Those semidiurnal tidal waves affecting the Taiwan Straits come respectively from the south and the north entrance of the channel, and the north tidal wave is stronger than the south one. The strongest tidal field is the area from the Meizhou Bay to the Xinhua Bay along the coast of Fujian Province, where the biggest amplitude of the M2 partial tide can reach 240 cm. The strongest tidal cur- rent fields lie in the Penghu watercourse, where the maximum velocity of the M2 partial tide can arrive at 196 m/s. In the horizontal structure of the tidal currets, we have found that there is a stream dot in the north of the channel, besides, there still exist four new ones. As for the vertical structure, it mainly is biassed to the right at the surface, and to the left near the bottom layer.

Numerical modeling of 3D DC resistivity method in the mixed space-wavenumber domain

Forward modeling is the basis of inversion imaging and quantitative interpretation for DC resistivity exploration.Currently,a numerical model of the DC resistivity method must be finely divided to obtain a highly accurate solution under complex conditions,resulting in a long calculation time and large storage.Therefore,we propose a 3D numerical simulation method in a mixed space-wavenumber domain to overcome this challenge.The partial differential equation about abnormal potential is transformed into many independent ordinary differential equations with different wavenumbers using a 2D Fourier transform along the x axis and y axis direction.In this way,a large-scale 3D numerical simulation problem is decomposed into several 1D numerical simulation problems,which significantly reduces the computational and storage requirements.In addition,these ordinary 1D differential equations with different wavenumbers are independent of each other and high parallelelism of the algorithm.They are solved using a finite-element algorithm combined with a chasing method,and the obtained solution is modified using a contraction operator.In this method,the vertical direction is reserved as the spatial domain,then grid size can be determined flexibly based on the underground current density distribution,which considers the solution accuracy and calculation efficiency.In addition,for the first time,we use the contraction operator in the integral equation method to iterate the algorithm.The algorithm takes advantage of the high efficiency of the standard Fourier transform and chasing method,as well as the fast convergence of the contraction operator.We verified the accuracy of the algorithm and the convergence of the contraction operator.Compared with a volume integral method and goal-oriented adaptive finite-element method,the proposed algorithm has lower memory requirements and high computational efficiency,making it suitable for calculating a model with large-scale nodes.Moreover,different examples are used to verify the high adaptability and parallelism of the proposed algorithm.The findings show that the 3D numerical simulation method of DC resistivity method in a mixed space-wavenumber domain is highly efficient,precise,and parallel.

Experimental study on the movement law of overlying rock non-pillar coal overhead mining

Pillarless coal mining technology is a new practical technology.Based on the compensating mechanical behavior of the Negative Poisson’s Ratio(NPR)anchor cable on the roof,the roadway was successfully retained by the top cutting and pressure relief technology.This study utilizes the Digital Speckle Monitoring(DIC monitoring),stress-strain monitoring,and infrared thermal imaging systems to conduct physical model experiment of similar materials from the displacement,stress-strain,and temperature fields to investigate in depth the fracture change law of the overlying rock.In addition,it uses FLAC3D numerical simulation to invert the surface displacement settlement.The results show that the non-pillar overhead mining under the 110 mining method has little influence on the rock crack in the middle of the coal seam,and the crack development area is mainly concentrated in the overlying rock mass of the upward coal seam.The compensatory mechanical behavior of NPR anchor cable and the dilatation characteristics of rock mass have a good effect of retaining roadway along goaf,and can also reduce surface settlement.The 110 mining method provides a scientific basis for ecological environment protection and the development of other kilometer deep soft rock high ground stress underground projects.

Dynamic response analysis of liquefiable ground due to sinusoidal waves of different frequencies of shield construction

Vibration induced by shield construction can lead to liquefaction of saturated sand.Based on FLAC3D software,a numerical model of tunnel excavation is established and sinusoidal velocity loads with different frequencies are applied to the excavation face.The pattern of the excess pore pressure ratio with frequency,as well as the dynamic response of soil mass under different frequency loads before excavation,is analyzed.When the velocity sinusoidal wave acts on the excavation surface of the shield tunnel with a single sand layer,soil liquefaction occurs.However,the ranges and locations of soil liquefaction are different at different frequencies,which proves that the vibration frequency influences the liquefaction location of the stratum.For sand-clay composite strata with liquefiable layers,the influence of frequency on the liquefaction range is different from that of a single stratum.In the frequency range of 5-30 Hz,the liquefaction area and surface subsidence decrease with an increase in vibration frequency.The research results in this study can be used as a reference in engineering practice for tunneling liquefiable strata with a shield tunneling machine.

Development and Properties of Glass Fiber Reinforced Plastics Geogrid

Glass fiber reinforced plastics geogrid has a wide application in the field of soil reinforcement because of its high strength, good toughness, and resistance to environmental stress, creep resistance and strong stability. In order to get high-powered glass fiber reinforced plastics geogrid and its mechanical characteristics, the properties and physical mechanical index of geogrid have been got through the study of its raw material, production process and important quality index. The analysis and study have been made to the geogrid＇s mechanical properties with loading speed, three-axial compression, temperature tensile test and FLAC3D numerical simulation, thus obtain the mechanical parameters of its displacement time curve, breaking strength and elongation at break. Some conclusions can be drawn as follows： （a） Using glass fiber materials, knurling and coated projection process, the f^acture strength and corrosion resistance of geogrid are greatly improved and the interlocking bite capability of soil is enhanced. （b） The fracture strength of geogrid is related to temperature and loading rate. When the surrounding rock pressure is fixed, the strength and anti-deformation ability of reinforced soil are significantly enhanced with increasing reinforced layers. （c） The pullout test shows the positive correlation between geogrid displacement and action time. （d） As a new reinforced material, the glass fiber reinforced plastics geogrid is not mature enough in theoretical research and practical experience, so it has become an urgent problem both in theoretical study and practical innovation.

Study on Wave Field Characteristics and Imaging of Collapse Column in Three-Dimensional Detection with Love Channel Wave Reflected outside the Working Face

The safety accidents caused by collapse column water diversion occur frequently, which has great hidden danger to the safety production of coal mine. Limited by the space of underground, the detection of collapse column on the outside of working face has been a difficult problem. Based on this, numerical simulation and imaging research were carried out in this paper. The results indicate that when a seismic source near the roadway is excited, a part of seismic wave propagates along the roadway direction, namely direct P-wave, direct S-wave and direct Love channel wave. When the body waves and Love channel wave propagating to the outside of working face meet the interface of collapse column, the reflected Love channel wave and reflected body waves are generated. Reflection body waves and direct waves are mixed in time domain, which is difficult to identify in seismic records, while reflected Love channel wave whose amplitude is relatively strong. The reflected Love channel wave which has a large interval from other wave trains in the time domain is easily recognizable in seismic record, which makes it suitable for advanced detection of collapse column. The signal-to-noise ratio of X component is higher than that of Y component and Z component. According to the seismic records, polarization filtering was carried out to enhance the effective wave, which removed the interference waves, and the signal was migrated to get the position parameters of collapse column interface, which was basically consistent with the model position.

Support technology for mine roadways in extreme weakly cemented strata and its application

For the engineering geology conditions of bad mine roadway roof and floor lithology in extremely weak cemented strata, the best section shape of the roadway is determined from the study of tunnel surrounding rock displacement, plastic zone and stress distribution in rectangular, circle arch and arch wall sections, respectively. Based on the mining depth and thickness of the coal seam, roadway support technology solutions with different buried depth and thickness of coal seam are proposed. Support schemes are amended and optimized in time through monitoring data of the deformation of roadway, roof separation, l-beam bracket, bolt and anchor cable force to ensure the long-term stability and security of the roadway surrounding rock and support structure. The monitoring results show that mine roadway support schemes for different buried depth and section can be adapted to the characteristics of ground pressure and deformation of the surrounding rock in different depth well, effectively control the roadway surrounding rock deformation and the floor heave and guarantee the safety of construction and basic stability of surrounding rock and support structure.

Beach stabilization by the laying of a drainage layer

In the present study, the stabilization effects of the laying of a drainage layer in a sandy beach were analyzed. The numerical model developed by Hur et al. （2011） was used to examine the characteristics of wave-sandy beach interaction with/without a drainage layer. This numericai model directly simulates Wave-Structure-Seabed/Sandy （WASS） beach interaction by consid- ering flow through a porous medium with inertia, laminar and turbulent resistances, and determines the eddy viscosity with the Large Eddy Simulation （LES） turbulent model in a 3-D wave field （LES-WASS-3D）. Based on the numerical results, it ap- pears that the installation of a drainage layer in the beach results in a decrease in both the mean groundwater level around the foreshore and the flow moving offshore on the beach surface. Spatial distributions of the mean groundwater level and wave setup around the foreshore with various arrangements of the drainage layer and conditions of the incident wave are also discussed in this paper.