Mekanismien dynamiikan simuloinnissa sovellettavia numeerisia- ja mallinnusmenetelmiä

Koneet voidaan usein jakaa osajärjestelmiin, joita ovat ohjaus- ja säätöjärjestelmät, voimaa tuottavat toimilaitteet ja voiman välittävät mekanismit. Eri osajärjestelmiä on simuloitu tietokoneavusteisesti jo usean vuosikymmenen ajan. Osajärjestelmien yhdistäminen on kuitenkin uudempi ilmiö. Usein esimerkiksi mekanismien mallinnuksessa toimilaitteen tuottama voimaon kuvattu vakiona, tai ajan funktiona muuttuvana voimana. Vastaavasti toimilaitteiden analysoinnissa mekanismin toimilaitteeseen välittämä kuormitus on kuvattu vakiovoimana, tai ajan funktiona työkiertoa kuvaavana kuormituksena. Kun osajärjestelmät on erotettu toisistaan, on niiden välistenvuorovaikutuksien tarkastelu erittäin epätarkkaa. Samoin osajärjestelmän vaikutuksen huomioiminen koko järjestelmän käyttäytymissä on hankalaa. Mekanismien dynamiikan mallinnukseen on kehitetty erityisesti tietokoneille soveltuvia numeerisia mallinnusmenetelmiä. Useimmat menetelmistä perustuvat Lagrangen menetelmään, joka mahdollistaa vapaasti valittaviin koordinaattimuuttujiin perustuvan mallinnuksen. Numeerista ratkaisun mahdollistamiseksi menetelmän avulla muodostettua differentiaali-algebraaliyhtälöryhmää joudutaan muokkaamaan esim. derivoimalla rajoiteyhtälöitä kahteen kertaan. Menetelmän alkuperäisessä numeerisissa ratkaisuissa kaikki mekanismia kuvaavat yleistetyt koordinaatit integroidaan jokaisella aika-askeleella. Tästä perusmenetelmästä johdetuissa menetelmissä riippumattomat yleistetyt koordinaatit joko integroidaan ja riippuvat koordinaatit ratkaistaan rajoiteyhtälöiden perusteella tai yhtälöryhmän kokoa pienennetään esim. käyttämällä nopeus- ja kiihtyvyysanalyyseissä eri kiertymäkoordinaatteja kuin asema-analyysissä. Useimmat integrointimenetelmät on alun perin tarkoitettu differentiaaliyhtälöiden (ODE) ratkaisuunjolloin yhtälöryhmään liitetyt niveliä kuvaavat algebraaliset rajoiteyhtälöt saattavat aiheuttaa ongelmia. Nivelrajoitteiden virheiden korjaus, stabilointi, on erittäin tärkeää mekanismien dynamiikan simuloinnin onnistumisen ja tulosten oikeellisuuden kannalta. Mallinnusmenetelmien johtamisessa käytetyn virtuaalisen työn periaatteen oletuksena nimittäin on, etteivät rajoitevoimat tee työtä, eli rajoitteiden vastaista siirtymää ei tapahdu. Varsinkaan monimutkaisten järjestelmien pidemmissä analyyseissä nivelrajoitteet eivät toteudu tarkasti. Tällöin järjestelmän energiatasapainoei toteudu ja järjestelmään muodostuu virtuaalista energiaa, joka rikkoo virtuaalisen työn periaatetta, Tästä syystä tulokset eivät enää pidäpaikkaansa. Tässä raportissa tarkastellaan erityyppisiä mallinnus- ja ratkaisumenetelmiä, ja vertaillaan niiden toimivuutta yksinkertaisten mekanismien numeerisessa ratkaisussa. Menetelmien toimivuutta tarkastellaan ratkaisun tehokkuuden, nivelrajoitteiden toteutumisen ja energiatasapainon säilymisen kannalta.

Dynamic Model of a Small Once-Through Boiler

