Entropy generation for microscale forced convection: effects of different thermal boundary conditions, velocity slip, temperature jump, viscous dissipation, and duct geometry

This work presents closed form solutions for fully developed temperature distribution and entropy generation due to forced convection in microelectromechanical systems (MEMS) in the Slip-flow regime, for which the Knudsen number lies within the range 0.001

Combined Conduction-Convection-Radiation Heat Transfer of Slip Flow inside a Micro-channel Filled with a Cellular Porous Material

Combined conduction–convection–radiation heat transfer is investigated numerically in a micro-channel filled with a saturated cellular porous medium, with the channel walls held at a constant heat flux. Invoking the velocity slip and temperature jump, the thermal behaviour of the porous–fluid system are studied by considering hydrodynamically fully developed flow and applying the Darcy–Brinkman flow model. One energy equation model based on the local thermal equilibrium condition is adopted to evaluate the temperature field within the porous medium. Combined conduction and radiation heat transfer is treated as an effective conduction process with a temperature-dependent effective thermal conductivity. Results are reported in terms of the average Nusselt number and dimensionless temperature distribution, as a function of velocity slip coefficient, temperature jump coefficient, porous medium shape parameter and radiation parameters. Results show that increasing the radiation parameter (Tr)(Tr) and the temperature jump coefficient flattens the dimensionless temperature profile. The Nusselt numbers are more sensitive to the variation in the temperature jump coefficient rather than to the velocity slip coefficient. Such that for high porous medium shape parameter, the Nusselt number is found to be independent of velocity slip. Furthermore, it is found that as the temperature jump coefficient increases, the Nusselt number decrease. In addition, for high temperature jump coefficients, the Nusselt number is found to be insensitive to the radiation parameters and porous medium shape parameter. It is also concluded that compared with the conventional macro-channels, wherein using a porous material enhances the rate of heat transfer (up to about 40 % compared to the clear channel), insertion of a porous material inside a micro-channel in slip regime does not effectively enhance the rate of heat transfer that is about 2 %.

Modeling of Heat Transfer in Microchannel Gas Flow

Due to the existence of a velocity slip and temperature jump on the solid walls, the heat transfer in microchannels significantly differs from the one in the macroscale. In our research, we have focused on the pressure driven gas flows in a simple finite microchannel geometry, with an entrance and an outlet, for low Reynolds (Re<200) and low Knudsen (Kn<0.01) numbers. For such a regime, the slip induced phenomena are strongly connected with the viscous effects. As a result, heat transfer is also significantly altered. For the optimization of flow conditions, we have investigated various temperature gradient configurations, additionally changing Reynolds and Knudsen numbers. The entrance effects, slip flow, and temperature jump lead to complex relations between flow behavior and heat transfer. We have shown that slip effects are generally insignificant for flow behavior. However, two configuration setups (hot wall cold gas and cold wall hot gas) are affected by slip in distinguishably different ways. For the first one, which concerns turbomachinery, the mass flow rate can increase by about 1% in relation to the no-slip case, depending on the wall-gas temperature difference. Heat transfer is more significantly altered. The Nusselt number between slip and no-slip cases at the outlet of the microchannel is increased by about 10%.

Laboratory shearing experiments at low sliding velocities of IODP Hole 343-C0019E

The 2011 Tohoku-Oki earthquake demonstrated that the shallowest reaches of plate boundary subduction megathrusts can host substantial coseismic slip that generates large and destructive tsunamis, contrary to the common assumption that the frictional properties of unconsolidated clay-rich sediments at depths less than View the MathML source should inhibit rupture. We report on laboratory shearing experiments at low sliding velocities (View the MathML source) using borehole samples recovered during IODP Expedition 343 (JFAST), spanning the plate-boundary décollement within the region of large coseismic slip during the Tohoku earthquake. We show that at sub-seismic slip rates the fault is weak (sliding friction µs=0.2-0.26), in contrast to the much stronger wall rocks (µs>~0.5). The fault is weak due to elevated smectite clay content and is frictionally similar to a pelagic clay layer of similar composition. The higher cohesion of intact wall rock samples coupled with their higher amorphous silica content suggests that the wall rock is stronger due to diagenetic cementation and low clay content. Our measurements also show that the strongly developed in-situ fabric in the fault zone does not contribute to its frictional weakness, but does lead to a near-cohesionless fault zone, which may facilitate rupture propagation by reducing shear strength and surface energy at the tip of the rupture front. We suggest that the shallow rupture and large coseismic slip during the 2011 Tohoku earthquake was facilitated by a weak and cohesionless fault combined with strong wall rocks that drive localized deformation within a narrow zone.

Constant wall heat flux boundary conditions in micro-channels filled with porous material with internal heat generation under local thermal non-equilibrium conditions

