• Kyungpook Mathematical Journal
• /
• v.61 no.2
• /
• pp.383-393
• /
• 2021

In this paper, numerical techniques are presented for solving initial value problems of fractional differential equations with variable coefficients. The method is derived by applying a Taylor vector approximation. Moreover, the operational matrix of fractional integration of a Taylor vector is provided in order to transform the continuous equations into a system of algebraic equations. Furthermore, numerical examples demonstrate that this method is applicable and accurate.

• International Journal of Highway Engineering
• /
• v.18 no.3
• /
• pp.21-32
• /
• 2016

PURPOSES : In this paper, the analytical solutions suggested to simulate the behavior of rheological fluids were rigorously re-derived and investigated for fixed conditions to evaluate the applicability for the solutions on a mini-cone slump test of cement paste. The selected solutions with proper boundary conditions can be used as reference solutions to evaluate the performance of numerical simulation approaches, such as the discrete element method. METHODS : The slump, height, and spread radius for the given boundary and yield stress conditions that are determined by five different analytical solutions are compared. RESULTS : The analytical solution based on fluid mechanics for pure shear flow shows similar results to that for intermediate flow at low yield stresses. The fluid mechanics-based analytical solution resulted in a very similar trend to the geometry-based analytical solution. However, it showed a higher slump at high yield stress and lower slump at low yield stress ranges than the geometry-based analytical model. The analytical solution based on the mini-cone geometry was not significantly affected by the yield criteria, such as von Mises and Tresca. CONCLUSIONS : Even though differences among the analytical solutions in terms of slump and spread radius existed, the difference can be considered insignificant when the solutions were used as reference to evaluate the appropriateness of numerical approaches, such as the discrete element method.

• 한국지구물리탐사학회:학술대회논문집
• /
• 2003.11a
• /
• pp.238-243
• /
• 2003

As there are many advantages on underground caverns, such as safety and operation, they can also be used for gas storage purpose. When liquefied gas is stored underground, the cryogenic temperature of the gas will affect the stability of the storage cavern. In order to store the liquefied gas successfully, it is essential to estimate the exact temperature distribution of the rock mass around the cavern. In this study, an analytic solution and a conceptual model that can estimate three-dimensional temperature distribution around the storage cavern are suggested. When calculating the heat transfer within a solid, it is likely to consider the solid as the intersection of two or more infinite or semi-infinite geometries. Therefore heat transfer solution for the solid is expressed by the product of the dimensionless temperatures of the geometries, which are used to form the combined solid. Based on the multi-dimensional transient heat transfer theory, the analytic solution is successfully derived by assuming the cavern shape to be of simplified geometry. Also, a conceptual model is developed by using the analytic solution of this study. By performing numerical experiments of this multi-dimensional model, the temperature distribution of the analytic solution is compared with that of numerical analysis and theoretical solutions.

• Advances in Computational Design
• /
• v.5 no.1
• /
• pp.55-89
• /
• 2020

For the first time, an accurate analytical solution, based on coupled three-dimensional (3D) piezoelasticity equations, is presented for free vibration analysis of the angle-ply elastic and piezoelectric flat laminated panels under arbitrary boundary conditions. The present analytical solution is applicable to composite, sandwich and hybrid panels having arbitrary angle-ply lay-up, material properties, and boundary conditions. The modified Hamiltons principle approach has been applied to derive the weak form of governing equations where stresses, displacements, electric potential, and electric displacement field variables are considered as primary variables. Thereafter, multi-term multi-field extended Kantorovich approach (MMEKM) is employed to transform the governing equation into two sets of algebraic-ordinary differential equations (ODEs), one along in-plane (x) and other along the thickness (z) direction, respectively. These ODEs are solved in closed-form manner, which ensures the same order of accuracy for all the variables (stresses, displacements, and electric variables) by satisfying the boundary and continuity equations in exact manners. A robust algorithm is developed for extracting the natural frequencies and mode shapes. The numerical results are reported for various configurations such as elastic panels, sandwich panels and piezoelectric panels under different sets of boundary conditions. The effect of ply-angle and thickness to span ratio (s) on the dynamic behavior of the panels are also investigated. The presented 3D analytical solution will be helpful in the assessment of various 1D theories and numerical methods.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.14 no.3
• /
• pp.175-187
• /
• 2010

This paper presents the numerical valuation of the two-asset step-down equitylinked securities (ELS) option by using the operator-splitting method (OSM). The ELS is one of the most popular financial options. The value of ELS option can be modeled by a modified Black-Scholes partial differential equation. However, regardless of whether there is a closedform solution, it is difficult and not efficient to evaluate the solution because such a solution would be represented by multiple integrations. Thus, a fast and accurate numerical algorithm is needed to value the price of the ELS option. This paper uses a finite difference method to discretize the governing equation and applies the OSM to solve the resulting discrete equations. The OSM is very robust and accurate in evaluating finite difference discretizations. We provide a detailed numerical algorithm and computational results showing the performance of the method for two underlying asset option pricing problems such as cash-or-nothing and stepdown ELS. Final option value of two-asset step-down ELS is obtained by a weighted average value using probability which is estimated by performing a MC simulation.

