Porous versus porthole fuel injection in a radical farming scramjet: A numerical analysis

Numerically computed engine performance of a nominally two-dimensional radical farming scramjet with porous (permeable C/C ceramic) and porthole fuel injection is presented. Inflow conditions with Mach number, stagnation pressure, and enthalpy of 6.44, 40.2MPa, and 4.31 MJ/kg respectively, and fuel/air equivalence ratio of 0.44 were maintained, along with engine geometry. Hydrogen fuel was injected at an axial location of 92.33mm downstream of the leading edge for each investigated injection method. Results from this study show that porous fuel injection results in enhanced mixing and combustion compared to porthole fuel injection. This is particularly evident within the first half of the combustion chamber where porous fuel injection resulted in mixing and combustion efficiencies of 76% and 63% respectively. At the same location, porthole fuel injection resulted in efficiencies respectively of 58% and 46%. Key mechanisms contributing to the observed improved performance were the formation of an attached oblique fuel injection shock and associated stronger shock-expansion train ingested by the engine, enhanced spreading of the fuel in all directions and a more rapidly growing mixing layer.

Numerical analysis for the time distributed-order and Riesz space fractional diffusions on bounded domains

Subdiffusion equations with distributed-order fractional derivatives describe some important physical phenomena. In this paper, we consider the time distributed-order and Riesz space fractional diffusions on bounded domains with Dirichlet boundary conditions. Here, the time derivative is defined as the distributed-order fractional derivative in the Caputo sense, and the space derivative is defined as the Riesz fractional derivative. First, we discretize the integral term in the time distributed-order and Riesz space fractional diffusions using numerical approximation. Then the given equation can be written as a multi-term time–space fractional diffusion. Secondly, we propose an implicit difference method for the multi-term time–space fractional diffusion. Thirdly, using mathematical induction, we prove the implicit difference method is unconditionally stable and convergent. Also, the solvability for our method is discussed. Finally, two numerical examples are given to show that the numerical results are in good agreement with our theoretical analysis.

Numerical analysis for the interlayer toughening of piezoelectric smart composite plate

The piezoelectric composite material could engender stress concentration resulting from small cracks during layers easily, as the cracks growth will lead to the failure of the whole structure. In this paper, a finite element model for piezoelectric composite materials by ABAQUS including interlayer crack was established, and the J integral and crack tip stress of different types PZT patches were calculated by using the equivalent integral method. Then, the J integral for adhesive layers with different thickness, elastic modulus considering and not considering piezoelectricity was investigated. The results show that the J integral of mode I, II reduces with thicker adhesive layer and lower elastic modules, and the J integral of mode II decreases more sharply than that of mode I.

Multigrid algorithms for symmetric discontinuous Galerkin methods on graded meshes

We study a class of symmetric discontinuous Galerkin methods on graded meshes. Optimal order error estimates are derived in both the energy norm and the L 2 norm, and we establish the uniform convergence of V-cycle, F-cycle and W-cycle multigrid algorithms for the resulting discrete problems. Numerical results that confirm the theoretical results are also presented.

Meshless Local Petrov-Galerkin Method for Rotating Euler-Bernoulli Beam

Free vibration problem of a rotating Euler-Bernoulli beam is solved with a truly meshless local Petrov-Galerkin method. Radial basis function and summation of two radial basis functions are used for interpolation. Radial basis function satisfies the Kronecker delta property and makes it simpler to apply the essential boundary conditions. Interpolation with summation of two radial basis functions increases the node carrying capacity within the sub-domain of the trial function and higher natural frequencies can be computed by selecting the complete domain as a sub-domain of the trial function. The mass and stiffness matrices are derived and numerical results for frequencies are obtained for a fixed-free beam and hinged-free beam simulating hingeless and articulated helicopter blades. Stiffness and mass distribution suitable for wind turbine blades are also considered. Results show an accurate match with existing literature.

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.

Numerical Analysis of Rainfall Infiltration in the Slope with a Fracture

With the finite volume method, a 2D numerical model for seepage in unsaturated soil has been established to study the rainfall infiltration in the fractured slope.The result shows that more rain may infiltrate into the slope due to existing fracture and then the pore pressure rises correspondingly. Very probably, it is one of the crucial factors accounting for slope failure. Furthermore a preliminary study has been conducted to investigate the influence of various fracture and rainfall factors such as the depth, width and location of a crack, surface condition, rainfall intensity and duration. Pore pressure and water volumetric content during the transient seepage are carefully examined to reveal the intrinsic mechanism.

Numerical Analysis of Theoretical Model of the Rf MEMs Switches

An improved electromechanical model of the RF MEMS (radio frequency microelectromechanical systems) switches is introduced, in which the effects of intrinsic residual stress from fabrication processes, axial stress due to stretching of beam, and fringing field are taken into account. Four dimensionless numbers are derived from the governing equation of the developed model. A semi-analytical method is developed to calculate the behavior of the RF MEMS switches. Subsequently the influence of the material and geometry parameters on the behavior of the structure is analyzed and compared, and the corresponding analysis with the dimensionless numbers is conducted too. The quantitative relationship between the presented parameters and the critical pull-in voltage is obtained, and the relative importance of those parameters is given.

