Fine grids with small spacing in boundary-fitted coordinates are normally used to treat the computation of fluid dynamics for estuarine areas and tidal flats. However, the adoption of Cartesian components of velocity ...Fine grids with small spacing in boundary-fitted coordinates are normally used to treat the computation of fluid dynamics for estuarine areas and tidal flats. However, the adoption of Cartesian components of velocity vectors in this kind of non-orthogonal coordinates will definitely result in a difficulty in solving implicitly the transformed momentum equations, and also complicate the wet-dry point judgement used for flood areas. To solve this problem, equations in terms of generalized contravariant velocity vectors in curvilinear coordinates are derived in the present study, by virtue of which, an Alternative-Direction-Implicit numerical scheme in non-orthogonal grids would then be easily obtained, and wet-dry point judgement would as well be largely simplified. A comparison is made between the explicit scheme and implicit scheme, showing that the present model is accurate and numerically stable for computations of fluid dynamics for estuarine areas and tidal flats.展开更多
In this study, the momentum equations describing an atmospheric flow over a NW Pacific region of Mexico are solved numerically. In order to capture the complex flow-topography interactions with detail, a combination o...In this study, the momentum equations describing an atmospheric flow over a NW Pacific region of Mexico are solved numerically. In order to capture the complex flow-topography interactions with detail, a combination of a numerical wind model in full 3D curvilinear coordinates, along with a high resolution boundary-fitted grid is used. Boundary conditions were obtained from ten years (2002-2012) of measured offshore wind data. Prevailing winds from April to September during that period of observations were selected for the simulations. For the cases analyzed, it was found that at the points of the study region (PSS, PSM, PM), wind speed increased about 10% to 20% of its offshore values, while inland they decreased about 86% to 96%. This spatial behavior agreed very well with the observed local winds. A coastal jet (CJ), 35 km long with speeds about 1.5 - 2 m/s, emanating from PSS was found for NNW winds. Modeled winds were also used to compute wind stresses, wind stress curl, and CUI fields. Wind stress values agreed very well to those reported in the literature. High values of wind stress curl, and CUI were found at the lee of the points (PSS, PSM, PM). Indirect estimations of sea surface currents were about 15 - 20 cm/s offshore and 5 - 10 cm/s at the coast.展开更多
In order to overcome the zigzag grids generated by conventional finite difference method on complicated casting boundaries in the simulation of casting process, the generation program for 2-D boundary-fitted coordinat...In order to overcome the zigzag grids generated by conventional finite difference method on complicated casting boundaries in the simulation of casting process, the generation program for 2-D boundary-fitted coordinate grid has been developed by solving a set of partial differential equations (PDE) numerically. The STL format files were treated as input data for 2-D physical regions. The equipartition method for boundary points was used to improve the self-adaptability of grid according to the characteristic of the STL format files. The program was demonstrated through some examples. The comparison between the conventional finite difference method and the proposed method shows that this program is effective and flexible for generation of boundary-fitted grid in any arbitrary 2-D complex domain, and the grid is in accordance with the variety of boundary curvature finely. The program also provides two types of boundary-fitted grids for double-connected region, O-type and C-type. The limitation of the step-like boundary with the rectangle grid could be avoided effectively. Therefore, the computational accuracy and efficiency would be improved and the computational time would be saved significantly by the application of boundary-fitted grids.展开更多
A type of mesh called a body-fi tted Cartesian mesh,very different from the traditional structured body-fi tted mesh,is established.At f irst,the right parallelepiped mesh is generated,then,a feature analysis is done ...A type of mesh called a body-fi tted Cartesian mesh,very different from the traditional structured body-fi tted mesh,is established.At f irst,the right parallelepiped mesh is generated,then,a feature analysis is done on the cross sections.These cross sections are the intersections of the casting shape with the right parallelepiped grids(under the Cartesian coordinate system).On the basis of the feature analysis,two sorts of body-f itted boundary grids,shape-keeping grids and shape-distortion grids,are def ined.Shape-distortion grids can be removed or weaken by increasing the number of grids or moving the coordinates of the mesh generation region,so actually the body-fitted Cartesian mesh generation is to get shape-keeping grids.A shape-keeping grid mainly consists of two sorts of surfaces(I type face and II type face),and each of them is joined by two types of points(I type point and II type point).If only these two types of points were given,the shape-keeping mesh would be constructed.In this paper,the cases of the above two boundary grids being generated were discussed.An algorithm was put forward to get the shape-keeping grids.Several body-fi tted Cartesian meshes generated on castings show the validity of the algorithm.The mesh generation examples show that the body-fi tted Cartesian mesh is more excellent than the right parallelepiped mesh in aspects of decreasing grids number and being closer to the shape of the casting solid.展开更多
Terrain characteristics can be accurately represented in spectrum space. Terrain spectra can quantitatively reflect the effect of topographic dynamic forcing on the atmosphere. In wavelength space, topographic spectra...Terrain characteristics can be accurately represented in spectrum space. Terrain spectra can quantitatively reflect the effect of topographic dynamic forcing on the atmosphere. In wavelength space, topographic spectral energy decreases with decreasing wavelength, in spite of several departures. This relationship is approximated by an exponential function. A power law relationship between the terrain height spectra and wavelength is fitted by the least-squares method, and the fitting slope is associated with grid-size selection for mesoscale models. The monotonicity of grid size is investigated, and it is strictly proved that grid size increases with increasing fitting exponent, indicating that the universal grid size is determined by the minimum fitting exponent. An example of landslide-prone areas in western Sichuan is given, and the universal grid spacing of 4.1 km is shown to be a requirement to resolve 90% of terrain height variance for mesoscale models, without resorting to the parameterization of subgrid-scale terrain variance. Comparison among results of different simulations shows that the simulations estimate the observed precipitation well when using a resolution of 4.1 km or finer. Although the main flow patterns are similar, finer grids produce more complex patterns that show divergence zones, convergence zones and vortices. Horizontal grid size significantly affects the vertical structure of the convective boundary layer. Stronger vertical wind components are simulated for finer grid resolutions. In particular, noticeable sinking airflows over mountains are captured for those model configurations.展开更多
A new method for constructing a fitting surface on a triangular grid is presented. Assuming images are obtained by sampling from the original scene. Conventional polynomial interpolation methods generally construct th...A new method for constructing a fitting surface on a triangular grid is presented. Assuming images are obtained by sampling from the original scene. Conventional polynomial interpolation methods generally construct the fitting surface on a square grid. Different from existing methods, the new method constructs the fitting surface on a triangular grid which can divide the original surface more detailed and improve approximation accuracy. As the quality of the image edges plays a key role in visual effects of image, the new method uses image edges as constraints to get a triangle grid. The new method constructs a cubic polynomial patch locally using image data to approximate the original surface. Experimental comparison results of the new method with other methods show that the new method can produce high-quality images and remove the zigzagging artifact.展开更多
基金National Natural Science Foundation of China and National Excellent Youth Foundation of China.(Grant No.49606069)
文摘Fine grids with small spacing in boundary-fitted coordinates are normally used to treat the computation of fluid dynamics for estuarine areas and tidal flats. However, the adoption of Cartesian components of velocity vectors in this kind of non-orthogonal coordinates will definitely result in a difficulty in solving implicitly the transformed momentum equations, and also complicate the wet-dry point judgement used for flood areas. To solve this problem, equations in terms of generalized contravariant velocity vectors in curvilinear coordinates are derived in the present study, by virtue of which, an Alternative-Direction-Implicit numerical scheme in non-orthogonal grids would then be easily obtained, and wet-dry point judgement would as well be largely simplified. A comparison is made between the explicit scheme and implicit scheme, showing that the present model is accurate and numerically stable for computations of fluid dynamics for estuarine areas and tidal flats.
文摘In this study, the momentum equations describing an atmospheric flow over a NW Pacific region of Mexico are solved numerically. In order to capture the complex flow-topography interactions with detail, a combination of a numerical wind model in full 3D curvilinear coordinates, along with a high resolution boundary-fitted grid is used. Boundary conditions were obtained from ten years (2002-2012) of measured offshore wind data. Prevailing winds from April to September during that period of observations were selected for the simulations. For the cases analyzed, it was found that at the points of the study region (PSS, PSM, PM), wind speed increased about 10% to 20% of its offshore values, while inland they decreased about 86% to 96%. This spatial behavior agreed very well with the observed local winds. A coastal jet (CJ), 35 km long with speeds about 1.5 - 2 m/s, emanating from PSS was found for NNW winds. Modeled winds were also used to compute wind stresses, wind stress curl, and CUI fields. Wind stress values agreed very well to those reported in the literature. High values of wind stress curl, and CUI were found at the lee of the points (PSS, PSM, PM). Indirect estimations of sea surface currents were about 15 - 20 cm/s offshore and 5 - 10 cm/s at the coast.
基金supported by State Key Laboratory of Materials Processing and Die&Mould Technology,Huazhong University of Science and Technology(09-04)National Natural Science Foundation of China(No.50775050)
文摘In order to overcome the zigzag grids generated by conventional finite difference method on complicated casting boundaries in the simulation of casting process, the generation program for 2-D boundary-fitted coordinate grid has been developed by solving a set of partial differential equations (PDE) numerically. The STL format files were treated as input data for 2-D physical regions. The equipartition method for boundary points was used to improve the self-adaptability of grid according to the characteristic of the STL format files. The program was demonstrated through some examples. The comparison between the conventional finite difference method and the proposed method shows that this program is effective and flexible for generation of boundary-fitted grid in any arbitrary 2-D complex domain, and the grid is in accordance with the variety of boundary curvature finely. The program also provides two types of boundary-fitted grids for double-connected region, O-type and C-type. The limitation of the step-like boundary with the rectangle grid could be avoided effectively. Therefore, the computational accuracy and efficiency would be improved and the computational time would be saved significantly by the application of boundary-fitted grids.
基金financially supported by the fund of the State Key Laboratory of Solidif ication Processing in NWPU(No:SKLSP201006)the fund of the National Basic Research Program of China(No:2011CB610402)
文摘A type of mesh called a body-fi tted Cartesian mesh,very different from the traditional structured body-fi tted mesh,is established.At f irst,the right parallelepiped mesh is generated,then,a feature analysis is done on the cross sections.These cross sections are the intersections of the casting shape with the right parallelepiped grids(under the Cartesian coordinate system).On the basis of the feature analysis,two sorts of body-f itted boundary grids,shape-keeping grids and shape-distortion grids,are def ined.Shape-distortion grids can be removed or weaken by increasing the number of grids or moving the coordinates of the mesh generation region,so actually the body-fitted Cartesian mesh generation is to get shape-keeping grids.A shape-keeping grid mainly consists of two sorts of surfaces(I type face and II type face),and each of them is joined by two types of points(I type point and II type point).If only these two types of points were given,the shape-keeping mesh would be constructed.In this paper,the cases of the above two boundary grids being generated were discussed.An algorithm was put forward to get the shape-keeping grids.Several body-fi tted Cartesian meshes generated on castings show the validity of the algorithm.The mesh generation examples show that the body-fi tted Cartesian mesh is more excellent than the right parallelepiped mesh in aspects of decreasing grids number and being closer to the shape of the casting solid.
基金supported by the Key Research Program of the Chinese Academy of Sciences (Grant No. KZZD-EW-05-01)the special grant (Grant No. 41375052) from the National Natural Science Foundation of Chinafunded by an open project of the State Key Laboratory of Severe Weather (Grant No. 2013LASW-A06)
文摘Terrain characteristics can be accurately represented in spectrum space. Terrain spectra can quantitatively reflect the effect of topographic dynamic forcing on the atmosphere. In wavelength space, topographic spectral energy decreases with decreasing wavelength, in spite of several departures. This relationship is approximated by an exponential function. A power law relationship between the terrain height spectra and wavelength is fitted by the least-squares method, and the fitting slope is associated with grid-size selection for mesoscale models. The monotonicity of grid size is investigated, and it is strictly proved that grid size increases with increasing fitting exponent, indicating that the universal grid size is determined by the minimum fitting exponent. An example of landslide-prone areas in western Sichuan is given, and the universal grid spacing of 4.1 km is shown to be a requirement to resolve 90% of terrain height variance for mesoscale models, without resorting to the parameterization of subgrid-scale terrain variance. Comparison among results of different simulations shows that the simulations estimate the observed precipitation well when using a resolution of 4.1 km or finer. Although the main flow patterns are similar, finer grids produce more complex patterns that show divergence zones, convergence zones and vortices. Horizontal grid size significantly affects the vertical structure of the convective boundary layer. Stronger vertical wind components are simulated for finer grid resolutions. In particular, noticeable sinking airflows over mountains are captured for those model configurations.
基金Supported by National Natural Science Foundation of China(61572292,61373078,61272430)NSFC Joint Fund with Guangdong under Key Project(U1201258)
文摘A new method for constructing a fitting surface on a triangular grid is presented. Assuming images are obtained by sampling from the original scene. Conventional polynomial interpolation methods generally construct the fitting surface on a square grid. Different from existing methods, the new method constructs the fitting surface on a triangular grid which can divide the original surface more detailed and improve approximation accuracy. As the quality of the image edges plays a key role in visual effects of image, the new method uses image edges as constraints to get a triangle grid. The new method constructs a cubic polynomial patch locally using image data to approximate the original surface. Experimental comparison results of the new method with other methods show that the new method can produce high-quality images and remove the zigzagging artifact.