• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• 제16권3호
• /
• pp.193-204
• /
• 2012

This paper is comparing numerical schemes for a differential equation with convection and fourth-order diffusion. Our model equation is $h_t+(h^2-h^3)_x=-(h^3h_{xxx})_x$, which arises in the context of thin film flow driven the competing effects of an induced surface tension gradient and gravity. These films arise in thin coating flows and are of great technical and scientific interest. Here we focus on the several numerical methods to apply the model equation and the comparison and analysis of the numerical results. The convection terms are treated with well known WENO methods and the diffusion term is treated implicitly. The diffusion and convection schemes are combined using a fractional step-splitting method.

• International Journal of Aerospace System Engineering
• /
• 제7권2호
• /
• pp.6-12
• /
• 2020

The precise prediction of future position of satellite depends on the accurate determination of orbit, which is also helpful in performing orbit maneuvers and trajectory correction maneuvers. For estimating the orbit of satellite many methods are being used. Some of the conventional methods are based on (i) Differential Correction (DC) (ii) Extended Kalman Filter (EKF). In this paper, Differential Evolution (DE) is used to determine the orbit. Orbit Determination using DC and EKF requires some initial guess of the state vector to initiate the algorithm, whereas DE does not require an initial guess since a wide range of bounds for the design unknown variables (orbital elements) is sufficient. This technique is uniformly valid for all orbits viz. circular, elliptic or hyperbolic. Simulated observations have been used to demonstrate the performance of the method. The observations are generated by including random noise. The simulation model that generates the observations includes the perturbation due to non-spherical earth up to second zonal harmonic term.

• Geomechanics and Engineering
• /
• 제24권4호
• /
• pp.337-348
• /
• 2021

Brillouin Optical Frequency Domain Analysis (BOFDA) is a distributed fiber optic sensing (DFOS) technique that has unique advantages for performance monitoring of piles. However, the complicated production process and harsh operating environment of offshore PHC pipe piles make it difficult to apply this method to pile load testing. In this study, sensing cables were successfully pre-installed into an offshore PHC pipe pile directly for the first time and the BOFDA technique was used for in-situ monitoring of the pile under axial load. High-resolution strain and internal force distributions along the pile were obtained by the BOFDA sensing system. A finite element analysis incorporating the Degradation and Hardening Hyperbolic Model (DHHM) was carried out to evaluate and predict the performance of the pile, which provides an improved insight into the offshore pile-soil interaction mechanism.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• 제26권4호
• /
• pp.333-342
• /
• 2022

In this paper, we perform a direct comparison study of real experimental data for domain rearrangement and the Cahn-Hilliard (CH) equation on the dynamics of morphological evolution. To validate a mathematical model for physical phenomena, we take initial conditions from experimental images by using an image segmentation technique. The image segmentation algorithm is based on the Mumford-Shah functional and the Allen-Cahn (AC) equation. The segmented phase-field profile is similar to the solution of the CH equation, that is, it has hyperbolic tangent profile across interfacial transition region. We use unconditionally stable schemes to solve the governing equations. As a test problem, we take domain rearrangement of lipid bilayers. Numerical results demonstrate that comparison of the evolutions with experimental data is a good benchmark test for validating a mathematical model.

• 대한수학회보
• /
• 제60권5호
• /
• pp.1237-1252
• /
• 2023

A fundamental problem in coding theory is to find n_q(k, d), the minimum length n for which an [n, k, d]_q code exists. We show that some q-divisible optimal linear codes of dimension 4 over 𝔽_q, which are not of Belov type, can be constructed geometrically using hyperbolic quadrics in PG(3, q). We also construct some new linear codes over 𝔽_q with q = 7, 8, which determine n₇(4, d) for 31 values of d and n₈(4, d) for 40 values of d.

• 천문학회보
• /
• 제45권1호
• /
• pp.55.1-55.1
• /
• 2020

