Theoretical study of a 16-µm CO2 downstream-mixing gasdynamic laser: A two-dimensional approach

A 16-µm CO2-N2 downstream-mixing gasdynamic laser, where a cold CO2 stream is mixed with a vibrationally excited N2 stream at the exit of the nozzle, is studied theoretically. The flow field is analyzed using a two-dimensional, unsteady, laminar and viscous flow model including appropriate finite-rate vibrational kinetic equations. The analysis showed that local small-signal gain up to 21.75 m−1 can be obtained for a N2 reservoir temperature of 2000 K and a velocity ratio of 1:1 between the CO2 and N2 mixing streams. Applied Physics Letters is copyrighted by The American Institute of Physics.

A Method of Characteristics for Solving Two-dimensional Unsteady Flow Problems

The aim of this investigation is to evolve a method of solving two-dimensional unsteady flow problems by the method of characteristics. This involves the reduction of the given system of equations to an equivalent system where only interior derivatives occur on a characteristic surface. From this system, four special bicharacteristic directional derivatives are chosen. A finite difference scheme is prescribed for solving the equations. General rectangular lattices are also considered. As an example, we investigate the propagation of an initial pressure distribution in a medium at rest.

A parametric differentiation version with finite-difference scheme applicable to a class of problems in boundary layer flow with massive blowing

A numerical procedure, based on the parametric differentiation and implicit finite difference scheme, has been developed for a class of problems in the boundary-layer theory for saddle-point regions. Here, the results are presented for the case of a three-dimensional stagnation-point flow with massive blowing. The method compares very well with other methods for particular cases (zero or small mass blowing). Results emphasize that the present numerical procedure is well suited for the solution of saddle-point flows with massive blowing, which could not be solved by other methods.

Finite element analysis of curved anisotropic laminated hollow bars

Curved hollow bars of laminated anisotropic construction are used as structural members in many industries. They are used in order to save weight without loss of stiffness in comparison with solid sections. In this paper are presented the details of the development of the stiffness matrices of laminated anisotropic curved hollow bars under line member assumptions for two typical sections, circular and square. They are 16dof elements which make use of one-dimensional first-order Hermite interpolation polynomials for the description of assumed displacement state. Problems for which analytical or other solutions are available are first solved using these elements. Good agreement was found between the results. In order to show the capability of the element, application is made to carbon fibre reinforced plastic layered anisotropic curved hollow bars.

A method for separation of stresses in Two- and three-dimensional photoelasticity

A method for separation of stresses in two and three-dimensional photo elasticity using the harmonisation ofjrst stress invariant along a straight section is deve- ,dped. For two-dimensions, the equations of equilibrium are reformulated in terms ojsum and difference of normal stresses and relations are obtained which can be used for harmonisation of the first invariant of stress along a straight section. A similar procedure is adopted for three-dimensions by making use of the Beltrmi-MicheN equations. The new relations are used in finite d~yerencefo rm to evaluate the sum of normal stresses along straight sections in a three-dimensional body. The method requires photoelastic data along the section as well ~rra djacent sections. This method could be used as an alternative to the shear d@erence method for separation of stresses in photoelasticity. 7he accuracy and reliability of the method is ver$ed by applying the method to problems whose solutions are known.

Matrices of finite Lorentz transformations in a noncompact basis. III. Completeness relation for O(2, 1)

A parametrization of the elements of the three-dimensional Lorentz group O(2, 1), suited to the use of a noncompact O(1, 1) basis in its unitary representations, is derived and used to set up the representation matrices for the entire group. The Plancherel formula for O(2, 1) is then expressed in this basis.

Unsteady Incompressible Three-Dimensional Asymmetric Stagnation-Point Boundary Layers

The unsteady laminar incompressible boundary-layer flow near the three-dimensional asymmetric stagnation point has been studied under the assumptions that the free-stream velocity, wall temperature, and surface mass transfer vary arbitrarily with time. The partial differential equations governing the flow have been solved numerically using an implicit finite-difference scheme. It is found that in contrast with the symmetric flow, the maximum heat transfer occurs away from the stagnation point due to the decrease in the boundary-layer thickness. The effect of the variation of the wall temperature with time on heat transfer is strong. The skin friction and heat transfer due to asymmetric flow only are comparatively less affected by the mass transfer as compared to those of symmetric flow.

Supersolid and solitonic phases in the one-dimensional extended Bose-Hubbard model

We report our findings on the quantum phase transitions in cold bosonic atoms in a one-dimensional optical lattice using the finite-size density-matrix renormalization-group method in the framework of the extended Bose-Hubbard model. We consider wide ranges of values for the filling factors and the nearest-neighbor interactions. At commensurate fillings, we obtain two different types of charge-density wave phases and a Mott insulator phase. However, departure from commensurate fillings yields the exotic supersolid phase where both the crystalline and the superfluid orders coexist. In addition, we obtain the signatures for the solitary waves and the superfluid phase.

Consistent quaternion interpolation for objective finite element approximation of geometrically exact beam

