Influence Of Free Surface Curvature Of A Liquid Layer On The Critical Marangoni Convection

The Pearson instability was suggested to discuss the onset of Marangoni convection in a liquid layer of large Prandtl number under an applied temperature difference perpendicular to the free surface in the microgravity environment. In this case, the temperature distribution on the curved free surface is nonuniform, and the thermocapillary convection is induced and coupled with the Marangoni convection. In the present paper the effect of volume ratio of the liquid layer on the critical Marangoni convection and the corresponding spatial variation of the convection structure in zero-gravity condition were numerically investigated by two-dimensional model. (C) 2008 Elsevier Ltd. All rights reserved.

Dispersion properties of surface waves in a two-layer plasma structure bounded by a metal

The effect of the nonuniformity of the electron density on the dispersion properties of surface waves propagating in a direction transverse to an external magnetic field is studied for the model of a two-layer plasma structure bounded by a metal. It is shown that the spectra of the waves can be effectively controlled by varying the degree of nonuniformity of the density and the dimensions of the layers.

On the Kelvin-Helmholtz instability of structured plasma layers in the magnetosphere

The hydromagnetic Kelvin-Helmholtz (K-H) instability problem is studied for a three-layered system analytically by arriving at the marginal instability condition. As the magnetic field directions are taken to vary in the three regions, both the angle and finite thickness effects are seen on the instability criterion. When the relative flow speed of the plasmas on the two sides of the interfaces separating the inner and the surrounding layers is U < Uc, where Uc is the critical speed, the system is stable both for symmetric and asymmetric perturbations. However, unlike the case of the interface bounded by two semiinfinite media, Uc is no longer the minimum critical speed above which the system will be unstable for all wavenumbers; another critical speed U* > Uc is introduced due to the finiteness of the system. When Uc < U < U*, the instability can set in either through the symmetric or asymmetric mode, depending on the ratio of the plasma parameters and angle between the magnetic field directions across the boundaries. The instability arises for a finite range of wavenumbers, thus giving rise to the upper and lower cut-off frequencies for the spectra of hydromagnetic surface waves generated by the K-H instability mechanism. When U > U*, both the modes are unstable for short wavelengths. The results are finally used to explain some observational features of the dependence of hydromagnetic energy spectra in the magnetosphere on the interplanetary parameters.

Resonance cones and mode conversion in a warm magnetized bounded plasma

The warm plasma resonance cone structure of the quasistatic field produced by a gap source in a bounded magnetized slab plasma is determined theoretically. This is initially determined for a homogeneous or mildly inhomogeneous plasma with source frequency lying between the lower hybrid frequency and the plasma frequency. It is then extended to the complicated case of an inhomogeneous plasma with two internal lower hybrid layers present, which is of interest to radio frequency heating of plasmas.

In the first case, the potential is obtained as a sum of multiply reflected warm plasma resonance cones, each of which has a similar structure, but a different size, amplitude, and position. An important interference between nearby multiply-reflected resonance cones is found. The cones are seen to spread out as they move away from the source, so that this interference increases and the individual resonance cones become obscured far away from the source.

In the second case, the potential is found to be expressible as a sum of multiply-reflected, multiply-tunnelled, and mode converted resonance cones, each of which has a unique but similar structure. The effects of both collisional and collisionless damping are included and their effects on the decay of the cone structure studied. Various properties of the cones such as how they move into and out of the hybrid layers, through the evanescent region, and transform at the hybrid layers are determined. It is found that cones can tunnel through the evanescent layer if the layer is thin, and the effect of the thin evanescent layer is to subdue the secondary maxima of cone relative to the main peak, while slightly broadening the main peak and shifting it closer to the cold plasma cone line.

Energy theorems for quasistatic fields are developed and applied to determine the power flow and absorption along the individual cones. This reveals the points of concentration of the flow and the various absorption mechanisms.

Scattering of electromagnetic waves by rough interfaces and inhomogeneous layers

