A comparison between experimental and numerical investigations of the motion of the water surface in a model surge tank

This paper describes an investigation into the motion of the water surface in a simple model surge tank, and the relevant factors governing its behaviour. The oscillation of the free water surface is an unsteady flow problem, which is amenable to a mathematical solution using a finite difference step-by-step integration procedure. For comparison, two such methods are presented: (i) a simple initial value method and (ii) a predictor-corrector technique. Computer programs have been developed linked to a graph plotter to give a visual presentation of the numerical solutions together with the experimental results of the damped oscillation of the water surface in the model surge tank.     

Coupled wave and surge modelling for the eastern Irish Sea and implications for model wind-stress

We revisit the surge of November 1977, a storm event which caused damage on the Sefton coast in NW England. A hindcast has been made with a coupled surge-tide-wave model, to investigate whether a wave-dependent surface drag is necessary for accurate surge prediction, and also if this can be represented by an optimised Charnock parameter. The Proudman Oceanographic Laboratory Coastal Modelling System-Wave Model (POLCOMS-WAM) has been used to model combined tides, surges, waves and wave-current interaction in the Irish Sea on a 1.85 km grid. This period has been previously thoroughly studied, e.g. Jones and Davies [Jones, J.E., Davies, A.M., 1998. Storm surge computations for the Irish Sea using a three-dimensional numerical model including wave-current interaction. Continental Shelf Research 18(2), 201–251] and we build upon this previous work to validate the POLCOMS-WAM model to test the accuracy of surge elevation predictions in the study area. A one-way nested approach has been set up from larger scale models to the Irish Sea model. It was demonstrated that (as expected) swell from the North Atlantic does not have a significant impact in the eastern Irish Sea. To capture the external surge generated outside of the Irish Sea a (1/9° by 1/6°) model extending beyond the continental shelf edge was run using the POLCOMS model for tide and surge.     

Testing a new sampler for measuring plot soil loss

In order to measure soil loss in equipped plots the estimate of the weight of solid material intercepted at their lower end is required. At the experimental area of Sparacia, Sicily, the runoff produced by an erosive event is collected within storage tanks with a capacity of about 1 m³. In this paper, the use of a new sampler is proposed to measure easily the weight of solid material eroded from an experimental plot and collected into a storage tank. The sampler is a cylinder having a closing valve at the bottom. Two different series of runs were carried out both to test the reliability of the sampler and to establish a sampling procedure, respectively. An analysis of various sampling configurations usable in the field differentiated by the number and location of sampling verticals in the tank cross‐section was finally carried out. The results of the present investigation are that the concentration measurement by the sampler was more accurate than that obtained by other methods involving a collection tank, agitation and sampling of the suspension. This sampler is cheap and usable in combination with a quick field sampling procedure which is particularly advisable when the number of plots equipped at an experimental area is large. The sampler was tested using a clay soil contained within cylinders and a cubic tank, but it appeared also to be usable with coarser sediment than clay and in combination with tanks having a different shape. Copyright © 2015 John Wiley & Sons, Ltd.     

The effect of seismic uplift on the shell stresses of liquid‐storage tanks

Previous theoretical studies have shown that tank uplift, that is, separation of the tank base from the foundation, generally reduces the base shear and the base moment. However, there is a paucity of experimental investigations concerning the effect of uplift on the tank wall stresses, which is the principal parameter that controls the seismic design of liquid‐storage tanks. This paper reports a series of shake table experiments on a polyvinyl chloride model tank containing water. A comparison of the seismic behaviour of the tank with and without anchorage is described. Stochastically generated ground motions, based on a Japanese design spectrum, and three tank aspect ratios (liquid‐height/radius) are considered. Measurements were made of the stresses at the outer shell of the tank, the tank wall acceleration and the horizontal displacement at the top of the tank. While the top displacement and the tank shell acceleration increased when uplift was allowed, axial compressive stresses decreased by between 35% and 64% with tank uplift. The effect of uplift on the hoop stresses was variable depending on the aspect ratio. A comparison of experimental values with a numerical model is provided. Copyright © 2015 John Wiley & Sons, Ltd.     

Application of an unstructured mesh model to storm surge propagation in the Mersey estuary region of the Irish Sea