In an attempt to investigate the nonlinear dynamics such as shock, shear, and turbulence associated with ultra-relativistic jets, we develop a new relativistic hydrodynamics (RHD) code based on the weighted essentially non-oscillatory (WENO) scheme. It is a 5th-order accurate, finite-difference scheme, which has been widely used for solving hyperbolic systems of conservation equations. The code is parallelized with MPI and OpenMP. Through an extensive set of tests, the accuracy and efficiency of different WENO reconstructions, and different time discretizations are assessed. Different implementations of the equation of state (EOS) for relativistic fluid are incorporated, As the fiducial setup for simulations of ultra-relativistic jets, we adopt the EOS in Ryu et al. (2006) to treat arbitrary adiabatic index of relativistic fluid, the WENO-Z reconstructions to minimize numerical dissipation without loss of stability, and the strong stability preserving Runge-Kutta (SSPRK) method to achieve stable time stepping with large CFL numbers. In addition, the code includes a high-order flux averaging along the transverse directions for multi-dimensional problems, and the modified eigenvalues for the acoustic modes to effectively control the carbuncle instability. We find that the new code performs satisfactorily simulations of ultra-relativistic jets.

• 한국지반공학회지:지반
• /
• 제12권4호
• /
• pp.157-178
• /
• 1996

기존의 보강토벽체에 주로 이용되어온 steel strict등 고강도 인장보강재는 주변 뒤채움흙에 비해 상대적으로 변형이 작기 때문에, 설계검토시 과강재 자체에서 유발되는 변형의 크기에 대해서는 크게 유의할 필요가 없었다. 그러나 비교적 저강도인 섬유보강재의 경우, 한계상태에서 예상되는 섬유보강재 자체의 변형량은 주변 뒤채움흙의 소성변형 유발에 필요시 되는 변형량을 종종 초과하게 되며, 이와같은 크기의 과도한 변형량은 보강토벽체 구조체 자체의 안정성 확보 측면에서 허용할 수 없는 경우가 대부분이다. 결국 보증토벽체 구조체의 전면부 발생변위에 대한 일반적인 허용조건을 충족하기 위해서는, 극한강도 보다 훨씬 작은 크기의 강도가 섬유보강재의 경우 발휘하는 것으로 보아야 할 것이며, 따라서 최종적인 구조체 안정검토를 위해서는 보강재 자체의 예상변형량에 대한 평가가 섬유보강재의 경우 특히 중요시 된다. 보강재의 인장응력 -변형률 관계는 강보강재의 경우 선형탄성거동으로 가정할 수 있으나, 섬 유보강재의 경우에는 일반적으로 비 선형거동을 나타낸다. 본 연구에서는 쌍곡선 함수를 이용하여 섬유보강재의 비선형 거동특성을 모델링하였으며,또한 뒤채움흙 다짐으로 인한 유발응력등을 고려하기 위해 Ehrlich SE Mitchell, Duncan등이 제안한 방법을 수정하여 섬유 보강토벽체의 안정 해석법을 제시하였다. 본 안정 해석법 에서는 침투수압의 영향 및 뒤채움흙의 구속효과에 따른 섬유보강재의 부분적인 상대강성 변화 등을 고려하였으며, 이를 토대로 깊이별 각 섬유보 강재의 최대인장력 및 변형량 등의 예측이 가능하다. 본 연구에서는 제시하리라 하는 안정해석법의 적용성을 위해, paraweb polyester fibre multicord, non-woven polyester 지오텍스타일 및 knitted polyester 지오그리드 등 3가지 종류 보강재의 인장응력-변형률 관계 실험결과를 회귀분석하여 쌍곡선 함수형태로 이와같은 섬유보 강재의 비선형거동을 모델링하였다. 또한 이를 토대로 한 븐 연구 해석법의 적합성 검토를 위해, Ho & Rowe가 제시한 유한요소해석결과 및 LCPC, FHWA등에서 시행한 시험결과와 깊이별 각 섬유보강재의 최대인장력,변형량 및 지점별 변형률 등에 대해서도 비교하였다. 아울러 섬유 보강재의 상대강성, 뒤채움흙의 깊이별 구속효과의 정도, 다짐정도 및 침투수압 등이 각 섬유보강재의 변형량 및 전체적인 변형형태 등에 미치는 영향을 종합적으로 분석하였다.

