Advances in discrete element modelling of underground excavations

The paper presents advances in the discrete element modelling of underground excavation processes extending modelling possibilities as well as increasing computational efficiency. Efficient numerical models have been obtained using techniques of parallel computing and coupling the discrete element method with finite element method. The discrete element algorithm has been applied to simulation of different excavation processes, using different tools, TBMs and roadheaders. Numerical examples of tunnelling process are included in the paper, showing results in the form of rock failure, damage in the material, cutting forces and tool wear. Efficiency of the code for solving large scale geomechanical problems is also shown.     

Ammonia Nitrogen at the Water–Sediment Interface in Puck Bay (Baltic Sea)

The magnitude of the exchange flux at the water–sediment interface was determined on the basis of the ammonia concentration gradient at the near-bottom water–interstitial interface and Fick's first law. It was established that in Puck Bay, ammonia almost always passes from the sediment to water. Ammonia flux varied from 5 to 1434 μmol NH₄-N m⁻² day⁻¹. In total,c. 138·2 tonneammonia year⁻¹pass from sediments of Internal Puck Bay to near-bottom water, the equivalent value for External Puck Bay being 686·9 tonne year⁻¹. In total, about 825 tonne ammonia year⁻¹passes from the sediment to near-bottom water of Puck Bay. In interstitial waters, ammonia occurred in concentrations varying over a wide range (3–1084 μmol NH₄-N dm⁻³).The basic factors affecting the magnitude of ammonia concentration in interstitial waters included: oxidation of organic matter, type of sediment, and inflow of fresh underground waters to the region examined.This paper involves preliminary studies only and constitutes a continuation of the studies on ionic macrocomponents and phosphorus in interstitial waters of Puck Bay undertaken previously.     

Spatial distribution of phosphorus forms in sediments in the Gulf of Gdańsk (southern Baltic Sea)

The paper analyses the concentrations of total phosphorus and its forms in sediments from the Gulf of Gdańsk on the basis of studies conducted at 25 sampling sites in 2001–05. The phosphorus speciation analysis was performed by sequential extraction. The extensive spatial variability of P_tot concentrations and speciation was found to be dependent on the physicochemical properties of the sediments, the oxygen conditions in the water and sediments, and the depth of the water column above the sediment surface. In the coastal zone, the sedimentation of riverine suspended matter and the sorption and chemisorption processes exert a considerable influence on P speciation. Over 70% of variation of total phosphorus concentration in sediments in the Gulf of Gdańsk could be explained by changes of proportion of fine fraction of sediments (grain size <0.0625 mm). Maximum P_tot concentrations were recorded in clays and silts in the deep water, stratified part of the Gulf of Gdańsk. In the coastal zone, where sandy sediments are dominant, phosphorus concentrations were much lower; this was due to the considerable dynamics of the bottom water and intensive sea floor transport. P_tot concentrations in the Gulf of Gdańsk sediments ranged from 1.75 to 957.17 μmol g⁻¹ d.w. Of all the forms of phosphorus, the highest concentrations were found for organic phosphorus (Org-P). Of its inorganic forms, the highest concentrations were of phosphorus bound to clay minerals and aluminium oxides (NaOH-P), the lowest ones were of loosely bound phosphorus (NaCl-P). On the basis of determinations of total phosphorus concentrations in sediments of a given type and the available data on the seabed areas covered by particular sediment types in the Gulf of Gdańsk, the mass of total phosphorus in the surficial sediment layer (0–2 cm) was estimated at ca. 15.6×10³ tonnes.     

Slope processes in the High Arctic: a 30‐year‐long record from Spitsbergen

