Determination of optimal trajectory under design constraints for a satellite launch vehicle

An optimal pitch steering programme of a solid-fuel satellite launch vehicle to maximize either (1) the injection velocity at a given altitude, or (2) the size of circular orbit, for a given payload is presented. The two-dimensional model includes the rotation of atmosphere with the Earth, the vehicle's lift and drag, variation of thrust with time and altitude, inverse-square gravitational field, and the specified initial vertical take-off. The inequality constraints on the aerodynamic load, control force, and turning rates are also imposed. Using the properties of the central force motion the terminal constraint conditions at coast apogee are transferred to the penultimate stage burnout. Such a transformation converts a time-free problem into a time-fixed one, reduces the number of terminal constraints, improves accuracy, besides demanding less computer memory and time. The adjoint equations are developed in a compact matrix form. The problem is solved on an IBM 360/44 computer using a steepest ascent algorithm. An illustrative analysis of a typical launch vehicle establishes the speed of convergence, and accuracy and applicability of the algorithm.

Interface between 2 dissimilar fluids beneath a sheetpile cofferdam

A closed form solution is presented for determining the shape and location of the interface between two dissimilar fluids (having different densities) when steady flow takes place through a homogeneous and isotropic porous medium, into a sheetpile cofferdam; the interface is assumed to be sharp and the lower fluid stationary. The solution is obtained using the inverse hodograph. Numerical results are presented in nondimensional form for various parametric conditions in the physical plane; the interface pattern, as also the seepage discharge and exit gradient distribution are shown. The critical conditions of the interface are studied.

Rank-Augmented LU-Algorithm for Computing Generalized Matrix Inverses

A rank-augmnented LU-algorithm is suggested for computing a generalized inverse of a matrix. Initially suitable diagonal corrections are introduced in (the symmetrized form of) the given matrix to facilitate decomposition; a backward-correction scheme then yields a desired generalized inverse.

Dynamic viscosity of ternary mixtures of dichlorodifluoromethane (R-12), chlorodifluoromethane (R-22), and dichlorotetrafluoroethane (R-114) vapors

The viscosities of ternary mixtures of R-12, R-22, and R-114 vapors were determined at ambient temperature and pressure within +-1% by using an oscillating disk viscometer. The empirical viscosity obtained by Wllke's equation compares very well with the experimental results obtained with this vlscometer. In the case of this ternary vapor mixture, as long as the molar fraction ratio of R-12 to R-114 Is maintained at approximately 2"' (=Inverse ratio of thelr molecular weights) the viscosity of the ternary mixture at ambient temperature and pressure remalns constant irrespective of the percentage of R-22 present in the mixture.

Modeling of Freezing Phenomena Induced by Chemical Reactions

Spatial variations in the concentration of a reactive solute in solution are often encountered in a catalyst particle, and this leads to variation in the freezing point of the solution. Depending on the operating temperature, this can result in freezing of the solvent oil a portion of catalyst, rendering that part of the active area ineffective Freezing call occur by formation of a sharp front or it mush that separates the solid and fluid phases. In this paper, we model the extent of reduction in the active area due to freezing. Assuming that the freezing point decreases linearly with solute concentration, conditions for freezing to occur have been derived. At steady state, the ineffective fraction of catalyst pellet is found to be the same irrespective of the mode of freezing. Progress of freezing is determined by both the heat of reaction and the latent heat of fusion Unlike in freezing of alloys where the latter plays a dominant role, the exothermicity of the reaction has a significant effect on freezing in the presence of chemical reactions. A dimensionless group analogous to the Stefan number could be defined to capture the combined effect of both of these.

Which One Is Preferred: Myers-Saito Cyclization of Ene-Yne-Allene or Garratt-Braverman Cyclization of Conjugated Bisallenic Sulfone? A Theoretical and Experimental Study

A competitive scenario between Myers-Saito (MS) and Garraff-Braverman (GB) cyclization has been created in a molecule. High-level computations indicate a preference for GB over MS cyclization. The activation energies for the rate-determining steps of the GB and MS cyclizations were found to be the same (24.4 kcal/mol) at the B3LYP/6-31G* level of theory; thus, from the kinetic point of view, both reactions are feasible. However, the main biradical intermediate GB2 of the GB reaction is 6.2 kcal/mol lower in energy than the biradical MS2, which is the main intermediate of MS reaction, so GB cyclization is thermodynamically favored over MS cyclization. To verify the prediction by computational techniques, bisenediynyl sulfones 1-4 and bisenediynyl sulfoxide 17 were synthesized. Under basic conditions, these molecules isomerized to a system possessing both the ene-yne-allene and the bisallenic sulfone. The isolation of only one product, identified as the corresponding naphthalene- or benzene-fused sulfone 8-11, indicated the occurrence of GB cyclization as the sole reaction pathway. No product corresponding to the MS cyclization pathway could be isolated. Though the theoretical prediction showed a preference for the GB pathway over the MS pathway, the exclusive preference for GB over MS cyclization is very striking. Further analysis showed that the intramolecular self-quenching nature of the GB pathway may play an important role in the complete preference for this reaction. Apart from the mechanistic studies, these sulfones showed DNA cleavage activity that had an inverse relation with the reactivity order. Our findings are important for the design of artificial DNA-cleaving agents.

