Numerical simulations of full-wave fi elds and analysis of channel wave characteristics in 3-D coal mine roadway models

Currently, numerical simulations of seismic channel waves for the advance detection of geological structures in coal mine roadways focus mainly on modeling two- dimensional wave fields and therefore cannot accurately simulate three-dimensional （3-D） full-wave fields or seismic records in a full-space observation system. In this study, we use the first-order velocity-stress staggered-grid finite difference algorithm to simulate 3-D full-wave fields with P-wave sources in front of coal mine roadways. We determine the three components of velocity Vx, Vy, and Vz for the same node in 3-D staggered-grid finite difference models by calculating the average value of Vy, and Vz of the nodes around the same node. We ascertain the wave patterns and their propagation characteristics in both symmetrical and asymmetric coal mine roadway models. Our simulation results indicate that the Rayleigh channel wave is stronger than the Love channel wave in front of the roadway face. The reflected Rayleigh waves from the roadway face are concentrated in the coal seam, release less energy to the roof and floor, and propagate for a longer distance. There are surface waves and refraction head waves around the roadway. In the seismic records, the Rayleigh wave energy is stronger than that of the Love channel wave along coal walls of the roadway, and the interference of the head waves and surface waves with the Rayleigh channel wave is weaker than with the Love channel wave. It is thus difficult to identify the Love channel wave in the seismic records. Increasing the depth of the receivers in the coal walls can effectively weaken the interference of surface waves with the Rayleigh channel wave, but cannot weaken the interference of surface waves with the Love channel wave. Our research results also suggest that the Love channel wave, which is often used to detect geological structures in coal mine stopes, is not suitable for detecting geological structures in front of coal mine roadways. Instead, the Rayleigh channel wave can be used for the advance detection of geological structures in coal mine roadways.

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.

Three-dimensional numerical simulation of geothermal fi eld in space-wavenumber domain

Large-scale,fine,and efficient numerical simulation of a geothermal field plays an important role in geothermal energy development.Confronted with the problem of large computation and high storage requirements for complex underground models in a three-dimensional(3-D)numerical simulation of a geothermal fi eld,a mixed space-wavenumber domain 3-D numerical simulation algorithm is proposed in this paper.According to the superposition principle of temperature field,the geothermal field is decomposed into background and abnormal temperature fi elds for calculation.The uniform layered model is used to solve the background field.When the abnormal field is solved,the horizontal two-dimensional(2-D)Fourier transform is used to transform the 3-D diff erential equation satisfi ed by an abnormal field into a series of one-dimensional ordinary differential equations with diff erent wavenumbers,which greatly reduces the calculation and storage.The unit division of an ordinary diff erential equation is fl exible,and the calculation amount is small.The algorithm fully takes advantage of the effi ciency of the Fourier transform and the quickness of the catch-up method to solve linear equations with a fixed bandwidth,which effectively improves the computational efficiency.Compared with the COMSOL Multiphysics professional simulation finite element software,the time consumption and memory requirements of the algorithm proposed in this paper are reduced by multiple orders of magnitude in terms of ensuring accuracy and the same mesh division.The more the number of calculated nodes is,the more obvious is the advantage.We design models to study the thermal conductivity,heat fl ux boundary,regional tectonic morphology,and topographic relief of the geothermal fi eld distribution.A 3-D geophysical model is developed based on topographic elevation data,geothermal geology,and geophysical exploration data in the Qiabuqia area of Gonghe Basin,Qinghai Province,China.Numerical simulation of the geothermal fi eld in this area is realized,which shows that the algorithm is suitable for precise and effi cient simulation of an arbitrary complex terrain and geological conditions.

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.

Numerical simulation and analysis of solid-liquid two-phase threedimensional unsteady flow in centrifugal slurry pump

Based on RNG k-ε turbulence model and sliding grid technique, solid-liquid two-phase three-dimensional(3-D) unsteady turbulence of full passage in slurry pump was simulated by means of Fluent software. The effects of unsteady flow characteristics on solid-liquid two-phase flow and pump performance were researched under design condition. The results show that clocking effect has a significant influence on the flow in pump, and the fluctuation of flow velocity and pressure is obvious, particularly near the volute tongue, at the position of small sections of volute and within diffuser. Clocking effect has a more influence on liquid-phase than on solid-phase, and the wake-jet structure of relative velocity of solid-phase is less obvious than liquid-phase near the volute tongue and the impeller passage outlet. The fluctuation of relative velocity of solid-phase flow is 7.6% smaller than liquid-phase flow at the impeller outlet on circular path. Head and radial forces of the impeller are 8.1% and 85.7% of fluctuation, respectively. The results provide a theoretical basis for further research for turbulence, improving efficient, reducing the hydraulic losses and wear. Finally, field tests were carried out to verify the operation and wear of slurry pump.

