Soft-spring wall based non-periodic boundary conditions for non-equilibrium molecular dynamics of dense fluids

Non-equilibrium molecular dynamics (MD) simulations require imposition of non-periodic boundary conditions (NPBCs) that seamlessly account for the effect of the truncated bulk region on the simulated MD region. Standard implementation of specular boundary conditions in such simulations results in spurious density and force fluctuations near the domain boundary and is therefore inappropriate for coupled atomistic-continuum calculations. In this work, we present a novel NPBC model that relies on boundary atoms attached to a simple cubic lattice with soft springs to account for interactions from particles which would have been present in an untruncated full domain treatment. We show that the proposed model suppresses the unphysical fluctuations in the density to less than 1% of the mean while simultaneously eliminating spurious oscillations in both mean and boundary forces. The model allows for an effective coupling of atomistic and continuum solvers as demonstrated through multiscale simulation of boundary driven singular flow in a cavity. The geometric flexibility of the model enables straightforward extension to nonplanar complex domains without any adverse effects on dynamic properties such as the diffusion coefficient. (c) 2015 AIP Publishing LLC.

Third post-Newtonian gravitational waveforms for compact binary systems in general orbits: Instantaneous terms

We compute the instantaneous contributions to the spherical harmonic modes of gravitational waveforms from compact binary systems in general orbits up to the third post-Newtonian (PN) order. We further extend these results for compact binaries in quasielliptical orbits using the 3PN quasi-Keplerian representation of the conserved dynamics of compact binaries in eccentric orbits. Using the multipolar post-Minkowskian formalism, starting from the different mass and current-type multipole moments, we compute the spin-weighted spherical harmonic decomposition of the instantaneous part of the gravitational waveform. These are terms which are functions of the retarded time and do not depend on the history of the binary evolution. Together with the hereditary part, which depends on the binary's dynamical history, these waveforms form the basis for construction of accurate templates for the detection of gravitational wave signals from binaries moving in quasielliptical orbits.

Continuous synthesis of high quality CdSe quantum dots in supercritical fluids

We demonstrate in here a powerful scalable technology to synthesize continuously high quality CdSe quantum dots (QDs) in supercritical hexane. Using a low cost, highly thermally stable Cd-precursor, cadmium deoxycholate, the continuous synthesis is performed in 400 mu m ID stainless steel capillaries resulting in CdSe QDs having sharp full-width-at-half-maxima (23 nm) and high photoluminescence quantum yields (45-55%). Transmission electron microscopy images show narrow particles sizes distribution (sigma <= 5%) with well-defined crystal lattices. Using two different synthesis temperatures (250 degrees C and 310 degrees C), it was possible to obtain zinc blende and wurtzite crystal structures of CdSe QDs, respectively. This synthetic approach allows achieving substantial production rates up to 200 mg of QDs per hour depending on the targeted size, and could be easily scaled to gram per hour.

Effect of viscosity contrast on plume formation in density stratified fluids

We perform two and three dimensional numerical simulations of plume formation in density and viscosity stratified fluid systems. We show that the ambient to plume fluid viscosity ratio strongly affects the near wall plume structures (line or sheet plumes) such as plume spacing and shape of plumes. We observe that where mushroom-like plumes are observed for lower viscosity ratios, taller plumes with bulbous heads form for high viscosity ratios. Plume structure and spacing are in good agreement with experimental results. By studying the geometry of the line plumes and the flow in the circulation cells, we discuss the mechanisms of their formation and the dynamics of merging. We show that an increase in the viscosity ratio decreases the total length of line plumes in the planform which indicates a decreased mixing at higher viscosity ratios. (C) 2015 Elsevier Ltd. All rights reserved.

ASYMPTOTIC PRESERVING SCHEME FOR A KINETIC MODEL DESCRIBING IN COMPRESSIBLE FLUIDS

