Computational modeling of stented coronary arteries

The application of computational fluid dynamics (CFD) and finite element analysis (FEA) has been growing rapidly in the various fields of science and technology. One of the areas of interest is in biomedical engineering. The altered hemodynamics inside the blood vessels plays a key role in the development of the arterial disease called atherosclerosis, which is the major cause of human death worldwide. Atherosclerosis is often treated with the stenting procedure to restore the normal blood flow. A stent is a tubular, flexible structure, usually made of metals, which is driven and expanded in the blocked arteries. Despite the success rate of the stenting procedure, it is often associated with the restenosis (re-narrowing of the artery) process. The presence of non-biological device in the artery causes inflammation or re-growth of atherosclerotic lesions in the treated vessels. Several factors including the design of stents, type of stent expansion, expansion pressure, morphology and composition of vessel wall influence the restenosis process. Therefore, the role of computational studies is crucial in the investigation and optimisation of the factors that influence post-stenting complications. This thesis focuses on the stent-vessel wall interactions followed by the blood flow in the post-stenting stage of stenosed human coronary artery. Hemodynamic and mechanical stresses were analysed in three separate stent-plaque-artery models. Plaque was modeled as a multi-layer (fibrous cap (FC), necrotic core (NC), and fibrosis (F)) and the arterial wall as a single layer domain. CFD/FEA simulations were performed using commercial software packages in several models mimicking the various stages and morphologies of atherosclerosis. The tissue prolapse (TP) of stented vessel wall, the distribution of von Mises stress (VMS) inside various layers of vessel wall, and the wall shear stress (WSS) along the luminal surface of the deformed vessel wall were measured and evaluated. The results revealed the role of the stenosis size, thickness of each layer of atherosclerotic wall, thickness of stent strut, pressure applied for stenosis expansion, and the flow condition in the distribution of stresses. The thicknesses of FC, and NC and the total thickness of plaque are critical in controlling the stresses inside the tissue. A small change in morphology of artery wall can significantly affect the distribution of stresses. In particular, FC is the most sensitive layer to TP and stresses, which could determine plaque’s vulnerability to rupture. The WSS is highly influenced by the deflection of artery, which in turn is dependent on the structural composition of arterial wall layers. Together with the stenosis size, their roles could play a decisive role in controlling the low values of WSS (<0.5 Pa) prone to restenosis. Moreover, the time dependent flow altered the percentage of luminal area with WSS values less than 0.5 Pa at different time instants. The non- Newtonian viscosity model of the blood properties significantly affects the prediction of WSS magnitude. The outcomes of this investigation will help to better understand the roles of the individual layers of atherosclerotic vessels and their risk to provoke restenosis at the post-stenting stage. As a consequence, the implementation of such an approach to assess the post-stented stresses will assist the engineers and clinicians in optimizing the stenting techniques to minimize the occurrence of restenosis.

Computational studies of deformation in stainless steel rings due to torch heating

A support ring of AISI 304L stainless steel that holds vertical, parallel wires arranged in a circle forming a cylinder is studied. The wires are attached to the ring with heat-induced shrinkage. When the ring is heated with a torch the heat affected zone tries to expand while the adjacent cool structure obstructs the expansion causing upsetting. During cooling, the ring shrinks smaller than its original size clamping the wires. The most important requirement for the ring is that it should be as round as possible and the deformations should occur as overall shrinkage in the ring diameter. A three-dimensional nonlinear transient sequential thermo-structural Abaqus model is used together with a Fortran code that enters the heat flux to each affected element. The local and overall deformations in one ring inflicted by the heating are studied with a small amount of inspection on residual stresses. A variety of different cases are chosen to be studied with the model constructed to provide directional knowledge; torch flux with the means of speed, location of the wires, heating location and structural factors. The decrease of heating speed increases heat flux that rises the temperature increasing shrinkage. In a single progressive heating uneven distribution of shrinkage appears to the start/end region that can be partially fixed with using speeded heating’s to strengthen the heating of that region. Location of the wires affect greatly to the caused shrinkage unlike heating location. The ring structure affects also greatly to the shrinkage; smaller diameter, bigger ring height, thinner thickness and greater number of wires increase shrinkage.

