Forebody Compressibility Research of Hypersonic Vehicle

Three kinds of forebody model of hypersonic vehicles were studied with numerical simulation method. It shows that the two- order compressive ramp model is the best selection among the three for its good evaluative parameters value at the cowl of the inlet . This model can provide higher value of flux coefficient and total pressure recovery coefficient and lower average Mach number compared with those of the other two models . Simultaneously different compressive angles may have different effects . The configuration which the firstorder of compressive angle is 4°and the second 5°is the optimum combination. Furthermore factors such as attack angle were concerned. Better result may be obtained with a range of attack angles . Based on the work above the integrated design for forebodyPinlet of a hypersonic vehicle was performed. The numerical result shows that this integrated model provides good flow field quality for inlet and engine work.

超声速大攻角旋成体迎风激波数值模拟

旋成弹体在大攻角超声速流动中其迎风面上有时会出现一种迎风激波.采用非定常雷诺平均的Navier-Stokes方程以及两方程kω-sst和代数B-L湍流模型成功地模拟了该现象,验证了该迎风激波是由于流场主旋涡导致超声速流动偏移造成的涡诱导激波的一部分的结论.

Elastic fields of interacting elliptic inhomogeneities

This paper presents a method for the calculation of two-dimensional elastic fields in a solid containing any number of inhomogeneities under arbitrary far field loadings. The method called 'pseudo-dislocations method', is illustrated for the solution of interacting elliptic inhomogeneities. It reduces the interacting inhomogeneities problem to a set of linear algebraic equations. Numerical results are presented for a variety of elliptic inhomogeneity arrangements, including the special cases of elliptic holes, cracks and circular inhomogeneities. All these complicated problems can be solved with high accuracy and efficiency.

The Effect of Conical Flowfields on the Performance of Waveriders at Mach 6

The performance of 23 kinds of waveriders, derived from different conical flowfields, is analyzed by the numerical computation under the conditions of fight speed of Mach 6, attack angle of 0° and flight altitude of 30 km. These results indicate that the performance is influenced by the shapes and the width to height ratios (W/H ) of generating cones. The geometrical parameter and the lift coefficient are proportional to W/H, while the drag coefficient and the lift to drag ratio (L/D ) have extreme values. Considering the base drag and the computation errors, the waverider with the highest L/D is cut from the elliptical cone’s flowfield (W/H = 1.5―1.618), and the configuration with the lowest drag can also be obtained at W/H = 1:1.5. Accordingly, good suggestions are proposed for practical design based on these computational results.

Modeling of Ammonothermal Growth of Nitrides

Bulk single crystals of GaN and AlN can be grown from supercritical fluids using the ammonothermal method, which utilizes ammonia as fluid rather than water as in the hydrothermal process. In this process, a mineralizer such as amide, imide or nitride is used to attack a bulk nitride feedstock at temperatures from 200°C to 500°C and pressures from 1 to 4 kbar. Ammonothermal systems have been modeled here using fluid dynamics, thermodynamics and heat transfer models. The nutrient is considered as a porous media bed and the fluid flow is simulated using the Darcy-Brinkman-Forchheimer model. The resulting governing equations are solved using the finite volume method. The effects of particle size on flow pattern and temperature distribution in an autoclave are analyzed.

Turbulence modeling for flows around convex features giving rapid eddy distortion

