• The Journal of Korean Institute of Electromagnetic Engineering and Science
• /
• v.17 no.2 s.105
• /
• pp.158-164
• /
• 2006

In this paper, we proposed a new algorithm of the inverse scattering for the reconstruction of unknown dielectric scatterers using the finite-difference time-domain method and the design sensitivity analysis. We introduced the design sensitivity analysis based on the gradient information for the fast convergence of the reconstruction. By introducing the adjoint variable method for the efficient calculation, we derived the adjoint variable equation. As an optimal algorithm, we used the steepest descent method and reconstructed the dielectric targets using the iterative estimation. To verify our algorithm, we will show the numerical examples for the two-dimensional $TM^2$ cases.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.37 no.8
• /
• pp.729-736
• /
• 2009

In order to understand flow characteristics and flight mechanism of hummingbird's flapping flight, two-dimensional numerical analysis is carried out on the flapping motion of hummingbird, Selasphorus rufus. Hummingbird's flapping wing motion is realistically modeled from wind tunnel experimental data to perform numerical analysis. Numerical simulation shows that, as freestream velocity changes, wing trajectory is also adjusted and it substantially affects lift and thrust generation mechanism. According to this tendency, flight domain is separated as "low speed" and "high speed" regime, and each flight domain is studied for physical understanding. As a result, the lift generation during downstroke can be explained by the well-known effects, such as leading edge vortex effect, delayed stall, wake capture and so on. In addition, the lift generation during upstroke, the unique character of hummingbird, is also examined by detailed flow analysis. The thrust generation mechanism is investigated by examining the hummingbird's wing bone structure, vortex generation pattern and the resulting pressure gradient.

• Ocean Systems Engineering
• /
• v.8 no.4
• /
• pp.381-407
• /
• 2018

This paper presents a hybrid numerical approach, which combines a two-phase Navier-Stokes model (NS) and the fully nonlinear potential theory (FNPT), for modelling wave-structure interaction. The former governs the computational domain near the structure, where the viscous and turbulent effects are significant, and is solved by OpenFOAM/InterDyMFoam which utilising the finite volume method (FVM) with a Volume of Fluid (VOF) for the phase identification. The latter covers the rest of the domain, where the fluid may be considered as incompressible, inviscid and irrotational, and solved by using the Quasi Arbitrary Lagrangian-Eulerian finite element method (QALE-FEM). These two models are weakly coupled using a zonal (spatially hierarchical) approach. Considering the inconsistence of the solutions at the boundaries between two different sub-domains governed by two fundamentally different models, a relaxation (transitional) zone is introduced, where the velocity, pressure and surface elevations are taken as the weighted summation of the solutions by two models. In order to tackle the challenges associated and maximise the computational efficiency, further developments of the QALE-FEM have been made. These include the derivation of an arbitrary Lagrangian-Eulerian FNPT and application of a robust gradient calculation scheme for estimating the velocity. The present hybrid model is applied to the numerical simulation of a fixed horizontal cylinder subjected to a unidirectional wave with or without following current. The convergence property, the optimisation of the relaxation zone, the accuracy and the computational efficiency are discussed. Although the idea of the weakly coupling using the zonal approach is not new, the present hybrid model is the first one to couple the QALE-FEM with OpenFOAM solver and/or to be applied to numerical simulate the wave-structure interaction with presence of current.

• ETRI Journal
• /
• v.15 no.2
• /
• pp.35-51
• /
• 1993

This paper presents function approximation based on nonparametric estimation. As an estimation model of function approximation, a three layered network composed of input, hidden and output layers is considered. The input and output layers have linear activation units while the hidden layer has nonlinear activation units or kernel functions which have the characteristics of bounds and locality. Using this type of network, a many-to-one function is synthesized over the domain of the input space by a number of kernel functions. In this network, we have to estimate the necessary number of kernel functions as well as the parameters associated with kernel functions. For this purpose, a new method of parameter estimation in which linear learning rule is applied between hidden and output layers while nonlinear (piecewise-linear) learning rule is applied between input and hidden layers, is considered. The linear learning rule updates the output weights between hidden and output layers based on the Linear Minimization of Mean Square Error (LMMSE) sense in the space of kernel functions while the nonlinear learning rule updates the parameters of kernel functions based on the gradient of the actual output of network with respect to the parameters (especially, the shape) of kernel functions. This approach of parameter adaptation provides near optimal values of the parameters associated with kernel functions in the sense of minimizing mean square error. As a result, the suggested nonparametric estimation provides an efficient way of function approximation from the view point of the number of kernel functions as well as learning speed.

• The Transactions of The Korean Institute of Electrical Engineers
• /
• v.57 no.6
• /
• pp.915-923
• /
• 2008

This paper focuses on optimal tuning of nonlinear parameters of a dual-input power system stabilizer(dual-input PSS), which can improve the system damping performance immediately following a large disturbance. Until recently, various PSS models have developed to bring stability and reliability to power systems, and some of these models are used in industry applications. However, due to non-smooth nonlinearities from the interaction between linear parameters(gains and time constants of linear controllers) and nonlinear parameters(saturation output limits), the output limit parameters cannot be determined by the conventional tuning methods based on linear analysis. Only ad hoc tuning procedures('trial and error' approach) have been used. Therefore, the steepest descent method is applied to implement the optimal tuning of the nonlinear parameters of the dual-input PSS. The gradient required in this optimization technique can be computed from trajectory sensitivities in hybrid system modeling with the differential-algebraic-impulsive-switched(DAIS) structure. The optimal output limits of the dual-input PSS are evaluated by time-domain simulation in both a single machine infinite bus(SMIB) system and a multi-machine power system in comparison with those of a single-input PSS.

