Gas-solid two-phase flow in the riser of circulating fluidized beds: mathematical modelling and numerical simulation

A mathematical model is developed for gas-solids flows in circulating fluidized beds. An Eulerian formulation is followed based on the two-fluids model approach where both the fluid and the particulate phases are treated as a continuum. The physical modelling is discussed, including the formulation of boundary conditions and the description of the numerical methodology. Results of numerical simulation are presented and discussed. The model is validated through comparison to experiment, and simulation is performed to investigate the effects on the flow hydrodynamics of the solids viscosity.

A Design Methodology for the Compensation of Positioning Deviation in Gantry Manipulators

This work presents recent results concerning a design methodology used to estimate the positioning deviation for a gantry (Cartesian) manipulator, related mainly to structural elastic deformation of components during operational conditions. The case-study manipulator is classified as gantry type and its basic dimensions are 1,53m x 0,97m x 1,38m. The dimensions used for the calculation of effective workspace due to end-effector path displacement are: 1m x 0,5m x 0,5m. The manipulator is composed by four basic modules defined as module X, module Y, module Z and terminal arm, where is connected the end-effector. Each module controlled axis performs a linear-parabolic positioning movement. The planning path algorithm has the maximum velocity and the total distance as input parameters for a given task. The acceleration and deceleration times are the same. Denavit-Hartemberg parameterization method is used in the manipulator kinematics model. The gantry manipulator can be modeled as four rigid bodies with three degrees-of-freedom in translational movements, connected as an open kinematics chain. Dynamic analysis were performed considering inertial parameters specification such as component mass, inertia and center of gravity position of each module. These parameters are essential for a correct manipulator dynamic modelling, due to multiple possibilities of motion and manipulation of objects with different masses. The dynamic analysis consists of a mathematical modelling of the static and dynamic interactions among the modules. The computation of the structural deformations uses the finite element method (FEM).

Petri Net approach for modelling system integration in intelligent buildings

In this paper, a Petri Net approach is introduced for modelling and simulation of control strategies in Intelligent Building. In this context, it is claimed that integration with other building systems can be achieved in a more systematic way considering a mechatronic approach (i.e. multidisciplinary concepts applied to the development of systems). The case study is the Ambulatory Building of Medical School Hospital of University of São Paulo. Particularly, the developed methodology is applied to the elevator system and to the HVAC (Heating, Ventilation and Air Conditioning) system. It is shown that using this approach, the control systems could be integrated, improving performance.

Modelling flexible asset management in industrial maintenance companies and networks

The tightening competition and increasing dynamism have created an emerging need for flexible asset management. This means that the changes of market demand should be responded to with adjustments in the amount of assets tied to the balance sheets of companies. On the other hand, industrial maintenance has recently experienced drastic changes, which have led to an increase in the number of maintenance networks (consisting of customer companies that buy maintenance services, as well as various supplier companies) and inter-organizational partnerships. However, the research on maintenance networks has not followed the changes in the industry. Instead, there is a growing need for new ways of collaboration between partnering companies to enhance the competitiveness of the whole maintenance network. In addition, it is more and more common for companies to pursue lean operations in their businesses. This thesis shows how flexible asset management can increase the profitability of maintenance companies and networks under dynamic operating conditions, and how the additional value can then be shared between the network partners. Firstly, I have conducted a systematic literature review to identify what kind of requirements for asset management models are set by the increasing dynamism. Then I have responded to these requirements by constructing an analytical model for flexible asset management, linking asset management to the profitability and financial state of a company. The thesis uses the model to show how flexible asset management can increase profitability in maintenance companies and networks, and how the created value can be shared in the networks to reach a win-win situation. The research indicates that the existing models for asset management are heterogeneous by nature due to the various definitions of ‘asset management’. I conclude that there is a need for practical asset management models which address assets comprehensively with an inter-organizational, strategic view. The comprehensive perspective, taking all kinds of asset types into account, is needed to integrate the research on asset management with the strategic management of companies and networks. I will show that maintenance companies can improve their profitability by increasing the flexibility of their assets. In maintenance networks, reorganizing the ownership of the assets among the different network partners can create additional value. Finally, I will introduce flexible asset management contracts for maintenance networks. These contracts address the value sharing related to reorganizing the ownership of assets according to the principles of win-win situations.