In this study, a model for the unsteady dynamic behaviour of a once-through counter flow boiler that uses an organic working fluid is presented. The boiler is a compact waste-heat boiler without a furnace and it has a preheater, a vaporiser and a superheater. The relative lengths of the boiler parts vary with the operating conditions since they are all parts of a single tube. The present research is a part of a study on the unsteady dynamics of an organic Rankine cycle power plant and it will be a part of a dynamic process model. The boiler model is presented using a selected example case that uses toluene as the process fluid and flue gas from natural gas combustion as the heat source. The dynamic behaviour of the boiler means transition from the steady initial state towards another steady state that corresponds to the changed process conditions. The solution method chosen was to find such a pressure of the process fluid that the mass of the process fluid in the boiler equals the mass calculated using the mass flows into and out of the boiler during a time step, using the finite difference method. A special method of fast calculation of the thermal properties has been used, because most of the calculation time is spent in calculating the fluid properties. The boiler was divided into elements. The values of the thermodynamic properties and mass flows were calculated in the nodes that connect the elements. Dynamic behaviour was limited to the process fluid and tube wall, and the heat source was regarded as to be steady. The elements that connect the preheater to thevaporiser and the vaporiser to the superheater were treated in a special way that takes into account a flexible change from one part to the other. The model consists of the calculation of the steady state initial distribution of the variables in the nodes, and the calculation of these nodal values in a dynamic state. The initial state of the boiler was received from a steady process model that isnot a part of the boiler model. The known boundary values that may vary during the dynamic calculation were the inlet temperature and mass flow rates of both the heat source and the process fluid. A brief examination of the oscillation around a steady state, the so-called Ledinegg instability, was done. This examination showed that the pressure drop in the boiler is a third degree polynomial of the mass flow rate, and the stability criterion is a second degree polynomial of the enthalpy change in the preheater. The numerical examination showed that oscillations did not exist in the example case. The dynamic boiler model was analysed for linear and step changes of the entering fluid temperatures and flow rates.The problem for verifying the correctness of the achieved results was that there was no possibility o compare them with measurements. This is why the only way was to determine whether the obtained results were intuitively reasonable and the results changed logically when the boundary conditions were changed. The numerical stability was checked in a test run in which there was no change in input values. The differences compared with the initial values were so small that the effects of numerical oscillations were negligible. The heat source side tests showed that the model gives results that are logical in the directions of the changes, and the order of magnitude of the timescale of changes is also as expected. The results of the tests on the process fluid side showed that the model gives reasonable results both on the temperature changes that cause small alterations in the process state and on mass flow rate changes causing very great alterations. The test runs showed that the dynamic model has no problems in calculating cases in which temperature of the entering heat source suddenly goes below that of the tube wall or the process fluid.

Välimiesmenettely, yrityksen vaihtoehtoinen riidanratkaisumenetelmä

Tämän tutkielman tavoitteena on tutkia yrityksen näkökulmasta niin kansallista kuin kansainvälistä ad hoc välimiesmenettelyä, välimiesmenettelylain ja YK:n kansainvälisen kauppaoikeuden toimikunnan UNCITRAL:in sääntöjen pohjalta. Tutkielmassa käsitellään välimiesmenettelyn käynnistymisen edellytyksiä ja itse menettelyä välimiesoikeudessa sekä ulkomaisen välitystuomion tunnustamista ja täytäntöönpanoa Suomessa. Tutkimuksessa käytettävä tutkimusmetodologia on käsiteanalyyttinen. Välimiesmenettely on lisännyt suosiotaan vaihtoehtoisena riidanratkaisumenetelmänä ja sitä käytetään varsinkin kaupan, teollisuuden ja rakennustoiminnan piirissä. Välimiesmenettelyn merkitys on suuri kansainvälisen kaupan alalla ja se on vakiinnuttanut asemansa yritysten välisten riitojen ratkaisukeinona sen tuomien etuuksien, kuten joustavuuden, nopeuden ja kansainvälisyyden johdosta.

Soodakattiloiden ilmakanavien hyödyntäminen jäykistävänä rakenteena

Työssä tutkittiin soodakattiloiden ilmakanavien hyödyntämistä jäykistävänä rakenteena. Työssä käsiteltiin yksittäisiä jäykistämättömiä ja jäykistettyjä levykenttiä ja niiden lommahduskestävyyttä Eurokoodi standardin mukaisesti ja elementtimenetelmän avulla. Lisäksi käsiteltiin lommahduksen teoriaa ja levykenttien käyttäytymistä yleisellä tasolla erilaisilla kuormituksilla ja reunaehdoilla. Työn tavoitteena oli selvittää kuinka lommahdus tutkitaan Eurokoodin mukaisesti ja elementtimenetelmää hyödyntäen, kun levykentän kuormituksena on poikittainen kuormitus tason suuntaisen kuormituksen lisäksi. Työssä tutkittiin kahden eri elementtimenetelmään pohjautuvan ratkaisuvaihtoehdon käyttöä lommahduslaskennassa. Työssä kehitettiin Eurokoodin sovellettu yhteisvaikutuskaavan käyttö lineaarisen ominaisarvotehtävän ratkaisun lisänä, jossa otetaan huomioon painekuorman vaikutus levykentän lommahduksessa. Kehitettyä menetelmää sovellettiin ilmakanavan esimerkkirakenteen mitoituksessa.

