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.

Novel methods for the encapsulation of meglumine antimoniate into liposomes

The antimonial drug, meglumine antimoniate, was successfully encapsulated in dehydration-rehydration vesicles and in freeze-dried empty liposomes (FDELs). High encapsulation efficiencies (from 28 to 58%) and low weight ratios of lipids to encapsulated antimony (from 1:0.15 to 1:0.3) were achieved. These formulations, contrary to those obtained by conventional methods, can be stored as intermediate lyophilized forms and reconstituted just before use. The efficacy of FDEL-encapsulated meglumine antimoniate was evaluated in hamsters experimentally infected with Leishmania chagasi. A significant reduction of liver parasite burdens was observed in animals treated with this preparation, when compared to control animals treated with empty liposomes. In contrast, free meglumine antimoniate was found to be inefficient when administered at a comparable dose of antimony. This novel liposome-based meglumine antimoniate formulation appears to be promising as a pharmaceutical product for the treatment of visceral leishmaniasis.

The estimation of tuberculosis transmission parameters using ABC and MCMC methods

The aim of this work is to apply approximate Bayesian computation in combination with Marcov chain Monte Carlo methods in order to estimate the parameters of tuberculosis transmission. The methods are applied to San Francisco data and the results are compared with the outcomes of previous works. Moreover, a methodological idea with the aim to reduce computational time is also described. Despite the fact that this approach is proved to work in an appropriate way, further analysis is needed to understand and test its behaviour in different cases. Some related suggestions to its further enhancement are described in the corresponding chapter.

Radio-frequency-based measurement methods for bearing current analysis in induction motors

The most common reason for a low-voltage induction motor breakdown is a bearing failure. Along with the increasing popularity of modern frequency converters, bearing failures have become the most important motor fault type. Conditions in which bearing currents are likely to occur are generated as a side effect of fast du/dt switching transients. Once present, different types of bearing currents can accelerate the mechanical wear of bearings by causing deformation of metal parts in the bearing and degradation of the lubricating oil properties.The bearing current phenomena are well known, and several bearing current measurement and mitigation methods have been proposed. Nevertheless, in order to develop more feasible methods to measure and mitigate bearing currents, better knowledge of the phenomena is required. When mechanical wear is caused by bearing currents, the resulting aging impact has to be monitored and dealt with. Moreover, because of the stepwise aging mechanism, periodically executed condition monitoring measurements have been found ineffective. Thus, there is a need for feasible bearing current measurement methods that can be applied in parallel with the normal operation of series production drive systems. In order to reach the objectives of feasibility and applicability, nonintrusive measurement methods are preferred. In this doctoral dissertation, the characteristics and conditions of bearings that are related to the occurrence of different kinds of bearing currents are studied. Further, the study introduces some nonintrusive radio-frequency-signal-based approaches to detect and measure parameters that are associated with the accelerated bearing wear caused by bearing currents.

Study of Control Methods of a High-Frequency Cycloconverter for Domestic Induction Heating Appliances

A high-frequency cyclonverter acts as a direct ac-to-ac power converter circuit that does not require a diode bidge rectifier. Bridgeless topology makes it possible to remove forward voltage drop losses that are present in a diode bridge. In addition, the on-state losses can be reduced to 1.5 times the on-state resistance of switches in half-bridge operation of the cycloconverter. A high-frequency cycloconverter is reviewed and the charging effect of the dc-capacitors in ``back-to-back'' or synchronous mode operation operation is analyzed. In addition, a control method is introduced for regulating dc-voltage of the ac-side capacitors in synchronous operation mode. The controller regulates the dc-capacitors and prevents switches from reaching overvoltage level. This can be accomplished by variating phase-shift between the upper and the lower gate signals. By adding phase-shift between the gate signal pairs, the charge stored in the energy storage capacitors can be discharged through the resonant load and substantially, the output resonant current amplitude can be improved. The above goals are analyzed and illustrated with simulation. Theory is supported with practical measurements where the proposed control method is implemented in an FPGA device and tested with a high-frequency cycloconverter using super-junction power MOSFETs as switching devices.

