Motion artifact reduction and vessel enhancement for free-breathing navigator-gated coronary MRA using 3D k-space reordering.

Breathing-induced bulk motion of the myocardium during data acquisition may cause severe image artifacts in coronary magnetic resonance angiography (MRA). Current motion compensation strategies include breath-holding or free-breathing MR navigator gating and tracking techniques. Navigator-based techniques have been further refined by the applications of sophisticated 2D k-space reordering techniques. A further improvement in image quality and a reduction of relative scanning duration may be expected from a 3D k-space reordering scheme. Therefore, a 3D k-space reordered acquisition scheme including a 3D navigator gated and corrected segmented k-space gradient echo imaging sequence for coronary MRA was implemented. This new zonal motion-adapted acquisition and reordering technique (ZMART) was developed on the basis of a numerical simulation of the Bloch equations. The technique was implemented on a commercial 1.5T MR system, and first phantom and in vivo experiments were performed. Consistent with the results of the theoretical findings, the results obtained in the phantom studies demonstrate a significant reduction of motion artifacts when compared to conventional (non-k-space reordered) gating techniques. Preliminary in vivo findings also compare favorably with the phantom experiments and theoretical considerations. Magn Reson Med 45:645-652, 2001.

Particle transport method for convection-diffusion-reaction problems

Convective transport, both pure and combined with diffusion and reaction, can be observed in a wide range of physical and industrial applications, such as heat and mass transfer, crystal growth or biomechanics. The numerical approximation of this class of problemscan present substantial difficulties clue to regions of high gradients (steep fronts) of the solution, where generation of spurious oscillations or smearing should be precluded. This work is devoted to the development of an efficient numerical technique to deal with pure linear convection and convection-dominated problems in the frame-work of convection-diffusion-reaction systems. The particle transport method, developed in this study, is based on using rneshless numerical particles which carry out the solution along the characteristics defining the convective transport. The resolution of steep fronts of the solution is controlled by a special spacial adaptivity procedure. The serni-Lagrangian particle transport method uses an Eulerian fixed grid to represent the solution. In the case of convection-diffusion-reaction problems, the method is combined with diffusion and reaction solvers within an operator splitting approach. To transfer the solution from the particle set onto the grid, a fast monotone projection technique is designed. Our numerical results confirm that the method has a spacial accuracy of the second order and can be faster than typical grid-based methods of the same order; for pure linear convection problems the method demonstrates optimal linear complexity. The method works on structured and unstructured meshes, demonstrating a high-resolution property in the regions of steep fronts of the solution. Moreover, the particle transport method can be successfully used for the numerical simulation of the real-life problems in, for example, chemical engineering.

Modeling dynamic behavior in tilting pad gas journal bearings

Thisresearch deals with the dynamic modeling of gas lubricated tilting pad journal bearings provided with spring supported pads, including experimental verification of the computation. On the basis of a mathematical model of a film bearing, a computer program has been developed, which can be used for the simulation of a special type of tilting pad gas journal bearing supported by a rotary spring under different loading conditions time dependently (transient running conditions due to geometry variations in time externally imposed). On the basis of literature, different transformations have been used in the model to achieve simpler calculation. The numerical simulation is used to solve a non-stationary case of a gasfilm. The simulation results were compared with literature results in a stationary case (steady running conditions) and they were found to be equal. In addition to this, comparisons were made with a number of stationary and non-stationary bearing tests, which were performed at Lappeenranta University of Technology using bearings designed with the simulation program. A study was also made using numerical simulation and literature to establish the influence of the different bearing parameters on the stability of the bearing. Comparison work was done with literature on tilting pad gas bearings. This bearing type is rarely used. One literature reference has studied the same bearing type as that used in LUT. A new design of tilting pad gas bearing is introduced. It is based on a stainless steel body and electron beam welding of the bearing parts. It has good operation characteristics and is easier to tune and faster to manufacture than traditional constructions. It is also suitable for large serial production.

Application of the N-alkane molecular alloys to thermally protected containers for catering

