Integrating strong and weak discontinuities without integration subcells and example applications in an XFEM/GFEM framework

Partition of unity methods, such as the extended finite element method, allows discontinuities to be simulated independently of the mesh (Int. J. Numer. Meth. Engng. 1999; 45:601-620). This eliminates the need for the mesh to be aligned with the discontinuity or cumbersome re-meshing, as the discontinuity evolves. However, to compute the stiffness matrix of the elements intersected by the discontinuity, a subdivision of the elements into quadrature subcells aligned with the discontinuity is commonly adopted. In this paper, we use a simple integration technique, proposed for polygonal domains (Int. J. Nuttier Meth. Engng 2009; 80(1):103-134. DOI: 10.1002/nme.2589) to suppress the need for element subdivision. Numerical results presented for a few benchmark problems in the context of linear elastic fracture mechanics and a multi-material problem show that the proposed method yields accurate results. Owing to its simplicity, the proposed integration technique can be easily integrated in any existing code. Copyright (C) 2010 John Wiley & Sons, Ltd.

On-board correction of systematic error of Earth sensors

Infrared Earth sensors are used in spacecraft for attitude sensing. Their accuracy is limited by systematic and random errors. Dominant sources of systematic errors are analyzed for a typical scanning infrared Earth sensor used in a remote-sensing satellite in a 900-km sun-synchronous orbit. The errors considered arise from 1) seasonable variation of infrared radiation, 2) oblate shape of the Earth, 3) ambient temperature of sensors, 4) changes in spin/scan period, and 5) misalignment of the axis of the sensors. Simple relations are derived using least-squares curve fitting for onboard correction of these errors. With these, it is possible to improve the accuracy of attitude determination by eight fold and achieve performance comparable to ground-based post-facto attitude computation.

A multi-station satellite radio beacon study of ionospheric variations during total solar eclipses

Faraday rotation data obtained at Delhi, Kurukshetra, Hyderabad, Bangalore, Waltair, Nagpur and Calcutta during the total solar eclipse of 16 February 1980 and at Delhi during the total solar eclipse of 31 July 1981 have been analysed to detect the gravity waves generated by a total solar eclipse as hypothesized by Chimonas and Hines (1970, J. geophys. Res. 75, 875). It has been found that gravity waves can be generated by a total solar eclipse but their detection at ionospheric heights is critically dependent on the location of the observing station in relation to the eclipse path geometry. The distance of the observing station from the eclipse path should be more than 500 km in order to detect such gravity waves.

Chemical Aspects of the Synergistic Hypergolic Ignition in Hybrid Systems with N2O4 as Oxidizer

Synergistic hypergolic ignition with nitrogen tetroxide ( N2O4) as oxidizer has been observed in hybrid systems comprising of a mixture of magnesium and Schiff bases as fuels. The ignition delays (IDs) measured using a modified device, have been compared with those of magnesium-Schiff base-WFNA systems under identical conditions. The ID has been found to vary with the nature of the substitution in both the benzene rings. A linear relationship emerges when the ignition delays are plotted against the Hammett substitution constants (σ). The preignition products of the reaction of N2O4 with magnesium and benzylidineaniline have been analysed to be Mg(NO3)2, benzenediazonium salt and benzaldehyde. Based on the preignition products isolated, a probable reaction mechanism has been proposed. The previously proposed preignition mechanism for the Schiff base-magnesium-WFNA system has been further supported from the present ignition delay data.

On the investigation of singular stress field around radial cracks in annulii subjected to diametrical tension

Results of photoelastic investigation conducted on annulii containing a radial crack at inner edge and subjected to diametrical tension are reported. The cracks are oriented at 90°, 60° and 45° to the loading direction. The Stress-Intensity Factors (SIFs) were determined by analysing the crack-tip stress fields. Smith and Smith's method [Engng Fracture Mech.4, 357–366 (1972)] and a modified method developed earlier by the authors (to be published) were adopted in the evaluation of SIFs.

Mixed-mode stress intensity factors for arbitrarily oriented cracks in cylindrical shells under axial pull and torsion

