ALBERT

Description: For planners, processes of complex spatial transformations today are comparable to uncharted land and an uncertain voyage. Many possible role images overlap and contrast to traditional and established ways of thinking and acting. The focus here is on navigating instead of controlling, about supporting instead of enforcing. Planning lacks tools to think and act when facing uncertainty. This paper proposes role-reflexive planning as an educational and experimental approach to thinking through different potentialities. It offers groundwork from the boundary between planning and transition studies, using role-based ideas as a bridge. It offers an overview about different roles that are relevant to working towards transformations as spatial planners. It develops an account of role-reflexive planning that connects between contexts, actions and back to individual modes of behaviour in planning processes. As a basis, this paper condenses experiences of a role-playing pilot workshop and discussions about potential elements of a transition towards 'post-growth planning'. It outlines how role-playing challenges the individual roles of actors beyond the game situations themselves. Conceptual ideas foster a renewed role-based debate on thinking and acting in the face of uncertainty and ways to navigate through the stormy waters of transformation.

Description: ZusammenfassungDie Leerstandsquote ist ein zentraler Indikator der Wohnungsmarktbeobachtung und wird zur Beurteilung von Wohnungsmarktsituationen herangezogen. In Deutschland variieren die Leerstandsquoten der Städte und Gemeinden stark, was unter anderem durch die Merkmale unterschiedlicher Bestände des Wohnungsmarktes bedingt ist, wie z. B. Baualter oder Wohnungsgröße. Für die Abbildung der gesamten dabei auftretenden Variationen wurde bislang kaum eine systematische Analyse von Leerstandsstrukturen erstellt. Vor diesem Hintergrund wird im vorliegenden Beitrag das Ziel verfolgt, die differenzierte Struktur des Leerstandes für deutsche Gemeinden aufzuzeigen. Dafür werden Daten der Gebäude- und Wohnungszählung (GWZ)2011 verwendet, die flächendeckend und kleinräumig Leerstand differenziert nach 50 Ausprägungen erfassten. Diese Ausprägungen werden auf ihren Erklärungsgehalt für die Leerstandsquote geprüft. Darauf aufbauend werden räumliche Leerstandstypen gebildet, die die Leerstandsstruktur in Deutschland auf kommunaler Ebene wiedergeben. Der Beitrag macht deutlich, dass in erster Linie Merkmale des Baualters sowie der Gebäude- und Wohnungsgröße die Leerstandsquote erklären, sich aber deutliche Unterschiede zwischen Ost- und Westdeutschland etwa in Bezug auf die Eigentumsform zeigen. Die kommunale Typisierung schärft und verräumlicht diese Erkenntnisse und unterstreicht, dass neben Ost-West- sowie Stadt-Land-Unterschieden auch Typen der Leerstandsstruktur existieren, die im gesamten Bundesgebiet vorkommen. Der Beitrag stellt eine Analyse vor, die mithilfe der entwickelten Typologie helfen soll, Leerstand zu differenzieren, zu bewerten und einzuordnen. Daraus lassen sich dann nicht nur Aussagen zum Einsatz von städtebaulichen und wohnungspolitischen Maßnahmen ableiten, auch die Raum- bzw. Wohnungsmarktbeobachtung auf Bundesebene lässt sich weiter qualifizieren.

Description: Coastal and marine areas represent an increasingly important and relevant action space for spatial planning. However, to a large extent marine (or maritime) spatial planning has emerged separately from terrestrial spatial planning, constituting its own epistemic community. In particular, previous studies indicate that Marine Spatial Planning often follows an expert-driven resource management rationale focused on sea-use regulation. This paper examines practices of Marine Spatial Planning and Integrated Coastal Zone Management at the German North Sea coast. The paper focuses in particular on the engagement of spatial planners with these practices and their perception of their role therein. We seek to understand what form spatial planning at the coast and at sea currently takes and how this might develop in the future in response to current and anticipated policy developments. We argue for the necessity of a communicative, cross-sectoral approach to spatial planning at sea, providing a spatial vision for the future that extends from the Exclusive Economic Zone to encompass both the coastal waters of the federal states and the land-sea interface in a substantive manner.

Description: ZusammenfassungWenn man räumliche Planung verstehen möchte, ist es wichtig, die Selbstverständnisse und beruflichen Identitäten der handelnden Akteure zu kennen. Bislang ist dazu in Deutschland kaum geforscht worden. Der Beitrag basiert auf einem Forschungsansatz, der berufliche Identitäten als Ergebnis des Zusammenwirkens gesellschaftlicher Anforderungen und eigener Handlungen in den Fokus rückt. Es wird gezeigt, welche Rollen Regionalplanern von außen zugewiesen werden und wie Regionalplaner selber ihre beruflichen Identitäten wahrnehmen und daran arbeiten. Anhand von Dokumentenanalysen und autobiographisch-narrativen Interviews wird untersucht, welche Subjektpositionen in schriftlichen Diskursfragmenten, die sich auf Regionalplanung und Windenergie beziehen, zu finden sind, welchen äußeren Anforderungen sich Regionalplaner ausgesetzt sehen, welche Ansprüche sie selber an ihr berufliches Handeln stellen, welche Techniken des Selbst sie praktizieren und welche diesbezüglichen Spannungen sie wahrnehmen. Die Ergebnisse knüpfen an internationale Untersuchungen über planerische Identitäten an und können als Angebot an die Praxis gelesen werden, individuelle Praktiken und Spielräume der Identitätsarbeit zu reflektieren.

Description: The process of gentrification in the peri-urban districts of French cities has scarcely been touched upon in recent research, which has hitherto seen the phenomenon as typically associated with core urban areas. The tendency has been to view the periphery through the lens of the social crisis of the banlieues. In contrast, the present article focuses on gentrification in the metropolitan area of Nancy (Grand Est region) as a development that also plays a role in municipalities around major cities and especially around regional metropolitan centres. Starting with a survey of current research approaches, the article first pinpoints some gaps and methodological imbalances that need to be tackled, before embarking on the case study of peri-urban Nancy. Statistical data and empirical surveys in the form of qualitative interviews indicate how Nancy's peri-urban districts have developed a logic of separation, exclusion and social decoupling – typical features of gentrification – particularly in connection with the construction of new single-family houses as a supplement to existing residential stock. Key questions here concern individual motives for choosing a particular residential location, and the creeping "segregation from above" that accompanies this process. The image of France's peri-urban spaces that arises from this study stands in explicit contrast to the received, markedly negative connotations of the "urban periphery".

Description: Kurzfassung Als Inbegriff des Urbanen erlebt der Stadttypus „Metropole“ seit nunmehr zwei bis drei Jahrzehnten eine fortwährende Konjunktur. Mit einer quantitativen Textanalyse der führenden deutschen überregionalen Printmedien werden die Muster des Sprachgebrauchs für den Begriff Metropole und für die beiden größten Städte Berlin und Hamburg herausgearbeitet. Sie werden als Bedeutungszuschreibungen in massenmedialen Kommunikationsprozessen interpretiert und erzeugen ein semantisches Feld, dessen Analyse Aussagen über Größe, Funktion und Symbolik einer Metropole auch im Zeitverlauf zulässt.

Description: Kurzfassung In der Wissensgesellschaft differenziert sich die Raumstruktur in eine Vielzahl von Raumtypen aus; auch der ländliche Raum stellt sich nicht als ein homogenes Gebilde dar. Die „reflexive Moderne“ der Wissensgesellschaft scheint keinen Determinismus in Bezug auf räumliche Entwicklungsunterschiede zu kennen, sondern den Menschen Spielräume für eine Aneignung und Inwertsetzung unterschiedlicher Räume zu öffnen. Diese erweiterten Gestaltungsmöglichkeiten könnten auch in dem bisher als zurückbleibend eingeschätzten ländlichen Raum neue Entwicklungsimpulse auslösen. Zur Überprüfung dieser Annahmen wird in diesem Beitrag aus einem theoretischen Blickwinkel gefragt, wie die kulturellen, sozialen, wirtschaftlichen und technischen Dimensionen der Wissensgesellschaft mit dem Raum interagieren und welche Herausforderungen sowie Strategien sich daraus für den ländlichen Raum ergeben können. Die empirische Überprüfung stützt sich auf Ergebnisse einer Arbeitsgruppe der Akademie für Raumforschung und Landesplanung (ARL) zur Wissensgesellschaft in ländlich-peripheren Regionen sowie ergänzend auf weitere quantitative und qualitative Studien zum Thema. Aus diesen Untersuchungen geht hervor, dass sich in vielen, wenn auch nicht in allen Regionen neuartige wissensgesellschaftlich geprägte Lebenskulturen in Abgrenzung zu städtischen Lebensmodellen entfalten. Sie belegen eine neue Attraktivität des ländlichen Raumes, die sich auch in wissensbasierten Projekten der Kultur, Bildung, Daseinsvorsorge und des unternehmerischen Engagements niederschlägt. Beispiele finden sich besonders in wirtschaftlich wettbewerbsfähigen Regionen, aber auch in den peripheren Regionen, die unter Bevölkerungsverlusten leiden. Es bleibt weiteren Untersuchungen vorbehalten zu eruieren, wie diese Ansätze von öffentlichen Akteuren strategisch genutzt werden können, um den wissensgesellschaftlichen Wandel voranzutreiben. Sie können helfen, in erfolgreichen Regionen gesammelte Erfahrungen für eine politische und planerische Unterstützung peripherer Städte und Regionen im ländlichen Raum zu nutzen.

Description: This paper explores different functions of co-presence for collaborative knowledge creation in the context of seed accelerator programs. Seed accelerators offer programs of three to six months to enhance the growth of early-stage start-ups through various means of training and organizational development. In this paper, seed accelerator programs are analysed as orchestrated sequences of different types of physical co-presence. By drawing on qualitative case study data, the paper identifies eight different ways in which physical co-presence is used during seed accelerator programs. Through these eight types of co-presence, the analysis reveals that physical co-presence unfolds unique social dynamics that are utilized in a carefully designed combination of presence and absence. It is shown that physical co-presence is a means to enact relational distance and to bridge this distance for the benefit of the entrepreneurial process. Therefore, this paper adds a new perspective on how co-presence is used to facilitate the generation of value through collaborative knowledge creation.

Description: The article presents problem of non-uniform foundation of structures in weak wet subsoil. The problem is illustrated with the case study of two-chamber-reinforced concrete water tank constructed in 1920s of 20th century, which cracked during construction. Under part of foundation, where the peat was found, the concrete piles were introduced.The results of five-year measurement of crack widths with crack gauges and geodesic measurements of vertical displacement of tank were presented. These results indicate that the tank is not stable and part of broken tank supported on piles is movable.On the basis of the presented data, the general conclusions concerning the non-uniform founding of tanks are formulated.

Description: The inclusions of geosynthetic materials (fibers, geomembranes and geotextiles) is a new improvement technique that ensures uniformity in the soil during construction. The use of tension resisting discreet inclusions like polypropylene fibers has attracted a significant amount of attention these past years in the improvement of soil performance in a cost-efficient manner. A series of direct shear box tests were conducted on unreinforced and reinforced Chlef sand with different contents of fibers (0, 0.25, 0.5 and0.75%) in order to study the mechanical behavior of sand reinforced with polypropylene fibers. Samples were prepared at three different relative densities 30%, 50% and 80% representing loose, medium dense and dense states,respectively, and performed at normal stresses of 50, 100 and 200 kPa. The experimental results show that the mechanical characteristics are improved with the addition of polypropylene fibers. The inclusion of randomly distributed fibers has a significant effect on the shear strength and dilation of sandy soil. The increase in strength is a function of fiber content, where it has been shown that the mechanical characteristics improve with the increase in fiber content up to 0.75%, this improvement is more significant at a higher normal stress and relative density.

Description: The present article examines the problem related to the axisymmetric torsion of an elastic layer by a circular rigid disc at the symmetry plane. The layer is sandwiched between two similar elastic half-spaces with two penny-shaped cracks symmetrically located at the interfaces between the two bonded dissimilar media. The mixed boundary-value problem is transformed, by means of the Hankel integral transformation, to dual integral equations, that are reduced, to a Fredholm integral equation of the second kind. The numerical methods are used to convert the resulting system to a system of infinite algebraic equations. Some physical quantities such as the stress intensity factor and the moment are calculated and presented numerically according to some relevant parameters. The numerical results show that the discontinuities around the crack and the inclusion cause a large increase in the stresses that decay with distance from the disc-loaded. Furthermore, the dependence of the stress intensity factor on the disc size, the distance between the crack and the disc, and the shear parameter is also observerd.

Description: A characteristic feature of soil-steel structures is that, unlike in typical bridges, the backfill and the carriageway pavement with its foundation play a major role in bearing loads. In the soil-steel structure model, one can distinguish two structural subsystems: the shell made of corrugated plates and the backfill with the pavement layers. The interactions between the subsystems are modelled as interfacial interactions, that is, forces normal and tangent to the surface of the shell. This is a static condition of the consistency of mutual interactions between the surrounding earth and the shell, considering that slip can arise at the interface between the subsystems. This paper presents an algorithm for determining the internal forces in the shell on the basis of the unit strains in the corrugated plates, and subsequently, the interfacial interactions. The effects of loads arising during the construction of a soil-steel bridge when, for example, construction machines drive over the structure, are taken into account in the analysis of the internal forces in the shell and in the surrounding earth. During construction, the forces in the shell are usually many times greater than the ones generated by service loads. Thus, the analytical results presented in this paper provide the basis for predicting the behaviour of the soil medium under operational loads.

