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.
A novel type of an axial, piston-driven high pressure hydraulic pump with variable capacity marks a significant improvement in the area of the hydraulic machinery design. Total discharge from hydrostatic forces eliminates a need for a servomechanism, thus simplifying operation, reducing weight and introducing the possibility of the pump displacement control by computer. PWK-type pumps, invented in the Gdansk University of Technology, offer high efficiency for pressure levels up to 55 MPa, ability to work self sucking even at high speed. However, the heart of the new invention, the commutation unit, creates harmful pressure peaks. Those peaks can be mitigated by the introduction of a compensation chamber with elastic walls. Owing to the dynamic character of events taking place in the pump, a need for computer simulation arouse in order to understand phenomena leading to the occurrence of pressure peaks and choose compensation chamber parameters accordingly. A CFD package alone would not be sufficient to reliably represent the interaction between the compensation chamber wall and the working fluid. This paper presents Fluid Structure Interaction approach comparing 3 different models: 2 simplified models of the pump and a full pump model.
The article presents an analysis of stresses in the current tool system of the die during the implementation of the third forging operation of the screw M12 class 10.9 with cylinder head and hexagonal socket. It was assumed that the level of negative cracking due to stress can be reduced by using a mounting interference between the die and the tube blank. Due to the design of the tool system value of the die, the interference value cannot be too large. Therefore, an analysis of the influence of the interference between the die and the tube blank in a die tool system on the value and distribution of stresses in the individual components. An analysis of the assembly stresses and the stresses occurring during the process of deformation of the shaped head of the screw was done. The calculations were performed using a commercial software package MARC / Mentat.
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.
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.
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 article presents the use of computer graphics methods and computational geometry for the analysis on changes of geometrical parameters for a mixed zone in resistance-heated samples. To perform the physical simulation series of resistance heating process, the Gleeble 3800 physical simulator, located in the Institute for Ferrous Metallurgy in Gliwice, was used. The paper presents a description of the test stand and the method for performing the experiment. The numerical model is based on the Fourier-Kirchoff differential equation for unsteady heat flow with an internal volumetric heat source. In the case of direct heating of the sample, geometrical parameters of the remelting zone change rapidly. The described methodology of using shape descriptors to characterise the studied zone during the process allows to parametrise the heat influence zones. The shape descriptors were used for the chosen for characteristic timing steps of the simulation, which allowed the authors to describe the changes of the studied parameters as a function of temperature. Additionally, to determine the impact of external factors, the remelting zone parameters were estimated for two types of grips holding the sample, so-called hot grips of a shorter contact area with the sample, and so-called cold grips. Based on the collected data, conclusions were drawn on the impact of the process parameters on the localisation and shape of the mushy zone.
Scientific research discussed in the present article is focused on the determination of the vertical conveyance capacity in the process of mining minerals, while applying a mathematical calculation and verification of the calculation results by simulation. Input parameters for the capacity calculation include the transport cycle time. The article presents the results of measuring a transport cycle during the operation and a calculation of the transport cycle while using known formulas. On the basis of the observed findings, two methods of increasing the hoisting machine capacity were proposed. The first method is increasing the velocity from the original value of 6 m.s–1 to the velocity of 7 m.s–1. In this case, we achieved the daily capacity increase in 2-2.5%. The second method consisted in changing the hoisting machine acceleration and deceleration modes by which we achieved as much as 9% increase in the daily capacity. The article also describes a transport cycle simulation model, with its output being the number of work cycles parameter. The obtained parameter was used again in the capacity calculation. The simulation model was used in experiments for both, the current status as well as proposed solutions. The simulation model serves also for calculation verification.
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.
Simplifications used in simulation program codes require the use of substitute parameters in the material databases (also called apparent or substitutive). On the one hand, they formally fit into the records used in the heat flow model, porosity, properties etc. and on the other hand they should be determined in conditions most similar to the real casting-mould system. The article presents results of a research on thermophysical parameters of gypsum mould used for precision casting moulds. Experiments were carried out on a cylindrical mould made of Plasticast gypsum, in which the heat source was a cylinder filled with liquid aluminium alloy of the temperature of 655°C. Energetic validation was carried out by using the NovaFlow&Solid ver. 6.3 simulation code. As a result of validation tests, substitute thermophysical parameters of gypsum were determined. For determined parameters, best-fit of solidification time from the experiment and simulation was obtained and the curves of gypsum mass heating were satisfactorily recreate.