Simulation of Steady-State Nonlinear Heat Transfer Problems Through the Minimization of Quadratic Functionals

In this work it is presented a systematic procedure for constructing the solution of a large class of nonlinear conduction heat transfer problems through the minimization of quadratic functionals like the ones usually employed for linear descriptions. The proposed procedure gives rise to an efficient and easy way for carrying out numerical simulations of nonlinear heat transfer problems by means of finite elements. To illustrate the procedure a particular problem is simulated by means of a finite element approximation.

Simulation of Active Control Using Fuzzy Logic Applied to a Pulse Combustor

This work analyzes an active fuzzy logic control system in a Rijke type pulse combustor. During the system development, a study of the existing types of control for pulse combustion was carried out and a simulation model was implemented to be used with the package Matlab and Simulink. Blocks which were not available in the simulator library were developed. A fuzzy controller was developed and its membership functions and inference rules were established. The obtained simulation showed that fuzzy logic is viable in the control of combustion instabilities. The obtained results indicated that the control system responded to pulses in an efficient and desirable way. It was verified that the system needed approximately 0.2 s to increase the tube internal pressure from 30 to 90 mbar, with an assumed total delay of 2 ms. The effects of delay variation were studied. Convergence was always obtained and general performance was not affected by the delay. The controller sends a pressure signal in phase with the Rijke tube internal pressure signal, through the speakers, when an increase the oscillations pressure amplitude is desired. On the other hand, when a decrease of the tube internal pressure amplitude is desired, the controller sends a signal 180º out of phase.

Wind tunnel simulation of atmospheric boundary layer flows

The present work shows how thick boundary layers can be produced in a short wind tunnel with a view to simulate atmospheric flows. Several types of thickening devices are analysed. The experimental assessment of the devices was conducted by considering integral properties of the flow and the spectra: skin-friction, mean velocity profiles in inner and outer co-ordinates and longitudinal turbulence. Designs based on screens, elliptic wedge generators, and cylindrical rod generators are analysed. The paper describes in detail the experimental arrangement, including the features of the wind tunnel and of the instrumentation. The results are compared with experimental data published by other authors and with naturally developed flows.

Seismic motion simulation based on Cassino Parallel Manipulator

In this paper we present a study of feasibility by using Cassino Parallel Manipulator (CaPaMan) as an earthquake simulator. We propose a suitable formulation to simulate the frequency, amplitude and acceleration magnitude of seismic motion by means of the movable platform motion by giving a suitable input motion. In this paper we have reported numerical simulations that simulate the three principal earthquake types for a seismic motion: one at the epicenter (having a vertical motion), another far from the epicenter (with the motion on a horizontal plane), and a combined general motion (with a vertical and horizontal motion).

Neutron tomography using projection data obtained by Monte Carlo simulation for nondestructive evaluation

This work present the application of a computer package for generating of projection data for neutron computerized tomography, and in second part, discusses an application of neutron tomography, using the projection data obtained by Monte Carlo technique, for the detection and localization of light materials such as those containing hydrogen, concealed by heavy materials such as iron and lead. For tomographic reconstructions of the samples simulated use was made of only six equal projection angles distributed between 0º and 180º, with reconstruction making use of an algorithm (ARIEM), based on the principle of maximum entropy. With the neutron tomography it was possible to detect and locate polyethylene and water hidden by lead and iron (with 1cm-thick). Thus, it is demonstrated that thermal neutrons tomography is a viable test method which can provide important interior information about test components, so, extremely useful in routine industrial applications.

Testing and simulation of fractionary electromechanical rotative drives

This paper concerns the development of drives that use electromechanical rotative motor systems. It is proposed an experimental drive test structure integrated to simulation softwares. The objective of this work is to show that an affordable model validation procedure can be obtained by combining a precision data acquisition with well tuned state-of-the-art simulation packages. This is required for fitting, in the best way, a drive to its load or, inversely, to adapt loads to given drive characteristics.

