Numerical study of laminar, standing hydraulic jumps in a planar geometry

We solve the two-dimensional, planar Navier-Stokes equations to simulate a laminar, standing hydraulic jump using a Volume-of-Fluid method. The geometry downstream of the jump has been designed to be similar to experimental conditions by including a pit at the edge of the platform over which liquid film flows. We obtain jumps with and without separation. Increasing the inlet Froude number pushes the jump downstream and makes the slope of the jump weaker, consistent with experimental observations of circular jumps, and decreasing the Reynolds number brings the jump upstream while making it steeper. We study the effect of the length of the domain and that of a downstream obstacle on the structure and location of the jump. The transient flow which leads to a final steady jump is described for the first time to our knowledge. In the moderate Reynolds number regime, we obtain steady undular jumps with a separated bubble underneath the first few undulations. Interestingly, surface tension leads to shortening of wavelength of these undulations. We show that the undulations can be explained using the inviscid theory of Benjamin and Lighthill (Proc. R. Soc. London, Ser. A, 1954). We hope this new finding will motivate experimental verification.

Integrating AlGaN/GaN high electron mobility transistor with Si: A comparative study of integration schemes

AlGaN/GaN high electron mobility transistor stacks deposited on a single growth platform are used to compare the most common transition, AlN to GaN, schemes used for integrating GaN with Si. The efficiency of these transitions based on linearly graded, step graded, interlayer, and superlattice schemes on dislocation density reduction, stress management, surface roughness, and eventually mobility of the 2D-gas are evaluated. In a 500 nm GaN probe layer deposited, all of these transitions result in total transmission electron microscopy measured dislocations densities of 1 to 3 x 10(9)/cm(2) and <1 nm surface roughness. The 2-D electron gas channels formed at an AlGaN-1 nm AlN/GaN interface deposited on this GaN probe layer all have mobilities of 1600-1900 cm(2)/V s at a carrier concentration of 0.7-0.9 x 10(13)/cm(2). Compressive stress and changes in composition in GaN rich regions of the AlN-GaN transition are the most effective at reducing dislocation density. Amongst all the transitions studied the step graded transition is the one that helps to implement this feature of GaN integration in the simplest and most consistent manner. (C) 2015 AIP Publishing LLC.

A modified approach to numerical solution of Fredholm integral equations of the second kind

A modified approach to obtain approximate numerical solutions of Fredholin integral equations of the second kind is presented. The error bound is explained by the aid of several illustrative examples. In each example, the approximate solution is compared with the exact solution, wherever possible, and an excellent agreement is observed. In addition, the error bound in each example is compared with the one obtained by the Nystrom method. It is found that the error bound of the present method is smaller than the ones obtained by the Nystrom method. Further, the present method is successfully applied to derive the solution of an integral equation arising in a special Dirichlet problem. (C) 2015 Elsevier Inc. All rights reserved.

ASYMPTOTIC PRESERVING SCHEME FOR A KINETIC MODEL DESCRIBING IN COMPRESSIBLE FLUIDS

The kinetic theory of fluid turbulence modeling developed by Degond and Lemou in 7] is considered for further study, analysis and simulation. Starting with the Boltzmann like equation representation for turbulence modeling, a relaxation type collision term is introduced for isotropic turbulence. In order to describe some important turbulence phenomenology, the relaxation time incorporates a dependency on the turbulent microscopic energy and this makes difficult the construction of efficient numerical methods. To investigate this problem, we focus here on a multi-dimensional prototype model and first propose an appropriate change of frame that makes the numerical study simpler. Then, a numerical strategy to tackle the stiff relaxation source term is introduced in the spirit of Asymptotic Preserving Schemes. Numerical tests are performed in a one-dimensional framework on the basis of the developed strategy to confirm its efficiency.

Weakly corrected numerical solutions to stochastically driven nonlinear dynamical systems