SIMULATION OF 3-D DEFORMATION AND MATERIAL FLOW DURING ROLL FORGING PROCESS USING SYSTEM OF OVAL-ROUND GROOVE

Basing on the analysis of the traits of the roll forging process, a system-model of computer simulation has been established. Three-dimensional rigid-plastic FEM has been used for the simulation of the deformation process in the oval and round pass rolling, including the entering, rolling, and separating stages. The analysis was conducted using the Deform-3D ver.5.0 code. The important information concerned with the deformation area characteristic, material flow, and velocity field has been presented. Otherwise, the location of the neutral plane in the deformation area was shown clearly.

Simulation of subcritical flow pattern in 180° uniform and convergent open-channel bends using SSIIM 3-D model

In meandering rivers, the flow pattern is highly complex, with specific characteristics at bends that are not observed along straight paths. A numerical model can be effectively used to predict such flow fields. Since river bends are not uniform-some are divergent and others convergent-in this study, after the SSIIM 3-D model was calibrated using the result of measurements along a uniform 180° bend with a width of 0.6 m, a similar but convergent 180v bend, 0.6 m to 0.45 m wide, was simulated using the SSI1M 3-D numerical model. Flow characteristics of the convergent 180° bend, including lengthwise and vertical velocity profiles, primary and secondary flows, lengthwise and widtbwise slopes of the water surface, and the helical flow strength, were compared with those of the uniform 180° bend. The verification results of the model show that the numerical model can effectively simulate the flow field in the uniform bend. In addition, this research indicates that, in a convergent channel, the maximum velocity path at a plane near the water surface crosses the channel＇s centerline at about a 30° to 40° cross-section, while in the uniform bend, this occurs at about the 50° cross-section. The varying range of the water surface elevation is wider in the convergent channel than in the uniform one, and the strength of the helical flow is generally greater in the uniform channel than in the convergent one. Also, unlike the uniform bend, the convergent bend exhibits no rotational cell against the main direction of secondary flow rotation at the 135° cross-section.

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 the 3D flow field in a ring die pellet mill considering viscous heating

Viscous heating has a substantial influence on the extrusion forming process and product quality of powder materials.This study selected the MUZL420 ring die pellet mill as the research object,from which a 3D flow physical model was established.The numerical simulation of 3D nonisothermal flow in the extrusion pelletizing process of granulated alfalfa was performed with POLYFLOW.The distribution laws of pressure,velocity,shear rate,viscosity,viscous heating and temperature in the flow field were revealed to thoroughly investigate the pelletizing process and provide a reference for structural optimization and process control.The results showed that two extrusion zones in the pelleting chamber were symmetrical with respect to the center,and the significant pressure gradient along the rotating direction of the ring die and the roller caused the material to flow back in the opposite direction.There were larger velocity gradients,shear rates and viscous heating levels in the deformation and compaction zone,the negative pressure zone behind the extrusion zone and the die holes.The distribution of viscosity was opposite to that of the shear rate.The temperature increase area caused by viscous heating gradually expanded from the material inlet to the bottom of the extrusion chamber along the Z-axis direction,and the temperature increased accordingly.The extrusion force and the forming temperature in the extrusion forming zone were captured in the numerical simulation.The extrusion forming density was calculated with the regression prediction model established through the simulation experiment of pelletizing with a ring die.Through a comparison with the results of mean alfalfa pellet density from the ring die pellet mill experiment,the relative error was less than 5%,which indicated that the numerical simulation method was reliable.

Numerical simulation on the evolution of cloud particles in 3-D convective cloud

A 3-D convective cloud model with compressible non-hydrostatic dynamics and the spectral bin microphysics of a 2-D slab-symmetric model has been used to simulate an observed supercell storm occurring on 29 June, 2000 near Bird City, Kansas, USA. The main

Simulation of Polymer Melt Flow Fields in Intermeshing Co-Rotating Three-Screw Extruders

