共查询到20条相似文献，搜索用时 31 毫秒
1.
In this paper, a feasibility of anelastic approach for numerical weather prediction (NWP) is examined. The study concerns
the anelastic nonhydrostatic model EULAG as a prospective candidate for the new dynamical core of a high-resolution NWP model.
Such an application requires a series of benchmark tests to be performed. The study presents the results of dry idealized
two-dimensional linear and non-linear tests. They include evolution of cold and warm density currents in neutrally stratified
atmosphere, inertia-gravity waves in short and long channels, as well as mountain gravity waves for a set of different flow
regimes. Detailed comparison of the results with the reference solutions, based mainly on the results of compressible models,
indicates a high level of conformity for all of the experiments. It verifies the anelastic approach as strongly consistent
with the compressible one for a broad class of atmospheric problems. It also corroborates the robustness of EULAG numerics,
an essential requirement of dynamical core of NWP model. 相似文献
2.
Zbigniew P. Piotrowski Andrzej A. Wyszogrodzki Piotr K. Smolarkiewicz 《Acta Geophysica》2011,59(6):1294-1311
This paper highlights progress with the development of a petascale implementation of general-purpose high-resolution (nonoscillatory)
hydrodynamical simulation code EULAG [Prusa et al. 2008, Comput. Fluids
37, 1193]. The applications addressed are anelastic atmospheric flows in the range of scales from micro to planetary. The new
modeldomain decomposition into a three dimensional processor array has been implemented to increase model performance and
scalability. The performance of the new code is demonstrated on the IBM BlueGene/L and Cray XT4/XT5 supercomputers. The results
show significant improvement of the model efficacy compared to the original decomposition into a two-dimensional processor
array in the horizontal — a standard in meteorological models. 相似文献
3.
The anelastic nonhydrostatic model EULAG is a candidate for the future dynamical core of a numerical weather prediction model.
Achieving such an objective requires a number of experiments focused on testing correctness of the solutions and robustness
of the solver. In the spirit of this idea, a set of tests related to standard atmospheric problems was performed, of which
the two regarding development and evolution of a supercell were employed as benchmarks of moist dynamics of the model. Their
results are discussed in this paper. Development and evolution of a stormsystem with a set of characteristic features such
as stormsplitting along with the generation of horizontal vorticity and cold pool formation is investigated. In addition,
the influence of domain geometry, boundary conditions and subgrid-scale mixing is examined. 相似文献
4.
EULAG is a computational model for simulating flows across a wide range of scales and physical scenarios. A standard option
employs an anelastic approximation to capture nonhydrostatic effects and simultaneously filter sound waves from the solution.
In this study, we examine a localized gravity wave packet generated by instabilities in Held-Suarez climates. Although still
simplified versus the Earth’s atmosphere, a rich set of planetary wave instabilities and ensuing radiated gravity waves can arise. Wave packets
are observed that have lifetimes ≤ 2 days, are negligibly impacted by Coriolis force, and do not show the rotational effects
of differential jet advection typical of inertia-gravity waves. Linear modal analysis shows that wavelength, period, and phase
speed fit the dispersion equation to within a mean difference of ∼ 4%, suggesting an excellent fit. However, the group velocities
match poorly even though a propagation of uncertainty analysis indicates that they should be predicted as well as the phase
velocities. Theoretical arguments suggest the discrepancy is due to nonlinearity — a strong southerly flow leads to a critical
surface forming to the southwest of the wave packet that prevents the expected propagation. 相似文献
5.
Herman W. J. Kernkamp Henri A. H. Petit Herman Gerritsen Erik D. de Goede 《Ocean Dynamics》2005,55(3-4):351-369
In this paper, the formulations of the primitive equations for shallow water flow in various horizontal co-ordinate systems
and the associated finite difference grid options used in shallow water flow modelling are reviewed. It is observed that horizontal
co-ordinate transformations do not affect the chosen co-ordinate system and representation in the vertical, and are the same
for the three- and two-dimensional cases. A systematic derivation of the equations in tensor notation is presented, resulting
in a unified formulation for the shallow water equations that covers all orthogonal horizontal grid types of practical interest.
This includes spherical curvilinear orthogonal co-ordinate systems on the globe. Computational efficiency can be achieved
in a single computer code. Furthermore, a single numerical algorithmic code implementation satisfies. All co-ordinate system
specific metrics are determined as part of a computer-aided model grid design, which supports all four orthogonal grid types.
Existing intuitive grid design and visual interpretation is conserved by appropriate conformal mappings, which conserve spherical
orthogonality in planar representation. A spherical curvilinear co-ordinate solution of wind driven steady channel flow applying
a strongly distorted grid is shown to give good agreement with a regular spherical co-ordinate model approach and the solution
based on a β-plane approximation. Especially designed spherical curvilinear boundary fitted model grids are shown for typhoon
surge propagation in the South China Sea and for ocean-driven flows through Malacca Straits. By using spherical curvilinear
grids the number of grid points in these single model grid applications is reduced by a factor of 50–100 in comparison with
regular spherical grids that have the same horizontal resolution in the area of interest. The spherical curvilinear approach
combines the advantages of the various grid approaches, while the overall computational effort remains acceptable for very
large model domains. 相似文献
6.
Three-dimensional grids representing a heterogeneous, ground water system are generated at 10 different resolutions in support of a site-scale flow and transport modeling effort. These grids represent hydrostratigraphy near Yucca Mountain, Nevada, consisting of 18 stratigraphic units with contrasting fluid flow and transport properties. The grid generation method allows the stratigraphy to be modeled by numerical grids of different resolution so that comparison studies can be performed to test for grid quality and determine the resolution required to resolve geologic structure and physical processes such as fluid flow and solute transport. The process of generating numerical grids with appropriate property distributions from geologic conceptual models is automated, thus making the entire process easy to implement with fewer user-induced errors. The series of grids of various resolutions are used to assess the level at which increasing resolution no longer influences the flow and solute transport results. Grid resolution is found to be a critical issue for ground water flow and solute transport. The resolution required in a particular instance is a function of the feature size of the model, the intrinsic properties of materials, the specific physics of the problem, and boundary conditions. The asymptotic nature of results related to flow and transport indicate that for a hydrologic model of the heterogeneous hydrostratigraphy under Yucca Mountain, a horizontal grid spacing of 600 m and vertical grid spacing of 40 m resolve the hydrostratigraphic model with sufficient precision to accurately model the hypothetical flow and solute transport to within 5% of the value that would be obtained with much higher resolution. 相似文献
7.
The nonhydrostatic pressure effects on the generation and propagation of wind-forced internal waves are studied with a two-dimensional
numerical ocean model. A one-way directed wind pulse over a stratified ocean initiates surface and internal waves in a closed
basin. The studies are performed with horizontal grid sizes in the range from 1 km to 62.5 m. The experiments are performed
with both a hydrostatic and a nonhydrostatic model, facilitating systematic studies of the sensitivity of the numerical model
results to the grid size and to the nonhydrostatic pressure adjustments. The results show that the nonhydrostatic pressure
effects are highly dependent on the grid size and grow with increased resolution. In the internal depression wave, the horizontal
nonhydrostatic pressure gradients reach the same order of magnitude as the hydrostatic gradients in the high-resolution nonhydrostatic
studies. In these studies, the nonhydrostatic pressure gradients approximately balance the corresponding hydrostatic pressure
gradients in the internal depression wave, and the wave degenerates into a train of soliton waves. The time for the soliton
form to develop agrees with the steepening timescale calculated from Korteweg-de Vries theory. In the high-resolution hydrostatic
model, the internal depression wave takes the form of a single wave front. When the internal waves are generated in the boundary
layers, the nonhydrostatic pressure gradients are much smaller than the hydrostatic gradients and the generation processes
are not effected by the nonhydrostatic pressure with the present range of grid sizes. 相似文献
8.
In this paper we explore the optimum assimilation of high‐resolution data into numerical models using the example of topographic data provision for flood inundation simulation. First, we explore problems with current assimilation methods in which numerical grids are generated independent of topography. These include possible loss of significant length scales of topographic information, poor representation of the original surface and data redundancy. These are resolved through the development of a processing chain consisting of: (i) assessment of significant length scales of variation in the input data sets; (ii) determination of significant points within the data set; (iii) translation of these into a conforming model discretization that preserves solution quality for a given numerical solver; and (iv) incorporation of otherwise redundant sub‐grid data into the model in a computationally efficient manner. This processing chain is used to develop an optimal finite element discretization for a 12 km reach of the River Stour in Dorset, UK, for which a high‐resolution topographic data set derived from airborne laser altimetry (LiDAR) was available. For this reach, three simulations of a 1 in 4 year flood event were conducted: a control simulation with a mesh developed independent of topography, a simulation with a topographically optimum mesh, and a further simulation with the topographically optimum mesh incorporating the sub‐grid topographic data within a correction algorithm for dynamic wetting and drying in fixed grid models. The topographically optimum model is shown to represent better the ‘raw’ topographic data set and that differences between this surface and the control are hydraulically significant. Incorporation of sub‐grid topographic data has a less marked impact than getting the explicit hydraulic calculation correct, but still leads to important differences in model behaviour. The paper highlights the need for better validation data capable of discriminating between these competing approaches and begins to indicate what the characteristics of such a data set should be. More generally, the techniques developed here should prove useful for any data set where the resolution exceeds that of the model in which it is to be used. Copyright © 2002 John Wiley & Sons, Ltd. 相似文献
9.
《Advances in water resources》1998,22(2):133-143
In modeling shocks in open channel flows, the traditional finite difference schemes become inefficient and warrant special numerical treatment for smooth computations. This paper provides a general introduction to the non-oscillatory high-resolution methodology, coupled with the advantages of using these conservative methods for open channel applications. Results of the numerical experiments are presented and compared to standard algorithm results, analytical solutions and experimental results to demonstrate the robustness of the high-resolution formulation. Finally results of flows with mixed flow conditions (as in hydraulic jump), indicate that the contribution of the Boussinesq pressure term is minimal. Besides showing robustness in capturing shocks in open channel, the results indicate that the effect of the grid spacing on the shock resolution is small. However, for jumps with high Froude numbers, the present formulation slightly underestimates the peak depth. 相似文献
10.
We propose an improvement of the overland‐flow parameterization in a distributed hydrological model, which uses a constant horizontal grid resolution and employs the kinematic wave approximation for both hillslope and river channel flow. The standard parameterization lacks any channel flow characteristics for rivers, which results in reduced river flow velocities for streams narrower than the horizontal grid resolution. Moreover, the surface areas, through which these wider model rivers may exchange water with the subsurface, are larger than the real river channels potentially leading to unrealistic vertical flows. We propose an approximation of the subscale channel flow by scaling Manning's roughness in the kinematic wave formulation via a relationship between river width and grid cell size, following a simplified version of the Barré de Saint‐Venant equations (Manning–Strickler equations). The too large exchange areas between model rivers and the subsurface are compensated by a grid resolution‐dependent scaling of the infiltration/exfiltration rate across river beds. We test both scaling approaches in the integrated hydrological model ParFlow. An empirical relation is used for estimating the true river width from the mean annual discharge. Our simulations show that the scaling of the roughness coefficient and the hydraulic conductivity effectively corrects overland flow velocities calculated on the coarse grid leading to a better representation of flood waves in the river channels. 相似文献
11.
Numerical studies of flow over a sill: sensitivity of the non-hydrostatic effects to the grid size 总被引：1，自引：1，他引：0
A non-hydrostatic terrain-following model in cross sectional form is applied to study the processes in the lee of a sill in
an idealized stratified fjord during super-critical tidal inflow. A sequence of numerical studies with horizontal grid sizes
in the range from 100 to 1.5625 m are performed. All experiments are repeated using both hydrostatic and non-hydrostatic versions
of the model, allowing a systematic study of possible non-hydrostatic pressure effects and also of the sensitivity of these
effects to the horizontal grid size. The length scales and periods of the internal waves in the lee of the sill are gradually
reduced and the amplitudes of these waves are increased as the grid size is reduced from 100 down to 12.5 m. With a further
reduction in grid size, more short time and space scale motions become superimposed on the internal waves. Associated with
the internal wave activity, there is a deeper separation point that is fairly robust to all parameters investigated. Another
separation point nearer to the top of the sill appears in the numerical results from the high-resolution studies with the
non-hydrostatic model. Associated with this shallower separation point, an overturning vortex appears in the same set of numerical
solutions. This vortex grows in strength with reduced grid size in the non-hydrostatic experiments. The effects of the non-hydrostatic
pressure on the velocity and temperature fields grow with reduced grid size. In the experiments with horizontal grid sizes
equal to 100 or 50 m, the non-hydrostatic pressure effects are small. For smaller grid sizes, the time mean velocity and temperature
fields are also clearly affected by the non-hydrostatic pressure adjustments. 相似文献
12.
A non-linear two-dimensional vertically stratified cross-sectional model of a constant depth basin without rotation is used
to investigate the influence of vertical and horizontal diffusion upon the wind-driven circulation in the basin and the associated
temperature field. The influence of horizontal grid resolution, in particular the application of an irregular grid with high
resolution in the coastal boundary layer is examined. The calculations show that the initial response to a wind impulse is
downwelling at the downwind end of the basin with upwelling and convective mixing at the opposite end. Results from a two-layer
analytical model show that the initial response is the excitation of an infinite number of internal seiche modes in order
to represent the initial response which is confined to a narrow near coastal region. As time progresses, at the downwind end
of the basin a density front propagates away from the boundary, with the intensity of its horizontal gradient and associated
vertical velocity determined by both horizontal and vertical viscosity values. Calculations demonstrate the importance of
high horizontal grid resolution in resolving this density gradient together with an accurate density advection scheme. The
application of an irregular grid in the horizontal with high grid resolution in the nearshore region enables the initial response
to be accurately reproduced although physically unrealistic short waves appear as the frontal region propagates onto the coarser
grid. Parameterization of horizontal viscosity using a Smagorinsky-type formulation acts as a selective grid size-dependent
filter, and removes the short-wave problem although enhanced smoothing can occur if the scaling coefficient in the formulation
is too large. Calculations clearly show the advantages of using an irregular grid but also the importance of using a grid
size-dependent filter to avoid numerical problems. 相似文献
13.
L.-Y. Oey 《Continental Shelf Research》1997,17(15):1929-1944
Observations indicate that, perhaps as a result of intense shears produced by the topographic steering of the slope current and merging of the warm North Atlantic water with cold Nordic waters, the Faroe-Shetland Channel is a fertile ground for energetic eddies and meanders. How these meso-scale ( 10 km) features are resolved in a three-dimensional numerical model is studied here with grid-nesting, in which the convergence of the simulated energetics is examined as a function of grid sizes Δ =10, 10/3, 10/5 and 10/7 km. Large-scale forcing to the fine-grid nest consists of steady slope-current transport (2 Sv) and north/south climatological density difference only. The mean fields were found to converge at Δ = 10/3 km, but the eddy fields required Δ = 2 km or less. At the finest grids, meso-scale features more in conformity with those observed were found: meanders with amplitudes of O(10 km) and along-slope propagation speeds of about 0.14 m s- 1 cross-channel vacillation of the slope current, with amplitude of about 10 km and periods of 5–10 days, and vertical excursions of isotherms of O(100 m). The meanders were highly ageostrophic (ζ/f> 1, ζ=relative vorticity, f = Coriolis parameter), and developed as a result of conditions favorable to flow instability originated from south of the channel. When adequately resolved, the meander energetics are comparable to those produced by (large-scale) winds, which suggests that a high-resolution grid (Δ < 2 km) will be required for hindcast and forecast studies. 相似文献
14.
This study presents an experimental analysis from aircraft measurements above the Pyrenees chain during the PYREX experiment. The Pyrenees chain, roughly WE oriented, is a major barrier for northerly and southerly airflows. We present a case of southerly flow (15 October 1990) and three successive cases of northerly flows above the Pyrenees (14, 15 and 16 November 1990) documented by two aircraft. The aircraft have described a vertical cross section perpendicular to the Pyrenean ridge. This area is described via the thermodynamical and dynamical fields which have a horizontal resolution of 10 km. Three methods for computing the vertical velocity of the air are presented. The horizontal advection terms which play a role in the budget equations are also evaluated. The altitude turbulence zone of 15 October are shown via turbulent fluxes, turbulent kinetic energy (TKE), dissipation rate of TKE and inertial length-scale. A comparison of results obtained by eddy-correlation and inertial-dissi-pation method is presented. The experimental results show a warm and dry downdraft for the southerly flow with large values for advection terms. All the mountain wave cases are also shown to present an important dynamical perturbation just above the Pyrenees at upper altitudes. 相似文献
15.
A high resolution (3–8 km grid), 3D numerical ocean model of the West Caribbean Sea (WCS) is used to investigate the variability
and the forcing of flows near the Meso-American Barrier Reef System (MBRS) which runs along the coasts of Mexico, Belize,
Guatemala and Honduras. Mesoscale variations in velocity and temperature along the reef were found in seasonal model simulations
and in observations; these variations are associated with meandering of the Caribbean current (CC) and the propagation of
Caribbean eddies. Diagnostic calculations and a simple assimilation technique are combined to infer the dynamically adjusted
flow associated with particular eddies. The results demonstrate that when a cyclonic eddy (negative sea surface height anomaly
(SSHA)) is found near the MBRS the CC shifts offshore, the cyclonic circulation in the Gulf of Honduras (GOH) intensifies,
and a strong southward flow results along the reef. However, when an anticyclonic eddy (positive SSHA) is found near the reef,
the CC moves onshore and the flow is predominantly westward across the reef. The model results help to explain how drifters
are able to propagate in a direction opposite to the mean circulation when eddies cause a reversal of the coastal circulation.
The effect of including the Meso-American Lagoon west of the Belize Reef in the model topography was also investigated, to
show the importance of having accurate coastal topography in determining the variations of transports across the MBRS. The
variations found in transports across the MBRS (on seasonal and mesoscale time scales) may have important consequences for
biological activities along the reef such as spawning aggregations; better understanding the nature of these variations will
help ongoing efforts in coral reef conservation and maintaining the health of the ecosystem in the region. 相似文献
16.
Seismic reflection methods provide continuous access both to stratigraphy (vertical) and to subsurface morphology (horizontal), for which the scales of interest may differ by orders of magnitude. Seismic surveys of Quaternary successions have generally sought to optimise vertical resolution, through the use of higher source frequency content. Here, I show that low-frequency bandwidth is not necessarily a limiting factor for the seismic resolution of glacigenic morpho-sedimentary features. Observations are presented from a buried network of large mid-Pleistocene (Elsterian) tunnel-valleys in the southern North Sea Basin, across a 100×130 km study area with water depths less than 30 m. Low-frequency 2D and 3D seismic multi-channel data, acquired for deeper hydrocarbon exploration, are compared with previously available high-frequency single- and multi-channel profiles (5–15 km grid spacing). The low-frequency data contribute to a new understanding of the basal morphology and fill stratigraphy of the tunnel-valleys, in part due to higher data densities (≥1 km grid spacing), but also to improved imaging of reflectors at depth. The tunnel-valleys are seen to be overdeepened troughs, shallow (≤0.5 km) relative to their widths (≤6 km). The basal unconformity defines a series of arborescent elements, convergent to the south; erosional overlap by younger elements to the north has resulted in anastomosing patterns in places. The fill is dominated by axially downlapping clinoforms, descending to the north, onlapped and overlain by subhorizontal reflectors. Well data show that sand-dominated glaciofluvial sediments are overlain by glaciolacustrine to marine muds. Better definition of the clinoforms on low- versus high-frequency multi-channel data is suggested to reflect the coarse spatial scale of the backset glaciofluvial strata. The results support a simple interpretation of time-transgressive tunnel-valley formation by coeval glaciofluvial erosion and backfill beneath the outer tens of kilometres of the northward receding Elsterian ice sheet margin. Comparable submarginal interpretations have been proposed for drainage features (tunnel-valleys and eskers) of the last deglaciation of both northern Europe and North America using integrated geomorphologic and stratigraphic methods. Commercial 2D and 3D seismic data are widely available from exploration areas such as the North Sea and are argued to constitute an underexploited resource for Quaternary research. 相似文献
17.
《Advances in water resources》1998,21(1):1-9
The unsteady free surface flow caused by sudden collapse of a dam produces discontinuities in the flow variables. As the flow surges downstream, it forms a moving bore front with steep gradients of water height and velocity. In the numerical simulation of this flow, proper grid distribution can play a crucial part in the prediction and resolution of the solutions. The use of presently available numerical schemes to solve this problem on a uniform course grid system fails to resolve the characteristic flow features and hence do a poor job in simulating this flow. In this paper, an adaptive grid which adjusts itself as the solution evolves is used for a better resolution of the flow properties. Rai and Anderson's12 method is used to determine the grid speed; however, a different partial differential equation based on the conservative principle of grid arc lengths for clustering grids in one-dimensional flow is used along with the St. Venant equations to numerically simulate the flow. Both the subcritical and the supercritical flows under extreme boundary conditions are solved using this technique. With a specified number of grid points, this provides better quality solutions as compared to those obtained with uniformly distributed grids. 相似文献
18.
Regionality and seasonality of submesoscale and mesoscale turbulence in the North Pacific Ocean 总被引：1，自引：0，他引：1
The kinetic energy (KE) seasonality has been revealed by satellite altimeters in many oceanic regions. Question about the mechanisms that trigger this seasonality is still challenging. We address this question through the comparison of two numerical simulations. The first one, with a 1/10° horizontal grid spacing, 54 vertical levels, represents dynamics of physical scales larger than 50 km. The second one, with a 1/30° grid spacing, 100 vertical levels, takes into account the dynamics of physical scales down to 16 km. Comparison clearly emphasizes in the whole North Pacific Ocean, not only a significant KE increase by a factor up to three, but also the emergence of seasonal variability when the scale range 16–50 km (called submesoscales in this study) is taken into account. But the mechanisms explaining these KE changes display strong regional contrasts. In high KE regions, such the Kuroshio Extension and the western and eastern subtropics, frontal mixed-layer instabilities appear to be the main mechanism for the emergence of submesoscales in winter. Subsequent inverse kinetic energy cascade leads to the KE seasonality of larger scales. In other regions, in particular in subarctic regions, results suggest that the KE seasonality is principally produced by larger-scale instabilities with typical scales of 100 km and not so much by smaller-scale mixed-layer instabilities. Using arguments from geostrophic turbulence, the submesoscale impact in these regions is assumed to strengthen mesoscale eddies that become more coherent and not quickly dissipated, leading to a KE increase. 相似文献
19.
This work presents a procedure for developing a high-resolution, regional climatology estimate, named RClimo, off the coast of central California. This high-resolution climatology may provide an alternative way to initialize numerical nowcast/forecast exercises in coastal regions. The methodology includes two primary steps: (1) averaging available data on a high-resolution grid and (2) objective interpolating the resulting average profiles onto a regular grid. The first step involves the computation of averages over density layers in the vertical and allowing for data gaps in the horizontal if data are unavailable at a high resolution. The OA in the second step uses anisotropic correlation length scales derived from the data themselves and an averaging radius to preserve the scales and variability of the synoptic fields. 相似文献
20.
The flow of dense water along continental slopes is considered. There is a large literature on the topic based on observations and laboratory experiments. In addition, there are many analytical and numerical studies of dense water flows. In particular, there is a sequence of numerical investigations using the dynamics of overflow mixing and entrainment (DOME) setup. In these papers, the sensitivity of the solutions to numerical parameters such as grid size and numerical viscosity coefficients and to the choices of methods and models is investigated. In earlier DOME studies, three different bottom boundary conditions and a range of vertical grid sizes are applied. In other parts of the literature on numerical studies of oceanic gravity currents, there are statements that appear to contradict choices made on bottom boundary conditions in some of the DOME papers. In the present study, we therefore address the effects of the bottom boundary condition and vertical resolution in numerical investigations of dense water cascading on a slope. The main finding of the present paper is that it is feasible to capture the bottom Ekman layer dynamics adequately and cost efficiently by using a terrain-following model system using a quadratic drag law with a drag coefficient computed to give near-bottom velocity profiles in agreement with the logarithmic law of the wall. Many studies of dense water flows are performed with a quadratic bottom drag law and a constant drag coefficient. It is shown that when using this bottom boundary condition, Ekman drainage will not be adequately represented. In other studies of gravity flow, a no-slip bottom boundary condition is applied. With no-slip and a very fine resolution near the seabed, the solutions are essentially equal to the solutions obtained with a quadratic drag law and a drag coefficient computed to produce velocity profiles matching the logarithmic law of the wall. However, with coarser resolution near the seabed, there may be a substantial artificial blocking effect when using no-slip. 相似文献