Reynolds averaged Navier-Stokes model performances in the stagnation and wake regions for turbulent flows with relatively large Lagrangian length scales (generally larger than the scale of geometrical features) approaching small cylinders (both square and circular) is explored. The effective cylinder (or wire) diameter based Reynolds number, ReW ≤ 2.5 × 103. The following turbulence models are considered: a mixing-length; standard Spalart and Allmaras (SA) and streamline curvature (and rotation) corrected SA (SARC); Secundov's νt-92; Secundov et al.'s two equation νt-L; Wolfshtein's k-l model; the Explicit Algebraic Stress Model (EASM) of Abid et al.; the cubic model of Craft et al.; various linear k-ε models including those with wall distance based damping functions; Menter SST, k-ω and Spalding's LVEL model. The use of differential equation distance functions (Poisson and Hamilton-Jacobi equation based) for palliative turbulence modeling purposes is explored. The performance of SA with these distance functions is also considered in the sharp convex geometry region of an airfoil trailing edge. For the cylinder, with ReW ≈ 2.5 × 103 the mixing length and k-l models give strong turbulence production in the wake region. However, in agreement with eddy viscosity estimates, the LVEL and Secundov νt-92 models show relatively little cylinder influence on turbulence. On the other hand, two equation models (as does the one equation SA) suggest the cylinder gives a strong turbulence deficit in the wake region. Also, for SA, an order or magnitude cylinder diameter decrease from ReW = 2500 to 250 surprisingly strengthens the cylinder's disruptive influence. Importantly, results for ReW ≪ 250 are virtually identical to those for ReW = 250 i.e. no matter how small the cylinder/wire its influence does not, as it should, vanish. Similar tests for the Launder-Sharma k-ε, Menter SST and k-ω show, in accordance with physical reality, the cylinder's influence diminishing albeit slowly with size. Results suggest distance functions palliate the SA model's erroneous trait and improve its predictive performance in wire wake regions. Also, results suggest that, along the stagnation line, such functions improve the SA, mixing length, k-l and LVEL results. For the airfoil, with SA, the larger Poisson distance function increases the wake region turbulence levels by just under 5%. © 2007 Elsevier Inc. All rights reserved.

Continuous spectrum growth and modal instability in swirling duct flow

We present results on the stability of compressible inviscid swirling flows in an annular duct. Such flows are present in aeroengines, for example in the by-pass duct, and there are also similar flows in many aeroacoustic or aeronautical applications. The linearised Euler equations have a ('critical layer') singularity associated with pure convection of the unsteady disturbance by the mean flow, and we focus our attention on this region of the spectrum. By considering the critical layer singularity, we identify the continuous spectrum of the problem and describe how it contributes to the unsteady field. We find a very generic family of instability modes near to the continuous spectrum, whose eigenvalue wavenumbers form an infinite set and accumulate to a point in the complex plane. We study this accumulation process asymptotically, and find conditions on the flow to support such instabilities. It is also found that the continuous spectrum can cause a new type of instability, leading to algebraic growth with an exponent determined by the mean flow, given in the analysis. The exponent of algebraic growth can be arbitrarily large. Numerical demonstrations of the continuous spectrum instability, and also the modal instabilities are presented.

Simulation of Coal Combustion by AUSM Turbulence-Chemistry Char Combustion Model and a Full Two-Fluid Model

An algebraic unified second-order moment (AUSM) turbulence-chemistry model of char combustion is introduced in this paper, to calculate the effect of particle temperature fluctuation on char combustion. The AUSM model is used to simulate gas-particle flows, in coal combustion in a pulverized coal combustor, together with a full two-fluid model for reacting gas-particle flows and coal combustion, including the sub-models as the k-epsilon-k(p) two-phase turbulence niodel, the EBU-Arrhenius volatile and CO combustion model, and the six-flux radiation model. A new method for calculating particle mass flow rate is also used in this model to correct particle outflow rate and mass flow rate for inside sections, which can obey the principle of mass conservation for the particle phase and can also speed up the iterating convergence of the computation procedure effectively. The simulation results indicate that, the AUSM char combustion model is more preferable to the old char combustion model, since the later totally eliminate the influence of particle temperature fluctuation on char combustion rate.

Weaknesses of Margulis and Ramanujan-Margulis low-density parity-check codes

We report weaknesses in two algebraic constructions of low-density parity-check codes based on expander graphs. The Margulis construction gives a code with near-codewords, which cause problems for the sum-product decoder; The Ramanujan-Margulis construction gives a code with low-weight codewords, which produce an error-floor. © 2004 Elsevier B.V.

Weaknesses of margulis and ramanujan-margulis low-density parity-check codes

We report weaknesses in two algebraic constructions of low-density parity-check codes based on expander graphs. The Margulis construction gives a code with near-codewords, which cause problems for the sum-product decoder; The Ramanujan-Margulis construction gives a code with low-weight codewords, which produce an error-floor. ©2003 Published by Elsevier Science B. V.

Control of Mean and Fluctuating Forces on a Circular Cylinder at High Reynolds Numbers

