Effective partitioning is crucial for enabling parallel restoration of power systems after blackouts.This paper proposes a novel partitioning method based on deep reinforcement learning.First,the partitioning decision...Effective partitioning is crucial for enabling parallel restoration of power systems after blackouts.This paper proposes a novel partitioning method based on deep reinforcement learning.First,the partitioning decision process is formulated as a Markov decision process(MDP)model to maximize the modularity.Corresponding key partitioning constraints on parallel restoration are considered.Second,based on the partitioning objective and constraints,the reward function of the partitioning MDP model is set by adopting a relative deviation normalization scheme to reduce mutual interference between the reward and penalty in the reward function.The soft bonus scaling mechanism is introduced to mitigate overestimation caused by abrupt jumps in the reward.Then,the deep Q network method is applied to solve the partitioning MDP model and generate partitioning schemes.Two experience replay buffers are employed to speed up the training process of the method.Finally,case studies on the IEEE 39-bus test system demonstrate that the proposed method can generate a high-modularity partitioning result that meets all key partitioning constraints,thereby improving the parallelism and reliability of the restoration process.Moreover,simulation results demonstrate that an appropriate discount factor is crucial for ensuring both the convergence speed and the stability of the partitioning training.展开更多
We propose a symplectic partitioned Runge-Kutta （SPRK） method with eighth-order spatial accuracy based on the extended Hamiltonian system of the acoustic waveequation. Known as the eighth-order NSPRK method, this te...We propose a symplectic partitioned Runge-Kutta （SPRK） method with eighth-order spatial accuracy based on the extended Hamiltonian system of the acoustic waveequation. Known as the eighth-order NSPRK method, this technique uses an eighth-orderaccurate nearly analytic discrete （NAD） operator to discretize high-order spatial differentialoperators and employs a second-order SPRK method to discretize temporal derivatives.The stability criteria and numerical dispersion relations of the eighth-order NSPRK methodare given by a semi-analytical method and are tested by numerical experiments. We alsoshow the differences of the numerical dispersions between the eighth-order NSPRK methodand conventional numerical methods such as the fourth-order NSPRK method, the eighth-order Lax-Wendroff correction （LWC） method and the eighth-order staggered-grid （SG）method. The result shows that the ability of the eighth-order NSPRK method to suppress thenumerical dispersion is obviously superior to that of the conventional numerical methods. Inthe same computational environment, to eliminate visible numerical dispersions, the eighth-order NSPRK is approximately 2.5 times faster than the fourth-order NSPRK and 3.4 timesfaster than the fourth-order SPRK, and the memory requirement is only approximately47.17% of the fourth-order NSPRK method and 49.41% of the fourth-order SPRK method,which indicates the highest computational efficiency. Modeling examples for the two-layermodels such as the heterogeneous and Marmousi models show that the wavefields generatedby the eighth-order NSPRK method are very clear with no visible numerical dispersion.These numerical experiments illustrate that the eighth-order NSPRK method can effectivelysuppress numerical dispersion when coarse grids are adopted. Therefore, this methodcan greatly decrease computer memory requirement and accelerate the forward modelingproductivity. In general, the eighth-order NSPRK method has tremendous potential value forseismic exploration and seismology research.展开更多
A partition of unity finite element method for numerical simulation of short wave propagation in solids is presented. The finite element spaces were constructed by multiplying the standard isoparametric finite element...A partition of unity finite element method for numerical simulation of short wave propagation in solids is presented. The finite element spaces were constructed by multiplying the standard isoparametric finite element shape functions, which form a partition of unity, with the local subspaces defined on the corresponding shape functions, which include a priori knowledge about the wave motion equation in trial spaces and approximately reproduce the highly oscillatory properties within a single element. Numerical examples demonstrate the performance of the proposed partition of unity finite element in both computational accuracy and efficiency.展开更多
The impact dynamics of a flexible multibody system is investigated. By using a partition method, the system is divided into two parts, the local impact region and the region away from the impact. The two parts are con...The impact dynamics of a flexible multibody system is investigated. By using a partition method, the system is divided into two parts, the local impact region and the region away from the impact. The two parts are connected by specific boundary conditions, and the system after partition is equivalent to the original system. According to the rigid-flexible coupling dynamic theory of multibody system, system＇s rigid-flexible coupling dynamic equations without impact are derived. A local impulse method for establishing the initial impact conditions is proposed. It satisfies the compatibility con- ditions for contact constraints and the actual physical situation of the impact process of flexible bodies. Based on the contact constraint method, system＇s impact dynamic equa- tions are derived in a differential-algebraic form. The contact/separation criterion and the algorithm are given. An impact dynamic simulation is given. The results show that system＇s dynamic behaviors including the energy, the deformations, the displacements, and the impact force during the impact process change dramatically. The impact makes great effects on the global dynamics of the system during and after impact.展开更多
In order to research the influence of coal-rock mass morphology of mining space on the flow law of gas,the laboratory physical model and numerical computation methods were adopted to simulate coal mining activities.Th...In order to research the influence of coal-rock mass morphology of mining space on the flow law of gas,the laboratory physical model and numerical computation methods were adopted to simulate coal mining activities.The simulation results indicate that,after coal seam mining,the loose rock accumulation body of free caving,ordered rock arrangement body of plate damage rich in longitudinal and transverse fractures and horizontal fissure body formed by rock mass deformation imbalance are formed from bottom to top in the mining space.For these three types of accumulation bodies,there are essential differences in the accumulation state,rock size and gas breakover characteristics.According to this,the coal-rock mass in the mining space is classified into gas turbulence channel area,gas transitional flow channel area and gas seepage channel area.In the turbulence channel area,the gas is distributed transversely and longitudinally and gas diffuses in the form of convection with Reynolds number R_e more than100;in the transitional flow channel area,one-way or two-way gas channels are crisscross and gas is of transitional flow regime with R,.between 10 and 100.In the seepage channel area,there are a few vertical gas channels with R,.less than 10.In this paper,the researches on the gas orientation method in different partitions were further carried out,gas orientation methods of low-level pipe burying,middle-level interception and high-level extraction were determined and an on-site industrial test was conducted,achieving the effective diversion of gas and verifying the reasonableness of gas channel partition.展开更多
The impact problem of a flexible multibody system is a non-smooth, high-transient, and strong-nonlinear dynamic process with variable boundary. How to model the contact/impact process accurately and efficiently is one...The impact problem of a flexible multibody system is a non-smooth, high-transient, and strong-nonlinear dynamic process with variable boundary. How to model the contact/impact process accurately and efficiently is one of the main difficulties in many engineering applications. The numerical approaches being used widely in impact analysis are mainly from two fields： multibody system dynamics （MBS） and computational solid mechanics （CSM）. Approaches based on MBS provide a more efficient yet less accurate analysis of the contact/impact problems, while approaches based on CSM are well suited for particularly high accuracy needs, yet require very high computational effort. To bridge the gap between accuracy and efficiency in the dynamic simulation of a flexible multibody system with contacts/impacts, a partition method is presented considering that the contact body is divided into two parts, an impact region and a non-impact region. The impact region is modeled using the finite element method to guarantee the local accuracy, while the non-impact region is modeled using the modal reduction approach to raise the global efficiency. A three-dimensional rod-plate impact experiment is designed and performed to validate the numerical results. The principle for how to partition the contact bodies is proposed： the maximum radius of the impact region can be estimated by an analytical method, and the modal truncation orders of the non-impact region can be estimated by the highest frequency of the signal measured. The simulation results using the presented method are in good agreement with the experimental results. It shows that this method is an effec-rive formulation considering both accuracy and efficiency. Moreover, a more complicated multibody impact problem of a crank slider mechanism is investigated to strengthen this conclusion.展开更多
The average risk indices,such as the loss of load expectation(LOLE)and expected demand not supplied(EDNS),have been widely used in risk assessment of power systems.However,the average indices can't distinguish bet...The average risk indices,such as the loss of load expectation(LOLE)and expected demand not supplied(EDNS),have been widely used in risk assessment of power systems.However,the average indices can't distinguish between the events of low probability but high damage and the events of high probability but low damage.In order to ov+rcome these shortcomings,this paper proposes an extended risk analysis framework for the power system based on the partitioned multi-objective risk method(PMRM).展开更多
A partition of unity finite element method for numerical simulation of short wave propagation in solids is presented. The finite element spaces were constructed by multiplying the standard isoparametric finite element...A partition of unity finite element method for numerical simulation of short wave propagation in solids is presented. The finite element spaces were constructed by multiplying the standard isoparametric finite element shape functions, which form a partition of unity, with the local subspaces defined on the corresponding shape functions, which include a priori knowledge about the wave motion equation in trial spaces and approximately reproduce the highly oscillatory properties within a single element. Numerical examples demonstrate the performance of the proposed partition of unity finite element in both computational accuracy and efficiency.展开更多
Optimized calculations of 75 PCDDs and their parent DD were carded out at the B3LYP/6-31G＊ level by density functional theory （DFT） method. The structural parameters were obtained and significant correlation betwee...Optimized calculations of 75 PCDDs and their parent DD were carded out at the B3LYP/6-31G＊ level by density functional theory （DFT） method. The structural parameters were obtained and significant correlation between the C1 substitution position and some structural parameters was found. Consequently, the number of C1 substitution positions was taken as theoretical descriptors to establish two novel QSPR models for predicting lgKow and -lgSw of all PCDD congeners. The two models achieved in this work contain two variables （Na and Nβ）, of which r = 0.9312, 0.9965 and SD = 0.27, 0.12 respectively, and t values are all large. The variation inflation factors （VIF） of variables in the two models herein are both less than 5.0, suggesting high accuracy of the lgKow and -lgSw predicting models, and the results of cross-validation test also show that the two models exhibit optimum stability and good predictive power. By comparison, the correlation and predictive ability of the present work are more advantageous than those obtained using semi-empirical AM1 and GC-RI methods.展开更多
Hardware/software partitioning is an important step in the design of embedded systems. In this paper, the hardware/software partitioning problem is modeled as a constrained binary integer programming problem, which is...Hardware/software partitioning is an important step in the design of embedded systems. In this paper, the hardware/software partitioning problem is modeled as a constrained binary integer programming problem, which is further converted equivalently to an unconstrained binary integer programming problem by a penalty method. A local search method, HSFM, is developed to obtain a discrete local minimizer of the unconstrained binary integer programming problem. Next, an auxiliary function, which has the same global optimal solutions as the unconstrained binary integer programming problem, is constructed, and its properties are studied. We show that applying HSFM to minimize the auxiliary function can escape from previous local optima by the increase of the parameter value successfully. Finally, a discrete dynamic convexized method is developed to solve the hardware/software partitioning problem. Computational results and comparisons indicate that the proposed algorithm can get high-quality solutions.展开更多
Numerical quadrature is an important ingredient of Galerkin meshless methods. A new numerical quadrature technique, partition of unity quadrature （PUQ）,for Galerkin meshless methods was presented. The technique is b...Numerical quadrature is an important ingredient of Galerkin meshless methods. A new numerical quadrature technique, partition of unity quadrature （PUQ）,for Galerkin meshless methods was presented. The technique is based on finite covering and partition of unity. There is no need to decompose the physical domain into small cell. It possesses remarkable integration accuracy. Using Element-free Galerkin methods as example, Galerkin meshless methods based on PUQ were studied in detail. Meshing is always not required in the procedure of constitution of approximate function or numerical quadrature, so Galerkin meshless methods based on PUQ are “truly” meshless methods.展开更多
This study proposed a partitioned method to analyze maneuvering of insects during flapping flight.This method decomposed the insect flapping flight into wing and body subsystems and then coupled them via boundary cond...This study proposed a partitioned method to analyze maneuvering of insects during flapping flight.This method decomposed the insect flapping flight into wing and body subsystems and then coupled them via boundary conditions imposed on the wing’s base using one-way coupling.In the wing subsystem,the strong coupling of the flexible wings and surrounding fluid was accurately analyzed using the finite element method to obtain the thrust forces acting on the insect’s body.The resulting thrust forces were passed from the wing subsystem to the body subsystem,and then rigid body motion was analyzed in the body subsystem.The rolling,yawing,and pitching motions were simulated using the proposed method as follows:In the rolling simulation,the difference of the stroke angle between the right and left wings caused a roll torque.In the yawing simulation,the initial feathering angle in the right wing only caused a yaw torque.In the pitching simulation,the difference between the front-and back-stroke angles in both the right and left wings caused a pitch torque.All three torques generated maneuvering motion comparable with that obtained in actual observations of insect flight.These results demonstrate that the proposed method can adequately simulate the fundamental maneuvers of insect flapping flight.In the present simulations,the maneuvering mechanisms were investigated at the governing equation level,which might be difficult using other approaches.Therefore,the proposed method will contribute to revealing the underlying insect flight mechanisms.展开更多
The implicit partition algorithm used to solve fluid–structure coupling problems has high accuracy,but it requires a long computation time.In this paper,a semi-implicit fluid–structure coupling algorithm based on mo...The implicit partition algorithm used to solve fluid–structure coupling problems has high accuracy,but it requires a long computation time.In this paper,a semi-implicit fluid–structure coupling algorithm based on modal force prediction-correction is proposed to improve the computational efficiency.In the pre-processing stage,the fluid domain is assumed to be a pseudo-elastic solid and merged with the solid domain to form a holistic system,and the normalized modal information of the holistic system is calculated and stored.During the sub-step cycle,the modal superposition method is used to obtain the response of the holistic system with the predicted modal force as the load,so that the deformation of the structure and the updating of the fluid mesh can be achieved simultaneously.After solving the Reynolds-averaged Navier-Stokes equations in the fluid domain,the predicted modal force is corrected and a new sub-step cycle is started until the converged result is obtained.In this method,the computation of the fluid equations and the updating of the dynamic mesh are done implicitly,while the deformation of the structure is done explicitly.Two numerical cases,vortex induced oscillation of an elastic beam and fluid–structure interaction of a final stage blade,are used to verify the efficiency and accuracy of the proposed algorithm.The results show that the proposed method achieves the same accuracy as the implicit method while the computational time is reduced.In the case of the vortex-induced oscillation problem,the computational time can be reduced to 18.6%.In the case of the final stage blade vibration,the computational time can be reduced to 53.8%.展开更多
Considering the defects of conventional optimization methods, a novel optimization algorithm is introduced in this paper. Target space partitioning method is used in this algorithm to solve multi-objective optimizatio...Considering the defects of conventional optimization methods, a novel optimization algorithm is introduced in this paper. Target space partitioning method is used in this algorithm to solve multi-objective optimization problem, thus achieve the coherent solution which can meet the requirements of all target functions, and improve the population's overall evolution level. The algorithm which guarantees diversity preservation and fast convergence to the Pareto set is applied to structural optimization problems. The empirical analysis supports the algorithm and gives an example with program.展开更多
Matching optimization of resonant parameters among the high power inverters,low power transformers and plasma reactors have significant effects on the performance and output of the reactor array when applying the part...Matching optimization of resonant parameters among the high power inverters,low power transformers and plasma reactors have significant effects on the performance and output of the reactor array when applying the partitioned operation method.In this paper,the Matlab/Simulink electrical model was established based on the method of partitioned operation.The matching relation between resonant parameters is analyzed on the basis of experimental result.As a consequence,transformer leakage inductance and working frequency are the important parameters influencing the operational efficiency of system,leakage inductance of transformer should be adjusted based on the equivalent capacitance of plasma reactor to realize the matching optimization of resonant parameters.展开更多
The development and rapid usage of numerical codes for fluid-structure interaction(FSI) problems are of great relevance to researchers in many engineering fields such as civil engineering and ocean engineering. This m...The development and rapid usage of numerical codes for fluid-structure interaction(FSI) problems are of great relevance to researchers in many engineering fields such as civil engineering and ocean engineering. This multidisciplinary field known as FSI has been expanded to engineering fields such as offshore structures, tall slender structures and other flexible structures applications. The motivation of this paper is to investigate the numerical model of two-way coupling FSI partitioned flexible plate structure under fluid flow. The adopted partitioned method and approach utilized the advantage of the existing numerical algorithms in solving the two-way coupling fluid and structural interactions. The flexible plate was subjected to a fluid flow which causes large deformation on the fluid domain from the oscillation of the flexible plate. Both fluid and flexible plate are subjected to the interaction of load transfer within two physics by using the strong and weak coupling methods of MFS and Load Transfer Physics Environment, respectively. The oscillation and deformation results have been validated which demonstrate the reliability of both strong and weak method in resolving the two-way coupling problem in contribution of knowledge to the feasibility field study of ocean engineering and civil engineering.展开更多
A gradient-based optimization method for producing a contoured beam by using a single-fed reflector antenna is presented. First, a quick and accurate pattern approximation formula based on physical optics(PO) is adopt...A gradient-based optimization method for producing a contoured beam by using a single-fed reflector antenna is presented. First, a quick and accurate pattern approximation formula based on physical optics(PO) is adopted to calculate the gradients of the directivity with respect to reflector's nodal displacements. Because the approximation formula is a linear function of nodal displacements, the gradient can be easily derived. Then, the method of the steepest descent is adopted, and an optimization iteration procedure is proposed. The iteration procedure includes two loops: an inner loop and an outer loop. In the inner loop, the gradient and pattern are calculated by matrix operation, which is very fast by using the pre-calculated data in the outer loop. In the outer loop, the ideal terms used in the inner loop to calculate the gradient and pattern are updated, and the real pattern is calculated by the PO method. Due to the high approximation accuracy, when the outer loop is performed once, the inner loop can be performed many times, which will save much time because the integration is replaced by matrix operation. In the end, a contoured beam covering the continental United States(CONUS) is designed, and simulation results show the effectiveness of the proposed algorithm.展开更多
A two-grid partition of unity method for second order elliptic problems is proposed and analyzed. The standard two-grid method is a local and parallel method usually leading to a discontinuous solution in the entire c...A two-grid partition of unity method for second order elliptic problems is proposed and analyzed. The standard two-grid method is a local and parallel method usually leading to a discontinuous solution in the entire computational domain. Partition of unity method is employed to glue all the local solutions together to get the global continuous one, which is optimal in HI-norm. Furthermore, it is shown that the L^2 error can be improved by using the coarse grid correction. Numerical experiments are reported to support the theoretical results.展开更多
The inverse and direct piezoelectric and circuit coupling are widely observed in advanced electro-mechanical systems such as piezoelectric energy harvesters.Existing strongly coupled analysis methods based on direct n...The inverse and direct piezoelectric and circuit coupling are widely observed in advanced electro-mechanical systems such as piezoelectric energy harvesters.Existing strongly coupled analysis methods based on direct numerical modeling for this phenomenon can be classified into partitioned or monolithic formulations.Each formulation has its advantages and disadvantages,and the choice depends on the characteristics of each coupled problem.This study proposes a new option:a coupled analysis strategy that combines the best features of the existing formulations,namely,the hybrid partitioned-monolithic method.The analysis of inverse piezoelectricity and the monolithic analysis of direct piezoelectric and circuit interaction are strongly coupled using a partitioned iterative hierarchical algorithm.In a typical benchmark problem of a piezoelectric energy harvester,this research compares the results from the proposed method to those from the conventional strongly coupled partitioned iterative method,discussing the accuracy,stability,and computational cost.The proposed hybrid concept is effective for coupled multi-physics problems,including various coupling conditions.展开更多
Advanced array processing approaches require accurate knowledge of the location of individual element in a sensor array.Most array shape estimation methods require the directions of sources.In this paper,an array shap...Advanced array processing approaches require accurate knowledge of the location of individual element in a sensor array.Most array shape estimation methods require the directions of sources.In this paper,an array shape estimation method based on eigen-decomposition is presented.The directions of sources do not need to be considered in advance and optimal array shape is generated through a series of iterations.To further improve the accuracy of this algorithm,a partitioned eigenstructure method is introduced.Numerical simulations using non-partitioned and partitioned method are conducted to verify the performance of the proposed technique.展开更多
基金funded by the Beijing Engineering Research Center of Electric Rail Transportation.
文摘Effective partitioning is crucial for enabling parallel restoration of power systems after blackouts.This paper proposes a novel partitioning method based on deep reinforcement learning.First,the partitioning decision process is formulated as a Markov decision process(MDP)model to maximize the modularity.Corresponding key partitioning constraints on parallel restoration are considered.Second,based on the partitioning objective and constraints,the reward function of the partitioning MDP model is set by adopting a relative deviation normalization scheme to reduce mutual interference between the reward and penalty in the reward function.The soft bonus scaling mechanism is introduced to mitigate overestimation caused by abrupt jumps in the reward.Then,the deep Q network method is applied to solve the partitioning MDP model and generate partitioning schemes.Two experience replay buffers are employed to speed up the training process of the method.Finally,case studies on the IEEE 39-bus test system demonstrate that the proposed method can generate a high-modularity partitioning result that meets all key partitioning constraints,thereby improving the parallelism and reliability of the restoration process.Moreover,simulation results demonstrate that an appropriate discount factor is crucial for ensuring both the convergence speed and the stability of the partitioning training.
基金This research was supported by the National Natural Science Foundation of China （Nos. 41230210 and 41204074）, the Science Foundation of the Education Department of Yunnan Province （No. 2013Z152）, and Statoil Company （Contract No. 4502502663）.
文摘We propose a symplectic partitioned Runge-Kutta （SPRK） method with eighth-order spatial accuracy based on the extended Hamiltonian system of the acoustic waveequation. Known as the eighth-order NSPRK method, this technique uses an eighth-orderaccurate nearly analytic discrete （NAD） operator to discretize high-order spatial differentialoperators and employs a second-order SPRK method to discretize temporal derivatives.The stability criteria and numerical dispersion relations of the eighth-order NSPRK methodare given by a semi-analytical method and are tested by numerical experiments. We alsoshow the differences of the numerical dispersions between the eighth-order NSPRK methodand conventional numerical methods such as the fourth-order NSPRK method, the eighth-order Lax-Wendroff correction （LWC） method and the eighth-order staggered-grid （SG）method. The result shows that the ability of the eighth-order NSPRK method to suppress thenumerical dispersion is obviously superior to that of the conventional numerical methods. Inthe same computational environment, to eliminate visible numerical dispersions, the eighth-order NSPRK is approximately 2.5 times faster than the fourth-order NSPRK and 3.4 timesfaster than the fourth-order SPRK, and the memory requirement is only approximately47.17% of the fourth-order NSPRK method and 49.41% of the fourth-order SPRK method,which indicates the highest computational efficiency. Modeling examples for the two-layermodels such as the heterogeneous and Marmousi models show that the wavefields generatedby the eighth-order NSPRK method are very clear with no visible numerical dispersion.These numerical experiments illustrate that the eighth-order NSPRK method can effectivelysuppress numerical dispersion when coarse grids are adopted. Therefore, this methodcan greatly decrease computer memory requirement and accelerate the forward modelingproductivity. In general, the eighth-order NSPRK method has tremendous potential value forseismic exploration and seismology research.
基金Project supported by the National Basic Research Program of China (973Project) (No.2002CB412709) and the National Natural Science Foundation of China (Nos.50278012,10272027,19832010)
文摘A partition of unity finite element method for numerical simulation of short wave propagation in solids is presented. The finite element spaces were constructed by multiplying the standard isoparametric finite element shape functions, which form a partition of unity, with the local subspaces defined on the corresponding shape functions, which include a priori knowledge about the wave motion equation in trial spaces and approximately reproduce the highly oscillatory properties within a single element. Numerical examples demonstrate the performance of the proposed partition of unity finite element in both computational accuracy and efficiency.
基金supported by the National Natural Science Foundation of China(Nos.11132007,11272155,and 10772085)the Fundamental Research Funds for the Central Universities(No.30920130112009)the 333 Project of Jiangsu Province of China(No.BRA2011172)
文摘The impact dynamics of a flexible multibody system is investigated. By using a partition method, the system is divided into two parts, the local impact region and the region away from the impact. The two parts are connected by specific boundary conditions, and the system after partition is equivalent to the original system. According to the rigid-flexible coupling dynamic theory of multibody system, system＇s rigid-flexible coupling dynamic equations without impact are derived. A local impulse method for establishing the initial impact conditions is proposed. It satisfies the compatibility con- ditions for contact constraints and the actual physical situation of the impact process of flexible bodies. Based on the contact constraint method, system＇s impact dynamic equa- tions are derived in a differential-algebraic form. The contact/separation criterion and the algorithm are given. An impact dynamic simulation is given. The results show that system＇s dynamic behaviors including the energy, the deformations, the displacements, and the impact force during the impact process change dramatically. The impact makes great effects on the global dynamics of the system during and after impact.
基金Financial supports for this work,provided by the State Key Basic Research Program of China(No.2011CB201204)
文摘In order to research the influence of coal-rock mass morphology of mining space on the flow law of gas,the laboratory physical model and numerical computation methods were adopted to simulate coal mining activities.The simulation results indicate that,after coal seam mining,the loose rock accumulation body of free caving,ordered rock arrangement body of plate damage rich in longitudinal and transverse fractures and horizontal fissure body formed by rock mass deformation imbalance are formed from bottom to top in the mining space.For these three types of accumulation bodies,there are essential differences in the accumulation state,rock size and gas breakover characteristics.According to this,the coal-rock mass in the mining space is classified into gas turbulence channel area,gas transitional flow channel area and gas seepage channel area.In the turbulence channel area,the gas is distributed transversely and longitudinally and gas diffuses in the form of convection with Reynolds number R_e more than100;in the transitional flow channel area,one-way or two-way gas channels are crisscross and gas is of transitional flow regime with R,.between 10 and 100.In the seepage channel area,there are a few vertical gas channels with R,.less than 10.In this paper,the researches on the gas orientation method in different partitions were further carried out,gas orientation methods of low-level pipe burying,middle-level interception and high-level extraction were determined and an on-site industrial test was conducted,achieving the effective diversion of gas and verifying the reasonableness of gas channel partition.
基金supported by the National Natural Science Foundation of China (Grants 11772188, 11132007)
文摘The impact problem of a flexible multibody system is a non-smooth, high-transient, and strong-nonlinear dynamic process with variable boundary. How to model the contact/impact process accurately and efficiently is one of the main difficulties in many engineering applications. The numerical approaches being used widely in impact analysis are mainly from two fields： multibody system dynamics （MBS） and computational solid mechanics （CSM）. Approaches based on MBS provide a more efficient yet less accurate analysis of the contact/impact problems, while approaches based on CSM are well suited for particularly high accuracy needs, yet require very high computational effort. To bridge the gap between accuracy and efficiency in the dynamic simulation of a flexible multibody system with contacts/impacts, a partition method is presented considering that the contact body is divided into two parts, an impact region and a non-impact region. The impact region is modeled using the finite element method to guarantee the local accuracy, while the non-impact region is modeled using the modal reduction approach to raise the global efficiency. A three-dimensional rod-plate impact experiment is designed and performed to validate the numerical results. The principle for how to partition the contact bodies is proposed： the maximum radius of the impact region can be estimated by an analytical method, and the modal truncation orders of the non-impact region can be estimated by the highest frequency of the signal measured. The simulation results using the presented method are in good agreement with the experimental results. It shows that this method is an effec-rive formulation considering both accuracy and efficiency. Moreover, a more complicated multibody impact problem of a crank slider mechanism is investigated to strengthen this conclusion.
文摘The average risk indices,such as the loss of load expectation(LOLE)and expected demand not supplied(EDNS),have been widely used in risk assessment of power systems.However,the average indices can't distinguish between the events of low probability but high damage and the events of high probability but low damage.In order to ov+rcome these shortcomings,this paper proposes an extended risk analysis framework for the power system based on the partitioned multi-objective risk method(PMRM).
文摘A partition of unity finite element method for numerical simulation of short wave propagation in solids is presented. The finite element spaces were constructed by multiplying the standard isoparametric finite element shape functions, which form a partition of unity, with the local subspaces defined on the corresponding shape functions, which include a priori knowledge about the wave motion equation in trial spaces and approximately reproduce the highly oscillatory properties within a single element. Numerical examples demonstrate the performance of the proposed partition of unity finite element in both computational accuracy and efficiency.
基金This work was supported by the National Natural Science Foundation of China (No. 20737001)(2003CB415002)China Postdoctoral Science Foundation (No. 2003033486)
文摘Optimized calculations of 75 PCDDs and their parent DD were carded out at the B3LYP/6-31G＊ level by density functional theory （DFT） method. The structural parameters were obtained and significant correlation between the C1 substitution position and some structural parameters was found. Consequently, the number of C1 substitution positions was taken as theoretical descriptors to establish two novel QSPR models for predicting lgKow and -lgSw of all PCDD congeners. The two models achieved in this work contain two variables （Na and Nβ）, of which r = 0.9312, 0.9965 and SD = 0.27, 0.12 respectively, and t values are all large. The variation inflation factors （VIF） of variables in the two models herein are both less than 5.0, suggesting high accuracy of the lgKow and -lgSw predicting models, and the results of cross-validation test also show that the two models exhibit optimum stability and good predictive power. By comparison, the correlation and predictive ability of the present work are more advantageous than those obtained using semi-empirical AM1 and GC-RI methods.
基金Supported by the National Natural Science Foundation of China(11301255)the Fund by Collaborative Innovation Center of IoT Industrialization and Intelligent Production,Minjiang University(IIC1703)+1 种基金Foundation of Minjiang University(MYK17032)the Program for New Century Excellent Talents in Fujian Province University
文摘Hardware/software partitioning is an important step in the design of embedded systems. In this paper, the hardware/software partitioning problem is modeled as a constrained binary integer programming problem, which is further converted equivalently to an unconstrained binary integer programming problem by a penalty method. A local search method, HSFM, is developed to obtain a discrete local minimizer of the unconstrained binary integer programming problem. Next, an auxiliary function, which has the same global optimal solutions as the unconstrained binary integer programming problem, is constructed, and its properties are studied. We show that applying HSFM to minimize the auxiliary function can escape from previous local optima by the increase of the parameter value successfully. Finally, a discrete dynamic convexized method is developed to solve the hardware/software partitioning problem. Computational results and comparisons indicate that the proposed algorithm can get high-quality solutions.
文摘Numerical quadrature is an important ingredient of Galerkin meshless methods. A new numerical quadrature technique, partition of unity quadrature （PUQ）,for Galerkin meshless methods was presented. The technique is based on finite covering and partition of unity. There is no need to decompose the physical domain into small cell. It possesses remarkable integration accuracy. Using Element-free Galerkin methods as example, Galerkin meshless methods based on PUQ were studied in detail. Meshing is always not required in the procedure of constitution of approximate function or numerical quadrature, so Galerkin meshless methods based on PUQ are “truly” meshless methods.
文摘This study proposed a partitioned method to analyze maneuvering of insects during flapping flight.This method decomposed the insect flapping flight into wing and body subsystems and then coupled them via boundary conditions imposed on the wing’s base using one-way coupling.In the wing subsystem,the strong coupling of the flexible wings and surrounding fluid was accurately analyzed using the finite element method to obtain the thrust forces acting on the insect’s body.The resulting thrust forces were passed from the wing subsystem to the body subsystem,and then rigid body motion was analyzed in the body subsystem.The rolling,yawing,and pitching motions were simulated using the proposed method as follows:In the rolling simulation,the difference of the stroke angle between the right and left wings caused a roll torque.In the yawing simulation,the initial feathering angle in the right wing only caused a yaw torque.In the pitching simulation,the difference between the front-and back-stroke angles in both the right and left wings caused a pitch torque.All three torques generated maneuvering motion comparable with that obtained in actual observations of insect flight.These results demonstrate that the proposed method can adequately simulate the fundamental maneuvers of insect flapping flight.In the present simulations,the maneuvering mechanisms were investigated at the governing equation level,which might be difficult using other approaches.Therefore,the proposed method will contribute to revealing the underlying insect flight mechanisms.
基金support of the National Natural Science Foundation of China(No.51675406)the Basic Research Project Group,China(No.514010106-205)。
文摘The implicit partition algorithm used to solve fluid–structure coupling problems has high accuracy,but it requires a long computation time.In this paper,a semi-implicit fluid–structure coupling algorithm based on modal force prediction-correction is proposed to improve the computational efficiency.In the pre-processing stage,the fluid domain is assumed to be a pseudo-elastic solid and merged with the solid domain to form a holistic system,and the normalized modal information of the holistic system is calculated and stored.During the sub-step cycle,the modal superposition method is used to obtain the response of the holistic system with the predicted modal force as the load,so that the deformation of the structure and the updating of the fluid mesh can be achieved simultaneously.After solving the Reynolds-averaged Navier-Stokes equations in the fluid domain,the predicted modal force is corrected and a new sub-step cycle is started until the converged result is obtained.In this method,the computation of the fluid equations and the updating of the dynamic mesh are done implicitly,while the deformation of the structure is done explicitly.Two numerical cases,vortex induced oscillation of an elastic beam and fluid–structure interaction of a final stage blade,are used to verify the efficiency and accuracy of the proposed algorithm.The results show that the proposed method achieves the same accuracy as the implicit method while the computational time is reduced.In the case of the vortex-induced oscillation problem,the computational time can be reduced to 18.6%.In the case of the final stage blade vibration,the computational time can be reduced to 53.8%.
基金National Natural Science Foundations of China (No. 60970004, No. 60743010)Natural Science Foundation of ShandongProvince, China (No. Z2008G02)
文摘Considering the defects of conventional optimization methods, a novel optimization algorithm is introduced in this paper. Target space partitioning method is used in this algorithm to solve multi-objective optimization problem, thus achieve the coherent solution which can meet the requirements of all target functions, and improve the population's overall evolution level. The algorithm which guarantees diversity preservation and fast convergence to the Pareto set is applied to structural optimization problems. The empirical analysis supports the algorithm and gives an example with program.
基金supported by National Key Technology Research and Development Program of the Ministry of Science and Technology of China(No.2013BAC06B02)Public Science and Technology Research Funds Projects of Ocean (No.201305027)+2 种基金National Natural Science Foundation of China(No.51877024)Liaoning Scientific Research Project of Department of Education of Liaoning Province(No.LZ2015007)High Level Talent Innovation Project of Dalian(No.2016RQ040)
文摘Matching optimization of resonant parameters among the high power inverters,low power transformers and plasma reactors have significant effects on the performance and output of the reactor array when applying the partitioned operation method.In this paper,the Matlab/Simulink electrical model was established based on the method of partitioned operation.The matching relation between resonant parameters is analyzed on the basis of experimental result.As a consequence,transformer leakage inductance and working frequency are the important parameters influencing the operational efficiency of system,leakage inductance of transformer should be adjusted based on the equivalent capacitance of plasma reactor to realize the matching optimization of resonant parameters.
文摘The development and rapid usage of numerical codes for fluid-structure interaction(FSI) problems are of great relevance to researchers in many engineering fields such as civil engineering and ocean engineering. This multidisciplinary field known as FSI has been expanded to engineering fields such as offshore structures, tall slender structures and other flexible structures applications. The motivation of this paper is to investigate the numerical model of two-way coupling FSI partitioned flexible plate structure under fluid flow. The adopted partitioned method and approach utilized the advantage of the existing numerical algorithms in solving the two-way coupling fluid and structural interactions. The flexible plate was subjected to a fluid flow which causes large deformation on the fluid domain from the oscillation of the flexible plate. Both fluid and flexible plate are subjected to the interaction of load transfer within two physics by using the strong and weak coupling methods of MFS and Load Transfer Physics Environment, respectively. The oscillation and deformation results have been validated which demonstrate the reliability of both strong and weak method in resolving the two-way coupling problem in contribution of knowledge to the feasibility field study of ocean engineering and civil engineering.
基金supported by the National Natural Science Foundation of China(51805399)the Fundamental Research Funds for the Central Universities(JB180403)+2 种基金the Chinese Academy of Sciences(CAS)"Light of West China" Program(2017-XBQNXZ-B-024)the National Basic Research Program of China(973 Program)(2015CB857100)the Operation,Maintenance and Upgrading Fund for Astronomical Telescopes and Facility Instruments,budgeted from the Ministry of Finance of China(MOF)and administrated by the CAS
文摘A gradient-based optimization method for producing a contoured beam by using a single-fed reflector antenna is presented. First, a quick and accurate pattern approximation formula based on physical optics(PO) is adopted to calculate the gradients of the directivity with respect to reflector's nodal displacements. Because the approximation formula is a linear function of nodal displacements, the gradient can be easily derived. Then, the method of the steepest descent is adopted, and an optimization iteration procedure is proposed. The iteration procedure includes two loops: an inner loop and an outer loop. In the inner loop, the gradient and pattern are calculated by matrix operation, which is very fast by using the pre-calculated data in the outer loop. In the outer loop, the ideal terms used in the inner loop to calculate the gradient and pattern are updated, and the real pattern is calculated by the PO method. Due to the high approximation accuracy, when the outer loop is performed once, the inner loop can be performed many times, which will save much time because the integration is replaced by matrix operation. In the end, a contoured beam covering the continental United States(CONUS) is designed, and simulation results show the effectiveness of the proposed algorithm.
基金Project supported by the National Natural Science Foundation of China(No.40074031)the Science Foundation of the Science and Technology Commission of Shanghai Municipalitythe Program for Young Excellent Talents in Tongji University(No.2007kj008)
文摘A two-grid partition of unity method for second order elliptic problems is proposed and analyzed. The standard two-grid method is a local and parallel method usually leading to a discontinuous solution in the entire computational domain. Partition of unity method is employed to glue all the local solutions together to get the global continuous one, which is optimal in HI-norm. Furthermore, it is shown that the L^2 error can be improved by using the coarse grid correction. Numerical experiments are reported to support the theoretical results.
基金supported by the Japan Society for the Promotion of Science,KAKENHI Grant No.23H00475.
文摘The inverse and direct piezoelectric and circuit coupling are widely observed in advanced electro-mechanical systems such as piezoelectric energy harvesters.Existing strongly coupled analysis methods based on direct numerical modeling for this phenomenon can be classified into partitioned or monolithic formulations.Each formulation has its advantages and disadvantages,and the choice depends on the characteristics of each coupled problem.This study proposes a new option:a coupled analysis strategy that combines the best features of the existing formulations,namely,the hybrid partitioned-monolithic method.The analysis of inverse piezoelectricity and the monolithic analysis of direct piezoelectric and circuit interaction are strongly coupled using a partitioned iterative hierarchical algorithm.In a typical benchmark problem of a piezoelectric energy harvester,this research compares the results from the proposed method to those from the conventional strongly coupled partitioned iterative method,discussing the accuracy,stability,and computational cost.The proposed hybrid concept is effective for coupled multi-physics problems,including various coupling conditions.
文摘Advanced array processing approaches require accurate knowledge of the location of individual element in a sensor array.Most array shape estimation methods require the directions of sources.In this paper,an array shape estimation method based on eigen-decomposition is presented.The directions of sources do not need to be considered in advance and optimal array shape is generated through a series of iterations.To further improve the accuracy of this algorithm,a partitioned eigenstructure method is introduced.Numerical simulations using non-partitioned and partitioned method are conducted to verify the performance of the proposed technique.
