• Journal of the Korean Mathematical Society
• /
• v.58 no.3
• /
• pp.553-569
• /
• 2021

In this paper, we present and analyze a fully discrete numerical method for solving the time-fractional diffusion wave equation: ∂^β_tu - div(a∇u) = f, 1 < β < 2. We first construct a difference formula to approximate ∂^β_tu by using an interpolation of derivative type. The truncation error of this formula is of O(△t^2+δ-β)-order if function u(t) ∈ C^2,δ[0, T] where 0 ≤ δ ≤ 1 is the Hölder continuity index. This error order can come up to O(△t^3-β) if u(t) ∈ C³ [0, T]. Then, in combinination with the linear finite volume discretization on spatial domain, we give a fully discrete scheme for the fractional wave equation. We prove that the fully discrete scheme is unconditionally stable and the discrete solution admits the optimal error estimates in the H¹-norm and L₂-norm, respectively. Numerical examples are provided to verify the effectiveness of the proposed numerical method.

• Nuclear Engineering and Technology
• /
• v.51 no.5
• /
• pp.1181-1194
• /
• 2019

A new moving-mesh Finite Volume Method (FVM) for the efficient solution of the two-dimensional neutron diffusion equation is introduced. Many other moving-mesh methods developed to solve the neutron diffusion problems use a relatively large number of sophisticated mathematical equations, and so suffer from a significant complexity of mathematical calculations. In this study, the proposed method is formulated based on simple mathematical algebraic equations that enable an efficient mesh movement and CV deformation for using in practical nuclear reactor applications. Accordingly, a computational framework relying on a new moving-mesh FVM is introduced to efficiently distribute the meshes and deform the CVs in regions with high gradient variations of reactor power. These regions of interest are very important in the neutronic assessment of the nuclear reactors and accordingly, a higher accuracy of the power densities is required to be obtained. The accuracy, execution time and finally visual comparison of the proposed method comprehensively investigated and discussed for three different benchmark problems. The results all indicated a higher accuracy of the proposed method in comparison with the conventional fixed-mesh FVM.

• Journal of computational fluids engineering
• /
• v.15 no.2
• /
• pp.71-78
• /
• 2010

As a step toward accurate prediction of droplet impingement and ice accretion on aircraft, an Eulerian-based droplet impingement and ice accretion code for air flows around an airfoil containing water droplets is developed. A CFD solver based on the finite volume method was also developed to solve the clean airflow. The finite-volume-based approach for simulating droplet impingement on an airfoil was employed owing to its compatibility with the CFD solver and robustness. For ice accretion module, a simple model based on the control volume is combined with the droplet impingement module that provides the collection efficiency. To validate the present code, it is compared with NASA Glenn IRT (Icing Research Tunnel) experimental data and other well-known icing codes such as LEWICE and FENSAP-ICE. It is shown that the collection efficiency and shape of ice accretion are in good agreement with previous experimental and simulation results.

• International Journal of Aerospace System Engineering
• /
• v.2 no.2
• /
• pp.40-43
• /
• 2015

A hybrid mesh based finite volume compressible flow code (PolySim) has been developed recently. Instead of the simple average method for the gradients of variables at each face, the volume average is applied for the calculation of the viscous flux. What is more, an improved Green Gauss method for the calculation of the gradient is also presented. These two techniques will improve both the accuracy and robustness of the code. The aerodynamic performance of this in-house cell centered code is examined by several widely-used bench-mark test cases. These cases include flows over flat plate and RAE 2822 etc. The comparisons on results between calculation and experiment are conducted. They show that the code can produce good numerical results which agree well with the corresponding experiment data.

• Proceedings of the KSME Conference
• /
• 2001.11b
• /
• pp.551-556
• /
• 2001

The filling pattern technique based on the finite element method and Eulerian mesh advancement approach has been developed to analyze incompressible transient viscous flow with free surfaces. The governing equation for flow analysis is Navier-Stokes equation including inertia and gravity effects. The penalty and predictor-corrector methods are used effectively for finite element formulation. The flow front surface and the volume inflow rate are calculated using the filling pattern technique to select an adequate pattern among four filling patterns at each triangular control volume. Using the proposed numerical technique, the collapse of a dam has been analyzed to predict flow phenomenon of fluid and the predicted front positions versus time have been compared with the reported experimental result.

• Journal of Korea Water Resources Association
• /
• v.31 no.3
• /
• pp.279-290
• /
• 1998

The height and speed of the shock wave are critical data in flood-control operations or in the design of channel walls and bridges along rivers with high flow velocities. Therefore, a numerical model is needed for simulating flow discontinuity over a wide range of conditions. In this study, a governing equation. As a Riemann solver Roe(1981)'s one is used. The model employs the modified MUSCL for handling the unstructured grids in this research. this model that adopts the explicit tradditional twl dimmensional dam break problems, two hydraulic dam break model is simulations, and a steady state simulation in a curved channel. Conclusions of this research are as follows : 1) the finite volume method can be combined with the Godonov-type method that is useful for modeling shocks. Hence, the finite volume method is suitable for modeling shocks. 2) The finite volume model combined with the modified MUSCL is successful in modeling shock. Therefore, modified MUSCL is proved to be valid.

• 한국전산유체공학회:학술대회논문집
• /
• 2003.10a
• /
• pp.6-6
• /
• 2003

• Proceedings of the Korean Society for Technology of Plasticity Conference
• /
• 1994.10a
• /
• pp.133-140
• /
• 1994

In this paper, we have proposed a new method to determine the initial billet for the forged products using a function approximation in neural networks. the architecture of neural network is a three-layer neural network and the back propagation algorithm is employed to train the network. By utilizing the ability of function approximation of neural network, an optimal billet is determined by applying nonlinear mathematical relationship between shape ratio in the initial billet and the final products. A volume of incomplete filling in the die is measured by the rigid-plastic finite element method. The neural network is trained with the initial billet shape ratio and that of the un-filled volume. After learning, the system is able to predict the filling region which are exactly the same or slightly different to results of finite element method. It is found that the prediction of the filling shape ratio region can be made successfully and the finite element method results are represented better by the neural network.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.18 no.4
• /
• pp.1009-1018
• /
• 1994

An attempt is made to develop a kind of hybrid numerical method for computations of the thermal stresses during a solidification process. In this algorithm, the phase-change heat transfer analysis is perrformed by a finite volume method(FVM) and the thermal stress analysis in a solidifying body by a finite element method(FEM). The temperatures at the grid points calculated in the heat transfer analysis are transferred to those of gauss points in elements by a bi-cubic surface patch technique for the thermal stress analysis. A hyperbolic-sine constitutive law is used to prescribe the inelastic strain rate of material. Results for the unidirectional solidification process of a pure aluminum are compared with those of others and shows good agreement.

• 한국전산유체공학회:학술대회논문집
• /
• 2006.10a
• /
• pp.57-60
• /
• 2006

This paper presents an assessment of vertex-centered and cell-centered finite-volume methods on unstructured meshes. The results indicate that the vertex-centered method is more reliable than the cell-centered method.