• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• 제26권2호
• /
• pp.108-120
• /
• 2022

We considered qualitative behaviour of the generalization of Einstein's model of Brownian motion when the key parameter of the time interval of free jump degenerates. Fluids will be characterised by number of particles per unit volume (density of fluid) at point of observation. Degeneration of the phenomenon manifests in two scenarios: a) flow of the fluid, which is highly dispersing like a non-dense gas and b) flow of fluid far away from the source of flow, when the velocity of the flow is incomparably smaller than the gradient of the density. First, we will show that both types of flows can be modeled using the Einstein paradigm. We will investigate the question: What features will particle flow exhibit if the time interval of the free jump is inverse proportional to the density and its gradient ? We will show that in this scenario, the flow exhibits localization property, namely: if at some moment of time t₀ in the region, the gradient of the density or density itself is equal to zero, then for some T during time interval [t₀, t₀ + T] there is no flow in the region. This directly links to Barenblatt's finite speed of propagation property for the degenerate equation. The method of the proof is very different from Barenblatt's method and based on the application of Ladyzhenskaya - De Giorgi iterative scheme and Vespri - Tedeev technique. From PDE point of view it assumed that solution exists in appropriate Sobolev type of space.

• Nonlinear Functional Analysis and Applications
• /
• 제29권1호
• /
• pp.131-164
• /
• 2024

In this research, we study a modified relaxed Tseng method with a single projection approach for solving common solution to a fixed point problem involving finite family of τ-demimetric operators and a quasi-monotone variational inequalities in real Hilbert spaces with alternating inertial extrapolation steps and adaptive non-monotonic step sizes. Under some appropriate conditions that are imposed on the parameters, the weak and linear convergence results of the proposed iterative scheme are established. Furthermore, we present some numerical examples and application of our proposed methods in comparison with other existing iterative methods. In order to show the practical applicability of our method to real word problems, we show that our algorithm has better restoration efficiency than many well known methods in image restoration problem. Our proposed iterative method generalizes and extends many existing methods in the literature.

• 대한조선학회지
• /
• 제22권1호
• /
• pp.21-30
• /
• 1985

This paper deals whth the buckling as well as postbuckling analysis of axisymmertric shells taking the initial deflection effects into account. Incremental equilibrium equations, based on the principle of virtual work, were derived by the finite element method, the successive step-by-step Newton-Raphson iterative technique was adopted. To define the transition pattern of postbuckling behavior from the prebuckling state more accurately, a simple solution method was developed, i.e. the critical load was calculated by the load extrapolation method with the determinant of tangent stiffness matrix and the equilibrium configuration in the immediate postbuckling stage was obtained by perturbation scheme and eigenvalue analysis. Degenerated isoparametric shell elements were used to analyse the axisymmetric shell of revolution. And by the method developed in this paper, the computer program applicable to the nonlinear analysis of both thin and moderately thick shells was constructed. To verify the capabilities and accuracies of the present solution method, the computed results were compared with the results of analytical solutions. These results coincided fairly well in both the small deflection and large deflection ranges. Various numerical analyses were done to show the effect of initial deflection and shape of shells on buckling load and postbuckling behavior. Futhermore, corrected directions of applied loads at every increment steps were used to determine the actual effects of large deflection in non-conservative load systems such as hydrostatic pressure load. The following conclusions can be obtained. (1) The method described in this paper was found to be both economic and effective in calculating buckling load and postbuckling behavior of shell structure. (2) Buckling and postbuckling behavior of spherical caps is critically dependent upon their geometric configuration, i.e. the shape of spherical cap and quantities of the initial deflection. (3) In the analysis of large deflection problems of shells by the incremental method, corrections of the applied load directions are needed at every incremental step to compensate the follower force effects.

• 한국수자원학회논문집
• /
• 제41권8호
• /
• pp.761-772
• /
• 2008

본 연구의 목적은 수공학 분야에서 수치해석이 난해한 문제를 해결하기 위한 모형을 개발하고, 해석해가 존재하는 다양한 수치실험, 즉 하상과 하폭이 함께 변하는 점변부정류 조건에서의 검증, 하상경사가 변화하는 세가지 정상상태 조건의 문제, 그리고 해석해가 있는 마찰하상에 적용함으로써 개발된 모형의 적용성을 검증하기 위한 것이다. 모형의 지배방정식은 보존 법칙을 만족하는 Saint-Venant 적분형 방정식이며, Riemann 해법에 의한 유한체적법이 사용되었다. 질량 및 운동량의 흐름율 계산에 HLL Riemann 근사해법이 사용되었고, 시간-공간에서 2차정확도를 위하여 MUSCL-Hancock 기법이 사용되었다. 본 연구에서는 비선형의 흐름율과 생성항과의 균형을 위하여, 중력과 흐름방향 하폭의 변화로 인한 정수압력에 의한 생성항을 차분하는 새롭고 간편한 기법을 소개하였다. 수치실험 모의결과는 개발된 모형이 생성항을 포함한 다양한 흐름조건에서 정확하고, 견고하며, 매우 안정적임을 보여주고, 또한 수공학 분야에서 일차원 적용에 적합한 모형임을 보여준다.