Description: This article presents the results of laboratory tests on soft, normally consolidated soils from the Vistula Marshlands. Samples of high-plasticity organic soils (muds) taken from 3.2–4.0 m and 9.5–10.0 m depth, as well as peat deposit at 14.0 m, are analysed. Presented case study confirms the applicability of the Norwegian Institute of Technology (NTH) method based on Cone Penetration Tests (CPTU) and allows for a conservative estimation of effective friction angle for muds. The plastification angle equal to 14.5° for organic silt, applied in the modified NTH method, fits well the triaxial test (TX) results. Moreover, the dilative-contractive behaviour according to the CPTU soil classification based on the Robertson’s proposal from 2016 corresponds well with volumetric changes observed in the consolidated drained triaxial compression tests. The internal friction angles of the Vistula Marshlands’ muds and peats are lower in comparison with the database of similar soft soils.

Description: Based on the response of small-scale model square footing, the present paper shows the results of an experimental bearing capacity of eccentrically loaded square footing, near a slope sand bed. To reach this aim, a steel model square footing of (150 mm × 150 mm) and a varied sand relative density of 30%, 50% and 70% are used. The bearing capacity-settlement relationship of footing located at the edge of a slope and the effect of various parameters such as eccentricity (e) and dimensions report (b/B) were studied. Test results indicate that ultimate bearing capacity decreases with increasing load eccentricity to the core boundary of footing and that as far as the footing is distant from the crest, the bearing capacity increases. Furthermore, the results also prove that there is a clear proportional relation between relative densities –bearing capacity. The model test provides qualitative information on parameters influencing the bearing capacity of square footing. These tests can be used to check the bearing capacity estimated by the conventional methods.

Description: Tunnel construction below or adjacent to piles will affect the performance and eventually the stability of piles due to ground deformation resulting in the movement of piles and changes in the axial force distribution along the piles. A three dimensional finite element analysis using PLAXIS 3D (2013) was performed to study the behaviour of a single pile and 3 x 3 piles group during the advancement of shield tunnelling in ground. The 10-node tetrahedral elements were used to model both the soil and the tunnel lining. The Hardening Soil (HS) model was used to simulate the soil structure interaction at the tunnel-soil interface. An isotropic elastic model was used for the pile, piles cap, tunnel lining and tunnel boring machine shield (TBM). Several parametric studies were attempted including the longitudinal, lateral, and vertical tunnel location relative to pile embedded in different types of soil (clay or sand). The results showed that the pile head settlement increases during the tunnelling advancement in larger values than that for ground surface settlement. A zone of influence was determined in the range of twice the tunnel diameter in the longitudinal direction (forward and backward of the pile), and transverse direction (left and right of the tunnel centreline). If the tunnel boring is kept off this zone then there is no fear of pile collapse.

Description: The paper evaluates the effectiveness of reinforcing a damaged earth structure with making counterfort drains in its slope. The system of counterfort drains changed the soil properties significantly over a long-term use. The evaluation was based on many years of field and laboratory tests and stability analysis. The field tests concerned the observation of N WST probing resistance change, and the laboratory tests concerned the change in soil consistency and water content. The paper presents the results of tests that were conducted over 13 years.

Description: The article describes a computer analysis of the pull-out test used to calculate the force needed to pull out a rock fragment and determine the shape of this broken fragment. The analyzed material is sandstone and porphyry. The analysis included the first approach to using own subroutine in the Simulia Abaqus system, that is, which task is undertaken to accurately determine the crack path of the Finite Element Method model. The work also contains a description of laboratory tests and analytical considerations.

Description: At present, the suspended monorail systems constitute a very common means of transportation in the Polish hard coal mines. The main advantages of the suspended monorail include the independence of the route from the working floor surface irregularities and the possibility to transport cargo of significant mass and size.The masses and dimensions of machines and devices transported via monorail have increased considerably in recent times. This particularly concerns the transport of longwall system elements. In Poland, the maximum speed of suspended monorail travel is 2 m/s. Due to the fact that preparations are currently underway to increase the maximum speed above 2 m/s, it is necessary to inspect what influence it will have on work safety and mining support stability.Current operational experience and tests have shown that dynamic loads induced by the suspended monorail transportation have a significant influence on the roadway support stability, working protection durability and on the monorail operators. This is particularly true during the emergency braking of a suspended monorail by means of a braking trolley, where the overloads reach 3g.Bench tests of the selected steel arch and rock bolt support elements utilised in the Polish hard coal mines were conducted in order to determine the resistance of steel arch and rock bolt supports to static and dynamic loads.The article presents the results of the tests conducted on a steel arch support in the form of the sliding joints of an ŁP/V29 yielding roadway support, which is commonly employed in the Polish hard coal mines. Tests of elements of the threaded bolts with trapezoidal threads over the entire rod length were conducted as well.The conducted strength tests of steel arch and rock bolt support elements under static and dynamic loading have shown that dynamic loading has decisive influence on the support’s retaining of its stability. Support element stability decreases along with the increase of the impact velocity. This concerns both the steel arch support and the rock bolt support.

Description: We show that the global non-linear stability threshold for convection in a double-diffusive couple-stress fluid saturating a porous medium is exactly the same as the linear instability boundary. The optimal result is important because it shows that linearized instability theory has captured completely the physics of the onset of convection. It is also found that couple-stress fluid saturating a porous medium is thermally more stable than the ordinary viscous fluid, and the effects of couple-stress parameter (F ) , solute gradient ( S f ) and Brinkman number ( D a ) on the onset of convection is also analyzed.

Description: This paper discusses the pull-out laboratory tests and the monitoring of expansion-shell bolts with a length of 1.82 m. The bolts comprised the KE-3W expansion shell, a rod with a diameter of 0.0183 m and a profiled, circular plate with a diameter of 0.14 m, and a gauge of 0.006 m. The bolts were installed in a concrete block with a compressive strength of 75 MPa. The tests were conducted on a state-of-the-art test stand owned by the Department of Underground Mining of the AGH University of Science and Technology. The test stand can be used to test roof bolts on a geometric scale of 1:1 under static and rapidly varying loads. Also, the stand is suitable for testing rods measuring 5.5 m in length. The stand has a special feature of providing the ongoing monitoring of bolt load, displacement and deformation. The primary aim of the study was to compare the results recorded by two different measurement systems with the innovative Self-Excited Acoustic System (SAS) for measuring stress variations in roof bolts. In order to use the SAS, a special handle equipped with an accelerometer and exciter mounted to the nut or the upset end of the rod was designed at the Faculties of Mining and Geoengineering and Mechanical Engineering and Robotics of the AGH University of Science and Technology. The SAS can be used for nondestructive evaluation of performance of bolts around mining workings and in tunnels. Through laboratory calibration tests, roof bolt loads can be assessed using the in-situ non-destructive method.

Description: This paper presents an experimental investigation of the discharge flow pressure in the vertical silo and the hopper due to the use of insert (top cone with trunk cone bottom). Using the Insert inside the silos is one of the proposed solutions to avoid the problems of having funnel flow pattern, which has a significant effect on the distribution of flow pressure exerted on the silo wall and the hopper. The experiments were performed on a metal cylinder prototype; corn was used as a granular material, and the wall and hopper pressure distribution was measured by a special pressure transducer. The experiments revealed an important result in the flow pressure due to the change in the location of the insert. The experiments were conducted in Damascus University laboratories.

Description: Water seepage is one of the most important features of embankment dams. To prevent and reduce seepage, it is necessary to seal the dam. Plastic concrete cutoff walls are one of the most efficient methods in waterproofing the foundation of embankment dams on permeable alluvial substrates. Sufficient resistance to loads, low permeability to maintain dam sealing, high ductility compatible with the foundation and deformation under load without cracking are the main requirements in plastic concrete cutoff walls. In this paper, the construction and implementation of the cutoff wall of Karkheh Dam, which is one the world’s largest water sealing projects, was studied. In addition, a numerical model using Seep-3D software was developed to evaluate the efficiency of the cut-off wall to decrease the seepage over the dam’s foundation. The numerical results validated by instrumentation statistics resulted from 17-years dam operation. According to the results, after the drainage of the reservoir, the cutoff wall optimally reduced the hydraulic gradient by 0.08 from 2.35 and the water leakage by 3.1 m/s from 18.3 m/s.

Description: The main issue of the paper is the estimation of soil hydraulic permeability based on the DMT test. DMTA, DMTC and SASK methods performed in the Nielisz dam, Stegny and the SGGW Campus of the Warsaw University of Life Sciences sites are described. The article presents the implementation of the dilatometer Marchetti test (DMT) in the determination of soil fraction and effects of its occurrence in the subsoil, tested in the Nielisz dam located in the Wieprz river valley in the Lublin province, and in various sites in Warsaw (Stegny site and SGGW Campus of the Warsaw University of Life Sciences). In order to acquire the needed data, the flat dilatometer test (DMT) method was used. A direct and indirect pressure methodology of interpreting soil swelling was characterized in the article. The paper shows the possibilities of determining sand, silt and clay soil fractions based on po and p1 pressures from dilatometer tests (DMT) and the effective (σ’vo) and total (σvo) vertical in situ overburden stress. Additionally, the main advantage of this paper is the proposal of use of a new chart to determine hydraulic permeability and soil fraction, based on DMT tests.

Description: This article deals with the vibrations of a nonprismatic thin-walled beam with an open cross section and any geometrical parameters. The thin-walled beam model presented in this article was described using the membrane shell theory, whilst the equations were derived based on the Vlasov theory assumptions. The model is a generalisation of the model presented by Wilde (1968) in ‘The torsion of thin-walled bars with variable cross-section’, Archives of Mechanics, 4, 20, pp. 431–443. The Hamilton principle was used to derive equations describing the vibrations of the beam. The equations were derived relative to an arbitrary rectilinear reference axis, taking into account the curving of the beam axis and the axis formed by the shear centres of the beam cross sections. In most works known to the present authors, the equations describing the nonprismatic thin-walled beam vibration problem do not take into account the effects stemming from the curving (the inclination of the walls of the thin-walledcross section towards to the beam axis) of the analysed systems. The recurrence algorithm described in Lewanowicz’s work (1976) ‘Construction of a recurrence relation of the lowest order for coefficients of the Gegenbauer series’, Applicationes Mathematicae, XV(3), pp. 345–396, was used to solve the derived equations with variable coefficients. The obtained solutions of the equations have the form of series relative to Legendre polynomials. A numerical example dealing with the free vibrations of the beam was solved to verify the model and the effectiveness of the presented solution method. The results were compared with the results yielded by finite elements method (FEM).

Description: This paper presents the results obtained from an experimental programme and numerical investigations conducted on model tests of strip footing resting on reinforced and unreinforced sand slopes. The study focused on the determination of ultimate bearing capacity of strip footing subjected to eccentric load located either towards or opposite to the slope facing. Strip footing models were tested under different eccentricities of vertical load. The obtained results from tests conducted on unreinforced sand slope showed that the increase in eccentricity of applied load towards the slope facing decreases the ultimate bearing capacity of footing. Predictions of the ultimate bearing capacity obtained by the effective width rule are in good agreement with those proposed from the consideration of total width of footing subjected to eccentric load. The ultimate bearing capacity of an eccentrically loaded footing on a reinforced sand slope can be derived from that of axially loaded footing resting on horizontal sand ground when adopting the effective width rule and the coefficient of reduction due to the slope. When increasing the distance between the footing border to the slope crest, for unreinforced and reinforced ground slope by geogrids, the ultimate bearing capacity of footing is no more affected by the slope ground.

Description: In urban areas, the control of ground surface settlement is an important issue during shield tunnel-boring machine (TBM) tunneling. These ground movements are affected by many machine control parameters. In this article, a finite difference (FD) model is developed using Itasca FLAC-3D to numerically simulate the whole process of shield TBM tunneling. The model simulates important components of the mechanized excavation process including slurry pressure on the excavation face, shield conicity, installation of segmental lining, grout injection in the annular void, and grout consolidation. The analysis results from the proposed method are compared and discussed in terms of ground movements (both vertical and horizontal) with field measurements data. The results reveal that the proposed 3D simulation is sufficient and can reasonably reproduce all the operations achieved by the TBM. In fact, the results show that the TBM parameters can be controlled to have acceptable levels of surface settlement. In particular, it seems that moderate face pressure can reduce ground movement significantly and, most importantly, can prevent the occurrence of face-expected instability when the shield crosses very weak soil layers. The shield conicity has also an important effect on ground surface settlement, which can be partly compensated by the grout pressure during tail grouting. Finally, the injection pressure at the rear of the shield significantly reduces the vertical displacements at the crown of the tunnel and, therefore, reduces the settlement at the ground surface.

Description: This paper presents an enhanced constitutive model integrating deviatoric hardening with a modified yield surface for overconsolidated clayey soils in a general framework of Cam-clay type models. Its performance was assessed with the simulation of drained and undrained triaxial tests on three clays at different consolidation states in comparison to two critical state models. The proposed model satisfactorily estimates the shear resistance, while capturing the smooth nonlinearity of the soil response.Shear triaxial tests at constant mean pressure were performed on an overconsolidated marl to study the shear response. Their simulation attests the importance of deviatoric hardening integration.

Description: The present study investigates the onset of penetrative convection in- duced by selective absorption of radiation in a magnetic nanofluid saturated porous medium. The influence of Brownian motion, thermophoresis, and magnetophoresis on magnetic nanofluid treatment is taken into consideration. The Darcy’s model is selected for the porous medium. We conduct a linear stability analysis to examine the onset of instability and evaluate the results for two different configurations, namely, when the layer is heated from below and when the layer is heated from above. The numerical investigations are carried out by applying the Chebyshev pseudospectral method. The effect of the porosity parameter E, parameter Y (represents the ratio of internal heating to boundary heating), Lewis number Le, concentration Rayleigh number Rn, Langevin parameter αL, width of nanofluid layer d, diffusivity ratio η, and modified diffusivity ratio NA is examined at the onset of convection. The results indicate that the convection commences easily with an increase in the value of Y, Le, and NA but opposite in the case with a decrease in the value of E, αL, η and d for both the two configurations. The parameter Rn advances the onset of convection when the layer is heated from below, while delays the onset of convection when the layer is heated from above.