Modelling Long-Run Relationship between Spot and Future Prices of Different Commodities

The purpose of this study is to investigate whether there exists any kind of relationship between the spot and future prices of the different commodities or not. Commodities like cocoa, coffee, crude oil, gold, natural gas and silver are considered from January 3, 2000 to December 31, 2012. For this purpose, ADF test and KPSS test are used in testing the stationarity whereas Johansen Cointegration test is used in testing the long-run relationship. Johansen co-integration test exhibits that there at least 5 co-integrating pairs out of 6 except crude oil. Moreover, the result of Granger Causality supports the fact that if two or more than two time series tend to be co-integrated there exists either uni-directional or bi-directional relationship. However, our results reveled that although there exists the co-integration between the variable, one might not granger causes another .VAR model is also used to measure the proportion of effects. These findings will help the derivative market and arbitragers in developing the strategies to gain the maximum profit in the financial market.

Probabilistic position and orientation estimation of rigid objects using an RGB-D sensor

The recent emergence of low-cost RGB-D sensors has brought new opportunities for robotics by providing affordable devices that can provide synchronized images with both color and depth information. In this thesis, recent work on pose estimation utilizing RGBD sensors is reviewed. Also, a pose recognition system for rigid objects using RGB-D data is implemented. The implementation uses half-edge primitives extracted from the RGB-D images for pose estimation. The system is based on the probabilistic object representation framework by Detry et al., which utilizes Nonparametric Belief Propagation for pose inference. Experiments are performed on household objects to evaluate the performance and robustness of the system.

Kattilan ja kanaviston eristykset 3D-laitossuunnitteluohjelmassa

Työn tavoitteena oli kehittää kanavistojen ja kattilan eristyksien mallintamista 3D-laitossuunitteluohjelmassa sekä siitä saatavia materiaalilistoja. Kanavistojen eristeitä kuvastavaa geometriaa muokattiin vastaamaan eristyksen todellista tilanvarausta. Kattilan eristyksen mallintaminen on aikaisemmin jouduttu tekemään manuaalisesti, työssä kehitettiin työkalu, jonka avulla eristysten materiaalitiedot ja oikeat tilanvaraukset saadaan lisättyä 3D-malliin. Kaikkia eristyksen mallintamista ei kuitenkaan pystytä mielekkäästi toteuttamaan automaattisesti työkalujen avulla, näitä kohteita ovat tukikehikkorakenteet kuten kattilan vinttikehikko. Lisäksi täytyi selvittää, mitä tietoja eristyksien materiaalimäärälistoille vaaditaan, ja kuinka nämä tiedot saataisiin listoihin mukaan. Jotta kaikki eristettävät kohteet saadaan tuotettaviin materiaalilistoihin, on ohjelmaan lisättävä käytettäväksi uusia komponentteja. Näitä ovat uudet kanavistojen osat sekä muun muassa venturit ja huoltoluukut kanavistoille. Kattilaan asennettavat laitteet kuten näkölasit ja polttimet vaativat tukirakenteen, joka samalla toimii kotelona putkistojen taivutuksille. Näitä kohteita on aikaisemmin tuotu vanhoista projekteista tai ne on tehty suunnittelijan toimesta uusiksi. Työssä esitettiin uusi tapa mallintaa kattilan pinnat, jolloin niihin voidaan liittää koteloita erillisestä komponenttikirjastosta. Kanavistojen eristyksien materiaalimäärälistojen tuottaminen suoraan ohjelmasta auto-matisoitiin. Kattilan ja siihen liittyvien eristyskehysrakenteiden listojen tuotanto kuitenkin vaatii yhteisten toimintamallien kehittämistä. Tämän työn uudet menetelmät mahdollistavat yhä uusien työkalujen kehittämisen, esimerkiksi automaattisen eristysten teknisten piirustusten kuvatuotannon tulevaisuudessa.