A thermally controlled transport device was designed and tested. As hot food needs to be transported at temperatures between 60 and 70ºC in order to avoid contamination by microorganisms, the use of Molecular Alloy Phase Change Materials (MAPCM) can lead to improvements in this field of application. A heat transfer numerical simulation of the box used for transporting the food was conducted. Despite obvious simplifications, a good agreement between numerical simulation and experimental results was obtained. Furthermore, we compared our experimental results with those from other experiments related to the transport of hot food. Here, pizza is taken as the example, and it is shown that delivering time can be increased three-fold.

Evolutionary Games in Networked Populations: Models and Experiments

Cooperation and coordination are desirable behaviors that are fundamental for the harmonious development of society. People need to rely on cooperation with other individuals in many aspects of everyday life, such as teamwork and economic exchange in anonymous markets. However, cooperation may easily fall prey to exploitation by selfish individuals who only care about short- term gain. For cooperation to evolve, specific conditions and mechanisms are required, such as kinship, direct and indirect reciprocity through repeated interactions, or external interventions such as punishment. In this dissertation we investigate the effect of the network structure of the population on the evolution of cooperation and coordination. We consider several kinds of static and dynamical network topologies, such as Baraba´si-Albert, social network models and spatial networks. We perform numerical simulations and laboratory experiments using the Prisoner's Dilemma and co- ordination games in order to contrast human behavior with theoretical results. We show by numerical simulations that even a moderate amount of random noise on the Baraba´si-Albert scale-free network links causes a significant loss of cooperation, to the point that cooperation almost vanishes altogether in the Prisoner's Dilemma when the noise rate is high enough. Moreover, when we consider fixed social-like networks we find that current models of social networks may allow cooperation to emerge and to be robust at least as much as in scale-free networks. In the framework of spatial networks, we investigate whether cooperation can evolve and be stable when agents move randomly or performing Le´vy flights in a continuous space. We also consider discrete space adopting purposeful mobility and binary birth-death process to dis- cover emergent cooperative patterns. The fundamental result is that cooperation may be enhanced when this migration is opportunistic or even when agents follow very simple heuristics. In the experimental laboratory, we investigate the issue of social coordination between indi- viduals located on networks of contacts. In contrast to simulations, we find that human players dynamics do not converge to the efficient outcome more often in a social-like network than in a random network. In another experiment, we study the behavior of people who play a pure co- ordination game in a spatial environment in which they can move around and when changing convention is costly. We find that each convention forms homogeneous clusters and is adopted by approximately half of the individuals. When we provide them with global information, i.e., the number of subjects currently adopting one of the conventions, global consensus is reached in most, but not all, cases. Our results allow us to extract the heuristics used by the participants and to build a numerical simulation model that agrees very well with the experiments. Our findings have important implications for policymakers intending to promote specific, desired behaviors in a mobile population. Furthermore, we carry out an experiment with human subjects playing the Prisoner's Dilemma game in a diluted grid where people are able to move around. In contrast to previous results on purposeful rewiring in relational networks, we find no noticeable effect of mobility in space on the level of cooperation. Clusters of cooperators form momentarily but in a few rounds they dissolve as cooperators at the boundaries stop tolerating being cheated upon. Our results highlight the difficulties that mobile agents have to establish a cooperative environment in a spatial setting without a device such as reputation or the possibility of retaliation. i.e. punishment. Finally, we test experimentally the evolution of cooperation in social networks taking into ac- count a setting where we allow people to make or break links at their will. In this work we give particular attention to whether information on an individual's actions is freely available to poten- tial partners or not. Studying the role of information is relevant as information on other people's actions is often not available for free: a recruiting firm may need to call a job candidate's refer- ences, a bank may need to find out about the credit history of a new client, etc. We find that people cooperate almost fully when information on their actions is freely available to their potential part- ners. Cooperation is less likely, however, if people have to pay about half of what they gain from cooperating with a cooperator. Cooperation declines even further if people have to pay a cost that is almost equivalent to the gain from cooperating with a cooperator. Thus, costly information on potential neighbors' actions can undermine the incentive to cooperate in dynamical networks.