Description: This article presents the results of tests on the energy properties of sedimentary rocks in the Upper Silesian Coal Basin. The rocks were tested both in an air-dry state and in a water saturation state. Samples of sedimentary rocks were collected from boreholes drilled in the underground workings of coal mines located within the area of the city of Jastrzębie, in the areas of the Chwałowice Trough and Rybnik Trough (south-western part of the Upper Silesian Coal Basin) and in the Main Trough. Influence of saturation condition on the values of the tested energy parameters was observed. The values of elastic energy and dissipated energy obtained for the samples tested in water saturation were lower compared to the values obtained for samples tested in air-dry state. As observed, an increase in the values of the given types of specific energy corresponds to an increase in the uniaxial compression strength in air-dry state and in water saturation state. Results of the tests are original and they can be applied while analysing the possibility of the occurrence of some dynamic phenomena and hazards in mine workings in Carboniferous rock mass in the Upper Silesian Coal Basin, caused by mining operations.

Description: Different types of foundations are used in steel, above-ground cylindrical storage tanks for liquids. If a sand-gravel foundation is used under the entire bottom of the tank or only in the central part of the tank, settlement can be expected, and it increases after many years of operation. The paper presents the typical kinds and types of soil settlements under the bottoms of the tanks, in which different types of foundations were used. Numerical analyses of the effect of the soil settlement on the state of deformations and stresses in steel sheets of the bottom under one of the real tanks, in which different types of foundations and different cases of settlement were assumed. The results of numerical analyses indicated the possibility of evaluating the state of the soil settlement and bottom sheet deformations on the basis of simple measurements of deformations of the lower part of the tank cylinder. These measurements can be very useful in assessing the possible risk of failure of the tank bottom during each period of its operation, as measurements of settlement of the bottom of a filled tank are not feasible in practice. It has been proposed that in each steel tank, the deformation of the cylinder’s sheets should be measured even before the beginning of exploitation, and that in subsequent periodical measurements, the influence of the soil settlement under the tank on the state of the cylinder deformation and bottom’s strain should be assessed more accurately.

Description: Sands reinforced by hydraulic binders (cement) have constituted in recent decades a major asset for the expansion of several areas of engineering. The mechanical behavior of sand-cement mixtures has undergone some controversies studied on the Chlef sand. In this paper, we present an experimental study to investigate the mechanical behavior of a sandy soil reinforced by a hydraulic binder (cement), using the direct shear apparatus emphasizing on the shear strength characteristics and the vertical deformation variation of cemented reinforced sand. The parameters used in this study are mainly: relative density (Dr = 80%), normal stress (σn = 100, 200, 400 kPa), water content (3, 7 and 10%), cement content (2.5, 5, 7.5 and 10 %) and cure time (7, 14 and 28 days). The experimental results show that the mechanical characteristics in terms of internal cohesion (C) and internal frication angle (φ) give a better mechanical performance with the binder inclusion, and the cure conditions play an effective role on the improvement of the shear strength. This result also showed that 10% of the cement content gave us a maximum value of shear strength and an optimal influence on the mechanical characteristics. The addition of cement not only improves the shear strength of soil, but also provides diversity in the resistance against the deformations imposed load, which can be established by a dilatant character.

Description: Extractive industries often use explosives to destroy rocks, and productivity requirements tend to increase the charges of the explosives. The blasts induce vibrations, which result in a potential damage of the surrounding structures. Therefore, the prediction of vibrations should be described with accuracy, in order to ensure the safety of engineered structures. However, the prediction of vibrations’ levels remain a complicated issue, because it involves numerous parameters correlated to the quarry site. In this paper, statistical analysis based on the peak particle velocity (PPV) and the attenuation law has been carried out to assess the safety charges (Q) for different distances (R) between the blast and the considered structure to secure. Moreover, the experimental investigations were conducted on the quarry site of “Sococim”, which is located on the south coast of Senegal. To ensure the safety of the “Conveyor belt” and “Panel 1 (Upper exploitation level)” sites, the PPV should be less than 10 mm/s. In fact, the attenuation model has been used to assess the safe charge weights of the explosive (Q) to be used at the “Conveyor belt” site and at the “Panel 1 (Upper exploitation level)” site. Therefore, the safe charge weights per delay (Q) were respectively 116 kg and 13.75 kg.

Description: This work presents the geotechnical problems occurring in the interaction area between road embankments and the bridge structures in case a subsoil characterised by complex and complicated geological and engineering conditions. These significant problems that occur during the design, performance and exploitation of the abutment structures, are illustrated on the example of engineering practice in Lower Silesia, concerning a road embankment that constitutes access to the bridge. The results of numerical analyses concerning the process of consolidation of low-strength soils and their impact on the settlements of road embankment indicate the need to carry out such analyses also in the cases, when the soft soils occur in the direct geotechnical layer under the designed embankment. The Authors included in this article a discussion regarding other effective actions and solutions that can be used in the design and performance phase, leading to the elimination or reduction of problems concerning the connection of engineering structures with road embankments, which have been recurring for years, ultimately resulting in the improvement of quality, comfort and safety of road exploitation.

Description: Understanding the quality of intact rock is one of the most important parts of any engineering projects in the field of rock mechanics. The expression of correlations between the engineering properties of intact rock has always been the scope of experimental research, driven by the need to depict the actual behaviour of rock and to calculate most accurately the design parameters. To determine the behaviour of intact rock, the value of important mechanical parameters such as Young’s modulus (E), Poisson’s ratio (ν) and the strength of rock (σcd) was calculated. Recently, for modelling the behaviour of intact rock, the crack initiation stress (σci) is another important parameter, together with the strain (σ). The ratio of Young’s modulus and the strength of rock is the modulus ratio (MR), which can be used for calculations. These parameters are extensively used in rock engineering when the deformation of different structural elements of underground storage, caverns, tunnels or mining opening must be computed. The objective of this paper is to investigate the relationship between these parameters for Hungarian granitic rock samples. To achieve this goal, the modulus ratio (MR = E/σc) of 50 granitic rocks collected from Bátaapáti radioactive waste repository was examined. Fifty high-precision uniaxial compressive tests were conducted on strong (σc 〉100 MPa) rock samples, exhibiting the wide range of elastic modulus (E = 57.425–88.937 GPa), uniaxial compressive strength (σc = 133.34–213.04 MPa) and Poisson’s ratio (ν = 0.18–0.32). The observed value (MR = 326–597) and mean value of MR = 439.4 are compared with the results of similar previous researches. Moreover, the statistical analysis for all studied rocks was performed and the relationshipbetween MR and other mechanical parameters such as maximum axial strain $left( {{varepsilon }_{ext{a,},ext{max}}} ight)$for studied rocks was discussed.

Description: Waste material such as used tires is increasing every year, which poses environmental problems. However, such material has been used in several geotechnical applications as alternative lightweight backfill in highway embankments and/or behind retaining walls, providing environmental, economic and technical benefits. These applications require knowledge of engineering properties of soil-tire rubber mixtures. The present study aims to show the possibility of tire rubber usage in sand by evaluating the shear strength and deformability of sand mixed with granulated rubber, in weight percentages between 0 and 50%. The tire rubber content was found to influence the stress-strain and deformation behavior of the mixtures. The shear strength of sand mixed with 10% or 20% tire rubber was higher than that measured for sand only. However, the trend for TRC = 30–50% was different. Samples with a rubber content of 30-50% exhibited a rapid decrease in the stress ratio compared with that of sand. The major principal strain at maximum stress ratio was found to increase with increasing tire rubber content. However, it was observed that the lateral strains (minor and intermediate principal strains) of samples reduced significantly with the addition of tire rubber to the sand.

Description: 〈div data-abstract-type="normal"〉〈p〉The effectiveness of streak modes in controlling the oblique-type breakdown in a supersonic boundary-layer at Mach 2.0 is investigated using direct numerical simulations. Investigations in the literature have shown the effectiveness of streak modes in delaying the onset of transition dominated by two-dimensional waves, but in oblique breakdown, three-dimensional waves and a strong streak mode dominate the transition process. Paredes 〈span〉et al.〈/span〉 (〈span〉J. Fluid Mech.〈/span〉, vol. 831, 2017, pp. 524–553) discussed the possible stabilization of supersonic boundary layers by optimally growing streaks using parabolized stability equations. However, no study has as yet been reported regarding direct nonlinear control of oblique breakdown. This study deals with the effects of large-amplitude decaying streak modes generated by a blowing–suction strip at the wall to control full breakdown in a reference case. Modes with four to five times the fundamental wavenumber are found to be beneficial for controlling the transition. In the first region after the control-mode forcing, the beneficial mean-flow distortion (MFD), generated by inducing the control mode, is solely responsible for hampering the growth of the fundamental-mode. On the whole, the MFD and the three-dimensional part of the control contribute equally towards controlling the oblique breakdown. The results show significant suppression of transition, and substantial improvements have been observed in the levels of the skin-friction coefficient and wall-temperature in comparison to the uncontrolled case. Moreover, refreshing the control using an additional downstream control strip increases the gain. However, the forcing amplitude must be carefully chosen in order not to introduce a generalized inflection point in the spanwise averaged mean flow invoking enhanced disturbance growth.〈/p〉〈/div〉

Description: 〈div data-abstract-type="normal"〉〈p〉The problem of interaction between disturbances and shock waves was solved by a theoretical approach called linear interaction analysis in the mid-twentieth century. More recently, great progress has been made in analysing shock–turbulence interactions by direct numerical simulation. However, an unsolved theoretical problem remains: What happens when no acoustic waves are stimulated behind the shock wave? The concept of a damped wave is introduced, which is a type of excited plane wave. Based on this, the dispersion and amplitude relationships between any incident plane wave and resulting stimulated waves are constructed analytically, systematically and comprehensively. The physical essence of damped waves and the existence of critical angles are clarified. It is demonstrated that a damped wave is a complex number space solution to the acoustic dispersion relationship under certain conditions. It acts as a bridge connecting fast and slow acoustic waves at the position where the 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190628130807083-0762:S0022112019004385:S0022112019004385_inline1.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 component of the group velocity is zero. There are two critical angles that can excite fast and slow acoustic waves, which determine the conditions that stimulate a damped wave. Our results show good agreement with theoretical and simulation results. The contribution of each excited wave to the transmission coefficient is evaluated, the distribution of the transmission coefficient is analysed and application to an engineering wedge model is performed.〈/p〉〈/div〉

Description: 〈div data-abstract-type="normal"〉〈p〉The onset of thermal convection in a rapidly rotating spherical shell is studied by linear stability analysis based on the fully compressible Navier–Stokes equations. Compressibility is quantified by the number of density scale heights 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190628092857401-0357:S0022112019004361:S0022112019004361_inline1.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉, which measures the intensity of density stratification of the motionless, polytropic base state. The nearly adiabatic flow with polytropic index 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190628092857401-0357:S0022112019004361:S0022112019004361_inline2.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 is considered, where 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190628092857401-0357:S0022112019004361:S0022112019004361_inline3.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 is the adiabatic polytropic index. By investigating the stability of the base state with respect to the disturbance of specified wavenumber, the instability process is found to be sensitive to the Prandtl number 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190628092857401-0357:S0022112019004361:S0022112019004361_inline4.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 and to 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190628092857401-0357:S0022112019004361:S0022112019004361_inline5.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉. For large 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190628092857401-0357:S0022112019004361:S0022112019004361_inline6.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 and small 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190628092857401-0357:S0022112019004361:S0022112019004361_inline7.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉, the quasi-geostrophic columnar mode loses stability first; while for relatively small 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190628092857401-0357:S0022112019004361:S0022112019004361_inline8.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 a new quasi-geostrophic compressible mode is identified, which becomes unstable first under strong density stratification. The inertial mode can also occur first for relatively small 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190628092857401-0357:S0022112019004361:S0022112019004361_inline9.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 and a certain intensity of density stratification in the parameter range considered. Although the Rayleigh numbers 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190628092857401-0357:S0022112019004361:S0022112019004361_inline10.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 for the onsets of the quasi-geostrophic compressible mode and columnar mode are different by several orders of magnitude, we find that they follow very similar scaling laws with the Taylor number. The critical 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190628092857401-0357:S0022112019004361:S0022112019004361_inline11.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 for convection onset is found to be always positive, in contrast with previous results based on the widely used anelastic model that convection can occur at negative 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190628092857401-0357:S0022112019004361:S0022112019004361_inline12.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉. By evaluating the relative magnitude of the time derivative of density perturbation in the continuity equation, we show that the anelastic approximation in the present system cannot be applied in the small-〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190628092857401-0357:S0022112019004361:S0022112019004361_inline13.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 and large-〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190628092857401-0357:S0022112019004361:S0022112019004361_inline14.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 regime.〈/p〉〈/div〉

Description: 〈div data-abstract-type="normal"〉〈p〉We study the breakup of confined fluid threads at low flow rates to understand instability mechanisms. To determine the critical conditions between the earlier quasi-stable necking stage and the later unstable collapse stage, simulations and experiments are designed to operate at an extremely low flow rate. The critical mean radii at the neck centres are identified by the stop-flow method for elementary microfluidic configurations. Two distinct origins of capillary instabilities are revealed for different confinement situations. One is the gradient of capillary pressure induced by the confinements of geometry and external flow, whereas the other is the competition between the capillary pressure and internal pressure determined by the confinements.〈/p〉〈/div〉

