Flux comparison of Eulerian and Lagrangian estimates of Agulhas leakage: A case study using a numerical model

Estimating the magnitude of Agulhas leakage, the volume flux of water from the Indian to the Atlantic Ocean, is difficult because of the presence of other circulation systems in the Agulhas region. Indian Ocean water in the Atlantic Ocean is vigorously mixed and diluted in the Cape Basin. Eulerian integration methods, where the velocity field perpendicular to a section is integrated to yield a flux, have to be calibrated so that only the flux by Agulhas leakage is sampled. Two Eulerian methods for estimating the magnitude of Agulhas leakage are tested within a high-resolution two-way nested model with the goal to devise a mooring-based measurement strategy. At the GoodHope line, a section halfway through the Cape Basin, the integrated velocity perpendicular to that line is compared to the magnitude of Agulhas leakage as determined from the transport carried by numerical Lagrangian floats. In the first method, integration is limited to the flux of water warmer and more saline than specific threshold values. These threshold values are determined by maximizing the correlation with the float-determined time series. By using the threshold values, approximately half of the leakage can directly be measured. The total amount of Agulhas leakage can be estimated using a linear regression, within a 90% confidence band of 12 Sv. In the second method, a subregion of the GoodHope line is sought so that integration over that subregion yields an Eulerian flux as close to the float-determined leakage as possible. It appears that when integration is limited within the model to the upper 300 m of the water column within 900 km of the African coast the time series have the smallest root-mean-square difference. This method yields a root-mean-square error of only 5.2 Sv but the 90% confidence band of the estimate is 20 Sv. It is concluded that the optimum thermohaline threshold method leads to more accurate estimates even though the directly measured transport is a factor of two lower than the actual magnitude of Agulhas leakage in this model.

(Table 1) Ice flux, polynya area and sea-ice drift speed in the Kara Sea during winter (Jan-Apr) 1997-2001

The Shelf Seas of the Arctic are known for their large sea-ice production. This paper presents a comprehensive view of the Kara Sea sea-ice cover from high-resolution numerical modeling and space-borne microwave radiometry. As given by the latter the average polynya area in the Kara Sea takes a value of 21.2 × 10**3 km**2 ± 9.1 × 10**3 km**2 for winters (Jan.-Apr.) 1996/97 to 2000/01, being as high as 32.0 × 10**3 km**2 in 1999/2000 and below 12 × 10**3 km**2 in 1998/99. Day-to-day variations of the Kara Sea polynya area can be as high as 50 × 10**3 km**2. For the seasons 1996/97 to 2000/01 the modeled cumulative winter ice-volume flux out of the Kara Sea varied between 100 km**3/a and 350 km**3/a. Modeled high (low) ice export coincides with a high (low) average and cumulative polynya area, and with a low (high) sea-ice compactness in the Kara Sea from remote sensing data, and with a high (low) sea-ice drift speed across its northern boundary derived from independent model data for the winters 1996/97 to 2000/01.

Pulsatile flow of a viscous fluid in elliptical pipe of variable cross-section

The pulsatile flow of an incompressible viscous fluid in an elliptical pipe of slowly varying cross-section is considered. Asymptotic series solutions for the velocity distribution and pressure gradient are obtained in terms of Mathieu functions for a low Reynold number flow in which the volume flux is prescribed. An expression for shear stress on the boundary is derived. The physically significant quantities governing the flow are computed numerically and analysed for different types of constrictions. The effect of eccentricity and Womerslay parameter on the flow is discussed.

Discrete multicomponent model for biodiesel spray combustion simulation

