• 한국해안해양공학회지
• /
• 제3권3호
• /
• pp.137-151
• /
• 1991

영역별로 속도포텐셜을 정의하고 부방파제 측면의 가상 경계에서 유체운동의 질량속과 운동량속의 연속조건으로부터 미지계수를 구하는 속도포텐셜 접속법을 이용하여 천해역에 기인장의 상태로 계유된 직사각형 부방파제의 운동을 해석하였다. 본 방법에서는 정도 등(1972)이 수중계유물제 운동해석시 사용한 방법을 부유체 계유로 수정하여 전체유체영역을 입사파 영역, 투과파 영역 및 부방파제 저면 영역의 3개 부영역으로 구분하였다. 또한, 전도 등(1972)이 무시한 기인장된 계유라인에 의한 복원력을 선형 스프링으로 모형화하였다. 본 방법을 수리모형실험 결과와 비교ㆍ분석한 결과 정확도가 향상되었으며, 이를 이용하여 부방파제의 흘수 및 폭의 변화에 따른 발파 성능을 조사하였다.

• 대한기계학회논문집B
• /
• 제35권10호
• /
• pp.1089-1095
• /
• 2011

열전달 향상을 위한 방법으로 많이 사용되고 있는 마이크로 핀을 포함한 표면 위에서의 핵비등을 액상과 기상에서 질량 및 운동량, 에너지에 대한 지배 방정식을 풀어 수치해석을 수행하였다. 핵비등에서의 기포거동을 계산하기 위해 sharp-interface 레벨셋(level-set) 방법을 상변화 효과와 핀과 캐비티와 같은 잠긴 고체에서의 점착 조건 및 접촉각, 마이크로 액체층에서의 증발 열유속을 포함하도록 수정하였다. 핀과 캐비티를 포함한 표면에서의 기포 생성, 성장, 이탈에 대한 해석을 통하여 핀-캐비티 배열, 핀-핀 간격이 핵비등에서의 기포거동에 중요한 역할을 하는 것을 확인하였다.

• 한국해안·해양공학회논문집
• /
• 제23권4호
• /
• pp.276-284
• /
• 2011

파랑 잉여응력의 영향이 자동적으로 고려되므로 파랑으로부터 발생되는 흐름을 수치모의할 수 있는 수평점성 및 난류항이 포함된 Boussinesq 방정식 모형인 FUNWAVE를 이용하여, 다방향 불규칙파 조건으로 해운대 해안에서 발생하는 이안류를 수치모의하였다. 수치모의는 다방향 불규칙파의 전파양상과 지형에 의한 비선형 파랑변형을 잘 보여주고 있으며, 이러한 파랑변형과 해운대 지형특성이 반영되어 시간에 따라 발달하는 파랑유도 연안흐름 양상을 잘 보여준다. 수치모의 결과로 부터 이안류는 연안방향으로 상대적으로 수심이 깊은 곳에서 그리고 파고가 낮은 곳에서 이안류가 돌발적으로 발생 혹은 증폭될 수 있음을 확인하였다.

• 한국추진공학회:학술대회논문집
• /
• 한국추진공학회 2017년도 제48회 춘계학술대회논문집
• /
• pp.787-790
• /
• 2017

본 연구는 이원추진제 추력기(thruster)에 사용되는 메탄-산소 연소특성규명의 선행연구로서 모델연소실 내 전단동축형인젝터를 통해 분사된 기체메탄-기체산소 비예혼합화염의 연소안정한계 및 화염형상을 도출하기 위한 실험적 연구가 수행되었다. DSLR 카메라를 이용하여 화염 직접이미지(direct image)를 촬영하였고, 이미지 후처리(post-processing)를 통해 연소특성파악 및 화염길이 정량화를 수행하였다. 그 결과, 산화제 레이놀즈 수($Re_o$)가 증가함에 따라 이론반응비(stoichiometric ratio)에서 안정된 화염이 발생하였고, 동일 인젝터직경 조건에서 난류화염의 길이가 늘어남을 확인하였다.

• 한국항공우주학회지
• /
• 제31권1호
• /
• pp.97-103
• /
• 2003

