This study presents a numerical analysis of three-dimensional steady laminar flow in a rectangular channel with a 180-degree sharp turn. The Navier-Stokes equations are solved by using finite difference method for Re ...This study presents a numerical analysis of three-dimensional steady laminar flow in a rectangular channel with a 180-degree sharp turn. The Navier-Stokes equations are solved by using finite difference method for Re = 900. Three-dimensional streamlines and limiting streamlines on wall surface are used to analyze the three-dimensional flow characteristics. Topological theory is applied to limiting streamlines on inner walls of the channel and two-dimensional streamlines at several cross sections. It is also shown that the flow impinges on the end wall of turn and the secondary flow is induced by the curvature in the sharp turn.展开更多
Fractal theory offers a powerful tool for the precise description and quantification of the complex pore structures in reservoir rocks,crucial for understanding the storage and migration characteristics of media withi...Fractal theory offers a powerful tool for the precise description and quantification of the complex pore structures in reservoir rocks,crucial for understanding the storage and migration characteristics of media within these rocks.Faced with the challenge of calculating the three-dimensional fractal dimensions of rock porosity,this study proposes an innovative computational process that directly calculates the three-dimensional fractal dimensions from a geometric perspective.By employing a composite denoising approach that integrates Fourier transform(FT)and wavelet transform(WT),coupled with multimodal pore extraction techniques such as threshold segmentation,top-hat transformation,and membrane enhancement,we successfully crafted accurate digital rock models.The improved box-counting method was then applied to analyze the voxel data of these digital rocks,accurately calculating the fractal dimensions of the rock pore distribution.Further numerical simulations of permeability experiments were conducted to explore the physical correlations between the rock pore fractal dimensions,porosity,and absolute permeability.The results reveal that rocks with higher fractal dimensions exhibit more complex pore connectivity pathways and a wider,more uneven pore distribution,suggesting that the ideal rock samples should possess lower fractal dimensions and higher effective porosity rates to achieve optimal fluid transmission properties.The methodology and conclusions of this study provide new tools and insights for the quantitative analysis of complex pores in rocks and contribute to the exploration of the fractal transport properties of media within rocks.展开更多
For solving two-dimensional incompressible flow in the vorticity form by the fourth-order compact finite difference scheme and explicit strong stability preserving temporal discretizations,we show that the simple boun...For solving two-dimensional incompressible flow in the vorticity form by the fourth-order compact finite difference scheme and explicit strong stability preserving temporal discretizations,we show that the simple bound-preserving limiter in Li et al.(SIAM J Numer Anal 56:3308–3345,2018)can enforce the strict bounds of the vorticity,if the velocity field satisfies a discrete divergence free constraint.For reducing oscillations,a modified TVB limiter adapted from Cockburn and Shu(SIAM J Numer Anal 31:607–627,1994)is constructed without affecting the bound-preserving property.This bound-preserving finite difference method can be used for any passive convection equation with a divergence free velocity field.展开更多
Energetic Semiconductor bridge(ESCB)based on reactive multilayered films(RMFs)has a promising application in the miniature and intelligence of initiator and pyrotechnics device.Understanding the ignition enhancement m...Energetic Semiconductor bridge(ESCB)based on reactive multilayered films(RMFs)has a promising application in the miniature and intelligence of initiator and pyrotechnics device.Understanding the ignition enhancement mechanism of RMFs on semiconductor bridge(SCB)during the ignition process is crucial for the engineering and practical application of advanced initiator and pyrotechnics devices.In this study,a one-dimensional(1D)gas-solid two-phase flow ignition model was established to study the ignition process of ESCB to charge particles based on the reactivity of Al/MoO_(3) RMFs.In order to fully consider the coupled exothermic between the RMFs and the SCB plasma during the ignition process,the heat release of chemical reaction in RMFs was used as an internal heat source in this model.It is found that the exothermal reaction in RMFs improved the ignition performance of SCB.In the process of plasma rapid condensation with heat release,the product of RMFs enhanced the heat transfer process between the gas phase and the solid charge particle,which accelerated the expansion of hot plasma,and heated the solid charge particle as well as gas phase region with low temperature.In addition,it made up for pressure loss in the gas phase.During the plasma dissipation process,the exothermal chemical reaction in RMFs acted as the main heating source to heat the charge particle,making the surface temperature of the charge particle,gas pressure,and gas temperature rise continuously.This result may yield significant advantages in providing a universal ignition model for miniaturized ignition devices.展开更多
In this paper, we construct a local supersonic flow in a 3-dimensional axis-symmetry nozzle when a uniform supersonic flow inserts the throat. We apply the local existence theory of boundary value problem for quasilin...In this paper, we construct a local supersonic flow in a 3-dimensional axis-symmetry nozzle when a uniform supersonic flow inserts the throat. We apply the local existence theory of boundary value problem for quasilinear hyperbolic system to solve this problem. The boundary value condition is set in particular to guarantee the character number condition. By this trick, the theory in quasilinear hyperbolic system can be employed to a large range of the boundary value problem.展开更多
Rainbow particle image velocimetry(PIV)can restore the three-dimensional velocity field of particles with a single camera;however,it requires a relatively long time to complete the reconstruction.This paper proposes a...Rainbow particle image velocimetry(PIV)can restore the three-dimensional velocity field of particles with a single camera;however,it requires a relatively long time to complete the reconstruction.This paper proposes a hybrid algorithm that combines the fast Fourier transform(FFT)based co-correlation algorithm and the Horn–Schunck(HS)optical flow pyramid iterative algorithm to increase the reconstruction speed.The Rankine vortex simulation experiment was performed,in which the particle velocity field was reconstructed using the proposed algorithm and the rainbow PIV method.The average endpoint error and average angular error of the proposed algorithm were roughly the same as those of the rainbow PIV algorithm;nevertheless,the reconstruction time was 20%shorter.Furthermore,the effect of velocity magnitude and particle density on the reconstruction results was analyzed.In the end,the performance of the proposed algorithm was verified using real experimental single-vortex and double-vortex datasets,from which a similar particle velocity field was obtained compared with the rainbow PIV algorithm.The results show that the reconstruction speed of the proposed hybrid algorithm is approximately 25%faster than that of the rainbow PIV algorithm.展开更多
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 images of granular ore media were captured by X-ray CT scanner. Combined with digital image processing and finite element techniques, the three-dimensional geometrical model, which represents the realistic pore st...The images of granular ore media were captured by X-ray CT scanner. Combined with digital image processing and finite element techniques, the three-dimensional geometrical model, which represents the realistic pore structure of the media, was constructed. With this model, three dimensional pore scale fluid flow among particles was simulated. Then the distributions of fluid flow velocity and pressure were analyzed and the hydraulic conductivity was calculated. The simulation results indicate the fluid flow behaviors are mainly dominated by the volume and topological structure of pore space. There exist obvious preferential flow and leaching blind zones simultaneously in the medium. The highest velocities generally occur in those narrow pores with high pressure drops. The hydraulic conductivity obtained by simulation is the same order of magnitude as the laboratory test result, which denotes the validity of the model. The pore-scale and macro-scale are combined and the established geometrical model can be used for the simulations of other phenomena during heap leaching process.展开更多
The hydrodynamic characteristics of a rigid, single, circular cylinder in a three dimensional, incompressible, uniform cross flow were calculated using the large-eddy simulation method of CFX5. Solutions to the three ...The hydrodynamic characteristics of a rigid, single, circular cylinder in a three dimensional, incompressible, uniform cross flow were calculated using the large-eddy simulation method of CFX5. Solutions to the three dimensional N-S equations were obtained by the finite volume method. The focus of this numerical simulation was to research the characteristics of pressure distribution (drag and litt forces) and vortex tubes at high Reynolds numbers. The results of the calculations showed that the forces at every section in the spanwise direction of the cylinder were symmetrical about the middle section and smaller than the forces calculated in two dimensional cases. Moreover, the flow around the cylinder obviously presents three dimensional characteristics.展开更多
Monte Carlo simulations are frequently utilized in radiation dose assessments.However,many researchers find the prevailing computing platforms to be intricate.This highlights a pressing need for a specialized framewor...Monte Carlo simulations are frequently utilized in radiation dose assessments.However,many researchers find the prevailing computing platforms to be intricate.This highlights a pressing need for a specialized framework for phantom dose evalua-tion.To address this gap,we developed a user-friendly radiation dose assessment platform using the Monte Carlo toolkit,Geant4.The Tsinghua University Phantom Dose(THUDosePD)augments the flexibility of Monte Carlo simulations in dosi-metric research.Originating from THUDose,a code with generic,functional,and application layers,THUDosePD focuses predominantly on anatomical phantom dose assessment.Additionally,it enables medical exposure simulation,intricate geometry creation,and supports both three-dimensional radiation dose analysis and phantom format transformations.The system operates on a multi-threaded parallel CPU architecture,with some modules enhanced for GPU parallel computing.Benchmark tests on the ICRP reference male illustrated the capabilities of THUDosePD in phantom dose assessment,covering the effective dose,three-dimensional dose distribution,and three-dimensional organ dose.We also conducted a voxelization conversion on the polygon mesh phantom,demonstrating the method’s efficiency and consistency.Extended applications based on THUDosePD further underline its broad adaptability.This intuitive,three-dimensional platform stands out as a valuable tool for phantom radiation dosimetry research.展开更多
Load transformation from the yielding part of the soil to the adjacent part is known as the soil arching effect,which plays an important role in the design of various geotechnical infrastructures.Terzaghi’s trapdoor ...Load transformation from the yielding part of the soil to the adjacent part is known as the soil arching effect,which plays an important role in the design of various geotechnical infrastructures.Terzaghi’s trapdoor test was an importantmilestone in the development of theories on soil arching.The research on earth pressure of the trapdoor problem is presented in this paper using the three-dimensional(3D)discrete element method(DEM).Five 3D trapdoor models with different heights are established by 3DDEMsoftware PFC 3D.The variation of earth pressure on the trapdoor with the downward movement of the trapdoor,the distribution of vertical earth pressure along the horizontal direction,the distribution of vertical earth pressure along the vertical direction,the distribution of lateral earth pressure coefficient along the depth direction,the magnitude and direction of contact force chain are studied,respectively.Related research results show that the earth pressure on the trapdoor decreases rapidly after the downward movement of the trapdoor,and then reaches the minimum earth pressure.After that,the earth’s pressure will rise slightly,and whether this phenomenon occurs depends on the depth ratio.For the bottom soil,due to the stress transfer caused by the soil arching effect,the ratio of earth pressure in the loose area decreases,while the ratio of earth pressure in the stable area increases.With the trapdoor moving down,the vertical earth pressure along the depth in the stable zone is basically consistent with the initial state,which shows an approximate linear distribution.After the trapdoor moves down,the distribution of earth pressure along with the depth in the loose area changes,which is far less than the theoretical value of vertical earth pressure of its self-weight.Because of the compression of the soil on both sides,the lateral earth pressure coefficient of most areas on the central axis of the loose zone is close to the passive earth pressure coefficient Kp.The existence of a‘soil arch’can be observed intuitively from the distribution diagram of the contact force chain in the loose zone.展开更多
Gaseous detonation propagating in a toroidal chamber was numerically studied for hydrogen/oxygen/nitrogen mixtures. The numerical method used is based on the three-dimensional Euler equations with detailed finiterate ...Gaseous detonation propagating in a toroidal chamber was numerically studied for hydrogen/oxygen/nitrogen mixtures. The numerical method used is based on the three-dimensional Euler equations with detailed finiterate chemistry. The results show that the calculated streak picture is in qualitative agreement with the picture recorded by a high speed streak camera from published literature. The three-dimensional flow field induced by a continuously rotating detonation was visualized and distinctive features of the rotating detonations were clearly depicted. Owing to the unconfined character of detonation wavelet, a deficit of detonation parameters was observed. Due to the effects of wall geometries, the strength of the outside detonation front is stronger than that of the inside portion. The detonation thus propagates with a constant circular velocity. Numerical simulation also shows three-dimensional rotating detonation structures, which display specific feature of the detonation- shock combined wave. Discrete burning gas pockets are formed due to instability of the discontinuity. It is believed that the present study could give an insight into the interest- ing properties of the continuously rotating detonation, and is thus beneficial to the design of continuous detonation propulsion systems.展开更多
In this communication a generalized three- dimensional steady flow of a viscous fluid between two infinite parallel plates is considered. The flow is generated due to uniform stretching of the lower plate in x- and y-...In this communication a generalized three- dimensional steady flow of a viscous fluid between two infinite parallel plates is considered. The flow is generated due to uniform stretching of the lower plate in x- and y-directions. It is assumed that the upper plate is uniformly porous and is subjected to constant injection. The governing system is fully coupled and nonlinear in nature. A complete analytic solution which is uniformly valid for all values of the dimensionless parameters β Re and λ is obtained by using a purely analytic technique, namely the homotopy analysis method. Also the effects of the parameters β Re and λ on the velocity field are discussed through graphs.展开更多
文摘This study presents a numerical analysis of three-dimensional steady laminar flow in a rectangular channel with a 180-degree sharp turn. The Navier-Stokes equations are solved by using finite difference method for Re = 900. Three-dimensional streamlines and limiting streamlines on wall surface are used to analyze the three-dimensional flow characteristics. Topological theory is applied to limiting streamlines on inner walls of the channel and two-dimensional streamlines at several cross sections. It is also shown that the flow impinges on the end wall of turn and the secondary flow is induced by the curvature in the sharp turn.
基金supported by the National Natural Science Foundation of China (Nos.52374078 and 52074043)the Fundamental Research Funds for the Central Universities (No.2023CDJKYJH021)。
文摘Fractal theory offers a powerful tool for the precise description and quantification of the complex pore structures in reservoir rocks,crucial for understanding the storage and migration characteristics of media within these rocks.Faced with the challenge of calculating the three-dimensional fractal dimensions of rock porosity,this study proposes an innovative computational process that directly calculates the three-dimensional fractal dimensions from a geometric perspective.By employing a composite denoising approach that integrates Fourier transform(FT)and wavelet transform(WT),coupled with multimodal pore extraction techniques such as threshold segmentation,top-hat transformation,and membrane enhancement,we successfully crafted accurate digital rock models.The improved box-counting method was then applied to analyze the voxel data of these digital rocks,accurately calculating the fractal dimensions of the rock pore distribution.Further numerical simulations of permeability experiments were conducted to explore the physical correlations between the rock pore fractal dimensions,porosity,and absolute permeability.The results reveal that rocks with higher fractal dimensions exhibit more complex pore connectivity pathways and a wider,more uneven pore distribution,suggesting that the ideal rock samples should possess lower fractal dimensions and higher effective porosity rates to achieve optimal fluid transmission properties.The methodology and conclusions of this study provide new tools and insights for the quantitative analysis of complex pores in rocks and contribute to the exploration of the fractal transport properties of media within rocks.
文摘For solving two-dimensional incompressible flow in the vorticity form by the fourth-order compact finite difference scheme and explicit strong stability preserving temporal discretizations,we show that the simple bound-preserving limiter in Li et al.(SIAM J Numer Anal 56:3308–3345,2018)can enforce the strict bounds of the vorticity,if the velocity field satisfies a discrete divergence free constraint.For reducing oscillations,a modified TVB limiter adapted from Cockburn and Shu(SIAM J Numer Anal 31:607–627,1994)is constructed without affecting the bound-preserving property.This bound-preserving finite difference method can be used for any passive convection equation with a divergence free velocity field.
基金supported by the National Natural Science Foundation of China(Grant Nos.22275092,52102107 and 52372084)the Fundamental Research Funds for the Central Universities(Grant No.30923010920)。
文摘Energetic Semiconductor bridge(ESCB)based on reactive multilayered films(RMFs)has a promising application in the miniature and intelligence of initiator and pyrotechnics device.Understanding the ignition enhancement mechanism of RMFs on semiconductor bridge(SCB)during the ignition process is crucial for the engineering and practical application of advanced initiator and pyrotechnics devices.In this study,a one-dimensional(1D)gas-solid two-phase flow ignition model was established to study the ignition process of ESCB to charge particles based on the reactivity of Al/MoO_(3) RMFs.In order to fully consider the coupled exothermic between the RMFs and the SCB plasma during the ignition process,the heat release of chemical reaction in RMFs was used as an internal heat source in this model.It is found that the exothermal reaction in RMFs improved the ignition performance of SCB.In the process of plasma rapid condensation with heat release,the product of RMFs enhanced the heat transfer process between the gas phase and the solid charge particle,which accelerated the expansion of hot plasma,and heated the solid charge particle as well as gas phase region with low temperature.In addition,it made up for pressure loss in the gas phase.During the plasma dissipation process,the exothermal chemical reaction in RMFs acted as the main heating source to heat the charge particle,making the surface temperature of the charge particle,gas pressure,and gas temperature rise continuously.This result may yield significant advantages in providing a universal ignition model for miniaturized ignition devices.
文摘In this paper, we construct a local supersonic flow in a 3-dimensional axis-symmetry nozzle when a uniform supersonic flow inserts the throat. We apply the local existence theory of boundary value problem for quasilinear hyperbolic system to solve this problem. The boundary value condition is set in particular to guarantee the character number condition. By this trick, the theory in quasilinear hyperbolic system can be employed to a large range of the boundary value problem.
基金the National Natural Science Foundation of China(Grant Nos.51874264 and 52076200)。
文摘Rainbow particle image velocimetry(PIV)can restore the three-dimensional velocity field of particles with a single camera;however,it requires a relatively long time to complete the reconstruction.This paper proposes a hybrid algorithm that combines the fast Fourier transform(FFT)based co-correlation algorithm and the Horn–Schunck(HS)optical flow pyramid iterative algorithm to increase the reconstruction speed.The Rankine vortex simulation experiment was performed,in which the particle velocity field was reconstructed using the proposed algorithm and the rainbow PIV method.The average endpoint error and average angular error of the proposed algorithm were roughly the same as those of the rainbow PIV algorithm;nevertheless,the reconstruction time was 20%shorter.Furthermore,the effect of velocity magnitude and particle density on the reconstruction results was analyzed.In the end,the performance of the proposed algorithm was verified using real experimental single-vortex and double-vortex datasets,from which a similar particle velocity field was obtained compared with the rainbow PIV algorithm.The results show that the reconstruction speed of the proposed hybrid algorithm is approximately 25%faster than that of the rainbow PIV algorithm.
文摘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. .
基金Projects (50934002, 51074013, 51104100) supported by the National Natural Science Foundation of China
文摘The images of granular ore media were captured by X-ray CT scanner. Combined with digital image processing and finite element techniques, the three-dimensional geometrical model, which represents the realistic pore structure of the media, was constructed. With this model, three dimensional pore scale fluid flow among particles was simulated. Then the distributions of fluid flow velocity and pressure were analyzed and the hydraulic conductivity was calculated. The simulation results indicate the fluid flow behaviors are mainly dominated by the volume and topological structure of pore space. There exist obvious preferential flow and leaching blind zones simultaneously in the medium. The highest velocities generally occur in those narrow pores with high pressure drops. The hydraulic conductivity obtained by simulation is the same order of magnitude as the laboratory test result, which denotes the validity of the model. The pore-scale and macro-scale are combined and the established geometrical model can be used for the simulations of other phenomena during heap leaching process.
文摘The hydrodynamic characteristics of a rigid, single, circular cylinder in a three dimensional, incompressible, uniform cross flow were calculated using the large-eddy simulation method of CFX5. Solutions to the three dimensional N-S equations were obtained by the finite volume method. The focus of this numerical simulation was to research the characteristics of pressure distribution (drag and litt forces) and vortex tubes at high Reynolds numbers. The results of the calculations showed that the forces at every section in the spanwise direction of the cylinder were symmetrical about the middle section and smaller than the forces calculated in two dimensional cases. Moreover, the flow around the cylinder obviously presents three dimensional characteristics.
基金This work was supported by the National Natural Science Foundation of China(General Program)(Nos.12175114,U2167209)the Foundation of Key Laboratory of Metrology and Calibration Technology(No.JLKG2022001C001)+2 种基金the Platform Development foundation of China Institute for Radiation Protection(No.YP21030101)the National Key R&D Program of China(No.2021YFF0603600)the Tsinghua University Initiative Scientific Research Program(No.20211080081).
文摘Monte Carlo simulations are frequently utilized in radiation dose assessments.However,many researchers find the prevailing computing platforms to be intricate.This highlights a pressing need for a specialized framework for phantom dose evalua-tion.To address this gap,we developed a user-friendly radiation dose assessment platform using the Monte Carlo toolkit,Geant4.The Tsinghua University Phantom Dose(THUDosePD)augments the flexibility of Monte Carlo simulations in dosi-metric research.Originating from THUDose,a code with generic,functional,and application layers,THUDosePD focuses predominantly on anatomical phantom dose assessment.Additionally,it enables medical exposure simulation,intricate geometry creation,and supports both three-dimensional radiation dose analysis and phantom format transformations.The system operates on a multi-threaded parallel CPU architecture,with some modules enhanced for GPU parallel computing.Benchmark tests on the ICRP reference male illustrated the capabilities of THUDosePD in phantom dose assessment,covering the effective dose,three-dimensional dose distribution,and three-dimensional organ dose.We also conducted a voxelization conversion on the polygon mesh phantom,demonstrating the method’s efficiency and consistency.Extended applications based on THUDosePD further underline its broad adaptability.This intuitive,three-dimensional platform stands out as a valuable tool for phantom radiation dosimetry research.
基金supports from National Natural Science Foundation of China (NSFC Grant No.52008373)Natural Science Foundation of Zhejiang Province of China (No.Q22E080445)are greatly acknowledged.
文摘Load transformation from the yielding part of the soil to the adjacent part is known as the soil arching effect,which plays an important role in the design of various geotechnical infrastructures.Terzaghi’s trapdoor test was an importantmilestone in the development of theories on soil arching.The research on earth pressure of the trapdoor problem is presented in this paper using the three-dimensional(3D)discrete element method(DEM).Five 3D trapdoor models with different heights are established by 3DDEMsoftware PFC 3D.The variation of earth pressure on the trapdoor with the downward movement of the trapdoor,the distribution of vertical earth pressure along the horizontal direction,the distribution of vertical earth pressure along the vertical direction,the distribution of lateral earth pressure coefficient along the depth direction,the magnitude and direction of contact force chain are studied,respectively.Related research results show that the earth pressure on the trapdoor decreases rapidly after the downward movement of the trapdoor,and then reaches the minimum earth pressure.After that,the earth’s pressure will rise slightly,and whether this phenomenon occurs depends on the depth ratio.For the bottom soil,due to the stress transfer caused by the soil arching effect,the ratio of earth pressure in the loose area decreases,while the ratio of earth pressure in the stable area increases.With the trapdoor moving down,the vertical earth pressure along the depth in the stable zone is basically consistent with the initial state,which shows an approximate linear distribution.After the trapdoor moves down,the distribution of earth pressure along with the depth in the loose area changes,which is far less than the theoretical value of vertical earth pressure of its self-weight.Because of the compression of the soil on both sides,the lateral earth pressure coefficient of most areas on the central axis of the loose zone is close to the passive earth pressure coefficient Kp.The existence of a‘soil arch’can be observed intuitively from the distribution diagram of the contact force chain in the loose zone.
基金supported by the National Natural Science Foundation of China (10872096)the Open Fund of State Key Laboratory of Explosion Science and Technology, Beijing University of Science and Technology (KFJJ09-13)
文摘Gaseous detonation propagating in a toroidal chamber was numerically studied for hydrogen/oxygen/nitrogen mixtures. The numerical method used is based on the three-dimensional Euler equations with detailed finiterate chemistry. The results show that the calculated streak picture is in qualitative agreement with the picture recorded by a high speed streak camera from published literature. The three-dimensional flow field induced by a continuously rotating detonation was visualized and distinctive features of the rotating detonations were clearly depicted. Owing to the unconfined character of detonation wavelet, a deficit of detonation parameters was observed. Due to the effects of wall geometries, the strength of the outside detonation front is stronger than that of the inside portion. The detonation thus propagates with a constant circular velocity. Numerical simulation also shows three-dimensional rotating detonation structures, which display specific feature of the detonation- shock combined wave. Discrete burning gas pockets are formed due to instability of the discontinuity. It is believed that the present study could give an insight into the interest- ing properties of the continuously rotating detonation, and is thus beneficial to the design of continuous detonation propulsion systems.
文摘In this communication a generalized three- dimensional steady flow of a viscous fluid between two infinite parallel plates is considered. The flow is generated due to uniform stretching of the lower plate in x- and y-directions. It is assumed that the upper plate is uniformly porous and is subjected to constant injection. The governing system is fully coupled and nonlinear in nature. A complete analytic solution which is uniformly valid for all values of the dimensionless parameters β Re and λ is obtained by using a purely analytic technique, namely the homotopy analysis method. Also the effects of the parameters β Re and λ on the velocity field are discussed through graphs.