Embracing Integrability in Stellar and Planetary Dynamics

Hamiltonian systems in stellar and planetary dynamics are typically near integrable. For example, Solar System planets are almost in two-body orbits, and in simulations of the Galaxy, the orbits of stars seem regular. For such systems, sophisticated numerical methods can be developed through integrable approximations. Following this theme, we discuss three distinct problems. We start by considering numerical integration techniques for planetary systems. Perturbation methods (that utilize the integrability of the two-body motion) are preferred over conventional "blind" integration schemes. We introduce perturbation methods formulated with Cartesian variables. In our numerical comparisons, these are superior to their conventional counterparts, but, by definition, lack the energy-preserving properties of symplectic integrators. However, they are exceptionally well suited for relatively short-term integrations in which moderately high positional accuracy is required. The next exercise falls into the category of stability questions in solar systems. Traditionally, the interest has been on the orbital stability of planets, which have been quantified, e.g., by Liapunov exponents. We offer a complementary aspect by considering the protective effect that massive gas giants, like Jupiter, can offer to Earth-like planets inside the habitable zone of a planetary system. Our method produces a single quantity, called the escape rate, which characterizes the system of giant planets. We obtain some interesting results by computing escape rates for the Solar System. Galaxy modelling is our third and final topic. Because of the sheer number of stars (about 10^11 in Milky Way) galaxies are often modelled as smooth potentials hosting distributions of stars. Unfortunately, only a handful of suitable potentials are integrable (harmonic oscillator, isochrone and Stäckel potential). This severely limits the possibilities of finding an integrable approximation for an observed galaxy. A solution to this problem is torus construction; a method for numerically creating a foliation of invariant phase-space tori corresponding to a given target Hamiltonian. Canonically, the invariant tori are constructed by deforming the tori of some existing integrable toy Hamiltonian. Our contribution is to demonstrate how this can be accomplished by using a Stäckel toy Hamiltonian in ellipsoidal coordinates.

Kinematics of pantograph masts

This paper deals with the kinematics of pantograph masts. Pantograph masts have widespread use in space application as deployable structures. They are over constrained mechanisms with degree-of-freedom, evaluated by the Grübler–Kutzback formula, as less than one. In this paper, a numerical algorithm is used to evaluate the degree-of-freedom of pantograph masts by obtaining the null space of a constraint Jacobian matrix. In the process redundant joints in the masts are obtained. A method based on symbolic computation, to obtain the closed-form kinematics equations of triangular and box shaped pantograph masts, is presented. In the process, the various configurations such masts can attain during deployment, are obtained. The closed-form solution also helps in identifying the redundant joints in the masts. The symbolic computations involving the Jacobian matrix also leads to a method to evaluate the global degree-of-freedom for these masts.

MPD arcjet cathode surface current density

The radial current density distribution on the cathode longitudinal surface of magnetoplasmadynamic arcjets for axisymmetric geometries has been obtained by simultaneous solution of the electromagnetic equations for a given uniform gas dynamic field. The problem formulation permits a parametric study of the effects of the Hall parameter and the magnetic Reynolds number. The solution for the current density distribution displays current concentrations at two locations, that is, at the upstream and downstream ends of the cathode. This result is in conformity with known experimental data. The parameters responsible for these current concentrations are identified. It is shown that the effect of the magnetic Reynolds number on the current density distribution is different depending on whether or not the Hall effect is included. This result is also found to be consistent with experimental data.

Organometallic derivatives of diphosphazanes. 2. Seven-coordinated Group 6 metal tricarbonyl complexes of diphosphazane ligands. X-ray crystal structure of [WI2(CO)3{P(OPh)2}2NPh]

