🌱 Sistema Integrado de Dimensionamento de Irrigação Localizada / Integrated Micro-Irrigation Design System

Acesse a Versão em Português (PT) Go to English Version (EN)

🌱 Sistema Integrado de Dimensionamento de Irrigação Localizada (PT)

Desenvolvido por: MSc. Valdemiro Pitoro; Prof. Rodrigo Sánchez-Román | 2026 | Contacto: valdemiro.pitoro@unesp.br

Este App compreende uma plataforma interativa unificada voltada ao Dimensionamento Agronômico e Hidráulico de Subunidades de Irrigação Localizada (gotejamento superficial). O repositório contém o código-fonte e o embasamento teórico/aplicado no desenvolvimento.

O software foi estruturado para operar como um fluxo contínuo de engenharia: os resultados de demanda hídrica e a geometria espacial validados no dimensionamento agronômico alimentam automaticamente o Solver numérico topográfico do bloco hidráulico, baseado nos critérios rigorosos de Jack Keller e Ron D. Bliesner (1990) (“Sprinkle and Trickle Irrigation”).


🌐 Acesso ao Aplicativo

O painel de engenharia e simulações do DripPro Designer está totalmente operacional e pode ser acessado em:

👉 DripPro Designer - Versão Hospedada

🔒 Segurança & Atualizações Futuras: Planeja-se uma atualização incremental em breve para implementar controle de acesso restrito. Após a migração, o acesso ao pipeline completo de dimensionamento e relatórios técnicos exigirá autenticação por credenciais exclusivas (nome de usuário e senha) para proteção das simulações de projetos agrícolas.


🚀 Funcionalidades Principais

  • Conectividade Total (Fluxo Contínuo): Elimina a necessidade de reinserção de dados. A vazão real (\(q_a\)) e o espaçamento entre emissores (\(S_e\)) calculados na etapa agronômica são herdados instantaneamente pelo motor hidráulico.
  • Modelagem Bidimensional Dinâmica do Bulbo Úmido: Curvas interativas de Profundidade vs. Raio geradas a partir de dados reais de ensaios de campo.
  • Filtro Crítico de Engenharia Hidropedológica: Validação automatizada dos bulbos com base na profundidade radicular efetiva da cultura (\(0.9z \le P \le 1.2z\)), prevenindo o estresse hídrico ou a lixiviação profunda.
  • Determinação da Demanda Hídrica da Cultura (ETc) e Correção: Cálculo automático da ETc com base nas directrizes da FAO e correção devido a localização (kl), advenção (ka) e condições climáticas (kaa), reduzindo a subjetividade do projetista.
  • Solver Numérico de Otimização Topográfica: Algoritmo de varredura fina que localiza a posição física ideal de conexão (ponto de alimentação) em linhas laterais e de derivação sob aclive ou declive, equilibrando as forças de atrito e gravidade.
  • Rigor na Curva do Gotejador Não-Compensante: Incorporação das variáveis de fabricação (Coeficiente de Descarga \(K\) e Expoente de Fluxo \(x\)) do emissor. A variação admissível de pressão (\(\Delta P_{adm}\)) é deduzida diretamente da física do emissor (\(q = K \cdot H^x\)).
  • Interface Monotela Ultra-Compacta: Organização em abas dinâmicas (st.tabs) e otimização de elementos visuais via CSS para análise imediata sem rolagem vertical excessiva.

📊 Módulo I: Engenharia e Dimensionamento Agronômico

O bloco agronômico mapeia a relação solo-água-planta-clima para estabelecer o manejo operacional e a geometria de campo.

1. Parâmetros de Área e Cobertura Foliar

A área total disponível ou ocupada por planta (\(S_p\), em \(\text{m}^2\)) é dada pelo compasso de plantio:

\[S_p = E_{\text{plantas}} \times E_{\text{linhas}}\]

A fração de Área Sombreada (\(A\)) ou projeção da copa da cultura é calculada por:

\[A = \frac{\frac{\pi}{4} \times D^2}{S_p}\]

Onde \(D\) representa o diâmetro médio da copa da planta (m).

2. Coeficiente de Localização (\(K_l\))

Para corrigir a evapotranspiração de referência de acordo com a área molhada e sombreada, mitigando superestimativas, utiliza-se a mediana de quatro metodologias clássicas (limitada ao teto de \(1.00\)):

  • Aljibury: \(K_l = 1.34 \times A\)
  • Decroix: \(K_l = A + 0.10\)
  • Horlye: \(K_l = 0.5 \times (A + 1.0)\)
  • Keller: \(K_l = 0.85 \times A + 0.15\)

