• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.9 no.3
• /
• pp.65-73
• /
• 1989

Continuum formulations for the expressions of dynamic energy release rates and computational methods for dynamic stress intensity factors are developed for the analysis of dynamic fracture problems subjected to stress wave loading. Explicit volume integral expressions for instantaneous dynamic energy release rates are derived by modeling virtual crack extensions with the dynamic Eulerian-Lagrangian kinematic description. In the finite element applications a finite region around a crack-tip is modeled by using quarter-point singular isoparametric elements, and the volume integrals are evaluated for each crack-tip element during virtual crack extensions while the singularity is maintained. It is shown that the use of the present method is more reliable and accurate for the dynamic fracture analysis than that of other path-independent integral methods when the effects of stress waves are significant.

• 한국전산유체공학회:학술대회논문집
• /
• 2011.05a
• /
• pp.483-487
• /
• 2011

In the present study, a CFD program based on a finite volume method was developed by using an unstructured polyhedral grid system for the accurate simulation with the complex geometry of computational domain. To simulate the transient flow induced by the moving solid object, the program used a fractional step method and a ALE (Algebric Lagrangian-Eulerian) method. The grid deformation for the moving of solid object were performed with a spring analogy based on the center coordinate of each computational grid. To verify the present program with these methodologies, the numerical results of the flow around the fixed and oscillating circular cylinder were compared with the previous numerical results.

• Wind and Structures
• /
• v.18 no.3
• /
• pp.215-233
• /
• 2014

This paper presents a method to extract flutter derivatives of bridge decks based on a combination of the computational fluid dynamics (CFD), system simulations and system identifications. The incompressible solver adopts an Arbitrary Lagrangian-Eulerian (ALE) formulation with the finite volume discretization in space. The imposed sectional motion in heaving or pitching relies on exponential time series as input, with aerodynamic forces time histories acting on the section evaluated as output. System identifications are carried out to fit coefficients of the inputs and outputs of ARMA models, as to establish discrete-time aerodynamic models. System simulations of the established models are then performed as to obtain the lift and moment exerting on the sections to a sinusoidal displacement. It follows that flutter derivatives are identified. The present approaches are applied to a hexagon thin plate and a real bridge deck. The results are compared to the Theodorsen closed-form solution and those from wind tunnel tests. Satisfactory agreements are observed.

• Journal of the Semiconductor & Display Technology
• /
• v.9 no.3
• /
• pp.29-34
• /
• 2010

Nanoimprint lithography (NIL) is an emerging technology enabling cost effective and high throughput nanofabrication. To successfully imprint a nanometer scale patterns, the understanding of the mechanism in nanoimprint forming is essential. In this paper, a numerical analysis of polymer flow in thermal NIL was performed. First, a finite element model of the periodic mold structure with prescribed boundary conditions was established. Then, the volume of fluid (VOF) and grid deformation method were utilized to calculate the free surfaces of the polymer flow based on an Eulerian grid system. From the simulation, the velocity fields and the imprinting pressure for constant imprinting velocity in thermal NIL were obtained. The velocity field is significant because it can directly describe the mode of the polymer deformation, which is the key role to determine the mechanism of nanoimprint forming. Effects of different mold shapes and various thicknesses of polymer resist were also investigated.

• Nuclear Engineering and Technology
• /
• v.42 no.2
• /
• pp.211-218
• /
• 2010

An analysis model on debris transport in the containment floor of pressurized water reactors is developed in which the flow field is calculated by Eulerian conservation equations of mass and momentum and the debris particles are traced by Lagrange equations of motion using the pre-determined flow field data. For the flow field calculation, two-dimensional Shallow Water Equations derived from Navier Stokes equations are solved using the Finite Volume Method, and the Harten-Lax-van Leer scheme is used for accuracy to capture the dry-to-wet interface. For the debris tracing, a simplified two-dimensional Lagrangian particle tracking model including drag force is developed. Advanced schemes to find the positions of particles over the containment floor and to determine the position of particles reflected from the solid wall are implemented. The present model is applied to calculate the transport fraction to the Hold-up Volume Tank in Advanced Power Reactors 1400. By the present model, the debris transport fraction is predicted, and the effect of particle density and particle size on transport is investigated.

