In this study, dynamic crack propagation in brittle materials has been studied using a regularized phase field approach.The phase field model used has been validated using specific experimental results of a dynamic in-plane fracture. The crack branching phenomena and existence of a limiting crack tip velocity has been validated using a mode I simulation set-up. A parametric study has also been performed so as to normalize the various numerical parameters that affect the velocity at the crack tip. Following the validation of the phase field model a stochastic analysis of a PBX microstructure has been performed. The microstructure has a high HMX volume fraction of 79\%. The energetic material is HMX and the binder used is Sylgard. Artificial defects are introduced in the system using phase field cracks. The analysis uses a finite element framework that accounts for various thermal-mechanical processes like deformation, heat generation, conduction, fracture and frictional heating at the crack surfaces. The effect on the temperature and damage field due to varying parameters like loading velocities and critical energy release rates is studied. Critical hotspot formation due to localized frictional heating is also studied. A concept of dirty binder is introduced to increase the grain volume fraction of the energetic in the composite. This amounts to a homogenized binder that accounts for the influence of the subsume particles that do not contribute to fracture but affect material properties of the binder.

Este trabalho apresenta um estudo da solidificação de metais puros utilizando o modelo de campo de fases. O modelo é utilizado para simular a solidificação com o intuito de obter a morfologia da interface sólido-líquido sob diversas condições de transferência de calor. Foram realizados testes de validação comparando as morfologias da interface sólido-líquido obtida com as morfologias apresentadas em trabalhos anteriores para os casos bi e tridimensionais. O modelo do campo de fases adotado consiste principalmente de duas equações diferenciais: uma para calcular a variável de campo de fases e outra para calcular o campo de temperaturas. As equações foram solucionadas numericamente para um oitavo do domínio devido a simetria do problema. Os cálculos do modelo indicam que um sólido esférico com um raio inicial menor que o raio crítico de nucleação refunde. Entretanto uma esfera de raio maior cresce. Quando o sólido inicial cresce em uma malha numérica relativamente grosseira, a forma do sólido desvia da forma esférica devido perturbações na interface sólido-líquido. Quando a malha é refinada, as perturbações não são detectadas; contudo, quando introduzidas artificialmente as perturbações crescem e distorcem o formato esférico. / This work presents a study of the solidification of pure metals using the phase field model. The model is used to simulate solidification in order to obtain the morphology of the solid-liquid interface under different heat transfer conditions. Validation tests were performed comparing the morphology of the solid-liquid interface with the morphologies obtained from previous works for two and three dimensional cases. The adopted phase-field model consisted mainly of two differential equations: one to calculate the field of phase variable and another for the temperature field. The equations were solved numerically in only one eighth of the domain owing to the symmetry of the problem. Model calculations show that a solid sphere with an initial radius smaller than the critical radius for nucleation shrinks, whereas a sphere with a larger radius grows. When it grows in a relatively coarse numerical mesh, the initial solid shape deviates from a sphere owing to perturbations at the solid-liquid interface. When the numerical mesh is refined, the growth of perturbations is not detected, but artificially introduced perturbations grow and distort the spherical shape.

