Instabilities of a liquid film flowing down an inclined porous plane

The problem of a film flowing down an inclined porous layer is considered. The fully developed basic flow is driven by gravitation. A careful linear instability analysis is carried out. We use Darcy's law to describe the porous layer and solve the coupling equations of the fluid and the porous medium rather than the decoupled equations of the one-sided model used in previous works. The eigenvalue problem is solved by means of a Chebyshev collocation method. We compare the instability of the two-sided model with the results of the one-sided model. The result reveals a porous mode instability which is completely neglected in previous works. For a falling film on an inclined porous plane there are three instability modes, i.e., the surface mode, the shear mode, and the porous mode. We also study the influences of the depth ratio d, the Darcy number delta, and the Beavers-Joseph coefficient alpha(BJ) on the instability of the system.

Improved performance in organic light emitting diodes with a mixed electron donor-acceptor film involved in hole injection

Organic light emitting diodes using a mixed layer of electron acceptor 3, 4, 9, 10 perylenetetracarboxylic dianhydride and electron donor copper phthalocyanine (PTCDA:CuPc) on indium tin oxide (ITO) anodes were fabricated. The device properties were found to be strongly dependent on the thickness of the PTCDA:CuPc film: both the power efficiency and the driving voltage of the device were optimized with a thickness of PTCDA:CuPc ranging from 10 to 20 nm. As compared to the conventional ITO/CuPc hole injection structure, the ITO/PTCDA:CuPc hole injection structure could remarkably enhance both the luminance and the power efficiencies of devices. A mechanism of static-induced, very efficient hole-electron pairs generation in mixed PTCDA:CuPc films was proposed to explain the experimental phenomena. The structural and optical properties of PTCDA:CuPc film were examined as well. (c) 2007 American Institute of Physics.

Instabilities and transient behaviors of a liquid film flowing down a porous inclined plane

The nonmodal linear stability of a falling film over a porous inclined plane has been investigated. The base flow is driven by gravity. We use Darcy's law to describe the flow in the porous medium. A simplified one-sided model is used to describe the fluid flow. In this model, the influence of the porous layer on the flow in the film can be identified by a parameter beta. The instabilities of a falling film have traditionally been investigated by linearizing the governing equations and testing for unstable eigenvalues of the linearized problem. However, the results of eigenvalue analysis agree poorly in many cases with experiments, especially for shear flows. In the present paper, we have studied the linear stability of three-dimensional disturbances using the nonmodal stability theory. Particular attentions are paid to the transient behavior rather than the long time behavior of eigenmodes predicted by traditional normal mode analysis. The transient behaviors of the response to external excitations and the response to initial conditions are studied by examining the pseudospectral structures and the energy growth function G(t) Before we study the nonmodal stability of the system, we extend the results of long-wave analysis in previous works by examining the linear stabilities for streamwise and spanwise disturbances. Results show that the critical conditions of both the surface mode and the shear mode instabilities are dependent on beta for streamwise disturbances. However, the spanwise disturbances have no unstable eigenvalue. 2010 American Institute of Physics. [doi:10.1063/1.3455503]

Physical mechanism and numerical simulations of surface layer temperature inversion in tropical ocean

The one-dimensional Kraus-Turner mixed layer model improved by Liu is developed to consider the effect of salinity and the equations of temperature and salinity under the mixed layer. On this basis, the processes of growth and death of surface layer temperature inversion is numerically simulated under different environmental parameters. At the same time, the physical mechanism is preliminarily discussed combining the observations at the station of TOGA-COARE 0 degrees N, 156 degrees E. The results indicate that temperature inversion sensitively depends on the mixed layer depth, sea surface wind speed and solar shortwave radiation, etc., and appropriately meteorological and hydrological conditions often lead to the similarly periodical occurrence of this inversion phenomenon.

Ocean Boundary Layer Dynamics and Air-Sea Interaction

Thesis (Ph.D.)--University of Washington, 2005-12

Studies on ocean surface layer responses to atmospheric forcing in the North Indian Ocean

