Made available in DSpace on 2016-12-23T14:04:30Z (GMT). No. of bitstreams: 1 Wesley Jose Nunes Pinto.pdf: 2090367 bytes, checksum: e676fde8423a3a2cfeac61da24020ea8 (MD5) Previous issue date: 2013-08-19 / Neste estudo foi aplicado do método numérico, sem malhas, baseado em partículas, denominado SPH (Smoothed Particles Hydrodynamics). E um código numérico na linguagem computacional FORTRAN foi utilizado para solucionar as equações de Navier-Stokes. O clássico problema da literatura da dinâmica dos fluidos Computacional, denotado como problema da cavidade quadrada bidimensional (Shear-Driven Cavity Flow) , foi estudado com a intenção de verificar o comportamento do código numérico em relação a resultados específicos já existentes do assunto. O citado problema físico das cavidades abertas é amplamente empregado como benchmark, visando a validação do método numérico utilizado no trabalho desenvolvido na pesquisa. O trabalho de análise e validação do código numérico foi dividido em três seções: a primeira lista as localizações dos centros dos vórtices principais gerados pelo escoamento na aresta superior das cavidades; a segunda plota os perfis das componentes das velocidades centrais das cavidades; e a terceira: lista os desvios absolutos dos perfis das velocidades centrais do presente trabalho, comparados com dados de outros estudos. Constata-se que o método SPH apresentou boa acurácia nas simulações realizadas, obtendo boa concordância entre os resultados das simulações dinâmicas com os dados de referências, validando-se o modelo numérico proposto, tendo melhores resultados para baixos números de Reynolds / In this study, it was applied the numerical method, grid-free, based on particles named SPH (Smoothed Particles Hydrodynamics). Also, a numerical code in the computer language FORTRAN was used to solve the Navier-Stokes Equations. This classic problem of the literature related to Computational Fluid Dynamics indicated as Shear-Driven Cavity Flow was studied to check the behavior of the numerical code regarding specific existing results. Such problem is highly used as Benchmark, aiming the validation of the numerical method used to develop the research. The analysis and validation of the numerical code was divided into three sections: the first one lists the location of the centre of the main vortex generated by the flow of the upper edge of the cavities; the second one plots the profiles of the components of the central speed of the cavities; the third one lists the absolute deviation of the profiles of the central speed of this study compared with other cases data. It is established that the SPH Method presented accuracy in the performed simulations, in a consonance between the results of the dynamic simulations and the reference data, thus the proposed numerical model was validated with better results for low Reynolds numbers

<p>A new model for simulating non-Brownian flexible fibres suspended in a Newtonian fluid has been developed. Special attention has been given to include realistic flow conditions found in the industrial papermaking process in the key features of the model; it is the intention of the author to employ the model in simulations of the forming section of the paper machine in future studies.</p><p>The model considers inert fibres of various shapes and finite stiffness, interacting with each other through normal, frictional and lubrication forces, and with the surrounding fluid medium through hydrodynamic forces. Fibre-fluid interactions in the non-creeping flow regime are taken into account, and the two-way coupling between the solids and the fluid phase is included by enforcing momentum conservation between phases. The incompressible three-dimensional Navier-Stokes equations are employed to model the motion of the fluid medium.</p><p>The validity of the model has been tested by comparing simulation results with experimental data from the literature. It was demonstrated that the model predicts the motion of isolated fibres in shear flow over a wide range of fibre flexibilities. It was also shown that the model predicts details of the orientation distribution of multiple straight, rigid fibres in a sheared suspension. Model predictions of the viscosity and first normal stress difference were in good agreement with experimental data found in the literature. Since the model is based solely on first-principles physics, quantitative predictions could be made without any parameter fitting.</p> / <p>En ny modell för simulering av rörelserna hos icke-Brownska böjliga fibrer dispergerade i en Newtonsk vätska har utvecklats. Eftersom det är författarens avsikt att modellen skall kunna tillämpas vid simulering av arkformning under de förhållanden som råder i en modern pappersmaskin, har särskilt omsorg givits till att inkludera motsvarande flödesvillkor i modellens giltighetsområde.</p><p>Modellen hanterar fibrer av varierande form, massa och styvhet, som växelverkar sinsemellan via normal-, friktions- och smörjkrafter. Deras växelverkan med den omgivande vätskan sker via hydrodynamiska krafter vid finita Reynolds-tal. Den så kallade tvåvägskopplingen mellan fibrerna och vätskefasen har tagits i beaktande genom att kräva att rörelsemängden bevaras vid interaktionen mellan faserna. Vidare har Navier-Stokes ekvationer för inkompressibla vätskor använts för att beskriva mediets rörelser.</p><p>Modellens giltighet har undersökts genom att jämföra resultat från simuleringar med experimentella data från litteraturen. Det har påvisats att modellen förutsäger rörelsen hos ensamma fibrer i ett skjuvflöde, för vitt skilda fiberflexibiliteter. Det visades också att modellen förutsäger detaljer hos fiberorienteringsdistributionen i suspensioner utsatta för skjuvflöde. Det kunde också konstateras att modellens förutsägelser av fibersuspensioners viskositet och första normalspänningsdifferens under skjuvning väl överensstämde med experimentella data i litteraturen. Kvantitativa förutsägelser har kunnat göras utan någon parameteranpassning, då modellen bygger uteslutande på väletablerade fysikaliska samband inom klassisk mekanik och strömningslära.</p>

