• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.19 no.2
• /
• pp.197-215
• /
• 2015

Modeling water flow in variably saturated, porous media is important in many branches of science and engineering. Highly nonlinear relationships between water content and hydraulic conductivity and soil-water pressure result in very steep wetting fronts causing numerical problems. These include poor efficiency when modeling water infiltration into very dry porous media, and numerical oscillation near a steep wetting front. A one-dimensional finite element formulation is developed for the numerical simulation of variably saturated flow systems. First order backward Euler implicit and second order Crank-Nicolson time discretization schemes are adopted as a solution strategy in this formulation based on Picard and Newton iterative techniques. Five examples are used to investigate the numerical performance of two approaches and the different factors are highlighted that can affect their convergence and efficiency. The first test case deals with sharp moisture front that infiltrates into the soil column. It shows the capability of providing a mass-conservative behavior. Saturated conditions are not developed in the second test case. Involving of dry initial condition and steep wetting front are the main numerical complexity of the third test example. Fourth test case is a rapid infiltration of water from the surface, followed by a period of redistribution of the water due to the dynamic boundary condition. The last one-dimensional test case involves flow into a layered soil with variable initial conditions. The numerical results indicate that the Crank-Nicolson scheme is inefficient compared to fully implicit backward Euler scheme for the layered soil problem but offers same accuracy for the other homogeneous soil cases.

• Journal of the Korean Mathematical Society
• /
• v.43 no.5
• /
• pp.1081-1098
• /
• 2006

Mathematical analysis is made on a mesh free method for the compressible Euler equations. In particular, the Moving Least Square Reproducing Kernel (MLSRK) method is employed for space approximation. With the backward-Euler method used for time discretization, existence of discrete solution and it's $L^2-error$ estimate are obtained under a regularity assumption of the continuous solution. The result of numerical experiment made on the biconvex airfoil is presented.

• Journal of computational fluids engineering
• /
• v.2 no.2
• /
• pp.26-34
• /
• 1997

An implicit viscous turbulent flow solver is developed for two-dimensional geon unstructured triangular meshes. The flux terms are discretized based on a cell-centered formulation with the Roe's flux-difference splitting. The solution is advanced in time us backward-Euler time-stepping scheme. At each time step, the linear system of equation approximately solved wi th the Gauss-Seidel relaxation scheme. The effect of turbulence is with a standard k-ε two-equation model which is solved separately from the mean flow equation the same backward-Euler time integration scheme. The triangular meshes are generated advancing-front/layer technique. Validations are made for flows over the NACA 0012 airfoil. Douglas 3-element airfoil. Good agreements are obtained between the numerical result experiment.

• East Asian mathematical journal
• /
• v.33 no.1
• /
• pp.53-65
• /
• 2017

Numerical solutions of the Rosenau-Burgers equation based on the cubic B-spline finite element method are introduced. The backward Euler method is used for discretization in time, and the obtained nonlinear algebraic system is changed to a linear system by the Newton's method. We show that those methods are unconditionally stable. Two test problems are studied to demonstrate the accuracy of the proposed method. The computational results indicate that numerical solutions are in good agreement with exact solutions.

• 한국전산유체공학회:학술대회논문집
• /
• 1997.10a
• /
• pp.29-37
• /
• 1997

An implicit viscous turbulent flow solver is developed for two-dimensional geometries on unstructured triangular meshes. The flux terms are discretized based on a cell-centered finite-volume formulation with the Roe's flux-difference splitting. The solution is advanced in time using an implicit backward-Euler time-stepping scheme. At each time step, the linear system of equations is approximately solved with the Gauss-Seidel relaxation scheme. The effect of turbulence effects is approximated with a standard $k-{\varepsilon}$ two-equation model which is solved separately from the mean flow equations using the same backward-Euler time integration scheme. The triangular meshes are generated using an advancing-front/layer technique. Validations are made for flows over the NACA0012 airfoil and the Douglas 3-element airfoil. Good agreements are obtained between the numerical results and the experiment.