L'objectiu de la tesi és estudiar la cinètica de les cristal·litzacions primàries en vidres metàl·lics mitjançant simulacions de tipus phase field. Una cristal·lització primària és una transició de fase sòlid-sòlid on la fase que cristal·litza (fase transformada o fase secundaria) té una composició química diferent de la fase precursora (fase no transformada o fase primària).Les dades experimentals obtingudes a partir de l'estudi calorimètric de cristal·litzacions primàries s'analitzen generalment en el marc del model KJMA (Kolmogorov, Johnson & Mehl, Avrami). Aquest model proporciona l'evolució temporal de la fracció transformada basant-se en tres hipòtesis: - Els nuclis de la fase secundaria estan distribuïts aleatòriament en tot l'espai.- El creixement d'aquests nuclis és isotròpic.- El creixement s'atura únicament per xoc directe (hard impingement).En la cristal·lizació de vidres metàl·lics s'ha observat experimentalment un alentiment de la cinètica respecte del comportament calculat emprant la citada cinètica KJMA. Aquest alentiment s'explica a la literatura en base a que en aquest tipus de transformacions, controlades per difusió, la interacció entre els cristalls no és directa sinó que es produeix a través dels perfils de concentració (soft impingement) i, a més, l'evolució d'aquests perfils de concentració causa canvis en la concentració de la matriu amorfa, estabilitzant la i per tant fent que la nucleació de nous cristalls esdevingui no aleatòria. Diversos autors han proposat modificacions del model KJMA per tal d'intentar superar aquestes limitacions, basats bé en consideracions geomètriques, bé en aproximacions de camp mitjà. A pesar de tot, cap d'aquests models és capaç d'explicar satisfactòriament la cinètica observada en cristal·litzacions primàries. L'objectiu d'aquest treball ha estat la simulació realista de la cinètica de les transformacions primàries per trobar una explicació consistent a les diferències observades entre les dades experimentals i els models teòrics disponibles.Per tal de poder descriure de forma realista el procés de cristal·lització primària s'ha d'estudiar el procés de nucleació i creixement de la fase secundaria alhora que es resol l'equació de difusió en la fase primària. En aquest treball s'ha emprat un model de simulació phase field que permet estudiar aquest sistema introduint una nova variable lligada al camp de concentració que pren dos valors diferents segons es tracti de fase transformada o no transformada. Amb aquest tipus de models també es poden introduir diferents protocols de nucleació i per tant estudiar independentment els efectes de la nucleació en la cinètica. D'aquesta manera s'han realitzat simulacions en 2 i 3 dimensions de cristal·litzacions primàries amb diferents graus de fracció transformada final). Els resultats de les simulacions s'ha comparat amb el model KJMA i, contra el que es preveia, s'ha obtingut un bon acord entre les fraccions transformades del model KJMA i de les simulacions. Donat que el model KJMA no reprodueix satisfactòriament el comportament experimental d'aquest resultat es dedueix que ni el soft impingement ni la nucleació no aleatòria son les responsables de l'alentiment de la cinètica obtingut en cristal·litzacions primàries. Per tal de trobar una explicació físicament convincent del comportament observat experimentalment s'ha aprofundit en l'estudi teòric de les cristali·litzaciones primàries, incloent-hi l'efecte dels canvis composicionals que tenen lloc en la matriu a mesura que la transformació es produeix. Aquest fet, tot i ser conegut a la bibliografia, ha estat sistemàticament ignorat en l'elaboració de models cinètics. En concret, s'ha fet palès que canvis en la composició química de la fase primària han d'afectar de forma radical a la viscositat, que varia fortament a prop de la transició vitrea, i han de produir canvis en les propietats de transport atòmic. Això s'ha modelat a través de l'assumpció d'un coeficient de difusió depenent de la concentració, en base a la relació modificada d'Stokes-Einstein entre la viscositat i el coeficient de difusió. Les simulacions phase-field amb un coeficient de difusió d'aquest tipus donen lloc a una cinètica més lenta i que mostra un acord excel·lent amb la cinètica experimentalment observada en cristal·litzacions primàries de vidres metàl·lics. Per tant, les simulacions phase field confirmen que la cinètica de les cristal·litzacions primàries està controlada fonamentalment pel canvi en les propietats de transport atòmic, mentre que els efectes de soft impingement i nucleació no aleatoria, tot i estar presents, son secundaris. / El objetivo de la tesi es estudiar la cinética de las cristalizaciones primarias en vidrios metálicos mediante simulaciones de tipo phase field. Una cristalización primaria es una transición de fase sólido-sólido donde la fase que cristaliza (fase transformada o fase secundaria) tiene una composición química diferente a la fase precursora (fase no transformada o fase primaria).Los datos experimentales obtenidos a partir del estudio calorimétrico de cristalizaciones primarias se analizan generalmente en el marco del modelo KJMA (Kolmogorov, Johnson & Mehl, Avrami). Este modelo proporciona la evolución temporal de la fracción transformada basándose en tres hipótesis: - Los núcleos de la fase secundaria están distribuidos aleatoriamente en todo el espacio- El crecimiento de estos núcleos es isotrópico- El crecimiento se detiene únicamente por choque directo (hard impingement).En la cristalización de vidrios metálicos se ha observado experimentalmente un retardo de la cinética respecto del comportamiento calculado usando la cinética KJMA. Este retardo se explica en la literatura en base a que en este tipo de transformaciones, controladas por difusión, la interacción entre los cristales no es directa sino que se produce a través de los perfiles de concentración (soft impingement) y, además, la evolución de estos perfiles de concentración causa cambios en la concentración de la matriz amorfa, estabilizándola y por tanto haciendo que la nucleación de nuevos cristales sea no aleatoria. Varios autores han propuesto modificaciones del modelo KJMA para intentar superar estas limitaciones, basados bien en consideraciones geométricas, bien en aproximaciones de campo medio. A pesar de todo, ninguno de estos modelos es capaz de explicar satisfactoriamente la cinética observada en cristalizaciones primarias. El objetivo de este trabajo ha sido la simulación realista de la cinética de las transformaciones primarias para hallar una explicación consistente a las diferencias entre los datos experimentales y los modelos teóricos disponibles.Para describir de manera realista el proceso de cristalización primaria se tiene que estudiar el proceso de nucleación y crecimiento de la fase secundaria a la vez que se resuelve la ecuación de difusión en la fase primaria. En este trabajo se ha usado un modelo de simulación phase-field que permite estudiar este sistema introduciendo una nueva variable ligada al campo de concentración que toma dos valores diferentes según se trate de fase transformada o no transformada. Con este tipo de modelos también se pueden introducir diferentes protocolos de nucleación y por tanto estudiar independientemente los efectos de la nucleación en la cinética. De esta manera se han realizado simulaciones en 2 y 3 dimensiones de cristalizaciones primarias con diferentes grados de fracción transformada final. Los resultados de la simulaciones se han comparado con el modelo KJMA y, en contra de lo que se preveía, se ha obtenido un buen acuerdo entre las fracciones transformadas del modelo KJMA y de las simulaciones. Dado que el modelo KJMA no reproduce satisfactoriamente el comportamiento experimental, de este resultado se deduce que ni el soft impingement ni la nucleación no aleatoria son las responsables del retardo en la cinética obtenido en cristalizaciones primarias.Para encontrar una explicación físicamente convincente del comportamiento observado experimentalmente se ha profundizado en el estudio teórico de las cristalizaciones primarias, incluyendo el efecto de los cambios composicionales que tienen lugar en la matriz a medida que la transformación se produce. Este hecho, aún y ser conocido en la bibliografía, ha sido sistemáticamente ignorado en la elaboración de modelos cinéticos. En concreto, se ha hecho patente que cambios en la composición química de la fase primaria tienen que afectar de forma radical a la viscosidad, que varía fuertemente cerca de la transición vítrea, y tienen que producirse cambios en las propiedades de transporte atómico. Esto se ha modelado a través de la asunción de un coeficiente de difusión dependiente de la concentración, en base a la relación de Stokes-Einstein modificada entre la viscosidad y el coeficiente de difusión. Las simulaciones phsae-field con un coeficiente de difusión de este tipo dan lugar a una cinética más lenta y que muestra un acuerdo excelente con la cinética experimentalmente observada en cristalizaciones primarias de vidrios metálicos. Por tanto, las simulaciones phase-field confirman que la cinética de las cristalizaciones primarias está controlada fundamentalmente por los cambios en las propiedades de transporte atómico, mientras que los efectos de soft-impingement y nucleación no aleatoria, aún y estar presentes, son secundarios. / The aim of this thesis is to study the kinetics of primary crystallization in metallic glasses by means of phase-field simulations. A primary crystallization is a solid-solid phase transformation where the crystallized phase (transformed phase or secondary phase) has a chemical composition different than the precursor phase (untransformed phase or primary phase).Experimental data from calorimetric studies of primary crystallization are usually studied in the framework of the KJMA model (Kolmogorov, Johnson & Mehl, Avrami). This model yields the temporal evolution of the transformed fraction on the basis of three main assumptions: - A random distribution of particle nuclei of the secondary phase- The growth of these nuclei is isotropic- The growth is only halted by direct collisions (hard impingement).In the crystallization of metallic glasses, a slowing down of the kinetics respect the behavior calculated with the KJMA kinetics has been observed. This delay is explained in the literature by the fact that in this kind of transformations, that are diffusion controlled, the interaction between the crystals is not direct but through the concentration profiles (soft impingement) and moreover, the evolution of these profiles causes changes in the concentration of the amorphous matrix, stabilizing it and thus, the nucleation of new nuclei become non random. Several authors had proposed modifications to the KJMA model to try to overcome these limitations, based either on geometrical considerations or in mean field approaches. However, none of these models is able to explain the observed kinetics in primary crystallizations. The aim of this work has been the realistic simulation of the kinetics of primary crystallization to find a explanation to the differences between the experimental data and the available theoretical models.In order to describe in a realistic way the process of a primary crystallization, the nucleation and growth process of the secondary phase has to be studied at the same time that the diffusion equation is solved in the primary phase. In this work, it has been used a phase field model for the simulations that allows to study this system introducing a new variable, coupled to the concentration field, that takes two different values in each of the existing phases. With these kinds of models, different nucleation protocols can also be introduced and thus, independently study the effects of the nucleation in the kinetics. Therefore, 2 and 3 dimensional simulations of primary crystallization have been performed with several degrees of final transformed fraction. The simulation results have been compared with the KJMA model and, unexpectedly, a good agreement between the simulations and the KJMA model has been obtained. As the KJMA model does not reproduce satisfactorily the experimental behavior, from this result can be deduced that neither the soft impingement nor the non random nucleation are the responsible of the slowing down observed in the kinetics of primary crystallization.In order to find a physical convincing explanation of the observed experimental behavior, the theoretical study of primary crystallization has been extended, including the effects of the compositional changes that take place in the matrix as the transformation proceed. This fact, notwithstanding being known in the literature, has been systematically ignored in the development of the kinetics models. In particular, it has become clear that changes in the chemical composition of the primary phase have to radically affect the viscosity, that strongly varies near the glass transition, and some changes in the atomic transport properties must occur. This has been modeled through the assumption of a compositional dependent diffusion coefficient, on the basis of a modified Stokes-Einstein relation between viscosity and diffusion coefficient. Phase field simulations with a diffusion coefficient of this type yield a slower kinetics and show an excellent agreement with the kinetics experimentally observed in primary crystallization of metallic glasses. Thus, phase field simulations confirm that the kinetics of primary crystallization is fundamentally controlled by the changes in the atomic transport properties, while the soft impingement and non random effects, although being present, are secondary.

