Nonlinear asymptotic theory of hypersonic flow past a circular cone

Summary The hypersonic small-disturbance theory is reexamined in this study. A systematic and rigorous approach is proposed to obtain the nonlinear asymptotic equation from the Taylor-Maccoll equation for hypersonic flow past a circular cone. Using this approach, consideration is made of a general asymptotic expansion of the unified supersonic-hypersonic similarity parameter together with the stretched coordinate. Moreover, the successive approximate solutions of the nonlinear hypersonic smalldisturbance equation are solved by iteration. Both of these approximations provide a closed-form solution, which is suitable for the analysis of various related flow problems. Besides the velocity components, the shock location and other thermodynamic properties are presented. Comparisons are also made of the zeroth-order with first-order approximations for shock location and pressure coefficient on the cone surface, respectively. The latter (including the nonlinear effects) demonstrates better correlation with exact solution than the zeroth-order approximation. This approach offers further insight into the fundamental features of hypersonic small-disturbance theory.Notation a speed of sound - H unified supersonic-hypersonic similarity parameter, - K hypersonic similarity parameter, M - M freestream Mach number - P pressure - T temperature - S entropy - u, v radial, polar velocities - V freestream velocity - shock angle - cone angle - density - density ratio, /() - ratio of specific heats - polar angle - stretched polar angle, / - (), (), () gage functions     

Thermal Conductivity of the Four-Leg Spin-Ladder System La2Cu2O5 Single Crystal

We have investigated the magnetic susceptibility, , and the thermal conductivity, , in magnetic fields for the four-leg spin-ladder system La₂Cu₂O₅ single crystal. The in a magnetic field parallel to the ladder exhibits a kink at 130 K in correspondence to the magnetic ordering. The along the ladder exhibits a peak at 25 K and a shoulder at 14 K, which are probably related to the thermal conductivity due to magnons, _magnon, and that due to phonons, _phonon, respectively. The perpendicular to the ladder, on the other hand, exhibits only one broad peak related to _phonon. The observed large anisotropy of has been explained based upon the anisotropy of _magnon.     

Simulation of the Process of Water Freezing in a Round Pipe

The problem of freezing of pure water in a round pipe is treated with due regard for convection under asymmetric thermal boundary conditions in the absence of motion along the pipe. The problem is solved numerically using the control volume approach, SIMPLER algorithm, and the enthalpy method. Results are obtained for three Grashof (Gr) and six Biot (Bi) numbers: Gr = 1.55 × 10⁶, Bi = 0.305 (0 < ), Bi = 0.044 ( < 2); Gr = 1.24 × 10⁷, Bi = 0.610 (0 < ), Bi = 0.087 ( < 2); Gr = 9.89 × 10⁷, Bi = 1.220 (0 < ), Bi = 0.174 ( < 2). The correctness of calculation of the problem disregarding free-convection flows is analyzed.     

Thermal transport properties of helium near the superfluid transition. II. Dilute3He-4He mixtures in the superfluid phase

Measurements of the average thermal conductivity _exp hQ/T and of the thermal relaxation time to reach steady-state equilibrium conditions are reported in the superfluid phase for dilute mixtures of³He in⁴He. Hereh is the cell height,Q is the heat flux, andT is the temperature difference across the fluid layer. The measurements were made over the impurity range 2×10^–9<X(³He)<3×10^–2 and with heat fluxes 0.3<Q<160 µW/cm². Assuming the boundary resistanceR _b , measured forX<10^–5, to be independent ofX over the whole range ofX, a calculation is given for _exp. ForQ smaller than a well-defined critical heat fluxQ _c (X) X ^0.9, _exp is independent of Q and can be identified with the local conductivity _eff, which is found to be independent of the reduced temperature = (T–T)/T for –10^–2. Its extrapolated value at T is found to depart forX10^–3 from the prediction X ^–1 , tending instead to a weaker divergence X ^–a witha0.08. A finite conductivity asX tends to zero is not excluded by the data, however. ForQ >Q _c (X), a nonlinear regime is entered. ForX10^–6, the measurements with the available temperature resolution are limited to the nonlinear conditions, but can be extrapolated into the linear regime forX2×10^–7. The results for _exp(Q),Q _c (X), and _eff(XX) are found to be internally consistent, as shown by comparison with a theory by Behringer based on Khalatnikov's transport equations. Furthermore, the observed relaxation times (X) in the linear regime are found to be consistent forX>10^–5 with the hydrodynamic calculations using the measured _eff(X). ForX<10^–5, a faster relaxation mechanism than predicted seems to dominate. The transport properties in the nonlinear regimes are presented and unexplained observations are discussed.     

