• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.25 no.4
• /
• pp.224-261
• /
• 2021

The interest in implicit-explicit (IMEX) integration methods has emerged as an alternative for dealing in a computationally cost-effective way with stiff ordinary differential equations arising from practical modeling problems. In this paper, we introduce implicit-explicit second derivative linear multi-step methods (IMEX SDLMM) with error control. The proposed IMEX SDLMM is based on second derivative backward differentiation formulas (SDBDF) and recursive SDBDF. The IMEX second derivative schemes are constructed with order p ranging from p = 1 to 8. The methods are numerically validated on well-known stiff equations.

• Journal of the Korean Society of Combustion
• /
• v.13 no.1
• /
• pp.17-22
• /
• 2008

A partly implicit/quasi-explicit method is introduced for the solution of detailed chemical kinetics with stiff source terms based on the standard fourth-order Runge-Kutta scheme. Present method solves implicitly only the stiff reaction rate equations, whereas the others explicitly. The stiff equations are selected based on the survey of the chemical Jaconian matrix and its Eigenvalues. As an application of the present method constant pressure combustion was analyzed by a detailed mechanism of hydrogen-air combustion with NOx chemistry. The sensitivity analysis reveals that only the 4 species in NOx chemistry has strong stiffness and should be solved implicitly among the 13 species. The implicit solution of the 4 species successfully predicts the entire process with same accuracy and efficiency at half the price.

• Steel and Composite Structures
• /
• v.31 no.3
• /
• pp.261-276
• /
• 2019

The study investigates the free vibration of a nano-scale sandwich beam by an extended high order approach, which has not been reported in the existing literature. First-order shear deformation theory for steel skins and so-called high-order sandwich panel theory for the core are applied. Next, the modified couple stress theory is used for both skins and cores. The Hamilton principle is utilized for deriving equations and corresponding boundary conditions. First, in the study the three-mode shapes natural frequencies for various material parameters are investigated. Also, obtained results are evaluated for two types of stiff and soft cores and isotropic, homogenous steel skins. In the research since the governing equations and also the boundary conditions are nonhomogeneous, therefore some closed-form solutions are not applicable. So, to obtain natural frequencies, the boundary conditions are converted to initial conditions called the shooting method as the numerical one. This method is one of the most robust approaches to solve complex equations and boundary conditions. Moreover, three types of simply supported on both sides of the beam (S-S), simply on one side and clamp supported on the other one (S-C) and clamped supported on both sides (C-C) are scrutinized. The parametric study is followed to evaluate the effect of nano-size scale, geometrical configurations for skins, core and material property change for cores as well. Results show that natural frequencies increase by an increase in skins thickness and core Young modulus and a decrease in beam length, core thickness as well. Furthermore, differences between obtained frequencies for soft and stiff cores increase in higher mode shapes; while, the more differences are evaluated for the stiff one.

• 한국전산유체공학회:학술대회논문집
• /
• 2001.05a
• /
• pp.78-84
• /
• 2001

A dual-time stepping algorithm combined with a parallelized multigrid DADI method is presented to predict the dynamic damping coefficients. The Basic Finner model is chosen to validate the prediction capability of the present unsteady Euler method. The linearity of the pitch- and roll-damping coefficients is shown in the low angular rates and the interesting large drop and stiff increment in transonic region for roll-damping coefficients are explained in detail. Through the analysis for the pressure distributions at Mach number 1.0 to 1.2, the sudden drop results from the normal shock and the stiff increment of roll-damping reflects the transition of the normal shock to the oblique shock. The results also show that the Euler equations can give the damping coefficients with a comparable accuracy.

• Journal of applied mathematics & informatics
• /
• v.8 no.1
• /
• pp.165-180
• /
• 2001

This paper presents a new difference scheme for numerical solution of stiff system of ODE’s. The present study is mainly motivated to develop an absolutely stable numerical method with a high order of approximation. In this work a double implicit A-stable difference scheme with the sixth order of approximation is suggested. Another purpose of this study is to introduce automatic choice of the integration step size of the difference scheme which is derived from the proposed scheme and the one step scheme of the fourth order of approximation. The algorithm was tested by means of solving the Kreiss problem and a chemical kinetics problem. The behavior of the gas explosive mixture (H₂+ O₂) in a closed space with a mobile piston is considered in test problem 2. It is our conclusion that a hydrogen-operated engine will permit to decrease the emitted levels of hazardous atmospheric pollutants.

