Methods and Models for Accelerating Dynamic Simulation of Fluid Power Circuits

The objective of this dissertation is to improve the dynamic simulation of fluid power circuits. A fluid power circuit is a typical way to implement power transmission in mobile working machines, e.g. cranes, excavators etc. Dynamic simulation is an essential tool in developing controllability and energy-efficient solutions for mobile machines. Efficient dynamic simulation is the basic requirement for the real-time simulation. In the real-time simulation of fluid power circuits there exist numerical problems due to the software and methods used for modelling and integration. A simulation model of a fluid power circuit is typically created using differential and algebraic equations. Efficient numerical methods are required since differential equations must be solved in real time. Unfortunately, simulation software packages offer only a limited selection of numerical solvers. Numerical problems cause noise to the results, which in many cases leads the simulation run to fail. Mathematically the fluid power circuit models are stiff systems of ordinary differential equations. Numerical solution of the stiff systems can be improved by two alternative approaches. The first is to develop numerical solvers suitable for solving stiff systems. The second is to decrease the model stiffness itself by introducing models and algorithms that either decrease the highest eigenvalues or neglect them by introducing steady-state solutions of the stiff parts of the models. The thesis proposes novel methods using the latter approach. The study aims to develop practical methods usable in dynamic simulation of fluid power circuits using explicit fixed-step integration algorithms. In this thesis, twomechanisms whichmake the systemstiff are studied. These are the pressure drop approaching zero in the turbulent orifice model and the volume approaching zero in the equation of pressure build-up. These are the critical areas to which alternative methods for modelling and numerical simulation are proposed. Generally, in hydraulic power transmission systems the orifice flow is clearly in the turbulent area. The flow becomes laminar as the pressure drop over the orifice approaches zero only in rare situations. These are e.g. when a valve is closed, or an actuator is driven against an end stopper, or external force makes actuator to switch its direction during operation. This means that in terms of accuracy, the description of laminar flow is not necessary. But, unfortunately, when a purely turbulent description of the orifice is used, numerical problems occur when the pressure drop comes close to zero since the first derivative of flow with respect to the pressure drop approaches infinity when the pressure drop approaches zero. Furthermore, the second derivative becomes discontinuous, which causes numerical noise and an infinitely small integration step when a variable step integrator is used. A numerically efficient model for the orifice flow is proposed using a cubic spline function to describe the flow in the laminar and transition areas. Parameters for the cubic spline function are selected such that its first derivative is equal to the first derivative of the pure turbulent orifice flow model in the boundary condition. In the dynamic simulation of fluid power circuits, a tradeoff exists between accuracy and calculation speed. This investigation is made for the two-regime flow orifice model. Especially inside of many types of valves, as well as between them, there exist very small volumes. The integration of pressures in small fluid volumes causes numerical problems in fluid power circuit simulation. Particularly in realtime simulation, these numerical problems are a great weakness. The system stiffness approaches infinity as the fluid volume approaches zero. If fixed step explicit algorithms for solving ordinary differential equations (ODE) are used, the system stability would easily be lost when integrating pressures in small volumes. To solve the problem caused by small fluid volumes, a pseudo-dynamic solver is proposed. Instead of integration of the pressure in a small volume, the pressure is solved as a steady-state pressure created in a separate cascade loop by numerical integration. The hydraulic capacitance V/Be of the parts of the circuit whose pressures are solved by the pseudo-dynamic method should be orders of magnitude smaller than that of those partswhose pressures are integrated. The key advantage of this novel method is that the numerical problems caused by the small volumes are completely avoided. Also, the method is freely applicable regardless of the integration routine applied. The superiority of both above-mentioned methods is that they are suited for use together with the semi-empirical modelling method which necessarily does not require any geometrical data of the valves and actuators to be modelled. In this modelling method, most of the needed component information can be taken from the manufacturer’s nominal graphs. This thesis introduces the methods and shows several numerical examples to demonstrate how the proposed methods improve the dynamic simulation of various hydraulic circuits.

