• Geophysics and Geophysical Exploration
• /
• v.11 no.2
• /
• pp.85-92
• /
• 2008

Magnetotelluric (MT) methods are widely applied as an effective exploration technique to geothermal surveys. Two-dimensional (2-D) analysis is frequently used to investigate a complicated subsurface structure in a geothermal region. A 2-D finite-element method (FEM) is usually applied to the MT analysis, but we must pay attention to the accuracy of so-called auxiliary fields. Rodi (1976) proposed an algorithm of improving the accuracy of auxiliary fields, and named it as the MOM method. Because it introduces zeros into the diagonal elements of coefficient matrix of the FEM total equation, a pivoting procedure applied to the symmetrical band matrix makes the numerical solution far less efficient. The MOM method was devised mainly for the inversion analysis, in which partial derivatives of both electric and magnetic fields with respect to model parameters are required. In the case of forward modeling, however, we do not have to resort to the MOM method; there is no need of modifying the coefficient matrix, and the auxiliary fields can be elicited from the regular FEM solution. The computational efficiency of the MOM method, however, can be greatly improved through a sophisticated rearrangement of the total equation.

• Ecology and Resilient Infrastructure
• /
• v.8 no.4
• /
• pp.165-178
• /
• 2021

We studied driftwood behaviors including generation and deposition in a tsunami using a numerical simulation. We used an integrated two-dimensional numerical model, which included a driftwood dynamics model. The study area was Sendai, Japan. Observation data collected by Inagaki et al. (2012) were used to verify the simulation results by comparing them with driftwood deposition patterns. A simplified model was developed to consider the threshold of driftwood generation by the drag force of water flows. To consider the volume of driftwood generated, we estimated the total wood number in the study area using Google Earth. Therefore, we simulated more than 13,000 pieces of driftwood that were generated and transported inland from approximately 300,000 trees that were growing in the forest. The final distribution of the driftwood was similar to the observation data. The reproducibility of the generation and deposition patterns of driftwood showed good agreement in terms of longitudinal deposition pattern. In the future, a sensitivity analysis on driftwood parameters, such as the size of the wood, boundary conditions, and grid size, will be implemented to predict the travel patterns of driftwood. Such modeling will be a useful methodology for disaster prediction based on water flow and driftwood.

• Transactions of the Korean Society of Automotive Engineers
• /
• v.18 no.5
• /
• pp.108-114
• /
• 2010

The injection nozzle of an electro-hydraulic injector is being opened and closed by movement of a injector's needle which is balanced by pressure at the nozzle seat and at the needle control chamber, at the opposite end of the needle. In this study, the effects of needle movement in a piezo-driven injector on unsteady cavitating flows behavior inside nozzle were investigated by cavitation numerical model based on the Eulerian-Lagrangian approach. Aimed at simulating the 3-D two-phase flow behavior, the three dimensional geometry model along the central cross-section regarding of one injection hole with real design data of a piezo-driven diesel injector has been used to simulate the cavitating flows for injection time by at fully transient simulation with cavitation model. The cavitation model incorporates many of the fundamental physical processes assumed to take place in cavitating flows. The simulations performed were both fully transient and 'pseudo' steady state, even if under steady state boundary conditions. As this research results, we found that it could analyze the effect the pressure drop to the sudden acceleration of fuel, which is due to the fastest response of needle, on the degree of cavitation existed in piezo-driven injector nozzle.

• Structural Engineering and Mechanics
• /
• v.74 no.6
• /
• pp.723-735
• /
• 2020

This study is reported the adhesion failure in adhesive bonded composite and specifically for the T-joint structure. Three-dimensional finite element analysis has been performed using a commercial tool and the necessary outcomes are obtained via an eight noded solid element (Solid 185-element) from the library of ANSYS. The structural analysis input has been incurred through ANSYS parametric design language (APDL) code. The normal and shear stress distributions along different layers of the joint structure have been evaluated as the final outcomes. Based on the stress distributions, failure location in the composite joint structure has been identified by using the Tsai-Wu stress failure criterion. It has been found that the failure index is maximum at the interface between flange and web part of the joint (top layer) which indicates the probable location of failure initiation. This kind of failures are considered as adhesion failure and the failure propagation is governed by strain energy release rate (SERR) of fracture mechanics. The different adhesion failure lengths are also considered at the failure location to calculate the SERR values i.e. mode I fracture (opening), mode II fracture (sliding) and mode III fracture (tearing) along the failure front. Also, virtual crack closure technique (VCCT) principle of fracture mechanics steps is used to calculate the above said SERRs. It is found that the mode I SERR is more dominating compared to other two modes of failure for the joint considered. Finally, the influences of various parametric (geometrical and material) effect on SERR of the joint structure are evaluated and discussed in details.

• Wind and Structures
• /
• v.31 no.4
• /
• pp.373-380
• /
• 2020