Three-dimension isothermal flows of polymer melt in the kneading blocks of triangularly-arranged and parallelly-arranged intermeshing co-rotating three-screw extruders are simulated using the finite element package POLYFLOW. Based on the velocity fields calculated, the particle trajectories in both machines are visualized using particle tracking technique. The numerical results indicate that the flow patterns in three-screw extruders are similar to those in twin-screw extruders. The triangularly-arranged three-screw extruder has the largest pumping capacity and also the highest extrusion stability in terms of flowrate fluctuation with screw rotation. The instantaneous mixing and cumulative residence time distribution （RTD） characteristics are also analyzed and compared with traditional intermeshing co-rotating twin-screw extruders. It is shown that the start section of the cumulative RTD curve for the triangularly-arranged machine has a small shoulder, which is attributed to the faster flow in the central region of this type of extruder.

The Application of 3-D Visible Technology to Reservoir Management

The paper deals with the application of 3 -D visible technology to reservoir management. Making use of this method for expanding - spread data point in reservoir management, can discard the false and retain the true during data recording. As a result, The quality of data recording is ensured. In reservoir description, the reservoir characteristics, such as space distribution,physical change and fluid distribution may be identified by restoring palaeostructures, building - up 3-D facics tract model and 3-D fracture system model. Seismic interpretation, geologic modeling and numerical simulation are well integrated so that they can be promote reservoir performance management to develop into the intensive management pattern.

A SIMPLIFIED NUMERICAL MODEL OF 3-D GROUNDWATER AND SOLUTE TRANSPORT AT LARGE SCALE AREA

A simplified numerical model of groundwater and solute transport is developed. At large scale area there exists a big spatial scale difference between horizontal and vertical length scales. In the resultant model, the seepage region is particularly divided into several virtual layers along the z direction and vertical 1-D columns covering x-y 2-D area according to stratum properties. The numerical algorithm is replacing the full 3-D water and mass balance analysis as the 2-D Galerkin finite element method in x- and y-directions and 1-D finite differential approach in the z direction. The reasonable method of giving minimum thickness is successfully used to handle transient change of water table, drying cells and problem of rewetting. The solution of the simplified model is preconditioned conjugate gradient and ORTHOMIN method. The validity of the developed 3-D groundwater model is tested with the typical pumping and backwater scenarios. Results of water balance of the computed example reveal the model computation reliability. Based on a representative 3-D pollution case, the solute transport module is tested against computing results using the MT3DMS. The capability and high efficiency to predict non-stationary situations of free groundwater surface and solute plume in regional scale problem is quantitatively investigated. It is shown that the proposed model is computationally effective.

Numerical Simulation of a Torrential Rain Event in the Northeast of Huaihe Basin.PartⅠ:Model Verification and Analysis of MCSs

A torrential rain event accompanying Typhoon Prapiroon occurred in 2000, with 24-h rainfall amount reaching 800 mm near the typhoon center. This event is simulated by the nonhydrostatic mesoscaie model ARPS （V5.2）, with thriple one-way nested-grids. Grid spacings of 45, 15, and 5 km are chosen for the three nested domains. The corresponding grid sizes are 75×75, 140× 140, and 180× 180, respectively. The NCEP/NCAR reanalysis data, radar echoes, and CMS-5 satellite images are assimilated with the ARPS model initially using a 3-D data assimilation system-ADAS. The simple ice phase scheme and the Kain- Fritsch cumulus parameterization scheme axe used. There are 35 layers in the vertical, with a vertical grid spacing of about 625 m. The integration is performed up to 48 h from 0800 BT 29 to 0800 BT 31 August 2000. Compared with radar echoes, GMS-5 satellite images, and intensive surface observations, the results show that the heavy rain area down between the 500-hPa trough and the subtropical high in the left-front of Prapiroon is well simulated by the model ARPS, and the simulated rainfall centers are consistent with observations. A comparison of the radar echoes with these retrieved from the simulated hydrometeors reveals that there are meso-β scale convective systems that exhibit distinctive characteristics, and there are four convective belts converging in the vicinity of Xiangshui, where the maximum rainfall is observed. A further comparision of skew T-lgp diagrams from simulated and observed data demonstrates significant instability in this torrential rain process. The persistent vertical wind shear provides kinetic energy for the development of the MCSs, hence promoting the baroclinic development of convective cells, and the concentration of heavy rain at the specific location. The consistency between model results and observations encourages a further study of the torrential rain event using the simulation data.

Ensemble Cloud Model Application in Simulating the Catastrophic Heavy Rainfall Event