Forced convection heat transfer in a micro-channel filled with a porous material saturated with rarefied gas with internal heat generation is studied analytically in this work. The study is performed by analysing the boundary conditions for constant wall heat flux under local thermal non-equilibrium (LTNE) conditions. Invoking the velocity slip and temperature jump, the thermal behaviour of the porous-fluid system is studied by considering thermally and hydrodynamically fully-developed conditions. The flow inside the porous material is modelled by the Darcy–Brinkman equation. Exact solutions are obtained for both the fluid and solid temperature distributions for two primary approaches models A and B using constant wall heat flux boundary conditions. The temperature distributions and Nusselt numbers for models A and B are compared, and the limiting cases resulting in the convergence or divergence of the two models are also discussed. The effects of pertinent parameters such as fluid to solid effective thermal conductivity ratio, Biot number, Darcy number, velocity slip and temperature jump coefficients, and fluid and solid internal heat generations are also discussed. The results indicate that the Nusselt number decreases with the increase of thermal conductivity ratio for both models. This contrasts results from previous studies which for model A reported that the Nusselt number increases with the increase of thermal conductivity ratio. The Biot number and thermal conductivity ratio are found to have substantial effects on the role of temperature jump coefficient in controlling the Nusselt number for models A and B. The Nusselt numbers calculated using model A change drastically with the variation of solid internal heat generation. In contrast, the Nusselt numbers obtained for model B show a weak dependency on the variation of internal heat generation. The velocity slip coefficient has no noticeable effect on the Nusselt numbers for both models. The difference between the Nusselt numbers calculated using the two models decreases with an increase of the temperature jump coefficient.

Measurements of velocity-dependent and slip-dependent frictional strength in laboratory friction experiments on natural samples from IODP Expedition 316 and ODP Hole 190-1174B

Slowslip forms part of the spectrum of fault behaviour between stable creep and destructive earthquakes. Slow slip occurs near the boundaries of large earthquake rupture zones and may sometimes trigger fast earthquakes. It is thought to occur in faults comprised of rocks that strengthen under fast slip rates, preventing rupture as a normal earthquake, or on faults that have elevated pore-fluid pressures. However, the processes that control slow rupture and the relationship between slow and normal earthquakes are enigmatic. Here we use laboratory experiments to simulate faulting in natural rock samples taken from shallow parts of the Nankai subduction zone, Japan, where very low-frequency earthquakes - a form of slow slip - have been observed.We find that the fault rocks exhibit decreasing strength over millimetre-scale slip distances rather than weakening due to increasing velocity. However, the sizes of the slip nucleation patches in our laboratory simulations are similar to those expected for the very lowfrequency earthquakes observed in Nankai. We therefore suggest that this type of fault-weakening behaviour may generate slow earthquakes. Owing to the similarity between the expected behaviour of slow earthquakes based on our data, and that of normal earthquakes during nucleation, we suggest that some types of slow slip may represent prematurely arrested earthquakes.

Velocity fluctuations of fracturelike disruptions of associating polymer solutions

Velocity has been measured as a function of time for propagating crack tips as water is injected into solutions of end-capped associating polymers in a rectanguar Hele-Shaw cell. Measurements were performed for flows with different values of cell gap, channel width, polymer molecular weight, and polymer concentration. The condition for the onset of fracturelike behavior is well described by a Deborah number which uses the shear-thinning shear rate of the polymer solution as a characteristic frequency for network relaxation. At low molecular weight, the onset of fracturelike pattern evolution is accompanied by an abrupt jump in tip velocity, followed by a lower and approximately constant acceleration. At high molecular weight, the transition to fracturelike behavior involves passing through a regime that may be understood in terms of stick-slip dynamics. The crack-tip wanders from side to side and fluctuates (in both speed and velocity along the channel) with a characteristic frequency which depends linearly on the invading fluid injection rate.

Algebraic decay of velocity fluctuations near a wall

Computer simulations of the dynamics of a colloidal particle suspended in a fluid confined by an interface show that the asymptotic decay of the velocity correlation functions is algebraic. The exponents of the long-time tails depend on the direction of motion of the particle relative to the surface, as well as on the specific nature of the boundary conditions. In particular, we find that for the angular velocity correlation function, the decay in the presence of a slip surface is faster than the one corresponding to a stick one. An intuitive picture is introduced to explain the various long-time tails, and the simulations are compared with theoretical expressions where available.

Velocity structure and seismic anisotropy in the crust and upper mantle from Receiver Function analysis: three case studies in Italy