A method to weakly correct the solutions of stochastically driven nonlinear dynamical systems, herein numerically approximated through the Eule-Maruyama (EM) time-marching map, is proposed. An essential feature of the method is a change of measures that aims at rendering the EM-approximated solution measurable with respect to the filtration generated by an appropriately defined error process. Using Ito's formula and adopting a Monte Carlo (MC) setup, it is shown that the correction term may be additively applied to the realizations of the numerically integrated trajectories. Numerical evidence, presently gathered via applications of the proposed method to a few nonlinear mechanical oscillators and a semi-discrete form of a 1-D Burger's equation, lends credence to the remarkably improved numerical accuracy of the corrected solutions even with relatively large time step sizes. (C) 2015 Elsevier Inc. All rights reserved.

Numerical Study of Counterflow Diffusion Flame and Water Spray Interaction

This paper reports numerical investigation concerning the interaction of a laminar methane-air counterflow diffusion flame with monodisperse and polydisperse water spray. Commercial code ANSYS FLUENT with reduced chemistry has been used for investigation. Effects of strain rate, Sauter mean diameter (SMD), and droplet size distribution on the temperature along stagnation streamline have been studied. Flame extinction using polydisperse water spray has also been explored. Comparison of monodisperse and polydisperse droplet distribution on flame properties reveals suitability of polydisperse spray in flame temperature reduction beyond a particular SMD. This study also provides a numerical framework to study flame-spray interaction and extinction.

Numerical analysis of segmented-electrode Orbitraps

Simple geometries which are possible alternatives for the Orbitrap are studied in this paper. We have taken up for numerical investigation two segmented-electrode structures, ORB1 and ORB2, to mimic the electric field of the Orbitrap. In the ORB1, the inner spindle-like electrode and the outer barrel-like electrode of the Orbitrap have been replaced by 35 rings and 35 discs of fixed radii, respectively. In this structure two segmented end cap electrodes have been added. In this geometry, different potentials are applied to the different electrodes keeping top-bottom symmetry intact. In the second geometry, ORB2, the inner and outer electrodes of the Orbitrap were replaced by an approximate step structure which follows the profile of the Orbitrap electrodes. In the present study 45 steps have been used. In the ORB2, like the Orbitrap, the inner electrode is held at a negative potential and the outer electrode is at ground potential. For the purpose of comparing the performance of ORB1 and ORB2 with that of the Orbitrap, the following studies have been undertaken: (1) variation of electric potential, (2) computation of ion trajectories, (3) simulation of image currents. These studies have been carried out using both 2D and 3D Boundary Element Method (BEM), the 3D BEM was developed specifically for this study. It has been seen in these investigations that ORB1 and ORB2 have performance similar to that of the Orbitrap, with the performance of the ORB1 being seen to be marginally superior to that of the ORB2. It has been shown that with proper optimization, geometries containing far fewer electrodes can be used as mass analyzers. A novel technique of optimization of the electric field has been proposed with the objective of minimizing the dependence of axial frequency of ion motion on the initial position of an ion. The results on the optimization of 9 and 15 segmented-electrode traps having the same design as ORB1 show that it can provide accurate mass analysis. (C) 2015 Elsevier B.V. All rights reserved.

The time finite element as a robust general scheme for solving nonlinear dynamic equations including chaotic systems

Schemes that can be proven to be unconditionally stable in the linear context can yield unstable solutions when used to solve nonlinear dynamical problems. Hence, the formulation of numerical strategies for nonlinear dynamical problems can be particularly challenging. In this work, we show that time finite element methods because of their inherent energy momentum conserving property (in the case of linear and nonlinear elastodynamics), provide a robust time-stepping method for nonlinear dynamic equations (including chaotic systems). We also show that most of the existing schemes that are known to be robust for parabolic or hyperbolic problems can be derived within the time finite element framework; thus, the time finite element provides a unification of time-stepping schemes used in diverse disciplines. We demonstrate the robust performance of the time finite element method on several challenging examples from the literature where the solution behavior is known to be chaotic. (C) 2015 Elsevier Inc. All rights reserved.