• Journal of the Semiconductor & Display Technology
• /
• v.9 no.4
• /
• pp.19-23
• /
• 2010

Nanoimprint lithography (NIL) is one of the most versatile and promising technology for micro/nano-patterning due to its simplicity, high throughput and low cost. Recently, one of the major trends of NIL is large-area patterning. Especially, the research of the application of NIL to TFT-LCD field has been increasing. Technical difficulties to keep the uniformity of the residual layer, however, become severer as the imprinting area increases. In this paper we performed a numerical study for a large area NIL (the $2^nd$ generation TFT-LCD glass substrate ($370{\times}470$ mm)) by using finite element method. First, a simple model considering the surrounding wall was established in order to simulate effectively and reduce the computing time. Then, the volume of fluid (VOF) and grid deformation method were utilized to calculate the free surfaces of the resist flow based on an Eulerian grid system. From the simulation, the velocity fields and the imprinting pressure during the filling process in the NIL were analyzed, and the effect of the surrounding wall and the uniformity of residual layer were investigated.

• Journal of the Korean Society of Safety
• /
• v.15 no.3
• /
• pp.30-39
• /
• 2000

The present paper deals with the numerical simulation for the spray characteristics with swirling turbulent flows and dilution flows from swirl injectors in a simplified can type of gas turbine combustor. The main objective is to investigate the characteristics of swirling turbulent flows with dilution flows and to provide the qualitative results for the spray characteristics such as the droplet distribution and Sauter Mean Diameter(SMD). The gas-phase equations based on Eulerian approach were discretized by Finite Volume Method, together with SIMPLE algorithm and the Reynolds -Stress-Model. The liquid-phase equations based on Lagrangian method were used to predict the droplet behavior. The results of preliminary test are generally in good agreement with experimental data, and show that the anisotropy exists in the primary zone due to swirl velocity and injected air from primary injector, and then gradually decays due to turbulent mixing and consequently near-isotropy occurs in the region between primary and dilution zones. For the spray characteristics, it is indicated that the swirling flows of primary jet region increase the droplet atomization. In addition, it is showed that the swirling flows at the inlet region lead the air-fuel mixture to be distributed near the igniter and can significantly affect the spray behavior in the primary jet region.

• Ocean Systems Engineering
• /
• v.7 no.4
• /
• pp.435-460
• /
• 2017

The main objective of this study is to investigate the turning and zig-zag maneuvering performance of the well-known naval surface combatant DTMB (David Taylor Model Basin) 5415 hull with URANS (Unsteady Reynolds-averaged Navier-Stokes) method. Numerical simulations of static drift tests have been performed by a commercial RANS solver based on a finite volume method (FVM) in an unsteady manner. The fluid flow is considered as 3-D, incompressible and fully turbulent. Hydrodynamic analyses have been carried out for a fixed Froude number 0.28. During the analyses, the free surface effects have been taken into account using VOF (Volume of Fluid) method and the hull is considered as fixed. First, the code has been validated with the available experimental data in literature. After validation, static drift, static rudder and drift and rudder tests have been simulated. The forces and moments acting on the hull have been computed with URANS approach. Numerical results have been applied to determine the hydrodynamic maneuvering coefficients, such as, velocity terms and rudder terms. The acceleration, angular velocity and cross-coupled terms have been taken from the available experimental data. A computer program has been developed to apply a fast maneuvering simulation technique. Abkowitz's non-linear mathematical model has been used to calculate the forces and moment acting on the hull during the maneuvering motion. Euler method on the other hand has been applied to solve the simultaneous differential equations. Turning and zig-zag maneuvering simulations have been carried out and the maneuvering characteristics have been determined and the numerical simulation results have been compared with the available data in literature. In addition, viscous effects have been investigated using Eulerian approach for several static drift cases.