Impact of regeneration method on stand structure prior to first thinning

Summary

ADVANCED NUMERICAL METHODS FOR SIMULATING MICROBUBBLE IN AN AERATED MIXING TANK

Fluid mixing in mechanically agitated tanks is one of the major unit operations in many industries. Bubbly flows have been of interest among researchers in physics, medicine, chemistry and technology over the centuries. The aim of this thesis is to use advanced numerical methods for simulating microbubble in an aerated mixing tank. Main components of the mixing tank are a cylindrical vessel, a rotating Rushton turbine and the air nozzle. The objective of Computational Fluid Dynamics (CFD) is to predict fluid flow, heat transfer, mass transfer and chemical reactions. The CFD simulations of a turbulent bubbly flow are carried out in a cylindrical mixing tank using large eddy simulation (LES) and volume of fluid (VOF) method. The Rushton turbine induced flow is modeled by using a sliding mesh method. Numerical results are used to describe the bubbly flows in highly complex liquid flow. Some of the experimental works related to turbulent bubbly flow in a mixing tank are briefly reported. Numerical simulations are needed to complete and interpret the results of the experimental work. Information given by numerical simulations has a major role in designing and scaling-up mixing tanks. The results of this work have been reported in the following scientific articles: ·Honkanen M., Koohestany A., Hatunen T., Saarenrinne P., Zamankhan P., Large eddy simulations and PIV experiments of a two-phase air-water mixer, in Proceedings of ASME Fluids Engineering Summer Conference (2005). ·Honkanen M., Koohestany A., Hatunen T., Saarenrinne P., Zamankhan P., Dynamical States of Bubbling in an Aerated Stirring Tank, submitted to J. Computational Physics.

Adaptation of TRIZ Method for Problem Solving in Process Engineering

TRIZ is one of the well-known tools, based on analytical methods for creative problem solving. This thesis suggests adapted version of contradiction matrix, a powerful tool of TRIZ and few principles based on concept of original TRIZ. It is believed that the proposed version would aid in problem solving, especially those encountered in chemical process industries with unit operations. In addition, this thesis would help fresh process engineers to recognize importance of various available methods for creative problem solving and learn TRIZ method of creative problem solving. This thesis work mainly provides idea on how to modify TRIZ based method according to ones requirements to fit in particular niche area and solve problems efficiently in creative way. Here in this case, the contradiction matrix developed is based on review of common problems encountered in chemical process industry, particularly in unit operations and resolutions are based on approaches used in past to handle those issues.

Analyzing hydraulic head of fluid flow in underground porous medium

Hydraulic head is distributed through a medium with porous aspect. The analysis of hydraulic head from one point to another is used by the Richard's equation. This equation is equivalent to the groundwater ow equation that predicts the volumetric water contents. COMSOL 3.5 is used for computation applying Richard's equation. A rectangle of 100 meters of length and 10 meters of large (depth) with 0,1 m/s fl ux of inlet as source of our fl uid is simulated. The domain have Richards' equation model in two dimension (2D). Hydraulic head increases proportional with moisture content.

A Numerical study of fluid flow in peristaltic pumps and polymeric elastic pipes

This thesis work deals with a mathematical description of flow in polymeric pipe and in a specific peristaltic pump. This study involves fluid-structure interaction analysis in presence of complex-turbulent flows treated in an arbitrary Lagrangian-Eulerian (ALE) framework. The flow simulations are performed in COMSOL 4.4, as 2D axial symmetric model, and ABAQUS 6.14.1, as 3D model with symmetric boundary conditions. In COMSOL, the fluid and structure problems are coupled by monolithic algorithm, while ABAQUS code links ABAQUS CFD and ABAQUS Standard solvers with single block-iterative partitioned algorithm. For the turbulent features of the flow, the fluid model in both codes is described by RNG k-ϵ. The structural model is described, on the basis of the pipe material, by Elastic models or Hyperelastic Neo-Hookean models with Rayleigh damping properties. In order to describe the pulsatile fluid flow after the pumping process, the available data are often defective for the fluid problem. Engineering measurements are normally able to provide average pressure or velocity at a cross-section. This problem has been analyzed by McDonald's and Womersley's work for average pressure at fixed cross section by Fourier analysis since '50, while nowadays sophisticated techniques including Finite Elements and Finite Volumes exist to study the flow. Finally, we set up peristaltic pipe simulations in ABAQUS code, by using the same model previously tested for the fl uid and the structure.

