Finite volume TVD formulation of lattice Boltzmann simulation on unstructured mesh

A numerical scheme is presented for accurate simulation of fluid flow using the lattice Boltzmann equation (LBE) on unstructured mesh. A finite volume approach is adopted to discretize the LBE on a cell-centered, arbitrary shaped, triangular tessellation. The formulation includes a formal, second order discretization using a Total Variation Diminishing (TVD) scheme for the terms representing advection of the distribution function in physical space, due to microscopic particle motion. The advantage of the LBE approach is exploited by implementing the scheme in a new computer code to run on a parallel computing system. Performance of the new formulation is systematically investigated by simulating four benchmark flows of increasing complexity, namely (1) flow in a plane channel, (2) unsteady Couette flow, (3) flow caused by a moving lid over a 2D square cavity and (4) flow over a circular cylinder. For each of these flows, the present scheme is validated with the results from Navier-Stokes computations as well as lattice Boltzmann simulations on regular mesh. It is shown that the scheme is robust and accurate for the different test problems studied.

Analysis of an orthotropic cylindrical shell having a transversely isotropic core subjected to axisymmetric load

A long two-layered circular cylinder having a thin orthotropic outer shell and a thick transversely isotropic core subjected to an axisymmetric radialv line load has been analysed. For analysis of the outer shell the classical thin shell theory was adopted and for analysis of the inner core the elasticity theory was used. The continuity of stresses and deformations at the interface has been satisfied by assumming perfect adhesion between the layers. Numerical results have been presented for two different ratios of outer shell thickness to inner radius and for three different ratios of modulus of elasticity in the radial direction of outer shell to inner core. The results have been compared with the elasticity solution of the same problem to bring out the reliability of this hybrid method. References

On the propagation of a weak shock front: Theory and application

In this paper the kinematics of a weak shock front governed by a hyperbolic system of conservation laws is studied. This is used to develop a method for solving problems, involving the propagation of nonlinear unimodal waves. It consists of first solving the nonlinear wave problem by moving along the bicharacteristics of the system and then fitting the shock into this solution field, so that it satisfies the necessary jump conditions. The kinematics of the shock leads in a natural way to the definition of ldquoshock-raysrdquo, which play the same role as the ldquoraysrdquo in a continuous flow. A special case of a circular cylinder introduced suddenly in a constant streaming flow is studied in detail. The shock fitted in the upstream region propagates with a velocity which is the mean of the velocities of the linear and the nonlinear wave fronts. In the downstream the solution is given by an expansion wave.

Simultaneous heat and mass transfer in unsteady free convection from two-dimensional and axisymmetric bodies

The unsteady laminar free convection boundary layer flows around two-dimensional and axisymmetric bodies placed in an ambient fluid of infinite extent have been studied when the flow is driven by thermal buoyancy forces and buoyancy forces from species diffusion. The unsteadiness in the flow field is caused by both temperature and concentration at the wall which vary arbitrarily with time. The coupled nonlinear partial differential equations with three independent variables governing the flow have been solved numerically using an implicit finite-difference scheme in combination with the quasilinearization technique. Computations have been performed for a circular cylinder and a sphere. The skin friction, heat transfer and mass transfer are strongly dependent on the variation of the wall temperature and concentration with time. Also the skin friction and heat transfer increase or decrease as the buoyancy forces from species diffusion assist and oppose, respectively, the thermal buoyancy force, whereas the mass transfer rate is higher for small values of the ratio of the buoyancy parameters than for large values. The local heat and mass transfer rates are maximum at the stagnation point and they decrease progressively with increase of the angular position from the stagnation point.

Surface-ray tracing on hybrid surfaces of revolution for UTD mutual coupling analysis

An angle invariance property based on Hertz's principle of particle dynamics is employed to facilitate the surface-ray tracing on nondevelopable hybrid quadric surfaces of revolution (h-QUASOR's). This property, when used in conjunction with a Geodesic Constant Method, yields analytical expressions for all the ray-parameters required in the UTD formulation. Differential geometrical considerations require that some of the ray-parameters (defined heuristically in the UTD for the canonical convex surfaces) be modified before the UTD can be applied to such hybrid surfaces. Mutual coupling results for finite-dimensional slots have been presented as an example on a satellite launch vehicle modeled by general paraboloid of revolution and right circular cylinder.

Unsteady separation and vortex shedding from a laminar separation bubble over a bluff body

Boundary layers are subject to favorable and adverse pressure gradients because of both the temporal and spatial components of the pressure gradient. The adverse pressure gradient may cause the flow to separate. In a closed loop unsteady tunnel we have studied the initiation of separation in unsteady flow past a constriction (bluff body) in a channel. We have proposed two important scalings for the time when boundary layer separates. One is based on the local pressure gradient and the other is a convective time scale based on boundary layer parameters. The flow visualization using a dye injection technique shows the flow structure past the body. Nondimensional shedding frequency (Strouhal number) is calculated based on boundary layer and momentum thicknesses. Strouhal number based on the momentum thickness shows a close agreement with that for flat plate and circular cylinder.

