A parallel virtual machine (PVM) protocol based parallel computation of 3-D hypersonic flows with chemical non-equilibrium on hybrid meshes is presented. The numerical simulation for hypersonic flows with chemical n...A parallel virtual machine (PVM) protocol based parallel computation of 3-D hypersonic flows with chemical non-equilibrium on hybrid meshes is presented. The numerical simulation for hypersonic flows with chemical non-equilibrium reactions encounters the stiffness problem, thus taking huge CPU time. Based on the domain decomposition method, a high efficient automatic domain decomposer for three-dimensional hybrid meshes is developed, and then implemented to the numerical simulation of hypersonic flows. Control equations are multicomponent N-S equations, and spatially discretized scheme is used by a cell-centered finite volume algorithm with a five-stage Runge-Kutta time step. The chemical kinetic model is a seven species model with weak ionization. A point-implicit method is used to solve the chemical source term. Numerical results on PC-Cluster are verified on a bi-ellipse model compared with references.展开更多
Efficient and robust solution strategies are developed for discontinuous Galerkin (DG) discretization of the Navier-Stokes (NS) and Reynolds-averaged NS (RANS) equations on structured/unstructured hybrid meshes....Efficient and robust solution strategies are developed for discontinuous Galerkin (DG) discretization of the Navier-Stokes (NS) and Reynolds-averaged NS (RANS) equations on structured/unstructured hybrid meshes. A novel line-implicit scheme is devised and implemented to reduce the memory gain and improve the computational eificiency for highly anisotropic meshes. A simple and effective technique to use the mod- ified Baldwin-Lomax (BL) model on the unstructured meshes for the DC methods is proposed. The compact Hermite weighted essentially non-oscillatory (HWENO) limiters are also investigated for the hybrid meshes to treat solution discontinuities. A variety of compressible viscous flows are performed to examine the capability of the present high- order DG solver. Numerical results indicate that the designed line-implicit algorithms exhibit weak dependence on the cell aspect-ratio as well as the discretization order. The accuracy and robustness of the proposed approaches are demonstrated by capturing com- plex flow structures and giving reliable predictions of benchmark turbulent problems.展开更多
In a Wireless Mesh Network(WMN),the convenience of a routing strategy strongly depends on the mobility of the intermediate nodes that compose the paths.Taking this behaviour into account,this paper presents a routing ...In a Wireless Mesh Network(WMN),the convenience of a routing strategy strongly depends on the mobility of the intermediate nodes that compose the paths.Taking this behaviour into account,this paper presents a routing scheme that works differently accordingly to the node mobility.In this sense,a proactive routing scheme is restricted to the backbone to promote the use of stable routes.Conversely,the reactive protocol is used for searching routes to or from a mobile destination.Both approaches are simultaneously implemented in the mesh nodes so that the routing protocols share routing information that optimises the network performance.Aimed at guaranteeing the IP compatibility,the combination of the two protocols in the core routers is carried out in the Medium Access Control(MAC)layer.In contrast to the operation in the IP layer where two routing protocols cannot work concurrently,the transfer of the routing tasks to the MAC layer enables the use of multiple independent forwarding tables.Simulation results show the advantage of the proposal in terms of packet losses and data delay.展开更多
The computational fluid dynamics (CFD) method is used to numerically simulate a propeller wake flow field in open water. A sub-domain hybrid mesh method was adopted in this paper. The computation domain was separate...The computational fluid dynamics (CFD) method is used to numerically simulate a propeller wake flow field in open water. A sub-domain hybrid mesh method was adopted in this paper. The computation domain was separated into two sub-domains, in which tetrahedral elements were used in the inner domain to match the complicated geometry of the propeller, while hexahedral elements were used in the outer domain. The mesh was locally refined on the propeller surface and near the wake flow field, and a size function was used to control the growth rate of the grid. Sections at different axial location were used to study the spatial evolution of the propeller wake in the region ranging from the disc to one propeller diameter (D) downstream. The numerical results show that the axial velocity fluctuates along the wake flow; radial velocity, which is closely related to vortices, attenuates strongly. The trailing vortices interact with the tip vortex at the blades' trailing edge and then separate. The strength of the vortex shrinks rapidly, and the radius decreases 20% at one diameter downstream.展开更多
As the maneuverability of a ship navigating close to a bank is influenced by the sidewall, the assessment of ship maneuvering stability is important. The hydrodynamic derivatives measured by the planar motion mechani...As the maneuverability of a ship navigating close to a bank is influenced by the sidewall, the assessment of ship maneuvering stability is important. The hydrodynamic derivatives measured by the planar motion mechanism (PMM) test provide a way to predict the change of ship maneuverability. This paper presents a numerical simulation of PMM model tests with variant distances to a vertical bank by using unsteady RANS equations. A hybrid dynamic mesh technique is developed to realize the mesh configuration and remeshing of dynamic PMM tests when the ship is close to the bank. The proposed method is validated by comparing numerical results with results of PMM tests in a circulating water channel. The first-order hydrodynamic derivatives of the ship are analyzed from the time history of lateral force and yaw moment according to the multiple-run simulating procedure and the variations of hydrodynamic derivatives with the ship-sidewall distance are given. The straight line stability and directional stability are also discussed and stable or unstable zone of proportional-derivative (PD) controller parameters for directional stability is shown, which can be a reference for course keeping operation when sailing near a bank.展开更多
A hybrid grid generation technique and a multigrid/parallel algorithm are presented in this paper for turbulence flow simulations over three-dimensional (3D) complex geometries. The hybrid grid generation technique ...A hybrid grid generation technique and a multigrid/parallel algorithm are presented in this paper for turbulence flow simulations over three-dimensional (3D) complex geometries. The hybrid grid generation technique is based on an agglomeration method of anisotropic tetrahedrons. Firstly, the complex computational domain is covered by pure tetrahedral grids, in which anisotropic tetrahedrons are adopted to discrete the boundary layer and isotropic tetrahedrons in the outer field. Then, the anisotropic tetrahedrons in the boundary layer are agglomerated to generate prismatic grids. The agglomeration method can improve the grid quality in boundary layer and reduce the grid quantity to enhance the numerical accuracy and efficiency. In order to accelerate the convergence history, a multigrid/parallel algorithm is developed also based on anisotropic agglomeration approach. The numerical results demonstrate the excellent accelerating capability of this multigrid method.展开更多
The lower-upper symmetric Gauss-Seidel (LU-SGS) implicit relaxation has been widely used because it has the merits of less dependency on grid topology, low numerical complexity and modest memory requirements. In ori...The lower-upper symmetric Gauss-Seidel (LU-SGS) implicit relaxation has been widely used because it has the merits of less dependency on grid topology, low numerical complexity and modest memory requirements. In original LU-SGS scheme, the implicit system matrix is constructed based on the splitting of convective flux Jacobian according to its spectral radius. Although this treatment has the merit of reducing computational complexity and helps to ensure the diagonally dominant property of the implicit system matrix, it can also cause serious distortions on the implicit system matrix because too many approximations are introduced by this splitting method if the contravariant velocity is small or close to sonic speed. To overcome this shortcoming, an improved LU-SGS scheme with a hybrid construction method for the implicit system matrix is developed in this paper. The hybrid way is that: on the cell faces having small contravariant velocity or transonic contravariant velocity, the accurate derivative of the convective flux term is used to construct more accurate implicit system matrix, while the original Jacobian splitting method is adopted on the other cell faces to reduce computational complexity and ensure the diagonally dominant property of the implicit system matrix. To investigate the convergence performance of the improved LU-SGS scheme, 2D and 3D turbulent flows around the NACA0012 airfoil, RAE2822 airfoil and LANN wing are simulated on hybrid unstructured meshes. The nu- merical results show that the improved LU-SGS scheme is significantly more efficient than the original LU-SGS scheme.展开更多
In this paper,we present a hybrid grid generation approach for viscous flow simulations by marching a surface triangulation on viscous walls along certain directions.Focuses are on the computing strategies used to det...In this paper,we present a hybrid grid generation approach for viscous flow simulations by marching a surface triangulation on viscous walls along certain directions.Focuses are on the computing strategies used to determine the marching directions and distances since these strategies determine the quality of the resulting elements and the reliability of the meshing procedure to a large extent.With respect to marching directions,three strategies featured with different levels of efficiencies and robustness performance are combined to compute the initial normals at front nodes to balance the trade-off between efficiency and robustness.A novel weighted strategy is used in the normal smoothing scheme,which evidently reduces the possibility of early stop of front generation at complex corners.With respect to marching distances,the distance settings at concave and/or convex corners are locally adjusted to smooth the front shape at first;a further adjustment is then conducted for front nodes in the neighbourhood of gaps between opposite viscous boundaries.These efforts,plus other special treatments such as multi-normal generation and fast detection of local/global intersection,as a whole enable the setup of a hybrid mesher that could generate qualitied viscous grids for geometries with industry-level complexities.展开更多
The main goal of this paper is to study the characteristics of regression rate of solid grain during thrust regulation process. For this purpose, an unsteady numerical model of regression rate is established. Gas–sol...The main goal of this paper is to study the characteristics of regression rate of solid grain during thrust regulation process. For this purpose, an unsteady numerical model of regression rate is established. Gas–solid coupling is considered between the solid grain surface and combustion gas.Dynamic mesh is used to simulate the regression process of the solid fuel surface. Based on this model, numerical simulations on a H2O2/HTPB(hydroxyl-terminated polybutadiene) hybrid motor have been performed in the flow control process. The simulation results show that under the step change of the oxidizer mass flow rate condition, the regression rate cannot reach a stable value instantly because the flow field requires a short time period to adjust. The regression rate increases with the linear gain of oxidizer mass flow rate, and has a higher slope than the relative inlet function of oxidizer flow rate. A shorter regulation time can cause a higher regression rate during regulation process. The results also show that transient calculation can better simulate the instantaneous regression rate in the operation process.展开更多
This paper proposes a novel algorithm, which can be used to model and analyzemesh tree hybrid power/ground distribution networks with multiple voltage supply in time domain.Not only this algorithm enhances common meth...This paper proposes a novel algorithm, which can be used to model and analyzemesh tree hybrid power/ground distribution networks with multiple voltage supply in time domain.Not only this algorithm enhances common method''s ability on analysis of power/ground network withirregular topology, but also very high accuracy it keeps. The accuracy and stability of thisalgorithm is proved using strict math method in this paper. Also, the usage of both preconditiontechnique based on Incomplete Choleskey Decomposition and fast variable elimination technique hasimproved the algorithm''s efficiency a lot. Experimental results show that it can finish the analysisof power/ground network with enormous, size within very short time. Also, this algorithm can beapplied to analyze the clock network, bus network, and signal network without buffer under highworking frequency because of the independence of the topology.展开更多
Network coding mechanisms, such as COPE, can improve network throughput effectively in Wireless Mesh Networks(WMN). While the Hybrid Wireless Mesh Protocol(HWMP) is suitable for WMN, its extension with COPE does n...Network coding mechanisms, such as COPE, can improve network throughput effectively in Wireless Mesh Networks(WMN). While the Hybrid Wireless Mesh Protocol(HWMP) is suitable for WMN, its extension with COPE does not provide any added benefits; specifically, HWMP cannot establish paths with more coding opportunities. As a result, the advantages of network coding cannot be exploited sufficiently. This paper proposes improvements upon HWMP with a new, network Coding-Aware routing protocol(CAHWMP) for WMN. In the CAHWMP protocol, we propose a coding criterion based on data streams to devise an algorithm for actively detecting coding opportunities during path discovery. CAHWMP subsequently establishes paths using the codingaware routing metric, which can balance channel resource consumption and the gain due to sharing resources introduced by network coding. Simulation results show that CAHWMP can establish paths with more coding opportunities; as a result, it improves network performance such as network throughput.展开更多
文摘A parallel virtual machine (PVM) protocol based parallel computation of 3-D hypersonic flows with chemical non-equilibrium on hybrid meshes is presented. The numerical simulation for hypersonic flows with chemical non-equilibrium reactions encounters the stiffness problem, thus taking huge CPU time. Based on the domain decomposition method, a high efficient automatic domain decomposer for three-dimensional hybrid meshes is developed, and then implemented to the numerical simulation of hypersonic flows. Control equations are multicomponent N-S equations, and spatially discretized scheme is used by a cell-centered finite volume algorithm with a five-stage Runge-Kutta time step. The chemical kinetic model is a seven species model with weak ionization. A point-implicit method is used to solve the chemical source term. Numerical results on PC-Cluster are verified on a bi-ellipse model compared with references.
基金Project supported by the National Basic Research Program of China(No.2009CB724104)
文摘Efficient and robust solution strategies are developed for discontinuous Galerkin (DG) discretization of the Navier-Stokes (NS) and Reynolds-averaged NS (RANS) equations on structured/unstructured hybrid meshes. A novel line-implicit scheme is devised and implemented to reduce the memory gain and improve the computational eificiency for highly anisotropic meshes. A simple and effective technique to use the mod- ified Baldwin-Lomax (BL) model on the unstructured meshes for the DC methods is proposed. The compact Hermite weighted essentially non-oscillatory (HWENO) limiters are also investigated for the hybrid meshes to treat solution discontinuities. A variety of compressible viscous flows are performed to examine the capability of the present high- order DG solver. Numerical results indicate that the designed line-implicit algorithms exhibit weak dependence on the cell aspect-ratio as well as the discretization order. The accuracy and robustness of the proposed approaches are demonstrated by capturing com- plex flow structures and giving reliable predictions of benchmark turbulent problems.
文摘In a Wireless Mesh Network(WMN),the convenience of a routing strategy strongly depends on the mobility of the intermediate nodes that compose the paths.Taking this behaviour into account,this paper presents a routing scheme that works differently accordingly to the node mobility.In this sense,a proactive routing scheme is restricted to the backbone to promote the use of stable routes.Conversely,the reactive protocol is used for searching routes to or from a mobile destination.Both approaches are simultaneously implemented in the mesh nodes so that the routing protocols share routing information that optimises the network performance.Aimed at guaranteeing the IP compatibility,the combination of the two protocols in the core routers is carried out in the Medium Access Control(MAC)layer.In contrast to the operation in the IP layer where two routing protocols cannot work concurrently,the transfer of the routing tasks to the MAC layer enables the use of multiple independent forwarding tables.Simulation results show the advantage of the proposal in terms of packet losses and data delay.
基金Supported by Fundamental Research Funds for the Central Universities(Grant No.HEUCFT1001)Ph.D Programs Foundation of Ministry of Education of China(Grant No.10702016)
文摘The computational fluid dynamics (CFD) method is used to numerically simulate a propeller wake flow field in open water. A sub-domain hybrid mesh method was adopted in this paper. The computation domain was separated into two sub-domains, in which tetrahedral elements were used in the inner domain to match the complicated geometry of the propeller, while hexahedral elements were used in the outer domain. The mesh was locally refined on the propeller surface and near the wake flow field, and a size function was used to control the growth rate of the grid. Sections at different axial location were used to study the spatial evolution of the propeller wake in the region ranging from the disc to one propeller diameter (D) downstream. The numerical results show that the axial velocity fluctuates along the wake flow; radial velocity, which is closely related to vortices, attenuates strongly. The trailing vortices interact with the tip vortex at the blades' trailing edge and then separate. The strength of the vortex shrinks rapidly, and the radius decreases 20% at one diameter downstream.
基金supported by the National Key Basic Research Program of China(Grant No.2014CB046804)
文摘As the maneuverability of a ship navigating close to a bank is influenced by the sidewall, the assessment of ship maneuvering stability is important. The hydrodynamic derivatives measured by the planar motion mechanism (PMM) test provide a way to predict the change of ship maneuverability. This paper presents a numerical simulation of PMM model tests with variant distances to a vertical bank by using unsteady RANS equations. A hybrid dynamic mesh technique is developed to realize the mesh configuration and remeshing of dynamic PMM tests when the ship is close to the bank. The proposed method is validated by comparing numerical results with results of PMM tests in a circulating water channel. The first-order hydrodynamic derivatives of the ship are analyzed from the time history of lateral force and yaw moment according to the multiple-run simulating procedure and the variations of hydrodynamic derivatives with the ship-sidewall distance are given. The straight line stability and directional stability are also discussed and stable or unstable zone of proportional-derivative (PD) controller parameters for directional stability is shown, which can be a reference for course keeping operation when sailing near a bank.
基金supported partially by National Basic Research Program of China (Grant No. 2009CB723800)National Natural Science Foundation of China (Grant Nos: 91016001 and 10872023)
文摘A hybrid grid generation technique and a multigrid/parallel algorithm are presented in this paper for turbulence flow simulations over three-dimensional (3D) complex geometries. The hybrid grid generation technique is based on an agglomeration method of anisotropic tetrahedrons. Firstly, the complex computational domain is covered by pure tetrahedral grids, in which anisotropic tetrahedrons are adopted to discrete the boundary layer and isotropic tetrahedrons in the outer field. Then, the anisotropic tetrahedrons in the boundary layer are agglomerated to generate prismatic grids. The agglomeration method can improve the grid quality in boundary layer and reduce the grid quantity to enhance the numerical accuracy and efficiency. In order to accelerate the convergence history, a multigrid/parallel algorithm is developed also based on anisotropic agglomeration approach. The numerical results demonstrate the excellent accelerating capability of this multigrid method.
基金Foundation item: National Natural Science Foundation of China (10802067)
文摘The lower-upper symmetric Gauss-Seidel (LU-SGS) implicit relaxation has been widely used because it has the merits of less dependency on grid topology, low numerical complexity and modest memory requirements. In original LU-SGS scheme, the implicit system matrix is constructed based on the splitting of convective flux Jacobian according to its spectral radius. Although this treatment has the merit of reducing computational complexity and helps to ensure the diagonally dominant property of the implicit system matrix, it can also cause serious distortions on the implicit system matrix because too many approximations are introduced by this splitting method if the contravariant velocity is small or close to sonic speed. To overcome this shortcoming, an improved LU-SGS scheme with a hybrid construction method for the implicit system matrix is developed in this paper. The hybrid way is that: on the cell faces having small contravariant velocity or transonic contravariant velocity, the accurate derivative of the convective flux term is used to construct more accurate implicit system matrix, while the original Jacobian splitting method is adopted on the other cell faces to reduce computational complexity and ensure the diagonally dominant property of the implicit system matrix. To investigate the convergence performance of the improved LU-SGS scheme, 2D and 3D turbulent flows around the NACA0012 airfoil, RAE2822 airfoil and LANN wing are simulated on hybrid unstructured meshes. The nu- merical results show that the improved LU-SGS scheme is significantly more efficient than the original LU-SGS scheme.
基金National Numerical Wind Tunnel Project of China(Grant No.NNW2018-ZT5B10).
文摘In this paper,we present a hybrid grid generation approach for viscous flow simulations by marching a surface triangulation on viscous walls along certain directions.Focuses are on the computing strategies used to determine the marching directions and distances since these strategies determine the quality of the resulting elements and the reliability of the meshing procedure to a large extent.With respect to marching directions,three strategies featured with different levels of efficiencies and robustness performance are combined to compute the initial normals at front nodes to balance the trade-off between efficiency and robustness.A novel weighted strategy is used in the normal smoothing scheme,which evidently reduces the possibility of early stop of front generation at complex corners.With respect to marching distances,the distance settings at concave and/or convex corners are locally adjusted to smooth the front shape at first;a further adjustment is then conducted for front nodes in the neighbourhood of gaps between opposite viscous boundaries.These efforts,plus other special treatments such as multi-normal generation and fast detection of local/global intersection,as a whole enable the setup of a hybrid mesher that could generate qualitied viscous grids for geometries with industry-level complexities.
基金co-supported by the Innovation Foundation of Beihang University for Ph.D. Graduatesthe National Natural Science Foundation of China (No. 51206007)
文摘The main goal of this paper is to study the characteristics of regression rate of solid grain during thrust regulation process. For this purpose, an unsteady numerical model of regression rate is established. Gas–solid coupling is considered between the solid grain surface and combustion gas.Dynamic mesh is used to simulate the regression process of the solid fuel surface. Based on this model, numerical simulations on a H2O2/HTPB(hydroxyl-terminated polybutadiene) hybrid motor have been performed in the flow control process. The simulation results show that under the step change of the oxidizer mass flow rate condition, the regression rate cannot reach a stable value instantly because the flow field requires a short time period to adjust. The regression rate increases with the linear gain of oxidizer mass flow rate, and has a higher slope than the relative inlet function of oxidizer flow rate. A shorter regulation time can cause a higher regression rate during regulation process. The results also show that transient calculation can better simulate the instantaneous regression rate in the operation process.
文摘This paper proposes a novel algorithm, which can be used to model and analyzemesh tree hybrid power/ground distribution networks with multiple voltage supply in time domain.Not only this algorithm enhances common method''s ability on analysis of power/ground network withirregular topology, but also very high accuracy it keeps. The accuracy and stability of thisalgorithm is proved using strict math method in this paper. Also, the usage of both preconditiontechnique based on Incomplete Choleskey Decomposition and fast variable elimination technique hasimproved the algorithm''s efficiency a lot. Experimental results show that it can finish the analysisof power/ground network with enormous, size within very short time. Also, this algorithm can beapplied to analyze the clock network, bus network, and signal network without buffer under highworking frequency because of the independence of the topology.
基金supported by the National Natural Science Foundation of China (No. 61301110)
文摘Network coding mechanisms, such as COPE, can improve network throughput effectively in Wireless Mesh Networks(WMN). While the Hybrid Wireless Mesh Protocol(HWMP) is suitable for WMN, its extension with COPE does not provide any added benefits; specifically, HWMP cannot establish paths with more coding opportunities. As a result, the advantages of network coding cannot be exploited sufficiently. This paper proposes improvements upon HWMP with a new, network Coding-Aware routing protocol(CAHWMP) for WMN. In the CAHWMP protocol, we propose a coding criterion based on data streams to devise an algorithm for actively detecting coding opportunities during path discovery. CAHWMP subsequently establishes paths using the codingaware routing metric, which can balance channel resource consumption and the gain due to sharing resources introduced by network coding. Simulation results show that CAHWMP can establish paths with more coding opportunities; as a result, it improves network performance such as network throughput.