• Journal of applied mathematics & informatics
• /
• v.26 no.5_6
• /
• pp.1273-1287
• /
• 2008

In this paper, the numerical integration method for general singularly perturbed two point boundary value problems with mixed boundary conditions of both left and right end boundary layer is presented. The original second order differential equation is replaced by an approximate first order differential equation with a small deviating argument. By using the trapezoidal formula we obtain a three term recurrence relation, which is solved using Thomas Algorithm. To demonstrate the applicability of the method, we have solved four linear (two left and two right end boundary layer) and one nonlinear problems. From the results, it is observed that the present method approximates the exact or the asymptotic expansion solution very well.

• Journal of applied mathematics & informatics
• /
• v.32 no.1_2
• /
• pp.39-53
• /
• 2014

In this paper, a non-asymptotic method is presented for solving singularly perturbed delay differential equations whose solution exhibits a boundary layer behavior. The second order singularly perturbed delay differential equation is replaced by an asymptotically equivalent first order neutral type delay differential equation. Then, Simpson's integration formula and linear interpolation are employed to get three term recurrence relation which is solved easily by Discrete Invariant Imbedding Algorithm. Some numerical examples are given to validate the computational efficiency of the proposed numerical scheme for various values of the delay and perturbation parameters.

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.2 no.3
• /
• pp.127-131
• /
• 2010

A two-layer fluid with free surface is simulated in the time domain by a two-dimensional potential-based Numerical Wave Tank (NWT). The developed NWT is based on the boundary element method and a leap-frog time integration scheme. A whole domain scheme including interaction terms between two layers is applied to solve the boundary integral equation. The time histories of surface elevations on both fluid layers in the respective wave modes are verified with analytic results. The amplitude ratios of upper to lower elevation for various density ratios and water depths are also compared.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.22 no.1
• /
• pp.36-44
• /
• 2018

The thermal response of carbon/phenolic used in a solid rocket nozzle liner was analyzed. In this paper, the numerical analysis of the thermal response of carbon/phenolic consists of (1) the integration equation of the boundary layer to obtain the convective heat transfer coefficient of the combustion gas on the rocket nozzle wall and (2) 1-D finite difference method for heat conduction of carbon/phenolic to calculate the ablation, char, and temperature. The calculated result was compared with the result of a blast-tube-type test motor. It is found that the calculated result shows good agreement with the thermal response of the test motor, except at the vicinity of the throat insert.

• Journal of Advanced Marine Engineering and Technology
• /
• v.19 no.2
• /
• pp.10-19
• /
• 1995

Heat trnasfer for absorbing and emitting media in laminar layer along the cylinders has been analyzed. Governing equation are transformed to local nonsimilarity equations by the dimensional analysis. The effects of the Stark number, Prandtl number, Optical radius and wall emissivity are mainly investigated. For the formal solution a numerical integration is performed and the results are compared with those obtained by P-1 and P-3 approximation. The results show that boundary layers consist of conduction-convection-radiation layer near the wall and convection-radiation layer far from the wall. As the Stark number of wall emissivity increases the local radiative heat flux is increased. The Pradtl number or curvature variations do not affect the radiative heat flux from the wall, but The Prandtl number or wall emissivity variations affect the conduction heat flux. Consequently the total heat flux from the wall are affected by the Prandtl number or wall emissivity variation.

• Structural Engineering and Mechanics
• /
• v.83 no.1
• /
• pp.67-77
• /
• 2022

In this study, frictionless continuous and discontinuous contact problems of a magneto-electro-elastic layer in the presence of the body force were discussed. The layer was indented by a rigid cylindrical insulating punch and supported by a rigid substrate without bond. Applying the Fourier integral transform technique, the general expressions of the problem were derived in the presence of body force. Thanks to the boundary conditions, the singular integral equations were obtained for both the continuous and the discontinuous contact cases. Gauss-Chebyshev integration formulas were used to transform the singular integral equations into a set of nonlinear equations. Contact width under the punch, initial separation distance, critical load, separation regions and contact stress under the punch and between the layer, and substrate were given as a result.

