Three-dimensional liquid sloshing in a tank with baffles

A numerical model has been developed to study three-dimensional (3D) liquid sloshing in a tank with baffles. The numerical model solves the spatially averaged Navier-Stokes equations, which are constructed on a non-inertial reference frame having six degree-of-freedom (DOF) of motions. The large-eddy-simulation (LES) approach is employed to model turbulence by using the Smagorinsky sub-grid scale (SGS) closure model. The two-step projection method is employed in the numerical solutions, aided by the Bi-CGSTAB technique to solve the pressure Poisson equation for the filtered pressure field. The second-order accurate volume-of-fluid (VOF) method is used to track the distorted and broken free surface. The baffles in the tank are modeled by the concept of virtual boundary force (VBF) method. The numerical model is first validated against the available analytical solution and experimental data for two-dimensional (2D) liquid sloshing in a tank without baffles. The 2D liquid sloshing in tanks with baffles is then investigated. The numerical results are compared with other results from available literatures. Good agreement is obtained. Finally, the model is used to study 3D liquid sloshing in a tank with vertical baffles. The effect of the baffle is investigated and discussed.     

Tectonic evolution of the Cape and Karoo basins of South Africa

The Cape and Karoo basins formed within the continental interior of Gondwana. Subsidence resulted from the vertical motion of rigid basement blocks and intervening crustal faults. Each basin episode records a three-stage evolution consisting of crustal uplift, fault-controlled subsidence, and long periods of regional subsidence largely unaccompanied by faulting or erosional truncation. The large-scale episodes of subsidence were probably the result of lithospheric deflection due to subduction-driven mantle flow. The early Paleozoic Cape basin records the combined effects of a north-dipping intra-crustal décollement (a late Neoproterozoic suture) and a right-stepping offset between thick Rio de la Plata craton and Namaqua basement. Following the Saldanian orogeny, a suite of small rift basins and their post-rift drape formed at this releasing stepover. Great thicknesses of quartz sandstone (Ordovician–Silurian) and mudstone (Devonian) accumulation are attributed to subsidence by rheological weakening and mantle flow. In contrast, the Karoo basin is a cratonic cover that mimics the underlying basement blocks. The Permian Ecca and lower Beaufort groups were deposited in a southward-deepening ramp syncline by extensional decoupling on the intra-crustal décollement. Reflection seismic and deep-burial diagenetic studies indicate that the Cape orogeny started in the Early Triassic. Deformation was partitioned into basement-involved strike-slip faults and thin-skinned thrusting. Uplift of the Namaqua basement resulted in erosion of the Beaufort cover. East of the Cape fold belt, contemporaneous subsidence and tilting of the Natal basement created a late Karoo transtensional foreland basin, the Stormberg depocentre. Early Jurassic tectonic resetting and continental flood basalts terminated the Karoo basin.     

A predictive numerical model for potential mapping of the gas hydrate stability zone in the Gulf of Cadiz

This paper presents a computational model for mapping the regional 3D distribution in which seafloor gas hydrates would be stable, that is carried out in a Geographical Information System (GIS) environment. The construction of the model is comprised of three primary steps, namely: (1) the construction of surfaces for the various variables based on available 3D data (seafloor temperature, geothermal gradient and depth-pressure); (2) the calculation of the gas function equilibrium functions for the various hydrocarbon compositions reported from hydrate and sediment samples; and (3) the calculation of the thickness of the hydrate stability zone. The solution is based on a transcendental function, which is solved iteratively in a GIS environment.The model has been applied in the northernmost continental slope of the Gulf of Cadiz, an area where an abundant supply for hydrate formation, such as extensive hydrocarbon seeps, diapirs and fault structures, is combined with deep undercurrents and a complex seafloor morphology. In the Gulf of Cadiz, the model depicts the distribution of the base of the gas hydrate stability zone for both biogenic and thermogenic gas compositions, and explains the geometry and distribution of geological structures derived from gas venting in the Tasyo Field (Gulf of Cadiz) and the generation of BSR levels on the upper continental slope.     

A control-volume finite-element method for three-dimensional multiphase basin modeling

In this paper we describe a 3D control-volume finite-element method to solve numerically the coupled partial differential equations (PDEs) governing geological processes involved in the evolution of sedimentary basins. These processes include sediment deposition and deformation, hydrocarbon generation, multiphase fluid flow, and heat transfer in deforming porous media.     

