Exploitation of hard coal seams by roadway system is applied by two coal mines in southern Poland in Upper Silesian Basin. It is a secondary mining exploitation carries out in safety pillars of urban areas and shafts within mining areas of closed coal mines. Roadway system is the excavation process of gateways which are made in parallel order leaving coal pillars between them. An optimal width of coal pillar makes roadway stable and reduces subsidence of terrain surface. The article presents results of subsidence simulation caused by partial extraction using empirical and numerical methods on the example of one exploitation field of “Siltech” coal mine. The asymptotic state of subsidence was considered after mining ceased in the study area. In order to simulate of subsidence, numerical model of rock mass and model of Knothe-Budryk theory were calibrated. Simulation of vertical displacements in numerical method was carried out using RS3 program by Rocscience based on finite element method. The assumption was made that model of rock mass is transversely isotropic medium, in which panels were designed according to order of extraction of coal seams. The results of empirical and numerical methods were compared with measured values of subsidence at benchmarks along drawn lines (subsidence profiles).
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 cohesion and internal friction angle were characterized as quadratic functions of strain and were assumed to follow the Mohr-Coulomb criterion after the yield of peak strength. These mechanical parameters and their variations in post-peak softening stage can be exactly ascertained through the simultaneous solution based on the data points of stress-strain curves of triaxial compression tests. Taking the influence of the fault into account, the variation of strata pressure and roadway convergence with coal advancement, the temporal and spatial distribution of axial bolt load were numerically simulated by FLAC3D (Fast Lagrangian Analysis of Continua) using the ascertained post-peak mechanical parameters according to the cohesion weakening and friction strengthening model. The change mechanism of axial load of single rock bolt as abutment pressure changes was analyzed, through the comparison analysis with the results of axial bolt load by field measurements at a coal mine face. The research results show that the simulated results such as the period of main roof weighting, temporal and spatial distribution of axial bolt load are in accordance with field measurement results, so the validity of the numerical model is testified. In front of the working face, the front abutment pressure increases first and then decreases, finally tends to be stable. A corresponding correlation exists between the variation of axial bolt load and rock deformation along the bolt body. When encountered by a fault, the maximum abutment pressure, the influential range of mining disturbance and the roadway convergence between roof and floor before the working face are all increased. In the roadways along the gob, axial bolt loads on the side of the working face decrease, while the other side one increases after the collapse of the roof. As superficial surrounding rock mass is damaged, the anchoring force of rock bolts will transfer to inner rock mass for balancing the tensile load of the bolts.
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.
Passive autocatalytic recombiners (PAR) is the only used method for hydrogen removal from the containment buildings in modern nuclear reactors. Numerical models of such devices, based on the CFD approach, are the subject of this paper. The models may be coupled with two types of computer codes: the lumped parameter codes, and the computational fluid dynamics codes. This work deals with 2D numerical model of PAR and its validation. Gaseous hydrogen may be generated in water nuclear reactor systems in a course of a severe accident with core overheating. Therefore, a risk of its uncontrolled combustion appears which may be destructive to the containment structure.
The main purpose of this article is to verify and validate the mathematical description of the airflow around a wind turbine with vertical axis of rotation, which could be considered as representative for this type of devices. Mathematical modeling of the airflow around wind turbines in particular those with the vertical axis is a problematic matter due to the complex nature of this highly swirled flow. Moreover, it is turbulent flow accompanied by a rotation of the rotor and the dynamic boundary layer separation. In such conditions, the key aspects of the mathematical model are accurate turbulence description, definition of circular motion as well as accompanying effects like centrifugal force or the Coriolis force and parameters of spatial and temporal discretization. The paper presents the impact of the different simulation parameters on the obtained results of the wind turbine simulation. Analysed models have been validated against experimental data published in the literature.
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.
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.
The paper presents a dynamic analysis of the damaged masonry building repaired with the Flexible Joint Method. Numerical analysis helped to determine the effect of the applied repairing method on natural frequencies as well as values of stresses and accelerations in the analyzed variants of numerical model. They confirmed efficiency of the proposed repair method.
The Dez dam was commissioned in 1963 and since sediments accumulated in the reservoir up to an elevation of approximately 15m below the intake of the power tunnel. One of the possible measures to improve operation of the reservoir is by heightening of the existing dam. This paper describes the conducted procedure for static and thermal calibration of this 203m dam in Iran based on micro geodesies measurements. Also the nonlinear response of existing dam is investigated under maximum credible earthquake ground motions considering joint behavior and mass concrete cracking and safety of dam is evaluated for possible heightening. For thermal calibration of provided numerical model, transient thermal analysis was conducted and results were compared with thermometers records installed in central block. In addition, for static calibration; thermal distribution within dam body, dam self weight, hydrostatic pressure and silt load applied on the 3D fi nite element model of dam-reservoir-foundation were considered. Results show that the distribution of stresses will be critical within dam for heightening case under seismic loads in MCL.
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.
In many therapeutic applications of a pulsed focused ultrasound with various intensities the finite- amplitude acoustic waves propagate in water before penetrating into tissues and their local heating. Water is used as the matching, cooling and harmonics generating medium. In order to design ultrasonic probes for various therapeutic applications based on the local tissue heating induced in selected organs as well as to plan ultrasonic regimes of treatment a knowledge of pressure variations in pulsed focused nonlinear acoustic beams produced in layered media is necessary. The main objective of this work was to verify experimentally the applicability of the recently developed numerical model based on the Time- Averaged Wave Envelope (TAWE) approach (Wójcik et al., 2006) as an effective research tool for predicting the pulsed focused nonlinear fields produced in two-layer media comprising of water and tested materials (with attenuation arbitrarily dependent on frequency) by clinically relevant axially-symmetric therapeutic sources. First, the model was verified in water as a reference medium with known linear and nonlinear acoustic properties. The measurements in water were carried out at a 25◦C temperature using a 2.25 MHz circular focused (f/3.0) transducer with an effective diameter of 29 mm. The measurement results obtained for 8-cycle tone bursts with three different initial pressure amplitudes varied between 37 kPa and 113 kPa were compared with the numerical predictions obtained for the source boundary condition parameters determined experimentally. The comparison of the experimental results with those simulated numerically has shown that the model based on the TAWE approach predicts well both the spatial-peak and spatial-spectral pressure variations in the pulsed focused nonlinear beams produced by the transducer used in water for all excitation levels complying with the condition corresponding to weak or moderate source-pressure levels. Quantitative analysis of the simulated nonlinear beams from circular transducers with ka ≫ 1 allowed to show that the axial distance at which sudden accretion of the 2nd or higher harmonics amplitude appears is specific for this transducer regardless of the excitation level providing weak to moderate nonlinear fields. For the transducer used, the axial distance at which the 2nd harmonics amplitude suddenly begins to grow was found to be equal to 60 mm. Then, the model was verified experimentally for two-layer parallel media comprising of a 60-mm water layer and a 60-mm layer of 1.3-butanediol (99%, Sigma-Aldrich Chemie GmbH, Steinheim, Germany). This medium was selected because of its tissue-mimicking acoustic properties and known nonlinearity parameter B/A. The measurements of both, the peak- and harmonic-pressure variations in the pulsed nonlinear acoustic beams produced in two-layer media (water/1.3-butanediol) were performed for the same source boundary conditions as in water. The measurement results were compared with those simulated numerically. The good agreement between the measured data and numerical calculations has shown that the model based on the TAWE approach is well suited to predict both the peak and harmonic pressure variations in the pulsed focused nonlinear sound beams produced in layered media by clinically relevant therapeutic sources. Finally, the pulsed focused nonlinear fields from the transducer used in two-layer media: water/castor oil, water/silicone oil (Dow Corning Ltd., Coventry, UK), water/human brain and water/pig liver were predicted for various values of the nonlinearity parameter of tested media.
Balanced distribution of air in coal-fired boiler is one of the most important factors in the combustion process and is strongly connected to the overall system efficiency. Reliable and continuous information about combustion airﬂow and fuel rate is essential for achieving optimal stoichiometric ratio as well as efficient and safe operation of a boiler. Imbalances in air distribution result in reduced boiler efficiency, increased gas pollutant emission and operating problems, such as corrosion, slagging or fouling. Monitoring of air ﬂow trends in boiler is an effective method for further analysis and can help to appoint important dependences and start optimization actions. Accurate real-time monitoring of the air distribution in boiler can bring economical, environmental and operational benefits. The paper presents a novel concept for online monitoring system of air distribution in coal-ﬁred boiler based on real-time numerical calculations. The proposed mathematical model allows for identification of mass ﬂow rates of secondary air to individual burners and to overﬁre air (OFA) nozzles. Numerical models of air and flue gas system were developed using software for power plant simulation. The correctness of the developed model was verified and validated with the reference measurement values. The presented numerical model for real-time monitoring of air distribution is capable of giving continuous determination of the complete air flows based on available digital communication system (DCS) data.
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.
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 paper presents the results of an analysis of gaseous sensors based on a surface acoustic wave (SAW) by means of the equivalent model theory. The applied theory analyzes the response of the SAW sensor in the steady state affected by carbon monoxide (CO) in air. A thin layer of WO3 has been used as a sensor layer. The acoustical replacing impedance of the sensor layer was used, which takes into account the profile of the concentration of gas molecules in the layer. Thanks to implementing the Ingebrigtsen equation, the authors determined analytical expressions for the relative changes of the velocity of the surface acoustic wave in the steady state. The results of the analysis have shown that there is an optimum thickness of the layer of CO sensor at which the acoustoelectric effect (manifested here as a change in the acoustic wave velocity) is at its highest. The theoretical results were verified and confirmed experimentally