The mathematical model and numerical simulations of the solidification of a cylindrical shaped casting, which take into account the process of filling the mould cavity by liquid metal and feeding the casting through the riser during its solidification, are presented in the paper. Mutual dependence of thermal and flow phenomena were taken into account because have an essential influence on solidification process. The effect of the riser shape on the effectiveness of feeding of the solidifying casting was determined. In order to obtain the casting without shrinkage defects, an appropriate selection of riser shape was made, which is important for foundry practice. Numerical calculations of the solidification process of system consisting of the casting and the conical or cylindrical riser were carried out. The velocity fields have been obtained from the solution of momentum equations and continuity equation, while temperature fields from solving the equation of heat conductivity containing the convection term. Changes in thermo-physical parameters as a function of temperature were considered. The finite element method (FEM) was used to solve the problem.
This study aims to design a novel air cleaning facility which conforms to the current situation in China, and moreover can satisfy our demand on air purification under the condition of poor air quality, as well as discuss the development means of a prototype product. Air conditions in the operating room of a hospital were measured as the research subject of this study. First, a suitable turbulence model and boundary conditions were selected and computational fluid dynamics (CFD) software was used to simulate indoor air distribution. The analysis and comparison of the simulation results suggested that increasing the area of air supply outlets and the number of return air inlets would not only increase the area of unidirectional flow region in main flow region, but also avoid an indoor vortex and turbulivity of the operating area. Based on the summary of heat and humidity management methods, the system operation mode and relevant parameter technologies as well as the characteristics of the thermal-humidity load of the operating room were analyzed and compiled. According to the load value and parameters of indoor design obtained after our calculations, the airflow distribution of purifying the air-conditioning system in a clean operating room was designed and checked. The research results suggested that the application of a secondary return air system in the summer could reduce energy consumption and be consistent with the concept of primary humidity control. This study analyzed the feasibility and energy conservation properties of cleaning air-conditioning technology in operating rooms, proposed some solutions to the problem, and performed a feasible simulation, which provides a reference for practical engineering.
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 aim of this paper is to present the procedure test for calibration and validation of the numerical model for X22CrMoV12-1 steel multilayer welding. On the real multilayer weld was described how to arrange the whole experiment in order to obtain not only relevant input data but also verification data. Tests on a specially prepared specimen, welded with 8 beads in 4 layers, allows to determine the actual geometry of the single welded beads, registration of welding thermal cycles and the hardness distribution in successively deposited beads together with determining the heat influence of subsequent layers. The results of the real welding tests were compared with the results obtained from the numerical simulations and extended by the calculated stresses and distortions distributions of the tested specimen. A new, improved hardness prediction algorithm for high-alloy martensitic and bainitic steels was also proposed.
This study, describing computer simulation of a glider crash against a non-deformable ground barrier, is a part of a larger glider crash modeling project. The studies were intended to develop a numerical model of the pilot - glider - environment system, whereby the dynamics of the human body and the composite cockpit structure during a crash would make it possible to analyze flight accidents with focus on the pilot's safety. Notwithstanding that accidents involving glider crash against a rigid barrier (a wall, for example) are not common, establishing a simulation model for such event may prove quite useful considering subsequent research projects. First, it is much easier to observe the process of composite cockpit structure destruction if the crash is against a rigid barrier. Furthermore, the use of a non-deformable barrier allows one to avoid the errors that are associated with the modeling of a deformable substrate, which in most cases is quite problematic. Crash test simulation, carried out using a MAYMO package, involved a glider crash against a wall positioned perpendicularly to the object moving at a speed of 77 km/h. Computations allowed for determination of time intervals of the signals that are required to assess the behavior of the cockpit and pilot's body - accelerations and displacements in selected points of the glider's structure and loads applied to the pilot's body: head and chest accelerations, forces at femur, lumbar spine and safety belts. Computational results were compared with the results of a previous experimental test that had been designed to verify the numerical model. The glider's cockpit was completely destroyed in the crash and the loads transferred to the pilot's body were very substantial - way over the permitted levels. Since modeling results are fairly consistent with the experimental test, the numerical model can be used for simulation of plane crashes in the future.
A verification study of ultrasound transmission numerical simulation results with experiment results is presented in this paper. The work considers a model of a transformer tank which is filled with electro insulating oil. In the experiment, performed under laboratory conditions, an ultrasound wave is generated by a piezoelectric transducer that is fixed in the centre of the tank and measured by another transducer mounted inside the tank at three distances: 10, 20 and 30 cm from the sound source. The transducer is able to measure and generate acoustic waves in the ultrasound frequency band up to 1 MHz. The simulation considers numerical calculation of acoustic pressure distribution inside the tank in which acoustic source emits waves with frequency equal to 100 kHz. Verification analysis has confirmed consistency of the numerically calculated values with the measurement results.
This paper presents the results of research focused on the lowering of ash flow temperature at semianthracite coal from Donbas district by means of additive (calcite) dosing. Ash fusion temperatures were set for two coal samples (A, B) and for five various states (samples of ash without any additives, with 1%, with 3%, with 5% and with 7% of the additive) in total. The macroscopicphotographic method was used for identifying all specific temperatures. Obtained outputs prove that A type coal has a lower value of sphere temperature than B type coal in the whole scope of percentage representation of the additive. The flow temperature dropped in total from 1489 °C to 1280 °C, i.e. by 14% during the test of coal of type A with 7% of the additive; while it was near 10% for coal of type B (from 1450 °C to 1308 °C). Numerical simulations of the process showed that it is not effective to add an additive with a grain size lower than 280 μm by means of wastevapour burners.
Widely used CFD codes enable modelling of PC boilers operation. One of the areas where these numerical simulations are especially promising is predicting deposition on heat transfer surfaces, mostly superheaters. The basic goal of all simulations is to determine trajectories of ash particles in the vicinity of superheater tubes. It results in finding where on the surface the tube will be hit by particles, and what diameter and mass flow of the particles are. This paper presents results of CFD simulations for a single tube and a bundle of in-line tubes as well. It has been shown that available parameters like ash particle density, shape factor, reflection coefficients affect the trajectories in a different way. All the simulations were carried out with Fluent code of Ansys software.
Snap-fit connections have been used for many years in various fields of technology and everyday objects. They often have complex shapes, which is allowed by the processing technology of the polymers from which they are made, but they are not designed to carry loads. Changing the material to a metal or fiber composite allows these types of joints to be used as replacements for rivets or screws, but there are problems with the closing technique – an increase in closing force due to the large Young’s modulus of these materials relative to polymers without reinforcement. One of the methods to solve this problem may be the use of a thermo-bimetallic effect consisting in heating both or one of the connection parts to the appropriate temperature. This kind of treatment results in deflection of the beam of the clip (Fig. 1), followed by assembly with zero force or less in relation to the case without heating. The paper presents the results of numerical simulations for the connection in which the beam of the clip consisted of two materials: (1) a fiber composite designed to carry loads, (2) thin metal layer tied with the composite and designed to create a thermo-bimetallic effect. In the case of this solution, the main parameter is the difference in coefficients of linear thermal expansion of both materials. The paper presents results for two cases of connection work: closing and opening. The calculations were carried out in the Abaqus/Standard solver using thermal-displacement steps.
The melt cleaning is an important aspect in the production of high-quality aluminum castings. Specifically inclusions within the melt and an excessively high hydrogen content lead to defects and undesired porosity in the castings. Although it is possible to reduce the amount of hydrogen and oxidic inclusions by purge gas treatment and the use of melting salts, it is impossible to remove oxides (bifilms) created during filling of gating system. Paper deals with the effects of melt quality and the placement of a filter in the filling system on Al-7%Si-Mg alloy mechanical properties. Three different filters were used: (a) rectangular ceramic pressed filter with 3 mm thickness (b) cubical pressed ceramic filter with thickness 10 mm (c) cubical pressed ceramic filter with thickness 22 mm. The results showed that the highest tensile strength values were obtained from the filter with thickness of 22 mm. Numerical simulation analysis of the filling process showed that velocity reduction by filter is the major phenomenon affecting the mechanical properties. Another evaluated aspect during experiments was capability of filters to retain old bifilms. For this purpose multiply remelted alloy was prepared and analyzed. Results showed that filter efficiency increases with decreasing melt quality as a result of possibility to retain “old” bifilms better than small and thin “new” bifilms.
Eutectoid growth, as the important reaction mechanism of the carbon steel heat treatment, is the basis to control the microstructure and performance. At present, most studies have focused on lamellar growth, and did not consider the nucleation process. Mainly due to the nucleation theory is inconclusive, a lot of research can support their own theory in a certain range. Based on the existing nucleation theory, this paper proposes a cooperative nucleation model to simulate the nucleation process of eutectoid growth. In order to ensure that the nucleation process is more suitable to the theoretical results, different correction methods were used to amend the model respectively. The results of numerical simulation show that when the model is unmodified, the lateral growth of single phase is faster than that of longitudinal growth, so the morphology is oval. Then, the effects of diffusion correction, mobility correction and ledges nucleation mechanism correction on the morphology of nucleation and the nucleation rate were studied respectively. It was found that the introduction of boundary diffusion and the nucleation mechanism of the ledges could lead to a more realistic pearlite.
Micro-channel heat sinks are used in a wide variety of applications, including microelectronic devices, computers and high-energy-laser mirrors. Due to the high power density that is encountered in these devices (the density of delivered electrical power up to a few kW/cm2) they require efficient cooling as their temperatures must generally not exceed 100 ◦C. In the paper a new design for micro-channel heat sink (MCHS) to be used for cooling laser diode arrays (LDA) is considered. It is made from copper and consisting of 37 micro-channels with length of 9.78 mm, width of 190 μm and depth of 180 μm with the deionized water as a cooling medium. Mathematical and numerical models of the proposed design of the heat sink were developed. A series of thermofluid numerical simulations were performed for various volumetric flow rates of the cooling medium, its inlet temperature and different thermal power released in the laser diode. The results show that the LDA temperature could be decreased from 14 to 17% in comparison with earlier proposed design of the heat sink with the further drop in temperature obtained by applying indium instead of gallium arsenide as the soldering material between the LDA and MCHS interface. Moreover, it was found that the maximum temperature, and therefore the thermal resistance of the considered heat sink, could be decreased by increasing the coolant flow rate.
A new concept of an electrostatic spray column for liquid-liquid extraction was investigated. An important problem for separation processes is the presence of azeotropic or close-boiling mixtures in their production, for example heptane with ethanol, since the separation is impossible by ordinary distillation. The use of ionic liquids (IL) as a dispersed solvent specially engineered for any specific organic mixture in terms of selectivity is a key factor to successful separation. As IL present particularly attractive combination of favorable characteristics for the separation of heptane and ethanol, in this work we use 1-butyl-3-methylimidazolium methyl sulfate [BMIM][MeSO4]. Because of high viscosity and relatively high cost of IL a new technique was introduced, consisting in the electrostatically spray generation to enhance the mass transport between the phases. In order to optimally design the geometry of the contactor a series of numerical simulation was performed. Especially multi-nozzle variants for better exploitation of contactor volume were investigated. Experiments showed excellent possibility of control of the dispersion characteristics by applied voltage and thus control of the rate of extraction. The preliminary simulations based on our mathematical model for a three nozzle variant exhibited visual agreement with the theory of electrostatics.
The aim of the study was the implementation of a numerical simulation of the air-water two-phase flow in the minichannel and comparing results obtained with the values obtained experimentally. To perform the numerical simulations commercial software ANSYS FLUENT 12 was used. The first step of the study was to reproduce the actual research installation as a three-dimensional model with appropriate and possible simplifications - future computational domain. The next step was discretisation of the computational domain and determination of the types of boundary conditions. ANSYS FLUENT 12 has three built-in basic models with which a two-phase flow can be described. However, in this work Volume-of-Fluid (VOF) model was selected as it meets the established requirements of research. Preliminary calculations were performed for a simplified geometry. The calculations were later verified whether or not the simplifications of geometry were chosen correctly and if they affected the calculation. The next stage was validation of the chosen model. After positive verification, a series of calculations was performed, in which the boundary conditions were the same as the starting conditions in laboratory experiments. A satisfactory description of the experimental data accuracy was attained.
Because the heat release of plutonium material, the composite structure is heated and the stress and strain of the composite structure will increase, which will affect the thermodynamic properties of the structure. The thermodynamic analysis of complex structures, which are composed of concentric structures of plutonium, beryllium, tungsten, explosives, and steel, was carried out. The results showed that when the structure is spherical, the temperature is higher than that of the ellipsoid structure. Stress of the elliptical structure is greater than the spherical structure. This study showed that the more flat the shell is, the greater the stress concentration point occurs at the long axis, and the maximum stress occurs at the beryllium layer. These conclusions provide theoretical support for the plutonium composite component testing.
Thermo-chemical treatments are known to increase the fatigue life of industrial parts. Due to the imprecise consideration of residual stresses in predicting the durability of components subjected to cyclic loading and their effect on the fatigue life, the authors developed a numerical model combining the influence of residual stresses with stresses caused by bending. The authors performed the numerical simulation with the use of Finite Element Method to analyse material behaviour during cyclic loading. The residual stress state developed during nitriding was introduced onto cross-section of the numerical specimen. The goal of this work was better understanding of the real conditions of the nitride steel fatigue processes and improving the knowledge about numerical predicting of the fatigue life for parts with residual stresses. The results of simulation were compared with plane bending fatigue tests. The presented method indicates the possibility of increasing the accuracy of the fatigue analysis of elements after surface treatment, increasing its certainty and the ability to perform better optimization of service life.
The aim of this work is to develop a numerical model capable of predicting the grain density in the Mg-based matrix phase of an AZ91/SiC composite, as a function of the total mass fraction of the embedded SiC particles. Based on earlier work in a range of alloy systems, we assume an exponential relationship between the grain density and the maximum supercooling during solidification. Analysis of data from cast samples with different thicknesses, and mass fractions of added SiCp, permits conclusions to be drawn on the role of SiCp in increasing grain density. By fitting the data, an empirical nucleation law is derived that can be used in a micro model. Numerical simulation based on the model can predict the grain density of magnesium alloys containing SiC particles, using the mass fraction of the particles as inputs. These predictions are compared with measured data.
Aluminum profiles play an important role in civil engineering (facades, walls with windows) as well as in mechanical engineering (production lines, constructions of 3D printers and plotters). To ensure quick assembly, disassembly or changed the dimensions of constructions it is not possible to use such methods as welding, adhesive or riveting joints. The solution may be to use the so-called “popular lock”. It is a mechanism, the closure of which is caused by tightening of the conical screw, joining the “T” profile in the node. In order to properly design using the presented type of connection, it is necessary to know its strength and stiffness both in simple and complex loads states, also including imperfections. In the literature there is no information about the operation of the construction node with the so-called “popular lock”. The paper presents the results of experimental tests for connections subjected to uniaxial tensile test, paying special attention to the defects that may appear during the assembly. In the next step, a 3D solid connection model was created. Numerical simulations were performed in the Abaqus / Explicite program for both uniaxial tensile test and bending tests in two planes. Limit values of loads above which there is a plastic deformation of the material were determined. Determination of stiffness and strength of a single node allowed to make a simplified connector model. Using the numerical model, the analysis was performed taking into account the influence of imperfections on the work of the entire connection.
Numerical simulation is an economical and effective method in the field of marine engineering. The dynamics of mooring cables has been analysed by a numerical simulation code that was created on a basis of a new element frame. This paper aims at verifying the accuracy of the numerical simulation code through comparisons with both the real experiments and a commercial simulation code. The real experiments are carried out with a catenary chain mooring in a water tank. The experimental results match the simulation results by the numerical simulation code well. Additionally, a virtual simulation of a large size chain mooring in ocean is carried out by both the numerical simulation code and a commercial simulation code. The simulation results by the numerical simulation code match those by the commercial simulation code well. Thus, the accuracy of the numerical simulation code for underwater chain mooring is verified by both the real experiments and commercial simulation code.
This paper is devoted to the application of ultrasonic wave propagation phenomena for the diagnostics of prestressed, concrete, bridge T-beams. A multi-point damage detection system is studied with use of numerically obtained data. The system is designed to detect the presence of the material discontinuities as well as their location.
While modeling water dynamics in dam reservoirs, it is usually assumed that the flow involves the whole water body. It is true for shallow reservoirs (up to several meters of depth) but may be false for deeper ones. The possible presence of a thermocline creates an inactive bottom layer that does not move, causing all the discharge to be carried by the upper strata. This study compares the results of hydrodynamic simulations performed for the whole reservoir to the ones carried out for the upper strata only. The validity of a non-stratified flow approximation is then discussed.
Lower Carboniferous limestone has been extracted in the “Czatkowice” open-pit hill-slope quarry in southern Poland since 1947, for the needs of metallurgical and building industries, as well as farming. We can distinguish two aquifers in the Czatkowice area: the Quaternary porous aquifer and the Carboniferous fissure-porous one. Two vertical zones representing different hydrodynamic characteristics can be indentified in the Carboniferous formations. One is a weathering zone and the other one the zone of fissures and interbedding planes. Groundwater inflows into the quarry workings have been observed at the lowest mining level (+315 m above the sea level (asl)) for over 30 years. This study concerns two hypotheses of the sources of such inflows originating either from (a) the aeration zone or from (b) the saturation zone. Inflows into the quarry combine into one stream flowing gravitationally to the doline under the pile in the western part of the quarry. This situation does not cause a dewatering need. Extending eastward mining and lowering of the exploitation level lead to increased inflows.
In this paper, the mathematical model and numerical simulations of the molten steel flow by the submerged entry nozzle and the filling process of the continuous casting mould cavity are presented. In the mathematical model, the temperature fields were obtained by solving the energy equation, while the velocity fields were calculated by solving the momentum equations and the continuity equation. These equations contain the turbulent viscosity which is found by solving two additional transport equations for the turbulent kinetic energy and its rate of dissipation. In the numerical simulations, coupling of the thermal and fluid flow phenomena by changes in the thermophysical parameters of alloy depending on the temperature has been taken into consideration. This problem (2D) was solved by using the finite element method. Numerical simulations of filling the continuous casting mould cavity were performed for two variants of liquid metal pouring. The effect of the cases of pouring the continuous casting mould on the velocity fields and the solid phase growth kinetics in the process of filling the continuous casting mould was evaluated as these magnitudes have an influence on the high quality of the continuous cast steel slab.
This paper deals with real-time (RT) simulators applied in power electronic applications and implemented in a real inverter. The process of preparing and starting up an active rectifier prototype (with an active filter function), using the real-time OPAL RT simulator is given. The control system of the converter and the results of simulation using the Matlab/Simulink suite are discussed.