The kinetic theory of fluid turbulence modeling developed by Degond and Lemou in 7] is considered for further study, analysis and simulation. Starting with the Boltzmann like equation representation for turbulence modeling, a relaxation type collision term is introduced for isotropic turbulence. In order to describe some important turbulence phenomenology, the relaxation time incorporates a dependency on the turbulent microscopic energy and this makes difficult the construction of efficient numerical methods. To investigate this problem, we focus here on a multi-dimensional prototype model and first propose an appropriate change of frame that makes the numerical study simpler. Then, a numerical strategy to tackle the stiff relaxation source term is introduced in the spirit of Asymptotic Preserving Schemes. Numerical tests are performed in a one-dimensional framework on the basis of the developed strategy to confirm its efficiency.

Screened and unscreened phases in sedimenting suspensions

A coarse-grained stochastic hydrodynamical description of velocity and concentration fluctuations in steadily sedimenting suspensions is constructed and analyzed using self-consistent and renormalization-group methods. We find a nonequilibrium phase transition from an "unscreened" phase in which we recover the Caflisch-Luke [Phys. Fluids 28, 759 (1985)] divergence of the velocity variance to a "screened" phase where the fluctuations have a finite correlation length depending on the volume fraction phi as phi(-1/3), in agreement with Segre et al. [Phys. Rev. Lett. 79, 2574 (1997)] (if their observation of a phi-independent diffusivity is used), and the velocity variance is independent of system size.

Time-temperature histories of spherical food products during bulk forced-air precooling

A general mathematical model for forced air precooling of spherical food products in bulk is developed. The food products are arranged inline to form a rectangular parallelepiped. Chilled air is blown along the height of the package. The governing equations for the transient two-dimensional conduction with internal heat generation in the product, simultaneous heat and mass transfer at the product-air interface and one-dimensional transient energy and species conservation equations for the moist air are solved numerically using finite difference methods. Results are presented in the form of time-temperature histories. Experiments are conducted with model foods in a laboratory scale air precooling tunnel. The agreement between the theoretical and experimental results is found to be good. In general, a single product analysis fails to predict the precooling characteristics of bulk loads of food products. In the range of values investigated, the respiration heat is found to have a negligible effect.

On the relationship between finger width, velocity, and fluxes in thermohaline convection

Double-diffusive finger convection occurs in many natural processes.The theories for double-diffusive phenomena that exist at present consider systems with linear stratification in temperature and salinity. The double-diffusive systems with step change in salinity and temperature are, however, not amenable to simple stability analysis. Hence factors that control the width of the finger, velocity, and fluxes in systems that have step change in temperature and salinity have not been understood so far. In this paper we provide new physical insight regarding factors that influence finger convection in two-layer double-diffusive system through two-dimensional numerical simulations. Simulations have been carried out for density stability ratios (R-rho) from 1.5 to 10. For each density stability ratio, the thermal Rayleigh number (Ra-T) has been systematically varied from 7x10(3) to 7x10(8). Results from these simulations show how finger width, velocity, and flux ratios in finger convection are interrelated and the influence of governing parameters such as density stability ratio and the thermal Rayleigh number. The width of the incipient fingers at the time of onset of instability has been shown to vary as Ra-T-1/3. Velocity in the finger varies as Ra(T)1/3/R-rho. Results from simulation agree with the scale analysis presented in the paper. Our results demonstrate that wide fingers have lower velocities and flux ratios compared to those in narrow fingers. This result contradicts present notions about the relation between finger width and flux ratio. A counterflow heat-exchanger analogy is used in understanding the dependence of flux ratio on finger width and velocity.

Implicit and multigrid procedures for steady-state computations with upwind algorithms

A computational study for the convergence acceleration of Euler and Navier-Stokes computations with upwind schemes has been conducted in a unified framework. It involves the flux-vector splitting algorithms due to Steger-Warming and Van Leer, the flux-difference splitting algorithms due to Roe and Osher and the hybrid algorithms, AUSM (Advection Upstream Splitting Method) and HUS (Hybrid Upwind Splitting). Implicit time integration with line Gauss-Seidel relaxation and multigrid are among the procedures which have been systematically investigated on an individual as well as cumulative basis. The upwind schemes have been tested in various implicit-explicit operator combinations such that the optimal among them can be determined based on extensive computations for two-dimensional flows in subsonic, transonic, supersonic and hypersonic flow regimes. In this study, the performance of these implicit time-integration procedures has been systematically compared with those corresponding to a multigrid accelerated explicit Runge-Kutta method. It has been demonstrated that a multigrid method employed in conjunction with an implicit time-integration scheme yields distinctly superior convergence as compared to those associated with either of the acceleration procedures provided that effective smoothers, which have been identified in this investigation, are prescribed in the implicit operator.

