Modelling of combustion and sorbent reactions in three-dimensional flow environment of a circulating fluidized bed furnace

This thesis presents a three-dimensional, semi-empirical, steady state model for simulating the combustion, gasification, and formation of emissions in circulating fluidized bed (CFB) processes. In a large-scale CFB furnace, the local feeding of fuel, air, and other input materials, as well as the limited mixing rate of different reactants produce inhomogeneous process conditions. To simulate the real conditions, the furnace should be modelled three-dimensionally or the three-dimensional effects should be taken into account. The only available methods for simulating the large CFB furnaces three-dimensionally are semi-empirical models, which apply a relatively coarse calculation mesh and a combination of fundamental conservation equations, theoretical models and empirical correlations. The number of such models is extremely small. The main objective of this work was to achieve a model which can be applied to calculating industrial scale CFB boilers and which can simulate all the essential sub-phenomena: fluid dynamics, reactions, the attrition of particles, and heat transfer. The core of the work was to develop the model frame and the required sub-models for determining the combustion and sorbent reactions. The objective was reached, and the developed model was successfully used for studying various industrial scale CFB boilers combusting different types of fuel. The model for sorbent reactions, which includes the main reactions for calcitic limestones, was applied for studying the new possible phenomena occurring in the oxygen-fired combustion. The presented combustion and sorbent models and principles can be utilized in other model approaches as well, including other empirical and semi-empirical model approaches, and CFD based simulations. The main achievement is the overall model frame which can be utilized for the further development and testing of new sub-models and theories, and for concentrating the knowledge gathered from the experimental work carried out at bench scale, pilot scale and industrial scale apparatus, and from the computational work performed by other modelling methods.

Design Aspects of Megawatt-Range Direct-Driven Permanent MagnetWind Generators

Direct-driven permanent magnet synchronous generator is one of the most promising topologies for megawatt-range wind power applications. The rotational speed of the direct-driven generator is very low compared with the traditional electrical machines. The low rotational speed requires high torque to produce megawatt-range power. The special features of the direct-driven generators caused by the low speed and high torque are discussed in this doctoral thesis. Low speed and high torque set high demands on the torque quality. The cogging torque and the load torque ripple must be as low as possible to prevent mechanical failures. In this doctoral thesis, various methods to improve the torque quality are compared with each other. The rotor surface shaping, magnet skew, magnet shaping, and the asymmetrical placement of magnets and stator slots are studied not only by means of torque quality, but also the effects on the electromagnetic performance and manufacturability of the machine are discussed. The heat transfer of the direct-driven generator must be designed to handle the copper losses of the stator winding carrying high current density and to keep the temperature of the magnets low enough. The cooling system of the direct-driven generator applying the doubly radial air cooling with numerous radial cooling ducts was modeled with a lumped-parameter-based thermal network. The performance of the cooling system was discussed during the steady and transient states. The effect of the number and width of radial cooling ducts was explored. The large number of radial cooling ducts drastically increases the impact of the stack end area effects, because the stator stack consists of numerous substacks. The effects of the radial cooling ducts on the effective axial length of the machine were studied by analyzing the crosssection of the machine in the axial direction. The method to compensate the magnet end area leakage was considered. The effect of the cooling ducts and the stack end area effects on the no-load voltages and inductances of the machine were explored by using numerical analysis tools based on the three-dimensional finite element method. The electrical efficiency of the permanent magnet machine with different control methods was estimated analytically over the whole speed and torque range. The electrical efficiencies achieved with the most common control methods were compared with each other. The stator voltage increase caused by the armature reaction was analyzed. The effect of inductance saturation as a function of load current was implemented to the analytical efficiency calculation.

CFD Modelling of Direct Contact Condensation in Suppression Pools by Applying Condensation Models of Separated Flow

