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.

Numerical and experimental investigations of three-dimensional container filling with Newtonian viscous fluids

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

Using urea dilution to standardise cellular and non-cellular components of pleural and bronchoalveolar lavage (BAL) fluids in the cat

A technique to standardise the analysis of cellular and non-cellular components in epithelial lining fluid (ELF) collected during saline lavage of pulmonary and pleural cavities was developed using the urea dilution method. Bronchoalveolar lavage (BAL) and pleural lavage (PL) fluids were collected from 12 clinically healthy cats. Total and differential cell counts in BAL fluid were within normal ranges for the cat, while cell Counts in PL fluid were assumed to be normal based on clinical health during examination, auscultation and lactate dehydrogenase (LDH) activities being comparable with other species. The major clinical implication of this study was that nucleated cell counts within feline ELF could not be predicted from analysis of lavage fluid which suggests that calculation of the proportion of ELF in lavage fluid by the urea dilution method may be necessary to avoid misdiagnosis of health or disease in pulmonary or pleural cavities. (C) 2005 ESFM and AAFP. Published by Elsevier Ltd. All rights reserved.

Conduction-radiation effect on natural convection flow in a fluid-saturated non-Darcy porous medium enclosed by non-isothermal walls

The effect of conduction-convection-radiation on natural convection flow of Newtonian optically thick gray fluid, confined in a non-Darcian porous media square cavity is numerically studied. For the gray fluid consideration is given to Rosseland diffusion approximation. Further assuming that (i) the temperature of the left vertical wall is varying linearly with height, (ii) cooled right vertical and top walls and (iii) the bottom wall is uniformly-heated. The governing equations are solved using the Alternate Direct Implicit method together with the Successive Over Relaxation technique. The investigation of the effect of governing parameters namely the Forschheimer resistance (Γ), the Planck constant (Rd), and the temperature difference (Δ), on flow pattern and heat transfer characteristics has been carried out. It was seen that the reduction of flow and heat transfer occurs as the Forschheimer resistance is increased. On the other hand both the strength of flow and heat transfer increases as the temperature ratio, Δ, is increased.

Delivering crushed paracetamol tablets using thickened fluids : does it alter the pharmacokinetics?

Dysphagia, often associated with conditions such as stroke, Parkinson’s disease, multiple sclerosis, and dementia, causes patients to have difficulty with swallowing food and/or liquids. These patients require their fluids to be thickened using gum-based thickening powders in order to facilitate safe swallowing. These thickened fluids are also used as a vehicle for delivery of crushed medicines. Our in vitro measurements suggest that thickened fluids can delay and reduce the dissolution of a number of medications. This study was conducted to assess the impact of the use of thickened fluids on the clinical pharmacokinetics of oral paracetamol. METHODS 20 Healthy volunteers were administered a single oral dose (1g) of paracetamol as either whole tablets, crushed with water, crushed with semi-solid jam, or crushed with thickened fluid according to a randomised, crossover design. Saliva samples were collected periodically over 8 hr and paracetamol concentration analysed by HPLC-UV. Non-compartmental pharmacokinetic analysis was conducted using Winnonlin®. RESULTS The mean peak concentration (Cmax) of paracetamol ranged between 5.62 – 8.00 μg/mL. Comparison between the crushed paracetamol with thickened water (Level 900) and other treatment options (whole, crushed with water, and crushed with jam) showed there was a significant difference in Cmax at 90% CI (p < 0.05). Also, whole tablet had a significant difference in Cmax between crushed with water and crushed with jam. There was no significant difference in AUC irrespective of the treatment. DISCUSSION The use of thickened water resulted in alteration in the absorption kinetics of paracetamol. Given this interaction, co-administration with thickened fluids may have important clinical implications for medications with a narrow therapeutic index.

Use of non-specific and specific interactions in the analysis of testosterone and related compounds by capillary electromigration techniques