Puumassapesurin suotautumisolojen kartoitus ja simulointi

Työn tavoitteena oli tutkia sellumassan suotautumiseen liittyviä tekijöitä Ahlstrom Machineryn Drum Displacer™ -puumassapesurissa (DD-pesuri). Teoriaosassa tarkasteltiin aluksi suotautumisen teoriaa kuitu-vesi-suspensiossa, minkä jälkeen esiteltiin suotautumisnopeuteen vaikuttavia fysikaalisia ja kemiallisia vaikutusmekanismeja. Seuraavaksi kuvattiin massan pesun yleisiä perusteita sekä teoriaa puumassapesureissa. Lopuksi tarkasteltiin pesurien kytkeytymistä muuhun kuitulinjaan sekä prosessista johtuvia pesun toiminnan ulkoisia häiriötekijöitä. Kokeellisen osan aluksi tarkasteltiin paine- ja lämpötilamittauksien avulla massapesurissa vallitsevia prosessioloja. Mittaustulosten perusteella pumppausolot pesurin suodoslinjoissa ovat vaikeahkot ja häiriötilanteita voi esiintyä, mutta käytäntö on osoittanut tästä olevan vain harvoin haittaa prosessin toiminnalle. Pesureissa toteutuneet syrjäytysnopeudet laskettiin ja niitä verrattiin syrjäytystestien antamiin tuloksiin. Kuitulinjasta riippuen testin vastaavuus tehdasprosessiin vaihteli suuresti. Syrjäytystesteillä kokeiltiin myös tehdasprosesseissa usein esiintyvien muuttujien vaikutusta sellukakun syrjäytettävyyteen. Kakun paksuus ja syrjäytyslämpötila vaikuttivat syrjäytysnopeuteen Darcyn lain mukaisesti. Alipaineen massakakun alapuolella havaittiin huonontavan syrjäytysnopeutta verrattuna tilanteeseen, jossa kakun alla vallitsi ilmanpaine. Tämä havainto on selvästi ristiriidassa suotautumisen teorian kanssa. Massakakun muodostumis-pH osoittautui ratkaisevaksi lopulliselle syrjäytysnopeudelle, sillä alkalisissa oloissa muodostetun kuitukakun syrjäytysnopeus ei enää parantunut happamalla syrjäytysnesteellä. Happamissa oloissa muodostetun kakun syrjäytysnopeus oli alkalisista parempi, mutta se alkoi hitaasti alentua, kun syrjäytysneste vaihtui alkaliseen. Massan laimentaminen ennen syrjäytystä alkalisella tehdassuodoksella puhtaan veden sijasta alensi ligniinipitoisella massalla lopullista syrjäytysnopeutta. Shirato-Tillerin mallilla ja Jönssonin staattisella mallilla simuloitiin numeerisesti syrjäytystestiä kahdessa eri pH:ssa, ja simulointituloksia verrattiin vastaavissa oloissa tehtyihin syrjäytystesteihin. Shirato-Tillerin mallin antamien syrjäytysnopeuksien havaittiin olevan lähellä syrjäytystestien nopeuksia, kun Jönssonin mallin antamat tulokset jäivät huomattavasti testituloksia alemmiksi. Herkkyystarkastelussa havaittiin mallien olevan varsin herkkiä parametrien virheille. Hajonta vaadittavien kuituparametrien määrityksissä ja menetelmien työläys rajoittavat numeerisen simuloinnin käytettävyyttä, sillä kuituparametrien määrityksen vaatima työmäärä on ainakin toistaiseksi syrjäytystestiä suurempi. Lopuksi todettiin, että oikeiden syrjäytysolosuhteiden käyttö on ensiarvoisen tärkeää oikeiden tulosten saamiseksi sekä kokeellisessa että numeerisessa simuloinnissa. Nykyinen syrjäytystestilaitteisto on pienin muutoksin käyttökelpoinen, kun massan testaus prosessioloissa tulee rutiininomaiseksi.

Synthesis of highly focused fields with circular polarization at any transverse plane