• Geomechanics and Engineering
• /
• 제23권4호
• /
• pp.313-325
• /
• 2020

Slope displacement and factor of safety (FOS) of a slope are two aspects that reflect the stability of a slope. However, the traditional limit equilibrium (LE) methods only give the result of the slope FOS and cannot be used to solve for the slope displacement. Therefore, developing a LE method to obtain the results of the slope FOS and slope displacement has significance for engineering applications. Based on a force-displacement coupled mode, this work expands the LE Morgenstern-Price (M-P) method. Except for the mechanical equilibrium conditions of a sliding body adopted in the traditional M-P method, the present method introduces a nonlinear model of the shear stress and shear displacement. Moreover, the energy equation satisfied by a sliding body under a small slope displacement is also applied. Therefore, the double solutions of the slope FOS and horizontal slope displacement are established. Furthermore, the flow chart for the expanded LE M-P method is given. By comparisons and analyses of slope examples, the present method has close results with previous research and numerical simulation methods, thus verifying the feasibility of the present method. Thereafter, from the parametric analysis, the following conclusions are obtained: (1) the shear displacement parameters of the soil affect the horizontal slope displacement but have little effect on the slope FOS; and (2) the curves of the horizontal slope displacement vs. the minimum slope FOS could be fitted by a hyperbolic model, which would be beneficial to obtain the horizontal slope displacement for the slope in the critical state.

• 한국지반공학회지:지반
• /
• 제5권2호
• /
• pp.19-32
• /
• 1989

매설흄관의 거동을 해석하기 위하여 흙의 쌍곡선 모델, 비선형 이력모델 및 흙-구조체 접촉면 모 델을 사용한 유한요소법에 의하여 흄관에 작용하는 연직(설계)하중을 계산하였으며, 이를 고전이론 인 Marston-Spangler이론 등과 비교·분석하였다. 수치해석을 통하여 고전이론에서 고려되지 못했 던 각종 토질상수와 터파기 형태별의 영향을 포함하여 흙덮개, 터파기 폭 등의 영향을 관찰하였다. 부가적으로 철장 실무자들이 쉽게 이용할 수 있도록 설계하중을 구하는 방법을 제안하였다.

• 한국지반공학회:학술대회논문집
• /
• 한국지반공학회 1999년도 봄 학술발표회 논문집
• /
• pp.3-46
• /
• 1999

A new method called SRICOS is proposed to predict the scour depth z versus time t around a cylindrical bridge pier of diameter D founded in clay. The steps involved are ; 1. taking samples at the bridge pier site, 2. testing them in an Erosion Function Apparatus called the EFA to obtain the scour rate z versus the hydraulic shear stress applied $\tau$, 3. predicting the maximum shear stress r max which will be induced around the pier by the water flowing at ν Ο before the scour hole starts to develop, 4. using the measured z versus r curve to obtain the initial scour rate zi corresponding to r max , 5. predicting the maximum depth of scour zmax for the pier, 6. using zi and zmarx to develop the hyperbolic function describing the scour depth z versus time t curve, and 7. reading the z vs. t curve at a time corresponding to the duration of the flood to find the scour depth which will develop around the pier. A new apparatus is developed to measure the z vs t curve of step 2, a series of advanced numerical simulations are performed to develop an equation for the $\tau$ max value of step 3, and a series of flume tests are performed to develop an equation for the zmax value of step 5. The method is evaluated by comparing predictions and measurements in 42 flume experiments.