• Korea-Australia Rheology Journal
• /
• v.13 no.2
• /
• pp.97-106
• /
• 2001

Three-dimensional flow analysis was performed by using the control volume finite-element method for design of a profile extrusion die. Because polymer melt behavior is complicated and cross-sectional shape of the profile extrusion die is changing continuously, the fluid flow within the die must be analyzed three-dimensionally. A commercially available polypropylene is used for theoretical and experimental investigations. Material properties are assumed to be constant except for the viscosity. The 5-constant modified Cross model is used for the numerical analysis. A test problem is examined in order to verify the accuracy of the numerical method. Simulations are performed for conditions of three different screw speeds and three different die temperatures. Predicted pressure distribution is compared with the experimental measurements and the results of the previous two-dimensional study. The computational results obtained by using three dimensional CVFEM agree with the experimental measurements and are more accurate than those obtained by using the two-dimensional cross-sectional method. The velocity profiles and the temperature distributions within several cross-sections of the die are given as contour plots.

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

Direct numerical simulation results of aeolian tones generated by a two-dimensional obstacle (circular cylinder, square cylinder, NACA0012 airfoil) in a uniform flow are presented and the generation and propagation mechanisms of the sound are discussed. The unsteady compressible Navier-Stokes equations are solved by a highly-accurate finite difference scheme over the entire region from near to far fields. The direct numerical simulation results are also compared with the results obtained by Curle's acoustic analogy.

• Journal of Electrical Engineering and Technology
• /
• v.1 no.3
• /
• pp.396-405
• /
• 2006

The visual analysis of buried channel (Be) devices such as buried channel MOSFETs and CCDs (Charge Coupled Devices) is investigated to give better understanding and insight for their electrical behaviours using a 3-dimensional (3-D) numerical simulation. This paper clearly demonstrates the capability of the numerical simulation of 'EVEREST' for characterising the analysis of a depletion mode MOSFET and BC CCD, which is a simulation software package of the semiconductor device. The inverse threshold and punch-through voltages obtained from the simulations showed an excellent agreement with those from the measurement involving errors of within approximately 1.8% and 6%, respectively, leading to the channel implanted doping profile of only approximately $4{\sim}5%$ error. For simulation of a buried channel CCD an advanced adaptive discretising technique was used to provide more accurate analysis for the potential barrier height between two channels and depletion depth of a deep depletion CCD, thereby reducing the CPU running time and computer storage requirements. The simulated result for the depletion depth also showed good agreement with the measurement. Thus, the results obtained from this simulation can be employed as the input data of a circuit simulator.

• 한국전산유체공학회:학술대회논문집
• /
• 1995.04a
• /
• pp.5-29
• /
• 1995

A three-dimensional flow field induced by two trains passing by each other inside a tunnel is studied based on the numerical simulation of the three-dimensional compressible Euler/Navier-Stokes equations formulated in the finite difference approximation. Domain decomposition method with the FSA(fortified solution algorithm) interface scheme is used to treat this moving-body problem. The computed resluts show basic characteristic of the flow field created when two trains passing by each other. History of the pressure distributions and the aerodynamic forces acting on the trains are mailnly discussed. The results indicate that the phenomenon is complicated due to the interaction of the flow induced by two trains. Strong side force occurs between the two trains when the front portion of the opposite train passes by. It fluctuates rapidly and maximum suction force occurs when two trains are aligned side by side. The results also indicate the effectiveness of the present numerical method for moving boundary problems.

• Nuclear Engineering and Technology
• /
• v.29 no.6
• /
• pp.459-467
• /
• 1997

The steady-state behavior of recycling neutral atoms in a tokamak edge region has been analyzed through a two-dimensional Monte Carlo simulation. A particle tracking algorithm used in earlier research on the neutral particle transport is applied to this Monte Carlo simulation in order to perform more accurate calculations with the EDGETRAN code which was previously developed for a two-dimensional edge plasma transport in the authors' laboratory. The physical model of neutral recycling includes charge-exchange and ionization interactions between plasmas and neutral atoms. The reflection processes of incident particles on the device wall are described by empirical formulas. Calculations for density, energy, and velocity distributions of neutral deuterium-tritium atoms have been carried out for a medium-sized tokamak with a double-null configuration based on the KT-2 conceptual design. The input plasma parameters such as plasma density, ion and electron temperatures, and ion fluid velocity are provided from the EDGETRAN calculations. As a result of the present numerical analysis, it is noticed that a significant drop of the neutral atom density appears in the region of high plasma density and that the similar distribution of neutral energy to that of plasma ions is present as frequently reported in other studies. Relations between edge plasma conditions and the neutral recycling behavior are discussed from the numerical results obtained herein.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2008.03a
• /
• pp.557-560
• /
• 2008