Gas, liquid and powder flow in the lower part of a blast furnace is complex phenomena. In order to understand the aerodynamics of the blast furnace properly, these phenomena must be included in their advanced form. Previous studies have shown that the conditions of blast furnace resemble the cold model experiments which have been done in decreasing gas velocities. Also, the recent studies have shown that liquid flow in a blast furnace can be represented more realistically considering it discrete in nature. In the current study, both the phenomena have been considered along with the injection of powder through a nozzle while studying the fluid flow behaviour in a packed bed. The situation resembles the lower part of an ironmaking blast furnace. In this study, gas flow has been modelled using k-ε turbulent model and has been coupled with previously developed stress model to calculate the raceway size. Coal powder is treated as continuum and has been modelled in the similar way as gas flow. After this gas and powder flow model were coupled with previously developed discrete liquid flow model. Liquid flow model has been considered for structured bed only. The governing equations for gas phase were discretized. Finite Volume method was used for the solution. Co-located grid is used for the simulation. Blending of upwind difference scheme and central difference scheme (deferred correction approach) is used to achieve the stability of upwind scheme and accuracy of central difference scheme. Similar treatment was employed for powder phase. For the solution of volume fraction of powder, powder phase continuity equation was used along with pseudo time step scheme. Results obtained from gas and powder models have been validated against published experimental data. Similarly, gas-liquid flow results have been validated against published experimental data on gas-powder flow. Results obtained by gas-powder-liquid model could not be validated against any experimental or theoretical data as they are not available in the literature. The effect of various parameters on the fluid flow (gas/liquid/powder) behaviour have been studied like the effect of increasing and decreasing gas velocities, flow rates of liquid, gas and powder, size of powder and packing etc. It is found that the above mentioned phenomena have significant effect on the fluid flow behaviour in a packed bed.