An important first step in spray combustion simulation is an accurate determination of the fuel properties which affects the modelling of spray formation and reaction. In a practical combustion simulation, the implementation of a multicomponent model is important in capturing the relative volatility of different fuel components. A Discrete Multicomponent (DM) model is deemed to be an appropriate candidate to model a composite fuel like biodiesel which consists of four components of fatty acid methyl esters (FAME). In this paper, the DM model is compared with the traditional Continuous Thermodynamics (CTM) model for both diesel and biodiesel. The CTM model is formulated based on mixing rules that incorporate the physical and thermophysical properties of pure components into a single continuous surrogate for the composite fuel. The models are implemented within the open-source CFD code OpenFOAM, and a semi-quantitative comparison is made between the predicted spray-combustion characteristics and optical measurements of a swirl-stabilised flame of diesel and biodiesel. The DM model performs better than the CTM model in predicting a higher magnitude of heat release rate in the top flame brush region of the biodiesel flame compared to that of the diesel flame. Using both the DM and CTM models, the simulation successfully reproduces the droplet size, volume flux, and droplet density profiles of diesel and biodiesel. The DM model predicts a longer spray penetration length for biodiesel compared to that of diesel, as seen in the experimental data. Also, the DM model reproduces a segregated biodiesel fuel vapour field and spray in which the most abundant FAME component has the longest vapour penetration. In the biodiesel flame, the relative abundance of each fuel component is found to dominate over the relative volatility in terms of the vapour species distribution and vice versa in the liquid species distribution. © 2014 Elsevier Ltd. All rights reserved.

Scaling arguments for the fluxes in turbulent miscible fountains

For established axisymmetric turbulent miscible Boussinesq fountains in quiescent uniform environments, expressions are developed for the fluxes of volume, momentum and buoyancy at the outflow from the fountain: the outflow referring to the counterflow at the horizontal plane of the source. The fluxes are expressed in terms of the fountain source conditions and two dimensionless functions of the source Froude number, Fr0: a radial function (relating a horizontal scale of the outflow to the source radius) and a volume flux function (relating the outflow and source volume fluxes). The forms taken by these two functions at low Fr0 and high Fr0 are deduced, thereby providing the outflow fluxes and outflow Froude number solely in terms of the source conditions. For high Fr0, the outflow Froude number, Frout, is shown to be invariant, indicating (by analogy with plumes for which the 'far-field' Froude number is invariant with source Froude number) that the outflow may be regarded as 'far-field' since the fluxes within the fountain have adjusted to attain a balance which is independent of the source conditions. Based on Frout, the fluxes in the plume that forms beyond the fountain outflow are deduced. Finally, from the results of previously published studies, we show that the scalings deduced for fountains are valid for 0.0025 ≲ Fr0 ≲ 1.0 for low Fr0 and Fr0≳ 3.0 for high Fr0. © 2014 Cambridge University Press.

Development of microvalves for gas flow control in micronozzles using PVDF piezoelectric polymer

This work describes a fabrication and test sequence of microvalves installed on micronozzles. The technique used to fabricate the micronozzles was powder blasting. The microvalves are actuators made from PVDF (polivinylidene fluoride), that is a piezoelectric polymer. The micronozzles have convergent-divergent shape with external diameter of 1mm and throat around 230μm. The polymer have low piezoelectric coefficient, for this reason a bimorph structure with dimensions of 2mm width and 4mm of length was build (two piezoelectric sheets were glued together with opposite polarization). Both sheets are recovered with a conductor thin film used as electrodes. Applying a voltage between the electrodes one sheet expands while the other contracts and this generate a vertical movement to the entire actuator. Appling +300V DC between the electrodes the volume flux rate, for a pressure ratio of 0.5, was 0.36 sccm. Applying -200V DC between the electrodes (that means it closed) the volume flux rate was 0.32 sccm, defining a possible range of flow between 0.32 and 0.36 sccm. The third measurement was performed using AC voltage (200V AC with frequency of 1Hz), where the actuator was oscillating. For pressure ratio of 0.5, the flow rate was 0.62 sccm. © 2008 IOP Publishing Ltd.

Scherverhalten kolloidaler Lagenphasen in Rohrströmungen