The reactions of the complexes [MI2(CO)3-(NCMe)2] (M = Mo, W) with the diphosphazane ligands RN{P(OPh)2}2 (R = Me, Ph) in CH2Cl2 at room temperature afford new seven-coordinated complexes of the type [MI2(CO)3{P(OPh)2}2NR]. The molybdenum complexes are sensitive to air oxidation even in the solid state, whereas the tungsten complexes are more stable in the solid state and in solution. The structure of the tungsten complex [WI2(CO)3{P(OPh)2}2NPh] has been determined by single-crystal X-ray diffraction. It crystallizes in the orthorhombic system with the space group Pna 2(1), a = 19.372 (2) angstrom, b = 11.511 (1) angstrom, c = 15.581 (1) angstrom, and Z = 4. Full-matrix least-squares refinement with 3548 reflections (I > 2.5-sigma-(I)) led to final R and R(w) values of 0.036 and 0.034, respectively. The complex adopts a slightly distorted pentagonal-bypyramidal geometry rarely observed for such a type of complexes; two phosphorus atoms of the diphosphazane ligand, two iodine atoms, and a carbonyl group occupy the equatorial plane, and the other two carbonyl groups, the apical positions.

Perturbation of critical solution temperatures by impurity doping

An expression is developed for the variation of the critical solution temperature of a binary liquid system when a third component (dopant) is added, using an extension of the regular solution theory. The model can be used for UCST, LCST and for closed loop systems and has the correct features in the limiting cases.

Reactivity Of Pph3 Toward Ruiiruiiicl(O2car)4 - Syntheses, Molecular-Structures, And Spectroscopic And Electrochemical Properties Of Ruiiruiii(Oh2)Cl(Mecn)(O2car)4(Pph3)2 And Ruii2(Oh2)(Mecn)2(O2car)4(Pph3)2

By the reaction of Ru2Cl(O2CAr)4 (1) and PPh3 in MeCN-H2O the diruthenium(II,III) and diruthenium(II) compounds of the type Ru2(OH2)Cl(MeCN)(O2CAr)4(PPh3)2 (2) and Ru2(OH2)(MeCN)2(O2CAr)4(PPh3)2 (3) were prepared and characterized by analytical, spectral, and electrochemical data (Ar is an aryl group, C6H4-p-X; X = H, OMe, Me, Cl, NO2). The molecular structure of Ru2(OH2)Cl(MeCN)(O2CC6H4-p-OMe)4(PPh3)2 was determined by X-ray crystallography. Crystal data are as follows: triclinic, P1BAR, a = 13.538 (5) angstrom, b = 15.650 (4) angstrom, c = 18.287 (7) angstrom, alpha = 101.39 (3)-degrees, beta = 105.99 (4)-degrees, gamma = 97.94 (3)-degrees, V = 3574 angstrom 3, Z = 2. The molecule is asymmetric, and the two ruthenium centers are clearly distinguishable. The Ru(III)-Ru(II), Ru(III)-(mu-OH2), and Ru(II)-(mu-OH2) distances and the Ru-(mu-OH2)-Ru angle in [{Ru(III)Cl(eta-1-O2CC6H4-p-OMe)(PPh3)}(mu-OH2)(mu-O2CC6H4-p-OMe)2{Ru(II)(MeCN)(eta-1-O2CC6H4-p-OMe)(PPh3)}] are 3.604 (1), 2.127 (8), and 2.141 (10) angstrom and 115.2 (5)-degrees, respectively. The compounds are paramagnetic and exhibit axial EPR spectra in the polycrystalline form. An intervalence transfer (IT) transition is observed in the range 900-960 nm in chloroform in these class II type trapped mixed-valence species 2. Compound 2 displays metal-centered one-electron reduction and oxidation processes near -0.4 and +0.6 V (vs SCE), respectively in CH2Cl2-TBAP. Compound 2 is unstable in solution phase and disproportionates to (mu-aquo)diruthenium(II) and (mu-oxo)diruthenium(III) complexes. The mechanistic aspects of the core conversion are discussed. The molecular structure of a diruthenium(II) compound, Ru2(OH2)(MeCN)2(O2CC6H4-p-NO2)4(PPh3)2.1.5CH2Cl2, was obtained by X-ray crystallography. The compound crystallizes in the space group P2(1)/c with a = 23.472 (6) angstrom, b = 14.303 (3) angstrom, c = 23.256 (7) angstrom, beta = 101.69 (2)-degrees, V = 7645 angstrom 3, and Z = 4. The Ru(II)-Ru(II) and two Ru(II)-(mu-OH2) distances and the Ru(II)-(mu-OH2)-Ru(II) angle in [{(PPh3)-(MeCN)(eta-1-O2CC6H4-p-NO2)Ru}2(mu-OH2)(mu-O2CC6H4-p-NO2)2] are 3.712 (1), 2.173 (9), and 2.162 (9) angstrom and 117.8 (4)-degrees, respectively. In both diruthenium(II,III) and diruthenium(II) compounds, each metal center has three facial ligands of varying pi-acidity and the aquo bridges are strongly hydrogen bonded with the eta-1-carboxylato facial ligands. The diruthenium(II) compounds are diamagnetic and exhibit characteristic H-1 NMR spectra in CDCl3. These compounds display two metal-centered one-electron oxidations near +0.3 and +1.0 V (vs SCE) in CH2Cl2-TBAP. The overall reaction between 1 and PPh3 in MeCN-H2O through the intermediacy of 2 is of the disproportionation type. The significant role of facial as well as bridging ligands in stabilizing the core structures is observed from electrochemical studies.

