• 제목/요약/키워드: Mechanism simulation

검색결과 2,645건 처리시간 0.035초

PSCAD/EMTDC를 이용한 ESS의 누설전류 모델링에 관한 연구 (A Study on Modeling of Leakage Current in ESS Using PSCAD/EMTDC)

  • 김지명;태동현;이일무;임건표;노대석
    • 한국산학기술학회논문지
    • /
    • 제22권2호
    • /
    • pp.810-818
    • /
    • 2021
  • ESS의 누설전류는 PCS(Power Control System)측 누설전류와 계통불평형 전류로 인한 누설전류로 구분되는데, PCS측의 누설전류는 정상 상태 운전 시, IGBT(Insulated Gate Bipolar Transistor) 스위칭의 전압 변화량과 IGBT와 방열판 사이에 존재하는 기생 커패시턴스에 의해 발생한다. 또한, 계통불평형 전류에 의한 누설전류는 불평형 부하로 인해 발생한 불평형 전류가 Yg-∆ 결선방식의 3각 철심이 적용된 태양광전원 연계형 변압기의 중성선을 통해 ESS로 유입된다. 따라서, 본 논문에서는 방열판 유도공식을 통해 산정한 기생 커패시턴스에 의하여 PCS측의 누설전류 발생 메커니즘을 제시하고 또한, 계통불평형에 의한 ESS측의 누설전류 발생 메커니즘을 제안한다. 이를 바탕으로, 배전계통 상용해석 프로그램인 PSCAD/EMTDC를 이용하여 배터리부, PCS부, AC전원부로 이루어진 PCS측의 누설전류 발생 메커니즘과 배전 계통부, 불평형 부하부, ESS부로 이루어진 계통불평형에 의한 ESS측의 누설전류 발생 메커니즘을 모델링하고, 누설전류의 특성을 평가한다. 상기의 모델링을 바탕으로 시뮬레이션을 수행한 결과, 외함의 저항과 접지저항의 크기에 따라 PCS측의 누설전류는 7[mA]에서 34[mA]로, 계통불평형에 의한 배터리 외함으로 흐르는 누설전류는 3.96[mA]에서 10.76[mA]로 증가하여 배터리측에 큰 영향을 미침을 알 수 있었다.

유역모델을 이용한 농업용수 신속회귀수량 산정 연구 (A study on estimating the quick return flow from irrigation canal of agricultural water using watershed model)

  • 이지완;정충길;김다예;맹승진;정현식;조영식;김성준
    • 한국수자원학회논문집
    • /
    • 제55권5호
    • /
    • pp.321-331
    • /
    • 2022
  • 본 연구에서는 유역단위의 물수지 분석과 농업용수의 수문학적 매커니즘을 모의 할 수 있는 유역 모델링 방법을 이용하여 회귀수량 산정기법을 제시하고자 하였다. SWAT 모델을 이용하여 영산강수계 대표적인 농업지역인 만봉천 표준유역 (97.34 km2)에 대해 담수 논 모의가 고려된 유역물수지 분석을 실시하였다. 회귀수량 산정에 앞서, 나주 유량관측소의 일 유량 자료를 이용하여 SWAT을 검·보정하였다. R2, Nash-Sutcliffe Efficiency (NSE), Root-Mean-Square Error (RMSE)는 각각 0.73, 0.70, 0.64 mm/day으로 분석되었다. 3년 동안(2015~2017) 모의 결과를 토대로 관개기간(4/1~9/30)에 대한 신속 회귀수량과 공급량 대비 회귀율을 산정하였고, 평균 53.4%로 분석되었다. 본 연구에서 제시한 유역 회귀수량 모델링 기법은 향후 합리적인 유역물관리를 위한 최적 농업용수 공급방안에 대한 기초자료 구축에 활용될 수 있다.