Description: 〈div data-abstract-type="normal"〉〈p〉We report on a combined experimental and numerical study of convective heat transfer along ratchet surfaces in vertical natural convection (VC). Due to the asymmetry of the convection system caused by the asymmetric ratchet-like wall roughness, two distinct states exist, with markedly different orientations of the large-scale circulation roll (LSCR) and different heat transport efficiencies. Statistical analysis shows that the heat transport efficiency depends on the strength of the LSCR. When a large-scale wind flows along the ratchets in the direction of their smaller slopes, the convection roll is stronger and the heat transport is larger than the case in which the large-scale wind is directed towards the steeper slope side of the ratchets. Further analysis of the time-averaged temperature profiles indicates that the stronger LSCR in the former case triggers the formation of a secondary vortex inside the roughness cavity, which promotes fluid mixing and results in a higher heat transport efficiency. Remarkably, this result differs from classical Rayleigh–Bénard convection (RBC) with asymmetric ratchets (Jiang 〈span〉et al.〈/span〉, 〈span〉Phys. Rev. Lett.〈/span〉, vol. 120, 2018, 044501), wherein the heat transfer is stronger when the large-scale wind faces the steeper side of the ratchets. We reveal that the reason for the reversed trend for VC as compared to RBC is that the flow is less turbulent in VC at the same 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190628072522266-0966:S0022112019004464:S0022112019004464_inline1.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉. Thus, in VC the heat transport is driven primarily by the coherent LSCR, while in RBC the ejected thermal plumes aided by gravity are the essential carrier of heat. The present work provides opportunities for control of heat transport in engineering and geophysical flows.〈/p〉〈/div〉

Description: 〈div data-abstract-type="normal"〉〈p〉We explore the dynamics of inclined temporal gravity currents using direct numerical simulation, and find that the current creates an environment in which the flux Richardson number 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190627132031350-0413:S0022112019004300:S0022112019004300_inline1.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉, gradient Richardson number 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190627132031350-0413:S0022112019004300:S0022112019004300_inline2.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 and turbulent flux coefficient 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190627132031350-0413:S0022112019004300:S0022112019004300_inline3.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 are constant across a large portion of the depth. Changing the slope angle 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190627132031350-0413:S0022112019004300:S0022112019004300_inline4.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 modifies these mixing parameters, and the flow approaches a maximum Richardson number 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190627132031350-0413:S0022112019004300:S0022112019004300_inline5.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 as 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190627132031350-0413:S0022112019004300:S0022112019004300_inline6.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 at which the entrainment coefficient 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190627132031350-0413:S0022112019004300:S0022112019004300_inline7.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉. The turbulent Prandtl number remains 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190627132031350-0413:S0022112019004300:S0022112019004300_inline8.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 for all slope angles, demonstrating that 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190627132031350-0413:S0022112019004300:S0022112019004300_inline9.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 is not caused by a switch-off of the turbulent buoyancy flux as conjectured by Ellison (〈span〉J. Fluid Mech.〈/span〉, vol. 2, 1957, pp. 456–466). Instead, 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190627132031350-0413:S0022112019004300:S0022112019004300_inline10.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 occurs as the result of the turbulence intensity going to zero as 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190627132031350-0413:S0022112019004300:S0022112019004300_inline11.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉, due to the flow requiring larger and larger shear to maintain the same level of turbulence. We develop an approximate model valid for small 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190627132031350-0413:S0022112019004300:S0022112019004300_inline12.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 which is able to predict accurately 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190627132031350-0413:S0022112019004300:S0022112019004300_inline13.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉, 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190627132031350-0413:S0022112019004300:S0022112019004300_inline14.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 and 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190627132031350-0413:S0022112019004300:S0022112019004300_inline15.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 as a function of 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190627132031350-0413:S0022112019004300:S0022112019004300_inline16.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 and their maximum attainable values. The model predicts an entrainment law of the form 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190627132031350-0413:S0022112019004300:S0022112019004300_inline17.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉, which is in good agreement with the simulation data. The simulations and model presented here contribute to a growing body of evidence that an approach to a marginally or critically stable, relatively weakly stratified equilibrium for stratified shear flows may well be a generic property of turbulent stratified flows.〈/p〉〈/div〉

Description: 〈div data-abstract-type="normal"〉〈p〉Ice scallops are a small-scale (5–20 cm) quasi-periodic ripple pattern that occurs at the ice–water interface. Previous work has suggested that scallops form due to a self-reinforcing interaction between an evolving ice-surface geometry, an adjacent turbulent flow field and the resulting differential melt rates that occur along the interface. In this study, we perform a series of laboratory experiments in a refrigerated flume to quantitatively investigate the mechanisms of scallop formation and evolution in high resolution. Using particle image velocimetry, we probe an evolving ice–water boundary layer at sub-millimetre scales and 15 Hz frequency. Our data reveal three distinct regimes of ice–water interface evolution: a transition from flat to scalloped ice; an equilibrium scallop geometry; and an adjusting scallop interface. We find that scalloped-ice geometry produces a clear modification to the ice–water boundary layer, characterized by a time-mean recirculating eddy feature that forms in the scallop trough. Our primary finding is that scallops form due to a self-reinforcing feedback between the ice-interface geometry and shear production of turbulent kinetic energy in the flow interior. The length of this shear production zone is therefore hypothesized to set the scallop wavelength.〈/p〉〈/div〉

Description: 〈div data-abstract-type="normal"〉〈p〉Direct numerical simulation (DNS) is performed for two wall-bounded flow configurations: laminar Couette flow at 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190628070320625-0031:S0022112019004191:S0022112019004191_inline1.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 and turbulent channel flow at 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190628070320625-0031:S0022112019004191:S0022112019004191_inline2.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉, where 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190628070320625-0031:S0022112019004191:S0022112019004191_inline3.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 is the shear stress at the wall. The top wall is smooth and the bottom wall is a realistically rough superhydrophobic surface (SHS), generated from a three-dimensional surface profile measurement. The air–water interface, which is assumed to be flat, is simulated using the volume-of-fluid (VOF) approach. The two flow cases are studied with varying interface heights 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190628070320625-0031:S0022112019004191:S0022112019004191_inline4.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 to understand its effect on slip and drag reduction (〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190628070320625-0031:S0022112019004191:S0022112019004191_inline5.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉). For the laminar Couette flow case, the presence of the surface roughness is felt up to 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190628070320625-0031:S0022112019004191:S0022112019004191_inline6.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 of the channel height in the wall-normal direction. Nonlinear dependence of 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190628070320625-0031:S0022112019004191:S0022112019004191_inline7.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 on 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190628070320625-0031:S0022112019004191:S0022112019004191_inline8.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 is observed with three distinct regions. A nonlinear curve fit is obtained for gas fraction 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190628070320625-0031:S0022112019004191:S0022112019004191_inline9.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 as a function of 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190628070320625-0031:S0022112019004191:S0022112019004191_inline10.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉, where 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190628070320625-0031:S0022112019004191:S0022112019004191_inline11.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 determines the amount of slip area exposed to the flow. A power law fit is obtained from the data for the effective slip length as a function of 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190628070320625-0031:S0022112019004191:S0022112019004191_inline12.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 and is compared to those derived for structured geometry. For the turbulent channel flow, statistics of the flow field are compared to that of a smooth wall to understand the effects of roughness and 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190628070320625-0031:S0022112019004191:S0022112019004191_inline13.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉. Four cases are simulated ranging from fully wetted to fully covered and two intermediate regions in between. Scaling laws for slip length, slip velocity, roughness function and 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190628070320625-0031:S0022112019004191:S0022112019004191_inline14.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 are obtained for different penetration depths and are compared to past work for structured geometry. 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190628070320625-0031:S0022112019004191:S0022112019004191_inline15.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 is shown to depend on a competing effect between slip velocity and turbulent losses due to the Reynolds shear stress contribution. Presence of trapped air in the cavities significantly alters near-wall flow physics where we examine near-wall structures and propose a physical mechanism for their behaviour. The fully wetted roughness increases the peak value of turbulent intensities, whereas the presence of the interface suppresses them. The pressure fluctuations have competing contributions between turbulent pressure fluctuations and stagnation due to asperities, the near-wall structure is altered and breaks down with increasing slip. Overall, there exists a competing effect between the interface and the asperities, the interface suppresses turbulence whereas the asperities enhance them. The present work demonstrates DNS over a realistic multiphase SHS for the first time, to the best of our knowledge.〈/p〉〈/div〉

Description: 〈div data-abstract-type="normal"〉〈p〉We developed a numerical method for the set of equations governing fully compressible convection in the limit of infinite Prandtl numbers. Reduced models have also been analysed, such as the anelastic approximation and the anelastic liquid approximation. The tests of our numerical schemes against self-consistent criteria have shown that our numerical simulations are consistent from the point of view of energy dissipation, heat transfer and entropy budget. The equation of state of an ideal gas has been considered in this work. Specific effects arising because of the compressibility of the fluid are studied, like the scaling of viscous dissipation and the scaling of the heat flux contribution due to the mechanical power exerted by viscous forces. We analysed the solutions obtained with each model (fully compressible model, anelastic and anelastic liquid approximations) in a wide range of dimensionless parameters and determined the errors induced by each approximation with respect to the fully compressible solutions. Based on a rationale on the development of the thermal boundary layers, we can explain reasonably well the differences between the fully compressible and anelastic models, in terms of both the heat transfer and viscous dissipation dependence on compressibility. This could be mostly an effect of density variations on thermal diffusivity. Based on the different forms of entropy balance between exact and anelastic models, we find that a necessary condition for convergence of the anelastic results to the exact solutions is that the product 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190625082702251-0920:S0022112019004208:S0022112019004208_inline1.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 must be small compared to unity, where 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190625082702251-0920:S0022112019004208:S0022112019004208_inline2.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 is the ratio of the superadiabatic temperature difference to the adiabatic difference, and 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190625082702251-0920:S0022112019004208:S0022112019004208_inline3.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 is the ratio of the superadiabatic heat flux to the heat flux conducted along the adiabat. The same condition seems also to be associated with a convergence of the computed heat fluxes. Concerning the anelastic liquid approximation, we confirm previous estimates by Anufriev 〈span〉et al.〈/span〉 (〈span〉Phys. Earth Planet. Inter.〈/span〉, vol. 152, 2005, pp. 163–190) and find that its results become generally close to those of the fully compressible model when 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190625082702251-0920:S0022112019004208:S0022112019004208_inline4.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 is small compared to unity, where 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190625082702251-0920:S0022112019004208:S0022112019004208_inline5.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 is the isobaric thermal expansion coefficient, 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190625082702251-0920:S0022112019004208:S0022112019004208_inline6.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 is the temperature (here 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190625082702251-0920:S0022112019004208:S0022112019004208_inline7.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 for an ideal gas) and 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190625082702251-0920:S0022112019004208:S0022112019004208_inline8.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 is the dissipation number.〈/p〉〈/div〉

Description: 〈div data-abstract-type="normal"〉〈p〉Progress in roughness research, mapping any given roughness geometry to its fluid dynamic behaviour, has been hampered by the lack of accurate and direct measurements of skin-friction drag, especially in open systems. The Taylor–Couette (TC) system has the benefit of being a closed system, but its potential for characterizing irregular, realistic, three-dimensional (3-D) roughness has not been previously considered in depth. Here, we present direct numerical simulations (DNSs) of TC turbulence with sand grain roughness mounted on the inner cylinder. The model proposed by Scotti (〈span〉Phys. Fluids〈/span〉, vol. 18, 031701, 2006) has been modified to simulate a random rough surface of monodisperse sand grains. Taylor numbers range from 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190621104521411-0582:S0022112019003768:S0022112019003768_inline1.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉(corresponding to 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190621104521411-0582:S0022112019003768:S0022112019003768_inline2.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉) to 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190621104521411-0582:S0022112019003768:S0022112019003768_inline3.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 (〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190621104521411-0582:S0022112019003768:S0022112019003768_inline4.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉). We focus on the influence of the roughness height 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190621104521411-0582:S0022112019003768:S0022112019003768_inline5.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 in the transitionally rough regime, through simulations of TC with rough surfaces, ranging from 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190621104521411-0582:S0022112019003768:S0022112019003768_inline6.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 up to 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190621104521411-0582:S0022112019003768:S0022112019003768_inline7.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉. We analyse the global response of the system, expressed both by the dimensionless angular velocity transport 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190621104521411-0582:S0022112019003768:S0022112019003768_inline8.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 and by the friction factor 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190621104521411-0582:S0022112019003768:S0022112019003768_inline9.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉. An increase in friction with increasing roughness height is accompanied with enhanced plume ejection from the inner cylinder. Subsequently, we investigate the local response of the fluid flow over the rough surface. The equivalent sand grain roughness 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190621104521411-0582:S0022112019003768:S0022112019003768_inline10.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 is calculated to be 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190621104521411-0582:S0022112019003768:S0022112019003768_inline11.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉, where 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190621104521411-0582:S0022112019003768:S0022112019003768_inline12.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 is the size of the sand grains. We find that the downwards shift of the logarithmic layer, due to transitionally rough sand grains exhibits remarkably similar behaviour to that of the Nikuradse (〈span〉VDI-Forsch.〈/span〉, vol. 361, 1933) data of sand grain roughness in pipe flow, regardless of the Taylor number dependent constants of the logarithmic layer. Furthermore, we find that the dynamical effects of the sand grains are contained to the roughness sublayer 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190621104521411-0582:S0022112019003768:S0022112019003768_inline13.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 with 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190621104521411-0582:S0022112019003768:S0022112019003768_inline14.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉.〈/p〉〈/div〉