Solution concepts in continuous-kernel multicriteria games

This paper considers nonzero-sum multicriteria games with continuous kernels. Solution concepts based on the notions of Pareto optimality, equilibrium, and security are extended to these games. Separate necessary and sufficient conditions and existence results are presented for equilibrium, Pareto-optimal response, and Pareto-optimal security strategies of the players.

Solution of a Boundary-Value Problem associated with Diffraction of Water Waves by a Nearly Vertical Barrier

An exact solution is derived for a boundary-value problem for Laplace's equation which is a generalization of the one occurring in the course of solution of the problem of diffraction of surface water waves by a nearly vertical submerged barrier. The method of solution involves the use of complex function theory, the Schwarz reflection principle, and reduction to a system of two uncoupled Riemann-Hilbert problems. Known results, representing the reflection and transmission coefficients of the water wave problem involving a nearly vertical barrier, are derived in terms of the shape function.

Steady incompressible flow of cohesionless granular materials through a wedge-shaped hopper: Frictional 14kinetic solution to the smooth wall, radial gravity problem

Hybrid frictional-kinetic equations are used to predict the velocity, grain temperature, and stress fields in hoppers. A suitable choice of dimensionless variables permits the pseudo-thermal energy balance to be decoupled from the momentum balance. These balances contain a small parameter, which is analogous to a reciprocal Reynolds number. Hence an approximate semi-analytical solution is constructed using perturbation methods. The energy balance is solved using the method of matched asymptotic expansions. The effect of heat conduction is confined to a very thin boundary layer near the exit, where it causes a marginal change in the temperature. Outside this layer, the temperature T increases rapidly as the radial coordinate r decreases. In particular, the conduction-free energy balance yields an asymptotic solution, valid for small values of r, of the form T proportional r-4. There is a corresponding increase in the kinetic stresses, which attain their maximum values at the hopper exit. The momentum balance is solved by a regular perturbation method. The contribution of the kinetic stresses is important only in a small region near the exit, where the frictional stresses tend to zero. Therefore, the discharge rate is only about 2.3% lower than the frictional value, for typical parameter values. As in the frictional case, the discharge rate for deep hoppers is found to be independent of the head of material.

Rapid synthesis of room temperature ferromagnetic Ag-doped LaMnO3 perovskite phases by the solution combustion method

We report the rapid solution combustion synthesis and characterization of Ag-substituted LaMnO3 phases at relatively low temperature using oxalyl dihydrazide, as fuel. Structural parameters were refined by the Rietveld method using powder X-ray diffraction data. While the parent LaMnO3 crystallizes in the orthorhombic structure, the Ag-substituted compounds crystallize in the rhombohedral symmetry. On increasing Ag-content, unit cell volume and Mn-O-Mn bond angle decreases. The Fourier transform infra red spectrum shows two absorption bands corresponding to Mn-O stretching vibration (v(s) mode) and Mn-O-Mn deformation vibration (v(b) mode) around 600 cm(-1) and 400 cm(-1) for the compositions x = 0.0, 0.05 and 0.10, respectively. Electrical resistivity measurements reveal that composition-controlled metal to insulator transition, with the maximum metal to insulator being 280 K for the composition La0.75Ag0.25MnO3. Increase in magnetic moment was observed with increase in Ag-content. The maximum magnetic moment of 35 emu/g was observed for the composition La0.80Ag0.20MnO3. (C) 2010 Elsevier Ltd. All rights reserved.

