Three-dimensional finite element simulation and reconstruction of jointed rock models using CT scanning and photogrammetry

The geometry of joints has a significant influence on the mechanical properties of rocks.To simplify the curved joint shapes in rocks,the joint shape is usually treated as straight lines or planes in most laboratory experiments and numerical simulations.In this study,the computerized tomography (CT) scanning and photogrammetry were employed to obtain the internal and surface joint structures of a limestone sample,respectively.To describe the joint geometry,the edge detection algorithms and a three-dimensional (3D) matrix mapping method were applied to reconstruct CT-based and photogrammetry-based jointed rock models.For comparison tests,the numerical uniaxial compression tests were conducted on an intact rock sample and a sample with a joint simplified to a plane using the parallel computing method.The results indicate that the mechanical characteristics and failure process of jointed rocks are significantly affected by the geometry of joints.The presence of joints reduces the uniaxial compressive strength (UCS),elastic modulus,and released acoustic emission (AE) energy of rocks by 37%–67%,21%–24%,and 52%–90%,respectively.Compared to the simplified joint sample,the proposed photogrammetry-based numerical model makes the most of the limited geometry information of joints.The UCS,accumulative released AE energy,and elastic modulus of the photogrammetry-based sample were found to be very close to those of the CT-based sample.The UCS value of the simplified joint sample (i.e.38.5 MPa) is much lower than that of the CT-based sample (i.e.72.3 MPa).Additionally,the accumulative released AE energy observed in the simplified joint sample is 3.899 times lower than that observed in the CT-based sample.CT scanning provides a reliable means to visualize the joints in rocks,which can be used to verify the reliability of photogrammetry techniques.The application of the photogrammetry-based sample enables detailed analysis for estimating the mechanical properties of jointed rocks.

Numerical simulation of three-dimensional fracturing fracture propagation in radial wells

A fracture propagation model of radial well fracturing is established based on the finite element-meshless method.The model considers the coupling effect of fracturing fluid flow and rock matrix deformation.The fracture geometries of radial well fracturing are simulated,the induction effect of radial well on the fracture is quantitatively characterized,and the influences of azimuth,horizontal principle stress difference,and reservoir matrix permeability on the fracture geometries are revealed.The radial wells can induce the fractures to extend parallel to their axes when two radial wells in the same layer are fractured.When the radial wells are symmetrically distributed along the direction of the minimum horizontal principle stress with the azimuth greater than 15,the extrusion effect reduces the fracture length of radial wells.When the radial wells are symmetrically distributed along the direction of the maximum horizontal principal stress,the extrusion increases the fracture length of the radial wells.The fracture geometries are controlled by the rectification of radial borehole,the extrusion between radial wells in the same layer,and the deflection of the maximum horizontal principal stress.When the radial wells are distributed along the minimum horizontal principal stress symmetrically,the fracture length induced by the radial well decreases with the increase of azimuth;in contrast,when the radial wells are distributed along the maximum horizontal principal stress symmetrically,the fracture length induced by the radial well first decreases and then increases with the increase of azimuth.The fracture length induced by the radial well decreases with the increase of horizontal principal stress difference.The increase of rock matrix permeability and pore pressure of the matrix around radial wells makes the inducing effect of the radial well on fractures increase.

Dynamic simulation insights into friction weakening effect on rapid long-runout landslides:A case study of the Yigong landslide in the Tibetan Plateau,China

This study proposed a novel friction law dependent on velocity,displacement and normal stress for kinematic analysis of runout process of rapid landslides.The well-known Yigong landslide occurring in the Tibetan Plateau of China was employed as the case,and the derived dynamic friction formula was included into the numerical simulation based on Particle Flow Code.Results showed that the friction decreased quickly from 0.64(the peak)to 0.1(the stead value)during the 5s-period after the sliding initiation,which explained the behavior of rapid movement of the landslide.The monitored balls set at different sections of the mass showed similar variation characteritics regarding the velocity,namely evident increase at the initial phase of the movement,followed by a fluctuation phase and then a stopping one.The peak velocity was more than 100 m/s and most particles had low velocities at 300s after the landslide initiation.The spreading distance of the landslide was calculated at the two-dimension(profile)and three-dimension scale,respectively.Compared with the simulation result without considering friction weakening effect,our results indicated a max distance of about 10 km from the initial unstable position,which fit better with the actual situation.

