A scaled unscented transformation based directed Gaussian sum filter for nonlinear dynamic system identification

Impoverishment of particles, i.e. the discretely simulated sample paths of the process dynamics, poses a major obstacle in employing the particle filters for large dimensional nonlinear system identification. A known route of alleviating this impoverishment, i.e. of using an exponentially increasing ensemble size vis-a-vis the system dimension, remains computationally infeasible in most cases of practical importance. In this work, we explore the possibility of unscented transformation on Gaussian random variables, as incorporated within a scaled Gaussian sum stochastic filter, as a means of applying the nonlinear stochastic filtering theory to higher dimensional structural system identification problems. As an additional strategy to reconcile the evolving process dynamics with the observation history, the proposed filtering scheme also modifies the process model via the incorporation of gain-weighted innovation terms. The reported numerical work on the identification of structural dynamic models of dimension up to 100 is indicative of the potential of the proposed filter in realizing the stated aim of successfully treating relatively larger dimensional filtering problems. (C) 2013 Elsevier Ltd. All rights reserved.

Studies on rheocasting using cooling slope

In the present work, a cooling channel is employed to produce semi-solid A356 alloy slurry. To understand the transport process involved, a 3D non-isothermal, multiphase volume averaging model has been developed for simulation of the semi-solid slurry generation process in the cooling channel. For simulation purpose, the three phases considered are the parent melt, the nearly spherical grains and air as separated but highly coupled interpenetrating continua. The conservation equations of mass, momentum, energy and species have been solved for each phase and the thermal and mechanical interactions (drag force) among the phases have been considered using appropriate model. The superheated liquid alloy is poured at the top of the cooling slope/channel, where specified velocity inlet boundary condition is used in the model, and allowed to flow along gravity through the channel. The melt loses its superheat and becomes semisolid up to the end of cooling channel due to the evolving -Al grains with decreasing temperature. The air phase forms a definable air/liquid melt interface, i.e. free surface, due its low density. The results obtained from the present model includes volume fractions of three different phases considered, grain evolution, grain growth rate, size and distribution of solid grains. The effect of key process variables such as pouring temperature, slope angle of the cooling channel and cooling channel wall temperature on temperature distribution, velocity distribution, grain formation and volume fraction of different phases are also studied. The results obtained from the simulations are validated by microstructure study using SEM and quantitative image analysis of the semi-solid slurry microstructure obtained from the experimental set-up.

Free turbulent shear layer in a point vortex gas as a problem in nonequilibrium statistical mechanics

This paper attempts to unravel any relations that may exist between turbulent shear flows and statistical mechanics through a detailed numerical investigation in the simplest case where both can be well defined. The flow considered for the purpose is the two-dimensional (2D) temporal free shear layer with a velocity difference Delta U across it, statistically homogeneous in the streamwise direction (x) and evolving from a plane vortex sheet in the direction normal to it (y) in a periodic-in-x domain L x +/-infinity. Extensive computer simulations of the flow are carried out through appropriate initial-value problems for a ``vortex gas'' comprising N point vortices of the same strength (gamma = L Delta U/N) and sign. Such a vortex gas is known to provide weak solutions of the Euler equation. More than ten different initial-condition classes are investigated using simulations involving up to 32 000 vortices, with ensemble averages evaluated over up to 10(3) realizations and integration over 10(4)L/Delta U. The temporal evolution of such a system is found to exhibit three distinct regimes. In Regime I the evolution is strongly influenced by the initial condition, sometimes lasting a significant fraction of L/Delta U. Regime III is a long-time domain-dependent evolution towards a statistically stationary state, via ``violent'' and ``slow'' relaxations P.-H. Chavanis, Physica A 391, 3657 (2012)], over flow time scales of order 10(2) and 10(4)L/Delta U, respectively (for N = 400). The final state involves a single structure that stochastically samples the domain, possibly constituting a ``relative equilibrium.'' The vortex distribution within the structure follows a nonisotropic truncated form of the Lundgren-Pointin (L-P) equilibrium distribution (with negatively high temperatures; L-P parameter lambda close to -1). The central finding is that, in the intermediate Regime II, the spreading rate of the layer is universal over the wide range of cases considered here. The value (in terms of momentum thickness) is 0.0166 +/- 0.0002 times Delta U. Regime II, extensively studied in the turbulent shear flow literature as a self-similar ``equilibrium'' state, is, however, a part of the rapid nonequilibrium evolution of the vortex-gas system, which we term ``explosive'' as it lasts less than one L/Delta U. Regime II also exhibits significant values of N-independent two-vortex correlations, indicating that current kinetic theories that neglect correlations or consider them as O(1/N) cannot describe this regime. The evolution of the layer thickness in present simulations in Regimes I and II agree with the experimental observations of spatially evolving (3D Navier-Stokes) shear layers. Further, the vorticity-stream-function relations in Regime III are close to those computed in 2D Navier-Stokes temporal shear layers J. Sommeria, C. Staquet, and R. Robert, J. Fluid Mech. 233, 661 (1991)]. These findings suggest the dominance of what may be called the Kelvin-Biot-Savart mechanism in determining the growth of the free shear layer through large-scale momentum and vorticity dispersal.