We develop a method for generating focused vector beams with circular polarization at any transverse plane. Based on the Richards-Wolf vector model, we derive analytical expressions to describe the propagation of these set of beams near the focal area. Since the polarization and the amplitude of the input beam are not uniform, an interferometric system capable of generating spatially-variant polarized beams has to be used. In particular, this wavefront is manipulated by means of spatial light modulators displaying computer generated holograms and subsequently focused using a high numerical aperture objective lens. Experimental results using a NA=0.85 system are provided: irradiance and Stokes images of the focused field at different planes near the focal plane are presented and compared with those obtained by numerical simulation.

Propagação da poluição atmosférica por flúor nas águas subterrâneas e solos de regiões próximas às indústrias de fertilizantes (Rio Grande, RS)

Fluoride concentration was determined in rainwater, ground water and soil in the zone of fertilizer industry in the city of Rio Grande. In contaminated rainwater fluorine concentration was registered up to the value of 4,4 mg.L-1. Fluorine concentration in the shallow ground water in general reflects its distribution in the atmosphere, but cannot be used as marker of atmosphere contamination in the urban area due to dissolving influence of residential effluents. The 0,01% HCl extracts from the set of surface soil samples demonstrates fluoride distribution in the zone of influence of industrial emissions, which coincides to the numerical simulation of fluorides dispersion in the air.

Chaotic dynamics of electric-field domains in periodically driven superlattices

Self-sustained time-dependent current oscillations under dc voltage bias have been observed in recent experiments on n-doped semiconductor superlattices with sequential resonant tunneling. The current oscillations are caused by the motion and recycling of the domain wall separating low- and high-electric-field regions of the superlattice, as the analysis of a discrete drift model shows and experimental evidence supports. Numerical simulation shows that different nonlinear dynamical regimes of the domain wall appear when an external microwave signal is superimposed on the dc bias and its driving frequency and driving amplitude vary. On the frequency-amplitude parameter plane, there are regions of entrainment and quasiperiodicity forming Arnold tongues. Chaos is demonstrated to appear at the boundaries of the tongues and in the regions where they overlap. Coexistence of up to four electric-field domains randomly nucleated in space is detected under ac+dc driving.

Application of a non isotropic turbulence model to stable atmospheric flows and dispersion over 3D topography

A non isotropic turbulence model is extended and applied to three dimensional stably stratified flows and dispersion calculations. The model is derived from the algebraic stress model (including wall proximity effects), but it retains the simplicity of the "eddy viscosity" concept of first order models. The "modified k-epsilon" is implemented in a three dimensional numerical code. Once the flow is resolved, the predicted velocity and turbulence fields are interpolated into a second grid and used to solve the concentration equation. To evaluate the model, various steady state numerical solutions are compared with small scale dispersion experiments which were conducted at the wind tunnel of Mitsubishi Heavy Industries, in Japan. Stably stratified flows and plume dispersion over three distinct idealized complex topographies (flat and hilly terrain) are studied. Vertical profiles of velocity and pollutant concentration are shown and discussed. Also, comparisons are made against the results obtained with the standard k-epsilon model.

Thermalhydraulic behavior of electrically heated rod during a critical heat flux transient

In nuclear reactors, the occurrence of critical heat flux leads to fuel rod overheating with clad fusion and radioactive products leakage. To predict the effects of such phenomenon, experiments are performed using electrically heated rods to simulate operational and accidental conditions of nuclear fuel rods. In the present work, it is performed a theoretical analysis of the drying and rewetting front propagation during a critical heat flux experiment, starting with the application of an electrical power step from steady state condition. After the occurrence of critical heat flux, the drying front propagation is predicted. After a few seconds, a power cut is considered and the rewetting front behavior is analytically observed. Studies performed with various values of coolant mass flow rate show that this variable has more influence on the drying front velocity than on the rewetting one.

Surface Tension Implementation for Gensmac2d