Modelling of hydrogen stratification and gas mixing in containment using the APROS code

Hydrogen stratification and atmosphere mixing is a very important phenomenon in nuclear reactor containments when severe accidents are studied and simulated. Hydrogen generation, distribution and accumulation in certain parts of containment may pose a great risk to pressure increase induced by hydrogen combustion, and thus, challenge the integrity of NPP containment. The accurate prediction of hydrogen distribution is important with respect to the safety design of a NPP. Modelling methods typically used for containment analyses include both lumped parameter and field codes. The lumped parameter method is universally used in the containment codes, because its versatility, flexibility and simplicity. The lumped parameter method allows fast, full-scale simulations, where different containment geometries with relevant engineering safety features can be modelled. Lumped parameter gas stratification and mixing modelling methods are presented and discussed in this master’s thesis. Experimental research is widely used in containment analyses. The HM-2 experiment related to hydrogen stratification and mixing conducted at the THAI facility in Germany is calculated with the APROS lump parameter containment package and the APROS 6-equation thermal hydraulic model. The main purpose was to study, whether the convection term included in the momentum conservation equation of the 6-equation modelling gives some remarkable advantages compared to the simplified lumped parameter approach. Finally, a simple containment test case (high steam release to a narrow steam generator room inside a large dry containment) was calculated with both APROS models. In this case, the aim was to determine the extreme containment conditions, where the effect of convection term was supposed to be possibly high. Calculation results showed that both the APROS containment and the 6-equation model could model the hydrogen stratification in the THAI test well, if the vertical nodalisation was dense enough. However, in more complicated cases, the numerical diffusion may distort the results. Calculation of light gas stratification could be probably improved by applying the second order discretisation scheme for the modelling of gas flows. If the gas flows are relatively high, the convection term of the momentum equation is necessary to model the pressure differences between the adjacent nodes reasonably.

Seed bank modelling of volunteer oil seed rape: from seeds fate in the soil to seedling emergence

Studies were conducted to estimate parameters and relationships associated with sub-processes in soil seed banks of oilseed rape in Gorgan, Iran. After one month of burial, seed viability decreased to 39%, with a slope of 2.03% per day, and subsequently decreased with a lower slope of 0.01 until 365 days following burial in the soil. Germinability remained at its highest value in autumn and winter and decreased from spring to the last month of summer. Non-dormant seeds of volunteer oilseed rape did not germinate at temperatures lower than 3.8 ºC and a water potential of -1.4 MPa ºd. The hydrothermal values were 36.2 and 42.9 MPa ºd for sub- and supra-optimal temperatures, respectively. Quantification of seed emergence as influenced by burial depth was performed satisfactorily (R² = 0.98 and RMSE = 5.03). The parameters and relationships estimated here can be used for modelling soil seed bank dynamics or establishing a new model for the environment.

LVDC power distribution system: computational modelling

In the doctoral dissertation, low-voltage direct current (LVDC) distribution system stability, supply security and power quality are evaluated by computational modelling and measurements on an LVDC research platform. Computational models for the LVDC network analysis are developed. Time-domain simulation models are implemented in the time-domain simulation environment PSCAD/EMTDC. The PSCAD/EMTDC models of the LVDC network are applied to the transient behaviour and power quality studies. The LVDC network power loss model is developed in a MATLAB environment and is capable of fast estimation of the network and component power losses. The model integrates analytical equations that describe the power loss mechanism of the network components with power flow calculations. For an LVDC network research platform, a monitoring and control software solution is developed. The solution is used to deliver measurement data for verification of the developed models and analysis of the modelling results. In the work, the power loss mechanism of the LVDC network components and its main dependencies are described. Energy loss distribution of the LVDC network components is presented. Power quality measurements and current spectra are provided and harmonic pollution on the DC network is analysed. The transient behaviour of the network is verified through time-domain simulations. DC capacitor guidelines for an LVDC power distribution network are introduced. The power loss analysis results show that one of the main optimisation targets for an LVDC power distribution network should be reduction of the no-load losses and efficiency improvement of converters at partial loads. Low-frequency spectra of the network voltages and currents are shown, and harmonic propagation is analysed. Power quality in the LVDC network point of common coupling (PCC) is discussed. Power quality standard requirements are shown to be met by the LVDC network. The network behaviour during transients is analysed by time-domain simulations. The network is shown to be transient stable during large-scale disturbances. Measurement results on the LVDC research platform proving this are presented in the work.