The condensation rate has to be high in the safety pressure suppression pool systems of Boiling Water Reactors (BWR) in order to fulfill their safety function. The phenomena due to such a high direct contact condensation (DCC) rate turn out to be very challenging to be analysed either with experiments or numerical simulations. In this thesis, the suppression pool experiments carried out in the POOLEX facility of Lappeenranta University of Technology were simulated. Two different condensation modes were modelled by using the 2-phase CFD codes NEPTUNE CFD and TransAT. The DCC models applied were the typical ones to be used for separated flows in channels, and their applicability to the rapidly condensing flow in the condensation pool context had not been tested earlier. A low Reynolds number case was the first to be simulated. The POOLEX experiment STB-31 was operated near the conditions between the ’quasi-steady oscillatory interface condensation’ mode and the ’condensation within the blowdown pipe’ mode. The condensation models of Lakehal et al. and Coste & Lavi´eville predicted the condensation rate quite accurately, while the other tested ones overestimated it. It was possible to get the direct phase change solution to settle near to the measured values, but a very high resolution of calculation grid was needed. Secondly, a high Reynolds number case corresponding to the ’chugging’ mode was simulated. The POOLEX experiment STB-28 was chosen, because various standard and highspeed video samples of bubbles were recorded during it. In order to extract numerical information from the video material, a pattern recognition procedure was programmed. The bubble size distributions and the frequencies of chugging were calculated with this procedure. With the statistical data of the bubble sizes and temporal data of the bubble/jet appearance, it was possible to compare the condensation rates between the experiment and the CFD simulations. In the chugging simulations, a spherically curvilinear calculation grid at the blowdown pipe exit improved the convergence and decreased the required cell count. The compressible flow solver with complete steam-tables was beneficial for the numerical success of the simulations. The Hughes-Duffey model and, to some extent, the Coste & Lavi´eville model produced realistic chugging behavior. The initial level of the steam/water interface was an important factor to determine the initiation of the chugging. If the interface was initialized with a water level high enough inside the blowdown pipe, the vigorous penetration of a water plug into the pool created a turbulent wake which invoked the chugging that was self-sustaining. A 3D simulation with a suitable DCC model produced qualitatively very realistic shapes of the chugging bubbles and jets. The comparative FFT analysis of the bubble size data and the pool bottom pressure data gave useful information to distinguish the eigenmodes of chugging, bubbling, and pool structure oscillations.

Electrically driven fertilizer applicator controlled by fuzzy logic

The present study shows the development, simulation and actual implementation of a closed-loop controller based on fuzzy logic that is able to regulate and standardize the mass flow of a helical fertilizer applicator. The control algorithm was developed using MATLAB's Fuzzy Logic Toolbox. Both open and closed-loop simulations of the controller were performed in MATLAB's Simulink environment. The instantaneous deviation of the mass flow from the set point (SP), its derivative, the equipment´s translation velocity and acceleration were all used as input signals for the controller, whereas the voltage of the applicator's DC electric motor (DCEM) was driven by the controller as output signal. Calibration and validation of the rules and membership functions of the fuzzy logic were accomplished in the computer simulation phase, taking into account the system's response to SP changes. The mass flow variation coefficient, measured in experimental tests, ranged from 6.32 to 13.18%. The steady state error fell between -0.72 and 0.13g s-1 and the recorded average rise time of the system was 0.38 s. The implemented controller was able to both damp the oscillations in mass flow that are characteristic of helical fertilizer applicators, and to effectively respond to SP variations.

Effect of companion ion on nitrate displacement through transport parameters analysis

Nitrate is the main form of nitrogen associated with water contamination; the high mobility of this species in soil justifies the concern regarding nitrogen management in agricultural soils. Therefore, the objective of this research was to assess the effect of companion cation on nitrate displacement, by analyzing nitrate transport parameters through Breakthrough Curves (BTCs) and their settings made by numerical model (STANMOD). The experiment was carried out in the Soil and Water Quality Laboratory of the Department of Biosystems Engineering, "Luiz de Queiroz" College of Agriculture in Piracicaba (SP), Brazil. It was performed using saturated soil columns in steady-state flow condition, in which two different sources of inorganic nitrate Ca(NO3)2 and NH4NO3 were applied at a concentration of 50 mg L-1 NO3-. Each column was filled with either a Red-Yellow Oxisol (S1) or an Alfisol (S2). Results are indicative that the companion ion had no effect on nitrate displacement. However, nitrate transport was influenced by soil texture, particle aggregation, solution speed in soil and organic matter presence. Nitrate mobility was higher in the Alfisol (S2).

Fundamentals of heat transfer

