H4 exponential finite-difference scheme for convective diffusion-equations

Perturbations are applied to the convective coefficients and source term of a convection-diffusion equation so that second-order corrections may be applied to a second-order exponential scheme. The basic Structure of the equations in the resulting fourth-order scheme is identical to that for the second order. Furthermore, the calculations are quite simple as the second-order corrections may be obtained in a single pass using a second-order scheme. For one to three dimensions, the fourth-order exponential scheme is unconditionally stable. As examples, the method is applied to Burgers' and other fluid mechanics problems. Compared with schemes normally used, the accuracies are found to be good and the method is applicable to regions with large gradients.

Character of simplified Navier-Stokes (SNS) equations near separation point for two-dimensional laminar flow over a flat plate

It is proved that the simplified Navier-Stokes (SNS) equations presented by Gao Zhi[1], Davis and Golowachof-Kuzbmin-Popof (GKP)[3] are respectively regular and singular near a separation point for a two-dimensional laminar flow over a flat plate. The order of the algebraic singularity of Davis and GKP equation[2,3] near the separation point is indicated. A comparison among the classical boundary layer (CBL) equations, Davis and GKP equations, Gao Zhi equations and the complete Navier-Stokes (NS) equations near the separation point is given.

Diffusion analogy and shock capturing for solving aerodynamic equations

Improving the resolution of the shock is one of the most important subjects in computational aerodynamics. In this paper the behaviour of the solutions near the shock is discussed and the reason of the oscillation production is investigated heuristically. According to the differential approximation of the difference scheme the so-called diffusion analogy equation and the diffusion analogy coefficient are defined. Four methods for improving the resolution of the shock are presented using the concept of diffusion analogy.

Some theoretical aspects of the simplified Navier-Stokes(SNS) equations

This study deals with the formulation, mathematical property and physical meaning of the simplified Navier-Stokes (SNS) equations. The tensorial SNS equations proposed is the simplest in form and is applicable to flow fields with arbitrary body boundaries. The zones of influence and dependence of the SNS equations, which are of primary importance to numerical solutions, are expounded for the first time from the viewpoint of subcharacteristics. Besides, a detailed analysis of the diffusion process in flow fields shows that the diffusion effect has an influence zone globally windward and an upwind propagation greatly depressed by convection. The maximum upwind influential distance of the viscous effect and the relative importance of the viscous effect in the flow direction to that in the direction normal to the flow are represented by the Reynolds number, which illustrates the conversion of the complete Navier-Stokes (NS) equations to the SNS equations for flows with large Reynolds number.

Basic equations for suspension flows with interphase mass-transfer

In the case of suspension flows, the rate of interphase momentum transfer M(k) and that of interphase energy transfer E(k), which were expressed as a sum of infinite discontinuities by Ishii, have been reduced to the sum of several terms which have concise physical significance. M(k) is composed of the following terms: (i) the momentum carried by the interphase mass transfer; (ii) the interphase drag force due to the relative motion between phases; (iii) the interphase force produced by the concentration gradient of the dispersed phase in a pressure field. And E(k) is composed of the following four terms, that is, the energy carried by the interphase mass transfer, the work produced by the interphase forces of the second and third parts above, and the heat transfer between phases. It is concluded from the results that (i) the term, (-alpha-k-nabla-p), which is related to the pressure gradient in the momentum equation, can be derived from the basic conservation laws without introducing the "shared-pressure presumption"; (ii) the mean velocity of the action point of the interphase drag is the mean velocity of the interface displacement, upsilonBAR-i. It is approximately equal to the mean velocity of the dispersed phase, upsilonBAR-d. Hence the work terms produced by the drag forces are f(dc) . upsilonBAR-d, and f(cd) . upsilonBAR-d, respectively, with upsilonBAR-i not being replaced by the mean velocity of the continuous phase, upsilonBAR-c; (iii) by analogy, the terms of the momentum transfer due to phase change are upsilonBAR-d-GAMMA-c, and upsilonBAR-d-GAMMA-d, respectively; (iv) since the transformation between explicit heat and latent heat occurs in the process of phase change, the algebraic sum of the heat transfer between phases is not equal to zero. Q(ic) and Q(id) are composed of the explicit heat and latent heat, so that the sum Q(ic) + Q(id)) is equal to zero.

Hierarchal structure of the simplified Navier-Stokes equations (SNSE) and its mechanical connotation and application

The hierarchial structure and mathematical property of the simplified Navier-Stokesequations (SNSE) are studied for viscous flow over a sphere and a jet of compressible flu-id. All kinds of the hierarchial SNSE can be divided into three types according to theirmathematical property and also into five groups according to their physical content. Amultilayers structure model for viscous shear flow with a main stream direction is pre-sented. For the example of viscous incompressible flow over a flat plate there existthree layers for both the separated flow and the attached flow; the character of thetransition from the three layers of attached flow to those of separated flow is elucidated.A concept of transition layer being situated between the viscous layer and inviscidlayer is introduced. The transition layer features the interaction between viscous flow andinviscid flow. The inner-outer-layers-matched SNSE proposed by the present author inthe past is developed into the layers matched (LsM)-SNSE.

Prediction of mixed-mode cracking direction in random, short-fibre composite-materials