For large-scale 5MW offshore wind turbines, the discrete equation of fluid domain and the motion equation of structural domain with geometric nonlinearity were built, the three-dimensional modeling of the blade considering fluid-structure interaction (FSI) was achieved by using Unigraphics (UG) and Geometry modules, and the numerical simulation and the analysis of the vibration characteristics for wind turbine structure under rotating effect were carried out based on ANSYS software. The results indicate that the rotating effect has an apparent effect on displacement and Von Mises stress, and the response and the distribution of displacement and Von Mises stress for the blade in direction of wingspan increase nonlinearly with the equal increase of rotational speeds. Compared with the single blade model, the blade vibration period of the whole machine model is much longer. The structural coupling effect reduces the response peak value of the blade displacement and Von Mises stress, and the increase of rotational speed enhances this coupling effect. The maximum displacement difference between two models decreases first and then increases along wingspan direction, the trend is more visible with the equal increase of rotational speed, and the boundary point with zero displacement difference moves towards the blade root. Furthermore, the Von Mises stress difference increases gradually with the increase of rotational speed and decreases nonlinearly from the blade middle to both sides. The results can provide technical reference for the safe operation and optimal design of offshore wind turbines.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.24 no.4
• /
• pp.287-294
• /
• 2012

Sedimentary environment in coastal zone has been changing due to a large number of coastal structures and continuous coastal development. As a result, the environment has been changing. In particular, the economic and environmental damage can occur due to cohesive sediment transport closely related with the fate of pollutants. Due to large sea wall construction the ebb dominance in the Mokpo coastal waters has been clearer. Cohesive sediment transport was simulated by the EFDC model. The simulated SS showed good agreements with the observed SS. From the sensitivity analysis of sediment parameters, we found out that the erosion rate, the critical shear stresses for erosion and deposition, and the settling velocity are important factors in cohesive sediment transport modeling.

• Geophysics and Geophysical Exploration
• /
• v.20 no.4
• /
• pp.226-231
• /
• 2017

Since microgravity survey aims to delineate subsurface density structures in small scale, it requires inversion method, which is able to resolve small scale structures. It can be achieved by adopting a stabilizing functional which separates density boundary distinctly, which is different concept from general inversion routines. We composed Matlab-based interactive two-dimensional microgravity data inversion package containing several kinds of inversion routines with different stabilizing functional, for handling various geologic conditions and survey purposes. Different kinds of inversion routines in the package were verified and examined with representative synthetic data sets generated by numerical modeling. In addition, we applied the developed package to a real microgravity survey data.

• The Journal of Korean Institute of Communications and Information Sciences
• /
• v.42 no.3
• /
• pp.577-584
• /
• 2017

Long Term Evolution device-to-device (LTE D2D) is a key technology to mitigate data traffic load in a cellular system. It facilitates direct data exchange between neighboring users, which is preceded by D2D discovery. Each device advertises its presence to neighboring devices by broadcasting its discovery message. In this paper, we develop a mathematical analysis to assess the probability that discovery messages are successfully transmitted at the D2D discovery stage. We make use of stochastic geometry for modeling spatial statistics of nodes in a two dimensional space. It reflects signal to noise plus interference ratio (SINR) degradation due to resource collision and in-band emission, which leads to the discovery message reception probability being modeled as a function of the distance between the transmitter and the receiver. Numerical results verify that the newly developed analysis accurately estimates discovery message reception probabilities of nodes at the D2D discovery stage.

• Journal of the Korean Society for Railway
• /
• v.16 no.3
• /
• pp.169-174
• /
• 2013

The effect of modifying the shape of a high-performance EMU train on the aerodynamic drag is studied here using Computational Fluid Dynamics(CFD) based on three dimensional Steady-state Navier-Stokes equation and two equation turbulence modeling. FLUENT 12 and Gambit 2.4.6 are employed for a numerical simulation of the aerodynamic drag of a streamlined-shape train as well as a proto type train. The characteristics of the aerodynamic drag of trains in tunnels are analyzed in a comparison with these characteristics in an open space. The contribution of the aerodynamic drag of each case is also investigated to establish principal pertaining to drag reduction for urban trains in tunnels. The aerodynamic drag of a streamlined train was reduced to 9.8% relative to a proto-type train with a blunt nose and a protruding roof facility and underbody shape: the running resistance is expected to be reduced by as much as 4% at a running speed of 80km/h.

• Advances in nano research
• /
• v.7 no.3
• /
• pp.209-221
• /
• 2019

Graphene, with impressive electronic properties, have high potential in the microelectronic field. However, graphene itself is a zero bandgap material which is not suitable for digital logic gates and its application. Thus, much focus is on graphene nanoribbons (GNRs) that are narrow strips of graphene. During GNRs fabrication process, the occurrence of defects that ultimately change electronic properties of graphene is difficult to avoid. The modelling of GNRs with defects is crucial to study the non-idealities effects. In this work, nearest-neighbor tight-binding (TB) model for GNRs is presented with three main simplifying assumptions. They are utilization of basis function, Hamiltonian operator discretization and plane wave approximation. Two major edges of GNRs, armchair-edged GNRs (AGNRs) and zigzag-edged GNRs (ZGNRs) are explored. With single vacancy (SV) defects, the components within the Hamiltonian operator are transformed due to the disappearance of tight-binding energies around the missing carbon atoms in GNRs. The size of the lattices namely width and length are varied and studied. Non-equilibrium Green's function (NEGF) formalism is employed to obtain the electronics structure namely band structure and density of states (DOS) and all simulation is implemented in MATLAB. The band structure and DOS plot are then compared between pristine and defected GNRs under varying length and width of GNRs. It is revealed that there are clear distinctions between band structure, numerical DOS and Green's function DOS of pristine and defective GNRs.