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.

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.

3-D fracture network dynamic simulation based on error analysis in rock mass of dam foundation

Accurate 3-D fracture network model for rock mass in dam foundation is of vital importance for stability,grouting and seepage analysis of dam foundation.With the aim of reducing deviation between fracture network model and measured data,a 3-D fracture network dynamic modeling method based on error analysis was proposed.Firstly,errors of four fracture volume density estimation methods(proposed by ODA,KULATILAKE,MAULDON,and SONG)and that of four fracture size estimation methods(proposed by EINSTEIN,SONG and TONON)were respectively compared,and the optimal methods were determined.Additionally,error index representing the deviation between fracture network model and measured data was established with integrated use of fractal dimension and relative absolute error(RAE).On this basis,the downhill simplex method was used to build the dynamic modeling method,which takes the minimum of error index as objective function and dynamically adjusts the fracture density and size parameters to correct the error index.Finally,the 3-D fracture network model could be obtained which meets the requirements.The proposed method was applied for 3-D fractures simulation in Miao Wei hydropower project in China for feasibility verification and the error index reduced from 2.618 to 0.337.

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.

3-D Simulation Study on Seismic Response of Bridge Piles in Landslide

The anti-slide support structure is widely used in the anti-seismic reinforcement of bridge foundations,but related experimental research was processing slowly. Based on the prototype of the Jiuzhaigou bridge at the Chengdu-Lanzhou Railway,a 3-D simulation model was established on the basis of the shaking table model test,and the rationality of the dynamic analysis model was verified by indicators such as the bending moment of the bridge piles,peak soil pressure,and PGA amplification factors. The results show that the inertia force of the bridge pier has an important influence on the deformation of the pile foundation. The bending moment and shearing force are larger in lateral bridge piles,and the maximum value is near the pile top. The PGA amplification factor is stronger in the back of the rear anti-slide piles and so is it in front of the bridge pier,and the soil is prone to slip and damage. The bedrock is rigid and the dynamic response is maintained at a low level. The anti-slide piles in the rear row play a major role in the anti-seismic reinforcement design,and the anti-slide piles in the front row can be used as an auxiliary support structure.

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.

Application of computer aided 3D simulation technique for complicated foot and ankle fractures

Objective To investigate the effect of computer aided 3D simulation technique for treating complicated foot and ankle fractures precisely.Methods From November 2007 to August 2009,255 patients with complicated foot and ankle fractures

Intrinsic background radiation of LaBr3(Ce)detector via coincidence measurements and simulations

The LaBr3(Ce)detector has attracted much attention in recent years because of its superior characteristics compared with other scintillating materials in terms of resolution and efficiency.However,it has a relatively high intrinsic background radiation because of the naturally occurring radioisotopes in lanthanum,actinium,and their daughter nuclei.This limits its applications in low-counting rate experiments.In this study,we identified the radioactive isotopes in theφ3"x 3"Saint-Gobain B380 detector by a coincidence measurement using a Clover detector in a low-background shielding system.Moreover,we carried out a Geant4 simulation of the experimental spectra to evaluate the activities of the main internal radiation components.The total activity of the background radiation of B380 is determined to be 1.523(34)Bq/cm^3.The main sources include 138La at 1.428(34)Bq/cm^3,207Tl at 0.0135(13)Bq/cm^3,211Bi at 0.0136(15)Bq/cm^3,215Po at 0.0135(3)Bq/cm^3,219Rn at 0.0125(12)Bq/cm^3,223Fr at 0.0019(11)Bq/cm^3,223Ra at 0.0127(10)Bq/cm^3,227Th at 0.0158(22)Bq/cm^3,and 227Ac at 0.0135(13)Bq/cm^3.Of these,the activities of 207Tl,211Po,215Po,223Fr,and 227Ac are deduced for the first time from the secular equilibrium established in the decay chain of 227Ac.

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.

Filling Simulation of Three-dimension Braided Composite in Resin Transfer Molding

This paper measured permeability of three-dimension braided preform by radial technology. The results show that principal permeability tensor coincided with their braiding axial direction. The software of one dimensional flow filling mold was designed using Visual C++ language. Filling time is predicted and validated. The result showed that the filling time of the mold centerline agrees with the prediction value. The filling time of the mould edge is shorter than that of the prediction. An actual plate of 3D braided preform/ modified polyarylacetylene composite is produced according to prediction value and validation analysis.

THE SYSTEM SIMULATION OF THREE-DIMENSIONAL RADAR

To provide a test platform for Electronic Warfare (EW) system, it is needed to simulate the radar received Intermediate Frequency (IF) signals and radar system functions.This letter gives a description of a radar system simulation software developed for frequencyphase scanning three-dimensional (3-D) radar. Experimental results prove that the software could be used for system evaluation and for training purposes as an attractive alternative to real EW system.