Description: 〈div data-abstract-type="normal"〉〈p〉Geometric, kinematic and dynamic properties of focusing deep-water surface gravity wave packets are examined in a simplified model with the intent of deriving a wave breaking threshold parameter. The model is based on the spatial modified nonlinear Schrödinger equation of Dysthe (〈span〉Proc. R. Soc. Lond.〈/span〉 A, vol. 369 (1736), 1979, pp. 105–114). The evolution of initially narrow-banded and weakly nonlinear chirped Gaussian wave packets are examined, by means of a trial function and a variational procedure, yielding analytic solutions describing the approximate evolution of the packet width, amplitude, asymmetry and phase during focusing. A model for the maximum free surface gradient, as a function of 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190621105221417-0357:S0022112019004282:S0022112019004282_inline1.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 and 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190621105221417-0357:S0022112019004282:S0022112019004282_inline2.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉, for 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190621105221417-0357:S0022112019004282:S0022112019004282_inline3.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 the linear prediction of the maximum slope at focusing and 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190621105221417-0357:S0022112019004282:S0022112019004282_inline4.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 the non-dimensional packet bandwidth, is proposed and numerically examined, indicating a quasi-self-similarity of these focusing events. The equations of motion for the fully nonlinear potential flow equations are then integrated to further investigate these predictions. It is found that a model of this form can characterize the bulk partitioning of 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190621105221417-0357:S0022112019004282:S0022112019004282_inline5.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 phase space, between non-breaking and breaking waves, serving as a breaking criterion. Application of this result to better understanding air–sea interaction processes is discussed.〈/p〉〈/div〉

Description: 〈div data-abstract-type="normal"〉〈p〉We conduct direct numerical simulations (DNS) of the Cahn–Hilliard–Navier–Stokes (CHNS) equations to investigate the statistical properties of a turbulent phase-separating symmetric binary-fluid mixture. Turbulence causes an arrest of the phase separation which leads to the formation of a statistically steady emulsion. We characterise turbulent velocity fluctuations in an emulsion for different values of the Reynolds number and the Weber number. Our scale-by-scale kinetic energy budget analysis shows that the interfacial terms in the CHNS equations provide an alternative route for the kinetic energy transfer. By studying the probability distribution function (p.d.f.) of the energy dissipation rate, the vorticity magnitude and the joint-p.d.f. of the velocity-gradient invariants we show that the statistics of the turbulent fluctuations do not change with the Weber number.〈/p〉〈/div〉

Description: 〈div data-abstract-type="normal"〉〈p〉With the aim to characterize the near-wall flow structures and their interaction with large-scale motions in the log-law region, time-resolved planar and volumetric flow field measurements were performed in the near-wall and log-law region of an adverse pressure gradient turbulent boundary layer following a zero pressure gradient turbulent boundary layer at a friction Reynolds number 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190621113954055-0112:S0022112019004087:S0022112019004087_inline1.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉. Due to the high spatial and temporal resolution of the measurements, it was possible to resolve and identify uniform-momentum zones in the region 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190621113954055-0112:S0022112019004087:S0022112019004087_inline2.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 or 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190621113954055-0112:S0022112019004087:S0022112019004087_inline3.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 and to relate them with well known coherent flow motions near the wall. The space–time results confirm that the turbulent superstructures have a strong impact even on the very near-wall flow motion and also their alternating appearance in time and intensity could be quantified over long time sequences. Using the time record of the velocity field, rare localized separation events appearing in the viscous sublayer were also analysed. By means of volumetric particle tracking velocimetry their three-dimensional topology and dynamics could be resolved. Based on the results, a conceptual model was deduced that explains their rare occurrence, topology and dynamics by means of a complex interaction process between low-momentum turbulent superstructures, near-wall low-speed streaks and tilted longitudinal and spanwise vortices located in the near-wall region.〈/p〉〈/div〉

Description: 〈div data-abstract-type="normal"〉〈p〉We experimentally investigate the extensional flow of a sheet – or curtain – of viscoelastic liquid falling freely from a slot at constant flow rate under gravity. Extruded liquids are aqueous solutions of flexible polyethylene oxide (PEO) and of semi-rigid partially hydrolysed polyacrylamide (HPAM) with low shear viscosities. Velocimetry measurements reveal that the mean velocity field 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190621102720940-0917:S0022112019003896:S0022112019003896_inline1.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 (where 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190621102720940-0917:S0022112019003896:S0022112019003896_inline2.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 is the distance from the slot exit) does not reduce to a free fall. More precisely, we show that the liquid falls initially with sub-gravitational accelerations up to a distance from the slot which scales as 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190621102720940-0917:S0022112019003896:S0022112019003896_inline3.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 (where 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190621102720940-0917:S0022112019003896:S0022112019003896_inline4.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 is gravity and 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190621102720940-0917:S0022112019003896:S0022112019003896_inline5.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 is the extensional relaxation time of the liquid) due to the stretching of polymer molecules. Beyond this elastic length, inertia dominates and the local acceleration reaches the asymptotic free-fall value 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190621102720940-0917:S0022112019003896:S0022112019003896_inline6.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉. The length of the sub-gravitational part of the curtain is shown to be much larger than the equivalent viscous length 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190621102720940-0917:S0022112019003896:S0022112019003896_inline7.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 for Newtonian liquids of density 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190621102720940-0917:S0022112019003896:S0022112019003896_inline8.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 and dynamic viscosity 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190621102720940-0917:S0022112019003896:S0022112019003896_inline9.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 which is usually small compared to the curtain length. By analogy with Newtonian curtains, we show that the velocity field 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190621102720940-0917:S0022112019003896:S0022112019003896_inline10.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 rescales on a master curve. Besides, the flow is shown to be only weakly affected by the history of polymer deformations in the die upstream of the curtain. Furthermore, investigations on the curtain stability reveal that polymer addition reduces the minimum flow rate required to maintain a continuous sheet of liquid.〈/p〉〈/div〉

Description: 〈div data-abstract-type="normal"〉〈p〉Previous experiments have revealed that shock waves driven through dissipative media may become unstable, for example, in granular gases, and in molecular gases undergoing strong relaxation effects. The current paper addresses this problem of shock stability at the Euler and Navier–Stokes continuum levels in a system of disks (two-dimensional) undergoing activated inelastic collisions. The dynamics of shock formation and stability is found to be in very good agreement with earlier molecular dynamic simulations (Sirmas & Radulescu, 〈span〉Phys. Rev.〈/span〉 E, vol. 91, 2015, 023003). It was found that the modelling of shock instability requires the introduction of molecular noise for its development and sustenance. This is confirmed in two stability problems. In the first, the evolution of shock formation dynamics is monitored without noise, with only initial noise and with continuous molecular noise. Only the latter reproduces the results of shock instability of molecular dynamics simulations. In the second problem, the steady travelling wave solution is obtained for the shock structure in the inviscid and viscous limits and its nonlinear stability is studied with and without molecular fluctuations, again showing that instability can be sustained only in the presence of fluctuations. The continuum results show that instability takes the form of a rippled front of a wavelength comparable with the relaxation thickness of the steady shock wave, at scales at which molecular fluctuations become important, in excellent agreement with the molecular dynamic simulations.〈/p〉〈/div〉

Description: 〈div data-abstract-type="normal"〉〈p〉We performed two-way coupled direct numerical simulations of turbulent channel flow with Lagrangian tracking of small, heavy spheres at a dimensionless gravitational acceleration of 0.077 in wall units, which is based on the flow condition in the experiment by Gerashchenko 〈span〉et al.〈/span〉 (〈span〉J. Fluid Mech.〈/span〉, vol. 617, 2008, pp. 255–281). We removed deposited particles after several collisions with the lower wall and then released new particles near the upper wall to observe direct interactions between particles and coherent structures of near-wall turbulence during gravitational settling through the mean shear. The results indicate that when the Stokes number is approximately 1 on the basis of the Kolmogorov time scale of the flow (〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190624134045701-0463:S0022112019004002:S0022112019004002_inline1.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉), the so-called preferential sweeping occurs in association with coherent streamwise vortices, while the effect of crossing trajectories becomes significant for 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190624134045701-0463:S0022112019004002:S0022112019004002_inline2.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉. Consequently, in either case, the settling particles deposit on the wall without strong accumulation in low-speed streaks in the viscous sublayer. When particles settle through near-wall turbulence from the upper wall, more small-scale vortical structures are generated in the outer layer as low-speed fluid is pulled farther in the direction of gravity, while the opposite is true near the lower wall.〈/p〉〈/div〉

Description: 〈div data-abstract-type="normal"〉〈p〉Streaks have been found to be an important part of wall-turbulence dynamics. In this paper, we extend the analysis for unbounded shear flows, in particular a Mach 0.4 round jet, using measurements taken using dual-plane, time-resolved, stereoscopic particle image velocimetry (PIV) taken at pairs of jet cross-sections, allowing the evaluation of the cross-spectral density of streamwise velocity fluctuations resolved into azimuthal Fourier modes. From the streamwise velocity results, two analyses are performed: the evaluation of wavenumber spectra (assuming Taylor’s hypothesis for the streamwise coordinate) and a spectral proper orthogonal decomposition (SPOD) of the velocity field using PIV planes in several axial stations. The methods complement each other, leading to the conclusion that large-scale streaky structures are also present in turbulent jets where they experience large growth in the streamwise direction, energetic structures extending up to eight diameters from the nozzle exit. Leading SPOD modes highlight the large-scale, streaky shape of the structures, whose aspect ratio (streamwise over azimuthal length) is approximately 15. The data were further analysed using SPOD, resolvent and transient growth analyses, good agreement being observed between the models and the leading SPOD mode for the wavenumbers considered. The models also indicate that the lift-up mechanism is active in turbulent jets, with streamwise vortices leading to streaks. The results show that large-scale streaks are a relevant part of the jet dynamics.〈/p〉〈/div〉

Description: 〈div data-abstract-type="normal"〉〈p〉We study single-phase and particle-laden turbulent channel flows bounded by two incompressible hyper-elastic walls with different deformability at bulk Reynolds number 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190621111752900-0025:S0022112019004130:S0022112019004130_inline1.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉. The solid volume fraction of finite-size neutrally buoyant rigid spherical particles considered is 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190621111752900-0025:S0022112019004130:S0022112019004130_inline2.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉. The elastic walls are assumed to be of a neo-Hookean material. A fully Eulerian formulation is employed to model the elastic walls together with a direct-forcing immersed boundary method for the coupling between the fluid and the particles. The data show a significant drag increase and the enhancement of the turbulence activity with growing wall elasticity for both the single-phase and particle-laden flows when compared with the single-phase flow over rigid walls. Drag reduction and turbulence attenuation is obtained, on the other hand, with highly elastic walls when comparing the particle-laden flow with the single-phase flow for the same wall properties; the opposite effect, drag increase, is observed upon adding particles to the flow over less elastic walls. This is explained by investigating the near-wall turbulence, where the strong asymmetry in the magnitude of the wall-normal velocity fluctuations (favouring positive 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190621111752900-0025:S0022112019004130:S0022112019004130_inline3.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉), is found to push the particles towards the channel centre. The particle layer close to the wall contributes to turbulence production by increasing the wall-normal velocity fluctuations, so that in the absence of this layer, smaller wall deformations and in turn turbulence attenuation is observed. For a moderate wall elasticity, we increase the particle volume fraction up to 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190621111752900-0025:S0022112019004130:S0022112019004130_inline4.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 and find that particle migration away from the wall is the cause of turbulence attenuation with respect to the flow over rigid walls. However, for this higher volume fractions, the particle induced stress compensates for the decreasing Reynolds shear stress, resulting in a higher overall drag for the case with elastic walls. The effect of the wall elasticity on the overall drag reduces significantly with increasing particle volume fraction.〈/p〉〈/div〉

Description: 〈div data-abstract-type="normal"〉〈p〉We present large-eddy simulation (LES) of flow past different airfoils with 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190621103821391-0158:S0022112019003604:S0022112019003604_inline1.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉, based on the free-stream velocity and airfoil chord length, ranging from 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190621103821391-0158:S0022112019003604:S0022112019003604_inline2.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 to 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190621103821391-0158:S0022112019003604:S0022112019003604_inline3.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉. To avoid the challenging resolution requirements of the near-wall region, we develop a virtual wall model in generalized curvilinear coordinates and incorporate the non-equilibrium effects via proper treatment of the momentum equations. It is demonstrated that the wall model dynamically captures the instantaneous skin-friction vector field on arbitrary curved surfaces at the resolved scale. By combining the present wall model with the stretched-vortex subgrid-scale model, we apply the wall-modelled LES approach to three different airfoil cases, spanning different geometrical parameters, different attack angles and low to high 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190621103821391-0158:S0022112019003604:S0022112019003604_inline4.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉. The numerical results are verified with direct numerical simulation (DNS) at low 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190621103821391-0158:S0022112019003604:S0022112019003604_inline5.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉, and validated with experiment data at higher 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190621103821391-0158:S0022112019003604:S0022112019003604_inline6.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉, including typical aerodynamic properties such as pressure coefficient distributions, velocity components and also more challenging measurements such as skin-friction coefficient and Reynolds stresses. All comparisons show reasonable agreement, providing a measure of validity that enables us to further probe simulation results into aspects of flow physics that are not available from experiments. Two techniques to quantify hitherto unexplored physics of flows past airfoils are employed: one is the construction of the anisotropy invariant map, and the second is skin-friction portraits with emphasis on flow transition and unsteady separation along the airfoil surface. The anisotropy maps for all three 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190621103821391-0158:S0022112019003604:S0022112019003604_inline7.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 cases, show clearly that a portion of the flow field is aligned along the axisymmetric expansion line, corresponding to the turbulent boundary layer log-law behaviour and the appearance of turbulent transition. The instantaneous skin-friction portraits reveal a monotonic shrinking of the near wall structure scale. At 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190621103821391-0158:S0022112019003604:S0022112019003604_inline8.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉, the interaction between the primary separation bubble and the secondary separation bubble contributes to turbulent transition, similar to the case of flow past a cylinder. At higher 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190621103821391-0158:S0022112019003604:S0022112019003604_inline9.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉, the primary separation breaks into several small separation bubbles. At even higher 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190621103821391-0158:S0022112019003604:S0022112019003604_inline10.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉, near the turbulent separation, the skin-friction lines show small-scale reversal flows that are similar to those observed in DNS of the flat plate turbulent separation. A notable feature of turbulent separation in flow past an airfoil is the appearance of turbulence structures and small-scale reversal flows in the spanwise direction due to the vortex shedding behaviour.〈/p〉〈/div〉

