Double diffusive Marangoni convection flow of electrically conducting fluid in a square cavity with chemical reaction

Double diffusive Marangoni convection flow of viscous incompressible electrically conducting fluid in a square cavity is studied in this paper by taking into consideration of the effect of applied magnetic field in arbitrary direction and the chemical reaction. The governing equations are solved numerically by using alternate direct implicit (ADI) method together with the successive over relaxation (SOR) technique. The flow pattern with the effect of governing parameters, namely the buoyancy ratio W, diffusocapillary ratio w, and the Hartmann number Ha, is investigated. It is revealed from the numerical simulations that the average Nusselt number decreases; whereas the average Sherwood number increases as the orientation of magnetic field is shifted from horizontal to vertical. Moreover, the effect of buoyancy due to species concentration on the flow is stronger than the one due to thermal buoyancy. The increase in diffusocapillary parameter, w caus

Long, vertically aligned single-walled carbon nanotubes from plasmas: Morpho-kinetic and alignment controls

The enhanced large-scale model and numerical simulations are used to clarify the growth mechanism and the differences between the plasma- and neutral gas-grown carbon nanotubes, and to reveal the underlying physics and the key growth parameters. The results show that the nanotubes grown by plasma can be longer due to the effects of hydrocarbon ions with velocities aligned with the nanotubes. We show that the low-temperature growth is possible when the hydrocarbon ion flux dominates over fluxes of other species. We have also analysed the dependencies of the nanotube growth rates on nanotube and process parameters. The results are verified by a direct comparison with the experimental data. The model is generic and can be used for other types of carbon nanostructures such as carbon nanowalls, vertical graphenes, etc.

Two-scale computational modelling of water flow in unsaturated soils containing irregular-shaped inclusions

The focus of this paper is two-dimensional computational modelling of water flow in unsaturated soils consisting of weakly conductive disconnected inclusions embedded in a highly conductive connected matrix. When the inclusions are small, a two-scale Richards’ equation-based model has been proposed in the literature taking the form of an equation with effective parameters governing the macroscopic flow coupled with a microscopic equation, defined at each point in the macroscopic domain, governing the flow in the inclusions. This paper is devoted to a number of advances in the numerical implementation of this model. Namely, by treating the micro-scale as a two-dimensional problem, our solution approach based on a control volume finite element method can be applied to irregular inclusion geometries, and, if necessary, modified to account for additional phenomena (e.g. imposing the macroscopic gradient on the micro-scale via a linear approximation of the macroscopic variable along the microscopic boundary). This is achieved with the help of an exponential integrator for advancing the solution in time. This time integration method completely avoids generation of the Jacobian matrix of the system and hence eases the computation when solving the two-scale model in a completely coupled manner. Numerical simulations are presented for a two-dimensional infiltration problem.

Natural convection in a triangular enclosure heated from below and non-uniformly cooled from top

This study is concerned with transient natural convection in an isosceles triangular enclosure subject to non-uniformly cooling at the inclined surfaces and uniformly heating at the base. The numerical simulations of the unsteady flows over a range of Rayleigh numbers and aspect ratios are carried out using Finite Volume Method. Since the upper inclined surfaces are linearly cooled and the bottom surface is heated, the flow is potentially unstable. It is revealed from the numerical simulations that the transient flow development in the enclosure can be classified into three distinct stages; an early stage, a transitional stage, and a steady stage. The flow inside the enclosure depends significantly on the governing parameters, Rayleigh number and aspect ratio. The effect of Rayleigh number and aspect ratio on the flow development and heat transfer rate are discussed. The key finding for this study is to analyze the pitchfork bifurcation of the flow about the geometric center line. The heat transfer through the roof and the ceiling as a form of Nusselt number is reported in this study.

On the relationship between finger width, velocity, and fluxes in thermohaline convection

