受冲击海上桥梁的非线性动力响应特征数值模拟

在受冲击海上桥梁的非线性动力响应特征数值模拟中,使用传统动力响应特征数值模拟方法模拟的准确度较低,因此提出一种受冲击海上桥梁的非线性动力响应特征数值模拟方法,利用有限元模型对海上桥梁进行基础非线性动力响应分析,通过自由度方程组建立动力响应特征数值模拟模型,结合等分能量法实现动力响应特征数值模拟,为了验证该方法的模拟准确度,利用该方法与传统动力响应特征数值模拟方法同时进行桥梁的非线性动力响应特征数值模拟,其模拟准确度分别为92.4%,72.3%,65.4%,69.1%,通过对比实验可知该方法的模拟准确度最高,验证了该方法的优越性。

离散裂缝网络模型与有限元方法在页岩气藏数值模拟中的研究

目前对砂岩单孔隙介质油藏的数值模拟理论和方法比较成熟,由于页岩基质和裂缝之间显著的物性差异以及页岩气的吸附规律使得我们无法继续采用单孔隙介质来对其进行气藏数值模拟研究,对于复杂而致密的页岩气藏,其理论基础还不完善,调研国内外页岩气数值模拟现状,搞清楚气水在页岩气藏中的渗流机理,建立此流动机制的页岩气藏的控制方程,并建立页岩气藏离散裂缝网络数学模型,探索页岩气藏离散裂缝模型的离散化方法,使用有限元方法求解离散裂缝网络模型,通过离散裂缝网络数值模拟方法优化页岩气藏地层参数和裂缝参数,并作出产能评价,这些方面对页岩气藏离散裂缝网络数值模拟研究的突破具有十分重大的实际意义。

2018年大连海域船舶大气污染物排放特征及影响分析

船舶废气已成为大气污染物的重要来源之一。为掌握大连周边海域船舶大气污染物排放特征,以及大连海域船舶排放对大连市大气环境的影响,基于船舶自动识别系统监测数据,采用基于功率的动力法估算出2018年大连海域船舶排放清单,并利用WRF、SMOKE、CMAQ空气质量数值模拟系统,研究了1月、4月、7月、10月船舶排放对大连市大气环境的影响。研究结果表明,2018年大连海域船舶共排放SO_(2)7606.23 t、NO_(x) 30990.13 t、PM_(10)1212.02 t、PM_(2.5)969.58 t、CO 3339.10 t、HC 1414.63 t、CO_(2)2546299.67 t。其中,客船、散货船、油船和多用途船的排放占比较大,排放区域主要集中在大连湾、大窑湾附近及周边航道。受气候因素影响,不同月份的船舶排放对大连市大气环境的影响不同,其中,1月、4月和10月受影响较大的区域主要是旅顺口区和金州区,而7月受影响较大的区域主要是甘井子区。

预测蒸汽驱后期转接替方式时机的新方法

蒸汽驱开发后期易汽窜,且原油产量呈递减趋势,当汽窜严重时,生产油汽比下降,蒸汽驱开发的效果越来越差。为了取得最佳的采出程度与开发的经济性,确定蒸汽驱后期转接替方式的时机非常重要。通过分析蒸汽驱后期产油量递减曲线,总结出蒸汽驱开发的递减规律;通过对蒸汽驱注采特征方程的推导得出蒸汽驱的采收率表达式,两式联立便可以求出蒸汽驱转接替方式的时机。通过对齐40块某井组的实例分析得出,该井组现阶段虽已出现汽窜,但并不是转接替方式的最佳时机。计算结果表明,2017年1月转汽水交替较为合适。进一步通过数值模拟验证,当蒸汽驱至2017年1月时,蒸汽腔在地层中发育较为成熟,汽窜程度更加严重,此时该区块从蒸汽驱转其他的接替方式较为合适。数值模拟很好的验证了计算结果,能够较好的指导现场开发。

副高边缘一次突发性暴雨过程的成因分析