3. Lâminas de Irrigação e Demanda Hídrica

A Lâmina Líquida Diária (\(L_l\), em \(\text{mm/dia}\)) quantifica a evapotranspiração real corrigida:

\[L_l = ETo \times K_c \times K_l \times K_a \times K_{cc}\]

A Lâmina Bruta (\(L_b\), em \(\text{mm/dia}\)) incorpora a eficiência de aplicação (\(E_{\text{ap}} = 0.90\)) e a necessidade de lavagem de sais (\(LR\), Fração de Lixiviação) para controle da salinidade do solo na zona radicular:

\[L_b = \frac{L_l}{E_{\text{ap}} \times (1 - LR)}\]

Onde a Necessidade de Lavagem de Sais (\(LR\)) é obtida diretamente das condutividades elétricas da água de irrigação (\(CE_i\)) e do extrato de saturação do solo (\(CE_e\)), conforme relação clássica:

\[LR = \frac{CE_i}{2 \times CE_e}\]


🔍 Critério Seleção do Bulbo e Tratamento de Exceções

O aplicativo confronta os dados geométricos do bulbo úmido obtidos em campo com o sistema radicular efetivo (\(z\)) informado:

  1. Piso Mínimo (\(0.9 \times z\)): Impede que a profundidade máxima do bulbo (\(P\)) fique abaixo de 90% da raiz, evitando o déficit na zona de maior absorção.
  2. Teto Máximo (\(1.2 \times z\)): Limita \(P\) a 120% de \(z\) para evitar que a água percole profundamente, gerando desperdício e lixiviação de nutrientes.

Tratamento de Incompatibilidade: Caso o solo apresente restrições físicas de infiltração vertical e nenhum ponto atenda à restrição (\(0.9z \le P \le 1.2z\)), o sistema identifica que o conjunto de dados está vazio (df_filtrado.empty) e interrompe o cálculo de forma controlada. Isso previne falhas matemáticas em cascata e alerta o projetista para revisar a profundidade informada ou alterar o padrão de emissores.


⚙️ Módulo II: Modelação e Otimização Hidráulica

O bloco hidráulico recebe as saídas operacionais e espaciais do módulo agronômico para calibrar as pressões e garantir a uniformidade física de aplicação.

1. Curva do Emissor e Variação Admissível de Pressão (\(\Delta P_{adm}\))

A tolerância de perda de carga dentro da subunidade varia de acordo com o mecanismo regulador do gotejador:

  • Emissores Autocompensantes: Baseia-se estritamente nos limites físicos da membrana de silicone fornecidos pelo fabricante: \[\Delta P_{adm} = P_{max} - P_{min}\]
  • Emissores Não-Compensantes: Calculada a partir da derivada da curva característica do emissor (\(q = K \cdot H^x\)) sob uma variação de vazão admissível (\(Y = 10\%\)): \[\Delta P_{adm} = \left[1 - (1 - Y)^{\frac{1}{x}}\right] \times P_{\text{projeto}}\]

2. Equação Universal de Darcy-Weisbach e Regimes de Fluxo

A perda de carga por atrito contínuo em condutos forçados (\(h_f\)) é equacionada por:

\[h_f = f \times \frac{L}{D} \times \frac{v^2}{2g}\]

O fator de atrito (\(f\)) é reavaliado a cada trecho com base no número de Reynolds (\(Re\)): * Laminar (\(Re \le 2000\)): \(f = \frac{64}{Re}\) * Transição (\(2000 < Re < 4000\)): Interpolação linear de transição física. * Turbulento (\(Re \ge 4000\)): Resolvido pela aproximação explícita de Swamee-Jain, adotando a rugosidade absoluta do polietileno (\(\epsilon = 1.5 \times 10^{-6}\text{ m}\)): \[f = \frac{0.25}{\left[\log_{10}\left(\frac{\epsilon}{3.7D} + \frac{5.74}{Re^{0.9}}\right)\right]^2}\]

3. Fator de Acoplamento de Múltiplas Saídas (\(F\))

Para compensar a saída gradual de vazão ao longo dos gotejadores ou derivações, a perda de carga bruta do trecho cego é corrigida pelo Fator de Christiansen:

\[h_{f,\text{total}} = h_f \times F \quad \text{onde} \quad F = \frac{1}{3} + \frac{1}{2N} + \frac{1}{6N^2}\]


🔄 Dinâmica do Solver Numérico de Balanço Topográfico