We consider a two-dimensional problem of scattering of a time-harmonic electromagnetic plane wave by an infinite inhomogeneous conducting or dielectric layer at the interface between semi-infinite homogeneous dielectric half-spaces. The magnetic permeability is assumed to be a fixed positive constant. The material properties of the media are characterized completely by an index of refraction, which is a bounded measurable function in the layer and takes positive constant values above and below the layer, corresponding to the homogeneous dielectric media. In this paper, we examine only the transverse magnetic (TM) polarization case. A radiation condition appropriate for scattering by infinite rough surfaces is introduced, a generalization of the Rayleigh expansion condition for diffraction gratings. With the help of the radiation condition the problem is reformulated as an equivalent mixed system of boundary and domain integral equations, consisting of second-kind integral equations over the layer and interfaces within the layer. Assumptions on the variation of the index of refraction in the layer are then imposed which prove to be sufficient, together with the radiation condition, to prove uniqueness of solution and nonexistence of guided wave modes. Recent, general results on the solvability of systems of second kind integral equations on unbounded domains establish existence of solution and continuous dependence in a weighted norm of the solution on the given data. The results obtained apply to the case of scattering by a rough interface between two dielectric media and to many other practical configurations.

“Unconformity-Bounded” Stratigraphic Units in the South Polar Layered Deposits (Promethei Lingula, Mars)

We conducted a stratigraphic analysis of the South Polar Layered Deposits (SPLDs) in Promethei Lingula (PL, Mars) based on the identification of regional unconformities at visible and radar wavelengths. According to the terrestrial classification, this approach constrains the stratigraphy of the region and remedies the ambiguous interpretation of stratigraphy through marker layers, bypassing the problem related to the morphologic and radiometric appearance of the layers. Thus, the approach does not exclude diverse classifications, but complements them, whereas other discriminant elements are doubtful or difficult/impossible to be defined. Using this approach, we defined two stratigraphic units (or synthems: PL1 and PL2) in PL, which are morphologically different and divided by a regional unconformity (AuR1). This stratigraphic architecture implies that the geological history of PL has been conditioned by periodic changes in climate, which in turn are related to orbital variations of Mars.

Convective instabilities in two superposed horizontal liquid layers heated laterally.

This work is devoted to the theoretical study of the stability of two superposed horizontal liquid layers bounded by two solid planes and subjected to a horizontal temperature gradient. The liquids are supposed to be immiscible with a nondeformable interface. The forces acting on the system are buoyancy and interfacial tension. Four different flow patterns and temperature profiles are found for the basic state. A linear perturbative analysis with respect to two- and three-dimensional perturbations reveals the existence of three kinds of patterns. Depending on the relative height of both liquids several situations are predicted: either wave propa- gation from cold to the hot regions, or waves propagating in the opposite direction or still stationary longitu- dinal rolls. The behavior of three different pairs of liquids which have been used in experiments on bilayers under vertical gradient by other authors have been examined. The instability mechanisms are discussed and a qualitative interpretation of the different behaviors exhibited by the system is provided. In some configurations it is possible to find a codimension-two point created by the interaction of two Hopf modes with different frequencies and wave numbers. These results suggest to consider two liquid layers as an interesting prototype ? nard-Marangoni problem.

High reynolds numbers turbulent boundary layers