Segmentation of overlapping convex objects

This thesis presents a framework for segmentation of clustered overlapping convex objects. The proposed approach is based on a three-step framework in which the tasks of seed point extraction, contour evidence extraction, and contour estimation are addressed. The state-of-art techniques for each step were studied and evaluated using synthetic and real microscopic image data. According to obtained evaluation results, a method combining the best performers in each step was presented. In the proposed method, Fast Radial Symmetry transform, edge-to-marker association algorithm and ellipse fitting are employed for seed point extraction, contour evidence extraction and contour estimation respectively. Using synthetic and real image data, the proposed method was evaluated and compared with two competing methods and the results showed a promising improvement over the competing methods, with high segmentation and size distribution estimation accuracy.

Modelling of ultrasound assisted mixing in Newtonian and non-Newtonian liquids

The purpose of this work is to obtain a better understanding of behaviour of possible ultrasound appliance on fluid media mixing. The research is done in the regard to Newtonian and non-Newtonian fluids. The process of ultrasound appliance on liquids is modelled in COMSOL Multiphysics software. The influence of ultrasound using is introduced as waveform equation. Turbulence modelling is fulfilled by the k-ε model in Newtonian fluid. The modeling of ultrasound assisted mixing in non-Newtonian fluids is based on the power law. To verify modelling results two practical methods are used: Particle Image Velocimetry and measurements of mixing time. Particle Image Velocimetry allows capturing of velocity flow field continuously and presents detailed depiction of liquid dynamics. The second way of verification is the comparison of mixing time of homogeneity. Experimentally achievement of mixing time is done by conductivity measurements. In modelling part mixing time is achieved by special module of COMSOL Multiphysics – the transport of diluted species. Both practical and modelling parts show similar radial mechanism of fluid flow under ultrasound appliance – from the horn tip fluid moves to the bottom and along the walls goes back. Velocity profiles are similar in modelling and experimental part in the case of Newtonian fluid. In the case of non-Newtonian fluid velocity profiles do not agree. The development track of ultrasound-assisted mixing modelling is presented in the thesis.

Diclofenac plasma protein binding: PK-PD modelling in cardiac patients submitted to cardiopulmonary bypass

Twenty-four surgical patients of both sexes without cardiac, hepatic, renal or endocrine dysfunctions were divided into two groups: 10 cardiac surgical patients submitted to myocardial revascularization and cardiopulmonary bypass (CPB), 3 females and 7 males aged 65 ± 11 years, 74 ± 16 kg body weight, 166 ± 9 cm height and 1.80 ± 0.21 m2 body surface area (BSA), and control, 14 surgical patients not submitted to CPB, 11 female and 3 males aged 41 ± 14 years, 66 ± 14 kg body weight, 159 ± 9 cm height and 1.65 ± 0.16 m2 BSA (mean ± SD). Sodium diclofenac (1 mg/kg, im Voltaren 75® twice a day) was administered to patients in the Recovery Unit 48 h after surgery. Venous blood samples were collected during a period of 0-12 h and analgesia was measured by the visual analogue scale (VAS) during the same period. Plasma diclofenac levels were measured by high performance liquid chromatography. A two-compartment open model was applied to obtain the plasma decay curve and to estimate kinetic parameters. Plasma diclofenac protein binding decreased whereas free plasma diclofenac levels were increased five-fold in CPB patients. Data obtained for analgesia reported as the maximum effect (EMAX) were: 25% VAS (CPB) vs 10% VAS (control), P<0.05, median measured by the visual analogue scale where 100% is equivalent to the highest level of pain. To correlate the effect versus plasma diclofenac levels, the EMAX sigmoid model was applied. A prolongation of the mean residence time for maximum effect (MRTEMAX) was observed without any change in lag-time in CPB in spite of the reduced analgesia reported for these patients, during the time-dose interval. In conclusion, the extent of plasma diclofenac protein binding was influenced by CPB with clinically relevant kinetic-dynamic consequences