Phase II Clinical Study of Three-Dimensional Printed Coplanar Template Combined with CT-Guided Percutaneous Core Needle Biopsy of Pulmonary Nodules in Elderly Patients

Background: As the population age structure gradually ages, more and more elderly people were found to have pulmonary nodules during physical examinations. Most elderly people had underlying diseases such as heart, lung, brain and blood vessels and cannot tolerate surgery. Computed tomography (CT)-guided percutaneous core needle biopsy (CNB) was the first choice for pathological diagnosis and subsequent targeted drugs, immune drugs or ablation treatment. CT-guided percutaneous CNB requires clinicians with rich CNB experience to ensure high CNB accuracy, but it was easy to cause complications such as pneumothorax and hemorrhage. Three-dimensional (3D) printing coplanar template (PCT) combined with CT-guided percutaneous pulmonary CNB biopsy has been used in clinical practice, but there was no prospective, randomized controlled study. Methods: Elderly patients with lung nodules admitted to the Department of Oncology of our hospital from January 2019 to January 2023 were selected. A total of 225 elderly patients were screened, and 30 patients were included after screening. They were randomly divided into experimental group (Group A: 30 cases) and control group (Group B: 30 cases). Group A was given 3D-PCT combined with CT-guided percutaneous pulmonary CNB biopsy, Group B underwent CT-guided percutaneous pulmonary CNB. The primary outcome measure of this study was the accuracy of diagnostic CNB, and the secondary outcome measures were CNB time, number of CNB needles, number of pathological tissues and complications. Results: The diagnostic accuracy of group A and group B was 96.67% and 76.67%, respectively (P = 0.026). There were statistical differences between group A and group B in average CNB time (P = 0.001), number of CNB (1 vs more than 1, P = 0.029), and pathological tissue obtained by CNB (3 vs 1, P = 0.040). There was no statistical difference in the incidence of pneumothorax and hemorrhage between the two groups (P > 0.05). Conclusions: 3D-PCT combined with CT-guided percutaneous CNB can improve the puncture accuracy of elderly patients, shorten the puncture time, reduce the number of punctures, and increase the amount of puncture pathological tissue, without increasing pneumothorax and hemorrhage complications. We look forward to verifying this in a phase III randomized controlled clinical study. .

Three-Dimensional Simulation of Hydrodynamic Mechanism of Fluidized Bed Methanation

Organic solid waste(OSW)contains many renewable materials.The pyrolysis and gasification of OSW can realize resource utilization,and its products can be used for methanation reaction to produce synthetic natural gas in the specific reactor.In order to understand the dynamic characteristics of the reactor,a three-dimensional numerical model has been established by the method of Computational Fluid Dynamics(CFD).Along the height of the reactor,the particle distribution in the bed becomes thinner and the mean solid volume fraction decreases from 4.18%to 0.37%.Meanwhile,the pressure fluctuation range decreased from 398.76 Pa at the entrance to a much lower value of 74.47 Pa at the exit.In this simulation,three parameters of gas inlet velocity,operating temperature and solid particle diameter are changed to explore their influences on gas-solid multiphase flow.The results show that gas velocity has a great influence on particle distribution.When the gas inlet velocity decreases from 6.51 to 1.98 m/s,the minimum height that particles can reach decreases from 169 to 100 mm.Additionally,as the operating temperature increases,the particle holdup inside the reactor changes from 0.843%to 0.700%.This indicates that the particle residence time reduces,which is not conducive to the follow-up reaction.Moreover,with the increase of particle size,the fluctuation range of the pressure at the bottom of the reactor increases,and its standard deviation increases from 55.34 to 1266.37 Pa.

Three-Dimensional and Cross-sectional Characteristics of Normal Grain Growth Based on Monte Carlo Simulation

An appropriate Monte Carlo method was developed to simulate the three-dimensional normal grain growth more completely. Comparative investigation on the three-dimensional and the cross-sectional characteristics of normal grain growth was done. It was found that the time exponent of grain growth determined from cross-section exhibits the same rule of increasing slowly with time and approaching the theoretical value n = 0.5 of steadygrain growth as the three-dimensional (3-D) system. From change of the number of grains per unit area with timemeasured in cross-section, the state of 3-D normal grain growth may be predicted. The gtain size distribution incross-section is different from that in 3-D system and can not express the evolution characteristic of the 3-D distribution. Furthermore, there exists statistical connection between the topological parameters in cross-section and thosein three-dimensions.

