Double diffusive unsteady free convection on two-dimensional and axisymmetric bodies in a porous medium

The unsteady laminar free convection flow of an incompressible electrically conducting fluid over two-dimensional and axisymmetric bodies embedded in a highly porous medium with an applied magnetic field has been studied. The unsteadiness in the flow field is caused by the variation of the wall temperature and concentration with time. The coupled nonlinear partial differential equations with three independent variables have been solved numerically using an implicit finite-difference scheme in combination with the quasilinearization technique. It is observed that the skin friction, heat transfer and mass transfer increase with the permeability parameter but decrease with the magnetic parameter. The results are strongly dependent on the variation of wall temperature and concentration with time. The skin friction and heat transfer increase or decrease as the buoyancy forces from species diffusion assist or oppose the thermal buoyancy force. However, the mass transfer is found to be higher for small values of the ratio of the buoyancy parameters than for large values

Second-order boundary-layer effects for the unsteady laminar incompressible three-dimensional stagnation-point flow

All the second-order boundary-layer effects on the unsteady laminar incompressible flow at the stagnation-point of a three-dimensional body for both nodal and saddle point regions have been studied. It has been assumed that the free-stream velocity, wall temperature and mass transfer vary arbitrarily with time. The effect of the Prandtl number has been taken into account. The partial differential equations governing the flow have been derived for the first time and then solved numerically unsteady free-stream velocity distributions, the nature of the using an implicit finite-difference scheme. It is found that the stagnation point and the mass transfer strongly affect the skin friction and heat transfer whereas the effects of the Prandtl number and the variation of the wall temperature with time are only on the heat transfer. The skin friction due to the combined effects of first- and second-order boundary layers is less than the skin friction due to, the first-order boundary layers whereas the heat transfer has the opposite behaviour. Suction increases the skin friction and heat transfer but injection does the opposite

Spatial survival probability for one-dimensional fluctuating interfaces in the steady state

We report numerical and analytic results for the spatial survival probability for fluctuating one-dimensional interfaces with Edwards-Wilkinson or Kardar-Parisi-Zhang dynamics in the steady state. Our numerical results are obtained from analysis of steady-state profiles generated by integrating a spatially discretized form of the Edwards-Wilkinson equation to long times. We show that the survival probability exhibits scaling behavior in its dependence on the system size and the "sampling interval" used in the measurement for both "steady-state" and "finite" initial conditions. Analytic results for the scaling functions are obtained from a path-integral treatment of a formulation of the problem in terms of one-dimensional Brownian motion. A "deterministic approximation" is used to obtain closed-form expressions for survival probabilities from the formally exact analytic treatment. The resulting approximate analytic results provide a fairly good description of the numerical data.

A two-dimensional model for crack closure effect in plates under bending

A two-dimensional model is proposed for taking into account the establishment of contact on the compression side of crack faces in plates under bending. An approximate but simple method is developed for evaluating reduction of stress intensity factor due to such ‘crack closure’. Analysis is first carried out permitting interference of the crack faces. Contact forces are then introduced on the crack faces and their magnitudes determined from the consideration that the interference is just eliminated. The method is based partly on finite element analysis and partly on a continuum analysis using Irwin's solution for point loads on the crack line.

A rectangular laminated anisotropic shallow thin shell finite element

This report contains the details of the development of the stiffness matrix for a rectangular laminated anisotropic shallow thin shell finite element. The derivation is done under linear thin shell assumptions. Expressing the assumed displacement state over the middle surface of the shell as products of one-dimensional first-order Hermite interpolation polynomials, it is possible to insure that the displacement state for the assembled set of such elements, to be geometrically admissible. Monotonic convergence of the total potential energy is therefore possible as the modelling is successively refined. The element is systematically evaluated for its performance considering various examples for which analytical or other solutions are available

Unsteady three-dimensional boundary layer flow over a flat plate with attached cylinder

