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.

Ferromagnetic properties of pyrolytic nanocarbon

The work reported in this thesis is dedicated to irreversible magnetic properties in pyrolytic nanocarbon samples. Based on atomic force microscope images, the samples consist of carbon clusters with radius 30..120 nm. These are treated as single-domain nanoparticles. Magnetic hysteresis, field cooled, zero field cooled and thermoremanent magnetization measurements were performed using an RF SQUID magnetometer and ferromagnetic behaviour was observed. Analysis suggests that the ferromagnetic ordering is associated with defects in a thin surface layer, whose thickness is independent of particle size. Critical radius for single-domain particles, critical radius for coherent rotation, magnetic layer thickness, distance between elementary magnetic moments, saturation magnetization, exchange stiffness constant and anisotropy energy density are also presented.

Leptogenesis as the origin of matter-antimatter asymmetry in extra dimensional and supersymmetric models

The Standard Model of particle physics is currently the best description of fundamental particles and their interactions. All particles save the Higgs boson have been observed in particle accelerator experiments over the years. Despite the predictive power the Standard Model there are many phenomena that the scenario does not predict or explain. Among the most prominent dilemmas is matter-antimatter asymmetry, and much effort has been made in formulating scenarios that accurately predict the correct amount of matter-antimatter asymmetry in the universe. One of the most appealing explanations is baryogenesis via leptogenesis which not only serves as a mechanism of producing excess matter over antimatter but can also explain why neutrinos have very small non-zero masses. Interesting leptogenesis scenarios arise when other possible candidates of theories beyond the Standard Model are brought into the picture. In this thesis, we have studied leptogenesis in an extra dimensional framework and in a modified version of supersymmetric Standard Model. The first chapters of this thesis introduce the standard cosmological model, observations made on the photon to baryon ratio and necessary preconditions for successful baryogenesis. Baryogenesis via leptogenesis is then introduced and its connection to neutrino physics is illuminated. The final chapters concentrate on extra dimensional theories and supersymmetric models and their ability to accommodate leptogenesis. There, the results of our research are also presented.

Multifrequency VLBI Observations of Selected Active Galactic Nuclei

In this dissertation, active galactic nuclei (AGN) are discussed, as they are seen with the high-resolution radio-astronomical technique called Very Long Baseline Interferometry (VLBI). This observational technique provides very high angular resolution (_ 10−300 = 1 milliarcsecond). VLBI observations, performed at different radio frequencies (multi-frequency VLBI), allow to penetrate deep into the core of an AGN to reveal an otherwise obscured inner part of the jet and the vicinity of the AGN’s central engine. Multi-frequency VLBI data are used to scrutinize the structure and evolution of the jet, as well as the distribution of the polarized emission. These data can help to derive the properties of the plasma and the magnetic field, and to provide constraints to the jet composition and the parameters of emission mechanisms. Also VLBI data can be used for testing the possible physical processes in the jet by comparing observational results with results of numerical simulations. The work presented in this thesis contributes to different aspects of AGN physics studies, as well as to the methodology of VLBI data reduction. In particular, Paper I reports evidence of optical and radio emission of AGN coming from the same region in the inner jet. This result was obtained via simultaneous observations of linear polarization in the optical and in radio using VLBI technique of a sample of AGN. Papers II and III describe, in detail, the jet kinematics of the blazar 0716+714, based on multi-frequency data, and reveal a peculiar kinematic pattern: plasma in the inner jet appears to move substantially faster that that in the large-scale jet. This peculiarity is explained by the jet bending, in Paper III. Also, Paper III presents a test of the new imaging technique for VLBI data, the Generalized Maximum Entropy Method (GMEM), with the observed (not simulated) data and compares its results with the conventional imaging. Papers IV and V report the results of observations of the circularly polarized (CP) emission in AGN at small spatial scales. In particular, Paper IV presents values of the core CP for 41 AGN at 15, 22 and 43 GHz, obtained with the help of the standard Gain transfer (GT) method, which was previously developed by D. Homan and J.Wardle for the calibration of multi-source VLBI observations. This method was developed for long multi-source observations, when many AGN are observed in a single VLBI run. In contrast, in Paper V, an attempt is made to apply the GT method to single-source VLBI observations. In such observations, the object list would include only a few sources: a target source and two or three calibrators, and it lasts much shorter than the multi-source experiment. For the CP calibration of a single-source observation, it is necessary to have a source with zero or known CP as one of the calibrators. If the archival observations included such a source to the list of calibrators, the GT could also be used for the archival data, increasing a list of known AGN with the CP at small spatial scale. Paper V contains also calculation of contributions of different sourced of errors to the uncertainty of the final result, and presents the first results for the blazar 0716+714.