Determination of testosterone and related compounds in body fluids is of utmost importance in doping control and the diagnosis of many diseases. Capillary electromigration techniques are a relatively new approach for steroid research. Owing to their electrical neutrality, however, separation of steroids by capillary electromigration techniques requires the use of charged electrolyte additives that interact with the steroids either specifically or non-specifically. The analysis of testosterone and related steroids by non-specific micellar electrokinetic chromatography (MEKC) was investigated in this study. The partial filling (PF) technique was employed, being suitable for detection by both ultraviolet spectrophotometry (UV) and electrospray ionization mass spectrometry (ESI-MS). Efficient, quantitative PF-MEKC UV methods for steroid standards were developed through the use of optimized pseudostationary phases comprising surfactants and cyclodextrins. PF-MEKC UV proved to be a more sensitive, efficient and repeatable method for the steroids than PF-MEKC ESI-MS. It was discovered that in PF-MEKC analyses of electrically neutral steroids, ESI-MS interfacing sets significant limitations not only on the chemistry affecting the ionization and detection processes, but also on the separation. The new PF-MEKC UV method was successfully employed in the determination of testosterone in male urine samples after microscale immunoaffinity solid-phase extraction (IA-SPE). The IA-SPE method, relying on specific interactions between testosterone and a recombinant anti-testosterone Fab fragment, is the first such method described for testosterone. Finally, new data for interactions between steroids and human and bovine serum albumins were obtained through the use of affinity capillary electrophoresis. A new algorithm for the calculation of association constants between proteins and neutral ligands is introduced.

Steady mixed convection stagnation-point flow of upper convected Maxwell fluids with magnetic field

The steady MHD mixed convection flow of a viscoelastic fluid in the vicinity of two-dimensional stagnation point with magnetic field has been investigated under the assumption that the fluid obeys the upper-convected Maxwell (UCM) model. Boundary layer theory is used to simplify the equations of motion. induced magnetic field and energy which results in three coupled non-linear ordinary differential equations which are well-posed. These equations have been solved by using finite difference method. The results indicate the reduction in the surface velocity gradient, surface heat transfer and displacement thickness with the increase in the elasticity number. These trends are opposite to those reported in the literature for a second-grade fluid. The surface velocity gradient and heat transfer are enhanced by the magnetic and buoyancy parameters. The surface heat transfer increases with the Prandtl number, but the surface velocity gradient decreases.

Stability of non-parabolic flow in a flexible tube

Flows with velocity profiles very different from the parabolic velocity profile can occur in the entrance region of a tube as well as in tubes with converging/diverging cross-sections. In this paper, asymptotic and numerical studies are undertaken to analyse the temporal stability of such 'non-parabolic' flows in a flexible tube in the limit of high Reynolds numbers. Two specific cases are considered: (i) developing flow in a flexible tube; (ii) flow in a slightly converging flexible tube. Though the mean velocity profile contains both axial and radial components, the flow is assumed to be locally parallel in the stability analysis. The fluid is Newtonian and incompressible, while the flexible wall is modelled as a viscoelastic solid. A high Reynolds number asymptotic analysis shows that the non-parabolic velocity profiles can become unstable in the inviscid limit. This inviscid instability is qualitatively different from that observed in previous studies on the stability of parabolic flow in a flexible tube, and from the instability of developing flow in a rigid tube. The results of the asymptotic analysis are extended numerically to the moderate Reynolds number regime. The numerical results reveal that the developing flow could be unstable at much lower Reynolds numbers than the parabolic flow, and hence this instability can be important in destabilizing the fluid flow through flexible tubes at moderate and high Reynolds number. For flow in a slightly converging tube, even small deviations from the parabolic profile are found to be sufficient for the present instability mechanism to be operative. The dominant non-parallel effects are incorporated using an asymptotic analysis, and this indicates that non-parallel effects do not significantly affect the neutral stability curves. The viscosity of the wall medium is found to have a stabilizing effect on this instability.

Effect of convection on domain growth during demixing transitions in fluids

