Robust Reentry Guidance of a Reusable Launch Vehicle Using Model Predictive Static Programming

A robust suboptimal reentry guidance scheme is presented for a reusable launch vehicle using the recently developed, computationally efficient model predictive static programming. The formulation uses the nonlinear vehicle dynamics with a spherical and rotating Earth, hard constraints for desired terminal conditions, and an innovative cost function having several components with associated weighting factors that can account for path and control constraints in a soft constraint manner, thereby leading to smooth solutions of the guidance parameters. The proposed guidance essentially shapes the trajectory of the vehicle by computing the necessary angle of attack and bank angle that the vehicle should execute. The path constraints are the structural load constraint, thermal load constraint, bounds on the angle of attack, and bounds on the bank angle. In addition, the terminal constraints include the three-dimensional position and velocity vector components at the end of the reentry. Whereas the angle-of-attack command is generated directly, the bank angle command is generated by first generating the required heading angle history and then using it in a dynamic inversion loop considering the heading angle dynamics. Such a two-loop synthesis of bank angle leads to better management of the vehicle trajectory and avoids mathematical complexity as well. Moreover, all bank angle maneuvers have been confined to the middle of the trajectory and the vehicle ends the reentry segment with near-zero bank angle, which is quite desirable. It has also been demonstrated that the proposed guidance has sufficient robustness for state perturbations as well as parametric uncertainties in the model.

Nanodiamond-Based Thermal Fluids

Dispersions of nanodiamond (average size similar to 6 nm) within dielectric insulator mineral oil are reported for their enhanced thermal conductivity properties and potential applications in thermal management. Dynamic and kinematic viscosities-very important parameters in thermal management by nanofluids-are investigated. The dependence of the dynamic viscosity is well-described by the theoretical predictions of Einstein's model. The temperature dependence of the dynamic viscosity obeys an Arrhenius-like behavior, where the activation energy and the pre-exponential factor have an exponential dependence on the filler fraction of nanodiamonds. An enhancement in thermal conductivity up to 70% is reported for nanodiamond based thermal fluids. Additional electron microscopy, Raman spectroscopy and X-ray diffraction analysis support the experimental data and their interpretation.

The generalized Onsager model for a binary gas mixture

The Onsager model for the secondary flow field in a high-speed rotating cylinder is extended to incorporate the difference in mass of the two species in a binary gas mixture. The base flow is an isothermal solid-body rotation in which there is a balance between the radial pressure gradient and the centrifugal force density for each species. Explicit expressions for the radial variation of the pressure, mass/mole fractions, and from these the radial variation of the viscosity, thermal conductivity and diffusion coefficient, are derived, and these are used in the computation of the secondary flow. For the secondary flow, the mass, momentum and energy equations in axisymmetric coordinates are expanded in an asymptotic series in a parameter epsilon = (Delta m/m(av)), where Delta m is the difference in the molecular masses of the two species, and the average molecular mass m(av) is defined as m(av) = (rho(w1)m(1) + rho(w2)m(2))/rho(w), where rho(w1) and rho(w2) are the mass densities of the two species at the wall, and rho(w) = rho(w1) + rho(w2). The equation for the master potential and the boundary conditions are derived correct to O(epsilon(2)). The leading-order equation for the master potential contains a self-adjoint sixth-order operator in the radial direction, which is different from the generalized Onsager model (Pradhan & Kumaran, J. Fluid Mech., vol. 686, 2011, pp. 109-159), since the species mass difference is included in the computation of the density, viscosity and thermal conductivity in the base state. This is solved, subject to boundary conditions, to obtain the leading approximation for the secondary flow, followed by a solution of the diffusion equation for the leading correction to the species mole fractions. The O(epsilon) and O(epsilon(2)) equations contain inhomogeneous terms that depend on the lower-order solutions, and these are solved in a hierarchical manner to obtain the O(epsilon) and O(epsilon(2)) corrections to the master potential. A similar hierarchical procedure is used for the Carrier-Maslen model for the end-cap secondary flow. The results of the Onsager hierarchy, up to O(epsilon(2)), are compared with the results of direct simulation Monte Carlo simulations for a binary hard-sphere gas mixture for secondary flow due to a wall temperature gradient, inflow/outflow of gas along the axis, as well as mass and momentum sources in the flow. There is excellent agreement between the solutions for the secondary flow correct to O(epsilon(2)) and the simulations, to within 15 %, even at a Reynolds number as low as 100, and length/diameter ratio as low as 2, for a low stratification parameter A of 0.707, and when the secondary flow velocity is as high as 0.2 times the maximum base flow velocity, and the ratio 2 Delta m/(m(1) + m(2)) is as high as 0.5. Here, the Reynolds number Re = rho(w)Omega R-2/mu, the stratification parameter A = root m Omega R-2(2)/(2k(B)T), R and Omega are the cylinder radius and angular velocity, m is the molecular mass, rho(w) is the wall density, mu is the viscosity and T is the temperature. The leading-order solutions do capture the qualitative trends, but are not in quantitative agreement.