Double-diffusive finger convection occurs in many natural processes.The theories for double-diffusive phenomena that exist at present consider systems with linear stratification in temperature and salinity. The double-diffusive systems with step change in salinity and temperature are, however, not amenable to simple stability analysis. Hence factors that control the width of the finger, velocity, and fluxes in systems that have step change in temperature and salinity have not been understood so far. In this paper we provide new physical insight regarding factors that influence finger convection in two-layer double-diffusive system through two-dimensional numerical simulations. Simulations have been carried out for density stability ratios (R-rho) from 1.5 to 10. For each density stability ratio, the thermal Rayleigh number (Ra-T) has been systematically varied from 7x10(3) to 7x10(8). Results from these simulations show how finger width, velocity, and flux ratios in finger convection are interrelated and the influence of governing parameters such as density stability ratio and the thermal Rayleigh number. The width of the incipient fingers at the time of onset of instability has been shown to vary as Ra-T-1/3. Velocity in the finger varies as Ra(T)1/3/R-rho. Results from simulation agree with the scale analysis presented in the paper. Our results demonstrate that wide fingers have lower velocities and flux ratios compared to those in narrow fingers. This result contradicts present notions about the relation between finger width and flux ratio. A counterflow heat-exchanger analogy is used in understanding the dependence of flux ratio on finger width and velocity.

Investigation of Missile-Shaped Body with Forward-Facing Cavity at Mach 8

A blunt-nosed hypersonic missile mounted with a forward-facing cavity is a good alternative to reduce the stagnation heating rates. The effects of a forward-racing cavity on heat transfer and aerodynamic coefficients are addressed in this paper. Tests were carried out in hypersonic shock tunnel HST2, at a hypersonic Mach number of 8 using a 41 deg apex-angle blunt cone. The aerodynamic forces on the test model with and without a forward-facing cavity at various angles of attack are measured by using an internally mountable accelerometer force balance system. Heat flux measurements have been carried out on the test model with and without a forward-facing cavity of the entire surface at zero degree angle of attack with platinum sensors. A numerical simulation was also carried out using the computational fluid dynamics code (CFX-Ansys 5.7). An important result of this study is that the smaller cavity diameter has the highest lift-to-drag ratio, whereas the medium cavity has the highest heat flux reduction. Theshock structure around the test model has also been visualized using the Schlieren flow visualization technique. The visualized shock structure and the measured aerodynamic forces on the missile-shaped body with cavity configurations agree well with the axisymmetric numerical simulations.

A novel approach to 6-DOF adaptive trajectory tracking control of an AUV in the presence of parameter uncertainties

In this paper, the trajectory tracking control of an autonomous underwater vehicle (AUVs) in six-degrees-of-freedom (6-DOFs) is addressed. It is assumed that the system parameters are unknown and the vehicle is underactuated. An adaptive controller is proposed, based on Lyapunov׳s direct method and the back-stepping technique, which interestingly guarantees robustness against parameter uncertainties. The desired trajectory can be any sufficiently smooth bounded curve parameterized by time even if consist of straight line. In contrast with the majority of research in this field, the likelihood of actuators׳ saturation is considered and another adaptive controller is designed to overcome this problem, in which control signals are bounded using saturation functions. The nonlinear adaptive control scheme yields asymptotic convergence of the vehicle to the reference trajectory, in the presence of parametric uncertainties. The stability of the presented control laws is proved in the sense of Lyapunov theory and Barbalat׳s lemma. Efficiency of presented controller using saturation functions is verified through comparing numerical simulations of both controllers.

Bearing capacity of square footings on geosynthetic reinforced sand

