Fast iterative solution of large undrained soil‐structure interaction problems

In view of rapid developments in iterative solvers, it is timely to re‐examine the merits of using mixed formulation for incompressible problems. This paper presents extensive numerical studies to compare the accuracy of undrained solutions resulting from the standard displacement formulation with a penalty term and the two‐field mixed formulation. The standard displacement and two‐field mixed formulations are solved using both direct and iterative approaches to assess if it is cost‐effective to achieve more accurate solutions. Numerical studies of a simple footing problem show that the mixed formulation is able to solve the incompressible problem ‘exactly’, does not create pressure and stress instabilities, and obviate the need for an ad hoc penalty number. In addition, for large‐scale problems where it is not possible to perform direct solutions entirely within available random access memory, it turns out that the larger system of equations from mixed formulation also can be solved much more efficiently than the smaller system of equations arising from standard formulation by using the symmetric quasi‐minimal residual (SQMR) method with the generalized Jacobi (GJ) preconditioner. Iterative solution by SQMR with GJ preconditioning also is more elegant, faster, and more accurate than the popular Uzawa method. Copyright © 2003 John Wiley & Sons, Ltd.     

Screening of polymers on selected Hawaii soils for erosion reduction and particle settling

In recent years, high‐molecular‐weight anionic polyacrylamides (PAMs) have been tested on a variety of soils, primarily in temperate climates. However, little information is available regarding the effectiveness of PAM for preventing soil loss through runoff in tropical settings. Screening tests were performed using three negatively charged PAMs and one positively charged PAM on five Hawaii soils (two Oxisols, one Vertisol, and two Aridisols) to determine erosion loss, sediment settling, and aggregate stability. A laboratory‐scale rainfall simulator was used to apply erosive rainfall at intensities from 5 to 8·5 cm h^?1 at various PAM doses applied in both dry and solution forms. Soil detachment due to splash and runoff, as well as the runoff and percolate water volumes, were measured for initial and successive storms. The impact of PAM on particle settling and aggregate stability was also evaluated for selected soil‐treatment combinations. Among the PAMs, Superfloc A‐836 was most effective, and significantly reduced runoff and splash sediment loss for the Wahiawa Oxisol and Pakini Andisol at rates varying between 10 and 50 kg ha^?1. Reduced runoff and splash sediment loss were also noted for PAM Aerotil‐D when applied in solution form to the Wahiawa Oxisol. Significant reductions in soil loss were not noted for either the Lualualei Vertisol or the Holomua Oxisol. It is believed that the high montmorillonite content of the Lualualei Vertisol and the low cation‐exchange capacity of the Holomua Oxisol diminished the effectiveness of the various PAMs tested. The polymers were also found to enhance sediment settling of all soils and helped improve their aggregate stability. This screening study shows the potential use of PAM for tropical soils for applications such as infiltration enhancement, runoff reduction, and enhanced sedimentation of detention ponds. Copyright © 2005 John Wiley & Sons, Ltd.     

Analytical solutions of one‐dimensional infiltration before and after ponding

This paper is a continuation of previous research, which obtained a convenient solution for arbitrary surface fluxes before ponding. By means of Fourier Transformation this has been extended to derive analytical solutions of a linearized Richards' equation for arbitrary input fluxes after surface saturation. Copyright © 2003 John Wiley & Sons, Ltd.     

Magnitude-frequency relationship of coastal sand delivery by a southern California stream

Coastal sand delivery by a stream in southern California is estimated based on a numerical model which stimulates unsteady flow, sediment transport, and the associated channel adjustments for a stream-delta system. An average annual sediment yield of 51,400 m³/yr is estimated for the San Dieguito River, which drains a semiarid watershed controlled by dams. Of the total sand delivery by the stream, 20.5 percent is contributed from floods greater than the 100-year flood; 17.6 percent from those between the 50- and 100-year events; 28.4 percent from those between the 25- and 50 year floods; and 33.5 percent from those smaller than the 25-year flood.     

Sedimentology of a subtidal,tide-dominated sand body in the Yellow Sea,southwest Korea

