• 한국토양비료학회지
• /
• 제39권5호
• /
• pp.259-267
• /
• 2006

수리 전도도는 수리구배에 대한 플럭스의 비율을 나타내며, 포화된 토양에서의 물의 이동이 포화수리전도도이고 불포화된 토양에서의 이동이 불포화수리전도도이다. 일반적인 밭 상태에서의 토양수분 조건은 불포화수리전도도로 표시하는 것이 적절하나 그 상태를 표현하기가 쉽지 않다. 토양의 불포화 상태를 나타내는데 가장 많이 쓰이고 있는 VGM(van Genuchten Mualem) 모형은 토양수분 포텐셜과 수분함량의 함수로 구성된 모형이며 몇 가지 매개변수가 필요하다. VGM 모형의 매개변수를 얻기 위해 본 연구에서는 VGM 모형의 매개변수를 계산해주는 프로그램인 Rosetta를 사용하였다. Rosetta 모형은 신경그물 얼개(neural network)를 이용하여 토양의 물리적 자료들인 토성이나 모래, 미사, 점토 함량 또는 용적밀도나 33 kPa, 1500 kPa에서의 토양수분 함량 자료를 가지고 VGM의 매개변수인 Ko(effedive saturated hydraulic conductivity), ${\theta}r$(residual soil water content), ${\theta}s$(saturated soil water content), L, n, m(=1-1/n)을 예측하는 모형으로 미국 농무성(USDA-ARS)에서 개발한 프로그램이다. Rosetta를 이용하여 10kPa에서의 불포화수리전도도를 예측하였다. 또한 Gardner와 Wooding의 모형을 기반으로 하여 만들어진 장력침투계의 포화수리전도도 값을 Gardner식에 적용하여 1, 3, 5, 7 kPa에서의 불포화수리전도도 값을 17개 토양통을 대상으로 하여 구했다. 토양수분 potential이 3 kPa에서는 물의 이동이 거의 없는 토양들이 있었는데 반해 남계통을 비롯한 학곡통, 회곡통, 백산통, 상주통, 석천통, 예산통 등 7개의 토양은 3 kPa에서도 약간의 물의 이동이 있었다. 또한, 1 kPa에서 물의 이동은 삼각통에서 $40.8{\times}10^{-5}cm{\cdot}sec^{-1}$로 이동 속도가 가장 컸으며 그 뒤로 예산통, 화봉통, 학곡통, 백산통 등이 토양에서 빠른 속도로 이동하였다. 가천통이나 석천통 및 우곡통은 1 kPa에서의 이동 속도가 아주 느린 토양으로 판단되었다. PTF와 VG모형에 의해 얻어진 10 kPa에서의 수분함량 예측 값을 VGM 모형에 적용해 불포화수리전도도를 구했을 때, VG모형에 의한 예측 값은 존재하는 반면 PTF에 의한 값은 결측 값이 존재해 그 적용에 한계가 있었다. 그리고 1 kPa에서 불포화 수리전도도를 VGM 모형으로 예측한 값과 측정된 값을 Gardner 모형으로 해석한 값을 비교했을 때 자갈이 없는 토양에서는 일정한 경향(exponential 함수)이 존재한 반면, 자갈이 있는 토양에서는 경향을 발견할 수가 없었다. 이상의 결과로 불포화 수리전도도 특성평가에 대한 VGM 모형의 적용성을 살펴보았을 때는 우리나라와 같이 경사지가 많고 토심이 깊지 않으면서 자갈함량이 많은 토양에서는 한계가 있을 것으로 판단되었다.

• 한국수자원학회논문집
• /
• 제52권5호
• /
• pp.361-372
• /
• 2019

