共查询到20条相似文献，搜索用时 31 毫秒
1.
A modified three-parameter model of turbulence for a thermally stratified atmospheric boundary layer (ABL) is presented. The model is based on tensor-invariant parametrizations for the pressure-strain and pressure-temperature correlations that are more complete than the parametrizations used in the Mellor-Yamada model of level 3.0. The turbulent momentum and heat fluxes are calculated with explicit algebraic models obtained with the aid of symbol algebra from the transport equations for momentum and heat fluxes in the approximation of weakly equilibrium turbulence. The turbulent transport of heat and momentum fluxes is assumed to be negligibly small in this approximation. The three-parameter E ? ε ? 2> model of thermally stratified turbulence is employed to obtain closed-form algebraic expressions for the fluxes. A computational test of a 24-h ABL evolution is implemented for an idealized two-dimensional region. Comparison of the computed results with the available observational data and other numerical models shows that the proposed model is able to reproduce both the most important structural features of the turbulence in an urban canopy layer near the urbanized ABL surface and the effect of urban roughness on a global structure of the fields of wind and temperature over a city. The results of the computational test for the new model indicate that the motion of air in the urban canopy layer is strongly influenced by mechanical factors (buildings) and thermal stratification. 相似文献
2.
A. S. Sarkisyan 《Oceanology》2016,56(5):615-620
In this paper, an integration method is proposed for a system of ocean dynamics equations (3DPEM) for creating a series of atlases with the joint title SOCRAT so that circulation in lower ocean layers can be adjusted to circulation in upper layers, avoiding smoothing of the intensity in upper layers caused by long integration or poor resolution. Specific simulations are presented, which were carried out by the author with a group of colleagues. These results conserve the realistic features of circulation thanks to the short integration period and high resolution. It is shown that due to long integration and/or simulation with rough resolution, simulations lose the main feature of ocean dynamics: intense alongshore currents; as a result, a nonexistant stable physical state uncharacteristic of an ocean is imposed. 相似文献
3.
《Coastal Engineering》2004,51(10):1021-1049
A numerical process-based model to forecast beach profile morphodynamics has been developed. In the present paper, an analysis of various modelling approaches and key parametrizations involved in the estimation of the wave-driven current and the suspended sediment concentration is carried out.Several resolution techniques for the 1DV horizontal (i.e., in the x-direction perpendicular to coastline) momentum equation governing the Mean Horizontal Velocity (MHV) are analysed. In the first kind of techniques, the mean horizontal velocity is computed from the momentum equation, whereas the Mean Water Level (MWL) is computed using a parametrization of the depth-averaged momentum equation. Two boundary or integral conditions are thus needed. In the second kind, both mean horizontal velocity and mean water level gradient in the x-direction are the unknowns of the momentum equation, thus, three boundary or integral conditions are needed. Various additional conditions are discussed. We show that using a technique of the first kind is equivalent to imposing the difference between the surface and the bottom shear stresses in the 1D vertical equation. Both techniques lead to results that are in good agreement with the Delta Flume experimental data, provided the Stokes drift flow discharge is imposed as an additional condition. The influence of the breaking roller model and of the turbulent viscosity parametrization are also analysed.Suspended sediment transport by the mean current and wave-induced bedload transport are taken into account in the sediment flux. Three turbulent diffusivity parametrizations are compared for suspended sediment concentration estimations. A linear profile for the turbulent diffusivity taking into account the wave bottom shear stress and the surface wave breaking turbulence production is shown to give the best results. Using experimental data, we put forward the poor estimation of the bottom sediment concentration given by the three implemented parametrizations. We thus propose a new parametrization relying on a Shields parameter based on the breaking roller induced surface shear stress. Using this new parametrization, the bottom profile used in the tests keeps its two bars which disappear otherwise. However, the morphodynamical model still overestimates the bars offshore motion, a bias already observed in other models. 相似文献
4.
A static adaptive grid approximates the topography and defines the vertical resolution in Vector-Ocean-Model (VOM). The adaptation to topography creates unstructured grids, which are organised in a one-dimensional vector by column-wise storage of only wet cells. The model’s name reflects this data structure. The intention of VOM is better resolving flow and stratification near topographic boundaries in Z-coordinates. This is the second part of a publication that describes the generation of adaptive grids (part I), and simulations with VOM in unstructured grids (this part). Adaptive grids generated for a synthetic topography in part I include shelf, continental slope, and ocean. Three of those grids are here utilised in upwelling simulations. Under the same forcing increased vertical resolution at seabed and slopes yields a significant increase in flow energy as compared to coarser grids. Results allow explaining the surface intensification of a continental slope jet by vertical displacements of water masses in the seabed Ekman layer. Results in unstructured grids are almost identical to reference simulations in equidistant grids where the respective smallest grid size of unstructured grids was used. Negative effects of grids on predicted flow and stratification are absent also over particularly rough topography, as demonstrated by using vertical velocity as most sensitive indicator. In a further simulation an overflow governed by the advection of water mass properties is presented to demonstrate the conservation properties of the model. After 5 months of simulation the predicted domain average temperature deviated by 10−8 from the initial temperature field. Compared to equidistant grids the advection/diffusion scheme looses about one order of magnitude in accuracy when used in an unstructured grid. The results of VOM, being defined in Z-coordinates, are void of coordinate transformation errors. In an arbitrary topography unforced zero-flow remains quiescent in a stratification that only varies in the vertical. VOM due to its depth-independent vertical resolution appears particularly suitable for simulations of ocean-shelf exchange. 相似文献
5.
MPAS-Ocean is used to simulate an idealized, density-driven overflow using the dynamics of overflow mixing and entrainment (DOME) setup. Numerical simulations are carried out using three of the vertical coordinate types available in MPAS-Ocean, including z-star with partial bottom cells, z-star with full cells, and sigma coordinates. The results are first benchmarked against other models, including the MITgcm’s z-coordinate model and HIM’s isopycnal coordinate model, which are used to set the base case used for this work. A full parameter study is presented that looks at how sensitive overflow simulations are to vertical grid type, resolution, and viscosity. Horizontal resolutions with 50 km grid cells are under-resolved and produce poor results, regardless of other parameter settings. Vertical grids ranging in thickness from 15 m to 120 m were tested. A horizontal resolution of 10 km and a vertical resolution of 60 m are sufficient to resolve the mesoscale dynamics of the DOME configuration, which mimics real-world overflow parameters. Mixing and final buoyancy are least sensitive to horizontal viscosity, but strongly sensitive to vertical viscosity. This suggests that vertical viscosity could be adjusted in overflow water formation regions to influence mixing and product water characteristics. Lastly, the study shows that sigma coordinates produce much less mixing than z-type coordinates, resulting in heavier plumes that go further down slope. Sigma coordinates are less sensitive to changes in resolution but as sensitive to vertical viscosity compared to z-coordinates. 相似文献
6.
Ｎｕｍｅｒｉｃａｌｓｉｍｕｌａｔｉｏｎｏｆｔｈｒｅｅ－ｄｉｍｅｎｓｉｏｎａｌｔｉｄａｌｃｕｒｒｅｎｔｉｎｔｈｅＢｏｂａｌＳｅａ￥ＤｅｎＺｈｅｎｘｉｎｇ;ＹａｎｇＬｉａｎｗｕ;ａｎｄＪｏｓｅＯｚｅｒ（ＲｅｃｅｉｖｅｄＭａｙ１,１９９３;ａｃｃｅｐｔｅｄ... 相似文献
7.
研究建立渤海海域三维斜压场的湍流能量模型。对计算网格无法分辩的湍流运动，引入湍流能量求得垂向动量、质量及能量交换系统，从而形成封闭的湍流运动方程组。定量地计算了湍流能量、湍流速度及有风、无风状况下湍流能量的垂直分布等，探讨了湍流运动对潮波运动的影响。 相似文献
8.
The impact of vertical resolution on the evolution and movement of tropical cyclones was studied using NCAR MM5 model with a horizontal resolution of 9 km. Four numerical experiments were performed with different vertical resolutions, that is, with 23 vertical levels as control experiment, and 36 vertical levels with high resolution in the lower troposphere, 33 vertical levels with high resolution in the upper troposphere and 46 vertical levels with increased vertical resolution throughout the troposphere as relative to base experiment. The results indicate that increased vertical resolution in the lower troposphere produces efficient intensification and better structure in terms of eye and eyewall. Increased vertical resolution at lower levels improves the prediction of vertical shear of horizontal wind. Experiments with high resolution in the lower troposphere and high resolution throughout the troposphere simulate better track up to 72 hours. 相似文献
9.
CHEN Yongping 《中国海洋工程》2001,(2):241-252
Based on the fluid motion equations, the physical meaning of eddy viscosity coefficient and the rationality of the Boussinesq hypothesis are discussed in this paper. The effect of the coefficient on numerical stability is analyzed briefly. A semi-enclosed rectangular sea area, with an orthogonal spur dike, is applied in a 2-D numerical model to study the effect of horizontal eddy viscosity coefficient (AH). The computed result shows that AH has little influence on the tidal level and averaged flow velocity, but has obvious influence on the intensity and the range of return flow around near the spur dike. Correspondingly, a wind-driven current pool and an annular current are applied in a 3-D numerical model respectively to study the effect of vertical eddy viscosity coefficient (Av). The computed result shows that the absolute value of Av is inversely proportional to that of horizontal velocity, and the vertical gradient value of Av determines the vertical distribution of horizontal velocity. The distrib 相似文献
10.
The numerical investigation of random wave slamming on superstructures of marine structures in the splash zone is presented in this paper. The impact pressures on the underside of the structure are computed based on the improved volume of fluid method (VOF). The governing equations are Reynolds time-averaged equations and the two equation k– model. The third order upwind difference scheme is applied to the convection term to reduce the effect of numerical viscosity. The numerical wave flume with random wave-maker suitable for VOF is established. Appropriate moving contact-line boundary conditions are introduced to the model wave in contact with and separated from the underside of structure. Parametric studies have been carried out for different incident waves, structure dimensions and structure clearance. The numerical results are verified by the experimental results. 相似文献
11.
S. N. Moshonkin V. B. Zalesny A. V. Gusev R. Tamsalu 《Izvestiya Atmospheric and Oceanic Physics》2014,50(1):49-60
A physical formulation of the problem is considered. A mathematical model and a numerical algorithm of the turbulence model as part of the ocean circulation model for simulations for decades are formulated. The model is based on the evolution equations for turbulent kinetic energy (TKE) and the frequency of its viscous dissipation. A numerical solution algorithm for both the circulation model and the turbulence model is based on implicit schemes of splitting with respect to physical processes and geometric coordinates. For the turbulence model, this provided analytical solutions at a splitting step related to TKE generation and dissipation. Numerical experiments have been performed with a model of the joint circulation of the North Atlantic, the Arctic Ocean, and the Bering Sea to reproduce the annual cycle and synoptic disturbances of ocean characteristics. The model has a resolution of 0.25° in latitude and longitude and 40 levels in the vertical, which are compressed toward the surface to reproduce the process of developed turbulence better. The results are compared with observations and with the results of simulations using traditional parameterizations of the upper ocean mixing. It is shown that the model reproduces ocean characteristics correctly, only slightly increasing the computation time in comparison with simple parameterizations. Spatial and temporal characteristics of turbulence are analyzed. 相似文献
12.
《Deep Sea Research Part II: Topical Studies in Oceanography》1999,46(6-7):1169-1198
Simulations from a coupled ice–ocean model that highlight the importance of synoptic forcing on sea-ice dynamics are described. The ocean model is a non-hydrostatic primitive equation model coupled to a dynamic thermodynamic sea ice model. The ice modelling sensitivity study presented here is part of an ongoing research programme to define the role played by sea ice in the energy balance of the Greenland Sea. The different categories of sea ice found in the subpolar regions are simulated through the use of equations for thin ice, thick ice and the Marginal Ice Zone. A basin scale numerical model of the Greenland, Iceland and Norwegian Seas has a horizontal resolution of 20 km and a vertical grid spacing of 50 m. This resolution is adequate for resolving the mesoscale topographic structures known to control the circulation in this region. The spin-up reproduces the main features of the circulation, including the cyclonic gyres in the Norwegian and Greenland Basins and Iceland Plateau. Topographic steering of the flow is evident. The baroclinic Rossby radius of deformation is between 5 and 10 km so that the model is not eddy-resolving. The coupled ice–ocean model was run for a period of two weeks. The influence of horizontal resolution of the atmospheric model was tested by comparing simulations using six hourly wind fields from the ECMWF with those generated using six hourly fields from a HIRLAM, with horizontal resolutions of 1° and 0.18° respectively. The simulations show reasonable agreement with satellite ice compactness data and data of ice transports across sections at 79°N, 75°N and Denmark Strait. 相似文献
13.
《Coastal Engineering》2006,53(4):363-379
Low crested structures for beach defence purposes are frequently configured in the form of cells. This paper investigates the hydrodynamic and morphological response of such defence schemes by means of laboratory experiments and numerical simulations.Experiments were carried out at the large mobile-bed basin of the “Laboratorio di Idraulica Costiera” of the “Politecnico di Bari”, IT. Tests peculiarity is the simultaneous simulation of three different types of lateral confinement: a channel, representative of an indefinitely long structure, and two areas with narrow and wide gaps, representative of a common scheme and of an almost isolated structure respectively. Submerged, zero freeboard and emerged conditions were tested. A typical storm of the Northern Adriatic sea was reproduced by 6 wave attacks.Water level and currents were measured in front and behind the structures, at gaps and roundheads. Results, focusing on rip current intensities and piling-up, are presented and discussed. An extensive plan view of wave intensities and flow patterns inside the basin is provided through numerical simulations with MIKE 21 numerical suite.The sand bed was surveyed along 30 profiles covering the narrow gap area and the isolated structure roundheads. For the submerged conditions tests, 15 of these profiles and the shoreline position were monitored after each step of the simulated storm.Average and maximum values of local erosion at the structure toe are presented. Time evolution of erosion at gaps is related to the main hydrodynamic variables through Bijker sediment transport formula. Evaluation of volume changes in the protected area shows that sand is in average trapped behind the barriers, being transported offshore from gaps and returning inshore over the crest. The observed and reconstructed transport mechanism is a combination of bed load and suspended transport, the latter being slightly more important.Bed level changes are well predicted by morphodynamic simulations performed with MIKE 21 CAMS. 相似文献
14.
基于微控制器(MCU)设计了一个通用的四维混沌系统数字硬件实验电路,由此实现了9×7网格涡卷的混沌和超混沌吸引子的生成.本文基于由Colpitts振荡器模型延伸出的四维多涡卷超混沌系统,通过引入单位锯齿波函数替换原系统中的三角波函数,构建了一个便于MCU数字硬件实现的新的网格涡卷超混沌系统,并对新系统网格涡卷吸引子的形成机理进行了分析和数值仿真.通过采用Euler算法对新系统进行离散化,在实验电路的有效动态范围内可以生成比原系统更多网格涡卷数量的吸引子.实验结果有效验证了本文基于MCU实现的网格涡卷超混沌 相似文献
15.
A mechanistic study is performed to examine the coupling between the in-line and the cross-flow motion of a cylindrical structure subjected to current forces. The structure represents a free span pipeline but concerns marine risers as well.A time domain model is formulated in which the in-line and cross-flow deflections are coupled through the axial tension which in turn is computed from the pipeline prolongation at any time. This formulation introduces time dependent tensions and non-linearity into the problem.Preliminary validation of the model simulations vs. physical test data are carried out for one specific case to ensure that the sag and the in-line deflection are correctly resolved by the model. Using this as the initial condition a series of calculations are carried out to examine cross-flow induced deflections induced by an in-line prescribed deflection and vice versa. Finally, an idealistic simulation of flow induced vibration is presented.The model simulations demonstrate that the coupling varies with the mode shape and with which component it is initially introduced into. However, it is evident that the coupling effects may be significant and not negligible. 相似文献
16.
ＡｎＡＤＩｔｅｃｈｎｉｑｕｅｆｏｒｓｏｌｖｉｎｇｔｈｒｅｅ－ｄｉｍｅｎｓｉｏｎａｌｃｏａｓｔａｌｃｉｒｃｕｌａｔｉｏｎＹａｎｇＬｉａｎｗｕ１，＊ａｎｄＪｏｓｅＯｚｅｒ２（ＲｅｃｅｉｖｅｄＮｏｖｅｍｂｅｒ２７，１９９６；ａｃｃｅｐｔｅｄＤｅｃｅｍｂｅｒ... 相似文献
17.
Benchmark computations of wave run-up on single cylinder and four cylinders by naoe-FOAM-SJTU solver
The benchmark simulations of wave run-up on a fixed single truncated circular cylinder and four circular cylinders are presented in this paper. Our in-house CFD solver naoe-FOAM-SJTU is adopted which is an unsteady two-phase CFD code based on the open source package OpenFOAM. The Navier-Stokes equations are employed as the governing equations, and the volume of fluid (VOF) method is applied for capturing the free surface. Monochromatic incident waves with the specified wave period and wave height are simulated and wave run-up heights around the cylinder are computed and recorded with numerical virtual wave probes. The relationship between the wave run-up heights and the incident wave parameters are analyzed. The numerical results indicate that the presented naoe-FOAM-SJTU solver can provide accurate predictions for the wave run-up on one fixed cylinder and four cylinders, which has been proved by the comparison of simulated results with experimental data. 相似文献
18.
This is the first part of a publication that describes the generation of adaptive grids (this part), and simulations with vector-ocean-model (VOM) in unstructured grids resulting from the adaptation (part II). A static vertical adaptive grid in z-coordinates allows improving the approximation of topography and vertical resolution at slopes. Adaptive grids use elements from a set of grid sizes by multiplying a basic smallest cell size with powers of two, as in cell division. Grids with locally isotropic vertical resolution at surface, seabed, and slopes can be generated whereby resolution decreases towards the ocean interior. The adaptation to topography yields unstructured grids that are organised in a one-dimensional vector by column-wise storage of cells, discarding land cells. The vector storage suggested the model’s name. Grids are generated by an iterative procedure that relies on rules, i.e. criteria and directives to control the grid structure in favour of a good representation of physics and smooth numerical operations. The directives govern vertical resolution at sea surface and seabed, and at slopes. For the latter vertical resolution is extended in the horizontal. In the ocean interior horizontal distances between changes in grid size can be controlled for the sake of smooth numerics. The use of a z-grid that avoids transformation errors, the depth-independence of vertical resolution, and the lateral extension of vertical resolution at slopes towards the ocean interior are the most significant differences of adaptive grids in comparison to vertical coordinate transformations. Unstructured grids do not rely on a smoothing of topography and can be used within any of the horizontal Arakawa-grids. For the same topography directives allow creating various grids as demonstrated for a shelf-ocean topography. The number of cells per column in two unstructured grids generated for the North Atlantic may locally well exceed typical layer numbers in conventional model matrices. But the domain average is similar to layer numbers of today’s ocean models. Thus, with the same investment of cells per domain a higher resolution in slope regions can be achieved by unstructured grids as compared to conventional z-grids. 相似文献
19.
利用热带太平洋TAO浮标阵列评估全球高分辨率浪-潮-流耦合模式同化结果 总被引：3，自引：0，他引：3
In order to evaluate the assimilation results from a global high resolution ocean model, the buoy observations from tropical atmosphere ocean(TAO) during August 2014 to July 2015 are employed. The horizontal resolution of wave-tide-circulation coupled ocean model developed by The First Institute of Oceanography(FIOCOM model) is 0.1°×0.1°, and ensemble adjustment Kalman filter is used to assimilate the sea surface temperature(SST), sea level anomaly(SLA) and Argo temperature/salinity profiles. The simulation results with and without data assimilation are examined. First, the overall statistic errors of model results are analyzed. The scatter diagrams of model simulations versus observations and corresponding error probability density distribution show that the errors of all the observed variables, including the temperature, isotherm depth of 20°C(D20), salinity and two horizontal component of velocity are reduced to some extent with a maximum improvement of 54% after assimilation. Second, time-averaged variables are used to investigate the horizontal and vertical structures of the model results. Owing to the data assimilation, the biases of the time-averaged distribution are reduced more than70% for the temperature and D20 especially in the eastern Pacific. The obvious improvement of D20 which represents the upper mixed layer depth indicates that the structure of the temperature after the data assimilation becomes more close to the reality and the vertical structure of the upper ocean becomes more reasonable. At last,the physical processes of time series are compared with observations. The time evolution processes of all variables after the data assimilation are more consistent with the observations. The temperature bias and RMSE of D20 are reduced by 76% and 56% respectively with the data assimilation. More events during this period are also reproduced after the data assimilation. Under the condition of strong 2014/2016 El Ni?o, the Equatorial Undercurrent(EUC) from the TAO is gradually increased during August to November in 2014, and followed by a decreasing process. Since the improvement of the structure in the upper ocean, these events of the EUC can be clearly found in the assimilation results. In conclusion, the data assimilation in this global high resolution model has successfully reduced the model biases and improved the structures of the upper ocean, and the physical processes in reality can be well produced. 相似文献
20.
迎风有限元法在三维潮流数值模拟中的应用 总被引：4，自引：0，他引：4
为研究和开发迎风有限元法在河口海岸三维波数值模拟中的应用，并为建立三维物质扩散模型，最终为解决实际工程问题打下基础，以非线性浅水波方程为基础，采用一咱平面迎风有奶元与垂向工有限差分相结合地数值计算方法建立三维潮流数学模型，基地实际的物理过程，在计算中将三维流动分成外重力波和内重力波，耦联求解潮位和流速的空间分布。应用本模式计算了天津新港附近渤海海域的三维潮流运动，以较高的分辩率揭示了潮流空间结构特 相似文献