Computational modelling of non-equilibrium condensing steam flows in low-pressure steam turbines

The steam turbines play a significant role in global power generation. Especially, research on low pressure (LP) steam turbine stages is of special importance for steam turbine man- ufactures, vendors, power plant owners and the scientific community due to their lower efficiency than the high pressure steam turbine stages. Because of condensation, the last stages of LP turbine experience irreversible thermodynamic losses, aerodynamic losses and erosion in turbine blades. Additionally, an LP steam turbine requires maintenance due to moisture generation, and therefore, it is also affecting on the turbine reliability. Therefore, the design of energy efficient LP steam turbines requires a comprehensive analysis of condensation phenomena and corresponding losses occurring in the steam tur- bine either by experiments or with numerical simulations. The aim of the present work is to apply computational fluid dynamics (CFD) to enhance the existing knowledge and understanding of condensing steam flows and loss mechanisms that occur due to the irre- versible heat and mass transfer during the condensation process in an LP steam turbine. Throughout this work, two commercial CFD codes were used to model non-equilibrium condensing steam flows. The Eulerian-Eulerian approach was utilised in which the mix- ture of vapour and liquid phases was solved by Reynolds-averaged Navier-Stokes equa- tions. The nucleation process was modelled with the classical nucleation theory, and two different droplet growth models were used to predict the droplet growth rate. The flow turbulence was solved by employing the standard k-ε and the shear stress transport k-ω turbulence models. Further, both models were modified and implemented in the CFD codes. The thermodynamic properties of vapour and liquid phases were evaluated with real gas models. In this thesis, various topics, namely the influence of real gas properties, turbulence mod- elling, unsteadiness and the blade trailing edge shape on wet-steam flows, are studied with different convergent-divergent nozzles, turbine stator cascade and 3D turbine stator-rotor stage. The simulated results of this study were evaluated and discussed together with the available experimental data in the literature. The grid independence study revealed that an adequate grid size is required to capture correct trends of condensation phenomena in LP turbine flows. The study shows that accurate real gas properties are important for the precise modelling of non-equilibrium condensing steam flows. The turbulence modelling revealed that the flow expansion and subsequently the rate of formation of liquid droplet nuclei and its growth process were affected by the turbulence modelling. The losses were rather sensitive to turbulence modelling as well. Based on the presented results, it could be observed that the correct computational prediction of wet-steam flows in the LP turbine requires the turbulence to be modelled accurately. The trailing edge shape of the LP turbine blades influenced the liquid droplet formulation, distribution and sizes, and loss generation. The study shows that the semicircular trailing edge shape predicted the smallest droplet sizes. The square trailing edge shape estimated greater losses. The analysis of steady and unsteady calculations of wet-steam flow exhibited that in unsteady simulations, the interaction of wakes in the rotor blade row affected the flow field. The flow unsteadiness influenced the nucleation and droplet growth processes due to the fluctuation in the Wilson point.