A support system for evaluating a suitable joining method in the production of sheet metal goods

This report illustrates a comparative study of various joining methods involved in sheet metal production. In this report it shows the selection of joining methods, which includes comparing the advantages and disadvantages of a method over the other ones and choosing the best method for joining. On the basis of various joining process from references, a table is generated containing set of criterion that helps in evaluation of various sheet metal joining processes and in selecting the most suitable process for a particular product. Three products are selected and a comprehensive study of the joining methods is analyzed with the help of various parameters. The table thus is the main part of the analysis process of this study and can be advanced with the beneficial results. It helps in a better and easy understanding and comparing the various methods, which provides the foundation of this study and analysis. The suitability of the joining method for various types of cases of different sheet metal products can be tested with the help of this table. The sections of the created table display the requirements of manufacturing. The important factor has been considered and given focus in the table, as how the usage of these parameters is important in percentages according to particular or individual case. The analysis of the methods can be extended or altered by changing the parameters according to the constraint. The use of this table is demonstrated by pertaining the cases from sheet metal production.

Discription of a New Calculation Method for Determination of Stoichiometric Dissociation Constants of WeakAcid in Salt Solutions from Potentiometric Titration Data

A recently developed calculation method to determine stoichiometric dissociation constants of weak acids from potentiometric titration data is described. The titration data from three different weak acids in aqueous salt solutions at 25 °C were used as examples of the use of the method. The salt alone determined the ionic strength of the solutions considered in this study, and salt molalities up to 0,5 mol kg -1 were used.

A novel approachfor assessing the fatigue strength of ultrasonic impact treated welded structures.

It is commonly observed that complex fabricated structures subject tofatigue loading fail at the welded joints. Some problems can be corrected by proper detail design but fatigue performance can also be improved using post-weld improvement methods. In general, improvement methods can be divided into two main groups: weld geometry modification methods and residual stress modification methods. The former remove weld toe defects and/or reduce the stress concentrationwhile the latter introduce compressive stress fields in the area where fatigue cracks are likely to initiate. Ultrasonic impact treatment (UIT) is a novel post-weld treatment method that influences both the residual stress distribution andimproves the local geometry of the weld. The structural fatigue strength of non-load carrying attachments in the as-welded condition has been experimentally compared to the structural fatigue strength of ultrasonic impact treated welds. Longitudinal attachment specimens made of two thicknesses of steel S355 J0 have been tested for determining the efficiency of ultrasonic impacttreatment. Treated welds were found to have about 50% greater structural fatigue strength, when the slope of the S-N-curve is three. High mean stress fatigue testing based on the Ohta-method decreased the degree of weld improvement only 19%. This indicated that the method could be also applied for large fabricated structures operating under high reactive residual stresses equilibrated within the volume of the structure. The thickness of specimens has no significant effect tothe structural fatigue strength. The fatigue class difference between 5 mm and 8 mm specimen was only 8%. It was hypothesized that the UIT method added a significant crack initiation period to the total fatigue life of the welded joints. Crack initiation life was estimated by a local strain approach. Material parameters were defined using a modified Uniform Material Law developed in Germany. Finite element analysis and X-ray diffraction were used to define, respectively, the stress concentration and mean stress. The theoretical fatigue life was found to have good accuracy comparing to experimental fatigue tests.The predictive behaviour of the local strain approach combined with the uniformmaterial law was excellent for the joint types and conditions studied in this work.

Mass transfer and particleinteractions in granular systems and suspension flows