Numerical Simulation of Stochastic Di erential Equations

Stochastic differential equation (SDE) is a differential equation in which some of the terms and its solution are stochastic processes. SDEs play a central role in modeling physical systems like finance, Biology, Engineering, to mention some. In modeling process, the computation of the trajectories (sample paths) of solutions to SDEs is very important. However, the exact solution to a SDE is generally difficult to obtain due to non-differentiability character of realizations of the Brownian motion. There exist approximation methods of solutions of SDE. The solutions will be continuous stochastic processes that represent diffusive dynamics, a common modeling assumption for financial, Biology, physical, environmental systems. This Masters' thesis is an introduction and survey of numerical solution methods for stochastic differential equations. Standard numerical methods, local linearization methods and filtering methods are well described. We compute the root mean square errors for each method from which we propose a better numerical scheme. Stochastic differential equations can be formulated from a given ordinary differential equations. In this thesis, we describe two kind of formulations: parametric and non-parametric techniques. The formulation is based on epidemiological SEIR model. This methods have a tendency of increasing parameters in the constructed SDEs, hence, it requires more data. We compare the two techniques numerically.

Development of a process for the direct synthesis of hydrogen peroxide in a novel microstructured reactor

Microreactors have proven to be versatile tools for process intensification. Over recent decades, they have increasingly been used for product and process development in chemical industries. Enhanced heat and mass transfer in the reactors due to the extremely high surfacearea- to-volume ratio and interfacial area allow chemical processes to be operated at extreme conditions. Safety is improved by the small holdup volume of the reactors and effective control of pressure and temperature. Hydrogen peroxide is a powerful green oxidant that is used in a wide range of industries. Reduction and auto-oxidation of anthraquinones is currently the main process for hydrogen peroxide production. Direct synthesis is a green alternative and has potential for on-site production. However, there are two limitations: safety concerns because of the explosive gas mixture produced and low selectivity of the process. The aim of this thesis was to develop a process for direct synthesis of hydrogen peroxide utilizing microreactor technology. Experimental and numerical approaches were applied for development of the microreactor. Development of a novel microreactor was commenced by studying the hydrodynamics and mass transfer in prototype microreactor plates. The prototypes were designed and fabricated with the assistance of CFD modeling to optimize the shape and size of the microstructure. Empirical correlations for the mass transfer coefficient were derived. The pressure drop in micro T-mixers was investigated experimentally and numerically. Correlations describing the friction factor for different flow regimes were developed and predicted values were in good agreement with experimental results. Experimental studies were conducted to develop a highly active and selective catalyst with a proper form for the microreactor. Pd catalysts supported on activated carbon cloths were prepared by different treatments during the catalyst preparation. A variety of characterization methods were used for catalyst investigation. The surface chemistry of the support and the oxidation state of the metallic phase in the catalyst play important roles in catalyst activity and selectivity for the direct synthesis. The direct synthesis of hydrogen peroxide was investigated in a bench-scale continuous process using the novel microreactor developed. The microreactor was fabricated based on the hydrodynamic and mass transfer studies and provided a high interfacial area and high mass transfer coefficient. The catalysts were prepared under optimum treatment conditions. The direct synthesis was conducted at various conditions. The thesis represents a step towards a commercially viable direct synthesis. The focus is on the two main challenges: mitigating the safety problem by utilization of microprocess technology and improving the selectivity by catalyst development.

Häiriöanalyysejä tehonkorotetun kiehutusvesivoimalaitoksen käytettävyyden varmistamiseksi