Esta tesis estudia las similitudes y diferencias entre los flujos turbulentos de pared de tipo externo e interno, en régimen incompresible, y a números de Reynolds moderada¬mente altos. Para ello consideramos tanto simulaciones numéricas como experimentos de capas límites con gradiente de presiones nulo y de flujos de canal, ambos a números de Reynolds en el rango δ+ ~ 500 - 2000. Estos flujos de cortadura son objeto de numerosas investigaciones debido a la gran importancia que tienen tanto a nivel tecnológico como a nivel de física fundamental. No obstante, todavía existen muchos interrogantes sobre aspectos básicos tales como la universalidad de los perfiles medios y de fluctuación de las velocidades o de la presión, tanto en la zona cercana a la pared como en la zona logarítmica, el escalado y el efecto del número de Reynolds, o las diferencias entre los flujos internos y externos en la zona exterior. En éste estudio hemos utilizado simulaciones numéricas ya existentes de canales y capas límites a números de Reynolds δ+ ~ 2000 y δ+ ~ 700, respectivamente. Para poder comparar ambos flujos a igual número de Reynolds hemos realizado una nueva simulación directa de capa límite en el rango δ+ ~ 1000-2000. Los resultados de la misma son presentados y analizados en detalle. Los datos sin postprocesar y las estadísticas ya postprocesadas están públicamente disponibles en nuestro sitio web.162 El análisis de las estadísticas usando un único punto confirma la existencia de perfiles logarítmicos para las fluctuaciones de la velocidad transversal w'2+ y de la presión p'2+ en ambos tipos de flujos, pero no para la velocidad normal v'2+ o la velocidad longitudinal u'2+. Para aceptar o rechazar la existencia de un rango logarítmico en u'2+ se requieren números de Reynolds más altos que los considerados en éste trabajo. Una de las conse¬cuencias más importantes de poseer tales perfiles es que el valor máximo de la intensidad, que se alcanza cerca de la pared, depende explícitamente del número de Reynolds. Esto ha sido confirmado tras analizar un gran número de datos experimentales y numéricos, cor¬roborando que el máximo de u'2+, p/2+, y w'2+ aumenta proporcionalmente con el log(δ+). Por otro lado, éste máximo es más intenso en los flujos externos que en los internos. La máxima diferencia ocurre en torno a y/δ ~ 0.3-0.5, siendo esta altura prácticamente independiente del número de Reynolds considerado. Estas diferencias se originan como consecuencia del carácter intermitente de las capas límites, que es inexistente en los flujos internos. La estructura de las fluctuaciones de velocidad y de presión, junto con la de los esfuer¬zos de Reynolds, se han investigado por medio de correlaciones espaciales tridimensionales considerando dos puntos de medida. Hemos obtenido que el tamaño de las mismas es gen¬eralmente mayor en canales que en capas límites, especialmente en el caso de la correlación longitudinal Cuu en la dirección del flujo. Para esta correlación se demuestra que las es¬tructuras débilmente correladas presentan longitudes de hasta 0(75), en el caso de capas límites, y de hasta 0(185) en el caso de canales. Estas longitudes se obtienen respecti-vamente en la zona logarítmica y en la zona exterior. Las longitudes correspondientes en la dirección transversal son significativamente menores en ambos flujos, 0(5 — 25). La organización espacial de las correlaciones es compatible con la de una pareja de rollos casi paralelos con dimensiones que escalan en unidades exteriores. Esta organización se mantiene al menos hasta y ~ 0.65, altura a la cual las capas límites comienzan a organi¬zarse en rollos transversales. Este comportamiento es sin embargo más débil en canales, pudiéndose observar parcialmente a partir de y ~ 0.85. Para estudiar si estas estructuras están onduladas a lo largo de la dirección transver¬sal, hemos calculado las correlaciones condicionadas a eventos intensos de la velocidad transversal w'. Estas correlaciones revelan que la ondulación de la velocidad longitudinal aumenta conforme nos alejamos de la pared, sugiriendo que las estructuras están más alineadas en la zona cercana a la pared que en la zona lejana a ella. El por qué de esta ondulación se encuentra posiblemente en la configuración a lo largo de diagonales que presenta w'. Estas estructuras no sólo están onduladas, sino que también están inclinadas respecto a la pared con ángulos que dependen de la variable considerada, de la altura, y de el contorno de correlación seleccionado. Por encima de la zona tampón e independien¬temente del número de Reynolds y tipo de flujo, Cuu presenta una inclinación máxima de unos 10°, las correlaciones Cvv y Cm son esencialmente verticales, y Cww está inclinada a unos 35°. Summary This thesis studies the similitudes and differences between external and internal in¬compressible wall-bounded turbulent flows at moderately-high Reynolds numbers. We consider numerical and experimental zero-pressure-gradient boundary layers and chan¬nels in the range of δ+ ~ 500 — 2000. These shear flows are subjects of intensive research because of their technological importance and fundamental physical interest. However, there are still open questions regarding basic aspects such as the universality of the mean and fluctuating velocity and pressure profiles at the near-wall and logarithmic regions, their scaling and the effect of the Reynolds numbers, or the differences between internal and external flows at the outer layer, to name but a few. For this study, we made use of available direct numerical simulations of channel and boundary layers reaching δ+ ~ 2000 and δ+ ~ 700, respectively. To fill the gap in the Reynolds number, a new boundary layer simulation in the range δ+ ~ 1000-2000 is presented and discussed. The original raw data and the post-processed statistics are publicly available on our website.162 The analysis of the one-point statistic confirms the existence of logarithmic profiles for the spanwise w'2+ and pressure p'2+ fluctuations for both type of flows, but not for the wall-normal v'2+ or the streamwise u'2+ velocities. To accept or reject the existence of a logarithmic range in u'2+ requires higher Reynolds numbers than the ones considered in this work. An important consequence of having such profiles is that the maximum value of the intensities, reached near the wall, depends on the Reynolds number. This was confirmed after surveying a wide number of experimental and numerical datasets, corrob¬orating that the maximum of ul2+, p'2+, and w'2+ increases proportionally to log(δ+). On the other hand, that maximum is more intense in external flows than in internal ones, differing the most around y/δ ~ 0.3-0.5, and essentially independent of the Reynolds number. We discuss that those differences are originated as a consequence of the inter¬mittent character of boundary layers that is absent in internal flows. The structure of the velocity and pressure fluctuations, together with those of the Reynolds shear stress, were investigated using three-dimensional two-point spatial correlations. We find that the correlations extend over longer distances in channels than in boundary layers, especially in the case of the streamwise correlation Cuu in the flow direc-tion. For weakly correlated structures, the maximum streamwise length of Cuu is O(78) for boundary layers and O(188) for channels, attained at the logarithmic and outer regions respectively. The corresponding lengths for the transverse velocities and for the pressure are shorter, 0(8 — 28), and of the same order for both flows. The spatial organization of the velocity correlations is shown to be consistent with a pair of quasi-streamwise rollers that scales in outer units. That organization is observed until y ~ 0.68, from which boundary layers start to organize into spanwise rollers. This effect is weaker in channels, and it appears at y ~ 0.88. We present correlations conditioned to intense events of the transversal velocity, w', to study if these structures meander along the spanwise direction. The results indicate that the streamwise velocity streaks increase their meandering proportionally to the distance to the wall, suggesting that the structures are more aligned close to the wall than far from it. The reason behind this meandering is probably due to the characteristic organization along diagonals of w'. These structures not only meander along the spanwise direction, but they are also inclined to the wall at angles that depend on the distance from the wall, on the variable being considered, and on the correlation level used to define them. Above the buffer layer and independent of the Reynolds numbers and type of flow, the maximum inclination of Cuu is about 10°, Cvv and Cpp are roughly vertical, and Cww is inclined by 35°.