Este trabalho apresenta um estudo da solidificação de metais puros utilizando o modelo de campo de fases. O modelo é utilizado para simular a solidificação com o intuito de obter a morfologia da interface sólido-líquido sob diversas condições de transferência de calor. Foram realizados testes de validação comparando as morfologias da interface sólido-líquido obtida com as morfologias apresentadas em trabalhos anteriores para os casos bi e tridimensionais. O modelo do campo de fases adotado consiste principalmente de duas equações diferenciais: uma para calcular a variável de campo de fases e outra para calcular o campo de temperaturas. As equações foram solucionadas numericamente para um oitavo do domínio devido a simetria do problema. Os cálculos do modelo indicam que um sólido esférico com um raio inicial menor que o raio crítico de nucleação refunde. Entretanto uma esfera de raio maior cresce. Quando o sólido inicial cresce em uma malha numérica relativamente grosseira, a forma do sólido desvia da forma esférica devido perturbações na interface sólido-líquido. Quando a malha é refinada, as perturbações não são detectadas; contudo, quando introduzidas artificialmente as perturbações crescem e distorcem o formato esférico. / This work presents a study of the solidification of pure metals using the phase field model. The model is used to simulate solidification in order to obtain the morphology of the solid-liquid interface under different heat transfer conditions. Validation tests were performed comparing the morphology of the solid-liquid interface with the morphologies obtained from previous works for two and three dimensional cases. The adopted phase-field model consisted mainly of two differential equations: one to calculate the field of phase variable and another for the temperature field. The equations were solved numerically in only one eighth of the domain owing to the symmetry of the problem. Model calculations show that a solid sphere with an initial radius smaller than the critical radius for nucleation shrinks, whereas a sphere with a larger radius grows. When it grows in a relatively coarse numerical mesh, the initial solid shape deviates from a sphere owing to perturbations at the solid-liquid interface. When the numerical mesh is refined, the growth of perturbations is not detected, but artificially introduced perturbations grow and distort the spherical shape.