본 연구의 목적은 운동량방정식에서 이송가속도항을 제외한 지배방정식을 이용하여 정형 사각 격자 기반의 2차원 지표면 침수해석 모형을 개발하는 것이다. 공간적 이산화는 유한체적법을 이용하였으며, 시간적 이산화는 음해법을 적용하였다. 모형의 실행시간을 단축하기 위해서 CPU를 이용한 병렬계산 기법을 적용하였다. 개발된 모형의 검증을 위해서 해석해와 비교하고, 가상 도메인에서 수치실험을 통해 모형의 거동을 평가하였다. 또한 국내의 장호원 지역과 모로코의 Sebou 강 지역에 대해서 각기 다른 공간해상도로 침수해석을 수행하고, 그 결과를 CAESER-LISFLOOD (CLF) 모형을 이용한 해석 결과와 비교하였다. 모형의 검증 결과 해석해와 잘 일치된 모의 결과를 나타내었고, 가상 도메인에서의 흐름 해석도 타당한 것으로 평가되었다. 장호원 지역과 Sebou 강 지역에 대한 본 연구와 CLF 모형의 침수모의 결과는 침수심과 침수범위에서 서로 유사하게 나타났으며, 장호원 지역의 경우 홍수위험지도의 침수범위와도 유사한 값을 보였다. 본 연구와 CLF 모형의 모의결과에서 상이한 부분에 대해서는 각각의 모의결과를 비교 평가하였다. 연구결과 본 연구에서 제시된 모형은 홍수터에서의 침수 양상을 잘 모의할 수 있는 것으로 평가되었다. 그러나 본 연구에서 제시된 모형을 이용하여 침수해석을 할 경우에는 도메인 구성 방법과 지배방정식 및 해석 방법에 의한 모형의 특징과 한계점을 충분히 고려해야 할 것이다.

• 지질공학
• /
• 제33권1호
• /
• pp.105-119
• /
• 2023

기후 변화에 따라 용수 공급량은 줄어들고 있으며, 해안 및 도서지역은 지하수의 의존도가 높아 활용 가능한 용수의 양이 줄어들고 있다. 지속 가능한 해안 도서 지역의 물 공급을 위해 현재의 상태를 정확히 진단하고 효율적인 용수의 배분 및 관리가 필요하다. 해안 도서 지역의 지하수 흐름의 정밀한 분석을 위해 제주 동부의 성산유역을 대상으로 해안 담지하수 유출량을 산정하였다. 중요 매개변수인 담수층의 두께의 산정은 다심도 관측공에서 측정되는 전기전도도를 이용하여 수직적인 보간을 통해 추정하는 방법과 Ghyben-Herzberg의 이론적인 비율을 이용하여 추정하는 방법을 사용하였다. Ghyben-Herzberg의 비율을 사용한 값은 변화하는 담-염수 경계면의 정밀한 추정이 불가능하며, 전기전도도로 분석한 값은 담-염수 경계면의 현재 상태를 나타낼 수 있다. 매개변수는 가상의 격자에 분포시키고, 각각의 격자에 대한 해안 담지하수 유출 유량의 평균을 유역 대표 유량으로 결정하였다. 유역 규모의 연도별 해안 담지하수 유출량 산정과 유량 분포를 나타냈으며, 2018년에서 2020년까지의 해안 담지하수 유출량은 전기전도도로 추정한 담수층의 두께를 적용한 방법은 6.27 × 10⁶ m³/year, Ghyben-Herzberg의 비율을 적용한 방법은 10.87 × 10⁶ m³/year로 나타났다. 본 연구에서 나타낸 결과는 해안 도서지역의 정밀한 물수지 분석 등을 이용하여 지속가능한 용수의 공급을 결정하는 정책의 근거 자료로 활용할 수 있다.

• 대한토목학회논문집
• /
• 제32권4B호
• /
• pp.243-251
• /
• 2012

