Search results

This paper introduces a novel algorithm for solving the two‐dimensional Euler and Navier‐Stokes compressible equations using a one‐step effective flux vector‐splitting implicit method. The new approach makes a contribution by deriving a simple and yet effective implicit scheme which has the features of an exact factorization and avoids the solving of block‐diagonal system of equations. This results in a significant improvement in computational efficiency as compared to the standard Beam‐Warming and Steger implicit factored schemes. The current work has advantageous characteristics in the creation of higher order numerical implicit terms. The scheme is stable if we could select the correct values of the scalars (λ±ξ and λ±η) for the respective split flux‐vectors (F± and G±) along the ξ− and η−directions. A simple solving procedure is suggested with the discussion of the implicit boundary conditions, stability analysis, time‐step length and convergence criteria. This method is spatially second‐order accurate, fully conservative and implemented with general co‐ordinate transformations for treating complex geometries. Also, the scheme shows a good convergence rate and acceptable accuracy in capturing the shock waves. Results calculated from the program developed include transonic flows through convergence‐divergence nozzle and turbine cascade. Comparisons with other well‐documented experimental data are presented and their agreements are very promising. The extension of the algorithm to 3D simulation is straightforward and under way.

The purpose of this paper is to design a simple and accurate a-posteriori Lagrangian-based error estimator is developed for the class of backward differentiation formula (BDF) algorithms with variable time step size, and the adaptive time-stepping in BDF algorithms is demonstrated for efficient time-dependent simulations in fluid flow and heat transfer.

The Lagrange interpolation polynomial is used to predict the time derivative, and then the accurate primary result is obtained by the Gauss integral, which is applied to evaluate the local error. Not only the generalized formula of the proposed error estimator is presented but also the specific expression for the widely applied BDF1/2/3 is illustrated. Two essential executable MATLAB functions to implement the proposed error estimator are appended for practical applications. Then, the adaptive time-stepping is demonstrated based on the newly proposed error estimator for BDF algorithms.

The validation tests show that the newly proposed error estimator is accurate such that the effectivity index is always close to unity for both linear and nonlinear problems, and it avoids under/overestimation of the exact local error. The applications for fluid dynamics and coupled fluid flow and heat transfer problems depict the advantage of adaptive time-stepping based on the proposed error estimator for time-dependent simulations.

In contrast to existing error estimators for BDF algorithms, the present work is more accurate for the local error estimation, and it can be readily extended to practical applications in engineering with a few changes to existing codes, contributing to efficient time-dependent simulations in fluid flow and heat transfer.

A comparison of the accuracy of the central discretization scheme with artificial dissipation and the upwind flux‐difference TVD scheme has been made for the compressible Navier‐Stokes equations for high Reynolds number flows. First, a comparison is made on two one‐dimensional model problems. Then the schemes are compared on flat plate boundary layer flow. It is shown that a central scheme basically has poor accuracy due to the isotropic nature of the artificial dissipation. An upwind scheme decomposes the flow into different components and adapts the dissipation to the velocity of the components. The associated ansitropic dissipation results in a good accuracy. It is further discussed how a central discretization scheme with artificial dissipation can be improved at the expense of the same complexity of an upwind scheme.

Isotropic artificial dissipation is added to the Navier‐Stokes equations along with a correction term which cancels the artificial dissipation term in the limit when the mesh size is zero. For a finite mesh size, the correction term replaces the artificial viscosity terms with hyperviscosity terms, i.e., with an artificial dissipation which depends on the fourth derivatives of the velocity. Hyperviscosity more effectively suppresses the higher wave number modes and has a smaller effect on the inertial modes of the flow field than does artificial viscosity. This scheme is implemented using the finite element method and therefore the required amount of dissipation is determined by analysing the discretization on a finite element. The scheme is used to simulate the flow in a driven cavity and over a backward facing step and the results are compared to existing results for these cases.

The purpose of this paper is to present two low diffusive convective-pressure flux split finite volume algorithms for solving incompressible flows in artificial compressibility framework.

The present method follows the framework similar to advection upwind splitting method of Liou and Steffen for compressible flows which is used by Vierendeels et al. to solve incompressible flow equations. Instead of discretizing the total inviscid flux using upwind scheme, the inviscid flux is first split into convective and pressure parts, and then discretized the two parts differently. The convective part is discretized using upwind method and the pressure part using central differencing. Since the Vierendeels type scheme may not be able to capture the divergence free velocity field due to the presence of artificial dissipation term, a strategy to progressively withdraw the dissipation with time step is proposed here that can ascertain the divergence free velocity condition to the level of residual error. This approach helps in reducing the amount of numerical dissipation due to upwind discretization, which is evident from the numerical test examples.

Upwind treatment of only the convective part of the inviscid flux terms, instead of the whole inviscid flux term, leads to more accurate solutions even at relatively coarse grids, which is substantiated by numerical test examples.

The method is presently applicable to Cartesian grid.

Although similar formulation is reported by Vierendeels et al., no detailed study of the accuracy is presented. Discretization and solution reconstructions used in the present approach differ from the approach reported by Vierendeels et al. A modification to Vierendeels type scheme is proposed that can help in achieving divergence free velocity condition. Finally the efficacy of the present approach to produce very accurate solutions even on coarse grids is successfully demonstrated using a few benchmark problems.