Three-dimensional model of bubbling fluidized bed combustion

 Diplomityön tarkoituksena on kehittää kolmiulotteinen malli kerrosleijupoltolle. Työn kirjallisuusosa sisältää seuraavat perusteet kerrosleijupolton tekniikasta: yleistiedot, leijutus- ja palamisilmiöt, kiinteän aineen ja kaasun sekoittuminen, päästöt ja lämmönsiirto. Lisäksi palamissysteemin mallinnuksen perusteet ja ratkaisumenetelmät ovat esitelty. Työn mallinnusosassa kehitetty koodi on ohjelmoitu Fortran-ohjelmointikielellä. Kehitetty malli perustuu olemassa olevaan malliin kiertoleijupoltosta. Yhtälö kiintoainekonsentraatioprofiilille on vaihdettu ja kiertovirta on poistettu koodista. Mallilla on tehty herkkyystarkasteluja polttoaineen ja kaasun sekoittumisen sekä reaktiokertoimen vaikutukselle. Visualisointi on tehty ohjelmassa Tecplot 360 ja mallinnustuloksia on vertailtu mitattuihin tuloksiin. Mallin laskemattulokset vastaavat hyvin mittaustuloksia ja kokemusperäisiä tietoja; monissa tapauksissa malli pystyy kvantitatiivisesti kuvaamaan parametrien variointia ja kaikissa tapauksissa malli antaa ainakin kvalitatiivisesti oikeita tuloksia. Työhön liittyvän kehityksen ja mallinnuskokemuksen perusteella on tehty ehdotukset mallin tulevaa kehitystä ja mittauksia varten.

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.

Torque vibration model of axial-flux surface-mounted permanent magnet synchronous machine

In order that the radius and thus ununiform structure of the teeth and otherelectrical and magnetic parts of the machine may be taken into consideration the calculation of an axial flux permanent magnet machine is, conventionally, doneby means of 3D FEM-methods. This calculation procedure, however, requires a lotof time and computer recourses. This study proves that also analytical methods can be applied to perform the calculation successfully. The procedure of the analytical calculation can be summarized into following steps: first the magnet is divided into slices, which makes the calculation for each section individually, and then the parts are submitted to calculation of the final results. It is obvious that using this method can save a lot of designing and calculating time. Thecalculation program is designed to model the magnetic and electrical circuits of surface mounted axial flux permanent magnet synchronous machines in such a way, that it takes into account possible magnetic saturation of the iron parts. Theresult of the calculation is the torque of the motor including the vibrations. The motor geometry and the materials and either the torque or pole angle are defined and the motor can be fed with an arbitrary shape and amplitude of three-phase currents. There are no limits for the size and number of the pole pairs nor for many other factors. The calculation steps and the number of different sections of the magnet are selectable, but the calculation time is strongly depending on this. The results are compared to the measurements of real prototypes. The permanent magnet creates part of the flux in the magnetic circuit. The form and amplitude of the flux density in the air-gap depends on the geometry and material of the magnetic circuit, on the length of the air-gap and remanence flux density of the magnet. Slotting is taken into account by using the Carter factor in the slot opening area. The calculation is simple and fast if the shape of the magnetis a square and has no skew in relation to the stator slots. With a more complicated magnet shape the calculation has to be done in several sections. It is clear that according to the increasing number of sections also the result will become more accurate. In a radial flux motor all sections of the magnets create force with a same radius. In the case of an axial flux motor, each radial section creates force with a different radius and the torque is the sum of these. The magnetic circuit of the motor, consisting of the stator iron, rotor iron, air-gap, magnet and the slot, is modelled with a reluctance net, which considers the saturation of the iron. This means, that several iterations, in which the permeability is updated, has to be done in order to get final results. The motor torque is calculated using the instantaneous linkage flux and stator currents. Flux linkage is called the part of the flux that is created by the permanent magnets and the stator currents passing through the coils in stator teeth. The angle between this flux and the phase currents define the torque created by the magnetic circuit. Due to the winding structure of the stator and in order to limit the leakage flux the slot openings of the stator are normally not made of ferromagnetic material even though, in some cases, semimagnetic slot wedges are used. In the slot opening faces the flux enters the iron almost normally (tangentially with respect to the rotor flux) creating tangential forces in the rotor. This phenomenon iscalled cogging. The flux in the slot opening area on the different sides of theopening and in the different slot openings is not equal and so these forces do not compensate each other. In the calculation it is assumed that the flux entering the left side of the opening is the component left from the geometrical centre of the slot. This torque component together with the torque component calculated using the Lorenz force make the total torque of the motor. It is easy to assume that when all the magnet edges, where the derivative component of the magnet flux density is at its highest, enter the slot openings at the same time, this will have as a result a considerable cogging torque. To reduce the cogging torquethe magnet edges can be shaped so that they are not parallel to the stator slots, which is the common way to solve the problem. In doing so, the edge may be spread along the whole slot pitch and thus also the high derivative component willbe spread to occur equally along the rotation. Besides forming the magnets theymay also be placed somewhat asymmetric on the rotor surface. The asymmetric distribution can be made in many different ways. All the magnets may have a different deflection of the symmetrical centre point or they can be for example shiftedin pairs. There are some factors that limit the deflection. The first is that the magnets cannot overlap. The magnet shape and the relative width compared to the pole define the deflection in this case. The other factor is that a shifting of the poles limits the maximum torque of the motor. If the edges of adjacent magnets are very close to each other the leakage flux from one pole to the other increases reducing thus the air-gap magnetization. The asymmetric model needs some assumptions and simplifications in order to limit the size of the model and calculation time. The reluctance net is made for symmetric distribution. If the magnets are distributed asymmetrically the flux in the different pole pairs will not be exactly the same. Therefore, the assumption that the flux flows from the edges of the model to the next pole pairs, in the calculation model from one edgeto the other, is not correct. If it were wished for that this fact should be considered in multi-pole pair machines, this would mean that all the poles, in other words the whole machine, should be modelled in reluctance net. The error resulting from this wrong assumption is, nevertheless, irrelevant.