• Wind and Structures
• /
• v.26 no.4
• /
• pp.205-214
• /
• 2018

In this paper, the thermo-mechanical buckling characteristics of functionally graded (FG) size-dependent Timoshenko nanobeams subjected to an in-plane thermal loading are investigated by presenting a Navier type solution for the first time. Material properties of FG nanobeam are supposed to vary continuously along the thickness according to the power-law form and the material properties are assumed to be temperature-dependent. The small scale effect is taken into consideration based on nonlocal elasticity theory of Eringen. The nonlocal governing equations are derived based on Timoshenko beam theory through Hamilton's principle and they are solved applying analytical solution. According to the numerical results, it is revealed that the proposed modeling can provide accurate critical buckling temperature results of the FG nanobeams as compared to some cases in the literature. The detailed mathematical derivations are presented and numerical investigations are performed while the emphasis is placed on investigating the effect of the several parameters such as material distribution profile, small scale effects and aspect ratio on the critical buckling temperature of the FG nanobeams in detail. It is explicitly shown that the thermal buckling of a FG nanobeams is significantly influenced by these effects. Numerical results are presented to serve as benchmarks for future analyses of FG nanobeams.

• Structural Engineering and Mechanics
• /
• v.72 no.5
• /
• pp.617-629
• /
• 2019

Based on the finite element method of traditional straight Euler-Bernoulli beams and the coupled relations between linear displacement and angular displacement of a pre-twisted Euler-Bernoulli beam, the shape functions and stiffness matrix are deduced. Firstly, the stiffness of pre-twisted Euler-Bernoulli beam is developed based on the traditional straight Euler-Bernoulli beam. Then, a new finite element model is proposed based on the displacement general solution of a pre-twisted Euler-Bernoulli beam. Finally, comparison analyses are made among the proposed Euler-Bernoulli model, the new numerical model based on displacement general solution and the ANSYS solution by Beam188 element based on infinite approach. The results show that developed numerical models are available for the pre-twisted Euler-Bernoulli beam, and which provide more accurate finite element model for the numerical analysis. The effects of pre-twisted angle and flexural stiffness ratio on the mechanical property are investigated.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.31 no.6 s.261
• /
• pp.499-506
• /
• 2007

This study aims to propose evaluation methods of chemical washing performance and estimate the washing capability by flowing detergent/water solution for application to home appliances such as dishwashers. Standard pollutant is stearic acid. A numerical study is also tried using a SIMPLER code. Preliminary experiments are performed by varying the concentration and temperature of the solution. From the pre-experiments, 10 minute pre-curing time is found to be necessary to remove the stearic acid. Stoichiometric ratio and detergent consumption coefficient of reaction between the detergent and stearic are estimated following a proposed method. Washing experiments of pollutant to compare with the numerical results are performed. The relative errors between the experimental and the numerical results with pre-curing time included are less than 7%. In conclusion, important mechanisms of chemical washing are revealed and methods of predicting washing performance are well established.

• Wind and Structures
• /
• v.18 no.3
• /
• pp.235-252
• /
• 2014

An effective numerical technique to calculate the reactions of a multi-spanned transmission line conductor system, under arbitrary loads varying along the spans, is developed. Such variable loads are generated by High Intensity Wind (HIW) events in the form of tornadoes and downburst. First, a semi-closed form solution is derived to obtain the displacements and the reactions at the ends of each conductor span. The solution accounts for the nonlinearity of the system and the flexibility of the insulators. Second, a numerical scheme to solve the derived closed-form solution is proposed. Two conductor systems are analyzed under loads resulting from HIW events for validation of the proposed technique. Non-linear Finite Element Analyses (FEA) are also conducted for the same two systems. The responses resulting from the technique are shown to be in a very good agreement with those resulting from the FEA, which confirms the technique accuracy. Meanwhile, the semi-closed form technique shows superior efficiency in terms of the required computational time. The saving in computational time has a great advantage in predicting the response of the conductors under HIW events, since this requires a large number of analyses to cover different potential locations and sizes of those localized events.

• Proceedings of the Korean Geotechical Society Conference
• /
• 2001.03a
• /
• pp.41-48
• /
• 2001

In this study, two factors are simultaneously considered for assessing tunnel face stability: one is the effective stress acting on the tunnel face calculated by upper bound solution; and the other is the seepage force calculated by numerical analysis under the condition of steady-state groundwater flow. The seepage forces calculated by numerical analysis are compared with the results of a model test. From the results of derivations of the upper bound solution with the consideration of seepage forces acting on the tunnel face, it could be found that the minimum support pressure for the face stability is equal to the sum of effective support pressure and seepage pressure acting on the tunnel face. Also it could be found that the average seepage pressure acting on the tunnel face is proportional to the hydrostatic pressure at the same elevation and the magnitude is about 22% of the hydrostatic pressure for the drainage type tunnel and about 28% for the water-proof type tunnel. The model tests performed with a tunnel model had a similar trend with the seepage pressure calculated by numerical analysis. From the model tests it could be also found that the collapse at the tunnel face occurs suddenly and leads to unlimited displacement.