The unsteady laminar incompressible three-dimensional boundary layer flow and heat transfer on a flat plate with an attached cylinder have been studied when the free stream velocity components and wall temperature vary inversely as linear and quadratic functions of time, respectively. The governing semisimilar partial differential equations with three independent variables have been solved numerically using a quasilinear finite-difference scheme. The results indicate that the skin friction increases with parameter ? which characterizes the unsteadiness in the free stream velocity and the streamwise distance Image , but the heat transfer decreases. However, the skin friction and heat transfer are found to change little along Image . The effect of the Prandtl number on the heat transfer is found to be more pronounced when ? is small, whereas the effect of the dissipation parameter is more pronounced when ? is comparatively large.

Unsteady mixed convection over two-dimensional and axisymmetric bodies

The unsteady laminar incompressible mixed convection flow over a two-dimensional body (cylinder) and an axisymmetric body (sphere) has been studied when the buboyancy forces arise from both thermal and mass diffusion and the unsteadiness in the flow field is introduced by the time dependent free stream velocity. The 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. The results indicate that for the thermally assisting flow the local skin friction, heat transfer and mass diffusion are enhanced when the buoyancy force from mass diffusion assists the thermal buoyancy force. But this trend is opposite for the thermally opposing flow. The point of zero skin friction moves upstream due to unsteadiness. No singularity is observed at the point of zero skin friction for unsteady flow unlike steady flow. The flow reversal is observed after a certain instant of time. The velocity overshoot occurs for assisting flows.

Nonadiabatic charge pumping by oscillating potentials in one dimension: Results for infinite system and finite ring

We study charge pumping when a combination of static potentials and potentials oscillating with a time period T is applied in a one-dimensional system of noninteracting electrons. We consider both an infinite system using the Dirac equation in the continuum approximation and a periodic ring with a finite number of sites using the tight-binding model. The infinite system is taken to be coupled to reservoirs on the two sides which are at the same chemical potential and temperature. We consider a model in which oscillating potentials help the electrons to access a transmission resonance produced by the static potentials and show that nonadiabatic pumping violates the simple sin phi rule which is obeyed by adiabatic two-site pumping. For the ring, we do not introduce any reservoirs, and we present a method for calculating the current averaged over an infinite time using the time evolution operator U(T) assuming a purely Hamiltonian evolution. We analytically show that the averaged current is zero if the Hamiltonian is real and time-reversal invariant. Numerical studies indicate another interesting result, namely, that the integrated current is zero for any time dependence of the potential if it is applied to only one site. Finally we study the effects of pumping at two sites on a ring at resonant and nonresonant frequencies, and show that the pumped current has different dependences on the pumping amplitude in the two cases.

An accurate numerical integration scheme for finite rotations using rotation vector parametrization