GAS EXCHANGE IN YOUNG PLANTS OF Tabebuia aurea(Bignoniaceae Juss.) SUBJECTED TO FLOODING STRESS1

ABSTRACT The Paratudo (Tabebuia aurea) is a species occurring in the Pantanal of Miranda, Mato Grosso do Sul, Brazil, an area characterized by seasonal flooding. To evaluate the tolerance of this plant to flooding, plants aged four months were grown in flooded soil and in non-flooded soil (control group). Stomatal conductance, transpiration and CO2 assimilation were measured during the stress (48 days) and recovery (11 days) period, totalling 59 days. The values of stomatal conductance of the control group and stressed plants at the beginning of the flooded were 0.33 mol m-2s-1 and reached 0.02 mol m-2 s-1 (46th day) at the end of this event. For the transpiration parameter, the initial rate was 3.1 mol m s-1, and the final rate reached 0.2 or 0.3 mol m-2 s-1 (47/48 th day). The initial photosynthesis rate was 8.9 mmol m-2s-1 and oscillated after the sixth day, and the rate reached zero on the 48th day. When the photosynthesis rate reached zero, the potted plants were dried, and the rate was analyzed (11th day). The following values were obtained for dried plants: stomatal conductance = 0.26 mol m-2 s-1, transpiration rate = 2.5 mol m-2 s-1 and photosynthesis rate = 7.8 mmol m-2 s-1. Flooded soil reduced photosynthesis and stomatal conductance, leading to the hypertrophy of the lenticels. These parameters recovered and after this period, and plants exhibited tolerance to flooding stress by reducing their physiological activities.

Modeling study of limestone reactions in fluidized bed conditions

Traditionally limestone has been used for the flue gas desulfurization in fluidized bed combustion. Recently, several studies have been carried out to examine the use of limestone in applications which enable the removal of carbon dioxide from the combustion gases, such as calcium looping technology and oxy-fuel combustion. In these processes interlinked limestone reactions occur but the reaction mechanisms and kinetics are not yet fully understood. To examine these phenomena, analytical and numerical models have been created. In this work, the limestone reactions were studied with aid of one-dimensional numerical particle model. The model describes a single limestone particle in the process as a function of time, the progress of the reactions and the mass and energy transfer in the particle. The model-based results were compared with experimental laboratory scale BFB results. It was observed that by increasing the temperature from 850 °C to 950 °C the calcination was enhanced but the sulfate conversion was no more improved. A higher sulfur dioxide concentration accelerated the sulfation reaction and based on the modeling, the sulfation is first order with respect to SO2. The reaction order of O2 seems to become zero at high oxygen concentrations.

The position effect of structural Eucalyptus round timber on the flexural modulus of elasticity

Round timber has great use in civil construction, performing the function of beams, columns, foundations, poles for power distribution among others, with the advantage of not being processed, such as lumber. The structural design of round timber requires determining the elastic properties, mainly the modulus of elasticity. The Brazilian standards responsible for the stiffness and strength determination of round timber are in effect for over twenty years with no technical review. Round timber, for generally present an axis with non-zero curvature according to the position of the element in the bending test, may exhibit different values of modulus of elasticity. This study aims to analyze the position effect of Eucalyptus grandis round timber on the flexural modulus of elasticity. The three-point bending test was evaluated in two different positions based on the longitudinal rotation of the round timber element. The results revealed that at least two different positions of the round timber element are desired to obtain significant modulus of elasticity.