Velocity potential formulations of highly accurate Boussinesq-type models

The highly accurate Boussinesq-type equations of Madsen et al. (Madsen, P.A., Bingham, H.B., Schäffer, H.A., 2003. Boussinesq-type formulations for fully nonlinear and extremely dispersive water waves: Derivation and analysis. Proc. R. Soc. Lond. A 459, 1075–1104; Madsen, P.A., Fuhrman, D.R., Wang, B., 2006. A Boussinesq-type method for fully nonlinear waves interacting with a rapidly varying bathymetry. Coast. Eng. 53, 487–504); Jamois et al. (Jamois, E., Fuhrman, D.R., Bingham, H.B., Molin, B., 2006. Wave-structure interactions and nonlinear wave processes on the weather side of reflective structures. Coast. Eng. 53, 929–945) are re-derived in a more general framework which establishes the correct relationship between the model in a velocity formulation and a velocity potential formulation. Although most work with this model has used the velocity formulation, the potential formulation is of interest because it reduces the computational effort by approximately a factor of two and facilitates a coupling to other potential flow solvers. A new shoaling enhancement operator is introduced to derive new models (in both formulations) with a velocity profile which is always consistent with the kinematic bottom boundary condition. The true behaviour of the velocity potential formulation with respect to linear shoaling is given for the first time, correcting errors made by Jamois et al. (Jamois, E., Fuhrman, D.R., Bingham, H.B., Molin, B., 2006. Wave-structure interactions and nonlinear wave processes on the weather side of reflective structures. Coast. Eng. 53, 929–945). An exact infinite series solution for the potential is obtained via a Taylor expansion about an arbitrary vertical position z = zˆ. For practical implementation however, the solution is expanded based on a slow variation of zˆ and terms are retained to first-order. With shoaling enhancement, the new models obtain a comparable accuracy in linear shoaling to the original velocity formulation. General consistency relations are also derived which are convenient for verifying that the differential operators satisfy a potential flow and/or conserve mass up to the order of truncation of the model. The performance of the new formulation is validated using computations of linear and nonlinear shoaling problems. The behaviour on a rapidly varying bathymetry is also checked using linear wave reflection from a shelf and Bragg scattering from an undulating bottom. Although the new models perform equally well for Bragg scattering they fail earlier than the existing model for reflection/transmission problems in very deep water.     

Two-dimensional,two-phase granular sediment transport model with applications to scouring downstream of an apron

We present a two-dimensional, two-phase model for non-cohesive sediment transport. This model solves concentration-weighted averaged equations of motion for both fluid and sediment phases. The model accounts for the interphase momentum transfer by considering drag forces. A collisional theory is used to compute the sediment stresses, while a two-equation (k–ε) fluid turbulence closure is implemented. A benchmark sediment transport problem concerning the scouring downstream of an apron is carried out as an example and numerical results agree with existing experimental data.     

Molecular ecological study of <Emphasis Type="Italic">Siganus spinus</Emphasis> and <Emphasis Type="Italic">S. guttatus</Emphasis> from Okinawan waters based on mitochondrial DNA control region sequences

The mitochondrial DNA (mtDNA) control region was sequenced to examine the genetic variability and gene flow of juvenile Siganus spinus and S. guttatus. In total, 461 nucleotide sequences were obtained from 69 specimens of juvenile S. spinus from Okinawa Island and Ishigaki Island, in which 17 variable sites and 22 haplotypes were identified. Haplotype diversity (h) was high (0.9244 in Okinawa and 0.8984 in Ishigaki), whereas nucleotide diversity (π) was low (0.0063 in Okinawa and 0.0059 in Ishigaki). The two populations were not genetically distinct. Siganus guttatus, which do not form large schools at recruitment, in contrast S. spinus, were also analyzed by studying 152 individuals collected off Okinawa Island, Miyako Island, and Ishigaki Island. Of 476 nucleotide sequences, 50 were variable, and 42 haplotypes were identified. Genetic variability values were high (h = 0.8766 and π = 0.0151 in Okinawa; h = 0.9640 and π = 0.0192 in Miyako; h = 0.9161 and π = 0.0199 in Ishigaki). The Okinawa population was genetically isolated from the Miyako and Ishigaki populations. As a result, genetic diversity was high for each of these siganid populations despite their being target species for fisheries; however, the degree of inter-population gene flow was higher for S. spinus than S. guttatus, suggesting that these species exhibit different dispersal strategies.     

