Numerical procedure for second order non-linear ordinary differential equations and application to heat transfer problem

In this paper, we have first given a numerical procedure for the solution of second order non-linear ordinary differential equations of the type y″ = f (x;y, y′) with given initial conditions. The method is based on geometrical interpretation of the equation, which suggests a simple geometrical construction of the integral curve. We then translate this geometrical method to the numerical procedure adaptable to desk calculators and digital computers. We have studied the efficacy of this method with the help of an illustrative example with known exact solution. We have also compared it with Runge-Kutta method. We have then applied this method to a physical problem, namely, the study of the temperature distribution in a semi-infinite solid homogeneous medium for temperature-dependent conductivity coefficient.

Dispersion of the stress-optic coefficient in glasses

The stress-optic coefficient (n3/2)(q11-q12) has been determined for a series of 18 optical glasses of different compositions in the wavelength range 5700-3200 Å. The coefficients are negative for all the glasses except for a high-lead-content glass of density 6·7 and refractive index 1·89. The numerical value of the coefficient decreases as one proceeds to the ultraviolet. This behaviour is just the opposite of what is observed in fused silica. By applying Mueller's theory, the strain polarizability constant and its dispersion have been evaluated.

Crystalline ethane-1,2-diol does not have intra-molecular hydrogen bonding: Experimental and theoretical charge density studies

The X-ray structure and electron density distribution of ethane-1,2-diol (ethylene glycol), obtained at a resolution extending to 1.00 Å−1 in sin θ/λ (data completion = 100% at 100 K) by in situ cryocrystallization technique is reported. The diol is in the gauche (g′Gt) conformation with the crystal structure stabilised by a network of inter-molecular hydrogen bonds. In addition to the well-recognized O–H···O hydrogen bonds there is topological evidence for C–H···O inter-molecular interactions. There is no experimental electron density based topological evidence for the occurrence of an intra-molecular hydrogen bond. The O···H spacing is not, vert, similar0.45 Å greater than in the gas-phase with an O–H···O angle close to 90°, calling into question the general assumption that the gauche conformation of ethane-1,2-diol is stabilised by the intra-molecular oxygen–hydrogen interaction.

Effect of relative density and confining pressure on Poisson ratio from bender and extender elements tests

By using the bender and extender elements tests, together with measurements of the travel times of shear (S) and primary (P) waves, the variation of Poisson ratio (nu) was determined for dry sands with respect to changes in relative densities and effective confining pressures (sigma(3)). The tests were performed for three different ranges of particle sizes. The magnitude of the Poisson ratio decreases invariably with an increase in both the relative density and the effective confining pressure. The effect of the confining pressure on the Poisson ratio was found to become relatively more significant for fine-grained sand as compared with the coarse-grained sand. For a given material, at a particular value of sigma(3), the magnitude of the Poisson ratio decreases, almost in a linear fashion, with an increase in the value of maximum shear modulus (G(max)). The two widely used correlations in literature, providing the relationships among G(max), void ratio (e) and effective confining pressure (sigma(3)), applicable for angular granular materials, were found to compare reasonably well with the present experimental data for the fine- and medium-grained sands. However, for the coarse-grained sand, these correlations tend to overestimate the values of G(max).

Nonlinear vibration analysis of composite laminated and sandwich plates with random material properties

Nonlinear vibration analysis is performed using a C-0 assumed strain interpolated finite element plate model based on Reddy's third order theory. An earlier model is modified to include the effect of transverse shear variation along the plate thickness and Von-Karman nonlinear strain terms. Monte Carlo Simulation with Latin Hypercube Sampling technique is used to obtain the variance of linear and nonlinear natural frequencies of the plate due to randomness in its material properties. Numerical results are obtained for composite plates with different aspect ratio, stacking sequence and oscillation amplitude ratio. The numerical results are validated with the available literature. It is found that the nonlinear frequencies show increasing non-Gaussian probability density function with increasing amplitude of vibration and show dual peaks at high amplitude ratios. This chaotic nature of the dispersion of nonlinear eigenvalues is also r

Fatigue crack propagation in rocker and roller–rocker bearings of railway steel bridges

In this work, a method is proposed for rolling contact fatigue crack propagation analysis using contact and fracture theories in conjunction with fatigue laws. The proposed method is used in the fatigue analysis of rocker and roller–rocker bearings of a railway open web girder bridge which is instrumented with strain gages. Using a contact algorithm based on the minimum energy principle for bodies in rolling contact with dry friction, the normal and tangential pressure distribution are computed. It is seen that the most critical location of a crack in bearings is at a point very close to the contact region, as expected.

