The aim of the presented paper is to show the results of shape optimization of railway polynomial transition curves (TCs) of 5th, 7th, and 9th degrees through the use of the full vehicle model and new criteria of assessement concerning the jerk value. The search for the proper shape of TCs means that in this work, the evaluation of TC properties is based on select quantities and the generation of such a shape through the use of mathematically understood optimization methods. The studies presented have got a character of the numerical tests. For this work, advanced vehicle models describing dynamical track-vehicle and vehicle-passenger interactions as well as optimization methods were exploited. In the software vehicle model of a 2-axle freight car, the track discrete model, non-linear descriptions of wheel-rail contact are applied. This part of the software, the vehicle simulation software, is combined with a library optimization procedure into the final computer program.
The paper presents the application of the newly developed method of the solution of nonlinear equations to the adaptive modelling and computer simulation. The approach is suitable when the system of equations can be arranged in such a way that it consists of a large number of linear equations and a smaller number of nonlinear equations. This situation occurs in the case of adaptive modelling of mechanical systems using finite elements or finite differences techniques. In this case the classical least square method becomes very effective. The paper presents several examples of the application of the method. A solution to the, so called, “black box” problem is also presented.
In this paper starch gelatinisation in Couette-Taylor flow (CTF) apparatus (equipped with a water heat jacket) has been investigated. CTF (characterised by the presence of Taylor vortices) provides good environment for gelatinisation, e.g. effective mixing, fast heat transfer, positive influence on starch rheological properties. During experiments starch gelatinisation degree and starch swelling has been studied. It was accompanied by temperature measurements performed along the apparatus. Additionally, starch gelatinisation was investigated by computer simulation. A complete starch gelatinisation was obtained for the shortest investigated residence time in the apparatus when the temperature in the heat jacket was above 85 °C. Nevertheless, it seems that it is still possible to reduce a residence time value, as well as, the value of Thj, but it may require some acceleration of rotor rotation. The swelling degree of gelatinised starch increased with growing values of residence time, rotor rotation and process temperature. Heat transferred could be affected by the structure of the Taylor vortex flow. No significant destruction of starch granules was observed during the treatment in Couette-Taylor flow apparatus. A quite satisfactory agreement between computer simulation and experiments results was achieved.
Steel is a versatile material that has found widespread use because of its mechanical properties, its relatively low cost, and the ease with which it can be used in manufacturing process such as forming, welding and machining. Regarding to mechanical properties are strongly affected by grain size and chemical composition variations. Many industrial developments have been carried out both from the point of view of composition variation and grain size in order to exploit the effect of these variables to improve the mechanical proprieties of steels. It is also evident that grain growth are relevant to the mechanical properties of steels, thus suggesting the necessity of mathematical models able to predict the microstructural evolution after thermo cycles. It is therefore of primary importance to study microstructural changes, such as grain size variations of steels during isothermal treatments through the application of a mathematical model, able in general to describe the grain growth in metals. This paper deals with the grain growth modelling of steels based on the statistical theory of grain growth originally developed by Lücke  and here integrated to take into account the Zener drag effect and is therefore focused on the process description for the determination of the kinetics of grain growth curves temperature dependence.
The main aim of the presented research was to check mechanical response of human body model under loads that can occur during airplane accidents and compare results of analysis with some results of experimental tests described in literature. In simulations, new multi-purpose human body model, the VIRTHUMAN, was used. The whole model, as well as its particular segments, was earlier validated based on experimental data, which proved its accuracy to simulate human body dynamic response under condition typical for car crashes, but it was not validated for loads with predominant vertical component (loads acting along spinal column), typical for airplane crashes. Due to limitation of available experimental data, the authors focused on conducting calculations for the case introduced in 14 CFR: Parts 23.562 and 25.562, paragraph (b)(1), knowing as the 60 pitch test. The analysis consists in comparison of compression load measured in lumbar section of spine of the FAA HIII Dummy (experimental model) and in the Virthuman (numerical model). The performed analyses show numerical stability of the model and satisfactory agreement between experimental data and simulated Virthuman responses. In that sense, the Virthuman model, although originally developed for automotive analyses, shows also great potential to become valuable tool for applications in aviation crashworthiness and safety analyses, as well.
There are certain well-known methods of diminishing concentrations of nitrogen compounds, but they are ineffective in case of nitrogen-rich wastewater with a low content of biodegradable carbon. Partial nitritation followed by anaerobic ammonium oxidation (Anammox) process appear to be an excellent alternative for traditional nitrification and denitrification. This paper presents the feasibility of successful start-up of Anammox process in a laboratory-scale membrane bioreactor (MBR). It was shown that the combination of membrane technology and Anammox process allowed to create a new highly efficient and compact system for nitrogen removal. It was possible to achieve average nitrogen removal efficiency equal to 76.7 ± 8.3%. It was shown that the start-up period of 6 months was needed to obtain high nitrogen removal efficiency. The applied biochemical model of the Anammox process was based on the state-of-the-art Activated Sludge Model No.1 (ASM 1) which was modified for accounting activity of autotrophs (nitrite-oxidising bacteria and nitrateoxidising bacteria) and anammox bacteria. In order to increase the predictive power of the simulation selected parameters of the model were adjusted during model calibration. Readjustment of the model parameters based on the critically evaluated data of the reactor resulted in a satisfactory match between the model predictions and the actual observations.
The purpose of this study is to identify relationships between the values of the fluidity obtained by computer simulation and by an experimental test in the horizontal three-channel mould designed in accordance with the Measurement Systems Analysis. Al-Si alloy was a model material. The factors affecting the fluidity varied in following ranges: Si content 5 wt.% – 12 wt.%, Fe content 0.15 wt.% – 0.3wt. %, the pouring temperature 605°C-830°C, and the pouring speed 100 g · s–1 – 400 g · s–1. The software NovaFlow&Solid was used for simulations. The statistically significant difference between the value of fluidity calculated by the equation and obtained by experiment was not found. This design simplifies the calculation of the capability of the measurement process of the fluidity with full replacement of experiments by calculation, using regression equation.
Metallographic investigations and a computer simulation of stresses in a gravity die-casting bushing were performed. Simulation of the casting process, solidification of the thick-walled bushing and calculations of the stress was performed using MAGMA5.3 software. The size variability of phases κIIaffecting the formation of phase stresses σf, depending on the location of the metallographic test area, was identified. The distribution of thermal σtand shrinkage stresses σs, depending on the location of the control point SC in the bushing's volume, was estimated. Probably the nature of these stresses will change slightly even after machining. This can cause variations in operating characteristics (friction coefficient, wear). Due to the strong inhomogeneity of the stress distribution in the bushing's casting, it is necessary to perform further tests of the possibility to conduct thermal treatment guaranteeing homogenization of the internal stresses in the casting, as well as to introduce changes in the bushing's construction and the casting technology. The paper presents the continuation of the results of research aimed at identifying the causes of defects in the thick-walled bushing, die-casting made of CuAl10Fe5Ni5Cr aluminium bronze.
The article is a case study of the steel milling ring casting of about 6 tonnes net weight. The casting has been cast in the steel foundry the authors have been cooperating with. The aim was to analyse the influence of the shape of the chills and the material which was used to make them on the casting crystallization process. To optimally design the chills the set of the computer simulation has been carried out with 3 chills’ shape versions and 3 material’s versions and the results have been compared with the technology being in use (no chills). The proposed chills were of different thermal conductivity from low to high. Their shapes were obviously dependant on the adjacent casting surface geometry but were the result of the attempt to optimise their effect with the minimum weight, too. The chills working efficiency was analysed jointly with the previously designed top feeders system. The following parameters have been chosen to compare their effectiveness and the crystallization process: time to complete solidification and so-called fed volume describing the casting feeding efficiency. The computer simulations have been carried out with use of MagmaSoft v. 5.2 software. Finally, the optimisation has led to 15% better steel yield thanks to 60% top feeders weight reduction and 40% shorter solidification time. The steel ring cast with use of such technology fulfil all quality criteria.
Experiments of filling the model moulds cavity of various inner shapes inserted in rectangular cavity of the casting die (dimensions: 280 mm (height) x 190 mm (width) x 10 mm (depth) by applying model liquids of various density and viscosity are presented in the paper. Influence of die venting as well as inlet system area and inlet velocity on the volumetric rate of filling of the model liquid – achieved by means of filming the process in the system of a cold-chamber casting die was tested. Experiments compared with the results of simulation performed by means of the calculation module Novacast (Novaflow&Solid) for the selected various casting conditions – are also presented in the paper.
Results of 3-D of discharge channel displacement simulation, acquired by means of the Fluent program during one current half-period of AC arc, indicate that the obtained images of the phenomenon are qualitatively similar to those, recorded with a high-speed digital camera, while the computer simulation enables much a more comprehensive analysis of the acquired data. In addition to selected arc simulation frames and corresponding distributions of mass velocity vectors and current density vectors on a plane, the distributions of temperature, current density and mass velocity values are presented on the axis of the electrode arrangement model. The composite motion (continuous and jumping) of discharge channels was analyzed, taking into account mass displacement and matter state changes.
The paper presents a multi-scale mathematical model dedicated to a comprehensive simulation of resistance heating combined with the melting and controlled cooling of steel samples. Experiments in order to verify the formulated numerical model were performed using a Gleeble 3800 thermo-mechanical simulator. The model for the macro scale was based upon the solution of Fourier-Kirchhoff equation as regards predicting the distribution of temperature fields within the volume of the sample. The macro scale solution is complemented by a functional model generating voluminal heat sources, resulting from the electric current flowing through the sample. The model for the micro-scale, concerning the grain growth simulation, is based upon the probabilistic Monte Carlo algorithm, and on the minimization of the system energy. The model takes into account the forming mushy zone, where grains degrade at the melting stage – it is a unique feature of the micro-solution. The solution domains are coupled by the interpolation of node temperatures of the finite element mesh (the macro model) onto the Monte Carlo cells (micro model). The paper is complemented with examples of resistance heating results and macro- and micro-structural tests, along with test computations concerning the estimation of the range of zones with diverse dynamics of grain growth.
The quantitative description of an airlift bioreactor, in which aerobic biodegradation limited by carbonaceous substrate and oxygen dissolved in a liquid takes place, is presented. This process is described by the double-substrate kinetics. Mathematical models based on the assumption of plug flow and dispersion flow of liquid through the riser and the downcomer in the reactor were proposed. Calculations were performed for two representative hydrodynamic regimes of reactor operation, i.e. with the presence of gas bubbles only within the riser and for complete gas circulation. The analysis aimed at how the choice of a mathematical model of the process would enable detecting the theoretical occurrence of oxygen deficiency in the airlift reactor. It was demonstrated that the simplification of numerical calculations by assuming the “plug flow” model instead of dispersion with high Péclet numbers posed a risk of improper evaluation of the presence of oxygen deficiency zones. Conclusions related to apparatusmodelling and process design were drawn on the basis of the results obtained. The paper is a continuation of an earlier publication (Grzywacz, 2012a) where an analysis of single-substrate models of the airlift reactor was presented.
CFD modelling of momentum and heat transfer using the Large Eddy Simulation (LES) approach has been presented for a Kenics static mixer. The simulations were performed with the commercial code ANSYS Fluent 15 for turbulent flow of three values of Reynolds number, Re = 5 000, 10 000 and 18 000. The numerical modelling began in the RANS model, where standard k−ε turbulence model and wall functions were used. Then the LES iterations started from the initial velocity and temperature fields obtained in RANS. In LES, the Smagorinsky–Lilly model was used for the sub-grid scale fluctuations along with wall functions for prediction of flow and heat transfer in the near-wall region. The performed numerical study in a Kenics static mixer resulted in highly fluctuating fields of both velocity and temperature. Simulation results were presented and analysed in the form of velocity and temperature contours. In addition, the surface-averaged heat transfer coefficient values for the whole insert length were computed and compared with the literature experimental data. Good compliance of the LES simulation results with the experimental correlation was obtained.
Rock excavation is a basic technological operation during tunnelling and drilling roadways in underground mines. Tunnels and roadways in underground mines are driven into a rock mass, which in the particular case of sedimentary rocks, often have a layered structure and complicated tectonics. For this reason, rock strata often have highly differentiated mechanical properties, diverse deposition patterns and varied thicknesses in the cross sections of such headings. In the field of roadheader technology applied to drilling headings, the structure of a rock mass is highly relevant when selecting the appropriate cutting method for the heading face. Decidedly differentiated values of the parameters which describe the mechanical properties of a particular rock layer deposited in the cross section of the drilled tunnel heading will influence the value and character of the load on the cutting system, generated by the cutting process, power demand, efficiency and energy consumption of the cutting process. The article presents a mathematical modelling process for cutting a layered structure rock mass with the transverse head of a boom-type roadheader. The assumption was made that the rock mass being cut consists of a certain number of rock layers with predefined mechanical properties, a specific thickness and deposition pattern. The mathematical model created was executed through a computer programme. It was used for analysing the impact deposition patterns of rock layers with varied mechanical properties, have on the amount of cutting power consumed and load placed on a roadheader cutting system. The article presents an example of the results attained from computer simulations. They indicate that variations in the properties of the rock cut – as cutting heads are moving along the surface of the heading face – may have, apart from multiple other factors, a significant impact on the value of the power consumed by the cutting process.
This article considers designing of a renewable electrical power generation system for self-contained homes away from conventional grids. A model based on a technique for the analysis and evaluation of two solar and wind energy sources, electrochemical storage and charging of a housing area is introduced into a simulation and calculation program that aims to decide, based on the optimized results, on electrical energy production system coupled or separated from the two sources mentioned above that must be able to ensure a continuous energy balance at any time of the day. Such system is the most cost-effective among the systems found. The wind system adopted in the study is of the low starting speed that meets the criteria of low winds in the selected region under study unlike the adequate solar resource, which will lead to an examination of its feasibility and profitability to compensate for the inactivity of photovoltaic panels in periods of no sunlight. That is a system with fewer photovoltaic panels and storage batteries whereby these should return a full day of autonomy. Two configurations are selected and discussed. The first is composed of photovoltaic panels and storage batteries and the other includes the addition of a wind system in combination with the photovoltaic system with storage but at a higher investment cost than the first. Consequently, this result proves that is preferable to opt for a purely photovoltaic system supported by the storage in this type of site and invalidates the interest of adding micro wind turbines adapted to sites with low wind resources.
According to the European Environment Agency (EEA 2018), air quality in Poland is one of the worst in Europe. There are several sources of air pollution, but the condition of the air in Poland is primarily the result of the so-called low-stack emissions from the household sector. The main reason for the emission of pollutants is the combustion of low-quality fuels (mainly low-quality coal) and waste, and the use of obsolete heating boilers with low efficiency and without appropriate filters. The aim of the study was to evaluate the impact of measures aimed at reducing low-stack emissions from the household sector (boiler replacement, change of fuel type, and thermal insulation of buildings), resulting from environmental regulations, on the improvement of energy efficiency and the emission of pollutants from the household sector in Poland. Stochastic energy and mass balance models for a hypothetical household, which were used to assess the impact of remedial actions on the energy efficiency and emission of pollutants, have been developed. The annual energy consumption and emissions of pollutants were estimated for hypothetical households before and after the implementation of a given remedial action. The calculations, using the Monte Carlo simulation, were carried out for several thousand hypothetical households, for which the values of the technical parameters (type of residential building, residential building area, unitary energy demand for heating, type of heat source) were randomly drawn from probability distributions developed on the basis of the analysis of the domestic structure of households. The model takes the coefficients of correlation between the explanatory variables in the model into account. The obtained results were multiplied so that the number of hypothetical households was equal to 14.1 million, i.e. the real number of households in Poland. The obtained results allowed for identifying the potential for reducing the emission of pollutants such as carbon dioxide, carbon monoxide, dust, and nitrogen oxides, and improving the energy efficiency as a result of the proposed and implemented measures, aimed at reducing low-stack emission, resulting from the policy. The potential for emissions of gaseous pollutants is 94% for CO, 49% for NOx, 90% for dust, and 87% for SO2. The potential for improving the energy efficiency in households is around 42%.
This work presents the methodology for analyzing the impact of ground vibrations induced during the drilling of gas/oil exploration wells on the surrounding constructions, as well as on humans and the natural environment. In the primary stage, this methodology is based on measurements of ground vibrations induced by a specific type of drilling system in the so-called reference site. In the next stage, ground vibrations are estimated in similar conditions to another design site, these conditions are assumed for a given drilling system, treated as a vibration source. In both sites, special seismic and geotechnical data are collected to construct numerical models for dynamic analyses. Finally, if it is required, a protection system is proposed with respect to the drilling technology and local conditions. The methodology presented has been tested on the terrain of an active natural gas mine used as the design site, and located in the southeastern part of Poland. The reference site was placed in the terrain of a working drilling system in similar conditions in the central part of Poland. Based on the results of numerical simulations, one may verify the different locations of the drilling rig in the design site with respect to the existing industrial structure. Due to the hazard from destructive ground vibrations, a certain vibroisolation system was proposed at the design site. Based on the results of numerical simulations one could rearrange the components of the drilling system in order to provide maximum security for the surrounding structures.
The engine simulations have become an integral part of engine design and development. They are based on approximations and assumptions. The precision of the results depends on the accuracy of these hypotheses. The simplified models of frozen composition, chemical equilibrium and chemical kinetics provide the compositions of combustion products for engine cycle simulations. This paper evaluates the effects of different operating conditions and hypotheses on the exergetic analysis of a spark-ignition engine. The Brazilian automotive market has the highest number of flex-fuel vehicles. Therefore, a flex-fuel engine is considered for simulations in order to demonstrate the effects of these different hypotheses. The stroke length and bore diameter have the same value of 80 mm. The in-cylinder irreversibility is calculated for each case at the closed part of the engine cycle. A comparative analysis of these hypotheses provides a comprehensive evaluation of their effects on exergetic analysis. Higher values of accumulated irreversibility are observed for the oversimplified hypothesis.
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.
By simulating the actual working conditions of a cable, the temperature variation rule of different measuring points under different load currents was analyzed. On this basis, a three-dimensional finite element model (FEM) was established, and the difference and influence factors between the simulation temperature and the experimental measured value were discussed, then the influence of thermal conductivity on the operating temperature of the conductor layer was studied. Finally, combined with the steady-state thermal conductivity model and the experimental measured data, the relation between thermal conductivity and load current was obtained.
The paper presents the evaluation of utility properties of a tube bend produced by bending process using local induction heating technique. Optimal process parameters were defined on the basis of numerical simulations. Mechanical testing procedures for the tube bend were carried out after normalising annealing. Heat treatment parameters have been chosen in compliance with the relevant standard. Geometrical measurements of the bend manufactured under industrial conditions indicate high accuracy of numerical simulations. Geometry and mechanical properties of the produced bend were consistent with the requirements of the applicable standards.
The paper presents an approach of numerical modelling of alloy solidification in permanent mold and transient heat transport between the casting and the mold in two-dimensional space. The gap of time-dependent width called "air gap", filled with heat conducting gaseous medium is included in the model. The coefficient of thermal conductivity of the gas filling the space between the casting and the mold is small enough to introduce significant thermal resistance into the heat transport process. The mathematical model of heat transport is based on the partial differential equation of heat conduction written independently for the solidifying region and the mold. Appropriate solidification model based on the latent heat of solidification is also included in the mathematical description. These equations are supplemented by appropriate initial and boundary conditions. The formation process of air gap depends on the thermal deformations of the mold and the casting. The numerical model is based on the finite element method (FEM) with independent spatial discretization of interacting regions. It results in multi-mesh problem because the considered regions are disconnected.
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.
Dynamic Mine disasters can be induced by the instability and failure of a composite structure of rock and coal layers during coal mining. Coal seam contains many native defects, severely affecting the instability and failure of the compound structure. In this study, the effects of coal persistent joint on the strength and failure characteristics of coal-rock composite samples were evaluated using PFC2D software. The results show that with the increase of included angle α between the loading direction and joint plane direction, the uniaxial compressive stress (UCS) and peak strain of composite samples first decrease and then gradually increase. The elastic moduli of composite samples do not change obviously with α. The peak strain at α of 45° is the lowest, and the UCS at α of 30° is the smallest. This is inconsistent with theoretical analysis of lowest UCS at α of 45°. This is because that the local stress concentration caused by the motion inconformity of composite samples may increase the average axial stress of upper wall in PFC2D software. Moreover, the coal persistent joint promotes the transformation from the unstable crack expansion to the macro-instability of composite samples, especially at α of 30° and 45°. The majority of failures for composite samples occur within the coal, and no obvious damage is observed in rock. Their failure modes are shear failure crossing or along the coal persistent joint. The failure of composite sample at α of 30° is a mixed failure, including the shear failure along the persistent joint in coal and tensile failure of rock induced by the propagation of coal persistent joint.