本文利用常规的探空资料、卫星云图和U-CAR/PSU的中尺度大气模式MM5,对2004年8月5日发生在山东半岛的一次中β尺度暴雨过程进行了天气学分析、数值模拟及动力特征分析。结果表明:西太平洋副高边缘加强西伸,高、低空西南风增大导致山东半岛纬向型切变线发展、辐合加强造成本次暴雨,暴雨发生在有利的水汽输送和水汽辐合的条件下,水汽的主要来源为长江中游和黄海。暴雨中心附近垂直动力结构特征为:450hPa以下为强烈的上升运动,低层正涡度、辐合,高层负涡度、辐散,这种结构有利于强对流天气的发生。

SIMULATION OF MICROSTRUCTURE AND RHEOGICAL PROPERTIES OF SURFACTANT SOLUTION

The Brownian dynamics （BD） simulation of a dilute surfactant solution is conducted in a steady shear flow. The rodlike micelle is assumed as a rigid rod composed of lined-up beads. A novel intercluster potential model is introduced for describing the interactions between, micelles. In the model, the Lennard-Jones and the soft-sphere potentials are used as inter-bead potentials for end-end and interior-interior beads, respectively. The micelles are combined at their ends to form a network structure at lower shear rates and are disconnected to become more and more parallel to the shear flow direction with increasing shear rate. The change of micellar microstructures with the variation of the shear rate results in shear thinning characteristics of the computed shear viscosities and first normal stress difference coefficients. The effects of surfactant solution concentration on the micellar structures and rheological properties are also investigated. Results show that the shear viscosities and the first normal stress difference coefficients increase with increasing the viscosity of the surfactant solution.

Experimental investigation into stress-relief characteristics with upward large height and upward mining under hard thick roof

According to geological conditions of No. 3 and No. 4 coal seams （namely A3 and B4） of the Pan＇er coal mine and the parameters of panels 11223, 11224, and 11124 with fully-mechanical coal mining, we built 2D similar material simulation and FLAC3D numerical simulation models to investigate the development of mining-induced stress and the extraction effect of pressure-relief gas with large height and upward mining. Based on a comprehensive analysis of experimental data and observations, we obtained the deformation and breakage characteristics of strata overlying the coal seam, the development patterns of the mining-induced stress and fracture, and the size of the stress-relief area. The stress-relief effect was investigated and analyzed in consideration with mining height and three thick hard strata. Because of the group of three hard thick strata located in the main roof and the residual stress of mined panel 11124, the deformation, breakage, mining-induced stress and fracture development, and the stress-relief coefficient were discontinuous and asymmetrical. The breakage angle of the overlying strata, and the compressive and expansive zones of coal deformation were mainly controlled by the number, thickness, and strength of the hard stratum. Compared with the value of breakage angle derived by the traditional empirical method, the experimental value was lower than the traditional results by 3°-4°below the hard thick strata group, and by 13°-19° above the hard thick strata group. The amount of gas extracted from floor drainage roadway of B4 over 17 months was variable and the amount of gas per month differed considerably, being much smaller when panel 11223 influenced the area of the three hard thick strata. Generally, the stress-relief zone of No. 4 coal seam was small under the influence of the hard thick strata located in the main roof, which played an important role in delaying the breakage time and increasing the breakage space. In this study we gained understanding of the stress-relief mechanism influenced by the hard thick roof. The research results and engineering practice show that the main roof of the multiple hard thick strata is a critical factor in the design of panel layout and roadways for integrated coal exploitation and gas extraction, provides a theoretical basis for safe and high-efficient mining of coal resources.

Forming condition of transient saturated zone and its distribution in residual slope under rainfall conditions

Rainfall, as one of the most significant factors triggering the residual soil slope failure, leads to not only the reduction of soil shear strength, but also the increase of soil weight and the decrease of matric suction as well. All these modifications in soil properties have important influence on the slope stability. The water infiltration and redistribution inside the slope are the preconditions of the slope stability under rainfall conditions. Based on the numerical simulation via finite element method, the water infiltration process under rainfall conditions was studied in the present work. The emphases are the formation, distribution and dissipation of transient saturated zone. As for the calculation parameters, the SWCC and the saturated permeability have been determined by pressure plate test and variable head test respectively. The entire process(formation, development, dissipation) of the transient saturated zone was studied in detail. The variations of volumetric water content, matric suction and hydraulic gradient inside the slope, and the eventually raise of groundwater table were characterized and discussed, too. The results show that the major cause of the formation of transient saturated zone is ascribed to the fact that the exudation velocity of rainwater on the wetting front is less than the infiltration velocity of rainfall; as a result, the water content of the soil increases. On the other hand, the formation and extension of transient saturated zone have a close relationship with rainfall intensity and duration. The results can help the geotechnical engineers for the deeper understanding of the failure of residual slope under rainfall condition. It is also suggested that the proper drainage system in the slope may be the cost-effective slope failure mitigation method.

