Search results

Extrapolation is a process used to accelerate the convergence of a sequence of approximations to the true value. Different stepsizes are used to obtain approximate solutions, which are combined to increase the order of the approximation by eliminating leading error terms. The smoothing technique is also applied to suppress order reduction and to dampen the oscillatory component in the numerical solution when solving stiff problems. The extrapolation and smoothing technique can be applied in either active, passive or the combination of both active and passive modes. In this paper, the authors investigate the best strategy of implementing extrapolation and smoothing technique and use this strategy to solve stiff ordinary differential equations. Based on the experiment, the authors suggest using passive smoothing in order to reduce the computation time.

The two-step smoothing is a composition of four steps of the symmetric method with different weights. It is used as the final two steps when combined with many steps of the symmetric method. The aim is to preserve symmetry and provide damping for stiff problem and to be more robust than the one-step smoothing. The two-step smoothing is L-stable. The new method is then applied with extrapolation process in passive and active modes to investigate the most efficient and accurate method of implementation.

In this paper, the authors constructed the two-step smoothing to be more robust than the one-step smoothing. The two-step smoothing is constructed to achieve as high order as possible and able to restore the classical order of particular method compared to the one-step active smoothing that is only able to achieve order-1 condition. The two-step smoothing for ITR is also superior in solving stiff case since it has the super-convergent order-4 behavior. In our experiments with extrapolation, it is proven that the two-step smoothing is more accurate and more efficient than the one-step smoothing, namely 1ASAX. It is also observed that the method with smoothing is comparable if not superior to the existing base method in certain cases. Based on the experiment, the authors would suggest using passive smoothing if the aim is to reduce computation time. It is of interest to conduct more experiment to validate the accuracy and efficiency of the smoothing formula with and without extrapolation.

The implementation of extrapolation on two-step symmetric Runge–Kutta method has not been tested on variety of other test problems yet. The two-step symmetrization is an extension of the one-step symmetrization and has not been constructed by other researchers yet. The method is constructed such that it preserves the asymptotic error expansion in even powers of stepsize, and when used with extrapolation the order might increase by 2 at a time. The method is also L-stable and eliminates the order reduction phenomenon when solving stiff ODEs. It is also of interest to observe other ways of implementing extrapolation using other sequences or with interpolation.

The purpose of this article is to derive an implicit-explicit local differential transform method (IELDTM) in dealing with the spatial approximation of the stiff advection-diffusion-reaction (ADR) equations.

A direction-free numerical approach based on local Taylor series representations is designed for the ADR equations. The differential equations are directly used for determining the local Taylor coefficients and the required degrees of freedom is minimized. The complete system of algebraic equations is constructed with explicit/implicit continuity relations with respect to direction parameter. Time integration of the ADR equations is continuously utilized with the Chebyshev spectral collocation method.

The IELDTM is proven to be a robust, high order, stability preserved and versatile numerical technique for spatial discretization of the stiff partial differential equations (PDEs). It is here theoretically and numerically shown that the order refinement (p-refinement) procedure of the IELDTM does not affect the degrees of freedom, and thus the IELDTM is an optimum numerical method. A priori error analysis of the proposed algorithm is done, and the order conditions are determined with respect to the direction parameter.

The IELDTM overcomes the known disadvantages of the differential transform-based methods by providing reliable convergence properties. The IELDTM is not only improving the existing Taylor series-based formulations but also provides several advantages over the finite element method (FEM) and finite difference method (FDM). The IELDTM offers better accuracy, even when using far less degrees of freedom, than the FEM and FDM. It is proven that the IELDTM produces solutions for the advection-dominated cases with the optimum degrees of freedom without producing an undesirable oscillation.