Observations of Industry Momentum Effect in European Equity Indices

This thesis examines the existence and nature of momentum effect in European equity indices. A set of predefined indicators is used to compose momentum portfolios and different holding periods are used to test the strategies over variable time periods as well as under different economical conditions. The data consists of daily closing prices of STOXX Europe 600 index and its 18 super sector indices. Over the study period we follow the performances of a long position in the Winner portfolio, a position in the market neutral zero-cost portfolio and also a position in the risk-controlled zero-cost portfolio. The investment ratio of the risk-controlled zero-cost portfolio is negatively correlated with the realized market volatility. The results show that momentum effect is present in European industries and is most prominent in the short-term. Indicators that are based on short-term performance tend predict the over- and underperformers for the 1-month holding period more reliably than any other indicator/holding period combination. The examination of the strategies under different economical conditions shows that the market neutral approach can create significant returns in times of recession but in times of economic boom the long position in Winner portfolio outperforms the market neutral portfolio by an extensive margin.

Effect of Hevea brasiliensis latex sap gel on healing of acute skin wounds induced on the back of rats

Objective : to evaluate the effect of topical delivery of latex cream-gel in acute cutaneous wounds induced on the back of rats. Methods : we subjected sixteen rats to dermo-epidermal excision of a round dorsal skin flap, with 2.5cm diameter. We divided the animals into two groups: Latex Group: application of cream-gel-based latex throughout the wound bed on postoperative days zero, three, six and nine; Control group: no treatment on the wound. Photographs of the lesions were taken on the procedure day and on the 6th and 14th postoperative days, for analyzing the area and the larger diameter of the wound. We carried out euthanasia of all animals on the 14th postoperative day, when we resected he dorsal skin and the underlying muscle layer supporting the wound for histopathological study. Results : there was no statistically significant difference in the percentage of wound closure, in the histopathological findings or in the reduction of the area and of the largest diameter of the wounds among the groups studied on the 14th postoperative day. Conclusion : according to the experimental conditions in which the study was conducted, latex cream-gel did not interfere in the healing of acute cutaneous wounds in rats.

Impact of Sacrospinous Colpopexy Associated with Anterior Colporrhaphy for the Treatment of Dome Prolapse on all Three Vaginal Compartments

Objective To evaluate the results of sacrospinous colpopexy surgery associated with anterior colporrhaphy for the treatment of women with post-hysterectomy vaginal vault prolapse. Methods This prospective study included 20women with vault prolapse, PelvicOrgan Prolapse Quantification System (POP-Q) stage≥2, treated between January 2003 and February 2006, and evaluated in a follow-up review (more than one year later). Genital prolapse was evaluated qualitatively in stages and quantitatively in centimeters. Prolapse stage < 2 was considered to be the cure criterion. Statistical analysis was performed using the Wilcoxon test (paired samples) to compare the points and stages of prolapse before and after surgery. Results Evaluation of the vaginal vault after one year revealed that 95% of subjects were in stage zero and that 5% were in stage 1. For cystocele, 50% were in stage 1, 10% were in stage 0 (cured) and 40% were in stage 2. For rectocele, three women were in stage 1 (15%), one was in stage 2 (5%) and 16 had no further prolapse. The most frequent complication was pain in the right buttock, with remission of symptoms in all three cases three months after surgery. Conclusions In this retrospective study, the surgical correction of vault prolapse using a sacrospinous ligament fixation technique associatedwith anterior colporrhaphy proved effective in resolving genital prolapse. Despite the low complication rates, there was a high rate of cystocele, which may be caused by posterior vaginal shifting due to either the technique or an overvaluation by the POP-Q system.

Hydrocortisone decreases apoptosis in jejunum of horses subjected to experimental ischemia and reperfusion