The fluid flow and tracer transport in a single rock fracture during shear processes has been an important issue in rock mechanics and is investigated in this thesis using Finite Element Method (FEM) and streamline particle tracking method, considering evolutions of aperture and transmissivity with shear displacement histories under different normal stresses, based on laboratory tests. The distributions of fracture aperture and its evolution during shear were calculated from the initial aperture fields, based on the laser-scanned surface roughness features of replicas of rock fracture specimens, and shear dilations measured during the coupled shear-flow-tracer tests in laboratory performed using a newly developed testing apparatus in Nagasaki University, Nagasaki, Japan. Three rock fractures of granite with different roughness characteristics were used as parent samples from which nine plaster replicas were made and coupled shear-flow tests was performed under three normal loading conditions (two levels of constant normal loading (CNL) and one constant normal stiffness (CNS) conditions). In order to visualize the tracer transport, transparent acrylic upper parts and plaster lower parts of the fracture specimens were manufactured from an artificially created tensile fracture of sandstone and the coupled shear-flow tests with fluid visualization was performed using a dye tracer injected from upstream and a CCD camera to record the dye movement. A special algorithm for treating the contact areas as zero-aperture elements was used to produce more accurate flow field simulations by using FEM, which is important for continued simulations of particle transport, but was often not properly treated in literature. The simulation results agreed well with the flow rate data obtained from the laboratory tests, showing that complex histories of fracture aperture and tortuous flow channels with changing normal stresses and increasing shear displacements, which were also captured by the coupled shear-flow tests of fracture specimens with visualization of the fluid flow. From the obtained flow velocity fields, the particle transport was predicted by the streamline particle tracking method with calculated flow velocity fields (vectors) from the flow simulations, obtaining results such as flow velocity profiles, total flow rates, particle travel time, breakthrough curves and the Péclet number, Pe, respectively. The fluid flow in the vertical 2-D cross-sections of a rock fracture was also simulated by solving both Navier-Stokes (NS) and Reynolds equations, and the particle transport was predicted by streamline particle tracking method. The results obtained using NS and Reynolds equations were compared to illustrate the degree of the validity of the Reynolds equation for general applications in practice since the later is mush more computationally efficient for large scale problems. The flow simulation results show that the total flow rate and the flow velocity predicted by NS equations are quite different from that as predicted by the Reynolds equation. The results show that a roughly 5-10 % overestimation on the flow rate is produced when the Reynolds equation is used, and the ideal parabolic velocity profiles defined by the local cubic law, when Reynolds equation is used, is no longer valid, especially when the roughness feature of the fracture surfaces changes with shear. These deviations of flow rate and flow velocity profiles across the fracture aperture have a significant impact on the particle transport behavior and the associated properties, such as the travel time and Péclet number. The deviations increase with increasing flow velocity and become more significant when fracture aperture geometry changes with shear. The scientific findings from these studies provided new insights to the physical behavior of fluid flow and mass transport in rock fractures which is the scientific basis for many rock mechanics problems at the fundamental level, and with special importance to rock engineering problems such as geothermal energy extraction (where flow rate in fractures dominates the productivity of a geothermal energy reservoir) and nuclear waste repositories (where radioactive nuclides transport through fractures dominates the final safety evaluations) in fractured rocks. / Vätskeflödet och spårämnestransporten i en enskild bergsspricka under skjuvningsprocesser har varit ett viktigt ämne inom bergmekanik. I denna avhandling undersöks ämnet med hjälp av finita element metoden (FEM) och en strömlinjebaserad partikelspårningsmetod. Hänsyn tas till utveckling av öppningar och transmissivitet med skjuvningens förflyttningshistoria under olika normala belastningar baserat på laboratorietester. Fördelningen av spricköppningar och deras utveckling under skjuvning beräknades från de initiala öppningsfälten baserat på det laserscannade provets ytas grovhetskännetecken sam tskjuvningsöppningar uppmätta under de kopplade skjuvning-flöde-spårämneslaboratorietesterna som utförts med nyutvecklad testapparatur i Nagasaki Universitet i Nagasaki, Japan. Tre bergssprickor i granit med olika grovhetskarakteristika användes som utgångsprover från vilka nio gipskopior gjordes. Kopplade skjuvning-flödes tester utfördes sedan under tre normala belastningstillstånd (två nivåer med konstant normal last (KNL) och en konstant normal styvhetstillstånd (KNS). För att visualisera spårämnestransporten tillverkades en transparent övre del av sprickproverna av akryl och en nedre del av gipsbaserat på en kostgjord spänningsspricka i sandsten och de kopplade skjuvning-flödes testerna med vätskevisualisering utfördes med färgspårämne injekterat uppströms och en CCD kamera monterad ovanför för att registrera färgens rörelse. En särskild algoritm användes för att behandla kontaktytorna som nollöppningsämnen användes för att åstadkomma mer exakta flödesfältssimuleringar med FEM. Detta är viktigt för kontinuerliga simuleringar av partikelflöden men uppmärksammas oftast inte tillräckligt i litteraturen. Simuleringsresultaten överensstämde väl med de flödesnivådata som erhölls från laboratorietesterna vilket visade att komplexa historier av spricköppningar och invecklade flöden överensstämde med ändrade normala belastningar och ökande skjuvningsförflyttningar, vilket även fångades av de kopplade skjuvning-flödestesterna av sprickproverna genom visualisering av vätskeflödet. Från de erhållna flödesfälten förutsågs partikeltransporten genom en strömlinjebaserad partikelspårningsmetod med kalkylerade flödeshastighetsfält (vektorer) från flödessimuleringarna genom vilka resultat som flödeshastighetsprofiler, totala flödesnivåer,partikeltransporttid, genombrottskurvor samt Pécletnumret, Pe, erhölls. Vätskeflödet i det vertikala tvådimensionella tvärsnittet av en bergsspricka simulerades även genom att både Navier-Stokes (NS) och Reynoldsekvationerna löstes och partikeltransporten förutsågs genom den strömlinjebaserade partikelspårningsmetoden. Resultaten som erhöllsmed NS och Reynoldsekvationerna jämfördes för att illustrera graden av tillförlitlighet för Reynoldsekvationen för allmänna tillämpningar i praktiken då den senare är betydligt mer beräkningseffektiv för storskaliga problem. Resultaten från flödessimuleringarna visar att den totala flödesnivån och den totala flödeshastigheten förutsedda med NS ekvationer är helt annorlunda motsvarande värden som förutsågs med Reynoldsekvationen. Resultaten visar att en ca 5-10 % för hög uppskattning av flödesnivån erhålls då Reynoldsekvationen används och de ideala parabola hastighetsprofilerna, som definieras av den lokala kubiklagen när Reynoldsekvationen används, inte längre är giltiga särskilt när sprickytornas grovhetskarakteristika ändras med skjuvning. De här avvikelserna i flödesnivå och flödeshastighetsprofiler längs med spricköppningen har en betydande påverkan på partikeltransportuppträdande och de tillhörande egenskaperna såsom rörelsetid och Pécletnummer. Avvikelserna ökar med ökande flödeshastighet och blir mer signifikanta när spricköppningarnas geometri ändras med skjuvning. Forskningsresultaten från dessa studier gav nya insikter i de fysiska uppträdandet av vätskeflöde och masstransporter i bergssprickor vilket är den vetenskapliga basen för många bergmekanikproblem på grundläggande nivå och som har särskild vikt för bergstekniksproblem såsom geotermisk energiutvinning (där flödesnivå i sprickor dominerar produktiviteten för en geotermisk energikälla) och kärnavfallsförvaringsplatser (där transporten av radioaktiva nuklider genom sprickor dominerar den slutgiltigasäkerhetsutvärderingen) i sprickigt berg. / QC 20100803

lthough extensive research work has been carried out on the drag reduction behavior of polymers and surfactants alone, little progress has been made on the synergistic effects of combined polymers and surfactants. A number of studies have demonstrated that certain types of polymers and surfactants interact with each other to form surfactant-polymer complexes. The formation of such complexes can cause changes in the solution properties and may result in better drag reduction characteristics as compared with pure additives. A series of drag-reducing surfactants and polymers were screened for the synergistic studies. The following two widely used polymeric drag reducing agents (DRA) were chosen: a copolymer of acrylamide and sodium acrylate (referred to as PAM) and polyethylene oxide (PEO). Among the different types of surfactants screened, a cationic surfactant octadecyltrimethylammonium chloride (OTAC) and an anionic surfactant Sodium dodecyl sulfate (SDS) were selected for the synergistic study. In the case of the cationic surfactant OTAC, sodium salicylate (NaSal) was used as a counterion. No counterion was used with anionic surfactant SDS. The physical properties such as viscosity, surface tension and electrical conductivity were measured in order to detect any interaction between the polymer and the surfactant. The drag reduction (DR) ability of both pure and mixed additives was investigated in a pipeline flow loop. The effects of different parameters such as additive concentration, type of water (deionized (DI) or tap), temperature, tube diameter, and mechanical degradation were investigated. The addition of OTAC to PAM solution has a significant effect on the properties of the system. The critical micelle concentration (CMC) of the mixed surfactant-polymer system is found to be different from that of the surfactant alone. The anionic PAM chains collapse upon the addition of cationic OTAC and a substantial decrease in the viscosity occurs. The pipeline flow behaviour of PAM/OTAC mixtures is found to be consistent with the bench scale results. The drag reduction ability of PAM is reduced upon the addition of OTAC. At low concentrations of PAM, the effect of OTAC on the drag reduction behavior is more pronounced. The drag reduction behavior of polymer solutions is strongly influenced by the nature of water (de-ionized or tap). The addition of OTAC to PEO solution exhibited a week interaction based on the viscosity and surface tension measurements. However, the pipeline results showed a considerable synergistic effect, that is, the mixed system gave a significantly higher drag reduction (lower friction factors) as compared with the pure additives (pure polymer or pure surfactant). The synergistic effect in the mixed system was stronger at low polymer concentrations and high surfactant concentrations. Also the resistance against mechanical degradation of the additive was improved upon the addition of OTAC to PEO. The mixed PEO/SDS system exhibited a strong interaction between the polymers (PEO) and the surfactant (SDS), Using electrical conductivity and surface tension measurements, the critical aggregation concentration (CAC) and the polymer saturation point (PSP) were determined. As the PEO concentration is increased, the CAC decreases and the PSP increase. The addition of SDS to the PEO solution exhibits a remarkable increase in the relative viscosity compared to the pure PEO solution. This increase is attributed to the changes in the hydrodynamic radius of the polymer coil. The pipeline flow exhibited a considerable increase in DR for the mixed system as compared to the pure PEO solution. The addition of surfactant always improves the extent of DR up to the PSP. Also the mixed PEO/ SDS system shows better resistance against shear degradation of the additive.

Well logs and formation testers are routinely used for detection and quantification of hydrocarbon reserves. Overbalanced drilling causes invasion of mud filtrate into permeable rocks, hence radial displacement of in-situ saturating fluids away from the wellbore. The spatial distribution of fluids in the near-wellbore region remains affected by a multitude of petrophysical and fluid factors originating from the process of mud-filtrate invasion. Consequently, depending on the type of drilling mud (e.g. water- and oil-base muds) and the influence of mud filtrate, well logs and formation-tester measurements are sensitive to a combination of in-situ (original) fluids and mud filtrate in addition to petrophysical properties of the invaded formations. This behavior can often impair the reliable assessment of hydrocarbon saturation and formation storage/mobility. The effect of mud-filtrate invasion on well logs and formation-tester measurements acquired in vertical wells has been extensively documented in the past. Much work is still needed to understand and quantify the influence of mud-filtrate invasion on well logs acquired in horizontal and deviated wells, where the spatial distribution of fluids in the near-wellbore region is not axial-symmetric in general, and can be appreciably affected by gravity segregation, permeability anisotropy, capillary pressure, and flow barriers. This dissertation develops a general algorithm to simulate the process of mud-filtrate invasion in vertical and deviated wells for drilling conditions that involve water- and oil-base mud. The algorithm is formulated in cylindrical coordinates to take advantage of the geometrical embedding imposed by the wellbore in the spatial distribution of fluids within invaded formations. In addition, the algorithm reproduces the formation of mudcake due to invasion in permeable formations and allows the simulation of pressure and fractional flow-rate measurements acquired with dual-packer and point-probe formation testers after the onset of invasion. An equation-of-state (EOS) formulation is invoked to simulate invasion with both water- and oil-base muds into rock formations saturated with water, oil, gas, or stable combinations of the three fluids. The algorithm also allows the simulation of physical dispersion, fluid miscibility, and wettability alteration. Discretized fluid flow equations are solved with an implicit pressure and explicit concentration (IMPEC) scheme. Thermodynamic equilibrium and mass balance, together with volume constraint equations govern the time-space evolution of molar and fluid-phase concentrations. Calculations of pressure-volume-temperature (PVT) properties of the hydrocarbon phase are performed with Peng-Robinson's equation of state. A full-tensor permeability formulation is implemented with mass balance equations to accurately model fluid flow behavior in horizontal and deviated wells. The simulator is rigorously and successfully verified with both analytical solutions and commercial simulators. Numerical simulations performed over a wide range of fluid and petrophysical conditions confirm the strong influence that well deviation angle can have on the spatial distribution of fluid saturation resulting from invasion, especially in the vicinity of flow barriers. Analysis on the effect of physical dispersion on the radial distribution of salt concentration shows that electrical resistivity logs could be greatly affected by salt dispersivity when the invading fluid has lower salinity than in-situ water. The effect of emulsifiers and oil-wetting agents present in oil-base mud was studied to quantify wettability alteration and changes in residual water saturation. It was found that wettability alteration releases a fraction of otherwise irreducible water during invasion and this causes electrical resistivity logs to exhibit an abnormal trend from shallow- to deep-sensing apparent resistivity. Simulation of formation-tester measurements acquired in deviated wells indicates that (i) invasion increases the pressure drop during both drawdown and buildup regimes, (ii) bed-boundary effects increase as the wellbore deviation angle increases, and (iii) a probe facing upward around the perimeter of the wellbore achieves the fastest fluid clean-up when the density of invading fluid is larger than that of in-situ fluid. / text

Στην παρούσα διατριβή, ερευνάται αν οι διακυμάνσεις της τριχοειδούς πίεσης είναι δυνατό να παράσχουν ποσοτική πληροφορία που αφορά τη χωρική κατανομή της διαβροχής σε πορώδη μέσα κλασματικής διαβροχής. Για το σκοπό αυτό πραγματοποιήθηκαν πειράματα εκτόπισης της ελαϊκής φάσης από υδατική υπό σταθερή παροχή σε πορώδη μέσα κλασματικής διαβροχής. Τα πορώδη μέσα κλασματικής διαβροχής κατασκευάστηκαν με τυχαία ανάμιξη γυάλινων και πλαστικών (PTFE) σφαιρών, απαράλλαχτης διαμέτρου, έτσι ώστε η γεωμετρία των πόρων να διατηρηθεί αναλλοίωτη. Η χωρική κατανομή των διαφορετικών ειδών σφαιρών διέφερε ακόμη και όταν τα σχετικά τους κλάσματα ήταν ίδια. Κατά τη διάρκεια των πειραμάτων λαμβάνονταν μετρήσεις της μεταβατικής απόκρισης της τριχοειδούς πίεσης και στιγμιότυπα της εκτόπισης ανά τακτά χρονικά διαστήματα. Αρχικά πραγματοποιήθηκαν πειράματα εκτόπισης σε πορώδη μέσα που ήταν ενδιάμεσα υδατοδιαβρεκτά και ισχυρά ελαιοδιαβρεκτά. Τα ληφθέντα σήματα τριχοειδούς πίεσης, αναλύθηκαν με τη μέθοδο των κυματιδίων (wavelet) και συσχετίσθηκαν με τα γεγονότα ροής που ελάμβαναν χώρα. Η διαβροχή των περιοχών των πορωδών μέσων που εισδύονταν από υδατική φάση περιγράφηκε μέσω ενός εμπειρικού δείκτη διαβροχής ο οποίος συσχετίστηκε με τα φάσματα λεπτομέρειας που προέκυψαν από την ανάλυση των σημάτων. Στην περίπτωση ομοιόμορφης και ενδιάμεσης υδατοδιαβροχής ευνοούνται γεγονότα ροής που συνοδεύονται από σημαντικές διακυμάνσεις της τριχοειδούς πίεσης ώστε να καθίσταται δύσκολη η παραπάνω συσχέτιση. Για το λόγο αυτό πραγματοποιήθηκαν πειράματα εκτόπισης με χρήση ενός δεύτερου ζεύγους ρευστών. Τα πορώδη μέσα κλασματικής διαβροχής ως προς το δεύτερο ζεύγος ρευστών ήταν υδατοδιαβρεκτά και ισχυρά ελαιοδιαβρεκτά. Στην περίπτωση του δεύτερου ζεύγους ρευστών, τα σήματα τριχοειδούς πίεσης χαρακτηρίζονταν από αργές και μικρές σε ύψος διακυμάνσεις. Αυτή η διαφορετικότητα επέβαλλε τη βελτίωση της μεθόδου ανάλυσης των σημάτων ώστε να μπορεί να εφαρμοστεί σε όλα τα σήματα που ελήφθησαν στα πειράματα εκτόπισης της παρούσας εργασίας. Υιοθετήθηκαν πλέον δύο εμπειρικοί δείκτες διαβροχής: ο πρώτος δείκτης, Λ, δηλώνει τη διαβροχή της ενεργούς διεπιφάνειας και ο δεύτερος δηλώνει τη διαβροχή της επιφάνειας του πορώδους μέσου που γεμίζει με υδατική φάση όταν κινηθεί η προαναφερθείσα διεπιφάνεια. Οι δύο δείκτες διαβροχής συσχετίστηκαν με την «ενέργεια» των σημάτων, τη «συχνότητα» των φασμάτων λεπτομερειών του βέλτιστου επιπέδου ανάλυσης των σημάτων και τις χρονομέσες τιμές τριχοειδούς πίεσης. / A transparent porous medium of controlled fractional wettability is fabricated by mixing water-wet glass microspheres with strongly oil-wet polytetrafluouroethylene microspheres and packing them between two transparent glass plates. The displacement experiments performed for two pairs of fluids: a) silicon oil and water and b) paraffin oil and water. The growth pattern is video recorded and the transient response of the pressure drop across the pore network is measured for various fractions of oil-wet particles. The measured global capillary pressure fluctuates as the result of the variation of the equilibrium curvature of menisci between local maxima and local minima. By using wavelet transform on the capillary pressure signal, the best level detail spectrum arises and its peaks indict the parts of the original signal with the most important fluctuations, each of those is analyzed again. The best level detail coefficients give the “energy” of the analyzed part of the capillary pressure signal while the number of the most important fluctuations of the best level detail spectrum to the time interval they took place gives the “frequency” of the best level detail spectrum. The time averaged capillary pressure, the “energy” and the “frequency” are associated with the frontal wettability of the interface separating the two fluids before the invasion step and one number which denotes the regional wettability of the invaded area.