Fatigue crack growth behaviour of tungsten carbide-cobalt hardmetals

Fatigue crack propagation studies have been carried out on a range of WC-Co hardmetals of varying cobalt content and grain size using a constant-stress intensity factor double torsion test specimen geometry. Results have confirmed the marked influence of mean stress (throughK _max), which is interpreted in terms of static modes of fracture occurring in conjunction with a true fatigue process, the existence of which can be rationalized through the absence of any frequency effect. Dramatic increases in fatigue crack growth rate are found asK _max approaches that value of stress intensity factor ( 0.9K_IC) for which static crack growth under monotonic load (or static fatigue) occurs in these materials. Lower crack growth rates, however, produce fractographic features indistinguishable from those resulting from fast fracture. These observations, and the important effect of increasing mean free path of the cobalt binder in reducing fatigue crack growth rate, can reasonably be explained through a consideration of the mechanism of fatigue crack advance through ligament rupture of the cobalt binder at the tip of a propagating crack.     

The electrical conductivity of some azine compounds

Correlation of the molecular structure and thermal study of the electrical conductivity of benzaldazine (BA) and its NO₂-derivatives were carried out. The values of the activation energies, obtained from the electrical conductivity measurements, as well as the energy gaps calculated from ultra violet (u.v.) and visible spectra obtained either in liquid or in solid forms indicate that all compounds behave like a semiconducting material. Both n*s and *s were found to be the main sources participatiog in the cooduction processes.     

On internal stress and activation volume in polymers

The two-site model is developed for the analysis of stress relaxation data. It is shown that the product of d In (– )/d and (- _i) is constant where is the applied stress, _i is the (deformation-induced) internal stress and = d/dt. The quantity d In ( )/d is often presented in the literature as the (experimental) activation volume, and there are many examples in which the above relationship with (- _i) holds true. This is in apparent contradiction to the arguments that lead to the association of the quantity d In (– )/d with the activation volume, since these normally start with the premise that the activation volume is independent of stress. In the modified theory presented here the source of this anomaly is apparent. Similar anomalies arise in the estimation of activation volume from creep or constant strain rate tests and these are also examined from the standpoint of the site model theory. In the derivation presented here full account is taken of the site population distribution and this is the major difference compared to most other analyses. The predicted behaviour is identical to that obtained with the standard linear solid. Consideration is also given to the orientation-dependence of stress-aided activation.     

Nonlinear axisymmetric response of orthotropic thin truncated conical and spherical caps

Summary An approximate analytical solution of the large deflection axisymmetric response of polar orthotropic thin truncated conical and spherical shallow caps is presented. Donnell type equations are employed. The deflection is approximated by a one term mode shape satisfying the boundary conditions. The Galerkin's method is used to get the governing equation for the deflection at the hole. Nonlinear free vibration response and the response under uniformly distributed static and step function loads are obtained. The effect of various parameters is investigated.Notations A, A ^* Inward and outward amplitudes - a, b, h Base radius, inner radius and thickness of the cap - D M h ³/[12(–v ² )] - E ,E Young's moduli - H ^*,H Apex height, dimensionless apex heght:H ^*/h - N _, Stress resultants - p ^1/2 - q Uniformly distributed load - Q,Q₀ Dimensionless load: , dimensionless step load - Q, Q ₀ Dimensionless load: , step load - t, Time, dimensionless time: t - T _A Ratio of nonlinear periodT for inward amplitudeA and the linear periodT _L - w ^* Normal displacement at middle surface - w Dimensionless displacement:w ^*/h - ₁ Linear parameter of static response - Orthotropic Parameter:E /E - Mass density - ₂,₃ Quadratic and cubic nonlinearity parameters - b/a - v ,v Poisson's ratios - Dimensionless radius:r/a - ^*, Stress function, dimensionless stress function: - ₀ ^* ,₀ Linnear frequency, dimensionless frequency: With 7 Figures     

Optical thermophysical characteristics of the condensed phase of combustion products and heatproof materials

