Scouring experiments were conducted using a three-dimensional laser scanning technology for angles of the jet spanning the interval from 0°to 30°,and the characteristics of the scour hole in equilibrium cond...Scouring experiments were conducted using a three-dimensional laser scanning technology for angles of the jet spanning the interval from 0°to 30°,and the characteristics of the scour hole in equilibrium conditions were investigated accordingly.The results indicate that the optimal scouring effects occur when the jet angle is in the ranges between 15°and 20°.Moreover,the dimensionless profiles of the scour hole exhibit a high degree of similarity at different jet angles.Numerical simulations conducted using the Flow-3D software to investigate the bed shear stress along the jet impingement surface have shown that this stress is influenced by both the resultant force and the jet impingement surface area.It reaches its maximum value when the jet is vertical,decreases rapidly as the jet starts to tilt,then increases slightly,and decreases again significantly when the angle exceeds 20°.展开更多
The geometry of joints has a significant influence on the mechanical properties of rocks.To simplify the curved joint shapes in rocks,the joint shape is usually treated as straight lines or planes in most laboratory e...The geometry of joints has a significant influence on the mechanical properties of rocks.To simplify the curved joint shapes in rocks,the joint shape is usually treated as straight lines or planes in most laboratory experiments and numerical simulations.In this study,the computerized tomography (CT) scanning and photogrammetry were employed to obtain the internal and surface joint structures of a limestone sample,respectively.To describe the joint geometry,the edge detection algorithms and a three-dimensional (3D) matrix mapping method were applied to reconstruct CT-based and photogrammetry-based jointed rock models.For comparison tests,the numerical uniaxial compression tests were conducted on an intact rock sample and a sample with a joint simplified to a plane using the parallel computing method.The results indicate that the mechanical characteristics and failure process of jointed rocks are significantly affected by the geometry of joints.The presence of joints reduces the uniaxial compressive strength (UCS),elastic modulus,and released acoustic emission (AE) energy of rocks by 37%–67%,21%–24%,and 52%–90%,respectively.Compared to the simplified joint sample,the proposed photogrammetry-based numerical model makes the most of the limited geometry information of joints.The UCS,accumulative released AE energy,and elastic modulus of the photogrammetry-based sample were found to be very close to those of the CT-based sample.The UCS value of the simplified joint sample (i.e.38.5 MPa) is much lower than that of the CT-based sample (i.e.72.3 MPa).Additionally,the accumulative released AE energy observed in the simplified joint sample is 3.899 times lower than that observed in the CT-based sample.CT scanning provides a reliable means to visualize the joints in rocks,which can be used to verify the reliability of photogrammetry techniques.The application of the photogrammetry-based sample enables detailed analysis for estimating the mechanical properties of jointed rocks.展开更多
Aiming at the problems that the simulation accuracy which is reduced due to the simplification of the model,a three-dimensional simulation method based on solid modeling is being proposed.By analyzing the motion relat...Aiming at the problems that the simulation accuracy which is reduced due to the simplification of the model,a three-dimensional simulation method based on solid modeling is being proposed.By analyzing the motion relationship and positional relationship between the caries knife and the workpiece,the coordinate system of the caries machining was established.With the MATLAB software,the cutting edge model and the blade sweeping surface model of the boring cutter are sequentially established.Boolean operation is performed on the blade swept surface formed by the tooth cutter teeth with time t and the workpiece tooth geometry as well as the undeformed three-dimensional chip geometry model and the instantaneous cogging geometry model are obtained at different times.Through the compare between gear end face simulation tooth profile and the theoretical inner arc tooth profile,we verified the accuracy and rationality of the proposed method.展开更多
A fracture propagation model of radial well fracturing is established based on the finite element-meshless method.The model considers the coupling effect of fracturing fluid flow and rock matrix deformation.The fractu...A fracture propagation model of radial well fracturing is established based on the finite element-meshless method.The model considers the coupling effect of fracturing fluid flow and rock matrix deformation.The fracture geometries of radial well fracturing are simulated,the induction effect of radial well on the fracture is quantitatively characterized,and the influences of azimuth,horizontal principle stress difference,and reservoir matrix permeability on the fracture geometries are revealed.The radial wells can induce the fractures to extend parallel to their axes when two radial wells in the same layer are fractured.When the radial wells are symmetrically distributed along the direction of the minimum horizontal principle stress with the azimuth greater than 15,the extrusion effect reduces the fracture length of radial wells.When the radial wells are symmetrically distributed along the direction of the maximum horizontal principal stress,the extrusion increases the fracture length of the radial wells.The fracture geometries are controlled by the rectification of radial borehole,the extrusion between radial wells in the same layer,and the deflection of the maximum horizontal principal stress.When the radial wells are distributed along the minimum horizontal principal stress symmetrically,the fracture length induced by the radial well decreases with the increase of azimuth;in contrast,when the radial wells are distributed along the maximum horizontal principal stress symmetrically,the fracture length induced by the radial well first decreases and then increases with the increase of azimuth.The fracture length induced by the radial well decreases with the increase of horizontal principal stress difference.The increase of rock matrix permeability and pore pressure of the matrix around radial wells makes the inducing effect of the radial well on fractures increase.展开更多
Organic solid waste(OSW)contains many renewable materials.The pyrolysis and gasification of OSW can realize resource utilization,and its products can be used for methanation reaction to produce synthetic natural gas i...Organic solid waste(OSW)contains many renewable materials.The pyrolysis and gasification of OSW can realize resource utilization,and its products can be used for methanation reaction to produce synthetic natural gas in the specific reactor.In order to understand the dynamic characteristics of the reactor,a three-dimensional numerical model has been established by the method of Computational Fluid Dynamics(CFD).Along the height of the reactor,the particle distribution in the bed becomes thinner and the mean solid volume fraction decreases from 4.18%to 0.37%.Meanwhile,the pressure fluctuation range decreased from 398.76 Pa at the entrance to a much lower value of 74.47 Pa at the exit.In this simulation,three parameters of gas inlet velocity,operating temperature and solid particle diameter are changed to explore their influences on gas-solid multiphase flow.The results show that gas velocity has a great influence on particle distribution.When the gas inlet velocity decreases from 6.51 to 1.98 m/s,the minimum height that particles can reach decreases from 169 to 100 mm.Additionally,as the operating temperature increases,the particle holdup inside the reactor changes from 0.843%to 0.700%.This indicates that the particle residence time reduces,which is not conducive to the follow-up reaction.Moreover,with the increase of particle size,the fluctuation range of the pressure at the bottom of the reactor increases,and its standard deviation increases from 55.34 to 1266.37 Pa.展开更多
Numerical simulation is the most powerful computational and analysis tool for a large variety of engineering and physical problems.For a complex problem relating to multi-field,multi-process and multi-scale,different ...Numerical simulation is the most powerful computational and analysis tool for a large variety of engineering and physical problems.For a complex problem relating to multi-field,multi-process and multi-scale,different computing tools have to be developed so as to solve particular fields at different scales and for different processes.Therefore,the integration of different types of software is inevitable.However,it is difficult to perform the transfer of the meshes and simulated results among software packages because of the lack of shared data formats or encrypted data formats.An image processing based method for three-dimensional model reconstruction for numerical simulation was proposed,which presents a solution to the integration problem by a series of slice or projection images obtained by the post-processing modules of the numerical simulation software.By means of mapping image pixels to meshes of either finite difference or finite element models,the geometry contour can be extracted to export the stereolithography model.The values of results,represented by color,can be deduced and assigned to the meshes.All the models with data can be directly or indirectly integrated into other software as a continued or new numerical simulation.The three-dimensional reconstruction method has been validated in numerical simulation of castings and case studies were provided in this study.展开更多
Three-dimensional direct numerical simulations of the wake flow downstream of a near-wall circular cylinder at different gap ratios and boundary layer thicknesses are carried out by using the iterative immersed bounda...Three-dimensional direct numerical simulations of the wake flow downstream of a near-wall circular cylinder at different gap ratios and boundary layer thicknesses are carried out by using the iterative immersed boundary method.The non-dimensional gap between the cylinder and the wall,G/D=0.2,0.6 and 1.0,the non-dimensional boundary layer thickness,δ/D=0.0,0.7 and 1.6,the Reynolds number,Re=350,and the aspect ratio of the cylinder,L/D=25are adopted.High-resolution visualizations of the complex vortex structures at differentδ/D and G/D are presented.The transition of the streamwise vortex mode,the combined effects ofδ/D and G/D on the flow statistics,the pressure and shear stress distribution and the hydrodynamic forces are analyzed.Results show that with decreasing G/D and increasingδ/D,the gap flow and its vortex-shedding are significantly weakened,together with an elongated wake and an enlarged low-velocity area near the wall,leading to the wake mode transition from the two-sided to one-sided vortex-shedding.Different relative positions of the cylinder regarding the boundary layer alter the flow features of the shear layers.With an increase inδ/D,the front stagnation point shifts to the upper surface,and the distance between the flow divergence point and the maximum pressure position increases.The mean drag coefficient and r.m.s.values of drag and lift coefficients decrease with a decrease in G/D and an increase inδ/D,while the mean lift coefficient increases with decreasing G/D but decreases with increasingδ/D.展开更多
Hypoxia is a typical feature of the tumor microenvironment,one of the most critical factors affecting cell behavior and tumor progression.However,the lack of tumor models able to precisely emulate natural brain tumor ...Hypoxia is a typical feature of the tumor microenvironment,one of the most critical factors affecting cell behavior and tumor progression.However,the lack of tumor models able to precisely emulate natural brain tumor tissue has impeded the study of the effects of hypoxia on the progression and growth of tumor cells.This study reports a three-dimensional(3D)brain tumor model obtained by encapsulating U87MG(U87)cells in a hydrogel containing type I collagen.It also documents the effect of various oxygen concentrations(1%,7%,and 21%)in the culture environment on U87 cell morphology,proliferation,viability,cell cycle,apoptosis rate,and migration.Finally,it compares two-dimensional(2D)and 3D cultures.For comparison purposes,cells cultured in flat culture dishes were used as the control(2D model).Cells cultured in the 3D model proliferated more slowly but had a higher apoptosis rate and proportion of cells in the resting phase(G0 phase)/gap I phase(G1 phase)than those cultured in the 2D model.Besides,the two models yielded significantly different cell morphologies.Finally,hypoxia(e.g.,1%O2)affected cell morphology,slowed cell growth,reduced cell viability,and increased the apoptosis rate in the 3D model.These results indicate that the constructed 3D model is effective for investigating the effects of biological and chemical factors on cell morphology and function,and can be more representative of the tumor microenvironment than 2D culture systems.The developed 3D glioblastoma tumor model is equally applicable to other studies in pharmacology and pathology.展开更多
The process of entrainment-mixing between cumulus clouds and the ambient air is important for the development of cumulus clouds.Accurately obtaining the entrainment rate(λ)is particularly important for its parameteri...The process of entrainment-mixing between cumulus clouds and the ambient air is important for the development of cumulus clouds.Accurately obtaining the entrainment rate(λ)is particularly important for its parameterization within the overall cumulus parameterization scheme.In this study,an improved bulk-plume method is proposed by solving the equations of two conserved variables simultaneously to calculateλof cumulus clouds in a large-eddy simulation.The results demonstrate that the improved bulk-plume method is more reliable than the traditional bulk-plume method,becauseλ,as calculated from the improved method,falls within the range ofλvalues obtained from the traditional method using different conserved variables.The probability density functions ofλfor all data,different times,and different heights can be well-fitted by a log-normal distribution,which supports the assumed stochastic entrainment process in previous studies.Further analysis demonstrate that the relationship betweenλand the vertical velocity is better than other thermodynamic/dynamical properties;thus,the vertical velocity is recommended as the primary influencing factor for the parameterization ofλin the future.The results of this study enhance the theoretical understanding ofλand its influencing factors and shed new light on the development ofλparameterization.展开更多
Earth’s magnetopause is a thin boundary separating the shocked solar wind plasma from the magnetospheric plasmas,and it is also the boundary of the solar wind energy transport to the magnetosphere.Soft X-ray imaging ...Earth’s magnetopause is a thin boundary separating the shocked solar wind plasma from the magnetospheric plasmas,and it is also the boundary of the solar wind energy transport to the magnetosphere.Soft X-ray imaging allows investigation of the large-scale magnetopause by providing a two-dimensional(2-D)global view from a satellite.By performing 3-D global hybrid-particle-in-cell(hybrid-PIC)simulations,we obtain soft X-ray images of Earth’s magnetopause under different solar wind conditions,such as different plasma densities and directions of the southward interplanetary magnetic field.In all cases,magnetic reconnection occurs at low latitude magnetopause.The soft X-ray images observed by a hypothetical satellite are shown,with all of the following identified:the boundary of the magnetopause,the cusps,and the magnetosheath.Local X-ray emissivity in the magnetosheath is characterized by large amplitude fluctuations(up to 160%);however,the maximum line-of-sight-integrated X-ray intensity matches the tangent directions of the magnetopause well,indicating that these fluctuations have limited impact on identifying the magnetopause boundary in the X-ray images.Moreover,the magnetopause boundary can be identified using multiple viewing geometries.We also find that solar wind conditions have little effect on the magnetopause identification.The Solar wind Magnetosphere Ionosphere Link Explorer(SMILE)mission will provide X-ray images of the magnetopause for the first time,and our global hybrid-PIC simulation results can help better understand the 2-D X-ray images of the magnetopause from a 3-D perspective,with particle kinetic effects considered.展开更多
We derived the properties of the terrestrial magnetopause(MP)from two modeling approaches,one global–fluid,the other local–kinetic,and compared the results with data collected in situ by the Magnetospheric Multiscal...We derived the properties of the terrestrial magnetopause(MP)from two modeling approaches,one global–fluid,the other local–kinetic,and compared the results with data collected in situ by the Magnetospheric Multiscale 2(MMS2)spacecraft.We used global magnetohydrodynamic(MHD)simulations of the Earth’s magnetosphere(publicly available from the NASA-CCMC[National Aeronautics and Space Administration–Community Coordinated Modeling Center])and local Vlasov equilibrium models(based on kinetic models for tangential discontinuities)to extract spatial profiles of the plasma and field variables at the Earth’s MP.The global MHD simulations used initial solar wind conditions extracted from the OMNI database at the time epoch when the MMS2 observes the MP.The kinetic Vlasov model used asymptotic boundary conditions derived from the same in situ MMS measurements upstream or downstream of the MP.The global MHD simulations provide a three-dimensional image of the magnetosphere at the time when the MMS2 crosses the MP.The Vlasov model provides a one-dimensional local view of the MP derived from first principles of kinetic theory.The MMS2 experimental data also serve as a reference for comparing and validating the numerical simulations and modeling.We found that the MP transition layer formed in global MHD simulations was generally localized closer to the Earth(roughly by one Earth radius)from the position of the real MP observed by the MMS.We also found that the global MHD simulations overestimated the thickness of the MP transition by one order of magnitude for three analyzed variables:magnetic field,density,and tangential speed.The MP thickness derived from the local Vlasov equilibrium was consistent with observations for all three of these variables.The overestimation of density in the Vlasov equilibrium was reduced compared with the global MHD solutions.We discuss our results in the context of future SMILE(Solar wind Magnetosphere Ionosphere Link Explorer)campaigns for observing the Earth’s MP.展开更多
Solar wind charge exchange produces emissions in the soft X-ray energy range which can enable the study of near-Earth space regions such as the magnetopause,the magnetosheath and the polar cusps by remote sensing tech...Solar wind charge exchange produces emissions in the soft X-ray energy range which can enable the study of near-Earth space regions such as the magnetopause,the magnetosheath and the polar cusps by remote sensing techniques.The Solar wind Magnetosphere Ionosphere Link Explorer(SMILE)and Lunar Environment heliospheric X-ray Imager(LEXI)missions aim to obtain soft Xray images of near-Earth space thanks to their Soft X-ray Imager(SXI)instruments.While earlier modeling works have already simulated soft X-ray images as might be obtained by SMILE SXI during its mission,the numerical models used so far are all based on the magnetohydrodynamics description of the space plasma.To investigate the possible signatures of ion-kinetic-scale processes in soft Xray images,we use for the first time a global hybrid-Vlasov simulation of the geospace from the Vlasiator model.The simulation is driven by fast and tenuous solar wind conditions and purely southward interplanetary magnetic field.We first produce global X-ray images of the dayside near-Earth space by placing a virtual imaging satellite at two different locations,providing meridional and equatorial views.We then analyze regional features present in the images and show that they correspond to signatures in soft X-ray emissions of mirrormode wave structures in the magnetosheath and flux transfer events(FTEs)at the magnetopause.Our results suggest that,although the time scales associated with the motion of those transient phenomena will likely be significantly smaller than the integration time of the SMILE and LEXI imagers,mirror-mode structures and FTEs can cumulatively produce detectable signatures in the soft X-ray images.For instance,a local increase by 30%in the proton density at the dayside magnetopause resulting from the transit of multiple FTEs leads to a 12%enhancement in the line-of-sight-and time-integrated soft X-ray emissivity originating from this region.Likewise,a proton density increase by 14%in the magnetosheath associated with mirror-mode structures can result in an enhancement in the soft X-ray signal by 4%.These are likely conservative estimates,given that the solar wind conditions used in the Vlasiator run can be expected to generate weaker soft X-ray emissions than the more common denser solar wind.These results will contribute to the preparatory work for the SMILE and LEXI missions by providing the community with quantitative estimates of the effects of small-scale,transient phenomena occurring on the dayside.展开更多
Spinal cord injury is considered one of the most difficult injuries to repair and has one of the worst prognoses for injuries to the nervous system.Following surgery,the poor regenerative capacity of nerve cells and t...Spinal cord injury is considered one of the most difficult injuries to repair and has one of the worst prognoses for injuries to the nervous system.Following surgery,the poor regenerative capacity of nerve cells and the generation of new scars can make it very difficult for the impaired nervous system to restore its neural functionality.Traditional treatments can only alleviate secondary injuries but cannot fundamentally repair the spinal cord.Consequently,there is a critical need to develop new treatments to promote functional repair after spinal cord injury.Over recent years,there have been seve ral developments in the use of stem cell therapy for the treatment of spinal cord injury.Alongside significant developments in the field of tissue engineering,three-dimensional bioprinting technology has become a hot research topic due to its ability to accurately print complex structures.This led to the loading of three-dimensional bioprinting scaffolds which provided precise cell localization.These three-dimensional bioprinting scaffolds co uld repair damaged neural circuits and had the potential to repair the damaged spinal cord.In this review,we discuss the mechanisms underlying simple stem cell therapy,the application of different types of stem cells for the treatment of spinal cord injury,and the different manufa cturing methods for three-dimensional bioprinting scaffolds.In particular,we focus on the development of three-dimensional bioprinting scaffolds for the treatment of spinal cord injury.展开更多
This study proposed a novel friction law dependent on velocity,displacement and normal stress for kinematic analysis of runout process of rapid landslides.The well-known Yigong landslide occurring in the Tibetan Plate...This study proposed a novel friction law dependent on velocity,displacement and normal stress for kinematic analysis of runout process of rapid landslides.The well-known Yigong landslide occurring in the Tibetan Plateau of China was employed as the case,and the derived dynamic friction formula was included into the numerical simulation based on Particle Flow Code.Results showed that the friction decreased quickly from 0.64(the peak)to 0.1(the stead value)during the 5s-period after the sliding initiation,which explained the behavior of rapid movement of the landslide.The monitored balls set at different sections of the mass showed similar variation characteritics regarding the velocity,namely evident increase at the initial phase of the movement,followed by a fluctuation phase and then a stopping one.The peak velocity was more than 100 m/s and most particles had low velocities at 300s after the landslide initiation.The spreading distance of the landslide was calculated at the two-dimension(profile)and three-dimension scale,respectively.Compared with the simulation result without considering friction weakening effect,our results indicated a max distance of about 10 km from the initial unstable position,which fit better with the actual situation.展开更多
The maintenance of safety and dependability in rail and road embankments is of utmost importance in order to facilitate the smooth operation of transportation networks.This study introduces a comprehensive methodology...The maintenance of safety and dependability in rail and road embankments is of utmost importance in order to facilitate the smooth operation of transportation networks.This study introduces a comprehensive methodology for soil slope stability evaluation,employing Monte Carlo Simulation(MCS)and Subset Simulation(SS)with the"UPSS 3.0 Add-in"in MS-Excel.Focused on an 11.693-meter embankment with a soil slope(inclination ratio of 2H:1V),the investigation considers earthquake coefficients(kh)and pore water pressure ratios(ru)following Indian zoning requirements.The chance of slope failure showed a considerable increase as the Coefficient of Variation(COV),seismic coefficients(kh),and pore water pressure ratios(ru)experienced an escalation.The SS approach showed exceptional efficacy in calculating odds of failure that are notably low.Within computational modeling,the study optimized the worst-case scenario using ANFIS-GA,ANFIS-GWO,ANFIS-PSO,and ANFIS-BBO models.The ANFIS-PSO model exhibits exceptional accuracy(training R2=0.9011,RMSE=0.0549;testing R2=0.8968,RMSE=0.0615),emerging as the most promising.This study highlights the significance of conducting thorough risk assessments and offers practical insights into evaluating and improving the stability of soil slopes in transportation infrastructure.These findings contribute to the enhancement of safety and reliability in real-world situations.展开更多
A toroidal soft x-ray imaging(T-SXRI)system has been developed to investigate threedimensional(3D)plasma physics on J-TEXT.This T-SXRI system consists of three sets of SXR arrays.Two sets are newly developed and locat...A toroidal soft x-ray imaging(T-SXRI)system has been developed to investigate threedimensional(3D)plasma physics on J-TEXT.This T-SXRI system consists of three sets of SXR arrays.Two sets are newly developed and located on the vacuum chamber wall at toroidal positionsφof 126.4°and 272.6°,respectively,while one set was established previously atφ=65.50.Each set of SXR arrays consists of three arrays viewing the plasma poloidally,and hence can be used separately to obtain SXR images via the tomographic method.The sawtooth precursor oscillations are measured by T-SXRI,and the corresponding images of perturbative SXR signals are successfully reconstructed at these three toroidal positions,hence providing measurement of the 3D structure of precursor oscillations.The observed 3D structure is consistent with the helical structure of the m/n=1/1 mode.The experimental observation confirms that the T-SXRI system is able to observe 3D structures in the J-TEXT plasma.展开更多
基金supported by Research on the Influence of Nozzle Structure on the Scouring Effect of Submerged Water Jet(2023R411045)Design and Control Strategy Research of PEM Fuel Cell Hybrid Propulsion System for Ships(2024R411015)+1 种基金Zhejiang Ocean University Outstanding Master’s Thesis Cultivation Project(ZJOUYJS20230018)General Program of Education Department of Zhejiang Province(Y202250817)which was gained by Chen.
文摘Scouring experiments were conducted using a three-dimensional laser scanning technology for angles of the jet spanning the interval from 0°to 30°,and the characteristics of the scour hole in equilibrium conditions were investigated accordingly.The results indicate that the optimal scouring effects occur when the jet angle is in the ranges between 15°and 20°.Moreover,the dimensionless profiles of the scour hole exhibit a high degree of similarity at different jet angles.Numerical simulations conducted using the Flow-3D software to investigate the bed shear stress along the jet impingement surface have shown that this stress is influenced by both the resultant force and the jet impingement surface area.It reaches its maximum value when the jet is vertical,decreases rapidly as the jet starts to tilt,then increases slightly,and decreases again significantly when the angle exceeds 20°.
基金supported by the National Natural Science Foundation of China(Grant Nos.42277150,41977219)Henan Provincial Science and Technology Research Project(Grant No.222102320271).
文摘The geometry of joints has a significant influence on the mechanical properties of rocks.To simplify the curved joint shapes in rocks,the joint shape is usually treated as straight lines or planes in most laboratory experiments and numerical simulations.In this study,the computerized tomography (CT) scanning and photogrammetry were employed to obtain the internal and surface joint structures of a limestone sample,respectively.To describe the joint geometry,the edge detection algorithms and a three-dimensional (3D) matrix mapping method were applied to reconstruct CT-based and photogrammetry-based jointed rock models.For comparison tests,the numerical uniaxial compression tests were conducted on an intact rock sample and a sample with a joint simplified to a plane using the parallel computing method.The results indicate that the mechanical characteristics and failure process of jointed rocks are significantly affected by the geometry of joints.The presence of joints reduces the uniaxial compressive strength (UCS),elastic modulus,and released acoustic emission (AE) energy of rocks by 37%–67%,21%–24%,and 52%–90%,respectively.Compared to the simplified joint sample,the proposed photogrammetry-based numerical model makes the most of the limited geometry information of joints.The UCS,accumulative released AE energy,and elastic modulus of the photogrammetry-based sample were found to be very close to those of the CT-based sample.The UCS value of the simplified joint sample (i.e.38.5 MPa) is much lower than that of the CT-based sample (i.e.72.3 MPa).Additionally,the accumulative released AE energy observed in the simplified joint sample is 3.899 times lower than that observed in the CT-based sample.CT scanning provides a reliable means to visualize the joints in rocks,which can be used to verify the reliability of photogrammetry techniques.The application of the photogrammetry-based sample enables detailed analysis for estimating the mechanical properties of jointed rocks.
基金The National Natural Science Foundation of China (No.52165060,12272189)Program for Young Talents of Science and Technology in Universities of Inner Mongolia Autonomous Region: (NJYT23022)+2 种基金Science and Technology Projects of Inner Mongolia Autonomous Region: (2021GG0432)Central Guiding Local Science and Technology Development Plan (2022ZY0013)Basic research business fee project for universities directly under Inner Mongolia Autonomous Region (GXKY22046).
文摘Aiming at the problems that the simulation accuracy which is reduced due to the simplification of the model,a three-dimensional simulation method based on solid modeling is being proposed.By analyzing the motion relationship and positional relationship between the caries knife and the workpiece,the coordinate system of the caries machining was established.With the MATLAB software,the cutting edge model and the blade sweeping surface model of the boring cutter are sequentially established.Boolean operation is performed on the blade swept surface formed by the tooth cutter teeth with time t and the workpiece tooth geometry as well as the undeformed three-dimensional chip geometry model and the instantaneous cogging geometry model are obtained at different times.Through the compare between gear end face simulation tooth profile and the theoretical inner arc tooth profile,we verified the accuracy and rationality of the proposed method.
基金Supported by the National Natural Science Foundation of China(51827804)CNPC Strategic Cooperation Science and Technology Major Project(ZLZX2020-01-05)Open Fund of State Key Laboratory of Rock Mechanics and Engineering(SKLGME021024).
文摘A fracture propagation model of radial well fracturing is established based on the finite element-meshless method.The model considers the coupling effect of fracturing fluid flow and rock matrix deformation.The fracture geometries of radial well fracturing are simulated,the induction effect of radial well on the fracture is quantitatively characterized,and the influences of azimuth,horizontal principle stress difference,and reservoir matrix permeability on the fracture geometries are revealed.The radial wells can induce the fractures to extend parallel to their axes when two radial wells in the same layer are fractured.When the radial wells are symmetrically distributed along the direction of the minimum horizontal principle stress with the azimuth greater than 15,the extrusion effect reduces the fracture length of radial wells.When the radial wells are symmetrically distributed along the direction of the maximum horizontal principal stress,the extrusion increases the fracture length of the radial wells.The fracture geometries are controlled by the rectification of radial borehole,the extrusion between radial wells in the same layer,and the deflection of the maximum horizontal principal stress.When the radial wells are distributed along the minimum horizontal principal stress symmetrically,the fracture length induced by the radial well decreases with the increase of azimuth;in contrast,when the radial wells are distributed along the maximum horizontal principal stress symmetrically,the fracture length induced by the radial well first decreases and then increases with the increase of azimuth.The fracture length induced by the radial well decreases with the increase of horizontal principal stress difference.The increase of rock matrix permeability and pore pressure of the matrix around radial wells makes the inducing effect of the radial well on fractures increase.
基金Funding Statement:This work was supported by the National Key Research and Development Program of China[Grant No.2019YFC1906802].
文摘Organic solid waste(OSW)contains many renewable materials.The pyrolysis and gasification of OSW can realize resource utilization,and its products can be used for methanation reaction to produce synthetic natural gas in the specific reactor.In order to understand the dynamic characteristics of the reactor,a three-dimensional numerical model has been established by the method of Computational Fluid Dynamics(CFD).Along the height of the reactor,the particle distribution in the bed becomes thinner and the mean solid volume fraction decreases from 4.18%to 0.37%.Meanwhile,the pressure fluctuation range decreased from 398.76 Pa at the entrance to a much lower value of 74.47 Pa at the exit.In this simulation,three parameters of gas inlet velocity,operating temperature and solid particle diameter are changed to explore their influences on gas-solid multiphase flow.The results show that gas velocity has a great influence on particle distribution.When the gas inlet velocity decreases from 6.51 to 1.98 m/s,the minimum height that particles can reach decreases from 169 to 100 mm.Additionally,as the operating temperature increases,the particle holdup inside the reactor changes from 0.843%to 0.700%.This indicates that the particle residence time reduces,which is not conducive to the follow-up reaction.Moreover,with the increase of particle size,the fluctuation range of the pressure at the bottom of the reactor increases,and its standard deviation increases from 55.34 to 1266.37 Pa.
基金funded by National Key R&D Program of China(No.2021YFB3401200)the National Natural Science Foundation of China(No.51875308)the Beijing Nature Sciences Fund-Haidian Originality Cooperation Project(L212002).
文摘Numerical simulation is the most powerful computational and analysis tool for a large variety of engineering and physical problems.For a complex problem relating to multi-field,multi-process and multi-scale,different computing tools have to be developed so as to solve particular fields at different scales and for different processes.Therefore,the integration of different types of software is inevitable.However,it is difficult to perform the transfer of the meshes and simulated results among software packages because of the lack of shared data formats or encrypted data formats.An image processing based method for three-dimensional model reconstruction for numerical simulation was proposed,which presents a solution to the integration problem by a series of slice or projection images obtained by the post-processing modules of the numerical simulation software.By means of mapping image pixels to meshes of either finite difference or finite element models,the geometry contour can be extracted to export the stereolithography model.The values of results,represented by color,can be deduced and assigned to the meshes.All the models with data can be directly or indirectly integrated into other software as a continued or new numerical simulation.The three-dimensional reconstruction method has been validated in numerical simulation of castings and case studies were provided in this study.
基金financially supported by the National Key R&D Program of China (Grant No.2022YFB2603000)the National Natural Science Foundation of China (Grant Nos.51779172 and 52179076)。
文摘Three-dimensional direct numerical simulations of the wake flow downstream of a near-wall circular cylinder at different gap ratios and boundary layer thicknesses are carried out by using the iterative immersed boundary method.The non-dimensional gap between the cylinder and the wall,G/D=0.2,0.6 and 1.0,the non-dimensional boundary layer thickness,δ/D=0.0,0.7 and 1.6,the Reynolds number,Re=350,and the aspect ratio of the cylinder,L/D=25are adopted.High-resolution visualizations of the complex vortex structures at differentδ/D and G/D are presented.The transition of the streamwise vortex mode,the combined effects ofδ/D and G/D on the flow statistics,the pressure and shear stress distribution and the hydrodynamic forces are analyzed.Results show that with decreasing G/D and increasingδ/D,the gap flow and its vortex-shedding are significantly weakened,together with an elongated wake and an enlarged low-velocity area near the wall,leading to the wake mode transition from the two-sided to one-sided vortex-shedding.Different relative positions of the cylinder regarding the boundary layer alter the flow features of the shear layers.With an increase inδ/D,the front stagnation point shifts to the upper surface,and the distance between the flow divergence point and the maximum pressure position increases.The mean drag coefficient and r.m.s.values of drag and lift coefficients decrease with a decrease in G/D and an increase inδ/D,while the mean lift coefficient increases with decreasing G/D but decreases with increasingδ/D.
基金supported by the National Natural Science Foundation of China (No. 52275291)the Fundamental Research Funds for the Central Universitiesthe Program for Innovation Team of Shaanxi Province,China (No. 2023-CX-TD-17)
文摘Hypoxia is a typical feature of the tumor microenvironment,one of the most critical factors affecting cell behavior and tumor progression.However,the lack of tumor models able to precisely emulate natural brain tumor tissue has impeded the study of the effects of hypoxia on the progression and growth of tumor cells.This study reports a three-dimensional(3D)brain tumor model obtained by encapsulating U87MG(U87)cells in a hydrogel containing type I collagen.It also documents the effect of various oxygen concentrations(1%,7%,and 21%)in the culture environment on U87 cell morphology,proliferation,viability,cell cycle,apoptosis rate,and migration.Finally,it compares two-dimensional(2D)and 3D cultures.For comparison purposes,cells cultured in flat culture dishes were used as the control(2D model).Cells cultured in the 3D model proliferated more slowly but had a higher apoptosis rate and proportion of cells in the resting phase(G0 phase)/gap I phase(G1 phase)than those cultured in the 2D model.Besides,the two models yielded significantly different cell morphologies.Finally,hypoxia(e.g.,1%O2)affected cell morphology,slowed cell growth,reduced cell viability,and increased the apoptosis rate in the 3D model.These results indicate that the constructed 3D model is effective for investigating the effects of biological and chemical factors on cell morphology and function,and can be more representative of the tumor microenvironment than 2D culture systems.The developed 3D glioblastoma tumor model is equally applicable to other studies in pharmacology and pathology.
基金supported by the National Natural Science Foundation of China(Grant Nos.42175099,42027804,42075073)the Innovative Project of Postgraduates in Jiangsu Province in 2023(Grant No.KYCX23_1319)+3 种基金supported by the National Natural Science Foundation of China(Grant No.42205080)the Natural Science Foundation of Sichuan(Grant No.2023YFS0442)the Research Fund of Civil Aviation Flight University of China(Grant No.J2022-037)supported by the National Key Scientific and Technological Infrastructure project“Earth System Science Numerical Simulator Facility”(Earth Lab)。
文摘The process of entrainment-mixing between cumulus clouds and the ambient air is important for the development of cumulus clouds.Accurately obtaining the entrainment rate(λ)is particularly important for its parameterization within the overall cumulus parameterization scheme.In this study,an improved bulk-plume method is proposed by solving the equations of two conserved variables simultaneously to calculateλof cumulus clouds in a large-eddy simulation.The results demonstrate that the improved bulk-plume method is more reliable than the traditional bulk-plume method,becauseλ,as calculated from the improved method,falls within the range ofλvalues obtained from the traditional method using different conserved variables.The probability density functions ofλfor all data,different times,and different heights can be well-fitted by a log-normal distribution,which supports the assumed stochastic entrainment process in previous studies.Further analysis demonstrate that the relationship betweenλand the vertical velocity is better than other thermodynamic/dynamical properties;thus,the vertical velocity is recommended as the primary influencing factor for the parameterization ofλin the future.The results of this study enhance the theoretical understanding ofλand its influencing factors and shed new light on the development ofλparameterization.
基金supported by the National Natural Science Foundation of China(NNSFC)grants 42074202,42274196Strategic Priority Research Program of Chinese Academy of Sciences grant XDB41000000ISSI-BJ International Team Interaction between magnetic reconnection and turbulence:From the Sun to the Earth。
文摘Earth’s magnetopause is a thin boundary separating the shocked solar wind plasma from the magnetospheric plasmas,and it is also the boundary of the solar wind energy transport to the magnetosphere.Soft X-ray imaging allows investigation of the large-scale magnetopause by providing a two-dimensional(2-D)global view from a satellite.By performing 3-D global hybrid-particle-in-cell(hybrid-PIC)simulations,we obtain soft X-ray images of Earth’s magnetopause under different solar wind conditions,such as different plasma densities and directions of the southward interplanetary magnetic field.In all cases,magnetic reconnection occurs at low latitude magnetopause.The soft X-ray images observed by a hypothetical satellite are shown,with all of the following identified:the boundary of the magnetopause,the cusps,and the magnetosheath.Local X-ray emissivity in the magnetosheath is characterized by large amplitude fluctuations(up to 160%);however,the maximum line-of-sight-integrated X-ray intensity matches the tangent directions of the magnetopause well,indicating that these fluctuations have limited impact on identifying the magnetopause boundary in the X-ray images.Moreover,the magnetopause boundary can be identified using multiple viewing geometries.We also find that solar wind conditions have little effect on the magnetopause identification.The Solar wind Magnetosphere Ionosphere Link Explorer(SMILE)mission will provide X-ray images of the magnetopause for the first time,and our global hybrid-PIC simulation results can help better understand the 2-D X-ray images of the magnetopause from a 3-D perspective,with particle kinetic effects considered.
基金support from the European Space Agency(ESA)PRODEX(PROgramme de Développement d’Expériences scientifiques)Project mission(No.PEA4000134960)Partial funding was provided by the Romanian Ministry of Research,Innovation and Digitalization under Romanian National Core Program LAPLAS VII(Contract No.30N/2023)+2 种基金the Belgian Solar-Terrestrial Centre of Excellencesupported by the project Belgian Research Action through Interdisciplinary Networks(BRAIN-BE)2.0(Grant No.B2/223/P1/PLATINUM)funded by the Belgian Office for Research(BELSPO)partially supported by a grant from the Romanian Ministry of Education and Research(CNCS-UEFISCDI,Project No.PN-III-P1-1.1TE-2021-0102)。
文摘We derived the properties of the terrestrial magnetopause(MP)from two modeling approaches,one global–fluid,the other local–kinetic,and compared the results with data collected in situ by the Magnetospheric Multiscale 2(MMS2)spacecraft.We used global magnetohydrodynamic(MHD)simulations of the Earth’s magnetosphere(publicly available from the NASA-CCMC[National Aeronautics and Space Administration–Community Coordinated Modeling Center])and local Vlasov equilibrium models(based on kinetic models for tangential discontinuities)to extract spatial profiles of the plasma and field variables at the Earth’s MP.The global MHD simulations used initial solar wind conditions extracted from the OMNI database at the time epoch when the MMS2 observes the MP.The kinetic Vlasov model used asymptotic boundary conditions derived from the same in situ MMS measurements upstream or downstream of the MP.The global MHD simulations provide a three-dimensional image of the magnetosphere at the time when the MMS2 crosses the MP.The Vlasov model provides a one-dimensional local view of the MP derived from first principles of kinetic theory.The MMS2 experimental data also serve as a reference for comparing and validating the numerical simulations and modeling.We found that the MP transition layer formed in global MHD simulations was generally localized closer to the Earth(roughly by one Earth radius)from the position of the real MP observed by the MMS.We also found that the global MHD simulations overestimated the thickness of the MP transition by one order of magnitude for three analyzed variables:magnetic field,density,and tangential speed.The MP thickness derived from the local Vlasov equilibrium was consistent with observations for all three of these variables.The overestimation of density in the Vlasov equilibrium was reduced compared with the global MHD solutions.We discuss our results in the context of future SMILE(Solar wind Magnetosphere Ionosphere Link Explorer)campaigns for observing the Earth’s MP.
基金the European Research Council for starting grant 200141-QuESpace,with which the Vlasiator model was developedconsolidator grant 682068-PRESTISSIMO awarded for further development of Vlasiator and its use in scientific investigations+4 种基金Academy of Finland grant numbers 338629-AERGELC’H,339756-KIMCHI,336805-FORESAIL,and 335554-ICT-SUNVACThe Academy of Finland also supported this work through the PROFI4 grant(grant number 3189131)support from the NASA grants,80NSSC20K1670 and 80MSFC20C0019the NASA GSFC FY23 IRADHIF funds。
文摘Solar wind charge exchange produces emissions in the soft X-ray energy range which can enable the study of near-Earth space regions such as the magnetopause,the magnetosheath and the polar cusps by remote sensing techniques.The Solar wind Magnetosphere Ionosphere Link Explorer(SMILE)and Lunar Environment heliospheric X-ray Imager(LEXI)missions aim to obtain soft Xray images of near-Earth space thanks to their Soft X-ray Imager(SXI)instruments.While earlier modeling works have already simulated soft X-ray images as might be obtained by SMILE SXI during its mission,the numerical models used so far are all based on the magnetohydrodynamics description of the space plasma.To investigate the possible signatures of ion-kinetic-scale processes in soft Xray images,we use for the first time a global hybrid-Vlasov simulation of the geospace from the Vlasiator model.The simulation is driven by fast and tenuous solar wind conditions and purely southward interplanetary magnetic field.We first produce global X-ray images of the dayside near-Earth space by placing a virtual imaging satellite at two different locations,providing meridional and equatorial views.We then analyze regional features present in the images and show that they correspond to signatures in soft X-ray emissions of mirrormode wave structures in the magnetosheath and flux transfer events(FTEs)at the magnetopause.Our results suggest that,although the time scales associated with the motion of those transient phenomena will likely be significantly smaller than the integration time of the SMILE and LEXI imagers,mirror-mode structures and FTEs can cumulatively produce detectable signatures in the soft X-ray images.For instance,a local increase by 30%in the proton density at the dayside magnetopause resulting from the transit of multiple FTEs leads to a 12%enhancement in the line-of-sight-and time-integrated soft X-ray emissivity originating from this region.Likewise,a proton density increase by 14%in the magnetosheath associated with mirror-mode structures can result in an enhancement in the soft X-ray signal by 4%.These are likely conservative estimates,given that the solar wind conditions used in the Vlasiator run can be expected to generate weaker soft X-ray emissions than the more common denser solar wind.These results will contribute to the preparatory work for the SMILE and LEXI missions by providing the community with quantitative estimates of the effects of small-scale,transient phenomena occurring on the dayside.
基金supported by the National Natural Science Foundation of China,No.82171380(to CD)Jiangsu Students’Platform for Innovation and Entrepreneurship Training Program,No.202110304098Y(to DJ)。
文摘Spinal cord injury is considered one of the most difficult injuries to repair and has one of the worst prognoses for injuries to the nervous system.Following surgery,the poor regenerative capacity of nerve cells and the generation of new scars can make it very difficult for the impaired nervous system to restore its neural functionality.Traditional treatments can only alleviate secondary injuries but cannot fundamentally repair the spinal cord.Consequently,there is a critical need to develop new treatments to promote functional repair after spinal cord injury.Over recent years,there have been seve ral developments in the use of stem cell therapy for the treatment of spinal cord injury.Alongside significant developments in the field of tissue engineering,three-dimensional bioprinting technology has become a hot research topic due to its ability to accurately print complex structures.This led to the loading of three-dimensional bioprinting scaffolds which provided precise cell localization.These three-dimensional bioprinting scaffolds co uld repair damaged neural circuits and had the potential to repair the damaged spinal cord.In this review,we discuss the mechanisms underlying simple stem cell therapy,the application of different types of stem cells for the treatment of spinal cord injury,and the different manufa cturing methods for three-dimensional bioprinting scaffolds.In particular,we focus on the development of three-dimensional bioprinting scaffolds for the treatment of spinal cord injury.
基金funded by the National Natural Science Foundation of China(42307248,U23A2047,42277187)Natural Science Foundation of Hebei Province(D2022202005)+1 种基金Planning and Natural Resources Research Project of Tianjin City(2022-40,KJ[2024]25)the support from the Graduated Student Innovation Funding Project of Hebei Province(CXZZSS2024007)。
文摘This study proposed a novel friction law dependent on velocity,displacement and normal stress for kinematic analysis of runout process of rapid landslides.The well-known Yigong landslide occurring in the Tibetan Plateau of China was employed as the case,and the derived dynamic friction formula was included into the numerical simulation based on Particle Flow Code.Results showed that the friction decreased quickly from 0.64(the peak)to 0.1(the stead value)during the 5s-period after the sliding initiation,which explained the behavior of rapid movement of the landslide.The monitored balls set at different sections of the mass showed similar variation characteritics regarding the velocity,namely evident increase at the initial phase of the movement,followed by a fluctuation phase and then a stopping one.The peak velocity was more than 100 m/s and most particles had low velocities at 300s after the landslide initiation.The spreading distance of the landslide was calculated at the two-dimension(profile)and three-dimension scale,respectively.Compared with the simulation result without considering friction weakening effect,our results indicated a max distance of about 10 km from the initial unstable position,which fit better with the actual situation.
文摘The maintenance of safety and dependability in rail and road embankments is of utmost importance in order to facilitate the smooth operation of transportation networks.This study introduces a comprehensive methodology for soil slope stability evaluation,employing Monte Carlo Simulation(MCS)and Subset Simulation(SS)with the"UPSS 3.0 Add-in"in MS-Excel.Focused on an 11.693-meter embankment with a soil slope(inclination ratio of 2H:1V),the investigation considers earthquake coefficients(kh)and pore water pressure ratios(ru)following Indian zoning requirements.The chance of slope failure showed a considerable increase as the Coefficient of Variation(COV),seismic coefficients(kh),and pore water pressure ratios(ru)experienced an escalation.The SS approach showed exceptional efficacy in calculating odds of failure that are notably low.Within computational modeling,the study optimized the worst-case scenario using ANFIS-GA,ANFIS-GWO,ANFIS-PSO,and ANFIS-BBO models.The ANFIS-PSO model exhibits exceptional accuracy(training R2=0.9011,RMSE=0.0549;testing R2=0.8968,RMSE=0.0615),emerging as the most promising.This study highlights the significance of conducting thorough risk assessments and offers practical insights into evaluating and improving the stability of soil slopes in transportation infrastructure.These findings contribute to the enhancement of safety and reliability in real-world situations.
基金supported by the National Magnetic Confinement Fusion Energy R&D Program of China(Nos.2018YFE0309100 and 2019YFE03010004)National Natural Science Foundation of China(No.51821005)。
文摘A toroidal soft x-ray imaging(T-SXRI)system has been developed to investigate threedimensional(3D)plasma physics on J-TEXT.This T-SXRI system consists of three sets of SXR arrays.Two sets are newly developed and located on the vacuum chamber wall at toroidal positionsφof 126.4°and 272.6°,respectively,while one set was established previously atφ=65.50.Each set of SXR arrays consists of three arrays viewing the plasma poloidally,and hence can be used separately to obtain SXR images via the tomographic method.The sawtooth precursor oscillations are measured by T-SXRI,and the corresponding images of perturbative SXR signals are successfully reconstructed at these three toroidal positions,hence providing measurement of the 3D structure of precursor oscillations.The observed 3D structure is consistent with the helical structure of the m/n=1/1 mode.The experimental observation confirms that the T-SXRI system is able to observe 3D structures in the J-TEXT plasma.