The interaction of riblets with wall-bounded turbulence

The purpose of this thesis is to give answer to the question: why do riblets stop working for a certain size? Riblets are small surface grooves aligned in the mean direction of an overlying turbulent flow, designed specifically to reduce the friction between the flow and the surface. They were inspired by biological surfaces, like the oriented denticles in the skin of fastswimming sharks, and were the focus of a significant amount of research in the late eighties and nineties. Although it was found that the drag reduction depends on the riblet size scaled in wall units, the physical mechanisms implicated have not been completely understood up to now. It has been explained how riblets of vanishing size interact with the turbulent flow, producing a change in the drag proportional to their size, but that is not the regime of practical interest. The optimum performance is achieved for larger sizes, once that linear behavior has broken down, but before riblets begin adopting the character of regular roughness and increasing drag. This regime, which is the most relevant from a technological perspective, was precisely the less understood, so we have focused on it. Our efforts have followed three basic directions. First, we have re-assessed the available experimental data, seeking to identify common characteristics in the optimum regime across the different existing riblet geometries. This study has led to the proposal of a new length scale, the square root of the groove crosssection, to substitute the traditional peak-to-peak spacing. Scaling the riblet dimension with this length, the size of breakdown of the linear behavior becomes roughly universal. This suggests that the onset of the breakdown is related to a certain, fixed value of the cross-section of the groove. Second, we have conducted a set of direct numerical simulations of the turbulent flow over riblets, for sizes spanning the full drag reduction range. We have thus been able to reproduce the gradual transition between the different regimes. The spectral analysis of the flows has proven particularly fruitful, since it has made possible to identify spanwise rollers immediately above the riblets, which begin to appear when the riblet size is close to the optimum. This is a quite surprising feature of the flow, not because of the uniqueness of the phenomenon, which had been reported before for other types of complex and porous surfaces, but because most previous studies had focused on the detail of the flow above each riblet as a unit. Our novel approach has provided the adequate tools to capture coherent structures with an extended spanwise support, which interact with the riblets not individually, but collectively. We have also proven that those spanwise structures are responsible for the increase in drag past the viscous breakdown. Finally, we have analyzed the stability of the flow with a simplified model that connects the appearance of rollers to a Kelvin–Helmholtz-like instability, as is the case also for the flow over plant canopies and porous surfaces. In spite of the model emulating the presence of riblets only in an averaged, general fashion, it succeeds to capture the essential attributes of the breakdown, and provides a theoretical justification for the scaling with the groove cross-section.

A coupled hydro-mechanical analysis for prediction of hydraulic fracture propagation in saturated porous media using EFG mesh-less method

The details of the Element Free Galerkin (EFG) method are presented with the method being applied to a study on hydraulic fracturing initiation and propagation process in a saturated porous medium using coupled hydro-mechanical numerical modelling. In this EFG method, interpolation (approximation) is based on nodes without using elements and hence an arbitrary discrete fracture path can be modelled.The numerical approach is based upon solving two governing partial differential equations of equilibrium and continuity of pore water simultaneously. Displacement increment and pore water pressure increment are discretized using the same EFG shape functions. An incremental constrained Galerkin weak form is used to create the discrete system of equations and a fully implicit scheme is used for discretization in the time domain. Implementation of essential boundary conditions is based on the penalty method. In order to model discrete fractures, the so-called diffraction method is used.Examples are presented and the results are compared to some closed-form solutions and FEM approximations in order to demonstrate the validity of the developed model and its capabilities. The model is able to take the anisotropy and inhomogeneity of the material into account. The applicability of the model is examined by simulating hydraulic fracture initiation and propagation process from a borehole by injection of fluid. The maximum tensile strength criterion and Mohr-Coulomb shear criterion are used for modelling tensile and shear fracture, respectively. The model successfully simulates the leak-off of fluid from the fracture into the surrounding material. The results indicate the importance of pore fluid pressure in the initiation and propagation pattern of fracture in saturated soils. © 2013 Elsevier Ltd.

A zero-stiffness elastic shell structure

A remarkable shell structure is described that, due to a particular combination of geometry and initial stress, has zero stiffness for any finite deformation along a twisting path; the shell is in a neutrally stable state of equilibrium. Initially the shell is straight in a longitudinal direction, but has a constant, nonzero curvature in the transverse direction. If residual stresses are induced in the shell by, for example, plastic deformation, to leave a particular resultant bending moment, then an analytical inextensional model of the shell shows it to have no change in energy along a path of twisted configurations. Real shells become closer to the inextensional idealization as their thickness is decreased; experimental thin-shell models have confirmed the neutrally stable configurations predicted by the inextensional theory. A simple model is described that shows that the resultant bending moment that leads to zero stiffness gives the shell a hidden symmetry, which explains this remarkable property.

Ultrafast collinear scattering and carrier multiplication in graphene.

Graphene is emerging as a viable alternative to conventional optoelectronic, plasmonic and nanophotonic materials. The interaction of light with charge carriers creates an out-of-equilibrium distribution, which relaxes on an ultrafast timescale to a hot Fermi-Dirac distribution, that subsequently cools emitting phonons. Although the slower relaxation mechanisms have been extensively investigated, the initial stages still pose a challenge. Experimentally, they defy the resolution of most pump-probe setups, due to the extremely fast sub-100 fs carrier dynamics. Theoretically, massless Dirac fermions represent a novel many-body problem, fundamentally different from Schrödinger fermions. Here we combine pump-probe spectroscopy with a microscopic theory to investigate electron-electron interactions during the early stages of relaxation. We identify the mechanisms controlling the ultrafast dynamics, in particular the role of collinear scattering. This gives rise to Auger processes, including charge multiplication, which is key in photovoltage generation and photodetectors.

The inverse problem in magnetic force microscopy--inferring sample magnetization from MFM images.

Nanomagnetic structures have the potential to surpass silicon's scaling limitations both as elements in hybrid CMOS logic and as novel computational elements. Magnetic force microscopy (MFM) offers a convenient characterization technique for use in the design of such nanomagnetic structures. MFM measures the magnetic field and not the sample's magnetization. As such the question of the uniqueness of the relationship between an external magnetic field and a magnetization distribution is a relevant one. To study this problem we present a simple algorithm which searches for magnetization distributions consistent with an external magnetic field and solutions to the micromagnetic equations' qualitative features. The algorithm is not computationally intensive and is found to be effective for our test cases. On the basis of our results we propose a systematic approach for interpreting MFM measurements.

Ultrafast non-thermal electron dynamics in single layer graphene

The impulsive optical excitation of carriers in graphene creates an out-of-equilibrium distribution, which thermalizes on an ultrafast timescale [1-4]. This hot Fermi-Dirac (FD) distribution subsequently cools via phonon emission within few hundreds of femtoseconds. While the relaxation mechanisms mediated by phonons have been extensively investigated, the initial stages, ruled by fundamental electron-electron (e-e) interactions still pose a challenge. © 2013 IEEE.