Numerical Analysis of LEC Growth of GaAs with an Axial Magnetic Field

A set of numerical analyses for momentum and heat transfer For a 3 in. (0.075 m) diameter Liquid Encapsulant Czochralski (LEC) growth of single-crystal GaAs with or without all axial magnetic field was carried Out using the finite-element method. The analyses assume a pseudosteady axisymmetric state with laminar floats. Convective and conductive heat transfers. radiative heat transfer between diffuse surfaces and the Navier-Stokes equations for both melt and encapsulant and electric current stream function equations Cor melt and crystal Lire considered together and solved simultaneously. The effect of the thickness of encapsulant. the imposed magnetic field strength as well as the rotation rate of crystal and crucible on the flow and heat transfer were investigated. (C) 2002 Published by Elsevier Science Ltd.

Numerical analysis and comparison of three metal-oxide-type super-resolution near field structures

Based on the Fresnel-Kirchkoff diffraction theory, we build up a Gaussian diffraction model of metal-oxide-type super-resolution near field structure (super-RENS), which can describe far field optical properties. The spectral contrast induced by refractive index and the structural changes in AgOx, PtOx and PdOx thin films, which are the key functional layers in super-RENS, are studied by using this model. Comparison results indicate that the spectral contrast intensively on laser-induced distribution and change of the refractive index in the metal-oxide films. The readout mechanism of the metal-oxide-type super-RENS optical disc is further clarified. This Gaussian diffraction model can be used as a simple and effective method for choosing proper active materials in super-RENS.

Numerical analysis of the coupled circuit and cooling holes for an electromagnetic shaker

This paper presents a time-stepping shaker modeling scheme. The new method improves the accuracy of analysis of armature-position-dependent inductances and force factors, analysis of axial variation of current density in copper plates (short-circuited turns), and analysis of cooling holes in the magnetic circuit. Linear movement modeling allows armature position to be precisely included in the shaker analysis. A more accurate calculation of eddy currents in the coupled circuit is in particular crucial for the shaker analysis in a mid-or high-frequency operation range. Large currents in a shaker, including eddy currents, incur large Joule losses, which in turn require the use of a cooling system to keep temperature at bay. Sizable cooling holes have influence on the saturation state of iron poles, and hence have to be properly taken into account.

The post-processed Galerkin method applied to non-linear shell vibrations

We present the results of a computational study of the post-processed Galerkin methods put forward by Garcia-Archilla et al. applied to the non-linear von Karman equations governing the dynamic response of a thin cylindrical panel periodically forced by a transverse point load. We spatially discretize the shell using finite differences to produce a large system of ordinary differential equations (ODEs). By analogy with spectral non-linear Galerkin methods we split this large system into a 'slowly' contracting subsystem and a 'quickly' contracting subsystem. We then compare the accuracy and efficiency of (i) ignoring the dynamics of the 'quick' system (analogous to a traditional spectral Galerkin truncation and sometimes referred to as 'subspace dynamics' in the finite element community when applied to numerical eigenvectors), (ii) slaving the dynamics of the quick system to the slow system during numerical integration (analogous to a non-linear Galerkin method), and (iii) ignoring the influence of the dynamics of the quick system on the evolution of the slow system until we require some output, when we 'lift' the variables from the slow system to the quick using the same slaving rule as in (ii). This corresponds to the post-processing of Garcia-Archilla et al. We find that method (iii) produces essentially the same accuracy as method (ii) but requires only the computational power of method (i) and is thus more efficient than either. In contrast with spectral methods, this type of finite-difference technique can be applied to irregularly shaped domains. We feel that post-processing of this form is a valuable method that can be implemented in computational schemes for a wide variety of partial differential equations (PDEs) of practical importance.

3D numerical analysis of tunnelling induced damage: the influence of the alignment of a masonry building with the tunnel axis.

The development of infrastructure in major cities often involves tunnelling, which can cause damage to existing structures. Therefore, these projects require a careful prediction of the risk of settlement induced damage. The simplified approach of current methods cannot account for three-dimensional structural aspects of buildings, which can result in an inaccurate evaluation of damage. This paper investigates the effect of the building alignment with the tunnel axis on structural damage. A three-dimensional, phased, fully coupled finite element model with non-linear material properties is used as a tool to perform a parametric study. The model includes the simulation of the tunnel construction process, with the tunnel located adjacent to a masonry building. Three different type of settlements are included (sagging, hogging and a combination of them), with seven different increasing angles of the building with respect to the tunnel axis. The alignment parameter is assessed, based on the maximum occurring crack width, measured in the building. Results show a significant dependency of the final damage on the building and tunnel alignment.

Finite element beam propagation method for analysis of plasmonic waveguide

The basic idea of the finite element beam propagation method (FE-BPM) is described. It is applied to calculate the fundamental mode of a channel plasmonic polariton (CPP) waveguide to confirm its validity. Both the field distribution and the effective index of the, fundamental mode are given by the method. The convergence speed shows the advantage and stability of this method. Then a plasmonic waveguide with a dielectric strip deposited on a metal substrate is investigated, and the group velocity is negative for the fundamental mode of this kind of waveguide. The numerical result shows that the power flow direction is reverse to that of phase velocity.