Current analytical work on the effect of convection on the late stages of spinodal decomposition in liquids is briefly described. The morphology formed during the spinodal decomposition process depends on the relative composition of the two species. Droplet spinodal decomposition occurs when the concentration of one of the species is small. Convective transport has a significant effect on the scaling laws in the late-stage coarsening of droplets in translational or shear flows. In addition, convective transport could result in an attractive interaction between non-Brownian droplets which could lead to coalescence. The effect of convective transport for the growth of random interfaces in a near-symmetric quench was analysed using an area distribution function, which gives the distribution of surface area of the interface in curvature space. It was found that the curvature of the interface decreases proportional to time t in the late stages of spinodal decomposition, and the surface area also decreases proportional to t.

Explicit filtering LES for turbulent non-premixed combustion

Two models for large eddy simulation of turbulent reacting flow in homogeneous turbulence were studied. The sub-grid stress arising out of non-linearities of the Navier-Stokes equations were modeled using an explicit filtering approach. A filtered mass density function (FMDF) approach was used for closure of the sub-grid scalar fluctuations. A posteriori calculations, when compared with the results from the direct numerical simulation, indicate that the explicit filtering is adequate in representing the effect of sub-grid stress on the filtered velocity field in the absence of reaction. Discrepancies arise when reactions occur, but the FMDF approach suffices to account for sub-grid scale fluctuations of the reacting scalars, accurately.

Effect of convection and viscoelasticity on spinodal decomposition in fluids

Current analytical work on the effect of convection and viscoelasticity on the early and late stages of spinodal decomposition is briefly described. In the early stages, the effect of viscoelastic stresses was analysed using a simple Maxwell model for the stress, which was incorporated in the Langevin equation for the momentum field. The viscoelastic stresses are found to enhance the rate of decomposition. In the late stages, the pattern formed depends on the relative composition of the two species. Droplet spinodal decomposition occurs when the concentration of one of the species is small. Convective transport does not have a significant effect on the growth of a single droplet, but it does result in an attractive interaction between non - Brownian droplets which could lead to coalescence. The effect of convective transport for the growth of random interfaces in a near symmetric quench was analysed using an 'area distribution function', which gives the distribution of surface area of the interface in curvature space. It was found that the curvature of the interface decreases proportional to t in the late stages of spinodal decomposition, and the surface area also decreases proportional to t.

Numerical simulations of bubble formation from submerged needles under non-uniform direct current electric field

In several chemical and space industries, small bubbles are desired for efficient interaction between the liquid and gas phases. In the present study, we show that non-uniform electric field with appropriate electrode configurations can reduce the volume of the bubbles forming at submerged needles by up to three orders of magnitude. We show that localized high electric stresses at the base of the bubbles result in slipping of the contact line on the inner surface of the needle and subsequent bubble formation occurs with contact line inside the needle. We also show that for bubble formation in the presence of highly non-uniform electric field, due to high detachment frequency, the bubbles go through multiple coalescences and thus increase the apparent volume of the detached bubbles. (C) 2013 AIP Publishing LLC.

Effect of potassium present in stratum corneum during non-invasive measurement of potassium in human subjects using reverse iontophoresis

Purpose: Reverse iontophoresis (RI) is one of the potential techniques used to monitor the concentration of various analytes in body fluids non -invasively. Transdermal extraction of potassium is investigated using RI. In the present work, the effect of potassium on stratum corneum (SC) during RI, feasibility of RI for continuous monitoring of potassium, and use of potassium as internal standard in RI, are investigated. Methods: Tape stripping experiment is carried out to find potassium concentration in SC. RI is carried out continuously for 180 min without passive diffusion and after passive diffusion for 60 min. Skin impedance measurements are done at 20 Hz and 20 kHz. Results: Potassium is found to be in the range 300-650 nmol/cm(2) on SC by tape stripping experiment. Correlation coefficient between blood potassium and extracted potassium through RI after passive diffusion (R-2 = 0.5870) is more than without passive diffusion (R-2 = 0.5117). The skin impedance measurement shows that RI has more effect on SC than superficial layer of SC during RI. Conclusion: The present investigations conclude that it is possible to monitor potassium continuously through RI and using potassium as internal standard in RI.