Vertical Gridding Research Articles

Simulating multi-scale oceanic processes around Taiwan on unstructured grids

We validate a 3D unstructured-grid (UG) model for simulating multi-scale processes as occurred in Northwestern Pacific around Taiwan using recently developed new techniques (Zhang et al., Ocean Modeling, 102, 64–81, 2016) that require no bathymetry smoothing even for this region with prevalent steep bottom slopes and many islands. The focus is on short-term forecast for several months instead of long-term variability. Compared with satellite products, the errors for the simulated Sea-surface Height (SSH) and Sea-surface Temperature (SST) are similar to a reference data-assimilated global model. In the nearshore region, comparison with 34 tide gauges located around Taiwan indicates an average RMSE of 13 cm for the tidal elevation. The average RMSE for SST at 6 coastal buoys is 1.2 °C. The mean transport and eddy kinetic energy compare reasonably with previously published values and the reference model used to provide boundary and initial conditions. The model suggests ∼2-day interruption of Kuroshio east of Taiwan during a typhoon period. The effect of tidal mixing is shown to be significant nearshore. The multi-scale model is easily extendable to target regions of interest due to its UG framework and a flexible vertical gridding system, which is shown to be superior to terrain-following coordinates.

3-D Numerical Study of Mechanical Behaviors of Pile-Anchor System

Mechanical behaviors of pile-soil effect and anchor-soil effect are significantly important in supporting engineering activities of foundation pit. In this paper, finite difference method (FDM) was utilized to perform the numerical simulation of pile-anchor system, composed of supporting piles and pre-stressed anchor cables. Numerical simulations were on the basis of the foundation pit of Jinan’s West Railway Station, and 3D simulation analysis of foundation pit has been prepared during the whole processes of excavation, supporting and construction. The paper also analyzed the changes of bending moments of piles and axial forces of cables, and discussed mechanical behaviors of pile-anchor system, through comparisons with field monitoring. The results show that the parameters concluding vertical gridding’s number, cohesion of pile and soil, and pile stiffness have robust influences on supporting elements’ behaviors. Mechanical behaviors of supporting pile and axial forces of anchor cable changed dramatically, indicating that the potential failure form was converted from toppling failure to sliding failure.

Validating the turbulence parameterization schemes of a numerical model using eddy dissipation rate and turbulent kinetic energy measurements in terrain-disrupted airflow

A number of turbulence parameterization schemes are available in the latest version (6.0) of the Regional Atmospheric Modelling System (RAMS). Chan in Meteorol Atmos Phys 103:145–157, (2009), studied the performance of these schemes by simulating the eddy dissipation rate (EDR) distribution in the vicinity of the Hong Kong International Airport (HKIA) and comparing with the EDR measurements of remote-sensing instruments at the airport. For the e-l (turbulent kinetic energy − mixing length) scheme considered in that study, the asymptotic mixing length was assumed to be a constant. This assumption is changed in the present paper, a variable asymptotic mixing length is chosen and simulations of EDR fields are repeated for terrain-disrupted airflow in the vicinity of HKIA. It is found that, with a variable asymptotic mixing length, the performance of the e-l scheme is greatly improved. With suitable choice of the empirical constants in the turbulence closure, the accuracy of the EDR profile (in comparison with LIDAR and wind profiler measurements) is found to be comparable with that predicted by the Deardorff scheme. A study on the sensitivity of the simulation results to these empirical constants has also been performed. Moreover, as a follow-up of the previous study of Chan in Meteorol Atmos Phys 103:145–157, (2009), case studies have been conducted on the following issues of the model simulation of turbulence for aviation application: (a) the effect of vertical gridding on the simulation results, (b) possibility of false alarm (such as over-forecasting of EDR value) in light turbulence cases, and (c) the performance in the simulation of other turbulence intensity metric for aviation purpose, e.g. TKE.

Characterization of atmospheric optical turbulence for laser propagation

