It has been evident that the theory and methods of dynamic derivatives are playing an increasingly important role in hybrid modeling and computations. Being constructed on various kinds of hybrid grids, that is, tim...It has been evident that the theory and methods of dynamic derivatives are playing an increasingly important role in hybrid modeling and computations. Being constructed on various kinds of hybrid grids, that is, time scales, dynamic derivatives offer superior accuracy and flexibility in approximating mathematically important natural processes with hard-to-predict singularities, such as the epidemic growth with unpredictable jump sizes and option market changes with high uncertainties, as compared with conventional derivatives. In this article, we shall review the novel new concepts, explore delicate relations between the most frequently used second-order dynamic derivatives and conventional derivatives. We shall investigate necessary conditions for guaranteeing the consistency between the two derivatives. We will show that such a consistency may never exist in general. This implies that the dynamic derivatives provide entirely different new tools for sensitive modeling and approximations on hybrid grids. Rigorous error analysis will be given via asymptotic expansions for further modeling and computational applications. Numerical experiments will also be given.展开更多
An unstructured Reynolds-averaged Navier-Stokes flow solver using the finite volume method is studied. The spatial discretisation is based on the Osher approximate Riemann solvers. A two-equation turbulence model (k-...An unstructured Reynolds-averaged Navier-Stokes flow solver using the finite volume method is studied. The spatial discretisation is based on the Osher approximate Riemann solvers. A two-equation turbulence model (k-ω model) is also developed for hybrid grids to compute the turbulence flow. The turbulence flow past NACA0012 airfoil and the double ellipsolids are computed, and the numerical results show that the above methods are very efficient.展开更多
Viscous CFD computations are performed to predict the aerodynamic coefficients and flowfield for a missile with grid fins by using hybrid grids. The computations arc made at Math numbvrs of 0.7 and 2.0. Full N-S equat...Viscous CFD computations are performed to predict the aerodynamic coefficients and flowfield for a missile with grid fins by using hybrid grids. The computations arc made at Math numbvrs of 0.7 and 2.0. Full N-S equations arc discretized into finite volume form and solved by an algorithm in LU-SGS. The comparisons between computation and experimental data are made, and the detailed flow structure near grid fin is shown and examined.展开更多
The influence of chemical nonequilibrium on the thermal characteristics is explored by using the 2Dhybrid grid direct simulation Monte Carlo(DSMC)parallel method.An improved molecule search algorithm is proposed,which...The influence of chemical nonequilibrium on the thermal characteristics is explored by using the 2Dhybrid grid direct simulation Monte Carlo(DSMC)parallel method.An improved molecule search algorithm is proposed,which can preserve the high efficiency of area search algorithm.This method can overcome the defects of area search algorithm,and give all information about molecules hitting surface.The heat flux calculation method for a rarefied hypersonic flow is established.In addition,the testing methods of chemical reaction probability for five species of mixed gas with limited speed chemical reactions are also selected.To validate the effectiveness of the present method,hypersonic flow around a cylinder is firstly simulated,and subsequently numerical simulations of the heat flux and flow field characteristics around the blunt body at different heights are carried out in two different cases:the thermal nonequilibrium condition and the thermochemical nonequilibrium condition.Numerical results demonstrate the validity and reliability of the proposed methods.展开更多
An efficient compressible Euler equation solver for vortex-dominated flows is presented based on the adaptive hybrid Cartesian mesh and vortex identifying method.For most traditional grid-based Euler solvers,the exces...An efficient compressible Euler equation solver for vortex-dominated flows is presented based on the adaptive hybrid Cartesian mesh and vortex identifying method.For most traditional grid-based Euler solvers,the excessive numerical dissipation is the great obstruction for vortex capturing or tracking problems.A vortex identifying method based on the curl of velocity is used to identify the vortex in flow field.Moreover,a dynamic adaptive mesh refinement(DAMR)process for hybrid Cartesian gird system is employed to track and preserve vortex.To validate the proposed method,a single compressible vortex convection flow is involved to test the accuracy and efficiency of DAMR process.Additionally,the vortex-dominated flow is investigated by the method.The obtained results are shown as a good agreement with the previous published data.展开更多
A hybrid Cartesian structured grid method is proposed for solving moving boundary unsteady problems. The near body region is discretized by using the body-fitted structured grids, while the remaining computational dom...A hybrid Cartesian structured grid method is proposed for solving moving boundary unsteady problems. The near body region is discretized by using the body-fitted structured grids, while the remaining computational domain is tessellated with the generated Cartesian grids. As the body moves, the structured grids move with the body and the outer boundaries of inside grids are used to generate new holes in the outside adaptive Cartesian grid to facilitate data communication. By using the alternating digital tree （ADT） algorithm, the computational time of hole-cutting and identification of donor cells can be reduced significantly. A compressible solver for unsteady flow problems is developed. A cell-centered, second-order accurate finite volume method is employed in spatial discreti- zation and an implicit dual-time stepping low-upper symmetric Gauss-Seidei （LU-SGS） approach is employed in temporal discretization. Geometry-based adaptation is used during unsteady simulation time steps when boundary moves and the flow solution is interpolated from the old Cartesian grids to the new one with inverse distance weigh- ting interpolation formula. Both laminar and turbulent unsteady cases are tested to demonstrate the accuracy and efficiency of the proposed method. Then, a 2-D store separation problem is simulated. The result shows that the hybrid Cartesian grid method can handle the unsteady flow problems involving large-scale moving boundaries.展开更多
Strong DC coupling with weak AC and large-scale renewable energy integration are the two significant characteristics of ultra-high-voltage AC/DC(UHVAC/DC)hybrid power grids in China.Strong coupling between AC and DC g...Strong DC coupling with weak AC and large-scale renewable energy integration are the two significant characteristics of ultra-high-voltage AC/DC(UHVAC/DC)hybrid power grids in China.Strong coupling between AC and DC grids and the different integration performance of renewable energy sources have profoundly changed the stability characteristics of the power system.The traditional stability control system is inadequate for the stability control of UHVAC/DC power grids.This paper analyzes the requirements for constructing an integrated defense system in a UHVAC/DC hybrid power grid(i.e.power system protection).The definition,connotation,and designing principles of power system protection are put forward.The relationship between the power system protection and the traditional three-defense lines is investigated.The design principles,general hardware structure and main functions of a power system protection are presented.Key problems and technologies are specified in the construction of the power system protection.展开更多
A projectile with exotic wraparound wings( WAW) configuration is designed to improve the finstabilized projectile shooting quality. Two fin-stabilized projectiles with the same body with and without exotic WAW configu...A projectile with exotic wraparound wings( WAW) configuration is designed to improve the finstabilized projectile shooting quality. Two fin-stabilized projectiles with the same body with and without exotic WAW configuration are simulated numerically by applying the Roe scheme. The shear-stress transport turbulence models and the lower-upper symmetric Gauss-Seidel implicit method are used to solve 3D Reynoldsaveraged Navier—Stokes equations. The differences in aerodynamic coefficients and aerodynamic characteristics of the projectiles when the Mach number varies from 0. 35 to 0. 95 are obtained,and the cause of these differences is analyzed. The calculation results indicate that the lift-to-drag ratio of the projectile significantly increases,the rolling moment decreases,and the position of the pressure center of the projectile shows relatively small changes when the exotic WAW configuration is used. Therefore,this projectile can obviously reduce rolling effect,enlarge range and improve flying stability.展开更多
An efficient MPI/OpenMP hybrid parallel Radial Basis Function (RBF) strategy for both continuous and discontinuous large-scale mesh deformation is proposed to reduce the computational cost and memory consumption.Unlik...An efficient MPI/OpenMP hybrid parallel Radial Basis Function (RBF) strategy for both continuous and discontinuous large-scale mesh deformation is proposed to reduce the computational cost and memory consumption.Unlike the conventional parallel methods in which all processors use the same surface displacement and implement the same operation,the present method employs different surface points sets and influence radius for each volume point movement,accompanied with efficient geometry searching strategy.The deformed surface points,also called Control Points (CPs),are stored in each processor.The displacement of spatial points is interpolated by using only 20-50 nearest control points,and the local influence radius is set to 5-20 times the maximum displacement of control points.To shorten the searching time for the nearest control point clouds,an Alternating Digital Tree (ADT) algorithm for 3D complex geometry is designed based on an iterative bisection technique.Besides,an MPI/OpenMP hybrid parallel approach is developed to reduce the memory cost in each High-Performance Computing (HPC) node for large-scale applications.Three 3D cases,including the ONERA-M6 wing and a commercial transport airplane standard model with up to 2.5 billion hybrid elements,are used to test the present mesh deformation method.The robustness and high parallel efficiency are demonstrated by a wing deflection case with a maximum bending angle of 450 and more than 80% parallel efficiency with 1024 MPI processors.In addition,the availability for both continuous and discontinuous surface deformation is verified by interpolating the projecting displacement with opposite directions surface points to the spatial points.展开更多
Considering the advantages and limitations of traditional identification method,combined with the strategy of active detection,the principle of DC grid pilot protection based on active detection is proposed to improve...Considering the advantages and limitations of traditional identification method,combined with the strategy of active detection,the principle of DC grid pilot protection based on active detection is proposed to improve the sensitivity and reliability of hybrid MMC DC grid protection,and to ensure the accurate identification of fault areas in DC grid.By using the DC fault ride-through control strategy of the hybrid sub-module MMC,the fault current at the converter station DC terminal is limited.Based on the high controllability of hybrid MMC,sinusoidal fault detection signals with the same frequency are injected into the line at each converter station.Based on model recognition,the capacitance model condition is satisfied by the detected signals at both ends during external faults whereas not satisfied during internal faults.The Spearman correlation coefficients is then introduced,and the correlation discriminant of capacitance model is constructed to realize fault area discrimination of DC grid.The simulation results show that the active detection protection scheme proposed in this paper can accurately identify the fault area of DC grid,and is not affected by fault impedance and has low sampling rate requirement.展开更多
The number of electric vehicles are increasing in the society as they are considered as zero emission vehicles and also because conventional fuels are becoming expensive. Additional electrical power should be produced...The number of electric vehicles are increasing in the society as they are considered as zero emission vehicles and also because conventional fuels are becoming expensive. Additional electrical power should be produced to meet the energy requirement of this increase in electric vehicle population. To use the existing grid infrastructure without any failure, installing distributed generator at secondary distribution network is essential. In this work, sizing of wind-driven permanent magnet synchronous generator photovoltaic hybrid distributed generating system has been attempted to meet the energy demand of electric vehicles of a particular residential area. Different feasible combinations for wind generator capacity and photovoltaic capacity are obtained to satisfy the additional energy requirement. Results are analyzed based on energy, financial payback periods and daily power profile of the hybrid system. Based on this analysis, the sizes of wind generator and photovoltaic array have been chosen to meet the energy demand of electric vehicles of that particular residential locality.展开更多
A new approach to high-order accuracy for the numerical solution of conservation laws introduced by Huynh and extended to simplexes by Wang and Gao is renamed CPR(correction procedure or collocation penalty via recons...A new approach to high-order accuracy for the numerical solution of conservation laws introduced by Huynh and extended to simplexes by Wang and Gao is renamed CPR(correction procedure or collocation penalty via reconstruction).The CPR approach employs the differential form of the equation and accounts for the jumps in flux values at the cell boundaries by a correction procedure.In addition to being simple and economical,it unifies several existing methods including discontinuous Galerkin,staggered grid,spectral volume,and spectral difference.To discretize the diffusion terms,we use the BR2(Bassi and Rebay),interior penalty,compact DG(CDG),and I-continuous approaches.The first three of these approaches,originally derived using the integral formulation,were recast here in the CPR framework,whereas the I-continuous scheme,originally derived for a quadrilateral mesh,was extended to a triangular mesh.Fourier stability and accuracy analyses for these schemes on quadrilateral and triangular meshes are carried out.Finally,results for the Navier-Stokes equations are shown to compare the various schemes as well as to demonstrate the capability of the CPR approach.展开更多
A CFD-based Numerical Virtual Flight(NVF)simulator is presented,which integrates an unsteady flow solver on moving hybrid grids,a Rigid-Body Dynamics(RBD)solver and a module of the Flight Control System(FCS).A techni...A CFD-based Numerical Virtual Flight(NVF)simulator is presented,which integrates an unsteady flow solver on moving hybrid grids,a Rigid-Body Dynamics(RBD)solver and a module of the Flight Control System(FCS).A technique of dynamic hybrid grids is developed to control the active control surfaces with body morphing,with a technique of parallel unstructured dynamic overlapping grids generating proper moving grids over the deflecting control surfaces(e.g.the afterbody rudders of a missile).For the flow/kinematic coupled problems,the 6 Degree-Of-Freedom(DOF)equations are solved by an explicit or implicit method coupled with the URANS CFD solver.The module of the control law is explicitly coupled into the NVF simulator and then improved by the simulation of the pitching maneuver process of a maneuverable missile model.A nonlinear dynamic inversion method is then implemented to design the control law for the pitching process of the maneuverable missile model.Simulations and analysis of the pitching maneuver process are carried out by the NVF simulator to improve the flight control law.Higher control response performance is obtained by adjusting the gain factors and adding an integrator into the control loop.展开更多
文摘It has been evident that the theory and methods of dynamic derivatives are playing an increasingly important role in hybrid modeling and computations. Being constructed on various kinds of hybrid grids, that is, time scales, dynamic derivatives offer superior accuracy and flexibility in approximating mathematically important natural processes with hard-to-predict singularities, such as the epidemic growth with unpredictable jump sizes and option market changes with high uncertainties, as compared with conventional derivatives. In this article, we shall review the novel new concepts, explore delicate relations between the most frequently used second-order dynamic derivatives and conventional derivatives. We shall investigate necessary conditions for guaranteeing the consistency between the two derivatives. We will show that such a consistency may never exist in general. This implies that the dynamic derivatives provide entirely different new tools for sensitive modeling and approximations on hybrid grids. Rigorous error analysis will be given via asymptotic expansions for further modeling and computational applications. Numerical experiments will also be given.
文摘An unstructured Reynolds-averaged Navier-Stokes flow solver using the finite volume method is studied. The spatial discretisation is based on the Osher approximate Riemann solvers. A two-equation turbulence model (k-ω model) is also developed for hybrid grids to compute the turbulence flow. The turbulence flow past NACA0012 airfoil and the double ellipsolids are computed, and the numerical results show that the above methods are very efficient.
文摘Viscous CFD computations are performed to predict the aerodynamic coefficients and flowfield for a missile with grid fins by using hybrid grids. The computations arc made at Math numbvrs of 0.7 and 2.0. Full N-S equations arc discretized into finite volume form and solved by an algorithm in LU-SGS. The comparisons between computation and experimental data are made, and the detailed flow structure near grid fin is shown and examined.
基金supported by the National Defense Basic Research Program during the Twelfth Five-Year Plan Period
文摘The influence of chemical nonequilibrium on the thermal characteristics is explored by using the 2Dhybrid grid direct simulation Monte Carlo(DSMC)parallel method.An improved molecule search algorithm is proposed,which can preserve the high efficiency of area search algorithm.This method can overcome the defects of area search algorithm,and give all information about molecules hitting surface.The heat flux calculation method for a rarefied hypersonic flow is established.In addition,the testing methods of chemical reaction probability for five species of mixed gas with limited speed chemical reactions are also selected.To validate the effectiveness of the present method,hypersonic flow around a cylinder is firstly simulated,and subsequently numerical simulations of the heat flux and flow field characteristics around the blunt body at different heights are carried out in two different cases:the thermal nonequilibrium condition and the thermochemical nonequilibrium condition.Numerical results demonstrate the validity and reliability of the proposed methods.
基金Supported by the National Natural Science Foundation of China(11102179)
文摘An efficient compressible Euler equation solver for vortex-dominated flows is presented based on the adaptive hybrid Cartesian mesh and vortex identifying method.For most traditional grid-based Euler solvers,the excessive numerical dissipation is the great obstruction for vortex capturing or tracking problems.A vortex identifying method based on the curl of velocity is used to identify the vortex in flow field.Moreover,a dynamic adaptive mesh refinement(DAMR)process for hybrid Cartesian gird system is employed to track and preserve vortex.To validate the proposed method,a single compressible vortex convection flow is involved to test the accuracy and efficiency of DAMR process.Additionally,the vortex-dominated flow is investigated by the method.The obtained results are shown as a good agreement with the previous published data.
基金supported partly by the National Basic Research Program of China(″973″Program)(No.2014CB046200)
文摘A hybrid Cartesian structured grid method is proposed for solving moving boundary unsteady problems. The near body region is discretized by using the body-fitted structured grids, while the remaining computational domain is tessellated with the generated Cartesian grids. As the body moves, the structured grids move with the body and the outer boundaries of inside grids are used to generate new holes in the outside adaptive Cartesian grid to facilitate data communication. By using the alternating digital tree （ADT） algorithm, the computational time of hole-cutting and identification of donor cells can be reduced significantly. A compressible solver for unsteady flow problems is developed. A cell-centered, second-order accurate finite volume method is employed in spatial discreti- zation and an implicit dual-time stepping low-upper symmetric Gauss-Seidei （LU-SGS） approach is employed in temporal discretization. Geometry-based adaptation is used during unsteady simulation time steps when boundary moves and the flow solution is interpolated from the old Cartesian grids to the new one with inverse distance weigh- ting interpolation formula. Both laminar and turbulent unsteady cases are tested to demonstrate the accuracy and efficiency of the proposed method. Then, a 2-D store separation problem is simulated. The result shows that the hybrid Cartesian grid method can handle the unsteady flow problems involving large-scale moving boundaries.
文摘Strong DC coupling with weak AC and large-scale renewable energy integration are the two significant characteristics of ultra-high-voltage AC/DC(UHVAC/DC)hybrid power grids in China.Strong coupling between AC and DC grids and the different integration performance of renewable energy sources have profoundly changed the stability characteristics of the power system.The traditional stability control system is inadequate for the stability control of UHVAC/DC power grids.This paper analyzes the requirements for constructing an integrated defense system in a UHVAC/DC hybrid power grid(i.e.power system protection).The definition,connotation,and designing principles of power system protection are put forward.The relationship between the power system protection and the traditional three-defense lines is investigated.The design principles,general hardware structure and main functions of a power system protection are presented.Key problems and technologies are specified in the construction of the power system protection.
基金Sponsored by the National Natural Science Foundation of China(Grant No.51076066)
文摘A projectile with exotic wraparound wings( WAW) configuration is designed to improve the finstabilized projectile shooting quality. Two fin-stabilized projectiles with the same body with and without exotic WAW configuration are simulated numerically by applying the Roe scheme. The shear-stress transport turbulence models and the lower-upper symmetric Gauss-Seidel implicit method are used to solve 3D Reynoldsaveraged Navier—Stokes equations. The differences in aerodynamic coefficients and aerodynamic characteristics of the projectiles when the Mach number varies from 0. 35 to 0. 95 are obtained,and the cause of these differences is analyzed. The calculation results indicate that the lift-to-drag ratio of the projectile significantly increases,the rolling moment decreases,and the position of the pressure center of the projectile shows relatively small changes when the exotic WAW configuration is used. Therefore,this projectile can obviously reduce rolling effect,enlarge range and improve flying stability.
基金supported by the National Key Research and Development Program of China (No.2016YFB0200701)the National Natural Science Foundation of China (Nos. 11532016 and 91530325)
文摘An efficient MPI/OpenMP hybrid parallel Radial Basis Function (RBF) strategy for both continuous and discontinuous large-scale mesh deformation is proposed to reduce the computational cost and memory consumption.Unlike the conventional parallel methods in which all processors use the same surface displacement and implement the same operation,the present method employs different surface points sets and influence radius for each volume point movement,accompanied with efficient geometry searching strategy.The deformed surface points,also called Control Points (CPs),are stored in each processor.The displacement of spatial points is interpolated by using only 20-50 nearest control points,and the local influence radius is set to 5-20 times the maximum displacement of control points.To shorten the searching time for the nearest control point clouds,an Alternating Digital Tree (ADT) algorithm for 3D complex geometry is designed based on an iterative bisection technique.Besides,an MPI/OpenMP hybrid parallel approach is developed to reduce the memory cost in each High-Performance Computing (HPC) node for large-scale applications.Three 3D cases,including the ONERA-M6 wing and a commercial transport airplane standard model with up to 2.5 billion hybrid elements,are used to test the present mesh deformation method.The robustness and high parallel efficiency are demonstrated by a wing deflection case with a maximum bending angle of 450 and more than 80% parallel efficiency with 1024 MPI processors.In addition,the availability for both continuous and discontinuous surface deformation is verified by interpolating the projecting displacement with opposite directions surface points to the spatial points.
基金supported by The National Natural Science Foundation key project(U1766209).
文摘Considering the advantages and limitations of traditional identification method,combined with the strategy of active detection,the principle of DC grid pilot protection based on active detection is proposed to improve the sensitivity and reliability of hybrid MMC DC grid protection,and to ensure the accurate identification of fault areas in DC grid.By using the DC fault ride-through control strategy of the hybrid sub-module MMC,the fault current at the converter station DC terminal is limited.Based on the high controllability of hybrid MMC,sinusoidal fault detection signals with the same frequency are injected into the line at each converter station.Based on model recognition,the capacitance model condition is satisfied by the detected signals at both ends during external faults whereas not satisfied during internal faults.The Spearman correlation coefficients is then introduced,and the correlation discriminant of capacitance model is constructed to realize fault area discrimination of DC grid.The simulation results show that the active detection protection scheme proposed in this paper can accurately identify the fault area of DC grid,and is not affected by fault impedance and has low sampling rate requirement.
文摘The number of electric vehicles are increasing in the society as they are considered as zero emission vehicles and also because conventional fuels are becoming expensive. Additional electrical power should be produced to meet the energy requirement of this increase in electric vehicle population. To use the existing grid infrastructure without any failure, installing distributed generator at secondary distribution network is essential. In this work, sizing of wind-driven permanent magnet synchronous generator photovoltaic hybrid distributed generating system has been attempted to meet the energy demand of electric vehicles of a particular residential area. Different feasible combinations for wind generator capacity and photovoltaic capacity are obtained to satisfy the additional energy requirement. Results are analyzed based on energy, financial payback periods and daily power profile of the hybrid system. Based on this analysis, the sizes of wind generator and photovoltaic array have been chosen to meet the energy demand of electric vehicles of that particular residential locality.
文摘A new approach to high-order accuracy for the numerical solution of conservation laws introduced by Huynh and extended to simplexes by Wang and Gao is renamed CPR(correction procedure or collocation penalty via reconstruction).The CPR approach employs the differential form of the equation and accounts for the jumps in flux values at the cell boundaries by a correction procedure.In addition to being simple and economical,it unifies several existing methods including discontinuous Galerkin,staggered grid,spectral volume,and spectral difference.To discretize the diffusion terms,we use the BR2(Bassi and Rebay),interior penalty,compact DG(CDG),and I-continuous approaches.The first three of these approaches,originally derived using the integral formulation,were recast here in the CPR framework,whereas the I-continuous scheme,originally derived for a quadrilateral mesh,was extended to a triangular mesh.Fourier stability and accuracy analyses for these schemes on quadrilateral and triangular meshes are carried out.Finally,results for the Navier-Stokes equations are shown to compare the various schemes as well as to demonstrate the capability of the CPR approach.
基金supported partially by National Key Research and Development Program (No. 2016YFB0200701)National Natural Science Foundation of China (Nos. 11532016 and 11672324)
文摘A CFD-based Numerical Virtual Flight(NVF)simulator is presented,which integrates an unsteady flow solver on moving hybrid grids,a Rigid-Body Dynamics(RBD)solver and a module of the Flight Control System(FCS).A technique of dynamic hybrid grids is developed to control the active control surfaces with body morphing,with a technique of parallel unstructured dynamic overlapping grids generating proper moving grids over the deflecting control surfaces(e.g.the afterbody rudders of a missile).For the flow/kinematic coupled problems,the 6 Degree-Of-Freedom(DOF)equations are solved by an explicit or implicit method coupled with the URANS CFD solver.The module of the control law is explicitly coupled into the NVF simulator and then improved by the simulation of the pitching maneuver process of a maneuverable missile model.A nonlinear dynamic inversion method is then implemented to design the control law for the pitching process of the maneuverable missile model.Simulations and analysis of the pitching maneuver process are carried out by the NVF simulator to improve the flight control law.Higher control response performance is obtained by adjusting the gain factors and adding an integrator into the control loop.
