Docstoc

Chen-Pey

Document Sample
Chen-Pey Powered By Docstoc
					  Numerical Computation of Crops Growing
       affected by Nitrogen in Land
Chen-Pey, Huang Chia-Chen, Kuo          Whey-Fone, Tsai
  National Center for High-performance Computing, Taiwan
                  Hwai-Ping, Cheng
          National Tsing Hua University, Taiwan
                     Ming-Hsu, Li
  Office of the National Science & Technology Program for
                   Hazards Mitigation, Taiwan
                    Ming-Tee, Hung
      Council of Agriculture, Executive Yuan, Taiwan
                 Index

•   Introduction
•   Methods description
    – Model and web-based post-processor
•   Case study
•   Results
•   Concluding remarks
•   Suggestions
              Introduction

• Background:
  – The heavy metals and high-
    nitrogen farming wastewater
    affect the growth of the crops in
    Taiwan
• Goals:
  – Pollution prevention
  – Maintain and Promote efficiently
    the soil fertility
       Methods Description
• Numerical experiment and
  analysis
  – 2DSOIL
  – A Modular Simulator of Soil and
    Root Processes
• Across-platform visualization
  – VWIS
    (Visual Water information system )
  – A Web-based post-processor
    (by JDK1.3)
                     2DSOIL
Five Classes of Modules
1) Management processes
2) Transport processes
3) Biotransformation
   processes
4) Interphase exchange
   processes
5) Root growth and water
   uptake processes
             2DSOIL
• Governing flow equation:
  Richard’s eqn.

    h      h  
    k    k   1   S  0
t x  x  z   z  
    :water content
   h :pressure head
   k :soil hydraulic conductivity
   S:extraction term (removal of water
             from the soil)
                    2DSOIL Calculation Flow Chart

  Initialize
variable default                SetSurface
    setting                  weather condition in
                           Soil-atmosphere boundary                                Solutemover

Get_Grid_and_
 Boundary                                                     watermover           Heatmover
Grid and boundary                SetTDB
                       time-dependent boundary condition
                                                                                   Gasmover

Synchronizer                                                                                        Output
Time Generation

                              ShootImitator                       RootWaterUptake
                            provide shoot variable            root growth and root water uptake
                                                                                                  ScreenOutput



                                                                    provide
               Mngm                     Soilnitrogen
                                                                   solute(N)      SoluteUptake
                                     nitrogen information in soil
      solute(N) input in specified time
                                                                  concentration
    Root_Water_Uptake Module
• Relationship between (Ec) and
  (φL) for each time step
  Ec :the potential transpiration rate
  φL :the leaf water potential
        MacroChem Module
• Chemical reactions in the
  porous media
  – Dissolution-precipitation of
    gypsum and calcite
  – Ca-Mg-Na cation exchange
  – Formation of ion pairs in solution
  – Dissolution of gaseous CO2
               SoilNitrogen Module
• Transport and reaction of the nitrogen in crop systems
                   lt s
                   an
                    r
                   Pto
                     o



                                    Qi tN
                                    1m
                                    5 oo 3
                                    - bn
                                     i
                                     m i
                                       z N
                                        a O
                                       l i -
                     ia g
                     nt n
                     e r
                      r
                      l
                     moi e
                       n
                                             R
                                             OA
                                              N
                                              GC
                                               I
        o
        i
        l
        S               --
                        Q
                        nr
                         t
                         ii        1m
                                  Qi t N
                                    m
                                    i
                                   4 oo 4
                                   - bn i
                                        z N
                                       l i -
                                         a H
            7 i
            -t
            d
            Qe
             nr-
               i
       t-
       m
       ao           O
                    N
                    3
                        6
                               H
                               N 4
                                             RT
                                             E R
                                             FZ
                                              I
                                              IE
                                              L
                                   Qlt
                                     1n i
                                     3 eo
                                      - r
                                      m n
                                       i az
                                          ia
       p
       h ft
       sr
        ec
        e ii
          aon            ft
                         ii
                         co
                          an                 PO
                                              O
                                              L




                               -i t
                               h i
                               Qa
                               2
                                mi
                                 c
                                 f n
                                uo
                                                U
                                                M
                                                HS
                                                 U
              I R
              T O
               T L
               E
              LPO                               O
                                                O
                                                PL
                              Java Web-based
                 Post-processor: VWIS
• Numerical data representation based on
  WAN (Wide Area Networking)
• Apply on soil-crop simulators combined
  with information technologies
                                                                            FORTRAN            FORTRAN
                                                                             connector          Program
                                    HTTP Interface




                                                                           OBSERVATION       OBSERVATION
                   VWIS                              Information Channel
                                                                             connector



                                                                           Other System


                    Client                                                   Server
       1. Information Channel provides the interface for VWIS connect to other systems, like FORTRAN, and OBSERVATION.
         Java Web-based
      Post-processor: VWIS