Odanam Satoe, a subtidal, tide-dominated sand body in the Yellow Sea, Korea, is linear in plan and asymmetrical in cross-section. It consists of fine- to medium-grained, well-sorted subangular sand. Bedforms consist of high-amplitude (1–2 m) sandwaves on the lower flanks of the gentler-sloping bar surface, and medium-amplitude (0.5-1 m) sandwaves on the sand body trough adjoining the steeper face, the bar crest and shallower parts of the gently sloping bar surface. Bedforms are absent on the relatively steeper bar surface, which is characterized by 2° slopes. Bedform orientation on the gentler slope is oblique by 30° to the bar crest, parallel to the sand-body crest on the crest itself, and opposite to the steeper sand-body face in the trough below the steeper slope of the bar.Bottom current velocity data show that tidal currents are semi-rotary with a flood time—velocity asymmetry over the gentler slope, and ebb time—velocity asymmetry over the steeper slope during most of the tidal cycle. Tidal-current flow parallels bar elongation over the steeper slope, whereas over the gentler slope, tidal-current flow is directed at 30° to the bar crest and changes to normal to the crest one hour prior to low tide. Bedform orientation mapped with side-scan sonar shows agreement with these flow directions.Sand dispersal around the sand body is controlled by time—velocity asymmetry and partial rotary flow directions of tidal currents. This circulation causes not only a trapezoidal mode of grain dispersal, but also westerly migration of the sand body documented from comparative bathymetric surveys in 1964 and 1980.     

1976年夏季东海陆架区浮游植物生态的研究

关于东海浮游植物,虽然国内、外都有人进行过一些调査研究,取得一些成果,但因调查时间较短或调查范围較小所获资料不全,都未能代表东海浮游植物的全貌。中国科学院海洋研究所继全国海洋综合调査之后,自1976年夏季又在东海陆架区进行了综合调查。这次调查，分析了东海浮游植物的个体数量和种数的分布状况,讨论了其主要种类的生态性质、群落组成和它们与水文环境的关系,对东海浮游植物的生态状况有了进一步了解。 1976年夏季在东海陆架区进行了两次(6月25日至7月23日和8月27日至9月21日)调查,调查范围南到北纬26°30′,东到东经127°,在沿纬度安排的6个断面共44个测站上(第一航次还在黄海南部增加一斜断面),分別进行了水文观测,并用网目为64平方微米的小型浮游生物网(Juday phytoplankton net),从水底到表面垂直拖网进行浮游植物采样,以个体计数法计数分析。     

A general method for calculating three-dimensional laminar and turbulent boundary layers on ship hulls

A general method for representing the flow properties in the three-dimensional boundary layers around ship hulls of arbitrary shape is described. It makes use of an efficient two-point finite-diffirence schem to solve the boundary-layer equations and includes an algebraic eddy-viscosity representaion of the Reynolds-stress ternsor. The numericzal method contains novel and desirable features and allows the calculation of flows in which the circumferential velocity component contains regions of flow reversal across the boundary layer. The inviscid pressure distribution is determined with the Douglas-Neumann method which, if necessary, can conveniently allow for the boundary-layer displacement surface. To allow its application to ships, and particularly to those with double-elliptic and flat-bottomed hulls, a non-orthogonal coordinate system has been developed and is shown to be economical, precise and comparatively easy to use. Present calculations relate to zero Froude number but they can be extended to include the effects of a water wave and local regions of flow separation which may stem from bulbous-bow geometries.     

Fully nonlinear numerical wave tank (NWT) simulations and wave run-up prediction around 3-D structures

A finite-difference scheme and a modified marker-and-cell (MAC) algorithm have been developed to investigate the interactions of fully nonlinear waves with two- or three-dimensional structures of arbitrary shape. The Navier–Stokes (NS) and continuity equations are solved in the computational domain and the boundary values are updated at each time step by the finite-difference time-marching scheme in the framework of a rectangular coordinate system. The fully nonlinear kinematic free-surface condition is implemented by the marker-density function (MDF) technique developed for two fluid layers.To demonstrate the capability and accuracy of the present method, the numerical simulation of backstep flows with free-surface, and the numerical tests of the MDF technique with limit functions are conducted. The 3D program was then applied to nonlinear wave interactions with conical gravity platforms of circular and octagonal cross-sections. The numerical prediction of maximum wave run-up on arctic structures is compared with the prediction of the Shore Protection Manual (SPM) method and those of linear and second-order diffraction analyses based on potential theory and boundary element method (BEM). Through this comparison, the effects of non-linearity and viscosity on wave loading and run-up are discussed.     

A higher-order σ-coordinate non-hydrostatic model for nonlinear surface waves

A higher-order non-hydrostatic model in a σ-coordinate system is developed. The model uses an implicit finite difference scheme on a staggered grid to simultaneously solve the unsteady Navier-Stokes equations (NSE) with the free-surface boundary conditions. An integral method is applied to resolve the top-layer non-hydrostatic pressure, allowing for accurately resolving free-surface wave propagation. In contrast to the previous work, a higher-order spatial discretization is utilized to approximate the large horizontal pressure gradient due to steep surface waves or rapidly varying topographies. An efficient direct solver is developed to solve the resulting block hepta-diagonal matrix system. Accuracy of the new model is validated by linear and nonlinear standing waves and progressive waves. The model is then used to examine freak (extreme) waves. Features of downshifting focusing location and wave asymmetry characteristics are predicted on the temporal and spatial domains of a freak wave.