Fatigue crack propagation in rocker and roller-rocker bearings of railway steel bridges

In this work, a method is proposed for rolling contact fatigue crack propagation analysis using contact and fracture theories in conjunction with fatigue laws. The proposed method is used in the fatigue analysis of rocker and roller-rocker bearings of a railway open web girder bridge which is instrumented with strain gages. Using a contact algorithm based on the minimum energy principle for bodies in rolling contact with dry friction, the normal and tangential pressure distribution are computed. It is seen that the most critical location of a crack in bearings is at a point very close to the contact region, as expected. (C) 2010 Elsevier Ltd. All rights reserved.

Effect of copper addition on the fracture and fatigue crack growth behavior of solution heat-treated SUS 304H austenitic steel

Small additions of Cu to the SUS 304H, a high temperature austenitic stainless steel, enhance its high temperature strength and creep resistance. As Cu is known to cause embrittlement, the effect of Cu on room temperature mechanical properties that include fracture toughness and fatigue crack threshold of as-solutionized SUS 304H steel were investigated in this work. Experimental results show a linear reduction in yield and ultimate strengths with Cu addition of up to 5 wt.% while ductility drops markedly for 5 wt.% Cu alloy. However, the fracture toughness and the threshold stress intensity factor range for fatigue crack initiation were found to be nearly invariant with Cu addition. This is because the fracture in this alloy is controlled by the debonding from the matrix of chromium carbide precipitates, as evident from fractography. Cu, on the other hand, remains either in solution or as nano-precipitates and hence does not influence the fracture characteristics. It is concluded that small additions of Cu to 304H will not have adverse effects on its fracture and fatigue behavior. (C) 2010 Elsevier B.V. All rights reserved.

Influence of propellant choice on MPD arcjet cathode surface current density distribution

The radial current density on an MPD arcjet cathode surface is theoretically investigated for five propellants. It is found that excessive current concentration at the upstream end of the cathode occurs in the case of hydrogen. This undesirable effect is traced to the higher electrical conductivity of hydrogen plasma.

A finite element analysis of quasistatic crack growth in a pressure sensitive constrained ductile layer

Polymeric adhesive layers are employed for bonding two components in a wide variety of technological applications, It has been observed that, unlike in metals, the yield behavior of polymers is affected by the state of hydrostatic stress. In this work, the effect of pressure sensitivity of yielding and layer thickness on quasistatic interfacial crack growth in a ductile adhesive layer is investigated. To this end, finite deformation, finite element analyses of a cracked sandwiched layer are carried out under plane strain, small-scale yielding conditions for a wide range of mode mixities. The Drucker-Prager constitutive equations are employed to represent the behavior of the layer. Crack propagation is simulated through a cohesive zone model, in which the interface is assumed to follow a prescribed traction-separation law. The results show that for a given mode mixity, the steady state Fracture toughness [K](ss) is enhanced as the degree of pressure sensitivity increases. Further, for a given level of pressure sensitivity, [K](ss) increases steeply as mode Il loading is approached. (C) 2000 Elsevier Science Ltd. All rights reserved.

Charge Density Analysis of Crystals of Nicotinamide with Salicylic Acid and Oxalic Acid: An Insight into the Salt to Cocrystal Continuum

Charge density analysis from both experimental and theoretical points of view on two molecular complexes: one is formed between nicotinamide and salicylic acid, and the other formed between nicotinamide and oxalic acid brings out the quantitative topological features to distinguish a cocrystal from a salt.

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.

Topological analysis off charge density distribution in concomitant polymorphs of 3-acetylcoumarin, a case off packing polymorphism

Detailed investigation of the charge density distribution in concomitant polymorphs of 3-acetylcoumarin in terms of experimental and theoretical densities shows significant differences in the intermolecular features when analyzed based on the topological properties via the quantum theory of atoms in molecules. The two forms, triclinic and monoclinic (Form A and Form B), pack in the crystal lattice via weak C-H---O and C-H---pi interactions. Form A results in a head-to-head molecular stack, while Form B generates a head-to-tail stack. Form A crystallizes in PI (Z' = 2) and Form B crystallizes in P2(1)/n (Z = 1). The electron density maps of the polymorphs demonstrate the differences in the nature of the charge density distribution in general. The charges derived from experimental and theoretical analysis show significant differences with respect to the polymorphic forms. The molecular dipole moments differ significantly for the two forms. The lattice energies evaluated at the HF and DFT (B3LYP) methods with 6-31G** basis set for the two forms clearly suggest that Form A is the thermodynamically stable form as compared to Form B. Mapping of electrostatic potential over the molecular surface shows dominant variations in the electronegative region, which bring out the differences between the two forms.