Application and development of numerical methods for the modelling of innovative gas cooled fission reactors

Innovative gas cooled reactors, such as the pebble bed reactor (PBR) and the gas cooled fast reactor (GFR) offer higher efficiency and new application areas for nuclear energy. Numerical methods were applied and developed to analyse the specific features of these reactor types with fully three dimensional calculation models. In the first part of this thesis, discrete element method (DEM) was used for a physically realistic modelling of the packing of fuel pebbles in PBR geometries and methods were developed for utilising the DEM results in subsequent reactor physics and thermal-hydraulics calculations. In the second part, the flow and heat transfer for a single gas cooled fuel rod of a GFR were investigated with computational fluid dynamics (CFD) methods. An in-house DEM implementation was validated and used for packing simulations, in which the effect of several parameters on the resulting average packing density was investigated. The restitution coefficient was found out to have the most significant effect. The results can be utilised in further work to obtain a pebble bed with a specific packing density. The packing structures of selected pebble beds were also analysed in detail and local variations in the packing density were observed, which should be taken into account especially in the reactor core thermal-hydraulic analyses. Two open source DEM codes were used to produce stochastic pebble bed configurations to add realism and improve the accuracy of criticality calculations performed with the Monte Carlo reactor physics code Serpent. Russian ASTRA criticality experiments were calculated. Pebble beds corresponding to the experimental specifications within measurement uncertainties were produced in DEM simulations and successfully exported into the subsequent reactor physics analysis. With the developed approach, two typical issues in Monte Carlo reactor physics calculations of pebble bed geometries were avoided. A novel method was developed and implemented as a MATLAB code to calculate porosities in the cells of a CFD calculation mesh constructed over a pebble bed obtained from DEM simulations. The code was further developed to distribute power and temperature data accurately between discrete based reactor physics and continuum based thermal-hydraulics models to enable coupled reactor core calculations. The developed method was also found useful for analysing sphere packings in general. CFD calculations were performed to investigate the pressure losses and heat transfer in three dimensional air cooled smooth and rib roughened rod geometries, housed inside a hexagonal flow channel representing a sub-channel of a single fuel rod of a GFR. The CFD geometry represented the test section of the L-STAR experimental facility at Karlsruhe Institute of Technology and the calculation results were compared to the corresponding experimental results. Knowledge was gained of the adequacy of various turbulence models and of the modelling requirements and issues related to the specific application. The obtained pressure loss results were in a relatively good agreement with the experimental data. Heat transfer in the smooth rod geometry was somewhat under predicted, which can partly be explained by unaccounted heat losses and uncertainties. In the rib roughened geometry heat transfer was severely under predicted by the used realisable k − epsilon turbulence model. An additional calculation with a v2 − f turbulence model showed significant improvement in the heat transfer results, which is most likely due to the better performance of the model in separated flow problems. Further investigations are suggested before using CFD to make conclusions of the heat transfer performance of rib roughened GFR fuel rod geometries. It is suggested that the viewpoints of numerical modelling are included in the planning of experiments to ease the challenging model construction and simulations and to avoid introducing additional sources of uncertainties. To facilitate the use of advanced calculation approaches, multi-physical aspects in experiments should also be considered and documented in a reasonable detail.