본 연구에서는 해수-담수-해안대수층의 비선형 상호작용을 직접 해석할 있는 PBM(Porous Body Model) 기반의 3차원 N-S Solver인 LES-WASS-3D ver 2.0을 적용하며, 해안대수층의 해수침투모의를 수행하였다. 이와 같은 N-S Solver를 적용한 해안대수층의 해수침투모의는 국내 최초 수행되는 것일 뿐만 아니라, 국외적으로도 찾아보기 어려운 새로운 수치해석방법이라고 할 수 있다. 먼저 적용하는 수치모델을 검증하기 위하여 해안대수층의 해수-담수 경계면에 관한 수리모형실험결과와 비교·검토 하여 수치모델의 타당성 및 유효성을 확인하였다. 그리고 해수위 및 지하수위 변화를 고려한 해안대수층 내의 해수침투모의를 수행하여 해수위-지하수위 차와 해수위의 비(△h/h)의 증가에 따른 해안대수층 내의 유동장 그리고 해수-담수 경계면 분포 특성에 관하여 논의하였다. 또한 기존의 비확산 수치모델에서 도출할 수 없었던, △h/h에 따른 해안대수층 내의 연직 염분농도로 부터 해수침투 특성을 파악하였으며, 최종적으로 지표화 할 수 있는 △h/h가 해안대수층 내의 해수침투거리에 미치는 영향에 대하여 분석하였다. 이 결과로부터 △h/h가 작을수록 해안대수층 내의 해수침투가 약해지는 메커니즘을 이해할 수 있었다.