연소기 내부에서 연소효율의 증대목적으로 사용되어지고 있는 후향계단후류의 측면제트 분사가 있는 초음속 난류유동장에 대한 수치적 연구를 수행하였다. 본 연구에서는 2차원 Navier_Stokes 방정식과 k-$\varepsilon$ 난류모델식을 고해상도 풍상차분법인 TVD 기법을 사용하여 비정상 유동장에 대한 수치계산을 수행하였다. 유동특성은 제트분사에 의한 마하면이 장애물 역할을 하여 후향계단과 마하면사이에 공동유동(cavity flow)과 비슷한 유동현상이 나타난다. 따라서, 공동 유동에서 발생하는 유동의 진동특성이 나타날 것을 예측할 수 있다. 수치계산 결과에서도 모맨텀플럭스와 분사위치에 따라, 공동유동에서 발생하는 진동특성을 갖는 비정상 유동 특성이 나타나는 것을 볼 수 있었다.

• 대한기계학회논문집
• /
• 제11권2호
• /
• pp.262-270
• /
• 1987

본 논문에서는 Eulerian 좌표축에서 온도 경계층과 혼합된 비정상 경계층 유 동을 효과적으로 다룰 수 있는 좌표변환 방법과 수치해석 방법을 도입하였다. 이전 의 수치적 방법들은 축차적이며, 뒷정체점 부근의 경계층 두께가 시간과 더불어 지수 적으로 증가하므로, 격자점의 수를 경계층 두께의 증가에 따라 상당히 늘려 주어야 된 다. 그러나 여기에서는, 격자점의 수를 무리하게 늘릴 필요가 없으며 일반적인 비선 형 차분방정식을 정확도를 떨어뜨리지 않고 선형화시킴으로써 비축차적인 수치해들이 허용된다. 이런 선형화 방법은 Beam & Warming에 의해 최초로 압축성 Navier-Stokes 방정식에 사용되었고, Orlandi & Ferziger, Kim & Chang에 의해 경계층 유동에 확장되 었다.또한 뒷정체점 부근에서 경계층 두께의 증가로 인한, 격자점 증가의 필요를 피하기 위하여, 몇 가지 서로 다른 종류의 변환변수들을 시간과 공간에 따라 선별적으 로 사용하여 수치적인 경계층 두께가 거의 일정하도록 만들었다. 이와 같은 변환변 수들은 또한 쌍-포물선형인 현 지배방정식의 초기조건들을 쉽게 구할 수 있도록 허용 해 준다.

• Nuclear Engineering and Technology
• /
• 제53권10호
• /
• pp.3344-3351
• /
• 2021

The Rare Isotope Science Project was launched in 2011 in Korea toward constructing the Rare isotope Accelerator complex for ON line experiments (RAON). RAON will house several experimental systems, including the Muon Spin Rotation/Relaxation/Resonance (μSR) facility in High Energy Experimental Building B. This facility will use 600-MeV protons with a maximum current of 660 pμA and beam power of 400 kW. The key μSR features will facilitate projects related to condensed-matter and nuclear physics. Typical experiments require a few million surface muons fully spin-polarized opposite to their momentum for application to small samples. Here, we describe the design of a muon transport beam line for delivering the requisite muon numbers and the electromagnetic-component specifications in the μSR facility. We determine the beam-line configuration via beam-optics calculations and the transmission efficiency via single-particle tracking simulations. The electromagnet properties, including fringe field effects, are applied for each component in the calculations. The designed surface-muon beamline is 17.3 m long, consisting of 2 solenoids, 2 dipoles affording 70° deflection, 9 quadrupoles, and a Wien filter to eliminate contaminant positrons. The average incident-muon flux and spin rotation angle are estimated as 5.2 × 10⁶ μ⁺/s and 45°, respectively.

• 천문학회보
• /
• 제40권1호
• /
• pp.73.2-74
• /
• 2015

The gas cloud G2 falling toward Sagittarius A* (Sgr A*), the supermassive black hole at the center of the Milky Way, is supposed to provide valuable information on the physics of accretion flows and the environment of the black hole. We observed Sgr A* with four European stations of the Global Millimeter Very Long Baseline Interferometry Array (GMVA) at 86 GHz on 1 October 2013 when parts of G2 had already passed the pericenter. We searched for possible transient asymmetric structure - such as jets or winds from hot accretion flows - around Sgr A* caused by accretion of material from G2. The interferometric closure phases (which are zero if the spatial brightness distribution of the target is symmetric, and deviate from zero otherwise) remained zero within errors during the observation time. We thus conclude that Sgr A* did not show significant asymmetric (in the observer frame) outflows in late 2013. Using simulations, we constrain the size of the outflows that we could have missed to ${\approx}2.5$ mas along the major axis, ${\approx}0.4$ mas along the minor axis of the beam, corresponding to approximately 232 and 35 Schwarzschild radii, respectively; we thus probe spatial scales on which the jets of radio galaxies are suspected to convert magnetic into kinetic energy. As probably less than 0.2 Jy of the flux from Sgr A* can be attributed to accretion from G2, one finds an effective accretion rate ${\eta}M{\leq}1.5{\times}10^9kg/s{\approx}7.7{\times}10^{-9}M_{earth}/yr$ for material from G2. Exploiting the kinetic jet power-accretion power relation of radio galaxies, one finds that the rate of accretion of matter that ends up in jets is limited to $M{\leq}10^{17}kg/s{\approx}0.5M_{Earth}/yr$ less than about 20% of the mass of G2. Accordingly, G2 appears to be largely stable against loss of angular momentum and subsequent (partial) accretion at least on time scales ${\leq}1$ year.