Results of photoelastic investigations conducted on cylindrical tubes (made of Araldite material) containing cracks oriented at 0°, 30°, 45°, 60° and 90° to the axis of the tube and subjected to axial and torsional loads are reported. The stress-intensity factors (SIFs) were determined by analysing the crack-tip stress fields. Smith and Smith's method [Engng Fracture Mech.4, 357–366 (1972)] and a new method developed by the authors by modifying Rakesh et al.'s method [Proc. 26th Congress of ISTAM, India (1981)] were employed to evaluate the mixed-mode SIFs.

Transient Momentum and Enthalpy Transfer in Packed Beds at High Reynolds Numbers

An experimental study for transient temperature response and pressure drop in a randomly packed bed at high Reynolds numbers is presented.The packed bed is used as a compact heat exchanger along with a solid-propellant gas generator, to generate room-temperature gases for use in control actuation, air bottle pressurization, etc. Packed beds of lengths 200 and 300 mm were characterized for packing-sphere-based Reynolds numbers ranging from 0.8 x 10(4) to 8.5 x 10(4).The solid packing used in the bed consisted of phi 9.5 mm steel spheres. The bed-to-particle diameter ratio was with the average packed-bed porosity around 0.43. The inlet flow temperature was unsteady and a mesh of spheres was used at either end to eliminate flow entrance and exit effects. Gas temperature and pressure were measured at the entry, exit,and at three axial locations along centerline in the packed beds. The solid packing temperature was measured at three axial locations in the packed bed. A correlation based on the ratio of pressure drop and inlet-flow momentum (Euler number) exhibited an asymptotically decreasing trend with increasing Reynolds number. Axial conduction across the packed bed was found to he negligible in the investigated Reynolds number range. The enthalpy absorption rate to solid packing from hot gases is plotted as a function of a nondimensional time constant for different Reynolds numbers. A longer packed bed had high enthalpy absorption rate at Reynolds number similar to 10(4), which decreased at Reynolds number similar to 10(5). The enthalpy absorption plots can be used for estimating enthalpy drop across packed bed with different material, but for a geometrically similar packing.

Helicopter blade flapping with and without small angle assumption in the presence of dynamic stall

The flapping equation for a rotating rigid helicopter blade is typically derived by considering (1)small flap angle, (2) small induced angle of attack and (3) linear aerodynamics. However, the use of nonlinear aerodynamics such as dynamic stall can make the assumptions of small angles suspect as shown in this paper. A general equation describing helicopter blade flap dynamics for large flap angle and large induced inflow angle of attack is derived. A semi-empirical dynamic stall aerodynamics model (ONERA model) is used. Numerical simulations are performed by solving the nonlinear flapping ordinary differential equation for steady state conditions and the validity of the small angle approximations are examined. It is shown that the small flapping assumption, and to a lesser extent, the small induced angle ofattack assumption, can lead to inaccurate predictions of the blade flap response in certain flight conditions for some rotors when nonlinear aerodynamics is considered. (C) 2010 Elsevier Inc. All rights reserved.

Evaluating ILS and MLS Sites without Flight Tests

Instrument landing systems (ILS) and the upcoming microwave landing systems (MLS) are (or are planned to be) very important navigational aids at most major airports of the world. However, their performance is directly affected by the features of the site in which they are located. Currently, validation of the ILS performance is through costly and time-consuming experimental methods. This paper outlines a powerful and versatile analytical approach for performing the site evaluation, as an alternative to the experimental methods. The approach combines a multi-plate model for the terrain with a powerful and exhaustive ray-tracing technique and a versatile and accurate formulation for estimating the electromagnetic fields due to the array antenna in the presence of the terrain. It can model the effects of the undulation, the roughness and the impedance (depending on the soil type) of the terrain at the site. The results computed from the analytical method are compared with the actual measurements and good agreement is shown. Considerations for site effects on MLS are also outlined.

Environmental effects on fibre-Polymer composites

Moisture absorption characteristics and its effects on the mechanical properties and failure process of polymers (neat epoxy and polyester resins) and composites with simple (glass, carbon and kevlar) and hybrid (glass-carbon, carbon-kevlar and kevlar-glass) fibres were experimentally determined before and after immersion in water at 343 K for 20 days. The maximum moisture content (Mm) and diffusion coefficient (Dx) of these composites were determined. The degradation in ultimate tensile strength and Young's modulus due to the moisture content were experimentally determined and found to be quite significant. Acoustic emissions, from specimens before and after exposure, were monitored during the load cycle, and revealed a significant change in the failure process of these composites. Scanning Electron Microscope (SEM) studies on failed exposed and unexposed specimens revealed resin leach out and fibre prominence.