• Structural Engineering and Mechanics
• /
• v.48 no.2
• /
• pp.241-255
• /
• 2013

The receding contact problem for two elastic layers whose elastic constants and heights are different supported by two elastic quarter planes is considered. The lower layer is supported by two elastic quarter planes and the upper elastic layer is subjected to symmetrical distributed load whose heights are 2a on its top surface. It is assumed that the contact between all surfaces is frictionless and the effect of gravity force is neglected. The problem is formulated and solved by using Theory of Elasticity and Integral Transform Technique. The problem is reduced to a system of singular integral equations in which contact pressures are the unknown functions by using integral transform technique and boundary conditions of the problem. Stresses and displacements are expressed depending on the contact pressures using Fourier and Mellin formula technique. The singular integral equation is solved numerically by using Gauss-Jacobi integration formulation. Numerical results are obtained for various dimensionless quantities for the contact pressures and the contact areas are presented in graphics and tables.

• International Journal of Aeronautical and Space Sciences
• /
• v.7 no.2
• /
• pp.155-174
• /
• 2006

Parabolized stability equations for compressible flows in general curvilinear coordinate system are derived to deal with a broad range of transition prediction problems on complex geometry. A highly accurate finite difference PSE code has been developed using an implicit marching procedure. Compressible and incompressible flat plate flow stability under two-dimensional and three¬dimensional disturbances has been investigated to test the present code. Results of the present computation are found to be in good agreement with the multiple scale analysis and DNS data. Stability calculation results by the present PSE code for compressible boundary layer at Mach numbers ranging from 0.02 to 1.5 are also presented and are again seen to be as accurate as the spectral method.

• Structural Engineering and Mechanics
• /
• v.78 no.5
• /
• pp.585-597
• /
• 2021

This paper presents a comparative study of analytical method, finite element method (FEM) and Multilayer Perceptron (MLP) for analysis of a contact problem. The problem consists of a functionally graded (FG) layer resting on a half plane and pressed with distributed load from the top. Firstly, analytical solution of the problem is obtained by using theory of elasticity and integral transform techniques. The problem is reduced a system of integral equation in which the contact pressure are unknown functions. The numerical solution of the integral equation was carried out with Gauss-Jacobi integration formulation. Secondly, finite element model of the problem is constituted using ANSYS software and the two-dimensional analysis of the problem is carried out. The results show that contact areas and the contact stresses obtained from FEM provide boundary conditions of the problem as well as analytical results. Thirdly, the contact problem has been extended based on the MLP. The MLP with three-layer was used to calculate the contact distances. Material properties and loading states were created by giving examples of different values were used at the training and test stages of MLP. Program code was rewritten in C++. As a result, average deviation values such as 0.375 and 1.465 was obtained for FEM and MLP respectively. The contact areas and contact stresses obtained from FEM and MLP are very close to results obtained from analytical method. Finally, this study provides evidence that there is a good agreement between three methods and the stiffness parameters has an important effect on the contact stresses and contact areas.

• Structural Engineering and Mechanics
• /
• v.51 no.4
• /
• pp.685-705
• /
• 2014

This paper presents a new formulation for forward scalar wave simulations in semi-infinite media. Perfectly-Matched-Layers (PMLs) are used as a wave absorbing boundary layer to surround a finite computational domain truncated from the semi-infinite domain. In this work, a hybrid formulation was developed for the simulation of scalar wave motion in two-dimensional PML-truncated domains. In this formulation, displacements and stresses are considered as unknowns in the PML domain, while only displacements are considered to be unknowns in the interior domain. This formulation reduces computational cost compared to fully-mixed formulations. To obtain governing wave equations in the PML region, complex coordinate stretching transformation was introduced to equilibrium, constitutive, and compatibility equations in the frequency domain. Then, equations were converted back to the time-domain using the inverse Fourier transform. The resulting equations are mixed (contain both displacements and stresses), and are coupled with the displacement-only equation in the regular domain. The Newmark method was used for the time integration of the semi-discrete equations.