Modeling of turbulent flow through radial diffuser

The work considers the modeling of turbulent flow in radial diffuser with axial feeding. Due to its claimed capability to predict flow including features such as separation, curvature and adverse pressure gradient, the RNG k-epsilon model of Orzag et al. (1993) is applied in the present analysis. The governing equations are numerically solved using the finite volume methodology. Experiments were conducted to assess the turbulence model. Numerical results of pressure distribution on the front disk surface for different flow conditions when compared to the experimental data indicated that the RNG k-epsilon model is adequate to predict this class of flow.

Solution of aerospace problems using structured and unstructured strategies

Products developed at industries, institutes and research centers are expected to have high level of quality and performance, having a minimum waste, which require efficient and robust tools to numerically simulate stringent project conditions with great reliability. In this context, Computational Fluid Dynamics (CFD) plays an important role and the present work shows two numerical algorithms that are used in the CFD community to solve the Euler and Navier-Stokes equations applied to typical aerospace and aeronautical problems. Particularly, unstructured discretization of the spatial domain has gained special attention by the international community due to its ease in discretizing complex spatial domains. This work has the main objective of illustrating some advantages and disadvantages of numerical algorithms using structured and unstructured spatial discretization of the flow governing equations. Numerical methods include a finite volume formulation and the Euler and Navier-Stokes equations are applied to solve a transonic nozzle problem, a low supersonic airfoil problem and a hypersonic inlet problem. In a structured context, these problems are solved using MacCormack’s implicit algorithm with Steger and Warming’s flux vector splitting technique, while, in an unstructured context, Jameson and Mavriplis’ explicit algorithm is used. Convergence acceleration is obtained using a spatially variable time stepping procedure.

Object-oriented Programming Applied to the Development of Structural Analysis and Optimization Software

This paper presents the development of a two-dimensional interactive software environment for structural analysis and optimization based on object-oriented programming using the C++ language. The main feature of the software is the effective integration of several computational tools into graphical user interfaces implemented in the Windows-98 and Windows-NT operating systems. The interfaces simplify data specification in the simulation and optimization of two-dimensional linear elastic problems. NURBS have been used in the software modules to represent geometric and graphical data. Extensions to the analysis of three-dimensional problems have been implemented and are also discussed in this paper.

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.

A Design Methodology for the Compensation of Positioning Deviation in Gantry Manipulators

This work presents recent results concerning a design methodology used to estimate the positioning deviation for a gantry (Cartesian) manipulator, related mainly to structural elastic deformation of components during operational conditions. The case-study manipulator is classified as gantry type and its basic dimensions are 1,53m x 0,97m x 1,38m. The dimensions used for the calculation of effective workspace due to end-effector path displacement are: 1m x 0,5m x 0,5m. The manipulator is composed by four basic modules defined as module X, module Y, module Z and terminal arm, where is connected the end-effector. Each module controlled axis performs a linear-parabolic positioning movement. The planning path algorithm has the maximum velocity and the total distance as input parameters for a given task. The acceleration and deceleration times are the same. Denavit-Hartemberg parameterization method is used in the manipulator kinematics model. The gantry manipulator can be modeled as four rigid bodies with three degrees-of-freedom in translational movements, connected as an open kinematics chain. Dynamic analysis were performed considering inertial parameters specification such as component mass, inertia and center of gravity position of each module. These parameters are essential for a correct manipulator dynamic modelling, due to multiple possibilities of motion and manipulation of objects with different masses. The dynamic analysis consists of a mathematical modelling of the static and dynamic interactions among the modules. The computation of the structural deformations uses the finite element method (FEM).

Numerical study of wedge supported oblique shock wave-oblique detonation wave transitions