Optimal sequential wireless relay placement on a random lattice path

Our work is motivated by impromptu (or ``as-you-go'') deployment of wireless relay nodes along a path, a need that arises in many situations. In this paper, the path is modeled as starting at the origin (where there is the data sink, e.g., the control center), and evolving randomly over a lattice in the positive quadrant. A person walks along the path deploying relay nodes as he goes. At each step, the path can, randomly, either continue in the same direction or take a turn, or come to an end, at which point a data source (e.g., a sensor) has to be placed, that will send packets to the data sink. A decision has to be made at each step whether or not to place a wireless relay node. Assuming that the packet generation rate by the source is very low, and simple link-by-link scheduling, we consider the problem of sequential relay placement so as to minimize the expectation of an end-to-end cost metric (a linear combination of the sum of convex hop costs and the number of relays placed). This impromptu relay placement problem is formulated as a total cost Markov decision process. First, we derive the optimal policy in terms of an optimal placement set and show that this set is characterized by a boundary (with respect to the position of the last placed relay) beyond which it is optimal to place the next relay. Next, based on a simpler one-step-look-ahead characterization of the optimal policy, we propose an algorithm which is proved to converge to the optimal placement set in a finite number of steps and which is faster than value iteration. We show by simulations that the distance threshold based heuristic, usually assumed in the literature, is close to the optimal, provided that the threshold distance is carefully chosen. (C) 2014 Elsevier B.V. All rights reserved.

Phase Field Simulation of Equiaxed Microstructure Formation during Semi-solid Processing of A380 Al Alloy

A phase field modelling approach is implemented in the present study towards simulation of microstructure evolution during cooling slope semi solid slurry generation process of A380 Aluminium alloy. First, experiments are performed to evaluate the number of seeds required within the simulation domain to simulate near spherical microstructure formation, occurs during cooling slope processing of the melt. Subsequently, microstructure evolution is studied employing a phase field method. Simulations are performed to understand the effect of cooling rate on the slurry microstructure. Encouraging results are obtained from the simulation studies which are validated by experimental observations. The results obtained from mesoscopic phase field simulations are grain size, grain density, degree of sphericity of the evolving primary Al phase and the amount of solid fraction present within the slurry at different time frames. Effect of grain refinement also has been studied with an aim of improving the slurry microstructure further. Insight into the process has been obtained from the numerical findings, which are found to be useful for process control.

An unusual demixing behavior in PS-PVME blends in the presence of nanoparticles