Description: 〈div data-abstract-type="normal"〉〈p〉In the framework of the complete formulation of the conjugate problem, the liquid–gas flow structure arising upon local heating using thermal sources is investigated numerically. The two-layer system is confined by solid impermeable walls. The Navier–Stokes equations in the Boussinesq approximation in the ‘streamfunction–vorticity’ variables are used to describe the media motion. The dynamic conditions at the interface are formulated in terms of the tangential and normal velocities, while the temperature conditions at the external boundaries of the system take into account the presence of local heaters. The influence of the number of heaters and heating modes on the dynamics and character of the appearing convective regimes is analysed. The steady and commutated heating modes for one and two heaters arranged at the lower boundary are investigated. The heating initiates convective and thermocapillary mechanisms causing the fluid motion. Transient regimes with the successive formation of two-vortex, quadruple-vortex and two-vortex flows are observed before the stabilization of the system in the uniform heating mode. A stable thermocapillary deflection appears at the interface above the heater. The commutated mode of heating entails oscillations of the interface with a change in the deflection form and the formation of travelling vortices in the fluids. The impact of particular mechanisms on the flow patterns is analysed. The paper presents typical distributions of the velocity and temperature fields in the system and the position of the interface for the considered cases.〈/p〉〈/div〉

Description: 〈div data-abstract-type="normal"〉〈p〉In the context of dynamic wetting, wall slip is often treated as a microscopic effect for removing viscous stress singularity at a moving contact line. In most drop spreading experiments, however, a considerable amount of slip may occur due to the use of polymer liquids such as silicone oils, which may cause significant deviations from the classical Tanner–de Gennes theory. Here we show that many classical results for complete wetting fluids may no longer hold due to wall slip, depending crucially on the extent of de Gennes’s slipping ‘foot’ to the relevant length scales at both the macroscopic and microscopic levels. At the macroscopic level, we find that for given liquid height 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190619140754267-0248:S0022112019003525:S0022112019003525_inline1.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 and slip length 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190619140754267-0248:S0022112019003525:S0022112019003525_inline2.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉, the apparent dynamic contact angle 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190619140754267-0248:S0022112019003525:S0022112019003525_inline3.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 can change from Tanner’s law 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190619140754267-0248:S0022112019003525:S0022112019003525_inline4.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 for 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190619140754267-0248:S0022112019003525:S0022112019003525_inline5.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 to the strong-slip law 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190619140754267-0248:S0022112019003525:S0022112019003525_inline6.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 for 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190619140754267-0248:S0022112019003525:S0022112019003525_inline7.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉, where 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190619140754267-0248:S0022112019003525:S0022112019003525_inline8.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 is the capillary number and 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190619140754267-0248:S0022112019003525:S0022112019003525_inline9.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 is the macroscopic length scale. Such a no-slip-to-slip transition occurs at the critical capillary number 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190619140754267-0248:S0022112019003525:S0022112019003525_inline10.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉, accompanied by the switch of the ‘foot’ of size 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190619140754267-0248:S0022112019003525:S0022112019003525_inline11.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 from the inner scale to the outer scale with respect to 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190619140754267-0248:S0022112019003525:S0022112019003525_inline12.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉. A more generalized dynamic contact angle relationship is also derived, capable of unifying Tanner’s law and the strong-slip law under 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190619140754267-0248:S0022112019003525:S0022112019003525_inline13.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉. We not only confirm the two distinct wetting laws using many-body dissipative particle dynamics simulations, but also provide a rational account for anomalous departures from Tanner’s law seen in experiments (Chen, 〈span〉J. Colloid Interface Sci〈/span〉., vol. 122, 1988, pp. 60–72; Albrecht 〈span〉et al.〈/span〉, 〈span〉Phys. Rev. Lett.〈/span〉, vol. 68, 1992, pp. 3192–3195). We also show that even for a common spreading drop with small macroscopic slip, slip effects can still be microscopically strong enough to change the microstructure of the contact line. The structure is identified to consist of a strongly slipping precursor film of length 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190619140754267-0248:S0022112019003525:S0022112019003525_inline14.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 followed by a mesoscopic ‘foot’ of width 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190619140754267-0248:S0022112019003525:S0022112019003525_inline15.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 ahead of the macroscopic wedge, where 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190619140754267-0248:S0022112019003525:S0022112019003525_inline16.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 is the molecular length. It thus turns out that it is the ‘foot’, rather than the film, contributing to the microscopic length in Tanner’s law, in accordance with the experimental data reported by Kavehpour 〈span〉et al.〈/span〉 (〈span〉Phys. Rev. Lett.〈/span〉, vol. 91, 2003, 196104) and Ueno 〈span〉et al.〈/span〉 (〈span〉Trans. ASME J. Heat Transfer〈/span〉, vol. 134, 2012, 051008). The advancement of the microscopic contact line is still led by the film whose length can grow as the 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190619140754267-0248:S0022112019003525:S0022112019003525_inline17.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 power of time due to 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190619140754267-0248:S0022112019003525:S0022112019003525_inline18.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉, as supported by the experiments of Ueno 〈span〉et al.〈/span〉 and Mate (〈span〉Langmuir〈/span〉, vol. 28, 2012, pp. 16821–16827). The present work demonstrates that the behaviour of a moving contact line can be strongly influenced by wall slip. Such slip-mediated dynamic wetting might also provide an alternative means for probing slippery surfaces.〈/p〉〈/div〉

Description: 〈div data-abstract-type="normal"〉〈p〉We numerically investigated the unsteady dynamics of a two-dimensional airfoil undergoing a continuous, prescribed pitch-up motion and freely translating as a response to aerodynamic forces and the gravity field. The pitch-up motion was applied about an axis located 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190619142124772-0727:S002211201900421X:S002211201900421X_inline1.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 chord away from the leading edge and was parameterized using the shape change number, with a Reynolds number set to 2000. It was shown that the minimum kinetic energy reached by the airfoil depends stochastically and asymptotically on shape change numbers for values below and above 1, respectively. Very low kinetic energy levels (close to zero) can be reached in both stochastic and asymptotic regions but high shape change numbers are accompanied by significant gain in altitude which may be undesirable from a practical perspective. Rather, shape change numbers in the range [0.1–0.3] allow us to reach relatively low levels of kinetic energy for close perching locations. We showed that highly nonlinear fluid–structure interactions induced by massive flow separations and strong vortices are conducive to low kinetic energy, but responsible for the stochastic dependence of kinetic energy to shape change number, which can make perching manoeuvres hardly controllable for flying vehicles.〈/p〉〈/div〉

Description: 〈div data-abstract-type="normal"〉〈p〉We present a dynamic decomposition analysis of the wake flow in fluid–structure interaction (FSI) systems under both laminar and turbulent flow conditions. Of particular interest is to provide the significance of low-dimensional wake flow features and their interaction dynamics to sustain the free vibration of a square cylinder at a relatively low mass ratio. To obtain the high-dimensional data, we employ a body-conforming variational FSI solver based on the recently developed partitioned iterative scheme and the dynamic subgrid-scale turbulence model for a moderate Reynolds number (〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190327055129472-0988:S002211201900140X:S002211201900140X_inline1.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉). The snapshot data from high-dimensional FSI simulations are projected to a low-dimensional subspace using the proper orthogonal decomposition (POD). We utilize each corresponding POD mode to detect features of the organized motions, namely, the vortex street, the shear layer and the near-wake bubble. We find that the vortex shedding modes contribute solely to the lift force, while the near-wake and shear layer modes play a dominant role in the drag force. We further examine the fundamental mechanism of this dynamical behaviour and propose a force decomposition technique via low-dimensional approximation. To elucidate the frequency lock-in, we systematically analyse the decomposed modes and their dynamical contributions to the force fluctuations for a range of reduced velocity at low Reynolds number laminar flow. These quantitative mode energy contributions demonstrate that the shear layer feeds the vorticity flux to the wake vortices and the near-wake bubble during the wake–body synchronization. Based on the decomposition of wake dynamics, we suggest an interaction cycle for the frequency lock-in during the wake–body interaction, which provides the interrelationship between the high-amplitude motion and the dominating wake features. Through our investigation of wake–body synchronization below critical 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190327055129472-0988:S002211201900140X:S002211201900140X_inline2.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 range, we discover that the bluff body can undergo a synchronized high-amplitude vibration due to flexibility-induced unsteadiness. Owing to the wake turbulence at a moderate Reynolds number of 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190327055129472-0988:S002211201900140X:S002211201900140X_inline3.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉, a distorted set of POD modes and the broadband energy distribution are observed, while the interaction cycle for the wake synchronization is found to be valid for the turbulent wake flow.〈/p〉〈/div〉

Description: 〈div data-abstract-type="normal"〉〈p〉We use resolvent analysis to design active control techniques for separated flows over a NACA 0012 airfoil. Spanwise-periodic flows over the airfoil at a chord-based Reynolds number of 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190325174522297-0677:S0022112019001630:S0022112019001630_inline1.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 and a free-stream Mach number of 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190325174522297-0677:S0022112019001630:S0022112019001630_inline2.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 are considered at two post-stall angles of attack of 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190325174522297-0677:S0022112019001630:S0022112019001630_inline3.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 and 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190325174522297-0677:S0022112019001630:S0022112019001630_inline4.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉. Near the leading edge, localized unsteady thermal actuation is introduced in an open-loop manner with two tunable parameters of actuation frequency and spanwise wavelength. To provide physics-based guidance for the effective choice of these control input parameters, we conduct global resolvent analysis on the baseline turbulent mean flows to identify the actuation frequency and wavenumber that provide large perturbation energy amplification. The present analysis also considers the use of a temporal filter to limit the time horizon for assessing the energy amplification to extend resolvent analysis to unstable base flows. We incorporate the amplification and response mode from resolvent analysis to provide a metric that quantifies momentum mixing associated with the modal structure. This metric is compared to the results from a large number of three-dimensional large-eddy simulations of open-loop controlled flows. With the agreement between the resolvent-based metric and the enhancement of aerodynamic performance found through large-eddy simulations, we demonstrate that resolvent analysis can predict the effective range of actuation frequency as well as the global response to the actuation input. We believe that the present resolvent-based approach provides a promising path towards mean flow modification by capitalizing on the dominant modal mixing.〈/p〉〈/div〉

Description: 〈div data-abstract-type="normal"〉〈p〉In this work we present and demonstrate the reliability of a theoretical framework for the study of thermally driven turbulence. It consists of scale-by-scale budget equations for the second-order velocity and temperature structure functions and their limiting cases, represented by the turbulent kinetic energy and temperature variance budgets. This framework represents an extension of the classical Kolmogorov and Yaglom equations to inhomogeneous and anisotropic flows, and allows for a novel assessment of the turbulent processes occurring at different scales and locations in the fluid domain. Two relevant characteristic scales, 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190327124217270-0217:S0022112019001198:S0022112019001198_inline1.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 for the velocity field and 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190327124217270-0217:S0022112019001198:S0022112019001198_inline2.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 for the temperature field, are identified. These variables separate the space of scales into a quasi-homogeneous range, characterized by turbulent kinetic energy and temperature variance cascades towards dissipation, and an inhomogeneity-dominated range, where the production and the transport in physical space are important. This theoretical framework is then extended to the context of large-eddy simulation to quantify the effect of a low-pass filtering operation on both resolved and subgrid dynamics of turbulent Rayleigh–Bénard convection. It consists of single-point and scale-by-scale budget equations for the filtered velocity and temperature fields. To evaluate the effect of the filter length 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190327124217270-0217:S0022112019001198:S0022112019001198_inline3.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 on the resolved and subgrid dynamics, the velocity and temperature fields obtained from a direct numerical simulation are split into filtered and residual components using a spectral cutoff filter. It is found that when 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190327124217270-0217:S0022112019001198:S0022112019001198_inline4.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 is smaller than the minimum values of the cross-over scales given by 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190327124217270-0217:S0022112019001198:S0022112019001198_inline5.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉, the resolved processes correspond to the exact ones, except for a depletion of viscous and thermal dissipations, and the only role of the subgrid scales is to drain turbulent kinetic energy and temperature variance to dissipate them. On the other hand, the resolved dynamics is much poorer in the near-wall region and the effects of the subgrid scales are more complex for filter lengths of the order of 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190327124217270-0217:S0022112019001198:S0022112019001198_inline6.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 or larger. This study suggests that classic eddy-viscosity/diffusivity models employed in large-eddy simulation may suffer from some limitations for large filter lengths, and that alternative closures should be considered to account for the inhomogeneous processes at subgrid level. Moreover, the theoretical framework based on the filtered Kolmogorov and Yaglom equations may represent a valuable tool for future assessments of the subgrid-scale models.〈/p〉〈/div〉