국립해양조사원의 해운대 부이로 부터 관측된 자료와 CCTV 영상자료를 분석하고 수치실험을 수행하여 해운대 해변에서 발생하는 이안류의 특성을 연구하였다. CCTV 영상자료를 분석한 결과로 부터 이안류가 발생할 때 규칙파에 가깝고 진행방향이 약간 다른 파랑이 중첩되어 입사되고 있는 것을 관찰할 수 있었으며, 부이 관측자료의 웨이브릿 분석으로부터 한 주기에 파랑에너지가 집중된 협대역 스펙트럼을 확인할 수 있었다. 이로부터 이미 알려진 이안류 발생 메커니즘가운데 하나인 벌집구조의 절점선 영역을 통한 이안류의 발생이 해운대 이안류의 주요한 발생 메커니즘임을 추론할 수 있었다. 이러한 이안류 발생 메커니즘은 다음과 같이 설명할 수 있다. 진행방향이 약간 다른 규칙파가 해안으로 전파될 때 서로 중첩되어 비선형 상호작용으로 파봉선이 벌집구조 형상으로 변형된다. 여기서, 파고가 영(zero)인 절점선이 해안선에서 외해방향으로 발달하게 된다. 파고가 상대적으로 작은 이 절점선 영역을 향하여 연안류의 질량수송이 집중되며 일시에 외해로 이안류를 발생시키게 된다. 벌집구조 파봉형상 입사파에 의한 이안류 발생을 재현하기 위해 완전 비선형 Boussinesq 방정식 모형을 이용해 수치모의를 수행하였다. 이 수치모의로부터 진행방향이 서로 다르므로, 그리고 일방향 파랑이 수중천퇴에 의하여 변형되므로 발생된 벌집구조 현상에 의하여 이안류가 잘 발생함을 확인하였다.

• 대한토목학회논문집
• /
• 제14권1호
• /
• pp.231-243
• /
• 1994

본 연구에서는 하천본류의 유량과 수질, 하수처리장 유입수량과 유입수질을 확률변수로 두고, 취수 문제와 수중보에 의한 재포기현상을 포함하는 추계학적 하천수질관리모형을 개발하였다. Streeter-Phelps식을 이용해서 각 구간 하천수질의 기대값과 분산을 계산하는 확률모형을 만들고, 최적화문제의 확률적 제약식은 chance constrained 방법을 이용해서 확정적 제약식으로 변환한다. 목적함수는 지역내 하수처리장의 년간처리비용으로 두었다. 건설비용함수와 유지관리비함수는 처리효율과 처리용량의 함수인 비선형의 단일식으로 유도되었다. 최적화문제는 비선형계획법으로 해를 구하였다. 본 모형을 한강하류부에 적용한 결과 서울시내 4개 하수처리장에서 2차처리를 하고, 지천유입수의 BOD부하량이 현재와 같을 때, 1996년 DO수질기준을 만족하는 신뢰도는 50% 정도였다. 그리고 탄천, 중량천, 안양천의 BOD부하량을 현재보다 50% 감소시켰을 때 1996년 DO수질기준을 만족하는 신뢰도는 60% 이상이었다. 따라서 한강하류부의 수질보전을 위해서는 유입지천의 수질개선과 2차처리 이상의 하수처리가 요구된다.

• 한국해안·해양공학회논문집
• /
• 제29권6호
• /
• pp.369-381
• /
• 2017

쇄파역에서 흔히 관측되는 왜도된 비선형 파랑은 정현파에 비해 상대적으로 큰 관성력을 지녀 투과율에 미치는 영향이 상당할 것으로 추정되어왔다. 본고에서는 이러한 추론을 확인하기 위해 왜도가 1.02에서 1.032에 분포하는 총 여섯 개의 파랑을 대상으로 한 수치모의를 수행하였다. 수치모의는 OpenFoam 기반 Tool box인 ihFoam을 사용하여 수행하였으며, 왜도된 비선형 파랑은 경사가 1:30인 천퇴를 이용하여 유도하였다. 투수층에서의 유동은 Navier Stokes Eq.에 추가 항력을 도입하여 해석하였으며, 모의결과 왜도된 비선형 파랑은 정현파에 비해 큰 관성력으로 인해 상대적으로 큰 투과율을 지니는 것으로 모의되었다. 또한 Porous층에서의 Damping으로 인한 파속변화 여부를 살펴 보기위해 Damping이 반영된 천수방정식으로부터 분산관계식을 새로이 유도하였으며, 그 결과 투과성 방파제에서의 파형과 마찰력 사이의 위상차가 파속에 영향을 미치는 것으로 보인다.

