• Ocean Systems Engineering
• /
• v.1 no.4
• /
• pp.317-336
• /
• 2011

A time domain finite element based method is employed to analyze wave radiation by multiple cylinders. The nonlinear free surface and body surface boundary conditions are satisfied based on the perturbation method up to the second order. The first- and second-order velocity potential problems at each time step are solved through a finite element method (FEM). The matrix equation of the FEM is solved through an iteration and the initial solution is obtained from the result at the previous time step. The three-dimensional (3D) mesh required is generated based on a two-dimensional (2D) hybrid mesh on a horizontal plane and its extension in the vertical direction. The hybrid mesh is generated by combining an unstructured grid away from cylinders and two structured grids near the cylinder and the artificial boundary, respectively. The fluid velocity on the free surface and the cylinder surface are calculated by using a differential method. Results for various configurations including two-cylinder and four-cylinder cases are provided to show the mutual influence due to cylinders on the first and second waves and forces.

• Interaction and multiscale mechanics
• /
• v.6 no.2
• /
• pp.173-195
• /
• 2013

Meshfree methods are known to have the capability to overcome the strict regularization requirements and numerical instabilities that encumber the finite element method (FEM) in large deformation problems. They are also more naturally suited for problems involving material perforation and fragmentation. To take advantage of the high efficiency of FEM and high accuracy of meshfree methods, a coupled finite element (FE) and reproducing kernel (RK, one of the meshfree approximations) formulation is described in this paper. The coupling of FE and RK approximation is implemented in an evolutionary fashion, where the extent and location of the evolution is dependent on a triggering criteria provided by the material constitutive laws. To enhance computational efficiency, Gauss quadrature is applied to integrate both FE and RK domains so that no state variable transfer is required when mesh conversion is performed. To control the hourglassing that might occur with 1-point integrated hexahedral grids, viscous type hourglass control is implemented. Meanwhile, the FEM version of the K&C concrete (KCC) model was modified to make it applicable in both FE and RK formulations. Results using this code and the KCC model are shown for the modeling of concrete responses under quasi-static, blast and impact loadings. These analyses demonstrate that fragmentation phenomena of the sort commonly observed under blast and impact loadings of concrete structures was able to be realistically captured by the coupled formulation.

• Nuclear Engineering and Technology
• /
• v.33 no.1
• /
• pp.93-101
• /
• 2001

The spacer grid is one of the main structural components in the fuel assembly, which supports the fuel rods, guides cooling water, and protects the system from an external impact load, such as earthquakes. Therefore, the mechanical and structural properties of the spacer grids must be extensively examined while designing them. In this paper, a numerical method for predicting the buckling strength of spacer grids is presented. Numerical analyses on the buckling behavior of the spacer grids are performed for a various array of sizes of the grids considering that the spacer grid is an assembled structure with thin-walled plates and imposing proper boundary conditions by nonlinear dynamic finite element method using ABAQUS/Explicit. Buckling tests on several numbers of specimens of the spacer grid were also carried out in order to compare the results between the test and the simulation result. The drop test is accomplished by dropping a carriage on the specimen at a pre-determined position. From this test, the specimens are buckled only at the uppermost and the lowermost layer among the multi-cells, which is similar to the local buckling at the weakest point of the grid structure. The simulated results also similarly predicted the local buckling phenomena and were found to give good correspondence with the experimental values for the thin-walled grid structures.

• Journal of the Korean Society of Fisheries and Ocean Technology
• /
• v.42 no.1
• /
• pp.38-43
• /
• 2006

The 3-D modeling of ocean finite element method(OFEM) using $k-{\varepsilon}$ turbulent model and tetrahedron grids has been used to investigate the internal wave generation during the expansion of the deep water from the open sea to the shelf with a simple shape, which can be widely used in the fields of submarine development, ocean environment and meteorology, etc. In this paper, the detailed configuration of internal wave with its length and height and also the distribution of salinity and turbulent kinematic energy, etc. were derived. It is hoped that this OFEM method can be successfully applied to the numerical calculation of internal wave for and the oceanographic problems (tidal flows around underwater hill, plateau, Georges Bank, etc.) and ocean engineering problems(flow past artificial sea reefs) in future.

• Journal of the Korean Mathematical Society
• /
• v.54 no.1
• /
• pp.193-213
• /
• 2017

