Boundary-layer effects on internal flow choking in dual-thrust solid rocket motors

Theoretical studies have been carried out to examine internal flow choking in the inert simulators of a dual-thrust motor. Using a two-dimensional k-omega turbulence model, detailed parametric studies have been carried out to examine aerodynamic choking and the existence of a fluid throat at the transition region during the startup transient of dual-thrust motors. This code solves standard k-omega turbulence equations with shear flow corrections using a coupled second-order-implicit unsteady formulation. In the numerical study, a fully implicit finite volume scheme of the compressible, Reynolds-averaged, Navier-Stokes equations is employed. It was observed that, at the subsonic inflow conditions, there is a possibility of the occurrence of internal flow choking in dual-thrust motors due to the formation of a fluid throat at the beginning of the transition region induced by area blockage caused by boundary-layer-displacement thickness. It has been observed that a 55% increase in the upstream port area of the dual-thrust motor contributes to a 25% reduction in blockage factor at the transition region, which could negate the internal How choking and supplement with an early choking of the dual-thrust motor nozzle. If the height of the upstream port relative to the motor length is too small, the developing boundary layers from either side of the port can interact, leading to a choked,flow. On the other hand, if the developing boundary layers are far enough apart, then choking does not occur. The blockage factor is greater in magnitude for the choked case than for the unchoked case. More tangible explanations are presented in this paper for the boundary-layer blockage and the internal flow choking in dual-thrust motors, which hitherto has been unexplored.

Fracture analysis of stiffened panels under combined tensile, bending, and shear loads

The objective is to present the formulation of numerically integrated modified virtual crack closure integral technique for concentrically and eccentrically stiffened panels for computation of strain-energy release rate and stress intensity factor based on linear elastic fracture mechanics principles. Fracture analysis of cracked stiffened panels under combined tensile, bending, and shear loads has been conducted by employing the stiffened plate/shell finite element model, MQL9S2. This model can be used to analyze plates with arbitrarily located concentric/eccentric stiffeners, without increasing the total number of degrees of freedom, of the plate element. Parametric studies on fracture analysis of stiffened plates under combined tensile and moment loads have been conducted. Based on the results of parametric,studies, polynomial curve fitting has been carried out to get best-fit equations corresponding to each of the stiffener positions. These equations can be used for computation of stress intensity factor for cracked stiffened plates subjected to tensile and moment loads for a given plate size, stiffener configuration, and stiffener position without conducting finite element analysis.

Adhesively-bonded Patch Repair with Composites

Adhesively-bonded composite patch repairs over cracked or corrosion-damaged metallic aircraft structures have shown great promise for extending life of ageing structures. This study presents the numerical investigation into the interface behaviour of adhesively-bonded cracked aluminum alloy substrate patched with fibre-reinforced composite material. The adhesive is modelled as an elasto-plastic bilinear material to characterise the debond behaviour, while the defective substrate is regarded as linear elastic continuum. Two typical patch shapes were selected based on information available in the literature. Geometric and material nonlinear analyses for square and octagonal patches were performed to capture peel and shear stresses developed between the substrate and the patch to examine the possibility of interface delamination/debonding. Parametric studies on adhesive thickness and patch thickness were carried out to predict their infuence on damage tolerance of repaired structures.

Behaviour of Bi-adhesive Joints

The use of relatively low modulus adhesive at the ends of overlap in a bi-adhesive bondline of a bonded joint can reduce the stress concentration significantly and, therefore, potentially lead to higher strength of the joint. This study presents the two-dimensional and three-dimensional nonlinear (geometric and material) finite element analyses of adhesively bonded single lap joints having modulus-graded bondline under monotonic loading conditions. The adhesives were modelled as an elasto-plastic multi-linear material, while the substrates were regarded as both linear elastic and bi-linear elasto-plastic material. The computational simulations have been performed to investigate the bondline behaviour by studying the stress and strain distributions both at the mid-plane as well as at the interface of the bondline. It has been observed that the static strength is higher for joints with bi-adhesive bondlines compared to those with single adhesives in bondline. Higher joint strength has also been observed for optimum bi-adhesive bondline ratio through parametric studies. Effects of load level, and bondline thickness on stress distribution in the bi-adhesive bondline have also been studied. 3D analysis results reveal the existence of complex multi-axial stress/strain state at the ends of the overlap in the bondline which cannot be observed in 2D plane strain analysis. About 1/3rd of the width of the joint from the free edge in the width direction has 3D stress state, especially in the compliant adhesive of the bondline. Magnitudes of longitudinal and lateral stress/strain components are comparable to peel stress/strain components. It has also been analytically shown that the in-plane global stiffness of the joint remains unaffected by modulus gradation of the bondline adhesive. (C) Koninklijke Brill NV, Leiden, 2010.

Robust prediction of residual velocities and ballistic limits of projectiles for impact on thin aluminium plates