An unconventional unit for experimentally studying the optical thermophysical properties of materials over a wide temperature range is described. Results are presented of studying the temperature function and dispersion of the absorption index for the condensed phase of the combustion products of a metalbearing fuel and of the emittance of fiber fireproof materials.Notation K spectral absorption coefficient - D spectral transmission coefficient - spectral absorption index - wavelength - d sample thickness - d_s surface layer thickness - I spectral intensity of radiation - , spectral and total emittances Kirov Polytechnical Institute. Kazan Chemical Engineering Institute, Russia. Translated from Inzhenerno-Fizicheskii Zhurnal, Vol. 64, No. 3, pp. 330–336. March, 1993     

Numerical simulation of an explosive plasma generator in the gasdynamic approximation

A two-dimensional gasdynamic model of a plasma generator is proposed. A numerical solution of the problem is obtained, and peculiarities of the gasdynamic flow are considered. Results are compared with experiment.Notation T temperature - p pressure - N number of gas particles per unit volume - w gas mass velocity vector with components w_z=u and w_r=v - z axial coordinate - r radial coordinate - t time - E total specific energy of gas - specific internal energy - density - V plate velocity - E_in initial plate energy - M plate mass - S plate area - p pressure difference between left and right sides of plate - t, z, r time and space steps - k<1 Courant number - c velocity of sound in the gas Indices 0 initial value - * characteristic dimensional quantities - i, j grid cell indices along z and r - n number of time step - () symbol denoting intermediate values of gasdynamic variables in a time layer Translated from Inzhenerno-Fizicheskii Zhurnal, Vol. 37, No. 5, pp. 859–867, November, 1979.     

A Monte Carlo Study of Disorder in the High Pressure Mixed Crystal N2-Ar

In order to obtain a better notion of the experimental results in our laboratory, Monte Carlo calculations have been performed of the N ₂-Ar crystal on the N ₂-rich side, in the p-T region where the and phases exist in pure N ₂. Considering the enthalpy, the system prefers the Ar atoms to be located on the sphere positions. The * phase is present for mixtures down to but is most likely metastable. The *-* transition shifts to lower temperatures with decreasing . The 2 ^nd order transition within the phase continues to exist to even smaller . In contrast to the * -* transition, the transition temperature for the 2 ^nd order transition does not shift to lower temperatures. For a mixture of it is within 5 K from the pure 2 ^nd order transition at a pressure of 7.0 GPa.     

Flux motion in anisotropic type-II superconductors nearH c 2 with arbitrary vortex orientation

Flux motion in anisotropic type-II superconductors near H_{c 2} is studied in the framework of the time-dependent Ginzburg-Landau (TDGL) theory for the case that the average flux densityB is oriented in an arbitrary direction relative to the principal axes of the sample. The linearized TDGL equation for a uniformly translating order parameter is solved and expressions for all elements of the flux-flow resistivity tensor _ij (including the off-diagonal Hall elements) are obtained. The diagonal elements _ii show the angular scaling property that _ii(B) = _ii(B/H_{c 2}(, ø)), whereas the Hall elements _ij (i j) have additional angular dependences that are not contained in H_c2. For the case that the normal state Hall elements _ij ⁽ⁿ⁾ B _k with i j k, the ratios _ij/ _ij ⁽ⁿ⁾ are functions of B/H_c2 (, ) only.     

Thermal conductivity of argon and argon-neon mixtures at high temperature

The thermal conductivity of neon and argon-neon mixtures is studied in the temperature range 400–1500°K. This is the first recording of such data above 793°K.Notation T_g true temperature drop in gas layer °K - ¯T temperature, °K, Q, effective thermal flux, W - Q_t, Qr thermal flux transmitted by thermal conductivity and radiation, respectively, W - T_sh correction for temperature shift, % - thermal conductivity of gas mixture, W/m·°K - x_i a molar concentration of neon - _i, _i, _ij, _ij potential function parameters for inter molecular interaction of homogeneous and inhomogeneous molecules - slope of exponential repulsion term Translated from Inzhenerno-Fizicheskii Zhurnal, Vol. 33, No. 5, pp. 848–856, November, 1977.     

Second-order wave force on a large cylinder