Deformation mechanism and excavation process of large span intersection within deep soft rock roadway

The FLAC3D software was used to simulate and analyze numerically the displacement, stress and plastic zone distribu-tion characteristics of a large span intersection in a deep soft rock roadway after the surrounding rock was excavated. Our simula-tion results show that there are two kinds of dominating factors affecting roadway stability at points of intersection, one is the angle between horizontal stress and axial direction of the roadway and the other are the angles at the points of intersection. These results are based on a study we carried out as follows: first, we analyzed the failure mechanism of a large span intersection and then we built a mechanical model of a rock pillar at one of the points of intersection. At the end of this analysis, we obtained the failure characteristics of the critical parts on the large span intersection. Given these failure characteristics, we proposed a new supporting method, i.e., a Double-Bolt Control Technology (DBCT). By way of numerical simulation, DBCT can effectively control the deformation of the surrounding rock at the points of intersection in roadways.

Mechanical response features and failure process of soft surrounding rock around deeply buried three-centered arch tunnel

Due to the extreme complexity of mechanical response of soft surrounding rock(SR) around a tunnel under high geostatic stress conditions, the integration of physical and numerical modeling techniques was adopted. Based on the similarity theory, new composite-similar material was developed, which showed good agreement with the similarity relation and successfully simulated physico-mechanical properties(PMP) of deep buried soft rock. And the 800 mm×800 mm×200 mm physical model(PM) was conducted, in which the endoscopic camera technique was adopted to track the entire process of failure of the model all the time. The experimental results indicate that the deformation of SR around a underground cavern possessed the characteristics of development by stages and in delay, and the initial damage of SR could induce rapid failure in the later stage, and the whole process could be divided into three stages, including the localized extension of crack(the horizontal load(HL) was in the range of 130 k N to 170 k N, the vertical load(VL) was in the range of 119 k N to 153.8 k N), rapid crack coalescence(the HL was in the range of 170 k N to 210 k N, the VL was in the range of 153.8 k N to 182.5 k N) and residual strength(the HL was greater than 210 k N, the VL was greater than 182.5 k N). Under the high stress conditions, the phenomenon of deformation localization in the SR became serious and different space positions show different deformation characteristics. In order to further explore the deformation localization and progressive failure phenomenon of soft SR around the deeply buried tunnel, applying the analysis software of FLAC3 D three-dimensional explicit finite-difference method, based on the composite strain-softening model of Mohr-Coulomb shear failure and tensile failure, the calculation method of large deformation was adopted. Then, the comparative analysis between the PM experiment and numerical simulation of the three centered arch tunnels was implemented and the relationship of deformation localization and progressive failure of SR around a tunnel under high stress conditions was discussed.

Analytical Solution of Boundary Integral Equations for 2-D Steady Linear Wave Problems

Based on the Fourier transform, the analytical solution of boundary integral equations formulated for the complex velocity of a 2-D steady linear surface flow is derived. It has been found that before the radiation condition is imposed,free waves appear both far upstream and downstream. In order to cancel the free waves in far upstream regions, the eigensolution of a specific eigenvalue, which satisfies the homogeneous boundary integral equation, is found and superposed to the analytical solution. An example, a submerged vortex, is used to demonstrate the derived analytical solution. Furthermore,an analytical approach to imposing the radiation condition in the numerical solution of boundary integral equations for 2-D steady linear wave problems is proposed.

Some issues on atomic force microscopy based surface characterization