A narrow strip is used to control mean and fluctuating forces on a circular cylinder at Reynolds numbers from 2.0 x 10(4) to 1.0 x 10(5). The axes of the strip and cylinder are parallel. The control parameters are strip width ratio and strip position characterized by angle of attack and distance from the cylinder. Wind tunnel tests show that the vortex shedding from both sides of the cylinder can be suppressed, and mean drag and fluctuating lift on the cylinder can be reduced if the strip is installed in an effective zone downstream of the cylinder. A phenomenon of mono-side vortex shedding is found. The strip-induced local changes of velocity profiles in the near wake of the cylinder are measured, and the relation between base suction and peak value in the power spectrum of fluctuating lift is studied. The control mechanism is then discussed from different points of view.

Numerical Simulation and Experimental Investigation about Internal and External Flows

In this paper, TASCflow3D is used to solve inner and outer 3D viscous incompressible turbulent flow (R-e = 5.6 X 10(6)) around axisymmetric body with duct. The governing equation is a RANS equation with standard k-epsilon turbulence model. The discrete method used is a finite volume method based on the finite element approach. In this method, the description of geometry is very flexible and at the same time important conservative properties are retained. The multi-block and algebraic multi-grid techniques are used for the convergence acceleration. Agreement between experimental results and calculation is good. It indicates that this novel approach can be used to simulate complex flow such as the interaction between rotor and stator or propulsion systems containing tip clearance and cavitation.

Studies of the Effect of a Coal Concentrator on NO Formation in Swirling Coal Combustion

To develop low-pollution burners, the effect of a coal concentrator on NO formation in swirling coal combustion is studied using both numerical simulation and experiments. The isothermal gas-particle two-phase velocities and particle concentration in a cold model of swirl burners with and without coal concentrators were measured using the phase Doppler particle anemometer (PDPA). A full two-fluid model of reacting gas-particle flows and coal combustion with an algebraic unified second-order moment (AUSM) turbulence-chemistry model for the turbulent reaction rate of NO formation are used to simulate swirling coal combustion and NO formation with different coal concentrators. The results give the turbulent kinetic energy, particle concentration, temperature and NO concentration in cases of with and without coal concentrators. The predicted results for cold two-phase flows are in good agreement with the PDPA measurement results, showing that the coal concentrator increases the turbulence and particle concentration in the recirculation zone. The combustion modeling results indicate that although the coal concentrator increases the turbulence and combustion temperature, but still can remarkably reduce the NO formation due to creating high coal concentration in the recirculation zone.

Equivalent Model Of Expansion Of Cement Mortar Under Sulphate Erosion

The expansion property of cement mortar under the attack of sulfate ions is studied by experimental and theoretical methods. First, cement mortars are fabricated with the ratio of water to cement of 0.4, 0.6, and 0.8. Secondly, the expansion of specimen immerged in sulphate solution is measured at different times. Thirdly, a theoretical model of expansion of cement mortar under sulphate erosion is suggested by virtue of represent volume element method. In this model, the damage evolution due to the interaction between delayed ettringite and cement mortar is taken into account. Finally, the numerical calculation is performed. The numerical and experimental results indicate that the model perfectly describes the expansion of the cement mortar.

Direct Numerical Simulation Of Hypersonic Boundary-Layer Transition Over A Blunt Cone

Direct numerical simulation of transition How over a blunt cone with a freestream Mach number of 6, Reynolds number of 10,000 based on the nose radius, and a 1-deg angle of attack is performed by using a seventh-order weighted essentially nonoscillatory scheme for the convection terms of the Navier-Stokes equations, together with an eighth-order central finite difference scheme for the viscous terms. The wall blow-and-suction perturbations, including random perturbation and multifrequency perturbation, are used to trigger the transition. The maximum amplitude of the wall-normal velocity disturbance is set to 1% of the freestream velocity. The obtained transition locations on the cone surface agree well with each other far both cases. Transition onset is located at about 500 times the nose radius in the leeward section and 750 times the nose radius in the windward section. The frequency spectrum of velocity and pressure fluctuations at different streamwise locations are analyzed and compared with the linear stability theory. The second-mode disturbance wave is deemed to be the dominating disturbance because the growth rate of the second mode is much higher than the first mode. The reason why transition in the leeward section occurs earlier than that in the windward section is analyzed. It is not because of higher local growth rate of disturbance waves in the leeward section, but because the growth start location of the dominating second-mode wave in the leeward section is much earlier than that in the windward section.