Quality potential of DIP and hardwood CTMP based raw materials in woodcontaining printing paper

The objective of the thesis was to define the quality potential of DIP and hardwood CTMP based raw material furnish for a printing paper production and to define the end product's pulp-based boundary conditions especially when thinking of Chinese markets. Although the Chinese paper industry expands rapidly, the production of some paper grades is still exiguous. Especially the softwood resources are limited in Asia, thus the purpose of the thesis was to find out the possibilities to produce printing paper in China from local raw materials. Bleached CTMP can be produced, for example, from fast-growing hardwood species like eucalyptus and poplar. Therefore in this thesis it was examined if good quality printing paper is possible to produce by using deinked pulp and hardwood CTMP based furnish. In the first section of experimental part, various deinked pulps and chemithermomechanical pulps were compared. The deinked pulps were from China, Central Europe and Finland. Central European was made for magazine papers, and the Chinese as well as the Finnish pulps were made for newsprints. Two of the BCTMPs were from China and those both were made from poplar, whereas one BCTMP was made from eucalyptus in a pilot plant. There were significant differences especially between BCTMPs and their paper properties. In the second section of experimental part, the deinked pulp and eucalyptus BCTMP were blended to produce handsheets. The results show that producing the highest quality printing paper would be difficult from these raw materials. Deinked pulp affected especially the strength and optical properties as well as calender blackening. The BCTMP was found to have effects mostly on the smoothness, strength and optical properties as well as calender blackening.

Seismic characterization of heterogeneous sedimentary and fractured rocks based on Biot's theory of poroelasticity