Kunnossapidon suorituskykymittariston kehittäminen

Fazerin kolmen suurimman suomalaisleipomon kesken on sovittu yhteisestä kunnossapidon suorituskykymittareista, mutta mittariston kehitykseen ei ollut käytetty mitään tiettyä menetelmää, joka varmistaisi, että käyttöönotetut mittarit olisivat linjassa yrityksen tai sen ydintoimintojen tavoitteiden kanssa. Tämän työn on tarkoitus tutkia, kuinka hyvä leipomoiden käyttämä nykyinen kunnossapidon suorituskykymittaristo on, ja voiko sitä parantaa entisestään. Ratkaisuksi arviointiongelmaan luotiin suorituskykymittariston kehitysmenetelmä ja sitä tukeva työkalu, jolla mittareita voidaan arvottaa ja vertailla keskenään. Kehitysmenetelmän luomisessa hyödynnettiin kirjallisuudesta löytyneitä teorioita suorituskyvyn mittaamisen periaatteista ja valmista menetelmää mittariston kehittämistä varten. Työn ohessa syntyi myös mittariston kehitysmenetelmää tukeva työkalu, joka helpottaa mittareiden arvottamista ja keskinäistä vertailua. Luodulla suorituskykymittariston kehitysmenetelmällä ja sitä tukevalla työkalulla saatiin arvioitua nykyisin käytössä olevan mittariston hyvyys suhteessa vaihtoehtoisiin, kilpaileviin mittareihin. Tulosten perusteella Fazerin Lahden leipomolla on päätetty pilotoida uusia kunnossapidon suorituskykymittareita.

DNA variation in a 13-Mb region including the F9 gene: inferring the genealogical history and causal role of a hemophilia B mutation (IVS 5+13 A -> G)

About 5.5% of all UK hemophilia B patients have the base substitution IVS 5+13 A-->G as the only change in their factor (F)IX gene (F9). This generates a novel donor splice site which fits the consensus better than the normal intron 5 donor splice. Use of the novel splice site should result in a missense mutation followed by the abnormal addition of four amino acids to the patients' FIX. In order to explain the prevalence of this mutation, its genealogical history is examined. Analysis of restriction fragment length polymorphism in the 21 reference UK individuals (from different families) with the above mutation showed identical haplotypes in 19 while two differed from the rest and from each other. In order to investigate the history of the mutation and to verify that it had occurred independently more than once, the sequence variation in 1.5-kb segments scattered over a 13-Mb region including F9 was examined in 18 patients and 15 controls. This variation was then analyzed with a recently developed Bayesian approach that reconstructs the genealogy of the gene investigated while providing evidence of independent mutations that contribute disconnected branches to the genealogical tree. The method also provides minimum estimates of the age of the mutation inherited by the members of coherent trees. This revealed that 17 or 18 mutant genes descend from a founder who probably lived 450 years ago, while one patient carries an independent mutation. The independent recurrence of the IVS5+13 A-->G mutation strongly supports the conclusion that it is the cause of these patients' mild hemophilia.

Non-steady state heat conduction in composite walls

The problem of heat conduction in one-dimensional piecewise homogeneous composite materials is examined by providing an explicit solution of the one-dimensional heat equation in each domain. The location of the interfaces is known, but neither temperature nor heat flux are prescribed there. Instead, the physical assumptions of their continuity at the interfaces are the only conditions imposed. The problem of two semi-infinite domains and that of two finite-sized domains are examined in detail. We indicate also how to extend the solution method to the setting of one finite-sized domain surrounded on both sides by semi-infinite domains, and on that of three finite-sized domains.

Orthogonal packing of identical rectangles within isotropic convex regions

A mixed integer continuous nonlinear model and a solution method for the problem of orthogonally packing identical rectangles within an arbitrary convex region are introduced in the present work. The convex region is assumed to be made of an isotropic material in such a way that arbitrary rotations of the items, preserving the orthogonality constraint, are allowed. The solution method is based on a combination of branch and bound and active-set strategies for bound-constrained minimization of smooth functions. Numerical results show the reliability of the presented approach. (C) 2010 Elsevier Ltd. All rights reserved.

Essays on decision sciences: exploring cognition, information processing, and complexity

