• 천문학회보
• /
• 제38권2호
• /
• pp.101.2-101.2
• /
• 2013

We present a new computational fluid dynamic (CFD) simulation code. The code employs the moving and polyhedral unstructured mesh scheme, which is known as a superior approach to the conventional SPH (smoothed particle hydrodynamics) and AMR (adaptive mesh refinement) schemes. The code first generates unstructured meshes by the Voronoi tessellation at every time step, and then solves the Riemann problem for surfaces of every Voronoi cell to update the hydrodynamic states as well as to move former generated meshes. For the second-order accuracy, the MUSCL-Hancock scheme is implemented. To increase efficiency for generating Voronoi tessellation we also develop the incremental expanding method, by which the CPU time is turned out to be just proportional to the number of particles, i.e., O(N). We will discuss the applications of our code in the context of cosmological simulations as well as numerical experiments for galaxy formation.

• 천문학회보
• /
• 제39권2호
• /
• pp.65.2-65.2
• /
• 2014

We present a new hydrodynamic simulation code based on the Voronoi tessellation for estimating the density precisely. The code employs both of Lagrangian and Eulerian description by adopting the movable mesh scheme, which is superior to the conventional SPH (smoothed particle hydrodynamics) and AMR (adaptive mesh refinement) schemes. The code first generates unstructured meshes by the Voronoi tessellation at every time step, and then solves the Riemann problem for all surfaces of each Voronoi cell so as to update the hydrodynamic states as well as to move current meshes. Besides, the IEM (incremental expanding method) is devised to compute the Voronoi tessellation to desired degree of speed, thereby the CPU time is turned out to be just proportional to the number of particles, i.e., O(N). We discuss the applications of our code in the context of cosmological simulations as well as numerical experiments for galaxy formation.

• 한국전산구조공학회논문집
• /
• 제27권2호
• /
• pp.71-78
• /
• 2014

회전익 항공기 중 군에서 운용하는 기동헬기는 전장상황에서 운용되기 때문에 연료셀 피탄 상황에 직면할 가능성이 높다. 연료셀 피탄에 따른 내부압력 증가로 내부폭발이나 화재가 발생할 수 있으며, 이는 승무원의 생존 가능성에 치명적인 영향을 주게 된다. 따라서, 승무원의 생존성을 극대화하기 위해서는 연료셀이 직면 가능한 극한 상황을 예측하여 설계에 반영해야 한다. 항공기 연료셀 설계시 고려해야 하는 데이타는 피탄에 의한 연료셀 내부압력, 수압램 영향에 의한 연료셀 자체 및 금속피팅부 응력, 탄환의 운동에너지 등이 포함될 수 있다. 이러한 설계 데이터 확보를 위해서는 실물 시험을 수행하는 것이 가장 바람직하지만, 시간과 비용의 부담과 더불어 시험실패와 같은 시행착오 위험성으로 많은 제약이 따른다. 따라서, 사전에 다양한 설계 데이터 예측과 시행착오의 최소화를 위해서는 피탄 상황에 대한 수치해석이 필요하다. 본 연구에서는 입자법을 사용하여 연료셀 피탄 조건에 대한 유체-구조 연성 수치해석을 수행하였다. 수치해석은 전용 충돌해석 프로그램인 LS-DYNA를 사용하였고, 결과로 얻어진 탄의 거동과 에너지, 연료셀 내부압력과 등가응력의 평가를 통해 연료셀 설계와 관련한 데이터 확보 가능성을 타진하였다.

• 천문학회보
• /
• 제37권2호
• /
• pp.101.2-101.2
• /
• 2012