lthough extensive research work has been carried out on the drag reduction behavior of polymers and surfactants alone, little progress has been made on the synergistic effects of combined polymers and surfactants. A number of studies have demonstrated that certain types of polymers and surfactants interact with each other to form surfactant-polymer complexes. The formation of such complexes can cause changes in the solution properties and may result in better drag reduction characteristics as compared with pure additives. A series of drag-reducing surfactants and polymers were screened for the synergistic studies. The following two widely used polymeric drag reducing agents (DRA) were chosen: a copolymer of acrylamide and sodium acrylate (referred to as PAM) and polyethylene oxide (PEO). Among the different types of surfactants screened, a cationic surfactant octadecyltrimethylammonium chloride (OTAC) and an anionic surfactant Sodium dodecyl sulfate (SDS) were selected for the synergistic study. In the case of the cationic surfactant OTAC, sodium salicylate (NaSal) was used as a counterion. No counterion was used with anionic surfactant SDS. The physical properties such as viscosity, surface tension and electrical conductivity were measured in order to detect any interaction between the polymer and the surfactant. The drag reduction (DR) ability of both pure and mixed additives was investigated in a pipeline flow loop. The effects of different parameters such as additive concentration, type of water (deionized (DI) or tap), temperature, tube diameter, and mechanical degradation were investigated. The addition of OTAC to PAM solution has a significant effect on the properties of the system. The critical micelle concentration (CMC) of the mixed surfactant-polymer system is found to be different from that of the surfactant alone. The anionic PAM chains collapse upon the addition of cationic OTAC and a substantial decrease in the viscosity occurs. The pipeline flow behaviour of PAM/OTAC mixtures is found to be consistent with the bench scale results. The drag reduction ability of PAM is reduced upon the addition of OTAC. At low concentrations of PAM, the effect of OTAC on the drag reduction behavior is more pronounced. The drag reduction behavior of polymer solutions is strongly influenced by the nature of water (de-ionized or tap). The addition of OTAC to PEO solution exhibited a week interaction based on the viscosity and surface tension measurements. However, the pipeline results showed a considerable synergistic effect, that is, the mixed system gave a significantly higher drag reduction (lower friction factors) as compared with the pure additives (pure polymer or pure surfactant). The synergistic effect in the mixed system was stronger at low polymer concentrations and high surfactant concentrations. Also the resistance against mechanical degradation of the additive was improved upon the addition of OTAC to PEO. The mixed PEO/SDS system exhibited a strong interaction between the polymers (PEO) and the surfactant (SDS), Using electrical conductivity and surface tension measurements, the critical aggregation concentration (CAC) and the polymer saturation point (PSP) were determined. As the PEO concentration is increased, the CAC decreases and the PSP increase. The addition of SDS to the PEO solution exhibits a remarkable increase in the relative viscosity compared to the pure PEO solution. This increase is attributed to the changes in the hydrodynamic radius of the polymer coil. The pipeline flow exhibited a considerable increase in DR for the mixed system as compared to the pure PEO solution. The addition of surfactant always improves the extent of DR up to the PSP. Also the mixed PEO/ SDS system shows better resistance against shear degradation of the additive.