Following over 170+ pages and additional appendixes are formed based on content of Course: Fundamentals of Heat Transfer. Mainly this summarizes relevant parts on Book of Fundamentals of Heat and Mass Transfer (Incropera), but also other references introducing the same concepts are included. Student’s point of view has been consideredwith following highlights: (1) Relevant topics are presented in a nutshell to provide fast digestion of principles of heat transfer. (2) Appendixes include terminology dictionary. (3) Totally 22 illustrating examples are connecting theory to practical applications and quantifying heat transfer to understandable forms as: temperatures, heat transfer rates, heat fluxes, resistances and etc. (4) Most important Learning outcomes are presented for each topic separately. The Book, Fundamentals of Heat and Mass Transfer (Incropera), is certainly recommended for those going beyond basic knowledge of heat transfer. Lecture Notes consists of four primary content-wise objectives: (1) Give understanding to physical mechanisms of heat transfer, (2)Present basic concepts and terminology relevant for conduction, convection and radiation (3) Introduce thermal performance analysis methods for steady state and transient conduction systems. (4) Provide fast-to-digest phenomenological understanding required for basic design of thermal models

PWR PACTEL -koelaitteiston ylätilan mallinnusvaihtoehtojen tarkastelu APROS-ohjelmistolla

Diplomityössä selvitetään PWR PACTEL -koelaitteiston APROS-mallin ylätilan toimi¬vuut¬ta ja tutkitaan tapoja parantaa mallin simulointitarkkuutta. Työssä on esitelty PWR PACTEL -koelaitteiston sekä APROS-simulointiohjelman pääpiirteet ja tutustuttu simu¬loin¬ti¬ohjelmien kelpoistukseen. Simulointiosiossa tarkastellaan yksinkertaisen ja muokatun ylätilamallin toimivuutta ver¬taa¬mal¬la simulointia koelaitteistolla tehtyjen kokeiden mittaustuloksiin. Lisäksi tutkitaan nel¬¬jää tapaa parantaa ylätilamallin toimivuutta: virtausvastusten säätö, lämpöhäviöiden sä䬬tö, rakenteelliset muutokset ja laskentakorrelaatioiden vaihtaminen. Simulointituloksia ver¬¬ra¬¬taan kahteen koelaitteistolla suoritettuun kokeeseen, joista toinen on tasapainotilankoe ja toinen dynaamisen tilan koe. Työssä havaittiin, että nykyisin käytössä olevaa mallia on kehitetty jo varsin paljon, ver¬rat¬tu¬na ns. yksinkertaiseen malliin. Testatuilla parannusvaihtoehdoilla saatiin hieman paran¬net¬tua ylätilan lämpötilojen simulointia tasapainotilassa. Dynaamisen tilan osalta havaittiin pien¬tä parannusta mallin toiminnassa ylätilamallin rakennemuutosten jälkeen. Tulosten poh¬jalta arvioitiin, että mallia on mahdollista kehittää muuttamalla ylätilamallin raken¬net¬ta.

Acidic dissolution of iron oxides and regeneration of a ceramic filter medium