The prediction of cracking direction in composite materials is of significance to the design of composite structures. This paper presents several methods for predicting the cracking direction in the double grooved tension-shear specimen which gives mixed-mode cracking. Five different criteria are used in this analysis: two of them have been used by other investigators and the others are proposed by the present authors. The strain energy density criterion proposed by G.C. Sih is modified to take account of the influence of the anisotropy of the strength on the direction of crack. The two failure criteria of Tsai-Hill and Norris are extended to predict the crack orientation. The stress distributions in the near-notch zone are calculated by using the 8-node quadrilateral isoparametric finite element method. The predictions of all the criteria except one are in good agreement with the experimental measurement. In addition, on the basis of the FEM results, the size of the zone in which the singular term is dominant is estimated.

Simplified Navier-Stokes equations (SNSE)

Ten kinds of the simplified Navier-Stokes equations (SNSE) are reviewed and also used to calculate the Jeffery-Hamel flow as well as to analyze briefly the seven kinds of flows to which the exact solutions of the complete Navier-Stokes equations (CNSE) have been found. Analysis shows that the actual differences among the solutions of the different SNSE can go beyond the range of the order of magnitude of Re-1/2 and result even in different flow patterns, therefore, how to choose the viscous terms included in the SNSE is worthy of notice where Re=S∞u∞ L/μ∞ is the Reynolds numbers. For the aforesaid eight kinds of flows, the solutions to the inner-outer-layer-matched SNSE and to the thin-layer-2-order SNSE agree completely with the exact solutions to CNSE. But the solutions to all the other SNSE are not completely consistent with the exact solutions to CNSE and not a few of them are actually the solutions of the classical boundary layer theory. The innerouter-layer-matched SNSE contains the shear stress causing angular displacement of the inormal axis with respect to the streamwise axis and the normal stress causing expansion-contraction in the direction of the normal axis and the viscous terms being of the order of magnitude of the normal stress; and it can also reasonably treat the inertial terms as well as the relation between the viscous and inertial terms. Therefore, it seems promising in respects of both mechanics and mathematics.

An asymptotic solution of hydrodynamic equations at the mid-plane of a plume in the earths mantle

From the partial differential equations of hydrodynamics governing the movements in the Earth's mantle of a Newtonian fluid with a pressure- and temperature-dependent viscosity, considering the bilateral symmetry of velocity and temperature distributions at the mid-plane of the plume, an analytical solution of the governing equations near the mid-plane of the plume was found by the method of asymptotic analysis. The vertical distribution of the upward velocity, viscosity and temperature at the mid-plane, and the temperature excess at the centre of the plume above the ambient mantle temperature were then calculated for two sets of Newtonian rheological parameters. The results obtained show that the temperature at the mid-plane and the temperature excess are nearly independent of the rheological parameters. The upward velocity at the mid-plane, however, is strongly dependent on the rheological parameters.

Perturbational finite volume scheme for the one-dimensional Navier-Stokes equations

Starting from the second-order finite volume scheme,though numerical value perturbation of the cell facial fluxes, the perturbational finite volume (PFV) scheme of the Navier-Stokes (NS) equations for compressible flow is developed in this paper. The central PFV scheme is used to compute the one-dimensional NS equations with shock wave.Numerical results show that the PFV scheme can obtain essentially non-oscillatory solution.

Basic equations for ferroelectrics

A set of new formula of energy functions for ferroelectrics was proposed, and then the new basic equations were derived in this paper. The finite element formulation based on the new basic equations was improved to avoid the equivalent nodal load produced by remnant polarization. With regard to the fundamentals of mathematics and physics, the new energy functions and basic equations are reasonable for the material element of ferroelectrics in finite element analysis.

Prediction of particle distribution in isotropic turbulence by large-eddy simulation

The recent application of large-eddy simulation (LES) to particle-laden turbulence requires that the LES with a subgrid scale (SGS) model could accurately predict particle distributions. Usually, a SGS particle model is used to recover the small-scale structures of velocity fields. In this study, we propose a rescaling technique to recover the effects of small-scale motions on the preferential concentration of inertial particles. The technique is used to simulate particle distribution in isotropic turbulence by LES and produce consistent results with direct numerical simulation (DNS). Key words: particle distribution, particle-laden turbulence, large-eddy simulation, subgrid scale model.

Energy Functions And Basic Equations For Ferroelectrics

Energy functions (or characteristic functions) and basic equations for ferroelectrics in use today are given by those for ordinary dielectrics in the physical and mechanical communications. Based on these basic equations and energy functions, the finite element computation of the nonlinear behavior of the ferroelectrics has been carried out by several research groups. However, it is difficult to process the finite element computation further after domain switching, and the computation results are remarkably deviating from the experimental results. For the crack problem, the iterative solution of the finite element calculation could not converge and the solutions for fields near the crack tip oscillate. In order to finish the calculation smoothly, the finite element formulation should be modified to neglect the equivalent nodal load produced by spontaneous polarization gradient. Meanwhile, certain energy functions for ferroelectrics in use today are not compatible with the constitutive equations of ferroelectrics and need to be modified. This paper proposes a set of new formulae of the energy functions for ferroelectrics. With regard to the new formulae of the energy functions, the new basic equations for ferroelectrics are derived and can reasonably explain the question in the current finite element analysis for ferroelectrics.

Prediction Of Shear-Band Thickness In Metallic Glasses

We derive an explicit expression for predicting the thicknesses of shear bands in metallic glasses. The model demonstrates that the shear-band thickness is mainly dominated by the activation size of the shear transformation zone (STZ) and its activation free volume concentration. The predicted thicknesses agree well with the results of measurements and simulations. The underlying physics is attributed to the local topological instability of the activated STZ. The result is of significance in understanding the origin of inhomogeneous flow in metallic glasses. (C) 2009 Acta Materialia Inc. Published by Elsevier Ltd. All rights reserved.