Τρεις περιοχές της Δυτικής Ελλάδας, στις οποίες είχαν καταγραφεί διαφυγές ρευστών από τον πυθμένα της θάλασσας μελετήθηκαν στα πλαίσια της παρούσας διδακτορικής διατριβής, το υποθαλάσσιο πεδίο κρατήρων διαφυγής ρευστών στον Πατραϊκό κόλπο, το υποθαλάσσιο πεδίο κρατήρων διαφυγής στον Όρμο του Ελαιώνα στο Δυτικό Κορινθιακό κόλπο και οι διαφυγές ρευστών στον κόλπο του Κατακόλου. Η διατριβή βασίζεται στη μελέτη θαλάσσιων γεωφυσικών δεδομένων, στη μελέτη των φυσικοχημικών παραμέτρων της υδάτινης στήλης, στη χημική και ισοτοπική ανάλυση των ρευστών, στη μελέτη των δεδομένων που προέκυψαν από την πρώτη παγκόσμια καταγραφή μακράς περιόδου διαφυγών ρευστών από τον πυθμένα, στη μελέτη των διαφυγών στην παράκτια ζώνη πλησίον των υποθαλάσσιων θέσεων και τέλος στην οπτική παρατήρηση των θέσεων διαφυγών ρευστών. Το πεδίο κρατήρων διαφυγής ρευστών του Πατραϊκού κόλπου σχηματίζεται σε Ολοκαινικές ιλυούχες αποθέσεις που υπέρκεινται Πλειστοκαινικών ιζημάτων, και ελέγχεται από ρηξιγενείς δομές. Το πεδίο κρατήρων διαφυγής ρευστών του Πατραϊκού κόλπου παρουσιάζει μία σχεδόν συνεχή μικροδραστηριότητα η οποία διακόπτεται από παροξυσμικά γεγονότα μεγάλων διαφυγών ρευστών που σχετίζονται με μεγάλους σεισμούς. Το καθεστώς ενεργοποίησης ελέγχεται από το μέγεθος, το επίκεντρο και το είδος του σεισμού. Τα ρευστά που διαφεύγουν είναι αέριοι υδρογονάνθρακες (μεθάνιο μικροβιακής προέλευσης) με ή χωρίς νερό των πόρων των ιζημάτων. Στην παράκτια ζώνη της Πάτρας, πλησίον του πεδίου εντοπίζονται τρεις περιοχές με έντονες διαφυγές αερίων από το έδαφος. Συνολικά η παράκτια ζώνη της Πάτρας, τόσο το χερσαίο όσο και το θαλάσσιο περιβάλλον της, προσφέρει τουλάχιστον 4.7tn μεθανίου ετησίως στην ατμόσφαιρα, με τη μέγιστη προσφορά μεθανίου να μπορεί να φτάσει τους 19 τόννους ετησίως, ενώ μετά από ένα ισχυρό σεισμικό γεγονός εκτιμήθηκε ότι διαφεύγουν 500 κιλά περίπου μεθανίου από το θαλάσσιο πυθμένα σε διάστημα 16 ημερών. Στην περιοχή του Όρμου του Ελαιώνα οι κρατήρες διαφυγής ρευστών σχηματίζονται σε Ολοκαινικές ιλύες και ο πυθμένας των κρατήρων φτάνει στην διεπιφάνεια Ολοκαίνου/Πλειστοκαίνου, με κύριο μηχανισμό σχηματισμού την εκφόρτιση του γλυκού νερού. Οι βενθικές βιοκοινωνίες που εντοπίστηκαν στο εσωτερικό των κρατήρων διαφυγής ρευστών οδηγεί στο συμπέρασμα ότι το υπόγειο νερό, εμπλουτισμένο με θρεπτικά συστατικά και διαλυμένο οξυγόνο δημιουργεί ένα εξαιρετικό οικότοπο για την ανάπτυξη βενθικών οργανισμών. Στην περιοχή του Κατακόλου εντοπίστηκαν και καταγράφηκαν πολύ μεγάλες διαφυγές ρευστών σε τρεις περιοχές, στον Λιμένα Κατακόλου, στη θέση «Φάρος» Κατακόλου και νότια του Λιμένα. Οι έντονες διαφυγές ρευστών, κυρίως θερμογενούς μεθανίου με σημαντικές ποσότητες υδροθείου, φαίνεται να συνδέονται με τις ενεργές διαπυρικές δομές στην περιοχή, που επηρεάζουν τις υπερκείμενες ακολουθίες σχηματίζοντας κανονικά ρήγματα, τα οποία συγκροτούν μία ιδανική δίοδο μετανάστευσης των αερίων από τους Μεσοζωικούς ταμιευτήρες υδρογονανθράκων προς την επιφάνεια. Στην υποθαλάσσια περιοχή του Λιμένα Κατακόλου εκτιμήθηκε ότι η ποσότητα μεθανίου που διαφεύγει από το θαλάσσιο πυθμένα είναι της τάξεως των 1260-1500 τόννων το χρόνο. επικινδυνότητα τόσο για τον άνθρωπο, όσο και τις κατασκευές στην περιοχή του Λιμένα. H υψηλή συγκέντρωση μεθανίου στην ατμόσφαιρα μπορεί να προκαλέσει έκρηξη ή πυρκαγιά, ενώ οι υποθαλάσσιες διαφυγές και η παρουσία αερίων υδρογονανθράκων στα επιφανειακά ιζήματα μπορούν να προκαλέσουν εκρήξεις κατά τη διάρκεια γεωτρήσεων ή να καταστρέψουν θαλάσσιες κατασκευές. Το υδρόθειο είναι το πιο επικίνδυνο και τοξικό γεωρευστό και μπορεί να προκαλέσει σημαντικές ιατρικές βλάβες στον άνθρωπο στις συγκεντρώσεις που μετρήθηκε στην περιοχή του λιμένα. / Three areas in Western Greece studied, where seabed fluid flow have been reported. The dissertation is based on the study of marine geophysical data, on physicochemical parameters of seawater, on chemical and isotopic analysis of fluids, on the study of the data collected during the first long-term monitoring of seabed fluid flow ever done and on the study of the fluid flow on the surroundings inshore areas. The Patras Gulf pockmark field is formed on Holocene muds which overlies Pleistocene coarser sediments and is controlled by normal faults. The Patras Gulf pockmark field is characterized by ‘normal’ micro-activity which is interrupting by high activity periods after strong earthquakes. The fluids that flows from the pockmarks is mainly methane of microbial origin with/or without porewater. An estimation of the total methane that flows from the pockmark field and the surrounding inshore area of the field shows that the methane that leaks to the atmosphere is about 4.7-19 tn/year. After strong earthquakes, when the pockmark field is activated, the estimation of the methane that leaks from the seabed is about 500kgr for a period of 16 days. The Eleonas Bay pockmark field is formed on Holocene muds whereas the base of the pockmarks reaches the Holocene/Pleistocene boundary. The main mechanism for the formation of these pockmarks is the groundwater discharge. Benthic macrofauna that observed on the base of the pockmarks shows that the groundwater which is enriched in nutrients and dissolved oxygen is creating an excellent environment for the growth of benthic fauna. In the area of Katakolo have been observed three areas of enhanced fluid flow offshore and onshore. The studies of these areas have shown that the seeps are caused by thermogenic methane that had accumulated in Mesozoic limestone and had migrated upward through faults, or zones of weakness, induced by salt diapirism. An estimation of the methane that flows from the offshore part of the Katakolo Harbour area have shown that the methane that reaches the atmosphere is about 1260-1500 tn/year. Methane seeps in potentially explosive amounts at the area of the Harbour, and hydrogen sulfide is over the levels necessary to induce toxicological diseases and lethal effects.

