Two-dimensional homogeneous isotropic fluid turbulence with polymer additives

We carry out an extensive and high-resolution direct numerical simulation of homogeneous, isotropic turbulence in two-dimensional fluid films with air-drag-induced friction and with polymer additives. Our study reveals that the polymers (a) reduce the total fluid energy, enstrophy, and palinstrophy; (b) modify the fluid energy spectrum in both inverse-and forward-cascade regimes; (c) reduce small-scale intermittency; (d) suppress regions of high vorticity and strain rate; and (e) stretch in strain-dominated regions. We compare our results with earlier experimental studies and propose new experiments.

CROSS-SECTIONAL ANALYSIS OF PRE-TWISTED THICK BEAMS USING VARIATIONAL ASYMPTOTIC METHOD

An asymptotically-exact methodology is presented for obtaining the cross-sectional stiffness matrix of a pre-twisted moderately-thick beam having rectangular cross sections and made of transversely isotropic materials. The anisotropic beam is modeled from 3-D elasticity, without any further assumptions. The beam is allowed to have large displacements and rotations, but small strain is assumed. The strain energy of the beam is computed making use of the constitutive law and the kinematical relations derived with the inclusion of geometrical nonlinearities and initial twist. Large displacements and rotations are allowed, but small strain is assumed. The Variational Asymptotic Method is used to minimize the energy functional, thereby reducing the cross section to a point on the reference line with appropriate properties, yielding a 1-D constitutive law. In this method as applied herein, the 2-D cross-sectional analysis is performed asymptotically by taking advantage of a material small parameter and two geometric small parameters. 3-D strain components are derived using kinematics and arranged as orders of the small parameters. Warping functions are obtained by the minimization of strain energy subject to certain set of constraints that renders the 1-D strain measures well-defined. Closed-form expressions are derived for the 3-D non-linear warping and stress fields. The model is capable of predicting interlaminar and transverse shear stresses accurately up to first order.

Stability of Poiseuille flow in a fluid overlying an anisotropic and inhomogeneous porous layer

We present the linear stability analysis of horizontal Poiseuille flow in a fluid overlying a porous medium with anisotropic and inhomogeneous permeability. The generalized Darcy model is used to describe the flow in the porous medium with the Beavers-Joseph condition at the interface of the two layers and the eigenvalue problem is solved numerically. The effect of major system parameters on the stability characteristics is addressed in detail. It is shown that the anisotropic and inhomogeneous modulation of the permeability of the underlying porous layer provides an effective means for passive control of the flow stability.

Linking monazite geochronology with fluid infiltration and metamorphic histories: Nature and experiment

Migmatised metapelites from the Kodaikanal region, central Madurai Block, southern India have undergone ultrahigh-temperature metamorphism (950-1000 degrees C; 7-8 kbar). In-situ electron microprobe Th-U-Pb isochron (CHIME) dating of monazites in a leucosome and surrounding silica-saturated and silica-poor restites from the same outcrop indicates three principal ages that can be linked to the evolutionary history of these rocks. Monazite grains from the silica-saturated restite have well-defined, inherited cores with thick rims that yield an age of ca. 1684 Ma. This either dates the metamorphism of the original metapelite or is a detrital age of inherited monazite. Monazite grains from the silica-poor restite, thick rims from the silica-saturated restite, and monazite cores from the leucosome have ages ranging from 520 to 540 Ma suggesting a mean age of 530 Ma within the error bars. In the leucosome the altered rim of the monazite gives an age of ca. 502 Ma. Alteration takes the form of Th-depleted lobes of monazite with sharp curvilinear boundaries extending from the monazite grain rim into the core. We have replicated experimentally these altered rims in a monazite-leucosome experiment at 800 degrees C and 2 kbar. This experiment, coupled with earlier published monazite-fluid experiments involving high pH alkali-bearing fluids at high P-T, helps to confirm the idea that alkali-bearing fluids, in the melt and along grain boundaries during crystallization, were responsible for the formation of the altered monazite grain rims via the process of coupled dissolution-reprecipitation. (C) 2015 Elsevier B.V. All rights reserved.

Effect of dark matter halo on global spiral modes in galaxies

Low surface brightness (LSB) galaxies form a major class of galaxies, and are characterized by low disc surface density and low star formation rate. These are known to be dominated by dark matter halo from the innermost regions. Here, we study the role of the dark matter halo on the grand-design, m = 2, spiral modes in a galactic disc by carrying out a global mode analysis in the WKB approximation. The Bohr-Sommerfeld quantization rule is used to determine how many discrete global spiral modes are permitted. First, a typical superthin, LSB galaxy UGC 7321 is studied by taking only the galactic disc, modelled as a fluid; and then the disc embedded in a dark matter halo. We find that both cases permit the existence of global spiral modes. This is in contrast to earlier results where the inclusion of dark matter halo was shown to nearly fully suppress local, swing-amplified spiral features. Although technically global modes are permitted in the fluid model as shown here, we argue that due to lack of tidal interactions, these are not triggered in LSB galaxies. For comparison, we carried out a similar analysis for the Galaxy, for which the dark matter halo does not dominate in the inner regions. We show that here too the dark matter halo has little effect, hence the disc embedded in a halo is also able to support global modes. The derived pattern speed of the global mode agrees fairly well with the observed value for the Galaxy.