We investigate an efficient approximation of solution to stochastic Burgers equation driven by an additive space-time noise. We discuss existence and uniqueness of a solution through the Orstein-Uhlenbeck (OU) process. To approximate the OU process, we introduce the Karhunen-$Lo{\grave{e}}ve$ expansion, and sparse grid stochastic collocation method. About spatial discretization of Burgers equation, two separate finite element approximations are presented: the conventional Galerkin method and Galerkin-conservation method. Numerical experiments are provided to demonstrate the efficacy of schemes mentioned above.

• Proceedings of the KIEE Conference
• /
• 1996.07c
• /
• pp.1536-1540
• /
• 1996

The Design of are quenching parts of molded case circuit breakers depends on the utilization of strength and distribution of the magnetic field by which the arc is forced. The magnetic field causes the are to move into a set of V-slotted iron grids, where the are is extinguished rapidly. This paper present the effective method 10 design V-slotted iron plates of the are breaking chamber of molded case circuit breakers. This magnetic force was calculated by using the flux densities in the arc which are obtained by three dimensional finite element method, as a result of that this paper verified by testing that a grid model which has biggest magnetic force is excellent in the are quenching ability.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.24 no.2
• /
• pp.143-159
• /
• 2020

We propose a consistent numerical method for elliptic interface problems with nonhomogeneous jumps. We modify the discontinuous bubble immersed finite element method (DB-IFEM) introduced in (Chang et al. 2011), by adding a consistency term to the bilinear form. We prove optimal error estimates in L² and energy like norm for this new scheme. One of the important technique in this proof is the Bramble-Hilbert type of interpolation error estimate for discontinuous functions. We believe this is a first time to deal with interpolation error estimate for discontinuous functions. Numerical examples with various interfaces are provided. We observe optimal convergence rates for all the examples, while the performance of early DB-IFEM deteriorates for some examples. Thus, the modification of the bilinear form is meaningful to enhance the performance.

• Journal of computational fluids engineering
• /
• v.1 no.1
• /
• pp.88-97
• /
• 1996

This paper treats an adaptive finite-element method for the viscous compressible flow governed by Navier-Stokes equations in two dimensions. The numerical algorithm is the two-step Taylor-Galerkin mettled using unstructured triangular grids. To increase accuracy and stability, combined moving node method and grid refinement method have been used for grid adaption. Validation of the present algorithm has been made by comparing the present computational results with the existing experimental data and other numerical solutions. Four benchmark problems are solved for demonstration of the present numerical approach. They include a subsonic flow over a flat plate, the Carter flat plate problem, a laminar shock-boundary layer interaction. and finally a laminar flow around NACA0012 airfoil at zero angle of attack and free stream Mach number of 0.85. The results indicates that the present adaptive triangular grid method is accurate and useful for laminar viscous flow calculations.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.21 no.11
• /
• pp.1538-1551
• /
• 1997

Kinetic energy conservation for fixed and moving grids is examined in time-accurate finite element computation of fully unsteady inviscid flows. As numerical algorithms, fractional step method (FSM) and modified SIMPLE are used. To simulate the flow in moving grid system, arbitrary Lagrangian-Eulerian (ALE) method is adopted. In the present study, the energy conserving time integration rule for finite element algorithm is proposed and discussed schematically. It is shown that the discretization by Crank-Nicolson in time and Galerkin (central difference) in space must be used to ensure energy conservation. The developed code has been tested for a standing vortex in fixed or moving grid system, sloshing in a tank and propagation of a solitary wave, and has been shown to be a completely energy conserving algorithm.

• Journal of Mechanical Science and Technology
• /
• v.16 no.5
• /
• pp.683-695
• /
• 2002

The objective of the present study is to analyze the fluid flow with moving boundary using a finite element method. The algorithm uses a fractional step approach that can be used to solve low-speed flow with large density changes due to intense temperature gradients. The explicit Lax-Wendroff scheme is applied to nonlinear convective terms in the momentum equations to prevent checkerboard pressure oscillations. The ALE (Arbitrary Lagrangian Eulerian) method is adopted for moving grids. The numerical algorithm in the present study is validated for two-dimensional unsteady flow in a driven cavity and a natural convection problem. To extend the present numerical method to engine simulations, a piston-driven intake flow with moving boundary is also simulated. The density, temperature and axial velocity profiles are calculated for the three-dimensional unsteady piston-driven intake flow with density changes due to high inlet fluid temperatures using the present algorithm. The calculated results are in good agreement with other numerical and experimental ones.