The purpose of this study was to investigate some important features of granular flows and suspension flows by computational simulation methods. Granular materials have been considered as an independent state ofmatter because of their complex behaviors. They sometimes behave like a solid, sometimes like a fluid, and sometimes can contain both phases in equilibrium. The computer simulation of dense shear granular flows of monodisperse, spherical particles shows that the collisional model of contacts yields the coexistence of solid and fluid phases while the frictional model represents a uniform flow of fluid phase. However, a comparison between the stress signals from the simulations and experiments revealed that the collisional model would result a proper match with the experimental evidences. Although the effect of gravity is found to beimportant in sedimentation of solid part, the stick-slip behavior associated with the collisional model looks more similar to that of experiments. The mathematical formulations based on the kinetic theory have been derived for the moderatesolid volume fractions with the assumption of the homogeneity of flow. In orderto make some simulations which can provide such an ideal flow, the simulation of unbounded granular shear flows was performed. Therefore, the homogeneous flow properties could be achieved in the moderate solid volume fractions. A new algorithm, namely the nonequilibrium approach was introduced to show the features of self-diffusion in the granular flows. Using this algorithm a one way flow can beextracted from the entire flow, which not only provides a straightforward calculation of self-diffusion coefficient but also can qualitatively determine the deviation of self-diffusion from the linear law at some regions nearby the wall inbounded flows. Anyhow, the average lateral self-diffusion coefficient, which was calculated by the aforementioned method, showed a desirable agreement with thepredictions of kinetic theory formulation. In the continuation of computer simulation of shear granular flows, some numerical and theoretical investigations were carried out on mass transfer and particle interactions in particulate flows. In this context, the boundary element method and its combination with the spectral method using the special capabilities of wavelets have been introduced as theefficient numerical methods to solve the governing equations of mass transfer in particulate flows. A theoretical formulation of fluid dispersivity in suspension flows revealed that the fluid dispersivity depends upon the fluid properties and particle parameters as well as the fluid-particle and particle-particle interactions.

Description and implementation of cut calculation-method in a sawmill

Työn tavoitteena oli kuvata ja ottaa käyttöön sahauseräkohtaisen kannattavuuden laskentamenetelmä sahalle, sekä tehdä laskentamalli menetelmän tueksi. Sahauksen peruskäsitteiden jälkeen työssä on esitelty sahan tuotantoprosessi. Tuotantoprosessi on kuvattu kirjallisuuden ja asiantuntijoiden haastattelujen perusteella. Seuraavaksi kartoitettiin hyötyjä ja vaikutuksia, mitä laskentamenetelmältä odotetaan.. Kustannuslaskennan teoriaa selvitettiin kirjallisuuslähteitä käyttäen silmälläpitäen juuri tätä kehitettävää laskentamenetelmää. Lisäksi esiteltiin Uimaharjun sahalla käytettävät ja laskentaan liittyvät laskenta- ja tietojärjestelmät.Nykyisin sahalla ei ole minkäänlaista menetelmää sahauseräkohtaisen tuloksen laskemiseksi. Pienillä muutoksilla sahan tietojärjestelmään ja prosessikoneisiin voidaan sahauserä kuljettaa prosessin läpi niin, että jokaisessa prosessin vaiheessa sille saadaan kohdistettua tuotantotietoa. Eri vaiheista saatua tietoa käyttämällä saadaan tarkasti määritettyä tuotteet, joita sahauserä tuotti ja paljonko tuotantoresursseja tuottamiseen kului. Laskentamalliin syötetään tuotantotietoja ja kustannustietoa ja saadaan vastaukseksi sahauserän taloudellinen tulos.Toimenpide ehdotuksena esitetään lisätutkimusta tuotantotietojen automaattisesta keräämisestä manuaalisen työn ja virheiden poistamiseksi. Suhteellisen pienillä panoksilla voidaan jokaiselle sahauserälle kerätä tuotantotiedot täysin automaattisesti. Lisäksi kehittämäni laskentamallin tilalle tulisi hankkia sovellus, joka käyttäisi paremmin hyväksi nykyisiä tietojärjestelmiä ja poistaisi manuaalisen työvaiheen laskennassa.

