1 / 66

Environmental System Research Laboratory Dept. of Environmental Engineering

KNU-ESRL. A Web Based Partial Differential Equation Solver as a General Modeling System for Environmental Problems. June 2007. Joon Hyun Kim, and Young Han Han|. Environmental System Research Laboratory Dept. of Environmental Engineering Kangwon National University, Korea. KNU-ESRL.

kruegerc
Download Presentation

Environmental System Research Laboratory Dept. of Environmental Engineering

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. KNU-ESRL A Web Based Partial Differential Equation Solver as a General Modeling System for Environmental Problems June 2007 Joon Hyun Kim, and Young Han Han| Environmental System Research Laboratory Dept. of Environmental Engineering Kangwon National University, Korea

  2. KNU-ESRL Introduction - Increase of importance for conservation and management of water environment - Development and application of computational model for environment management : Use of analytical solution and numerical analysis technique : Difficulty and troublesomeness of program development for each problem - Commercial partial differential equation solver (CFD) : PHOENICS, FLUENT, ABAQUS, ANSYS, etc. : Mechanical, aerospace, atomic, chemical, civil engineering , etc. : Unsuitableness to construct of main variables, parameters, and reaction mechanism of analyzed problem - Need the partial differential equation solver which can easily analyze the environmental problems

  3. KNU-ESRL 1. Research Objective • - Migration of fluids and pollutants • : Representation into general partial differential equation based upon • mass, momentum, and energy equation • - Development of partial differential equation analysis • program using multidimensional finite element model • : Numerical analysis modules for each differential term • : Independently integrated and object oriented system • : Construction of main variables, parameters, and reaction mechanism • of analyzed problems • Program can easily analyze the environmental problems • : basin, surface and subsurface water, sewer, and sea • - Operation on the basis of world wide web

  4. KNU-ESRL 2. Survey on Past Models ㅇTechniques of phenomena analysis

  5. KNU-ESRL ㅇGeneral numerical analysis methods

  6. KNU-ESRL ㅇCommercial CFD program - PHOENICS, ANSYS, ABAQUS, NASTRAN, FLUENT, etc. ㅇWatershed management model - SWMM, STORM, RRL, ILLUDAS, HSPF, UCURM, etc. ㅇSurface water hydraulic model - TIDE4, HEC-2, FESWMS-2DH, RMA2, DYNHYD2, etc. ㅇSurface water quality model - MFEMWASP, QUAL2E, WASP5, WQRRS, CE-QUAL-W2, RMA-4, etc.

  7. KNU-ESRL ㅇHydraulics and water quality model of subsurface - PM5 : MODFLOW, MOC 3D, MT3D : MFEMGAS ㅇUnderground gas migration model : CMGM ㅇComposite multiphase groundwater model : TRANSPORT, EXTRAN ㅇSewer management model : ASIM, PHOENICS ㅇPredictive control model of WWTP : FEMCL ㅇWaterworks management model

  8. KNU-ESRL ㅇWeb based operation - Terminal Client Service of Window2000 - Operation on internet using ActiveX Document of MS-VB ㅇIntegrated operation with GIS • - Combination of GIS and computational model • - Construction of image and contribute data for study area • (on ArcView) • Analysis of input data and modeling result of model • Web-GIS using ArcExplorer, MapCafe and MapObject • of ESRI, and MS-VB

  9. KNU-ESRL 3. Development of Comprehensive Partial Differential Equation Solver(WPDES) - Development of multidimensional finite element model : Minimization of numerical error and instability : Combination of module for each differential term of PDE : Analysis of boundary condition is easy : Easy combination with GIS using a square grid networks : Object oriented program based on independently numerical analysis logic - Integral modules to estimate the differential term(space, time) - Discretization of governing equation using integral modules - Algorithm of comprehensive PDE analysis model

  10. KNU-ESRL 3.1 Multidimensional Finite Element Modules - Construction of weighting residual eq. for the governing eq. - Discretization of main variables and parameters using basis function - Estimation of node coordinates for each elements using basis function

  11. KNU-ESRL ㅇTransformation of global coordinates into local ones

  12. KNU-ESRL ㅇDefinition of element integral matrices using FEM - Integral matrix without spatial differential term - Differential term for weighting function - Differential term for main variables - Integral matrix with 2nd order differential term for main variables

  13. KNU-ESRL ㅇAnalysis of boundary condition using FEM - 1st boundary condition - 2nd boundary condition 1-dim. : 2-dim. : 3-dim. : - 3rd boundary condition

  14. 3.2 Application of developed method to environmental problems KNU-ESRL ㅇ Rainfall runoff and pollutants migration model in basin - 2-dim. runoff by rainfall - Flow in waterway

  15. KNU-ESRL ㅇ Surface water hydraulic model - 1-dim. flow of surface water (continuity eq.) - 1-dim. flow of surface water (momentum eq.)

  16. KNU-ESRL - 2-dim. flow of shallow sea (continuity eq.) - 2-dim. flow of shallow sea (momentum eq.)

  17. KNU-ESRL - 2-dim. flow of lake and river (continuity eq.) - 2-dim. flow of lake and river (momentum eq.)

  18. KNU-ESRL - Assembling of continuity and momentum eq. into finite element matrix

  19. KNU-ESRL ㅇ Surface water quality model - Water quality model - Thermal change model of waterbody

  20. KNU-ESRL ㅇSubsurface water hydraulics and pollutants migration model - Subsurface water hydraulic model - Pollutants migration model

  21. KNU-ESRL ㅇMultidimensional gas flow model in soil - Gas flow model in soil - Composite migration model in gas phase

  22. KNU-ESRL ㅇComposite multiphase model for NAPL analysis

  23. KNU-ESRL ㅇSewer flow and pollution analysis model - Continuity equation - Momentum equation ㅇWaterworks management model

  24. KNU-ESRL 3.3 Modules of WPDES model - Input and output - Modules related to space domain : Basis and Weighting function : Estimation of nodes and coefficients for each element - Calculation of element matrices - Assembling of element matrices - Analysis of nonlinear system - Analysis of reaction and constitution equation - Analysis of boundary condition

  25. KNU-ESRL 3.4 Example of modules in WPDES model c c********************************************** subroutine femc c********************************************** c include ’mfemwasp.cmn' c c... form the r.h.s global matrix [sk] and c the l.h.s load vector f ... c if(node.eq.2) call asem1 if(node.eq.4) call asem2 if(node.eq.8) call asem3 c c... include boundary conditions for each system ... c if(indb1(nis).eq.1) call bound1 if(node.eq.2.and.indb2(nis).eq.1) call bound12 if(node.eq.4.and.indb2(nis).eq.1) call bound22 if(node.eq.8.and.indb2(nis).eq.1) call bound32 if(node.eq.2.and.indb3(nis).eq.1) call bound13 if(node.eq.4.and.indb3(nis).eq.1) call bound23 if(node.eq.8.and.indb3(nis).eq.1) call bound33 c c c... solve the system of linear equations ... c c if(time.eq.dt) then call solve(1,sk,f,np,ihalfb,nis) c endif call solve(2,sk,f,np,ihalfb,nis) c c... store the solution in vector ce(i,nis) ... c... reset the variable for time domain ... c do ip=1,np ce(ip,nis)=f(ip,nis) if(ce(ip,nis).le.0) ce(ip,nis)=-0.6*ce(ip,nis) oldce(ip,nis)=ce(ip,nis) end do c c return end ㅇ Concentration calculation module

  26. KNU-ESRL ㅇAssembling module c c********************************************** subroutine asem3 c********************************************** c include 'femwasp.cmn' include 'elemen.cmn' c eps=ep-1.0 c c... initialize global matrix [sk] and load vector f c c if(time.eq.dt) then do i=1,np f(i,nis)=0.d0 do j=1,iband sk(i,j,nis)=0.d0 end do end do c endif c c... evaluate parameters ... c call reaction call source c c... begin assembling process ... c do nel=1,ne c c... set local coordinates and point load in each element ... c do ia=1,node ijkia=ijk(ia,nel) xe(ia)=x(ijkia) ye(ia)=y(ijkia) ze(ia)=z(ijkia) rke(ia)=rk(ijkia,nis) soure(ia)=sour(ijkia,nis) end do c c... set parameters in each element ... c c... velocity ... c v1x=vx(nel) v1y=vy(nel) v1z=vz(nel) c c... dispersion coefficient tensor ... c d1xx=dxx(nel)*v1x d1xy=dxy(nel)*v1x d1xz=dxz(nel)*v1x d1yx=dyx(nel)*v1y d1yy=dyy(nel)*v1y d1yz=dyz(nel)*v1y d1zx=dzx(nel)*v1z d1zy=dzy(nel)*v1z d1zz=dzz(nel)*v1z c c... compute the [et], [ev], [ed], [etk], [ets] matrix ... c call element3(et,evx,evy,evz,edxx,edxy,edxz,edyx,edyy,edyz, & edzx,edzy,edzz,etk,ets,v1x,v1y,v1z,d1xx,d1xy,d1xz,d1yx,d1yy,d1yz, & d1zx,d1zy,d1zz,xe,ye,ze,rke,soure) c c c... assemble [et], [ev], [ed], [etk], [ets] and form [sk] ... c c if(time.eq.dt) then do i=1,node ijki=ijk(i,nel) do j=1,node ijkj=ijk(j,nel) ji=ijkj+(ihalfb+1)-ijki if(ji.le.0) goto 99 if(ji.gt.iband) goto 99 sk(ijki,ji,nis)=sk(ijki,ji,nis)+et(i,j)+ & ep*(evx(i,j)+evy(i,j)+evz(i,j)+edxx(i,j)+edxy(i,j)+edxz(i,j) & +edyx(i,j)+edyy(i,j)+edyz(i,j)+edzx(i,j)+edzy(i,j)+edzz(i,j) & -etk(i,j)) end do end do c endif c c... assemble [et],[ev],[ed],[etk],[ets] c and form r.h.s. vector f ... c do i=1,node fe(i,nis)=0.d0 ijki=ijk(i,nel) do j=1,node ijkj=ijk(j,nel) fe(i,nis)=fe(i,nis)+(et(i,j) &+eps*(evx(i,j)+evy(i,j)+evz(i,j)+edxx(i,j)+edxy(i,j)+edxz(i,j) & +edyx(i,j)+edyy(i,j)+edyz(i,j)+edzx(i,j)+edzy(i,j)+edzz(i,j) & -etk(i,j)))*oldce(ijkj,nis)+ets(i,j) end do f(ijki,nis)=f(ijki,nis)+fe(i,nis) end do c end do goto 999 c 99 write(*,*) 'computation stop due to band width' c 999 return end

  27. KNU-ESRL 4. Verification and Application of WPDES model 4.1 Verification of WPDES model - Analytical solution - Comparison of 1,2,3-dim modeling results for 1-dim. problem - Comparison with experimental data - Comparison with modeling result of another model

  28. KNU-ESRL ㅇModel verification using analytical solution - Element network used for model verification - Convective dominant transport (Velocity V=0.369m/day, Diffusion coeff. D=0.0001725m2/day) - Dispersive dominant transport (Velocity V=0.369m/day, Diffusion coeff. D=0.01725m2/day) - Middle case (Velocity V=0.369m/day, Diffusion coeff. D=0.001725m2/day)

  29. KNU-ESRL - Verification result using analytical solution

  30. KNU-ESRL ㅇComparison of 1,2,3-dim modeling - Element network used for model verification

  31. KNU-ESRL - 1,2, and 3-dim. modeling result

  32. KNU-ESRL - 1,2, and 3-dim. modeling result

  33. KNU-ESRL ㅇComparison with experimental data - Comparison result for gas migration problem

  34. KNU-ESRL - Comparison result for gas migration problem

  35. KNU-ESRL ㅇComparison of two modeling results - Model and input data used for verification

  36. KNU-ESRL - Comparison of modeling results

  37. KNU-ESRL 4.2 Combination of WPDES and GIS ㅇCombination of ArcView and Fortran compiler

  38. KNU-ESRL ㅇOperation of WPDES model on ArcView

  39. KNU-ESRL ㅇAnalysis of WPDES modeling result

  40. KNU-ESRL 4.3 Operation of Web based model ㅇScheme of web based operation

  41. KNU-ESRL ㅇRemote access using TSC of Window 2000

  42. KNU-ESRL ㅇRemote operation using TSC

  43. KNU-ESRL ㅇScheme of Web-GIS operation - Web based operation of GIS in combination with WPDES - Extension of application using MapObject - Construction of GUI using MS-VB - Web based operation using ActiveX Control - Preparation of web page using ASP(Active Server Page) - Application and control using IIS(Internet Information Service) and Active Directory of Window2000

  44. KNU-ESRL ㅇ Web GIS using ActiveX control

  45. KNU-ESRL ㅇ Web GIS using ArcView IMS and ArcIMS

  46. KNU-ESRL ㅇ Web based model operation using ASP

  47. KNU-ESRL ㅇConstruction of interface using VB

  48. KNU-ESRL 4.4 Application of WPDES model ㅇSurface water environment management - Calculation of pollutant loading for each sub-baisin (Youngwol Lake)

  49. KNU-ESRL - Construction of image data (site information)

  50. KNU-ESRL - Construction of image data (non-point source)

More Related