In the present work we describe a method which allows the incorporation of surface tension into the GENSMAC2D code. This is achieved on two scales. First on the scale of a cell, the surface tension effects are incorporated into the free surface boundary conditions through the computation of the capillary pressure. The required curvature is estimated by fitting a least square circle to the free surface using the tracking particles in the cell and in its close neighbors. On a sub-cell scale, short wavelength perturbations are filtered out using a local 4-point stencil which is mass conservative. An efficient implementation is obtained through a dual representation of the cell data, using both a matrix representation, for ease at identifying neighbouring cells, and also a tree data structure, which permits the representation of specific groups of cells with additional information pertaining to that group. The resulting code is shown to be robust, and to produce accurate results when compared with exact solutions of selected fluid dynamic problems involving surface tension.

Dynamic Analysis Model of Spherical Roller Bearings with Defects

Rolling element bearings are essential components of rotating machinery. The spherical roller bearing (SRB) is one variant seeing increasing use, because it is self-aligning and can support high loads. It is becoming increasingly important to understand how the SRB responds dynamically under a variety of conditions. This doctoral dissertation introduces a computationally efficient, three-degree-of-freedom, SRB model that was developed to predict the transient dynamic behaviors of a rotor-SRB system. In the model, bearing forces and deflections were calculated as a function of contact deformation and bearing geometry parameters according to nonlinear Hertzian contact theory. The results reveal how some of the more important parameters; such as diametral clearance, the number of rollers, and osculation number; influence ultimate bearing performance. Distributed defects, such as the waviness of the inner and outer ring, and localized defects, such as inner and outer ring defects, are taken into consideration in the proposed model. Simulation results were verified with results obtained by applying the formula for the spherical roller bearing radial deflection and the commercial bearing analysis software. Following model verification, a numerical simulation was carried out successfully for a full rotor-bearing system to demonstrate the application of this newly developed SRB model in a typical real world analysis. Accuracy of the model was verified by comparing measured to predicted behaviors for equivalent systems.

SEPARATION OF SATURED AND UNSATURATED FATTY ACIDS FROM PALM FATTY ACIDS DISTILLATES IN CONTINUOUS MULTISTAGE COUNTERCURRENT COLUMNS WITH SUPERCRITICAL CARBON DIOXIDE AS SOLVENT: A PROCESS DESIGN METHODOLOGY

In this work the separation of multicomponent mixtures in counter-current columns with supercritical carbon dioxide has been investigated using a process design methodology. First the separation task must be defined, then phase equilibria experiments are carried out, and the data obtained are correlated with thermodynamic models or empirical functions. Mutual solubilities, Ki-values, and separation factors aij are determined. Based on this data possible operating conditions for further extraction experiments can be determined. Separation analysis using graphical methods are performed to optimize the process parameters. Hydrodynamic experiments are carried out to determine the flow capacity diagram. Extraction experiments in laboratory scale are planned and carried out in order to determine HETP values, to validate the simulation results, and to provide new materials for additional phase equilibria experiments, needed to determine the dependence of separation factors on concetration. Numerical simulation of the separation process and auxiliary systems is carried out to optimize the number of stages, solvent-to-feed ratio, product purity, yield, and energy consumption. Scale-up and cost analysis close the process design. The separation of palmitic acid and (oleic+linoleic) acids from PFAD-Palm Fatty Acids Distillates was used as a case study.

Effects of Fabrication Processes on the Behavior of S960 Ultra High Strength Steel (UHSS)

Strenx® 960 MC is a direct quenched type of Ultra High Strength Steel (UHSS) with low carbon content. Although this material combines high strength and good ductility, it is highly sensitive towards fabrication processes. The presence of stress concentration due to structural discontinuity or notch will highlight the role of these fabrication effects on the deformation capacity of the material. Due to this, a series of tensile tests are done on both pure base material (BM) and when it has been subjected to Heat Input (HI) and Cold Forming (CF). The surface of the material was dressed by laser beam with a certain speed to study the effect of HI while the CF is done by bending the specimen to a certain angle prior to tensile test. The generated results illustrate the impact of these processes on the deformation capacity of the material, specially, when the material has HI experience due to welding or similar processes. In order to compare the results with those of numerical simulation, LS-DYNA explicit commercial package has been utilized. The generated results show an acceptable agreement between experimental and numerical simulation outcomes.