On Commutation and Conjugacy of Rational Languages and the Fixed Point Method

The research on language equations has been active during last decades. Compared to the equations on words the equations on languages are much more difficult to solve. Even very simple equations that are easy to solve for words can be very hard for languages. In this thesis we study two of such equations, namely commutation and conjugacy equations. We study these equations on some limited special cases and compare some of these results to the solutions of corresponding equations on words. For both equations we study the maximal solutions, the centralizer and the conjugator. We present a fixed point method that we can use to search these maximal solutions and analyze the reasons why this method is not successful for all languages. We give also several examples to illustrate the behaviour of this method.

Objective Evaluation of Nasal Dimensions in Children with Acoustic Rhinometry

Aims: This study was carried out to investigate the usefulness of acoustic rhinometry in the evaluation of intranasal dimensions in children. The aim was to define reference values for school children. In addition, the role of the VAS scale in the subjective evaluation of nasal obstruction in children was studied. Materials and methods: Measurements were done with Acoustic Rhinometry A1. The values of special interest were the minimal cross-sectional area (MCA) and the anterior volume of the nose (VOL). The data for reference values included 124 voluntary school children with no permanent nasal symptoms, aged between 7 and 14 years. Data were collected at baseline and after decongestion of the nose; the VAS scale was filled in before measurements. The subjects in the follow-up study (n=74, age between 1 and 12 years) were receiving intranasal spray of insulin or placebo. The nasal symptoms were recorded and acoustic rhinometry was measured at each control visit. Results: In school children, the mean total MCA was 0.752 cm2 (SD 0.165), and the mean total VOL was 4.00 cm3 (SD 0.63) at baseline. After decongestion, a significant increase in the mean TMCA and in the mean TVOL was found. A correlation was found between TMCA and age, and between TVOL and height of a child. There was no difference between boys and girls. A correlation was found between unilateral acoustic values and VAS at baseline, but not after decongestion. No difference wasfound in acoustic values or symptoms between the insulin and placebo group in the follow-up study of two years. Conclusions: Acoustic rhinometry is a suitable objective method to examine intranasal dimensions in children. It is easy to perform and well tolerated. Reference values for children between 7 and 14 years were established.

Future prospects of short distance rail based passenger transport in Northern Europe

The worlds’ population is increasing and cities have become more crowded with people and vehicles. Communities in the fringe of metropolitans’ increase the traffic done with private cars, but also increase the need for public transportation. People have typically needs traveling to work located in city centers during the morning time, and return to suburbs in the afternoon or evening. Rail based passenger transport is environmentally friendly transport mode with high capacity to transport large volume of people. Railways have been regulated markets with national incumbent having monopoly position. Opening the market for competition is believed to have a positive effect by increasing the efficiency of the industry. National passenger railway market is opened for competition only in few countries, where as international traffic in EU countries was deregulated in 2010. The objective of this study is to examine the passenger railway market of three North European countries, Sweden, Denmark and Estonia. The interest was also to get an understanding of the current situation and how the deregulation has proceeded. Theory of deregulation is unfolded with literature analyses and empirical part of the study is constructed from two parts. Customer satisfaction survey was chosen as a method to collect real life experiences from the passengers and measure their knowledge of the market situation and possible changes appeared. Interviews of experts from the industry and labor unions give more insights and able better understanding for example of social consequences caused from opening the market for competition. Expert interviews were conducted by using semi-structured theme interview. Based on the results of this study, deregulation has proceeded quite differently in the three countries researched. Sweden is the most advanced country, where the passenger railway market is open for new entrants. Denmark and Estonia are lagging behind. Opening the market is considered positive among passengers and most of the experts interviewed. Common for the interviews were the labour unions negative perspective concerning deregulation. Despite the fact deregulation is considered positive among the respondents of the customer satisfaction survey, they could not name railway undertakings operating in their country. Generally respondents were satisfied with the commuter trains. Ticket price, punctuality of trains and itinerary affect the most to customer satisfaction.