The dewatering of iron ore concentrates requires large capacity in addition to producing a cake with low moisture content. Such large processes are commonly energy intensive and means to lower the specific energy consumption are needed. Ceramic capillary action disc filters incorporate a novel filter medium enabling the harnessing of capillary action, which results in decreased energy consumption in comparison to traditional filtration technologies. As another benefit, the filter medium is mechanically and chemically more durable than, for example, filter cloths and can, thus, withstand harsh operating conditions and possible regeneration better than other types of filter media. In iron ore dewatering, the regeneration of the filter medium is done through a combination of several techniques: (1) backwashing, (2) ultrasonic cleaning, and (3) acid regeneration. Although it is commonly acknowledged that the filter medium is affected by slurry particles and extraneous compounds, published research, especially in the field of dewatering of mineral concentrates, is scarce. Whereas the regenerative effect of backwashing and ultrasound are more or less mechanical, regeneration with acids is based on chemistry. The chemistry behind the acid regeneration is, naturally, dissolution. The dissolution of iron oxide particles has been extensively studied over several decades but those studies may not necessarily be directly applicable in the regeneration of the filter medium which has undergone interactions with the slurry components. The aim of this thesis was to investigate if free particle dissolution indeed correlates with the regeneration of the filter medium. For this purpose, both free particle dissolution and dissolution of surface adhered particles were studied. The focus was on acidic dissolution of iron oxide particles and on the study of the ceramic filter medium used in the dewatering of iron ore concentrates. The free particle dissolution experiments show that the solubility of synthetic fine grained iron oxide particles in oxalic acid could be explained through linear models accounting for the effects of temperature and acid concentration, whereas the dissolution of a natural magnetite is not so easily explained by such models. In addition, the kinetic experiments performed both support and contradict the work of previous authors: the suitable kinetic model here supports previous research suggesting solid state reduction to be the reaction mechanism of hematite dissolution but the formation of a stable iron oxalate is not supported by the results of this research. Several other dissolution mechanisms have also been suggested for iron oxide dissolution in oxalic acid, indicating that the details of oxalate promoted reductive dissolution are not yet agreed and, in this respect, this research offers added value to the community. The results of the regeneration experiments with the ceramic filter media show that oxalic acid is highly effective in removing iron oxide particles from the surface of the filter medium. The dissolution of those particles did not, however, exhibit the expected behaviour, i.e. complete dissolution. The results of this thesis show that although the regeneration of the ceramic filter medium with acids incorporates the dissolution of slurry particles from the surface of the filter medium, the regeneration cannot be assessed purely based upon free particle dissolution. A steady state, dependent on temperature and on the acid concentration, was observed in the dissolution of particles from the surface even though the limit of solubility of free iron oxide particles had not been reached. Both the regeneration capacity and efficiency, with regards to the removal of iron oxide particles, was found to be temperature dependent, but was not affected by the acid concentration. This observation further suggests that the removal of the surface adhered particles does not follow the dissolution of free particles, which do exhibit a dependency on the acid concentration. In addition, changes in the permeability and in the pore structure of the filter medium were still observed after the bulk concentration of dissolved iron had reached a steady state. Consequently, the regeneration of the filter medium continued after the dissolution of particles from the surface had ceased. This observation suggests that internal changes take place at the final stages of regeneration. The regeneration process could, in theory, be divided into two, possibly overlapping, stages: (1) dissolution of surface-adhered particles, and (2) dissolution of extraneous compounds from within the pore structure. In addition to the fundamental knowledge generated during this thesis, tools to assess the effects of parameters on the regeneration of the ceramic filter medium are needed. It has become clear that the same tools used to estimate the dissolution of free particles cannot be used to estimate the regeneration of a filter medium unless only a robust characterisation of the order of regeneration efficiency is needed.

Accuracy Analysis of Frequency Converter Estimates in Diagnostic Applications

Transportation of fluids is one of the most common and energy intensive processes in the industrial and HVAC sectors. Pumping systems are frequently subject to engineering malpractice when dimensioned, which can lead to poor operational efficiency. Moreover, pump monitoring requires dedicated measuring equipment, which imply costly investments. Inefficient pump operation and improper maintenance can increase energy costs substantially and even lead to pump failure. A centrifugal pump is commonly driven by an induction motor. Driving the induction motor with a frequency converter can diminish energy consumption in pump drives and provide better control of a process. In addition, induction machine signals can also be estimated by modern frequency converters, dispensing with the use of sensors. If the estimates are accurate enough, a pump can be modelled and integrated into the frequency converter control scheme. This can open the possibility of joint motor and pump monitoring and diagnostics, thereby allowing the detection of reliability-reducing operating states that can lead to additional maintenance costs. The goal of this work is to study the accuracy of rotational speed, torque and shaft power estimates calculated by a frequency converter. Laboratory tests were performed in order to observe estimate behaviour in both steady-state and transient operation. An induction machine driven by a vector-controlled frequency converter, coupled with another induction machine acting as load was used in the tests. The estimated quantities were obtained through the frequency converter’s Trend Recorder software. A high-precision, HBM T12 torque-speed transducer was used to measure the actual values of the aforementioned variables. The effect of the flux optimization energy saving feature on the estimate quality was also studied. A processing function was developed in MATLAB for comparison of the obtained data. The obtained results confirm the suitability of this particular converter to provide accurate enough estimates for pumping applications.

The Effects of Anesthetic Induced Loss of Consciousness on Quantitative Electroen Cephalogram, and Bispectral and Spectral Entropy Indices. Studies on Healthy Male