The time finite element as a robust general scheme for solving nonlinear dynamic equations including chaotic systems

Schemes that can be proven to be unconditionally stable in the linear context can yield unstable solutions when used to solve nonlinear dynamical problems. Hence, the formulation of numerical strategies for nonlinear dynamical problems can be particularly challenging. In this work, we show that time finite element methods because of their inherent energy momentum conserving property (in the case of linear and nonlinear elastodynamics), provide a robust time-stepping method for nonlinear dynamic equations (including chaotic systems). We also show that most of the existing schemes that are known to be robust for parabolic or hyperbolic problems can be derived within the time finite element framework; thus, the time finite element provides a unification of time-stepping schemes used in diverse disciplines. We demonstrate the robust performance of the time finite element method on several challenging examples from the literature where the solution behavior is known to be chaotic. (C) 2015 Elsevier Inc. All rights reserved.

Scaling analysis and numerical studies on water vapour adsorption in a columnar porous silica gel bed

This article presents a theoretical analysis of heat and mass transfer in a silica gel + water adsorption process using scaling principles. A two-dimensional columnar packed adsorber domain is chosen for the study, with side and bottom walls cooled and vapour inlet from the top. The adsorption process is initiated from the cold walls with a temperature jump of 15 K, whereas the water vapour supply is maintained at a constant inlet pressure of 1 kPa. The first part of the study is dedicated to deriving relevant scales for the adsorption process by an order of magnitude analysis of energy, continuity and momentum equations. In the latter part, the derived scales are compared with the outcome of numerical studies performed for various domain widths and aspect ratio of bed. A good correlation between scaling and simulation results is observed, thereby validating the scaling approach. (C) 2015 Elsevier Ltd. All rights reserved.

Numerical analysis of segmented-electrode Orbitraps

Simple geometries which are possible alternatives for the Orbitrap are studied in this paper. We have taken up for numerical investigation two segmented-electrode structures, ORB1 and ORB2, to mimic the electric field of the Orbitrap. In the ORB1, the inner spindle-like electrode and the outer barrel-like electrode of the Orbitrap have been replaced by 35 rings and 35 discs of fixed radii, respectively. In this structure two segmented end cap electrodes have been added. In this geometry, different potentials are applied to the different electrodes keeping top-bottom symmetry intact. In the second geometry, ORB2, the inner and outer electrodes of the Orbitrap were replaced by an approximate step structure which follows the profile of the Orbitrap electrodes. In the present study 45 steps have been used. In the ORB2, like the Orbitrap, the inner electrode is held at a negative potential and the outer electrode is at ground potential. For the purpose of comparing the performance of ORB1 and ORB2 with that of the Orbitrap, the following studies have been undertaken: (1) variation of electric potential, (2) computation of ion trajectories, (3) simulation of image currents. These studies have been carried out using both 2D and 3D Boundary Element Method (BEM), the 3D BEM was developed specifically for this study. It has been seen in these investigations that ORB1 and ORB2 have performance similar to that of the Orbitrap, with the performance of the ORB1 being seen to be marginally superior to that of the ORB2. It has been shown that with proper optimization, geometries containing far fewer electrodes can be used as mass analyzers. A novel technique of optimization of the electric field has been proposed with the objective of minimizing the dependence of axial frequency of ion motion on the initial position of an ion. The results on the optimization of 9 and 15 segmented-electrode traps having the same design as ORB1 show that it can provide accurate mass analysis. (C) 2015 Elsevier B.V. All rights reserved.

Numerical Investigations for Leakage and Windage Heating in Straight-Through Labyrinth Seals

