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).
The behaviour of porous sinters, during compression and compression with reverse cyclic torsion tests is investigated in the article based on the combination of experimental and numerical techniques. The sinters manufactured from the Distaloy AB powder are examined. First, series of simple uniaxial compression tests were performed on samples with three different porosity volume fractions: 15, 20 and 25%. Obtained data were then used during identification procedure of the Gurson-Tvergaard-Needleman finite element based model, which can capture influence of porosity evolution on plasticity. Finally, the identified Gurson-Tvergaard- Needleman model was validated under complex compression with reverse cyclic torsion conditions and proved its good predictive capabilities. Details on both experimental and numerical investigations are presented within the paper.
Detailed studies of the movement of liquid steel (hydrodynamics) on a real object are practically impossible. The solution to this problem are physical modelling carried out on water models and numerical modelling using appropriate programs. The method of numerical modelling thanks to the considerable computing power of modern computers gives the possibility of solving very complex problems. The paper presents the results of model tests of liquid flow through tundish. The examined object was model of the twonozzle tundish model. The ANSYS Fluent program was used to describe the behavior of liquid in the working area of the tundish model. Numerical simulations were carried out using two numerical methods of turbulence description: RANS (Reynolds-Averaged Navier-Stokes) – model k-ε and LES (Large Eddy Simulation). The results obtained from CFD calculations were compared with the results obtained using the water model.
The presented results of investigations are part of a larger study focused on the optimization of the flow and mixing of liquid steel in the industrial tundish of continuous casting machine. The numerical simulations were carried out concern the analysis of hydrodynamic conditions of liquid steel flow in a tundish operating in one of the national steelworks. Numerical simulations were performed using the commercial code ANSYS Fluent. The research concerns two different speeds of steel casting. In real conditions, these speeds are the most commonly used in the technological process when casting two different groups of steel. As a result of computational fluid dynamics (CFD) calculations, predicted spatial distributions of velocity and liquid steel turbulence fields and residence time distribution (RTD) curves were obtained. The volume fractions of different flows occurring in the tundish were also calculated. The results of the research allowed a detailed analysis of the influence of casting speed on the formation of hydrodynamic conditions prevailing in the reactor.
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.
New York Bay is one of the most important transition regions of ships trading to east America. The region plays an important role in the commerce of the New York metropolitan area. The area is surrounded with the coasts that have various levels of environmental sensitivity. The area accommodates high diversity of native ecosystems and species that are rather vulnerable in case of oil spill. Thus getting well informed about the likelihood, or fate, of oil spills around this region is of great importance so that proactive measures can be taken. The purpose of this study is to investigate the oil spill and predict the future accidents likely to be encountered around the Bay of New York. Two trajectory models have been conducted for the study. ADIOS (Automated Data Inquiry for Oil Spills), has been conducted for natural degradation calculations, and, GNOME (General NOAA Operational Modeling Environment), has been conducted for surface spread simulation. The results gained through these efforts are hoped to be useful for many organizations dealing with oil spill response operations and contribute to an effective and efficient coordination among the relevant institutions.
The paper evaluates two approaches of numerical modelling of solidification of continuously cast steel billets by finite element method, namely by the numerical modelling under the Steady-State Thermal Conditions, and by the numerical modelling with the Traveling Boundary Conditions. In the paper, the 3D drawing of the geometry, the preparation of computational mesh, the definition of boundary conditions and also the definition of thermo-physical properties of materials in relation to the expected results are discussed. The effect of thermo-physical properties on the computation of central porosity in billet is also mentioned. In conclusion, the advantages and disadvantages of two described approaches are listed and the direction of the next research in the prediction of temperature field in continuously cast billets is also outlined.
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 this paper, a new simple method for determination of flow parameters, axial dispersion coefficients DL and Péclet numbers Pe was presented. This method is based on an accurate measurement model considering pulse tracer response. Our method makes it possible to test the character of gas flow motion and precisely measure flow parameters for different pressures and temperatures. The idea of combining the transfer function, numerical inversion of the Laplace transform and optimisation method gives many benefits like a simple and effective way of finding solution of inverse problem and model coefficients. The calculated values of flow parameters (DL and/or Pe) suggest that in the considered case the gas flow is neither plug flow nor perfect mixing under operation condition. The obtained outcomes agree with the gas flow theory. Calculations were performed using the CAS program type, Maple®.
In the paper the thermal processes proceeding in the solidifying metal are analyzed. The basic energy equation determining the course of solidification contains the component (source function) controlling the phase change. This component is proportional to the solidification rate ¶ fS/¶ t (fS Î [0, 1], is a temporary and local volumetric fraction of solid state). The value of fS can be found, among others, on the basic of laws determining the nucleation and nuclei growth. This approach leads to the so called micro/macro models (the second generation models). The capacity of internal heat source appearing in the equation concerning the macro scale (solidification and cooling of domain considered) results from the phenomena proceeding in the micro scale (nuclei growth). The function fS can be defined as a product of nuclei density N and single grain volume V (a linear model of crystallization) and this approach is applied in the paper presented. The problem discussed consists in the simultaneous identification of two parameters determining a course of solidification. In particular it is assumed that nuclei density N (micro scale) and volumetric specific heat of metal (macro scale) are unknown. Formulated in this way inverse problem is solved using the least squares criterion and gradient methods. The additional information which allows to identify the unknown parameters results from knowledge of cooling curves at the selected set of points from solidifying metal domain. On the stage of numerical realization the boundary element method is used. In the final part of the paper the examples of computations are presented.
In the paper the use of the artificial neural network to the control of the work of heat treating equipment for the long axisymmetric steel elements with variable diameters is presented. It is assumed that the velocity of the heat source is modified in the process and is in real time updated according to the current diameter. The measurement of the diameter is performed at a constant distance from the heat source (∆z = 0). The main task of the model is control the assumed values of temperature at constant parameters of the heat source such as radius and power. Therefore the parameter of the process controlled by the artificial neural network is the velocity of the heat source. The input data of the network are the values of temperature and the radius of the heated element. The learning, testing and validation sets were determined by using the equation of steady heat transfer process with a convective term. To verify the possibilities of the presented algorithm, based on the solve of the unsteady heat conduction with finite element method, a numerical simulation is performed. The calculations confirm the effectiveness of use of the presented solution, in order to obtain for example the constant depth of the heat affected zone for the geometrically variable hardened axisymmetric objects.