The present study on upper ocean responses to atmospheric forcing (associated with cyclone passage) in North Indian Ocean revealed significant variability between AS and BoB. The analysis of cyclone frequency during 1947 to 2006 exhibited lesser frequency of cyclones in AS than that of BoB. The analysis also revealed significant reduction in cyclone frequency after the year 1976 with substantial reduction during monsoon season. The long term SST data at selected points in AS and BoB could not reveal any relation with reduction in cyclone frequency. However the SLP at same locations exhibited considerable increase during mid 1970’s, which could have contributed to the observed reduction in cyclone frequency after the year 1976.The response in waves during cyclone passage exhibited significant asymmetry on either side of the track in AS and BoB and the response is observed at 100’s of kilometers away from the track. The significant clockwise rotation in wave direction is observed on the right side of the track starting from near the track to far away locations, which existed for a longer duration. However, the anticlockwise rotation in wave direction is observed over a shorter distance on the left side of the track and dissipated immediately.Inertial oscillation is observed in surface current and in the mixed layer temperature associated with cyclone passage, which revealed the role of relative location(s) on either side of the track. The inertial peak closer to the local inertial period indicates maximum transfer of energy during the cyclone passage in both AS and BoB. The absence of strong inertial oscillation even with clockwise rotation in surface current and wind indicates the dominant role of duration of strong wind in generating inertial oscillation.The oceanic response associated with cyclone passage reveal the variable response(s) which depends on cyclone intensity, the proximity to track and cyclone translation speed. It is observed that resonance with wind generates higher response in surface current, wave and SST on the right side of the track and it lasts for a longer duration. The maximum oceanic response is observed at a few kilometers away on right side of the track. However lesser rightward bias in the location of maximum cooling is observed for cyclones with low cyclone translation speed. The response on the left side of the track is less and is limited over a shorter distance and dissipates immediately. It is observed that the ocean response, in general, increases with intensity of cyclones. However the differential cooling produced by the same intensity cyclones in AS and in BoB indicates the dominant role of low cyclone translation speed in oceanic response.The surface cooling exhibited strikingly differential responses between AS and BoB. The TMI-SST and buoy observations exhibited significant cooling for a longer duration in AS compared to that of BoB. The spatial extent of cooling is also much higher in AS than that of BoB. The wide spread cooling associated with cyclone passage in AS indicates the dominant role of thermal structure in oceanic response in AS than that of BoB.

Dynamics and mixing in the Upper Ocean Layer

S'estudia la resposta de la capa de barreja oceànica al forçament atmosfèric considerant dades obtingudes durant 12 dies d'abril del 2001 a 42 estacions a través de l'Atlàntic nord seguint aproximadament la latitud de 53ºN. Aquestes dades inclouen, a més de variables atmosfèriques, mesures de CTD, velocitats amb ADCP i dades de microestructura obtingudes amb un perfilador de caiguda lliure. En aquest últim cas, s'han desenvolupat tècniques de processament de les dades que també es presenten aquí. El transsecte estudiat segueix la posició climatològica del rotacional mitjà anual del vent igual a zero i travessa el corrent del Labrador i algunes branques i meandres del Corrent Atlàntic Nord. El forçament atmosfèric es va caracteritzar per vents intensos i fluxos superficials de calor negatius, tot i que, tal com es dedueix de la comparació del gruix de la capa de barreja amb la longitud de Monin-Obukov, la barreja induïda pel vent va dominar sobre la convectiva durant tot el transsecte.

Formation and circulation of the cold intermediate layer in the Gulf of Saint Lawrence