Influences of tip radius and sampling interval on applying atomic force microscopy （AFM） in quantitative surface evaluations are investigated by numerical simulations and experiments. Several evaluation parameters of surfaces ranging from amplitude to functional parameters are studied. Numerical and experimental results are in good agreements. The accuracy of estimating tip radius on random rough surface with Ganssian distribution of heights using a blind reconstruction method is also discussed theoretically. It is found that the accuracy is greatly depending on the ratio of actual tip radius to root-mean-square （rms） radius of curvature. To obtain an accurate estimation of tip radius under Gaussian rough surface, the ratio has to be larger than 3/2.

Application of high temperature heat pipe in hypersonic vehicles thermal protection

In order to develop further the application of high temperature heat pipe in hypersonic vehicles thermal protection, the principles and characteristics of high temperature heat pipe used in hypersonic vehicles thermal protection were introduced. The methods of numerical simulation, theory analysis and experiment research were utilized to analyze the frozen start-up and steady state characteristic of the heat pipe as well as the machining improvement for fabricating irregularly shaped heat pipe which is suitable for leading edge of hypersonic vehicles. The results indicate that the frozen start-up time of heat pipe is long （10 min） and there exists large temperature difference along the heat pipe （47 ℃/cm）, but the heat pipe can reduce the temperature in stagnation area of hypersonic vehicles from 1 926 to 982 ℃ and work normally during 1 000-1 200℃. How to improve the maximum heat transfer capability and reduce the time needed for start-up from frozen state of the heat pipe by optimizing thermostructure such as designing of a novel wick with high performance is the key point in hypersonic vehicles thermal protection of heat pipe.

THE INFLUENCE OF THE ZHOUSHAN ARCHIPELAGO TERRAIN ON THE TROPICAL CYCLONE TODD (NO.9806) BY NUMERICAL SIMULATION

Taking Tropical Cyclone (TC) No.9806 (Todd) as an example, the effects of Zhoushanarchipelago terrain on landfall TC are investigated by use of numerical simulation. Results show that, undertopographic influences of Zhoushan Islands, the westward-moving landfall TC deflects. And, small orographichighs and enhanced rainfall caused by climbing airflow on the windward slope of main mountains of theseislands are a result of effects of Zhoushan Islands. These results display some particular laws of effects ofsmall-sized islands on the landfall of TC.

CLIMATIC CHARACTERISTICS OF THE ONSET OF SOUTH CHINASEA SUMMERMONSOON II.INTER-DECADAL VARIATION

By using the 40-year NCEP (1958-1997) grid point reanalysis meteorological data, we analyzed the inter-decadal variation on the climatic characteristics of the onset of South China Sea summer monsoon. The results are as follows. (1) There was great difference on the onset date of the SCS summer monsoon between the first two decades and the last two decades. It was late on the 6th pentad of May for the first two decades and was on the 4th and 5th pentad of May for the next two decades. (2) Except for the third decade (1978-1987), the establishment of the monsoon rainfall was one to two pentads earlier than the onset of the summer monsoon in all other three decades. (3) The onset of the SCS monsoon is the result of the abrupt development and eastward advancement of the southwesterly monsoon over the Bay of Bengal. The four-decade analysis shows that there were abrupt development of the southwesterly monsoon over the Bay of Bengal between the 3rd and 4th pentad of May, but there was great difference between its eastward movement and its onset intensity. These may have important effect to the earlier or later onset of the SCS summer monsoon. (4) During the onset of the SCS summer monsoon, there were great difference in the upper and lower circulation feature between the first two and the next two decades. At the lower troposphere of the first two decades, the Indian-Burma trough was stronger and the center of the subtropical high was located more eastward. At the upper troposphere, the northward movement of the center of subtropical high was large and located more northward after it landed on the Indo-China Peninsula. After comparison, we can see that the circulation feature of the last two decades was favorable to the establishment and development of the SCS summer monsoon.

Analysis of Physical and Chemical Characteristics and Influence Factors of UCG

Based on the UCG(underground coal gasification) theory, the “three zones” which are oxidization zone, reduction zone, and drying zone, were divided; physical and chemical properties of each zone were analyzed. Factors, such as temperature, rate of water pouring, quantity of air blast, thickness of coal seam, and the operation pressure were discussed. Among the influencing factors, the temperature is the most important one.

Analysis on characteristics of damaged field of surrounding rock and face of FMTC