• 한국해안·해양공학회논문집
• /
• 제29권3호
• /
• pp.139-146
• /
• 2017

해안선 변동을 제어하기 위해 설치되는 잠제의 효과를 이해하기 위해 해안선과 평행한 단일 잠제 주변 연안의 흐름변화 특성을 수리실험을 통해 연구하였다. 잠제와 입사파 파라미터에 따른 잠제 후면의 흐름패턴을 관측하기 위해 LSPIV(Large-Scale Particle Image Velocimetry) 이미지 분석 기법을 적용하였고, 잠제 주변의 역학적 특성을 파악하기 위해 평균수면 및 파고분포를 관측하였다. 수리실험을 통해 흐름패턴은 잠제에 의한 파고(파랑응력) 분포변화와 그에 따라 발달하는 잠제 후면과 측면 쇄파유도류 및 평균수면상승효과의 상호작용에 의해 변화됨을 알 수 있었다. 잠제 주변의 분기 또는 수렴하는 흐름패턴에 따라 각각 해안선 변동을 침식 또는 퇴적반응으로 구분하는 Ranasinghe et al.(2010)의 경험식과 수리실험 결과를 비교하여, 전반적으로 경험식에 부합하는 흐름패턴을 확인할 수 있었다. 그러나 일부 실험결과에서 기존의 경험식을 적용하여 구분하기 어려운 침식과 퇴적이 혼재된 반응을 유도하는 흐름패턴도 관찰할 수 있었다.

• 한국진공학회:학술대회논문집
• /
• 한국진공학회 2014년도 제46회 동계 정기학술대회 초록집
• /
• pp.221-221
• /
• 2014

The present work proposes an improved numerical simulator for design and modification of large area capacitively coupled plasma (CCP) processing chamber. CCP, as notoriously well-known, demands the tremendously huge computational cost for carrying out transient analyses in realistic multi-dimensional models, because electron dissociations take place in a much smaller time scale (${\Delta}t{\approx}10-8{\sim}10-10$) than time scale of those happened between neutrals (${\Delta}t{\approx}10-1{\sim}10-3$), due to the rf drive frequencies of external electric field. And also, for spatial discretization of electron flux (Je), exponential scheme such as Scharfetter-Gummel method needs to be used in order to alleviate the numerical stiffness and resolve exponential change of spatial distribution of electron temperature (Te) and electron number density (Ne) in the vicinity of electrodes. Due to such computational intractability, it is prohibited to simulate CCP deposition in a three-dimension within acceptable calculation runtimes (<24 h). Under the situation where process conditions require thickness non-uniformity below 5%, however, detailed flow features of reactive gases induced from three-dimensional geometric effects such as gas distribution through the perforated plates (showerhead) should be considered. Without considering plasma chemistry, we therefore simulated flow, temperature and species fields in three-dimensional geometry first, and then, based on that data, boundary conditions of two-dimensional plasma discharge model are set. In the particular case of SiH4-NH3-N2-He CCP discharge to produce deposition of SiNxHy thin film, a cylindrical showerhead electrode reactor was studied by numerical modeling of mass, momentum and energy transports for charged particles in an axi-symmetric geometry. By solving transport equations of electron and radicals simultaneously, we observed that the way how source gases are consumed in the non-isothermal flow field and such consequences on active species production were outlined as playing the leading parts in the processes. As an example of application of the model for the prediction of the deposited thickness uniformity in a 300 mm wafer plasma processing chamber, the results were compared with the experimentally measured deposition profiles along the radius of the wafer varying inter-electrode gap. The simulation results were in good agreement with experimental data.