Effect of shrinkage induced flow on binary alloy dendrite growth: An equivalent undercooling model

This paper presents a theoretical model for studying the effects of shrinkage induced flow on the growth rate of binary alloy dendrites. An equivalent undercooling of the melt is defined in terms of ratio of the phase densities to represent the change in dendrite growth rate due to variation in solutal and thermal transport resulting from shrinkage induced flow. Subsequently, results for dendrite growth rate predicted by the equivalent undercooling model is compared with the corresponding predictions obtained using an enthalpy based numerical method for dendrite growth with shrinkage. The agreement is found to be good. Published by Elsevier Ltd.

The cold mode: a phenomenological model for the evolution of density perturbations in the intracluster medium

Cool cluster cores are in global thermal equilibrium but are locally thermally unstable. We study a non-linear phenomenological model for the evolution of density perturbations in the intracluster medium (ICM) due to local thermal instability and gravity. We have analysed and extended a model for the evolution of an overdense blob in the ICM. We find two regimes in which the overdense blobs can cool to thermally stable low temperatures. One for large t(cool)/t(ff) (t(cool) is the cooling time and t(ff) is the free-fall time), where a large initial overdensity is required for thermal runaway to occur; this is the regime which was previously analysed in detail. We discover a second regime for t(cool)/t(ff) less than or similar to 1 (in agreement with Cartesian simulations of local thermal instability in an external gravitational field), where runaway cooling happens for arbitrarily small amplitudes. Numerical simulations have shown that cold gas condenses out more easily in a spherical geometry. We extend the analysis to include geometrical compression in weakly stratified atmospheres such as the ICM. With a single parameter, analogous to the mixing length, we are able to reproduce the results from numerical simulations; namely, small density perturbations lead to the condensation of extended cold filaments only if t(cool)/t(ff) less than or similar to 10.

Electron transfer statistics and thermal fluctuations in molecular junctions

We derive analytical expressions for probability distribution function (PDF) for electron transport in a simple model of quantum junction in presence of thermal fluctuations. Our approach is based on the large deviation theory combined with the generating function method. For large number of electrons transferred, the PDF is found to decay exponentially in the tails with different rates due to applied bias. This asymmetry in the PDF is related to the fluctuation theorem. Statistics of fluctuations are analyzed in terms of the Fano factor. Thermal fluctuations play a quantitative role in determining the statistics of electron transfer; they tend to suppress the average current while enhancing the fluctuations in particle transfer. This gives rise to both bunching and antibunching phenomena as determined by the Fano factor. The thermal fluctuations and shot noise compete with each other and determine the net (effective) statistics of particle transfer. Exact analytical expression is obtained for delay time distribution. The optimal values of the delay time between successive electron transfers can be lowered below the corresponding shot noise values by tuning the thermal effects. (C) 2015 AIP Publishing LLC.

Phase diagram of the half-filled ionic Hubbard model