In order to evaluate the effect of hydrocortisone on apoptosis in the jejunum of horses subjected to ischemia and reperfusion, ten horses were paired and grouped into two groups - treated (n=5) and non treated (n=5). Segments of the jejunum were used as controls (C), or as venous ischemia (VIsc), which were subjected to 2h of ischemia followed by 2 or 12h of reperfusion. C samples were collected at time zero (prior to ischemia) and VIsc samples were collected at 2h of ischemia and at 2 and 12h of reperfusion. TUNEL positive apoptotic cells were counted in 10 microscopical fields in deep mucosa from each horse throughout the time course. After 12h of reperfusion, the number of apoptotic cells in treated group were significantly lower than in untreated animals, indicating that hydrocortisone inhibits apoptosis. These results indicate that hydrocortisone has a beneficial effects favoring the maintenance of jejunal integrity in horses with ischemia and reperfusion injuries by preventing apoptotic cell death.

Oxidation rates of carbon and nitrogen in char residues from solid fuels

Computational fluid dynamics (CFD) modeling is an important tool in designing new combustion systems. By using CFD modeling, entire combustion systems can be modeled and the emissions and the performance can be predicted. CFD modeling can also be used to develop new and better combustion systems from an economical and environmental point of view. In CFD modeling of solid fuel combustion, the combustible fuel is generally treated as single fuel particles. One of the limitations with the CFD modeling concerns the sub-models describing the combustion of single fuel particles. Available models in the scientific literature are in many cases not suitable as submodels for CFD modeling since they depend on a large number of input parameters and are computationally heavy. In this thesis CFD-applicable models are developed for the combustion of single fuel particles. The single particle models can be used to improve the combustion performance in various combustion devices or develop completely new technologies. The investigated fields are oxidation of carbon (C) and nitrogen (N) in char residues from solid fuels. Modeled char-C oxidation rates are compared to experimental oxidation rates for a large number of pulverized solid fuel chars under relevant combustion conditions. The experiments have been performed in an isothermal plug flow reactor operating at 1123-1673 K and 3-15 vol.% O2. In the single particle model, the char oxidation is based on apparent kinetics and depends on three fuel specific parameters: apparent pre-exponential factor, apparent activation energy, and apparent reaction order. The single particle model can be incorporated as a sub-model into a CFD code. The results show that the modeled char oxidation rates are in good agreement with experimental char oxidation rates up to around 70% of burnout. Moreover, the results show that the activation energy and the reaction order can be assumed to be constant for a large number of bituminous coal chars under conditions limited by the combined effects of chemical kinetics and pore diffusion. Based on this, a new model based on only one fuel specific parameter is developed (Paper III). The results also show that reaction orders of bituminous coal chars and anthracite chars differ under similar conditions (Paper I and Paper II); reaction orders of bituminous coal chars were found to be one, while reaction orders of anthracite chars were determined to be zero. This difference in reaction orders has not previously been observed in the literature and should be considered in future char oxidation models. One of the most frequently used comprehensive char oxidation models could not explain the difference in the reaction orders. In the thesis (Paper II), a modification to the model is suggested in order to explain the difference in reaction orders between anthracite chars and bituminous coal chars. Two single particle models are also developed for the NO formation and reduction during the oxidation of single biomass char particles. In the models the char-N is assumed to be oxidized to NO and the NO is partly reduced inside the particle. The first model (Paper IV) is based on the concentration gradients of NO inside and outside the particle and the second model is simplified to such an extent that it is based on apparent kinetics and can be incorporated as a sub-model into a CFD code (Paper V). Modeled NO release rates from both models were in good agreement with experimental measurements from a single particle reactor of quartz glass operating at 1173-1323 K and 3-19 vol.% O2. In the future, the models can be used to reduce NO emissions in new combustion systems.

Identification of Some Additional Loss Components in High-Power Low-Voltage Permanent Magnet Generators