[ 1] The local heat content and formation rate of the cold intermediate layer (CIL) in the Gulf of Saint Lawrence are examined using a combination of new in situ wintertime observations and a three-dimensional numerical model. The field observations consist of five moorings located throughout the gulf over the period of November 2002 to June 2003. The observations demonstrate a substantially deeper surface mixed layer in the central and northeast gulf than in regions downstream of the buoyant surface outflow from the Saint Lawrence Estuary. The mixed-layer depth in the estuary remains shallow (< 60 m) throughout winter, with the arrival of a layer of near-freezing waters between 40 and 100 m depth in April. An eddy-permitting ice-ocean model with realistic forcing is used to hindcast the period of observation. The model simulates well the seasonal evolution of mixed-layer depth and CIL heat content. Although the greatest heat losses occur in the northeast, the most significant change in CIL heat content over winter occurs in the Anticosti Trough. The observed renewal of CIL in the estuary in spring is captured by the model. The simulation highlights the role of the northwest gulf, and in particular, the separation of the Gaspe Current, in controlling the exchange of CIL between the estuary and the gulf. In order to isolate the effects of inflow through the Strait of Belle Isle on the CIL heat content, we examine a sensitivity experiment in which the strait is closed. This simulation shows that the inflow has a less important effect on the CIL than was suggested by previous studies.

The formation of ice in a long-lived supercooled layer cloud

This article focuses on the characteristics of persistent thin single-layer mixed-phase clouds. We seek to answer two important questions: (i) how does ice continually nucleate and precipitate from these clouds, without the available ice nuclei becoming depleted? (ii) how do the supercooled liquid droplets persist in spite of the net flux of water vapour to the growing ice crystals? These questions are answered quantitatively using in situ and radar observations of a long-lived mixed-phase cloud layer over the Chilbolton Observatory. Doppler radar measurements show that the top 500 m of cloud (the top 250 m of which is mixed-phase, with ice virga beneath) is turbulent and well-mixed, and the liquid water content is adiabatic. This well-mixed layer is bounded above and below by stable layers. This inhibits entrainment of fresh ice nuclei into the cloud layer, yet our in situ and radar observations show that a steady flux of ≈100 m−2s−1 ice crystals fell from the cloud over the course of ∼1 day. Comparing this flux to the concentration of conventional ice nuclei expected to be present within the well-mixed layer, we find that these nuclei would be depleted within less than 1 h. We therefore argue that nucleation in these persistent supercooled clouds is strongly time-dependent in nature, with droplets freezing slowly over many hours, significantly longer than the few seconds residence time of an ice nucleus counter. Once nucleated, the ice crystals are observed to grow primarily by vapour deposition, because of the low liquid water path (21 g m−2) yet vapour-rich environment. Evidence for this comes from high differential reflectivity in the radar observations, and in situ imaging of the crystals. The flux of vapour from liquid to ice is quantified from in situ measurements, and we show that this modest flux (3.3 g m−2h−1) can be readily offset by slow radiative cooling of the layer to space.

Impact of Laptev Sea flaw polynyas on the atmospheric boundary layer and ice production using idealized mesoscale simulations

The interaction between polynyas and the atmospheric boundary layer is examined in the Laptev Sea using the regional, non-hydrostatic Consortium for Small-scale Modelling (COSMO) atmosphere model. A thermodynamic sea-ice model is used to consider the response of sea-ice surface temperature to idealized atmospheric forcing. The idealized regimes represent atmospheric conditions that are typical for the Laptev Sea region. Cold wintertime conditions are investigated with sea-ice–ocean temperature differences of up to 40 K. The Laptev Sea flaw polynyas strongly modify the atmospheric boundary layer. Convectively mixed layers reach heights of up to 1200 m above the polynyas with temperature anomalies of more than 5 K. Horizontal transport of heat expands to areas more than 500 km downstream of the polynyas. Strong wind regimes lead to a more shallow mixed layer with strong near-surface modifications, while weaker wind regimes show a deeper, well-mixed convective boundary layer. Shallow mesoscale circulations occur in the vicinity of ice-free and thin-ice covered polynyas. They are forced by large turbulent and radiative heat fluxes from the surface of up to 789 W m−2, strong low-level thermally induced convergence and cold air flow from the orographic structure of the Taimyr Peninsula in the western Laptev Sea region. Based on the surface energy balance we derive potential sea-ice production rates between 8 and 25 cm d−1. These production rates are mainly determined by whether the polynyas are ice-free or covered by thin ice and by the wind strength.

Langmuir turbulence and surface heating in the ocean surface boundary layer