The research for this PhD project consisted in the application of the RFs analysis technique to different data-sets of teleseismic events recorded at temporary and permanent stations located in three distinct study regions: Colli Albani area, Northern Apennines and Southern Apennines. We found some velocity models to interpret the structures in these regions, which possess very different geologic and tectonics characteristics and therefore offer interesting case study to face. In the Colli Albani some of the features evidenced in the RFs are shared by all the analyzed stations: the Moho is almost flat and is located at about 23 km depth, and the presence of a relatively shallow limestone layer is a stable feature; contrariwise there are features which vary from station to station, indicating local complexities. Three seismic stations, close to the central part of the former volcanic edifice, display relevant anisotropic signatures­­­ with symmetry axes consistent with the emplacement of the magmatic chamber. Two further anisotropic layers are present at greater depth, in the lower crust and the upper mantle, respectively, with symmetry axes directions related to the evolution of the volcano complex. In Northern Apennines we defined the isotropic structure of the area, finding the depth of the Tyrrhenian (almost 25 km and flat) and Adriatic (40 km and dipping underneath the Apennines crests) Mohos. We determined a zone in which the two Mohos overlap, and identified an anisotropic body in between, involved in the subduction and going down with the Adiratic Moho. We interpreted the downgoing anisotropic layer as generated by post-subduction delamination of the top-slab layer, probably made of metamorphosed crustal rocks caught in the subduction channel and buoyantly rising toward the surface. In the Southern Apennines, we found the Moho depth for 16 seismic stations, and highlighted the presence of an anisotropic layer underneath each station, at about 15-20 km below the whole study area. The moho displays a dome-like geometry, as it is shallow (29 km) in the central part of the study area, whereas it deepens peripherally (down to 45 km); the symmetry axes of anisotropic layer, interpreted as a layer separating the upper and the lower crust, show a moho-related pattern, indicated by the foliation of the layer which is parallel to the Moho trend. Moreover, due to the exceptional seismic event occurred on April 6th next to L’Aquila town, we determined the Vs model for two station located next to the epicenter. An extremely high velocity body is found underneath AQU station at 4-10 km depth, reaching Vs of about 4 km/s, while this body is lacking underneath FAGN station. We compared the presence of this body with other recent works and found an anti-correlation between the high Vs body, the max slip patches and earthquakes distribution. The nature of this body is speculative since such high velocities are consistent with deep crust or upper mantle, but can be interpreted as a as high strength barrier of which the high Vs is a typical connotation.

Total internal reflection fluorescence cross-correlation spectroscopy: theory and application for studying boundary slip phenomenon

I present a new experimental method called Total Internal Reflection Fluorescence Cross-Correlation Spectroscopy (TIR-FCCS). It is a method that can probe hydrodynamic flows near solid surfaces, on length scales of tens of nanometres. Fluorescent tracers flowing with the liquid are excited by evanescent light, produced by epi-illumination through the periphery of a high NA oil-immersion objective. Due to the fast decay of the evanescent wave, fluorescence only occurs for tracers in the ~100 nm proximity of the surface, thus resulting in very high normal resolution. The time-resolved fluorescence intensity signals from two laterally shifted (in flow direction) observation volumes, created by two confocal pinholes are independently measured and recorded. The cross-correlation of these signals provides important information for the tracers’ motion and thus their flow velocity. Due to the high sensitivity of the method, fluorescent species with different size, down to single dye molecules can be used as tracers. The aim of my work was to build an experimental setup for TIR-FCCS and use it to experimentally measure the shear rate and slip length of water flowing on hydrophilic and hydrophobic surfaces. However, in order to extract these parameters from the measured correlation curves a quantitative data analysis is needed. This is not straightforward task due to the complexity of the problem, which makes the derivation of analytical expressions for the correlation functions needed to fit the experimental data, impossible. Therefore in order to process and interpret the experimental results I also describe a new numerical method of data analysis of the acquired auto- and cross-correlation curves – Brownian Dynamics techniques are used to produce simulated auto- and cross-correlation functions and to fit the corresponding experimental data. I show how to combine detailed and fairly realistic theoretical modelling of the phenomena with accurate measurements of the correlation functions, in order to establish a fully quantitative method to retrieve the flow properties from the experiments. An importance-sampling Monte Carlo procedure is employed in order to fit the experiments. This provides the optimum parameter values together with their statistical error bars. The approach is well suited for both modern desktop PC machines and massively parallel computers. The latter allows making the data analysis within short computing times. I applied this method to study flow of aqueous electrolyte solution near smooth hydrophilic and hydrophobic surfaces. Generally on hydrophilic surface slip is not expected, while on hydrophobic surface some slippage may exists. Our results show that on both hydrophilic and moderately hydrophobic (contact angle ~85°) surfaces the slip length is ~10-15nm or lower, and within the limitations of the experiments and the model, indistinguishable from zero.