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 irregularity profiles of steel samples after vapour blasting were measured. A correlation analysis of profile parameters was then carried out. As the result, the following parameters were selected: Pq, Pt, PDq, Pp/Pt and Pku. Surface profiles after vapour blasting were modeled. The modeled surfaces were correctly matched to measured surfaces in 78% of all analyzed cases. The vapour blasting experiment was then carried out using an orthogonal selective research plan. The distance between the nozzle and sample d and the pressure of feed system p were input parameters; selected surface texture coefficients were output parameters. As the result of the experiment, regression equations connecting vapour blasting process parameters p and d with selected profile parameters were obtained. Finally, 2D profiles of steel samples were forecasted for various values of vapour blasting parameters. Proper matching accuracy of modeled to measured profiles was assured in 75% of analyzed cases.
A hybrid artificial boundary condition (HABC) that combines the volume-based acoustic damping layer (ADL) and the local face-based characteristic boundary condition (CBC) is presented to enhance the absorption of acoustic waves near the computational boundaries. This method is applied to the prediction of aerodynamic noise from a circular cylinder immersed in uniform compressible viscous flow. Different ADLs are designed to assess their effectiveness whereby the effect of the mesh-stretch direction on wave absorption in the ADL is analysed. Large eddy simulation (LES) and FW-H acoustic analogy method are implemented to predict the far-field noise, and the sensitivities of each approach to the HABC are compared. In the LES computed propagation field of the fluctuation pressure and the frequency-domain results, the spurious reflections at edges are found to be significantly eliminated by the HABC through the effective dissipation of incident waves along the wave-front direction in the ADL. Thereby, the LES results are found to be in a good agreement with the acoustic pressure predicted using FW-H method, which is observed to be just affected slightly by reflected waves.
Technological development offers a wide range of new possibilities for implementation of production processes. Continual production development is the main key to success and competitiveness improvement, labour productivity and image-building for all manufacturing companies. The article deals with designing of new workplace with implementation and utilization of automated robot for faster and safer handling of cast stock. The new layout of workplace is created in software Process Simulate.
The paper is focused on properties testing of materials used in form of iso-exo sleeves for risers in ferrous alloys foundry. They are grainyfibrous materials, containing components which initiate and upkeep exothermic reaction. Thermo-physical parameters characterizing such sleeves are necessary also to fill in reliable databases for computer simulation of processes in the casting-mould layout. Studies with use of a liquid alloy, especially regarding different sleeves bring valuable results, but are also relatively expensive and require longer test preparation time. A simplified method of study in laboratory conditions was proposed, in a furnace heated to a temperature above ignition temperature of sleeve material (initiation of exothermic reaction). This method allows to determine the basic parameters of each new sleeve supplied to foundries and assures relatively quick evaluation of sleeve quality, by comparison with previous sleeve supplies or with sleeves brought by new providers.
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.
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 problem of the design of a perfect reduced-order unknown-input observer for standard systems is formulated and solved. The procedure of designing the observer using well-known canonical form is proposed and illustrated with a numerical example. Necessary and sufficient conditions for the solvability of the procedure are given.
The goal of the project is to investigate the influence of elastic mechanisms on technical, bipedal locomotion. In particular, the paper presents the parameter identification for a biologically inspired two-legged robot model. The simulation model consists of a rigid body model equipped with rubber straps. The arrangement of the rubber straps is based on the arrangement of certain muscle groups in a human being. The parameters of the elastic elements are identified applying numerical optimisation. Thus two optimisation algorithms are investigated and compared with respect to robustness and computing time. Moreover, different objective functions are defined and discussed. The behaviour of the resulting configuration of the system is explored in terms of biomechanics.
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.
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.
Measurements of the absorption rate of carbon dioxide into aqueous solutions of N-methyldiethanoloamine (MDEA) and 2-ethylaminoethanol (EAE) have been carried out. On this basis a mathematical model of the performance of an absorption column operated with aqueous solution of a blend of the above amines at elevated temperatures and pressures have been proposed. The results of simulations obtained by means of this model are described. The work is a part of a wider program, aimed at the development of a new process.
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.
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.
In order to study the effects of various gating systems on the casting of a complex aluminum alloyed multi-way valve body, both software simulation analysis and optimization were carried out. Following, the aluminum alloyed multi-way valve body was cast to check the pouring of the aluminum alloy valve body. The computer simulation results demonstrated that compared to the single side casting mode, the casting method of both sides of the gating system would reduce the filling of the external gas, while the air contact time would be lower. Adversely, due to the pouring on both sides, the melt cannot reach at the same time, leading to the liquid metal speed into the cavity to differ, which affected the liquid metal filling stability. The riser unreasonable setting led to the solidification time extension, resulting in a high amount of casting defects during solidification. Also, both gating systems led the entire casting inconsequential solidification. To overcome the latter problems, a straight gate was set at the middle pouring and the horizontal gate diversion occurred on both sides of pouring, which could provide better casting results for the aluminum alloyed multi-valve body.
We propose the time slot routing, a novel routing scheme that allows for a simple design of interconnection networks. The simulative results show that the proposed scheme demonstrates optimal performance at the maximal uniform network load, and for uniform loads the network throughput is greater than for deflection routing.
Foundry technologists use their own style of gating system designing. Most of their patterns are caused by experience. The designs differ from plant to plant and give better or worse results. This shows that the theory of gating systems is not brought into general use sufficiently and therefore not applied in practise very often. Hence, this paper describes the theory and practical development of one part of gating systems - sprue base for automated horizontal moulding lines used for iron castings. Different geometries of sprue bases with gating system and casting were drawn in Solid Edge ST9. The metal flow through the gating systems was then simulated with use of MAGMA Express 188.8.131.52, and the results were achieved. The quality of flow was considered in a few categories: splashes, air entrapment, vortex generation and air contact. The economical aspect (weight of runner) was also taken under consideration. After quantitative evaluation, the best shape was chosen and optimised in other simulations with special attention on its impact on filling velocity and mould erosion. This design (a sprue base with notch placed in drag and cope) is recommended to be used in mass production iron foundries to reduce oxide creation in liquid metal and especially to still metal stream to improve filtration.
Aluminum 6082-T6 panels were joined by friction stir welding utilizing a bobbin tool. A thermal simulation of the process was developed based upon machine torque and the temperature dependent yield stress utilizing a slip factor and an assumed coefficient of friction. The torque-based approach was compared to another simulation established on the shear layer methodology (SLM), which does not require the slip factor or coefficient of friction as model inputs. The SLM simulation, however, only models heat generation from the leading edges of the tool. Ultimately, the two approaches yielded matching temperature predictions as both methodologies predicted the same overall total heat generation from the tool. A modified shear layer approach is proposed that adopts the flexibility and convenience of the shear layer method, yet models heat generation from all tool/workpiece interfaces.