• 대한조선학회논문집
• /
• 제32권1호
• /
• pp.42-57
• /
• 1995

자유표면의 유동문제는 저항추진성능과 내항성능이 우수한 선박과 파랑중 작업성능이 우수한 해양구조물의 설계와 관련되어 조선해양공학분야에서 지속적으로 관심의 대상이 되어온 연구분야이다. 본 논문에서는 선체주위 유동을 정확하고 효율적으로 해석하기 위한 3차원 수치해법의 개발을 목적으로 하였다. 수치해법으로 경계요소법을 사용하였으며, 그린함수는 간단한 랜킨소오스를 사용하였다. 전 경계요소면은 8점 경계요소로 표시하여 기하학적 특성을 정밀하게 반영하고자 하였다. 자유표면에서 속도포텐셜의 변화를 정규화된 8점 경계요소에서 이중 2차 스플라인함수(bi-quadratic spline function)로 표시함으로써 자유표면에서의 수치감쇠 및 분산오차를 개선하였다. 한편 물체표면에서의 물리량은 8점 경계요소의 특성을 살려 이중 2차 다항식(bi-quadratic function)으로 근사하였다. 이와같이 계산영역에 따라 해의 특성에 부합하는 수치방법을 채택함으로써 수치해의 정확성과 효율성이 향상되도록 하였다. 개발한 수치해법의 효능을 검증하기 위해 계산예로서 정상유동 및 비정상유동의 경우 Neumann-Kelvin문제를 다루었다. 본 방법에 의한 몰수 타원체 및 Series 60선에 대한 조파저항 계산결과는 적은 파넬수를 사용하고도 기존의 계산치는 물론 실험치와 좋은 일치를 보였다. 변형된 Wigley선형에 대한 동유체력 계산결과도 기존의 실험치 및 계산치와 비교적 잘 일치하였다. 비정상 유동의 경우 랜킨소오스법에서 일반적으로 적용하는 상류방사조건은 무차원주파수가 1/4보다 큰 경우에만 유효하므로, 본 논문에서는 파동방정식 연산자를 이용하여 무차원주파수가 1/4보다 작은 경우에 적용할 수 있는 상류방사조건을 유도하였다. 수면하에서 전진하며 동요하는 소오스에 대하여 적용한 결과 본 논문에서 유도한 방사조건이 유효함을 입증하였다.

• Applied Biological Chemistry
• /
• 제27권2호
• /
• pp.55-63
• /
• 1984

Pilot scale 연속식(連速式) 김치순간살균장치(瞬間殺菌裝置)를 사용한 무우김치살균시(殺菌時) 살균효과(殺菌效果)를 분석(分析)하기위해 살균시간(殺菌時間)에 따른 생균수(生菌數)의 변화(變化)를 측정(測定)하였으며 simulation model을 설정(設定)하여 김치내(內) 생균수(生菌數)의 감소(減少)를 예측(豫測)하고 실제 실험치(實驗値)와 비교(比較) 분석(分析)함으로써 simulation model의 적용타당성(適用妥當性)을 검정(檢定)하였다. 1. 김치숙성균(熟成菌)의 D값은 가열온도(加熱溫度) 60, 64, 70, 75, $80^{\circ}C$에서 각각(各各) 2.21, 1.62, 0.73, 0.32, 0.21분(分)이었으며 z값은 $19^{\circ}C$였다. 2. 김치시료(試料)를 4 l/min의 유속(流速)으로 1cycle 살균처리(殺菌處理)하는데 소요(所要)되는 시간(時間)은 0.99분(分)으로 약 1분(分)이 소요(所要)되었다. 3. 총살균효과(總殺菌效果)에서 예열부(豫熱部)에서의 살균효과(殺菌效果)가 차지하는 비율(比率)은 30%, 살균(殺菌) 예냉(豫冷) 냉각부(冷却部)에서의 살균효과(殺菌效果)가 차지하는 비율(比率)은 70%이었다. 4. 살균부(殺菌部)의 온도(溫度)가 65, $70^{\circ}C$일때 관내(管內)를 통과하는 김치액(液)의 시간(時間)에 대한 직선적(直線的) 온도변화(溫度變化)와 지수적(指數的) 온도변화(溫度變化)로부터 유도(誘導)된 두개의 simulation model equation에 의한 예측치(豫測値)와 실제 실험치(實驗値)가 비교적(比較的) 잘 일치(一致)하였다.