Changes in the electroencephalography (EEG) signal have been used to study the effects of anesthetic agents on the brain function. Several commercial EEG based anesthesia depth monitors have been developed to measure the level of the hypnotic component of anesthesia. Specific anesthetic related changes can be seen in the EEG, but still it remains difficult to determine whether the subject is consciousness or not during anesthesia. EEG reactivity to external stimuli may be seen in unconsciousness subjects, in anesthesia or even in coma. Changes in regional cerebral blood flow, which can be measured with positron emission tomography (PET), can be used as a surrogate for changes in neuronal activity. The aim of this study was to investigate the effects of dexmedetomidine, propofol, sevoflurane and xenon on the EEG and the behavior of two commercial anesthesia depth monitors, Bispectral Index (BIS) and Entropy. Slowly escalating drug concentrations were used with dexmedetomidine, propofol and sevoflurane. EEG reactivity at clinically determined similar level of consciousness was studied and the performance of BIS and Entropy in differentiating consciousness form unconsciousness was evaluated. Changes in brain activity during emergence from dexmedetomidine and propofol induced unconsciousness were studied using PET imaging. Additionally, the effects of normobaric hyperoxia, induced during denitrogenation prior to xenon anesthesia induction, on the EEG were studied. Dexmedetomidine and propofol caused increases in the low frequency, high amplitude (delta 0.5-4 Hz and theta 4.1-8 Hz) EEG activity during stepwise increased drug concentrations from the awake state to unconsciousness. With sevoflurane, an increase in delta activity was also seen, and an increase in alpha- slow beta (8.1-15 Hz) band power was seen in both propofol and sevoflurane. EEG reactivity to a verbal command in the unconsciousness state was best retained with propofol, and almost disappeared with sevoflurane. The ability of BIS and Entropy to differentiate consciousness from unconsciousness was poor. At the emergence from dexmedetomidine and propofol induced unconsciousness, activation was detected in deep brain structures, but not within the cortex. In xenon anesthesia, EEG band powers increased in delta, theta and alpha (8-12Hz) frequencies. In steady state xenon anesthesia, BIS and Entropy indices were low and these monitors seemed to work well in xenon anesthesia. Normobaric hyperoxia alone did not cause changes in the EEG. All of these results are based on studies in healthy volunteers and their application to clinical practice should be considered carefully.

Effect of exposure to sublethal concentrations of sodium cyanide on the carbohydrate metabolism of the Indian Major Carp Labeo rohita (Hamilton, 1822)

Experiments were designed to study in-vivo effects of sodium cyanide on biochemical endpoints in the freshwater fish Labeo rohita. Fish were exposed to two sublethal concentrations (0.106 and 0.064mg/L) for a period of 15 days. Levels of glycogen, pyruvate, lactate and the enzymatic activities of lactate dehydrogenase (LDH), succinate dehydrogenase (SDH), glucose-6-phosphate dehydrogenase (G6PDH), phosphorylase, alkaline phosphatase (ALP), acid phosphatase (AcP) were assessed in different tissues (liver, muscle and gills). Result indicated a steady decrease in glycogen, pyruvate, SDH, ALP and AcP activity with a concomitant increase in the lactate, phosphorylase, LDH and G6PD activity in all selected tissues. The alterations in all the above biochemical parameters were significantly (p<0.05) time and dose dependent. In all the above parameters, liver pointing out the intensity of cyanide intoxication compare to muscle and gills. Study revealed change in the metabolic energy by means of altered metabolic profile of the fish. Further, these observations indicated that even sublethal concentrations of sodium cyanide might not be fully devoid of deleterious influence on metabolism in L. rohita.

Microbusiness: A new social business model to alleviate poverty in developing countries?

Poverty alleviation views have shifted from seeing the poor as victims or as potential consumers, to seeing them as gainers. Social businesses include microfinancing and microfranchising, which engage people at the bottom of the pyramid using business instead of charity. There are, however, social business firms that do not fit to the existing social business model theory. These firms provide markets to poor producers and mix traditional, local craftsmanship with western design. Social business models evolve faster than the academic literature can study them. This study contributes to filling this gap. The purpose of this Master’s thesis is to develop the concept of social business as poverty alleviation method in developing countries. It also aims; 1) to describe the means for poverty alleviation in developing countries; 2) to introduce microbusiness as a social business model; and 3) to examine the challenges of microbusinesses. Qualitative case study is used as a research strategy and theme interviews as a data collecting method. The empirical data is gathered from four interviews of Finnish or Finnish-owned firms that employ microbusiness – Mifuko, Tensira, Mangomaa and Tikau – and this is supported with secondary data including articles on case companies. The results show that microbusiness is a valid new social business model that aims at poverty alleviation by engaging the poor at the bottom of the pyramid. It is possible to map the value proposition, value constellation, and economic and social profit equations of the case firms. Two major types of firms emerge from the results; the first consists of design-oriented firms that emphasize the quality and design of the products, and the second consists of bazaar-like firms whose product portfolio is less sophisticated and who promote more the stories of the products – not the design. All microbusiness firms provide markets, promote traditional handicrafts, form close relationships to their producers, and aim at enhancing lives through their businesses. The attitudes towards social businesses are sometimes negative, but this is changing for the better. In conclusion, microbusiness answers to two different needs at the same time – consumers’ needs for ethical products and the social needs of the producers – but the social need is the ultimate reason why the entrepreneurs started business. Microbusiness continues as a poverty alleviation tool that sees the poor as gainers; by providing them steady employment, microbusiness increases the poor’s self-esteem and enables them for a better living. Academic literature has not been able to offer enough alternative business models to cover all social businesses; the current study contributes to this by concluding that microbusiness is another social business model.