The results of a numerical study of premixed Hydrogen-air flows ignition by an oblique shock wave (OSW) stabilized by a wedge are presented, in situations when initial and boundary conditions are such that transition between the initial OSW and an oblique detonation wave (ODW) is observed. More precisely, the objectives of the paper are: (i) to identify the different possible structures of the transition region that exist between the initial OSW and the resulting ODW and (ii) to evidence the effect on the ODW of an abrupt decrease of the wedge angle in such a way that the final part of the wedge surface becomes parallel to the initial flow. For such a geometrical configuration and for the initial and boundary conditions considered, the overdriven detonation supported by the initial wedge angle is found to relax towards a Chapman-Jouguet detonation in the region where the wedge surface is parallel to the initial flow. Computations are performed using an adaptive, unstructured grid, finite volume computer code previously developed for the sake of the computations of high speed, compressible flows of reactive gas mixtures. Physico-chemical properties are functions of the local mixture composition, temperature and pressure, and they are computed using the CHEMKIN-II subroutines.

Voimalaitoksen kanavien kannakkeiden rakennesuunnittelu

Tässä diplomityössä esitetään voimalaitoksen kanavien kannakkeiden rakennesuunnittelussa tarvittavat laskentamenetelmät. Työssä rakenteiden suunnitteluun ja mitoitukseen käytetään pääasiassa Eurokoodi 3 teräsrakenteiden suunnittelustandardin mukaista rajatilamitoitusta. Lisäksi kehitetään mitoitustyökaluja tärkeimpien kanavakannakkeiden suunnitteluun. Toteutettujen mitoitustyökalujen toiminta verifioidaan lujuusopin elementtimenetelmällä tehtävin tarkistuslaskelmin. Laskentatyökalujen analyyttisen ratkaisun verifioitiin olevan varmalla puolella kaikissa tutkituissa ilmiöissä. Työssä verifioituja menetelmiä voidaan soveltaa myös muiden vastaavien rakenteiden mitoittamiseen. Työssä luotujen laskentatyökalujen sisältämät laskentamenetelmät mahdollistavat monenlaisten rakenteiden vaatimustenmukaisen suunnittelun.

Turveperävaunun väsymiskestävyyden parantaminen

Työssä tutkittiin turveperävaunun renkaan kiinnitysrakenteen väsymiskestävyyttä. Väsymiskestävyyden parantamiseksi renkaiden kiinnityskohdan muotoilu ja mitoitus suunniteltiin uudelleen. Suunnittelussa keskityttiin ensisijaisesti nimellisten ja rakenteellisten jännitysheilahdusten pienentämiseen. Sekä vanhan, että uuden rakenneratkaisun väsymiskestoikiä tarkasteltiin käyttämällä hot spot – menetelmää.

Value engineering of InnoTrackTM

Middle section module of InnoTrackTM moving walk was re-engineered according to value analysis process. Self-supporting steel structure for moving walk was created as a result of this process. Designed structure was verified and validated by prototype tests and finite element method calculations. Self-supporting steel structure replaces the original design of middle section module in InnoTrackTM. Designed structure provides higher satisfaction to customers’ needs and at the same time, it uses less resources. The redesigned middle section module provides higher value to the customer.

Weed control via intercropping with gliricidia: II. corn crop

A reduction in herbicide use is one of modern agriculture's main interests and several alternatives are being investigated with this objective, including intercropping. Gliricídia (Gliricidia sepium) mulch has no allelopathic effect on corn or beans but significantly decreased the population of some weed species. The objective of this study was to evaluate green ear and grain yield in corn cultivars as a response to weed control achieved via intercropping with gliricidia. A completely randomized block design with five replicates and split-plots was used. Cultivars AG 1051, AG 2060, BRS 2020, and PL 6880 (assigned to plots) were submitted to the following treatments: no hoeing, hoeing (performed at 20 and 40 days after sowing the corn), and corn intercropped with gliricidia. Gliricidia was grown in a transplanting system to ensure uniform germination and fast establishment in the field. Seeding was made in 200-cell trays with one seed per cell (35 mL volume). The plants emerged two to three days after sowing and were transplanted to a permanent site two to three days after emergence. Corn was sown on the same day gliricidia was transplanted. Sixteen weed species occurred at different frequencies, with uneven distribution in the experimental area. Cultivars AG 1051 and AG 2060 were the best with reference to most characteristics employed to evaluate green corn yield. Cultivar AG 1051 provided the highest grain yield. The highest green ear yield and grain yield values were obtained with hoeing. However, the fact that intercropped plots showed intermediate yield between the values obtained for hoed and non-hoed plots indicates that gliricidia was beneficial to corn, and exerted a certain level of weed control.