An Experimental Model of Prostatic Inflammation for Drug Discovery

There is increasing evidence to support a significant role for chronic non-bacterial, prostatic inflammation in the development of human voiding dysfunction and prostate cancer. Their increased prevalence with age suggests that the decrease of testosterone concentration and/or the ratio of testosterone-to-estradiol in serum may have a role in their development. The main objective of this study was to explore prostatic inflammation and its relationship with voiding dysfunction and prostate carcinogenesis by developing an experimental model. A novel selective estrogen receptor modulator (SERM), fispemifene, was tested for the prevention and treatment of prostatic inflammation in this model. Combined treatment of adult Noble rats with testosterone and estradiol for 3 to 6 weeks induced gradually developing prostatic inflammation in the dorsolateral prostatic lobes. Inflammatory cells, mainly T-lymphocytes, were first seen around capillaries. Thereafter, the lymphocytes migrated into the stroma and into periglandular space. When the treatment time was extended to 13 weeks, the number of inflamed acini increased. Urodynamical recordings indicated voiding dysfunction. When the animals had an above normal testosterone and estradiol concentrations but still had a decreased testosterone-to-estradiol ratio in serum, they developed obstructive voiding. Furthermore, they developed precancerous lesions and prostate cancers in the ducts of the dorsolateral prostatic lobes. Interestingly, inflammatory infiltrates were observed adjacent to precancerous lesions but not in the adjacency of adenocarcinomas suggesting that inflammation has a role in the early stages of prostate carcinogenesis. Fispemifene, a novel SERM tested in this experimental model, showed anti-inflammatory action by attenuating the number of inflamed acini in the dorsolateral prostate. Fispemifene exhibited also antiestrogenic properties by decreasing expression of estrogen-induced biomarkers in the acinar epithelium. These findings suggest that SERMs could be considered as a new therapeutic possibility in the prevention and in the treatment of chronic prostatic inflammation

Transgenic mice overexpressing neuropeptide Y: An experimental model of metabolic and cardiovascular diseases