Index properties of illite-bentonite mixtures in electrolyte solutions

Bentonite, commonly used for liner constructions in waste containment systems, possesses many limitations. Illite or illite containing bentonite has been proposed as an alternative material for liner construction. Their properties in different types of pore fluids are important to assess the long-term performance of the liner. Further, the illite-bentonite interaction occurs and changes their properties. The effect of these interactions is known when the pore fluid is only water. How their properties are modified in electrolyte solutions has been brought out in this paper. The index properties have been studied since they give an indication of their engineering properties. Due to reduction in the thickness of the diffused double layer and consequent particle aggregation in bentonite, the effect of clay-clay interaction reduces in electrolyte solutions. In electrolyte solutions, the liquid limit, the plasticity index, and free swell index of bentonite are lower than illite. The plasticity index of bentonite is further reduced in KCI solution. Clays with a higher plasticity index perform better to retain pollutants and reduce permeability. Hence, the presence of both illite and bentonite ensures better performance of the liner in different fluids.

Transitional intermittency in boundary layers subjected to pressure gradient

Results are reported from an extensive series of experiments on boundary layers in which the location of pressure gradient and transition onset could be varied almost independently, by judicious use of tunnel wall liners and transition-fixing devices. The experiments show that the transition zone is sensitive to the pressure gradient especially near onset, and can be significantly asymmetric; no universal similarity appears valid in general. Observed intermittency distributions cannot be explained on the basis of the hypothesis, often made, that the spot propagates at speeds proportional to the local free-stream velocity but is otherwise unaffected by the pressure gradient.

Microcontinuum approach to the pulsatile flow in tubes with and without longitudinal vibration

A fully developed pulsatile flow in a circular rigid tube is analysed by a microcontinuum approach. Solutions for radial variation of axial velocity and cell rotational velocity across the tube are obtained using the momentum integral method. Simplified forms of the solutions are presented for the relevant physiological data. Marked deviations in the results are observed when compared to a Newtonian fluid model. It is interesting to see that there is sufficient reduction in the mass flow rate, phase lag and friction due to the micropolar character of the fluid.

Effect of a magnetic field on the flow and blood oxygenation in channels of variable cross-section

The effect of a magnetic field on the flow and oxygenation of an incompressible Newtonian conducting fluid in channels with irregular boundaries has been investigated. The geometric parameter δ, which is a ratio of the mean half width of the channel d to the characteristic length λ along the channel over which the significant changes in the flow quantities occur, has been used for perturbing the governing equations. Closed form solutions of the various order equations are presented for the stream function. The equations for oxygen partial pressure remain nonlinear even after perturbation, therefore a numerical solution is presented. The expressions for shear stress at a wall and pressure distributions are derived. Here the separation in the flow occurs at a higher Reynolds number than the corresponding non-magnetic case. It is found that the magnetic field has an effect on local oxygen concentration but has a little effect on the saturation length.

Flow of micropolar fluid through a tube with stenosis

mathematical model for the steady flow of non-Newtonian fluid through a stenotic region is presented. The results indicate that the general shape and size of the stenosis together with rheological properties of blood are important in understanding the flow characteristics and the presence of flow separation.

Instabilities of a compressible stratified fluid in horizontal sheared motion

The stability characteristics of a Helmholtz velocity profile in a stably stratified, compressible atmosphere in the presence of a lower boundary are studied. A jump in the Brunt–Väisälä frequency is introduced and the level at which this jump occurs is assumed to be different from the shear zone, to simulate sharp temperature discontinuities in the atmosphere. The results are compared with those of Pellacani, Tebaldi, and Tosi and Lindzen and Rosenthal. In the present configuration, new unstable modes with larger growth rates are found. The wavelengths of the most unstable gravity waves for the parameters pertaining to observed cases of clear air turbulence agree quite closely with the experimental values. Physics of Fluids is copyrighted by The American Institute of Physics