The storm surge period of 13–16 November 1977 when there was a major positive surge followed by a negative surge in the Irish Sea is investigated using a two-dimensional unstructured mesh model of the west coast of Britain. The model accounts for tidal and external surge forcing across its open boundaries which are situated in the Celtic Sea and off the west coast of Scotland. Although this period has been examined previously using a uniform finite-difference model, and a finite element model, neither of these could resolve the Mersey estuary which is the focus of the present study. By using a finite element model with very high mesh resolution within the Mersey, the spatial variability of surge elevations and currents within the Mersey to rapidly changing surge dynamics can be examined. The mesh in the model varies from about 7 km in deep water, to the order of 100 m in the Mersey, with the largest mesh length reaching 17 km in deep offshore regions, and smallest of order 26 m occurring in shallow coastal regions of the Mersey estuary. The model accounts for wetting/drying which occurs in shallow water coastal areas. Calculations showed that during the positive surge period, the amplitude and speed of propagation of the surge was largest in the deep water channels. This gave rise to significant spatial variability of surge elevations and currents within the estuary. As wind stresses decreased over the Irish Sea, a negative surge occurred over Liverpool Bay and at the entrance to the Mersey. However, within the Mersey there was a local positive surge which continued to propagate down the estuary. This clearly showed that although the large scale response of the Irish Sea to changing wind fields occurred rapidly, the response in the Mersey was much slower. These calculations with a west coast variable mesh model that included a high-resolution representation of the Mersey revealed for the first time how elevations and currents within the Mersey responded to Irish Sea surges that rapidly changed from positive to negative.     

Storm surge computations for the west coast of Britain using a finite element model (TELEMAC)

An unstructured mesh finite element model of the sea region off the west coast of Britain is used to examine the storm surge event of November 1977. This period is chosen because accurate meteorological data to drive the model and coastal observations for validation purposes are available. In addition, previous published results from a coarse-grid (resolution 7 km) finite difference model of the region and high-resolution (1 km) limited area (namely eastern Irish Sea) model are available for comparison purposes. To enable a “like with like” comparison to be made, the finite element model covers the same domain and has the same meteorological forcing as these earlier finite difference models. In addition, the mesh is based on an identical set of water depths. Calculations show that the finite element model can reproduce both the “external” and “internal” components of the surge in the region. This shows that the “far field” (external) component of the surge can accurately propagate through the irregular mesh, and the model responds accurately, without over- or under-damping, to local wind forcing. Calculations show significant temporal and spatial variability in the surge in close agreement with that found in earlier finite difference calculations. In addition, root mean square errors between computed and observed surge are comparable to those found in previous finite different calculations. The ability to vary the mesh in nearshore regions reveals appreciable small-scale variability that was not found in the previous finite difference solutions. However, the requirement to perform a “like with like” comparison using the same water depths means that the full potential of the unstructured grid model to improve resolution in the nearshore region is inhibited. This is clearly evident in the Mersey estuary region where a higher resolution unstructured mesh model, forced with uniform winds, had shown high topographic variability due to small-scale variations in topography that are not resolved here. Despite the lack of high resolution in the nearshore region, the model showed results that were consistent with the previous storm surge models of the region. Calculations suggest that to improve on these earlier results, a finer nearshore mesh is required based upon accurate nearshore topography.     

液压阻尼系统控制浮放立式储罐提离失稳

立式圆柱储罐在油田、石油化工等企业用途广泛，多用作存贮易燃、易爆介质，一旦遭遇震害，后果将十分严重。本文针对浮放立式储罐地震响应的提离失稳问题，采用液压阻尼系统（HDS）用以控制储罐的提离反应，建立了安装HDS储罐的力学分析模型并进行了数值计算，对安装HDS后储罐的提离反应进行了分析。仿真计算结果表明：HDS可以有效地减小提离应力，从而地震烈度可以降低1度进行储罐的提离设计。     

Simulation of the 1953 storm surge in the North Sea