Three-Dimensional Finite Element Numerical Simulation and Physical Experiment for Magnetism-Stress Detecting in Oil Casing

The casing damage has been a big problem in oilfield production. The current detection methods mostly are used after casing damage, which is not very effective. With the rapid development of China's offshore oil industry, the number of offshore oil wells is becoming larger and larger. Because the cost of offshore oil well is very high, the casing damage will cause huge economic losses. What's more, it can also bring serious pollution to marine environment. So the effective methods of detecting casing damage are required badly. The accumulation of stress is the main reason for the casing damage. Magnetic anisotropy technique based on counter magnetostriction effect can detect the stress of casing in real time and help us to find out the hidden dangers in time. It is essential for us to prevent the casing damage from occurring. However, such technique is still in the development stage. Previous studies mostly got the relationship between stress and magnetic signals by physical experiment, and the study of physical mechanism in relative magnetic permeability connecting the stress and magnetic signals is rarely reported. The present paper uses the ANSYS to do the three-dimensional finite element numerical simulation to study how the relative magnetic permeability works for the oil casing model. We find that the quantitative relationship between the stress' s variation and magnetic induction intensity's variation is: Δδ =K* ΔB, K = 8.04×109, which is proved correct by physical experiment.

Three-dimensional discrete element numerical simulation of Paleogene salt structures in the western Kuqa foreland thrust belt

Taking the Paleogene salt strata in the west of Kuqa foreland thrust belt as study object, the deformation features of salt structure in the compression direction and perpendicular to the compression direction were examined to find out the control factors and formation mechanisms of the salt structures. By using the three-dimensional discrete element numerical simulation method, the formation mechanisms of typical salt structures of western Kuqa foreland thrust belt in Keshen and Dabei work areas were comprehensively analyzed. The simulation results show that the salt deformation in Keshen and Dabei work areas is of forward spread type, with deformation concentrated in the piedmont zone;the salt deformation is affected by the early uplift near the compression end, pre-existing basement faults, synsedimentary process and the initial salt depocenter;in the direction perpendicular to the compression direction, salt rocks near the compression end have strong lateral mobility with the velocity component moving towards the middle part, and the closer to the middle, the larger the velocity will be, so that salt rocks will aggregate towards the middle and deform intensely, forming complex folds and separation of salt structures from salt source, and local outcrop with thrust faults. Compared with 2 D simulation, 3 D simulation can analyze salt structures in the principal stress direction and direction perpendicular to the principal stress, give us a full view of the formation mechanisms of salt structures, and guide the exploration of oil and gas reservoirs related to salt structures.

Three-Dimensional Rigid-Plastic FEM Simulation of Bulk Forming Processes with New Contact and Remeshing Techniques

Some techniques such as die surface description, contact judgement algorithm and remeshing are proposed to improve the robustness of the numerical solution. Based on these techniques, a three-dimensional rigid-plastic FEM code has been developed. Isothermal forging process of a cylindrical housing has been simulated. The simulation results show that the given techniques and the FEM code are reasonable and feasible for three-dimensional bulk forming processes.

Numerical simulations and comparative analysis of two- and three-dimensional circulating fluidized bed reactors for CO2 capture

Carbon dioxide(CO2),the main gas emitted from fossil burning,is the primary contributor to global warming.Circulating fluidized bed reactor(CFBR)is proved as an energy-efficient method for post-combustion CO2 capture.The numerical simulation by computational fluid dynamics(CFD)is believed as a promising tool to study CO2 adsorption process in CFBR.Although three-dimensional(3D)simulations were proved to have better predicting performance with the experimental results,two-dimensional(2D)simulations have been widely reported for qualitative and quantitative studies on gas-solid behavior in CFBR for its higher computational efficiency recently.However,the discrepancies between 2D and 3D simulations have rarely been evaluated by detailed study.Considering that the differences between the 2D and 3D simulations will vary substantially with the changes of independent operating conditions,it is beneficial to lower computational costs to clarify the effects of dimensionality on the numerical CO2 adsorption runs under various operating conditions.In this work,the comparative analysis for CO2 adsorption in 2D and 3D simulations was conducted to enlighten the effects of dimensionality on the hydrodynamics and reaction behaviors,in which the separation rate,species distribution and hydrodynamic characteristics were comparatively studied for both model frames.With both accuracy and computational costs considered,the viable suggestions were provided in selecting appropriate model frame for the studies on optimization of operating conditions,which directly affect the capture and energy efficiencies of cyclic CO2 capture process in CFBR.

