Numerical analysis on thermoacoustic prime movers for development of pulse tube cryocoolers

Thermo Acoustic Prime Movers (TAPMs) are being considered as the ideal choice for driving the Pulse Tube Cryocoolers replacing the conventional compressors. The advantages are the absence of moving components and they can be driven by low grade energy as such as fuel, gas, solar energy, waste heat etc. While the development of such TAPMs is in progress in our laboratory, their design and fabrication should be guided by numerical modeling and this may be done by several methods such as solving the energy equation 1], enthalpy flow model 2], CFD 3], etc. We have used CFD technique, since it provides a better insight into the velocity and temperature profiles. The analysis is carried out by varying parameters such as (a) temperature difference across the stack, (b) stack and resonator lengths and (c) different working fluids such as air, nitrogen, argon etc. The theoretical results are compared with the experimental data wherever possible and they are in reasonably good agreement with each other. The analysis indicate that (i) larger temperature difference across the stack leads to increased acoustic amplitude, (ii) longer resonator leads to decrease in frequency with lesser amplitude and (iii) there exists an optimal stack length for the best performance of TAPM. These results are presented here.

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.

Numerical analysis of geocell reinforced ballast overlying soft clay subgrade

Geotextiles and geogrids have been in use for several decades in variety of geo-structure applications including foundation of embankments, retaining walls, pavements. Geocells is one such variant in geosynthetic reinforcement of recent years, which provides a three dimensional confinement to the infill material. Although extensive research has been carried on geocell reinforced sand, clay and layered soil subgrades, limited research has been reported on the aggregates/ballast reinforced with geocells. This paper presents the behavior of a railway sleeper subjected to monotonic loading on geocell reinforced aggregates, of size ranging from 20 to 75 mm, overlying soft clay subgrades. Series of tests were conducted in a steel test tank of dimensions 700 mm x 300 mm x 700 mm. In addition to the laboratory model tests, numerical simulations were performed using a finite difference code to predict the behavior of geocell reinforced ballast. The results from numerical simulations were compared with the experimental data. The numerical and experimental results manifested the importance that the geocell reinforcement has a significant effect on the ballast behaviour. The results depicted that the stiffness of underlying soft clay subgrade has a significant influence on the behavior of the geocell-aggregate composite material in redistributing the loading system.

Effect of Inlet Jet injection angle on a calandria based reactor-An investigation with numerical analysis of heat and mass transfer for an optimum reactor design

Heat and mass transfer studies in a calandria based reactor is quite complex both due to geometry and due to the complex mixing flow. It is challenging to devise optimum operating conditions with efficient but safe working range for such a complex configuration. Numerical study known to be very effective is taken up for investigation. In the present study a 3D RANS code with turbulence model has been used to compute the flow fields and to get the heat transfer characteristics to understand certain design parameters of engineering importance. The angle of injection and of the coolant liquid has a large effect on the heat transfer within the reactor.

Three-dimensional numerical analysis of geocell-reinforced soft clay beds by considering the actual geometry of geocell pockets

Due to its complex honeycomb structure, the numerical modeling of the geocell has always been a big challenge. Generally, the equivalent composite approach is used to model the geocells. In the equivalent composite approach, the geocellsoil composite is treated as the soil layer with improved strength and stiffness values. Though this approach is very simple, it is unrealistic to model the geocells as the soil layer. This paper presents a more realistic approach of modeling the geocells in three-dimensional (3D) framework by considering the actual curvature of the geocell pocket. A square footing resting on geocell reinforced soft clay bed was modeled using the ``fast Lagrangian analysis of continua in 3D'' (FLAC(3D)) finite difference package. Three different material models, namely modified Cam-clay, Mohr-Coulomb, and linear elastic were used to simulate the behaviour of foundation soil, infill soil and the geocell, respectively. It was found that the geocells distribute the load laterally to the wider area below the footing as compared to the unreinforced case. More than 50% reduction in the stress was observed in the clay bed in the presence of geocells. In addition to geocells, two other cases, namely, only geogrid and geocell with additional basal geogrid cases were also simulated. The numerical model was systematically validated with the results of the physical model tests. Using the validated numerical model, parametric studies were conducted to evaluate the influence of various geocell properties on the performance of reinforced clay beds.

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 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.

Error analysis of discontinuous Galerkin methods for the Stokes problem under minimal regularity