This thesis provides three original contributions to the field of Decision Sciences. The first contribution explores the field of heuristics and biases. New variations of the Cognitive Reflection Test (CRT--a test to measure "the ability or disposition to resist reporting the response that first comes to mind"), are provided. The original CRT (S. Frederick [2005] Journal of Economic Perspectives, v. 19:4, pp.24-42) has items in which the response is immediate--and erroneous. It is shown that by merely varying the numerical parameters of the problems, large deviations in response are found. Not only the final results are affected by the proposed variations, but so is processing fluency. It seems that numbers' magnitudes serve as a cue to activate system-2 type reasoning. The second contribution explores Managerial Algorithmics Theory (M. Moldoveanu [2009] Strategic Management Journal, v. 30, pp. 737-763); an ambitious research program that states that managers display cognitive choices with a "preference towards solving problems of low computational complexity". An empirical test of this hypothesis is conducted, with results showing that this premise is not supported. A number of problems are designed with the intent of testing the predictions from managerial algorithmics against the predictions of cognitive psychology. The results demonstrate (once again) that framing effects profoundly affect choice, and (an original insight) that managers are unable to distinguish computational complexity problem classes. The third contribution explores a new approach to a computationally complex problem in marketing: the shelf space allocation problem (M-H Yang [2001] European Journal of Operational Research, v. 131, pp.107--118). A new representation for a genetic algorithm is developed, and computational experiments demonstrate its feasibility as a practical solution method. These studies lie at the interface of psychology and economics (with bounded rationality and the heuristics and biases programme), psychology, strategy, and computational complexity, and heuristics for computationally hard problems in management science.

Envelhecimento acelerado em sementes de alface

Objetivou-se adequar a metodologia do teste de envelhecimento acelerado para avaliação do potencial fisiológico de sementes de quatro lotes dos genótipos de alface 'Babá de verão', 'Tainá' e 'Vera', aplicando-se o método tradicional e com o uso de solução saturada de NaCl. Para tal, as sementes foram submetidas à determinação do teor de água, testes de germinação e vigor (primeira contagem, condutividade elétrica, emergência em bandeja, índice de velocidade de emergência) incluindo o teste de envelhecimento acelerado, conduzido nos períodos de 24, 48, e 72 horas, a 41°C com e sem solução saturada de NaCl. O experimento foi conduzido em delineamento inteiramente casualizado com quatro repetições. Os procedimentos tradicional e com o uso da solução saturada de NaCl por 72 horas foram sensíveis em estratificar os lotes em níveis de vigor, além de estarem de acordo com os demais testes utilizados. Portanto, o envelhecimento acelerado de sementes de alface deve ser feito utilizando o período de 72 horas, podendo ser realizado tanto pelo método tradicional quanto pelo método da solução saturada de NaCl.

Estudo da cristalização de parafinas em sistemas solventes/tensoativos/água

The WAT is the temperature at the beginning of the appearance of wax crystals. At this temperature the first wax crystals are formed by the cooling systems paraffin / solvents. Paraffins are composed of a mixture of saturated hydrocarbons of high molecular weight. The removal of petroleum from wells and the production lines means a surcharge on produced oil, thus solubilize these deposits formed due to modifications of thermodynamics has been a constant challenge for companies of oil exploration. This study combines the paraffin solubilization by microemulsion systems, the determination of WAT systems paraffin / solvent and performance of surfactant in reducing the crystallization. We used the methods: rheological and the photoelectric signal, validating the latter which was developed to optimize the data obtained due to sensitivity of the equipment used. Methods developed for description of wax precipitation are often in poor agreement with the experimental data, they tend to underestimate the amount of wax at temperatures below the turbidity point. The Won method and the Ideal solution method were applied to the WAT data obtained in solvent systems, best represented by the second interaction of Won method using the solvents naphtha, hexane and LCO. It was observed that the results obtained by WAT photoelectric signal when compared with the viscosity occur in advance, demonstrating the greatest sensitivity of the method developed. The ionic surfactant reduced the viscosity of the solvent systems as it acted modifying the crystalline structure and, consequently, the pour point. The curves show that the WAT experimental data is, in general, closer to the modeling performed by the method of Won than to the one performed by the ideal solution method, because this method underestimates the curve predicting the onset of paraffin hydrocarbons crystallization temperature. This occurs because the actual temperature measured was the crystallization temperature and the method proposes the fusion temperature measurement.

Estudo da digestão de sedimento de curtume visando a extração do cromo por microemulsão