This paper deals with the simulation-driven study of the impact of hardened steel projectiles on thin aluminium target plates using explicit finite element analysis as implemented in LS-DYNA. The evaluation of finite element modelling includes a comprehensive mesh convergence study using shell elements for representing target plates and the solid element-based representation of ogivalnosed projectiles. A user-friendly automatic contact detection algorithm is used for capturing interaction between the projectile and the target plate. It is shown that the proper choice of mesh density and strain rate-dependent material properties is crucial as these parameters significantly affect the computed residual velocity. The efficacy of correlation with experimental data is adjudged in terms of a 'correlation index' defined in the present study for which values close to unity are desirable.By simulating laboratory impact tests on thin aluminium plates carried out by earlier investigators, extremely good prediction of experimental ballistic limits has been observed with correlation indices approaching unity. Additional simulation-based parametric studies have been carried out and results consistent with test data have been obtained. The simulation procedures followed in the present study can be applied with confidence in designing thin aluminium armour plates for protection against low calibre projectiles.

Adhesively-bonded Patch Repair with Composites

Adhesively-bonded composite patch repairs over cracked or corrosion-damaged metallic aircraft structures have shown great promise for extending life of ageing structures. This study presents the numerical investigation into the interface behaviour of adhesively-bonded cracked aluminum alloy substrate patched with fibre-reinforced composite material. The adhesive is modelled as an elasto-plastic bilinear material to characterise the debond behaviour, while the defective substrate is regarded as linear elastic continuum. Two typical patch shapes were selected based on information available in the literature. Geometric and material nonlinear analyses for square and octagonal patches were performed to capture peel and shear stresses developed between the substrate and the patch to examine the possibility of interface delamination/debonding. Parametric studies on adhesive thickness and patch thickness were carried out to predict their infuence on damage tolerance of repaired structures.

A dragonfly inspired flapping wing actuated by electro active polymers

An energy-based variational approach is used for structural dynamic modeling of the IPMC (Ionic Polymer Metal Composites) flapping wing. Dynamic characteristics of the wing are analyzed using numerical simulations. Starting with the initial design, critical parameters which have influence on the performance of the wing are identified through parametric studies. An optimization study is performed to obtain improved flapping actuation of the IPMC wing. It is shown that the optimization algorithm leads to a flapping wing with dimensions similar to the dragonfly Aeshna Multicolor wing. An unsteady aerodynamic model based on modified strip theory is used to obtain the aerodynamic forces. It is found that the IPMC wing generates sufficient lift to support its own weight and carry a small payload. It is therefore a potential candidate for flapping wing of micro air vehicles.

Flow-Acoustic Performance of Perforated Element Mufflers With Application to Design

A fairly comprehensive computer program incorporating explicit expressions for the four-pole parameters of concentric-tube resonators, plug mufflers, and three-duct cross-flow perforated elements has been used for parametric studies. The parameters considered are hole diameter, the center-to-center distance between consecutive holes (which decides porosity), the incoming mean flow Mach number, the area expansion ratio, the number of partitions of chambers within a given overall shell length, and the relative lengths of these partitions or chambers, all normalized with respect to the exhaust pipe diameter. Transmission loss has been plotted as a function of a normalized frequency parameter. Additionally, the effect of the tail pipe length on insertion loss for an anechoic source has also been studied. These studies have been supplemented by empirical expressions for the normalized static pressure drop for different types of perforated-element mufflers developed from experimental observations.

Elliptic slot antenna for broadband wireless communications

This paper presents the design of a broadband antenna suitable for wireless communications operating over the frequency range of 3.1-10.6 GHz. Parametric studies on the effect of stub and elliptic slot have been carried out to arrive at optimum dimensions to achieve enhanced bandwidth of the proposed antenna. An experimental antenna has been designed and tested to validate the proposed design. Measured return loss characteristics have been compared against the simulation results. Simulated radiation patterns at 3.1 GHz, 6.85 GHz and 10.6 GHz have also been presented in this paper.

Oscillatory flow and temperature fields in an open tube with temperature difference across the ends

The oscillating flow and temperature field in an open tube subjected to cryogenic temperature at the cold end and ambient temperature at the hot end is studied numerically. The flow is driven by a time-wise sinusoidally varying pressure at the cold end. The conjugate problem takes into account the interaction of oscillatory flow with the heat conduction in the tube wall. The full set of compressible flow equations with axisymmetry assumption are solved with a pressure correction algorithm. Parametric studies are conducted with frequencies of 5-15 Hz, with one end maintained at 100 K and other end at 300 K. The flow and temperature distributions and the cooldown characteristics are obtained. The frequency and pressure amplitude have negligible effect on the time averaged Nusselt number. Pressure amplitude is an important factor determining the enthalpy flow through the solid wall. The frequency of operation has considerable effect on penetration of temperature into the tube. The density variation has strong influence on property profiles during cooldown. The present study is expected to be of interest in applications such as pulse tube refrigerators and other cryocoolers, where oscillatory flows occur in open tubes. (C) 2011 Elsevier Ltd. All rights reserved.

Effective thermal conductivity of a heat generating rod bundle dissipating heat by natural convection and radiation

