Experimental and Numerical Studies in a Compact Trapped Vortex Combustor: Stability Assessment and Augmentation

Fundamental studies on a compact trapped vortex combustor indicate that cavity injection strategies play a major role on flame stability. Detailed experiments indicate that blow-out occurs for a certain range of cavity air flow velocities. An unsteady RANS-based reacting flow simulation tool has been utilized to study the basic dynamics of cavity vortex for various flow conditions. The phenomenon of flame blow-out at certain intermediate cavity air velocities is explained on the basis of transition from a cavity-stabilized mode to an opposed flow stagnation mode. A novel strategy is proposed for achieving flame stability at all conditions. This involves using a flow guide vane in the path of the main flow to direct a portion of the main flow into the cavity. This seems to result in a desirable dual vortex structure, i.e., a small clockwise vortex behind the vane and large counterclockwise vortex in the cavity. Experimental results show stable flame at all flow conditions with the flow guide vane, and pressure drop is estimated to be within acceptable limits. Cold flow simulations show self-similar velocity profiles for a range of main inlet velocities, and high reverse velocity ratios (-0.3) are observed. Such a high-velocity ratio in the reverse flow shear layer profile leads to enhanced production of turbulence imperative to compact combustors. Reacting flow simulations show even higher reverse velocity ratios (above -0.7) due to flow acceleration. The flame is observed to be stable, even though minor shear layer oscillations are present in the form of vortex shedding. Self-similarity is also observed in reacting flow temperature profiles at combustor exit over the entire range of the mainstream velocity. This indicates that the present configuration holds a promise of delivering robust performance invariant of the flow operating conditions.

A simple numerical method for three-dimensional analysis of simple expansion chamber mufflers of rectangular as well as circular cross-section with a stationary medium

Three-dimensional effects are a primary source of discrepancy between the measured values of automotive muffler performance and those predicted by the plane wave theory at higher frequencies. The basically exact method of (truncated) eigenfunction expansions for simple expansion chambers involves very complicated algebra, and the numerical finite element method requires large computation time and core storage. A simple numerical method is presented in this paper. It makes use of compatibility conditions for acoustic pressure and particle velocity at a number of equally spaced points in the planes of the junctions (or area discontinuities) to generate the required number of algebraic equations for evaluation of the relative amplitudes of the various modes (eigenfunctions), the total number of which is proportional to the area ratio. The method is demonstrated for evaluation of the four-pole parameters of rigid-walled, simple expansion chambers of rectangular as well as circular cross-section for the case of a stationary medium. Computed values of transmission loss are compared with those computed by means of the plane wave theory, in order to highlight the onset (cutting-on) of various higher order modes and the effect thereof on transmission loss of the muffler. These are also compared with predictions of the finite element methods (FEM) and the exact methods involving eigenfunction expansions, in order to demonstrate the accuracy of the simple method presented here.

On numerical prediction of the acoustic source characteristics of an engine exhaust system

The contemporary methods for source characterization rely mainly on experiments. These methods produce inaccurate results in the low‐frequency band, where the characteristics are all the more important. Moreover, the experimental methods cannot be used at the design stage. Hence, a numerical technique to obtain the source characteristics is desirable. In this paper, the pressure‐time history and the mass‐flux‐time history obtained by means of the time‐domain analysis have been used, along with the two‐load method to compute the source characteristics. Two new computational methods for obtaining the source characteristics have been described. These are much simpler, and computationally more economical than the complete time‐domain simulation, which makes use of the method of characteristics.

Accelerating Numerical Linear Algebra Kernels on a Scalable Run Time Reconfigurable Platform

