In the era of big data, government, business and personal digital information will be possible for data mining. Data mining requires massive data as a support. However, the direct release of the original mass data, wh...In the era of big data, government, business and personal digital information will be possible for data mining. Data mining requires massive data as a support. However, the direct release of the original mass data, which usually contain some sensitive information of personal or analysis, will result in leakage of user privacy. Therefore, it is becoming increasingly important to protect privacy information in data publishing. In this paper, we focus on the multi-type self-identified format-preserving encryption. First, we introduce a multi-type self-identified format-preserving encryption system and discuss the encryption of various types of data in this system. Then, for the format preserving encryption(FPE) about Chinese name, we study from the encryption model construction and basic encryption scheme. The format-preserving encryption model about Chinese name is constructed and the concept of the name library is presented. Based on this, it is used to not only limit the message space to reduce complexity, but also ensure the cipher in accordance with the Chinese naming habits. In addition, according to the encryption and decryption model, format-preserving encryption process of Chinese name is designed. In order to add new names, the algorithm of name space expansion is proposed. Based on the Prefix, this paper put forward an algorithm named Cycle-Prefix, which enhances the security and dynamics of FPE by using two adjustment factors and the circular encryption. Compared with the traditional Prefix algorithm, experiments show that Cycle-Prefix can not only complete the task of FPE for Chinese name, but also encrypt same plain text into different ciphers under the premise of similar efficiency with Prefix.展开更多
This approach represents the relative susceptibility of the topography of the earth to active deformation by means of geometrical distinctiveness of the river networks. This investigation employs the fractal analysis ...This approach represents the relative susceptibility of the topography of the earth to active deformation by means of geometrical distinctiveness of the river networks. This investigation employs the fractal analysis of drainage system extracted from ASTER Global Digital Elevation Model (GDEM-30m resolution). The objective is to mark active structures and to pinpoint the areas robustly influenced by neotectonics. This approach was examined in the Hindukush, NE-Afghanistan. This region is frequently affected by deadly earthquakes and the modern fault activities and deformation are driven by the collision between the northward-moving Indian subcontinent and Eurasia. This attempt is based on the fact that drainage system is strained to linearize due to neotectonic deformation. Hence, the low fractal dimensions of the Kabul, Panjsher, Laghman, Andarab, Alingar and Kocha Rivers are credited to active tectonics. A comprehensive textural examination is conducted to probe the linearization, heterogeneity and connectivity of the drainage patterns. The aspects for these natural textures are computed by using the fractal dimension (FD), lacunarity (LA) and succolarity (SA) approach. All these methods are naturally interrelated, i.e. objects with similar FD can be further differentiated with LA and/or SA analysis. The maps of FD, LA and SA values are generated by using a sliding window of 50 arc seconds by 50 arc seconds (50" × 50"). Afterwards, the maps are interpreted in terms of regional susceptibility to neotectonics. This method is useful to pinpoint numerous zones where the drainage system is highly controlled by Hindukush active structures. In the North-Northeast of the Kabul block, we recognized active tectonic blocks. The region comprising, Kabul, Panjsher, Andrab, Alingar and Badakhshan is more susceptible to damaging events. This investigation concludes that the fractal analysis of the river networks is a bonus tool to localize areas vulnerable to deadly incidents influencing the Earth’s topography and consequently intimidate human lives.展开更多
In this paper, the authors used the Princeton Ocean Model (POM) to simulate the seasonal evolu- tions of circulation and thermal structure in the Yellow Sea. The simulated circulation showed that the Yellow Sea Warm C...In this paper, the authors used the Princeton Ocean Model (POM) to simulate the seasonal evolu- tions of circulation and thermal structure in the Yellow Sea. The simulated circulation showed that the Yellow Sea Warm Current (YSWC) was a compensation current of monsoon-driven current, and that in winter, the YSWC became stronger with depth, and could flow across the Bohai Strait in the north. Sensitivity and control- ling tests led to the following conclusions. In winter, the direction of the Yellow Sea Coastal Current in the sur- face layer was controlled partly by tide instead of wind. In summer, a cyclonic horizontal gyre existed in the middle and eastern parts of the Yellow Sea below 10 m. The downwelling in upper layer and upwelling in lower layer were somehow similar to Hu et al. (1991) conceptual model. The calculated thermal structure showed an obvious northward extending YSWC tongue in winter, its position and coverage of the Yellow Sea Cold Water Mass in summer.展开更多
The fight-bank slope of the Dagangshan hydropower station located in Southwest China is a highly unloaded rock slope. Moreover, large-scale natural faults were detected in the slope body; some excavation-induced unloa...The fight-bank slope of the Dagangshan hydropower station located in Southwest China is a highly unloaded rock slope. Moreover, large-scale natural faults were detected in the slope body; some excavation-induced unloading fractures were discovered at elevations between lo75m and 1146m. Because of poor tectonic stability, the excavation work was suspended in September 2009, and six large- scale anti-shear galleries were employed to replace the weak zone in the slope body to reinforce the fight- bank slope. In this study, based on microseismic- monitoring technology and a numerical-simulation method, the stabilities of the slope with and without the reinforcement are analysed. An in-situ microseismic-monitofing system is used to obtain quantitative information about the damage location, extent, energy, and magnitude of the rocks. Thus, any potential sliding block in the fight-bank slope can be identified. By incorporating the numerical results along with the microseismic-monitoring data, the stress concentration is found to largely occur aroundthe anti-shear galleries, and the seismic deformation near the anti-shear galleries is apparent, particularly at elevations of 121o, 118o, 115o, and 112om. To understand the interaction mechanism between the anti-shear gallery and the surrounding rock, a 2D simulation of the potential damage process occurring in an anti-shear gallery is performed. The numerical simulation helps in obtaining additional information about the stress distribution and failure-induced stress re-distribution in the vicinity of the anti-shear galleries that cannot be directly observed in the field. Finally, the potential sliding surface of the right-bank slope is numerically obtained, which generally agrees with the spatial distribution of the in-situ monitored microseismic events. The safety factor of the slope reinforced with the anti-shear gallery increases by approximately 36.2%. Both the numerical results and microseismic data show that the anti-shear galleries have a good reinforcement effect.展开更多
The paper intends to analyze economic factors that influence electricity consumption in the OECD economies. A special interest in this context is given to spillover effects of trade on electricity consumption. For thi...The paper intends to analyze economic factors that influence electricity consumption in the OECD economies. A special interest in this context is given to spillover effects of trade on electricity consumption. For this purpose, a model is constructed that using a dynamic panel study approach. The model is estimated in a GMM framework in which a dynamic procedure is conducted along the balanced growth path for electricity consumption in each economy. In advance, the long run dynamic behavior of prices, GDP, and trade induced spillover variables is determined. In a further step, the short run dynamic mechanism is pursued by estimating the partial adjustment dynamic coefficient on the target level of electricity consumption. The analysis is conducted for industrial, as well as residential electricity consumption. Alternatively, the same procedure is estimated by the application of a fixed period model. The model provides a benchmark tool for electricity policy decisions and for electricity consumption projections.展开更多
In this paper, we revisit foundations of the applications of physical measurement and Lindenmayer system to the modeling of plants. The measurement is proposed to a formal procedure and measuring the mass of leaves on...In this paper, we revisit foundations of the applications of physical measurement and Lindenmayer system to the modeling of plants. The measurement is proposed to a formal procedure and measuring the mass of leaves on a tree, tailored to branching plant structures with Simpson' s rule and Monte Carlo Methods. L-system is possible to visualize mathematical models of biological structures and processes. The formalism is illustrated using theoretical branching systems, and applied to analyze total leaves number as well as total weight of them.展开更多
In order to improve the automatic retrieval ability of English vocabulary, for the distribution of semantic attributes in English vocabulary, an automatic classification method of English vocabulary is proposed based ...In order to improve the automatic retrieval ability of English vocabulary, for the distribution of semantic attributes in English vocabulary, an automatic classification method of English vocabulary is proposed based on association rules, English vocabulary data storage model is constructed, a two element linguistic feature function is constructed for describing the directionality of English lexical retrieval scheduling, English vocabulary classification decision making model is constructed based on contextual relations of English vocabulary, the features of the association rules of English vocabulary are extracted, the adaptive learning method is used to realize the automatic classification of English vocabulary. The simulation results show that the method of English vocabulary classification has good performance, the classification error rate is low, the retrieval precision is high, and the computational overhead is small.展开更多
Flip buckets are commonly used to discharge flow away from a hydraulic structure into the downstream to dissipate energy.A new leak-floor flip bucket is presented,making the ski-jump water jet a typical long-narrow na...Flip buckets are commonly used to discharge flow away from a hydraulic structure into the downstream to dissipate energy.A new leak-floor flip bucket is presented,making the ski-jump water jet a typical long-narrow nappe.Based on the model experiments and numerical simulation,the flow pattern,formation process and mechanism of the leak-floor flip bucket are studied.The results show that cross section flow shape develops from the"Y-type"to"|-type",and this is because the natural pressure difference is generated when water flows through the leak-floor area and moves transversely from both sides to the center.Different from the slit-type flip bucket with sidewall contraction,the leak-floor flit bucket makes the water jet narrow and long without high pressure on the side walls of the flip bucket.Under the same jet length condition,the maximum sidewall pressure of the slit-type is 4.67 times that of the leak-floor flip bucket.The effects of flow discharge on the jet length are less significant for the leak-floor bucket than for the slit-type bucket.展开更多
Shale gas has been discovered in the Upper Triassic Yanchang Formation, Ordos Basin, China. Due to the weak tectonic activities in which the shale plays, core observations indicate abundant random non-tectonic micro- ...Shale gas has been discovered in the Upper Triassic Yanchang Formation, Ordos Basin, China. Due to the weak tectonic activities in which the shale plays, core observations indicate abundant random non-tectonic micro- fractures in the producing shales. The non-tectonic micro-fractures are different from tectonic fractures and are characterized by being irregular, curved, discontinuous, and randomly distributed. The role of micro-fractures in hydraulic fracturing for shale gas development is currently poorly understood yet potentially critical. Two-dimensional computational modeling studies have been used in an initial attempt toward understanding how naturally random fractured reservoirs respond during hydraulic fracturing. The aim of the paper is to investigate the effect of random non-tectonic fractures on hydraulic fracturing. The numerical models with random non-tectonic micro-fractures are built by extracting the fractures of rock blocks after repeated heating and cooling, using a digital image process. Simulations were conducted as a function of:(1) the in-situ stress ratio;(2) internal friction angle of random fractures;(3) cohesion of random fractures;(4) operational variables such as injection rate; and(5) variable injection rate technology. A sensitivity study reveals a number of interesting observations resulting from these parameters on the shear stimulation in a natural fracture system. Three types of fracturing networks were observed from the studied simulations, and the results also show that variable injection rate technology is most promising for producing complex fracturing networks. This work strongly links the production technology and geomechanical evaluation. It can aid in the understanding and optimization of hydraulic fracturing simulations in naturally random fractured reservoirs.展开更多
The present work proposes a novel methodology for constructing coarse-grained (CG) models, which aims at minimizing the difference between CG model and the corresponding original system. The difference is defined as...The present work proposes a novel methodology for constructing coarse-grained (CG) models, which aims at minimizing the difference between CG model and the corresponding original system. The difference is defined as a functional of their equilibrium conformationaJ probability densities, then is estimated from equilibrium averages of many independent physical quantities denoted as basis functions. An orthonormalization strategy is adopted to get the independent basis functions from su^ciently preselected interesting physical quantities of the system. Thus the current method is named as probability density matching coarse-graining (PMCG) scheme, which effectively takes into account the overall cha,~acteristics of the original systems to construct CG model, and it is a rtatural improvement of the usual CG scheme wherein some physical quantities are intuitively chosen without considering their correlations. We verify the general PMCG framework in constructing a one-site CG water model from TIP3P model. Both structure of liquids and pressure of the TIP3P water system are found to be well reproduced at the same time in the constructed CG model.展开更多
Complementing our previous publications, this paper presents the information schema constructs (ISCs) that underpin the programming of specific system manifestation feature (SMF) orientated information management ...Complementing our previous publications, this paper presents the information schema constructs (ISCs) that underpin the programming of specific system manifestation feature (SMF) orientated information management and composing system models. First, we briefly present (1) the general process of pre-embodiment design with SMFs, (2) the procedures of creating genotypes and phenotypes of SMFs, (3) the specific procedure of instantiation of phenotypes of SMFs, and (4) the procedure of system model management and processing. Then, the chunks of information needed for instantiation of phenotypes of SMFs are discussed, and the ISCs designed for instantiation presented. Afterwards, the information management aspects of system modeling are addressed. Methodologically, system modeling involves (1) placement of phenotypes of SMF in the modeling space, (2) combining them towards the desired architecture and operation, (3) assigning values to the parameters and checking the satisfac- tion of constraints, and (4) storing the system model in the SMFs-based warehouse database. The final objective of the reported research is to develop an SMFs-based toolbox to support modeling of cyber-physical systems (CPSs).展开更多
基金supported by the National Natural Science Foundation of China under Grant(No.61772085),(No.61672109),(No.1472024)and.(No.61532012)
文摘In the era of big data, government, business and personal digital information will be possible for data mining. Data mining requires massive data as a support. However, the direct release of the original mass data, which usually contain some sensitive information of personal or analysis, will result in leakage of user privacy. Therefore, it is becoming increasingly important to protect privacy information in data publishing. In this paper, we focus on the multi-type self-identified format-preserving encryption. First, we introduce a multi-type self-identified format-preserving encryption system and discuss the encryption of various types of data in this system. Then, for the format preserving encryption(FPE) about Chinese name, we study from the encryption model construction and basic encryption scheme. The format-preserving encryption model about Chinese name is constructed and the concept of the name library is presented. Based on this, it is used to not only limit the message space to reduce complexity, but also ensure the cipher in accordance with the Chinese naming habits. In addition, according to the encryption and decryption model, format-preserving encryption process of Chinese name is designed. In order to add new names, the algorithm of name space expansion is proposed. Based on the Prefix, this paper put forward an algorithm named Cycle-Prefix, which enhances the security and dynamics of FPE by using two adjustment factors and the circular encryption. Compared with the traditional Prefix algorithm, experiments show that Cycle-Prefix can not only complete the task of FPE for Chinese name, but also encrypt same plain text into different ciphers under the premise of similar efficiency with Prefix.
文摘This approach represents the relative susceptibility of the topography of the earth to active deformation by means of geometrical distinctiveness of the river networks. This investigation employs the fractal analysis of drainage system extracted from ASTER Global Digital Elevation Model (GDEM-30m resolution). The objective is to mark active structures and to pinpoint the areas robustly influenced by neotectonics. This approach was examined in the Hindukush, NE-Afghanistan. This region is frequently affected by deadly earthquakes and the modern fault activities and deformation are driven by the collision between the northward-moving Indian subcontinent and Eurasia. This attempt is based on the fact that drainage system is strained to linearize due to neotectonic deformation. Hence, the low fractal dimensions of the Kabul, Panjsher, Laghman, Andarab, Alingar and Kocha Rivers are credited to active tectonics. A comprehensive textural examination is conducted to probe the linearization, heterogeneity and connectivity of the drainage patterns. The aspects for these natural textures are computed by using the fractal dimension (FD), lacunarity (LA) and succolarity (SA) approach. All these methods are naturally interrelated, i.e. objects with similar FD can be further differentiated with LA and/or SA analysis. The maps of FD, LA and SA values are generated by using a sliding window of 50 arc seconds by 50 arc seconds (50" × 50"). Afterwards, the maps are interpreted in terms of regional susceptibility to neotectonics. This method is useful to pinpoint numerous zones where the drainage system is highly controlled by Hindukush active structures. In the North-Northeast of the Kabul block, we recognized active tectonic blocks. The region comprising, Kabul, Panjsher, Andrab, Alingar and Badakhshan is more susceptible to damaging events. This investigation concludes that the fractal analysis of the river networks is a bonus tool to localize areas vulnerable to deadly incidents influencing the Earth’s topography and consequently intimidate human lives.
基金Contribution No. 4616 from Institute of Oceanology, CAS. Project No.40406025 supported by NSFC.
文摘In this paper, the authors used the Princeton Ocean Model (POM) to simulate the seasonal evolu- tions of circulation and thermal structure in the Yellow Sea. The simulated circulation showed that the Yellow Sea Warm Current (YSWC) was a compensation current of monsoon-driven current, and that in winter, the YSWC became stronger with depth, and could flow across the Bohai Strait in the north. Sensitivity and control- ling tests led to the following conclusions. In winter, the direction of the Yellow Sea Coastal Current in the sur- face layer was controlled partly by tide instead of wind. In summer, a cyclonic horizontal gyre existed in the middle and eastern parts of the Yellow Sea below 10 m. The downwelling in upper layer and upwelling in lower layer were somehow similar to Hu et al. (1991) conceptual model. The calculated thermal structure showed an obvious northward extending YSWC tongue in winter, its position and coverage of the Yellow Sea Cold Water Mass in summer.
基金jointly supported by grants from the National Key Research and Development Program(Grant No.2016YFC0801607,2016YFC0801602)the National Natural Science Foundation of China(Grant No.51279024)the National Basic Research Program of China(Grant No.2014CB047103)
文摘The fight-bank slope of the Dagangshan hydropower station located in Southwest China is a highly unloaded rock slope. Moreover, large-scale natural faults were detected in the slope body; some excavation-induced unloading fractures were discovered at elevations between lo75m and 1146m. Because of poor tectonic stability, the excavation work was suspended in September 2009, and six large- scale anti-shear galleries were employed to replace the weak zone in the slope body to reinforce the fight- bank slope. In this study, based on microseismic- monitoring technology and a numerical-simulation method, the stabilities of the slope with and without the reinforcement are analysed. An in-situ microseismic-monitofing system is used to obtain quantitative information about the damage location, extent, energy, and magnitude of the rocks. Thus, any potential sliding block in the fight-bank slope can be identified. By incorporating the numerical results along with the microseismic-monitoring data, the stress concentration is found to largely occur aroundthe anti-shear galleries, and the seismic deformation near the anti-shear galleries is apparent, particularly at elevations of 121o, 118o, 115o, and 112om. To understand the interaction mechanism between the anti-shear gallery and the surrounding rock, a 2D simulation of the potential damage process occurring in an anti-shear gallery is performed. The numerical simulation helps in obtaining additional information about the stress distribution and failure-induced stress re-distribution in the vicinity of the anti-shear galleries that cannot be directly observed in the field. Finally, the potential sliding surface of the right-bank slope is numerically obtained, which generally agrees with the spatial distribution of the in-situ monitored microseismic events. The safety factor of the slope reinforced with the anti-shear gallery increases by approximately 36.2%. Both the numerical results and microseismic data show that the anti-shear galleries have a good reinforcement effect.
文摘The paper intends to analyze economic factors that influence electricity consumption in the OECD economies. A special interest in this context is given to spillover effects of trade on electricity consumption. For this purpose, a model is constructed that using a dynamic panel study approach. The model is estimated in a GMM framework in which a dynamic procedure is conducted along the balanced growth path for electricity consumption in each economy. In advance, the long run dynamic behavior of prices, GDP, and trade induced spillover variables is determined. In a further step, the short run dynamic mechanism is pursued by estimating the partial adjustment dynamic coefficient on the target level of electricity consumption. The analysis is conducted for industrial, as well as residential electricity consumption. Alternatively, the same procedure is estimated by the application of a fixed period model. The model provides a benchmark tool for electricity policy decisions and for electricity consumption projections.
文摘In this paper, we revisit foundations of the applications of physical measurement and Lindenmayer system to the modeling of plants. The measurement is proposed to a formal procedure and measuring the mass of leaves on a tree, tailored to branching plant structures with Simpson' s rule and Monte Carlo Methods. L-system is possible to visualize mathematical models of biological structures and processes. The formalism is illustrated using theoretical branching systems, and applied to analyze total leaves number as well as total weight of them.
文摘In order to improve the automatic retrieval ability of English vocabulary, for the distribution of semantic attributes in English vocabulary, an automatic classification method of English vocabulary is proposed based on association rules, English vocabulary data storage model is constructed, a two element linguistic feature function is constructed for describing the directionality of English lexical retrieval scheduling, English vocabulary classification decision making model is constructed based on contextual relations of English vocabulary, the features of the association rules of English vocabulary are extracted, the adaptive learning method is used to realize the automatic classification of English vocabulary. The simulation results show that the method of English vocabulary classification has good performance, the classification error rate is low, the retrieval precision is high, and the computational overhead is small.
基金supported by the National Natural Science Foundation of China(Grant Nos.51379138&51179113)Sichuan Province Funds for Distinguished Young Scientists(Grant No.2013JQ0007)
文摘Flip buckets are commonly used to discharge flow away from a hydraulic structure into the downstream to dissipate energy.A new leak-floor flip bucket is presented,making the ski-jump water jet a typical long-narrow nappe.Based on the model experiments and numerical simulation,the flow pattern,formation process and mechanism of the leak-floor flip bucket are studied.The results show that cross section flow shape develops from the"Y-type"to"|-type",and this is because the natural pressure difference is generated when water flows through the leak-floor area and moves transversely from both sides to the center.Different from the slit-type flip bucket with sidewall contraction,the leak-floor flit bucket makes the water jet narrow and long without high pressure on the side walls of the flip bucket.Under the same jet length condition,the maximum sidewall pressure of the slit-type is 4.67 times that of the leak-floor flip bucket.The effects of flow discharge on the jet length are less significant for the leak-floor bucket than for the slit-type bucket.
基金supported by the National Natural Science Foundation of China(Grant Nos.4122790141330643&41502294)+2 种基金China Postdoctoral Science Foundation Funded Project(Grants No.2015M571118)the Strategic Priority Research Program of the Chinese Academy of Sciences(Grants Nos.XDB10030000XDB10030300&XDB10050400)
文摘Shale gas has been discovered in the Upper Triassic Yanchang Formation, Ordos Basin, China. Due to the weak tectonic activities in which the shale plays, core observations indicate abundant random non-tectonic micro- fractures in the producing shales. The non-tectonic micro-fractures are different from tectonic fractures and are characterized by being irregular, curved, discontinuous, and randomly distributed. The role of micro-fractures in hydraulic fracturing for shale gas development is currently poorly understood yet potentially critical. Two-dimensional computational modeling studies have been used in an initial attempt toward understanding how naturally random fractured reservoirs respond during hydraulic fracturing. The aim of the paper is to investigate the effect of random non-tectonic fractures on hydraulic fracturing. The numerical models with random non-tectonic micro-fractures are built by extracting the fractures of rock blocks after repeated heating and cooling, using a digital image process. Simulations were conducted as a function of:(1) the in-situ stress ratio;(2) internal friction angle of random fractures;(3) cohesion of random fractures;(4) operational variables such as injection rate; and(5) variable injection rate technology. A sensitivity study reveals a number of interesting observations resulting from these parameters on the shear stimulation in a natural fracture system. Three types of fracturing networks were observed from the studied simulations, and the results also show that variable injection rate technology is most promising for producing complex fracturing networks. This work strongly links the production technology and geomechanical evaluation. It can aid in the understanding and optimization of hydraulic fracturing simulations in naturally random fractured reservoirs.
基金Supported by National Natural Science Foundation of China under Grant No.11175250
文摘The present work proposes a novel methodology for constructing coarse-grained (CG) models, which aims at minimizing the difference between CG model and the corresponding original system. The difference is defined as a functional of their equilibrium conformationaJ probability densities, then is estimated from equilibrium averages of many independent physical quantities denoted as basis functions. An orthonormalization strategy is adopted to get the independent basis functions from su^ciently preselected interesting physical quantities of the system. Thus the current method is named as probability density matching coarse-graining (PMCG) scheme, which effectively takes into account the overall cha,~acteristics of the original systems to construct CG model, and it is a rtatural improvement of the usual CG scheme wherein some physical quantities are intuitively chosen without considering their correlations. We verify the general PMCG framework in constructing a one-site CG water model from TIP3P model. Both structure of liquids and pressure of the TIP3P water system are found to be well reproduced at the same time in the constructed CG model.
文摘Complementing our previous publications, this paper presents the information schema constructs (ISCs) that underpin the programming of specific system manifestation feature (SMF) orientated information management and composing system models. First, we briefly present (1) the general process of pre-embodiment design with SMFs, (2) the procedures of creating genotypes and phenotypes of SMFs, (3) the specific procedure of instantiation of phenotypes of SMFs, and (4) the procedure of system model management and processing. Then, the chunks of information needed for instantiation of phenotypes of SMFs are discussed, and the ISCs designed for instantiation presented. Afterwards, the information management aspects of system modeling are addressed. Methodologically, system modeling involves (1) placement of phenotypes of SMF in the modeling space, (2) combining them towards the desired architecture and operation, (3) assigning values to the parameters and checking the satisfac- tion of constraints, and (4) storing the system model in the SMFs-based warehouse database. The final objective of the reported research is to develop an SMFs-based toolbox to support modeling of cyber-physical systems (CPSs).