Kolloidale Suspensionen aus identischen kugelförmigen, geladenen Partikeln in wässrigen Medien stellen ein ideales Modellsystem zur Untersuchung des Gleichgewichtsverhaltens, aber auch des Nicht-Gleichgewichtsverhaltens Weicher Materie dar. So bilden derartige Systeme bei hinreichend starker und langreichweitiger elektrostatischer Repulsion fluid und kristallin geordnete Strukturen aus, die wegen der weitreichenden Analogie zu atomar kondensierter Materie als kolloidale Fluide und Kristalle bezeichnet werden. Von großem Vorteil ist dabei die Möglichkeit zur kontrollierten Einstellung der Wechselwirkung und die gute optische Zugänglichkeit für Mikroskopie und Lichtstreuung sowie die Weichheit der Materialien, aufgrund derer sich auch Zustände fernab des mechanischen Gleichgewichts gezielt präparieren lassen. Themenstellung der vorliegenden Arbeit ist die Untersuchung des Phasenverhaltens und der Fließmechanismen kolloidaler Kristalle in einer Rohrströmung. Im ersten Teil der Arbeit wird gezeigt, dass beim Fluss durch eine zylindrische Röhre Mehrphasenkoexistenz auftritt, wobei ein polykristalliner Kern von einer isotropen Scherschmelze umgeben ist. Zusätzlich treten an der Grenze zwischen diesen Phasen und an der Rohrwand Phasen hexagonal geordneter übereinander hinweggleitender Lagen auf. Der Vergleich zwischen auf der Basis der Navier-Stokes-Gleichung theoretisch berechneten und gemessenen Geschwindigkeitsprofilen zeigt, dass jede dieser Phasen für sich Newtonsches Fließverhalten aufweist. Die Gesamtviskosität ist hingegen durch die mit dem Durchsatz veränderliche Phasenzusammensetzung Nicht-Newtonsch. Damit gelang es, die erstmalig von Würth beschriebene Scherverdünnung auf eine Veränderung der Phasenzusammensetzung zurückzuführen. Im zweiten Teil der Arbeit wurde erstmals das Fließverhalten der Lagenphasen mittels Lichtstreuung und Korrelationsanalyse untersucht. Dafür wurde ein im Prinzip einfacher, aber leistungsstarker Aufbau realisiert, der es erlaubt, die zeitliche Veränderung der Bragg-Reflexe der Lagenphase in radialer und azimutaler Richtung zu verfolgen und mittels Fourieranalyse zu analysieren. In Abhängigkeit vom Durchsatz geht die zunächst rastend gleitende Lagenphase in eine frei gleitende Lagenphase über, wobei charakteristische Veränderungen der Spektren sowie der Korrelationsfunktionen auftreten, die detailliert diskutiert werden. Der Übergang im Gleitmechanismus ist mit einem Verlust der Autokorrelation der Rotationskomponente der periodischen Intra-Lagenverzerrung verbunden, während die Kompressionskomponente erhalten bleibt. Bei hohen Durchflüssen lassen die Reflexbewegungen auf das Auftreten einer Eigenschwingung der frei gleitenden Lagen schließen. Diese Schwingung lässt sich als Rotationsbewegung, gekoppelt mit einer transversalen Auslenkung in Vortexrichtung, beschreiben. Die Ergebnisse erlauben eine detaillierte Diskussion von verschiedenen Modellvorstellungen anderer Autoren.

Kuroshio Intrusion and the Circulation in the South China Sea

The Princeton Ocean Model is used to study the circulation in the South China Sea (SCS) and its seasonal transition. Kuroshio enters ( leaves) the SCS through the southern ( northern) portion of the Luzon Strait. The annually averaged net volume flux through the Luzon Strait is similar to2 Sv into the SCS with seasonal reversals. The inflow season is from May to January with the maximum intrusion of Kuroshio water reaching the western SCS during fall in compensation of summertime surface offshore transport associated with coastal upwelling. From February to April the net transport reverses from the SCS to the Pacific. The intruded Kuroshio often forms an anticyclonic current loop west of the Luzon Strait. The current loop separates near the Dongsha Islands with the northward branch continuously feeding the South China Sea Warm Current (SCSWC) near the shelf break and the westward branch becoming the South China Sea Branch of Kuroshio on the slope, which is most apparent in the fall. The SCSWC appears from December to February on the seaward side of the shelf break, flowing eastward against the prevailing wind. Diagnosis shows that the onshore Ekman transport due to northeasterly monsoon generates upwelling when moving upslope, and the particular distributions of the density and sea level associated with the cross shelf motion supports the SCSWC.

Viscous-Flow Approach to In Situ Infiltration and In Vitro Saturated Hydraulic Conductivity Determination

