The three-dimensional （3D） processing maps considering strain based on the two-dimensional （2D） processing maps proposed by PRASAD can describe the distribution of the efficiency of power dissipation and flow inst...The three-dimensional （3D） processing maps considering strain based on the two-dimensional （2D） processing maps proposed by PRASAD can describe the distribution of the efficiency of power dissipation and flow instability regions at various temperatures, strain rates and strains, which exhibit intrinsic workability related to material itself. Finite element （FE） simulation can obtain the distribution of strain, strain rate, temperature and die filling status, which indicates state-of-stress （SOS） workability decided by die shape and different processing conditions. On the basis of this, a new material driven analysis method for hot deformation was put forward by the combination of FE simulation with 3D processing maps, which can demonstrate material workability of the entire hot deformation process including SOS workability and intrinsic workability. The hot forging process for hard-to-work metal magnesium alloy was studied, and the 3D thermomechanical FE simulation including 3D processing maps of complex hot forging spur bevel gear was first conducted. The hot forging experiments were carried out. The results show that the new method is reasonable and suitable to determine the aoorooriate nrocess narameters.展开更多
Simulations of photoresist etching,aerial image,exposure,and post-bake processes are integrated to obtain a photolithography process simulation for microelectromechanical system（MEMS） and integrated circuit（IC） fa...Simulations of photoresist etching,aerial image,exposure,and post-bake processes are integrated to obtain a photolithography process simulation for microelectromechanical system（MEMS） and integrated circuit（IC） fabrication based on three-dimensional （3D） cellular automata（CA）. The simulation results agree well with available experimental results. This indicates that the 3D dynamic CA model for the photoresist etching simulation and the 3D CA model for the post-bake simulation could be useful for the monolithic simulation of various lithography processes. This is determined to be useful for the device-sized fabrication process simulation of IC and MEMS.展开更多
Sequential-modular-based process flowsheeting software remains an indispensable tool for process design,control,and optimization.Yet,as the process industry advances in intelligent operation and maintenance,convention...Sequential-modular-based process flowsheeting software remains an indispensable tool for process design,control,and optimization.Yet,as the process industry advances in intelligent operation and maintenance,conventional sequential-modular-based process-simulation techniques present challenges regarding computationally intensive calculations and significant central processing unit(CPU)time requirements,particularly in large-scale design and optimization tasks.To address these challenges,this paper proposes a novel process-simulation parallel computing framework(PSPCF).This framework achieves layered parallelism in recycling processes at the unit operation level.Notably,PSPCF introduces a groundbreaking concept of formulating simulation problems as task graphs and utilizes Taskflow,an advanced task graph computing system,for hierarchical parallel scheduling and the execution of unit operation tasks.PSPCF also integrates an advanced work-stealing scheme to automatically balance thread resources with the demanding workload of unit operation tasks.For evaluation,both a simpler parallel column process and a more complex cracked gas separation process were simulated on a flowsheeting platform using PSPCF.The framework demonstrates significant time savings,achieving over 60%reduction in processing time for the simpler process and a 35%–40%speed-up for the more complex separation process.展开更多
In the ultra-deep strata of the Tarim Basin,the vertical growth process of strike-slip faults remains unclear,and the vertical distribution of fractured-cavity carbonate reservoirs is complex.This paper investigates t...In the ultra-deep strata of the Tarim Basin,the vertical growth process of strike-slip faults remains unclear,and the vertical distribution of fractured-cavity carbonate reservoirs is complex.This paper investigates the vertical growth process of strike-slip faults through field outcrop observations in the Keping area,interpretation of seismic data from the Fuman Oilfield,Tarim Basim,NW China,and structural physical simulation experiments.The results are obtained mainly in four aspects.First,field outcrops and ultra-deep seismic profiles indicate a three-layer structure within the strike-slip fault,consisting of fault core,fracture zone and primary rock.The fault core can be classified into three parts vertically:fracture-cavity unit,fault clay and breccia zone.The distribution of fracture-cavity units demonstrates a distinct pattern of vertical stratification,owing to the structural characteristics and growth process of the slip-strike fault.Second,the ultra-deep seismic profiles show multiple fracture-cavity units in the strike-slip fault zone.These units can be classified into four types:top fractured,middle connected,deep terminated,and intra-layer fractured.Third,structural physical simulation experiments and ultra-deep seismic data interpretation reveal that the strike-slip faults have evolved vertically in three stages:segmental rupture,vertical growth,and connection and extension.The particle image velocimetry detection demonstrates that the initial fracture of the fault zone occurred at the top or bottom and then evolved into cavities gradually along with the fault growth,accompanied by the emergence of new fractures in the middle part of the strata,which subsequently connected with the deep and shallow cavities to form a complete fault zone.Fourth,the ultra-deep carbonate strata primarily develop three types of fractured-cavity reservoirs:flower-shaped fracture,large and deep fault and staggered overlap.The first two types are larger in size with better reservoir conditions,suggesting a significant exploration potential.展开更多
There is little low-and-high frequency information on seismic data in seismic exploration,resulting in narrower bandwidth and lower seismic resolution.It considerably restricts the prediction accuracy of thin reservoi...There is little low-and-high frequency information on seismic data in seismic exploration,resulting in narrower bandwidth and lower seismic resolution.It considerably restricts the prediction accuracy of thin reservoirs and thin interbeds.This study proposes a novel method to constrain improving seismic resolution in the time and frequency domain.The expected wavelet spectrum is used in the frequency domain to broaden the seismic spectrum range and increase the octave.In the time domain,the Frobenius vector regularization of the Hessian matrix is used to constrain the horizontal continuity of the seismic data.It eff ectively protects the signal-to-noise ratio of seismic data while the longitudinal seismic resolution is improved.This method is applied to actual post-stack seismic data and pre-stack gathers dividedly.Without abolishing the phase characteristics of the original seismic data,the time resolution is signifi cantly improved,and the structural features are clearer.Compared with the traditional spectral simulation and deconvolution methods,the frequency distribution is more reasonable,and seismic data has higher resolution.展开更多
Due to the restricted satellite payloads in LEO mega-constellation networks(LMCNs),remote sensing image analysis,online learning and other big data services desirably need onboard distributed processing(OBDP).In exist...Due to the restricted satellite payloads in LEO mega-constellation networks(LMCNs),remote sensing image analysis,online learning and other big data services desirably need onboard distributed processing(OBDP).In existing technologies,the efficiency of big data applications(BDAs)in distributed systems hinges on the stable-state and low-latency links between worker nodes.However,LMCNs with high-dynamic nodes and long-distance links can not provide the above conditions,which makes the performance of OBDP hard to be intuitively measured.To bridge this gap,a multidimensional simulation platform is indispensable that can simulate the network environment of LMCNs and put BDAs in it for performance testing.Using STK's APIs and parallel computing framework,we achieve real-time simulation for thousands of satellite nodes,which are mapped as application nodes through software defined network(SDN)and container technologies.We elaborate the architecture and mechanism of the simulation platform,and take the Starlink and Hadoop as realistic examples for simulations.The results indicate that LMCNs have dynamic end-to-end latency which fluctuates periodically with the constellation movement.Compared to ground data center networks(GDCNs),LMCNs deteriorate the computing and storage job throughput,which can be alleviated by the utilization of erasure codes and data flow scheduling of worker nodes.展开更多
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.展开更多
When castings become complicated and the demands for precision of numerical simulation become higher,the numerical data of casting numerical simulation become more massive.On a general personal computer,these massive ...When castings become complicated and the demands for precision of numerical simulation become higher,the numerical data of casting numerical simulation become more massive.On a general personal computer,these massive numerical data may probably exceed the capacity of available memory,resulting in failure of rendering.Based on the out-of-core technique,this paper proposes a method to effectively utilize external storage and reduce memory usage dramatically,so as to solve the problem of insufficient memory for massive data rendering on general personal computers.Based on this method,a new postprocessor is developed.It is capable to illustrate filling and solidification processes of casting,as well as thermal stess.The new post-processor also provides fast interaction to simulation results.Theoretical analysis as well as several practical examples prove that the memory usage and loading time of the post-processor are independent of the size of the relevant files,but the proportion of the number of cells on surface.Meanwhile,the speed of rendering and fetching of value from the mouse is appreciable,and the demands of real-time and interaction are satisfied.展开更多
Large eddy simulation （LES） using the Smagorinsky eddy viscosity model is added to the two-dimensional nine velocity components （D2Q9） lattice Boltzmann equation （LBE） with multi-relaxation-time （MRT） to simul...Large eddy simulation （LES） using the Smagorinsky eddy viscosity model is added to the two-dimensional nine velocity components （D2Q9） lattice Boltzmann equation （LBE） with multi-relaxation-time （MRT） to simulate incompressible turbulent cavity flows with the Reynolds numbers up to 1 × 10^7. To improve the computation efficiency of LBM on the numerical simulations of turbulent flows, the massively parallel computing power from a graphic processing unit （GPU） with a computing unified device architecture （CUDA） is introduced into the MRT-LBE-LES model. The model performs well, compared with the results from others, with an increase of 76 times in computation efficiency. It appears that the higher the Reynolds numbers is, the smaller the Smagorinsky constant should be, if the lattice number is fixed. Also, for a selected high Reynolds number and a selected proper Smagorinsky constant, there is a minimum requirement for the lattice number so that the Smagorinsky eddy viscosity will not be excessively large.展开更多
The absorbing process in isolating and coating process of α-olefin drag reducing polymer was studied by molecular dynamic simulation method, on basis of coating theory of α-olefin drag reducing polymer particles wit...The absorbing process in isolating and coating process of α-olefin drag reducing polymer was studied by molecular dynamic simulation method, on basis of coating theory of α-olefin drag reducing polymer particles with polyurethane as coating material. The distributions of sodium laurate, sodium dodeeyl sulfate, and sodium dodeeyl benzene sulfonate on the surface of α-olefin drag reducing polymer particles were almost the same, but the bending degrees of them were obviously different. The bending degree of SLA molecules was greater than those of the other two surfactant molecules. Simulation results of absorbing and accumulating structure showed that, though hydrophobie properties of surfactant molecules were almost the same, water density around long chain sulfonate sodium was bigger than that around alkyl sulfate sodium. This property goes against useful absorbing and accumulating on the surface of α-olefin drag reducing polymer particles; simulation results of interactions of different surfactant and multiple hydroxyl compounds on surface of particles showed that, interactions of different surfaetant and one kind of multiple hydroxyl compound were similar to those of one kind of surfaetant and different multiple hydroxyl compounds. These two contrast types of interactions also exhibited the differences of absorbing distribution and closing degrees to surface of particles. The sequence of closing degrees was derived from simulation; control step of addition polymerization interaction in coating process was absorbing mass transfer process, so the more closed to surface of particle the multiple hydroxyl compounds were, the easier interactions With isoeyanate were. Simulation results represented the compatibility relationship between surfactant and multiple hydroxyl compounds. The isolating and coating processes of α-olefin drag reducing polymer were further understood on molecule and atom level through above simulation research, and based on the simulation, a referenced theoretical basis was provided for practical optimal selection and experimental preparation of α-olefin drag reducing polymer particles suspension isolation agent.展开更多
This article presents a multiscale simulation approach starting at the molecular level for the adsorption process development. A grand canonical Monte Carlo method is used for the prediction of adsorption isotherms of...This article presents a multiscale simulation approach starting at the molecular level for the adsorption process development. A grand canonical Monte Carlo method is used for the prediction of adsorption isotherms of methanol on an activated carbon at the molecular level. The adsorption isotherms obtained in the linear region （or adsorption constant） are exploited as a model parameter required for the adsorption process simulation. The adsorption process model described by a set of partial differential equations （PDEs） is solved by using the conservation element and solution element method, which produces a fast and an accurate numerical solution to PDEs. The simulation results obtained from the adsorption constant estimated at the molecular level are in good agreement with the experimental results of the pulse response. The systematical multiscale simulation approach addressed in this study may be useful to accelerate the adsorption process development by reducing the number of experiments.展开更多
According to the characteristic of Team Software Process （TSP）, it adopts a hierarchy-based model combined discrete event model with system dynamics model. This model represents TSP as form of two levels, the inner ...According to the characteristic of Team Software Process （TSP）, it adopts a hierarchy-based model combined discrete event model with system dynamics model. This model represents TSP as form of two levels, the inner level embodies the continuity of the software process, the outer embodies the software development process by phases, and the structure and principle of the model is explained in detail, then formalization description of the model is offered. At last, an example is presented to demonstrate the simulation process and result. This model can simulate team software process from various angles, supervise and predict the software process. Also it can make the management of software development become more scientific and improve the quality of software.展开更多
Physical geography and human geography are the principal branches of the geographical sciences. Physical process simulation and human process simulation in geography are both quantitative methods used to recover past ...Physical geography and human geography are the principal branches of the geographical sciences. Physical process simulation and human process simulation in geography are both quantitative methods used to recover past events and even to forecast events based on precisely determined parameters. There are four differences between physical process simulation and human process simulation in geography, which we summarize with two specific cases, one of which is about a typhoon＇s development and its precipitation, and the other of which is regarding the evolution of three industrial structures in China. The differences focus on four aspects： the main factors of the research framework; the knowledge back-ground of the systematic analysis framework; the simulation data sources and quantitative method; and the core of the study object and the method of forecast application. As the human-land relationship is the key ideology of the man-land system, the relationship between the physical and human factors is becoming increasingly close at present. Physical process simulation and human process simulation in geography will exhibit crossing and blending in the future to reflect the various geographical phenomena better.展开更多
Since processing parameters have always been assumed to be stable in the current finite element numerical simulation of dieless drawing process, the simulation results for the product dimension tend to stabilize gradu...Since processing parameters have always been assumed to be stable in the current finite element numerical simulation of dieless drawing process, the simulation results for the product dimension tend to stabilize gradually. In fact, the dimension fluctuation exists in the forming process all the while. A mathematical model of Gauss distribution for processing parameters was employed and a finite element numerical model of dieless drawing process with non-steady processing parameters was established. Dieless drawing processing of Ni-Ti alloy wire was conducted for verifying the proposed model. The results indicated that the non-steady processing parameters model had higher simulation accuracy of the wire diameter than that given by the steady parameters model. Furthermore, the model could also be used to analyze the fluctuation characteristics in the whole dieless drawing process.展开更多
A molecular dynamics simulation study has been performed for the microstructure evolution in a liquid metal Ni system during crystallization process at two cooling rates by adopting the embedded atom method （EAM） mo...A molecular dynamics simulation study has been performed for the microstructure evolution in a liquid metal Ni system during crystallization process at two cooling rates by adopting the embedded atom method （EAM） model potential. The bond-type index method of Honeycutt-Andersen （HA） and a new cluster-type index method （CTIM-2） have been used to detect and analyse the microstructures in this system. It is demonstrated that the cooling rate plays a critical role in the microstructure evolution： below the crystallization temperature Tc, the effects of cooling rate are very remarkable and can be fully displayed. At different cooling rates of 2.0 × 10^13 K·s^-1 and 1.0 × 10^12 K·s^-1, two different kinds of crystal structures are obtained in the system. The first one is the coexistence of the hcp （expressed by （12 0 0 0 6 6） in CTIM-2） and the fcc （12 0 0 0 12 0） basic clusters consisting of 1421 and 1422 bond-types, and the hcp basic cluster becomes the dominant one with decreasing temperature, the second one is mainly the fcc （12 0 0 0 12 0） basic clusters consisting of 1421 bond-type, and their crystallization temperatures Tc would be 1073 and 1173 K, respectively.展开更多
The western coast of Hainan Island exhibits a savanna landscape. Many types of sand dunes, including transverse dune ridges, longitudinal dune ridges, elliptical dunes, coppice dunes, and climbing dunes, are widely di...The western coast of Hainan Island exhibits a savanna landscape. Many types of sand dunes, including transverse dune ridges, longitudinal dune ridges, elliptical dunes, coppice dunes, and climbing dunes, are widely distributed in the coastal zone. In winter, high-frequency and high-energy NE winds （dominant winds） are prevalent, with a resultant drift direction （RDD） of S35.6°W. In spring, low-frequency and low-energy SW secondary winds prevail, with a RDD of N25.1°E. Wind tunnel simulations revealed that the airflow over the dune surface is the main factor controlling the erosion and deposition patterns of dune surfaces and the morphological development of dunes. In the region＇s bidirectional wind environment, with two seasonally distinct energy levels, the airflow over the surface of elliptical dunes, barchan dunes, and transverse dune ridges will exhibit a transverse pattern, whereas the airflow over longitudinal dunes ridges exhibits a lateral pattern and that over climbing dunes exhibits a climbing-circumfluent pattern. These patterns represent different dynamic processes. The coastal dunes on the western coast of Hainan Island are influenced by factors such as onshore winds, sand sources, coastal slopes, rivers, and forest shelter belts. The source of the sand that supplements these dunes particularly influences the development pattern： when there is more sand, the pattern shows positive equilibrium deposition between dune ridges and dunes; otherwise, it shows negative equilibrium deposition. The presence or absence of forest shelter belts also influences deposition and dune development patterns and transformation of dune forms. Coastal dunes and inland desert dunes experience similar dynamic processes, but the former have more diversified shapes and more complex formation mechanisms.展开更多
Turbidity currents represent a major agent for sediment transport in lakes, seas and oceans. In particu-lar, they formulate the most significant clastic accumulations in the deep sea, which become many of the world＇s...Turbidity currents represent a major agent for sediment transport in lakes, seas and oceans. In particu-lar, they formulate the most significant clastic accumulations in the deep sea, which become many of the world＇s most important hydrocarbon reservoirs. Several boreholes in the Qiongdongnan Basin, the north-western South China Sea, have recently revealed turbidity current deposits as significant hydrocarbon res-ervoirs. However, there are some arguments for the potential provenances. To solve this problem, it is es-sential to delineate their sedimentary processes as well as to evaluate their qualities as reservoir. Numerical simulations have been developed rapidly over the last several years, offering insights into turbidity current behaviors, as geologically significant turbidity currents are difficult to directly investigate due to their large scale and often destructive nature. Combined with the interpretation of the turbidity system based on high-resolution 3D seismic data, the paleotophography is acquired via a back-stripping seismic profile integrated with a borehole, i.e., Well A, in the western Qiongdongnan Basin; then a numerical model is built on the basis of this back-stripped profile. After defining the various turbidity current initial boundary conditions, includ-ing grain size, velocity and sediment concentration, the structures and behaviors of turbidity currents are investigated via numerical simulation software ANSYS FLUENT. Finally, the simulated turbidity deposits are compared with the interpreted sedimentary bodies based on 3D seismic data and the potential provenances of the revealed turbidites by Well A are discussed in details. The simulation results indicate that a sedimen-tary body develops far away from its source with an average grain size of 0.1 mm, i.e., sand-size sediment. Taking into account the location and orientation of the simulated seismic line, the consistence between normal forward simulation results and the revealed cores in Well A indicates that the turbidites should have been transported from Vietnam instead of Hainan Island. This interpretation has also been verified by the planar maps of sedimentary systems based on integration of boreholes and seismic data. The identification of the turbidity provenance will benefit the evaluation of extensively distributed submarine fans for hydro-carbon exploration in the deepwater areas.展开更多
A new 18-lump kinetic model for naphtha catalytic reforming reactions is discussed. By developing this model as a user module, a whole industrial continuous catalytic reforming process is simulated on Aspen plus plat-...A new 18-lump kinetic model for naphtha catalytic reforming reactions is discussed. By developing this model as a user module, a whole industrial continuous catalytic reforming process is simulated on Aspen plus plat-form. The technique utilizes the strong databases, complete sets of modules, and flexible simulation tools of the Aspen plus system and retains the characteristics of the proposed kinetic model. The calculated results are in fair agreement with the actual operating data. Based on the model of the whole reforming process, the process is opti-mized and the optimization results are tested in the actual industrial unit for about two months. The test shows that the process profit increases about 1000yuan·h-1 averagely, which is close to the calculated result.展开更多
Numerical simulation of casting＇s mold filling process is the main and the most important aspect of the foundry CAE technology. But it is time-consuming; it may take dozens of hours or several days. While with the de...Numerical simulation of casting＇s mold filling process is the main and the most important aspect of the foundry CAE technology. But it is time-consuming; it may take dozens of hours or several days. While with the development of computer hardware, numerical simulation of casting＇ s mold filling process has made rapid progress. The simulation results, therefore, have become more and more practical. This study tries to find some clues of the computational time of mold filling process. Firstly, this paper introduces mathematic model and the basic route of numerical simulation of casting＇s mold filling process. Then the computational time of mold filling process has been carefully studied, and some new and useful results have been gained from the study of the computational time. Finally, this paper has given some real applications of numerical simulation of casting＇s mold filling process.展开更多
A 3-D mathematical model considering turbulence phenomena has been established based on a computational fluid dynamics technique, so called 3-D SOLA-VOF (Solution Algorithm-Volume of Fluid), to simulate the fluid flow...A 3-D mathematical model considering turbulence phenomena has been established based on a computational fluid dynamics technique, so called 3-D SOLA-VOF (Solution Algorithm-Volume of Fluid), to simulate the fluid flow of mold filling process of die casting. In addition, the mathematical model for simulating the heat transfer in die casting process has also been established. The computation program has been developed by the authors with the finite difference method (FDM) recently. As verification, the mold filling process of a S-shaped die casting has been simulated and the simulation results coincide with that of the benchmark test. Finally, as a practical application, the gating design of a motorcycle component was modified by the mold filling simulation and the dies design of another motorcycle component was optimized by the heat transfer simulation. All the optimized designs were verified by the production practice.展开更多
基金Project(2011ZX04014-051)supported by the Key Scientific and Technical Project of ChinaProjects(51375306,50905110)supported by the National Natural Science Foundation of China
文摘The three-dimensional （3D） processing maps considering strain based on the two-dimensional （2D） processing maps proposed by PRASAD can describe the distribution of the efficiency of power dissipation and flow instability regions at various temperatures, strain rates and strains, which exhibit intrinsic workability related to material itself. Finite element （FE） simulation can obtain the distribution of strain, strain rate, temperature and die filling status, which indicates state-of-stress （SOS） workability decided by die shape and different processing conditions. On the basis of this, a new material driven analysis method for hot deformation was put forward by the combination of FE simulation with 3D processing maps, which can demonstrate material workability of the entire hot deformation process including SOS workability and intrinsic workability. The hot forging process for hard-to-work metal magnesium alloy was studied, and the 3D thermomechanical FE simulation including 3D processing maps of complex hot forging spur bevel gear was first conducted. The hot forging experiments were carried out. The results show that the new method is reasonable and suitable to determine the aoorooriate nrocess narameters.
文摘Simulations of photoresist etching,aerial image,exposure,and post-bake processes are integrated to obtain a photolithography process simulation for microelectromechanical system（MEMS） and integrated circuit（IC） fabrication based on three-dimensional （3D） cellular automata（CA）. The simulation results agree well with available experimental results. This indicates that the 3D dynamic CA model for the photoresist etching simulation and the 3D CA model for the post-bake simulation could be useful for the monolithic simulation of various lithography processes. This is determined to be useful for the device-sized fabrication process simulation of IC and MEMS.
基金supported by the National Key Research and Development Program of China(2022YFB3305900)the National Natural Science Foundation of China(Key Program)(62136003)+2 种基金the National Natural Science Foundation of China(62394345)the Major Science and Technology Projects of Longmen Laboratory(LMZDXM202206)the Fundamental Research Funds for the Central Universities.
文摘Sequential-modular-based process flowsheeting software remains an indispensable tool for process design,control,and optimization.Yet,as the process industry advances in intelligent operation and maintenance,conventional sequential-modular-based process-simulation techniques present challenges regarding computationally intensive calculations and significant central processing unit(CPU)time requirements,particularly in large-scale design and optimization tasks.To address these challenges,this paper proposes a novel process-simulation parallel computing framework(PSPCF).This framework achieves layered parallelism in recycling processes at the unit operation level.Notably,PSPCF introduces a groundbreaking concept of formulating simulation problems as task graphs and utilizes Taskflow,an advanced task graph computing system,for hierarchical parallel scheduling and the execution of unit operation tasks.PSPCF also integrates an advanced work-stealing scheme to automatically balance thread resources with the demanding workload of unit operation tasks.For evaluation,both a simpler parallel column process and a more complex cracked gas separation process were simulated on a flowsheeting platform using PSPCF.The framework demonstrates significant time savings,achieving over 60%reduction in processing time for the simpler process and a 35%–40%speed-up for the more complex separation process.
基金Supported by the National Natural Science Foundation of China(42362026)Key R&D Project of Xinjiang Uygur Autonomous Region(2024B01015).
文摘In the ultra-deep strata of the Tarim Basin,the vertical growth process of strike-slip faults remains unclear,and the vertical distribution of fractured-cavity carbonate reservoirs is complex.This paper investigates the vertical growth process of strike-slip faults through field outcrop observations in the Keping area,interpretation of seismic data from the Fuman Oilfield,Tarim Basim,NW China,and structural physical simulation experiments.The results are obtained mainly in four aspects.First,field outcrops and ultra-deep seismic profiles indicate a three-layer structure within the strike-slip fault,consisting of fault core,fracture zone and primary rock.The fault core can be classified into three parts vertically:fracture-cavity unit,fault clay and breccia zone.The distribution of fracture-cavity units demonstrates a distinct pattern of vertical stratification,owing to the structural characteristics and growth process of the slip-strike fault.Second,the ultra-deep seismic profiles show multiple fracture-cavity units in the strike-slip fault zone.These units can be classified into four types:top fractured,middle connected,deep terminated,and intra-layer fractured.Third,structural physical simulation experiments and ultra-deep seismic data interpretation reveal that the strike-slip faults have evolved vertically in three stages:segmental rupture,vertical growth,and connection and extension.The particle image velocimetry detection demonstrates that the initial fracture of the fault zone occurred at the top or bottom and then evolved into cavities gradually along with the fault growth,accompanied by the emergence of new fractures in the middle part of the strata,which subsequently connected with the deep and shallow cavities to form a complete fault zone.Fourth,the ultra-deep carbonate strata primarily develop three types of fractured-cavity reservoirs:flower-shaped fracture,large and deep fault and staggered overlap.The first two types are larger in size with better reservoir conditions,suggesting a significant exploration potential.
基金supported by the PetroChina Prospective,Basic,and Strategic Technology Research Project(No.2021DJ0606).
文摘There is little low-and-high frequency information on seismic data in seismic exploration,resulting in narrower bandwidth and lower seismic resolution.It considerably restricts the prediction accuracy of thin reservoirs and thin interbeds.This study proposes a novel method to constrain improving seismic resolution in the time and frequency domain.The expected wavelet spectrum is used in the frequency domain to broaden the seismic spectrum range and increase the octave.In the time domain,the Frobenius vector regularization of the Hessian matrix is used to constrain the horizontal continuity of the seismic data.It eff ectively protects the signal-to-noise ratio of seismic data while the longitudinal seismic resolution is improved.This method is applied to actual post-stack seismic data and pre-stack gathers dividedly.Without abolishing the phase characteristics of the original seismic data,the time resolution is signifi cantly improved,and the structural features are clearer.Compared with the traditional spectral simulation and deconvolution methods,the frequency distribution is more reasonable,and seismic data has higher resolution.
基金supported by National Natural Sciences Foundation of China(No.62271165,62027802,62201307)the Guangdong Basic and Applied Basic Research Foundation(No.2023A1515030297)+2 种基金the Shenzhen Science and Technology Program ZDSYS20210623091808025Stable Support Plan Program GXWD20231129102638002the Major Key Project of PCL(No.PCL2024A01)。
文摘Due to the restricted satellite payloads in LEO mega-constellation networks(LMCNs),remote sensing image analysis,online learning and other big data services desirably need onboard distributed processing(OBDP).In existing technologies,the efficiency of big data applications(BDAs)in distributed systems hinges on the stable-state and low-latency links between worker nodes.However,LMCNs with high-dynamic nodes and long-distance links can not provide the above conditions,which makes the performance of OBDP hard to be intuitively measured.To bridge this gap,a multidimensional simulation platform is indispensable that can simulate the network environment of LMCNs and put BDAs in it for performance testing.Using STK's APIs and parallel computing framework,we achieve real-time simulation for thousands of satellite nodes,which are mapped as application nodes through software defined network(SDN)and container technologies.We elaborate the architecture and mechanism of the simulation platform,and take the Starlink and Hadoop as realistic examples for simulations.The results indicate that LMCNs have dynamic end-to-end latency which fluctuates periodically with the constellation movement.Compared to ground data center networks(GDCNs),LMCNs deteriorate the computing and storage job throughput,which can be alleviated by the utilization of erasure codes and data flow scheduling of worker nodes.
基金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 the New Century Excellent Talents in University(NCET-09-0396)the National Science&Technology Key Projects of Numerical Control(2012ZX04014-031)+1 种基金the Natural Science Foundation of Hubei Province(2011CDB279)the Foundation for Innovative Research Groups of the Natural Science Foundation of Hubei Province,China(2010CDA067)
文摘When castings become complicated and the demands for precision of numerical simulation become higher,the numerical data of casting numerical simulation become more massive.On a general personal computer,these massive numerical data may probably exceed the capacity of available memory,resulting in failure of rendering.Based on the out-of-core technique,this paper proposes a method to effectively utilize external storage and reduce memory usage dramatically,so as to solve the problem of insufficient memory for massive data rendering on general personal computers.Based on this method,a new postprocessor is developed.It is capable to illustrate filling and solidification processes of casting,as well as thermal stess.The new post-processor also provides fast interaction to simulation results.Theoretical analysis as well as several practical examples prove that the memory usage and loading time of the post-processor are independent of the size of the relevant files,but the proportion of the number of cells on surface.Meanwhile,the speed of rendering and fetching of value from the mouse is appreciable,and the demands of real-time and interaction are satisfied.
基金supported by College of William and Mary,Virginia Institute of Marine Science for the study environment
文摘Large eddy simulation （LES） using the Smagorinsky eddy viscosity model is added to the two-dimensional nine velocity components （D2Q9） lattice Boltzmann equation （LBE） with multi-relaxation-time （MRT） to simulate incompressible turbulent cavity flows with the Reynolds numbers up to 1 × 10^7. To improve the computation efficiency of LBM on the numerical simulations of turbulent flows, the massively parallel computing power from a graphic processing unit （GPU） with a computing unified device architecture （CUDA） is introduced into the MRT-LBE-LES model. The model performs well, compared with the results from others, with an increase of 76 times in computation efficiency. It appears that the higher the Reynolds numbers is, the smaller the Smagorinsky constant should be, if the lattice number is fixed. Also, for a selected high Reynolds number and a selected proper Smagorinsky constant, there is a minimum requirement for the lattice number so that the Smagorinsky eddy viscosity will not be excessively large.
文摘The absorbing process in isolating and coating process of α-olefin drag reducing polymer was studied by molecular dynamic simulation method, on basis of coating theory of α-olefin drag reducing polymer particles with polyurethane as coating material. The distributions of sodium laurate, sodium dodeeyl sulfate, and sodium dodeeyl benzene sulfonate on the surface of α-olefin drag reducing polymer particles were almost the same, but the bending degrees of them were obviously different. The bending degree of SLA molecules was greater than those of the other two surfactant molecules. Simulation results of absorbing and accumulating structure showed that, though hydrophobie properties of surfactant molecules were almost the same, water density around long chain sulfonate sodium was bigger than that around alkyl sulfate sodium. This property goes against useful absorbing and accumulating on the surface of α-olefin drag reducing polymer particles; simulation results of interactions of different surfactant and multiple hydroxyl compounds on surface of particles showed that, interactions of different surfaetant and one kind of multiple hydroxyl compound were similar to those of one kind of surfaetant and different multiple hydroxyl compounds. These two contrast types of interactions also exhibited the differences of absorbing distribution and closing degrees to surface of particles. The sequence of closing degrees was derived from simulation; control step of addition polymerization interaction in coating process was absorbing mass transfer process, so the more closed to surface of particle the multiple hydroxyl compounds were, the easier interactions With isoeyanate were. Simulation results represented the compatibility relationship between surfactant and multiple hydroxyl compounds. The isolating and coating processes of α-olefin drag reducing polymer were further understood on molecule and atom level through above simulation research, and based on the simulation, a referenced theoretical basis was provided for practical optimal selection and experimental preparation of α-olefin drag reducing polymer particles suspension isolation agent.
基金the Basic Research Program of the Korea Science & Engineering Foundation (KoSEF, No. R01-2006-000-10786-0).
文摘This article presents a multiscale simulation approach starting at the molecular level for the adsorption process development. A grand canonical Monte Carlo method is used for the prediction of adsorption isotherms of methanol on an activated carbon at the molecular level. The adsorption isotherms obtained in the linear region （or adsorption constant） are exploited as a model parameter required for the adsorption process simulation. The adsorption process model described by a set of partial differential equations （PDEs） is solved by using the conservation element and solution element method, which produces a fast and an accurate numerical solution to PDEs. The simulation results obtained from the adsorption constant estimated at the molecular level are in good agreement with the experimental results of the pulse response. The systematical multiscale simulation approach addressed in this study may be useful to accelerate the adsorption process development by reducing the number of experiments.
基金Supported by the National Defense Basic ResearchFoundation (K1503063165)
文摘According to the characteristic of Team Software Process （TSP）, it adopts a hierarchy-based model combined discrete event model with system dynamics model. This model represents TSP as form of two levels, the inner level embodies the continuity of the software process, the outer embodies the software development process by phases, and the structure and principle of the model is explained in detail, then formalization description of the model is offered. At last, an example is presented to demonstrate the simulation process and result. This model can simulate team software process from various angles, supervise and predict the software process. Also it can make the management of software development become more scientific and improve the quality of software.
基金National Natural Science Foundation of China, No.41125005 Knowledge Innovation Program of the Chinese Academy of Sciences, No.KACX1-YW-1001
文摘Physical geography and human geography are the principal branches of the geographical sciences. Physical process simulation and human process simulation in geography are both quantitative methods used to recover past events and even to forecast events based on precisely determined parameters. There are four differences between physical process simulation and human process simulation in geography, which we summarize with two specific cases, one of which is about a typhoon＇s development and its precipitation, and the other of which is regarding the evolution of three industrial structures in China. The differences focus on four aspects： the main factors of the research framework; the knowledge back-ground of the systematic analysis framework; the simulation data sources and quantitative method; and the core of the study object and the method of forecast application. As the human-land relationship is the key ideology of the man-land system, the relationship between the physical and human factors is becoming increasingly close at present. Physical process simulation and human process simulation in geography will exhibit crossing and blending in the future to reflect the various geographical phenomena better.
基金National Basic Research Program of China (2006CB605200)National Natural Science Foundation of China (50634010, 50674008)Program for New Century Excellent Talents in University (NCET-06-0083)
文摘Since processing parameters have always been assumed to be stable in the current finite element numerical simulation of dieless drawing process, the simulation results for the product dimension tend to stabilize gradually. In fact, the dimension fluctuation exists in the forming process all the while. A mathematical model of Gauss distribution for processing parameters was employed and a finite element numerical model of dieless drawing process with non-steady processing parameters was established. Dieless drawing processing of Ni-Ti alloy wire was conducted for verifying the proposed model. The results indicated that the non-steady processing parameters model had higher simulation accuracy of the wire diameter than that given by the steady parameters model. Furthermore, the model could also be used to analyze the fluctuation characteristics in the whole dieless drawing process.
基金Project supported by the National Natural Science Foundation of China （Grant Nos 50271026 and 50571037）.
文摘A molecular dynamics simulation study has been performed for the microstructure evolution in a liquid metal Ni system during crystallization process at two cooling rates by adopting the embedded atom method （EAM） model potential. The bond-type index method of Honeycutt-Andersen （HA） and a new cluster-type index method （CTIM-2） have been used to detect and analyse the microstructures in this system. It is demonstrated that the cooling rate plays a critical role in the microstructure evolution： below the crystallization temperature Tc, the effects of cooling rate are very remarkable and can be fully displayed. At different cooling rates of 2.0 × 10^13 K·s^-1 and 1.0 × 10^12 K·s^-1, two different kinds of crystal structures are obtained in the system. The first one is the coexistence of the hcp （expressed by （12 0 0 0 6 6） in CTIM-2） and the fcc （12 0 0 0 12 0） basic clusters consisting of 1421 and 1422 bond-types, and the hcp basic cluster becomes the dominant one with decreasing temperature, the second one is mainly the fcc （12 0 0 0 12 0） basic clusters consisting of 1421 bond-type, and their crystallization temperatures Tc would be 1073 and 1173 K, respectively.
基金National Natural Science Foundation of China, No.40671186 No.40271012
文摘The western coast of Hainan Island exhibits a savanna landscape. Many types of sand dunes, including transverse dune ridges, longitudinal dune ridges, elliptical dunes, coppice dunes, and climbing dunes, are widely distributed in the coastal zone. In winter, high-frequency and high-energy NE winds （dominant winds） are prevalent, with a resultant drift direction （RDD） of S35.6°W. In spring, low-frequency and low-energy SW secondary winds prevail, with a RDD of N25.1°E. Wind tunnel simulations revealed that the airflow over the dune surface is the main factor controlling the erosion and deposition patterns of dune surfaces and the morphological development of dunes. In the region＇s bidirectional wind environment, with two seasonally distinct energy levels, the airflow over the surface of elliptical dunes, barchan dunes, and transverse dune ridges will exhibit a transverse pattern, whereas the airflow over longitudinal dunes ridges exhibits a lateral pattern and that over climbing dunes exhibits a climbing-circumfluent pattern. These patterns represent different dynamic processes. The coastal dunes on the western coast of Hainan Island are influenced by factors such as onshore winds, sand sources, coastal slopes, rivers, and forest shelter belts. The source of the sand that supplements these dunes particularly influences the development pattern： when there is more sand, the pattern shows positive equilibrium deposition between dune ridges and dunes; otherwise, it shows negative equilibrium deposition. The presence or absence of forest shelter belts also influences deposition and dune development patterns and transformation of dune forms. Coastal dunes and inland desert dunes experience similar dynamic processes, but the former have more diversified shapes and more complex formation mechanisms.
基金The National Science and Technology Major Project of China under contract No.2011ZX05025-002-02the National Natural Science Foundation of China under contract Nos 41476032,91028009 and 40806019
文摘Turbidity currents represent a major agent for sediment transport in lakes, seas and oceans. In particu-lar, they formulate the most significant clastic accumulations in the deep sea, which become many of the world＇s most important hydrocarbon reservoirs. Several boreholes in the Qiongdongnan Basin, the north-western South China Sea, have recently revealed turbidity current deposits as significant hydrocarbon res-ervoirs. However, there are some arguments for the potential provenances. To solve this problem, it is es-sential to delineate their sedimentary processes as well as to evaluate their qualities as reservoir. Numerical simulations have been developed rapidly over the last several years, offering insights into turbidity current behaviors, as geologically significant turbidity currents are difficult to directly investigate due to their large scale and often destructive nature. Combined with the interpretation of the turbidity system based on high-resolution 3D seismic data, the paleotophography is acquired via a back-stripping seismic profile integrated with a borehole, i.e., Well A, in the western Qiongdongnan Basin; then a numerical model is built on the basis of this back-stripped profile. After defining the various turbidity current initial boundary conditions, includ-ing grain size, velocity and sediment concentration, the structures and behaviors of turbidity currents are investigated via numerical simulation software ANSYS FLUENT. Finally, the simulated turbidity deposits are compared with the interpreted sedimentary bodies based on 3D seismic data and the potential provenances of the revealed turbidites by Well A are discussed in details. The simulation results indicate that a sedimen-tary body develops far away from its source with an average grain size of 0.1 mm, i.e., sand-size sediment. Taking into account the location and orientation of the simulated seismic line, the consistence between normal forward simulation results and the revealed cores in Well A indicates that the turbidites should have been transported from Vietnam instead of Hainan Island. This interpretation has also been verified by the planar maps of sedimentary systems based on integration of boreholes and seismic data. The identification of the turbidity provenance will benefit the evaluation of extensively distributed submarine fans for hydro-carbon exploration in the deepwater areas.
基金Supported by the National Natural Science Foundation of China (No.60421002).
文摘A new 18-lump kinetic model for naphtha catalytic reforming reactions is discussed. By developing this model as a user module, a whole industrial continuous catalytic reforming process is simulated on Aspen plus plat-form. The technique utilizes the strong databases, complete sets of modules, and flexible simulation tools of the Aspen plus system and retains the characteristics of the proposed kinetic model. The calculated results are in fair agreement with the actual operating data. Based on the model of the whole reforming process, the process is opti-mized and the optimization results are tested in the actual industrial unit for about two months. The test shows that the process profit increases about 1000yuan·h-1 averagely, which is close to the calculated result.
文摘Numerical simulation of casting＇s mold filling process is the main and the most important aspect of the foundry CAE technology. But it is time-consuming; it may take dozens of hours or several days. While with the development of computer hardware, numerical simulation of casting＇ s mold filling process has made rapid progress. The simulation results, therefore, have become more and more practical. This study tries to find some clues of the computational time of mold filling process. Firstly, this paper introduces mathematic model and the basic route of numerical simulation of casting＇s mold filling process. Then the computational time of mold filling process has been carefully studied, and some new and useful results have been gained from the study of the computational time. Finally, this paper has given some real applications of numerical simulation of casting＇s mold filling process.
文摘A 3-D mathematical model considering turbulence phenomena has been established based on a computational fluid dynamics technique, so called 3-D SOLA-VOF (Solution Algorithm-Volume of Fluid), to simulate the fluid flow of mold filling process of die casting. In addition, the mathematical model for simulating the heat transfer in die casting process has also been established. The computation program has been developed by the authors with the finite difference method (FDM) recently. As verification, the mold filling process of a S-shaped die casting has been simulated and the simulation results coincide with that of the benchmark test. Finally, as a practical application, the gating design of a motorcycle component was modified by the mold filling simulation and the dies design of another motorcycle component was optimized by the heat transfer simulation. All the optimized designs were verified by the production practice.