Summary A new technique has been developed for estimating the wave loadings on large circular cylinders. The theory, mainly due to Lighthill, has been applied to the case of large cylinders. Comparison with the previously reported experimental results shows favourable agreement for the range of parameters indicated in the graphs.Notations The following symbols are used in this paper a Wave amplitude - b Radius of the cylinder - D Diameter of the cylinder - F _l Linear force - F _d Dynamic force - F _q Quadratic force - F _w Water line force - g Acceleration due to gravity - h Depth of water - H=2a Total wave height - k Wave number,k=2/L - L Wave length - n Outward normal to body surface - n _x Direction cosine between the normaln and the given force - S Body surface - t Time - T Wave period - x, y, z Cartesian coordinates - r, ,z Cylindrical coordinates - W Vertical velocity of fluid - Density of fluid - Wave frequency, =2/T - Total velocity potential - _r , _z Partial derivatives with respect tor andz, respectively - ⁽¹⁾, _l Linear theory velocity potential - ⁽²⁾, _q Second-order velocity potential - Water surface elevation or wave height - _i Linear theory wave height - _q Quadratic theory wave height With 3 Figures     

Effect of fibre misalignment on fracture behaviour of fibre-reinforced composites

When a matrix crack encounters a fibre that is inclined relative to the direction of crack opening, geometry requires that the fibre flex is bridging between the crack faces. Conversely, the degree of flexing is a function of the crack face separation, as well as of (1) the compliance of the supporting matrix, (2) the crossing angle, (3) the bundle size, and (4) the shear coupling of the fibre to the matrix. At some crack face separation the stress level in the fibre bundle will cause it to fail. Other bundles, differing in size and orientation, will fail at other values of the crack separation. Such bridging contributes significantly to the resistance of the composite to crack propagation and to ultimate failure. The stress on the composite needed to produce a given crack face separation is inferred by analysing the forces and displacements involved. The resulting model computes stress versus crack-opening behaviour, ultimate strengths, and works of failure. Although the crack is assumed to be planar and to extend indefinitely, the model should also be applicable to finite cracks.Glossary of Symbols a radius of fibre bundle - C 2 _f /aE _f - ^* critical failure strain of fibre bundle - _b bending strain in outer fibre of a bundle - _c background strain in composite - _f axial strain in fibre - _s strain in fibre bundle due to fibre stretching = _f - () strain in composite far from crack - E Young's modulus of fibre bundle - E _c Young's modulus of composite - E _f Young's modulus of fibre - E _m Young's modulus of matrix - f() number density per unit area of fibres crossing crack plane in interval to + d - F total force exerted by fibre bundle normal to crack plane - F _s component of fibre stretching force normal to crack plane - F _b component of bending force normal to crack plane - G _m shear modulus of matrix - h crack face opening relative to crack mid-point - h _m matrix contraction contribution to h - h _f fibre deformation contribution to h - h _max crack opening at which bridging stress is a maximum - I moment of inertia of fibre bundle - k fibre stress decay constant in non-slip region - k ₀ force constant characterizing an elastic foundation (see Equation 7) - L exposed length of bridging fibre bundle (see Equation 1a) - L _f half-length of a discontinuous fibre - m, n parameters characterizing degree of misalignment - N number of bundles intersecting a unit area of crack plane - P _b bending force normal to bundle axis at crack midpoint - P _s stretching force parallel to bundle axis in crack opening - Q() distribution function describing the degree of misalignment - s _f fibre axial tensile stress - s _f ^* fibre tensile failure stress - S stress supported by totality of bridging fibre bundles - S _max maximum value of bridging stress - v fibre displacement relative to matrix - v elongation of fibre in crack bridging region - u _coh non-slip contribution to fibre elongation - U fibre elongation due to crack bridging - v overall volume fraction of fibres - v _f volume fraction of bundles - v _m volume fraction matrix between bundles - w transverse deflection of bundle at the crack mid-point - x distance along fibre axis, origin defined by context - X distance between the end of discontinuous fibre and the crack face - X ^* threshold (minimum) value of X that results in fibre failure instead of complete fibre pullout - y displacement of fibre normal to its undeflected axis - Z() area fraction angular weighting function - tensile strain in fibre relative to applied background strain - ^* critical value of to cause fibre/matrix debonding - angle at which a fibre bundle crosses the crack plane - (k ₀/4EI)^1/4, a parameter in cantilever beam analysis - v_m Poisson's ratio of matrix - L (see Equation 9) - shear stress - _* interlaminar shear strength of bundle - _d fibre/matrix interfacial shear strength - _f frictional shear slippage stress at bundle/matrix interface - angular deviation of fibre bundle from mean orientation of all bundles - angle between symmetry axis and crack plane     