Infiltration is dominantly gravity driven, and a viscous-flow approach was developed. Laminar film flow equilibrates gravity with the viscous force and a constant flow velocity evolves during a period lasting 3/2 times the duration of a constant input rate, qS. Film thickness F and the specific contact area L of the film per unit soil volume are the key parameters. Sprinkler irrigation produced in situ time series of volumetric water contents, θ(z,t), as determined with TDR probes. The wetting front velocity v and the time series of the mobile water content, w(z,t) were deduced from θ(z,t). In vitro steady flow in a core of saturated soil provided volume flux density, q(z,t), and flow velocity, v, as determined from a heat front velocity. The F and L parameters of the in situ and the in vitro experiments were compared. The macropore-flow restriction states that, for a particular permeable medium, the specific contact area L must be independent from qS i.e., dL/dqS = 0. If true, then the relationship of qS ∝ v3/2 could scale a wide range of input rates 0 ≤ qS ≤ saturated hydraulic conductivity, Ksat, into a permeable medium, and kinematic-wave theory would become a versatile tool to deal with non-equilibrium flow. The viscous-flow approach is based on hydromechanical principles similar to Darcy’s law, but currently it is not suited to deduce flow properties from specified individual spatial structures of permeable media.

Accumulation, velocities, fluxes and mass balance of glacial tributaries to the Lambert Glacier-Amery Ice Shelf system, Antarctica

By incorporating recently available remote sensing data, we investigated the mass balance for all individual tributary glacial basins of the Lambert Glacier-Amery Ice Shelf system, East Antarctica. On the basis of the ice flow information derived from SAR interferometry and ICESat laser altimetry, we have determined the spatial configuration of eight tributary drainage basins of the Lambert-Amery glacial system. By combining the coherence information from SAR interferometry and the texture information from SAR and MODIS images, we have interpreted and refined the grounding line position. We calculated ice volume flux of each tributary glacial basin based on the ice velocity field derived from Radarsat three-pass interferometry together with ice thickness data interpolated from Australian and Russian airborne radio echo sounding (RES) surveys and inferred from ICESat laser altimetry data. Our analysis reveals that three tributary basins have a significant net positive imbalance, while five other subbasins are slightly positive or close to zero balance. Overall, in contrast to previous studies, we find that the grounded ice in Lambert Glacier-Amery Ice Shelf system has a positive mass imbalance of 22.9 ± 4.4 Gt/a. The net basal melting for the entire Amery Ice Shelf is estimated to be 27.0 ± 7.0 Gt/a. The melting rate decreases rapidly from the grounding zone to the ice shelf front. Significant basal refreezing is detected in the downstream section of the ice shelf. The mass balance estimates for both the grounded ice sheet and the ice shelf mass differ substantially from other recent estimates.

Sedimentology and Cd/Ca ratios of benthic foraminifera from ODP Site 172-1060 in the western Nort Atlantic

The Western Boundary Undercurrent (WBUC), off eastern America, is an important component of the Atlantic Meridional Overturning circulation and is the principal route for southward transport of North Atlantic waters and southward return of Southern Source Water (SSW). Here a direct flow speed proxy (mean grain size of the sortable silt) is used to infer the vigour of flow of the palaeo-WBUC at Blake Outer Ridge, (ODP Site 1060, depth 3481 m) during Marine Isotope Stage (MIS) 3. The overall shape of the flow speed proxy record shows a complex pattern of variability, with generally more vigorous flow and larger-scale flow variations between 35 and 60 ka than in the younger part of MIS 3 and MIS 2 (b35 ka). Six events of reduced bottom flow vigour (Slow Events, SEs) occur. These appear uncorrelated with Heinrich events, but are instead synchronous with the warming phases of Antarctic Warm Events A-1 to A-4 (with one new one, A-1a and one poorly defined, 'A-0'). This indicates that Antarctic climate exerts a stronger control on deep flow vigour in the North Atlantic during MIS 3 than Northern Hemisphere climate. The correspondence of SEs with Antarctic warming suggests a weaker WBUC flow due to reduced volume flux at SSW source or reduced SSW density. Because the variability of the lower limb of the WBUC was not connected to sharp North Atlantic changes in temperature, it is unlikely that the Dansgaard/Oeschger cycles were associated with a mode of MOC variation involving wholeocean overturn, but more likely with perturbations of only the shallow Glacial Gulf Stream-Glacial Northern Source Intermediate Water cell. Nutrient proxies (benthic carbon isotopes and Cd/Ca of Uvigerina peregrina) at this site show similar trends to the GRIP delta18O record. This correlation has previously been attributed mainly to hydrographic and flow changes but is here shown to be better explained by variations in surface ocean productivity and subsequent decomposition of 12C rich organic material on the sea floor.