Neuropeptide Y (NPY) is an abundant neurotransmitter in the brain and sympathetic nervous system (SNS). Hypothalamic NPY is known to be a key player in food intake and energy expenditure. NPY’s role in cardiovascular regulation has also been shown. In humans, a Leucine 7 to Proline 7 single nucleotide polymorphism (p.L7P) in the signal peptide of the NPY gene has been associated with traits of metabolic syndrome. The p.L7P subjects also show increased stress-related release of NPY, which suggests that more NPY is produced and released from SNS. The main objective of this study was to create a novel mouse model with noradrenergic cell-targeted overexpression of NPY, and to characterize the metabolic and vascular phenotype of this model. The mouse model was named OE-NPY^DBH mouse. Overexpression of NPY in SNS and brain noradrenergic neurons led to increased adiposity without significant weight gain or increased food intake. The mice showed lipid accumulation in the liver at young age, which together with adiposity led to impaired glucose tolerance and hyperinsulinemia with age. The mice displayed stress-related increased mean arterial blood pressure, increased plasma levels of catecholamines and enhanced SNS activity measured by GDP binding activity to brown adipose tissue mitochondria. Sexual dimorphism in NPY secretion pattern in response to stress was also seen. In an experimental model of vascular injury, the OE-NPY^DBH mice developed more pronounced neointima formation compared with wildtype controls. These results together with the clinical data indicate that NPY in noradrenergic cells plays an important role in the pathogenesis of metabolic syndrome and related diseases. Furthermore, new insights on the role of the extrahypothalamic NPY in the process have been obtained. The OE-NPY^DBH model provides an important tool for further stress and metabolic syndrome-related studies.

Borrelia burgdorferi evades the effects of ceftriaxone treatment in a mouse model of Lyme borreliosis

Borrelia burgdorferi infektoitujen hiirten antibioottihoidon jälkeinen oireilu Lymen borrelioosi on puutiaisten välittämä monimuotoinen infektiotauti, jonka tunnetuin oire on ns. vaeltava ihottuma eli erythema migrans. Muita tavallisia ilmentymiä ovat erityisesti nivel- ja hermosto-oireet sekä harvemmin sydän- ja silmäoireet. Suurin osa potilaista paranee täysin terveeksi antibioottihoidon avulla, mutta jopa 10 % borrelioosiin sairastuneista oireilee suositusten mukaisesta hoidosta huolimatta. Pitkittyneen oireilun on ajateltu johtuvan mm. infektion laukaisemasta autoimmuunitaudista tai kroonisesta infektiosta, mutta teorioiden tueksi ei ole kyetty esittämään kiistattomia todisteita. Onkin todennäköistä, että antibioottihoidon jälkeisen oireilun takana on useampia mekanismeja eikä yksi teoria selitä kaikkien potilaiden oireilua. Tässä väitöskirjatyössä on tutkittu hoidonjälkeistä borrelioosia hiirimallin avulla. Varhaisvaiheessa (2 viikkoa infektoinnin jälkeen) annettu antibiootti vähensi hiirten nivelturvotusta ja esti B. burgdorferi – bakteerin kasvun kudoksista otetuista näytteissä. Hoidettu¬jen hiirten B. burgdorferi -spesifiset IgG-luokan vasta-aineet pysyivät kuitenkin koholla ja osasta kudosnäytteistä löytyi B. burgdorferi:n DNA:ta PCR-tutkimuksen avulla. Mikäli hiiret hoidettiin myöhäisessä vaiheessa (yli 18 viikkoa infektoinnista) tulokset olivat muuten samanlaiset, mutta keftriaksoni ei vaikuttanut nivelturvotukseen. Näin hiirissä oli aikaansaatu tilanne, joka on hyvin samankaltainen ihmisen hoitoresistentin borrelia-artriitin kanssa: oireet jatkuvat, mutta taudinaiheuttajaa ei saada esiin. Inflammaatiota vaimentavaa anti-TNF-alphaa on käytetty nivelreuman hoidossa menestyksekkäästi huonosti muuhun hoitoon reagoivilla potilailla ja siitä syystä sen ajateltiin voivan vaikuttaa suotuisasti myös B. burgdorferi -infektoitujen hiirten hoidonjälkeiseen nivelturvotukseen. Sillä ei kuitenkaan ollut vaikutusta nivelturvotukseen, mutta yllättäen hoidon jälkeen osa hiirten kudosnäytteistä osoittautui viljelypositiivisiksi. On siis ilmeistä, että hiirimallissamme osa B. burgdorferi spirokeetoista pystyy välttämään keftriaksonihoidon vaikutuksen joko hakeutumalla elimistössä kudokseen, jossa antibiootin pitoisuus ei nouse riittävän korkeaksi, tai ne kykenevät muuntautumaan metabolisesti inaktiiviin tilaan eikä mikrobilääke yhdessä immuunipuolustuksen kanssa onnistu tappamaan niitä. Jatkotutkimuksissa selvitimme B. burgdorferi - spirokeetan mahdollista piilopaikkaa tutkimalla antibioottihoidon jälkeen useita eri kudoksia PCR-menetelmällä. Tulosten perusteella spirokeetta näyttää suosivan nivelkudosta tai soluja, joita esiintyy nivelessä runsaasti. On kuitenkin edelleen epäselvää, missä muodossa B. burgdorferi –spirokeetat säilyvät kudoksessa antibioottihoidon jälkeen.