Advanced analysis and design methods for preparative chromatographic separation processes

Preparative liquid chromatography is one of the most selective separation techniques in the fine chemical, pharmaceutical, and food industries. Several process concepts have been developed and applied for improving the performance of classical batch chromatography. The most powerful approaches include various single-column recycling schemes, counter-current and cross-current multi-column setups, and hybrid processes where chromatography is coupled with other unit operations such as crystallization, chemical reactor, and/or solvent removal unit. To fully utilize the potential of stand-alone and integrated chromatographic processes, efficient methods for selecting the best process alternative as well as optimal operating conditions are needed. In this thesis, a unified method is developed for analysis and design of the following singlecolumn fixed bed processes and corresponding cross-current schemes: (1) batch chromatography, (2) batch chromatography with an integrated solvent removal unit, (3) mixed-recycle steady state recycling chromatography (SSR), and (4) mixed-recycle steady state recycling chromatography with solvent removal from fresh feed, recycle fraction, or column feed (SSR–SR). The method is based on the equilibrium theory of chromatography with an assumption of negligible mass transfer resistance and axial dispersion. The design criteria are given in general, dimensionless form that is formally analogous to that applied widely in the so called triangle theory of counter-current multi-column chromatography. Analytical design equations are derived for binary systems that follow competitive Langmuir adsorption isotherm model. For this purpose, the existing analytic solution of the ideal model of chromatography for binary Langmuir mixtures is completed by deriving missing explicit equations for the height and location of the pure first component shock in the case of a small feed pulse. It is thus shown that the entire chromatographic cycle at the column outlet can be expressed in closed-form. The developed design method allows predicting the feasible range of operating parameters that lead to desired product purities. It can be applied for the calculation of first estimates of optimal operating conditions, the analysis of process robustness, and the early-stage evaluation of different process alternatives. The design method is utilized to analyse the possibility to enhance the performance of conventional SSR chromatography by integrating it with a solvent removal unit. It is shown that the amount of fresh feed processed during a chromatographic cycle and thus the productivity of SSR process can be improved by removing solvent. The maximum solvent removal capacity depends on the location of the solvent removal unit and the physical solvent removal constraints, such as solubility, viscosity, and/or osmotic pressure limits. Usually, the most flexible option is to remove solvent from the column feed. Applicability of the equilibrium design for real, non-ideal separation problems is evaluated by means of numerical simulations. Due to assumption of infinite column efficiency, the developed design method is most applicable for high performance systems where thermodynamic effects are predominant, while significant deviations are observed under highly non-ideal conditions. The findings based on the equilibrium theory are applied to develop a shortcut approach for the design of chromatographic separation processes under strongly non-ideal conditions with significant dispersive effects. The method is based on a simple procedure applied to a single conventional chromatogram. Applicability of the approach for the design of batch and counter-current simulated moving bed processes is evaluated with case studies. It is shown that the shortcut approach works the better the higher the column efficiency and the lower the purity constraints are.

Comparison of the quantitative competitive and semiquantitative RT-PCR methods for the determination of interferon-gamma mRNA levels in AIDS-free HIV-infected individuals

IFN-gamma mRNA expression was evaluated in nonstimulated peripheral blood mononuclear cells (PBMC) of HIV-infected and seronegative individuals using quantitative competitive and semiquantitative RT-PCR and the sensitivity of these methods was compared. A significant correlation was found between quantitative competitive and semiquantitative RT-PCR in samples of both HIV-seronegative (P = 0.004) and HIV-infected individuals (P = 0.0004). PBMC from HIV-infected individuals presented a remarkable increase of IFN-gamma mRNA expression, as determined by both types of RT-PCR methods. Semiquantitative RT-PCR even without an internal standard is also acceptable for measuring cytokine mRNA expression, but less reliable if small amounts are quantified. Moreover, we found that increased IFN-gammamRNA expression is independent of CD4+ cell count in AIDS-free HIV-infected patients.

Reporting on methods of subgroup analysis in clinical trials: a survey of four scientific journals