The prediction of fracture is of utmost importance regarding the design of modern, specifically tailored engineering materials. These materials are often heterogeneous, \ie their properties vary in space. This can either be achieved purposefully by combining two or more constituents to profit from a more resilient composite material, or happen due to unavoidable imperfections. In any case, purely experimental assessment of failure is tedious and circuitous as soon as the structure of interest gets more complex, and the involvement of numerical models is inevitable. In this work, the phase-field approach to fracture is applied which is able to capture manifold crack phenomena inherently and sidesteps the need for remeshing by describing the crack as a continuous field. The phase-field model is extended to a fully diffuse incorporation of heterogeneities: A static order parameter smoothly transitions from one to the other bulk material constituent, while the weak, brittle interface is incorporated via a continuous fracture toughness distribution. The sharp interface jump conditions still hold for the diffuse representation since a partial rank-1 relaxation is employed in accordance with the unilateral contact condition of the phase-field model. Moreover, a compensation procedure ensures the independence of the interface fracture toughness from interface and phase-field length scales. The model is validated by a comparison to analytical results and the predictive power is demonstrated by the deduction of direction-dependent, effective fracture properties of heterogeneous microstructures. / Die Vorhersage des Bruchverhaltens ist für die Entwicklung moderner, speziell zugeschnittener technischer Werkstoffe von größter Bedeutung. Diese Materialien sind oft heterogen, d.h. ihre Eigenschaften variieren im Raum. Dies kann entweder absichtlich durch die Kombination zweier oder mehrerer Bestandteile erreicht werden, um von einem widerstandsfähigeren Verbundwerkstoff zu profitieren, oder durch unvermeidbare Imperfektionen geschehen. In jedem Fall ist eine rein experimentelle Versagensbewertung komplexer Strukturen mühsam und umständlich, und die Nutzung numerischer Modelle unvermeidlich. In dieser Arbeit werden Brüche mittels Phasenfeldmethode modelliert, wodurch vielfältige Rissphänomene erfasst werden können und die Notwendigkeit einer Neuvernetzung durch die Beschreibung des Risses als kontinuierliches Feld entfällt. Das Phasenfeldmodell wird um eine vollständig diffuse Einbindung von Heterogenitäten erweitert: Ein statischer Orderparameter beschreibt den glatten Übergang zwischen zwei Bestandteilen des Materials, während die geschwächte, spröde Grenzfläche durch eine kontinuierliche Bruchzähigkeitsverteilung eingebunden wird. Die scharfen Grenzflächensprungbedingungen gelten auch für die diffuse Darstellung, da eine partielle Rang-1 Relaxation in Übereinstimmung mit der unilateralen Rissflächenkontaktbedingung genutzt wird. Darüber hinaus gewährleistet ein Kompensationsverfahren die Unabhängigkeit der Grenzflächenbruchzähigkeit von inhärenten Längenskalen der Grenzfläche und des Rissphasenfelds. Das Modell wird durch einen Vergleich mit analytischen Ergebnissen validiert und die Vorhersagekraft wird durch die Ableitung richtungsabhängiger, effektiver Brucheigenschaften heterogener Mikrostrukturen demonstriert.