Effects of particle size and loading rate on the tensile failure of asphalt specimens based on a direct tensile test and particle flow code simulation

  • Q. Wang;D.C. Wang;J.W. Fu;Vahab Sarfarazi;Hadi Haeri;C.L. Guo;L.J. Sun;Mohammad Fatehi Marji
    • Structural Engineering and Mechanics
    • /
    • 제86권5호
    • /
    • pp.607-619
    • /
    • 2023
  • This study, it was tried to evaluate the asphalt behavior under tensile loading conditions through indirect Brazilian and direct tensile tests, experimentally and numerically. This paper is important from two points of view. The first one, a new test method was developed for the determination of the direct tensile strength of asphalt and its difference was obtained from the indirect test method. The second one, the effects of particle size and loading rate have been cleared on the tensile fracture mechanism. The experimental direct tensile strength of the asphalt specimens was measured in the laboratory using the compression-to-tensile load converting (CTLC) device. Some special types of asphalt specimens were prepared in the form of slabs with a central hole. The CTLC device is then equipped with this specimen and placed in the universal testing machine. Then, the direct tensile strength of asphalt specimens with different sizes of ingredients can be measured at different loading rates in the laboratory. The particle flow code (PFC) was used to numerically simulate the direct tensile strength test of asphalt samples. This numerical modeling technique is based on the versatile discrete element method (DEM). Three different particle diameters were chosen and were tested under three different loading rates. The results show that when the loading rate was 0.016 mm/sec, two tensile cracks were initiated from the left and right of the hole and propagated perpendicular to the loading axis till coalescence to the model boundary. When the loading rate was 0.032 mm/sec, two tensile cracks were initiated from the left and right of the hole and propagated perpendicular to the loading axis. The branching occurs in these cracks. This shows that the crack propagation is under quasi-static conditions. When the loading rate was 0.064 mm/sec, mixed tensile and shear cracks were initiated below the loading walls and branching occurred in these cracks. This shows that the crack propagation is under dynamic conditions. The loading rate increases and the tensile strength increases. Because all defects mobilized under a low loading rate and this led to decreasing the tensile strength. The experimental results for the direct tensile strengths of asphalt specimens of different ingredients were in good accordance with their corresponding results approximated by DEM software.

방사선 조사시 산소가 세포에 미치는 영향의 이론적 분석 (A Theoretical Study for Estimation of Oxygen Effect in Radiation Therapy)

  • Rena J. Lee;HyunSuk Suh
    • 한국의학물리학회지:의학물리
    • /
    • 제11권2호
    • /
    • pp.157-165
    • /
    • 2000
  • 목적: 방사선조사에 의해 형성된 DNA 손상정도가 산소에 의해 변화되는 추세을 이론적 모델을 이용하여 평가 및 분석하였다. 방법 및 대상: Water radiolysis(물의 방사성분해)시 생성되는 유리기 들간의 반응, 손상된 DNA의 형성, 손상의 복구, 및 산소에 의한 손상의 고정들을 시간 미분 방정식을 이용하여 표현하였다. 문헌에 나와 있는 반응상수(rate constants)들은 그대로 사용되었고 알려지지 않은 상수들은 실험에서 얻어진 데이터에서 얻은 곡선을 대입하여 얻었다. 화학반응상수 변화와 산소농도변화에 따른 세포 생존도를 구하였다. 모델을 통해 얻어진 세포 생존도와 방사선량의 상관관계를 세포 생존곡선으로 표시하였다. 산소에 의한 손상의 fixation(고착화)과정과 산소농도가 세포생존에 미치는 영향을 분석하여 보았다. 산소가 세포생존에 미치는 정도를 정량화 하기 위하여 산소증가효과비(Oxygen Ehhancement Ratio, OER)를 계산하여 실험치와 이론치를 비교하였다. 결과: 산소에 대한 민감성 연구에서 DNA생존도는 산소의 농도와 화학반응속도상수에 영향을 받았다. 산소의 농도가 0에서 1.0 $\times$ $10^{7}$ 으로 변화할 때 세포 생존도에는 변화가 없었다. 그러나 1.0 x $10^{7}$ 에서 1.0 $\times$ $10^{10}$ 변화할 경우 세포의 생존률이 감소하였다. 모의연구에서 얻은 OER치는 10% 세포생존시에는 2.32, 45% 세포생존시에는 1.9 이었다. 결론: 산소의 감수성 연구에서 보여주듯이 방정식을 이용하여 얻어진 새포생존곡선은 실험에서 얻어진 세포생존곡선을 재연할 수 있었다 또한 방사선조사시 생성되는 손상된 세포의 양은 산소 반응상수가 가장 크고 산소농도가 증가되는 경우에 증가됨을 알수있었다. 모의연구에서 얻은 산소증가효과비는 세포생존이 low level인 경우 높은 일치성을 보여 주었다 따라서 본 연구는 세포살해시 산소의 효과를 semi-empirical 모델을 사용하여 예측할수 있다는 것을 보여주고 있다.

  • PDF

