Mountainous regions have disadvantages in economic development because of harsh physical and climatic conditions.However,winter tourism activities are one of the key components for supporting economic development in t...Mountainous regions have disadvantages in economic development because of harsh physical and climatic conditions.However,winter tourism activities are one of the key components for supporting economic development in the highlands.Establishing a ski resort area supports direct and indirect employment in a region,and it stops immigration from mountainous regions to other places.This research aimed to assess the potential ski areas using a multi criteria evaluation technique in the Van region which is located in the eastern part of Türkiye.In this context,snow cover duration,sun effect,slope,slope length,elevation,population density,distance from main roads and lake visibility were used as input factors in the decision making process.Each factor was standardized using a fuzzy technique based on existing well-known ski centers in Türkiye.The weight of inputs was defined by applying a survey to the professional skiers.The most important factors were detected as transportation opportunities and snow covers whereas,the least important factors were elevation and population density.Additionally,lake visibility was very important to make a difference from other existing facilities in the region.Therefore,it was included as constraints and lake visible areas were extracted at the final stage of the research.Potential ski areas were mapped in three levels as professional,intermediate and beginner skiers.One of the suitable areas was selected as a sample projection and for the 3D simulation of the ski investment area.Potential costs and benefits were discussed.It was found that a ski tourism area investment can be amortized in 3 years in the region.展开更多
Airblasts,as one common phenomenon accompanied by rapid movements of landslides or rock/snow avalanches,commonly result in catastrophic damages and are attracting more and more scientific attention.To quantitatively a...Airblasts,as one common phenomenon accompanied by rapid movements of landslides or rock/snow avalanches,commonly result in catastrophic damages and are attracting more and more scientific attention.To quantitatively analyze the intensity of airblast initiated by landslides,the Wangjiayan landslide,occurred in the Wenchuan earthquake,is selected here with the landslide propagation and airblast evolution being studied using FLUENT by introducing the Voellmy rheological law.The results reveal that:(1)For the Wangjiayan landslide,its whole travelling duration is only 12 s with its maximum velocity reaching 36 m/s at t=10 s;(2)corresponding to the landslide propagation,the maximum velocity,28 m/s,of the airblast initiated by the landslide also appears at t=10 s with its maximum pressure reaching594.8 Pa,which is equivalent to violent storm;(3)under the attack of airblast,the load suffered by buildings in the airblast zone increases to 1300 Pa at t=9.4 s and sharply decreased to-7000 Pa as the rapid decrease of the velocity of the sliding mass at t=10 s,which is seriously unfavorable for buildings and might be the key reason for the destructive collapse of buildings in the airblast zone of the Wangjiayan landslide.展开更多
Single event transient of a real p-n junction in a 0.18μm bulk process is studied by 3D TCAD simulation. The impact of voltage, temperature, substrate concentration, and LET on SET is studied. Our simulation results ...Single event transient of a real p-n junction in a 0.18μm bulk process is studied by 3D TCAD simulation. The impact of voltage, temperature, substrate concentration, and LET on SET is studied. Our simulation results demonstrate that biases in the range 1.62 to 1.98V influence DSET current shape greatly and total collected charge weakly. Peak current and charge collection within 2ns decreases as temperature increases,and temperature has a stronger influence on SET currents than on total charge. Typical variation of substrate concentration in modern VDSM processes has a negligible effect on SEEs. Both peak current and total collection charge increases as LET increases.展开更多
The standard three dimensional(3D) k-ε turbulence model was applied to simulate the flow field of a small scale combined oxidation ditch. The moving mesh approach was used to model the rotor of the ditch. Compariso...The standard three dimensional(3D) k-ε turbulence model was applied to simulate the flow field of a small scale combined oxidation ditch. The moving mesh approach was used to model the rotor of the ditch. Comparison of the computed and the measured data is acceptable. A vertical reverse flow zone in the ditch was found, and it played a very important role in the ditch flow behavior. The flow pattern in the ditch is discussed in detail, and approaches are suggested to improve the hydrodynamic performance in the ditch.展开更多
This paper presents a 3D simulator used for studying the motion control and autonomous navigation of robotic fish. The simulator’s system structure and computation flow are presented. Simplified kinematics and hydrod...This paper presents a 3D simulator used for studying the motion control and autonomous navigation of robotic fish. The simulator’s system structure and computation flow are presented. Simplified kinematics and hydrodynamics models for a virtual robotic fish are proposed. Many other object models are created for water, obstacles, sonar sensors and a swimming pool. Experimental results show that the simulator provides a realistic and convenient way to develop autonomous navigation algorithms for robotic fish.展开更多
In this paper, near-fault strong ground motions caused by a surface rupture fault (SRF) and a buried fault (BF) are numerically simulated and compared by using a time-space-decoupled, explicit finite element metho...In this paper, near-fault strong ground motions caused by a surface rupture fault (SRF) and a buried fault (BF) are numerically simulated and compared by using a time-space-decoupled, explicit finite element method combined with a multi-transmitting formula (MTF) of an artificial boundary. Prior to the comparison, verification of the explicit element method and the MTF is conducted. The comparison results show that the final dislocation of the SRF is larger than the BF for the same stress drop on the fault plane. The maximum final dislocation occurs on the fault upper line for the SRF; however, for the BE the maximum final dislocation is located on the fault central part. Meanwhile, the PGA, PGV and PGD of long period ground motions (≤ 1 Hz) generated by the SRF are much higher than those of the BF in the near-fault region. The peak value of the velocity pulse generated by the SRF is also higher than the BE Furthermore, it is found that in a very narrow region along the fault trace, ground motions caused by the SRF are much higher than by the BF. These results may explain why SRFs almost always cause heavy damage in near-fault regions compared to buried faults.展开更多
Weft knitted fancy fabrics are widely used in knitted garment design. Due to the complexity of the structures, their modeling and simulation needs to be solved in three-dimensional (3D) CAD developments. In this paper...Weft knitted fancy fabrics are widely used in knitted garment design. Due to the complexity of the structures, their modeling and simulation needs to be solved in three-dimensional (3D) CAD developments. In this paper, 3D loop geometrical models of weft knitted fancy structures, including tuck stitch, jacquard stitch, transfer stitch and fleecy stitch, were developed based on an improved model of plain loop, and their central axes as some 3D space curves were achieved by using Non-Uniform Rational B-Splines (NURBS). The 3D visual simulation programme was written in C++ programming language using OpenGL, which was a function library of 3D graphics. Some examples of weft knitted fancy fabrics were generated and practical application of 3D simulation was discussed.展开更多
When heavy machines and large scaled receiver system of communication equipment are manufactured, it always needs to produce large-sized steel castings, aluminum castings and etc. Some defects of hot cracking by therm...When heavy machines and large scaled receiver system of communication equipment are manufactured, it always needs to produce large-sized steel castings, aluminum castings and etc. Some defects of hot cracking by thermal stress often appear during solidification process as these castings are produced, which results in failure of castings. Therefore predicting the effects of technological parameters for production of castings on the thermal stress during solidification process becomes an important means. In this paper, the mathematical models have been established and numerical calculation of temperature fields by using finite difference method (FDM) and then thermal stress fields by using finite element method (FEM) during solidification process of castings have been carried out. The technological parameters of production have been optimized by the results of calculation and the defects of hot cracking have been eliminated. Modeling and simulation of 3D thermal stress during solidification processes of large-sized castings provided a scientific basis, which promoted further development of advanced manufacturing technique.展开更多
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 accu...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.展开更多
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 p...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.展开更多
Minimally invasive surgery is an important technique used for cytopathological examination.Recently,multiple studies have been conducted on a three-dimensional(3D)puncture simulation model as it can reveal the interna...Minimally invasive surgery is an important technique used for cytopathological examination.Recently,multiple studies have been conducted on a three-dimensional(3D)puncture simulation model as it can reveal the internal deformation state of the tissue at the micro level.In this study,a viscoelastic constitutive equation suitable for muscle tissue was derived.Additionally,a method was developed to define the fracture characteristics of muscle tissue material during the simulation process.The fracture of the muscle tissue in contact with the puncture needle was simulated using the cohesive zone model and a 3D puncture finite element model was established to analyze the deformation of the muscle tissue.The stress nephogram and reaction force under different parameters were compared and analyzed to study the deformation of the biological soft tissue and guide the actual operation process and reduce pain.展开更多
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 ...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.展开更多
This paper presents a new approach to parallelize 3D lattice Monte Carlo algorithms used in the numerical simulation of polymer on ZiQiang 2000 a cluster of symmetric multiprocessors (SMPs). The combined load for cel...This paper presents a new approach to parallelize 3D lattice Monte Carlo algorithms used in the numerical simulation of polymer on ZiQiang 2000 a cluster of symmetric multiprocessors (SMPs). The combined load for cell and energy calculations over the time step is balanced together to form a single spatial decomposition. Basic aspects and strategies of running Monte Carlo calculations on parallel computers are studied. Different steps involved in porting the software on a parallel architecture based on ZiQiang 2000 running under Linux and MPI are described briefly. It is found that parallelization becomes more advantageous when either the lattice is very large or the model contains many cells and chains.展开更多
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 o...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.展开更多
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 an...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展开更多
Objective To study the effect of using improved 2D computer-assisted fluoroscopic navigation through simulating 3D vertebrae image to guide pedicle screw internal fixation.Methods Posterior pedicle screw internal fixa...Objective To study the effect of using improved 2D computer-assisted fluoroscopic navigation through simulating 3D vertebrae image to guide pedicle screw internal fixation.Methods Posterior pedicle screw internal fixation,distraction展开更多
Reliable prediction of the shale fracturing process is a challenging problem in exploiting deep shale oil and gas resources.Complex fracture networks need to be artificially created to employ deep shale oil and gas re...Reliable prediction of the shale fracturing process is a challenging problem in exploiting deep shale oil and gas resources.Complex fracture networks need to be artificially created to employ deep shale oil and gas reserves.Randomly distributed minerals and heterogeneities in shales significantly affect mechanical properties and fracturing behaviors in oil and gas exploitation.Describing the actual microstructure and associated heterogeneities in shales constitutes a significant challenge.The RFPA3D(rock failure process analysis parallel computing program)-based modeling approach is a promising numerical technique due to its unique capability to simulate the fracturing behavior of rocks.To improve traditional numerical technology and study crack propagation in shale on the microscopic scale,a combination of high-precision internal structure detection technology with the RFPA^(3D) numerical simulation method was developed to construct a real mineral structure-based modeling method.First,an improved digital image processing technique was developed to incorporate actual shale microstructures(focused ion beam scanning electron microscopy was used to capture shale microstructure images that reflect the distri-butions of different minerals)into the numerical model.Second,the effect of mineral inhomogeneity was considered by integrating the mineral statistical model obtained from the mineral nanoindentation experiments into the numerical model.By simulating a shale numerical model in which pyrite particles are wrapped by organic matter,the effects of shale microstructure and applied stress state on microcrack behavior and mechanical properties were investigated and analyzed.In this study,the effect of pyrite particles on fracture propagation was systematically analyzed and summarized for the first time.The results indicate that the distribution of minerals and initial defects dominated the fracture evolution and the failure mode.Cracks are generally initiated and propagated along the boundaries of hard mineral particles such as pyrite or in soft minerals such as organic matter.Locations with collections of hard minerals are more likely to produce complex fractures.This study provides a valuable method for un-derstanding the microfracture behavior of shales.展开更多
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 ...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.展开更多
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 frequenci...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.展开更多
文摘Mountainous regions have disadvantages in economic development because of harsh physical and climatic conditions.However,winter tourism activities are one of the key components for supporting economic development in the highlands.Establishing a ski resort area supports direct and indirect employment in a region,and it stops immigration from mountainous regions to other places.This research aimed to assess the potential ski areas using a multi criteria evaluation technique in the Van region which is located in the eastern part of Türkiye.In this context,snow cover duration,sun effect,slope,slope length,elevation,population density,distance from main roads and lake visibility were used as input factors in the decision making process.Each factor was standardized using a fuzzy technique based on existing well-known ski centers in Türkiye.The weight of inputs was defined by applying a survey to the professional skiers.The most important factors were detected as transportation opportunities and snow covers whereas,the least important factors were elevation and population density.Additionally,lake visibility was very important to make a difference from other existing facilities in the region.Therefore,it was included as constraints and lake visible areas were extracted at the final stage of the research.Potential ski areas were mapped in three levels as professional,intermediate and beginner skiers.One of the suitable areas was selected as a sample projection and for the 3D simulation of the ski investment area.Potential costs and benefits were discussed.It was found that a ski tourism area investment can be amortized in 3 years in the region.
基金supported by the National Natural Science Foundation of China(42322702,42177131)。
文摘Airblasts,as one common phenomenon accompanied by rapid movements of landslides or rock/snow avalanches,commonly result in catastrophic damages and are attracting more and more scientific attention.To quantitatively analyze the intensity of airblast initiated by landslides,the Wangjiayan landslide,occurred in the Wenchuan earthquake,is selected here with the landslide propagation and airblast evolution being studied using FLUENT by introducing the Voellmy rheological law.The results reveal that:(1)For the Wangjiayan landslide,its whole travelling duration is only 12 s with its maximum velocity reaching 36 m/s at t=10 s;(2)corresponding to the landslide propagation,the maximum velocity,28 m/s,of the airblast initiated by the landslide also appears at t=10 s with its maximum pressure reaching594.8 Pa,which is equivalent to violent storm;(3)under the attack of airblast,the load suffered by buildings in the airblast zone increases to 1300 Pa at t=9.4 s and sharply decreased to-7000 Pa as the rapid decrease of the velocity of the sliding mass at t=10 s,which is seriously unfavorable for buildings and might be the key reason for the destructive collapse of buildings in the airblast zone of the Wangjiayan landslide.
文摘Single event transient of a real p-n junction in a 0.18μm bulk process is studied by 3D TCAD simulation. The impact of voltage, temperature, substrate concentration, and LET on SET is studied. Our simulation results demonstrate that biases in the range 1.62 to 1.98V influence DSET current shape greatly and total collected charge weakly. Peak current and charge collection within 2ns decreases as temperature increases,and temperature has a stronger influence on SET currents than on total charge. Typical variation of substrate concentration in modern VDSM processes has a negligible effect on SEEs. Both peak current and total collection charge increases as LET increases.
基金The Specialized Research Fund for the Doctoral Programof Higher Education(No.20010610023) and the Sino-Finnish Scientific and TechnologicalCooperation Program
文摘The standard three dimensional(3D) k-ε turbulence model was applied to simulate the flow field of a small scale combined oxidation ditch. The moving mesh approach was used to model the rotor of the ditch. Comparison of the computed and the measured data is acceptable. A vertical reverse flow zone in the ditch was found, and it played a very important role in the ditch flow behavior. The flow pattern in the ditch is discussed in detail, and approaches are suggested to improve the hydrodynamic performance in the ditch.
文摘This paper presents a 3D simulator used for studying the motion control and autonomous navigation of robotic fish. The simulator’s system structure and computation flow are presented. Simplified kinematics and hydrodynamics models for a virtual robotic fish are proposed. Many other object models are created for water, obstacles, sonar sensors and a swimming pool. Experimental results show that the simulator provides a realistic and convenient way to develop autonomous navigation algorithms for robotic fish.
基金National Natural Science Foundation of China Under Grant No. 50408003National Scientifi c and Technical Supporting Programs Funded by Ministry of Science & Technology of China Under Grant No. 2006BAC13B01
文摘In this paper, near-fault strong ground motions caused by a surface rupture fault (SRF) and a buried fault (BF) are numerically simulated and compared by using a time-space-decoupled, explicit finite element method combined with a multi-transmitting formula (MTF) of an artificial boundary. Prior to the comparison, verification of the explicit element method and the MTF is conducted. The comparison results show that the final dislocation of the SRF is larger than the BF for the same stress drop on the fault plane. The maximum final dislocation occurs on the fault upper line for the SRF; however, for the BE the maximum final dislocation is located on the fault central part. Meanwhile, the PGA, PGV and PGD of long period ground motions (≤ 1 Hz) generated by the SRF are much higher than those of the BF in the near-fault region. The peak value of the velocity pulse generated by the SRF is also higher than the BE Furthermore, it is found that in a very narrow region along the fault trace, ground motions caused by the SRF are much higher than by the BF. These results may explain why SRFs almost always cause heavy damage in near-fault regions compared to buried faults.
基金Natural Science Foundation of Tianjin,China( No. 11JCYBJC26400) Tianjin High School Scientific and Technology Fund Planning Project,China( No. 20100310)
文摘Weft knitted fancy fabrics are widely used in knitted garment design. Due to the complexity of the structures, their modeling and simulation needs to be solved in three-dimensional (3D) CAD developments. In this paper, 3D loop geometrical models of weft knitted fancy structures, including tuck stitch, jacquard stitch, transfer stitch and fleecy stitch, were developed based on an improved model of plain loop, and their central axes as some 3D space curves were achieved by using Non-Uniform Rational B-Splines (NURBS). The 3D visual simulation programme was written in C++ programming language using OpenGL, which was a function library of 3D graphics. Some examples of weft knitted fancy fabrics were generated and practical application of 3D simulation was discussed.
文摘When heavy machines and large scaled receiver system of communication equipment are manufactured, it always needs to produce large-sized steel castings, aluminum castings and etc. Some defects of hot cracking by thermal stress often appear during solidification process as these castings are produced, which results in failure of castings. Therefore predicting the effects of technological parameters for production of castings on the thermal stress during solidification process becomes an important means. In this paper, the mathematical models have been established and numerical calculation of temperature fields by using finite difference method (FDM) and then thermal stress fields by using finite element method (FEM) during solidification process of castings have been carried out. The technological parameters of production have been optimized by the results of calculation and the defects of hot cracking have been eliminated. Modeling and simulation of 3D thermal stress during solidification processes of large-sized castings provided a scientific basis, which promoted further development of advanced manufacturing technique.
文摘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.
基金financially supported by the State Oceanic Administration of China(Grant No.GHME2011CL01)the Program of State Key Laboratory of Coastal and Offshore Engineering(Grant No.LP1102)
文摘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.
基金Natural Science Foundation of Shandong Province(Grant No.ZR2019JQ19)Interdisciplinary Research Project of Shandong University(Grant No.2017JC027)China Scholarship Council(CSC).
文摘Minimally invasive surgery is an important technique used for cytopathological examination.Recently,multiple studies have been conducted on a three-dimensional(3D)puncture simulation model as it can reveal the internal deformation state of the tissue at the micro level.In this study,a viscoelastic constitutive equation suitable for muscle tissue was derived.Additionally,a method was developed to define the fracture characteristics of muscle tissue material during the simulation process.The fracture of the muscle tissue in contact with the puncture needle was simulated using the cohesive zone model and a 3D puncture finite element model was established to analyze the deformation of the muscle tissue.The stress nephogram and reaction force under different parameters were compared and analyzed to study the deformation of the biological soft tissue and guide the actual operation process and reduce pain.
基金Project supported partly by the National Natural Science Foundation of China (Grant Nos 10301034 and 40574069), The authors thank Professor Du Q very much for his important discussions.
文摘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.
文摘This paper presents a new approach to parallelize 3D lattice Monte Carlo algorithms used in the numerical simulation of polymer on ZiQiang 2000 a cluster of symmetric multiprocessors (SMPs). The combined load for cell and energy calculations over the time step is balanced together to form a single spatial decomposition. Basic aspects and strategies of running Monte Carlo calculations on parallel computers are studied. Different steps involved in porting the software on a parallel architecture based on ZiQiang 2000 running under Linux and MPI are described briefly. It is found that parallelization becomes more advantageous when either the lattice is very large or the model contains many cells and chains.
文摘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.
文摘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
文摘Objective To study the effect of using improved 2D computer-assisted fluoroscopic navigation through simulating 3D vertebrae image to guide pedicle screw internal fixation.Methods Posterior pedicle screw internal fixation,distraction
基金supported by the Central Program of Basic Science of the National Natural Science Foundation of China(No.72088101)"The theory and application of resource and environment management in the digital economy era"+1 种基金The National Natural Science Foundation of China(No.41941018)Scientific research and technological development program of RIPED,"major research of basic geologic and synergy research of engineering practice on Gulong shale oil"(No.2021ycq01).
文摘Reliable prediction of the shale fracturing process is a challenging problem in exploiting deep shale oil and gas resources.Complex fracture networks need to be artificially created to employ deep shale oil and gas reserves.Randomly distributed minerals and heterogeneities in shales significantly affect mechanical properties and fracturing behaviors in oil and gas exploitation.Describing the actual microstructure and associated heterogeneities in shales constitutes a significant challenge.The RFPA3D(rock failure process analysis parallel computing program)-based modeling approach is a promising numerical technique due to its unique capability to simulate the fracturing behavior of rocks.To improve traditional numerical technology and study crack propagation in shale on the microscopic scale,a combination of high-precision internal structure detection technology with the RFPA^(3D) numerical simulation method was developed to construct a real mineral structure-based modeling method.First,an improved digital image processing technique was developed to incorporate actual shale microstructures(focused ion beam scanning electron microscopy was used to capture shale microstructure images that reflect the distri-butions of different minerals)into the numerical model.Second,the effect of mineral inhomogeneity was considered by integrating the mineral statistical model obtained from the mineral nanoindentation experiments into the numerical model.By simulating a shale numerical model in which pyrite particles are wrapped by organic matter,the effects of shale microstructure and applied stress state on microcrack behavior and mechanical properties were investigated and analyzed.In this study,the effect of pyrite particles on fracture propagation was systematically analyzed and summarized for the first time.The results indicate that the distribution of minerals and initial defects dominated the fracture evolution and the failure mode.Cracks are generally initiated and propagated along the boundaries of hard mineral particles such as pyrite or in soft minerals such as organic matter.Locations with collections of hard minerals are more likely to produce complex fractures.This study provides a valuable method for un-derstanding the microfracture behavior of shales.
基金the National Natural Science Foundation of China(No.42272204)the Fundamental Research Funds for the Central Universities(Grant No.2021JCCXDC02)+3 种基金the Gansu Province Science and Technology Major Project(19ZD2GA005)for their supportfinancially supported by the State Key Laboratory for Geomechanics and Deep Underground Engineering(SKLGDUEK2020)Huaneng Group headquarters science and technology project(HNKJ21-H07)the Coal Burst Research Center of Jiangsu,China。
文摘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.
基金Research Grants for Returned Students of China under Grant No.2020-038the National Natural Science Foundation of China under Grant No.51408392。
文摘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.