Decomposition analysis on direct material input and dematerialization of mining cities in Northeast China

Material dematerialization is a basic approach to reduce the pressure on the resources and environment and to realize the sustainable development. The material flow analysis and decomposition method are used to calculate the direct material input (DMI) of 14 typical mining cities in Northeast China in 1995-2004 and to analyze the demateri-alization and its driving factors in the different types of mining cities oriented by coal, petroleum, metallurgy and multi-resources. The results are as follows: 1) from 1995 to 2006, the increase rates of the DMI and the material input intensity of mining cities declined following the order of multi-resources, metallurgy, coal, and petroleum cities, and the material utilizing efficiency did following the order of petroleum, coal, metallurgy, and multi-resources cities; 2) during the research period, all the kinds of mining cities were in the situation of weak sustainable development in most years; 3) the pressure on resources and environment in the multi-resources cities was the most serious; 4) the petro-leum cities showed the strong trend of sustainable development; and 5) in recent years, the driving function of eco-nomic development for material consuming has continuously strengthened and the controlling function of material utilizing efficiency for it has weakened. The key approaches to promote the development of circular economy of min-ing cities in Northeast China are put forward in the following aspects: 1) to strengthen the research and development of the technique of resources' cycling utilization, 2) to improve the utilizing efficiency of resources, and 3) to carry out the auditing system of resources utilization.     

Soil evolution features of debris flow waste-shoal land

The reclamation and utilization of debris flow waste-shoal land plays an important role in the mitigation and control of debris flow hazards, which thus contributes a lot to the exploitation of insufficient land resources in mountainous areas and the reduction of losses caused by debris flow. The aim of this paper is to discuss the features and mechanism of soil evolution of debris flow waste-shoal land so as to search for the available modes of its reclamation and utilization. The Jiangjiagou Ravine, a typical debris flow ravine, was selected to study soil evolution features of debris flow waste-shoal land based on the analysis of soil physieochemical properties and soil microstructure. It was found that the soil evolution rates of debris flow waste-shoal land varied with different modes of reclamation. For the land which had been reclaimed for less than lO years, soil evolved most rapidly in paddy fields, and more rapidly in dry farmland than in naturally restored waste-shoal land. For the land which had been used for more than lo years, the soil evolution rates of dry farmland, naturally restored waste-shoal land and paddy farmland decreased in the file. For the same utilization period of time, significant differences were recognized in soil evolution features under different modes of reclamation. Analysis data showed that soil clay content, soil thickness, the psephicity of skeleton particles and contents of microaggregates （〈0.02 mm） in paddy farmland were all highest. Soil nutrients and porosity of dry farmland were better than those of paddy farmland and naturally restored waste-shoal land, and those of paddy farmland were superior to those of naturally restored waste-shoal land. Paddy farmland characterized by rapid pedogenesis, stable evolution and high utilizability was the priority candidate for the reclamation and utilization of debris flow waste -shoal land.     

Regional hydrologic mapping of flows in stream networks

River low flow (LF) is an important hydrological characteristic used in management of the quantity and quality of water resources. A common way to transpose at-station low flow quantiles to ungaged locations involves development of regional multivariate regression models. In site selection studies of water control facilities, need arises to automate the mapping of LF characteristics onto the stream network in a continuous manner. This would allow estimation of LF at any desired ungaged location. A mapping algorithm is introduced that determines the value of regression model's independent variables (input components) over the drainage area of each stream pixel and calculates the flow (output). The mapping algorithm relies mainly on the digital elevation model (DEM) and its derivatives, such as flow direction and flow accumulation. Furthermore, the contribution of each independent variable of the regional model to the total flow may be plotted to represent the flow profiles along the streams. A case study involving LF mapping in rivers of Gilan province, Iran, where LF quantiles are required for issuing water withdrawal permits as well as maintaining water quality standards, is also described in the paper. Overall, the maps and profiles of flow statistics in the region of interest provide convenient visualization and assessment tools for water resource and environmental engineers. The transposition algorithm may also be applied in regional mapping of other flow characteristics such as flood or average flows.