Low-Earth-Orbit satellite constellation networks(LEO-SCN)can provide low-cost,largescale,flexible coverage wireless communication services.High dynamics and large topological sizes characterize LEO-SCN.Protocol develo...Low-Earth-Orbit satellite constellation networks(LEO-SCN)can provide low-cost,largescale,flexible coverage wireless communication services.High dynamics and large topological sizes characterize LEO-SCN.Protocol development and application testing of LEO-SCN are challenging to carry out in a natural environment.Simulation platforms are a more effective means of technology demonstration.Currently available simulators have a single function and limited simulation scale.There needs to be a simulator for full-featured simulation.In this paper,we apply the parallel discrete-event simulation technique to the simulation of LEO-SCN to support large-scale complex system simulation at the packet level.To solve the problem that single-process programs cannot cope with complex simulations containing numerous entities,we propose a parallel mechanism and algorithms LP-NM and LP-YAWNS for synchronization.In the experiment,we use ns-3 to verify the acceleration ratio and efficiency of the above algorithms.The results show that our proposed mechanism can provide parallel simulation engine support for the LEO-SCN.展开更多
With wider applications of power electronic devices in modern power systems,simulation using traditional electro-mechanical and electromagnetic simulation tools suffer from low speed and imprecision.Multi-rate technol...With wider applications of power electronic devices in modern power systems,simulation using traditional electro-mechanical and electromagnetic simulation tools suffer from low speed and imprecision.Multi-rate technologies can greatly improve simulation efficiency by avoiding simulating the entire system using a small time-step.However,the drawbacks of the current synchronization mechanisms is that they introduce numerical errors and numerical instabilities in multi-rate parallel simulations.An improved multi-rate parallel technology,node splitting interface(NSI),is proposed to reduce errors and enhance simulation stability.A new synchronization mechanism is used to avoid prediction and signal delays.Theoretical analyses are carried out to prove the convergence and absolute stability of the proposed NSI algorithm.This algorithm is particularly suitable for simultaneously investigating long term dynamics of DC grids and fast transients of power electronic converters.展开更多
The geometry of joints has a significant influence on the mechanical properties of rocks.To simplify the curved joint shapes in rocks,the joint shape is usually treated as straight lines or planes in most laboratory e...The geometry of joints has a significant influence on the mechanical properties of rocks.To simplify the curved joint shapes in rocks,the joint shape is usually treated as straight lines or planes in most laboratory experiments and numerical simulations.In this study,the computerized tomography (CT) scanning and photogrammetry were employed to obtain the internal and surface joint structures of a limestone sample,respectively.To describe the joint geometry,the edge detection algorithms and a three-dimensional (3D) matrix mapping method were applied to reconstruct CT-based and photogrammetry-based jointed rock models.For comparison tests,the numerical uniaxial compression tests were conducted on an intact rock sample and a sample with a joint simplified to a plane using the parallel computing method.The results indicate that the mechanical characteristics and failure process of jointed rocks are significantly affected by the geometry of joints.The presence of joints reduces the uniaxial compressive strength (UCS),elastic modulus,and released acoustic emission (AE) energy of rocks by 37%–67%,21%–24%,and 52%–90%,respectively.Compared to the simplified joint sample,the proposed photogrammetry-based numerical model makes the most of the limited geometry information of joints.The UCS,accumulative released AE energy,and elastic modulus of the photogrammetry-based sample were found to be very close to those of the CT-based sample.The UCS value of the simplified joint sample (i.e.38.5 MPa) is much lower than that of the CT-based sample (i.e.72.3 MPa).Additionally,the accumulative released AE energy observed in the simplified joint sample is 3.899 times lower than that observed in the CT-based sample.CT scanning provides a reliable means to visualize the joints in rocks,which can be used to verify the reliability of photogrammetry techniques.The application of the photogrammetry-based sample enables detailed analysis for estimating the mechanical properties of jointed rocks.展开更多
A new three dimensional simulation method is introduced to study the workspace of a 6 PSS (P denotes a prismatic kinematic pair, S denotes a spherical kinematic pair) parallel machine tool. This algorithm adopts the...A new three dimensional simulation method is introduced to study the workspace of a 6 PSS (P denotes a prismatic kinematic pair, S denotes a spherical kinematic pair) parallel machine tool. This algorithm adopts the method of numerical analysis to investigate the boundary points in a series of sections which form the surface of the workspace. That is, to study such points that have the largest polar radius on a certain section in a system of polar coordinates according to conditions of constraint. The constraint conditions considered in the article include the maximum and minimum displacements of each dieblock, the maximum and minimum angles of oscillation in each hinge. By converting the constraint inequalities into constraint equations, the largest polar radius corresponding to every constraint condition can be evaluated and the minimum one is used to decide the boundary point. This algorithm greatly simplifies the computational process and can be used to analyze any section of the workspace. It provides a theoretical basis for the structural design of such a machine tool.展开更多
Traditional simulation methods are unable to meet the requirements of lunar takeo simulations, such as high force output precision, low cost, and repeated use. Considering that cable-driven parallel mechanisms have th...Traditional simulation methods are unable to meet the requirements of lunar takeo simulations, such as high force output precision, low cost, and repeated use. Considering that cable-driven parallel mechanisms have the advantages of high payload to weight ratio, potentially large workspace, and high-speed motion, these mechanisms have the potential to be used for lunar takeo simulations. Thus, this paper presents a parallel mechanism driven by nine cables. The purpose of this study is to optimize the dimensions of the cable-driven parallel mechanism to meet dynamic workspace requirements under cable tension constraints. The dynamic workspace requirements are derived from the kinematical function requests of the lunar takeo simulation equipment. Experimental design and response surface methods are adopted for building the surrogate mathematical model linking the optimal variables and the optimization indices. A set of dimensional parameters are determined by analyzing the surrogate mathematical model. The volume of the dynamic workspace increased by 46% after optimization. Besides, a force control method is proposed for calculating output vector and sinusoidal forces. A force control loop is introduced into the traditional position control loop to adjust the cable force precisely, while controlling the cable length. The e ectiveness of the proposed control method is verified through experiments. A 5% vector output accuracy and 12 Hz undulation force output can be realized. This paper proposes a cable-driven parallel mechanism which can be used for lunar takeo simulation.展开更多
In this paper a class of real-time parallel modified Rosenbrock methods of numerical simulation is constructed for stiff dynamic systems on a multiprocessor system, and convergence and numerical stability of these met...In this paper a class of real-time parallel modified Rosenbrock methods of numerical simulation is constructed for stiff dynamic systems on a multiprocessor system, and convergence and numerical stability of these methods are discussed. A-stable real-time parallel formula of two-stage third-order and A(α)-stable real-time parallel formula with o ≈ 89.96° of three-stage fourth-order are particularly given. The numerical simulation experiments in parallel environment show that the class of algorithms is efficient and applicable, with greater speedup.展开更多
This paper is concerned with three-dimensional numerical simulation of a plunging liquid jet. The transient processes of forming an air cavity around the jet, capturing an initially large air bubble, and the break-up ...This paper is concerned with three-dimensional numerical simulation of a plunging liquid jet. The transient processes of forming an air cavity around the jet, capturing an initially large air bubble, and the break-up of this large toroidal-shaped bubble into smaller bubbles were analyzed. A stabilized finite element method (FEM) was employed under parallel numerical simulations based on adaptive, unstructured grid and coupled with a level-set method to track the interface between air and liquid. These simulations show that the inertia of the liquid jet initially depresses the pool's surface, forming an annular air cavity which surrounds the liquid jet. A toroidal liquid eddy which is subse- quently formed in the liquid pool results in air cavity collapse, and in turn entrains air into the liquid pool from the unstable annular air gap region around the liquid jet.展开更多
The cavity formation and propagation process of stress wave from parallel hole cut blasting was simulated with ANSYS/LS-DYNA 3D nonlinear dynamic finite element software. The distribution of element plastic strain, no...The cavity formation and propagation process of stress wave from parallel hole cut blasting was simulated with ANSYS/LS-DYNA 3D nonlinear dynamic finite element software. The distribution of element plastic strain, node velocity, node time-acceleration history and the blasting cartridge volume ratio during the process were analyzed. It was found that the detonation of charged holes would cause the interaction of stress wave with the wall of uncharged holes. Initial rock cracking and displacement to neighboring uncharged holes become the main mechanism of cavity formation in early stage.2008 University of Science and Technology Beijing. All rights reserved.展开更多
In this paper a dynamical model of a two degree of freedom,tendon based,parallel manipulator (TBPM) system is proposed. The motion control methods of the TBPM system were designed. Using MATLAB,the motion control simu...In this paper a dynamical model of a two degree of freedom,tendon based,parallel manipulator (TBPM) system is proposed. The motion control methods of the TBPM system were designed. Using MATLAB,the motion control simulation of this model TBPM system was implemented in preparation for actual experiments. The results of the simulation demonstrated that the response time of the system was in a reasonable range,the motion behavior of the platform was stable and the tension forces acting on the tendons were in a safe range and acceptable. Furthermore,the parameters of the controllers were optimized using MATLAB and better results for the time response were obtained.展开更多
From such actual conditions as the effects of characteristics of miltilayer petroleum geology and permeation fluid mechanics, a new numerical model is put forward and coupling splitting-up implicit interactive scheme ...From such actual conditions as the effects of characteristics of miltilayer petroleum geology and permeation fluid mechanics, a new numerical model is put forward and coupling splitting-up implicit interactive scheme is formed. For the actual situation of Dongying hollow (four-layer) and Tanhai region (three-layer) of Shengli Petroleum Field, the numerical simulation test results and the actual conditions are coincident.展开更多
The desire to increase spatial and temporal resolution in modeling groundwater system has led to the requirement for intensive computational ability and large memory space. In the course of satisfying such requirement...The desire to increase spatial and temporal resolution in modeling groundwater system has led to the requirement for intensive computational ability and large memory space. In the course of satisfying such requirement, parallel computing has played a core role over the past several decades. This paper reviews the parallel algebraic linear solution methods and the parallel implementation technologies for groundwater simulation. This work is carried out to provide guidance to enable modelers of groundwater systems to make sensible choices when developing solution methods based upon the current state of knowledge in parallel computing.展开更多
In this paper, a parallel simulation algorithm for the control problem in differential algebraic system is presented. The error of the algorithm is estimated. The stability analysis is made for a model problem and the...In this paper, a parallel simulation algorithm for the control problem in differential algebraic system is presented. The error of the algorithm is estimated. The stability analysis is made for a model problem and the stability region is given. The numerical example demonstrates that the method is efficient.展开更多
Background Petrochemical products possess a high risk of flammability,explosivity,and toxicity,making petrochemical accidents exceedingly destructive.Therefore,disaster analysis,prediction,and real-time simulations ha...Background Petrochemical products possess a high risk of flammability,explosivity,and toxicity,making petrochemical accidents exceedingly destructive.Therefore,disaster analysis,prediction,and real-time simulations have become important means of controlling and reducing accident hazards.Methods This study proposes a complete real-time simulation solution of gas diffusion with coordinate and concentration data,which was mainly aimed at simulating the types of harmful gas leakage and diffusion accidents in the petrochemical industry.The rendering effect was more continuous and accurate through grid homogenization and trilinear interpolation.This study presents a data processing and rendering parallelization process to enhance simulation efficiency.Gas concentration and fragment transparency were combined to synthesize transparent pixels in a scene.To ensure the approximate accuracy of the rendering effect,improve the efficiency of real-time rendering,and meet the requirement of intuitive perception using concentration data,a weighted blended order-independent transparency(OIT)with enhanced alpha weight is presented,which can provide a more intuitive perception of the hierarchical information of concentration data while preserving depth information.This study compares and analyzes three OIT algorithms-depth peeling,weighted blended OIT,and weighted blended OIT with enhanced alpha weight-in terms of rendering image quality,rendering time,required memory,and hierarchical information.Results Using weighted blended OIT with an enhanced alpha weight technique,the rendering time was shortened by 53.2%compared with that of the depth peeling algorithm,and the texture memory required was significantly smaller than that of the depth peeling algorithm.The rendering results of weighted blended OIT with an enhanced alpha weight were approximately accurate compared with those of the depth peeling algorithm as the ground truth,and there was no popping when surfaces passed through one another.Simultaneously,compared with weighted blended OIT,weighted blended OIT with an enhanced alpha weight achieved an intuitive perception of the hierarchical information of concentration data.展开更多
The primary goal of cloth simulation is to express object behavior in a realistic manner and achieve real-time performance by following the fundamental concept of physic.In general,the mass–spring system is applied t...The primary goal of cloth simulation is to express object behavior in a realistic manner and achieve real-time performance by following the fundamental concept of physic.In general,the mass–spring system is applied to real-time cloth simulation with three types of springs.However,hard spring cloth simulation using the mass–spring system requires a small integration time-step in order to use a large stiffness coefficient.Furthermore,to obtain stable behavior,constraint enforcement is used instead of maintenance of the force of each spring.Constraint force computation involves a large sparse linear solving operation.Due to the large computation,we implement a cloth simulation using adaptive constraint activation and deactivation techniques that involve the mass-spring system and constraint enforcement method to prevent excessive elongation of cloth.At the same time,when the length of the spring is stretched or compressed over a defined threshold,adaptive constraint activation and deactivation method deactivates the spring and generate the implicit constraint.Traditional method that uses a serial process of the Central Processing Unit(CPU)to solve the system in every frame cannot handle the complex structure of cloth model in real-time.Our simulation utilizes the Graphic Processing Unit(GPU)parallel processing with compute shader in OpenGL Shading Language(GLSL)to solve the system effectively.In this paper,we design and implement parallel method for cloth simulation,and experiment on the performance and behavior comparison of the mass-spring system,constraint enforcement,and adaptive constraint activation and deactivation techniques the using GPU-based parallel method.展开更多
Two kinds of 2-dof parallel mechanisms are proposed in this paper which can be used as the actuator for the plane sprayer. The direct and inverse kinematics solutions of the two kinds of mechanisms are derived on the ...Two kinds of 2-dof parallel mechanisms are proposed in this paper which can be used as the actuator for the plane sprayer. The direct and inverse kinematics solutions of the two kinds of mechanisms are derived on the end operating point and two workspaces are analyzed and compared. The kinematics models of the end operating point of two mechanisms are simulated by Matlab examples obtaining variation of kinematics parameters of these two mechanisms. The research of this paper provides the basis for the selection of mechanism, trajectory planning of the end operating point on the sprayer and often some practical value for trajectory analysis and structure design of the plane sprayer.展开更多
A parallel arithmetic program for the molecular dynamics (MD) simulation study of a large sized system consisting of 50 000100 000 atoms of liquid metals is reformed, based on the cascade arithmetic program used for t...A parallel arithmetic program for the molecular dynamics (MD) simulation study of a large sized system consisting of 50 000100 000 atoms of liquid metals is reformed, based on the cascade arithmetic program used for the molecular dynamics simulation study of a small sized system consisting of 5001 000 atoms. The program is used to simulate the rapid solidification processes of liquid metal Al system. Some new results, such as larger clusters composed of more than 36 smaller clusters (icosahedra or defect icosahedra) obtained in the system of 50 000 atoms, however, the larger clusters can not be seen in the small sized system of 5001 000 atoms. On the other hand, the results from this simulation study would be more closed to the real situation of the system under consideration because the influence of boundary conditions is decreased remarkably. It can be expected that from the parallel algorithm combined with the higher performance super computer, the total number of atoms in simulation system can be enlarged again up to tens, even hundreds times in the near future.展开更多
This paper presents partially asynchronous parallel simulation of continuous-system (PAPSoCS) and some approaches to the issues of its implementation on a multicomputer system. To guarantee the simulation results cor...This paper presents partially asynchronous parallel simulation of continuous-system (PAPSoCS) and some approaches to the issues of its implementation on a multicomputer system. To guarantee the simulation results correct and speedup the simulation, the scheme for efficient PAPSoCS is proposed and the virtual topology star is constructed to match the path of message passing for solving algorithm-architecture adequation problem. Under the circumstances that messages frequently passed inter-processor are much shorter, typically within several 4 bytes, asynchronous communication mode is employed to reduce the communication ratio. Experiment results show that asynchronous parallel simulation has much higher efficiency than its synchronous counterpart.展开更多
In this paper, a homogenous parallel simulation system is presented in detail for continuous--system simulation. The system is collstructed by a host computer and I I transputers connected into a topologyof 'Super...In this paper, a homogenous parallel simulation system is presented in detail for continuous--system simulation. The system is collstructed by a host computer and I I transputers connected into a topologyof 'Super--Node' which is very suitable for simulation of stiff systems. An automatic software interface runin the host is developed to partition simulation model, either equations or block diagrams, into several equitable segments and then pack them into parallel simulation program to be executed in the parallel system.This interface frees simulation users from parallel programming to focus on their simulation experiments.展开更多
We employ the parallel computing technology to study numerically the three-dimensional structure of quantized vortices of Bose-Einstein condensates, For anisotropic cases, the bending process of vortices is described ...We employ the parallel computing technology to study numerically the three-dimensional structure of quantized vortices of Bose-Einstein condensates, For anisotropic cases, the bending process of vortices is described in detail by the decrease of Gross-Pitaevskii energy. A completely straight vortex and the steady and symmetrical multiple-vortex configurations are obtained. We analyse the effect of initial conditions and angular velocity on the number and shape of vortices.展开更多
In this paper, the stability analysis for parallel real-time digital simulation models is discussed. The coupling coefficient perturbation method and the simulation stepsize perturbation method are established. For tw...In this paper, the stability analysis for parallel real-time digital simulation models is discussed. The coupling coefficient perturbation method and the simulation stepsize perturbation method are established. For two classes of systems of test equations, we construct the parallel simulation models and prove that they have the stability behaviour which is similar to the original continuous systems.展开更多
基金supported by Jiangsu Provincial Key Research and Development Program (No.BE20210132)the Zhejiang Provincial Key Research and Development Program (No.2021C01040)the team of S-SET
文摘Low-Earth-Orbit satellite constellation networks(LEO-SCN)can provide low-cost,largescale,flexible coverage wireless communication services.High dynamics and large topological sizes characterize LEO-SCN.Protocol development and application testing of LEO-SCN are challenging to carry out in a natural environment.Simulation platforms are a more effective means of technology demonstration.Currently available simulators have a single function and limited simulation scale.There needs to be a simulator for full-featured simulation.In this paper,we apply the parallel discrete-event simulation technique to the simulation of LEO-SCN to support large-scale complex system simulation at the packet level.To solve the problem that single-process programs cannot cope with complex simulations containing numerous entities,we propose a parallel mechanism and algorithms LP-NM and LP-YAWNS for synchronization.In the experiment,we use ns-3 to verify the acceleration ratio and efficiency of the above algorithms.The results show that our proposed mechanism can provide parallel simulation engine support for the LEO-SCN.
基金This work was supported in part by the People Programme(Marie Curie Actions)of the European Union’s Seventh Framework Programme FP7/2007-2013/under REA grant agreement(No.317221)project title MEDOW,in part supported by the Project of National Science Foundation for Distinguished Young Scholars of China(No.51407164).
文摘With wider applications of power electronic devices in modern power systems,simulation using traditional electro-mechanical and electromagnetic simulation tools suffer from low speed and imprecision.Multi-rate technologies can greatly improve simulation efficiency by avoiding simulating the entire system using a small time-step.However,the drawbacks of the current synchronization mechanisms is that they introduce numerical errors and numerical instabilities in multi-rate parallel simulations.An improved multi-rate parallel technology,node splitting interface(NSI),is proposed to reduce errors and enhance simulation stability.A new synchronization mechanism is used to avoid prediction and signal delays.Theoretical analyses are carried out to prove the convergence and absolute stability of the proposed NSI algorithm.This algorithm is particularly suitable for simultaneously investigating long term dynamics of DC grids and fast transients of power electronic converters.
基金supported by the National Natural Science Foundation of China(Grant Nos.42277150,41977219)Henan Provincial Science and Technology Research Project(Grant No.222102320271).
文摘The geometry of joints has a significant influence on the mechanical properties of rocks.To simplify the curved joint shapes in rocks,the joint shape is usually treated as straight lines or planes in most laboratory experiments and numerical simulations.In this study,the computerized tomography (CT) scanning and photogrammetry were employed to obtain the internal and surface joint structures of a limestone sample,respectively.To describe the joint geometry,the edge detection algorithms and a three-dimensional (3D) matrix mapping method were applied to reconstruct CT-based and photogrammetry-based jointed rock models.For comparison tests,the numerical uniaxial compression tests were conducted on an intact rock sample and a sample with a joint simplified to a plane using the parallel computing method.The results indicate that the mechanical characteristics and failure process of jointed rocks are significantly affected by the geometry of joints.The presence of joints reduces the uniaxial compressive strength (UCS),elastic modulus,and released acoustic emission (AE) energy of rocks by 37%–67%,21%–24%,and 52%–90%,respectively.Compared to the simplified joint sample,the proposed photogrammetry-based numerical model makes the most of the limited geometry information of joints.The UCS,accumulative released AE energy,and elastic modulus of the photogrammetry-based sample were found to be very close to those of the CT-based sample.The UCS value of the simplified joint sample (i.e.38.5 MPa) is much lower than that of the CT-based sample (i.e.72.3 MPa).Additionally,the accumulative released AE energy observed in the simplified joint sample is 3.899 times lower than that observed in the CT-based sample.CT scanning provides a reliable means to visualize the joints in rocks,which can be used to verify the reliability of photogrammetry techniques.The application of the photogrammetry-based sample enables detailed analysis for estimating the mechanical properties of jointed rocks.
基金Ministerial Level Foundation(96J185 .1BQ0150) Fund for Reasearch on Doctoral Programs in Institutions of Higher Learning(1997000716)
文摘A new three dimensional simulation method is introduced to study the workspace of a 6 PSS (P denotes a prismatic kinematic pair, S denotes a spherical kinematic pair) parallel machine tool. This algorithm adopts the method of numerical analysis to investigate the boundary points in a series of sections which form the surface of the workspace. That is, to study such points that have the largest polar radius on a certain section in a system of polar coordinates according to conditions of constraint. The constraint conditions considered in the article include the maximum and minimum displacements of each dieblock, the maximum and minimum angles of oscillation in each hinge. By converting the constraint inequalities into constraint equations, the largest polar radius corresponding to every constraint condition can be evaluated and the minimum one is used to decide the boundary point. This algorithm greatly simplifies the computational process and can be used to analyze any section of the workspace. It provides a theoretical basis for the structural design of such a machine tool.
基金Supported by National Natural Science Foundation of China(Grant No.51405024)
文摘Traditional simulation methods are unable to meet the requirements of lunar takeo simulations, such as high force output precision, low cost, and repeated use. Considering that cable-driven parallel mechanisms have the advantages of high payload to weight ratio, potentially large workspace, and high-speed motion, these mechanisms have the potential to be used for lunar takeo simulations. Thus, this paper presents a parallel mechanism driven by nine cables. The purpose of this study is to optimize the dimensions of the cable-driven parallel mechanism to meet dynamic workspace requirements under cable tension constraints. The dynamic workspace requirements are derived from the kinematical function requests of the lunar takeo simulation equipment. Experimental design and response surface methods are adopted for building the surrogate mathematical model linking the optimal variables and the optimization indices. A set of dimensional parameters are determined by analyzing the surrogate mathematical model. The volume of the dynamic workspace increased by 46% after optimization. Besides, a force control method is proposed for calculating output vector and sinusoidal forces. A force control loop is introduced into the traditional position control loop to adjust the cable force precisely, while controlling the cable length. The e ectiveness of the proposed control method is verified through experiments. A 5% vector output accuracy and 12 Hz undulation force output can be realized. This paper proposes a cable-driven parallel mechanism which can be used for lunar takeo simulation.
基金This project was supported by the National Natural Science Foundation of China (No. 19871080).
文摘In this paper a class of real-time parallel modified Rosenbrock methods of numerical simulation is constructed for stiff dynamic systems on a multiprocessor system, and convergence and numerical stability of these methods are discussed. A-stable real-time parallel formula of two-stage third-order and A(α)-stable real-time parallel formula with o ≈ 89.96° of three-stage fourth-order are particularly given. The numerical simulation experiments in parallel environment show that the class of algorithms is efficient and applicable, with greater speedup.
基金supported by the Office of Naval Research(Grant ONRDC14292111)
文摘This paper is concerned with three-dimensional numerical simulation of a plunging liquid jet. The transient processes of forming an air cavity around the jet, capturing an initially large air bubble, and the break-up of this large toroidal-shaped bubble into smaller bubbles were analyzed. A stabilized finite element method (FEM) was employed under parallel numerical simulations based on adaptive, unstructured grid and coupled with a level-set method to track the interface between air and liquid. These simulations show that the inertia of the liquid jet initially depresses the pool's surface, forming an annular air cavity which surrounds the liquid jet. A toroidal liquid eddy which is subse- quently formed in the liquid pool results in air cavity collapse, and in turn entrains air into the liquid pool from the unstable annular air gap region around the liquid jet.
文摘The cavity formation and propagation process of stress wave from parallel hole cut blasting was simulated with ANSYS/LS-DYNA 3D nonlinear dynamic finite element software. The distribution of element plastic strain, node velocity, node time-acceleration history and the blasting cartridge volume ratio during the process were analyzed. It was found that the detonation of charged holes would cause the interaction of stress wave with the wall of uncharged holes. Initial rock cracking and displacement to neighboring uncharged holes become the main mechanism of cavity formation in early stage.2008 University of Science and Technology Beijing. All rights reserved.
基金Project HI 370/24-1 supported by the German Research Community (DFG)
文摘In this paper a dynamical model of a two degree of freedom,tendon based,parallel manipulator (TBPM) system is proposed. The motion control methods of the TBPM system were designed. Using MATLAB,the motion control simulation of this model TBPM system was implemented in preparation for actual experiments. The results of the simulation demonstrated that the response time of the system was in a reasonable range,the motion behavior of the platform was stable and the tension forces acting on the tendons were in a safe range and acceptable. Furthermore,the parameters of the controllers were optimized using MATLAB and better results for the time response were obtained.
文摘From such actual conditions as the effects of characteristics of miltilayer petroleum geology and permeation fluid mechanics, a new numerical model is put forward and coupling splitting-up implicit interactive scheme is formed. For the actual situation of Dongying hollow (four-layer) and Tanhai region (three-layer) of Shengli Petroleum Field, the numerical simulation test results and the actual conditions are coincident.
基金supported by the National Basic Research Program (973 Program) of China under Grant No.2010CB428804 and 2011CB 309702
文摘The desire to increase spatial and temporal resolution in modeling groundwater system has led to the requirement for intensive computational ability and large memory space. In the course of satisfying such requirement, parallel computing has played a core role over the past several decades. This paper reviews the parallel algebraic linear solution methods and the parallel implementation technologies for groundwater simulation. This work is carried out to provide guidance to enable modelers of groundwater systems to make sensible choices when developing solution methods based upon the current state of knowledge in parallel computing.
文摘In this paper, a parallel simulation algorithm for the control problem in differential algebraic system is presented. The error of the algorithm is estimated. The stability analysis is made for a model problem and the stability region is given. The numerical example demonstrates that the method is efficient.
基金Supported by National Key R&D Program of China (2020YFB1710400)。
文摘Background Petrochemical products possess a high risk of flammability,explosivity,and toxicity,making petrochemical accidents exceedingly destructive.Therefore,disaster analysis,prediction,and real-time simulations have become important means of controlling and reducing accident hazards.Methods This study proposes a complete real-time simulation solution of gas diffusion with coordinate and concentration data,which was mainly aimed at simulating the types of harmful gas leakage and diffusion accidents in the petrochemical industry.The rendering effect was more continuous and accurate through grid homogenization and trilinear interpolation.This study presents a data processing and rendering parallelization process to enhance simulation efficiency.Gas concentration and fragment transparency were combined to synthesize transparent pixels in a scene.To ensure the approximate accuracy of the rendering effect,improve the efficiency of real-time rendering,and meet the requirement of intuitive perception using concentration data,a weighted blended order-independent transparency(OIT)with enhanced alpha weight is presented,which can provide a more intuitive perception of the hierarchical information of concentration data while preserving depth information.This study compares and analyzes three OIT algorithms-depth peeling,weighted blended OIT,and weighted blended OIT with enhanced alpha weight-in terms of rendering image quality,rendering time,required memory,and hierarchical information.Results Using weighted blended OIT with an enhanced alpha weight technique,the rendering time was shortened by 53.2%compared with that of the depth peeling algorithm,and the texture memory required was significantly smaller than that of the depth peeling algorithm.The rendering results of weighted blended OIT with an enhanced alpha weight were approximately accurate compared with those of the depth peeling algorithm as the ground truth,and there was no popping when surfaces passed through one another.Simultaneously,compared with weighted blended OIT,weighted blended OIT with an enhanced alpha weight achieved an intuitive perception of the hierarchical information of concentration data.
基金supported by the Basic Science Research Program through the National Research Foundation of Korea(NRF-2019R1F1A1062752)funded by the Ministry of Education+1 种基金funded by BK21 FOUR(Fostering Outstanding Universities for Research)(No.:5199990914048)supported by the Soonchunhyang University Research Fund.
文摘The primary goal of cloth simulation is to express object behavior in a realistic manner and achieve real-time performance by following the fundamental concept of physic.In general,the mass–spring system is applied to real-time cloth simulation with three types of springs.However,hard spring cloth simulation using the mass–spring system requires a small integration time-step in order to use a large stiffness coefficient.Furthermore,to obtain stable behavior,constraint enforcement is used instead of maintenance of the force of each spring.Constraint force computation involves a large sparse linear solving operation.Due to the large computation,we implement a cloth simulation using adaptive constraint activation and deactivation techniques that involve the mass-spring system and constraint enforcement method to prevent excessive elongation of cloth.At the same time,when the length of the spring is stretched or compressed over a defined threshold,adaptive constraint activation and deactivation method deactivates the spring and generate the implicit constraint.Traditional method that uses a serial process of the Central Processing Unit(CPU)to solve the system in every frame cannot handle the complex structure of cloth model in real-time.Our simulation utilizes the Graphic Processing Unit(GPU)parallel processing with compute shader in OpenGL Shading Language(GLSL)to solve the system effectively.In this paper,we design and implement parallel method for cloth simulation,and experiment on the performance and behavior comparison of the mass-spring system,constraint enforcement,and adaptive constraint activation and deactivation techniques the using GPU-based parallel method.
文摘Two kinds of 2-dof parallel mechanisms are proposed in this paper which can be used as the actuator for the plane sprayer. The direct and inverse kinematics solutions of the two kinds of mechanisms are derived on the end operating point and two workspaces are analyzed and compared. The kinematics models of the end operating point of two mechanisms are simulated by Matlab examples obtaining variation of kinematics parameters of these two mechanisms. The research of this paper provides the basis for the selection of mechanism, trajectory planning of the end operating point on the sprayer and often some practical value for trajectory analysis and structure design of the plane sprayer.
文摘A parallel arithmetic program for the molecular dynamics (MD) simulation study of a large sized system consisting of 50 000100 000 atoms of liquid metals is reformed, based on the cascade arithmetic program used for the molecular dynamics simulation study of a small sized system consisting of 5001 000 atoms. The program is used to simulate the rapid solidification processes of liquid metal Al system. Some new results, such as larger clusters composed of more than 36 smaller clusters (icosahedra or defect icosahedra) obtained in the system of 50 000 atoms, however, the larger clusters can not be seen in the small sized system of 5001 000 atoms. On the other hand, the results from this simulation study would be more closed to the real situation of the system under consideration because the influence of boundary conditions is decreased remarkably. It can be expected that from the parallel algorithm combined with the higher performance super computer, the total number of atoms in simulation system can be enlarged again up to tens, even hundreds times in the near future.
文摘This paper presents partially asynchronous parallel simulation of continuous-system (PAPSoCS) and some approaches to the issues of its implementation on a multicomputer system. To guarantee the simulation results correct and speedup the simulation, the scheme for efficient PAPSoCS is proposed and the virtual topology star is constructed to match the path of message passing for solving algorithm-architecture adequation problem. Under the circumstances that messages frequently passed inter-processor are much shorter, typically within several 4 bytes, asynchronous communication mode is employed to reduce the communication ratio. Experiment results show that asynchronous parallel simulation has much higher efficiency than its synchronous counterpart.
文摘In this paper, a homogenous parallel simulation system is presented in detail for continuous--system simulation. The system is collstructed by a host computer and I I transputers connected into a topologyof 'Super--Node' which is very suitable for simulation of stiff systems. An automatic software interface runin the host is developed to partition simulation model, either equations or block diagrams, into several equitable segments and then pack them into parallel simulation program to be executed in the parallel system.This interface frees simulation users from parallel programming to focus on their simulation experiments.
基金Project supported partly by the National Natural Science Foundation of China (Grant Nos 10301034 and 40574069), The authors thank Professor Du Q very much for his important discussions.
文摘We employ the parallel computing technology to study numerically the three-dimensional structure of quantized vortices of Bose-Einstein condensates, For anisotropic cases, the bending process of vortices is described in detail by the decrease of Gross-Pitaevskii energy. A completely straight vortex and the steady and symmetrical multiple-vortex configurations are obtained. We analyse the effect of initial conditions and angular velocity on the number and shape of vortices.
基金This work is supported partly by the National Natural Science Foundation of China
文摘In this paper, the stability analysis for parallel real-time digital simulation models is discussed. The coupling coefficient perturbation method and the simulation stepsize perturbation method are established. For two classes of systems of test equations, we construct the parallel simulation models and prove that they have the stability behaviour which is similar to the original continuous systems.