The 1953 North Sea floods, the Big Flood, was one of the worst natural disasters in Europe in modern times and is probably one of the most studied severe coastal floods. Several factors led to the devastating storm surge along the southern North Sea coast in combination of strong and sustained northerly winds, invert barometric effect, high spring tide, and an accumulation of the large surge in the Strait of Dover. However, the storm waves and their roles during the 1953 North Sea storm surge are not well investigated. Therefore, the effect of wave setup due to breaking waves in the storm surge processes is investigated through numerical experiments. A coupled process-based tide-wave-surge model was used to investigate and simulate the storm surge in the North Sea during January 31–February 1, 1953 and validated by comparing with historical water level records at tide gauges and wave observations at light vessels in the North Sea. Meteorological forcing inputs for the period, January 27–February 3, 1953 are reproduced from ERA-20C reanalysis data with a constant correction factor for winds. From the simulation results, it is found that, in addition to the high water due to wind setup, wave setup due to breaking waves nearshore play a role of approximately 10% of the storm surge peaks with approximately 0.2 m. The resulting modeling system can be used extensively for the preparedness of the storm surge and wave of extreme condition, and usual barotropic forecast.     

Design and application of a direct-push vadose zone gravel permeameter

A borehole permeameter is well suited for testing saturated hydraulic conductivity (K(sat)) at specific depths in the vadose zone. Most applications of the method involve fine-grained soils that allow hand auguring of test holes and require a small water reservoir to maintain a constant head. In non-cohesive gravels, hand-dug test holes are difficult to excavate, holes are prone to collapse, and large volumes of water are necessary to maintain a constant head for the duration of the test. For coarse alluvial gravels, a direct-push steel permeameter was designed to place a slotted pipe at a specific sampling depth. Measurements can be made at successive depths at the same location. A 3790 L (1000 gallons) trailer-mounted water tank maintained a constant head in the permeameter. Head in the portable tank was measured with a pressure transducer and flow was calculated based on a volumetric rating curve. A U.S. Bureau of Reclamation analytical method was utilized to calculate K(sat). Measurements with the permeameter at a field site were similar to those reported from falling-head tests.     

Soil-structure interaction effects for laterally excited liquid storage tanks

A comprehensive study is made of the effects of soil-structure interaction on the response of liquid containing, upright, circular cylindrical tanks subjected to a horizontal component of ground shaking. A simple, physically motivated method of analysis is employed which elucidates the effects and relative importance of the principal actions involved. Both the impulsive and convective actions of the liquid are examined. The interrelationship of the tank responses to horizontal and rocking actions of the foundation is established, and the well known mechanical model for laterally excited, rigid tanks supported on a non-deformable medium is generalized to permit consideration of the effects of tank and ground flexibilities and base rocking. Critical responses are evaluated for harmonic and seismic excitations over wide ranges of tank proportions and soil stiffnesses, and the results are presented in a form convenient for use in practical applications. In addition to a precise method of analysis, an approximate, hand-computation method is presented with which the effects of the primary parameters may be evaluated readily. The soil-structure interaction effects in the latter approach are provided for by modifying the natural frequency and damping of the tank-liquid system and evaluating its response to the prescribed free-field ground motion considering the tank to be rigidly supported at the base. The requisite modifications may be determined from information presented herein. It is shown that soil-structure interaction may reduce significantly the impulsive components of response but that it has a negligible effect on the convective components.     

Interrelations among pyroclastic surge,pyroclastic flow,and lahars in Smith Creek valley during first minutes of 18 May 1980 eruption of Mount St. Helens,USA

