This study proposes a combined hybrid energy storage system(HESS) and transmission grid(TG) model, and a corresponding time series operation simulation(TSOS) model is established to relieve the peak-shaving pressure o...This study proposes a combined hybrid energy storage system(HESS) and transmission grid(TG) model, and a corresponding time series operation simulation(TSOS) model is established to relieve the peak-shaving pressure of power systems under the integration of renewable energy. First, a linear model for the optimal operation of the HESS is established, which considers the different power-efficiency characteristics of the pumped storage system, electrochemical storage system, and a new type of liquid compressed air energy storage. Second, a TSOS simulation model for peak shaving is built to maximize the power entering the grid from the wind farms and HESS. Based on the proposed model, this study considers the transmission capacity of a TG. By adding the power-flow constraints of the TG, a TSOS-based HESS and TG combination model for peak shaving is established. Finally, the improved IEEE-39 and IEEE-118 bus systems were considered as examples to verify the effectiveness and feasibility of the proposed model.展开更多
Calculation grid and turbulence model for numerical simulating pressure fluctuations in a high-speed train tunnel are studied through the comparison analysis of numerical simulation and moving model test.Compared the ...Calculation grid and turbulence model for numerical simulating pressure fluctuations in a high-speed train tunnel are studied through the comparison analysis of numerical simulation and moving model test.Compared the waveforms and peak-peak values of pressure fluctuations between numerical simulation and moving model test,the structured grid and the SST k-ωturbulence model are selected for numerical simulating the process of high-speed train passing through the tunnel.The largest value of pressure wave amplitudes of numerical simulation and moving model test meet each other.And the locations of the largest value of the initial compression and expansion wave amplitude of numerical simulation are in agreement with that of moving model test.The calculated pressure at the measurement point fully conforms to the propagation law of compression and expansion waves in the tunnel.展开更多
Process algorithm, numerical model and techno-economic assessment of charge calculation and furnace bath optimization for target alloy for induction furnace-based steelmaking is presented in this study. The developed ...Process algorithm, numerical model and techno-economic assessment of charge calculation and furnace bath optimization for target alloy for induction furnace-based steelmaking is presented in this study. The developed algorithm combines the make-to-order (MTO) and charge optimization planning (COP) of the steel melting shop in the production of target steel composition. Using a system-level approach, the unit operations involved in the melting process were analyzed with the purpose of initial charge calculation, prevailing alloy charge prediction and optimizing the sequence of melt chemistry modification. The model performance was established using real-time production data from a cast iron-based foundry with a 1- and 2-ton induction furnace capacity and a medium carbon-based foundry with a 10- and 15-ton induction furnace capacity. A simulation engine (CastMELT) was developed in Java IDE with a MySQL database for continuous interaction with changing process parameters to run the model for validation. The comparison between the model prediction and production results was analyzed for charge prediction, melt modification and ferroalloy optimization and possible cost savings. The model performance for elemental charge prediction and calculation purpose with respect to the charge input (at overall scrap meltdown) gave R-squared, Standard Error, Pearson correlation and Significance value of (0.934, 0.06, 0.97, 0.0003) for Carbon prediction, (0.962, 0.06, 0.98, 0.00009) for Silicon prediction, (0.999, 0.048, 0.999, 9E -11) for Manganese Prediction, and (0.997, 0.076, 0.999, 6E -7) for Chromium prediction respectively. Correlation analysis for melt modification (after charging of ferroalloy) using the model for after-alloying spark analysis compared with the target chemistry is at 99.82%. The results validate the suitability of the developed model as a functional system of induction furnace melting for combined charge calculation and melt optimization Techno-economic evaluation results showed that 0.98% - 0.25% ferroalloy saving per ton of melt is possible using the model. This brings about an annual production cost savings of 100,000 $/y in foundry A (medium carbon steel) and 20,000 $/y in foundry B (cast iron) on the use of different ferroalloy materials.展开更多
In this paper, two sub-grid scale (SGS) models are introduced into the Lattice Boltzmann Method (LBM), i.e., the dynamics SGS model and the dynamical system SGS model, and applied to numerically solving three-dimensio...In this paper, two sub-grid scale (SGS) models are introduced into the Lattice Boltzmann Method (LBM), i.e., the dynamics SGS model and the dynamical system SGS model, and applied to numerically solving three-dimensional high Re turbulent cavity flows. Results are compared with those obtained from the Smagorinsky model and direct numerical simulation for the same cases. It is shown that the method with LBM dynamics SGS model has advantages of fast computation speed, suitable to simulate high Re turbulent flows. In addition, it can capture detailed fine structures of turbulent flow fields. The method with LBM dynamical system SGS model dose not contain any adjustable parameters, and can be used in simulations of various complicated turbulent flows to obtain correct information of sub-grid flow field, such as the backscatter of energy transportation between large and small scales. A new average method of eliminating the inherent unphysical oscillation of LBM is also given in the paper.展开更多
Computer grids are infrastructures in which heterogeneous and distributed resources offer very high computing or storage performance. If they offer extreme computing performance, they are also subject to the appearanc...Computer grids are infrastructures in which heterogeneous and distributed resources offer very high computing or storage performance. If they offer extreme computing performance, they are also subject to the appearance of many failures related to this type of architecture. While performing tasks, if the response time of a node in the system incomprehensibly exceeds the requirements of the specifications, the node experiences an omission failure. The task running in the failed node will be unavailable until the node resumes normal activity. Waiting not being a possible solution, many fault tolerance methods have been proposed. Despite this large number of fault tolerance methods on offer, computer grids are still prone to many failures by omission. In this work, a numerical study of the failures by omission which occur in the calculation grids during the execution of the tasks was carried out and a model allowing anticipating its failures was proposed with the formalism PDEVS (Parallel Discret EVent system Specification).展开更多
Biot-flow and squirt-flow are the two most important fluid flow mechanisms in porous media containing fluids. Based on the BISQ (Biot-Squirt) model where the two mechanisms are treated simultaneously, the elastic wa...Biot-flow and squirt-flow are the two most important fluid flow mechanisms in porous media containing fluids. Based on the BISQ (Biot-Squirt) model where the two mechanisms are treated simultaneously, the elastic wave-field simulation in the porous medium is limited to two-dimensions and two-components (2D2C) or two-dimensions and three-components (2D3C). There is no previous report on wave simulation in three- dimensions and three-components. Only through three dimensional numerical simulations can we have an overall understanding of wave field coupling relations and the spatial distribution characteristics between the solid and fluid phases in the dual-phase anisotropic medium. In this paper, based on the BISQ equation, we present elastic wave propagation in a three dimensional dual-phase anisotropic medium simulated by the staggered-grid high-order finite-difference method. We analyze the resulting wave fields and show that the results are an improvement.展开更多
The flows behind the base of a generic rocket, at both hypersonic and subsonic flow conditions, are numerically studied. The main concerns are addressed to the evaluation of turbulence models and the using of grid ada...The flows behind the base of a generic rocket, at both hypersonic and subsonic flow conditions, are numerically studied. The main concerns are addressed to the evaluation of turbulence models and the using of grid adaptation techniques. The investigation focuses on two configurations, related to hypersonic and subsonic experiments. The applicability tests of different turbu- lence models are conducted on the level of two-equation models calculating the steady state solution of the Reynolds-averaged Navier-Stokes(RANS) equations. All used models, the original Wilcox k-co, the Menter shear-stress transport (SST) and the ex- plicit algebraic Reynolds stress model(EARSM) formulation, predict an asymmetric base flow in both cases caused by the support of the models. A comparison with preliminary experimental results indicates a preference for the SST and EARSM results over the results from the older k-co model. Sensitivity studies show no significant influence of the grid topology or the location of the laminar to turbulent transition on the base flow field, but a strong influence of even small angles of attack is reported from the related experiments.展开更多
In order to study the extraction and acceleration mechanism of the dual-stage grid,a three-dimensional model based on the Particle-In-Cell/Monte Carlo Collision(PIC/MCC)method is performed.Dual-stage grid ion thruster...In order to study the extraction and acceleration mechanism of the dual-stage grid,a three-dimensional model based on the Particle-In-Cell/Monte Carlo Collision(PIC/MCC)method is performed.Dual-stage grid ion thruster is a new type of electrostatic ion thruster,which can break through the limitations of traditional gridded ion thrusters,and greatly improve the specific impulse.The high performance also makes the grid sensitive to operating parameters.In this paper,the influence of grid parameters on xenon ion thruster’s performance in a wide range is systematically simulated,and the optimal operating condition is given.Both the over-focusing of the plume,and the transparency of the screen grid are improved,and the grid corrosion is reduced through simulation optimization.The specific impulse under the given working conditions is 9877.24 s and the thrust is 7.28 mN.Based on the simulation optimization,the limitation of the dual-stage grid is discussed.The grid performs well under high voltage conditions(>3000 V)but not well under low voltage conditions(<2000 V).Finally,since argon is cheaper and more advantageous in future engineering applications,the plasma distribution and grid extraction ability under xenon and argon are analyzed and compared to study the flexibility of the dual-stage grid ion thruster.The simulation results show that a set of optimal parameters is only applicable to the corresponding propellant,which needs to be optimized for different propellant types.展开更多
This paper proposes to use the power system simulation software CYME to plan, model and simulate for an actual distribution network for improving the reliability and efficiency, enhancing the efficiency and capacity, ...This paper proposes to use the power system simulation software CYME to plan, model and simulate for an actual distribution network for improving the reliability and efficiency, enhancing the efficiency and capacity, simulating the abnormal condition of distribution network, and presenting operation program of safe, reliable and having simulation record statements. The modeling simulation results show that the software module has lots of advantages including high accuracy, ideal reliability, powerful practicality in simulation and analysis of distribution network, it only need to create once model, the model can sufficiently satisfy multifarious types of simulation analysis required for the distribution network planning.展开更多
Grid direction selection and grid size design are two important elements that need to be considered in the grid direction design in reservoir numerical simulation. Reservoir engineers normally utilize geological data ...Grid direction selection and grid size design are two important elements that need to be considered in the grid direction design in reservoir numerical simulation. Reservoir engineers normally utilize geological data (such as the distribution of fractures, low permeability zones, faults and major stress) and simulation experiences to design the grid direction of simulation model qualitatively. The research of the paper indicates that the key to determine the grid direction is to determine the principal permeability direction. Under the circumstances of few static materials, a new grid direction determination method has been developed by using field data (well location map and inter-well permeability) on the bases of Darcy’s law and tensor analysis theory. The grid direction of WZ11-7 Oilfield simulation model has been determined using four production wells and two production zones (L1 and L3) in WZ11-7-2 well group, the results are in conformity with the geological studied major stress. Therefore, this method can give insights into the numerical simulation study.展开更多
As the smart grid concepts are emphasized lately, the need to modernize the power engineering education is also well recognized. This paper presents a set of newly developed modeling, simulation and testing tools aime...As the smart grid concepts are emphasized lately, the need to modernize the power engineering education is also well recognized. This paper presents a set of newly developed modeling, simulation and testing tools aimed at better understanding of the design concept and related applications for protective relaying and substation automation solutions for the smart grid. Since the smart grid applications require integration of data from multiple IEDs (intelligent electronic devices), understanding properties of each IED type in detail, as well as their responses to the power system events is needed. In addition, understanding the communication requirements to perform data integration is also important. To illustrate how the mentioned smart grid issues may be taught, the following advanced teaching approaches are presented: (1) Use of modeling and simulation means to better understand interaction between the relays and power system; (2) Use of IED test facilities to better understand performance of physical devices used for protection, monitoring and control; (3) Utilization of communication network modeling tools to simulate the communication network within SAS (substation automation system). Examples showing the use of proposed techniques for teaching the fundamentals and applications are presented. The examples demonstrate the adequacy and efficiency of the proposed techniques.展开更多
The interconnection of urban critical infrastructure poses new challenges to the secure operation of power grid‐centric urban lifeline systems.The interdependencies among infrastructure systems increase the risk of c...The interconnection of urban critical infrastructure poses new challenges to the secure operation of power grid‐centric urban lifeline systems.The interdependencies among infrastructure systems increase the risk of cascading fault propagation,thereby threatening urban public safety.Coupling modelling and fault propagation simulation of urban lifeline systems provide a foundation for analysing disaster evolution mechanisms and support research on resilience assessment and enhancement technologies.In this paper,a coupling modelling and simulation method for urban lifeline systems under extreme disaster scenarios is proposed.First,based on fault propagation analysis requirements,the basic principles for coupling modelling are established,and simulation models suitable for lifeline systems under different time scales are analysed.Second,a co-simulation method combining time-driven and event-driven mechanisms is designed.Customised interfaces for information exchange between systems are developed based on mature simulation software,enabling the simulation of the disaster impact propagation process.Finally,the effectiveness of the proposed method is verified using a constructed urban lifeline coupling test system.展开更多
In this paper, the numerical model of the net cage with the grid mooring system in waves is set up by the lumped mass method and rigid kinematics theory, and then the motion equations of floating system, net system, m...In this paper, the numerical model of the net cage with the grid mooring system in waves is set up by the lumped mass method and rigid kinematics theory, and then the motion equations of floating system, net system, mooring system, and floaters are solved by the Runge-Kutta fifth-order method. For the verification of the numerical model, a series of physical model tests have been carried out. According to the comparisons between the simulated and experimental results, it can be found that the simulated and experimental results agree well in each condition. Then, the effects of submerged depth of grid and direction of incident wave propagation on hydrodynamic behaviors of the net cage are analyzed. According to the simulated results, it can be found that with the increase of submerged depth of grid, the forces acting on mooring lines and bridle lines increase, while the forces on grid lines decrease; the horizontal motion amplitudes of floating collar decrease obviously, while the vertical motion amplitudes of floating collar change little. When the direction of incident wave propagation changes, forces on mooring lines and motion of net cage also change accordingly. When the propagation direction of incident wave changes from 0° to 45°, forces on the main ropes and bridle ropes increase, while the forces on the grid ropes decrease. With the increasing propagation direction of incident wave, the horizontal amplitude of the forces collar decreases, while the vertical amplitude of the floating collar has little variation.展开更多
The injection of fuel-generated CO2 into oil reservoirs will lead to two benefits in both enhanced oil recovery (EOR) and the reduction in atmospheric emission of CO2. To get an insight into CO2 miscible flooding pe...The injection of fuel-generated CO2 into oil reservoirs will lead to two benefits in both enhanced oil recovery (EOR) and the reduction in atmospheric emission of CO2. To get an insight into CO2 miscible flooding performance in oil reservoirs, a multi-compositional non-isothermal CO2 miscible flooding mathematical model is developed. The convection and diffusion of CO2-hydrocarbon mixtures in multiphase fluids in reservoirs, mass transfer between CO2 and crude, and formation damages caused by asphaltene precipitation are fully considered in the model. The governing equations are discretized in space using the integral finite difference method. The Newton-Raphson iterative technique was used to solve the nonlinear equation systems of mass and energy conservation. A numerical simulator, in which regular grids and irregular grids are optional, was developed for predicting CO2 miscible flooding processes. Two examples of one-dimensional (1D) regular and three-dimensional (3D) rectangle and polygonal grids are designed to demonstrate the functions of the simulator. Experimental data validate the developed simulator by comparison with 1D simulation results. The applications of the simulator indicate that it is feasible for predicting CO2 flooding in oil reservoirs for EOR.展开更多
Based on twin-roll casting, a cast-rolling force model was proposed to predict the rolling force in the bimetal solid-liquid cast-rolling bonding(SLCRB) process. The solid-liquid bonding zone was assumed to be below t...Based on twin-roll casting, a cast-rolling force model was proposed to predict the rolling force in the bimetal solid-liquid cast-rolling bonding(SLCRB) process. The solid-liquid bonding zone was assumed to be below the kiss point(KP). The deformation resistance of the liquid zone was ignored. Then, the calculation model was derived. A 2D thermal-flow coupled simulation was established to provide a basis for the parameters in the model, and then the rolling forces of the Cu/Al clad strip at different rolling speeds were calculated. Meanwhile, through measurement experiments, the accuracy of the model was verified. The influence of the rolling speed, the substrate strip thickness, and the material on the rolling force was obtained. The results indicate that the rolling force decreases with the increase of the rolling speed and increases with the increase of the thickness and thermal conductivity of the substrate strip. The rolling force is closely related to the KP height. Therefore, the formulation of reasonable process parameters to control the KP height is of great significance to the stability of cast-rolling forming.展开更多
The characteristics of three-dimensional (3-D) tidal current in the Oujiang Estuary are investigated according to in situ observations. The Oujiang Estuary has features of irregular coastline, complex topography, ma...The characteristics of three-dimensional (3-D) tidal current in the Oujiang Estuary are investigated according to in situ observations. The Oujiang Estuary has features of irregular coastline, complex topography, many islands, moveable boundary, and submerged dyke, therefore, σ 3-D numerical model oil an unstructured triangular grid has been degeloped. The σ coordinate transforination, the moveable boundary and submerged dyke treatment techniques were employed in the model so it is suitable for the tidal simulations in the Oujing Estuary with submerged dyke and moveable boundary problems. The model is evaluated with the in situ data, and the results show that the calculated water elevations at 19 stations and currents at 19 profiler stations are in good agreement with measured data both in magnitude and phase. This numerical model is applied to the 3-D tidal circulation simulations of experiments in stopping flow transport through the South Branch of the Oujiang Estuary, and the feasibility to cutoff the flow in the South Branch of the Oujiang Estuary is demonstrated by numerical simulation experiments. The developed numerical model simulated the 3-D tidal current circulations in complicated coastal and estuarine waters very well.展开更多
The development of new engineering alloy chemistries and heat treatments is a time-consuming and iterative process.Here,a hybrid approach of the high-throughput precipitation simulations and decisive experiments is de...The development of new engineering alloy chemistries and heat treatments is a time-consuming and iterative process.Here,a hybrid approach of the high-throughput precipitation simulations and decisive experiments is developed to optimize the composition and manipulate the microstructure of Al-Zn-Mg-Cu alloys to achieve the expected yield strength and elongation.For that purpose,a multi-class Kampmann-Wagner numerical(KWN)framework is established and the contributions to precipitation kinetics and strength from primary phases and precipitates formed before age hardening are introduced for the first time.The composition/process-structure-property relationship of Al-Zn-Mg-Cu alloys is pre-sented and discussed in detail.Coupled with thermodynamic calculations,two concentration-optimized Al-Zn-Mg-Cu alloys with expected high yield strength and long elongation are designed,prepared,and characterized.The excellent strength and elongation of the designed alloys and the good agreement between the measured and model-predicted mechanical properties for these two alloys underscores the remarkable predictive power of the presently developed material design strategy.This work establishes a novel material design strategy for rapidly exploring the compositional space and investigating the effects of composition and heat treatment on the microstructure and performance of ultrahigh strength Al alloys and other materials.展开更多
基金supported by the State Grid Science and Technology Project (No.52999821N004)。
文摘This study proposes a combined hybrid energy storage system(HESS) and transmission grid(TG) model, and a corresponding time series operation simulation(TSOS) model is established to relieve the peak-shaving pressure of power systems under the integration of renewable energy. First, a linear model for the optimal operation of the HESS is established, which considers the different power-efficiency characteristics of the pumped storage system, electrochemical storage system, and a new type of liquid compressed air energy storage. Second, a TSOS simulation model for peak shaving is built to maximize the power entering the grid from the wind farms and HESS. Based on the proposed model, this study considers the transmission capacity of a TG. By adding the power-flow constraints of the TG, a TSOS-based HESS and TG combination model for peak shaving is established. Finally, the improved IEEE-39 and IEEE-118 bus systems were considered as examples to verify the effectiveness and feasibility of the proposed model.
文摘Calculation grid and turbulence model for numerical simulating pressure fluctuations in a high-speed train tunnel are studied through the comparison analysis of numerical simulation and moving model test.Compared the waveforms and peak-peak values of pressure fluctuations between numerical simulation and moving model test,the structured grid and the SST k-ωturbulence model are selected for numerical simulating the process of high-speed train passing through the tunnel.The largest value of pressure wave amplitudes of numerical simulation and moving model test meet each other.And the locations of the largest value of the initial compression and expansion wave amplitude of numerical simulation are in agreement with that of moving model test.The calculated pressure at the measurement point fully conforms to the propagation law of compression and expansion waves in the tunnel.
文摘Process algorithm, numerical model and techno-economic assessment of charge calculation and furnace bath optimization for target alloy for induction furnace-based steelmaking is presented in this study. The developed algorithm combines the make-to-order (MTO) and charge optimization planning (COP) of the steel melting shop in the production of target steel composition. Using a system-level approach, the unit operations involved in the melting process were analyzed with the purpose of initial charge calculation, prevailing alloy charge prediction and optimizing the sequence of melt chemistry modification. The model performance was established using real-time production data from a cast iron-based foundry with a 1- and 2-ton induction furnace capacity and a medium carbon-based foundry with a 10- and 15-ton induction furnace capacity. A simulation engine (CastMELT) was developed in Java IDE with a MySQL database for continuous interaction with changing process parameters to run the model for validation. The comparison between the model prediction and production results was analyzed for charge prediction, melt modification and ferroalloy optimization and possible cost savings. The model performance for elemental charge prediction and calculation purpose with respect to the charge input (at overall scrap meltdown) gave R-squared, Standard Error, Pearson correlation and Significance value of (0.934, 0.06, 0.97, 0.0003) for Carbon prediction, (0.962, 0.06, 0.98, 0.00009) for Silicon prediction, (0.999, 0.048, 0.999, 9E -11) for Manganese Prediction, and (0.997, 0.076, 0.999, 6E -7) for Chromium prediction respectively. Correlation analysis for melt modification (after charging of ferroalloy) using the model for after-alloying spark analysis compared with the target chemistry is at 99.82%. The results validate the suitability of the developed model as a functional system of induction furnace melting for combined charge calculation and melt optimization Techno-economic evaluation results showed that 0.98% - 0.25% ferroalloy saving per ton of melt is possible using the model. This brings about an annual production cost savings of 100,000 $/y in foundry A (medium carbon steel) and 20,000 $/y in foundry B (cast iron) on the use of different ferroalloy materials.
基金Supported by the Key Project of National Natural Science Foundation of China (Grant No. 10532030)
文摘In this paper, two sub-grid scale (SGS) models are introduced into the Lattice Boltzmann Method (LBM), i.e., the dynamics SGS model and the dynamical system SGS model, and applied to numerically solving three-dimensional high Re turbulent cavity flows. Results are compared with those obtained from the Smagorinsky model and direct numerical simulation for the same cases. It is shown that the method with LBM dynamics SGS model has advantages of fast computation speed, suitable to simulate high Re turbulent flows. In addition, it can capture detailed fine structures of turbulent flow fields. The method with LBM dynamical system SGS model dose not contain any adjustable parameters, and can be used in simulations of various complicated turbulent flows to obtain correct information of sub-grid flow field, such as the backscatter of energy transportation between large and small scales. A new average method of eliminating the inherent unphysical oscillation of LBM is also given in the paper.
文摘Computer grids are infrastructures in which heterogeneous and distributed resources offer very high computing or storage performance. If they offer extreme computing performance, they are also subject to the appearance of many failures related to this type of architecture. While performing tasks, if the response time of a node in the system incomprehensibly exceeds the requirements of the specifications, the node experiences an omission failure. The task running in the failed node will be unavailable until the node resumes normal activity. Waiting not being a possible solution, many fault tolerance methods have been proposed. Despite this large number of fault tolerance methods on offer, computer grids are still prone to many failures by omission. In this work, a numerical study of the failures by omission which occur in the calculation grids during the execution of the tasks was carried out and a model allowing anticipating its failures was proposed with the formalism PDEVS (Parallel Discret EVent system Specification).
基金National Natural Science Foundation (Project number 40604013).
文摘Biot-flow and squirt-flow are the two most important fluid flow mechanisms in porous media containing fluids. Based on the BISQ (Biot-Squirt) model where the two mechanisms are treated simultaneously, the elastic wave-field simulation in the porous medium is limited to two-dimensions and two-components (2D2C) or two-dimensions and three-components (2D3C). There is no previous report on wave simulation in three- dimensions and three-components. Only through three dimensional numerical simulations can we have an overall understanding of wave field coupling relations and the spatial distribution characteristics between the solid and fluid phases in the dual-phase anisotropic medium. In this paper, based on the BISQ equation, we present elastic wave propagation in a three dimensional dual-phase anisotropic medium simulated by the staggered-grid high-order finite-difference method. We analyze the resulting wave fields and show that the results are an improvement.
基金German Research Foundation (Deutsche Forschungs-gemeinschaft-DFG) Sonderforschungsbereich Transregio 40
文摘The flows behind the base of a generic rocket, at both hypersonic and subsonic flow conditions, are numerically studied. The main concerns are addressed to the evaluation of turbulence models and the using of grid adaptation techniques. The investigation focuses on two configurations, related to hypersonic and subsonic experiments. The applicability tests of different turbu- lence models are conducted on the level of two-equation models calculating the steady state solution of the Reynolds-averaged Navier-Stokes(RANS) equations. All used models, the original Wilcox k-co, the Menter shear-stress transport (SST) and the ex- plicit algebraic Reynolds stress model(EARSM) formulation, predict an asymmetric base flow in both cases caused by the support of the models. A comparison with preliminary experimental results indicates a preference for the SST and EARSM results over the results from the older k-co model. Sensitivity studies show no significant influence of the grid topology or the location of the laminar to turbulent transition on the base flow field, but a strong influence of even small angles of attack is reported from the related experiments.
基金co-supported by the National Key R & D Program for Intergovernmental International Scientific and Technological Innovation Cooperation,China(No. 2021YFE0116000)the National Natural Science Foundation of China (Nos. 12175032, 12102082 and 12275044),the National Natural Science Foundation of China+5 种基金the Belarusian Republican Foundation for Fundamental Research (No. 12211530449)the Fundamental Research Funds for the Central Universities of China (No. DUT22QN232)the S&T Program of Hebei, China (Nos. YCYZ202201 and 216Z1901G)the S&T Innovation Program of Hebei, China (Nos. SJMYF2022X18 and SJMYF2022X06)the Science and Technology Project of Hebei Education Department, China (No. ZC2023144)S&T Program of Langfang, China (No. 2022011039)
文摘In order to study the extraction and acceleration mechanism of the dual-stage grid,a three-dimensional model based on the Particle-In-Cell/Monte Carlo Collision(PIC/MCC)method is performed.Dual-stage grid ion thruster is a new type of electrostatic ion thruster,which can break through the limitations of traditional gridded ion thrusters,and greatly improve the specific impulse.The high performance also makes the grid sensitive to operating parameters.In this paper,the influence of grid parameters on xenon ion thruster’s performance in a wide range is systematically simulated,and the optimal operating condition is given.Both the over-focusing of the plume,and the transparency of the screen grid are improved,and the grid corrosion is reduced through simulation optimization.The specific impulse under the given working conditions is 9877.24 s and the thrust is 7.28 mN.Based on the simulation optimization,the limitation of the dual-stage grid is discussed.The grid performs well under high voltage conditions(>3000 V)but not well under low voltage conditions(<2000 V).Finally,since argon is cheaper and more advantageous in future engineering applications,the plasma distribution and grid extraction ability under xenon and argon are analyzed and compared to study the flexibility of the dual-stage grid ion thruster.The simulation results show that a set of optimal parameters is only applicable to the corresponding propellant,which needs to be optimized for different propellant types.
文摘This paper proposes to use the power system simulation software CYME to plan, model and simulate for an actual distribution network for improving the reliability and efficiency, enhancing the efficiency and capacity, simulating the abnormal condition of distribution network, and presenting operation program of safe, reliable and having simulation record statements. The modeling simulation results show that the software module has lots of advantages including high accuracy, ideal reliability, powerful practicality in simulation and analysis of distribution network, it only need to create once model, the model can sufficiently satisfy multifarious types of simulation analysis required for the distribution network planning.
文摘Grid direction selection and grid size design are two important elements that need to be considered in the grid direction design in reservoir numerical simulation. Reservoir engineers normally utilize geological data (such as the distribution of fractures, low permeability zones, faults and major stress) and simulation experiences to design the grid direction of simulation model qualitatively. The research of the paper indicates that the key to determine the grid direction is to determine the principal permeability direction. Under the circumstances of few static materials, a new grid direction determination method has been developed by using field data (well location map and inter-well permeability) on the bases of Darcy’s law and tensor analysis theory. The grid direction of WZ11-7 Oilfield simulation model has been determined using four production wells and two production zones (L1 and L3) in WZ11-7-2 well group, the results are in conformity with the geological studied major stress. Therefore, this method can give insights into the numerical simulation study.
文摘As the smart grid concepts are emphasized lately, the need to modernize the power engineering education is also well recognized. This paper presents a set of newly developed modeling, simulation and testing tools aimed at better understanding of the design concept and related applications for protective relaying and substation automation solutions for the smart grid. Since the smart grid applications require integration of data from multiple IEDs (intelligent electronic devices), understanding properties of each IED type in detail, as well as their responses to the power system events is needed. In addition, understanding the communication requirements to perform data integration is also important. To illustrate how the mentioned smart grid issues may be taught, the following advanced teaching approaches are presented: (1) Use of modeling and simulation means to better understand interaction between the relays and power system; (2) Use of IED test facilities to better understand performance of physical devices used for protection, monitoring and control; (3) Utilization of communication network modeling tools to simulate the communication network within SAS (substation automation system). Examples showing the use of proposed techniques for teaching the fundamentals and applications are presented. The examples demonstrate the adequacy and efficiency of the proposed techniques.
基金supported by the National Key Research and Development Program of China,Grant:2022YFB2405500The Science and Technology Project of State Grid Corporation of China,Grant:52094023001HShanghai Electric Power Research Institute,State Grid Corporation of China,Grant:52094023003J.
文摘The interconnection of urban critical infrastructure poses new challenges to the secure operation of power grid‐centric urban lifeline systems.The interdependencies among infrastructure systems increase the risk of cascading fault propagation,thereby threatening urban public safety.Coupling modelling and fault propagation simulation of urban lifeline systems provide a foundation for analysing disaster evolution mechanisms and support research on resilience assessment and enhancement technologies.In this paper,a coupling modelling and simulation method for urban lifeline systems under extreme disaster scenarios is proposed.First,based on fault propagation analysis requirements,the basic principles for coupling modelling are established,and simulation models suitable for lifeline systems under different time scales are analysed.Second,a co-simulation method combining time-driven and event-driven mechanisms is designed.Customised interfaces for information exchange between systems are developed based on mature simulation software,enabling the simulation of the disaster impact propagation process.Finally,the effectiveness of the proposed method is verified using a constructed urban lifeline coupling test system.
基金supported by the National Natural Science Foundation of China(Grant No.50809014)the National High Technology Research and Development Program of China(863 Program,Grant No.2006AA100301)the Specialized Research Fund for the Doctoral Program of Higher Education(Grant No.200801411094)
文摘In this paper, the numerical model of the net cage with the grid mooring system in waves is set up by the lumped mass method and rigid kinematics theory, and then the motion equations of floating system, net system, mooring system, and floaters are solved by the Runge-Kutta fifth-order method. For the verification of the numerical model, a series of physical model tests have been carried out. According to the comparisons between the simulated and experimental results, it can be found that the simulated and experimental results agree well in each condition. Then, the effects of submerged depth of grid and direction of incident wave propagation on hydrodynamic behaviors of the net cage are analyzed. According to the simulated results, it can be found that with the increase of submerged depth of grid, the forces acting on mooring lines and bridle lines increase, while the forces on grid lines decrease; the horizontal motion amplitudes of floating collar decrease obviously, while the vertical motion amplitudes of floating collar change little. When the direction of incident wave propagation changes, forces on mooring lines and motion of net cage also change accordingly. When the propagation direction of incident wave changes from 0° to 45°, forces on the main ropes and bridle ropes increase, while the forces on the grid ropes decrease. With the increasing propagation direction of incident wave, the horizontal amplitude of the forces collar decreases, while the vertical amplitude of the floating collar has little variation.
基金Parts of this work were supported by the National Science and Technology Major Projects (2011ZX05009-002, 2011ZX05009–006)the Fundamental Research Funds for the Central Universities, the Project-sponsored by SRF for ROCS, SEM, and the joint research on "Investigation of Mathematical Models and Their Applications for Oil, Water and CO2 Flow in Reservoirs" between Colorado School of Mines, U.S.A and PetroChina Research Institute of Petroleum Exploration & Development (RIPED), CNPC, China
文摘The injection of fuel-generated CO2 into oil reservoirs will lead to two benefits in both enhanced oil recovery (EOR) and the reduction in atmospheric emission of CO2. To get an insight into CO2 miscible flooding performance in oil reservoirs, a multi-compositional non-isothermal CO2 miscible flooding mathematical model is developed. The convection and diffusion of CO2-hydrocarbon mixtures in multiphase fluids in reservoirs, mass transfer between CO2 and crude, and formation damages caused by asphaltene precipitation are fully considered in the model. The governing equations are discretized in space using the integral finite difference method. The Newton-Raphson iterative technique was used to solve the nonlinear equation systems of mass and energy conservation. A numerical simulator, in which regular grids and irregular grids are optional, was developed for predicting CO2 miscible flooding processes. Two examples of one-dimensional (1D) regular and three-dimensional (3D) rectangle and polygonal grids are designed to demonstrate the functions of the simulator. Experimental data validate the developed simulator by comparison with 1D simulation results. The applications of the simulator indicate that it is feasible for predicting CO2 flooding in oil reservoirs for EOR.
基金The authors are grateful for the financial supports from the National Natural Science Foundation of China(51974278)the Distinguished Young Fund of Natural Science Foundation of Hebei Province,China(E2018203446).
文摘Based on twin-roll casting, a cast-rolling force model was proposed to predict the rolling force in the bimetal solid-liquid cast-rolling bonding(SLCRB) process. The solid-liquid bonding zone was assumed to be below the kiss point(KP). The deformation resistance of the liquid zone was ignored. Then, the calculation model was derived. A 2D thermal-flow coupled simulation was established to provide a basis for the parameters in the model, and then the rolling forces of the Cu/Al clad strip at different rolling speeds were calculated. Meanwhile, through measurement experiments, the accuracy of the model was verified. The influence of the rolling speed, the substrate strip thickness, and the material on the rolling force was obtained. The results indicate that the rolling force decreases with the increase of the rolling speed and increases with the increase of the thickness and thermal conductivity of the substrate strip. The rolling force is closely related to the KP height. Therefore, the formulation of reasonable process parameters to control the KP height is of great significance to the stability of cast-rolling forming.
基金The Natural Science Foundation of Tianjin, China under contract No.08JCZDZT00200
文摘The characteristics of three-dimensional (3-D) tidal current in the Oujiang Estuary are investigated according to in situ observations. The Oujiang Estuary has features of irregular coastline, complex topography, many islands, moveable boundary, and submerged dyke, therefore, σ 3-D numerical model oil an unstructured triangular grid has been degeloped. The σ coordinate transforination, the moveable boundary and submerged dyke treatment techniques were employed in the model so it is suitable for the tidal simulations in the Oujing Estuary with submerged dyke and moveable boundary problems. The model is evaluated with the in situ data, and the results show that the calculated water elevations at 19 stations and currents at 19 profiler stations are in good agreement with measured data both in magnitude and phase. This numerical model is applied to the 3-D tidal circulation simulations of experiments in stopping flow transport through the South Branch of the Oujiang Estuary, and the feasibility to cutoff the flow in the South Branch of the Oujiang Estuary is demonstrated by numerical simulation experiments. The developed numerical model simulated the 3-D tidal current circulations in complicated coastal and estuarine waters very well.
基金supported by the National Key Research and Development Program of China(No.2018YFB0704003)the Funds for International Cooperation and Exchange of the National Natural Science Foundation of China(No.51820105001).
文摘The development of new engineering alloy chemistries and heat treatments is a time-consuming and iterative process.Here,a hybrid approach of the high-throughput precipitation simulations and decisive experiments is developed to optimize the composition and manipulate the microstructure of Al-Zn-Mg-Cu alloys to achieve the expected yield strength and elongation.For that purpose,a multi-class Kampmann-Wagner numerical(KWN)framework is established and the contributions to precipitation kinetics and strength from primary phases and precipitates formed before age hardening are introduced for the first time.The composition/process-structure-property relationship of Al-Zn-Mg-Cu alloys is pre-sented and discussed in detail.Coupled with thermodynamic calculations,two concentration-optimized Al-Zn-Mg-Cu alloys with expected high yield strength and long elongation are designed,prepared,and characterized.The excellent strength and elongation of the designed alloys and the good agreement between the measured and model-predicted mechanical properties for these two alloys underscores the remarkable predictive power of the presently developed material design strategy.This work establishes a novel material design strategy for rapidly exploring the compositional space and investigating the effects of composition and heat treatment on the microstructure and performance of ultrahigh strength Al alloys and other materials.