We study the phase diagram of the ionic Hubbard model (IHM) at half filling on a Bethe lattice of infinite connectivity using dynamical mean-field theory (DMFT), with two impurity solvers, namely, iterated perturbation theory (IPT) and continuous time quantum Monte Carlo (CTQMC). The physics of the IHM is governed by the competition between the staggered ionic potential Delta and the on-site Hubbard U. We find that for a finite Delta and at zero temperature, long-range antiferromagnetic (AFM) order sets in beyond a threshold U = U-AF via a first-order phase transition. For U smaller than U-AF the system is a correlated band insulator. Both methods show a clear evidence for a quantum transition to a half-metal (HM) phase just after the AFM order is turned on, followed by the formation of an AFM insulator on further increasing U. We show that the results obtained within both methods have good qualitative and quantitative consistency in the intermediate-to-strong-coupling regime at zero temperature as well as at finite temperature. On increasing the temperature, the AFM order is lost via a first-order phase transition at a transition temperature T-AF(U,Delta) or, equivalently, on decreasing U below U-AF(T,Delta)], within both methods, for weak to intermediate values of U/t. In the strongly correlated regime, where the effective low-energy Hamiltonian is the Heisenberg model, IPT is unable to capture the thermal (Neel) transition from the AFM phase to the paramagnetic phase, but the CTQMC does. At a finite temperature T, DMFT + CTQMC shows a second phase transition (not seen within DMFT + IPT) on increasing U beyond U-AF. At U-N > U-AF, when the Neel temperature T-N for the effective Heisenberg model becomes lower than T, the AFM order is lost via a second-order transition. For U >> Delta, T-N similar to t(2)/U(1 - x(2)), where x = 2 Delta/U and thus T-N increases with increase in Delta/U. In the three-dimensional parameter space of (U/t, T/t, and Delta/t), as T increases, the surface of first-order transition at U-AF(T,Delta) and that of the second-order transition at U-N(T,Delta) approach each other, shrinking the range over which the AFM order is stable. There is a line of tricritical points that separates the surfaces of first- and second-order phase transitions.

Active dendrites mediate stratified gamma-range coincidence detection in hippocampal model neurons

Hippocampal pyramidal neurons exhibit gamma-phase preference in their spikes, selectively route inputs through gamma frequency multiplexing and are considered part of gamma-bound cell assemblies. How do these neurons exhibit gamma-frequency coincidence detection capabilities, a feature that is essential for the expression of these physiological observations, despite their slow membrane time constant? In this conductance-based modelling study, we developed quantitative metrics for the temporal window of integration/coincidence detection based on the spike-triggered average (STA) of the neuronal compartment. We employed these metrics in conjunction with quantitative measures for spike initiation dynamics to assess the emergence and dependence of coincidence detection and STA spectral selectivity on various ion channel combinations. We found that the presence of resonating conductances (hyperpolarization-activated cyclic nucleotide-gated or T-type calcium), either independently or synergistically when expressed together, led to the emergence of spectral selectivity in the spike initiation dynamics and a significant reduction in the coincidence detection window (CDW). The presence of A-type potassium channels, along with resonating conductances, reduced the STA characteristic frequency and broadened the CDW, but persistent sodium channels sharpened the CDW by strengthening the spectral selectivity in the STA. Finally, in a morphologically precise model endowed with experimentally constrained channel gradients, we found that somatodendritic compartments expressed functional maps of strong theta-frequency selectivity in spike initiation dynamics and gamma-range CDW. Our results reveal the heavy expression of resonating and spike-generating conductances as the mechanism underlying the robust emergence of stratified gamma-range coincidence detection in the dendrites of hippocampal and cortical pyramidal neurons.

A novel interferometric technique to estimate thermal diffusivity of `optically transparent solid' using `Isothermal surface velocimetry'