A devastating pyroclastic surge and resultant lahars at Mount St. Helens on 18 May 1980 produced several catastrophic flowages into tributaries on the northeast volcano flank. The tributaries channeled the flows to Smith Creek valley, which lies within the area devastated by the surge but was unaffected by the great debris avalanche on the north flank. Stratigraphy shows that the pyroclastic surge preceded the lahars; there is no notable “wet” character to the surge deposits. Therefore the lahars must have originated as snowmelt, not as ejected water-saturated debris that segregated from the pyroclastic surge as has been inferred for other flanks of the volcano. In stratigraphic order the Smith Creek valley-floor materials comprise (1) a complex valley-bottom facies of the pyroclastic surge and a related pyroclastic flow, (2) an unusual hummocky diamict caused by complex mixing of lahars with the dry pyroclastic debris, and (3) deposits of secondary pyroclastic flows. These units are capped by silt containing accretionary lapilli, which began falling from a rapidly expanding mushroom-shaped cloud 20 minutes after the eruption's onset. The Smith Creek valley-bottom pyroclastic facies consists of (a) a weakly graded basal bed of fines-poor granular sand, the deposit of a low-concentration lithic pyroclastic surge, and (b) a bed of very poorly sorted pebble to cobble gravel inversely graded near its base, the deposit of a high-concentration lithic pyroclastic flow. The surge apparently segregated while crossing the steep headwater tributaries of Smith Creek; large fragments that settled from the turbulent surge formed a dense pyroclastic flow along the valley floor that lagged behind the front of the overland surge. The unusual hummocky diamict as thick as 15 m contains large lithic clasts supported by a tough, brown muddy sand matrix like that of lahar deposits upvalley. This unit contains irregular friable lenses and pods meters in diameter, blocks incorporated from the underlying dry and hot pyroclastic material that had been deposited only moments earlier. The hummocky unit is the deposit of a high-viscosity debris flow which formed when lahars mingled with the pyroclastic materials on Smith Creek valley floor. Overlying the debris flow are voluminous pyroclastic deposits of pebbly sand cut by fines-poor gas-escape pipes and containing charred wood. The deposits are thickest in topographic lows along margins of the hummocky diamict. Emplaced several minutes after the hot surge had passed, this is the deposit of numerous secondary pyroclastic flows derived from surge material deposited unstably on steep valley sides.     

A 3D ERT study of solute transport in a large experimental tank

A high resolution, cross-borehole, 3D electrical resistivity tomography (ERT) study of solute transport was conducted in a large experimental tank. ERT voxels comprising the time sequence of electrical images were converted into a 3D array of ERT estimated fluid conductivity breakthrough curves and compared with direct measurements of fluid conductivity breakthrough made in wells. The 3D ERT images of solute transport behaviour were also compared with predictions based on a 3D finite-element, coupled flow and transport model, accounting for gravity induced flow caused by concentration differences.The tank (dimensions 185×245×186 cm) was filled with medium sand, with a gravel channel and a fine sand layer installed. This heterogeneous system was designed to complicate solute transport behaviour relative to a homogeneous sand tank, and to thus provide a challenging but insightful analysis of the ability of 3D ERT to resolve transport phenomena. Four ERT arrays and 20 piezometers were installed during filling. A NaCl tracer (conductivity 1.34 S/m) was injected and intensively monitored with 3D ERT and direct sampling of fluid chemistry in piezometers.We converted the bulk conductivity estimate for 250 voxels in the ERT imaged volume into ERT estimated voxel fluid conductivity by assuming that matrix conduction in the tank is negligible. In general, the ERT voxel response is in reasonable agreement with the shape of fluid conductivity breakthrough observed in six wells in which direct measurements of fluid conductivity were made. However, discrepancies occur, particularly at early times, which we attribute to differences between the scale of the image voxels and the fluid conductivity measurement, measurement errors mapped into the electrical inversion and artificial image roughness resulting from the inversion.ERT images revealed the 3D tracer distribution at 15 times after tracer injection. The general pattern and timing of solute breakthrough observed with ERT agreed with that predicted from the flow/transport modelling. However, the ERT images indicate a vertical component of tracer transport and preferential flow paths in the medium sand. We attribute this to transient vertical gradients established during tracer injection, and heterogeneity caused by sorting of the sand resulting from the filling procedure. In this study, ERT provided a unique dataset of 250 voxel breakthrough curves in 1.04 m³. The use of 3D ERT to generate an array of densely sampled estimated fluid conductivity breakthrough curves is a potentially powerful tool for quantifying solute transport processes.     

Improved water-level forecasting for the Northwest European Shelf and North Sea through direct modelling of tide, surge and non-linear interaction