Regulation of photosynthetic electron transfer in plant chloroplasts

This thesis focuses on the molecular mechanisms regulating the photosynthetic electron transfer reactions upon changes in light intensity. To investigate these mechanisms, I used mutants of the model plant Arabidopsis thaliana impaired in various aspects of regulation of the photosynthetic light reactions. These included mutants of photosystem II (PSII) and light harvesting complex II (LHCII) phosphorylation (stn7 and stn8), mutants of energy-dependent non-photochemical quenching (NPQ) (npq1 and npq4) and of regulation of photosynthetic electron transfer (pgr5). All of these processes have been extensively investigated during the past decades, mainly on plants growing under steady-state conditions, and therefore many aspects of acclimation processes may have been neglected. In this study, plants were grown under fluctuating light, i.e. the alternation of low and high intensities of light, in order to maximally challenge the photosynthetic regulatory mechanisms. In pgr5 and stn7 mutants, the growth in fluctuating light condition mainly damaged PSI while PSII was rather unaffected. It is shown that the PGR5 protein regulates the linear electron transfer: it is essential for the induction of transthylakoid ΔpH that, in turn, activates energy-dependent NPQ and downregulates the activity of cytochrome b6f. This regulation was shown to be essential for the photoprotection of PSI under fluctuations in light intensity. The stn7 mutants were able to acclimate under constant growth light conditions by modulating the PSII/PSI ratio, while under fluctuating growth light they failed in implementing this acclimation strategy. LHCII phosphorylation ensures the balance of the excitation energy distribution between PSII and PSI by increasing the probability for excitons to be trapped by PSI. LHCII can be phosphorylated over all of the thylakoid membrane (grana cores as well as stroma lamellae) and when phosphorylated it constitutes a common antenna for PSII and PSI. Moreover, LHCII was shown to work as a functional bridge that allows the energy transfer between PSII units in grana cores and between PSII and PSI centers in grana margins. Consequently, PSI can function as a quencher of excitation energy. Eventually, the LHCII phosphorylation, NPQ and the photosynthetic control of linear electron transfer via cytochrome b6f work in concert to maintain the redox poise of the electron transfer chain. This is a prerequisite for successful plant growth upon changing natural light conditions, both in short- and long-term.

Experimental analysis of the flow between stay and guide vanes of a pump-turbine in pumping mode

An experimental study was conducted in a pump-turbine model in pumping mode, in order to characterize the flow field structure in the region between stay and guide vanes, using mainly the laser-Doppler anemometry in a two-color and back-scattered light-based system. The structure of the steady and unsteady flow was analyzed. The measurements were performed at three operation points. The obtained data provide appropriate boundary conditions and a good base of validation for numerical codes, and for the understanding of main loss mechanisms of this complex flow.

Application of a non isotropic turbulence model to stable atmospheric flows and dispersion over 3D topography

A non isotropic turbulence model is extended and applied to three dimensional stably stratified flows and dispersion calculations. The model is derived from the algebraic stress model (including wall proximity effects), but it retains the simplicity of the "eddy viscosity" concept of first order models. The "modified k-epsilon" is implemented in a three dimensional numerical code. Once the flow is resolved, the predicted velocity and turbulence fields are interpolated into a second grid and used to solve the concentration equation. To evaluate the model, various steady state numerical solutions are compared with small scale dispersion experiments which were conducted at the wind tunnel of Mitsubishi Heavy Industries, in Japan. Stably stratified flows and plume dispersion over three distinct idealized complex topographies (flat and hilly terrain) are studied. Vertical profiles of velocity and pollutant concentration are shown and discussed. Also, comparisons are made against the results obtained with the standard k-epsilon model.