A computational study of gas flow in a De‐Laval micronozzle at different throat diameters

A numerical study has been carried out to investigate the gas flows in a micronozzle using a continuum model under both slip and no‐slip boundary conditions. The governing equations were solved with a finite volume method. The numerical model was validated with available experimental data. Numerical results of exit thrust showed good agreement with experimental data except at very low Reynolds numbers. For parametric studies on the effect of geometric scaling, the nozzle throat diameter was varied from 10 to 0.1 mm, whereas throat Reynolds number was varied from 5 to 2000. A correlation has also been developed to calculate the specific impulse at specified throat diameter and Reynolds number. The effect of different gases on the specific impulse of the nozzle, such as helium, nitrogen, argon and carbon dioxide, was also examined. Copyright © 2008 John Wiley & Sons, Ltd.     

An accurate finite element scheme with moving meshes for computing 3D‐axisymmetric interface flows

An accurate finite element scheme for computing 3D‐axisymmetric incompressible free surface and interface flows is proposed. It is based on the arbitrary Lagrangian Eulerian (ALE) approach using free surface/interface‐resolved moving meshes. Key features like the surface force, consisting of surface tension and the local curvature, and jumps in the density and viscosity over different fluid phases are precisely incorporated in the finite element formulation. The local curvature is approximated by using the Laplace–Beltrami operator technique combined with a boundary approximation by isoparametric finite elements. A new approach is used to derive the 3D‐axisymmetric form from the variational form in 3D‐Cartesian coordinates. Several test examples show the high accuracy and the robustness of the proposed scheme. Copyright © 2007 John Wiley & Sons, Ltd.     

Nonstandard perturbation approximation and travelling wave solutions of nonlinear reaction diffusion equations

This paper deals with the construction of a nonstandard numerical method to compute the travelling wave solutions of nonlinear reaction diffusion equations at high wave speeds. Related general properties are studied using the perturbation approximation. At high wave speed the perturbation parameter approaches to zero and the problem exhibits a multiscale character. That is, there are thin layers where the solution varies rapidly, while away from these layers the solution behaves regularly and varies slowly. Most of the conventional methods fail to capture this layer behavior. Thus, the quest for some new numerical techniques that may handle the travelling wave solutions at high wave speeds earns relevance. In this paper, one such parameter robust nonstandard numerical scheme is constructed, in the sense that its numerical solution converges in the maximum norm to the exact solution uniformly well for all finite wave speeds. To overcome the difficulty due to the nonlinearity, the problem is linearized using the quasilinearization process followed by nonstandard finite difference discretization. An extensive amount of analysis is carried out which uses a suitable decomposition of the error into smooth and singular component and a comparison principle combined with appropriate barrier functions. The error estimates are obtained, which ensures uniform convergence of the method. A set of numerical experiment is carried out in support of the predicted theory that validates computationally the theoretical results. © 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2007     

The locally conservative Galerkin (LCG) method for solving the incompressible Navier–Stokes equations

In this paper, the locally conservative Galerkin (LCG) method (Numer. Heat Transfer B Fundam. 2004; 46 :357–370; Int. J. Numer. Methods Eng. 2007) has been extended to solve the incompressible Navier–Stokes equations. A new correction term is also incorporated to make the formulation to give identical results to that of the continuous Galerkin (CG) method. In addition to ensuring element‐by‐element conservation, the method also allows solution of the governing equations over individual elements, independent of the neighbouring elements. This is achieved within the CG framework by breaking the domain into elemental sub‐domains. Although this allows discontinuous trial function field, we have carried out the formulation using the continuous trial function space as the basis. Thus, the changes in the existing CFD codes are kept to a minimum. The edge fluxes, establishing the continuity between neighbouring elements, are calculated via a post‐processing step during the time‐stepping operation. Therefore, the employed formulation needs to be carried out using either a time‐stepping or an equivalent iterative scheme that allows post‐processing of fluxes. The time‐stepping algorithm employed in this paper is based on the characteristic‐based split (CBS) scheme. Both steady‐ and unsteady‐state examples presented show that the element‐by‐element formulation employed is accurate and robust. Copyright © 2007 John Wiley & Sons, Ltd.     

On a theorem of Tate

We study applications of divisibility properties of recurrence sequences to Tate’s theory of abelian varieties over finite fields.      

A comparison of boundary element and finite element methods for modeling axisymmetric polymeric drop deformation

A modified boundary element method (BEM) and the DEVSS‐G finite element method (FEM) are applied to model the deformation of a polymeric drop suspended in another fluid subjected to start‐up uniaxial extensional flow. The effects of viscoelasticity, via the Oldroyd‐B differential model, are considered for the drop phase using both FEM and BEM and for both the drop and matrix phases using FEM. Where possible, results are compared with the linear deformation theory. Consistent predictions are obtained among the BEM, FEM, and linear theory for purely Newtonian systems and between FEM and linear theory for fully viscoelastic systems. FEM and BEM predictions for viscoelastic drops in a Newtonian matrix agree very well at short times but differ at longer times, with worst agreement occurring as critical flow strength is approached. This suggests that the dominant computational advantages held by the BEM over the FEM for this and similar problems may diminish or even disappear when the issue of accuracy is appropriately considered. Fully viscoelastic problems, which are only feasible using the FEM formulation, shed new insight on the role of viscoelasticity of the matrix fluid in drop deformation. Copyright © 2001 John Wiley & Sons, Ltd.     

An assessment of a parallel,finite element method for three‐dimensional,moving‐boundary flows driven by capillarity for simulation of viscous sintering

A parallel, finite element method is presented for the computation of three‐dimensional, free‐surface flows where surface tension effects are significant. The method employs an unstructured tetrahedral mesh, a front‐tracking arbitrary Lagrangian–Eulerian formulation, and fully implicit time integration. Interior mesh motion is accomplished via pseudo‐solid mesh deformation. Surface tension effects are incorporated directly into the momentum equation boundary conditions using surface identities that circumvent the need to compute second derivatives of the surface shape, resulting in a robust representation of capillary phenomena. Sample results are shown for the viscous sintering of glassy ceramic particles. The most serious performance issue is error arising from mesh distortion when boundary motion is significant. This effect can be severe enough to stop the calculations; some simple strategies for improving performance are tested. Copyright © 2001 John Wiley & Sons, Ltd.     

蓄冰球内固定融化问题的理论与实验研究

本文针对一种自行开发的具有较高融点的新型有机混合蓄冷工质，首先建立了描述融化过程的数学模型，并应用边界元方法给出了适用于该过程的数值解，最终得到了在不同工况下相界面和热流随时间的变化规律，为进一步优化蓄冷系统提供了理论依据；同时，在小型蓄冷实验台上进行了相应工况的实验研究，验证了理论解。     

Riemann问题的间断有限元

用时空全离散间断零次有限元对Riemarm问题进行了数值求解，没有出现振荡，很好的模拟了稀疏波的逐渐稀疏化和激波的剧烈变化。