Vääntökuormitetun koteloprofiilipuomin väsymiskestävyys

Tässä työssä tutkittiin vääntökuormitettua hitsattua koteloprofiilipuomia. Vääntökuormitus aiheuttaa koteloprofiiliin vinouttavan voimasysteemin, joka aiheuttaa kotelopalkkiin sekä poikittaisia taivutusjännityksiä että pitkittäisiä jännityksiä. Vinoutumisen aiheuttamia lisära-situksia on tutkittu analyyttisesti BEF-analogian avulla sekä elementtimenetelmää apuna käyttäen. Lisäksi vääntökuormitus aiheuttaa puomiin estetyn väännön jännityksiä. Väsymiskestoikää tutkittiin laboratoriossa suoritettujen väsytyskokeiden avulla sekä las-kennallisesti. Kestoikälaskennassa käytettiin tehollisen lovijännityksen menetelmää hitsin juuren puolen väsymisen arvioinnissa sekä hot spot- jännityksen menetelmää hitsin rajavii-van väsymisen arvioinnissa. Teholliset lovijännitykset sekä hot spot jännitykset ovat määri-tetty elementtimenetelmän avulla. Laboratoriokokeiden ja elementtimenetelmä laskennan perusteella saatiin rakenteen kes-toiän kannalta kriittiset detaljit määritettyä. Kriittiset detaljit sijaitsevat puomin päädyn rakenteissa. Tutkimuksessa saatiin selville, että kriittisten detaljien rakenteellisilla ja valmis-tusteknisillä ratkaisuilla on merkittävä vaikutus lopullisen tuotteen väsymiskestoikään.

Vibrations in Rotating Machinery Arising from Minor Imperfections in Component Geometries

Vibrations in machines can cause noise, decrease the performance, or even damage the machine. Vibrations appear if there is a source of vibration that excites the system. In the worst case scenario, the excitation frequency coincides with the natural frequency of the machine causing resonance. Rotating machines are a machine type, where the excitation arises from the machine itself. The excitation originates from the mass imbalance in the rotating shaft, which always exists in machines that are manufactured using conventional methods. The excitation has a frequency that is dependent on the rotational speed of the machine. The rotating machines in industrial use are usually designed to rotate at a constant rotational speed, the case where the resonances can be easily avoided. However, the machines that have a varying operational speed are more problematic due to a wider range of frequencies that have to be avoided. Vibrations, which frequencies equal to rotational speed frequency of the machine are widely studied and considered in the typical machine design process. This study concentrates on vibrations, which arise from the excitations having frequencies that are multiples of the rotational speed frequency. These vibrations take place when there are two or more excitation components in a revolution of a rotating shaft. The dissertation introduces four studies where three kinds of machines are experiencing vibrations caused by different excitations. The first studied case is a directly driven permanent magnet generator used in a wind power plant. The electromagnetic properties of the generator cause harmonic excitations in the system. The dynamic responses of the generator are studied using the multibody dynamics formulation. In another study, the finite element method is used to study the vibrations of a magnetic gear due to excitations, which frequencies equal to the rotational speed frequency. The objective is to study the effects of manufacturing and assembling inaccuracies. Particularly, the eccentricity of the rotating part with respect to non-rotating part is studied since the eccentric operation causes a force component in the direction of the shortest air gap. The third machine type is a tube roll of a paper machine, which is studied while the tube roll is supported using two different structures. These cases are studied using different formulations. In the first case, the tube roll is supported by spherical roller bearings, which have some wavinesses on the rolling surfaces. Wavinesses cause excitations to the tube roll, which starts to resonate at the frequency that is a half of the first natural frequency. The frequency is in the range where the machine normally operates. The tube roll is modeled using the finite element method and the bearings are modeled as nonlinear forces between the tube roll and the pedestals. In the second case studied, the tube roll is supported by freely rotating discs, which wavinesses are also measured. The above described phenomenon is captured as well in this case, but the simulation methodology is based on the flexible multibody dynamics formulation. The simulation models that are used in both of the last two cases studied are verified by measuring the actual devices and comparing the simulated and measured results. The results show good agreement.