In this article, we analyse several discontinuous Galerkin (DG) methods for the Stokes problem under minimal regularity on the solution. We assume that the velocity u belongs to H-0(1)(Omega)](d) and the pressure p is an element of L-0(2)(Omega). First, we analyse standard DG methods assuming that the right-hand side f belongs to H-1(Omega) boolean AND L-1(Omega)](d). A DG method that is well defined for f belonging to H-1(Omega)](d) is then investigated. The methods under study include stabilized DG methods using equal-order spaces and inf-sup stable ones where the pressure space is one polynomial degree less than the velocity space.

Vertical uplift capacity of two interfering horizontal anchors in sand using an upper bound limit analysis

The vertical uplift resistance of two interfering rigid rough strip anchors embedded horizontally in sand at shallow depths has been examined. The analysis is performed by using an upper bound theorem o limit analysis in combination with finite elements and linear programming. It is specified that both the anchors are loaded to failure simultaneously at the same magnitude of the failure load. For different clear spacing (S) between the anchors, the magnitude of the efficiency factor (xi(gamma)) is determined. On account of interference, the magnitude of xi(gamma) is found to reduce continuously with a decrease in the spacing between the anchors. The results from the numerical analysis were found to compare reasonably well with the available theoretical data from the literature.

An Analysis of a Boundary Layer Diffusion Flame Over a Porous Flat Plate in a Confined Flow

A numerical analysis of the gas dynamic structure of a two-dimensional laminar boundary layer diffusion flame over a porous flat plate in a confined flow is made on the basis of the familiar boundary layer and flame sheet approximations neglecting buoyancy effects. The governing equations of aerothermochemistry with the appropriate boundary conditions are solved using the Patankar-Spalding method. The analysis predicts the flame shape, profiles of temperature, concentrations of variousspecies, and the density of the mixture across the boundary layer. In addition, it also predicts the pressure gradient in the flow direction arising from the confinement ofthe flow and the consequent velocity overshoot near the flame surface. The results of thecomputation performed for an n-pentane-air system are compared with experimental data andthe agreement is found to be satisfactory.

Reliability Analysis of Earth Dams

In the present study, results of reliability analyses of four selected rehabilitated earth dam sections, i.e., Chang, Tapar, Rudramata, and Kaswati, under pseudostatic loading conditions, are presented. Using the response surface methodology, in combination with first order reliability method and numerical analysis, the reliability index (beta) values are obtained and results are interpreted in conjunction with conventional factor of safety values. The influence of considering variability in the input soil shear strength parameters, horizontal seismic coefficient (alpha(h)), and location of reservoir full level on the stability assessment of the earth dam sections is discussed in the probabilistic framework. A comparison of results with those obtained from other method of reliability analysis, viz., Monte Carlo simulations combined with limit equilibrium approach, provided a basis for discussing the stability of earth dams in probabilistic terms, and the results of the analysis suggest that the considered earth dam sections are reliable and are expected to perform satisfactorily.

Total risk rating and stability analysis of embankment dams in the Kachchh Region, Gujarat, India

A method for total risk analysis of embankment dams under earthquake conditions is discussed and applied to the selected embankment dams, i.e., Chang, Tapar, Rudramata, and Kaswati located in the Kachchh region of Gujarat, India, to obtain the seismic hazard rating of the dam site and the risk rating of the structures. Based on the results of the total risk analysis of the dams, coupled non-linear dynamic numerical analyses of the dam sections are performed using acceleration time history record of the Bhuj (India) earthquake as well as five other major earthquakes recorded worldwide. The objective of doing so is to perform the numerical analysis of the dams for the range of amplitude, frequency content and time duration of input motions. The deformations calculated from the numerical analyses are also compared with other approaches available in literature, viz, Makdisi and Seed (1978) approach, Jansen's approach (1990), Swaisgood's method (1995), Bureau's method (1997). Singh et al. approach (2007), and Saygili and Rathje approach (2008) and the results are utilized to foresee the stability of dams in future earthquake scenario. (C) 2010 Elsevier B.V. All rights reserved.

Some nonstandard error analysis of discontinuous Galerkin methods for elliptic problems

An a priori error analysis of discontinuous Galerkin methods for a general elliptic problem is derived under a mild elliptic regularity assumption on the solution. This is accomplished by using some techniques from a posteriori error analysis. The model problem is assumed to satisfy a GAyenrding type inequality. Optimal order L (2) norm a priori error estimates are derived for an adjoint consistent interior penalty method.