Permittivity of composites used for Luneburg lens antennas by drilling holes based on 3-D printing technique

Due to the attractive performances such as the ability of beam focus,broadband,multi-beam scanning and other features,Luneburg lens antennas are applied in multi-beam antenna,which overcomes the problem of gain loss produced by multi-beam parabolic antenna.Based on 3-D printing technique,Luneburg lens antennas by drilling holes are studied.Permittivity and loss tangent of the equivalent lens materials can be influenced by original materials,hole shapes,hole directions,and porosity.After tests,polystyrene with waxes may be the most appropriate materials for Luneburg lens with high strength.Permittivity with the shape of triangle is the lowest due to the homogeneity.Relative permittivities with the direction at a range of 15°-45°are lower while loss tangent at a range of 0°-30°.Radial directional holes are more appropriate for Luneburg lens.The relative permittivity is decreased with the increment of porosity.After calculations,the forecasts calculated by Looyenga and A-BG theory are more precise.Finally,Luneburg lens with two layers is fabricated by 3-D printing.

Magnetic field of mathematical modeling and simulation of 3D magnetic pole array spherical actuator

In this paper a new spherical actuator is designed and its advantages are compared to an existing spherical actuator, which function is limited by several design bottlenecks. First the output torque is too small. Second, the attitude is difficult to be accurately detected. The new three-dimen- sional magnetic pole array can solve these major problems. The new actuator features an outer rotor with multiple permanent magnet （PM） poles. Using an analytical solution and the finite element so- lution simulation, the feasibility of the approach is verified. A prototype was developed, tested, and experiments were conducted to obtain the practical value of the magnetic flux density.

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.

Three-dimensional Building Damage Simulation Based on GIS

GIS technology has been applied to building damage analysis around the world. However, most previous studies focused on the application of 2-D GIS technology, and the results from traditional earthquake damage prediction are displayed in 2-D figures and charts, which is incapable of demonstrating the 3-D spatial characteristics of buildings. Taking brick-concrete building as an example, we study the characteristics of building damage, and effectively combine the information of building textures and earthquake damage. Then, we apply Google SketchUp techniques to create building models and display them with seismic damage texture in the ArcGIS Engine software development environment. In this paper we propose a solid idea for 3-D simulation of earthquake damage, which is helpful in earthquake damage prediction, virtual emergency rescue practice and earthquake knowledge education.

Reconstruction of 3-D digital cores using a hybrid method

A 3-D digital core describes the pore space microstructure of rocks. An X-ray micro CT scan is the most accurate and direct but costly method to obtain a 3-D digital core. In this study, we propose a hybrid method which combines sedimentation simulation and simulated annealing （SA） method to generate 3-D digital cores based on 2-D images of rocks. The method starts with the sedimentation simulation to build a 3-D digital core, which is the initial configuration for the SA method. We update the initial digital core using the SA method to match the auto-correlation function of the 2-D rock image and eventually build the final 3-D digital core. Compared with the typical SA method, the hybrid method has significantly reduced the computation time. Local porosity theory is applied to quantitatively compare the reconstructed 3-D digital cores with the X-ray micro CT 3-D images. The results indicate that the 3-D digital cores reconstructed by the hybrid method have homogeneity and geometric connectivity similar to those of the X-ray micro CT image. The formation factors and permeabilities of the reconstructed 3-D digital cores are estimated using the finite element method （FEM） and lattice Boltzmann method （LBM）, respectively. The simulated results are in good agreement with the experimental measurements. Comparison of the simulation results suggests that the digital cores reconstructed by the hybrid method more closely reflect the true transport properties than the typical SA method alone.

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 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.

Visualization of 3-D temperature distribution in a 300 MW twin-furnace coal-fired boiler

Until now, it has been difficult to obtain on-line three-dimensional (3-D) temperature distribution information which can reflect the overall combustion condition in the furnace of a coal-fired power plant boiler. A combustion monitoring system is introduced which can solve the problem efficiently. Through this system, the 3-D temperature distribution in a coal-fired boiler furnace can be obtained using a novel flame image processing technique. Briefly, we first outline the visualization principle. Then, the hardware and software design of the system in a 300 MW twin-furnace coal-fired boiler are introduced in detail. The visualization of the 3-D temperature distribution in the twin-furnace boiler is realized with an industrial computer and the Distributed Control System (DCS) of the boiler. The practical operation of the system shows that it can provide valuable combustion information of a furnace and is useful for the combustion diagnosis and adjustment in coal-fired power plants.

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.