The purpose of this paper is to focus on modeling buoyancy driven viscous flow and heat transfer through saturated packed pebble‐beds via a set of homogeneous volume‐averaged conservation equations in which local thermal disequilibrium is accounted for.

The local thermal disequilibrium accounted for refers to the solid and liquid phases differing in temperature in a volume‐averaged sense, which is modeled by describing each phase with its own governing equation. The partial differential equations are discretized and solved via a vertex‐centered edge‐based dual‐mesh finite volume algorithm. A compact stencil is used for viscous terms, as this offers improved accuracy compared to the standard finite volume formulation. A locally preconditioned artificial compressibility solution strategy is employed to deal with pressure incompressibility, whilst stabilisation is achieved via a scalar‐valued artificial dissipation scheme.

The developed technology is demonstrated via the solution of natural convective flow inside a heated porous axisymmetric cavity. Predicted results were in general within 10 per cent of experimental measurements.

This is the first instance in which both artificial compressibility and artificial dissipation is employed to model flow through saturated porous materials.

The purpose of this paper is to present numerical study on the behaviour of 2D unsteady incompressible laminar wakes behind square cylinders.

The numerical method that has been developed is based on a finite point formulation characterised by its weak connectivity requirements. This formulation allows for a patched unstructured approach to computational domain modelling that is of interest for industrial applications. Time evolution of pressure is computed by using a pseudo‐compressibility relaxation model that is based on physical considerations.

This model is characterised by the fact that no sub‐iterations on a numerical pseudo‐time are required so that computational efficiency is increased. Algorithm stability requires the use of second and fourth order artificial viscosity operators that effectively change the order of the equations. A discussion is included regarding the boundary conditions for these operators that do not influence vortex shedding behaviour.

Bearing in mind the industrial drive (MEMS design) that the authors have in mind, solver validation has been addressed at two levels: global coefficients (lift, drag and Strouhal number) were compared with those published in the specialised literature, while local velocity and rms profiles were compared with those obtained after performing a specific low velocity wind tunnel testing campaign (Reynolds numbers in the range from 110 to 268).

A sensitivity analysis of the results obtained is presented and it shows that the solver numerical robustness makes it amenable for project oriented applications.

The formulation being presented is competitive and could be considered as a potential alternative to other approaches.

The purpose of this paper is to experimentally and theoretically investigate slippers, which have an important role on power dissipation in the swash plate axial piston pumps.

The slipper geometry and working conditions affected on the slipper performance have been analyzed experimentally. The model of the slipper system has been established by original neural network (NN) method.

First, the effects of the slipper geometry with smooth and conical sliding surfaces on the slipper performance were experimentally analyzed. Smooth sliding surface slippers showed a better performance then the conical surface ones. According to the results, the neural predictor would be used as a predictor for possible experimental applications on modeling this type of system.

This paper discusses a new modeling scheme known as artificial NNs an experimental and a NN approach have been employed for analyzing axial piston pumps. The simulation results suggest that the neural predictor would be used as a predictor for possible experimental applications on modeling bearing system.

This study aims to understand the difference between irreversibility in heat and work transfer processes. It also aims to explain that Helmholtz or Gibbs energy does not represent “free” energy but is a measure of loss of Carnot (reversible) work opportunity.

The entropy of mass is described as the net temperature-standardised heat transfer to mass under ideal conditions measured from a datum value. An expression for the “irreversibility” is derived in terms of work loss (W_loss) in a work transfer process, unaccounted heat dissipation (Q_loss) in a heat transfer process and loss of net Carnot work (CW_net) opportunity resulting from spontaneous heat transfer across a finite temperature difference during the process. The thermal irreversibility is attributed to not exploiting the capability for extracting work by interposing a combination of Carnot engine(s) and/or Carnot heat pump(s) that exchanges heat with the surrounding and operates across the finite temperature difference.

It is shown, with an example, how the contribution of thermal irreversibility, in estimating reversible input work, amounts to a loss of an opportunity to generate the net work output. The opportunity is created by exchanging heat with surroundings whilst transferring the same amount of heat across finite temperature difference. An entropy change is determined with a numerical simulation, including calculation of local entropy generation values, and results are compared with estimates based on an analytical expression.

A new interpretation of entropy combined with an enhanced mental image of a combination of Carnot engine(s) and/or Carnot heat pump(s) is used to quantify thermal irreversibility.

The purpose of this study is to use the Taguchi Method for parametric design in the early stages of product development. electromagnetic compatibility (EMC) issues can be considered in the early stages of product design to reduce counter-measure components, product cost and labor consumption increases due to a number of design changes in the R&D cycle and to accelerate the R&D process.

The three EMC characteristics, including radiated emission, conducted emission and fast transient impulse immunity of power, are considered response values; control factors are determined with respect to the relevant parameters for printed circuit board and mechanical design of the product and peripheral devices used in conjunction with the product are considered as noise factors. The optimal parameter set is determined by using the principal component gray relational analysis in conjunction with both response surface methodology and artificial neural network.

Market specifications and cost of components are considered to propose an optimal parameter design set with the number of grounded screw holes being 14, the size of the shell heat dissipation holes being 3 mm and the arrangement angle of shell heat dissipation holes being 45 degrees, to dispose of 390 O filters on the noise source.

The optimal parameter set can improve EMC effectively to accommodate the design specifications required by customers and pass test regulations.