The effect of silver nanoparticles (sNP) on the demixing and the evolution of morphology in off-critical blends of 90/10 and 10/90 (wt/wt) PS/PVME polystyrene/poly(vinyl methyl ether)] was probed here using shear rheology and optical microscopy. The faster component (PVME) has a higher molecular weight (80 kDa) than the slower component (PS, 35 kDa), which makes this system quite interesting to study with respect to the evolving morphology, as the blends transit through the binodal and the spinodal envelopes. An unusual demixing behavior was observed in both PVME rich and PS rich blends. Temperature modulated differential scanning calorimetry measurements showed that the T-g value for the blends with sNP was slightly lower than that of the neat blends. A decreased volume of cooperativity at T-g suggests confined segmental dynamics in the presence of sNP. Although, the addition of sNP had no influence on the thermodynamic demixing temperature, it significantly altered the elasticity of the minor component during the transition of the blend from the homogeneous to the heterogeneous state. This is manifested from energetically driven localization of the sNP in the PVME phase during demixing. As a direct consequence of this, the formation of the microstructures upon demixing was observed to be delayed in the presence of sNP. Interestingly, in the intermediate quench depth, the higher viscoelastic phase evolved as an interconnected network, which subsequently coarsened into discrete droplets in the late stages for the 90/10 PS/PVME blends. Similar observations were made for 10/90 PS/PVME blends where threads of PVME appeared at deeper quench depths in the presence of sNP. The interconnected network formation of the minor phase (here PVME), which is also the faster component in the blend, was different from the usual demixing behavior.

Formation Flying with Nonlinear Partial Integrated Guidance and Control

Using the dynamic inversion philosophy, a nonlinear partial integrated guidance and control approach is presented in this paper for formation flying. It is based on the evolving philosophy of integrated guidance and control. However, it also retains the advantages of the conventional guidance then control philosophy by retaining the timescale separation between translational and rotational dynamics explicitly. Simulation studies demonstrate that the proposed technique is effective in bringing the vehicles into formation quickly and maintaining the formation.

The Evolution of Animal Welfare and the 3Rs in Brazil, China, and India

Increasingly, scientific collaborations and contracts cross country borders. The need for assurance that the quality of animal welfare and the caliber of animal research conducted are equivalent among research partners around the globe is of concern to the scientific and laboratory animal medicine communities, the general public, and other key stakeholders. Therefore, global harmonization of animal care and use standards and practices, with the welfare of the animals as a cornerstone, is essential. In the evolving global landscape of enhanced attention to animal welfare, a widely accepted path to achieving this goal is the successful integration of the 3Rs in animal care and use programs. Currently, awareness of the 3Rs, their implementation, and the resulting animal care and use standards and practices vary across countries. This variability has direct effects on the animals used in research and potentially the data generated and may also have secondary effects on the country's ability to be viewed as a global research partner. Here we review the status of implementation of the 3Rs worldwide and focus on 3 countries-Brazil, China and India-with increasing economic influence and an increasing footprint in the biomedical research enterprise.

Thermalization of Green functions and quasinormal modes

We develop a new method to study the thermalization of time dependent retarded Green function in conformal field theories holographically dual to thin shell AdS Vaidya space times. The method relies on using the information of all time derivatives of the Green function at the shell and then evolving it for later times. The time derivatives of the Green function at the shell is given in terms of a recursion formula. Using this method we obtain analytic results for short time thermalization of the Green function. We show that the late time behaviour of the Green function is determined by the first quasinormal mode. We then implement the method numerically. As applications of this method we study the thermalization of the retarded time dependent Green function corresponding to a minimally coupled scalar in the AdS 3 and AdS 5 thin Vaidya shells. We see that as expected the late time behaviour is determined by the first quasinormal mode. We apply the method to study the late time behaviour of the shear vector mode in AdS 5 Vaidya shell. At small momentum the corresponding time dependent Green function is expected to relax to equilibrium by the shear hydrodynamic mode. Using this we obtain the universal ratio of the shear viscosity to entropy density from a time dependent process.

Intertwined vorticity and elastodynamics in flapping wing propulsion