Chiral molecular self-assembly of phospholipid tubules: A circular dichroism study

We report on spectroscopic studies of the chiral structure in phospholipid tubules formed in mixtures of alcohol and water. Synthetic phospholipids containing diacetylenic moieties in the acyl chains self-assemble into hollow, cylindrical tubules in appropriate conditions. Circular dichroism provides a direct measure of chirality of the molecular structure. We find that the CD spectra of tubules formed in mixtures of alcohol and water depends strongly on the alcohol used and the lipid concentration. The relative spectral intensity of different circular dichroism bands correlates with the number of bilayers observed using microscopy. The results provide experimental evidence that tubule formation is based on chiral packing of the lipid molecules and that interbilayer interactions are important to the tubule structure

Unsteady Compressible Boundary-Layer Flow On The Stagnation Line Of An Infinite Swept Cylinder

Numerical solutions of flow and heat transfer process on the unsteady flow of a compressible viscous fluid with variable gas properties in the vicinity of the stagnation line of an infinite swept cylinder are presented. Results are given for the case where the unsteady temperature field is produced by (i) a sudden change in the wall temperature (enthalpy) as the impulsive motion is started and (ii) a sudden change in the free-stream velocity. Solutions for the simultaneous development of the thermal and momentum boundary layers are obtained by using quasilinearization technique with an implicit finite difference scheme. Attention is given to the transient phenomenon from the initial flow to the final steady-state distribution. Results are presented for the skin friction and heat transfer coefficients as well as for the velocity and enthalpy profiles. The effects of wail enthalpy parameter, sweep parameter, fluid properties and transpiration cooling on the heat transfer and skin friction are considered.

Bearing capacity of circular footing on geocell-sand mattress overlying clay bed with void

The potential benefits of providing geocell reinforced sand mattress over clay subgrade with void have been investigated through a series of laboratory scale model tests. The parameters varied in the test programme include, thickness of unreinforced sand layer above clay bed, width and height of geocell mattress, relative density of the sand fill in the geocells, and influence of an additional layer of planar geogrid placed at the base of the geocell mattress. The test results indicate that substantial improvement in performance can be obtained with the provision of geocell mattress, of adequate size, over the clay subgrade with void. In order to have beneficial effect, the geocell mattress must spread beyond the void at least a distance equal to the diameter of the void. The influence of the void over the performance of the footing reduces for height of geocell mattress greater than 1.8 times the diameter of the footing. Better improvement in performance is obtained for geocells filled with dense soil. (C) 2008 Elsevier Ltd. All rights reserved.

Broadband Suspended Microstrip Antenna For Circular Polarization

In this paper we propose a circularly polarized (CP) microstrip antenna on a suspended substrate with a coplanar capacitive feed and a slot within the rectangular patch. The antenna has an axial ratio bandwidth (< 3 dB) of 7.1%. The proposed antenna exhibits a much higher impedance bandwidth of about 49% (S11 < -10 dB) and also yields return loss better than -15 dB in the useful range of circular polarization. Measured characteristics of the antenna are in good agreement with the simulated results. The radiation patterns indicate good cross polarization rejection and low back lobe radiations. The design proposed here can be scaled to any frequency of interest.

Simplified dispersion curves for circular cylindrical shells using shallow shell theory.

An alternative derivation of the dispersion relation for the transverse vibration of a circular cylindrical shell is presented. The use of the shallow shell theory model leads to a simpler derivation of the same result. Further, the applicability of the dispersion relation is extended to the axisymmetric mode and the high frequency beam mode.

Forced convection along a longitudinal cylinder embedded in a saturated porous medium

The thermal boundary layer along an isothermal cylinder in a porous 3edium is studied numerically by a finite difference scheme and also using the method of extended perturbation series. The series in terms of the transverse curvature parameter ξ extended to seven terms and is subsequently improved by applying the Shanks transformation twice and thrice, respectively. Results for heat transfer characteristics are found in very good agreement.

Unsteady Attachment Line Flow On A Flat-Plate With Attached Cylinder

The unsteady laminar incompressible boundary-layer attachment-line flow on a flat plate with attached cylinder with heat and mass transfer has been studied when the free stream velocity, mass transfer and surface wall temperature vary arbitrarily with time. The governing partial differential equations with three independent variables have been solved numerically using an implicit finite-difference scheme. The heat transfer was found to be strongly dependent on the Prandtl number, variation of wall temperature with time and dissipation parameter (for large times). However, the free stream velocity distribution and mass transfer affect both the heat transfer and skin friction.

Bending Of Circular Plate On Elastic-Foundation

The governing differential equation of linear, elastic, thin, circular plate of uniform thickness, subjected to uniformly distributed load and resting on Winkler-Pasternak type foundation is solved using ``Chebyshev Polynomials''. Analysis is carried out using Lenczos' technique, both for simply supported and clamped plates. Numerical results thus obtained by perturbing the differential equation for plates without foundation are compared and are found to be in good agreement with the available results. The effect of foundation on central deflection of the plate is shown in the form of graphs.