유무선랜 환경에서 TCP Flow의 성능향상을 위한 MAC 계층과 TCP 계층의 연동기법 (Interaction Between TCP and MAC-layer to Improve TCP Flow Performance over WLANs)

  • 김재훈;정광수
    • 한국정보과학회논문지:정보통신
    • /
    • 제35권2호
    • /
    • pp.99-111
    • /
    • 2008
  • 최근 무선 랜은 SOHO (Small Office Home Office) 및 Hot Spot과 같은 환경에서 공간의 제약에 구애받지 않고, 인터넷에 접속할 수 있는 기술로서 사용자의 요구가 크게 증가하였다. 하지만, 무선 랜 환경에서의 통신은 유선망과 달리 불안정한 무선 채널의 특성으로 인해 연집적인 패킷 손실이 발생하여 통신상의 제약이 많은 특징을 가진다. 연집적인 패킷 손실은 AP(Access Point) 와 무선 단말의 거리가 증가하거나, AP와 무선 단말사이에 장애물 등이 일시적으로 지나갈 때 주로 발생하는 현상이다 결국, 현재 인터넷상에서 가장 광범위하게 사용되고 있는 무선 랜 기술인 IEEE 802.11은 이러한 특성으로 인해 사용자의 요구에 만족할만한 전송 성능을 나타내지 못하며, 특히 전송 계층에 TCP가 사용될 경우 불필요한 혼잡 제어 기법을 사용하게 함으로써 심각한 성능저하를 야기한다. 이러한 무선 랜 환경의 문제점을 해결하기 위해 MAC-layer LDA(Loss Differentiation Algorithm)가 제안되었다. MAC-layer LDA는 MAC 계풍의 Retry limit을 기반으로 CRD(Consecutive Retry Duration)를 무선 구간의 연집된 패킷손실 기간 이상 증가시켜, TCP의 불필요한 Timeout 발생 이전에 손실된 패킷을 효율적으로 복구하는 기법이다. 하지만, MAC-layer LDA 기법은 한정된 Retry limit의 증가로 인해 CRD가 연집된 패킷 손실 구간 보다 적은 경우가 발생하여 심각한 전송성능 저하를 가져온다. 또한, CRD의 증가는 무선 구간의 패킷 처리 시간을 증가시켜 대역폭과 무선 단말의 한정된 에너지 자원을 불필요하게 낭비하는 문제를 초래한다. 본 논문에서는 이러한 문제점을 개선하기 위해 Cross-layer 기법을 적용한 재전송 기법인 BLD(Burst Loss Detection) 모듈을 제안한다. BLD 모듈의 알고리즘은 현재 무선 랜 환경에서 가장 널리 사용되는 IEEE 802.11 MAC 프로토콜 기반의 재전송 기법으로서, MAC 계층과 TCP에서 사용되는 재전송 기법의 효율적인 연동을 통해 손실된 패킷을 복구한다. ns-2(Network Simulator) 시뮬레이터를 이용한 실험을 통해 BLD 모듈은 무선 구간의 연집적인 패킷 손실에 대해 효율적인 보상을 수행하여 전송 성능과 에너지 효율성을 향상시킬 수 있음을 확인하였다.