Direct numerical simulation of autoignition in a non-premixed, turbulent medium

In this paper, direct numerical simulation of autoignition in an initially non-premixed medium under isotropic, homogeneous, and decaying turbulence is presented. The pressure-based method developed herein is a spectral implementation of the sequential steps followed in the predictor-corrector type of algorithms; it includes the effects of density fluctuations caused by spatial inhomogeneities ill temperature and species. The velocity and pressure field are solved in the spectral space while the scalars and density field are solved in the physical space. The presented results reveal that the autoignition spots originate and evolve at locations where (1) the composition corresponds to a small range around a specific mixture fraction, and (2) the conditional scaler dissipation rate is low. A careful examination of the data obtained indicates that the autoignition spots originate in the vortex cores, and the hot gases travel outward as combustion progresses. Hence, the applicability of the transient laminar flamelet model for this problem is questioned. The dependence of autoignition characteristics on parameters such as (1) die initial eddy-turnover time and (2) the initial ratio of length scale of scalars to that of velocities are investigated. Certain implications of new results on the conditional moment closure modeling are discussed.

Thermal decomposition studies on copolyurethanes based on hydroxyl terminated polybutadiene and poly(12-hydroxy stearic acid-co-TMP) ester polyol

Thermal degradation of copolyurethanes based on hydroxyl terminated polybutadiene (HTPB) and poly(12-hydroxy stearic acid-co-TMP) ester polyol (PEP) with varying compositions has been studied by thermo-gravimetric and pyrolysis-GC techniques. The copolyurethanes were found to decompose in multiple stages and the kinetic parameters were found to be dependent on the method of their evaluation. The activation energy for the initial stage of decomposition was found to increase, and for the main stage decreases with the increase in PEP content. The pyrolysis-GC studies on the ammonium perchlorate filled copolyurethanes (solid propellants) showed that the major products during the pyrolysis were C-2, C-3 hydrocarbons and butadiene. The amount of C-2 fraction in the pyrolyslate increased with solid loading, as well as with the HTPB content in the copolyurethanes. A linear relationship apparently exists between the amount of C-2 fraction and the burn rates of the solid propellants. (C) 2000 Elsevier Science Ltd. All rights reserved.

A numerical study of the role of the vertical structure of vorticity during tropical cyclone genesis

An eight-level axisymmetric model with simple parameterizations for clouds and the atmospheric boundary layer was developed to examine the evolution of vortices that are precursors to tropical cyclones. The effect of vertical distributions of vorticity, especially that arising from a merger of mid-level vortices, was studied by us to provide support for a new vortex-merger theory of tropical cyclone genesis. The basic model was validated with the analytical results available for the spin-down of axisymmetric vortices. With the inclusion of the cloud and boundary layer parameterizations, the evolution of deep vortices into hurricanes and the subsequent decay are simulated quite well. The effects of several parameters such as the initial vortex strength, radius of maximum winds, sea-surface temperature and latitude (Coriolis parameter) on the evolution were examined. A new finding is the manner in which mid-level vortices of the same strength decay and how, on simulated merger of these mid-level vortices, the resulting vortex amplifies to hurricane strength in a realistic time frame. The importance of sea-surface temperature on the evolution of full vortices was studied and explained. Also it was found that the strength of the surface vortex determines the time taken by the deep vortex to amplify to hurricane strength.

Three-Dimensional Unsteady Flow With Heat and Mass Transfer Over a Continuous Stretching Surface

A numerical solution of the unsteady boundary layer equations under similarity assumptions is obtained. The solution represents the three-dimensional unsteady fluid motion caused by the time-dependent stretching of a flat boundary. It has been shown that a self-similar solution exists when either the rate of stretching is decreasing with time or it is constant. Three different numerical techniques are applied and a comparison is made among them as well as with earlier results. Analysis is made for various situations like deceleration in stretching of the boundary, mass transfer at the surface, saddle and nodal point flows, and the effect of a magnetic field. Both the constant temperature and constant heat flux conditions at the wall have been studied.