Leather tanneries generate effluents with high content of heavy metals, especially chromium, which is used in the mineral tanning process. Microemulsions have been studied in the extraction of heavy metals from aqueous solutions. Considering the problems related with the sediment resulting from the tanning process, due to its high content in chromium, in this work this sediment was characterized and microemulsion systems were applied for chromium removal. The extraction process consists in the removal of heavy metal ions present in an aqueous feeding solution (acid digestion solution) by a microemulsion system. First three different solid sludge digestion methods were evaluated, being chosen the method with higher digestion capacity. For this digestion method, seeking its optimization, was evaluated the influence of granule size, temperature and digestion time. Experimental results showed that the method proposed by USEPA (Method A) was the most efficient one, being obtained 95.77% of sample digestion. Regarding to the evaluated parameters, the best results were achieved at 95°C, 14 Mesh granule size, and 60 minutes digestion time. For chromium removal, three microemulsion extraction methods were evaluated: Method 1, in a Winsor II region, using as aqueous phase the acid digestion solution; Method 2, in a Winsor IV region, being obtained by the addition of the acid digestion solution to a microemulsion phase, whose aqueous phase is distilled water, until the formation of Winsor II system; and Method 3, in a Winsor III region, consisting in the formation of a Winsor III region using as aqueous phase the acid digestion solution, diluted in NaOH 0.01N. Seeking to optimize the extraction process only Method 1 (Systems I, II, and VIII) and Method 2 (System IX) were evaluated, being chosen points inside the interest regions (studied domains) to study the influence of contact time and pH in the extraction percentiles. The studied systems present the following compositions: System I: Surfactant Saponified coconut oil, Cosurfactant 1-Butanol, Oil phase Kerosene, Aqueous phase 2% NaCl solution; System II: Aqueous phase Acid digestion solution with pH adjusted using KOH (pH 3.5); System VIII: Aqueous phase - Acid digestion solution (pH 0.06); and System IX Aqueous phase Distilled water (pH 10.24), the other phases of Systems II, VIII and IX are similar to System I. Method 2 showed to be the more efficient one regarding chromium extraction percentile (up to 96.59% - pH 3.5). Considering that with Method 2 the microemulsion region only appears in the Winsor II region, it was studied Method 3 (System X) for the evaluation and characterization of a triphasic system, seeking to compare with a biphases system. System X is composed by: Surfactant Saponified coconut oil, Cosurfactant 1-Butanol, Oil phase Kerosene, Aqueous phase Acid digestion solution diluted with water and with its pH adjusted using 0.01N NaOH solution. The biphasic and triphasic microemulsion systems were analyzed regarding its viscosity, extraction efficiency and drop effective diameter. The experimental results showed that for viscosity studies the obtained values were low for all studied systems, the diameter of the drop is smaller in the Winsor II region, with 15.5 nm, reaching 46.0 nm in Winsor III region, being this difference attributed to variations in system compositions and micelle geometry. In chromium extraction, these points showed similar results, being achieved 99.76% for Winsor II system and 99.62% for Winsor III system. Winsor III system showed to be more efficient due to the obtaining of a icroemulsion with smaller volume, with the possibility to recover the oil phase in excess, and the use of a smaller proportion of surfactant and cosurfactant (C/S)

C programs for solving the time-dependent Gross-Pitaevskii equation in a fully anisotropic trap

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Inferência fuzzy para o problema de corte de estoque com sobras aproveitáveis de material

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Doping of polyaniline by corona discharge

In the present article it is shown that a corona discharge can be employed to dope thin films of polyaniline (PANI) coated on poly(ethylene terephthalate) films, allowing the electrical conductivity to be tuned within the range 10(-10) to 0.3 S cm(-1). A study of the effect of different corona conditions, namely corona treatment for positive and negative polarities, air humidity, treatment time, corona current, and the geometry of the corona triode, on the electrical conductivity of the polyaniline is presented. The results indicate that the corona discharge leads to protonic doping of polyaniline similar to that which occurs in conventional protonic acid solution doping. Atomic force microscopic analysis shows that, as the PANI is exposed to the corona discharge, its globular morphology is disrupted leading to the appearance of droplet-like features and a significant decrease in the average height and surface roughness. Doping by corona discharge presents several advantages over the conventional solution method namely that it is a dry process which does not require use of chemicals reagents, and which is both rapid and avoids dopant migration. The latter can be important for applications of PANI in microelectronic devices. (C) 2000 American Institute of Physics. [S0021-8979(00)01608-X].