Pattern competition in homogeneously heated fluid layers

Simulations examining pattern competition have been performed on a horizontal homogeneously heated layer that is bounded by an isothermal plane above an adiabatic plane. Several different circulation patterns arose as the heating regime applied to the horizontal layer was modified. The sequence of the patterns formed as the Grashof number was increased had the following order: laminar, z-axis rolls, squares, hexagons and pentagons, pentagons and then two different square modes of differing orientations. Fourier analysis was used to determine how the key modes interact in the presence of different patterns.

Pattern competition in homogeneously heated fluid layers

Simulations examining pattern competition have been performed on a horizontal homogeneously heated layer that is bounded by an isothermal plane above an adiabatic plane. Several different circulation patterns arose as the heating regime applied to the horizontal layer was modified. The sequence of the patterns formed as the Grashof number was increased had the following order: laminar layer, rolls, squares, hexagons and pentagons, and then two square modes of differing orientations. Fourier analysis was used to determine how the key modes interact with each pattern.

Flow of a micropolar fluid bounded by a stretching sheet

We consider boundary layer flow of a micropolar fluid driven by a porous stretching sheet. A similarity solution is defined, and numerical solutions using Runge-Kutta and quasilinearisation schemes are obtained. A perturbation analysis is also used to derive analytic solutions to first order in the perturbing parameter. The resulting closed form solutions involve relatively complex expressions, and the analysis is made more tractable by a combination of offline and online work using a computational algebra system (CAS). For this combined numerical and analytic approach, the perturbation analysis yields a number of benefits with regard to the numerical work. The existence of a closed form solution helps to discriminate between acceptable and spurious numerical solutions. Also, the expressions obtained from the perturbation work can provide an accurate description of the solution for ranges of parameters where the numerical approaches considered here prove computationally more difficult.