3차원 게임에서 객체들의 상호 작용을 디자인하기 위한 제어 기법 (A Control Method for designing Object Interactions in 3D Game)

  • 김기현;김상욱
    • 한국정보과학회논문지:컴퓨팅의 실제 및 레터
    • /
    • 제9권3호
    • /
    • pp.322-331
    • /
    • 2003
  • 3차원 게임은 게임 시나리오의 다양한 요소에 의한 복잡도가 증가함에 따라 게임 객체들의 상호 관계를 제어하기 위한 문제점을 가진다. 그러므로, 게임 시스템은 각 게임 객체들의 응답을 조정하는 방법의 필요성을 가진다. 또한, 게임 시나리오의 결과에 따라 게임 객체들의 행동 애니메이션을 제어하기 위한 개념들도 필요하다. 사실적 게임 시뮬레이션을 생성하기 위해 시스템은 게임 객체들의 상호작용을 디자인 할 수 있는 구조를 포함해야 한다. 본 논문에서 게임 시나리오상에 게임 객체들의 상호작용 설계를 위해 동적 제어를 디자인하는 기법을 소개한다. 이 방법을 위해 특정 규칙을 이용한 의사결정이 가능한 지능적 에이전트 기반 구조로써 게임 에이전트 시스템을 제안한다. 게임 에이전트 시스템은 환경 데이터 처리, 게임 객체 시뮬레이션, 게임 객체들간의 상호작용 제어, 게임 객체들의 다양한 상호 관계를 정의할수 있는 시각 저작 인터페이스를 제공하기 위해 이용되어진다. 이들 기술들은 게임 객체의 자율성과 연관된 충돌 회피 기법 등을 처리한다. 또한, 장면의 변경으로부터 게임 객체들의 일관된 의사 결정력을 가능하게 한다. 본 논문에서는 규칙기반 행동 제어가 게임 객체의 시뮬레이션을 안내하기 위해 디자인되어졌다. 시각적 요소들로 구성된 에이전트 상태 결정 네트워크는 정보전달과 게임 객체들 사이의 현상태를 추론할 수 있다. 이들 기법들은 실시간으로 게임 객체들간의 동작 상태 변이를 체크하고 모니터링 할 수 있다. 마지막으로 간단한 사례 연구 예와 함께 제어 기법의 타당성을 제시한다.

실규모 현장시험 및 유한요소해석을 통한 강관매입말뚝의 공학적 거동에 대한 연구 (A Study on the Engineering Behaviour of Prebored and Precast Steel Pipe Piles from Full-Scale Field Tests and Finite Element Analysis)

  • 김정섭;정경자;정상섬;전영진;이철주
    • 한국지반환경공학회 논문집
    • /
    • 제19권4호
    • /
    • pp.5-16
    • /
    • 2018
  • 본 연구에서는 현장에서 실규모로 시공된 8본의 강관매입말뚝의 공학적 거동을 정재하시험, 동재하시험(EOID 및 restrike 시험) 및 Class-A 및 C1 type의 수치해석을 수행하여 상세히 고찰하였다. 이를 통해 말뚝의 하중-침하, 설계지지력 및 전단응력 전이 특성 등을 분석하였다. 정재하시험을 통해 평가된 설계지지력에 비해 동재하시험 및 수치해석은 이를 매우 상이하게 평가하는 것으로 분석되었으며, 이 가운데 restrike 시험이 전반적으로 최상의 결과를 도출하였고, 수치해석 결과의 신뢰성은 동재하시험에 비해 낮은 것으로 분석되었다. 각 시험에서 산정된 설계지지력을 정재하시험의 결과와 비교할 때 EOID는 정재하시험값의 20.0%-181.0%(평균: 69.3%)의 범위를 보였고, restrike 시험의 경우 48.2%-181.1%(평균: 92.1%)의 범위를 보였다. 또한 Class-A type에서는 37.1-210.5%(평균: 121.2%)로 가장 큰 분산을 보였다. EOID 시험에서는 선단지지력이 말뚝의 전체지지력의 대부분을 차지하는데 비해, restrike 시험에서는 주면마찰력 및 선단지지력이 비슷한 정도로 발현되었다. 이때 restrike 시험에서 측정된 선단지지력은 EOID 시험에서 측정된 크기보다 작은 것으로 분석되었다. 즉 restrike 시험에서 타격에너지가 충분하지 않은 경우 말뚝의 선단지지력이 과소평가될 가능성이 있는 것으로 나타났다. 말뚝의 축력분포로부터 측정된 전단응력은 말뚝의 심도가 깊어질수록 증가하는 양상을 보였으며, 수치해석 결과는 정재하시험과는 상당한 정도의 차이를 보였다. 말뚝선단 인근에 슬라임이 존재하는 경우 말뚝의 거동에 매우 큰 영향을 주는 것으로 분석되었는데, 슬라임의 탄성계수가 작을수록 그리고 슬라임의 두께가 두꺼울수록 말뚝의 침하량이 큰 폭으로 증가하는 것을 확인하였다.