In real-time operational coastal forecasting systems for the northwest European shelf, the representation accuracy of tide–surge models commonly suffers from insufficiently accurate tidal representation, especially in shallow near-shore areas with complex bathymetry and geometry. Therefore, in conventional operational systems, the surge component from numerical model simulations is used, while the harmonically predicted tide, accurately known from harmonic analysis of tide gauge measurements, is added to forecast the full water-level signal at tide gauge locations. Although there are errors associated with this so-called astronomical correction (e.g. because of the assumption of linearity of tide and surge), for current operational models, astronomical correction has nevertheless been shown to increase the representation accuracy of the full water-level signal. The simulated modulation of the surge through non-linear tide–surge interaction is affected by the poor representation of the tide signal in the tide–surge model, which astronomical correction does not improve. Furthermore, astronomical correction can only be applied to locations where the astronomic tide is known through a harmonic analysis of in situ measurements at tide gauge stations. This provides a strong motivation to improve both tide and surge representation of numerical models used in forecasting. In the present paper, we propose a new generation tide–surge model for the northwest European Shelf (DCSMv6). This is the first application on this scale in which the tidal representation is such that astronomical correction no longer improves the accuracy of the total water-level representation and where, consequently, the straightforward direct model forecasting of total water levels is better. The methodology applied to improve both tide and surge representation of the model is discussed, with emphasis on the use of satellite altimeter data and data assimilation techniques for reducing parameter uncertainty. Historic DCSMv6 model simulations are compared against shelf wide observations for a full calendar year. For a selection of stations, these results are compared to those with astronomical correction, which confirms that the tide representation in coastal regions has sufficient accuracy, and that forecasting total water levels directly yields superior results.     

Storm surge computations in estuarine and near-coastal regions: the Mersey estuary and Irish Sea area

An unstructured grid storm surge model of the west coast of Britain, incorporating a high-resolution representation of the Mersey estuary is used to examine storm surge dynamics in the region. The focus of the study is the major surge that occurred during the period 11–14 November 1977, which has been investigated previously using uniform grid finite difference models and a finite element model of the west coast of Britain. However, none of these models included the Mersey estuary. Comparison of solutions in the eastern Irish Sea with those computed with these earlier models showed that, away from the Liverpool Bay region, the inclusion of the Mersey estuary had little effect. However, at the entrance to the Mersey, its inclusion did influence the solution. By including a detailed representation of the Mersey estuary within the model, it was possible to conduct a detailed study of storm surge propagation in the Mersey, which had never previously been performed. This detailed study showed for the first time that the surge’s temporal variability within the estuary is influenced by surge elevation at its entrance. This varies with time as a function of spatial and temporal variations of wind stress over the west coast of Britain. Within the Mersey, calculations show that the spatial variability is mainly determined by changes in bottom topography, which had not been included in earlier finite difference models. However, since water depth is influenced by variations in tidal elevation, this, together with tide surge interaction through bottom friction and momentum advection, influences the surge. The ability of the finite element model to vary the mesh in near-shore regions to such an extent that it can resolve the Mersey and hence the impact of the Mersey estuary upon the Liverpool Bay circulation shows that it has distinct advantages over earlier finite difference models. In the absence of detailed measurements within the Mersey at the time of the surge, it was not possible to validate predicted surge elevations within the Mersey. However, significant insight into physical processes influencing the surge propagation down the estuary, its reflection and spatial/temporal variability could be gained.     

A simplified seismic analysis of rigid base liquid storage tanks under vertical exicitation with soil-structure interaction

A study is carried out to evaluate dynamic response of an elastic circular cylindrical tank having a rigid base under a vertical excitation taking into consideration the interaction with the foundation soil. At first, the soil is represented by frequency-independent parameters. Two coupled differential equations, governing the motion of the shell and the base, are solved using a step by step integration technique. The hydrodynamic pressures, acting on the shell and on the base, are derived from a velocity potential function which satisfies the Laplace equation and the appropriate boundary conditions. The response of the simplified model of a tank having a rigid base on a stiff foundation soil is compated to that obtained elsewhere to check the accuracy of the present model. Reasonable agreement is found between the maximum wall displacement and the associated stresses with those found by a more elaborate model. The interaction of the tank and the soil reduces the response than that calculated under the assumption of a rigid foundation soil. A parametric study to examine the effects of the height-to-radius ratio of the tank, and the effects of the shear wave velocity of the soil on the response is conducted. Varieties of foundation models are used to assess the sensitivity of the response to the variation in the soil parameters. Finally, a more representative solution for the problem in the frequency domain is obtained where the soil is appropriately modelled by frequency-dependent parameters. The transfer functions of the response of the tank wall and of the relative base motion are evaluated, and a comparison between the frequency-dependent and the frequency-independent solutions is made.     

The 1982 eruptions of El Chichon volcano,Mexico (3): Physical properties of pyroclastic surges