The ability to quantify leakage flow and windage heating for labyrinth seals with honeycomb lands is critical in understanding gas turbine engine system performance and predicting its component life. Variety of labyrinth seal configurations (number of teeth, stepped or straight, honeycomb cell size) are in use in gas turbines, and for each configuration, there are many geometric factors that can impact a seal's leakage and windage characteristics. This paper describes the development of a numerical methodology aimed at studying the effect of honeycomb lands on leakage and windage heating. Specifically, a three-dimensional computational fluid dynamics (CFD) model is developed utilizing commercial finite volume-based software incorporating the renormalization group (RNG) k-epsilon turbulence model with modified Schmidt number. The modified turbulence model is benchmarked and fine-tuned based on several experiments. Using this model, a broad parametric study is conducted by varying honeycomb cell size, pressure ratio (PR), and radial clearance for a four-tooth straight-through labyrinth seal. The results show good agreement with available experimental data. They further indicate that larger honeycomb cells predict higher seal leakage and windage heating at tighter clearances compared to smaller honeycomb cells and smooth lands. However, at open seal clearances larger honeycomb cells have lower leakage compared to smaller honeycomb cells.

Numerical prediction of nozzle flow separation: Issue of turbulence modeling

Numerical simulation of separated flows in rocket nozzles is challenging because existing turbulence models are unable to predict it correctly. This paper addresses this issue with the Spalart-Allmaras and Shear Stress Transport (SST) eddy-viscosity models, which predict flow separation with moderate success. Their performances have been compared against experimental data for a conical and two contoured subscale nozzles. It is found that they fail to predict the separation location correctly, exhibiting sensitivity to the nozzle pressure ratio (NPR) and nozzle type. A careful assessment indicated how the model had to be tuned for better, consistent prediction. It is learnt that SST model's failure is caused by limiting of the shear stress inside boundary layer according to Bradshaw's assumption, and by over prediction of jet spreading rate. Accordingly, SST's coefficients were empirically modified to match the experimental wall pressure data. Results confirm that accurate RANS prediction of separation depends on the correct capture of the jet spreading rate, and that it is feasible over a wide range of NPRs by modified values of the diffusion coefficients in the turbulence model. (C) 2015 Elsevier Masson SAS. All rights reserved.

Analysis of sine-triangle and zero-sequence injection modulation schemes for split-phase induction motor drive

A split-phase induction motor is fed from two three-phase voltage source inverters for speed control. This study analyses carrier-comparison based pulse width modulation (PWM) schemes for a split-phase motor drive, from a space-vector perspective. Sine-triangle PWM, one zero-sequence injection PWM where the same zero-sequence signal is used for both the inverters, and another zero-sequence injection PWM where different zero-sequence signals are employed for the two inverters are considered. The set of voltage vectors applied, the sequence in which the voltage vectors are applied, and the resulting current ripple vector are analysed for all the PWM methods. Besides all the PWM methods are compared in terms of dc bus utilisation. For the same three-phase sine reference, the PWM method with different zero-sequence signals for the two inverters is found to employ a set of vectors different from the other methods. Both analysis and experimental results show that this method results in lower total harmonic distortion and higher dc bus utilisation than the other two PWM methods.

Numerical study of driven flows of shear thinning viscoelastic fluids in rectangular cavities

In this work, we present a numerical study of flow of shear thinning viscoelastic fluids in rectangular lid driven cavities for a wide range of aspect ratios (depth to width ratio) varying from 1/16 to 4. In particular, the effect of elasticity, inertia, model parameters and polymer concentration on flow features in rectangular driven cavity has been studied for two shear thinning viscoelastic fluids, namely, Giesekus and linear PTT. We perform numerical simulations using the symmetric square root representation of the conformation tensor to stabilize the numerical scheme against the high Weissenberg number problem. The variation in flow structures associated with merging and splitting of elongated vortices in shallow cavities and coalescence of corner eddies to yield a second primary vortex in deep cavities with respect to the variation in flow parameters is discussed. We discuss the effect of the dominant eigenvalues and the corresponding eigenvectors on the location of the primary eddy in the cavity. We also demonstrate, by performing numerical simulations for shallow and deep cavities, that where the Deborah number (based on convective time scale) characterizes the elastic behaviour of the fluid in deep cavities, Weissenberg number (based on shear rate) should be used for shallow cavities. (C) 2016 Elsevier B.V. All rights reserved.