• 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.

• Structural Engineering and Mechanics
• /
• v.91 no.5
• /
• pp.459-468
• /
• 2024

This study presents a frictionless receding contact problem for an orthotropic elastic layer. It is assumed that the layer is supported by two rigid quarter planes and the material of the layer is orthotropic. The layer of thickness h is indented by a rigid cylindrical punch of radius R. The problem is modeled by using the singular integral equation method with the help of the Fourier transform technique. Applying the boundary conditions of the problem the system of singular integral equations is obtained. In this system, the unknowns are the contact stresses and contact widths under the punch and between the layer and rigid quarter planes. The Gauss-Chebyshev integration method is applied to the obtained system of singular integral equations of Cauchy type. Five different orthotropic materials are considered during the analysis. Numerical results are presented to interpret the effect of the material property and the other parameters on the contact stress and the contact width.

• 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.