Simulation-Based Construction of Three-Dimensional Process Model for Punching Cartridge Cases

A method of constructing three-dimensional process model for the punching cartridge cases is presented based on DEFORM simulation analysis. Using DEFORM software,the finite element simulation models for the punching and forming process of cartridge cases are established,and the corresponding simulation result model of each intermediate procedure is obtained by continuously performing the forming process simulation. The simulation model cannot annotate size and process information due to poor interface between DEFORM software and CAD software. Thus,a 3D annotation module is developed with secondary development technology of UG NX software. Consequently,the final process model with dimension and process information is obtained. Then,with the current 3D process management system,the 3D punching and forming process design of cartridge cases can be completed further. An example is also provided to illustrate that the relative error between the simulation process model and the physical model is less than 2%,which proves the validity and reliability of the proposed method in this study.

Numerical Simulation of the Three-Dimensional Wave-Induced Currents on Unstructured Grid

By coupling the three-dimensional hydrodynamic model with the wave model, numerical simulations of the three- dimensional wave-induced current are carried out in this study. The wave model is based on the numerical solution of the modified wave action equation and eikonal equation, which can describe the wave refraction and diffraction. The hydrodynamic model is driven by the wave-induced radiation stresses and affected by the wave turbulence. The numerical implementation of the module has used the finite-volume schemes on unstructured grid, which provides great flexibility for modeling the waves and currents in the complex actual nearshore, and ensures the conservation of energy propagation. The applicability of the proposed model is evaluated in calculating the cases of wave set-up, longshore currents, undertow on a sloping beach, rip currents and meandering longshore currents on a tri-cuspate beach. The results indicate that it is necessary to introduce the depth-dependent radiation stresses into the numerical simulation of wave-induced currents, and comparisons show that the present model makes better prediction on the wave procedure as well as both horizontal and vertical structures in the wave-induced current field.

A three-dimensional Eulerian method for the numerical simulation of high-velocity impact problems

In the present paper, a three-dimensional （3D） Eulerian technique for the 3D numerical simulation of high-velocity impact problems is proposed. In the Eulerian framework, a complete 3D conservation element and solution element scheme for conservative hyperbolic governing equations with source terms is given. A modified ghost fluid method is proposed for the treatment of the boundary conditions. Numerical simulations of the Taylor bar problem and the ricochet phenomenon of a sphere impacting a plate target at an angle of 60~ are carried out. The numerical results are in good agreement with the corresponding experimental observations. It is proved that our computational technique is feasible for analyzing 3D high-velocity impact problems.

Application and Analysis of Computer Virtual Simulation Technology in Three-dimensional Animation Production

With the continuous promotion of computer technology, the application system of virtual simulation technology has been further optimized and improved, and has been widely used in various fields of social development, such as urban construction, interior design, industrial simulation and tourism teaching. China's three-dimensional animation production started relatively late, but has achieved good results with the support of related advanced technology in the process of development. Computer virtual simulation technology is an important technical support in the production of three-dimensional animation. In this paper, firstly, the related content of computer virtual simulation technology was introduced. Then, the specific application of this technology in the production of three-dimensional animation was further elaborated, so as to provide some reference for the improvement of the production effect of three-dimensional animation in the future.

Three-dimensional Monte Carlo simulation of bulk fin field effect transistor

In this paper, we investigate the performance of the bulk fin field effect transistor （FinFET） through a three- dimensional （3D） full band Monte Carlo simulator with quantum correction. Several scattering mechanisms, such as the acoustic and optical phonon scattering, the ionized impurity scattering, the impact ionization scattering and the surface roughness scattering are considered in our simulator. The effects of the substrate bias and the surface roughness scattering near the Si/SiO2 interface on the performance of bulk FinFET are mainly discussed in our work. Our results show that the on-current of bulk FinFET is sensitive to the surface roughness and that we can reduce the substrate leakage current by modulating the substrate bias voltage.