충돌하중을 받는 이방향 비부착 프리스트레스트 콘크리트 패널부재의 충돌저항성능에 대한 실험적 거동 평가 (Experimental Evaluation of Bi-directionally Unbonded Prestressed Concrete Panel Impact-Resistance Behavior under Impact Loading)

  • 이나현;이상원;이승재;김장호
    • 콘크리트학회논문집
    • /
    • 제25권5호
    • /
    • pp.485-496
    • /
    • 2013
  • 2001년 9.11테러로 인한 미국 세계무역센터 및 미국 국방성 펜타곤이 붕괴된 이후, 전 세계적으로 충돌, 폭발 등의 극한하중에 의한 테러가 빈번하게 발생하고 있으며, 극한하중에 의한 구조물의 거동에 대한 사회적 불안감은 더욱 증가되고 있다. 그러나 지금까지의 사회기반시설구조물에는 폭발 및 충돌 등과 같은 극한하중에 대한 방호 및 방재개념을 설계에 고려하지 못하고 있는 실정이며, 원전격납건물, 가스탱크, 교량, 터널 등에 널리 사용되는 프리스트레스트 콘크리트 구조물에 대한 극한하중 연구는 전 세계적으로 미흡하다. 충돌과 같은 극한하중은 집약된 에너지의 급작스런 방출로 인한 높은 충돌압력을 형성하므로, 극한하중의 특성 및 전파 메커니즘을 이해하는 것이 필요하다. 그러므로 이 연구에서는 이방향 비부착 프리스트레스트 콘크리트 패널의 충돌저항성능을 비교검토하기 위하여, $1400mm{\times}1000mm{\times}300mm$의 철근콘크리트(RC), 프리스트레스 텐던으로만 보강된 콘크리트(PS), 프리스트레스 텐던과 철근으로 보강된 콘크리트(PSR, 일반적인 PSC) 시편을 제작하였다. 실험 조건에 맞춰 14 kN의 추를 10 m, 5 m, 4 m 높이에서 낙하하는 예비 실험과 3.5 m 높이의 본 실험으로 구성하여 충돌하중에 대한 프리스트레스트 콘크리트 구조물의 실험적 평가를 수행하였다. 또한, 충돌실험을 위한 기본적인 실험 구성 및 계측시스템을 구축하였다. 충돌 저항성능은 균열형상, 손상면적, 에너지 흡수, 처짐, 변형률, 가속도 등의 충돌에 의한 계측데이터를 이용한 거동분석 뿐만 아니라, 충돌 후 잔류구조성능 실험을 수행하여 이방향 비부착 프리스트레스트 콘크리트 패널의 충돌저항성능을 검토하였다. 본 실험은 향후 국내외 프리스트레스트 콘크리트에 대한 충돌 방호설계 및 충돌해석 등 관련 연구분야의 기초자료가 될 것이라고 판단되는 바이다.