Shocks are ubiquitous in astrophysical plasmas: bow shocks are formed by the interaction of solar wind with planetary magnetic fields, and supernova explosions and jets produce shocks in interstellar and intergalactic spaces. The global morphologies of these shocks are usually described by a set of magnetohydrodynamic (MHD) equations which tacitly assumes local thermal equilibrium, and the resulting Rankine-Hugoniot shock jump conditions are applied to obtain the relationship between the upstream and downstream physical quantities. While thermal equilibrium can be achieved easily in collisional fluids, it is generally believed that collisions are infrequent in astrophysical settings. In fact, shock widths are much smaller than collisional mean free paths and a variety of kinetic phenomena are seen at the shock fronts according to in situ observations of planetary shocks. Hence, both the MHD and kinetic equations have been adopted in theoretical and numerical studies to describe different aspects of the physical phenomena associated with astrophysical shocks. In this paper, we present the results of 3D relativistic particle-in-cell (PIC) simulations for ion-electron plasmas, with focus on the shock structures: when a jet propagates into an unmagnetized ambient plasma, a shock forms in the nonlinear stage of the Weibel instability. As the shock shows the structures that resemble those predicted in MHD systems, we compare the results with those predicted in the MHD shocks. We also discuss the thermalization processes of the upstream flows based on the time evolutions of the phase space and the velocity distribution, as well as the wave spectra analyses.

• 천문학회보
• /
• 제44권2호
• /
• pp.55.2-55.2
• /
• 2019

The firehose instability is driven by a pressure anisotropy in a magnetized plasma when the temperature along the magnetic field is higher than the perpendicular temperature. Such condition occurs commonly in astrophysical and space environments, for instance, when there are beams aligned with the background magnetic field. Recently, it was argued that, in weak quasi-perpendicular shocks in the high-β intracluster medium (ICM), shock-reflected electrons propagating upstream cause the temperature anisotropy. This electron temperature anisotropy can trigger the electron firehose instability (EFI), which excites oblique waves in the shock foot. Scattering of electrons by these waves enables multiple cycles of shock drift acceleration (SDA) in the preshock region, leading to the electron injection to diffusive shock acceleration (DSA). In the study, the kinetic properties of the EFI are examined by the linear stability analysis based on the kinetic Vlasov-Maxwell theory and then further investigated by 2D Particle-in-Cell (PIC) simulations, especially focusing on those in high-β (β~100) plasmas. We then discuss the basic properties of the firehose instability, and the implication of our work on electron acceleration in ICM shock.

• Smart Structures and Systems
• /
• 제12권1호
• /
• pp.1-22
• /
• 2013

Magnetic nanoparticle based bioseparation in microfluidics is a multiphysics phenomenon that involves interplay of various parameters. The ability to understand the dynamics of these parameters is a prerequisite for designing and developing more efficient magnetic cell/bio-particle separation systems. Therefore, in this work proof-of-concept experiments are combined with advanced numerical simulation to design and optimize the capturing process of magnetic nanoparticles responsible for efficient microfluidic bioseparation. A low cost generic microfluidic platform was developed using a novel micromolding method that can be done without a clean room techniques and at much lower cost and time. Parametric analysis using both experiments and theoretical predictions were performed. It was found that flow rate and magnetic field strength greatly influence the transport of magnetic nanoparticles in the microchannel and control the capturing efficiency. The results from mathematical model agree very well with experiments. The model further demonstrated that a 12% increase in capturing efficiency can be achieved by introducing of iron-grooved bar in the microfluidic setup that resulted in increase in magnetic field gradient. The numerical simulations were helpful in testing and optimizing key design parameters. Overall, this work demonstrated that a simple low cost experimental proof-of-concept setup can be synchronized with advanced numerical simulation not only to enhance the functional performance of magneto-fluidic capturing systems but also to efficiently design and develop microfluidic bioseparation systems for biomedical applications.

• 천문학회보
• /
• 제46권1호
• /
• pp.36.2-36.2
• /
• 2021