No description available.

Multiphase flow is a familiar phenomenon from daily life and occupies an important role in physics, engineering, and medicine. The understanding of multiphase flows relies largely on the theory of interfaces, which is not well understood in many cases. To date, the Navier-Stokes-Korteweg equations and the Cahn-Hilliard equation have represented two major branches of phase-field modeling. The Navier-Stokes-Korteweg equations describe a single component fluid material with multiple states of matter, e.g., water and water vapor; the Cahn-Hilliard type models describe multi-component materials with immiscible interfaces, e.g., air and water. In this dissertation, a unified multiphase fluid modeling framework is developed based on rigorous mathematical and thermodynamic principles. This framework does not assume any ad hoc modeling procedures and is capable of formulating meaningful new models with an arbitrary number of different types of interfaces. In addition to the modeling, novel numerical technologies are developed in this dissertation focusing on the Navier-Stokes-Korteweg equations. First, the notion of entropy variables is properly generalized to the functional setting, which results in an entropy-dissipative semi-discrete formulation. Second, a family of quadrature rules is developed and applied to generate fully discrete schemes. The resulting schemes are featured with two main properties: they are provably dissipative in entropy and second-order accurate in time. In the presence of complex geometries and high-order differential terms, isogeometric analysis is invoked to provide accurate representations of computational geometries and robust numerical tools. A novel periodic transformation operator technology is also developed within the isogeometric context. It significantly simplifies the procedure of the strong imposition of periodic boundary conditions. These attributes make the proposed technologies an ideal candidate for credible numerical simulation of multiphase flows. A general-purpose parallel computing software, named PERIGEE, is developed in this work to provide an implementation framework for the above numerical methods. A comprehensive set of numerical examples has been studied to corroborate the aforementioned theories. Additionally, a variety of application examples have been investigated, culminating with the boiling simulation. Importantly, the boiling model overcomes several challenges for traditional boiling models, owing to its thermodynamically consistent nature. The numerical results indicate the promising potential of the proposed methodology for a wide range of multiphase flow problems. / text

