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.展开更多
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.展开更多
Tree rings provide long-term records of tree growth and climate changes,which makes them ideal benchmarks for forest modeling.Tree-ring information has greatly improved the reliability of 3-PG,which is one of the most...Tree rings provide long-term records of tree growth and climate changes,which makes them ideal benchmarks for forest modeling.Tree-ring information has greatly improved the reliability of 3-PG,which is one of the most commonly used process-based forest growth models.Here,we strengthen 3-PG's ability to simulate tree-ring width and stable carbon isotopes(δ^(13)C)by enhancing its descriptions of tree physiology.The major upgrade was adding a carbon storage pool for tree-ring formation using stored carbohydrates.We also incorporated previous modifications(replacing the age modifier with a height modifier)of 3-PG and tested their efficacy in improving tree-ring simulations.We ran the model based on two grand fir(Abies grandis)stands.The updated model greatly improved the simulations for both tree-ring widths andδ^(13)C.The results represent one of the best tree-ringδ^(13)C simulations,which accurately captured the amplitude in annual variations ofδ^(13)C.The correlations(R^(2))between simulations and observations reached 0.50 and 0.73 at two stands,respectively.The new model also greatly improved the simulations of raw tree-ring widths and detrended ring-width index(RWI).Because of better descriptions of tree physiology and more accurate simulations of tree rings than the previous model version,the updated 3-PG should provide more reliable simulations than previous 3-PG versions when tree-ring information is used as a benchmark in future studies.展开更多
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.展开更多
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.展开更多
Neutron well logging,using instruments equipped with neutron source and detectors(e.g.,^(3)He-tubes,Nal,BGO),plays a key role in lithological differentiation,porosity determination,and fluid property evaluation in the...Neutron well logging,using instruments equipped with neutron source and detectors(e.g.,^(3)He-tubes,Nal,BGO),plays a key role in lithological differentiation,porosity determination,and fluid property evaluation in the petroleum industry.The growing trend of multifu nctional neutron well logging,which enables simultaneous extraction of multiple reservoir characteristics,requiring high-performance detectors capable of withstanding high-temperature downhole conditions,limited space,and instrument vibrations,while also detecting multiple particle types.The Cs_(2)LiYCl_(6):Ce^(3+)(CLYC)elpasolite scintillator demonstrates excellent temperature resistance and detection efficiency,making it become a promising candidate for leading the development of the novel neutron-based double-particle logging technology.This study employed Monte Carlo simulations to generate equivalent gamma spectra and proposed a pulse shape discrimination simulation method based on theoretical analysis and probabilistic iteration.The performance of CLYC was compared to that of common detectors in terms of physical properties and detection efficiency.A double-particle pulsed neutron detection system for porosity determination was established,based on the count ratio of equivalent gamma rays from the range of 2.95-3.42 MeVee energy bins.Results showed that CLYC can effectively replace ^(3)He-tubes for porosity measurement,providing consistent responses.This study offers numerical simulation support for the design of future neutron well logging tools and the application of double-particle detectors in logging systems.展开更多
The Secondary Air System(SAS)plays an important role in the safe operation and performance of aeroengines.The traditional 1D-3D coupling method loses information when used for secondary air systems,which affects the c...The Secondary Air System(SAS)plays an important role in the safe operation and performance of aeroengines.The traditional 1D-3D coupling method loses information when used for secondary air systems,which affects the calculation accuracy.In this paper,a Cross-dimensional Data Transmission method(CDT)from 3D to 1D is proposed by introducing flow field uniformity into the data transmission.First,a uniformity index was established to quantify the flow field parameter distribution characteristics,and a uniformity index prediction model based on the locally weighted regression method(Lowess)was established to quickly obtain the flow field information.Then,an information selection criterion in 3D to 1D data transmission was established based on the Spearman rank correlation coefficient between the uniformity index and the accuracy of coupling calculation,and the calculation method was automatically determined according to the established criterion.Finally,a modified function was obtained by fitting the ratio of the 3D mass-average parameters to the analytical solution,which are then used to modify the selected parameters at the 1D-3D interface.Taking a typical disk cavity air system as an example,the results show that the calculation accuracy of the CDT method is greatly improved by a relative 53.88%compared with the traditional 1D-3D coupling method.Furthermore,the CDT method achieves a speedup of 2 to 3 orders of magnitude compared to the 3D calculation.展开更多
The fabrication process dependent effects on single event effects （SEEs） are investigated in a commercial silicon- germanium heterojunction bipolar transistor （SiGe HBT） using three-dimensional （3D） TCAD simulat...The fabrication process dependent effects on single event effects （SEEs） are investigated in a commercial silicon- germanium heterojunction bipolar transistor （SiGe HBT） using three-dimensional （3D） TCAD simulations. The influences of device structure and doping concentration on SEEs are discussed via analysis of current transient and charge collection induced by ions strike. The results show that the SEEs representation of current transient is different from representation of the charge collection for the same process parameters. To be specific, the area of C/S junction is the key parameter that affects charge collection of SEE. Both current transient and charge collection are dependent on the doping of collector and substrate. The base doping slightly influences transient currents of base, emitter, and collector terminals. However, the SEEs of SiGe HBT are hardly affected by the doping of epitaxial base and the content of Ge.展开更多
To provide a test platform for Electronic Warfare (EW) system, it is needed to simulate the radar received Intermediate Frequency (IF) signals and radar system functions.This letter gives a description of a radar syst...To provide a test platform for Electronic Warfare (EW) system, it is needed to simulate the radar received Intermediate Frequency (IF) signals and radar system functions.This letter gives a description of a radar system simulation software developed for frequencyphase scanning three-dimensional (3-D) radar. Experimental results prove that the software could be used for system evaluation and for training purposes as an attractive alternative to real EW system.展开更多
Traditional clothing design models based on adaptive meshes cannot reflect.To solve this problem,a clothing simulation design model based on 3D image analysis technology is established.The model uses feature extractio...Traditional clothing design models based on adaptive meshes cannot reflect.To solve this problem,a clothing simulation design model based on 3D image analysis technology is established.The model uses feature extraction and description of image evaluation parameters,and establishes the mapping relationship between image features and simulation results by using the optimal parameter values,thereby obtaining a three-dimensional image simulation analysis environment.On the basis of this model,by obtaining the response results of clothing collision detection and the results of local adaptive processing of clothing meshes,the cutting form and actual cutting effect of clothing are determined to construct a design model.The simulation results show that compared with traditional clothing design models,clothing simulation design based on 3D image analysis technology has a better effect,with the definition of fabric folds increasing by 40%.More striking contrast between light and dark,the resolution increasing by 30%,and clothing details getting a more real manifestation.展开更多
A program MVFT3D of large-eddy simulation is developed and performed to solve the multi compressible Navier- Stokes equations. The SGS dissipation and molecular viscosity dissipation have been analyzed, and the former...A program MVFT3D of large-eddy simulation is developed and performed to solve the multi compressible Navier- Stokes equations. The SGS dissipation and molecular viscosity dissipation have been analyzed, and the former is much larger than the later. Our test shows that the SGS dissipation of Vreman model is smaller than the Smagorinsky model. We mainly simulate the experiment of fluid instability of shock-accelerated interface by Poggi in this paper. The decay of the turbulent kinetic energy before the first reflected shock wave–mixing zone interaction and its strong enhancement by re-shocks are presented in our numerical simulations. The computational mixing zone width under double re-shock agreement well with the experiment, and the decaying law of the turbulent kinetic energy is consistent with Mohamed and Larue’s investigation. Also, by using MVFT3D we give some simulation results of the inverse Chevron model from AWE. The numerical simulations presented in this paper allow us to characterize and better understand the Richtmyer-Meshkov instability induced turbulence, and the code MVFT3D is validated.展开更多
3D (three-dimensional) process simulation is currently one of the most challenging fields of research on debris flow.Large scale terrain rendering is the most basic task of 3D scenery construction in debris flow simul...3D (three-dimensional) process simulation is currently one of the most challenging fields of research on debris flow.Large scale terrain rendering is the most basic task of 3D scenery construction in debris flow simulation.As the major trigger for debris flow,rainfall will substantially enhance the realistic sense.Terrain and rainfall rendering in 3D debris flow simulations poses great challenges for numerical computation and graphical processing capability.In this paper,we propose to integrate GPU technology,LoD algorithms,and particle systems to realize 3D scenery modeling and rendering.The real-time LoD-based terrain modeling and rendering algorithm is presented first,and then a particle system-based rainfall scenery rendering method is implemented.Experimental results demonstrate that the 3D scenery rendered with the proposed approach exhibits sound performance and fair visual effects,which lays a solid foundation for the whole process simulation of debris flow disasters.展开更多
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.展开更多
Large-scale,fine,and efficient numerical simulation of a geothermal field plays an important role in geothermal energy development.Confronted with the problem of large computation and high storage requirements for com...Large-scale,fine,and efficient numerical simulation of a geothermal field plays an important role in geothermal energy development.Confronted with the problem of large computation and high storage requirements for complex underground models in a three-dimensional(3-D)numerical simulation of a geothermal fi eld,a mixed space-wavenumber domain 3-D numerical simulation algorithm is proposed in this paper.According to the superposition principle of temperature field,the geothermal field is decomposed into background and abnormal temperature fi elds for calculation.The uniform layered model is used to solve the background field.When the abnormal field is solved,the horizontal two-dimensional(2-D)Fourier transform is used to transform the 3-D diff erential equation satisfi ed by an abnormal field into a series of one-dimensional ordinary differential equations with diff erent wavenumbers,which greatly reduces the calculation and storage.The unit division of an ordinary diff erential equation is fl exible,and the calculation amount is small.The algorithm fully takes advantage of the effi ciency of the Fourier transform and the quickness of the catch-up method to solve linear equations with a fixed bandwidth,which effectively improves the computational efficiency.Compared with the COMSOL Multiphysics professional simulation finite element software,the time consumption and memory requirements of the algorithm proposed in this paper are reduced by multiple orders of magnitude in terms of ensuring accuracy and the same mesh division.The more the number of calculated nodes is,the more obvious is the advantage.We design models to study the thermal conductivity,heat fl ux boundary,regional tectonic morphology,and topographic relief of the geothermal fi eld distribution.A 3-D geophysical model is developed based on topographic elevation data,geothermal geology,and geophysical exploration data in the Qiabuqia area of Gonghe Basin,Qinghai Province,China.Numerical simulation of the geothermal fi eld in this area is realized,which shows that the algorithm is suitable for precise and effi cient simulation of an arbitrary complex terrain and geological conditions.展开更多
基金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.
基金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.
基金supported by the National Natural Science Foundation of China(Nos.42271048,42430503,and 31971492).
文摘Tree rings provide long-term records of tree growth and climate changes,which makes them ideal benchmarks for forest modeling.Tree-ring information has greatly improved the reliability of 3-PG,which is one of the most commonly used process-based forest growth models.Here,we strengthen 3-PG's ability to simulate tree-ring width and stable carbon isotopes(δ^(13)C)by enhancing its descriptions of tree physiology.The major upgrade was adding a carbon storage pool for tree-ring formation using stored carbohydrates.We also incorporated previous modifications(replacing the age modifier with a height modifier)of 3-PG and tested their efficacy in improving tree-ring simulations.We ran the model based on two grand fir(Abies grandis)stands.The updated model greatly improved the simulations for both tree-ring widths andδ^(13)C.The results represent one of the best tree-ringδ^(13)C simulations,which accurately captured the amplitude in annual variations ofδ^(13)C.The correlations(R^(2))between simulations and observations reached 0.50 and 0.73 at two stands,respectively.The new model also greatly improved the simulations of raw tree-ring widths and detrended ring-width index(RWI).Because of better descriptions of tree physiology and more accurate simulations of tree rings than the previous model version,the updated 3-PG should provide more reliable simulations than previous 3-PG versions when tree-ring information is used as a benchmark in future studies.
基金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.
文摘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.
基金the support of the National Natural Science Foundation of China(42174147,42474155)the Scientific and Technological Innovation Projects of Laoshan Laboratory(LSKJ20220347)。
文摘Neutron well logging,using instruments equipped with neutron source and detectors(e.g.,^(3)He-tubes,Nal,BGO),plays a key role in lithological differentiation,porosity determination,and fluid property evaluation in the petroleum industry.The growing trend of multifu nctional neutron well logging,which enables simultaneous extraction of multiple reservoir characteristics,requiring high-performance detectors capable of withstanding high-temperature downhole conditions,limited space,and instrument vibrations,while also detecting multiple particle types.The Cs_(2)LiYCl_(6):Ce^(3+)(CLYC)elpasolite scintillator demonstrates excellent temperature resistance and detection efficiency,making it become a promising candidate for leading the development of the novel neutron-based double-particle logging technology.This study employed Monte Carlo simulations to generate equivalent gamma spectra and proposed a pulse shape discrimination simulation method based on theoretical analysis and probabilistic iteration.The performance of CLYC was compared to that of common detectors in terms of physical properties and detection efficiency.A double-particle pulsed neutron detection system for porosity determination was established,based on the count ratio of equivalent gamma rays from the range of 2.95-3.42 MeVee energy bins.Results showed that CLYC can effectively replace ^(3)He-tubes for porosity measurement,providing consistent responses.This study offers numerical simulation support for the design of future neutron well logging tools and the application of double-particle detectors in logging systems.
基金supported by the National Science and Technology Major Project,China(No.2017-III-0010-0036).
文摘The Secondary Air System(SAS)plays an important role in the safe operation and performance of aeroengines.The traditional 1D-3D coupling method loses information when used for secondary air systems,which affects the calculation accuracy.In this paper,a Cross-dimensional Data Transmission method(CDT)from 3D to 1D is proposed by introducing flow field uniformity into the data transmission.First,a uniformity index was established to quantify the flow field parameter distribution characteristics,and a uniformity index prediction model based on the locally weighted regression method(Lowess)was established to quickly obtain the flow field information.Then,an information selection criterion in 3D to 1D data transmission was established based on the Spearman rank correlation coefficient between the uniformity index and the accuracy of coupling calculation,and the calculation method was automatically determined according to the established criterion.Finally,a modified function was obtained by fitting the ratio of the 3D mass-average parameters to the analytical solution,which are then used to modify the selected parameters at the 1D-3D interface.Taking a typical disk cavity air system as an example,the results show that the calculation accuracy of the CDT method is greatly improved by a relative 53.88%compared with the traditional 1D-3D coupling method.Furthermore,the CDT method achieves a speedup of 2 to 3 orders of magnitude compared to the 3D calculation.
基金supported by the National Natural Science Foundation of China(Grant Nos.61274106,11175138,and 61601352)
文摘The fabrication process dependent effects on single event effects （SEEs） are investigated in a commercial silicon- germanium heterojunction bipolar transistor （SiGe HBT） using three-dimensional （3D） TCAD simulations. The influences of device structure and doping concentration on SEEs are discussed via analysis of current transient and charge collection induced by ions strike. The results show that the SEEs representation of current transient is different from representation of the charge collection for the same process parameters. To be specific, the area of C/S junction is the key parameter that affects charge collection of SEE. Both current transient and charge collection are dependent on the doping of collector and substrate. The base doping slightly influences transient currents of base, emitter, and collector terminals. However, the SEEs of SiGe HBT are hardly affected by the doping of epitaxial base and the content of Ge.
文摘To provide a test platform for Electronic Warfare (EW) system, it is needed to simulate the radar received Intermediate Frequency (IF) signals and radar system functions.This letter gives a description of a radar system simulation software developed for frequencyphase scanning three-dimensional (3-D) radar. Experimental results prove that the software could be used for system evaluation and for training purposes as an attractive alternative to real EW system.
文摘Traditional clothing design models based on adaptive meshes cannot reflect.To solve this problem,a clothing simulation design model based on 3D image analysis technology is established.The model uses feature extraction and description of image evaluation parameters,and establishes the mapping relationship between image features and simulation results by using the optimal parameter values,thereby obtaining a three-dimensional image simulation analysis environment.On the basis of this model,by obtaining the response results of clothing collision detection and the results of local adaptive processing of clothing meshes,the cutting form and actual cutting effect of clothing are determined to construct a design model.The simulation results show that compared with traditional clothing design models,clothing simulation design based on 3D image analysis technology has a better effect,with the definition of fabric folds increasing by 40%.More striking contrast between light and dark,the resolution increasing by 30%,and clothing details getting a more real manifestation.
文摘A program MVFT3D of large-eddy simulation is developed and performed to solve the multi compressible Navier- Stokes equations. The SGS dissipation and molecular viscosity dissipation have been analyzed, and the former is much larger than the later. Our test shows that the SGS dissipation of Vreman model is smaller than the Smagorinsky model. We mainly simulate the experiment of fluid instability of shock-accelerated interface by Poggi in this paper. The decay of the turbulent kinetic energy before the first reflected shock wave–mixing zone interaction and its strong enhancement by re-shocks are presented in our numerical simulations. The computational mixing zone width under double re-shock agreement well with the experiment, and the decaying law of the turbulent kinetic energy is consistent with Mohamed and Larue’s investigation. Also, by using MVFT3D we give some simulation results of the inverse Chevron model from AWE. The numerical simulations presented in this paper allow us to characterize and better understand the Richtmyer-Meshkov instability induced turbulence, and the code MVFT3D is validated.
基金the National Basic Research Program of China (973 Program) (No. 2010CB731504)the One Hundred Person Program of Chinese Academy of Sciences (No.110900K242)+1 种基金the West Light Foundation of Chinese Academy of Sciences (No. Y1R2030030)the Youth Talent Team Program of IMHE,Chinese Academy of Sciences (No.SDSQB-2010-03)
文摘3D (three-dimensional) process simulation is currently one of the most challenging fields of research on debris flow.Large scale terrain rendering is the most basic task of 3D scenery construction in debris flow simulation.As the major trigger for debris flow,rainfall will substantially enhance the realistic sense.Terrain and rainfall rendering in 3D debris flow simulations poses great challenges for numerical computation and graphical processing capability.In this paper,we propose to integrate GPU technology,LoD algorithms,and particle systems to realize 3D scenery modeling and rendering.The real-time LoD-based terrain modeling and rendering algorithm is presented first,and then a particle system-based rainfall scenery rendering method is implemented.Experimental results demonstrate that the 3D scenery rendered with the proposed approach exhibits sound performance and fair visual effects,which lays a solid foundation for the whole process simulation of debris flow disasters.
基金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.
基金supported by National Natural Science Foundation of China(No.41574127,42174080)Innovation research team project of Guangxi Natural Science Foundation(No.GXNSFGA380004)Central South University independent exploration and innovation project for Postgraduates(Nos.2021zzts0831,2021zzts0271)
文摘Large-scale,fine,and efficient numerical simulation of a geothermal field plays an important role in geothermal energy development.Confronted with the problem of large computation and high storage requirements for complex underground models in a three-dimensional(3-D)numerical simulation of a geothermal fi eld,a mixed space-wavenumber domain 3-D numerical simulation algorithm is proposed in this paper.According to the superposition principle of temperature field,the geothermal field is decomposed into background and abnormal temperature fi elds for calculation.The uniform layered model is used to solve the background field.When the abnormal field is solved,the horizontal two-dimensional(2-D)Fourier transform is used to transform the 3-D diff erential equation satisfi ed by an abnormal field into a series of one-dimensional ordinary differential equations with diff erent wavenumbers,which greatly reduces the calculation and storage.The unit division of an ordinary diff erential equation is fl exible,and the calculation amount is small.The algorithm fully takes advantage of the effi ciency of the Fourier transform and the quickness of the catch-up method to solve linear equations with a fixed bandwidth,which effectively improves the computational efficiency.Compared with the COMSOL Multiphysics professional simulation finite element software,the time consumption and memory requirements of the algorithm proposed in this paper are reduced by multiple orders of magnitude in terms of ensuring accuracy and the same mesh division.The more the number of calculated nodes is,the more obvious is the advantage.We design models to study the thermal conductivity,heat fl ux boundary,regional tectonic morphology,and topographic relief of the geothermal fi eld distribution.A 3-D geophysical model is developed based on topographic elevation data,geothermal geology,and geophysical exploration data in the Qiabuqia area of Gonghe Basin,Qinghai Province,China.Numerical simulation of the geothermal fi eld in this area is realized,which shows that the algorithm is suitable for precise and effi cient simulation of an arbitrary complex terrain and geological conditions.