We explore an isoparametric interpolation of total quaternion for geometrically consistent, strain-objective and path-independent finite element solutions of the geometrically exact beam. This interpolation is a variant of the broader class known as slerp. The equivalence between the proposed interpolation and that of relative rotation is shown without any recourse to local bijection between quaternions and rotations. We show that, for a two-noded beam element, the use of relative rotation is not mandatory for attaining consistency cum objectivity and an appropriate interpolation of total rotation variables is sufficient. The interpolation of total quaternion, which is computationally more efficient than the one based on local rotations, converts nodal rotation vectors to quaternions and interpolates them in a manner consistent with the character of the rotation manifold. This interpolation, unlike the additive interpolation of total rotation, corresponds to a geodesic on the rotation manifold. For beam elements with more than two nodes, however, a consistent extension of the proposed quaternion interpolation is difficult. Alternatively, a quaternion-based procedure involving interpolation of relative rotations is proposed for such higher order elements. We also briefly discuss a strategy for the removal of possible singularity in the interpolation of quaternions, proposed in [I. Romero, The interpolation of rotations and its application to finite element models of geometrically exact rods, Comput. Mech. 34 (2004) 121–133]. The strain-objectivity and path-independence of solutions are justified theoretically and then demonstrated through numerical experiments. This study, being focused only on the interpolation of rotations, uses a standard finite element discretization, as adopted by Simo and Vu-Quoc [J.C. Simo, L. Vu-Quoc, A three-dimensional finite rod model part II: computational aspects, Comput. Methods Appl. Mech. Engrg. 58 (1986) 79–116]. The rotation update is achieved via quaternion multiplication followed by the extraction of the rotation vector. Nodal rotations are stored in terms of rotation vectors and no secondary storages are required.

The partition of unity finite element method for elastic wave propagation in Reissner-Mindlin plates

This paper reports a numerical method for modelling the elastic wave propagation in plates. The method is based on the partition of unity approach, in which the approximate spectral properties of the infinite dimensional system are embedded within the space of a conventional finite element method through a consistent technique of waveform enrichment. The technique is general, such that it can be applied to the Lagrangian family of finite elements with specific waveform enrichment schemes, depending on the dominant modes of wave propagation in the physical system. A four-noded element for the Reissner-indlin plate is derived in this paper, which is free of shear locking. Such a locking-free property is achieved by removing the transverse displacement degrees of freedom from the element nodal variables and by recovering the same through a line integral and a weak constraint in the frequency domain. As a result, the frequency-dependent stiffness matrix and the mass matrix are obtained, which capture the higher frequency response with even coarse meshes, accurately. The steps involved in the numerical implementation of such element are discussed in details. Numerical studies on the performance of the proposed element are reported by considering a number of cases, which show very good accuracy and low computational cost. Copyright (C)006 John Wiley & Sons, Ltd.

Thermal stresses in a finite hollow cylinder due to an axisymmetric temperature field at the end surface

A three-dimensional exact solution for determining the thermal stresses in a finite hollow cylinder subject to a steady state axisymmetric temperature field over one of its end surfaces has been given. Numerical results for a hollow cylinder, having lenght to outer diameter ratio equal to one and inner to outer diameter ratio equal to 0.75, subjected to a symmetric temperature variation over the end surfaces of the cylinder have been given.

Thermal stresses in a finite solid cylinder due to an axisymmetric temperature field at the end surface.

A three-dimensional rigorous solution for determining thermal stresses in a finite solid cylinder due to a steady state axisymmetric temperature field over one of its end surfaces is given. Numerical results for a solid cylinder having a length to diameter ratio equal to one and subjected to a symmetric temperature variation over half the radius of the cylinder at the end surfaces are included. These results have been compared with the results of the approximate solution given by W. Nowacki.

Comparison of elasticity and shell-theory solutions for finite circular cylindrical shells

A finite circular cylindrical shell subjected to a band of uniform pressure on its outer rim was investigated, using three-dimensional elasticity theory and the classical shell theories of Timoshenko (or Donnell) and Flügge. Detailed comparison of the resulting stresses and displacements was carried out for shells with ratios of inner to outer shell radii equal to 0.80, 0.85, 0.90 and 0.93 and for ratios of outer shell diameter to length of the shell equal to 0.5, 1 and 2. The ratio of band width to length of the shell was 0.2 and Poisson's ratio used was equal to 0.3. An Elliot 803 digital computer was used for numerical computations.

Three-Dimensional simulation of deformation fields underneath vickers indenter:Effects of power-law Plasticity

The effects of power-law plasticity (yield strength and strain hardening exponent) on the plastic strain distribution underneath a Vickers indenter was systematically investigated by recourse to three-dimensional finite element analysis, motivated by the experimental macro-and micro-indentation on heat-treated Al-Zn-Mg alloy. For meaningful comparison between simulated and experimental results, the experimental heat treatment was carefully designed such that Al alloy achieve similar yield strength with different strain hardening exponent, and vice versa. On the other hand, full 3D simulation of Vickers indentation was conducted to capture subsurface strain distribution. Subtle differences and similarities were discussed based on the strain field shape, size and magnitude for the isolated effect of yield strength and strain hardening exponent.

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.