We performed numerical experiments on a one-dimensional elastic solid oscillating in a two-dimensional viscous incompressible fluid with the intent of discerning the interplay of vorticity and elastodynamics in flapping wing propulsion. Perhaps for the first time, we have established the role of foil deflection topology and its influence on vorticity generation, through spatially and temporally evolving foil slope and curvature. Though the frequency of oscillation of the foil has a definite role, it is the phase relation between foil slope and pressure that determines thrust or drag. Similarly, the phase difference between flapping velocity, and pressure and inertial forces, determine the power input to the foil, and in turn drives propulsive efficiency. At low frequencies of oscillation, the sympathetic slope and curvature of deformation of the foil allow generation of leading-edge vortices that do not separate; they cause substantial rise in pressure between the leading edge and mid-chord. The circulatory component of pressure is determined primarily by the leading-edge vortex and therefore thrust too is predominantly circulatory in origin at low frequencies. In the intermediate and high-frequency range, thrust and drag on the foil spatially alternate and non-circulatory forces dominate over circulatory and viscous forces. For the mass ratios we simulated, thrust due to flapping varies quadratically as a function of Strouhal number or trailing-edge flapping velocity; further, the trailing edge flapping velocities peak at the same set of frequencies where the thrust is also a maximum. Propulsive efficiency, on the other hand, is roughly a mirror image of the thrust variation with respect to Strouhal number. Given that most instances of flapping propulsion in nature are primarily through distributed muscular actuation that enables precise control of deformation shape, leading to high thrust and efficiency, the results presented here are pointers towards understanding some of the mechanisms that drive thrust and propulsive efficiency.

A flame particle tracking analysis of turbulence-chemistry interaction in hydrogen-air premixed flames

Interactions of turbulence, molecular transport, and energy transport, coupled with chemistry play a crucial role in the evolution of flame surface geometry, propagation, annihilation, and local extinction/re-ignition characteristics of intensely turbulent premixed flames. This study seeks to understand how these interactions affect flame surface annihilation of lean hydrogen-air premixed turbulent flames. Direct numerical simulations (DNSs) are conducted at different parametric conditions with a detailed reaction mechanism and transport properties for hydrogen-air flames. Flame particle tracking (FPT) technique is used to follow specific flame surface segments. An analytical expression for the local displacement flame speed (S-d) of a temperature isosurface is considered, and the contributions of transport, chemistry, and kinematics on the displacement flame speed at different turbulence-flame interaction conditions are identified. In general, the displacement flame speed for the flame particles is found to increase with time for all conditions considered. This is because, eventually all flame surfaces and their resident flame particles approach annihilation by reactant island formation at the end of stretching and folding processes induced by turbulence. Statistics of principal curvature evolving in time, obtained using FPT, suggest that these islands are ellipsoidal on average enclosing fresh reactants. Further examinations show that the increase in S-d is caused by the increased negative curvature of the flame surface and eventual homogenization of temperature gradients as these reactant islands shrink due to flame propagation and turbulent mixing. Finally, the evolution of the normalized, averaged, displacement flame speed vs. stretch Karlovitz number are found to collapse on a narrow band, suggesting that a unified description of flame speed dependence on stretch rate may be possible in the Lagrangian description. (C) 2015 The Combustion Institute. Published by Elsevier Inc. All rights reserved.

Effect of applied pressure on densification of monolithic ZrCx ceramic by reactive hot pressing

The effect of applied pressure on reactive hot pressing (RHP) of zirconium (Zr):graphite (C) in molar ratios of 1:0.5, 1:0.67, 1:0.8, and 1:1 was studied at 1200 degrees C for 60 min. The relative density achievable increased with increasing pressure and ranged from 99% at 4 MPa for ZrC0.5 to 93% for stoichiometric ZrC at 100 MPa. The diminishing influence of pressure on the final density with increasing stoichiometry is attributed to two causes: the decreasing initial volume fraction of the plastically deforming Zr metal which leads to the earlier formation of a contiguous, stress shielding carbide skeleton and the larger molar volume shrinkage during reaction which leads to pore formation in the final stages. A numerical model of the creep densification of a dynamically evolving microstructure predicts densities that are consistent with observations and confirm that the availability of a soft metal is primarily responsible for the achievement of such elevated densification during RHP. The ability to densify nonstoichiometric compositions like ZrC0.5 at pressures as low as 4 MPa offers an alternate route to fabricating dense nonstoichiometric carbides.