Ao simular tubulações sob declividades ou aclives lineares, o aplicativo executa um algoritmo de Solver convergente. O Solver divide o comprimento total (\(L_{\text{total}}\)) em dois subtrechos opostos (\(L_1\) e \(L_2 = L_{\text{total}} - L_1\)) alimentados por uma única conexão comum, calculando os desníveis geométricos (\(\Delta Z = L \times S\), onde \(S\) é a declividade em decimal) e as perdas por atrito de cada lado.


🎓 Validação Científica e Metodologias Homologadas

O motor de cálculos deste aplicativo foi rigorosamente validado de forma integrada (blocos agronômico e hidráulico) aplicando resoluções manuais, exercícios práticos e metodologias propostas na literatura consagrada de engenharia de irrigação. As referências fundamentais de validação incluem:

  1. Jack Keller & Ron D. Bliesner (1990) - Sprinkle and Trickle Irrigation
    • Embasamento fundamental para a modulação hidráulica de linhas de gotejadores não-compensantes e compensantes, perda de carga contínua por Darcy-Weisbach e distribuição de pressões ao longo das linhas laterais e secundárias.
  2. Fernando Pizarro Cabello (1996) - Riegos Localizados de Alta Frecuencia (RLAF)
    • Formulação clássica dos balanços agronômicos e modelagem dimensional de irrigação localizada, definindo frações de área molhada recomendadas e dinâmica de bulbos úmidos.

Principais Variáveis de Dimensionamento Validadas:

O processo de validação em malha unificada cobriu com alta precisão e conformidade as seguintes variáveis de projeto: 1. Necessidade de água a aplicar (Lâmina Bruta / Volume de Rega): Correção da evapotranspiração potencial por localização multiautores, incorporando simultaneamente a eficiência de aplicação e a lixiviação de sais. 2. Fração de área molhada (\(P\) em % ou \(S_w\)): Percentual correspondente à superfície horizontal umedecida acumulada por planta em relação ao espaçamento geométrico, em total conformidade com as restrições literárias para o tipo de solo e arranjo de emissores. 3. Número de emissores por planta (\(n\)): Distribuição espacial ideal de pontos de emissão para cobrir a fração de área molhada alvo sem gerar sobreposição excessiva de bulbos ou percolação profunda. 4. Tempo de irrigação (\(T_i\)): Dimensionamento exato da janela operacional necessária para repor a necessidade de água calculada com base na vazão acumulada dos gotejadores. 5. Variação de pressão máxima admissível (\(\Delta P_{\text{adm}}\)): Determinação exata da tolerância de perda de carga do setor que garanta vazões uniformes, com variação máxima de vazão de 10% no emissor (\(\Delta q \le\) 10%). 6. Variação de pressão real (\(\Delta P_{\text{real}}\)): Solução numérica das perdas de carga contínuas (atrito nos condutos de PVC/PE) e localizadas (gotejadores inseridos em linha) calculadas ao longo de toda a subunidade de irrigação. 7. Vazão da subunidade de irrigação (\(Q_{\text{setor}}\)): Vazão acumulada total exigida para o suprimento hidráulico unificado do setor de campo, servindo de base sólida para o dimensionamento subsequente de redes principais e seleção de motobomba.


🌱 Integrated Micro-Irrigation Design System (EN)

Developed by: MSc. Valdemiro Pitoro; Prof. Rodrigo Sánchez-Román | 2026 | Contact: valdemiro.pitoro@unesp.br

This App comprises an interactive, unified platform designed for the Agronomic and Hydraulic Sizing of Micro-Irrigation Subunits (surface drip irrigation). This repository contains the source code and the theoretical/applied foundations utilized during development.

The software is structured as a continuous engineering workflow: the crop water requirements and spatial geometry validated in the agronomic design phase automatically feed the topographic numerical solver of the hydraulic block, based on the rigorous criteria of Jack Keller and Ron D. Bliesner (1990) (“Sprinkle and Trickle Irrigation”).


🌐 Web-App Access

The engineering design and simulation dashboard for DripPro Designer is fully operational and can be accessed at:

👉 DripPro Designer - Hosted Version

🔒 Security & Future Updates: An incremental update is planned in the near future to implement restricted access control. After migration, access to the complete design pipeline and technical reports will require authentication via exclusive credentials (username and password) to protect agricultural project simulations.