• Korea-Australia Rheology Journal
• /
• 제15권3호
• /
• pp.131-150
• /
• 2003

A new numerical algorithm of finite element methods is presented to solve high Deborah number flow problems with geometric singularities. The steady inertialess planar 4 : 1 contraction flow is chosen for its test. As a viscoelastic constitutive equation, we have applied the globally stable (dissipative and Hadamard stable) Leonov model that can also properly accommodate important nonlinear viscoelastic phenomena. The streamline upwinding method with discrete elastic-viscous stress splitting is incorporated. New interpolation functions classified as rational interpolation, an alternative formalism to enhance numerical convergence at high Deborah number, are implemented not for the whole set of finite elements but for a few elements attached to the entrance comer, where stress singularity seems to exist. The rational interpolation scheme contains one arbitrary parameter b that controls the singular behavior of the rational functions, and its value is specified to yield the best stabilization effect. The new interpolation method raises the limit of Deborah number by 2∼5 times. Therefore on average, we can obtain convergent solution up to the Deborah number of 200 for which the comer vortex size reaches 1.6 times of the half width of the upstream reservoir. Examining spatial violation of the positive definiteness of the elastic strain tensor, we conjecture that the stabilization effect results from the peculiar behavior of rational functions identified as steep gradient on one domain boundary and linear slope on the other. Whereas the rational interpolation of both elastic strain and velocity distorts solutions significantly, it is shown that the variation of solutions incurred by rational interpolation only of the elastic strain is almost negligible. It is also verified that the rational interpolation deteriorates speed of convergence with respect to mesh refinement.

• Structural Engineering and Mechanics
• /
• 제74권5호
• /
• pp.679-698
• /
• 2020

In this paper, the meshless local Petrov-Galerkin (MLPG) method is developed for dynamic analysis of non-symmetric nanocomposite cylindrical shell equations of elastic wave motion with nonlinear grading patterns under shock loading. The mechanical properties of the nanocomposite cylinder are obtained based on a micro-mechanical model. In this study, four kinds of grading patterns are assumed for carbon nanotube mechanical properties. The displacements can be approximated using shape function so, the multiquadrics (MQ) Radial Basis Functions (RBF) are used as the shape function. In order to discretize the derived equations in time domains, the Newmark time approximation scheme with suitable time step is used. To demonstrate the accuracy of the present method for dynamic analysis, at the first a problem verifies with analytical solution and then the present method compares with the finite element method (FEM), finally, the present method verifies by using the element free Galerkin (EFG) method. The comparison shows the high capacity and accuracy of the present method in the dynamic analysis of cylindrical shells. The capability of the present method to dynamic analysis of non-symmetric nanocomposite cylindrical shell is demonstrated by dynamic analysis of the cylinder with different kinds of grading patterns and angle of nanocomposite reinforcements. The present method shows high accuracy, efficiency and capability to dynamic analysis of non-symmetric nanocomposite cylindrical shell, which it furnishes a ground for a more flexible design.

• 한국지반공학회논문집
• /
• 제27권3호
• /
• pp.75-83
• /
• 2011

일반적인 불포화지반의 거동뿐만 아니라, 열과 관련된 분야, 지반환경 분야 등 다양한 분야에서 불포화 다공질 재료의 열적-수리적-역학적으로 결합된 문제들이 대두되면서 이러한 문제들을 해석하기 위한 수치도구 개발의 필요성이 대두되고 있다. 본 논문에서는 거시적 접근법(macroscopic approach)에 근거하여 불포화지반에 대한 열-수리-역학거동의 수식화를 하였다. 흙, 물, 공기에 대한 질량보존의 법칙, 에너지 보존법칙, 그리고 하중평형 조건식으로부터 결합된(coupled) 4개의 지배방정식을 유도하였다. Galerkin 간략화와 시간적분으로부터 주 변수인 변위(u), 가스압($P_g$), 유체압($P_1$), 온도(T)를 Newton의 반복과정을 이용하여 구하는 유한요소 프로그램(FEM)을 작성하였다. 개발된 프로그램을 이용하여 다공질재료에서 2상 흐름 문제 중 일차원 배수실험(u-$P_g-P_1$), 온도 압밀(u-$P_1$-T), 그리고 지표면 온도변화에 의한 말뚝의 주변지반에 대한 영향(u-$P_1$-T)에 대하여 수치해석을 수행하고 논의하였다.