An online physical simulation method for enhanced oil recovery by air injection in shale oil

In order to understand the mechanism of air flooding shale oil, an online physical simulation method for enhanced shale oil recovery by air injection was established by integrating CT scanning and nuclear magnetic resonance(NMR). The development effect of shale oil by air flooding under different depletion pressures, the micro-production characteristics of pore throats with different sizes and the mechanism of shale oil recovery by air flooding were analyzed. The effects of air oxygen content, permeability, gas injection pressure, and fractures on the air flooding effect in shale and crude oil production in pores with different sizes were analyzed. The recovery of shale oil can be greatly improved by injecting air into the depleted shale reservoir, but the oil displacement efficiency and the production degree of different levels of pore throats vary with the injection timing. The higher the air oxygen content and the stronger the low-temperature oxidation, the higher the production degree of pores with different sizes and the higher the shale oil recovery. The higher the permeability and the better the pore throat connectivity, the stronger the fluid flow capacity and the higher the shale oil recovery. As the injection pressure increases, the lower limit of the production degree of pore throats decreases, but gas channeling may occur to cause a premature breakthrough;as a result, the recovery increases and then decreases. Fractures can effectively increase the contact area between gas and crude oil, and increase the air sweep coefficient and matrix oil drainage area by supplying oil to fractures through the matrix, which means that a proper fracturing before air injection can help to improve the oil displacement effect under a reasonable production pressure difference.

Numerical simulation of asphalt mixture based on three-dimensional heterogeneous specimen

In order to verify the validity of finite element numerical simulation method for asphalt mixture, which consists of aggregates, mastic （where mastic is a kind of fine mixture composed of asphalt binder mixed with fines and fine aggregates） and air voids, based on three-dimensional （3D） heterogeneous specimen, X-ray computerized tomography （X-ray CT） was used to scan the asphalt specimens to obtain the real internal microstrnctures of asphalt mixture. CT images were reconstructed to build up 3D digital specimen, and the viscoelastic properties of mastic were described with Burgers model The uniaxial creep numerical simulations of three different levels of aggregate gradation were conducted. The simulation results agree well with the experimental results. The numerical simulation of asphalt mixture incorporated with real 3D microstructure based on finite element method is a promising application to conduct research of asphalt concrete. Additionally, this method can increase the mechanistic understanding of global viscoelastic properties of asphalt mixtures by linking the real 3D microstructure.

Three-dimensional visualization and virtual reality simulation role in hepatic surgery:Further research warranted

Artificial intelligence(AI)is the study of algorithms that enable machines to analyze and execute cognitive activities including problem solving,object and word recognition,reduce the inevitable errors to improve the diagnostic accuracy,and decision-making.Hepatobiliary procedures are technically complex and the use of AI in perioperative management can improve patient outcomes as discussed below.Three-dimensional(3D)reconstruction of images obtained via ultrasound,computed tomography scan or magnetic resonance imaging,can help surgeons better visualize the surgical sites with added depth perception.Preoperative 3D planning is associated with lesser operative time and intraoperative complications.Also,a more accurate assessment is noted,which leads to fewer operative complications.Images can be converted into physical models with 3D printing technology,which can be of educational value to students and trainees.3D images can be combined to provide 3D visualization,which is used for preoperative navigation,allowing for more precise localization of tumors and vessels.Nevertheless,AI enables surgeons to provide better,personalized care for each patient.

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.

Three-Dimensional Simulations of RESET Operation in Phase-Change Random Access Memory with Blade-Type Like Phase Change Layer by Finite Element Modeling

An optimized device structure for reducing the RESET current of phase-change random access memory （PCRAM） with blade-type like （BTL） phase change layer is proposed. The electrical thermal analysis of the BTL cell and the blade heater contactor structure by three-dimensional finite element modeling are compared with each other during RESET operation. The simulation results show that the programming region of the phase change layer in the BTL cell is much smaller, and thermal electrical distributions of the BTL cell are more concentrated on the TiN/GST interface. The results indicate that the BTL cell has the superiorities of increasing the heating efficiency, decreasing the power consumption and reducing the RESET current from 0.67mA to 0.32mA. Therefore, the BTL cell will be appropriate for high performance PCRAM device with lower power consumption and lower RESET current.