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.

Numerical Simulation of the Cavitation in the Hydrodynamic Lubrication of Journal Bearings: a Parallel Algorithm

In this paper we present an algorithm for the numerical simulation of the cavitation in the hydrodynamic lubrication of journal bearings. Despite the fact that this physical process is usually modelled as a free boundary problem, we adopted the equivalent variational inequality formulation. We propose a two-level iterative algorithm, where the outer iteration is associated to the penalty method, used to transform the variational inequality into a variational equation, and the inner iteration is associated to the conjugate gradient method, used to solve the linear system generated by applying the finite element method to the variational equation. This inner part was implemented using the element by element strategy, which is easily parallelized. We analyse the behavior of two physical parameters and discuss some numerical results. Also, we analyse some results related to the performance of a parallel implementation of the algorithm.

Direct Simulation Monte Carlo Method Applied to Aerothermodynamics

An axisymmetric supersonic flow of rarefied gas past a finite cylinder was calculated applying the direct simulation Monte Carlo method. The drag force, the coefficients of pressure, of skin friction, and of heat transfer, the fields of density, of temperature, and of velocity were calculated as function of the Reynolds number for a fixed Mach number. The variation of the Reynolds number is related to the variation of the Knudsen number, which characterizes the gas rarefaction. The present results show that all quantities in the transition regime (Knudsen number is about the unity) are significantly different from those in the hydrodynamic regime, when the Knudsen number is small.

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.

Seulonnan erotustehokkuuden määrittäminen - tarkastelussa Reichertin kartioseula

Seulonta on vanha ja teollisuudessa yleisesti käytetty erotusmenetelmä. Huolimat-ta yksinkertaisesta prosessista seulontaa ilmiönä on tutkittu vähän. Prosessin yk-sinkertaisuus on voinut osaltaan vaikuttaa tutkijoiden mielenkiintoon aihetta koh-taan. Tutkimuksen vähyys näkyy etenkin painovoimaiseen seulontaan liittyvän kirjallisuuden vähyytenä. Kirjallisuusosassa käsitellään seulonnan taustaa ja märkä- sekä kuivaseulontaa. Lisäksi määritetään seulontaan vaikuttavat tekijät ja kar-toitetaan seulontaan liittyvää matematiikkaa. Työn kokeellisessa osassa tutkitaan Reichertin kartioseulan erotustehokkuutta 10, 12,5 ja 15 asteen kallistuskulmissa. Tutkimukset tehdään ideaalisille pallomaisille lasihelmille. Mittaustulosten avulla tutkitaan myös partikkelien liikettä seulapin-nalla ja syötetyn massan vaikutusta erotustehokkuuteen. Kokeellisia mittaustulok-sia käytetään MATLAB-mallinnuksessa, jonka avulla optimoidaan kartioseulan kallistuskulmaa ja massapanosta suhteessa erotustehokkuuteen. Tutkimuksen kannalta olennaisia kysymyksiä ovat: ”Onko painovoimaan perustuva erottelu seuloilla riittävä takaamaan hyvän erotustehokkuuden?” ja ”Mikä on seulan ero-tustehokkuuden kannalta optimaalinen kallistuskulma ja syötettävä massa?” Reichertin kartioseulan etuja ovat sen energiatehokkuus ja yksinkertainen helposti säädettävä prosessi. Koemittauksista havaitaan, että erotustehokkuus vastaa opti-moinnin tuloksia. Seulan liian pieni kallistuskulma ja suuri massapanos pienentävät seulan erotustehokkuutta. Lisätutkimuksia tarvitaan eri kallistuskulmilla ja materiaaleilla, jotta tuloksia voidaan verrata todellisiin seulontaprosesseihin.

Rainfall effect on dissipation and movement of diuron and simazine in a vineyard soil

From 2003 to 2007, a field study was performed in a vineyard in Chile to investigate diuron and simazine soil behavior and the effect of additional rainfall. Both herbicides were applied once a year at a rate of 2.0 kg ha-1 a.i. Herbicide concentrations in soil were measured at 0, 10, 20, 40, 90 and 340 days after application, under two pluviometric conditions, natural rainfall and natural rainfall plus irrigation with 180 mm of simulated rainfall during the first 90 days after application. Soil partition coefficient (Kd) varied in the soil profile (0 to 90 cm deep) from 6.75 to 2.04 mL g-1 and from 1.4 to 0.66 mL g-1 and the maximum soil adsorption capacity was approximately 18.3 mg g-1 and 8.3 mg g-1 for diuron and simazine, respectively. Diuron and simazine reached up to 90 and 120 cm of soil depth, with an average of 8.3% and 62.4% of herbicide moved below 15 cm in the soil, respectively. Simazine soil half-life (DT50) was 38.1 days and 7.5 days, whereas the half life for diuron varied from 68.0 and 24.6 for natural rainfall and irrigated, respectively. The average of residual simazine remaining in the whole soil profile after 90 DAA was 25.4% and 39.9% for diuron, with no effect of additional rainfall amount. At 340 DAA the amount of simazine in the whole soil profile corresponded to 13.2% of the initial amount applied, being diuron more persistent with 21.5% of the initial herbicide applied. The high movement in soil of both herbicides could be due to a non-equilibrium sorption process explained by preferential flow, low Kd and high desorption.

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.