• 한국전산구조공학회논문집
• /
• 제34권3호
• /
• pp.151-157
• /
• 2021

본 연구에서는 비국부 적분 연산기로 표현되는 페리다이나믹 방정식의 수렴성을 검토한다. 정적/준정적 손상 해석 문제를 효율적으로 해석하기 위해 페리다이나믹 방정식의 implicit 정식화가 필요하다. 이 과정에서 페리다이나믹 비국부 적분 방정식으로부터 대수방정식 형태가 나타나게 되어 시스템 행렬 계산을 위해 많은 시간이 소요되기 때문에, 효율적인 계산을 위해 수렴성이 중요한 요소가 된다. 특히 radial influence 함수를 적분 kernel로 사용하는 경우 fractional Laplacian 적분 방정식이 유도된다. 비국부 적분 연산기의 교윳값 성질에 의해 대수방정식의 condition number가 radial influence 함수의 차수 및 비국부 영역의 크기에 영향을 받는 것이 수학적으로 확인되었다. 본 연구에서는 이를 토대로 균열이 있는 페리다이나믹 정적 해석 문제를 Newton-Raphson 방법으로 해석할 때 적분 커널의 차수, 비국부 영역의 크기 등이 대수방정식의 condition number와 preconditioned conjugate gradient (PCG) 방법으로 계산 시 수렴성 및 계산 시간에 미치는 영향을 수치적으로 분석한다.

• Structural Engineering and Mechanics
• /
• 제85권1호
• /
• pp.105-118
• /
• 2023

The present work is an attempt to develop a simple and accurate finite element formulation for the assessment of thermal shock/thermally induced vibrations in pretwisted and tapered functionally graded material thin (FGM) blades obtained from Voigt and local representative volume elements (LRVE) homogenization models, based on neutral surface approach. The neutral surface of the FGM blade does not coincide with its mid-surface. A finite element model (FEM) is developed using first-order shear deformation theory (FSDT) and the FGM turbine blade is modelled according to the shallow shell theory. The top and the bottom layers of the FGM blade are made of pure ceramic and pure metal, respectively and temperature-dependent material properties are functionally graded in the thickness direction, the position of the neutral surface also depends on the temperature. The material properties are estimated according to two different homogenization models viz., Voigt or LRVE. The top layer of the FGM blade is subjected to high temperature and the bottom surface is either thermally insulated or kept at room temperature. The solution of the nonlinear profile of the temperature in the thickness direction is obtained from the Fourier law of heat conduction in the unsteady state. The results obtained from the present FEM are compared with the benchmark examples. Next, the effect of angle of twist, intensity of thermal shock, variable chord and span and volume fraction index on the transient response due to thermal shock obtained from the two homogenization models viz., Voigt and LRVE scheme is investigated. It is shown that there can be a significant difference in the transient response calculated by the two homogenization models for a particular set of material and geometric parameters.

• 지구물리와물리탐사
• /
• 제9권3호
• /
• pp.241-249
• /
• 2006

한국의 지질 환경은 암석 분포가 매우 다양하고 복잡한 구조 활동의 영향을 받아 지하매질의 이방성 특성이 국부적으로 심하게 변화한다. 기존의 이방성 주시 모델링의 경우 지질 모델을 2차원으로 단순화시킴으로써 이러한 복잡한 지질 환경을 제대로 고려할 수 없었다. 또한 약 이방성 가정을 사용하여 실제로 나타날 수 있는 지하 매질의 심각한 이방성 영향을 주시 모델링에서 고려할 수 없었다. 이에 이 연구에서는 보다 실제적이고 복잡한 3차원 횡등방성 매질(transversely isotropic media)에서 q-P파의 초동 주시 양상을 모사할 수 있는 주시 모델링 알고리듬을 개발하였다. 이 알고리듬에서는 2차원 비선형 주시 내삽(2D nonlinear traveltime interpolation) 기법과 주시의 3차원 격자 채움법(mapping)을 이용한 직접 전파법(direct calculation)을 통해 급격한 물성의 변화에도 주시 계산이 가능하도록 하였다. 또한, 최소 주시 계산과정에서 수치 미분을 통한 최대 경사법(steepest descent method)을 사용하여 약 이방성 가정을 극복하였다. 개발된 알고리듬은 해석해와 비교하여 그 타당성을 검증하였고 3차원 2층구조에 대한 주시 계산을 수행하여 물성이 급격히 변화하는 모델에 대해서도 안정적으로 주시 계산이 이루어짐을 확인하였다. 이 연구에서 개발한 3차원 주시 모델링 알고리듬은 향후 구조보정이나 토모그래피 알고리듬 개발에 사용될 수 있을 것으로 기대한다.