The results from laboratory model tests and numerical simulations on square footings resting on sand are presented. Bearing capacity of footings on geosynthetic reinforced sand is evaluated and the effect of various reinforcement parameters like the type and tensile strength of geosynthetic material, amount of reinforcement, layout and configuration of geosynthetic layers below the footing on the bearing capacity improvement of the footings is studied through systemati model studies. A steel tank of size 900 x 900 x 600 mm is used for conducting model tests. Four types of grids, namely strong biaxial geogrid, weak biaxial geogrid, uniaxial geogrid and a geonet, each with different tensile strength, are used in the tests. Geosynthetic reinforcement is provided in the form of planar layers, varying the depth of reinforced zone below the footing, number of geosynthetic layers within the reinforced zone and the width of geosynthetic layers in different tests. Influence of all these parameters on the bearing capacity improvement of square footing and its settlement is studied by comparing with the test on unreinforced sand. Results show that the effective depth of reinforcement is twice the width of the footing and optimum spacing of geosynthetic layers is half the width of the footing. It is observed that the layout and configuration of reinforcement play a vital role in bearing capacity improvement rather than the tensile strength of the geosynthetic material. Experimental observations are supported by the findings from numerical analyses.

Generation of rotationally dominated galaxies by mergers of pressure-supported progenitors

Through the analysis of a set of numerical simulations of major mergers between initially non-rotating, pressure-supported progenitor galaxies with a range of central mass concentrations, we have shown that: (1) it is possible to generate elliptical-like galaxies, with outside one effective radius, as a result of the conversion of orbital- into internal-angular momentum; (2) the outer regions acquire part of the angular momentum first; (3) both the baryonic and the dark matter components of the remnant galaxy acquire part of the angular momentum, the relative fractions depending on the initial concentration of the merging galaxies. For this conversion to occur the initial baryonic component must be sufficiently dense and/or the encounter should take place on an orbit with high angular momentum. Systems with these hybrid properties have recently been observed through a combination of stellar absorption lines and planetary nebulae for kinematic studies of early-type galaxies. Our results are in qualitative agreement with these observations and demonstrate that even mergers composed of non rotating, pressure-supported progenitor galaxies can produce early-type galaxies with significant rotation at large radii.

Inverse Problem for Fractional Brownian Motion with Discrete Data

The problem of recovering information from measurement data has already been studied for a long time. In the beginning, the methods were mostly empirical, but already towards the end of the sixties Backus and Gilbert started the development of mathematical methods for the interpretation of geophysical data. The problem of recovering information about a physical phenomenon from measurement data is an inverse problem. Throughout this work, the statistical inversion method is used to obtain a solution. Assuming that the measurement vector is a realization of fractional Brownian motion, the goal is to retrieve the amplitude and the Hurst parameter. We prove that under some conditions, the solution of the discretized problem coincides with the solution of the corresponding continuous problem as the number of observations tends to infinity. The measurement data is usually noisy, and we assume the data to be the sum of two vectors: the trend and the noise. Both vectors are supposed to be realizations of fractional Brownian motions, and the goal is to retrieve their parameters using the statistical inversion method. We prove a partial uniqueness of the solution. Moreover, with the support of numerical simulations, we show that in certain cases the solution is reliable and the reconstruction of the trend vector is quite accurate.

Dynamic Instabilities in Slender Space Launch Vehicles under Propulsive Thrust and Aerodynamic Forces

A mechanics based linear analysis of the problem of dynamic instabilities in slender space launch vehicles is undertaken. The flexible body dynamics of the moving vehicle is studied in an inertial frame of reference, including velocity induced curvature effects, which have not been considered so far in the published literature. Coupling among the rigid-body modes, the longitudinal vibrational modes and the transverse vibrational modes due to asymmetric lifting-body cross-section are considered. The model also incorporates the effects of aerodynamic forces and the propulsive thrust of the vehicle. The effects of the coupling between the combustion process (mass variation, developed thrust etc.) and the variables involved in the flexible body dynamics (displacements and velocities) are clearly brought out. The model is one-dimensional, and it can be employed to idealised slender vehicles with complex shapes. Computer simulations are carried out using a standard eigenvalue problem within h-p finite element modelling framework. Stability regimes for a vehicle subjected to propulsive thrust are validated by comparing the results from published literature. Numerical simulations are carried out for a representative vehicle to determine the instability regimes with vehicle speed and propulsive thrust as the parameters. The phenomena of static instability (divergence) and dynamic instability (flutter) are observed. The results at low Mach number match closely with the results obtained from previous models published in the literature.

