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.

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.

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.

Thermomechanical generalization of pin joints in finite plates

The understanding of thermoelastic behaviour of joints is significant in order to ensure the integrity of large and complex structures exposed to a thermal environment, particularly in fields such as aerospace and nuclear engineering. Thermomechanical generalization of partial contact behaviour of a pin joint under combined in-plane mechanical loading and on-axis unidirectional heat flow has already been established by the authors for the analytically simpler domains of large plates. This paper successfully extends the on-going investigation to a single pin in a finite rectangular isotropic plate as a two-dimensional abstraction from a practical situation of a multipin fastener joint. The finite element method is used to analyse the joint problem under on-axis thermomechanical loading and unified load-contact relationships are established for a class of loading conditions.

Time-domain-finite-wave analysis of the engine exhaust system by means of the stationary-frame method of characteristics. Part I. Theory

Time-domain-finite-wave analysis of the engine exhaust system is usually done using the method of characteristics. This makes use of either the moving frame method, or the stationary frame method. The stationary frame method is more convenient than its counterpart inasmuch as it avoids the tedium of graphical computations. In this paper (part I), the stationary-frame computational scheme along with the boundary conditions has been implemented. The analysis of a uniform tube, cavity-pipe junction including the engine and the radiation ends, and also the simple area discontinuities has been presented. The analysis has been done accounting for wall friction and heat-transfer for a one-dimensional unsteady flow. In the process, a few inconsistencies in the formulations reported in the literature have been pointed out and corrected. In the accompanying paper (part II) results obtained from the simulation are shown to be in good agreement with the experimental observations.

Unsteady three-dimensional boundary layer flow due to a stretching surface

In this numerical study, the unsteady laminar incompressible boundary-layer flow over a continuously stretching surface has been investigated when the velocity of the stretching surface varies arbitrarily with time. Both the nodal and the saddle point regions of flow have been considered for the analysis. Also, constant wall temperature/concentration and constant heat/mass flux at the stretching surface have been taken into account. The quasilinearisation method with an implicit finite-difference scheme is used in the nodal point region (0 less-than-or-equal-to c less-than-or-equal-to 1) where c denotes the stretching ratio. This method fails in the saddle point region (-1 less-than-or-equal-to c less-than-or-equal-to 0) due to the occurrence of reverse flow in the y-component of velocity. In order to overcome this difficulty, the method of parametric differentiation with an implicit finite-difference scheme is used, where the values at c = 0 are taken as starting values. Results have been obtained for the stretching velocities which are accelerating and decelerating with time. Results show that the skin friction, the heat transfer and the mass transfer parameters respond significantly to the time dependent stretching velocities. Suction (A > 0) is found to be an important parameter in obtaining convergent solution in the case of the saddle point region of flow. The Prandtl number and the Schmidt number strongly affect the heat and mass transfer of the diffusing species, respectively.

Solid finite elements through three decades

conventionally, solid finite elements have been looked upon as just generalizations of two-dimensional finite elements. In this article we trace their development starting from the days of their inception. Keeping in tune with our perceptions on developing finite elements, without taking recourse to any extra variational techniques, we discuss a few of the techniques which have been applied to solid finite elements. Finally we critically examine our own work on formulating solid finite elements based on the solutions to the Navier equations.

Unsteady nonsimilar two-dimensional and axisymmetric water boundary layers with variable viscosity and prandtl number

The influence of temperature-dependent viscosity and Prandtl number on the unsteady laminar nonsimilar forced convection flow over two-dimensional and axisymmetric bodies has been examined where the unsteadiness and (or) nonsimilarity are (is) due to the free stream velocity, mass transfer, and transverse curvature. The partial differential equations governing the flow which involve three independent variables have been solved numerically using an implicit finite-difference scheme along with a quasilinearization technique. It is found that both the skin friction and heat transfer strongly respond to the unsteady free stream velocity distributions. The unsteadiness and injection cause the location of zero skin friction to move upstream. However, the effect of variable viscosity and Prandtl number is to move it downstream. The heat transfer is found to depend strongly on viscous dissipation, but the skin friction is little affected by it. In general, the results pertaining to variable fluid properties differ significantly, from those of constant fluid properties.

Flap-lag-torsion stability in hover and forward flight with three-dimensional wake

The contributions of full-wake dynamics in trim analysis are demonstrated for finding the control inputs and periodic responses simultaneously, as well as in Floquet eigenanalysis for finding the damping levels. The equations of flap bending, lag bending, and torsion are coupled with a three-dimensional, finite state wake, and low-frequency (<1/rev) to high frequency (>1/rev) multiblade modes are considered. Full blade-wake dynamics is used in trim analysis and Floquet eigenanalysis. A uniform cantilever blade in trimmed flight is investigated over a range of thrust levels, advance ratios, number of blades, and blade torsional frequencies. The investigation includes the convergence characteristics of control inputs, periodic responses, and damping levels with respect to the number of spatial azimuthal harmonics and radial shape functions in the wake representation. It also includes correlation with the measured lag damping of a three-bladed untrimmed rotor. The parametric study shows the dominant influence of wake dynamics on control inputs, periodic responses, and damping levels, and wake theory generally improves the correlation.