• 한국지반공학회지:지반
• /
• 제14권3호
• /
• pp.73-94
• /
• 1998

사면보강 또는 굴착면의 안정성 확보를 위해 쏘일례일링 공법이 종종 적용되고 있다. 그러나 오목형태 또는 볼록형태 모서리부와 같은 특수한 지역에 쏘일네일링 공법이 적용되어질 경우, 편기각보강형태, 즉 skew 쏘일네일 형태로 주로 시공쥐고 있다. 하지만, 지금까지 skew 쏘일례일링 공법이 적용된 굴착 모서리부에 대한 3차원 안정해석 및 거동분석 등에 대한 실험이나 연구결과는 미흡한 실정이며, 따라서 보강재의 배치형태, 삽입각도 및 길이 등 관련 설계변수값 결정에 관하여 주로 경험에 의존하고 있는 실정이다. 따글윽 본 연구의 주된 목적은, skew 쏘일네일링 공법이 오목형태 굴착 모서리부에 적용되는 경우 이에 대한 3차원 한계평형 안정성 평가기법을 제시하는 데 있다. 3차원 예상 파괴흙쐐기의 형상은 FLAC 프로그램 모델링 및 해석을 통해 결정하였으며, 모서리부에 대한 3차원 침투수압 산정식의 제시 및 해석시 다층지반조건의 고려 등이 포함되었다. 또한 제시된 3차원 안정해석법 을 이용해, 관련 설계변수들의 모서리부 안정성에 미치는 영향 정도를 분석하였다. 아울러 기 제시된 볼록형태 굴착 모서리부의 3차원 안정해석 법을 이용해 skew 쏘일네일 보강패턴의 효율성을 분석하였으며, 또한 굴착과정을 통해 전면부 벽체변위 및 인접지반의 침하 등이 상대적으로 문 제시되는 볼록형태 굴착 모서리부에 대한 변위예측을 위해 준 3차원 유한요소 해석기법 및 중첩기법 등의 적용을 시도하였다.

• 한국항공우주학회지
• /
• 제48권5호
• /
• pp.323-334
• /
• 2020

경계층 천이는 극초음속 비행체의 열공력 성능에 직접적인 영향을 미치기 때문에, 성능 해석과 설계를 위해 천이지점의 정확한 예측 능력이 필수적으로 요구된다. 본 연구에서는 극초음속 경계층 천이 예측을 위한 γ-Re_θt 모델을 개선하고 검증하였다. 천이 시작 위치 결정에 영향을 주는 간헐도 수송방정식 생성항의 계수를 압축성 경계층 상사해를 기반으로 마하수와 벽면온도, 자유류 정체온도 등에 대한 함수 형태로 구성하고 적용하였다. 기본 모델의 Re_θc와 F_length 상관관계식에 마하수에 따른 천이 시작 운동량두께 레이놀즈수와 천이구간의 길이 변화를 반영할 수 있도록 마하수에 대한 관계식을 추가로 결정하여 적용하였다. 실용적인 사용을 고려하여 제안 모델을 상용 CFD 코드에 적용하였으며, 검증을 위해 모델을 사용하여 극초음속 평판과 원뿔 경계층 해석을 수행하였다. 실험 결과와의 비교를 통해 마하수와 단위 레이놀즈수 변화에 대한 개선된 예측성능을 확인하였다.