Understanding and quantifying seismic energy dissipation, which manifests itself in terms of velocity dispersion and attenuation, in fluid-saturated porous rocks is of considerable interest, since it offers the perspective of extracting information with regard to the elastic and hydraulic rock properties. There is increasing evidence to suggest that wave-induced fluid flow, or simply WIFF, is the dominant underlying physical mechanism governing these phenomena throughout the seismic, sonic, and ultrasonic frequency ranges. This mechanism, which can prevail at the microscopic, mesoscopic, and macroscopic scale ranges, operates through viscous energy dissipation in response to fluid pressure gradients and inertial effects induced by the passing wavefield. In the first part of this thesis, we present an analysis of broad-band multi-frequency sonic log data from a borehole penetrating water-saturated unconsolidated glacio-fluvial sediments. An inherent complication arising in the interpretation of the observed P-wave attenuation and velocity dispersion is, however, that the relative importance of WIFF at the various scales is unknown and difficult to unravel. An important generic result of our work is that the levels of attenuation and velocity dispersion due to the presence of mesoscopic heterogeneities in water-saturated unconsolidated clastic sediments are expected to be largely negligible. Conversely, WIFF at the macroscopic scale allows for explaining most of the considered data while refinements provided by including WIFF at the microscopic scale in the analysis are locally meaningful. Using a Monte-Carlo-type inversion approach, we compare the capability of the different models describing WIFF at the macroscopic and microscopic scales with regard to their ability to constrain the dry frame elastic moduli and the permeability as well as their local probability distribution. In the second part of this thesis, we explore the issue of determining the size of a representative elementary volume (REV) arising in the numerical upscaling procedures of effective seismic velocity dispersion and attenuation of heterogeneous media. To this end, we focus on a set of idealized synthetic rock samples characterized by the presence of layers, fractures or patchy saturation in the mesocopic scale range. These scenarios are highly pertinent because they tend to be associated with very high levels of velocity dispersion and attenuation caused by WIFF in the mesoscopic scale range. The problem of determining the REV size for generic heterogeneous rocks is extremely complex and entirely unexplored in the given context. In this pilot study, we have therefore focused on periodic media, which assures the inherent self- similarity of the considered samples regardless of their size and thus simplifies the problem to a systematic analysis of the dependence of the REV size on the applied boundary conditions in the numerical simulations. Our results demonstrate that boundary condition effects are absent for layered media and negligible in the presence of patchy saturation, thus resulting in minimum REV sizes. Conversely, strong boundary condition effects arise in the presence of a periodic distribution of finite-length fractures, thus leading to large REV sizes. In the third part of the thesis, we propose a novel effective poroelastic model for periodic media characterized by mesoscopic layering, which accounts for WIFF at both the macroscopic and mesoscopic scales as well as for the anisotropy associated with the layering. Correspondingly, this model correctly predicts the existence of the fast and slow P-waves as well as quasi and pure S-waves for any direction of wave propagation as long as the corresponding wavelengths are much larger than the layer thicknesses. The primary motivation for this work is that, for formations of intermediate to high permeability, such as, for example, unconsolidated sediments, clean sandstones, or fractured rocks, these two WIFF mechanisms may prevail at similar frequencies. This scenario, which can be expected rather common, cannot be accounted for by existing models for layered porous media. Comparisons of analytical solutions of the P- and S-wave phase velocities and inverse quality factors for wave propagation perpendicular to the layering with those obtained from numerical simulations based on a ID finite-element solution of the poroelastic equations of motion show very good agreement as long as the assumption of long wavelengths remains valid. A limitation of the proposed model is its inability to account for inertial effects in mesoscopic WIFF when both WIFF mechanisms prevail at similar frequencies. Our results do, however, also indicate that the associated error is likely to be relatively small, as, even at frequencies at which both inertial and scattering effects are expected to be at play, the proposed model provides a solution that is remarkably close to its numerical benchmark. -- Comprendre et pouvoir quantifier la dissipation d'énergie sismique qui se traduit par la dispersion et l'atténuation des vitesses dans les roches poreuses et saturées en fluide est un intérêt primordial pour obtenir des informations à propos des propriétés élastique et hydraulique des roches en question. De plus en plus d'études montrent que le déplacement relatif du fluide par rapport au solide induit par le passage de l'onde (wave induced fluid flow en anglais, dont on gardera ici l'abréviation largement utilisée, WIFF), représente le principal mécanisme physique qui régit ces phénomènes, pour la gamme des fréquences sismiques, sonique et jusqu'à l'ultrasonique. Ce mécanisme, qui prédomine aux échelles microscopique, mésoscopique et macroscopique, est lié à la dissipation d'énergie visqueuse résultant des gradients de pression de fluide et des effets inertiels induits par le passage du champ d'onde. Dans la première partie de cette thèse, nous présentons une analyse de données de diagraphie acoustique à large bande et multifréquences, issues d'un forage réalisé dans des sédiments glaciaux-fluviaux, non-consolidés et saturés en eau. La difficulté inhérente à l'interprétation de l'atténuation et de la dispersion des vitesses des ondes P observées, est que l'importance des WIFF aux différentes échelles est inconnue et difficile à quantifier. Notre étude montre que l'on peut négliger le taux d'atténuation et de dispersion des vitesses dû à la présence d'hétérogénéités à l'échelle mésoscopique dans des sédiments clastiques, non- consolidés et saturés en eau. A l'inverse, les WIFF à l'échelle macroscopique expliquent la plupart des données, tandis que les précisions apportées par les WIFF à l'échelle microscopique sont localement significatives. En utilisant une méthode d'inversion du type Monte-Carlo, nous avons comparé, pour les deux modèles WIFF aux échelles macroscopique et microscopique, leur capacité à contraindre les modules élastiques de la matrice sèche et la perméabilité ainsi que leur distribution de probabilité locale. Dans une seconde partie de cette thèse, nous cherchons une solution pour déterminer la dimension d'un volume élémentaire représentatif (noté VER). Cette problématique se pose dans les procédures numériques de changement d'échelle pour déterminer l'atténuation effective et la dispersion effective de la vitesse sismique dans un milieu hétérogène. Pour ce faire, nous nous concentrons sur un ensemble d'échantillons de roches synthétiques idéalisés incluant des strates, des fissures, ou une saturation partielle à l'échelle mésoscopique. Ces scénarios sont hautement pertinents, car ils sont associés à un taux très élevé d'atténuation et de dispersion des vitesses causé par les WIFF à l'échelle mésoscopique. L'enjeu de déterminer la dimension d'un VER pour une roche hétérogène est très complexe et encore inexploré dans le contexte actuel. Dans cette étude-pilote, nous nous focalisons sur des milieux périodiques, qui assurent l'autosimilarité des échantillons considérés indépendamment de leur taille. Ainsi, nous simplifions le problème à une analyse systématique de la dépendance de la dimension des VER aux conditions aux limites appliquées. Nos résultats indiquent que les effets des conditions aux limites sont absents pour un milieu stratifié, et négligeables pour un milieu à saturation partielle : cela résultant à des dimensions petites des VER. Au contraire, de forts effets des conditions aux limites apparaissent dans les milieux présentant une distribution périodique de fissures de taille finie : cela conduisant à de grandes dimensions des VER. Dans la troisième partie de cette thèse, nous proposons un nouveau modèle poro- élastique effectif, pour les milieux périodiques caractérisés par une stratification mésoscopique, qui prendra en compte les WIFF à la fois aux échelles mésoscopique et macroscopique, ainsi que l'anisotropie associée à ces strates. Ce modèle prédit alors avec exactitude l'existence des ondes P rapides et lentes ainsi que les quasis et pures ondes S, pour toutes les directions de propagation de l'onde, tant que la longueur d'onde correspondante est bien plus grande que l'épaisseur de la strate. L'intérêt principal de ce travail est que, pour les formations à perméabilité moyenne à élevée, comme, par exemple, les sédiments non- consolidés, les grès ou encore les roches fissurées, ces deux mécanismes d'WIFF peuvent avoir lieu à des fréquences similaires. Or, ce scénario, qui est assez commun, n'est pas décrit par les modèles existants pour les milieux poreux stratifiés. Les comparaisons des solutions analytiques des vitesses des ondes P et S et de l'atténuation de la propagation des ondes perpendiculaires à la stratification, avec les solutions obtenues à partir de simulations numériques en éléments finis, fondées sur une solution obtenue en 1D des équations poro- élastiques, montrent un très bon accord, tant que l'hypothèse des grandes longueurs d'onde reste valable. Il y a cependant une limitation de ce modèle qui est liée à son incapacité à prendre en compte les effets inertiels dans les WIFF mésoscopiques quand les deux mécanismes d'WIFF prédominent à des fréquences similaires. Néanmoins, nos résultats montrent aussi que l'erreur associée est relativement faible, même à des fréquences à laquelle sont attendus les deux effets d'inertie et de diffusion, indiquant que le modèle proposé fournit une solution qui est remarquablement proche de sa référence numérique.

The ESRg matrix for strong field d5 systems

This review has tried to collect and correlate all the various equations for the g matrix of strong field d5 systems obtained from different basis sets using full electron and hole formalism calculations. It has corrected mistakes found in the literature and shown how the failure to properly take in symmetry boundary conditions has produced a variety of apparently inconsistent equations in the literature. The review has reexamined the problem of spin-orbit interaction with excited t4e states and finds that the earlier reports that it is zero in octahedral symmetry is not correct. It has shown how redefining what x, y, and z are in the principal coordinate system simplifies, compared to previous methods, the analysis of experimental g values with the equations.

Modelo de comportamiento en sevicio para estructuras reticulares de hormigón. Parte I: Descripción del modelo y ejemplos de aplicación

A general method for instantaneous and time-dependent serviceability analysis of plane concrete frames is presented. The methodology is based in an extension of the classic matrix formulation for bars. The main aspects influencing the behaviour of the structural concrete are considered: cracking, creep, shrinkage or prestress losses. To simulate the effect of cracking a smeared model (developed in Part II) based on the modification of the tensile branch of the concrete stress-strain relationship is adopted. The general approach considered permits the application to different materials and constitutive laws. Sequential construction (sectional and structural), incorporation of reinforcement, consideration of the loads history; placing and removing shores, and restraining or releasing in boundary conditions are considered. Some examples are included to highlight the capabilities of the model

Vaasan Sähköverkko Oy:n keskijänniteverkon nykytilan määritys sekä kehittämissuunnitelma käyttövarmuuden näkökulmasta

Tämän työn tavoitteena on analysoida Vaasan Sähköverkko Oy:n keskijänniteverkon nykytila sekä selvittää verkon kehittämisen suuntaviivoja. Nykytilan määrityksessä selvitetään kuormitettavuusasteet, jännitteenalenemat, häviöt oiko- ja maasulkusuojauksen toimivuus sekä keskeytyskustannukset. Nykytilan määrityksen apuna käytettiin verkkotietojärjestelmä Xpowerin tehojako-, oikosulku- ja maasulkulaskentoja. Nykytila on määritetty sekä normaalissa kytkentätilanteessa, että tilanteessa jossa yksi sähköasema kerrallaan on poissa käytöstä. Nykytilan määrityksen perusteella voidaan todeta, että normaalissa kytkentätilanteessa keskijänniteverkko täyttää suurilta osin sille asetetut tekniset reunaehdot, jotka ovat jännitteenalenema, kuormitettavuus, oikosulkukestoisuus sekä oiko- ja maasulkusuojauksen havahtuvaisuus. Jotta reunaehdot toteutuisivat kokonaisuudessaan, tulee oikosulkukestottomien lähtöjen releasetteluja muuttaa tai vaihtaa oikosulkukestottomat johto-osuudet suurempi poikkipintaiseen johtimeen. Nykytilan määrityksen perusteella neljän päämuuntajan kuormitettavuusaste on yli 90 %, joten näiden asemien alueiden kehitystä tulee tarkkailla. Sähköasemien korvaustarkastelun perusteella viiden sähköaseman korvaus tuottaa ongelmia jo nykyisillä huippukuormilla. Työssä tehtyjen selvitystöiden perusteella Vaasan Sähköverkko Oy:n keskijänniteverkkoon lisättiin 10 maastokatkaisijaa sekä selvitettiin verkon mahdollisia tulevaisuuden kehittämisen suuntaviivoja.

Leaving the Dark ages with AMIGA

We present an Analytic Model of Intergalactic-medium and GAlaxy (AMIGA) evolution since the dark ages. AMIGA is in the spirit of the popular semi-analytic models of galaxy formation, although it does not use halo merger trees but interpolates halo properties in grids that are progressively built. This strategy is less memory-demanding and allows one to start modeling at sufficiently high redshifts and low halo masses to have trivial boundary conditions. The number of free parameters is minimized by making a causal connection between physical processes usually treated as independent of each other, which leads to more reliable predictions. However, the strongest points of AMIGA are the following: (1) the inclusion of molecular cooling and metal-poor, population III (Pop III) stars with the most dramatic feedback and (2) accurate follow up of the temperature and volume filling factor of neutral, singly ionized, and doubly ionized regions, taking into account the distinct halo mass functions in those environments. We find the following general results. Massive Pop III stars determine the intergalactic medium metallicity and temperature, and the growth of spheroids and disks is self-regulated by that of massive black holes (MBHs) developed from the remnants of those stars. However, the properties of normal galaxies and active galactic nuclei appear to be quite insensitive to Pop III star properties due to the much higher yield of ordinary stars compared to Pop III stars and the dramatic growth of MBHs when normal galaxies begin to develop, which cause the memory loss of the initial conditions.

Jäykistetyn laattarakenteen mitoittaminen suurten taipumien teorialla

Tässä diplomityössä tutkitaan epälineaarisen teorian hyödyntämistä laattarakenteiden analysoinnissa ja pyritään muodostamaan yksinkertaisia mitoitusohjeita laattarakenteille. Analysointia varten laattarakenne jaetaan lohkoihin, joille mitoitussäännöt määritetään. Tarkastellaan myös lineaarisen ja epälineaarisen teorian antamien tulosten eroavai-suutta ja suurten taipumien huomioimisen aiheuttamaa hyötyä. Työn lähtökohtana on standardi SFS-EN 1993-1-7, joka antaa yksinkertaistettuja mitoitusohjeita laattarakenteiden mitoitukseen lineaarisella ja epälineaarisella teorialla. Työssä selvitetään mistä nämä mitoitusohjeet ovat peräisin. Lisäksi tarkastellaan laattalohkojen erilaisia liitoksia, jotta voidaan tehdä oletuksia laatan reunaehdoista. Tulokseksi saadaan suuntaa-antavia mitoitusohjeita, joiden soveltaminen onnistuu vain muutamille yksinkertaisille rakenteille. Mitoitusohjeiden perusteella saadaan laskettua laattalohkon jännitykset, taipuma sekä mitoitusohjeiden vaatiman jäykisteen koko. Jotta mitoitusohjeista saataisiin tarkemmat, tulisi mitoitettava rakenne pystyä määrittelemään yksityiskohtaisemmin. Vertailtaessa lineaarista ja epälineaarista laskentaa havaitaan, että huomioimalla suurten taipumien aikaansaamat kalvovoimat, laskevat laattarakenteeseen vaikuttavat taipumat ja jännitykset merkittävästikin.

Estudo teórico das transições eletrônicas usando métodos simples e sofisticados

In this paper, the use of both simple and sophisticated models in the study of electronic transitions was explored for a set of molecular systems: C2H4, C4H4, C4H6, C6H6, C6H8, "C8", C60, and [H2NCHCH(CHCH)kCHNH2]+, where k = 0 to 4. The simple model of the free particle (1D, 2D, and 3D boxes, rings or spherical surfaces), considering the boundary conditions, was found to yield similar results to the sophisticated theoretical methods such as EOM-CCSD/6-311++G** or TD(NStates=5,Root=1)-M06-2X/6-311++G**.

Fatigue analysis of brackets with and without scallop for bridge and other deck structure using effective notch stress.

Due to functional requirement of a structural detail brackets with and without scallop are frequently used in bridges, decks, ships and offshore structure. Scallops are designed to serve as passage way for fluids, to reduce weld length and plate distortions. Moreover, scallops are used to avoid intersection of two or more welds for the fact that there is the presence of inventible inherent initial crack except for full penetrated weld and the formation of multi-axial stress state at the weld intersection. Welding all around the scallop corner increase the possibility of brittle fracture even for the case the bracket is not loaded by primary load. Avoiding of scallop will establish an initial crack in the corner if bracket is welded by fillet welds. If the two weld run pass had crossed, this would have given a 3D residual stress situation. Therefore the presences and absence of scallop necessitates the 3D FEA fatigue resistance of both types of brackets using effective notch stress approach ( ). FEMAP 10.1 with NX NASTRAN was used for the 3D FEA. The first and main objective of this research was to investigate and compare the fatigue resistance of brackets with and without scallop. The secondary goal was the fatigue design of scallops in case they cannot be avoided for some reason. The fatigue resistance for both types of brackets was determined based on approach using 1 mm fictitiously rounded radius based on IIW recommendation. Identical geometrical, boundary and loading conditions were used for the determination and comparison of fatigue resistance of both types of brackets using linear 3D FEA. Moreover the size effect of bracket length was also studied using 2D SHELL element FEA. In the case of brackets with scallop the flange plate weld toe at the corner of the scallop was found to exhibit the highest and made the flange plate weld toe critical for fatigue failure. Whereas weld root and weld toe at the weld intersections were the highly stressed location for brackets without scallop. Thus weld toe for brackets with scallop, and weld root and weld toe for brackets without scallop were found to be the critical area for fatigue failure. Employing identical parameters on both types of brackets, brackets without scallop had the highest except for full penetrated weld. Furthermore the fatigue resistance of brackets without scallop was highly affected by the lack of weld penetration length and it was found out that decreased as the weld penetration was increased. Despite the fact that the very presence of scallop reduces the stiffness and also same time induce stress concentration, based on the 3D FEA it is worth concluding that using scallop provided better fatigue resistance when both types of brackets were fillet welded. However brackets without scallop had the highest fatigue resistance when full penetration weld was used. This thesis also showed that weld toe for brackets with scallop was the only highly stressed area unlike brackets without scallop in which both weld toe and weld root were the critical locations for fatigue failure when different types of boundary conditions were used. Weld throat thickness, plate thickness, scallop radius, lack of weld penetration length, boundary condition and weld quality affected the fatigue resistance of both types of brackets. And as a result, bracket design procedure, especially welding quality and post weld treatment techniques significantly affect the fatigue resistance of both type of brackets.

Ultrakevyen rahtialuksen kansiluukkujen sekä luukkujen siirto- ja varastointijärjestelmän suunnittelu hyödyntäen kevytrakennekonsepteja

Tämän työn tarkoituksena oli tutkia miten rahtialuksen kansiluukut voitaisiin valmistaa mahdollisimman kevyiksi. Katettavan ruuman pinta-ala on n. 10 m x 40 m. Luukkujen suuresta jännevälistä johtuen, rakenteelta vaaditaan suurta jäykkyyttä. Erilaisina vaihtoehtoina tutkittiin vaahtomaista alumiinia, alumiinisia kennorakenteita ja polyuretaanisia sandwich-rakenteita. Työssä vertailtiin myös erilaisia geometrisia ratkaisuja, joilla kansiluukkujen jäykkyyttä pyrittiin lisäämään ja sitä kautta pääsemään pienempään materiaalin tarpeeseen. Geometriaa suunniteltaessa huomioitiin myös vaikutukset ruuman tilavuuteen ja lainsäädännön asettamat reunaehdot. Lainsäädännöstä saatiin esimerkiksi turvakaiteiden minimikorkeus, joka vaikuttaa suoraan ruuman tilavuuteen, kun aluksen korkeimmaksi kohdaksi on valittu laivan keskilinja ja tämä korkeus on annettu suunnittelun lähtötietona. Tietokoneavusteisen lujuuslaskennan avulla eri vaihtoehdoista muodostettiin elementtimallit. Malleja varioimalla ja tuloksia vertailemalla saatiin selville kevyin mahdollinen rakenne ja geometria. Malleista saatiin selville myös luukkujen tukireaktiovoimat, eli voimat, jotka luukut kohdistavat aluksen muihin rakenteisiin. Lisäksi työssä mietittiin erilaisia tapoja ruuman avaamiseen ja avaamistavan vaikutusta kansiluukkujen painoon, geometriaan ja ruuman tilavuuteen.

Permanent Magnet Synchronous Machine for Parallel Hybrid Vehicle

The aim of this thesis is to describe hybrid drive design problems, the advantages and difficulties related to the drive. A review of possible hybrid constructions, benefits of parallel, series and series-parallel hybrids is done. In the thesis analytical and finite element calculations of permanent magnet synchronous machines with embedded magnets were done. The finite element calculations were done using Cedrat’s Flux 2D software. This machine is planned to be used as a motor-generator in a low power parallel hybrid vehicle. The boundary conditions for the design were found from Lucas-TVS Ltd., India. Design Requirements, briefly: • The system DC voltage level is 120 V, which implies Uphase = 49 V (RMS) in a three phase system. • The power output of 10 kW at base speed 1500 rpm (Torque of 65 Nm) is desired. • The maximum outer diameter should not be more than 250 mm, and the maximum core length should not exceed 40 mm. The main difficulties which the author met were the dimensional restrictions. After having designed and analyzed several possible constructions they were compared and the final design selected. Dimensioned and detailed design is performed. Effects of different parameters, such as the number of poles, number of turns and magnetic geometry are discussed. The best modification offers considerable reduction of volume.

Design of marine generators for alternative diesel-electric power systems

This study compares different electric propulsion systems. Results of the analysis of all the advantages and disadvantages of the different propulsion systems are given. This thesis estimates possibilities to apply different diesel-electric propulsion concepts for different vessel types. Small and medium size vessel’s power ranges are studied. The optimal delivery system is chosen. This choice is made on the base of detailed study of the concepts, electrical equipment market and comparison of mass, volume and efficiency parameters. In this thesis three marine generators are designed. They are: salient pole synchronous generator and two permanent magnet synchronous generators. Their electrical, dimensional, cost and efficiency parameters are compared. To understand all the benefits diagrams with these parameters are prepared. Possible benefits and money savings are estimated. As the result the advantages, disadvantages and boundary conditions for the permanent magnet synchronous generator application in marine electric-power systems are found out.

Miehistönkuljetusajoneuvon rungon suunnittelu

Tämän diplomityön tavoitteena oli suunnitella miehistönkuljetusajoneuvon runko. Rungosta suunniteltiin mahdollisimman hyvin energiaa absorboiva. Rakenne toteutettiin kennora-kenteena. Suunnittelussa sovellettiin koneensuunnittelun periaatteiden lisäksi energiaa ab-sorboivien rakenteiden suunnittelun periaatteita. Myös valmistustekniset näkökohdat otet-tiin huomioon. Rakenteessa hyödynnettiin Ruukki Oy:n Ramor 500 suojausterästä sekä OPTIM 500 MC terästä. Lisäksi erilaisten täyteaineiden käyttöä tutkittiin. Suunnittelun työkaluna käytettiin epälineaarista elementtimenetelmää, koska energiaa ab-sorboivien rakenteiden suunnittelussa on otettava huomioon materiaalien epälineaarinen käyttäytyminen. Rakenteen suunnittelu jakaantui viiteen vaiheeseen. Aluksi rakenteeseen kohdistuvat kuormitukset laskettiin elementtimenetelmän avulla. Esisuunnittelussa lasket-tiin plastisuusteorian avulla alustavasti tarvittavat materiaalipaksuudet. Tämän jälkeen ra-kenteen ydingeometria optimoitiin mahdollisimman hyvin energiaa absorboivaksi. Opti-moinnissa hyödynnettiin elementtimenetelmää. Seuraavassa vaiheessa varmistettiin raken-teen globaalit ominaisuudet. Lopuksi rakenteen kestävyyttä tarkasteltiin elementtimene-telmällä. Runko ei mallien mukaan kestänyt siltä vaadittuja kuormitustapauksia. Mallin kaikki ole-tukset pidettiin varmalla puolella. Reunaehdot oletettiin todellisuutta jäykemmiksi. Myös-kään materiaalin venymänopeudesta johtuvaa lujittumista ei otettu huomioon. Koska mii-naräjähdys on monimutkainen tapahtuma, rungon todellinen kestävyys joudutaan ar-viomaan räjähdystesteillä. Elementtimallien perusteella voidaan kuitenkin sanoa, että ener-giaa absorboiva ajoneuvon runko on mahdollista toteuttaa kennorakenteena. Lisäksi voi-daan todeta, että elementtimenetelmää sopii työvälineeksi tämän tyyppisten rakenteiden suunnitteluun.

Application of Computational Fluid Dynamics in Modeling Blood Flow in Human Thoracic Aorta

The aim of this study was to simulate blood flow in thoracic human aorta and understand the role of flow dynamics in the initialization and localization of atherosclerotic plaque in human thoracic aorta. The blood flow dynamics in idealized and realistic models of human thoracic aorta were numerically simulated in three idealized and two realistic thoracic aorta models. The idealized models of thoracic aorta were reconstructed with measurements available from literature, and the realistic models of thoracic aorta were constructed by image processing Computed Tomographic (CT) images. The CT images were made available by South Karelia Central Hospital in Lappeenranta. The reconstruction of thoracic aorta consisted of operations, such as contrast adjustment, image segmentations, and 3D surface rendering. Additional design operations were performed to make the aorta model compatible for the numerical method based computer code. The image processing and design operations were performed with specialized medical image processing software. Pulsatile pressure and velocity boundary conditions were deployed as inlet boundary conditions. The blood flow was assumed homogeneous and incompressible. The blood was assumed to be a Newtonian fluid. The simulations with idealized models of thoracic aorta were carried out with Finite Element Method based computer code, while the simulations with realistic models of thoracic aorta were carried out with Finite Volume Method based computer code. Simulations were carried out for four cardiac cycles. The distribution of flow, pressure and Wall Shear Stress (WSS) observed during the fourth cardiac cycle were extensively analyzed. The aim of carrying out the simulations with idealized model was to get an estimate of flow dynamics in a realistic aorta model. The motive behind the choice of three aorta models with distinct features was to understand the dependence of flow dynamics on aorta anatomy. Highly disturbed and nonuniform distribution of velocity and WSS was observed in aortic arch, near brachiocephalic, left common artery, and left subclavian artery. On the other hand, the WSS profiles at the roots of branches show significant differences with geometry variation of aorta and branches. The comparison of instantaneous WSS profiles revealed that the model with straight branching arteries had relatively lower WSS compared to that in the aorta model with curved branches. In addition to this, significant differences were observed in the spatial and temporal profiles of WSS, flow, and pressure. The study with idealized model was extended to study blood flow in thoracic aorta under the effects of hypertension and hypotension. One of the idealized aorta models was modified along with the boundary conditions to mimic the thoracic aorta under the effects of hypertension and hypotension. The results of simulations with realistic models extracted from CT scans demonstrated more realistic flow dynamics than that in the idealized models. During systole, the velocity in ascending aorta was skewed towards the outer wall of aortic arch. The flow develops secondary flow patterns as it moves downstream towards aortic arch. Unlike idealized models, the distribution of flow was nonplanar and heavily guided by the artery anatomy. Flow cavitation was observed in the aorta model which was imaged giving longer branches. This could not be properly observed in the model with imaging containing a shorter length for aortic branches. The flow circulation was also observed in the inner wall of the aortic arch. However, during the diastole, the flow profiles were almost flat and regular due the acceleration of flow at the inlet. The flow profiles were weakly turbulent during the flow reversal. The complex flow patterns caused a non-uniform distribution of WSS. High WSS was distributed at the junction of branches and aortic arch. Low WSS was distributed at the proximal part of the junction, while intermedium WSS was distributed in the distal part of the junction. The pulsatile nature of the inflow caused oscillating WSS at the branch entry region and inner curvature of aortic arch. Based on the WSS distribution in the realistic model, one of the aorta models was altered to induce artificial atherosclerotic plaque at the branch entry region and inner curvature of aortic arch. Atherosclerotic plaque causing 50% blockage of lumen was introduced in brachiocephalic artery, common carotid artery, left subclavian artery, and aortic arch. The aim of this part of the study was first to study the effect of stenosis on flow and WSS distribution, understand the effect of shape of atherosclerotic plaque on flow and WSS distribution, and finally to investigate the effect of lumen blockage severity on flow and WSS distributions. The results revealed that the distribution of WSS is significantly affected by plaque with mere 50% stenosis. The asymmetric shape of stenosis causes higher WSS in branching arteries than in the cases with symmetric plaque. The flow dynamics within thoracic aorta models has been extensively studied and reported here. The effects of pressure and arterial anatomy on the flow dynamic were investigated. The distribution of complex flow and WSS is correlated with the localization of atherosclerosis. With the available results we can conclude that the thoracic aorta, with complex anatomy is the most vulnerable artery for the localization and development of atherosclerosis. The flow dynamics and arterial anatomy play a role in the localization of atherosclerosis. The patient specific image based models can be used to diagnose the locations in the aorta vulnerable to the development of arterial diseases such as atherosclerosis.

Comparison of diffusion models for description of osmotic dehydration of radish slices dipped in brine

ABSTRACT This paper aims at describing the osmotic dehydration of radish cut into cylindrical pieces, using one- and two-dimensional analytical solutions of diffusion equation with boundary conditions of the first and third kind. These solutions were coupled with an optimizer to determine the process parameters, using experimental data. Three models were proposed to describe the osmotic dehydration of radish slices in brine at low temperature. The two-dimensional model with boundary condition of the third kind well described the kinetics of mass transfers, and it enabled prediction of moisture and solid distributions at any given time.