Soft gluon k(t)-resummation and the Froissart bound

We study soft gluon k(t)-resurnmation and the relevance of InfraRed (IR) gluons for the energy dependence of total hadronic cross-sections. In our model, consistency with the Froissart bound is directly related to the ansatz that the IR behaviour of the QCD coupling constant follows an inverse power law.

Transition from radiatively inefficient to cooling dominated phase in two temperature accretion disks around black holes

We investigate the transition of a radiatively inefficient phase of a viscous two temperature accreting flow to a cooling dominated phase and vice versa around black holes. Based on a global sub-Keplerian accretion disk model in steady state, including explicit cooling processes self-consistently, we show that general advective accretion flow passes through various phases during its infall towards a black hole. Bremsstrahlung, synchrotron and inverse Comptonization of soft photons are considered as possible cooling mechanisms. Hence the flow governs a much lower electron temperature similar to 10(8) - 10(9.5) K compared to the hot protons of temperature similar to 10(10.2) - 10(11.8) K in the range of the accretion rate in Eddington units 0.01 less than or simiar to (M) over dot less than or similar to 100. Therefore, the solutions may potentially explain the hard X-rays and the gamma-rays emitted from AGNs and X-ray binaries. We finally compare the solutions for two different regimes of viscosity and conclude that a weakly viscous flow is expected to be cooling dominated compared to its highly viscous counterpart which is radiatively inefficient. The flow is successfully able to reproduce the observed minosities of the under-fed AGNs and quasars (e.g. Sgr A*), ultra-luminous X-ray sources (e.g. SS433), as well as the highly luminous AGNs and ultra-luminous quasars (e.g. PKS 0743-67) at different combinations of the mass accretion rate and ratio of specific heats.

Two temperature viscous accretion flows around rotating black holes: Description of under-fed systems to ultra-luminous X-ray sources