Trukkitraktorin rungon suunnittelu ja laskenta

Tässä työssä suunniteltiin lappeenrantalaisen Astex Oy:n tilauksesta liikkuvan työkoneen runkorakenne. Työ tehtiin Lappeenrannan teknillisen yliopiston teräsrakenteiden laboratoriossa. Suunniteltava rakenne kuului linkkuohjattuun, noin 5000 kg painoiseen trukkitraktoriin. Lähtökohtana rakenteen suunnittelulle olivat tilaajan asettamat rakenteen geometriaan ja suorituskykyyn liittyvät reunaehdot ja rajoitteet. Uuden rakenteen suunnittelussa hyödynnettiin myös tilaajatahon kehittelemää vastaavan tyyppistä prototyyppirakennetta. Rakenteen suunnitteluprosessi koostui neljästä eri vaiheesta. Suunniteltavan rakenteen lähtökohtana olleelle prototyyppirakenteelle suoritettiin koneen käytönaikaisia venymäliuskamittauksia rakenteen kuormitushistorian selvittämiseksi. Mittauksista saatujen tulosten perusteella määritettiin kuormitukset uudelle suunniteltavalle rakenteelle. Määritettyjä kuormituksia hyödyntäen ideoitiin, suunniteltiin ja mallinnettiin uusi tilaajan vaatimuksia vastaava rakenne. Uudelle rakenteelle suoritettiin lujuusanalyysit FE-analyysiä hyödyntäen. Uuden rakenteen suunnittelussa kiinnitettiin huomiota rakenteen hyvään valmistettavuuteen ja suunniteltiin rakenneratkaisut tilaajataholle mahdollisimman optimaalisiksi. Suunnittelu- ja mallinnustyö tehtiin Solidworks 2014 ohjelmistolla. Rakenteen lujuustekniset tarkastelut sisälsivät rakennedetaljien analyyttistä mitoitusta ja laskentaa. FE-laskennalla selvitettiin rakenteen ääri- ja väsymiskestävyys. Laskenta sisälsi koko rakenteen globaaleja tarkasteluja, sekä eri kriittisten rakennedetaljien yksityiskohtaisempia analyysejä. FE-laskennan pääpaino oli rakenteen väsymisanalyyseissä, jotka toteutettiin Hot-Spot- ja särönkasvumenetelmillä. Rakenteen FE-analyysien suorittamisessa käytettiin Femap, NxNastran ja Abaqus-ohjelmistoja.

A priori estimation of accuracy and of the number of wells to be employed in limiting dilution assays

The use of limiting dilution assay (LDA) for assessing the frequency of responders in a cell population is a method extensively used by immunologists. A series of studies addressing the statistical method of choice in an LDA have been published. However, none of these studies has addressed the point of how many wells should be employed in a given assay. The objective of this study was to demonstrate how a researcher can predict the number of wells that should be employed in order to obtain results with a given accuracy, and, therefore, to help in choosing a better experimental design to fulfill one's expectations. We present the rationale underlying the expected relative error computation based on simple binomial distributions. A series of simulated in machina experiments were performed to test the validity of the a priori computation of expected errors, thus confirming the predictions. The step-by-step procedure of the relative error estimation is given. We also discuss the constraints under which an LDA must be performed.