Numerical simulation-FLAC^(3D) and equivalent material simulation were carried out to analyze the damaged patterns and lows, distribution of plastic width of face, and surrounding rock of FMTC during the advance of work face with different thick coal seams based on engineering geology and exploitation technology of 1151(3) fully mechanized top-coal caving (FMTC) face in Xieqiao Colliery.The results show that there is damage,and the destruction characteristics of surrounding rock and coal mass are different obviously in asymmetric exploitation layout.The damaged zone in surrounding rock and the coal of the return airway is larger than that of the intake airway.Moreover, the retained coal pillars are all damaged by tension and shear fracture, and plastic zone in coal mass in the dip direction ahead of Face is nonuniform.There are large damage zones in roof and floor strata, surrounding rock, and coal of return and intake airways near work face.The damaged zone in the upper part of Face is larger than that in middle and lower parts.The fruits of this paper are of guiding significance for engineering practices, such as support design and choice, roadway supporting and maintaining, rock pressure control of FMTC face, etc.

Rules of distribution in a plastic zone of rocks surrounding a roadway affected by tectonic stress

In order to study the rules of distribution in a plastic zone of rocks, surrounding a roadway, affected by tectonic stress, we first analyzed the mechanics of a roadway affected by tectonic stress and derived a theoretical formula for the plastic zone of rocks surrounding a roadway. We also analyzed the distribution characteristics of the plastic zone under different levels of tectonic stress, vertical pressure, cohesion and friction angle of the surrounding rock. Secondly, we used numerical simulation to analyze the range and shape features of the plastic zone of rocks surrounding the roadway, given different tectonic stress levels. Finally we used a rock drilling detector to carry out field measurements on the broken state of rock surrounding the roadway at the –700 substation and channels in the Xinzhuang mine of the Shenhuo mining area. Given the measured ground stress, we analyzed the relationship between tectonic stress and the distribution of this plastic zone. Our results show that the range of the plastic zone at the top and bottom of the roadway increases with an increase in tectonic stress and this increase is especially obvious at the roadway corner.

The conditions for onset of beta oscillations in an extended subthalamic nucleus-globus pallidus network

The persistence of beta oscillations generated in the basal ganglia is typical pathological characteristic of Parkin-son's disease.In this paper,we construct a new network model,which is mainly composed of the subthalamic nu-cleus and the globus pallidus.Based on theory analysis,we obtained the beta shock conditions in the proposed model.It is shown that numerical simulations are in good agreement with analytical results.At last,the effects of cortex and striatum inputs on the generation of oscillations are also discussed.We hope that the obtained results may find applications to treatment of Parkinson’s disease and study of neural oscillations produced in other parts of the brain neural network.

Characteristics of nonlinear evolution of wavepackets in boundary layers

The nonlinear evolution of a finite-amplitude disturbance in a 3-D supersonic boundary layer over a cone was investigated recently by Liu et al. using direct numerical simulation （DNS）. It was found that certain small-scale 3-D disturbances amplified rapidly. These disturbances exhibit the characteristics of second modes, and the most amplified components have a well- defined spanwise wavelength, indicating a clear selectivity of the amplification. In the case of a cone, the three-dimensionality of the base flow and the disturbances themselves may be responsible for the rapid amplification. In order to ascertain which of these two effects are essential, in this study we carried out DNS of the nonlinear evolution of a spanwise localized disturbance （wavepacket） in a flat-plate boundary layer. A similar amplification of small-scale disturbances was observed, suggesting that the direct reason for the rapid amplification is the three-dimensionality of the disturbances rather than the three-dimensional nature of the base flow, even though the latter does alter the spanwise distribution of the disturbance. The rapid growth of 3-D waves may be attributed to the secondary instability mechanism. Further simulations were performed for a wavepacket of first modes in a supersonic boundary layer and of Tollmien-Schlichting （T-S） waves in an incompressible boundary layer. The re- suits show that the amplifying components are in the band centered at zero spanwise wavenumber rather than at a finite spanwise wavenumber. It is therefore concluded that the rapid growth of 3-D disturbances in a band centered at a preferred large spanwise wavenumber is the main characteristic of nonlinear evolution of second mode disturbances in supersonic boundary layers.