🚀 Key Features

  • Full Connectivity (Continuous Workflow): Eliminates the need for data re-entry. The actual emitter discharge (\(q_a\)) and emitter spacing (\(S_e\)) calculated in the agronomic phase are instantly inherited by the hydraulic engine.
  • Dynamic Two-Dimensional Wetting Pattern Modeling: Interactive Depth vs. Radius curves generated using real field trial data.
  • Critical Pedohydraulic Engineering Filter: Automated bulb geometry validation based on the effective crop root depth (\(0.9z \le P \le 1.2z\)), in order to prevent crop water stress or deep soil percolation.
  • Crop Water Requirement (ETc) Sizing & Correction: Automatic determination of ETc using FAO guidelines, adjusted for localization (\(k_l\)), advection (\(k_a\)), and microclimate conditions (\(k_{cc}\)), reducing designer subjectivity.
  • Topographic Sizing Numerical Solver: A fine-sweep algorithm that finds the optimal physical manifold inlet position (feed point) for lateral and manifold lines on upward or downward slopes, balancing friction and gravitational head.
  • Precision in Non-Pressure-Compensating Emitter Curves: Incorporates emission manufacturing variables (discharge coefficient \(K\) and flow exponent \(x\)). The allowable pressure variation (\(\Delta P_{adm}\)) is deduced directly from the emitter physics (\(q = K \cdot H^x\)).
  • Compact Single-Screen Interface: Organized into dynamic tabs (st.tabs) and styled carefully via CSS for instant analysis without excessive vertical scrolling.

📊 Module I: Agronomic Sizing and Engineering

The agronomic module maps the soil-water-plant-atmosphere continuum to establish operational management and field layout.

1. Area and Canopy Cover Parameters

The total area available or occupied per plant (\(S_p\), in \(\text{m}^2\)) is calculated from the planting spacing:

\[S_p = E_{\text{plantas}} \times E_{\text{linhas}}\]

The Shaded Area Fraction (\(A\)) (canopy projection) is given by:

\[A = \frac{\frac{\pi}{4} \times D^2}{S_p}\]

Where \(D\) represents the average plant canopy diameter (m).

2. Sizing Coefficient (\(K_l\))

To correct reference evapotranspiration based on the wet and shaded area, avoiding overestimation, the median of four classical methodologies is taken (capped at a maximum of \(1.00\)):

  • Aljibury: \(K_l = 1.34 \times A\)
  • Decroix: \(K_l = A + 0.10\)
  • Horlye: \(K_l = 0.5 \times (A + 1.0)\)
  • Keller: \(K_l = 0.85 \times A + 0.15\)

3. Irrigation Depths and Water Demand

The Daily Net Irrigation Depth (\(L_l\), in \(\text{mm/day}\)) quantifies the corrected actual evapotranspiration:

\[L_l = ETo \times K_c \times K_l \times K_a \times K_{cc}\]

The Gross Irrigation Depth (\(L_b\), in \(\text{mm/day}\)) incorporates the application efficiency (\(E_{\text{ap}} = 0.90\)) and the salt leaching requirement (\(LR\), Leaching Fraction) for soil salinity control in the active root zone:

\[L_b = \frac{L_l}{E_{\text{ap}} \times (1 - LR)}\]

Where the Leaching Requirement (\(LR\)) is derived directly from the electrical conductivities of the irrigation water (\(CE_i\)) and the soil saturation extract (\(CE_e\)), using the classical relationship:

\[LR = \frac{CE_i}{2 \times CE_e}\]


🔍 Wetting Pattern Selection is and Exception Handling

The application evaluates the field-obtained wetting pattern geometry against the specified effective root depth (\(z\)):

  1. Floor Minimum (\(0.9 \times z\)): Prevents the maximum wetting pattern depth (\(P\)) from dropping below 90% of the root zone depth, avoiding deficit in the highest-absorption zone.
  2. Ceiling Maximum (\(1.2 \times z\)): Caps \(P\) at 120% of \(z\) to prevent water from percolating deep down, conserving water and nutrients.

Incompatibility Rules: If the soil presents vertical water movement restrictions and no data point satisfies the boundary constraint (\(0.9z \le P \le 1.2z\)), the system registers that the filtered dataset is empty (df_filtrado.empty) and halts the calculation gracefully. This prevents cascading mathematical failures and warns the designer to review the input root depth or choose a different emitter series.


⚙️ Module II: Hydraulic Modeling and Optimization

The hydraulic block takes the operational and physical outputs from the agronomic module to calibrate pressures and guarantee uniform emission.

1. Emitter Performance Curve and Allowable Pressure Variation (\(\Delta P_{adm}\))