AbstractThe performance of imaging or transmission systems is limited by atmospheric inhomogeneities and fluctuations. Temperature and pressure variations associated with turbulent eddies cause random fluctuations of the refractive index structure constant (Cn2), which distort optical waves. The distortions may lead to significant blurring, scintillations, broadening and wander of the laser beam. High resolution nested numerical simulations are used to predict and characterize Optical Turbulence (OT) induced by jet streams and gravity waves under various local atmospheric conditions. Non‐homogeneous, anisotropic, non‐Kolmogorov patchy shear‐stratified turbulence in the Upper Troposphere and Lower Stratosphere (UTLS) requires that a fine mesh be used to resolve stiff velocity and temperature gradient profiles. Computational approach is based on nesting and refined vertical gridding in nested mesoscale/microscale codes. This methodology is applied to the analysis of field data from observations. Local distributions of simulated optical turbulence in the UTLS are obtained using very high resolution explicit simulations and parametrization formula. They show strongly laminated structures with thin layers of high values of the refractive index. These optical layers are characterized by steep vertical gradients of potential temperature and are located at the edges of relatively well mixed regions produced by shear instabilities and wave breaking.

Non-uniform adaptive vertical grids for 3D numerical ocean models

A new strategy for the vertical gridding in terrain-following 3D ocean models is presented here. The vertical grid adaptivity is partially given by a vertical diffusion equation for the vertical layer positions, with diffusivities being proportional to shear, stratification and distance from the boundaries. In the horizontal, the grid can be smoothed with respect to z-levels, grid layer slope and density. Lagrangian tendency of the grid movement is supported. The adaptive terrain-following grid can be set to be an Eulerian–Lagrangian grid, a hybrid σ – ρ or σ – z grid and combinations of these with great flexibility. With this, internal flow structures such as thermoclines can be well resolved and followed by the grid. A set of idealised examples is presented in the paper, which show that the introduced adaptive grid strategy reduces pressure gradient errors and numerical mixing significantly. The grid adaption strategy is easy to implement in various types of terrain-following ocean models. The idealised examples give evidence that the adaptive grids can improve realistic, long-term simulations of stratified seas while keeping the advantages of terrain-following coordinates.

Falcon: A Production Quality Distributed Memory Reservoir Simulator

Summary We describe a new production model, Falcon, that has achieved speeds on parallel computers that are 100 times faster on real world problems than current production models on a vector computer. Falcon has been used to conduct the largest, geostatistical reservoir study ever conducted within Amoco. In this paper we discuss the following: Falcon's data parallel paradigm with FORTRAN 90 and high performance FORTRAN (HPF); its single program, multiple data (SPMD) paradigm with message passing; efficient memory management that enables simulation of enormous studies; a numerical formulation that reconciles the generalized compositional approach (based on component masses and pressure) with earlier approaches (based on pressures and saturations), in a more general and more efficient approach. We also discuss Falcon's scalability up to 512 processor nodes and performance (timings and memory) achieved on a number of parallel platforms, including Cray Research's T3D and T3E, SGI's Power Challenge and Origin 2000, Thinking Machines' CM5, and IBM's SP2. Falcon also runs on single processor computers such as PC's and IBM's RS6000. We discuss a new parallel linear solver technology based on a fully parallel scalable implementation of incomplete lower-upper (ILU) preconditioning coupled with a GMRES or Orthomin iteration process. This naturally ordered global ILU preconditioner is scalable to hundreds of processors, efficiently solving the matrix problems arising from large scale simulations The use of the techniques described in this paper has enabled us to run problem sizes of up to 16.5 million gridblocks. Falcon was used to simulate fifty geostatistically derived realizations of a large, black oil waterflood system. The realizations, each with 2.3 million cells and 1,039 wells, took an average of 4.2 hours to execute on a 128-node CM5 computer, thus enabling the simulation study to finish in less than a month. In this field study, we bypassed upscaling through the use of fine vertical resolution gridding. Our focus has been on the applicability of Falcon to real world problems. Falcon can be used for modeling both small and very large reservoirs, including reservoirs characterized by geostatistics. It can be used to simulate black oil, gas/water, and dry gas reservoirs. And, a fully compositional feature is being developed.