Fluid mechanical and electrostatic study on the osmotic flow through cylindrical pores

When two solutions differing in solute concentration are separated by a porous membrane, the osmotic pressure will generate a net volume flux of the suspending fluid across the membrane; this is termed osmotic flow. We consider the osmotic flow across a membrane with circular cylindrical pores when the solute and the pore walls are electrically charged, and the suspending fluid is an electrolytic solution containing small cations and anions. Under the condition in which the radius of the pores and that of the solute molecules greatly exceed those of the solvent as well as the ions, a fluid mechanical and electrostatic theory is introduced to describe the osmotic flow in the presence of electric charge. The interaction energy, including the electrostatic interaction between the solute and the pore wall, plays a key role in determining the osmotic flow. We examine the electrostatic effect on the osmotic flow and discuss the difference in the interaction energy determined from the nonlinear Poisson-Boltzmann equation and from its linearized equation (the Debye-Hückel equation).

Finite volume evolution Galerkin method for hyperbolic conservation laws with spatially varying flux functions

We present a generalization of the finite volume evolution Galerkin scheme [M. Lukacova-Medvid'ova,J. Saibertov'a, G. Warnecke, Finite volume evolution Galerkin methods for nonlinear hyperbolic systems, J. Comp. Phys. (2002) 183 533-562; M. Luacova-Medvid'ova, K.W. Morton, G. Warnecke, Finite volume evolution Galerkin (FVEG) methods for hyperbolic problems, SIAM J. Sci. Comput. (2004) 26 1-30] for hyperbolic systems with spatially varying flux functions. Our goal is to develop a genuinely multi-dimensional numerical scheme for wave propagation problems in a heterogeneous media. We illustrate our methodology for acoustic waves in a heterogeneous medium but the results can be generalized to more complex systems. The finite volume evolution Galerkin (FVEG) method is a predictor-corrector method combining the finite volume corrector step with the evolutionary predictor step. In order to evolve fluxes along the cell interfaces we use multi-dimensional approximate evolution operator. The latter is constructed using the theory of bicharacteristics under the assumption of spatially dependent wave speeds. To approximate heterogeneous medium a staggered grid approach is used. Several numerical experiments for wave propagation with continuous as well as discontinuous wave speeds confirm the robustness and reliability of the new FVEG scheme.

Pulsed-field magnetization of drilled bulk high-temperature superconductors: flux front propagation in the volume and on the surface

We present a method for characterizing the propagation of the magnetic flux in an artificially drilled bulk high-temperature superconductor (HTS) during a pulsed-field magnetization. As the magnetic pulse penetrates the cylindrical sample, the magnetic flux density is measured simultaneously in 16 holes by means of microcoils that are placed across the median plane, i.e. at an equal distance from the top and bottom surfaces, and close to the surface of the sample. We discuss the time evolution of the magnetic flux density in the holes during a pulse and measure the time taken by the external magnetic flux to reach each hole. Our data show that the flux front moves faster in the median plane than on the surface when penetrating the sample edge; it then proceeds faster along the surface than in the bulk as it penetrates the sample further. Once the pulse is over, the trapped flux density inside the central hole is found to be about twice as large in the median plane than on the surface. This ratio is confirmed by modelling.

Study by Hall probe mapping of the trapped flux modification produced by local heating in YBCO HTS bulks for different surface/volume ratios

The aim of this report is to compare the trapped field distribution under a local heating created at the sample edge for different sample morphologies. Hall probe mappings of the magnetic induction trapped in YBCO bulk samples maintained out of thermal equilibrium were performed on YBCO bulk single domains, YBCO single domains with regularly spaced hole arrays, and YBCO superconducting foams. The capability of heat draining was quantified by two criteria: the average induction decay and the size of the thermally affected zone caused by a local heating of the sample. Among the three investigated sample shapes, the drilled single domain displays a trapped induction which is weakly affected by the local heating while displaying a high trapped field. Finally, a simple numerical modelling of the heat flux spreading into a drilled sample is used to suggest some design rules about the hole configuration and their size. © 2005 IOP Publishing Ltd.