We discuss two temperature accretion disk flows around rotating black holes. As we know that to explain observed hard X-rays the choice of Keplerian angular momentum profile is not unique, we consider the sub-Keplerian regime of the disk. Without any strict knowledge of the magnetic field structure, we assume the cooling mechanism is dominated by bremsstrahlung process. We show that in a range of Shakura-Sunyaev viscosity parameter 0.2 greater than or similar to alpha greater than or similar to 0.0005, flow behavior varies widely, particularly by means of the size of disk, efficiency of cooling and corresponding temperatures of ions and electrons. We also show that the disk around a rotating black hole is hotter compared to that around a Schwarzschild black hole, rendering a larger difference between ion and electron temperatures in the former case. With all the theoretical solutions in hand, finally we reproduce the observed luminosities (L) of two extreme cases-the under-fed AGNs and quasars (e.g. Sgr A') with L greater than or similar to 10(33) erg/s to ultra-luminous X-ray sources with L similar to 10(41) erg/s, at different combinations of mass accretion rate, ratio of specific heats, Shakura-Sunyaev viscosity parameter and Kerr parameter, and conclude that Sgr A' may be an intermediate spinning black hole.

A real-time algorithm for simulation of flexible objects and hyper-redundant manipulators

Flexible objects such as a rope or snake move in a way such that their axial length remains almost constant. To simulate the motion of such an object, one strategy is to discretize the object into large number of small rigid links connected by joints. However, the resulting discretised system is highly redundant and the joint rotations for a desired Cartesian motion of any point on the object cannot be solved uniquely. In this paper, we revisit an algorithm, based on the classical tractrix curve, to resolve the redundancy in such hyper-redundant systems. For a desired motion of the `head' of a link, the `tail' is moved along a tractrix, and recursively all links of the discretised objects are moved along different tractrix curves. The algorithm is illustrated by simulations of a moving snake, tying of knots with a rope and a solution of the inverse kinematics of a planar hyper-redundant manipulator. The simulations show that the tractrix based algorithm leads to a more `natural' motion since the motion is distributed uniformly along the entire object with the displacements diminishing from the `head' to the `tail'.

The virial expansion of a classical interacting system

We consider N particles interacting pairwise by an inverse square potential in one dimension (Calogero-Sutherland-Moser model). For a system placed in a harmonic trap, its classical partition function for the repulsive regime is recognised in the literature. We start by presenting a concise re-derivation of this result. The equation of state is then calculated both for the trapped and the homogeneous gas. Finally, the classical limit of Wu's distribution function for fractional exclusion statistics is obtained and we re-derive the classical virial expansion of the homogeneous gas using this distribution function.

Two-temperature accretion around rotating black holes: a description of the general advective flow paradigm in the presence of various cooling processes to explain low to high luminous sources

We investigate viscous two-temperature accretion disc flows around rotating black holes. We describe the global solution of accretion flows with a sub-Keplerian angular momentum profile, by solving the underlying conservation equations including explicit cooling processes self-consistently. Bremsstrahlung, synchrotron and inverse Comptonization of soft photons are considered as possible cooling mechanisms. We focus on the set of solutions for sub-Eddington, Eddington and super-Eddington mass accretion rates around Schwarzschild and Kerr black holes with a Kerr parameter of 0.998. It is found that the flow, during its infall from the Keplerian to sub-Kepleria transition region to the black hole event horizon, passes through various phases of advection: the general advective paradigm to the radiatively inefficient phase, and vice versa. Hence, the flow governs a much lower electron temperature similar to 10(8)-10(9.5) K, in the range of accretion rate in Eddington units 0.01 less than or similar to (M) over dot less than or similar to 100, compared to the hot protons of temperature similar to 10(10.2)-10(11.8) K. Therefore, the solution may potentially explain the hard X-rays and gamma-rays emitted from active galactic nuclei (AGNs) and X-ray binaries. We then compare the solutions for two different regimes of viscosity. We conclude that a weakly viscous flow is expected to be cooling dominated, particularly at the inner region of the disc, compared to its highly viscous counterpart, which is radiatively inefficient. With all the solutions in hand, we finally reproduce the observed luminosities of the underfed AGNs and quasars (e. g. Sgr A*) to ultraluminous X-ray sources (e. g. SS433), at different combinations of input parameters, such as the mass accretion rate and the ratio of specific heats. The set of solutions also predicts appropriately the luminosity observed in highly luminous AGNs and ultraluminous quasars (e. g. PKS 0743-67).

Convergence analysis of the Newton algorithm and a pseudo-time marching scheme for diffuse correlation tomography

We propose a self-regularized pseudo-time marching scheme to solve the ill-posed, nonlinear inverse problem associated with diffuse propagation of coherent light in a tissuelike object. In particular, in the context of diffuse correlation tomography (DCT), we consider the recovery of mechanical property distributions from partial and noisy boundary measurements of light intensity autocorrelation. We prove the existence of a minimizer for the Newton algorithm after establishing the existence of weak solutions for the forward equation of light amplitude autocorrelation and its Frechet derivative and adjoint. The asymptotic stability of the solution of the ordinary differential equation obtained through the introduction of the pseudo-time is also analyzed. We show that the asymptotic solution obtained through the pseudo-time marching converges to that optimal solution provided the Hessian of the forward equation is positive definite in the neighborhood of optimal solution. The superior noise tolerance and regularization-insensitive nature of pseudo-dynamic strategy are proved through numerical simulations in the context of both DCT and diffuse optical tomography. (C) 2010 Optical Society of America.

Benchmark scenarios for the NMSSM

We discuss constrained and semi--constrained versions of the next--to--minimal supersymmetric extension of the Standard Model (NMSSM) in which a singlet Higgs superfield is added to the two doublet superfields that are present in the minimal extension (MSSM). This leads to a richer Higgs and neutralino spectrum and allows for many interesting phenomena that are not present in the MSSM. In particular, light Higgs particles are still allowed by current constraints and could appear as decay products of the heavier Higgs states, rendering their search rather difficult at the LHC. We propose benchmark scenarios which address the new phenomenological features, consistent with present constraints from colliders and with the dark matter relic density, and with (semi--)universal soft terms at the GUT scale. We present the corresponding spectra for the Higgs particles, their couplings to gauge bosons and fermions and their most important decay branching ratios. A brief survey of the search strategies for these states at the LHC is given.

Oxide Particle Surface Chemistry and Ion Transport in "Soggy Sand" Electrolytes

The crucial role of oxide surface chemical composition on ion transport in "soggy sand" electrolytes is discussed in a systematic manner. A prototype soggy sand electrolytic system comprising aerosil silica functionalized with various hydrophilic and hydrophobic moieties dispersed in lithium perchlorate-ethylene glycol solution was used for the study. Detailed rheology studies show that the attractive particle network in the case of the composite with unmodified aerosil silica (with surface silanol groups) is most favorable for percolation in ionic conductivity, as well as rendering the composite with beneficial elastic mechanical properties: Though weaker in strength compared to the composite with unmodified aerosil particles, attractive particle networks are also observed in composites of aerosil particles with surfaces partially substituted with hydrophobic groups. The percolation in ionic conductivity is, however, dependent on the size of the hydrophobic moiety. No spanning attractive particle network was formed for aerosil particles with surfaces modified with stronger hydrophilic groups (than silanol), and as a result, no percolation in ionic conductivity was observed. The composite with hydrophilic particles was a sol, contrary to gels obtained in the case of unmodified aerosil, and partially substituted with hydrophobic groups.