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 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.
This paper proposes a new approach for calculating the impulse response of room acoustics. Impulse response provides unique characterization of any discrete lineartime invariant (LTI) systems. We assume that the room is a linear time-invariant system and the impulse response is calculated simply by sending a Dirac Impulse into the system as input and getting the response from the output. Then, the output of the system is represented as a sum of time-shifted weighted impulse responses. Both mathematical justifications for the proposed method and results from simulation software developed to evaluate the proposed approach are presented in detail.
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.
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.
This paper presents the beam tracing with refraction method, developed to examine the possibility of creating the beam tracing simulation of sound propagation in environments with piecewise non- homogenous media. The beam tracing with refraction method (BTR) is developed as an adaptive beam tracing method that simulates not only the reflection but also the refraction of sound. The scattering and the diffraction of sound are not simulated. The BTR employs 2D and 3D topology in order to efficiently simulate scenes containing non-convex media. After the beam tracing is done all beams are stored in a beam tree and kept in the computer memory. The level of sound intensity at the beginning of each beam is also memorized. This beam data structure enables fast recalculation of results for stationary source and geometry. The BTR was compared with two commercial ray tracing simulations, to check the speed of BTR algorithms. This comparison demonstrated that the BTR has a performance similar to state-of- the-art room-acoustics simulations. To check the ability to simulate refraction, the BTR was compared with a commercial Finite Elements Method (FEM) simulation. In this comparison the BTR simulated the focusing of the ultrasound with an acoustic lens, with good accuracy and excellent performance.
Virtual Reality (VR) systems are used in engineering, architecture, design and in applications of biomedical research. The component of acoustics in such VR systems enables the creation of audio-visual stimuli for applications in room acoustics, building acoustics, automotive acoustics, environmental noise control, machinery noise control, and hearing research. The basis is an appropriate acoustic simulation and auralization technique together with signal processing tools. Auralization is based on time-domain modelling of the components of sound source characterization, sound propagation, and on spatial audio technology. Whether the virtual environment is considered sufficiently accurate or not, depends on many perceptual factors, and on the pre-conditioning and immersion of the user in the virtual environment. In this paper the processing steps for creation of Virtual Acoustic Environments and the achievable degree of realism are briefly reviewed. Applications are discussed in examples of room acoustics, archeological acoustics, aircraft noise, and audiology.
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.
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.
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 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.
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.
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.
n the paper, the simulation PROP5 program with the road model defined as a noise source and road surroundings model, is used to predict the efficiency of noise protecting means for the chosen building. The appropriate models of verified accuracy have been chosen by comparison of the simulation results with field measurements (Walerian et al., 2010). Here, using the pre-tested simulation program, the possibility of acoustic climate improvement has been analyzed in the ranges of practical variations of the input parameters. The road parameters: its geometry (number of lanes and their positions) and traffic structure over lanes (vehicle flow rates and their average speeds) have been taken under consideration as changeable parameters, that could be corrected to obtain acoustical climate improvement. Moreover, an acoustical screen designing has been considered. The screen efficiency has been evaluated under conditions defined by the input parameters of the road and its surroundings.
In the paper, the simulation PROP5 program is used to predict the sound level in proximity of a road with defined surroundings. The simulation involves road geometry (number of lanes and their positions) and traffic structure (vehicle flow rates and their average speeds), with equivalent omnidirectional point sources representing vehicles. In Part I of the paper, the agreement between measurement and simulation results is tested to verify the accuracy degree of the applied models of a road, as a noise source and propagation throughout surrounding space. In Part II, using the pre-tested simulation program, the possibility of acoustic climate improvement has been analyzed.
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.
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.
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.