A numerical study of conjugate natural convection and surface radiation in a horizontal hexagonal sheath housing 19 solid heat generating rods with cladding and argon as the fill gas, is performed. The natural convection in the sheath is driven by the volumetric heat generation in the solid rods. The problem is solved using the FLUENT CFD code. A correlation is obtained to predict the maximum temperature in the rod bundle for different pitch-to-diameter ratios and heat generating rates. The effective thermal conductivity is related to the heat generation rate, maximum temperature and the sheath temperature. Results are presented for the dimensionless maximum temperature, Rayleigh number and the contribution of radiation with changing emissivity, total wattage and the pitch-to-diameter ratio. In the simulation of a larger system that contains a rod bundle, the effective thermal conductivity facilitates simplified modelling of the rod bundle by treating it as a solid of effective thermal conductivity. The parametric studies revealed that the contribution of radiation can be 38-65% of the total heat generation, for the parameter ranges chosen. Data for critical Rayleigh number above which natural convection comes into effect is also presented. (C) 2011 Elsevier B.V. All rights reserved.

Effect of natural convection on oscillating flow in a pipe with cryogenic temperature difference across the ends

The effect of natural convection on the oscillatory flow in an open-ended pipe driven by a timewise sinusoidally varying pressure at one end and subjected to an ambient-to-cryogenic temperature difference across the ends, is numerically studied. Conjugate effects arising out of the interaction of oscillatory flow with heat conduction in the pipe wall are taken into account by considering a finite thickness wall with an insulated exterior surface. Two cases, namely, one with natural convection acting downwards and the other, with natural convection acting upwards, are considered. The full set of compressible flow equations with axissymmetry are solved using a pressure correction algorithm. Parametric studies are conducted with frequencies in the range 5-15 Hz for an end-to-end temperature difference of 200 and 50 K. Results are obtained for the variation of velocity, temperature. Nusselt number and the phase relationship between mass flow rate and temperature. It is found that the Rayleigh number has a minimal effect on the time averaged Nusselt number and phase angle. However, it does influence the local variation of velocity and Nusselt number over one cycle. The natural convection and pressure amplitude have influence on the energy flow through the gas and solid. (C) 2011 Elsevier Ltd. All rights reserved.

Acoustical behavior of single inlet and multiple outlet elliptical cylindrical chamber muffler

Transmission loss (TL) of an elliptical cylindrical chamber muffler having a single side/end inlet and multiple side/end outlet is analyzed by means of the 3-D semi-analytical formulation based upon the modal expansion (in terms of the angular and radial Mathieu functions) and the Green's function. The acoustic pressure response obtained in terms of Green's function is integrated over surface area of the side/end ports (modeled as rigid pistons) and upon subsequent division by the port area, yields the acoustic pressure response or impedance Z] matrix parameters due to the uniform piston-driven model. The 3-D semi-analytical results are found to be in excellent agreement with the results obtained by means of 3-D FEA (SYSNOISE) simulations, thereby validating the semi-analytical procedure suggested in this work. Parametric studies such as the effect of chamber length (L), angular and axial locations of the ports, interchanging the locations of inlet and outlet ports as well as the addition of an outlet port for double outlet mufflers on the TL performance are reported, thereby leading to the formulation of design guidelines for obtaining muffler configurations exhibiting a broad-band TL spectrum. One such configuration is an axially long chamber having side-inlet and side-outlet ports such that one of the side ports is located at half the axial length on themajor/minor axis and the other side port is located at three-quarters (or one-quarter) of the axial length on the minor/major axis. (C) 2012 Institute of Noise Control Engineering.

Wave propagation in stiffened structures using spectrally formulated finite element

In this work, the wave propagation analysis of built-up composite structures is performed using frequency domain spectral finite elements, to study the high frequency wave responses. The paper discusses basically two methods for modeling stiffened structures. In the first method, the concept of assembly of 2D spectral plate elements is used to model a built-up structure. In the second approach, spectral finite element method (SFEM) model is developed to model skin-stiffener structures, where the skin is considered as plate element and the stiffener as beam element. The SFEM model developed using the plate-beam coupling approach is then used to model wave propagation in a multiple stiffened structure and also extended to model the stiffened structures with different cross sections such as T-section, I-section and hat section. A number of parametric studies are performed to capture the mode coupling, that is, the flexural-axial coupling present in the wave responses.

Modeling of Piezocomposites using Variational Asymptotic Method

This work presents the development of piezocomposites made up of Macro Fiber Composites (MFCs) for aerospace applications and specifically involves, their computational analysis, material characterization and certain parametric studies. MFC was developed by NASA Langley Research Center in 1996 and currently is being distributed by Smart Material Co. 1] worldwide and finds applications both as an actuator as well as for sensor in various engineering applications. In this work, MFC is being modeled as an actuator and a theoretical formulation based on Variational Asymptotic Method (VAM) 2] is presented to analyse the laminates made up of MFCs. VAM minimizes the total electro-mechanical energy for the MFC laminate and approaches the exact solution asymptotically by making use of certain small parameters inherent to the problem through dimensional reduction. VAM provides closed form solutions for 1D constitutive law, recovery relations of warpings, 3D stress/strain fields and displacements and hence an ideal tool for carrying out parametric and design studies in such applications. VAM is geometrically exact and offers rigorous material characterization through cross-sectional analysis and dimensional reduction.