Two singular points finite elements in the analysis of kinked cracks

A formulation for incorporating two singular points (TSP) of variable orders in a single finite element is presented. Though the element does not satisfy any of the convergence criterion, its performance is found to be good, which has been demonstrated by considering number of examples on kinked cracks. In each case only one such element is incorporated in the whole discretization. These examples illustrate the usefulness of the element to analyse kinked cracks of various sizes and shapes and subjected to different loading and boundary conditions. Computed J-integrals are compared with analytical solutions, wherever possible, and the accuracy appears quite good. Effect of size of the element on, and the path independence of, J are also examined.List of symbols , conventional natural coordinate system - ₁, ₂ another elemental natural coordinate system - L _ij equation of side ij of an element - A _ij , B _ij constants - N _i shape function associated node i - ₁, ₂ constants associated with the order of singularities - , an elemental natural coordinate system - u _i , v _i displacement components in the Cartesian directions     

Phase decomposition in extruded Zn-Al based alloy

Ageing characteristics of an extruded eutectoid Zn-Al based alloy were investigated using X-ray diffraction and scanning electron microscopy techniques. The extruded alloy consisted of Al rich phase and Zn rich _E and phases. The original cast eutectoid Zn-Al alloy was extruded at 250 °C. Both supersaturated _s and _s phase decomposed during extrusion and appeared as fine and coarse lamellar structures. The _E and phases particles formed in the original interdendritic region. It was found that two Zn rich phases _E and decomposed sequentially during ageing at 170, 140 °C. The decomposition of the _E phase occurred as a discontinuous precipitation in the early stage of ageing and the decomposition of the phase took place in a four phase transformation: + T + in the prolonged ageing. Two typical morphologies of the decomposition of the Zn rich phases _E and were distinctive in back-scattered scanning electron microscopy.     

Thermal conductivity of methane

This paper reports thermal conductivity data for methane measured in the temperature range 120–400 K and pressure range 25–700 bar with a maximum uncertainty of ± 1%. A simple correlation of these data accurate to within about 3% is obtained and used to prepare a table of recommended values.Nomenclature a _k ,b _ij ,b _k Parameters of the regression model, k= 0 to n; i =0 to m; j =0 to n - P Pressure (MPa or bar) - Q _kl Heat flux per unit length (mW · m^–1) - t time (s) - T Temperature (K) - T _cr Critical temperature (K) - T _r reduced temperature (= T/T _cr) - T _w Temperature rise of wire between times t ₁ and t ₂ (deg K) - T ^* Reduced temperature difference (TT _cr)/T _cr - Thermal conductivity (mW · m^–1 · K^–1) - ₁ Thermal conductivity at 1 bar (mW · m^–1 · K^–1) - _bg Background thermal conductivity (mW · m^–1 · K^–1) - _cr Anomalous thermal conductivity (mW · m^–1 · K^–1) - _e Excess thermal conductivity (mW · m^–1 · K^–1) - Density (g · cm^–3) - _cr Critical density (g · cm^–3) - _r Reduced density (= / _cr) - ^* Reduced density difference (– _cr )/ _cr      

Phase decomposition of liquid-quenched β-type Ti-Cr alloys

Phase decomposition behaviour of liquid-quenched (bcc) type Ti-Cr alloys was investigated by means of transmission electron microscopy and hardness measurements. It was found that decomposition of to ₁ (Ti-rich, bcc) + ₂ (Ti-lean, bcc) takes place in the intermediate composition range of the Ti-Cr system. This experimental result proves the theoretical prediction made by Menon and Aaronson, but the observed ₁ + ₂ two-phase field expands towards higher temperatures than the predicted binodal line. The coherent ₁ + ₂ two-phase state exhibits the so-called 100 modulated structure and it was concluded that the formation of such a structure is a result of spinodal decomposition of the -phase. We obtained time-temperature-transformation (TTT) diagrams of -type Ti-30, 40 and 50 at % Cr alloys. A typical sequence of structural change is coherent ₁ + ₂ incoherent ₁ + ₂ incoherent ₁ + ₂ + grain boundary precipitates stable state of + TiCr₂ or + TiCr₂. Not all the states in the above sequence appear, depending on alloy composition, liquid-quenching rate and ageing temperature.