With respect to oceanic fluid dynamics,certain models have appeared,e.g.,an extended time-dependent(3+1)-dimensional shallow water wave equation in an ocean or a river,which we investigate in this paper.Using symbolic...With respect to oceanic fluid dynamics,certain models have appeared,e.g.,an extended time-dependent(3+1)-dimensional shallow water wave equation in an ocean or a river,which we investigate in this paper.Using symbolic computation,we find out,on one hand,a set of bilinear auto-Backlund transformations,which could connect certain solutions of that equation with other solutions of that equation itself,and on the other hand,a set of similarity reductions,which could go from that equation to a known ordinary differential equation.The results in this paper depend on all the oceanic variable coefficients in that equation.展开更多
This paper proposes a novel cargo loading algorithm applicable to automated conveyor-type loading systems.The algorithm offers improvements in computational efficiency and robustness by utilizing the concept of discre...This paper proposes a novel cargo loading algorithm applicable to automated conveyor-type loading systems.The algorithm offers improvements in computational efficiency and robustness by utilizing the concept of discrete derivatives and introducing logistics-related constraints.Optional consideration of the rotation of the cargoes was made to further enhance the optimality of the solutions,if possible to be physically implemented.Evaluation metrics were developed for accurate evaluation and enhancement of the algorithm’s ability to efficiently utilize the loading space and provide a high level of dynamic stability.Experimental results demonstrate the extensive robustness of the proposed algorithm to the diversity of cargoes present in Business-to-Consumer environments.This study contributes practical advancements in both cargo loading optimization and automation of the logistics industry,with potential applications in last-mile delivery services,warehousing,and supply chain management.展开更多
Lie symmetry analysis is applied to a(3+1)-dimensional combined potential Kadomtsev-Petviashvili equation with B-type Kadomtsev-Petviashvili equation(pKP-BKP equation)and the corresponding similarity reduction equatio...Lie symmetry analysis is applied to a(3+1)-dimensional combined potential Kadomtsev-Petviashvili equation with B-type Kadomtsev-Petviashvili equation(pKP-BKP equation)and the corresponding similarity reduction equations are obtained with the different infinitesimal generators.Invariant solutions with arbitrary functions and constants for the(3+1)-dimensional pKP-BKP equation,including the lump solution,the periodic-lump solution,the two-kink solution,the breather solution and the lump-two-kink solution,have been studied analytically and graphically.展开更多
Contact detection is the most time-consuming stage in 3D discontinuous deformation analysis(3D-DDA)computation.Improving the efficiency of 3D-DDA is beneficial for its application in large-scale computing.In this stud...Contact detection is the most time-consuming stage in 3D discontinuous deformation analysis(3D-DDA)computation.Improving the efficiency of 3D-DDA is beneficial for its application in large-scale computing.In this study,aiming at the continuous-discontinuous simulation of 3D-DDA,a highly efficient contact detection strategy is proposed.Firstly,the global direct search(GDS)method is integrated into the 3D-DDA framework to address intricate contact scenarios.Subsequently,all geometric elements,including blocks,faces,edges,and vertices are divided into searchable and unsearchable parts.Contacts between unsearchable geometric elements would be directly inherited,while only searchable geometric elements are involved in contact detection.This strategy significantly reduces the number of geometric elements involved in contact detection,thereby markedly enhancing the computation efficiency.Several examples are adopted to demonstrate the accuracy and efficiency of the improved 3D-DDA method.The rock pillars with different mesh sizes are simulated under self-weight.The deformation and stress are consistent with the analytical results,and the smaller the mesh size,the higher the accuracy.The maximum speedup ratio is 38.46 for this case.Furthermore,the Brazilian splitting test on the discs with different flaws is conducted.The results show that the failure pattern of the samples is consistent with the results obtained by other methods and experiments,and the maximum speedup ratio is 266.73.Finally,a large-scale impact test is performed,and approximately 3.2 times enhanced efficiency is obtained.The proposed contact detection strategy significantly improves efficiency when the rock has not completely failed,which is more suitable for continuous-discontinuous simulation.展开更多
Current transient analysis predominantly relies on zero-dimensional/one-dimensional tools,proficient at capturing aerothermodynamic variations across critical engine stations but insufficient for analyzing the interna...Current transient analysis predominantly relies on zero-dimensional/one-dimensional tools,proficient at capturing aerothermodynamic variations across critical engine stations but insufficient for analyzing the internal flow field evolution during transients.Addressing this gap,the study presents an enhanced quasi-three dimensional(quasi-3D)transient simulation technique that integrates component volume effects,offering a significant leap from the preceding quasi-3D transient simulation method based on quasi-steady assumption.By embedding the component volume effects on density,momentum,and energy within the physical temporal dimension of the Navier-Stokes equations,the refined quasi-3D transient model achieves a closer representation of physical phenomena.Validation against a single-shaft turbofan engine’s experimental data confirms the model’s accuracy.Average errors for key performance indicators,including shaft speed,thrust,mass flow rate,and critical component exit temperature and pressure,remain below 0.41%,5.69%,2.55%,3.18%and 0.67%,respectively.Crucially,the model exposes a discernible temporal lag in the compressor outlet pressure and temperature response due to volume effects—previously unquantified in quasi-3D transient simulations.And further exploration of the meridional flow field emphasizes the consequential role of volumes in transient flow field evolution.Incorporating volume effects within quasi-3D transient simulations enhances engine modeling and is pivotal for precise transient analysis in engine design and optimization.展开更多
A accelerated arithmetic algorithm of the dynamic computing regions was designed,and 3-dimensional numerical simulation of isothermal solidification for a binary alloy was implemented.The dendritic growth and the reca...A accelerated arithmetic algorithm of the dynamic computing regions was designed,and 3-dimensional numerical simulation of isothermal solidification for a binary alloy was implemented.The dendritic growth and the recalescence of Ni-Cu binary alloy during the solidification at different cooling rates were investigated.The effects of cooling rate on dendritic patterns and microsegregation patterns were studied.The computed results indicate that,with the increment of the cooling rate,the dendritic growth velocity increases,both the main branch and side-branches become slender,the secondary dendrite arm spacing becomes smaller,the inadequate solute diffusion in solid aggravates,and the severity of microsegregation ahead of interface aggravates.At a higher cooling rate,the binary alloy presents recalescence;while the cooling rate is small,no recalescence occurs.展开更多
A new quasi-three-dimensional （quasi-3D） numeric simulation method for a high-voltage level-shifting circuit structure is proposed. The performances of the 3D structure are analyzed by combining some 2D device struc...A new quasi-three-dimensional （quasi-3D） numeric simulation method for a high-voltage level-shifting circuit structure is proposed. The performances of the 3D structure are analyzed by combining some 2D device structures; the 2D devices are in two planes perpendicular to each other and to the surface of the semiconductor. In comparison with Davinci, the full 3D device simulation tool, the quasi-3D simulation method can give results for the potential and current distribution of the 3D high-voltage level-shifting circuit structure with appropriate accuracy and the total CPU time for simulation is significantly reduced. The quasi-3D simulation technique can be used in many cases with advantages such as saving computing time, making no demands on the high-end computer terminals, and being easy to operate.展开更多
Upon the conservation of mass, momentum and energy, volume fraction and surface penetrative rate were employed to modify the conservative equations to simulate the effect of blockages on fluid flows and heat transfer....Upon the conservation of mass, momentum and energy, volume fraction and surface penetrative rate were employed to modify the conservative equations to simulate the effect of blockages on fluid flows and heat transfer. These equations were solved numerically with the finite differential method and the primitive variable approach. This method uses staggered grid and pressure correction schemes. A computer code FASTOR3D integrated the aforementioned algorithm. The preliminary results have been compared with conventional benchmark solutions. With auxiliary software DV, the numerical results were visualized in colorful images to demonstrate the variation of flow patterns and temperature profiles during the transient process. The results of the simulation code for the fluid flows and heat transfer in the sodium pool of a fast breeder reactor are acceptable.展开更多
Three-dimensional(3D)braided composites have significant potential for use in engineering structural materials.However,conventional 3D braiding machines are insufficient for designing composites with complex geometrie...Three-dimensional(3D)braided composites have significant potential for use in engineering structural materials.However,conventional 3D braiding machines are insufficient for designing composites with complex geometries.This paper proposes a programmable design methodology for 3D rotary braiding machines using circle-cutting and combination strategies.By introducing varying numbers of incisions on the circle,a diverse range of horn gears can be designed.Different combinations of these cut-circles allow the horn gears to be assembled into various 3D rotary braiders.The parametric equation for the braider plate is derived,showing that a combination strategy involving two cut-circles is feasible for braider design,whereas integrating three cut-circles simultaneously is impossible for a single machine.The construction of an automatic 6-3 type 3D braiding machine demonstrates the effectiveness of the proposed design strategy.This flexible braider design approach provides a practical solution for producing 3D braided composites with complex geometries.展开更多
Convectively unstable processes caused by dense water subsidence are common occurrences in high-latitude oceanic regions,and significantly modulate mass and heat transport and mixing processes in the ocean.An idealize...Convectively unstable processes caused by dense water subsidence are common occurrences in high-latitude oceanic regions,and significantly modulate mass and heat transport and mixing processes in the ocean.An idealized numerical experiment using the large eddy simulation method was conducted to analyze the three-dimensional flow field structure and the mechanism for dense water subsidence.Specifically,a negative salt flux is set at the sea surface,in which salt flux enters the sea surface to simulate the icing and salting-out phenomena that occur at high latitudes.Results show that the mean-state 3D flow field of dense water subsidence exhibits a hollow conical distribution.The horizontal flow field is characterized by a cyclonic vortex that driven primarily by the pressure gradient and influenced by the Coriolis effect.Moreover,the inverse vertical pressure gradient generated by this vortex inhibits the sinking of the plume,leading to its off-axis deflection and the development of an anticyclonic precession.In addition,the impact of rotation on the structure of a sinking plume within a stratified environment is discussed.Both horizontal vortex intensity and cone angle of the hollow cone flow field are increased with increasing rotation rate,resulting in a decrease in the plume’s maximum sinking depth.Variances in rotation direction cause the horizontal vortex and sinking plumes of dense water in the northern and southern hemispheres to rotate in opposite directions.展开更多
Wire arc additive manufacturing(WAAM)has emerged as a promising technique for producing large-scale metal components,favoured by high deposition rates,flexibility and low cost.Despite its potential,the complexity of W...Wire arc additive manufacturing(WAAM)has emerged as a promising technique for producing large-scale metal components,favoured by high deposition rates,flexibility and low cost.Despite its potential,the complexity of WAAM processes,which involves intricate thermal dynamics,phase transitions,and metallurgical,mechanical,and chemical interactions,presents considerable challenges in final product qualities.Simulation technologies in WAAM have proven invaluable,providing accurate predictions in key areas such as material properties,defect identification,deposit morphology,and residual stress.These predictions play a critical role in optimising manufacturing strategies for the final product.This paper provides a comprehensive review of the simulation techniques applied in WAAM,tracing developments from 2013 to 2023.Initially,it analyses the current challenges faced by simulation methods in three main areas.Subsequently,the review explores the current modelling approaches and the applications of these simulations.Following this,the paper discusses the present state of WAAM simulation,identifying specific issues inherent to WAAM simulation itself.Finally,through a thorough review of existing literature and related analysis,the paper offers future perspectives on potential advancements in WAAM simulation strategies.展开更多
Drilling and blasting,characterized by their efficiency,ubiquity,and cost-effectiveness,have emerged as predominant techniques in rock excavation;however,they are accompanied by enormous destructive power.Accurately c...Drilling and blasting,characterized by their efficiency,ubiquity,and cost-effectiveness,have emerged as predominant techniques in rock excavation;however,they are accompanied by enormous destructive power.Accurately controlling the blasting energy and achieving the directional fracture of a rock mass have become common problems in the field.A two-dimensional blasting(2D blasting)technique was proposed that utilizes the characteristic that the tensile strength of a rock mass is significantly lower than its compressive strength.After blasting,only a 2D crack surface is generated along the predetermined direction,eliminating the damage to the reserved rock mass caused by conventional blasting.However,the interior of a natural rock mass is a"black box",and the process of crack propagation is difficult to capture,resulting in an unclear 2D blasting mechanism.To this end,a single-hole polymethyl methacrylate(PMMA)test piece was used to conduct a 2D blasting experiment with the help of a high-speed camera to capture the dynamic crack propagation process and the digital image correlation(DIC)method to analyze the evolution law of surface strain on the test piece.On this basis,a three-dimensional(3D)finite element model was established based on the progressive failure theory to simulate the stress,strain,damage,and displacement evolution process of the model under 2D blasting.The simulation results were consistent with the experimental results.The research results reveal the 2D blasting mechanism and provide theoretical support for the application of 2D blasting technology in the field of rock excavation.展开更多
To further research on high-parameter plasma,we plan to develop a two-dimensional hard X-ray(HXR)imaging system at the HL-3 tokamak to measure HXRs with energies ranging from 20 to 300 keV.The application of an array-...To further research on high-parameter plasma,we plan to develop a two-dimensional hard X-ray(HXR)imaging system at the HL-3 tokamak to measure HXRs with energies ranging from 20 to 300 keV.The application of an array-structured detector ensures that this system can measure HXR-radiation spectra from the entire plasma cross section.Therefore,it is suitable for the study of fast-electron physics,such as radio-frequency wave current drives,fast electrons driving instabilities,and plasma disruptions in fusion research.In this study,we develop a simulation for calculating fast-electron bremsstrahlung in the HL-3 tokamak based on the Monte Carlo simulation code Geant4,in which the plasma geometry and forward scattering of fast-electron bremsstrahlung are considered.The preliminary calculation results indicate that the HXR energy deposi-tion on the detector is symmetrically distributed,even though the plasma distribution is asymmetric owing to the toroidal effect.These simulation results are helpful in constructing the relationship between the energy deposition on the detector and parameter distribution on the plasma cross section during HL-3 experiments.This is beneficial for the reconstruction of the fast-electron-distribution function and for optimizing the design of the HXR-imaging system.展开更多
Lie group analysis method is applied to the extended(3+1)-dimensional Kadomtsev–Petviashvili–Boussinesq equation and the corresponding similarity reduction equations are obtained with various infinitesimal generator...Lie group analysis method is applied to the extended(3+1)-dimensional Kadomtsev–Petviashvili–Boussinesq equation and the corresponding similarity reduction equations are obtained with various infinitesimal generators.By selecting suitable arbitrary functions in the similarity reduction solutions,we obtain abundant invariant solutions,including the trigonometric solution,the kink-lump interaction solution,the interaction solution between lump wave and triangular periodic wave,the two-kink solution,the lump solution,the interaction between a lump and two-kink and the periodic lump solution in different planes.These exact solutions are also given graphically to show the detailed structures of this high dimensional integrable system.展开更多
In this paper,we mainly focus on proving the existence of lump solutions to a generalized(3+1)-dimensional nonlinear differential equation.Hirota’s bilinear method and a quadratic function method are employed to deri...In this paper,we mainly focus on proving the existence of lump solutions to a generalized(3+1)-dimensional nonlinear differential equation.Hirota’s bilinear method and a quadratic function method are employed to derive the lump solutions localized in the whole plane for a(3+1)-dimensional nonlinear differential equation.Three examples of such a nonlinear equation are presented to investigate the exact expressions of the lump solutions.Moreover,the 3d plots and corresponding density plots of the solutions are given to show the space structures of the lump waves.In addition,the breath-wave solutions and several interaction solutions of the(3+1)-dimensional nonlinear differential equation are obtained and their dynamics are analyzed.展开更多
Radiative cooling materials have gained prominence as a zero-energy solution for mitigating global warming.However,a comprehensive understanding of the atomic-scale optical properties and macroscopic optical performan...Radiative cooling materials have gained prominence as a zero-energy solution for mitigating global warming.However,a comprehensive understanding of the atomic-scale optical properties and macroscopic optical performance of radiative cooling materials remains elusive,limiting insight into the underlying physics of their optical response and cooling efficacy.La_(2)O_(3)and HfO_(2),which represent rare earth and third/fourth subgroup inorganic oxides,respectively,show promise for radiative cooling applications.In this study,we used multiscale simulations to investigate the optical properties of La_(2)O_(3)and HfO_(2)across a broad spectrum.First-principles calculations revealed their dielectric functions and intrinsic refractive indices,and the results indicated that the slightly smaller bandgap of La_(2)O_(3)compared to HfO_(2)induces a higher refractive index in the solar band.Additionally,three-phonon scattering was found to provide more accurate infrared optical properties than two-phonon scattering,which enhanced the emissivity in the sky window.Monte Carlo simulations were also used to determine the macroscopic optical properties of La_(2)O_(3)and HfO_(2)coatings.Based on the simulated results,we identified that the particle size and particle volume fraction play a dominant role in the optical properties.Our findings underscore the potential of La_(2)O_(3)and HfO_(2)nanocomposites for environment-friendly cooling and offer a new approach for high-throughput screening of optical materials through multiscale simulations.展开更多
The perfluorosulfonic acid(PFSA) membrane doped with two-dimensional conductive filler Ti_(3)C_(2)T_(x) is a fuel cell proton exchange membrane with high application potential. Experimental studies showed that the pro...The perfluorosulfonic acid(PFSA) membrane doped with two-dimensional conductive filler Ti_(3)C_(2)T_(x) is a fuel cell proton exchange membrane with high application potential. Experimental studies showed that the proton conductivity of Nafion/Ti_(3)C_(2)T_(x) composite membrane is improved significantly compared with that in pure Nafion. However, the microscopic mechanism of doping on the enhancement of membrane performance is remain unclear now. In this work, molecular dynamics simulation was used to investigate the microscopic morphology and proton transport behaviors of Nafion/Ti_(3)C_(2)T_(x) composite membrane at the molecular level. The results shown that there were significant differences about the diffusion kinetics of water molecules and hydroxium ions in Nafion/Ti_(3)C_(2)T_(x) at low and high hydration levels in the nanoscale region.With the increase of water content, Ti_(3)C_(2)T_(x) in membrane was gradually surrounded by ambient water molecules to form a hydration layer, and forming a relatively continuous proton transport channel between Nafion polymer and Ti_(3)C_(2)T_(x) monomer. The continuous proton transport channel could increase the number of binding sites of proton and thus achieving high proton conductivity and high mobility of water molecules at higher hydration level. The current work can provide a theoretical guidance for designing new type of Nafion composite membranes.展开更多
The numerical simulation is based on the authors' high-order models with a dissipative term for nonlinear and dispersive wave in water of varying depth. Corresponding finite-difference equations and general condit...The numerical simulation is based on the authors' high-order models with a dissipative term for nonlinear and dispersive wave in water of varying depth. Corresponding finite-difference equations and general conditions for open and fixed natural boundaries with an arbitrary reflection coefficient and phase shift are also given in this paper. The systematical tests of numerical simulation show that the theoretical models, the finite-difference algorithms and the boundary conditions can give good calculation results for the wave propagating in shallow and deep water with an arbitrary slope varying from gentle to steep.展开更多
In this paper, the three-dimensional variational data assimilation scheme (3DVAR) in the mesoscale model version 5 (MM5) of the US Pennsylvania State University/National Center for Atmospheric Research is used to stud...In this paper, the three-dimensional variational data assimilation scheme (3DVAR) in the mesoscale model version 5 (MM5) of the US Pennsylvania State University/National Center for Atmospheric Research is used to study the effect of assimilating the sea-wind data from QuikSCAT on the prediction of typhoon track and intensity. The case of Typhoon Dujuan (2003) is first tested and the results show appreciable improvements. Twelve other cases in 2003 are then evaluated. The assimilation of the QuikSCAT data produces significant impacts on the structure of Dujuan in terms of the horizontal and vertical winds, sea-level pressure and temperature at the initial time. With the assimilation, the 24-h (48-h) track prediction of 11 (10) out of the 12 typhoons is improved. The 24-h (48-h) prediction of typhoon intensity is also improved in 10 (9) of the 12 cases. These experiments therefore demonstrate that assimilation of the QuikSCAT sea-wind data can increase the accuracy of typhoon track and intensity predictions through modification of the initial fields associated with the typhoon.展开更多
In order to verify the validity of finite element numerical simulation method for asphalt mixture, which consists of aggregates, mastic （where mastic is a kind of fine mixture composed of asphalt binder mixed with fi...In order to verify the validity of finite element numerical simulation method for asphalt mixture, which consists of aggregates, mastic （where mastic is a kind of fine mixture composed of asphalt binder mixed with fines and fine aggregates） and air voids, based on three-dimensional （3D） heterogeneous specimen, X-ray computerized tomography （X-ray CT） was used to scan the asphalt specimens to obtain the real internal microstrnctures of asphalt mixture. CT images were reconstructed to build up 3D digital specimen, and the viscoelastic properties of mastic were described with Burgers model The uniaxial creep numerical simulations of three different levels of aggregate gradation were conducted. The simulation results agree well with the experimental results. The numerical simulation of asphalt mixture incorporated with real 3D microstructure based on finite element method is a promising application to conduct research of asphalt concrete. Additionally, this method can increase the mechanistic understanding of global viscoelastic properties of asphalt mixtures by linking the real 3D microstructure.展开更多
基金financially supported by the Scientific Research Foundation of North China University of Technology(Grant Nos.11005136024XN147-87 and 110051360024XN151-86).
文摘With respect to oceanic fluid dynamics,certain models have appeared,e.g.,an extended time-dependent(3+1)-dimensional shallow water wave equation in an ocean or a river,which we investigate in this paper.Using symbolic computation,we find out,on one hand,a set of bilinear auto-Backlund transformations,which could connect certain solutions of that equation with other solutions of that equation itself,and on the other hand,a set of similarity reductions,which could go from that equation to a known ordinary differential equation.The results in this paper depend on all the oceanic variable coefficients in that equation.
基金supported by the BK21 FOUR funded by the Ministry of Education of Korea and National Research Foundation of Korea,a Korea Agency for Infrastructure Technology Advancement(KAIA)grant funded by the Ministry of Land,Infrastructure,and Transport(Grant 1615013176)IITP(Institute of Information&Coummunications Technology Planning&Evaluation)-ICAN(ICT Challenge and Advanced Network of HRD)grant funded by the Korea government(Ministry of Science and ICT)(RS-2024-00438411).
文摘This paper proposes a novel cargo loading algorithm applicable to automated conveyor-type loading systems.The algorithm offers improvements in computational efficiency and robustness by utilizing the concept of discrete derivatives and introducing logistics-related constraints.Optional consideration of the rotation of the cargoes was made to further enhance the optimality of the solutions,if possible to be physically implemented.Evaluation metrics were developed for accurate evaluation and enhancement of the algorithm’s ability to efficiently utilize the loading space and provide a high level of dynamic stability.Experimental results demonstrate the extensive robustness of the proposed algorithm to the diversity of cargoes present in Business-to-Consumer environments.This study contributes practical advancements in both cargo loading optimization and automation of the logistics industry,with potential applications in last-mile delivery services,warehousing,and supply chain management.
文摘Lie symmetry analysis is applied to a(3+1)-dimensional combined potential Kadomtsev-Petviashvili equation with B-type Kadomtsev-Petviashvili equation(pKP-BKP equation)and the corresponding similarity reduction equations are obtained with the different infinitesimal generators.Invariant solutions with arbitrary functions and constants for the(3+1)-dimensional pKP-BKP equation,including the lump solution,the periodic-lump solution,the two-kink solution,the breather solution and the lump-two-kink solution,have been studied analytically and graphically.
基金financially supported by the National Key R&D Program of China(Grant No.2023YFC3081200)the National Natural Science Foundation of China(Grant Nos.U21A20159 and 52179117).
文摘Contact detection is the most time-consuming stage in 3D discontinuous deformation analysis(3D-DDA)computation.Improving the efficiency of 3D-DDA is beneficial for its application in large-scale computing.In this study,aiming at the continuous-discontinuous simulation of 3D-DDA,a highly efficient contact detection strategy is proposed.Firstly,the global direct search(GDS)method is integrated into the 3D-DDA framework to address intricate contact scenarios.Subsequently,all geometric elements,including blocks,faces,edges,and vertices are divided into searchable and unsearchable parts.Contacts between unsearchable geometric elements would be directly inherited,while only searchable geometric elements are involved in contact detection.This strategy significantly reduces the number of geometric elements involved in contact detection,thereby markedly enhancing the computation efficiency.Several examples are adopted to demonstrate the accuracy and efficiency of the improved 3D-DDA method.The rock pillars with different mesh sizes are simulated under self-weight.The deformation and stress are consistent with the analytical results,and the smaller the mesh size,the higher the accuracy.The maximum speedup ratio is 38.46 for this case.Furthermore,the Brazilian splitting test on the discs with different flaws is conducted.The results show that the failure pattern of the samples is consistent with the results obtained by other methods and experiments,and the maximum speedup ratio is 266.73.Finally,a large-scale impact test is performed,and approximately 3.2 times enhanced efficiency is obtained.The proposed contact detection strategy significantly improves efficiency when the rock has not completely failed,which is more suitable for continuous-discontinuous simulation.
基金supported by the National Natural Science Foundation of China(No.52376021).
文摘Current transient analysis predominantly relies on zero-dimensional/one-dimensional tools,proficient at capturing aerothermodynamic variations across critical engine stations but insufficient for analyzing the internal flow field evolution during transients.Addressing this gap,the study presents an enhanced quasi-three dimensional(quasi-3D)transient simulation technique that integrates component volume effects,offering a significant leap from the preceding quasi-3D transient simulation method based on quasi-steady assumption.By embedding the component volume effects on density,momentum,and energy within the physical temporal dimension of the Navier-Stokes equations,the refined quasi-3D transient model achieves a closer representation of physical phenomena.Validation against a single-shaft turbofan engine’s experimental data confirms the model’s accuracy.Average errors for key performance indicators,including shaft speed,thrust,mass flow rate,and critical component exit temperature and pressure,remain below 0.41%,5.69%,2.55%,3.18%and 0.67%,respectively.Crucially,the model exposes a discernible temporal lag in the compressor outlet pressure and temperature response due to volume effects—previously unquantified in quasi-3D transient simulations.And further exploration of the meridional flow field emphasizes the consequential role of volumes in transient flow field evolution.Incorporating volume effects within quasi-3D transient simulations enhances engine modeling and is pivotal for precise transient analysis in engine design and optimization.
基金Project(10964004)supported by the National Natural Science Foundation of ChinaProject(20070231001)supported by Research Fund for the Doctoral Program of Higher Education of China+1 种基金Project(096RJZA104)supported by the Natural Science Foundation of Gansu Province,ChinaProject(SB14200801)supported by the Doctoral Fund of Lanzhou University of Technology
文摘A accelerated arithmetic algorithm of the dynamic computing regions was designed,and 3-dimensional numerical simulation of isothermal solidification for a binary alloy was implemented.The dendritic growth and the recalescence of Ni-Cu binary alloy during the solidification at different cooling rates were investigated.The effects of cooling rate on dendritic patterns and microsegregation patterns were studied.The computed results indicate that,with the increment of the cooling rate,the dendritic growth velocity increases,both the main branch and side-branches become slender,the secondary dendrite arm spacing becomes smaller,the inadequate solute diffusion in solid aggravates,and the severity of microsegregation ahead of interface aggravates.At a higher cooling rate,the binary alloy presents recalescence;while the cooling rate is small,no recalescence occurs.
基金supported by the National Natural Science Foundation of China(No.50777005)the Young Foundation of University of Electronic Science and Technology of China(No.JX0832)
文摘A new quasi-three-dimensional （quasi-3D） numeric simulation method for a high-voltage level-shifting circuit structure is proposed. The performances of the 3D structure are analyzed by combining some 2D device structures; the 2D devices are in two planes perpendicular to each other and to the surface of the semiconductor. In comparison with Davinci, the full 3D device simulation tool, the quasi-3D simulation method can give results for the potential and current distribution of the 3D high-voltage level-shifting circuit structure with appropriate accuracy and the total CPU time for simulation is significantly reduced. The quasi-3D simulation technique can be used in many cases with advantages such as saving computing time, making no demands on the high-end computer terminals, and being easy to operate.
文摘Upon the conservation of mass, momentum and energy, volume fraction and surface penetrative rate were employed to modify the conservative equations to simulate the effect of blockages on fluid flows and heat transfer. These equations were solved numerically with the finite differential method and the primitive variable approach. This method uses staggered grid and pressure correction schemes. A computer code FASTOR3D integrated the aforementioned algorithm. The preliminary results have been compared with conventional benchmark solutions. With auxiliary software DV, the numerical results were visualized in colorful images to demonstrate the variation of flow patterns and temperature profiles during the transient process. The results of the simulation code for the fluid flows and heat transfer in the sodium pool of a fast breeder reactor are acceptable.
基金funded by the Shanghai Natural Science Foundation of Shanghai Municipal Science and Technology Commission(20ZR1400600)the Fundamental Research Funds for the Central Universities(2232023G-06)through collaborative research with the Advanced Fibrous Materials Lab(AFML)at the University of British Columbia.
文摘Three-dimensional(3D)braided composites have significant potential for use in engineering structural materials.However,conventional 3D braiding machines are insufficient for designing composites with complex geometries.This paper proposes a programmable design methodology for 3D rotary braiding machines using circle-cutting and combination strategies.By introducing varying numbers of incisions on the circle,a diverse range of horn gears can be designed.Different combinations of these cut-circles allow the horn gears to be assembled into various 3D rotary braiders.The parametric equation for the braider plate is derived,showing that a combination strategy involving two cut-circles is feasible for braider design,whereas integrating three cut-circles simultaneously is impossible for a single machine.The construction of an automatic 6-3 type 3D braiding machine demonstrates the effectiveness of the proposed design strategy.This flexible braider design approach provides a practical solution for producing 3D braided composites with complex geometries.
基金Supported by the National Natural Science Foundation of China(Nos.42250710152,42192562)the Southern Laboratory of Ocean Science and Engineering(Guangdong Zhuhai)(No.SML 2020 SP 007)。
文摘Convectively unstable processes caused by dense water subsidence are common occurrences in high-latitude oceanic regions,and significantly modulate mass and heat transport and mixing processes in the ocean.An idealized numerical experiment using the large eddy simulation method was conducted to analyze the three-dimensional flow field structure and the mechanism for dense water subsidence.Specifically,a negative salt flux is set at the sea surface,in which salt flux enters the sea surface to simulate the icing and salting-out phenomena that occur at high latitudes.Results show that the mean-state 3D flow field of dense water subsidence exhibits a hollow conical distribution.The horizontal flow field is characterized by a cyclonic vortex that driven primarily by the pressure gradient and influenced by the Coriolis effect.Moreover,the inverse vertical pressure gradient generated by this vortex inhibits the sinking of the plume,leading to its off-axis deflection and the development of an anticyclonic precession.In addition,the impact of rotation on the structure of a sinking plume within a stratified environment is discussed.Both horizontal vortex intensity and cone angle of the hollow cone flow field are increased with increasing rotation rate,resulting in a decrease in the plume’s maximum sinking depth.Variances in rotation direction cause the horizontal vortex and sinking plumes of dense water in the northern and southern hemispheres to rotate in opposite directions.
基金supported in part by China Scholarship Council under Grant 202208200010。
文摘Wire arc additive manufacturing(WAAM)has emerged as a promising technique for producing large-scale metal components,favoured by high deposition rates,flexibility and low cost.Despite its potential,the complexity of WAAM processes,which involves intricate thermal dynamics,phase transitions,and metallurgical,mechanical,and chemical interactions,presents considerable challenges in final product qualities.Simulation technologies in WAAM have proven invaluable,providing accurate predictions in key areas such as material properties,defect identification,deposit morphology,and residual stress.These predictions play a critical role in optimising manufacturing strategies for the final product.This paper provides a comprehensive review of the simulation techniques applied in WAAM,tracing developments from 2013 to 2023.Initially,it analyses the current challenges faced by simulation methods in three main areas.Subsequently,the review explores the current modelling approaches and the applications of these simulations.Following this,the paper discusses the present state of WAAM simulation,identifying specific issues inherent to WAAM simulation itself.Finally,through a thorough review of existing literature and related analysis,the paper offers future perspectives on potential advancements in WAAM simulation strategies.
基金supported by the National Natural Science Foundation of China(Grant Nos.52404155 and 52304111)State Key Laboratory for Geomechanics and Deep Underground Engineering,China University of Mining&Technology,Beijing(Grant No.XD2024006).
文摘Drilling and blasting,characterized by their efficiency,ubiquity,and cost-effectiveness,have emerged as predominant techniques in rock excavation;however,they are accompanied by enormous destructive power.Accurately controlling the blasting energy and achieving the directional fracture of a rock mass have become common problems in the field.A two-dimensional blasting(2D blasting)technique was proposed that utilizes the characteristic that the tensile strength of a rock mass is significantly lower than its compressive strength.After blasting,only a 2D crack surface is generated along the predetermined direction,eliminating the damage to the reserved rock mass caused by conventional blasting.However,the interior of a natural rock mass is a"black box",and the process of crack propagation is difficult to capture,resulting in an unclear 2D blasting mechanism.To this end,a single-hole polymethyl methacrylate(PMMA)test piece was used to conduct a 2D blasting experiment with the help of a high-speed camera to capture the dynamic crack propagation process and the digital image correlation(DIC)method to analyze the evolution law of surface strain on the test piece.On this basis,a three-dimensional(3D)finite element model was established based on the progressive failure theory to simulate the stress,strain,damage,and displacement evolution process of the model under 2D blasting.The simulation results were consistent with the experimental results.The research results reveal the 2D blasting mechanism and provide theoretical support for the application of 2D blasting technology in the field of rock excavation.
基金supported by the National Natural Science Foundation of China(No.12305239)Scientific Research Foundation of Chongqing University of Technology(No.2023ZDZ053)National Key R&D Program of China(No.2019YFE03010001).
文摘To further research on high-parameter plasma,we plan to develop a two-dimensional hard X-ray(HXR)imaging system at the HL-3 tokamak to measure HXRs with energies ranging from 20 to 300 keV.The application of an array-structured detector ensures that this system can measure HXR-radiation spectra from the entire plasma cross section.Therefore,it is suitable for the study of fast-electron physics,such as radio-frequency wave current drives,fast electrons driving instabilities,and plasma disruptions in fusion research.In this study,we develop a simulation for calculating fast-electron bremsstrahlung in the HL-3 tokamak based on the Monte Carlo simulation code Geant4,in which the plasma geometry and forward scattering of fast-electron bremsstrahlung are considered.The preliminary calculation results indicate that the HXR energy deposi-tion on the detector is symmetrically distributed,even though the plasma distribution is asymmetric owing to the toroidal effect.These simulation results are helpful in constructing the relationship between the energy deposition on the detector and parameter distribution on the plasma cross section during HL-3 experiments.This is beneficial for the reconstruction of the fast-electron-distribution function and for optimizing the design of the HXR-imaging system.
文摘Lie group analysis method is applied to the extended(3+1)-dimensional Kadomtsev–Petviashvili–Boussinesq equation and the corresponding similarity reduction equations are obtained with various infinitesimal generators.By selecting suitable arbitrary functions in the similarity reduction solutions,we obtain abundant invariant solutions,including the trigonometric solution,the kink-lump interaction solution,the interaction solution between lump wave and triangular periodic wave,the two-kink solution,the lump solution,the interaction between a lump and two-kink and the periodic lump solution in different planes.These exact solutions are also given graphically to show the detailed structures of this high dimensional integrable system.
基金supported by the National Natural Science Foundation of China(Nos.12101572,12371256)2023 Shanxi Province Graduate Innovation Project(No.2023KY614)the 19th Graduate Science and Technology Project of North University of China(No.20231943)。
文摘In this paper,we mainly focus on proving the existence of lump solutions to a generalized(3+1)-dimensional nonlinear differential equation.Hirota’s bilinear method and a quadratic function method are employed to derive the lump solutions localized in the whole plane for a(3+1)-dimensional nonlinear differential equation.Three examples of such a nonlinear equation are presented to investigate the exact expressions of the lump solutions.Moreover,the 3d plots and corresponding density plots of the solutions are given to show the space structures of the lump waves.In addition,the breath-wave solutions and several interaction solutions of the(3+1)-dimensional nonlinear differential equation are obtained and their dynamics are analyzed.
基金the National Natural Science Foundation of China(Grant Nos.U23A20565,52301194,and 52101178)the Shanghai Science and Technology Commission(Grant No.22511100400)+1 种基金the startup funding from Shanghai Jiao Tong University(Grant No.WH220405009)Innovation Program of Shanghai Municipal Education Commission(Grant No.2023ZKZD15)for providing funding support for this research。
文摘Radiative cooling materials have gained prominence as a zero-energy solution for mitigating global warming.However,a comprehensive understanding of the atomic-scale optical properties and macroscopic optical performance of radiative cooling materials remains elusive,limiting insight into the underlying physics of their optical response and cooling efficacy.La_(2)O_(3)and HfO_(2),which represent rare earth and third/fourth subgroup inorganic oxides,respectively,show promise for radiative cooling applications.In this study,we used multiscale simulations to investigate the optical properties of La_(2)O_(3)and HfO_(2)across a broad spectrum.First-principles calculations revealed their dielectric functions and intrinsic refractive indices,and the results indicated that the slightly smaller bandgap of La_(2)O_(3)compared to HfO_(2)induces a higher refractive index in the solar band.Additionally,three-phonon scattering was found to provide more accurate infrared optical properties than two-phonon scattering,which enhanced the emissivity in the sky window.Monte Carlo simulations were also used to determine the macroscopic optical properties of La_(2)O_(3)and HfO_(2)coatings.Based on the simulated results,we identified that the particle size and particle volume fraction play a dominant role in the optical properties.Our findings underscore the potential of La_(2)O_(3)and HfO_(2)nanocomposites for environment-friendly cooling and offer a new approach for high-throughput screening of optical materials through multiscale simulations.
基金financially supported by the National Key R&D Program of China (Nos.2020YFB1505500 and 2020YFB1505503)。
文摘The perfluorosulfonic acid(PFSA) membrane doped with two-dimensional conductive filler Ti_(3)C_(2)T_(x) is a fuel cell proton exchange membrane with high application potential. Experimental studies showed that the proton conductivity of Nafion/Ti_(3)C_(2)T_(x) composite membrane is improved significantly compared with that in pure Nafion. However, the microscopic mechanism of doping on the enhancement of membrane performance is remain unclear now. In this work, molecular dynamics simulation was used to investigate the microscopic morphology and proton transport behaviors of Nafion/Ti_(3)C_(2)T_(x) composite membrane at the molecular level. The results shown that there were significant differences about the diffusion kinetics of water molecules and hydroxium ions in Nafion/Ti_(3)C_(2)T_(x) at low and high hydration levels in the nanoscale region.With the increase of water content, Ti_(3)C_(2)T_(x) in membrane was gradually surrounded by ambient water molecules to form a hydration layer, and forming a relatively continuous proton transport channel between Nafion polymer and Ti_(3)C_(2)T_(x) monomer. The continuous proton transport channel could increase the number of binding sites of proton and thus achieving high proton conductivity and high mobility of water molecules at higher hydration level. The current work can provide a theoretical guidance for designing new type of Nafion composite membranes.
文摘The numerical simulation is based on the authors' high-order models with a dissipative term for nonlinear and dispersive wave in water of varying depth. Corresponding finite-difference equations and general conditions for open and fixed natural boundaries with an arbitrary reflection coefficient and phase shift are also given in this paper. The systematical tests of numerical simulation show that the theoretical models, the finite-difference algorithms and the boundary conditions can give good calculation results for the wave propagating in shallow and deep water with an arbitrary slope varying from gentle to steep.
基金This research was supported by the National Natural Science Foundation of China under Grant No.40333025.
文摘In this paper, the three-dimensional variational data assimilation scheme (3DVAR) in the mesoscale model version 5 (MM5) of the US Pennsylvania State University/National Center for Atmospheric Research is used to study the effect of assimilating the sea-wind data from QuikSCAT on the prediction of typhoon track and intensity. The case of Typhoon Dujuan (2003) is first tested and the results show appreciable improvements. Twelve other cases in 2003 are then evaluated. The assimilation of the QuikSCAT data produces significant impacts on the structure of Dujuan in terms of the horizontal and vertical winds, sea-level pressure and temperature at the initial time. With the assimilation, the 24-h (48-h) track prediction of 11 (10) out of the 12 typhoons is improved. The 24-h (48-h) prediction of typhoon intensity is also improved in 10 (9) of the 12 cases. These experiments therefore demonstrate that assimilation of the QuikSCAT sea-wind data can increase the accuracy of typhoon track and intensity predictions through modification of the initial fields associated with the typhoon.
基金Project(51038004) supported by the National Natural Science Foundation of China
文摘In order to verify the validity of finite element numerical simulation method for asphalt mixture, which consists of aggregates, mastic （where mastic is a kind of fine mixture composed of asphalt binder mixed with fines and fine aggregates） and air voids, based on three-dimensional （3D） heterogeneous specimen, X-ray computerized tomography （X-ray CT） was used to scan the asphalt specimens to obtain the real internal microstrnctures of asphalt mixture. CT images were reconstructed to build up 3D digital specimen, and the viscoelastic properties of mastic were described with Burgers model The uniaxial creep numerical simulations of three different levels of aggregate gradation were conducted. The simulation results agree well with the experimental results. The numerical simulation of asphalt mixture incorporated with real 3D microstructure based on finite element method is a promising application to conduct research of asphalt concrete. Additionally, this method can increase the mechanistic understanding of global viscoelastic properties of asphalt mixtures by linking the real 3D microstructure.