Friction head tolerance within the irrigation subunit varies according to the emitter regulation mechanism:

  • Pressure-Compensating (PC) Emitters: Based strictly on the physical limits of the silicone diaphragm provided by the manufacturer: \[\Delta P_{adm} = P_{max} - P_{min}\]
  • Non-PC Emitters: Calculated using the derivative of the characteristic emitter curve (\(q = K \cdot H^x\)) based on an allowable discharge variation (\(Y = 10\%\)): \[\Delta P_{adm} = \left[1 - (1 - Y)^{\frac{1}{x}}\right] \times P_{\text{design}}\]

2. Darcy-Weisbach Universal Equation and Flow Regimes

The head loss due to continuous friction in pressurized conduits (\(h_f\)) is modeled as:

\[h_f = f \times \frac{L}{D} \times \frac{v^2}{2g}\]

The friction factor (\(f\)) is updated at each segment according to the Reynolds number (\(Re\)): * Laminar Flow (\(Re \le 2000\)): \(f = \frac{64}{Re}\) * Transition Flow (\(2000 < Re < 4000\)): Linear physical transition interpolation. * Turbulent Flow (\(Re \ge 4000\)): Solved using the explicit Swamee-Jain approximation, adopting the absolute roughness of polyethylene (\(\epsilon = 1.5 \times 10^{-6}\text{ m}\)): \[f = \frac{0.25}{\left[\log_{10}\left(\frac{\epsilon}{3.7D} + \frac{5.74}{Re^{0.9}}\right)\right]^2}\]

3. Multiple Outlet Friction Adjustment Factor (\(F\))

To account for the progressive reduction of discharge along lateral or manifold lines, the continuous head loss is adjusted using the Christiansen Factor:

\[h_{f,\text{total}} = h_f \times F \quad \text{where} \quad F = \frac{1}{3} + \frac{1}{2N} + \frac{1}{6N^2}\]


🔄 Topographic Balance Numerical Solver Dynamics

When simulating pipes on linear slopes, the application runs a convergent Solver algorithm. The Solver divides the total pipeline length (\(L_{\text{total}}\)) into two opposing sub-segments (\(L_1\) and \(L_2 = L_{\text{total}} - L_1\)) fed by a single common connection, calculating geometric elevation differences (\(\Delta Z = L \times S\), where \(S\) is the slope in decimal format) and friction head losses on each side.


🎓 Scientific Validation and Standardized Methodologies

The mathematical engine of this application was rigorously validated in an integrated manner (agronomic and hydraulic blocks) using hand calculations, practical exercises, and methodologies proposed in classic irrigation engineering literature. The core validation references include:

  1. Jack Keller & Ron D. Bliesner (1990) - Sprinkle and Trickle Irrigation
    • Fundamental guidelines for the hydraulic modeling of pressure-compensating and non-compensating lateral lines, continuous head loss via Darcy-Weisbach, and pressure distribution across lateral and manifold lines.
  2. Fernando Pizarro Cabello (1996) - Riegos Localizados de Alta Frecuencia (RLAF)
    • Classical formulations for agronomic water balances, wetting pattern characteristics, recommended wetted area percentages, and wet bulb dynamics.

Key Design Variables Validated:

The unified validation process verified with high precision and conformity the following design parameters: 1. Required Water Application (Gross Irrigation Depth / Water Volume): Reference evapotranspiration correction by multi-author localized coefficients, incorporating application efficiency and soil salinity leaching requirements. 2. Wetted Area Fraction (\(P\) in % or \(S_w\)): Ratio of total soil surface horizontal wetted area per plant to total planting area, matching literature guidelines for specific soil textures and emitter layouts. 3. Number of Emitters per Plant (\(n\)): Optimal spatial distribution of emission points to cover the target wetted area fraction without creating excessive overlapping or deep percolation. 4. Irrigation Duration (\(T_i\)): Exact runtime window needed to replenish the crop water requirement based on cumulative emitter discharge. 5. Allowable Pressure Variation (\(\Delta P_{\text{adm}}\)): Precise allowable head loss for the subunit to guarantee high uniformity coefficients, with a maximum emitter flow variation of 10% (\(\Delta q \le\) 10%). 6. Real Pressure Variation (\(\Delta P_{\text{real}}\)): Numerical solver of continuous friction losses (PVC/PE pipes) and localized losses (inserted emitters) calculated along the entire irrigation subunit. 7. Irrigation Subunit Discharge (\(Q_{\text{sector}}\)): Cumulative system flow demand for the field sector, forming a solid basis for sub-main sizing, main lines, and pumping station selection.