• Address crop productivity and
  environmental issues in Taiwan
               Case Study
• Use 2DSOIL model
  –   Root water uptake
  –   Solute transport
  –   Soil nitrogen transformation
  –   Heat transport
  –   2-dimensional vertical simulation
  –   Crop growing during fertilization
             Case Study
• Simulation time :
   – From 130th to 270th in one year
• 2-Step fertilizing:
   – Planting and 30 days after
     planting
• Nitrogen adding rate :
   – 150kg/ha for the first
     fertilization
   – 100 kg/ha for the second
     fertilization
            Case Study
• Initial amount of nitrate-N:
  – 169 kg N/ha
• The soil organic matter
  content :
  – 1.6-kg/kg soil
• Carbon content of the organic
  matter :
  – 60 ﹪
• C/N ratio :
  – 10:1
Simulation grid
Weather data
              o t r cn t
                      e
             Misue o t n                       y r uc o d civ ,
                                                             t
                                              Hda li c n ut i y

0.6                                                                0
                                                                  10
.5
05
                                                                  10
0.5
05
.4                                                                1
0.4                                                               0.1
.3
05
                                                                  .0
                                                                  01
0.3
.2
05                                                                . 1
                                                                   0
                                                                  00
0.2
                                                                   . 0
                                                                    0
                                                                  0 01
.1
05
0.1                                                               E0
                                                                  1- 5

       20
      - 00             10
                      - 50     - 00
                                10            50
                                              -0           0
                               r s r edc
                                e
                              P s ueh a , m

                        ht o y r
                          a
                       Te f rlae 1         f r yr
                                          Ko lae 1

                        ht o y r
                          a
                       Te f rlae 2         f r yr
                                          Ko lae 2



                   Soil property information
Velocity field for 170th day         Velocity field for 210th day




                   Velocity field for 250th day
      200.00                                                                           200.00

      180.00                                                                           180.00

      160.00                                                                           160.00

      140.00                                                                           140.00

      120.00                                                                           120.00

      100.00                                                                           100.00
                                                    200.00
       80.00
                                                                                        80.00
                                                    180.00
       60.00
                                                                                        60.00
                                                    160.00
       40.00
                                                                                        40.00
                                                    140.00
       20.00
                                                                                        20.00
                                                    120.00
        0.00
           0.00 20.00                                                                    0.00
                                                    100.00                                  0.00 20.00

      temperature contour when t = 170 day                                             temperature contour when t = 210 day
                                                     80.00


Temperature distribution                             60.00
                                                                                Temperature distribution
                                                     40.00

for 170th day                                        20.00
                                                                                for 210th day
                                                       0.00
                                                          0.00 20.00

                                             temperature contour w       e
                                                                  hen tim = 250 days



                                      Temperature distribution
                                      for 250th day
                                                                                               200.00
       200.00


                                                                                               180.00
       180.00


                                                                                               160.00
       160.00


                                                                                               140.00
       140.00


                                                                                               120.00
       120.00



       100.00                                                                                  100.00



        80.00                                                                                   80.00
                                                   200.00

        60.00                                                                                   60.00
                                                   180.00

        40.00                                                                                   40.00
                                                   160.00

        20.00                                                                                   20.00
                                                   140.00

         0.00                                                                                    0.00
            0.00      20.00                                                                         0.00      20.00
                                                   120.00

       t = 170 day, contour of solute N
                                                                                               t = 210 day, contour of solute N
                                                   100.00



Concentration distribution                          80.00
                                                                                     Concentration distribution
of nitrogen for 170th day                           60.00



                                                    40.00
                                                                                     of nitrogen for 210th day
                                                    20.00



                                                      0.00
                                                         0.00      20.00



                                                  t = 250 day, contour of solute N



                                          Concentration distribution
                                          of nitrogen for 250th day
                                                   plant height

                        80.00




                        60.00
     plan height (cm)




                        40.00
         t




                        20.00




                        0.00
                            120.00   160.00        200.00         240.00   280.00
                                         time ( from 130 day to 270 day)


Plant height during growing ( day 130 to day 270)
Network computing result using VWIS
       Concluding Remarks

•   2DSOIL is introduced and
    is imported into Taiwan by
    setting at NCHC
•   Java is applied to develop
    the network computational
    environment named VWIS
               Suggestions
•   Support the development of
    the module
    –   solve the nitrogen and the
        heavy metal pollution problem
        in Taiwan
•   To support the management
    system, try to analyse
    –   the indexes of irrigation water
        quality
    –   the quantity of fertilizing in
        Taiwan

				
DOCUMENT INFO
Shared By:
Categories:
Tags:
Stats:
views:14
posted:8/19/2012
language:
pages:25