In the outskirts of galaxy clusters, weak shocks with M_s < ~3 appear as radio relics where the synchrotron radiation is emitted from cosmic-ray (CR) electrons. To understand the production of CR electrons through the so-called diffusive shock acceleration (DSA), the electron injection into the DSA process at shocks in the hot intracluster medium (ICM) has to be described. However, the injection remains as an unsolved, outstanding problem. To explore this problem, 2D Particle-in-Cell (PIC) simulations were performed. In this talk, we present the electron preacceleration mechanism mediated by multi-scale plasma waves in the shock transition zone. In particular, we find that the electron preacceleration is effective only in the supercritical shocks, which have the sonic Mach number M_s > M_crit ≈ 2.3 in the high-beta (β~100) plasma of the ICM, because the Alfven ion cyclotron instability operates and hence multi-scale plasma waves are induced only in such supercritical shocks. Our findings will help to understand the nature of radio relics in galaxy clusters.

• 한국지하수토양환경학회:학술대회논문집
• /
• 한국지하수토양환경학회 2001년도 추계학술발표회
• /
• pp.27-30
• /
• 2001

Flow and transport at fracture intersections, and their effects on network scale transport, are investigated in three-dimensional random fracture networks. Fracture intersection mixing rules complete mixing and streamline routing are defined in terms of fluxes normal to the intersection line between two fractures. By analyzing flow statistics and particle transfer probabilities distributed along fracture intersections, it is shown that for various network structures with power law size distributions of fractures, the choice of intersection mixing rule makes comparatively little difference in the overall simulated solute migration patterns. The occurrence and effects of local flows around an intersection (local flow cells) are emphasized. Transport simulations at fracture intersections indicate that local flow circulations can arise from variability within the hydraulic head distribution along intersections, and from the internal no flow condition along fracture boundaries. These local flow cells act as an effective mechanism to enhance the nondiffusive breakthrough tailing often observed in discrete fracture networks. It is shown that such non-Fickian (anomalous) solute transport can be accounted for by considering only advective transport, in the framework of a continuous time random walk model. To clarify the effect of forest environmental changes (forest type difference and clearcut) on water storage capacity in soil and stream flow, watershed had been investigated.

• International Journal of Fluid Machinery and Systems
• /
• 제2권1호
• /
• pp.31-39
• /
• 2009

The flow instability in a low specific speed mixed-flow pump, having a positive slope of head-flow characteristics was investigated. Based on the static pressure measurements, it was found that a rotating stall in the vaned diffuser occurs at about 65% flow rate of best efficiency point (BEP). A dynamic Particle Image Velocimetry (DPIV) measurement and the numerical simulations were conducted in order to investigate the flow fields. As a result, the diffuser rotating stall was simulated even by Computational Fluid Dynamics (CFD) and the calculated periodic flow patterns agree well with the measured ones by DPIV. It is clarified that a periodical large scaled backflow, generated at the leading edge of the suction surface of the diffuser vane, causes the instability. Furthermore, the growth of the strong vortex at the leading edge of the diffuser vane induces the strong backflow from the diffuser outlet to the inlet. The scale of one stall cell is covered over four-passages in total thirteen vane-passages.

• Nuclear Engineering and Technology
• /
• 제37권5호
• /
• pp.447-456
• /
• 2005

Laser driven accelerators promise to provide an alternative to conventional accelerator technology. They rely on the excitation of large amplitude density waves in a plasma by the photon pressure of an intense laser. The density oscillations in which electrons and ions are separated, result in extremely large longitudinal electric fields that can be several orders of magnitude larger than those that are used in today's radio-frequency accelerators. Whereas this principle had been demonstrated experimentally for nearly two decades, it was not until 2004 that the production of high quality electron beams around 100 MeV was demonstrated. Analysis, aided by particle-in-cell simulations, as well as experiments with various plasma lengths and densities, indicate that tailoring the length of the accelerator, together with loading of the accelerating structure with beam, are the keys to production of mono-energetic electron beams. Increasing the energy towards a GeV and beyond will require reducing the plasma density and design criteria are discussed for an optimized accelerator module. The current progress and future directions are summarized through comparison with conventional accelerators, highlighting the unique short and long term prospects for intense radiation sources and high energy accelerators based on laser-drivenplasma accelerators.