Partially‐decoupled upwind‐based total‐variation‐diminishing (TVD) finite‐difference schemes for the solution of the conservation laws governing two‐dimensional non‐equilibrium vibrationally relaxing and chemically reacting flows of thermally‐perfect gaseous mixtures are presented. In these methods, a novel partially‐decoupled flux‐difference splitting approach is adopted. The fluid conservation laws and species concentration and vibrational energy equations are decoupled by means of a frozen flow approximation. The resulting partially‐decoupled gas‐dynamic and thermodynamic subsystems are then solved alternately in a lagged manner within a time marching procedure, thereby providing explicit coupling between the two equation sets. Both time‐split semi‐implicit and factored implicit flux‐limited TVD upwind schemes are described. The semi‐implicit formulation is more appropriate for unsteady applications whereas the factored implicit form is useful for obtaining steady‐state solutions. Extensions of Roe's approximate Riemann solvers, giving the eigenvalues and eigenvectors of the fully coupled systems, are used to evaluate the numerical flux functions. Additional modifications to the Riemann solutions are also described which ensure that the approximate solutions are not aphysical. The proposed partially‐decoupled methods are shown to have several computational advantages over chemistry‐split and fully coupled techniques. Furthermore, numerical results for single, complex, and double Mach reflection flows, as well as corner‐expansion and blunt‐body flows, using a five‐species four‐temperature model for air demonstrate the capabilities of the methods.

This paper aims to develop a meshfree algorithm based on local radial basis functions (RBFs) combined with the differential quadrature (DQ) method to provide numerical approximations of the solutions of time-dependent, nonlinear and spatially one-dimensional reaction-diffusion systems and to capture their evolving patterns. The combination of local RBFs and the DQ method is applied to discretize the system in space; implicit multistep methods are subsequently used to discretize in time.

In a method of lines setting, a meshless method for their discretization in space is proposed. This discretization is based on a DQ approach, and RBFs are used as test functions. A local approach is followed where only selected RBFs feature in the computation of a particular DQ weight.

The proposed method is applied on four reaction-diffusion models: Huxley’s equation, a linear reaction-diffusion system, the Gray–Scott model and the two-dimensional Brusselator model. The method captured the various patterns of the models similar to available in literature. The method shows second order of convergence in space variables and works reliably and efficiently for the problems.

The originality lies in the following facts: A meshless method is proposed for reaction-diffusion models based on local RBFs; the proposed scheme is able to capture patterns of the models for big time T; the scheme has second order of convergence in both time and space variables and Nuemann boundary conditions are easy to implement in this scheme.

For a partial differential equation with a fourth-order derivative such as the Cahn-Hilliard equation, it is always a challenge to design numerical schemes that can handle the restrictive time step introduced by this higher order term. The purpose of this paper is to employ a fractional splitting method to isolate the convective, the nonlinear second-order and the fourth-order differential terms.

The full equation is then solved by consistent schemes for each differential term independently. In addition to validating the second-order accuracy, the authors will demonstrate the efficiency of the proposed method by validating the dissipation of the Ginzberg-Lindau energy and the coarsening properties of the solution.

The scheme is second-order accuracy, the authors will demonstrate the efficiency of the proposed method by validating the dissipation of the Ginzberg-Lindau energy and the coarsening properties of the solution.

The authors believe that this is the first time the equation is handled numerically using the fractional step method. Apart from the fact that the fractional step method substantially reduces computational time, it has the advantage of simplifying a complex process efficiently. This method permits the treatment of each segment of the original equation separately and piece them together, in a way that will be explained shortly, without destroying the properties of the equation.

Due to the strong reliance on element quality, there exist some inherent shortcomings of the traditional finite element method (FEM). The model of FEM behaves overly stiff, and the solutions of automated generated linear elements are generally of poor accuracy about especially gradient results. The proposed cell-based smoothed point interpolation method (CS-PIM) aims to improve the results accuracy of the thermoelastic problems via properly softening the overly-stiff stiffness.

This novel approach is based on the newly developed G space and weakened weak (w²) formulation, and of which shape functions are created using the point interpolation method and the cell-based gradient smoothing operation is conducted based on the linear triangular background cells.

Owing to the property of softened stiffness, the present method can generally achieve better accuracy and higher convergence results (especially for the temperature gradient and thermal stress solutions) than the FEM does by using the simplest linear triangular background cells, which has been examined by extensive numerical studies.

The CS-PIM is capable of producing more accurate results of temperature gradients as well as thermal stresses with the automated generated and unstructured background cells, which make it a better candidate for solving practical thermoelastic problems.

It is the first time that the novel CS-PIM was further developed for solving thermoelastic problems, which shows its tremendous potential for practical implications.