Ten monitoring sites were established in 1986 near the Scottbreen and Renardbreen glaciers (Bellsund, Spitsbergen) in order to determine the rates of gelifluction on raised marine terraces and on talus cones, accumulation on talus cones, and movement of the protalus rampart. The measurements were performed in 1987, 2007 and 2016. The results obtained revealed a very slow rate of changes on the surface of the analysed area. During 30 years, only one out of 84 points where the gelifluction rate was measured was displaced by 15 cm down‐slope. The maximal annual and average gelifluction rates were calculated at 0.5 cm and 0.18 cm a^?1, respectively. Several measurement points did not change their position over the observation period. As shown by the observations, the gelifluction rate largely does not depend on slope inclination. Studies of talus cones allowed the level of cone increment to be determined as 314.3 kg of scree per square metre. Analysis of factors influencing the rate of mass movement allowed for correlation of the small‐scale gelifluction movement at the analysed sites with the general trends of climate warming. In this part of the Arctic, these trends are manifested by soil desiccation. It appears that slope processes depend on very local topoclimatic factors. Analysis of the obtained data with regard to palaeogeographical and climatic investigations suggests that the ice segregation forms on raised marine terraces developed during the last Holocene cooling, i.e. the Little Ice Age.     

Comparison of continuous and discontinuous constitutive models to simulate concrete behaviour under mixed‐mode failure conditions

The paper presents detailed FE simulation results of concrete elements under mixed‐mode failure conditions according to the so‐called shear‐tension test by Nooru‐Mohamed, characterized by curved cracks. A continuous and discontinuous numerical two‐dimensional approach was used. In order to describe the concrete's behaviour within continuum mechanics, two different constitutive models were used. First, an elasto‐plastic model with isotropic hardening and softening was assumed. In a compression regime, a Drucker–Prager criterion with a non‐associated flow rule was used. In turn, in a tensile regime, a Rankine criterion with an associated flow rule was adopted. Second, an isotropic damage constitutive model was applied with a single scalar damage parameter and different definitions of the equivalent strain. Both constitutive laws were enriched by a characteristic length of micro‐structure to capture properly strain localization. As an alternative approach, the extended finite element method was used. Our results were compared with the experimental ones and with results of other FE simulations reported in the literature. Copyright © 2015 John Wiley & Sons, Ltd.     

The assessment of environmental pollution caused by mining and metallurgy wastes from highly polluted post-industrial regions in Southern Poland

Stored metallurgy and mining wastes contain relatively high amounts of potentially toxic elements. To monitor the distribution of contaminants originating from dumps, the chemical and physical properties of the wastes must be characterised. In this study, the chemical properties of wastes deposited in two different locations in Southern Poland (Szklary and Zloty Stok) were evaluated. Heaps located in Zloty Stok contain wastes from gold mineralisation comprising arsenic while wastes in Szklary originate from a factory that produced an iron-nickel alloy. In Szklary the total concentrations of Ca, Mg, Fe, Zn, Mn, Cr, Co, Cu, Ni, Tl, Ag, Cd and Pb were determined, while in Zloty Stok also As is an important contaminant. To assess the risk of contamination of the surrounding environment and to select the proper method for removing the contaminants, information on the distribution of elements between operationally defined phases must be obtained. For this purpose, a six-step sequential extraction was used. The mobility of most elements in the wastes from Szklary and Zloty Stok was relatively low; however, the large amount of As, 40–180 mg L^?1, that could be released to environment in case of Zloty Stok was high. The results of fractionation studies indicated that Ag, Cd, Tl, Co, Cu, Cr, Pb, Zn and Mn may be released into environment under low pH and low redox potential conditions, which can be induced by bacterial activity.     

Investigation of X-Ray Flares with Long Rising Phases

We have used Yohkoh and GOES X-ray observations to investigate flares with a long rising phase. We have found that a characteristic feature of such flares is a long time interval, Δ t ≥ 20 min, between the temperature maximum and the maximum of the emission measure. We have carried out detailed analysis for 10 limb flares of this type. Time variation of the heating function, E_H(t), has been determined for their loop-top X-ray kernels. The time variation of E_H(t), together with the temperature–density diagnostic diagrams, have been used to explain the large value of the time interval, Δ t. The main point is that for these flares the heating function E_H(t) decreases so slowly after the temperature maximum, that for the long time, Δ t, the energy flux reaching flare foot points is sufficient to maintain significant chromospheric evaporation. Investigation of the flare evolution in the temperature–density diagnostic diagrams allowed us to work out a new method of determination of the density for flare kernels. This method can be applied to all the kernels for which their altitudes can be estimated. The advantage of this method is that for the density determination it is not necessary to assume what is the extension of the emitting plasma along the line of sight.