An attempt has been made in the present research to simulate a deadly flash-flood event over the City of Skopje,Macedonia on 6 August 2016.A cloud model ensemble forecast method is developed to simulate a super-cell storm’s initiation and evolutionary features.Sounding data are generated using an ensemble approach,that utilizes a triple-nested WRF model.A three-dimensional(3-D)convective cloud model(CCM)with a very fine horizontal grid resolution of 250-m is initialized,using the initial representative sounding data,derived from the WRF 1-km forecast outputs.CCM is configured and run with an open lateral boundary conditions LBC,allowing explicit simulation of convective scale processes.This preliminary study showed that the ensemble approach has some advantages in the generation of the initial data and the model initialization.The applied method minimizes the uncertainties and provides a more qualitative-quantitative assessment of super-cell storm initiation,cell structure,evolutionary properties,and intensity.A high-resolution 3-D run is capable to resolve detailed aspects of convection,including high-intensity convective precipitation.The results are significant not only from the aspect of the cloud model’s ability to provide a qualitative-quantitative assessment of intense precipitation but also for a deeper understanding of the essence of storm development,its vortex dynamics,and the meaning of micro-physical processes for the production and release of large amounts of precipitation that were the cause of the catastrophic flood in an urban area.After a series of experiments and verification,such a system could be a reliable tool in weather services for very short-range forecasting(now-casting)and early warning of weather disasters.

Numerical study on structure of meso-β scale rainstorm system during IOP523 of HUAMEX

One of the scientific goals of HUAMEX (Torrential Rainfall Experiment over Both Sides of the Taiwan Strait and Adjacent Area), a special project in the National '863' Scientific Project, is to reveal meso-β scale structure of heavy rainfall. However, limited by the observational instruments, the resolution of the data attained from the outfield scientific observation network cannot reveal the meso-β scale structure. A 36-h numerical simulation of heavy rainfall case from May 23 to May 24, 1998 during IOP523 of HUAMEX was performed by the mesoscale numerical model of PSU/NCAR MM5, and the meso-p scale structure of con-vective system was examined based on the model output in this note. The results point out that the meso-β scale features of the heavy rainfall are clear considerably at low, middle and upper levels, respectively. There is a strong ageostrophic feature between the geopotential height field and the wind field on the upper level. The vertical cross sections of 54 km coarse grid and 18 km

离心压缩机级内三维粘性流动数值分析

本文介绍了用时间推进法进行离心压缩机级内全三维粘性流场的计算方法．该方法是作者在已完成的求解孤立叶排Ｎ－Ｓ方程的计算程序的基础上发展的，通过使用周向非均匀的叶排间隙混合平面概念来考虑离心叶轮与扩压器之间的相互影响，并依据质量守恒原则，对下游的进口气流角进行了修正，以保证了上下游叶排在相容的流动特性下工作．完成了离心压缩机级（包括叶轮与扩压器）三维粘性定常流场的计算及性能预测．应用本文所发展的程序，对某高炉鼓风机的叶轮及其叶片扩压器进行了计算．

离心风机现代设计方法研究

提出一种新的离心风机现代设计方法，它包含3部分内容：①以现有的工程设计为基础;②利用正在发展的技术来进行离心风机气动力优化设计局用现场性能试验进行考核。其中，关键和难点是三维粘性流场数值模拟。按此方法研制出多种样机，气动和噪声性能均有明显提高，证实这种方法是正确的。

轴流压气机级内三维粘性流动的数值模拟

考虑轴流压气机中实际气流沿径向和周向都不均匀的特点,采用亚松弛方法改进传统的混合界面模型,发展了轴流式压气机级内全三维粘性流场的数值计算程序,并以某轴流式双排对转压气机的进口导叶/转子为研究对象,将本文发展的程序与商用NUMECA软件的结果进行对比,对本文发展的程序进行了验证。对比结果表明:由于混合界面法引入了"人工掺混",商业软件计算得出的总压在通过交界面时存在一定的总压损失,而自编程序计算结果保证了总压在交界面处的一致性。流场分析结果表明:该转子叶片在设计工况下处于负攻角状态,在导叶出口轮毂处、压力面50%弦长处、吸力面靠近尾缘处均出现了轻微的流动分离现象。

三级风扇的三维流场计算

本文运用三维粘性流体计算软件Fine/NUMECA对某三级风扇多个转速下的三维流场进行了计算,并对其第一级三个典型转速进行了单独计算。分析了此风扇在不同转速下的特性变化,重点关注了设计点和近喘振边界工作点的流场特性,并分析了在设计转速时发生喘振的原因。结果表明,此风扇设计为典型的跨声速压气机设计,具有较高的效率,但仍有改进的余地。