Exact solution for the unsteady motion of a viscous-fluid in a porous annulus

An exact solution to the problem of time-dependent motion of a viscous fluid in an annulus with porous walls is obtained under the assumption that the rate of suction at one wall is equal to the rate of injection at the other. Finite Hankel transform is used to obtain a closed-form solution for the axial velocity. The average axial velocity profiles are depicted graphically.

Dilute solution properties of methyl methacrylate—acrylonitrile copolymer (MA2)

This article deals with studies of the dilute solution properties of methyl methacrylate-acrylonitrile (MMA-AN) copolymer of 0.415 mole fraction (mf) of acrylonitrile composition. Mark—Houwink parameters for this copolymer have been evaluated in acetonitrile (MeCN), 2-butanone (MEK), dimethyl formamide (DMF), and γ-butyrolactone (γ-BL). The Mark-Houwink exponent a in all four solvents at all temperatures is larger than the corresponding values of the parent homopolymers. The solvent power is in the order of DMF < γ–BL < MEK < MeCN; [η] decreases with an increase in temperature, which is behavior characteristic of polymers in good solvent. The unperturbed dimensions (K0) values, obtained by the Stockmayer–Fixman method, are lower than those for the parent homopolymers and depend on solvent as well as temperature. The solute—solvent interaction parameter X1 values are close to 0.5; X1 is independent of temperature. The excess interaction parameter XABvalues are negative. The results for this copolymer system in regard to low second virial coefficient A2, large X1, and high a values suggest that the large extension of these copolymer chains is due to the unusual short-range interactions.

Comparative study of different formulations and solution methods for two phase flow in saturated porous media

Six models (Simulators) are formulated and developed with all possible combinations of pressure and saturation of the phases as primary variables. A comparative study between six simulators with two numerical methods, conventional simultaneous and modified sequential methods are carried out. The results of the numerical models are compared with the laboratory experimental results to study the accuracy of the model especially in heterogeneous porous media. From the study it is observed that the simulator using pressure and saturation of the wetting fluid (PW, SW formulation) is the best among the models tested. Many simulators with nonwetting phase as one of the primary variables did not converge when used along with simultaneous method. Based on simulator 1 (PW, SW formulation), a comparison of different solution methods such as simultaneous method, modified sequential and adaptive solution modified sequential method are carried out on 4 test problems including heterogeneous and randomly heterogeneous problems. It is found that the modified sequential and adaptive solution modified sequential methods could save the memory by half and as also the CPU time required by these methods is very less when compared with that using simultaneous method. It is also found that the simulator with PNW and PW as the primary variable which had problem of convergence using the simultaneous method, converged using both the modified sequential method and also using adaptive solution modified sequential method. The present study indicates that pressure and saturation formulation along with adaptive solution modified sequential method is the best among the different simulators and methods tested.

Precursor Dependent Microstructure Evolution and Nonstoichiometry in Nanostructured Cerium Oxide Coatings Using the Solution Precursor Plasma Spray Technique

A solution precursor plasma spray (SPPS) technique has been used for direct deposition of cerium oxide nanoparticles (CNPs) from various cerium salt solutions as precursors. Solution precursors were injected into the hot zone of a plasma plume to deposit CNP coatings. A numerical study of the droplet injection model has been employed for microstructure development during SPPS. The decomposition of each precursor to cerium oxide was analyzed by thermogravimetric-differential thermal analysis and validated by thermodynamic calculations. The presence of the cerium oxide phase in the coatings was confirmed by X-ray diffraction studies. Transmission electron microscopy studies confirmed nanocrystalline (grain size <14 nm) characteristic of the coatings. X-ray photoelectron spectroscopy studies indicated the presence of a high concentration of Ce3+ (up to 0.32) in the coating prepared by SPPS. The processing and microstructure evolution of cerium oxide coatings with high nonstoichiometry are reported.