Two major pyroclastic surges generated during the 4 April 1982 eruption of El Chichon devastated an area of 153 km² with a quasi-radial distribution around the volcano. The hot surge clouds carbonized wood throughout their extent and were too hot to allow accretionary lapilli formation by vapor condensation. Field evidence indicates voidage fraction of 0.99 in the surge cloud with extensive entrainment of air. Thermal calculations indicate that heat content of pyroclasts can heat entrained air and maintain high temperatures in the surge cloud. The dominant bed form of the surge deposits are sand waves shaped in dune forms with vertical form index of 10–20, characterized by stoss-side erosion and lee-side deposition of 1–10 cm reversely graded laminae. A systematic decrease in maximum lithic diameter with distance from source is accompanied by decrease in wavelength and amplitude. Modal analysis indicates fractionation of glass and pumice from the surge cloud relative to crystals, resulting in loss of at least 10%–25% of the cloud mass due to winnowing out of fines during surge emplacement. Greatest fractionation from the –1.0–0.0– grain sizes reflects relatively lower pumice particle density in this range and segregation in the formative stages of the surge cloud. Extensive pumice rounding indicates abrasion during bed-load transport. Flow of pyroclastic debris in the turbulent surge cloud was by combination of bed-load and suspended-load transport. The surges are viewed as expanding pyroclastic gravity flows, which entrain and mix with air during transport. The balance between sedimentation at the base of the surge cloud and expansion due to entrainment of air contributed to low cloud density and internal turbulence, which persisted to the distal edge of the surge zone.     

Tank damage during the may 1983 coalinga earthquake

The magnitude M = 6-5 Coalinga earthquake of 2 May 1983 caused intense ground shaking throughout the epicentral region. Unanchored cylindrical ground supported tanks located at six sites within this oil producing area were damaged; damages included elephant's foot buckling at the base of three moderate sized tanks, joint rupture and top shell buckling in one large old rivetted tank, bottom plate rupture of a relatively new welded tank and damage to the floating roofs of 11 tanks. Also oil spilled over the top of many tanks and secondary damages occurred in pipe connections, ladders, etc. In this paper an estimate is made of the intensity of ground motion at each of the tank sites, based on strong motion records made during the main shock and the strongest aftershock. Then response parameters specified by current codes are correlated with the damages observed at each tank site. Based on this comparison, it is concluded that current U.S. practice under-estimates the sloshing response of tanks with floating roofs and does not adequately address the uplifting mechanism of unanchored ground supported tanks.     

Rapidly installed temporary gauging for hurricane waves and surge,and application to Hurricane Gustav

Hurricanes can produce extreme nearshore waves and surge, but permanent gauging stations are often much sparser than is desired. This paper describes the rationale behind and outline for rapidly installed temporary coastal gauges, and presents results during Hurricane Gustav (2008). Within 48 h prior to landfall, twenty self-recording pressure gauges were deployed in depths of 1.4–23 m over more than 700 km of coastline, using helicopters to cover the large distances. Results showed a complex picture that was strongly dependent on location. East of the Mississippi Delta, open coast waves were large, and surge reached 3.8 m NAVD88 in marshes. West of the delta but near landfall, waves and surge were generally smaller as the river levees blocked flow from East to West. West of landfall, both waves and surge were very small and the most prominent feature was a water level drawdown that reached 1.5 m. Wave spectra varied strongly depending both on location and time from landfall.     

Automatic calibration of the tank model / L'étalonnage automatique d'un modèle à cisterne

Abstract

The automatic calibration is done not by a hill-top climbing method but by a trial and error method carried out automatically by a computer program. The feedback procedure is made by comparing some criteria obtained from the observed hydrograph and the calculated hydrograph output from the working tank model. The two criteria are discharge volume and the shape of the hydrograph. The feedbacks of these two criteria correspond to dispacement feedback and velocity feedback in automatic control. The output of the working tank model is composed of components, the outputs from each of the tanks. Correspondingly, the whole period is divided into subperiods, in each of which each of the components plays the main part. The volume and shape are calculated in each subperiod and are used for the adjustment of the respective tanks. The feedback procedure starts from some initial model and converges very quickly after several (usually less than 15) iterations, and the result obtained is very good.