Perturbational Finite Volume Method For The Solution of 2-D Navier-Stokes Equations on Unstructured and Structured Colocated Meshes

Based on the first-order upwind and second-order central type of finite volume( UFV and CFV) scheme, upwind and central type of perturbation finite volume ( UPFV and CPFV) schemes of the Navier-Stokes equations were developed. In PFV method, the mass fluxes of across the cell faces of the control volume (CV) were expanded into power series of the grid spacing and the coefficients of the power series were determined by means of the conservation equation itself. The UPFV and CPFV scheme respectively uses the same nodes and expressions as those of the normal first-order upwind and second-order central scheme, which is apt to programming. The results of numerical experiments about the flow in a lid-driven cavity and the problem of transport of a scalar quantity in a known velocity field show that compared to the first-order UFV and second-order CFV schemes, upwind PFV scheme is higher accuracy and resolution, especially better robustness. The numerical computation to flow in a lid-driven cavity shows that the under-relaxation factor can be arbitrarily selected ranging from 0.3 to 0. 8 and convergence perform excellent with Reynolds number variation from 102 to 104.

Numerical-Value Perturbation Method with Application of Solving Convective-Diffusion Integral Equation

The convective--diffusion equation is of primary importance in such fields as fluid dynamics and heat transfer hi the numerical methods solving the convective-diffusion equation, the finite volume method can use conveniently diversified grids (structured and unstructured grids) and is suitable for very complex geometry The disadvantage of FV methods compared to the finite difference method is that FV-methods of order higher than second are more difficult to develop in three-dimensional cases. The second-order central scheme (2cs) offers a good compromise among accuracy, simplicity and efficiency, however, it will produce oscillatory solutions when the grid Reynolds numbers are large and then very fine grids are required to obtain accurate solution. The simplest first-order upwind (IUW) scheme satisfies the convective boundedness criteria, however. Its numerical diffusion is large. The power-law scheme, QMCK and second-order upwind (2UW) schemes are also often used in some commercial codes. Their numerical accurate are roughly consistent with that of ZCS. Therefore, it is meaningful to offer higher-accurate three point FV scheme. In this paper, the numerical-value perturbational method suggested by Zhi Gao is used to develop an upwind and mixed FV scheme using any higher-order interpolation and second-order integration approximations, which is called perturbational finite volume (PFV) scheme. The PFV scheme uses the least nodes similar to the standard three-point schemes, namely, the number of the nodes needed equals to unity plus the face-number of the control volume. For instanc6, in the two-dimensional (2-D) case, only four nodes for the triangle grids and five nodes for the Cartesian grids are utilized, respectively. The PFV scheme is applied on a number of 1-D problems, 2~Dand 3-D flow model equations. Comparing with other standard three-point schemes, The PFV scheme has much smaller numerical diffusion than the first-order upwind (IUW) scheme, its numerical accuracy are also higher than the second-order central scheme (2CS), the power-law scheme (PLS), the QUICK scheme and the second-order upwind(ZUW) scheme.

Numerical simulation of a single fluid particle rising in non-Newtonian fluids

In this work, a level set method is developed for simulating the motion of a fluid particle rising in non-Newtonian fluids described by generalized Newtonian as well as viscoelastic model fluids. As the shear-thinning model we use a Carreau-Yasuda model, and the viscoelastic effect can be modeled with Oldroyd-B constitutive equations. The control volume formulation with the SIMPLEC algorithm incorporated is used to solve the governing equations on a staggered Eulerian grid. The level set method is implemented to compute the motion of a bubble in a Newtonian fluid as one of typical examples for validation, and the computational results are in good agreement with the reported experimental data．The level set method is also applied for simulating a Newtonian drop rising in Carreau-Yasuda and Oldroyd-B fluids．Numerical results including noticeably negative wake behind the drop and viscosity field are obtained, and compare satisfactorily with the known literature data.

Numerical simulation of a bubble rising in shear-thinning fluids

The motion of a single bubble rising freely in quiescent non-Newtonian viscous fluids was investigated experimentally and computationally. The non-Newtonian effects in the flow of viscous inelastic fluids are modeled by the Carreau theological model. An improved level set approach for computing the incompressible two-phase flow with deformable free interface is used. The control volume formulation with the SIMPLEC algorithm incorporated is used to solve the governing equations on a staggered Eulerian grid. The simulation results demonstrate that the algorithm is robust for shear-thinning liquids with large density (rho(1)/rho(g) up to 10(3)) and high viscosity (eta(1)/eta(g) up to 10(4)). The comparison of the experimental measurements of terminal bubble shape and velocity with the computational results is satisfactory. It is shown that the local change in viscosity around a bubble greatly depends on the bubble shape and the zero-shear viscosity of non-Newtonian shear-thinning liquids. The shear-rate distribution and velocity fields are used to elucidate the formation of a region of large viscosity at the rear of a bubble as a result of the rather stagnant flow behind the bubble. The numerical results provide the basis for further investigations, such as the numerical simulation of viscoelastic fluids. (C) 2010 Elsevier B.V. All rights reserved.

Design of spectrometer based on volume phase grating for near infrared range

In order to design and fabricate a spectrometer for the infrared range widely used in the different applications, Volume Phase Grating (VPG) with. low Polarization Dependence Loss (PDL) and high efficiency has been adopted as the dispersion element. VPG is constructed by coating an optical substrate with a thin film of dichromated. gelatin and exposing the film to two mutually coherent laser beams to form index modulation. The diffraction efficiency for a VPG is governed by Bragg effects. The depth (d) and index modulation contrast of the grating structure control the efficiency at which the light is diffracted when the Bragg condition is satisfied. Gradient index lens with high performance and low aberration are used as collimating system instead of standard lens. The spot diagrams and MTF curve of the collimating lens are shown in the paper. The receive system is InCaAs photodiode (PD) array including 512 pixels with 25 mum pitch. The spectrum resolution of the spectrometer reaches to 0.2nm and wavelength accuracy is 40pm.

Droplet Image Feedback Control System in Evaporation Experiment

In order to realize the steady-state droplet evaporation, image feedback control system is designed based on DSP. The system has three main functions: to capture and store droplet images during the experiment; to calculate droplet geometrical and physical parameters such as volume, surface area, surface tension and evaporation velocity at a high-precision level; to keep the droplet volume constant. The DSP can drive an injection controller with the PID control to inject liquid so as to keep the droplet volume constant. The evaporation velocity of droplet can be calculated by measuring the injected volume during the evaporation. The structure of hardware and software of the control system, key processing methods such as contour fitting and experimental results are described.

Control of self-organized low-dimensional morphology in Poly(styrene-b-4vinylpyridine)/polystyrene blend thin films

The surface morphologies of poly(styrene-b-4vinylpyridine) (PS-b-P4VP) diblock copolymer and homopolystyrene (hPS) binary blend thin films were investigated by atomic force microscopy as a function of total volume fraction of PS (phi(PS)) in the mixture. It was found that when hPS was added into symmetric PS-b-P4VP diblock copolymers, the surface morphology of this diblock copolymer was changed to a certain degree. With phi(PS) increasing at first, hPS was solubilized into the corresponding domains of block copolymer and formed cylinders. Moreover, the more solubilized the hPS, the more cylinders exist. However, when the limit was reached, excessive hPS tended to separate from the domains independently instead of solubilizing into the corresponding domains any longer, that is, a macrophase separation occurred. A model describing transitions of these morphologies with an increase in phi(PS) is proposed. The effect of composition on the phase morphology of blend films when graphite is used as a substrate is also investigated.