Structure-rheology relationship in a sheared lamellar fluid

The structure-rheology relationship in the shear alignment of a lamellar fluid is studied using a mesoscale model which provides access to the lamellar configurations and the rheology. Based on the equations and free energy functional, the complete set of dimensionless groups that characterize the system are the Reynolds number (rho gamma L-2/mu), the Schmidt number (mu/rho D), the Ericksen number (mu(gamma)/B), the interface sharpness parameter r, the ratio of the viscosities of the hydrophilic and hydrophobic parts mu(r), and the ratio of the system size and layer spacing (L/lambda). Here, rho and mu are the fluid density and average viscosity, (gamma) over dot is the applied strain rate, D is the coefficient of diffusion, B is the compression modulus, mu(r) is the maximum difference in the viscosity of the hydrophilic and hydrophobic parts divided by the average viscosity, and L is the system size in the cross-stream direction. The lattice Boltzmann method is used to solve the concentration and momentum equations for a two dimensional system of moderate size (L/lambda = 32) and for a low Reynolds number, and the other parameters are systematically varied to examine the qualitative features of the structure and viscosity evolution in different regimes. At low Schmidt numbers where mass diffusion is faster than momentum diffusion, there is fast local formation of randomly aligned domains with ``grain boundaries,'' which are rotated by the shear flow to align along the extensional axis as time increases. This configuration offers a high resistance to flow, and the layers do not align in the flow direction even after 1000 strain units, resulting in a viscosity higher than that for an aligned lamellar phase. At high Schmidt numbers where momentum diffusion is fast, the shear flow disrupts layers before they are fully formed by diffusion, and alignment takes place by the breakage and reformation of layers by shear, resulting in defects (edge dislocations) embedded in a background of nearly aligned layers. At high Ericksen number where the viscous forces are large compared to the restoring forces due to layer compression and bending, shear tends to homogenize the concentration field, and the viscosity decreases significantly. At very high Ericksen number, shear even disrupts the layering of the lamellar phase. At low Ericksen number, shear results in the formation of well aligned layers with edge dislocations. However, these edge dislocations take a long time to anneal; the relatively small misalignment due to the defects results in a large increase in viscosity due to high layer stiffness and due to shear localization, because the layers between defects get pinned and move as a plug with no shear. An increase in the viscosity contrast between the hydrophilic and hydrophobic parts does not alter the structural characteristics during alignment. However, there is a significant increase in the viscosity, due to pinning of the layers between defects, which results in a plug flow between defects and a localization of the shear to a part of the domain.

A systematic asymptotic approach to determine the dispersion characteristics of structural-acoustic waveguides with arbitrary fluid loading

Structural-acoustic waveguides of two different geometries are considered: a 2-D rectangular and a circular cylindrical geometry. The objective is to obtain asymptotic expansions of the fluid-structure coupled wavenumbers. The required asymptotic parameters are derived in a systematic way, in contrast to the usual intuitive methods used in such problems. The systematic way involves analyzing the phase change of a wave incident on a single boundary of the waveguide. Then, the coupled wavenumber expansions are derived using these asymptotic parameters. The phase change is also used to qualitatively demarcate the dispersion diagram as dominantly structure-originated, fluid originated or fully coupled. In contrast to intuitively obtained asymptotic parameters, this approach does not involve any restriction on the material and geometry of the structure. The derived closed-form solutions are compared with the numerical solutions and a good match is obtained. (C) 2016 Elsevier Ltd. All rights reserved.

COUPLED OPTICAL-THERMAL-FLUID MODELING OF A DIRECTLY HEATED TUBULAR SOLAR RECEIVER FOR SUPERCRITICAL CO2 BRAYTON CYCLE

Recent studies have evaluated closed-loop supercritical carbon dioxide (s-CO2) Brayton cycles to be a higher energy density system in comparison to conventional superheated steam Rankine systems. At turbine inlet conditions of 923K and 25 MPa, high thermal efficiency (similar to 50%) can be achieved. Achieving these high efficiencies will make concentrating solar power (CSP) technologies a competitive alternative to current power generation methods. To incorporate a s-CO2 Brayton power cycle in a solar power tower system, the development of a solar receiver capable of providing an outlet temperature of 923 K (at 25 MPa) is necessary. The s-CO2 will need to increase in temperature by similar to 200 K as it passes through the solar receiver to satisfy the temperature requirements of a s-CO2 Brayton cycle with recuperation and recompression. In this study, an optical-thermal-fluid model was developed to design and evaluate a tubular receiver that will receive a heat input similar to 2 MWth from a heliostat field. The ray-tracing tool SolTrace was used to obtain the heat-flux distribution on the surfaces of the receiver. Computational fluid dynamics (CFD) modeling using the Discrete Ordinates (DO) radiation model was used to predict the temperature distribution and the resulting receiver efficiency. The effect of flow parameters, receiver geometry and radiation absorption by s-CO2 were studied. The receiver surface temperatures were found to be within the safe operational limit while exhibiting a receiver efficiency of similar to 85%.