Некоторые наблюдения по поводу построения компьютерного морфологического анализатора марийского языка. [Problems in writing a two-level model of Mari morphology]

The article describes some concrete problems that were encountered when writing a two-level model of Mari morphology. Mari is an agglutinative Finno-Ugric language spoken in Russia by about 600 000 people. The work was begun in the 1980s on the basis of K. Koskenniemi’s Two-Level Morphology (1983), but in the latest stage R. Beesley’s and L. Karttunen’s Finite State Morphology (2003) was used. Many of the problems described in the article concern the inexplicitness of the rules in Mari grammars and the lack of information about the exact distribution of some suffixes, e.g. enclitics. The Mari grammars usually give complete paradigms for a few unproblematic verb stems, whereas the difficult or unclear forms of certain verbs are only superficially discussed. Another example of phenomena that are poorly described in grammars is the way suffixes with an initial sibilant combine to stems ending in a sibilant. The help of informants and searches from electronic corpora were used to overcome such difficulties in the development of the two-level model of Mari. The variation of the order of plural markers, case suffixes and possessive suffixes is a typical feature of Mari. The morphotactic rules constructed for Mari declensional forms tend to be recursive and their productivity must be limited by some technical device, such as filters. In the present model, certain plural markers were treated like nouns. The positional and functional versatility of the possessive suffixes can be regarded as the most challenging phenomenon in attempts to formalize the Mari morphology. Cyrillic orthography, which was used in the model, also caused problems. For instance, a Cyrillic letter may represent a sequence of two sounds, the first being part of the word stem while the other belongs to a suffix. In some cases, letters for voiced consonants are also generalized to represent voiceless consonants. Such orthographical conventions distance a morphological model based on orthography from the actual (morpho)phonological processes in the language.

Computational Evaluation of Print Unevenness According to Human Vision

Print quality and the printability of paper are very important attributes when modern printing applications are considered. In prints containing images, high print quality is a basic requirement. Tone unevenness and non uniform glossiness of printed products are the most disturbing factors influencing overall print quality. These defects are caused by non ideal interactions of paper, ink and printing devices in high speed printing processes. Since print quality is a perceptive characteristic, the measurement of unevenness according to human vision is a significant problem. In this thesis, the mottling phenomenon is studied. Mottling is a printing defect characterized by a spotty, non uniform appearance in solid printed areas. Print mottle is usually the result of uneven ink lay down or non uniform ink absorption across the paper surface, especially visible in mid tone imagery or areas of uniform color, such as solids and continuous tone screen builds. By using existing knowledge on visual perception and known methods to quantify print tone variation, a new method for print unevenness evaluation is introduced. The method is compared to previous results in the field and is supported by psychometric experiments. Pilot studies are made to estimate the effect of optical paper characteristics prior to printing, on the unevenness of the printed area after printing. Instrumental methods for print unevenness evaluation have been compared and the results of the comparison indicate that the proposed method produces better results in terms of visual evaluation correspondence. The method has been successfully implemented as ail industrial application and is proved to be a reliable substitute to visual expertise.