TVO suunnittelee reaktoritehon 10 %:n korotusta Olkiluoto 1 ja 2 -voimalaitoksille. Reaktoriteho nostetaan 2500 MW:sta 2750 MW:iin polttoaineen rikastusastetta nostamalla ja pääkiertovirtausta kasvattamalla. Samalla syöttövesivirtaus reaktoriin ja tuorehöyryvirtaus turpiineille kasvaa. Lauhteenpuhdistusjärjestelmän kapasiteettia ei voida kuitenkaan kasvattaa, joten massavirran lisäys toteutetaan ottamalla käyttöön korkeapainesivulauhteen eteenpäinpumppaus. Lauhteen esilämmityslinjojen, lauhteenpuhdistuksen ja syöttövesipumppujen massavirta säilyy siten nykyisellään. Muita merkittäviä tehonkorotukseen liittyviä laitosmuutoksia ovat pääkiertopumppujen uusinta ja korkeapaineturpiinin muutokset. Tehonkorotetun prosessin käytettävyyden varmistamiseksi tehdään häiriöanalyysejä Apros-prosessisimulointiohjelmistoa käyttäen. OL1 ja OL2 -laitoksista on olemassa validoitu 2500 MW:n laitosmalli, josta muokatulla 2750 MW:n laitosmallilla simuloinnit tehdään. Häiriöanalyysien avulla selvitetään säätöjärjestelmien kyky pitää prosessin tila hallinnassa ilman suojausautomaation laukeamista. Simuloituihin tapauksiin kuuluu pumppujen ja venttiilien vikaantumistapauksia sekä turpiini- ja reaktoripuolen pikasulku- ja osittaispikasulkutapauksia. Myös meriveden lämpötilan vaikutusta häiriötilanteisiin tarkastellaan. Analyysien perusteella voimalaitosten ohjaus- ja suojausautomaatio toimivat hyvin myös korotetulla teholla. Tehonkorotuksen jälkeiset suuremmat massavirrat aiheuttavat kuitenkin voimakkaampia reaktoripaineen ja -tehon vaihteluita varsinkin venttiilien sulkeutumistapauksissa. Simuloinnit osoittivat, että tehonkorotus 2500 MW:sta 2750 MW:iin on mahdollinen, mutta aiheuttaa pieniä muutoksia laitoksen suojausjärjestelmien laukaisurajoihin.

Integration of Torrefaction with Steam Power Plant

Torrefaction is one of the pretreatment technologies to enhance the fuel characteristics of biomass. The efficient and continuous operation of a torrefaction reactor, in the commercial scale, demands a secure biomass supply, in addition to adequate source of heat. Biorefinery plants or biomass-fuelled steam power plants have the potential to integrate with the torrefaction reactor to exchange heat and mass, using available infrastructure and energy sources. The technical feasibility of this integration is examined in this study. A new model for the torrefaction process is introduced and verified by the available experimental data. The torrefaction model is then integrated in different steam power plants to simulate possible mass and energy exchange between the reactor and the plants. The performance of the integrated plant is investigated for different configurations and the results are compared.

OL1/OL2-suojarakennuksen jäähdytteenmenetysonnettomuusanalyysit APROS 5.09 -ohjelmistolla

Tässä työssä on tutkittu OL1/OL2-suojarakennuksen käyttäytymistä jäähdytteenmenetysonnettomuuden eli LOCA:n aikana. Onnettomuuden simulointiin on kehitetty suojarakennusmalli suomalaiseen APROS 5.09 - ohjelmistoon sisältyvällä LP-koodilla (Lumped Parameter Code). Työssä on keskitytty suojarakennuksen kannalta oleellisimpien suureiden: kaasutilavuuksien paineen sekä lämpötilan ja lauhdutusaltaan lämpötilan ja pinnankorkeuden ajalliseen käyttäytymiseen. Mallinnetut LOCA:t ovat päähöyrylinjan ja sammutetun reaktorin jäähdytysjärjestelmän putkikatkoksia. Simulointeja on tehty laitoksen täyden tehon ja kuumavalmiuden lähtötiloissa ja tarkasteltavien suureiden käyttäytymistä on tutkittu erikseen valituilla konservatiivisilla oletuksilla. Laskentatapaukset rajoittuvat lyhyeen aikaväliin 27.78 tuntia kuvitellusta putkirikosta eteenpäin. Tuloksia on verrattu OL1/OL2- laitostoimittajan, Westinghouse Electric Sweden AB:n, tekemiin lisensiointianalyyseihin. Työssä on myös kuvattu APROS-laskennan epävarmuustekijöitä. Suojarakennuksen on todettu käyttäytyvän fysikaalisesti yhtenevästi lisensiointianalyyseissä ja APROS:lla tehdyissä vertailulaskuissa.