보안 서비스를 고려한 이동 에이전트 모델과 클라이언트-서버 모델의 성능 비교 (A Performance Comparison of the Mobile Agent Model with the Client-Server Model under Security Conditions)

  • 한승완;정기문;박승배;임형석
    • 한국정보과학회논문지:정보통신
    • /
    • 제29권3호
    • /
    • pp.286-298
    • /
    • 2002
  • 분산 컴퓨팅 환경에서 프로세스 사이의 상호 협력을 위한 통신으로 인격 프로시져 호출이 전통적으로 사용되고 있다. 분산 응용이 더욱 복잡해짐에 따라 최근 이동 에이전트 패러다임이 등장하였다. 이처럼 다양한 상호 협력을 위한 통신 패러다임이 등장함에 따라 각 패러다임의 성능에 대한 평가와 비교 연구가 이루어지고 있다. 그러나 기존의 연구에서 성능 평가를 위해 사용한 성능 모델들은 보안 서비스를 위한 평가 요소를 고려하고 있지 않기 때문에 실제 분산 환경을 제대로 반영하지 못한다. 분산 환경은 개방되어 있으므로 정보의 노출이나 도청과 같은 공격에 있어서 상당히 취약하다. 이러한 분산 환경에서 안전하게 작업을 수행하기 위해서는 여러 가지 공격으로부터 응용 프로그램이나 정보를 보호하기 위한 보안 서비스가 고려되어야 한다. 본 논문에서는 상호 협력을 위한 통신 패러다임 중 인격 프로시져 호출과 이동 에이전트의 성능을 평가하고 비교한다. 분산 응용 프로그램을 안전하게 수행하기 위해 고려해야 하는 보안 서비스에 관하여 알아보고, 이러한 보안 서비스를 적용한 새로운 성능 모델을 제시한다. N개의 데이타베이스 서버에서 사용자가 필요한 정보를 검색하는 작업을 Petri Net으로 모델링하고, 각 파라미터에 수치 값을 할당해서 수행속도를 측정하여 두 패러다임의 성능을 비교한다. 본 논문에서 안전한 통신을 위하여 보안 서비스를 적용한 두 성능 모델의 비교 결과는 다음과 같다. 원격 프로시져 호출은 연산 비용이 높은 암호화 메커니즘을 포함하는 통신 횟수와 통신량이 많기 때문에 실행 시간이 급격하게 증가하지만, 이동 에이전트 패러다임은 통신 횟수와 통신량을 줄인 수 있으므로 실행시간이 완만하게 증가하는 것을 살펴볼 수 있다.멀티미디어 제작환경을 구축하는 것이 디지털 방송 시대의 방송 사업자에게 가장 중요한 과제중의 하나가 되었다. 멀티미디어 제작환경을 구축함으로써 영상, 음성 및 다양한 부가 데이터를 포함하는 멀티미디어 프로그램을 편리하게 제작할 수 있으며, 데이터베이스로부터 필요한 영상 이미지를 자유롭게 합성, 조작하는 등, 매우 다양하고 편리한 제작기법을 활용할 수 있다. 또한 멀티미디어를 응용한 제작 분야로서 컴퓨터 그래픽스 기술은 방송의 사전제작에 커다란 기여를 하고 있으며, 이미 선거방송을 비롯한 여러가지 프로그램은 가상스튜디오와 가상캐릭터 기술을 활용하여 제작하고 있다. 방송사업자는 이러한 멀티미디어 제작시스템을 근간으로 영상검색, 영상 합성, 스크립트 편집, 가상현실 응용 등 고도의 제작 기법을 활용함으로써 사용자 친화성, 다이나믹한 표현, 실시간, 대화성을 특징으로 하는 다양한 멀티미디어 서비스를 시청자에게 제공할 수 있을 것이다.is. Using the results, we can identify dominant overestimation sources that should be analyzed more accurately to get tighter WCET estimations. To make our method independent of any existing analysis techniques, we use simulation based methodology. We have implemented a MIPS R3000 simulator equipped with several switches, each of which determines the accuracy level of the timing analysis