• Journal of Mechanical Science and Technology
• /
• v.20 no.9
• /
• pp.1514-1524
• /
• 2006

This work presents the three-dimensional analysis of flow and heat transfer performed for a wire-wrapped fuel assembly of liquid metal reactor using Reynolds-averaged Wavier-Stokes analysis in conjunction with 557 model as a turbulence closure. The whole fuel assembly has been analyzed for one period of the wire-spacer using periodic boundary conditions at inlet and outlet of the calculation domain. Three different assemblies, two 7-pin wire-spacer fuel assemblies and one bare rod bundle, apart from the pressure drop calculations for a 19-pin case, have been analyzed. Individual as well as a comparative analysis of the flow field and heat transfer have been discussed. Also, discussed is the position of hot spots observed in the wire-spacer fuel assembly. The flow field in the subchannels of a bare rod bundle and a wire-spacer fuel assembly is found to be different. A directional temperature gradient is found to exist in the subchannels of a wire-spacer fuel assembly Local Nusselt number in the subchannels of wire-spacer fuel assemblies is found to vary according to the wire-wrap position while in case of bare rod bundle, it's found to be constant.

• Advanced Composite Materials
• /
• v.18 no.4
• /
• pp.297-314
• /
• 2009

A design method is proposed to optimize the stacking sequence of laminated composite plates for desired vibration characteristics. The objective functions are the natural frequencies of the laminated plates, and three types of optimization problems are studied where the fundamental frequency and the difference of two adjacent frequencies are maximized, and the difference between the target and actual frequencies is minimized. The design variables are a set of discrete values of fiber orientation angles with prescribed increment in the layers of the plates. The four lamination parameters are used to describe the bending property of a symmetrically laminated plate, and are optimized by a gradient method in the first stage. A new technique is introduced in the second stage to convert from the optimum four lamination parameters into the stacking sequence that is composed of the optimum fiber orientation angles of all the layers. Plates are divided into sub-domains composed of the small number of layers and designed sequentially from outer domains. For each domain, the optimum angles are determined by minimizing the errors between the optimum lamination parameters obtained in the first step and the parameters for all possible discrete stacking sequence designs. It is shown in numerical examples that this design method can provide with accurate optimum solutions for the stacking sequence of vibrating composite plates with various boundary conditions.

• Korea-Australia Rheology Journal
• /
• v.19 no.4
• /
• pp.201-209
• /
• 2007

The present paper is concerned with non-isothermal spherical Couette flow of Oldroyd-B fluid in the annular region between two concentric spheres. The inner sphere rotates with a uniform angular velocity while the outer sphere is kept at rest. Moreover, the two spherical boundaries are maintained at fixed temperature values. Hence, the fluid is effect by two heat sources; namely, the viscous heating and the temperature gradient between the two spheres. The viscoelasticity of the fluid is assumed to dominate the inertia such that the latter can be neglected. An approximate analytical solution of the energy and momentum equations is obtained through the expansion of the dynamical fields in power series of Nahme number. The analysis show that, the temperature variation due to the external source appears in the zero order solution and its effect extends to the fluid velocity distribution up to present second order. Viscous heating contributes in the first and second order solutions. In contrast to isothermal case, a first order axial velocity and a second order stream function fields has been appeared. Moreover, at higher orders the temperature distribution depends on the gap width between the two spheres. Finally, there exist a thermal distribution of positive and negative values depend on their positions in the domain region between the two spheres.

• Journal of Korean Society on Water Environment
• /
• v.28 no.1
• /
• pp.26-37
• /
• 2012

Multi-dimensional hydrodynamic and water quality models are widely used to simulate the physical and biogeochemical processes in the surface water systems such as reservoirs and estuaries. Most of the models have adopted the Eulerian grid modeling framework, mainly because it can reasonably simulate physical dynamics and chemical species concentrations throughout the entire model domain. Determining the optimum grid cell size is important when using the Eulerian grid-based three-dimensional water quality models because the characteristics of species are assumed uniform in each of the grid cells and chemical species are represented by concentration (mass per volume). The objective of this study was to examine the effect of grid-size of a three dimensional hydrodynamic and water quality model (EFDC) on hydrodynamics and mass transport in the Saemangeum Reservoir. Three grid resolutions, respectively representing coarse (CG), medium (MG), and fine (FG) grid cell sizes, were used for a sensitivity analysis. The simulation results of numerical tracer showed that the grid resolution affects on the flow path, mass transport, and mixing zone of upstream inflow, and results in a bias of temporal and spatial distribution of the tracer. With the CG, in particular, the model overestimates diffusion in the mixing zone, and fails to identify the gradient of concentrations between the inflow and the ambient water.

• The Transactions of The Korean Institute of Electrical Engineers
• /
• v.58 no.7
• /
• pp.1294-1303
• /
• 2009

This paper presents a new numerical algorithm for fault discrimination and fault location estimation when occur to arcing ground and arcing line to line on transmission lines. The object of this paper is developed from new numerical algorithm to calculate the fault distance and simultaneously to make a distinction between transient and permanent faults. so the first of object for propose algorithm would be distinguish the permanent from the transient faults. This arcing fault discrimination algorithm is used if calculated value of arc voltage amplitude is greater than product of arc voltage gradient and the length of the arc path, which is equal or greater than the flashover length of a suspension insulator string[1-3]. Also, each algorithm is separated from short distance and long distance. This is difference to with/without capacitance between short to long distance. To test the validity of the proposed algorithms, the results of algorithm testing through various computer simulations are given. The test was simulated in EMTP/ATP simulator under a number of scenarios and calculate of algorithm was used to MATLAB.