Gobi spans a large area of China,surpassing the combined expanse of mobile dunes and semi-fixed dunes.Its presence significantly influences the movement of sand and dust.However,the complex origins and diverse materia...Gobi spans a large area of China,surpassing the combined expanse of mobile dunes and semi-fixed dunes.Its presence significantly influences the movement of sand and dust.However,the complex origins and diverse materials constituting the Gobi result in notable differences in saltation processes across various Gobi surfaces.It is challenging to describe these processes according to a uniform morphology.Therefore,it becomes imperative to articulate surface characteristics through parameters such as the three-dimensional(3D)size and shape of gravel.Collecting morphology information for Gobi gravels is essential for studying its genesis and sand saltation.To enhance the efficiency and information yield of gravel parameter measurements,this study conducted field experiments in the Gobi region across Dunhuang City,Guazhou County,and Yumen City(administrated by Jiuquan City),Gansu Province,China in March 2023.A research framework and methodology for measuring 3D parameters of gravel using point cloud were developed,alongside improved calculation formulas for 3D parameters including gravel grain size,volume,flatness,roundness,sphericity,and equivalent grain size.Leveraging multi-view geometry technology for 3D reconstruction allowed for establishing an optimal data acquisition scheme characterized by high point cloud reconstruction efficiency and clear quality.Additionally,the proposed methodology incorporated point cloud clustering,segmentation,and filtering techniques to isolate individual gravel point clouds.Advanced point cloud algorithms,including the Oriented Bounding Box(OBB),point cloud slicing method,and point cloud triangulation,were then deployed to calculate the 3D parameters of individual gravels.These systematic processes allow precise and detailed characterization of individual gravels.For gravel grain size and volume,the correlation coefficients between point cloud and manual measurements all exceeded 0.9000,confirming the feasibility of the proposed methodology for measuring 3D parameters of individual gravels.The proposed workflow yields accurate calculations of relevant parameters for Gobi gravels,providing essential data support for subsequent studies on Gobi environments.展开更多
Understanding the pore water pressure distribution in unsaturated soil is crucial in predicting shallow landslides triggered by rainfall,mainly when dealing with different temporal patterns of rainfall intensity.Howev...Understanding the pore water pressure distribution in unsaturated soil is crucial in predicting shallow landslides triggered by rainfall,mainly when dealing with different temporal patterns of rainfall intensity.However,the hydrological response of vegetated slopes,especially three-dimensional(3D)slopes covered with shrubs,under different rainfall patterns remains unclear and requires further investigation.To address this issue,this study adopts a novel 3D numerical model for simulating hydraulic interactions between the root system of the shrub and the surrounding soil.Three series of numerical parametric studies are conducted to investigate the influences of slope inclination,rainfall pattern and rainfall duration.Four rainfall patterns(advanced,bimodal,delayed,and uniform)and two rainfall durations(4-h intense and 168-h mild rainfall)are considered to study the hydrological response of the slope.The computed results show that 17%higher transpiration-induced suction is found for a steeper slope,which remains even after a short,intense rainfall with a 100-year return period.The extreme rainfalls with advanced(PA),bimodal(PB)and uniform(PU)rainfall patterns need to be considered for the short rainfall duration(4 h),while the delayed(PD)and uniform(PU)rainfall patterns are highly recommended for long rainfall durations(168 h).The presence of plants can improve slope stability markedly under extreme rainfall with a short duration(4 h).For the long duration(168 h),the benefit of the plant in preserving pore-water pressure(PWP)and slope stability may not be sufficient.展开更多
The 2D limit equilibrium method is widely used for slope stability analysis.However,with the advancement of dump engineering,composite slopes often exhibit significant 3D mechanical effects.Consequently,it is of signi...The 2D limit equilibrium method is widely used for slope stability analysis.However,with the advancement of dump engineering,composite slopes often exhibit significant 3D mechanical effects.Consequently,it is of significant importance to develop an effective 3D stability calculation method for composite slopes to enhance the design and stability control of open-pit slope engineering.Using the composite slope formed by the mining stope and inner dump in Baiyinhua No.1 and No.2 open-pit coal mine as a case study,this research investigates the failure mode of composite slopes and establishes spatial shape equations for the sliding mass.By integrating the shear resistance and sliding force of each row of microstrip columns onto the bottom surface of the strip corresponding to the main sliding surface,a novel 2D equivalent physical and mechanical parameters analysis method for the strips on the main sliding surface of 3D sliding masses is proposed.Subsequently,a comprehensive 3D stability calculation method for composite slopes is developed,and the quantitative relationship between the coordinated development distance and its 3D stability coefficients is examined.The analysis reveals that the failure mode of the composite slope is characterized by cutting-bedding sliding,with the arc serving as the side interface and the weak layer as the bottom interface,while the destabilization mechanism primarily involves shear failure.The spatial form equation of the sliding mass comprises an ellipsoid and weak plane equation.The analysis revealed that when the coordinated development distance is 1500 m,the error rate between the 3D stability calculation result and the 2D stability calculation result of the composite slope is less than 8%,thereby verifying the proposed analytical method of equivalent physical and mechanical parameters and the 3D stability calculation method for composite slopes.Furthermore,the3D stability coefficient of the composite slope exhibits an exponential correlation with the coordinated development distance,with the coefficient gradually decreasing as the coordinated development distance increases.These findings provide a theoretical guideline for designing similar slope shape parameters and conducting stability analysis.展开更多
The important supporting component in a gas turbine is the casing,which has the characteristics of large size,complex structure,and thin wall.In the context of existing 3DP sand casting processes,casting crack defects...The important supporting component in a gas turbine is the casing,which has the characteristics of large size,complex structure,and thin wall.In the context of existing 3DP sand casting processes,casting crack defects are prone to occur.This leads to an increase in the scrap rate of casings,causing significant resource wastage.Additionally,the presence of cracks poses a significant safety hazard after the casings are put into service.The generation of different types of crack defects in stainless steel casings is closely related to casting stress and the high-temperature concession of the sand mold.Therefore,the types and causes of cracks in stainless steel casing products,based on their structural characteristics,were systematically analyzed.Various sand molds with different internal topology designs were printed using the 3DP technology to investigate the impact of sand mold structures on high-temperature concession.The optimal sand mold structure was used to cast casings,and the crack suppression effect was verified by analyzing its eddy current testing results.The experimental results indicate that the skeleton structure has an excellent effect on suppressing cracks in the casing.This research holds important theoretical and engineering significance in improving the quality of casing castings and reducing production costs.展开更多
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, lu...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. .展开更多
The heterogeneity of unconventional reservoir rock tremendously affects its hydrofracturing behavior. A visual representation and accurate characterization of the three-dimensional (3D) growth and distribution of hy...The heterogeneity of unconventional reservoir rock tremendously affects its hydrofracturing behavior. A visual representation and accurate characterization of the three-dimensional (3D) growth and distribution of hydrofracturing cracks within heterogeneous rocks is of particular use to the design and implementation of hydrofracturing stimulation of unconventional reservoirs. However, because of the difficulties involved in visually representing and quantitatively characterizing a 3D hydrofracturing crack-network, this issue remains a challenge. In this paper, a novel method is proposed for physically visualizing and quantitatively characterizing the 3D hydrofracturing crack-network distributed through a heterogeneous structure based on a natural glutenite sample. This method incorporates X-ray microfocus computed tomography (μCT), 3D printing models and hydrofracturing triaxial tests to represent visually the heterogeneous structure, and the 3D crack growth and distribution within a transparent rock model during hydrofracturing. The coupled effects of material heterogeneity and confining geostress on the 3D crack initiation and propagation were analyzed. The results indicate that the breakdown pressure of a heterogeneous rock model is significantly affected by material heterogeneity and confining geostress. The measured breakdown pressures of heterogeneous models are apparently different from those predicted by traditional theories. This study helps to elucidate the quantitative visualization and characterization of the mechanism and influencing factors that determine the hydrofracturing crack initiation and propagation in heterogeneous reservoir rocks.展开更多
To investigate the long-term stability of deep rocks,a three-dimensional(3D)time-dependent model that accounts for excavation-induced damage and complex stress state is developed.This model comprises three main compon...To investigate the long-term stability of deep rocks,a three-dimensional(3D)time-dependent model that accounts for excavation-induced damage and complex stress state is developed.This model comprises three main components:a 3D viscoplastic isotropic constitutive relation that considers excavation damage and complex stress state,a quantitative relationship between critical irreversible deformation and complex stress state,and evolution characteristics of strength parameters.The proposed model is implemented in a self-developed numerical code,i.e.CASRock.The reliability of the model is validated through experiments.It is indicated that the time-dependent fracturing potential index(xTFPI)at a given time during the attenuation creep stage shows a negative correlation with the extent of excavationinduced damage.The time-dependent fracturing process of rock demonstrates a distinct interval effect of the intermediate principal stress,thereby highlighting the 3D stress-dependent characteristic of the model.Finally,the influence of excavation-induced damage and intermediate principal stress on the time-dependent fracturing characteristics of the surrounding rocks around the tunnel is discussed.展开更多
A technique for modelling of three-dimensional(3D)quasi-statically propagating cracks in elastic bodies by the displacement discontinuity method(DDM)was described.When the crack is closed,the Mohr-coulomb rule on the ...A technique for modelling of three-dimensional(3D)quasi-statically propagating cracks in elastic bodies by the displacement discontinuity method(DDM)was described.When the crack is closed,the Mohr-coulomb rule on the two contacted surfaces of the crack must be satisfied.A simple iterative method was adopted in order to consider three different states of cracks.Under the assumption that the advance of the point on the crack front would occur only in the normal plane which is through this edge point,the maximum energy release rate criterion is modified to be used as the criterion for the crack growth.With discretization,the process of crack propagation can be seen as the advance of the vertices of the crack front.The program MCP3D was developed based on these theories to simulate the 3D quasi-static crack propagation.A numerical example of a penny-shaped crack subject to tension and compression in an infinite elastic media was analyzed with MCP3D,and the results in comparison with others' show that the present method for 3D crack propagation is effective.展开更多
To achieve the loading of the stress path of hard rock,the spherical discrete element model(DEM)and the new flexible membrane technology were utilized to realize the transient loading of three principal stresses with ...To achieve the loading of the stress path of hard rock,the spherical discrete element model(DEM)and the new flexible membrane technology were utilized to realize the transient loading of three principal stresses with arbitrary magnitudes and orientations.Furthermore,based on the deep tunnel of China Jinping Underground Laboratory II(CJPL-II),the deformation and fracture evolution characteristics of deep hard rock induced by excavation stress path were analyzed,and the mechanisms of transient loading-unloading and stress rotation-induced fractures were revealed from a mesoscopic perspective.The results indicated that the stressestrain curve exhibits different trends and degrees of sudden changes when subjected to transient changes in principal stress,accompanied by sudden changes in strain rate.Stress rotation induces spatially directional deformation,resulting in fractures of different degrees and orientations,and increasing the degree of deformation anisotropy.The correlation between the degree of induced fracture and the unloading magnitude of minimum principal stress,as well as its initial level is significant and positive.The process of mechanical response during transient unloading exhibits clear nonlinearity and directivity.After transient unloading,both the minimum principal stress and minimum principal strain rate decrease sharply and then tend to stabilize.This occurs from the edge to the interior and from the direction of the minimum principal stress to the direction of the maximum principal stress on theε1-ε3 plane.Transient unloading will induce a tensile stress wave.The ability to induce fractures due to changes in principal stress magnitude,orientation and rotation paths gradually increases.The analysis indicates a positive correlation between the abrupt change amplitude of strain rate and the maximum unloading magnitude,which is determined by the magnitude and rotation of principal stress.A high tensile strain rate is more likely to induce fractures under low minimum principal stress.展开更多
Weak structural plane deformation is responsible for the non-uniform large deformation disasters in layered rock tunnels,resulting in steel arch distortion and secondary lining cracking.In this study,a servo biaxial t...Weak structural plane deformation is responsible for the non-uniform large deformation disasters in layered rock tunnels,resulting in steel arch distortion and secondary lining cracking.In this study,a servo biaxial testing system was employed to conduct physical modeling tests on layered rock tunnels with bedding planes of varying dip angles.The influence of structural anisotropy in layered rocks on the micro displacement and strain field of surrounding rocks was analyzed using digital image correlation(DIC)technology.The spatiotemporal evolution of non-uniform deformation of surrounding rocks was investigated,and numerical simulation was performed to verify the experimental results.The findings indicate that the displacement and strain field of the surrounding layered rocks are all maximized at the horizontal bedding planes and decrease linearly with the increasing dip angle.The failure of the layered surrounding rock with different dip angles occurs and extends along the bedding planes.Compressive strain failure occurs after excavation under high horizontal stress.This study provides significant theoretical support for the analysis,prediction,and control of non-uniform deformation of tunnel surrounding rocks.展开更多
Three-dimensional(3D) synthetic aperture radar(SAR)extends the conventional 2D images into 3D features by several acquisitions in different aspects. Compared with 3D techniques via multiple observations in elevation, ...Three-dimensional(3D) synthetic aperture radar(SAR)extends the conventional 2D images into 3D features by several acquisitions in different aspects. Compared with 3D techniques via multiple observations in elevation, e.g. SAR interferometry(InSAR) and SAR tomography(TomoSAR), holographic SAR can retrieve 3D structure by observations in azimuth. This paper focuses on designing a novel type of orbit to achieve SAR regional all-azimuth observation(AAO) for embedded targets detection and holographic 3D reconstruction. The ground tracks of the AAO orbit separate the earth surface into grids. Target in these grids can be accessed with an azimuth angle span of360°, which is similar to the flight path of airborne circular SAR(CSAR). Inspired from the successive coverage orbits of optical sensors, several optimizations are made in the proposed method to ensure favorable grazing angles, the performance of 3D reconstruction, and long-term supervision for SAR sensors. Simulation experiments show the regional AAO can be completed within five hours. In addition, a second AAO of the same area can be duplicated in two days. Finally, an airborne SAR data process result is presented to illustrate the significance of AAO in 3D reconstruction.展开更多
The shear failure of intact rock under thermo-mechanical(TM)coupling conditions is common,such as in enhanced geothermal mining and deep mine construction.Under the effect of a continuous engineering disturbance,shear...The shear failure of intact rock under thermo-mechanical(TM)coupling conditions is common,such as in enhanced geothermal mining and deep mine construction.Under the effect of a continuous engineering disturbance,shear-formed fractures are prone to secondary instability,posing a severe threat to deep engineering.Although numerous studies regarding three-dimensional(3D)morphologies of fracture surfaces have been conducted,the understanding of shear-formed fractures under TM coupling conditions is limited.In this study,direct shear tests of intact granite under various TM coupling conditions were conducted,followed by 3D laser scanning tests of shear-formed fractures.Test results demonstrated that the peak shear strength of intact granite is positively correlated with the normal stress,whereas it is negatively correlated with the temperature.The internal friction angle and cohesion of intact granite significantly decrease with an increase in the temperature.The anisotropy,roughness value,and height of the asperities on the fracture surfaces are reduced as the normal stress increases,whereas their variation trends are the opposite as the temperature increases.The macroscopic failure mode of intact granite under TM coupling conditions is dominated by mixed tensileeshear and shear failures.As the normal stress increases,intragranular fractures are developed ranging from a local to a global distribution,and the macroscopic failure mode of intact granite changes from mixed tensileeshear to shear failure.Finally,3D morphological characteristics of the asperities on the shear-formed fracture surfaces were analyzed,and a quadrangular pyramid conceptual model representing these asperities was proposed and sufficiently verified.展开更多
This study is concerned with the three-dimensional(3D)stagnation-point for the mixed convection flow past a vertical surface considering the first-order and secondorder velocity slips.To the authors’knowledge,this is...This study is concerned with the three-dimensional(3D)stagnation-point for the mixed convection flow past a vertical surface considering the first-order and secondorder velocity slips.To the authors’knowledge,this is the first study presenting this very interesting analysis.Nonlinear partial differential equations for the flow problem are transformed into nonlinear ordinary differential equations(ODEs)by using appropriate similarity transformation.These ODEs with the corresponding boundary conditions are numerically solved by utilizing the bvp4c solver in MATLAB programming language.The effects of the governing parameters on the non-dimensional velocity profiles,temperature profiles,skin friction coefficients,and the local Nusselt number are presented in detail through a series of graphs and tables.Interestingly,it is reported that the reduced skin friction coefficient decreases for the assisting flow situation and increases for the opposing flow situation.The numerical computations of the present work are compared with those from other research available in specific situations,and an excellent consensus is observed.Another exciting feature for this work is the existence of dual solutions.An important remark is that the dual solutions exist for both assisting and opposing flows.A linear stability analysis is performed showing that one solution is stable and the other solution is not stable.We notice that the mixed convection and velocity slip parameters have strong effects on the flow characteristics.These effects are depicted in graphs and discussed in this paper.The obtained results show that the first-order and second-order slip parameters have a considerable effect on the flow,as well as on the heat transfer characteristics.展开更多
Similar to hydraulic fracturing(HF), the coalescence and fracture of cracks are induced within a rock under the action of an ultrasonic field, known as ultrasonic fracturing(UF). Investigating UF is important in both ...Similar to hydraulic fracturing(HF), the coalescence and fracture of cracks are induced within a rock under the action of an ultrasonic field, known as ultrasonic fracturing(UF). Investigating UF is important in both hard rock drilling and oil and gas recovery. A three-dimensional internal laser-engraved crack(3D-ILC) method was introduced to prefabricate two parallel internal cracks within the samples without any damage to the surface. The samples were subjected to UF. The mechanism of UF was elucidated by analyzing the characteristics of fracture surfaces. The crack propagation path under different ultrasonic parameters was obtained by numerical simulation based on the Paris fatigue model and compared to the experimental results of UF. The results show that the 3D-ILC method is a powerful tool for UF research.Under the action of an ultrasonic field, the fracture surface shows the characteristics of beach marks and contains powder locally, indicating that the UF mechanism includes high-cycle fatigue fracture, shear and friction, and temperature load. The two internal cracks become close under UF. The numerical result obtained by the Paris fatigue model also shows the attraction of the two cracks, consistent with the test results. The 3D-ILC method provides a new tool for the experimental study of UF. Compared to the conventional numerical methods based on the analysis of stress-strain and plastic zone, numerical simulation can be a good alternative method to obtain the crack path under UF.展开更多
The gravity inversion is to restore genetic density distribution of the underground target to be explored for explaining the internal structure and distribution of the Earth.In this paper,we propose a new 3D gravity i...The gravity inversion is to restore genetic density distribution of the underground target to be explored for explaining the internal structure and distribution of the Earth.In this paper,we propose a new 3D gravity inversion method based on 3D U-Net++.Compared with two-dimensional gravity inversion,three-dimensional(3D)gravity inversion can more precisely describe the density distribution of underground space.However,conventional 3D gravity inversion method input is two-dimensional,the input and output of the network proposed in our method are three-dimensional.In the training stage,we design a large number of diversifi ed simulation model-data pairs by using the random walk method to improve the generalization ability of the network.In the test phase,we verify the network performance by using the model-data pairs generated by the simulation.To further illustrate the eff ectiveness of the algorithm,we apply the method to the inversion of the San Nicolas mining area,and the inversion results are basically consistent with the borehole measurement results.Moreover,the results of the 3D U-Net++inversion and the 3D U-Net inversion are compared.The density models of the 3D U-Net++inversion have higher resolution,more concentrated inversion results,and a clearer boundary of the density model.展开更多
Due to high cost of full-scale experimental setup, this study presents a numerical model on fatigue behaviours of offshore pipeline with multiple coplanar cracks under cyclic tensile loadings. The validation on numeri...Due to high cost of full-scale experimental setup, this study presents a numerical model on fatigue behaviours of offshore pipeline with multiple coplanar cracks under cyclic tensile loadings. The validation on numerical results is made by other researchers' experimental results, and significant parameters affecting fatigue crack growth are studied.展开更多
The three-dimensional (3D) crack propagation is a hot issue in rock mechanics. To properly simulate 3D crack propagation, a modified maximum tangential tensile stress criterion is proposed. In this modified criterio...The three-dimensional (3D) crack propagation is a hot issue in rock mechanics. To properly simulate 3D crack propagation, a modified maximum tangential tensile stress criterion is proposed. In this modified criterion, it is supposed that cracks propagate only at crack front in the principal normal plane. The tangential tensile stress at crack front in the principal normal plane in local coordinates is employed to determine crack propagation, which is calculated through coordinate transformation from global to local coordinates. New cracks will propagate when the maximum tangential tensile stress at crack front in the principal normal plane reaches the tensile strength of rock-like materials. Compared with the previous crack propagation criteria, the modified crack propagation criterion is helpful in calculating 3D crack stress intensity factor, and can overcome the limitations of propagation step determined by individual experiences in previous studies. Finally, the 3D crack propagation process is traced by element-free Galerkin method. The numerical results agree well with the experimental ones for a frozen resin sample with prefabricated 3D cracks.展开更多
Adhesions between different cells and extracellular matrix have been studied extensively in vitro, but little is known about their functions in testicular tissue counterparts. Spermatogonia and their companion somatic...Adhesions between different cells and extracellular matrix have been studied extensively in vitro, but little is known about their functions in testicular tissue counterparts. Spermatogonia and their companion somatic cells maintain a close association throughout spermatogenesis and this association is necessary for normal spermatogenesis. In order to keep the relative integrity of the testicular tissues, and to detect the development in vitro, culture testicular tissues in a three- dimensional (3D) agarose matrix was examined. Testicular tissues isolated from 6.5 d postpartum (dpp) mouse were cultured on the top of the matrix for 26 d with a medium height up to 4/5 of the 3D agarose matrix. The results showed that in this 3D culture environment, each type of testicular cells kept the same structure, localization and function as in vivo and might be more biologically relevant to living organisms. After culture, germ cell marker VASA and meiosis markers DAZL and SCP3 showed typical positive analysed by immunofluorescence staining and RT-PCR. It demonstrated that this 3D culture system was able to maintain the number of germ cells and promote the meiosis initiation of male germ cells.展开更多
The three-dimensional finite element method of lines is presented, and the basic processing description of 3D FEMOL in cracking questions is given in detail. Applications to 3D bodies with cracks indicate that good ac...The three-dimensional finite element method of lines is presented, and the basic processing description of 3D FEMOL in cracking questions is given in detail. Applications to 3D bodies with cracks indicate that good accuracy can be obtained with relatively coarse girds. In particular, application to the tension specimen shows very good agreement with the evaluation of stress intensity factors, which is better than the results of other methods. This implies a considerable potential for using this method in the 3D analysis of finite geometry solids and suggests a possible extension of this technique to nonlinear material behavior.展开更多
Upon the conservation of mass, momentum and energy, volume fraction and surface penetrative rate were employed to modify the conservative equations to simulate the effect of blockages on fluid flows and heat transfer....Upon the conservation of mass, momentum and energy, volume fraction and surface penetrative rate were employed to modify the conservative equations to simulate the effect of blockages on fluid flows and heat transfer. These equations were solved numerically with the finite differential method and the primitive variable approach. This method uses staggered grid and pressure correction schemes. A computer code FASTOR3D integrated the aforementioned algorithm. The preliminary results have been compared with conventional benchmark solutions. With auxiliary software DV, the numerical results were visualized in colorful images to demonstrate the variation of flow patterns and temperature profiles during the transient process. The results of the simulation code for the fluid flows and heat transfer in the sodium pool of a fast breeder reactor are acceptable.展开更多
基金funded by the National Natural Science Foundation of China(42071014).
文摘Gobi spans a large area of China,surpassing the combined expanse of mobile dunes and semi-fixed dunes.Its presence significantly influences the movement of sand and dust.However,the complex origins and diverse materials constituting the Gobi result in notable differences in saltation processes across various Gobi surfaces.It is challenging to describe these processes according to a uniform morphology.Therefore,it becomes imperative to articulate surface characteristics through parameters such as the three-dimensional(3D)size and shape of gravel.Collecting morphology information for Gobi gravels is essential for studying its genesis and sand saltation.To enhance the efficiency and information yield of gravel parameter measurements,this study conducted field experiments in the Gobi region across Dunhuang City,Guazhou County,and Yumen City(administrated by Jiuquan City),Gansu Province,China in March 2023.A research framework and methodology for measuring 3D parameters of gravel using point cloud were developed,alongside improved calculation formulas for 3D parameters including gravel grain size,volume,flatness,roundness,sphericity,and equivalent grain size.Leveraging multi-view geometry technology for 3D reconstruction allowed for establishing an optimal data acquisition scheme characterized by high point cloud reconstruction efficiency and clear quality.Additionally,the proposed methodology incorporated point cloud clustering,segmentation,and filtering techniques to isolate individual gravel point clouds.Advanced point cloud algorithms,including the Oriented Bounding Box(OBB),point cloud slicing method,and point cloud triangulation,were then deployed to calculate the 3D parameters of individual gravels.These systematic processes allow precise and detailed characterization of individual gravels.For gravel grain size and volume,the correlation coefficients between point cloud and manual measurements all exceeded 0.9000,confirming the feasibility of the proposed methodology for measuring 3D parameters of individual gravels.The proposed workflow yields accurate calculations of relevant parameters for Gobi gravels,providing essential data support for subsequent studies on Gobi environments.
文摘Understanding the pore water pressure distribution in unsaturated soil is crucial in predicting shallow landslides triggered by rainfall,mainly when dealing with different temporal patterns of rainfall intensity.However,the hydrological response of vegetated slopes,especially three-dimensional(3D)slopes covered with shrubs,under different rainfall patterns remains unclear and requires further investigation.To address this issue,this study adopts a novel 3D numerical model for simulating hydraulic interactions between the root system of the shrub and the surrounding soil.Three series of numerical parametric studies are conducted to investigate the influences of slope inclination,rainfall pattern and rainfall duration.Four rainfall patterns(advanced,bimodal,delayed,and uniform)and two rainfall durations(4-h intense and 168-h mild rainfall)are considered to study the hydrological response of the slope.The computed results show that 17%higher transpiration-induced suction is found for a steeper slope,which remains even after a short,intense rainfall with a 100-year return period.The extreme rainfalls with advanced(PA),bimodal(PB)and uniform(PU)rainfall patterns need to be considered for the short rainfall duration(4 h),while the delayed(PD)and uniform(PU)rainfall patterns are highly recommended for long rainfall durations(168 h).The presence of plants can improve slope stability markedly under extreme rainfall with a short duration(4 h).For the long duration(168 h),the benefit of the plant in preserving pore-water pressure(PWP)and slope stability may not be sufficient.
基金supported by the National Natural Science Foundation of China (No.52374124)National Youth Science Foundation of China (No.52204135)+3 种基金Xing Liao Talent Plan (No.XLYC2202004)Young Elite Scientists Sponsorship Program by CAST (No.2023QNRC001)Liaoning Province International Science and Technology Cooperation Plan (No.2022JH2/1070004)Liaoning Natural Science Foundation Program (No.2022-BS-327)。
文摘The 2D limit equilibrium method is widely used for slope stability analysis.However,with the advancement of dump engineering,composite slopes often exhibit significant 3D mechanical effects.Consequently,it is of significant importance to develop an effective 3D stability calculation method for composite slopes to enhance the design and stability control of open-pit slope engineering.Using the composite slope formed by the mining stope and inner dump in Baiyinhua No.1 and No.2 open-pit coal mine as a case study,this research investigates the failure mode of composite slopes and establishes spatial shape equations for the sliding mass.By integrating the shear resistance and sliding force of each row of microstrip columns onto the bottom surface of the strip corresponding to the main sliding surface,a novel 2D equivalent physical and mechanical parameters analysis method for the strips on the main sliding surface of 3D sliding masses is proposed.Subsequently,a comprehensive 3D stability calculation method for composite slopes is developed,and the quantitative relationship between the coordinated development distance and its 3D stability coefficients is examined.The analysis reveals that the failure mode of the composite slope is characterized by cutting-bedding sliding,with the arc serving as the side interface and the weak layer as the bottom interface,while the destabilization mechanism primarily involves shear failure.The spatial form equation of the sliding mass comprises an ellipsoid and weak plane equation.The analysis revealed that when the coordinated development distance is 1500 m,the error rate between the 3D stability calculation result and the 2D stability calculation result of the composite slope is less than 8%,thereby verifying the proposed analytical method of equivalent physical and mechanical parameters and the 3D stability calculation method for composite slopes.Furthermore,the3D stability coefficient of the composite slope exhibits an exponential correlation with the coordinated development distance,with the coefficient gradually decreasing as the coordinated development distance increases.These findings provide a theoretical guideline for designing similar slope shape parameters and conducting stability analysis.
基金financially supported by the National Natural Science Foundation of China(No.52175352)the Xing Liao Ying Cai Project of Liaoning Province(No.XLYC2008036)the Shenyang Youth Innovation Talent Support Program(No.RC220429)。
文摘The important supporting component in a gas turbine is the casing,which has the characteristics of large size,complex structure,and thin wall.In the context of existing 3DP sand casting processes,casting crack defects are prone to occur.This leads to an increase in the scrap rate of casings,causing significant resource wastage.Additionally,the presence of cracks poses a significant safety hazard after the casings are put into service.The generation of different types of crack defects in stainless steel casings is closely related to casting stress and the high-temperature concession of the sand mold.Therefore,the types and causes of cracks in stainless steel casing products,based on their structural characteristics,were systematically analyzed.Various sand molds with different internal topology designs were printed using the 3DP technology to investigate the impact of sand mold structures on high-temperature concession.The optimal sand mold structure was used to cast casings,and the crack suppression effect was verified by analyzing its eddy current testing results.The experimental results indicate that the skeleton structure has an excellent effect on suppressing cracks in the casing.This research holds important theoretical and engineering significance in improving the quality of casing castings and reducing production costs.
文摘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. .
基金We gratefully acknowledge the financial support of the National Natural Science Foundation of China (Grants 51374213 and 51674251), National Natural Science Fund for Distinguished Young Scholars of China (Grant 51125017), Science Fund for Creative Research Groups of the National Natural Science Foundation of China (Grant 51421003), Fund for Innovative Research and Development Group Program of Jiangsu Province (Grant 2014-27), and the Priority Academic Program Development of Jiangsu Higher Education Institutions (Grant PAPD 2014).
文摘The heterogeneity of unconventional reservoir rock tremendously affects its hydrofracturing behavior. A visual representation and accurate characterization of the three-dimensional (3D) growth and distribution of hydrofracturing cracks within heterogeneous rocks is of particular use to the design and implementation of hydrofracturing stimulation of unconventional reservoirs. However, because of the difficulties involved in visually representing and quantitatively characterizing a 3D hydrofracturing crack-network, this issue remains a challenge. In this paper, a novel method is proposed for physically visualizing and quantitatively characterizing the 3D hydrofracturing crack-network distributed through a heterogeneous structure based on a natural glutenite sample. This method incorporates X-ray microfocus computed tomography (μCT), 3D printing models and hydrofracturing triaxial tests to represent visually the heterogeneous structure, and the 3D crack growth and distribution within a transparent rock model during hydrofracturing. The coupled effects of material heterogeneity and confining geostress on the 3D crack initiation and propagation were analyzed. The results indicate that the breakdown pressure of a heterogeneous rock model is significantly affected by material heterogeneity and confining geostress. The measured breakdown pressures of heterogeneous models are apparently different from those predicted by traditional theories. This study helps to elucidate the quantitative visualization and characterization of the mechanism and influencing factors that determine the hydrofracturing crack initiation and propagation in heterogeneous reservoir rocks.
基金supported by the National Natural Science Foundation of China(Grant No.52125903)the China Postdoctoral Science Foundation(Grant No.2023M730367)the Fundamental Research Funds for Central Public Welfare Research Institutes of China(Grant No.CKSF2023323/YT).
文摘To investigate the long-term stability of deep rocks,a three-dimensional(3D)time-dependent model that accounts for excavation-induced damage and complex stress state is developed.This model comprises three main components:a 3D viscoplastic isotropic constitutive relation that considers excavation damage and complex stress state,a quantitative relationship between critical irreversible deformation and complex stress state,and evolution characteristics of strength parameters.The proposed model is implemented in a self-developed numerical code,i.e.CASRock.The reliability of the model is validated through experiments.It is indicated that the time-dependent fracturing potential index(xTFPI)at a given time during the attenuation creep stage shows a negative correlation with the extent of excavationinduced damage.The time-dependent fracturing process of rock demonstrates a distinct interval effect of the intermediate principal stress,thereby highlighting the 3D stress-dependent characteristic of the model.Finally,the influence of excavation-induced damage and intermediate principal stress on the time-dependent fracturing characteristics of the surrounding rocks around the tunnel is discussed.
文摘A technique for modelling of three-dimensional(3D)quasi-statically propagating cracks in elastic bodies by the displacement discontinuity method(DDM)was described.When the crack is closed,the Mohr-coulomb rule on the two contacted surfaces of the crack must be satisfied.A simple iterative method was adopted in order to consider three different states of cracks.Under the assumption that the advance of the point on the crack front would occur only in the normal plane which is through this edge point,the maximum energy release rate criterion is modified to be used as the criterion for the crack growth.With discretization,the process of crack propagation can be seen as the advance of the vertices of the crack front.The program MCP3D was developed based on these theories to simulate the 3D quasi-static crack propagation.A numerical example of a penny-shaped crack subject to tension and compression in an infinite elastic media was analyzed with MCP3D,and the results in comparison with others' show that the present method for 3D crack propagation is effective.
基金the financial support from the National Natural Science Foundation of China(Grant No.51839003)Liaoning Revitalization Talents Program(Grant No.XLYCYSZX 1902)Hubei Key Laboratory for Efficient Utilization and Agglomeration of Metallurgic Mineral Resources(Grant No.2023zy002).
文摘To achieve the loading of the stress path of hard rock,the spherical discrete element model(DEM)and the new flexible membrane technology were utilized to realize the transient loading of three principal stresses with arbitrary magnitudes and orientations.Furthermore,based on the deep tunnel of China Jinping Underground Laboratory II(CJPL-II),the deformation and fracture evolution characteristics of deep hard rock induced by excavation stress path were analyzed,and the mechanisms of transient loading-unloading and stress rotation-induced fractures were revealed from a mesoscopic perspective.The results indicated that the stressestrain curve exhibits different trends and degrees of sudden changes when subjected to transient changes in principal stress,accompanied by sudden changes in strain rate.Stress rotation induces spatially directional deformation,resulting in fractures of different degrees and orientations,and increasing the degree of deformation anisotropy.The correlation between the degree of induced fracture and the unloading magnitude of minimum principal stress,as well as its initial level is significant and positive.The process of mechanical response during transient unloading exhibits clear nonlinearity and directivity.After transient unloading,both the minimum principal stress and minimum principal strain rate decrease sharply and then tend to stabilize.This occurs from the edge to the interior and from the direction of the minimum principal stress to the direction of the maximum principal stress on theε1-ε3 plane.Transient unloading will induce a tensile stress wave.The ability to induce fractures due to changes in principal stress magnitude,orientation and rotation paths gradually increases.The analysis indicates a positive correlation between the abrupt change amplitude of strain rate and the maximum unloading magnitude,which is determined by the magnitude and rotation of principal stress.A high tensile strain rate is more likely to induce fractures under low minimum principal stress.
基金support from the National Natural Science Foundation of China (Grant No.42207199)Zhejiang Provincial Postdoctoral Science Foundation (Grant Nos.ZJ2022155 and ZJ2022156).
文摘Weak structural plane deformation is responsible for the non-uniform large deformation disasters in layered rock tunnels,resulting in steel arch distortion and secondary lining cracking.In this study,a servo biaxial testing system was employed to conduct physical modeling tests on layered rock tunnels with bedding planes of varying dip angles.The influence of structural anisotropy in layered rocks on the micro displacement and strain field of surrounding rocks was analyzed using digital image correlation(DIC)technology.The spatiotemporal evolution of non-uniform deformation of surrounding rocks was investigated,and numerical simulation was performed to verify the experimental results.The findings indicate that the displacement and strain field of the surrounding layered rocks are all maximized at the horizontal bedding planes and decrease linearly with the increasing dip angle.The failure of the layered surrounding rock with different dip angles occurs and extends along the bedding planes.Compressive strain failure occurs after excavation under high horizontal stress.This study provides significant theoretical support for the analysis,prediction,and control of non-uniform deformation of tunnel surrounding rocks.
基金supported by the National Natural Science Foundation of China (62001436)the Natural Science Foundation of Jiangsu Province under (BK 20190143,JSGG20190823094603691)。
文摘Three-dimensional(3D) synthetic aperture radar(SAR)extends the conventional 2D images into 3D features by several acquisitions in different aspects. Compared with 3D techniques via multiple observations in elevation, e.g. SAR interferometry(InSAR) and SAR tomography(TomoSAR), holographic SAR can retrieve 3D structure by observations in azimuth. This paper focuses on designing a novel type of orbit to achieve SAR regional all-azimuth observation(AAO) for embedded targets detection and holographic 3D reconstruction. The ground tracks of the AAO orbit separate the earth surface into grids. Target in these grids can be accessed with an azimuth angle span of360°, which is similar to the flight path of airborne circular SAR(CSAR). Inspired from the successive coverage orbits of optical sensors, several optimizations are made in the proposed method to ensure favorable grazing angles, the performance of 3D reconstruction, and long-term supervision for SAR sensors. Simulation experiments show the regional AAO can be completed within five hours. In addition, a second AAO of the same area can be duplicated in two days. Finally, an airborne SAR data process result is presented to illustrate the significance of AAO in 3D reconstruction.
基金supported by the National Natural Science Foundation of China(Grant No.51974173)the Natural Science Foundation of Shandong Province,China(Grant No.ZR2020QD122).
文摘The shear failure of intact rock under thermo-mechanical(TM)coupling conditions is common,such as in enhanced geothermal mining and deep mine construction.Under the effect of a continuous engineering disturbance,shear-formed fractures are prone to secondary instability,posing a severe threat to deep engineering.Although numerous studies regarding three-dimensional(3D)morphologies of fracture surfaces have been conducted,the understanding of shear-formed fractures under TM coupling conditions is limited.In this study,direct shear tests of intact granite under various TM coupling conditions were conducted,followed by 3D laser scanning tests of shear-formed fractures.Test results demonstrated that the peak shear strength of intact granite is positively correlated with the normal stress,whereas it is negatively correlated with the temperature.The internal friction angle and cohesion of intact granite significantly decrease with an increase in the temperature.The anisotropy,roughness value,and height of the asperities on the fracture surfaces are reduced as the normal stress increases,whereas their variation trends are the opposite as the temperature increases.The macroscopic failure mode of intact granite under TM coupling conditions is dominated by mixed tensileeshear and shear failures.As the normal stress increases,intragranular fractures are developed ranging from a local to a global distribution,and the macroscopic failure mode of intact granite changes from mixed tensileeshear to shear failure.Finally,3D morphological characteristics of the asperities on the shear-formed fracture surfaces were analyzed,and a quadrangular pyramid conceptual model representing these asperities was proposed and sufficiently verified.
基金Project supported by the Executive Agency for Higher Education Research Development and Innovation Funding of Romania(No.PN-III-P4-PCE-2021-0993)。
文摘This study is concerned with the three-dimensional(3D)stagnation-point for the mixed convection flow past a vertical surface considering the first-order and secondorder velocity slips.To the authors’knowledge,this is the first study presenting this very interesting analysis.Nonlinear partial differential equations for the flow problem are transformed into nonlinear ordinary differential equations(ODEs)by using appropriate similarity transformation.These ODEs with the corresponding boundary conditions are numerically solved by utilizing the bvp4c solver in MATLAB programming language.The effects of the governing parameters on the non-dimensional velocity profiles,temperature profiles,skin friction coefficients,and the local Nusselt number are presented in detail through a series of graphs and tables.Interestingly,it is reported that the reduced skin friction coefficient decreases for the assisting flow situation and increases for the opposing flow situation.The numerical computations of the present work are compared with those from other research available in specific situations,and an excellent consensus is observed.Another exciting feature for this work is the existence of dual solutions.An important remark is that the dual solutions exist for both assisting and opposing flows.A linear stability analysis is performed showing that one solution is stable and the other solution is not stable.We notice that the mixed convection and velocity slip parameters have strong effects on the flow characteristics.These effects are depicted in graphs and discussed in this paper.The obtained results show that the first-order and second-order slip parameters have a considerable effect on the flow,as well as on the heat transfer characteristics.
基金supported by the National Natural Science Foundation of China (Grant Nos. 52104125, U1765204 and 51739008)
文摘Similar to hydraulic fracturing(HF), the coalescence and fracture of cracks are induced within a rock under the action of an ultrasonic field, known as ultrasonic fracturing(UF). Investigating UF is important in both hard rock drilling and oil and gas recovery. A three-dimensional internal laser-engraved crack(3D-ILC) method was introduced to prefabricate two parallel internal cracks within the samples without any damage to the surface. The samples were subjected to UF. The mechanism of UF was elucidated by analyzing the characteristics of fracture surfaces. The crack propagation path under different ultrasonic parameters was obtained by numerical simulation based on the Paris fatigue model and compared to the experimental results of UF. The results show that the 3D-ILC method is a powerful tool for UF research.Under the action of an ultrasonic field, the fracture surface shows the characteristics of beach marks and contains powder locally, indicating that the UF mechanism includes high-cycle fatigue fracture, shear and friction, and temperature load. The two internal cracks become close under UF. The numerical result obtained by the Paris fatigue model also shows the attraction of the two cracks, consistent with the test results. The 3D-ILC method provides a new tool for the experimental study of UF. Compared to the conventional numerical methods based on the analysis of stress-strain and plastic zone, numerical simulation can be a good alternative method to obtain the crack path under UF.
基金supported by the Key Laboratory of Geological Survey and Evaluation of Ministry of Education (China University of Geosciences)(No. GLAB2020ZR13)
文摘The gravity inversion is to restore genetic density distribution of the underground target to be explored for explaining the internal structure and distribution of the Earth.In this paper,we propose a new 3D gravity inversion method based on 3D U-Net++.Compared with two-dimensional gravity inversion,three-dimensional(3D)gravity inversion can more precisely describe the density distribution of underground space.However,conventional 3D gravity inversion method input is two-dimensional,the input and output of the network proposed in our method are three-dimensional.In the training stage,we design a large number of diversifi ed simulation model-data pairs by using the random walk method to improve the generalization ability of the network.In the test phase,we verify the network performance by using the model-data pairs generated by the simulation.To further illustrate the eff ectiveness of the algorithm,we apply the method to the inversion of the San Nicolas mining area,and the inversion results are basically consistent with the borehole measurement results.Moreover,the results of the 3D U-Net++inversion and the 3D U-Net inversion are compared.The density models of the 3D U-Net++inversion have higher resolution,more concentrated inversion results,and a clearer boundary of the density model.
文摘Due to high cost of full-scale experimental setup, this study presents a numerical model on fatigue behaviours of offshore pipeline with multiple coplanar cracks under cyclic tensile loadings. The validation on numerical results is made by other researchers' experimental results, and significant parameters affecting fatigue crack growth are studied.
基金Supported by the National Natural Science Foundation of China (50979052,40872203, 41072234)the Provincial Natural Science Foundation of Shandong (ZR2009FM041,ZR2010EM032,ZR2009AZ001)
文摘The three-dimensional (3D) crack propagation is a hot issue in rock mechanics. To properly simulate 3D crack propagation, a modified maximum tangential tensile stress criterion is proposed. In this modified criterion, it is supposed that cracks propagate only at crack front in the principal normal plane. The tangential tensile stress at crack front in the principal normal plane in local coordinates is employed to determine crack propagation, which is calculated through coordinate transformation from global to local coordinates. New cracks will propagate when the maximum tangential tensile stress at crack front in the principal normal plane reaches the tensile strength of rock-like materials. Compared with the previous crack propagation criteria, the modified crack propagation criterion is helpful in calculating 3D crack stress intensity factor, and can overcome the limitations of propagation step determined by individual experiences in previous studies. Finally, the 3D crack propagation process is traced by element-free Galerkin method. The numerical results agree well with the experimental ones for a frozen resin sample with prefabricated 3D cracks.
基金supported by the National Natural Science Foundation of China(31272518)the program for the New Century Excellent Talents of Ministry of Education of China(NCET-09-0654)+1 种基金the Doctoral Fund of Ministry of Education of P.R.China(RFDP,20120204110030)the Fundamental Research Funds for the Central Universities,China(QN2011012)
文摘Adhesions between different cells and extracellular matrix have been studied extensively in vitro, but little is known about their functions in testicular tissue counterparts. Spermatogonia and their companion somatic cells maintain a close association throughout spermatogenesis and this association is necessary for normal spermatogenesis. In order to keep the relative integrity of the testicular tissues, and to detect the development in vitro, culture testicular tissues in a three- dimensional (3D) agarose matrix was examined. Testicular tissues isolated from 6.5 d postpartum (dpp) mouse were cultured on the top of the matrix for 26 d with a medium height up to 4/5 of the 3D agarose matrix. The results showed that in this 3D culture environment, each type of testicular cells kept the same structure, localization and function as in vivo and might be more biologically relevant to living organisms. After culture, germ cell marker VASA and meiosis markers DAZL and SCP3 showed typical positive analysed by immunofluorescence staining and RT-PCR. It demonstrated that this 3D culture system was able to maintain the number of germ cells and promote the meiosis initiation of male germ cells.
文摘The three-dimensional finite element method of lines is presented, and the basic processing description of 3D FEMOL in cracking questions is given in detail. Applications to 3D bodies with cracks indicate that good accuracy can be obtained with relatively coarse girds. In particular, application to the tension specimen shows very good agreement with the evaluation of stress intensity factors, which is better than the results of other methods. This implies a considerable potential for using this method in the 3D analysis of finite geometry solids and suggests a possible extension of this technique to nonlinear material behavior.
文摘Upon the conservation of mass, momentum and energy, volume fraction and surface penetrative rate were employed to modify the conservative equations to simulate the effect of blockages on fluid flows and heat transfer. These equations were solved numerically with the finite differential method and the primitive variable approach. This method uses staggered grid and pressure correction schemes. A computer code FASTOR3D integrated the aforementioned algorithm. The preliminary results have been compared with conventional benchmark solutions. With auxiliary software DV, the numerical results were visualized in colorful images to demonstrate the variation of flow patterns and temperature profiles during the transient process. The results of the simulation code for the fluid flows and heat transfer in the sodium pool of a fast breeder reactor are acceptable.