La plupart des modules photovolta¨ıques produits sont `a base de silicium.L’efficacit´e de ces modules d´epend fortement de la qualit´e cristalline du siliciumutilis´ee ainsi que de la quantit´e d’impuret´es pr´esente dans le lingotd’origine d’o`u sont issus les modules. Les lingots de silicium sont obtenus aucours d’un proc´ed´e de solidification, durant lequel les impuret´es sont extraitespar ph´enom`ene de s´egr´egation. Le processus de s´egr´egation est influenc´e parl’´ecoulement dans le liquide durant l’´etape de solidification. L’utilisation d’unchamp magn´etique externe permet le contrˆole de l’´ecoulement dans le bainliquide. Dans cette ´etude, l’effet d’un ´ecoulement forc´e sur le processus des´egr´egation est ´etudi´e. Pour cela un dispositif exp´erimental (VB2) et unmod`ele num´erique ont ´et´e d´evelopp´es dans le but de comprendre le rˆole del’´ecoulement sur la forme de l’interface et sur la s´egr´egation des impuret´es. / The photovoltaic modules are generally produced using silicon wafers. Theirelectrical efficiency is related to the crystal quality, which is influenced bythe presence of pollutants in the ingots from which the wafers are cut. Siliconingots are obtained as a result of solidification processes, which implygrowing a crystal from melt. During this solidification process, impurities areseparated from the silicon. The segregation process is greatly influenced bythe melt velocity during the solidification process. The control of the meltflow during the crystallization process can be achieved using external magneticfields. This thesis presents the results of the study on the influence ofthe forced convection induced by a travelling magnetic field (TMF) duringthe solidification process, using both an experimental set-up (VB2) and anumerical model.