Results of subgroup analysis (SA) reported in randomized clinical trials (RCT) cannot be adequately interpreted without information about the methods used in the study design and the data analysis. Our aim was to show how often inaccurate or incomplete reports occur. First, we selected eight methodological aspects of SA on the basis of their importance to a reader in determining the confidence that should be placed in the author's conclusions regarding such analysis. Then, we reviewed the current practice of reporting these methodological aspects of SA in clinical trials in four leading journals, i.e., the New England Journal of Medicine, the Journal of the American Medical Association, the Lancet, and the American Journal of Public Health. Eight consecutive reports from each journal published after July 1, 1998 were included. Of the 32 trials surveyed, 17 (53%) had at least one SA. Overall, the proportion of RCT reporting a particular methodological aspect ranged from 23 to 94%. Information on whether the SA preceded/followed the analysis was reported in only 7 (41%) of the studies. Of the total possible number of items to be reported, NEJM, JAMA, Lancet and AJPH clearly mentioned 59, 67, 58 and 72%, respectively. We conclude that current reporting of SA in RCT is incomplete and inaccurate. The results of such SA may have harmful effects on treatment recommendations if accepted without judicious scrutiny. We recommend that editors improve the reporting of SA in RCT by giving authors a list of the important items to be reported.

Comparison of methods for urinary albumin determination in patients with type 1 diabetes

We tested the correlation of the albumin-to-creatinine ratio (A/C) in an early-morning urine sample, measured with a commercial kit (DCA 2000®), with the conventional immunoturbidimetric determination in the laboratory and with overnight albumin excretion rate (reference method). Fifty-five type 1 diabetic adolescents had their first-morning urine collected on the 1st and 8th day of the period. Urinary albumin and creatinine were determined immediately using the DCA 2000® kit. Samples were also stored for laboratory analysis. To evaluate the correlation between early-morning urinary A/C ratio and overnight albumin excretion rate, 16 subjects had a timed overnight urine collection. A/C ratios determined with the DCA 2000® kit and by the laboratory method were 13.1 ± 20.5 and 20.4 ± 46.3 mg/g, respectively. A/C results by both methods proved to be strongly correlated (r = 0.98, P<0.001). DCA 2000®-determined A/C showed 50% sensitivity and 100% specificity when compared to the reference method. Spot urinary A/C of the subset of 16 subjects significantly correlated with their overnight albumin excretion rate (r = 0.98, P<0.001). Intraindividual variation ranged from 17 to 32% and from 9 to 63% for A/C and overnight albumin excretion rate, respectively. In conclusion, an early-morning specimen should be used instead of timed overnight urine and the A/C ratio is an accurate, reliable and easily determined parameter for the screening of diabetic nephropathy. Immediate measurement of the A/C ratio is feasible using the DCA 2000® kit. Intraindividual variability indicates the need for repeated determinations to confirm microalbuminuria and the diagnosis of incipient diabetic nephropathy.

Clinical features of panic patients sensitive to hyperventilation or breath-holding methods for inducing panic attacks

Our aim was to compare the clinical features of panic disorder (PD) patients sensitive to hyperventilation or breath-holding methods of inducing panic attacks. Eighty-five PD patients were submitted to both a hyperventilation challenge test and a breath-holding test. They were asked to hyperventilate (30 breaths/min) for 4 min and a week later to hold their breath for as long as possible, four times with a 2-min interval. Anxiety scales were applied before and after the tests. We selected the patients who responded with a panic attack to just one of the tests, i.e., those who had a panic attack after hyperventilating (HPA, N = 24, 16 females, 8 males, mean age ± SD = 38.5 ± 12.7 years) and those who had a panic attack after breath holding (BHPA, N = 20, 11 females, 9 males, mean age ± SD = 42.1 ± 10.6 years). Both groups had similar (chi² = 1.28, d.f. = 1, P = 0.672) respiratory symptoms (fear of dying, chest/pain disconfort, shortness of breath, paresthesias, and feelings of choking) during a panic attack. The criteria of Briggs et al. [British Journal of Psychiatry, 1993; 163: 201-209] for respiratory PD subtype were fulfilled by 18 (75.0%) HPA patients and by 14 (70.0%) BHPA patients. The HPA group had a later onset of the disease compared to BHPA patients (37.9 ± 11.0 vs 21.3 ± 12.9 years old, Mann-Whitney, P < 0.001), and had a higher family prevalence of PD (70.8 vs 25.0%, chi² = 19.65, d.f. = 1, P = 0.041). Our data suggest that these two groups - HPA and BHPA patients - may be specific subtypes of PD.