The purpose of this study is to present and demonstrate a numerical method for solving chemically reacting flows. These are important for energy conversion devices, which rely on chemical reactions as their operational mechanism, with heat generated from the combustion of the fuel, often gases, being converted to work.

The numerical study of such flows requires the set of Navier-Stokes equations to be extended to include multiple species and the chemical reactions between them. The numerical method implemented in this study also accounts for changes in the material properties because of temperature variations and the process to handle steep spatial fronts and stiff source terms without incurring any numerical instabilities. An all-speed numerical framework is used through simple low-dissipation advection upwind splitting (SLAU) convective scheme, and it has been extended in a multi-component species framework on the in-house density-based flow solver. The capability of solving turbulent combustion is also implemented using the Eddy Dissipation Concept (EDC) framework and the recent k-kl turbulence model.

The numerical implementation has been demonstrated for several stiff problems in laminar and turbulent combustion. The laminar combustion results are compared from the corresponding results from the Cantera library, and the turbulent combustion computations are found to be consistent with the experimental results.

This paper has extended the single gas density-based framework to handle multi-component gaseous mixtures. This paper has demonstrated the capability of the numerical framework for solving non-reacting/reacting laminar and turbulent flow problems. The all-speed SLAU convective scheme has been extended in the multi-component species framework, and the turbulent model k-kl is used for turbulent combustion, which has not been done previously. While the former method provides the capability of solving for low-speed flows using the density-based method, the later is a length-scale-based method that includes scale-adaptive simulation characteristics in the turbulence modeling. The SLAU scheme has proven to work well for unsteady flows while the k-kL model works well in non-stationary turbulent flows. As both these flow features are commonly found in industrially important reacting flows, the convection scheme and the turbulence model together will enhance the numerical predictions of such flows.

Self‐reported health (and the extent to which this was associated with partner abuse or psychosocial variables) was investigated in 132 women recruited from a domestic violence service. The survey instrument included abuse disability, life event and daily stress exposure, social support, anger expression style, and perceived health status. The prevalence estimates for this sample were significantly higher than standard estimates across a range of health problems. Regression models demonstrated that whilst the extent of partner abuse predicted the prevalence of three conditions, psychosocial factors were more substantial predictors of health and well‐being in domestic violence victims. Of these, life event frequency and anger expression were the most significant. These findings provide important information about the health of domestic violence victims as they seek support from domestic violence agencies, with relevance for practitioners working with victims who have terminated a violent relationship and for those supporting victims who remain with a violent partner.

The purpose of this paper is to discuss published research in rotorcraft which has taken place in India during the last ten years. The helicopter research is divided into the following parts: health monitoring, smart rotor, design optimization, control, helicopter rotor dynamics, active control of structural response (ACSR) and helicopter design and development. Aspects of health monitoring and smart rotor are discussed in detail. Further work needed and areas for international collaboration are pointed out.

The archival journal papers on helicopter engineering published from India are obtained from databases and are studied and discussed. The contribution of the basic research to the state‐of‐the‐art in helicopter engineering science is brought out.

It is found that strong research capabilities have developed in rotor system health and usage monitoring, rotor blade design optimization, ACSR, composite rotor blades and smart rotor development. Furthermore, rotorcraft modeling and analysis aspects are highly developed with considerable manpower available and being generated in these areas.

Two helicopter projects leading to the “advanced light helicopter” and “light combat helicopter” have been completed by Hindustan Aeronautics Ltd These helicopter programs have benefited from the basic research and also provide platforms for further basic research and deeper industry academic collaborations. The development of well‐trained helicopter engineers is also attractive for international helicopter design and manufacturing companies. The basic research done needs to be further developed for practical and commercial applications.

This is the first comprehensive research on rotorcraft research in India, an important emerging market, manufacturing and sourcing destination for the industry.

A fast assembly scheme for FEM sparse matrices is given. For the solution we compare the conjugate gradient and two preconditionings for this method in the case of plane strain elasticity and orthotropic materials with very stiff coefficients. The influence of fibre orientation on the number of iterations is tested. It is suggested to use the simple incomplete Crout scheme even when the resulting preconditioning matrix is not positive definite.