Spiral-Wave Turbulence and Its Control in the Presence of Inhomogeneities in Four Mathematical Models of Cardiac Tissue

Regular electrical activation waves in cardiac tissue lead to the rhythmic contraction and expansion of the heart that ensures blood supply to the whole body. Irregularities in the propagation of these activation waves can result in cardiac arrhythmias, like ventricular tachycardia (VT) and ventricular fibrillation (VF), which are major causes of death in the industrialised world. Indeed there is growing consensus that spiral or scroll waves of electrical activation in cardiac tissue are associated with VT, whereas, when these waves break to yield spiral- or scroll-wave turbulence, VT develops into life-threatening VF: in the absence of medical intervention, this makes the heart incapable of pumping blood and a patient dies in roughly two-and-a-half minutes after the initiation of VF. Thus studies of spiral- and scroll-wave dynamics in cardiac tissue pose important challenges for in vivo and in vitro experimental studies and for in silico numerical studies of mathematical models for cardiac tissue. A major goal here is to develop low-amplitude defibrillation schemes for the elimination of VT and VF, especially in the presence of inhomogeneities that occur commonly in cardiac tissue. We present a detailed and systematic study of spiral- and scroll-wave turbulence and spatiotemporal chaos in four mathematical models for cardiac tissue, namely, the Panfilov, Luo-Rudy phase 1 (LRI), reduced Priebe-Beuckelmann (RPB) models, and the model of ten Tusscher, Noble, Noble, and Panfilov (TNNP). In particular, we use extensive numerical simulations to elucidate the interaction of spiral and scroll waves in these models with conduction and ionic inhomogeneities; we also examine the suppression of spiral- and scroll-wave turbulence by low-amplitude control pulses. Our central qualitative result is that, in all these models, the dynamics of such spiral waves depends very sensitively on such inhomogeneities. We also study two types of control chemes that have been suggested for the control of spiral turbulence, via low amplitude current pulses, in such mathematical models for cardiac tissue; our investigations here are designed to examine the efficacy of such control schemes in the presence of inhomogeneities. We find that a local pulsing scheme does not suppress spiral turbulence in the presence of inhomogeneities; but a scheme that uses control pulses on a spatially extended mesh is more successful in the elimination of spiral turbulence. We discuss the theoretical and experimental implications of our study that have a direct bearing on defibrillation, the control of life-threatening cardiac arrhythmias such as ventricular fibrillation.

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.

Critical lattice size limit for synchronized chaotic state in one- and two-dimensional diffusively coupled map lattices

We consider diffusively coupled map lattices with P neighbors (where P is arbitrary) and study the stability of the synchronized state. We show that there exists a critical lattice size beyond which the synchronized state is unstable. This generalizes earlier results for nearest neighbor coupling. We confirm the analytical results by performing numerical simulations on coupled map lattices with logistic map at each node. The above analysis is also extended to two-dimensional P-neighbor diffusively coupled map lattices.

Generalized survival in step fluctuations

The properties of the generalized survival probability, that is, the probability of not crossing an arbitrary location R during relaxation, have been investigated experimentally (via scanning tunneling microscope observations) and numerically. The results confirm that the generalized survival probability decays exponentially with a time constant tau(s)(R). The distance dependence of the time constant is shown to be tau(s)(R)=tau(s0)exp[-R/w(T)], where w(2)(T) is the material-dependent mean-squared width of the step fluctuations. The result reveals the dependence on the physical parameters of the system inherent in the prior prediction of the time constant scaling with R/L-alpha, with L the system size and alpha the roughness exponent. The survival behavior is also analyzed using a contrasting concept, the generalized inside survival S-in(t,R), which involves fluctuations to an arbitrary location R further from the average. Numerical simulations of the inside survival probability also show an exponential time dependence, and the extracted time constant empirically shows (R/w)(lambda) behavior, with lambda varying over 0.6 to 0.8 as the sampling conditions are changed. The experimental data show similar behavior, and can be well fit with lambda=1.0 for T=300 K, and 0.5