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.

Dynamiikan simulointi poralaitteiden tuotekehityksessä

Diplomityö käsittelee kallioporalaitteen eliniän kuormitusten arviointiin tarkoitetun simulointimallin rakentamista. Työssä luotiin modulaarinen malli, jolla voidaan simuloida dynaamisesti esteen yliajo. Esteen yliajo on perinteinen tehdastesti uusien laitteiden prototyypeille. Sillä pyritään saamaan selville laitteen maksimikuormitustapaukset. Pintaporalaitteen simulointimalli tehtiin ADAMS-ohjelmistolla laitteen suunnittelussa tehtyjä CAD-malleja hyödyntäen. ADAMS-ohjelmistolla mallinnettiin laitteen yksinkertaistettu mekaniikka siten, että laite oli jaettu viiteen eri moduuliin. Runko ja telasto olivat omana osionaan. Puomi ja syöttölaite oli oma kokoonpanonsa ja ohjaamo sekä katteet olivat kaksi erillistä modulia. Viides moduuli oli hydrauliikka, joka simuloitiin rinnakkaissimulointina EASY5-ohjelmistolla. Simuloituja esteen yliajon tuloksia verrattiin prototyyppilaitteen tehdastestien mitattuihin kuormituksiin. Vertailu tehtiin mallinnettujen neljän sylinterin paineita tarkastelemalla. Tuloksia tarkastellessa havaittiin, että simulaatiomalli antaa varsinkin staattisissa tarkasteluissa oikean suuntaisia tuloksia. Dynaamisissa tilanteissa koneen maksimikuormituksilla mallinnustarkkuus ei tämän työn laajuudessa ole riittävä varsinaiseen eliniän analysointiin. Sen sijaan mallinnustekniikka todettiin periaatteessa toimivaksi ja jatkokehityskelpoiseksi.

Introduction of the T-10 static line parachuting capability to the NH90 helicopter

The purpose of this study was to investigate the suitability of the Finnish Defence Forces’ NH90 helicopter for parachuting operations with the T-10 static line parachute system. The work was based on the Army Command’s need to compensate for the reduction in the outsourced flight hours for the military static line parachuting training. The aim of the research was to find out the procedures and limitations with which the NH90 IOC+ or FOC version helicopter could be used for static line parachutist training with the T- 10B/MC1-1C parachutes. The research area was highly complicated and non-linear. Thus analytical methods could not be applied with sufficient confidence, even with present-day computing power. Therefore an empirical research method was selected, concentrating on flight testing supported with literature study and some calculated estimations. During three flights and 4.5 flight hours in Utti, Finland on 17−20 September 2012, a total of 44 parachute drops were made. These consisted of 16 dummy drops and 28 paratrooper jumps. The test results showed that when equipped with the floor mounted PASI-1 anchor line, the deflector bar of the NHIndustries’ Parachuting Kit and Patria’s floor protection panels the Finnish NH90 variant could be safely used for T-10B/MC1-1C static line parachuting operations from the right cabin door at airspeed range of 50−80 KIAS (∼90–150 km/h). The ceiling mounted anchor lines of the NHI’s Parachuting Kit were not usable with the T-10 system. This was due to the static lines’ unsafe behaviour in slipstream when connected to the cabin ceiling level. In conclusion, the NH90 helicopter can be used to meet the Army Command’s requirement for an additional platform for T-10 static line parachutist training. Material dropping, the effect of additional equipment and jumping from the rear ramp should be further studied.

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.