Description: 〈div data-abstract-type="normal"〉〈p〉Motivated by recent studies that have revealed the existence of trapped acoustic waves in subsonic jets (Towne 〈span〉et al.〈/span〉, 〈span〉J. Fluid Mech.〈/span〉, vol. 825, 2017, pp. 1113–1152), we undertake a more general exploration of the physics associated with acoustic modes in jets and wakes, using a double vortex-sheet model. These acoustic modes are associated with eigenvalues of the vortex-sheet dispersion relation; they are discrete modes, guided by the vortex sheet; they may be either propagative or evanescent; and under certain conditions they behave in the manner of acoustic-duct modes. By analysing these modes we show how jets and wakes may both behave as waveguides under certain conditions, emulating ducts with soft or hard walls, with the vortex-sheet impedance providing effective ‘wall’ conditions. We consider, in particular, the role that upstream-travelling acoustic modes play in the dispersion-relation saddle points that underpin the onset of absolute instability. The analysis illustrates how departure from duct-like behaviour is a necessary condition for absolute instability, and this provides a new perspective on the stabilising and destabilising effects of reverse flow, temperature ratio and compressibility; it also clarifies the differing symmetries of jet (symmetric) and wake (antisymmetric) instabilities. An energy balance, based on the vortex-sheet impedance, is used to determine stability conditions for the acoustic modes: these may become unstable in supersonic flow due to an energy influx through the shear layers. Finally, we construct the impulse response of flows with zero and finite shear-layer thickness. This allows us to show how the long-time wavepacket behaviour is indeed determined by interaction between Kelvin–Helmholtz and acoustic modes.〈/p〉〈/div〉

Description: 〈div data-abstract-type="normal"〉〈p〉Following the recent observation that turbulent pipe flow can be relaminarised by a relatively simple modification of the mean velocity profile, we here carry out a quantitative experimental investigation of this phenomenon. Our study confirms that a flat velocity profile leads to a collapse of turbulence and in order to achieve the blunted profile shape, we employ a moving pipe segment that is briefly and rapidly shifted in the streamwise direction. The relaminarisation threshold and the minimum shift length and speeds are determined as a function of Reynolds number. Although turbulence is still active after the acceleration phase, the modulated profile possesses a severely decreased lift-up potential as measured by transient growth. As shown, this results in an exponential decay of fluctuations and the flow relaminarises. While this method can be easily applied at low to moderate flow speeds, the minimum streamwise length over which the acceleration needs to act increases linearly with the Reynolds number.〈/p〉〈/div〉

Description: 〈div data-abstract-type="normal"〉〈p〉We investigate the gravity-driven flow of a thin film of liquid metal on a conducting conical substrate in the presence of a strong toroidal magnetic field (transverse to the flow and parallel to the substrate). We solve the leading-order governing equations in a physically relevant asymptotic limit to find the free-surface profile. We find that the leading-order fluid flow rate is a non-monotonic bounded function of the film height, and this can lead to singularities in the free-surface profile. We perform a detailed stability analysis and identify values of the relevant geometric, hydrodynamic and magnetic parameters such that the flow is stable.〈/p〉〈/div〉

Description: 〈div data-abstract-type="normal"〉〈p〉We present a local stability analysis to investigate the effects of differential diffusion between momentum and density (quantified by the Schmidt number 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190327085207695-0702:S0022112019001472:S0022112019001472_inline1.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉) on the three-dimensional, short-wavelength instabilities in planar vortices with a uniform stable stratification along the vorticity axis. Assuming small diffusion in both momentum and density, but arbitrary values for 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190327085207695-0702:S0022112019001472:S0022112019001472_inline2.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉, we present a detailed analytical/numerical analysis for three different classes of base flows: (i) an axisymmetric vortex, (ii) an elliptical vortex and (iii) the flow in the neighbourhood of a hyperbolic stagnation point. While a centrifugally stable axisymmetric vortex remains stable for any 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190327085207695-0702:S0022112019001472:S0022112019001472_inline3.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉, it is shown that 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190327085207695-0702:S0022112019001472:S0022112019001472_inline4.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 can have significant effects in a centrifugally unstable axisymmetric vortex: the range of unstable perturbations increases when 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190327085207695-0702:S0022112019001472:S0022112019001472_inline5.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 is taken away from unity, with the extent of increase being larger for 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190327085207695-0702:S0022112019001472:S0022112019001472_inline6.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 than for 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190327085207695-0702:S0022112019001472:S0022112019001472_inline7.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉. Additionally, for 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190327085207695-0702:S0022112019001472:S0022112019001472_inline8.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉, we report the possibility of oscillatory instability. In an elliptical vortex with a stable stratification, 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190327085207695-0702:S0022112019001472:S0022112019001472_inline9.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 is shown to non-trivially influence the three different inviscid instabilities (subharmonic, fundamental and superharmonic) that have been previously reported, and also introduce a new branch of oscillatory instability that is not present at 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190327085207695-0702:S0022112019001472:S0022112019001472_inline10.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉. The unstable parameter space for the subharmonic (instability IA) and fundamental (instability IB) inviscid instabilities are shown to be significantly increased for 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190327085207695-0702:S0022112019001472:S0022112019001472_inline11.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 and 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190327085207695-0702:S0022112019001472:S0022112019001472_inline12.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉, respectively. Importantly, for sufficiently small and large 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190327085207695-0702:S0022112019001472:S0022112019001472_inline13.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉, respectively, the maximum growth rate for instabilities IA and IB occurs away from the inviscid limit. The new oscillatory instability (instability III) is shown to occur only for sufficiently small 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190327085207695-0702:S0022112019001472:S0022112019001472_inline14.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉, the signature of which is nevertheless present with zero growth rate in the inviscid limit. The Schmidt number is then shown to play no role in the evolution of transverse perturbations on the flow around a hyperbolic stagnation point with a stable stratification. We conclude by discussing the physical length scales associated with the 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190327085207695-0702:S0022112019001472:S0022112019001472_inline15.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 instabilities, and their potential relevance in various realistic settings.〈/p〉〈/div〉

Description: 〈div data-abstract-type="normal"〉〈p〉Instabilities and flow characteristics in the far wake of a circular cylinder are examined through direct numerical simulations. The transitions to the two-layered and secondary vortex streets are quantified by a new method based on the time-averaged transverse velocity field. Two processes for the transition to the secondary vortex street are observed: (i) the merging of two same-sign vortices over a range of low Reynolds numbers (〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190326132244592-0183:S0022112019001678:S0022112019001678_inline1.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉) between 200 and 300, and (ii) the pairing of two opposite-sign vortices, followed by the merging of the paired vortices into subsequent vortices, over a range of 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190326132244592-0183:S0022112019001678:S0022112019001678_inline2.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 between 400 and 1000. Single vortices may be generated between the merging cycles due to mismatch of the vortices. The irregular merging process results in flow irregularity and an additional frequency signal 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190326132244592-0183:S0022112019001678:S0022112019001678_inline3.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 (in addition to the primary vortex shedding frequency 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190326132244592-0183:S0022112019001678:S0022112019001678_inline4.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉) in the two-layered and secondary vortex streets. In particular, a gradual energy transfer from 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190326132244592-0183:S0022112019001678:S0022112019001678_inline5.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 to 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190326132244592-0183:S0022112019001678:S0022112019001678_inline6.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 with distance downstream is observed in the two-layered vortex street prior to the merging. The frequency spectra of 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190326132244592-0183:S0022112019001678:S0022112019001678_inline7.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 are broad-band for 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190326132244592-0183:S0022112019001678:S0022112019001678_inline8.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉–300 but become increasingly sharp-peaked with increasing 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190326132244592-0183:S0022112019001678:S0022112019001678_inline9.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 because the vortex merging process becomes increasingly regular. The ratio of the sharp-peaked frequencies 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190326132244592-0183:S0022112019001678:S0022112019001678_inline10.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 and 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190326132244592-0183:S0022112019001678:S0022112019001678_inline11.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 is equal to the ratio of the numbers of vortices observed after and before the merging. The general conclusions drawn from a circular cylinder are expected to be applicable to other bluff bodies.〈/p〉〈/div〉

Description: 〈div data-abstract-type="normal"〉〈p〉Results are presented from a laboratory study on the free-surface signal generated over an array of submerged circular cylinders, representative of submerged aquatic vegetation. We aim to understand whether aquatic ecosystems generate a surface signature that is indicative of both what is beneath the water surface as well as how it is altering the flow. A shear layer forms over the canopy, generating coherent vortex structures which eventually manifest in the free-surface slope field. We connect the vortex properties measured at the surface with measurements of the bulk flow, and show that correlations between these quantities are adequate to create a parameterized model in which the interior velocity profile can be predicted solely from measurements taken at the free surface. Experimental surface observations yield a Strouhal number that is twice the most amplified mode predicted by linear stability theory, suggesting that vortices may evolve between generation at the canopy height and their manifestation at the water surface. Additionally, the surface signal continues evolving with distance downstream, with vortices becoming spread farther apart and the passage frequency gradually decreasing. By the trailing edge of the canopy, surface-impacting boils emerge for canopies with higher submergence ratios, suggesting a transition from coherent two-dimensional rollers to transversely varying structures.〈/p〉〈/div〉

Description: 〈div data-abstract-type="normal"〉〈p〉Results of experimental and numerical investigations of a supersonic flow around a cylinder with a frontal gas-permeable insert made of a high-porosity cellular material are presented. The measurements are performed in a T-327 supersonic blowdown wind tunnel at the free-stream Mach numbers 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190325181423329-0720:S0022112019001654:S0022112019001654_inline1.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉, 7 and 21 in the range of the unit Reynolds numbers 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190325181423329-0720:S0022112019001654:S0022112019001654_inline2.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉. The drag coefficients for a cylinder with an aerospike and a cylinder with a frontal gas-permeable porous insert are obtained. For the cylinder with the frontal gas-permeable porous insert, variations of the insert length, cylinder diameter and pore size are considered, and the mechanism of drag reduction is found, which includes two supplementary processes: attenuation of the bow shock wave in a system of weaker shock waves, and formation of an effective pointed body. The experiments are accompanied by numerical simulations of the flow around the cylinder with the frontal high-porosity insert: the fields of parameters of the external flow and the flow inside the porous insert are obtained, the drag coefficients are calculated, and the shape of the effective body for the examined model is found. The structure of the high-porosity material is modelled by a system of staggered rings of different diameters aligned in the radial and longitudinal directions (skeleton model of a porous medium). Numerical simulations of the problem are performed by means of solving two-dimensional Reynolds-averaged Navier–Stokes equations written in an axisymmetric form. The experimental and numerical data reveal significant drag reduction in a wide range of supersonic flow conditions. The results obtained on the drag coefficient for the cylinder are generalized with the use of a parameter which includes the ratio of the cylinder diameter to the pore diameter in the insert and the Mach number. This parameter is proposed as a similarity criterion for the problem of a supersonic flow around a cylinder with a frontal high-porosity insert.〈/p〉〈/div〉

Description: 〈div data-abstract-type="normal"〉〈p〉Decay of honeycomb-generated turbulence in a duct with a static transverse magnetic field is studied via direct numerical simulations. The simulations follow the revealing experimental study of Sukoriansky 〈span〉et al.〈/span〉 (〈span〉Exp. Fluids〈/span〉, vol. 4 (1), 1986, pp. 11–16), in particular the paradoxical observation of high-amplitude velocity fluctuations, which exist in the downstream portion of the flow when the strong transverse magnetic field is imposed in the entire duct including the honeycomb exit, but not in other configurations. It is shown that the fluctuations are caused by the large-scale quasi-two-dimensional structures forming in the flow at the initial stages of the decay and surviving the magnetic suppression. Statistical turbulence properties, such as the energy decay curves, two-point correlations and typical length scales are computed. The study demonstrates that turbulence decay in the presence of a magnetic field is a complex phenomenon critically depending on the state of the flow at the moment the field is introduced.〈/p〉〈/div〉

Description: 〈div data-abstract-type="normal"〉〈p〉Faults in the Roer Valley Rift System (RVRS) act as barriers to horizontal groundwater flow. This causes steep cross-fault groundwater level steps (hydraulic head differences). An overview of the size and distribution of fault-related groundwater level steps and associated fault zone permeabilities is thus far lacking. Such an overview would provide useful insights for nature restoration projects in areas with shallow groundwater levels (〈span〉wijstgronden〈/span〉) on the foot wall of fault zones. In this review study, data on fault zone permeabilities and cross-fault hydraulic head differences were compiled from 39 sources of information, consisting of literature (starting from 1948), internal reports (e.g. from research institutes and drinking water companies), groundwater models, a geological database and new fieldwork. The data are unevenly distributed across the RVRS. Three-quarters of the data sources are related to the Peel Boundary Fault zone (PBFZ). This bias is probably caused by the visibility of fault scarps and fault-adjacent wet areas for the PBFZ, with the characteristic iron-rich groundwater seepage. Most data demonstrate a cross-fault phreatic groundwater level step of 1.0 to 2.5 m. Data for the Feldbiss Fault zone (FFZ) show phreatic cross-fault hydraulic head differences of 1.0 to 2.0 m. 〈span〉In situ〈/span〉 measured hydraulic conductivity data (K) are scarce. Values vary over three orders of magnitude, from 0.013 to 22.1 m d〈span〉−1〈/span〉, are non-directional and do not take into account heterogeneity caused by fault zones. The hydraulic conductivity (and hydraulic resistance) values used in three different groundwater models are obtained by calibration using field measurements. They also cover a large range, from 0.001 to 32 m d〈span〉−1〈/span〉 and from 5 to 100,000 days. Heterogeneity is also not taken into account in these models. The overview only revealed locations with a clear cross-fault groundwater level step, and at many locations the faults are visible on aerial photographs as cropmarks or as soil moisture contrasts at the surface. Therefore, it seems likely that all faults have a reduced permeability, which determines the size of the groundwater level steps. In addition, our results show that cross-fault hydraulic head gradients also correlate with topographic, fault-induced offsets, for both the Peel Boundary and the Feldbiss fault zone.〈/p〉〈/div〉

Description: 〈div data-abstract-type="normal"〉〈p〉〈img orientation="portrait" mimesubtype="gif" mimetype="image" position="float" type="simple" href="S0022112019008000_figAb" src="http://static.cambridge.org/content/id/urn%3Acambridge.org%3Aid%3Aarticle%3AS0022112019008000/resource/name/S0022112019008000_figAb.gif?pub-status=live"〉〈/p〉〈/div〉 〈div data-abstract-type="normal"〉〈p〉We consider a Stokeslet applied to a viscous fluid next to an infinite, flat wall, or in between two parallel walls. We calculate the forces exerted by the resulting flow on the confining boundaries, and use the results obtained to estimate the hydrodynamic contribution to the pressure exerted on boundaries by force-free self-propelled particles.〈/p〉〈/div〉