This thesis is devoted to several aspects of geometry and morphology in wetting problems and hard sphere packings. First, we propose a new method to simulate wetting and slip on nanostructured substrates: a phase field model associated with a dynamical density theory approach. We showed omniphobicity, meaning repellency, no matter the chemical properties of the liquid on monovalued surfaces, i.e. surfaces without overhangs, which is in contradiction with the macroscopic Cassie-Baxter-Wenzel theory, can produce so-called We checked systematically the impact of the surface parameters on omniphobic repellency, and we show that the key ingredient are line tensions, which emerge from needle shaped surface structures. Geometrical effects have also an important influence on glassy or jammed systems, for example amorphous hard sphere systems in infinite pressure limit. Such hard sphere packings got stuck in a so-called jammed phase, and we shall demonstrate that the local structure in such systems is universal, i.e. independent of the protocol of the generation. For this, robust order parameters - so-called Minkowski tensors - are developed, which overcome robustness deficiencies of widely used order parameters. This leads to a unifying picture of local order parameters, based on geometrical principles. Furthermore, we find with the Minkowski tensor analysis crystallization in jammed sphere packs at the random closed packing point

Phase field models have played an important role in shaping our understanding of a variety of micro structural phenomena in materials. Their attractive features include (a) their ability to capture instabilities in microstructures, and (b) their ability to handle topological transitions { such as splitting or coalescence { gracefully. Therefore, we have chosen to use a phase field model in our study of instabilities in cylindrical pores in nanoporous membranes which eventually lead to their failure. Our study is motivated by recent studies on thermal stability of nanoporous membranes of alumina, titania and zirconia. The key feature in our model is its ability to incorporate surface discussion as the mechanism for mass transport. We first benchmark the model through a critical comparison of our results on early stages of surface evolution during Rayleigh instability and grain boundary grooving with those from linear theories of these phenomena. We have then used longer simulations (which go beyond early stages, and therefore, can incorporate non-lineare effects) to study instabilities in a hollow cylinder in three different systems: single crystal or amorphous solid (which fails through Rayleigh instability), a model sys-tem with parallel grain boundaries (which fails through grain boundary grooving), and a polycrystal (whose failure depends on a combination of grain growth and grooving). In all the cases, the surface energy is assumed to be isotropic, and the operative mechanism for mass transport is assumed to be surface discussion.

We present a new diffuse interface model for the dynamics of inextensible vesicles in a viscous fluid with inertial forces. A new feature of this work is the implementation of the local inextensibility condition in the diffuse interface context. Local inextensibility is enforced by using a local Lagrange multiplier, which provides the necessary tension force at the interface. We introduce a new equation for the local Lagrange multiplier whose solution essentially provides a harmonic extension of the multiplier off the interface while maintaining the local inextensibility constraint near the interface. We also develop a local relaxation scheme that dynamically corrects local stretching/compression errors thereby preventing their accumulation. Asymptotic analysis is presented that shows that our new system converges to a relaxed version of the inextensible sharp interface model. This is also verified numerically. To solve the equations, we use an adaptive finite element method with implicit coupling between the Navier-Stokes and the diffuse interface inextensibility equations. Numerical simulations of a single vesicle in a shear flow at different Reynolds numbers demonstrate that errors in enforcing local inextensibility may accumulate and lead to large differences in the dynamics in the tumbling regime and smaller differences in the inclination angle of vesicles in the tank-treading regime. The local relaxation algorithm is shown to prevent the accumulation of stretching and compression errors very effectively. Simulations of two vesicles in an extensional flow show that local inextensibility plays an important role when vesicles are in close proximity by inhibiting fluid drainage in the near contact region.