• Journal of the Korean Society of Combustion
• /
• v.6 no.2
• /
• pp.15-22
• /
• 2001

A pressure-based, unstructured finite volume method has been applied to couple the chemical kinetics and fluid dynamics and to capture effectively and accurately the steep gradient flame field. The pressure-velocity coupling is handled by two methodologies including the pressure-correction algorithm and the projection scheme. A stiff, operator-split projection scheme for the detailed nonequilibrium chemistry has been employed to treat the stiff reaction source terms. The conservative form of the governing equations are integrated over a cell-centered control volume with collocated storage for all transport variables. Computations using detailed chemistry and variable transport properties were performed for two laminar reacting flows: a counterflow hydrogen-air diffusion flame and a lifted methane-air triple flame. Numerical results favorably agree with measurements in terms of the detailed flame structure.

• 한국연소학회:학술대회논문집
• /
• 1999.10a
• /
• pp.219-230
• /
• 1999

The Equations of Chemical kinetics are very stiff, which forces the use of an implicit scheme. The problem of implicit scheme, however, is that the jacobian must be solved at each time step. In this paper, we examined the methodology that can be stable without full chemical jacobian, This method is derived by applying the different time steps to the chemical source term. And the lower triangular chemical jacobian is derived. This is called the preconditioned time differencing method and represents partial implicit method. We show that this method is more stable in chemical kinetics than the full implicit method and that this is more efficient in supersonic combustion problem than the full jacobian method with same accuracy.

• Journal of applied mathematics & informatics
• /
• v.35 no.1_2
• /
• pp.165-180
• /
• 2017

In this study, Taylor matrix algorithm is designed for the approximate solution of linear and non-linear differential equation systems. The algorithm is essentially based on the expansion of the functions in differential equation systems to Taylor series and substituting the matrix forms of these expansions into the given equation systems. Using the Mathematica program, the matrix equations are solved and the unknown Taylor coefficients are found approximately. The presented numerical approach is discussed on samples from various linear and non-linear differential equation systems as well as stiff systems. The computational data are then compared with those of some earlier numerical or exact results. As a result, this comparison demonstrates that the proposed method is accurate and reliable.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.24 no.3
• /
• pp.243-291
• /
• 2020

In this paper, multi-block generalized backward differentiation methods for numerical solutions of ordinary differential and differential algebraic equations are introduced. This class of linear multi-block methods is implemented as multi-block boundary value methods (MB₂ VMs). The root distribution of the stability polynomial of the new class of methods are determined using the Wiener-Hopf factorization of a matrix polynomial for the purpose of their correct implementation. Numerical tests, showing the potential of such methods for output of multi-block of solutions of the ordinary differential equations in the new approach are also reported herein. The methods which output multi-block of solutions of the ordinary differential equations on application, are unlike the conventional linear multistep methods which output a solution at a point or the conventional boundary value methods and multi-block methods which output only a block of solutions per step. The MB₂ VMs introduced herein is a novel approach at developing very large scale integration methods (VLSIM) in the numerical solution of differential equations.

• Journal of computational fluids engineering
• /
• v.5 no.2
• /
• pp.47-54
• /
• 2000

A robust Navier-Stokes code has been developed to efficiently predict hypersonic flows in chemical nonequilibrium. The HLLE+ flux discretization scheme is used to improve accuracy and robustness of hypersonic flow analysis. An efficient LU approximate factorization method is also used to solve the flow equations and species continuity equations in fully coupled fashion to implicitly treat stiff source terms of chemical reactions. The HLLE+ scheme shows lower grid dependency for the wall heating rates than other schemes. The developed code has been used to compute chemical nonequilibrium air flow through expanding hypersonic nozzle and past two and three dimensional blunt-nosed bodies. The results are in good agreement with existing numerical and experimental results.