Permanent magnet generators (PMG) represent the cutting edge technology in modern wind mills. The efficiency remains high (over 90%) at partial loads. To improve the machine efficiency even further, every aspect of machine losses has to be analyzed. Additional losses are often given as a certain percentage without providing any detailed information about the actual calculation process; meanwhile, there are many design-dependent losses that have an effect on the total amount of additional losses and that have to be taken into consideration. Additional losses are most often eddy current losses in different parts of the machine. These losses are usually difficult to calculate in the design process. In this doctoral thesis, some additional losses are identified and modeled. Further, suggestions on how to minimize the losses are given. Iron losses can differ significantly between the measured no-load values and the loss values under load. In addition, with embedded magnet rotors, the quadrature-axis armature reaction adds losses to the stator iron by manipulating the harmonic content of the flux. It was, therefore, re-evaluated that in salient pole machines, to minimize the losses and the loss difference between the no-load and load operation, the flux density has to be kept below 1.5 T in the stator yoke, which is the traditional guideline for machine designers. Eddy current losses may occur in the end-winding area and in the support structure of the machine, that is, in the finger plate and the clamping ring. With construction steel, these losses account for 0.08% of the input power of the machine. These losses can be reduced almost to zero by using nonmagnetic stainless steel. In addition, the machine housing may be subjected to eddy current losses if the flux density exceeds 1.5 T in the stator yoke. Winding losses can rise rapidly when high frequencies and 10–15 mm high conductors are used. In general, minimizing the winding losses is simple. For example, it can be done by dividing the conductor into transposed subconductors. However, this comes with the expense of an increase in the DC resistance. In the doctoral thesis, a new method is presented to minimize the winding losses by applying a litz wire with noninsulated strands. The construction is the same as in a normal litz wire but the insulation between the subconductors has been left out. The idea is that the connection is kept weak to prevent harmful eddy currents from flowing. Moreover, the analytical solution for calculating the AC resistance factor of the litz-wire is supplemented by including an end-winding resistance in the analytical solution. A simple measurement device is developed to measure the AC resistance in the windings. In the case of a litz-wire with originally noninsulated strands, vacuum pressure impregnation (VPI) is used to insulate the subconductors. In one of the two cases studied, the VPI affected the AC resistance factor, but in the other case, it did not have any effect. However, more research is needed to determine the effect of the VPI on litz-wire with noninsulated strands. An empirical model is developed to calculate the AC resistance factor of a single-layer formwound winding. The model includes the end-winding length and the number of strands and turns. The end winding includes the circulating current (eddy currents that are traveling through the whole winding between parallel strands) and the main current. The end-winding length also affects the total AC resistance factor.

A simulation model of the optimal maintenance strategy: conceptual agent-based approach

In this Master’s thesis agent-based modeling has been used to analyze maintenance strategy related phenomena. The main research question that has been answered was: what does the agent-based model made for this study tell us about how different maintenance strategy decisions affect profitability of equipment owners and maintenance service providers? Thus, the main outcome of this study is an analysis of how profitability can be increased in industrial maintenance context. To answer that question, first, a literature review of maintenance strategy, agent-based modeling and maintenance modeling and optimization was conducted. This review provided the basis for making the agent-based model. Making the model followed a standard simulation modeling procedure. With the simulation results from the agent-based model the research question was answered. Specifically, the results of the modeling and this study are: (1) optimizing the point in which a machine is maintained increases profitability for the owner of the machine and also the maintainer with certain conditions; (2) time-based pricing of maintenance services leads to a zero-sum game between the parties; (3) value-based pricing of maintenance services leads to a win-win game between the parties, if the owners of the machines share a substantial amount of their value to the maintainers; and (4) error in machine condition measurement is a critical parameter to optimizing maintenance strategy, and there is real systemic value in having more accurate machine condition measurement systems.

Investigation on diamond turning of silicon crystal - generation mechanism of surface cut with worn tool

This paper discusses the effect of tool wear on surface finish in single-point diamond turning of single crystal silicon. The morphology and topography of the machined surface clearly show the type of cutting edge wear reproduced onto the cutting grooves. Scanning electron microscopy is used in order to correlate the cutting edge damage and microtopography features observed through atomic force microscopy. The possible wear mechanisms affecting tool performance and surface generation during cutting are also discussed. The zero degree rake angle single point diamond tool presented small nicks on the cutting edge. The negative rake angle tools presented more a type of crater wear on the rake face. No wear was detected on flank face of the diamond tools.