A numerical integration procedure for rotational motion using a rotation vector parametrization is explored from an engineering perspective by using rudimentary vector analysis. The incremental rotation vector, angular velocity and acceleration correspond to different tangent spaces of the rotation manifold at different times and have a non-vectorial character. We rewrite the equation of motion in terms of vectors lying in the same tangent space, facilitating vector space operations consistent with the underlying geometric structure. While any integration algorithm (that works within a vector space setting) may be used, we presently employ a family of explicit Runge-Kutta algorithms to solve this equation. While this work is primarily motivated out of a need for highly accurate numerical solutions of dissipative rotational systems of engineering interest, we also compare the numerical performance of the present scheme with some of the invariant preserving schemes, namely ALGO-C1, STW, LIEMIDEA] and SUBCYC-M. Numerical results show better local accuracy via the present approach vis-a-vis the preserving algorithms. It is also noted that the preserving algorithms do not simultaneously preserve all constants of motion. We incorporate adaptive time-stepping within the present scheme and this in turn enables still higher accuracy and a `near preservation' of constants of motion over significantly longer intervals. (C) 2010 The Franklin Institute. Published by Elsevier Ltd. All rights reserved.

Unsteady laminar compressible boundary-layer flow at a three-dimensional stagnation point

Unsteady laminar compressible boundary-layer flow with variable properties at a three-dimensional stagnation point for both cold and hot walls has been studied for the case when the velocity of the incident stream varies arbitrarily with time. The partial differential equations governing the flow have been solved numerically using an implicit finite-difference scheme. Computations have been carried out for two particular unsteady free-stream velocity distributions: (i) an accelerating stream and (ii) a fluctuating stream. The results indicate that the variation of the density-viscosity product across the boundary layer, the wall temperature and the nature of stagnation point significantly affect the skin friction and heat transfer.

Analysis of edge delaminations in laminates through combined use of quasi-three-dimensional, eight-noded, two-noded and transition elements

The use of appropriate finite elements in different regions of a stressed solid can be expected to be economical in computing its stress response. This concept is exploited here in studying stresses near free edges in laminated coupons. The well known free edge problem of [0/90], symmetric laminate is considered to illustrate the application of the concept. The laminate is modelled as a combination of three distinct regions. Quasi-three-dimensional eight-noded quadrilateral isoparametric elements (Q3D8) are used at and near the free edge of the laminate and two-noded line elements (Q3D2) are used in the region away from the free edge. A transition element (Q3DT) provides a smooth inter-phase zone between the two regions. Significant reduction in the problem size and hence in the computational time and cost have been achieved at almost no loss of accuracy.

On the relation between rotation increments in different tangent spaces

In computational mechanics, finite rotations are often represented by rotation vectors. Rotation vector increments corresponding to different tangent: spaces are generally related by a linear operator, known as the tangential transformation T. In this note, we derive the higher order terms that are usually left out in linear relation. The exact nonlinear relation is also presented. Errors via the linearized T are numerically estimated. While the concept of T arises out of the nonlinear characteristics of the rotation manifold, it has been derived via tensor analysis in the context of computational mechanics (Cardona and Geradin, 1988). We investigate the operator T from a Lie group perspective, which provides a better insight and a 1-1 correspondence between approaches based on tensor analysis and the standard matrix Lie group theory. (C) 2010 Elsevier Ltd. All rights reserved.

2-Dimensional Linear-Stability of Premixed Laminar Flames Under Zero-Gravity

This paper reports on the numerical study of the linear stability of laminar premixed flames under zero gravity. The study specifically addresses the dependence of stability on finite rate chemistry with low activation energy and variable thermodynamic and transport properties. The calculations show that activation energy and details of chemistry play a minor role in altering the linear neutral stability results from asymptotic analysis. Variable specific heat makes a marginal change to the stability. Variable transport properties on the other hand tend to substantially enhance the stability from critical wave number of about 0.5 to 0.20. Also, it appears that the effects of variable properties tend to nullify the effects of non-unity Lewis number. When the Lewis number of a single species is different from unity, as will happen in a hydrogen-air premixed flame, the stability results remain close to that of unity Lewis number.

Analysis of the evanescent field coupling between a fiber and a planar waveguide using the explicit finite difference method with a nonuniform mesh (Journal Paper)

We investigate an optical waveguide system consisting of an unclad fiber core suspended at a constant distance parallel to the surface of a planar waveguide. The coupling and propagation of light in the combined system is studied using the three-dimensional explicit finite difference beam propagation method with a nonuniform mesh configuration. The power loss in the fiber and the field distribution in the waveguide are studied as a function of various parameters, such as index changes, index profile, and propagation distance, for the combined system.

Flap‐Lag Damping in Hover and Forward Flight with a Three‐Dimensional Wake

Prediction of lag damping is difficult owing to the delicate balance of drag, induced drag and Coriolis forces in the in‐plane direction. Moreover, induced drag” is sensitive to dynamic wake, bath shed and trailing components, and thus its prediction requires adequate unsteady‐wake representation. Accordingly, rigid‐blade flap‐lag equations are coupled with a three‐dimensional finite‐state wake model; three isolatcd rotor canfigurations with three, four and five blades are treated over a range of thrust levels, tack numbers, lag frequencies and advance ratios. The investigation includes convergence characteristics of damping with respect to the number of radial shape functions and harmonics of the wake model for multiblade modes of low frequency (< 1/ rev.) to high frequency (> 1/rev.). Predicted flap and lag damping levels are then compared with similar predictions with 1) rigid wake (no unsteady induced now), 2) Loewy lift deficiency and 3) dynamic inflow. The coverage also includes correlations with the measured lag regressive‐mode damping in hover and forward flight and comparisons with similar correlations with dynamic inflow. Lag‐damping predictions with the dynamic wake model are consistently higher than the predictions with the dynamic inflow model; even for the low frequency lag regressive mode, the number of wake harmonics should at least be equal to twice the number of blades.