A two-dimensional particle-in-cell(PIC) simulation with a Monte-Carlo Collision(MCC) has been developed to investigate the discharge characteristics of the acceleration channel of a HET. The dynamics of electrons and ions are treated with PIC method at the time scale of electrons in order to investigate the particle transport. The densities of charged particles are coupled with Poisson's equation. Xenon neutrals are injected from the anode and experience elastic, excitation, and ionization collisions with electrons, and are scattered by ions. These collisions are simulated by using an MCC model. The effects of control parameters such as magnetic field profile, electron current density, and the applied voltage have been investigated. The secondary electron emission on the dielectric surface is also considered.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2008.03a
• /
• pp.589-596
• /
• 2008

Performance of a pulsed heat source high temperature inert gas plasma MHD electrical power generator, which can be one of the candidates of space-based laser-to-electrical power converter, is examined by a time dependent two dimensional numerical simulation. In the present MHD generator, the inert gas is assumed to be ideally heated to about $10^4K$ pulsed-likely within short time(${\sim}1{\mu}s$) in a stagnant energy input volume, and the energy of high temperature inert gas is converted to the electricity with the medium of pure inert gas plasma without seeding. The numerical simulation results show that an enthalpy extraction ratio(=electrical output energy/pulsed heat energy) of several tens of % can be achieved, which is the same level as the conventional seeded non-equilibrium plasma MHD generator. Although there still exist many phenomena to be clarified and many problems to be overcome in order to realize the system, the pulsed heat source high temperature inert gas MHD generator is surely worth examining in more detail.

• Journal of Fashion Business
• /
• v.19 no.5
• /
• pp.1-16
• /
• 2015

The appearances and geometry structures of knitted fabrics have important effects on their functions as textile fabrics. Structural design of the woven fabric, prior to the manufacturing processes in the weaving mill, often leads to a similar predictable appearance in the final outcome with the corresponding weave design. The increase of the employment of elastic textile yarns in knitting fabrics for comfort stretch or outdoor sports wear knit products has, however, resulted in difficulties in predicting the final appearance of the knit structure design. Due to the stretchability and exceptional recovery behavior of the elastic yarns such as polyurethane elastomeric yarns, the appearance of the final product often differs from the initial knit design. At textile CAD program for preparing tricot knit designs has been employed in this study to predict the two dimensional appearance of the design. The similarities between the designs and corresponding knit products seem to be acceptable for the two-dimensional textile CAD program in this study. However, when elastomeric yarns are partially employed in the polyester filament tricot product, a considerable amount of departure from the design is apparent due to the constriction and/or deformation of property differences in the elastomeric yarns and polyester filament yarns. Therefore, another purpose of this study is to measure the departure of the final tricot product from the initial tricot design, especially in the case employing elastomeric yarns in the knit structure together with regular polyester filament yarns. For measuring the three-dimensional departure, a 3D scanning system has been used for the mesh reconstruction of the fabric specimen. Hopefully, the result from this study will be used as a guide to modify and improve the current textile CAD program proposed for the two-dimensional simulation of the tricot.

• Journal of Astronomy and Space Sciences
• /
• v.6 no.1
• /
• pp.1-15
• /
• 1989

Magnetospheic substorm in the magnetotail region is studied numerically by means of a three dimensional MHD code. The analytic solution for the quiet magnetotail is employed as an initial configuration. The localized solar wind is modeled to enter the simulation domain through the boundaries located in the magnetotail lobe region. As a result of the interaction between the solar wind and the magnetosphere, the magnetic field lines are stretched, and the plasma sheet becomes thinner and thinner. When the current driven resistivity is generated, magnetic reconnection is triggered by this resistivity. The resulting plasma jetting is found to be super-magnetosonic. Although the plasmoid formation and its tailward motion is not quite clear as in the two dimensional simulation, which is mainly because of the numerical model chosen for the present simulation, the rarification of the plasmas near the x-point is observed. Field aligned currents are observed in the late expansive stage of the magnetospheric substorm. These field aligned currents flow from the tail toward the ionosphere on the dawn side from the ionosphere to ward the tail on the dusk side, namely in the same sense of the region 1 current. As the field aligned currents develop, it is found that the cross tail current in the earth side midnight section of the magnetic x-point is reduced.

• Structural Engineering and Mechanics
• /
• v.23 no.1
• /
• pp.29-42
• /
• 2006

This paper describes the design scheme of the three-dimensional structures based on the concept of the cellular automata simulation. The cellular automata simulation is performed according to the local rule. In this paper, the local rule is derived in the mathematical formulation from the optimization problem. The cell density is taken as the design variable. Two objective functions are defined for reducing the total weight of the structure and obtaining the fully stressed structure. The constraint condition is defined for defining the local rule. The penalty function is defined from the objective functions and the constraint condition. Minimization of the penalty function with respect to the design parameter leads to the local rule. The derived rule is applied to the design of the three-dimensional structure first. The final structure can be obtained successfully. However, the computational cost is expensive. So, in order to reduce the computational cost, the material parameters $c_1$ and $c_2$ and the value of the cell rejection criterion (CRC) are changed. The results show that the computational cost depends on the parameters and the CRC value.