Viscous coupling of shear-free turbulence across nearly flat fluid interfaces

The interactions between shear-free turbulence in two regions (denoted as + and − on either side of a nearly flat horizontal interface are shown here to be controlled by several mechanisms, which depend on the magnitudes of the ratios of the densities, ρ+/ρ−, and kinematic viscosities of the fluids, μ+/μ−, and the root mean square (r.m.s.) velocities of the turbulence, u0+/u0−, above and below the interface. This study focuses on gas–liquid interfaces so that ρ+/ρ− ≪ 1 and also on where turbulence is generated either above or below the interface so that u0+/u0− is either very large or very small. It is assumed that vertical buoyancy forces across the interface are much larger than internal forces so that the interface is nearly flat, and coupling between turbulence on either side of the interface is determined by viscous stresses. A formal linearized rapid-distortion analysis with viscous effects is developed by extending the previous study by Hunt & Graham (J. Fluid Mech., vol. 84, 1978, pp. 209–235) of shear-free turbulence near rigid plane boundaries. The physical processes accounted for in our model include both the blocking effect of the interface on normal components of the turbulence and the viscous coupling of the horizontal field across thin interfacial viscous boundary layers. The horizontal divergence in the perturbation velocity field in the viscous layer drives weak inviscid irrotational velocity fluctuations outside the viscous boundary layers in a mechanism analogous to Ekman pumping. The analysis shows the following. (i) The blocking effects are similar to those near rigid boundaries on each side of the interface, but through the action of the thin viscous layers above and below the interface, the horizontal and vertical velocity components differ from those near a rigid surface and are correlated or anti-correlated respectively. (ii) Because of the growth of the viscous layers on either side of the interface, the ratio uI/u0, where uI is the r.m.s. of the interfacial velocity fluctuations and u0 the r.m.s. of the homogeneous turbulence far from the interface, does not vary with time. If the turbulence is driven in the lower layer with ρ+/ρ− ≪ 1 and u0+/u0− ≪ 1, then uI/u0− ~ 1 when Re (=u0−L−/ν−) ≫ 1 and R = (ρ−/ρ+)(v−/v+)1/2 ≫ 1. If the turbulence is driven in the upper layer with ρ+/ρ− ≪ 1 and u0+/u0− ≫ 1, then uI/u0+ ~ 1/(1 + R). (iii) Nonlinear effects become significant over periods greater than Lagrangian time scales. When turbulence is generated in the lower layer, and the Reynolds number is high enough, motions in the upper viscous layer are turbulent. The horizontal vorticity tends to decrease, and the vertical vorticity of the eddies dominates their asymptotic structure. When turbulence is generated in the upper layer, and the Reynolds number is less than about 106–107, the fluctuations in the viscous layer do not become turbulent. Nonlinear processes at the interface increase the ratio uI/u0+ for sheared or shear-free turbulence in the gas above its linear value of uI/u0+ ~ 1/(1 + R) to (ρ+/ρ−)1/2 ~ 1/30 for air–water interfaces. This estimate agrees with the direct numerical simulation results from Lombardi, De Angelis & Bannerjee (Phys. Fluids, vol. 8, no. 6, 1996, pp. 1643–1665). Because the linear viscous–inertial coupling mechanism is still significant, the eddy motions on either side of the interface have a similar horizontal structure, although their vertical structure differs.

The influence of Langmuir turbulence on the scaling for the dissipation rate in the oceanic boundary layer

A model for estimating the turbulent kinetic energy dissipation rate in the oceanic boundary layer, based on insights from rapid-distortion theory, is presented and tested. This model provides a possible explanation for the very high dissipation levels found by numerous authors near the surface. It is conceived that turbulence, injected into the water by breaking waves, is subsequently amplified due to its distortion by the mean shear of the wind-induced current and straining by the Stokes drift of surface waves. The partition of the turbulent shear stress into a shear-induced part and a wave-induced part is taken into account. In this picture, dissipation enhancement results from the same mechanism responsible for Langmuir circulations. Apart from a dimensionless depth and an eddy turn-over time, the dimensionless dissipation rate depends on the wave slope and wave age, which may be encapsulated in the turbulent Langmuir number La_t. For large La_t, or any Lat but large depth, the dissipation rate tends to the usual surface layer scaling, whereas when Lat is small, it is strongly enhanced near the surface, growing asymptotically as ɛ ∝ La_t^{-2} when La_t → 0. Results from this model are compared with observations from the WAVES and SWADE data sets, assuming that this is the dominant dissipation mechanism acting in the ocean surface layer and statistical measures of the corresponding fit indicate a substantial improvement over previous theoretical models. Comparisons are also carried out against more recent measurements, showing good order-of-magnitude agreement, even when shallow-water effects are important.

Dissipation of shear-free turbulence near boundaries

The rapid-distortion model of Hunt & Graham (1978) for the initial distortion of turbulence by a flat boundary is extended to account fully for viscous processes. Two types of boundary are considered: a solid wall and a free surface. The model is shown to be formally valid provided two conditions are satisfied. The first condition is that time is short compared with the decorrelation time of the energy-containing eddies, so that nonlinear processes can be neglected. The second condition is that the viscous layer near the boundary, where tangential motions adjust to the boundary condition, is thin compared with the scales of the smallest eddies. The viscous layer can then be treated using thin-boundary-layer methods. Given these conditions, the distorted turbulence near the boundary is related to the undistorted turbulence, and thence profiles of turbulence dissipation rate near the two types of boundary are calculated and shown to agree extremely well with profiles obtained by Perot & Moin (1993) by direct numerical simulation. The dissipation rates are higher near a solid wall than in the bulk of the flow because the no-slip boundary condition leads to large velocity gradients across the viscous layer. In contrast, the weaker constraint of no stress at a free surface leads to the dissipation rate close to a free surface actually being smaller than in the bulk of the flow. This explains why tangential velocity fluctuations parallel to a free surface are so large. In addition we show that it is the adjustment of the large energy-containing eddies across the viscous layer that controls the dissipation rate, which explains why rapid-distortion theory can give quantitatively accurate values for the dissipation rate. We also find that the dissipation rate obtained from the model evaluated at the time when the model is expected to fail actually yields useful estimates of the dissipation obtained from the direct numerical simulation at times when the nonlinear processes are significant. We conclude that the main role of nonlinear processes is to arrest growth by linear processes of the viscous layer after about one large-eddy turnover time.

Effects of agricultural management practices on earthworm populations and crop yield: validation and application of a mechanistic modelling approach

There is little consensus on how agriculture will meet future food demands sustainably. Soils and their biota play a crucial role by mediating ecosystem services that support agricultural productivity. However, a multitude of site-specific environmental factors and management practices interact to affect the ability of soil biota to perform vital functions, confounding the interpretation of results from experimental approaches. Insights can be gained through models, which integrate the physiological, biological and ecological mechanisms underpinning soil functions. We present a powerful modelling approach for predicting how agricultural management practices (pesticide applications and tillage) affect soil functioning through earthworm populations. By combining energy budgets and individual-based simulation models, and integrating key behavioural and ecological drivers, we accurately predict population responses to pesticide applications in different climatic conditions. We use the model to analyse the ecological consequences of different weed management practices. Our results demonstrate that an important link between agricultural management (herbicide applications and zero, reduced and conventional tillage) and earthworms is the maintenance of soil organic matter (SOM). We show how zero and reduced tillage practices can increase crop yields while preserving natural ecosystem functions. This demonstrates how management practices which aim to sustain agricultural productivity should account for their effects on earthworm populations, as their proliferation stimulates agricultural productivity. Synthesis and applications. Our results indicate that conventional tillage practices have longer term effects on soil biota than pesticide control, if the pesticide has a short dissipation time. The risk of earthworm populations becoming exposed to toxic pesticides will be reduced under dry soil conditions. Similarly, an increase in soil organic matter could increase the recovery rate of earthworm populations. However, effects are not necessarily additive and the impact of different management practices on earthworms depends on their timing and the prevailing environmental conditions. Our model can be used to determine which combinations of crop management practices and climatic conditions pose least overall risk to earthworm populations. Linking our model mechanistically to crop yield models would aid the optimization of crop management systems by exploring the trade-off between different ecosystem services.

Energy consistent discontinuous Galerkin methods for the Navier–Stokes–Korteweg system

We design consistent discontinuous Galerkin finite element schemes for the approximation of the Euler-Korteweg and the Navier-Stokes-Korteweg systems. We show that the scheme for the Euler-Korteweg system is energy and mass conservative and that the scheme for the Navier-Stokes-Korteweg system is mass conservative and monotonically energy dissipative. In this case the dissipation is isolated to viscous effects, that is, there is no numerical dissipation. In this sense the methods are consistent with the energy dissipation of the continuous PDE systems. - See more at: http://www.ams.org/journals/mcom/2014-83-289/S0025-5718-2014-02792-0/home.html#sthash.rwTIhNWi.dpuf

The non-conservation of potential vorticity by a dynamical core compared with the effects of parametrized physical processes

Numerical models of the atmosphere combine a dynamical core, which approximates solutions to the adiabatic, frictionless governing equations for fluid dynamics, with tendencies arising from the parametrization of other physical processes. Since potential vorticity (PV) is conserved following fluid flow in adiabatic, frictionless circumstances, it is possible to isolate the effects of non-conservative processes by accumulating PV changes in an air-mass relative framework. This “PV tracer technique” is used to accumulate separately the effects on PV of each of the different non-conservative processes represented in a numerical model of the atmosphere. Dynamical cores are not exactly conservative because they introduce, explicitly or implicitly, some level of dissipation and adjustment of prognostic model variables which acts to modify PV. Here, the PV tracers technique is extended to diagnose the cumulative effect of the non-conservation of PV by a dynamical core and its characteristics relative to the PV modification by parametrized physical processes. Quantification using the Met Office Unified Model reveals that the magnitude of the non-conservation of PV by the dynamical core is comparable to those from physical processes. Moreover, the residual of the PV budget, when tracing the effects of the dynamical core and physical processes, is at least an order of magnitude smaller than the PV tracers associated with the most active physical processes. The implication of this work is that the non-conservation of PV by a dynamical core can be assessed in case studies with a full suite of physics parametrizations and directly compared with the PV modification by parametrized physical processes. The nonconservation of PV by the dynamical core is shown to move the position of the extratropical tropopause while the parametrized physical processes have a lesser effect at the tropopause level.

Non-LTE Monte Carlo radiative transfer. II. Nonisothermal solutions for viscous Keplerian disks

We discuss the basic hydrodynamics that determines the density structure of the disks around hot stars. Observational evidence supports the idea that these disks are Keplerian (rotationally supported) gaseous disks. A popular scenario in the literature, which naturally leads to the formation of Keplerian disks, is the viscous decretion model. According to this scenario, the disks are hydrostatically supported in the vertical direction, while the radial structure is governed by the viscous transport. This suggests that the temperature is one primary factor that governs the disk density structure. In a previous study we demonstrated, using three-dimensional non-LTE Monte Carlo simulations, that viscous Keplerian disks can be highly nonisothermal. In this paper we build on our previous work and solve the full problem of the steady state nonisothermal viscous diffusion and vertical hydrostatic equilibrium. We find that the self-consistent solution departs significantly from the analytic isothermal density, with potentially large effects on the emergent spectrum. This implies that nonisothermal disk models must be used for a detailed modeling of Be star disks.

Effects of N supply on the accumulation of photosynthetic pigments and photoprotectors in Gracilaria tenuistipitata (Rhodophyta) cultured under UV radiation

We have studied the effects of nitrate supply under photosynthetic active radiation (PAR) plus ultraviolet radiation (UVR) exposure on photosynthetic pigments (chlorophyll a and carotenoids), photoprotective UV screen mycosporine-like amino acids (MAAs), and photosynthetic parameters, including the maximum quantum yield (F(v)/F(m)) and electron transport rate (ETR) on the red agarophyte Gracilaria tenuistipitata. Apical tips of G. tenuistipitata were cultivated under ten different concentrations of NO(3)(-) for 7 days. It has been shown that G. tenuistipitata cultured under laboratory conditions has the ability to accumulate high amounts of MAAs following a nitrate concentration-dependent manner under PAR+UVR. Two MAAs were identified, shinorine and porphyra-334. The relative concentration of the first increased under high concentrations of nitrate, while the second one decreased. The presence of antheraxanthin is reported for the first time in this macro-algae, which also contains zeaxanthin, lutein, and beta-carotene. The accumulation of pigments, photoprotective compounds, and photosynthetic parameters of G. tenuistipitata is directly related to N availability. All variables decreased under low N supplies and reached constant maximum values with supplements higher than 0.5 mM NO(3)(-). Our results suggest a high potential to acclimation and photoprotection against stress factors (including high PAR and UVR) directly related to N availability for G. tenuistipitata.

Effects of starch gelatinization and oxidation on the rheological behavior of chitosan/starch blends

Chitosan/starchblends represent an interesting alternative for the preparation of biocompatible drug delivery systems, packing materials and edible films. This paper reports on the effects of starch gelatinization and oxidation on the rheological behavior of chitosan/starch blends. The results show that the modifications in the starch structure cause changes in G` (storage modulus) and G `` (lossmodulus) as a function of frequency. For chitosan/starch, G `` is higher than G`, showing a viscous behavior. However, for chitosan/gelatinized starch and chitosan/oxidized starch, an increase in the angular frequency promotes a modulus crossover at omega = 0.02 and 0.04 rad s(-1), respectively. The viscosity curves as a function of shear rate show that both modifications cause an increase in viscosity, and all blends show a non-Newtonian behavior. (C) 2011 Society of Chemical Industry

Effects of crystallization conditions on dielectric and ferroelectric properties of PZT thin films

This paper reports studies on dielectric and ferroelectric properties of lead zirconate titanate (PZT) thin films crystallized by conventional thermal annealing (CTA) and rapid thermal annealing (RTA) in air, oxygen and nitrogen atmospheres to better understand, control and optimize these properties. The dielectric constant (epsilon) and dissipation factor (tan delta) values, at a frequency of 100 kHz; for film crystallized in air by CTA process, were 358 and 0.039, respectively. Considering the same frequency for film crystallized in air by RTA, these values were 611 and 0.026, respectively. The different dielectric values were justified by a space-charge or interfacial polarization in films, often characterized as Maxwell-Wagner type. This effect was also responsible to dispersion at frequencies above 1 MHz in film crystallized in air by CTA process and film crystallized by RTA in oxygen atmosphere. The film crystallized by RTA under nitrogen atmosphere presented an evident dispersion at frequencies around 100 Hz, characterized by an increase in both epsilon and tan delta. This dispersion was attributed to conductivity effects. The remanent polarization (P-r) and coercive field (E-c) were also obtained for all films. Films obtained from RTA in air presented higher P-r (17.8 muC cm(-2)) than film crystallized from CTA (7.8 muC cm(-2)). As a function of the crystallization atmospheres, films crystallized by RTA in air and nitrogen presented essentially the same P-r values (around 18 muC cm(-2)) but the P-r (3.9 muC cm(-2)) obtained from film crystallized under oxygen atmosphere was profoundly influenced.

PRODUCTION of THERMAL PHOTONS IN VISCOUS FLUID DYNAMICS WITH TEMPERATURE-DEPENDENT SHEAR VISCOSITY

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Dissipation and its consequences in the scaling exponents for a family of two-dimensional mappings

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Effects of milling condition on the surface integrity of hot forged steel

This paper presents a study on the influence of milling condition on workpiece surface integrity focusing on hardness and roughness. The experimental work was carried out on a CNC machining center considering roughing and finishing operations. A 25 mm diameter endmill with two cemented carbide inserts coated with TiN layer were used for end milling operation. Low carbon alloyed steel Cr-Mo forged at 1200 degrees C was used as workpiece on the tests. Two kinds of workpiece conditions were considered, i.e. cur cooled after hot forging and normalized at 950 degrees C for 2 h. The results showed that finishing operation was able to significantly decrease the roughness by at least 46% without changing the hardness. on the other hand, roughing operation caused an increase in hardness statistically significant by about 6%. The machined surface presented deformed regions within feed marks, which directly affected the roughness. Surface finish behavior seems to correlate to the chip ratio given the decrease of 25% for roughing condition, which damaged the chip formation. The material removal rate for finishing operation 41% greater than roughing condition demonstrated to be favorable to the heat dissipation and minimized the effect on material hardness.

Effects of the Compaction Pressure on the Magnetic Properties of a Sintered Fe-Based Alloy

The structural and magnetic properties of a Fe-based alloy before and after sintering have been analyzed. X ray diffraction measurements confirm the deformation of the magnetic particles in the compacted samples. After sintering, hysteresis energy dissipation, remanence and intrinsic coercivity differ by less than 10% as porosity changes from 15 to 7%.

Effects of photobioreactor configuration, nitrogen source and light intensity on the fed-batch cultivation of Arthrospira (Spirulina) platensis. Bioenergetic aspects

Bioenergetic analysis may be applied in order to predict microbial growth yields, based on the Gibbs energy dissipation and mass conservation principles of the overall growth reaction. The bioenergetics of the photoautotrophic growth of the cyanobacterium Arthrospira (Spirulina) platensis was investigated in different bioreactor configurations (tubular photobioreactor and open ponds) using different nitrogen sources (nitrate and urea) and under different light intensity conditions to determine the best growing conditions in terms of Gibbs energy dissipation, number of photons to sustain cell growth and phototrophic energy yields distribution in relation to the ATP and NADPH formation, and release of heat. Although an increase in the light intensity increased the Gibbs energy dissipated for cell growth and maintenance with both nitrogen sources, it did not exert any appreciable influence on the moles of photons absorbed by the system to produce one C-mol biomass. On the other hand, both bioenergetic parameters were higher in cultures with nitrate than with urea, likely because of the higher energy requirements needed to reduce the former nitrogen source to ammonia. They appreciably increased also when open ponds were substituted by the tubular photobioreactor, where a more efficient light distribution ensured a remarkably higher cell mass concentration. The estimated percentages of the energy absorbed by the cell showed that, compared with nitrate, the use of urea as nitrogen source allowed the system to address higher energy fractions to ATP production and light fixation by the photosynthetic apparatus, as well as a lower fraction released as heat. The best energy yields values on Gibbs energy necessary for cell growth and maintenance were achieved in up to 4-5 days of cultivation, indicating that it would be the optimum range to maintain cell growth. Thanks to this better bioenergetic situation, urea appears to be a quite promising low-cost, alternative nitrogen source for Arthrospira platensis cultures in photobioreactors. (C) 2011 Elsevier Ltd. All rights reserved.