• Journal of Institute of Control, Robotics and Systems
• /
• v.6 no.9
• /
• pp.836-843
• /
• 2000

A sequence planning method is proposed to reduce the assembly time of gantey-type chip mounters with single head. The overall path of the chip mounter is divided into forward and backward path, and formulate the optimization problem is formulated as an transpoetation problem and an Euler's tour problem. The transportation alforithm is applied to find optimal backward path, and Euler's tour algorithm used to generate an assembly sequence. Simulation results are presented to verify the usefulness of the proposed method.

• Journal of applied mathematics & informatics
• /
• v.40 no.3_4
• /
• pp.657-668
• /
• 2022

In the present paper we introduce and study Euler sequence spaces of fractional difference and backward difference operators. We make an effort to prove that these spaces are BK-spaces and linearly isomorphic. Further, Schauder basis for Euler fractional difference sequence spaces $e^{\varsigma}_{0,p}({\Delta}^{(\tilde{\beta})},\;{\nabla}^m)$ and $e^{\varsigma}_{c,p}({\Delta}^{(\tilde{\beta})},\;{\nabla}^m)$ are also elaborate. In addition to this, we determine the 𝛼-, 𝛽- and 𝛾- duals of these spaces.

• Steel and Composite Structures
• /
• v.49 no.1
• /
• pp.1-17
• /
• 2023

Based on Ahmadzadeh-Varvani hardening rule (A-V model), multiaxial ratcheting effect of Z2CND18.12N austenitic stainless steel is simulated by ABAQUS with user subroutine UMAT. The results show that the predicted results of the origin multiaxial A-V model are lower than the experimental data, and it is difficult to control ratcheting strain rate. In order to improve the predicted capability of A-V model, the A-V model is modified. In this study. Moreover, under the assumption of the von Mises yield criterion and normal plasticity flow rule, we develop a numerical algorithm of plastic strain with the improved model to implement the finite element calculation of the model. Internal iteration in the numerical algorithm was implemented with the Euler backward method, which calculated the trial strain for each equilibrium iteration using the consistent tangent matrix. With a user subroutine, the proposed model is programmed into ABAQUS for a user - executable version. By simulating the uniaxial ratcheting of a round bar made of Z2CND18.12N austenitic stainless steel, we observe that the predicted results simulated by ABAQUS with UMAT are compared with the experimental data. The predicted results of the improved multiaxial A-V model are consistent well with the experimental data.

• Advances in nano research
• /
• v.8 no.3
• /
• pp.245-254
• /
• 2020

This work examines the fundamental vibrational characteristics of a spinning CNT-based nano-rotor assuming a nonlocal elasticity Euler-Bernoulli beam theory. The rotary inertia, gyroscopic, and rotor mass unbalance effects are all taken into consideration in the beam model. Assuming a nonlocal theory, two coupled 6th-order partial differential equations governing the vibration of the rotating SWCNT are first derived. In order to acquire the natural frequencies and dynamic response of the nano-rotor system, the nonlinear equations of motion are numerically solved. The nano-rotor system frequency spectrum is shown to exhibit two distinct frequencies: one positive and one negative. The positive frequency is known as to represent the forward whirling mode, whereas the negative characterizes the backward mode. First, the results obtained within the framework of this numerical study are compared with few existing data (i.e., molecular dynamics) and showed an overall acceptable agreement. Then, a thorough and detailed parametric study is carried out to study the effect of several parameters on the nano-rotor frequencies such as: the nanotube radius, the input angular velocity and the small scale parameters. It is shown that the vibration characteristics of a spinning SWCNT are significantly influenced when these parameters are changed.

• Journal of the Korean Mathematical Society
• /
• v.34 no.3
• /
• pp.515-531
• /
• 1997

In this paper, finite difference method is applied to approximate the generalized solutions of Sobolev equations. Using the Steklov mollifier and Bramble-Hilbert Lemma, a priori error estimates in discrete $L^2$ as well as in discrete $H^1$ norms are derived frist for the semidiscrete methods. For the fully discrete schemes, both backward Euler and Crank-Nicolson methods are discussed and related error analyses are also presented.