One-dimensional transient heat flow is interpreted as a procession of `macro-scale translatory motion of indexed isothermal surfaces'. A new analytical model is proposed by introducing velocity of isothermal surface in Fourier heat diffusion equation. The velocity dependent function is extracted by revisiting `the concept of thermal layer of heat conduction in solid' and `exact solution' to estimate thermal diffusivity. The experimental approach involves establishment of 1 D unsteady heat flow inside the sample through Step-temperature excitation. A novel self-reference interferometer is utilized to separate a `unique isothermal surface' in time-varying temperature field. The translatory motion of the said isothermal surface is recorded using digital camera to estimate its velocity. From the knowledge of thermo-optic coefficient, temperature of the said isothermal surface is predicted. The performance of proposed method is evaluated for Quartz sample and compared with literature.

Meridional gradients in aerosol vertical distribution over Indian Mainland: Observations and model simulations

Multi-year observations from the network of ground-based observatories (ARFINET), established under the project `Aerosol Radiative Forcing over India' (ARFI) of Indian Space Research Organization and space-borne lidar `Cloud Aerosol Lidar with Orthogonal Polarization' (CALIOP) along with simulations from the chemical transport model `Goddard Chemistry Aerosol Radiation and Transport' (GOCART), are used to characterize the vertical distribution of atmospheric aerosols over the Indian landmass and its spatial structure. While the vertical distribution of aerosol extinction showed higher values close to the surface followed by a gradual decrease at increasing altitudes, a strong meridional increase is observed in the vertical spread of aerosols across the Indian region in all seasons. It emerges that the strong thermal convections cause deepening of the atmospheric boundary layer, which although reduces the aerosol concentration at lower altitudes, enhances the concentration at higher elevations by pumping up more aerosols from below and also helping the lofted particles to reach higher levels in the atmosphere. Aerosol depolarization ratios derived from CALIPSO as well as the GOCART simulations indicate the dominance of mineral dust aerosols during spring and summer and anthropogenic aerosols in winter. During summer monsoon, though heavy rainfall associated with the Indian monsoon removes large amounts of aerosols, the prevailing southwesterly winds advect more marine aerosols over to landmass (from the adjoining oceans) leading to increase in aerosol loading at lower altitudes than in spring. During spring and summer months, aerosol loading is found to be significant, even at altitudes as high as 4 km, and this is proposed to have significant impacts on the regional climate systems such as Indian monsoon. (C) 2015 Elsevier Ltd. All rights reserved.

Rayleigh-Taylor instability in a thermocline based thermal storage tank

Numerical simulations are performed to study the stability characteristics of a molten salt thermocline storage unit. Perturbations are introduced into a stable flow field in such a way as to make the top-fluid heavier than the fluid at the bottom, thereby causing a possible instability in the system. The evolution pattern of the various disturbances are examined in detail. Disturbances applied for short duration get decayed before they could reach the thermocline, whereas medium and long duration disturbances evolve into a ``falling spike'' or ``stalactite-like'' structure and destabilize the thermocline. Rayleigh Taylor instability is observed inside the storage tank. The effect of the duration, velocity and temperature of the disturbance on thermocline thickness and penetration length are studied. A quadratic time dependence of penetration length was observed. New perspectives on thermocline breakdown phenomena are obtained from the numerical flow field. (C) 2015 Elsevier Masson SAS. All rights reserved.

Fabrication and Nonlinear Thermomechanical Analysis of SU8 Thermal Actuator

SU8-based micromechanical structures are widely used as thermal actuators in the development of compliant micromanipulation tools. This paper reports the design, nonlinear thermomechanical analysis, fabrication, and thermal actuation of SU8 actuators. The thermomechanical analysis of the actuator incorporates nonlinear temperature-dependent properties of SU8 polymer to accurately model its thermal response during actuation. The designed SU8 thermal actuators are fabricated using surface micromachining techniques and the electrical interconnects are made to them using flip-chip bonding. The issues due to thermal stress during fabrication are discussed and a novel strategy is proposed to release the thermal stress in the fabricated actuators. Subsequent characterization of the actuator using an optical profilometer reveals excellent thermal response, good repeatability, and low hysteresis. The average deflection is similar to 8.5 mu m for an actuation current of similar to 5 mA. The experimentally obtained deflection profile and the tip deflection at different currents are both shown to be in good agreement with the predictions of the nonlinear thermomechanical model. This underscores the need to consider nonlinearities when modeling the response of SU8 thermal actuators. 2015-0087]