공동이 있는 수직 분사 초음속 연소기 내의 불안정 연소유동 해석 (Numerical Analysis of Unstable Combustion Flows in Normal Injection Supersonic Combustor with a Cavity)

  • Jeong-Yeol Choi;Vigor Yang
    • 한국추진공학회:학술대회논문집
    • /
    • 한국추진공학회 2003년도 제20회 춘계학술대회 논문집
    • /
    • pp.91-93
    • /
    • 2003
  • A comprehensive numerical study is carried out to investigate for the understanding of the flow evolution and flame development in a supersonic combustor with normal injection of ncumally injecting hydrogen in airsupersonic flows. The formulation treats the complete conservation equations of mass, momentum, energy, and species concentration for a multi-component chemically reacting system. For the numerical simulation of supersonic combustion, multi-species Navier-Stokes equations and detailed chemistry of H2-Air is considered. It also accommodates a finite-rate chemical kinetics mechanism of hydrogen-air combustion GRI-Mech. 2.11[1], which consists of nine species and twenty-five reaction steps. Turbulence closure is achieved by means of a k-two-equation model (2). The governing equations are spatially discretized using a finite-volume approach, and temporally integrated by means of a second-order accurate implicit scheme (3-5).The supersonic combustor consists of a flat channel of 10 cm height and a fuel-injection slit of 0.1 cm width located at 10 cm downstream of the inlet. A cavity of 5 cm height and 20 cm width is installed at 15 cm downstream of the injection slit. A total of 936160 grids are used for the main-combustor flow passage, and 159161 grids for the cavity. The grids are clustered in the flow direction near the fuel injector and cavity, as well as in the vertical direction near the bottom wall. The no-slip and adiabatic conditions are assumed throughout the entire wall boundary. As a specific example, the inflow Mach number is assumed to be 3, and the temperature and pressure are 600 K and 0.1 MPa, respectively. Gaseous hydrogen at a temperature of 151.5 K is injected normal to the wall from a choked injector.A series of calculations were carried out by varying the fuel injection pressure from 0.5 to 1.5MPa. This amounts to changing the fuel mass flow rate or the overall equivalence ratio for different operating regimes. Figure 1 shows the instantaneous temperature fields in the supersonic combustor at four different conditions. The dark blue region represents the hot burned gases. At the fuel injection pressure of 0.5 MPa, the flame is stably anchored, but the flow field exhibits a high-amplitude oscillation. At the fuel injection pressure of 1.0 MPa, the Mach reflection occurs ahead of the injector. The interaction between the incoming air and the injection flow becomes much more complex, and the fuel/air mixing is strongly enhanced. The Mach reflection oscillates and results in a strong fluctuation in the combustor wall pressure. At the fuel injection pressure of 1.5MPa, the flow inside the combustor becomes nearly choked and the Mach reflection is displaced forward. The leading shock wave moves slowly toward the inlet, and eventually causes the combustor-upstart due to the thermal choking. The cavity appears to play a secondary role in driving the flow unsteadiness, in spite of its influence on the fuel/air mixing and flame evolution. Further investigation is necessary on this issue. The present study features detailed resolution of the flow and flame dynamics in the combustor, which was not typically available in most of the previous works. In particular, the oscillatory flow characteristics are captured at a scale sufficient to identify the underlying physical mechanisms. Much of the flow unsteadiness is not related to the cavity, but rather to the intrinsic unsteadiness in the flowfield, as also shown experimentally by Ben-Yakar et al. [6], The interactions between the unsteady flow and flame evolution may cause a large excursion of flow oscillation. The work appears to be the first of its kind in the numerical study of combustion oscillations in a supersonic combustor, although a similar phenomenon was previously reported experimentally. A more comprehensive discussion will be given in the final paper presented at the colloquium.

  • PDF