Computational studies for the design of process equipment with complex geometries

This study combines several projects related to the flows in vessels with complex shapes representing different chemical apparata. Three major cases were studied. The first one is a two-phase plate reactor with a complex structure of intersecting micro channels engraved on one plate which is covered by another plain plate. The second case is a tubular microreactor, consisting of two subcases. The first subcase is a multi-channel two-component commercial micromixer (slit interdigital) used to mix two liquid reagents before they enter the reactor. The second subcase is a micro-tube, where the distribution of the heat generated by the reaction was studied. The third case is a conventionally packed column. However, flow, reactions or mass transfer were not modeled. Instead, the research focused on how to describe mathematically the realistic geometry of the column packing, which is rather random and can not be created using conventional computeraided design or engineering (CAD/CAE) methods. Several modeling approaches were used to describe the performance of the processes in the considered vessels. Computational fluid dynamics (CFD) was used to describe the details of the flow in the plate microreactor and micromixer. A space-averaged mass transfer model based on Fick’s law was used to describe the exchange of the species through the gas-liquid interface in the microreactor. This model utilized data, namely the values of the interfacial area, obtained by the corresponding CFD model. A common heat transfer model was used to find the heat distribution in the micro-tube. To generate the column packing, an additional multibody dynamic model was implemented. Auxiliary simulation was carried out to determine the position and orientation of every packing element in the column. This data was then exported into a CAD system to generate desirable geometry, which could further be used for CFD simulations. The results demonstrated that the CFD model of the microreactor could predict the flow pattern well enough and agreed with experiments. The mass transfer model allowed to estimate the mass transfer coefficient. Modeling for the second case showed that the flow in the micromixer and the heat transfer in the tube could be excluded from the larger model which describes the chemical kinetics in the reactor. Results of the third case demonstrated that the auxiliary simulation could successfully generate complex random packing not only for the column but also for other similar cases.

Simulation of new p-type pixel strip detector with enhanced multiplication effect

This Master’s Thesis is dedicated to the simulation of new p-type pixel strip detector with enhanced multiplication effect. It is done for high-energy physics experiments upgrade such as Super Large Hadron Collider especially for Compact Muon Solenoid particle track silicon detectors. These detectors are used in very harsh radiation environment and should have good radiation hardness. The device engineering technology for developing more radiation hard particle detectors is used for minimizing the radiation degradation. New detector structure with enhanced multiplication effect is proposed in this work. There are studies of electric field and electric charge distribution of conventional and new p-type detector under reverse voltage bias and irradiation. Finally, the dependence of the anode current from the applied cathode reverse voltage bias under irradiation is obtained in this Thesis. For simulation Silvaco Technology Computer Aided Design software was used. Athena was used for creation of doping profiles and device structures and Atlas was used for getting electrical characteristics of the studied devices. The program codes for this software are represented in Appendixes.

Transient current technique simulation of semiconductor detectors

Nowadays advanced simulation technologies of semiconductor devices occupies an important place in microelectronics production process. Simulation helps to understand devices internal processes physics, detect new effects and find directions for optimization. Computer calculation reduces manufacturing costs and time. Modern simulation suits such as Silcaco TCAD allow simulating not only individual semiconductor structures, but also these structures in the circuit. For that purpose TCAD include MixedMode tool. That tool can simulate circuits using compact circuit models including semiconductor structures with their physical models. In this work, MixedMode is used for simulating transient current technique setup, which include detector and supporting electrical circuit. This technique was developed by RD39 collaboration project for investigation radiation detectors radiation hard properties.