Description: 〈div data-abstract-type="normal"〉〈p〉〈img orientation="portrait" mimesubtype="gif" mimetype="image" position="float" type="simple" href="S0022112019007973_figAb" src="http://static.cambridge.org/content/id/urn%3Acambridge.org%3Aid%3Aarticle%3AS0022112019007973/resource/name/S0022112019007973_figAb.gif?pub-status=live"〉〈/p〉〈/div〉 〈div data-abstract-type="normal"〉〈p〉By means of three-dimensional direct numerical simulations, we investigate the influence of the regular roughness of heated and cooled plates on the mean heat transport in a cylindrical Rayleigh–Bénard convection cell of aspect ratio one. The roughness is introduced by a set of isothermal obstacles, which are attached to the plates and have a form of concentric rings of the same width. The considered Prandtl number 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20191105072504660-0909:S0022112019007973:S0022112019007973_inline1.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 equals 1, the Rayleigh number 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20191105072504660-0909:S0022112019007973:S0022112019007973_inline2.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 varies from 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20191105072504660-0909:S0022112019007973:S0022112019007973_inline3.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 to 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20191105072504660-0909:S0022112019007973:S0022112019007973_inline4.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉, the number of rings on each plate is 1, 2, 4, 8 or 10, the height of the rings is varied from 1.5 % to 49 % of the cylinder height and the gap between the rings is varied from 1.5 % to 18.8 % of the cell diameter. Totally, 135 different cases are analysed. Direct numerical simulations show that with small 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20191105072504660-0909:S0022112019007973:S0022112019007973_inline5.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 and wide roughness rings, a small reduction of the mean heat transport (the Nusselt number 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20191105072504660-0909:S0022112019007973:S0022112019007973_inline6.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉) is possible, but, in most cases, the presence of the heated and cooled obstacles generally leads to an increase of 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20191105072504660-0909:S0022112019007973:S0022112019007973_inline7.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉, compared to the case of classical Rayleigh–Bénard convection with smooth plates. When the rings are very tall and the gaps between them are sufficiently wide, the effective mean heat flux can be several times larger than in the smooth case. For a fixed geometry of the obstacles, the scaling exponent in the 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20191105072504660-0909:S0022112019007973:S0022112019007973_inline8.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 versus 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20191105072504660-0909:S0022112019007973:S0022112019007973_inline9.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 scaling first increases with growing 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20191105072504660-0909:S0022112019007973:S0022112019007973_inline10.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 up to approximately 0.5, but then smoothly decreases back towards the exponent in the no-obstacle case.〈/p〉〈/div〉

Description: 〈div data-abstract-type="normal"〉〈p〉A stacked aeolian sequence with intercalated soils is presented from the southern Netherlands, which fully covers the Late Weichselian and Holocene periods. An integrated sedimentological (sedimentary structures, grain size), palynological (pollen) and dating approach (radiocarbon, optically stimulated luminescence (OSL)) was applied to unravel climatic and human forcing factors. The dating results of soils and sediments are compatible, and no large hiatuses between the radiocarbon-dated top of the soils and OSL-dated overlying sands were observed. It is argued that the peaty top of wet-type podzols can be used for reliable radiocarbon dating. This study reveals more phases than previously known of landscape stability (Usselo Soil and two podzol soils) and instability (Younger Coversand I and II, two drift-sand units) that are related to Late Weichselian climate change and Holocene human occupation. Regional aeolian deposition in source-bordering (river) dunes (Younger Coversand II) took place in the second part of the Younger Dryas, after 12.3 ka cal. BP, implying a delayed response to Younger Dryas cooling, vegetation cover decline and river pattern change of the Scheldt. The onset of podzolisation and development of ericaceous vegetation occurred prior to the introduction of Neolithic farming, which is earlier than previously assumed. Early podzolisation was followed by a short phase of local drift-sand deposition, at 〈span〉c〈/span〉.5500 cal. BP, that possibly relates to agriculture. Strong human impact on the landscape by deforestation and agriculture resulted in a second phase of widespread drift-sand deposition covering the younger podzol soil after AD 1000.〈/p〉〈/div〉

Description: 〈div data-abstract-type="normal"〉〈p〉Early-maturing provitamin A (PVA) quality protein maize (QPM) hybrids with combined drought and low soil nitrogen (low-N) tolerance are needed to address malnutrition and food security problems in sub-Saharan Africa (SSA). The current study's objectives were to (i) examine combining ability of selected early maturing PVA-QPM inbreds for grain yield and other agronomic traits under drought, low-N, optimal environments and across environments, (ii) determine gene action conditioning PVA accumulation under optimal environments, (iii) classify inbreds into heterotic groups and identify testers and (iv) assess yield and stability of hybrids across environments. Ninety-six hybrids generated from 24 inbred lines using the North Carolina Design II together with four commercial hybrid controls were evaluated under drought, low-N and optimal environments in Nigeria in 2016 and 2017. Fifty-four selected hybrids were assayed for PVA carotenoid and tryptophan content. Additive genetic effects were greater than non-additive effects for grain yield and most agronomic traits under each and across environments. The gene action conditioning accumulation of PVA carotenoids under optimal growing conditions followed a pattern similar to that of grain yield and other yield-related traits. The inbred lines were categorized into four heterotic groups consistent with the pedigree records and with TZEIORQ 29 identified as the best male and female tester for heterotic group IV. No tester was found for the other groups. Hybrid TZEIORQ 24 × TZEIORQ 41 was the highest yielding and most stable across environments and should be further tested for consistent performance for commercialization in SSA.〈/p〉〈/div〉

Description: 〈div data-abstract-type="normal"〉〈p〉〈span〉Hordeum marinum〈/span〉 commonly known as sea barley is a salinity-tolerant species of grass. In the current study, 150 lines from ten populations of 〈span〉H. marinum〈/span〉 ssp. 〈span〉marinum〈/span〉 collected from five Tunisian bioclimatic sites were screened for polymorphism with 13 selected random amplified polymorphic DNA primers. Results exhibited a high level of polymorphism (160 polymorphic bands with an average of 12.46 per primer) and a high level of genetic diversity in all the studied populations (on average 〈span〉UHe〈/span〉 = 0.247 and 〈span〉I〈/span〉 = 0.358). High discrimination capacity was found for the 13 primers and a combination of three allowed assignation of a unique profile for each of the 150 lines. The partition of genetic diversity with Analysis of Molecular Variance suggested that the majority of genetic variation (67%) was within populations. The components between-populations within ecoregions and between-ecoregions explained 21 and 12%, respectively, of the total genetic variance. There was no significant association of population differentiation (Ф〈span〉PT〈/span〉) with geographical distance or altitudinal difference. Results also showed that the 150 lines grouped into three clusters with no respect to geographic origin. A sub-set of 13 lines was identified, which captured the maximum genetic diversity of the entire collection. The genetic variation found in this collection of 〈span〉H. marinum〈/span〉 is deemed to be useful in formulating conservation strategies for this species.〈/p〉〈/div〉

Description: 〈div data-abstract-type="normal"〉〈p〉Potential planting area for tuber mustard was simulated using the Maxent model under current and future conditions based on 591 coordinates and 22 environmental layers. Model accuracy was excellent, with area under the receiving operator curve values of 0.967 and 0.958 for model training and testing, respectively. Dominant factors were mean diurnal range, mean temperature of the coldest quarter, annual mean temperature and minimum temperature of the coldest month, with thresholds of 6.5–7.5, 5.5–9, 16–19 and 2.0–6.5 °C, respectively. Under current conditions, suitable habitat areas (2.16% of total land in China) were concentrated mainly in Central, Southwest and East China, which can be defined as three occurrence and diffusion centres. In the 2050s and 2070s, suitable habitat areas are predicted to change to 3.72 and 3.92%, and 3.60 and 3.73% under scenarios RCP4.5 and RCP6.0, respectively, indicating that suitable habitat areas will increase slightly. However, future distribution of tuber mustard was predicted to differ among provinces or cities, i.e. predicted suitable habitat areas in Sichuan Province increased up to the 2050s but remained relatively unchanged between the 2050s and 2070s; in Chongqing city they first increased and then decreased; in Hunan, Anhui, Jiangsu, Zhejiang and Fujian Provinces they increased continuously; and in Guizhou, Hubei, Jiangxi Provinces and Shanghai city they first decreased, and then increased. The results from the current study provide useful information for management decisions of tuber mustard.〈/p〉〈/div〉

Description: 〈div data-abstract-type="normal"〉〈p〉A completely randomized experiment was designed to evaluate the effects of 〈span〉α〈/span〉-amylase (AMY) and glucoamylase (GLU) on total losses, fermentative profile, chemical composition and amylolytic activity of rehydrated maize. Eighty-four experimental silos of rehydrated maize [0.33 litres/kg ground maize, 4-mm theoretical particle size, and 625 g/kg dry matter (DM)] were assigned to the following treatments: (1) control (CON), no enzyme addition; (2) GLU added at 300 µl/kg of ground maize (as-fed); and (3) AMY added at 300 µl/kg of ground maize. Seven silos from each treatment were opened after 7, 14, 21 and 28 days. Differences among treatments were evaluated through orthogonal contrasts (CON 〈span〉v.〈/span〉 enzymes, and AMY 〈span〉v.〈/span〉 GLU). Time effects were decomposed using polynomial regression. Glucoamylase silage exhibited greater total losses than AMY. Enzymes increased acetate and lactic acid concentrations and decreased ethanol concentration. Regardless of treatment, gas, effluent and total fermentative losses linearly increased, whereas DM recovery linearly decreased with higher storage length. Glucoamylase silage had lower ammonia nitrogen and higher lactic acid concentrations than AMY. Enzyme treatments decreased silage neutral detergent fibre content and increased 〈span〉in vitro〈/span〉 DM degradation. Glucoamylase silage exhibited a more moderate starch content and greater 〈span〉in vitro〈/span〉 DM degradation than AMY. Storage time linearly decreased DM, starch and fibre content of rehydrated maize. 〈span〉In vitro〈/span〉 degradation of DM linearly increased as the storage length increased. This study showed evidence that enzymes with amylolytic activity, particularly GLU, improve the fermentative profile and DM degradation of rehydrated maize silage.〈/p〉〈/div〉

Description: 〈div data-abstract-type="normal"〉〈p〉We argue in this article that the social and economic conditions in the Angkorian society of the tenth century or earlier contributed to the decline in status of some middle-level officials, as is evident from the mid-eleventh century. Many Angkorian inscriptions written between the late ninth and late twelfth centuries record purchases and donations of lands acquired for religious foundations. The texts often contain details of transactions and disputes seeking to validate title to these holdings. The buyers include middle-ranking 〈span〉loñ〈/span〉 and 〈span〉vāp〈/span〉, and increasingly, higher-ranking officials. An analysis of the roles and activities of the officials reveals something of their relative status and helps explain the disappearance of 〈span〉vāp〈/span〉 from the inscriptions in the eleventh century, and the relegation of 〈span〉loñ〈/span〉 to temple roles by the twelfth century. The transfer of communal lands and lands owned by these officials to elites is attributed to hierarchical restrictions on land purchases, a reduction in fiscal immunities, and the need for taxes to be paid to the centre with high-value goods in Angkor's moneyless economy.〈/p〉〈/div〉

Description: 〈div data-abstract-type="normal"〉〈p〉The importance of spring grass in the diet of early lactation dairy cows has been widely acknowledged. Numerous studies completed on a plot/paddock basis have identified methods of increasing herbage availability in spring, but little focus was placed on how this impacts animal production. The aim of the current study was to investigate the impact of opening farm cover (OFC; grass availability on farm at turnout in spring; 1036 (high), 748 (medium) and 544 (low) kg DM/ha) and spring rotation length (fast – 56 days and normal – 63 days) on animal performance in early lactation and herbage production and quality. Spring rotation length had little effect on animal performance, herbage variables or sward composition. High OFC increased pre-grazing herbage mass, allowing for increased daily herbage allowance (DHA) compared to medium and low OFC. There was a reduced proportion of leaf in the sward of the high OFC compared to the low OFC, resulting in lower organic matter digestibility. Despite the reduction in sward quality observed as a consequence of achieving high OFC in spring, the greater DHA available to animals increased milk production (+1.4 kg milk/cow/day). Additionally, animals grazing a medium or low OFC had a greater requirement for silage supplementation in spring (+1.3 kg DM/cow/day). The benefits of the higher DHA highlighted in the current study suggest that autumn grazing management must be adapted to increase herbage availability in spring. However, the benefits observed in milk production did not persist beyond the first grazing rotation.〈/p〉〈/div〉

Description: 〈div data-abstract-type="normal"〉〈p〉Lord Reid played a vital role as chairman in a Commonwealth commission in framing the Malayan Independence constitution between 1956 and 1957. The Scottish Lord of Appeal sought to ensure the commission's impartiality and to achieve a fair balance between the demands of the various interest groups. The Federation of Malaya was a complex emerging nation-state with a diverse population and the framers had to manage competing interests and demands. This article, through a close examination of the primary constitutional documents, considers Reid's influence on the framing of the Malayan (and hence, Malaysian) federal constitution. The article begins with a brief discussion of Lord Reid's appointment to head the commission and then considers in some depth areas where his influence on the framing of the draft constitution is evident. The article argues that Reid was the main playmaker and moderator during the constitution-framing process and played a critical role in ensuring a balance was achieved between the competing demands of the federal government and the states, safeguarding the fundamental rights of the citizens against the state, and in moderating the various communal demands.〈/p〉〈/div〉