Statistical segmentation methods and color variance analysis of retinal images

In this research, the effectiveness of Naive Bayes and Gaussian Mixture Models classifiers on segmenting exudates in retinal images is studied and the results are evaluated with metrics commonly used in medical imaging. Also, a color variation analysis of retinal images is carried out to find how effectively can retinal images be segmented using only the color information of the pixels.

Novel virtual environment and real-time simulation based methods for improving life-cycle efficiency of non-road mobile machinery

Virtual environments and real-time simulators (VERS) are becoming more and more important tools in research and development (R&D) process of non-road mobile machinery (NRMM). The virtual prototyping techniques enable faster and more cost-efficient development of machines compared to use of real life prototypes. High energy efficiency has become an important topic in the world of NRMM because of environmental and economic demands. The objective of this thesis is to develop VERS based methods for research and development of NRMM. A process using VERS for assessing effects of human operators on the life-cycle efficiency of NRMM was developed. Human in the loop simulations are ran using an underground mining loader to study the developed process. The simulations were ran in the virtual environment of the Laboratory of Intelligent Machines of Lappeenranta University of Technology. A physically adequate real-time simulation model of NRMM was shown to be reliable and cost effective in testing of hardware components by the means of hardware-in-the-loop (HIL) simulations. A control interface connecting integrated electro-hydraulic energy converter (IEHEC) with virtual simulation model of log crane was developed. IEHEC consists of a hydraulic pump-motor and an integrated electrical permanent magnet synchronous motorgenerator. The results show that state of the art real-time NRMM simulators are capable to solve factors related to energy consumption and productivity of the NRMM. A significant variation between the test drivers is found. The results show that VERS can be used for assessing human effects on the life-cycle efficiency of NRMM. HIL simulation responses compared to that achieved with conventional simulation method demonstrate the advances and drawbacks of various possible interfaces between the simulator and hardware part of the system under study. Novel ideas for arranging the interface are successfully tested and compared with the more traditional one. The proposed process for assessing the effects of operators on the life-cycle efficiency will be applied for wider group of operators in the future. Driving styles of the operators can be analysed statistically from sufficient large result data. The statistical analysis can find the most life-cycle efficient driving style for the specific environment and machinery. The proposed control interface for HIL simulation need to be further studied. The robustness and the adaptation of the interface in different situations must be verified. The future work will also include studying the suitability of the IEHEC for different working machines using the proposed HIL simulation method.

Balancing the quotation process by LEAN methods

This thesis studies the possibility to use lean tools and methods in a quotation process which is carried out in an office environment. The aim of the study was to find out and test the relevant lean tools and methods which can help to balance and standardize the quotation process, and reduce the variance in quotation lead times and in quality. Seminal works, researches and guide books related to the topic were used as the basis for the theory development. Based on the literature review and the case company’s own lean experience, the applicable lean tools and methods were selected to be tested by a sales support team. Leveling production, by product categorization and value stream mapping, was a key method to be used to balance the quotation process. 5S method was started concurrently for standardizing the work. Results of the testing period showed that lean tools and methods are applicable in office process and selected tools and methods helped to balance and standardize the quotation process. Case company’s sales support team decided to implement new lean based quotation process model.

Lean Methods and Digital Technologies in Healtchare

The experiences of several healthcare organizations were considered to distinguish the most frequently used lean tools, the success and failure factors, and the obstacles that may appear while implementing lean. As a result, two approaches to “go lean” were defined, and analyzed from the prospective of the applicability to healthcare processes. Industrialization of healthcare was studied, and the most promising digital technology tools to improve healthcare process were highlighted. Finally, the analysis of healthcare challenges and feasible ways to address them was conducted and presented as the main result of this work. The possible ways of implementation of the findings and limitations were described in the conclusion.