This study uses large-eddy simulation to investigate the structure of the ocean surface boundary layer (OSBL) in the presence of Langmuir turbulence and stabilizing surface heat fluxes. The OSBL consists of a weakly stratified layer, despite a surface heat flux, above a stratified thermocline. The weakly stratified (mixed) layer is maintained by a combination of a turbulent heat flux produced by the wave-driven Stokes drift and downgradient turbulent diffusion. The scaling of turbulence statistics, such as dissipation and vertical velocity variance, is only affected by the surface heat flux through changes in the mixed layer depth. Diagnostic models are proposed for the equilibrium boundary layer and mixed layer depths in the presence of surface heating. The models are a function of the initial mixed layer depth before heating is imposed and the Langmuir stability length. In the presence of radiative heating, the models are extended to account for the depth profile of the heating.

STUDY OF THE EQUATORIAL ATLANTIC OCEAN MIXING LAYER USING A ONE-DIMENSIONAL TURBULENCE MODEL

The General Ocean Turbulence Model (GOTM) is applied to the diagnostic turbulence field of the mixing layer (ML) over the equatorial region of the Atlantic Ocean. Two situations were investigated: rainy and dry seasons, defined, respectively, by the presence of the intertropical convergence zone and by its northward displacement. Simulations were carried out using data from a PIRATA buoy located on the equator at 23 degrees W to compute surface turbulent fluxes and from the NASA/GEWEX Surface Radiation Budget Project to close the surface radiation balance. A data assimilation scheme was used as a surrogate for the physical effects not present in the one-dimensional model. In the rainy season, results show that the ML is shallower due to the weaker surface stress and stronger stable stratification; the maximum ML depth reached during this season is around 15 m, with an averaged diurnal variation of 7 m depth. In the dry season, the stronger surface stress and the enhanced surface heat balance components enable higher mechanical production of turbulent kinetic energy and, at night, the buoyancy acts also enhancing turbulence in the first meters of depth, characterizing a deeper ML, reaching around 60 m and presenting an average diurnal variation of 30 m.

Submesoscale activity over the Argentinian shelf

Submesoscale activity over the Argentinian shelf is investigated by means of high resolution primitive equation numerical solutions. These reveal energetic turbulent activity (visually similar to the one occasionally seen in satellite images) at scales O(5 km) in fall and winter that is linked to mixed layer baroclinic instability. The air-sea heat flux responsible for (i) deepening the upper ocean boundary layer (at these seasons) and (ii) maintaining a cross-shelf background density gradient is the key environmental parameter controlling submesoscale activity. Implications of submesoscale turbulence are investigated. Its mixing efficiency estimated by computing a diffusivity coefficient is above 30 m(2) s(-1) away from the shallowest regions. Aggregation of surface buoyant material by submesoscale currents occurs within hours and is presumably important to the ecosystem.

Coupling between upper ocean layer variability and size-fractionated phytoplankton in a non-nutrient-limited environment

[EN] We describe the coupling between upper ocean layer variability and size-fractionated phytoplankton distribution in the non-nutrient-limited Bransfield Strait region (BS) of Antarctica. For this purpose we use hydrographic and size-fractionated chlorophyll a data from a transect that crossed 2 fronts and an eddy, together with data from 3 stations located in a deeply mixed region, the Antarctic Sound (AS). In the BS transect, small phytoplankton (<20 μm equivalent spherical diameter [ESD]) accounted for 80% of total chl a and their distribution appeared to be linked to cross-frontal variability. On the deepening upper mixed layer (UML) sides of both fronts we observed a deep subducting column-like structure of small phytoplankton biomass. On the shoaling UML sides of both fronts, where there were signs of restratification, we observed a local shallow maximum of small phytoplankton biomass. We propose that this observed phytoplankton distribution may be a response to the development of frontal vertical circulation cells. In the deep, turbulent environment of the AS, larger phytoplankton (>20 μm ESD) accounted for 80% of total chl a. The proportion of large phytoplankton increases as the depth of the upper mixed layer (ZUML), and the corresponding rate of vertical mixing, increases. We hypothesize that this change in phytoplankton composition with varying ZUML is related to the competition for light, and results from modification of the light regime caused by vertical mixing.