This paper presents a complex study of anhydrite interbeds influence on the cavern stability in the Mechelinki salt deposit. The impact of interbeds on the cavern shape and the stress concentrations were also considered. The stability analysis was based on the 3D numerical modelling. Numerical simulations were performed with use of the Finite Difference Method (FDM) and the FLAC3D v. 6.00 software. The numerical model in a cuboidal shape and the following dimensions: length 1400, width 1400, height 1400 m, comprised the part of the Mechelinki salt deposit. Three (K-6, K-8, K-9) caverns were projected inside this model. The mesh of the numerical model contained about 15 million tetrahedral elements. The occurrence of anhydrite interbeds within the rock salt beds had contributed to the reduction in a diameter and irregular shape of the analysed caverns. The results of the 3D numerical modelling had indicated that the contact area between the rock salt beds and the anhydrite interbeds is likely to the occurrence of displacements. Irregularities in a shape of the analysed caverns are prone to the stress concentration. However, the stability of the analysed caverns are not expected to be affected in the assumed operation conditions and time period (9.5 years).
This article describes stability issues of main excavations in deep copper mines in Poland, from the perspective of mining work safety. To protect main transportation and ventilation routes, parts of rock are left untaken to form so-called protective pillars. The problem was to determine the size of main excavations protective pillars in deep underground copper mines in which provide stability of main excavations. The results of numerical simulations of the stability of protective pillars under specific geological and mining conditions are presented, covering: underground depth and width of protective pillar, number, size and layout geometry of protected excavations, as well as the impact of parameters of surrounding gob areas. Problem was solved applying numerical simulations based on the finite element method which were performed in a plane state of strain by means of Phase2 v. 8.0 software. The behavior of the rock mass under load was described by an elastic-plastic model. The Mohr-Coulomb criterion was used to assess the stability of the rock mass. The results of numerical modeling have practical applications in the designing of protective pillars primarily in determining their width. These results were used to prepare new guidelines for protective pillars in Polish copper mines in the Legnica-Glogow Copper District.
The paper presents a numerical model of the novel design of the axial magnetic bearing with six cylindrical poles. The motivation behind this idea was to eliminate vibrations in rotating machinery due to the axial load. Common conception of such a bearing provides a single component of the electromagnetic force, which is not enough to reduce transverse and lateral vibrations of the armature. The proposed design allows for avoiding wobbling of the disc with the use of a few axial force components that are able to actively compensate the axial load and stabilise the disc in a balanced position. Before a real device is manufactured, a virtual prototype should be prepared. The accurate numerical model will provide essential knowledge about the performance of the axial magnetic bearing.
This paper constitutes the sensitivity study of application the Polar WRF model to the Svalbard area with testing selected parameterizations, including planetary boundary layer, radiation and microphysics schemes. The model was configured, using three one-way nested domains with 27 km, 9 km and 3 km grid cell resolutions. Results from the innermost domain were presented and compared against measured wind speed and air temperature at 10 meteorological stations. The study period covers two months: June 2008 and January 2009. Significant differences between simulations results occurred for planetary boundary layer (PBL) schemes in January 2009. The Mellor-Yamada-Janjic (MYJ) planetary boundary layer (PBL) scheme resulted in the lowest errors for air temperature, according to mean error (ME), mean absolute error (MAE) and correlation coefficient values, where for wind speed this scheme was the worst from all the PBL schemes tested. In the case of June 2008, shortwave and longwave radiation schemes influenced the results the most. Generally, higher correlations were obtained for January, both for air temperature and wind speed. However, the model performs better for June in terms of ME and MAE error statistics. The results were also analyzed spatially, to summarize the uncertainty of the model results related to the analyzed parameterization schemes groups. Significant variability among simulations was calculated for January 2009 over the northern part of Spitsbergen and fjords for the PBL schemes. Standard deviations for monthly average simulated values were up to 3.5°C for air temperature and around 1 m s-1 for wind speed.
A numerical model of binary alloy crystallization, based on the cellular automaton technique, is presented. The model allows to follow the crystallization front movement and to generate the images of evolution of the dendritic structures during the solidification of a binary alloy. The mathematic description of the model takes into account the proceeding thermal, diffusive, and surface phenomena. There are presented the results of numerical simulations concerning the multi-dendritic growth of solid phase along with the accompanying changes in the alloying element concentration field during the solidification of Al + 5% wt. Mg alloy. The model structure of the solidified casting was achieved and compared with the actual structure of a die casting. The dendrite interaction was studied with respect to its influence on the generation and growth of the primary and secondary dendrite arms and on the evolution of solute segregation both in the liquid and in the solid state during the crystallization of the examined alloy. The morphology of a single, free-growing dendritic crystal was also modelled. The performed investigations and analyses allowed to state e.g. that the developed numerical model correctly describes the actual evolution of the dendritic structure under the non-equilibrium conditions and provides for obtaining the qualitatively correct results of simulation of the crystallization process.
This work presents the results of numerical modeling of Karman vortex street generation performed with ANSYS/FLUENT package application. The influence of the mechanical elements located downstream of the bluff body on the vortex frequency has been found during earlier laboratory investigations. Five various geometrical configurations have been tested. Considerable differences in pictures of distributions of pressure, horizontal and vertical velocities have appeared for various configurations. Qualitative as well as quantitative results are presented in the paper. They confirm the significant dependence of the Karman vortex street parameters on the meter configuration.
The generic mathematical model and computational algorithm considering hydrodynamics, heat and mass transfer processes during casting and forming steel ingots and castings are offered. Usage domains for turbulent, convective and non-convective models are determined depending on ingot geometry and thermal overheating of the poured melt. The expert system is developed, enabling to choose a mathematical model depending on the physical statement of a problem.
The main purpose of the present work was to validate the numerical model for the pulse-step liquid steel alloying method using a physical simulator that enables the observation and recording of phenomena occurring during the continuous steel casting process. The facility under investigation was a single-nozzle tundish equipped with a dam. To physical trials the glass water model was made on a scale of 2:5. For the mathematical description of turbulence during liquid steel alloying process, the k-ε and k-ω models were employed in the simulations. Based on the computer simulations and physical trials carried out, alloy addition behaviour and mixing curves for different tundish alloy addition feeding positions were obtained. The change in the location of alloy addition feeding to the liquid steel had an effect on the process of alloy addition spread in the liquid steel bulk and on the mixing time.