Numerical Linear Algebra (NLA) kernels are at the heart of all computational problems. These kernels require hardware acceleration for increased throughput. NLA Solvers for dense and sparse matrices differ in the way the matrices are stored and operated upon although they exhibit similar computational properties. While ASIC solutions for NLA Solvers can deliver high performance, they are not scalable, and hence are not commercially viable. In this paper, we show how NLA kernels can be accelerated on REDEFINE, a scalable runtime reconfigurable hardware platform. Compared to a software implementation, Direct Solver (Modified Faddeev's algorithm) on REDEFINE shows a 29X improvement on an average and Iterative Solver (Conjugate Gradient algorithm) shows a 15-20% improvement. We further show that solution on REDEFINE is scalable over larger problem sizes without any notable degradation in performance.

Numerical analysis of transpiration cooling in carbon dioxide atmosphere at hypersonic Mach numbers

The numerical solutions are obtained for skin friction, heat transfer to the wall and growth of boundary layer along the flat plate by employing two dimensional Navier-Stokes equations governing the hypersonic flow coupled with species continuity equations. Flow fields have been computed along the flat plate in CO2 atmosphere in the presence of transpiration cooling using air and carbon dioxide.

Effect of capillarity and heterogeneity in the numerical modelling of multiphase flow of fluids in unsaturated porous medium

The multiphase flow of fluids in the unsaturated porous medium is considered as a three phase flow of water, NAPL, and air simultaneously in the porous medium. The adaptive solution fully implicit modified sequential method is used for the numerical modelling. The effect of capillarity and heterogeneity effect at the interface between the media is studied and it is observed that the interface criteria has to be taken into account for the correct prediction of NAPL migration especially in heterogeneous media. The modified Newton Raphson method is used for the linearization and Hestines and Steifel Conjugate Gradient method is used as the solver.

Numerical simulation of nano-indentation on rough surfaces

Nano-indentation is a technique used to measure various mechanical properties like hardness, Young's modulus and the adherence of thin films and surface layers. It can be used as a quality control tool for various surface modification techniques like ion-implantation, film deposition processes etc. It is important to characterise the increasing scatter in the data measured at lower penetration depths observed in the nano-indentation, for the technique to be effectively applied. Surface roughness is one of the parameters contributing for the scatter. This paper is aimed at quantifying the nature and the amount of scatter that will be introduced in the measurement due to the roughness of the surface on which the indentation is carried out. For this the surface is simulated using the Weierstrass-Mandelbrot function which gives a self-affine fractal. The contact area of this surface with a conical indenter with a spherical cap at the tip is measured numerically. The indentation process is simulated using the spherical cavity model. This eliminates the indentation size effect observed at the micron and sub-micron scales. It has been observed that there exists a definite penetration depth in relation to the surface roughness beyond which the scatter is reduced such that reliable data could be obtained.

A numerical study of variable coefficient elliptic Cauchy problem via projection method

In this paper, we investigate a numerical method for the solution of an inverse problem of recovering lacking data on some part of the boundary of a domain from the Cauchy data on other part for a variable coefficient elliptic Cauchy problem. In the process, the Cauchy problem is transformed into the problem of solving a compact linear operator equation. As a remedy to the ill-posedness of the problem, we use a projection method which allows regularization solely by discretization. The discretization level plays the role of regularization parameter in the case of projection method. The balancing principle is used for the choice of an appropriate discretization level. Several numerical examples show that the method produces a stable good approximate solution.

Studying the resistivity imaging of chicken tissue phantoms with different current patterns in Electrical Impedance Tomography (EIT)

A current injection pattern in Electrical Impedance Tomography (EIT) has its own current distribution profile within the domain under test. Hence, different current patterns have different sensitivity, spatial resolution and distinguishability. Image reconstruction studies with practical phantoms are essential to assess the performance of EIT systems for their validation, calibration and comparison purposes. Impedance imaging of real tissue phantoms with different current injection methods is also essential for better assessment of the biomedical EIT systems. Chicken tissue paste phantoms and chicken tissue block phantoms are developed and the resistivity image reconstruction is studied with different current injection methods. A 16-electrode array is placed inside the phantom tank and the tank is filled with chicken muscle tissue paste or chicken tissue blocks as the background mediums. Chicken fat tissue, chicken bone, air hole and nylon cylinders are used as the inhomogeneity to obtained different phantom configurations. A low magnitude low frequency constant sinusoidal current is injected at the phantom boundary with opposite and neighboring current patterns and the boundary potentials are measured. Resistivity images are reconstructed from the boundary data using EIDORS and the reconstructed images are analyzed with the contrast parameters calculated from their elemental resistivity profiles. Results show that the resistivity profiles of all the phantom domains are successfully reconstructed with a proper background resistivity and high inhomogeneity resistivity for both the current injection methods. Reconstructed images show that, for all the chicken tissue phantoms, the inhomogeneities are suitably reconstructed with both the current injection protocols though the chicken tissue block phantom and opposite method are found more suitable. It is observed that the boundary potentials of the chicken tissue block phantoms are higher than the chicken tissue paste phantom. SNR of the chicken tissue block phantoms are found comparatively more and hence the chicken tissue block phantom is found more suitable for its lower noise performance. The background noise is found less in opposite method for all the phantom configurations which yields the better resistivity images with high PCR and COC and proper IRMean and IRMax neighboring method showed higher noise level for both the chicken tissue paste phantoms and chicken tissue block phantoms with all the inhomogeneities. Opposite method is found more suitable for both the chicken tissue phantoms, and also, chicken tissue block phantoms are found more suitable compared to the chicken tissue paste phantom. (C) 2012 Elsevier Ltd. All rights reserved.

A pseudo-time EnKF incorporating shape based reconstruction for diffuse optical tomography

Purpose: The authors aim at developing a pseudo-time, sub-optimal stochastic filtering approach based on a derivative free variant of the ensemble Kalman filter (EnKF) for solving the inverse problem of diffuse optical tomography (DOT) while making use of a shape based reconstruction strategy that enables representing a cross section of an inhomogeneous tumor boundary by a general closed curve. Methods: The optical parameter fields to be recovered are approximated via an expansion based on the circular harmonics (CH) (Fourier basis functions) and the EnKF is used to recover the coefficients in the expansion with both simulated and experimentally obtained photon fluence data on phantoms with inhomogeneous inclusions. The process and measurement equations in the pseudo-dynamic EnKF (PD-EnKF) presently yield a parsimonious representation of the filter variables, which consist of only the Fourier coefficients and the constant scalar parameter value within the inclusion. Using fictitious, low-intensity Wiener noise processes in suitably constructed ``measurement'' equations, the filter variables are treated as pseudo-stochastic processes so that their recovery within a stochastic filtering framework is made possible. Results: In our numerical simulations, we have considered both elliptical inclusions (two inhomogeneities) and those with more complex shapes (such as an annular ring and a dumbbell) in 2-D objects which are cross-sections of a cylinder with background absorption and (reduced) scattering coefficient chosen as mu(b)(a)=0.01mm(-1) and mu('b)(s)=1.0mm(-1), respectively. We also assume mu(a) = 0.02 mm(-1) within the inhomogeneity (for the single inhomogeneity case) and mu(a) = 0.02 and 0.03 mm(-1) (for the two inhomogeneities case). The reconstruction results by the PD-EnKF are shown to be consistently superior to those through a deterministic and explicitly regularized Gauss-Newton algorithm. We have also estimated the unknown mu(a) from experimentally gathered fluence data and verified the reconstruction by matching the experimental data with the computed one. Conclusions: The PD-EnKF, which exhibits little sensitivity against variations in the fictitiously introduced noise processes, is also proven to be accurate and robust in recovering a spatial map of the absorption coefficient from DOT data. With the help of shape based representation of the inhomogeneities and an appropriate scaling of the CH expansion coefficients representing the boundary, we have been able to recover inhomogeneities representative of the shape of malignancies in medical diagnostic imaging. (C) 2012 American Association of Physicists in Medicine. [DOI: 10.1118/1.3679855]

Numerical Investigation of Heat Transfer from a Two-Dimensional Sudden Expansion Flow Using Nanofluids

Laminar forced convection heat transfer from two-dimensional sudden expansion flow of different nanofluids is studied numerically. The governing equations are solved using the unsteady stream function-vorticity method. The effect of volume fraction of the nanoparticles and type of nanoparticles on heat transfer is examined and found to have a significant impact. Local and average Nusselt numbers are reported in connection with various nanoparticle, volume fraction, and Reynolds number for expansion ratio 2. The Nusselt number reaches peak values near the reattachment point and reaches asymptotic value in the downstream. Bottom wall eddy and volume fraction shows a significant impact on the average Nusselt number.

A Numerical Scheme for Three-Dimensional Front Propagation and Control of Jordan Mode

As an example of a front propagation, we study the propagation of a three-dimensional nonlinear wavefront into a polytropic gas in a uniform state and at rest. The successive positions and geometry of the wavefront are obtained by solving the conservation form of equations of a weakly nonlinear ray theory. The proposed set of equations forms a weakly hyperbolic system of seven conservation laws with an additional vector constraint, each of whose components is a divergence-free condition. This constraint is an involution for the system of conservation laws, and it is termed a geometric solenoidal constraint. The analysis of a Cauchy problem for the linearized system shows that when this constraint is satisfied initially, the solution does not exhibit any Jordan mode. For the numerical simulation of the conservation laws we employ a high resolution central scheme. The second order accuracy of the scheme is achieved by using MUSCL-type reconstructions and Runge-Kutta time discretizations. A constrained transport-type technique is used to enforce the geometric solenoidal constraint. The results of several numerical experiments are presented, which confirm the efficiency and robustness of the proposed numerical method and the control of the Jordan mode.