Modeling of seepage flow through layered solis

A two-dimensional finite difference model, which solves mixed type of Richards' equation, whose non-linearity is dealt with modified Picard's iteration and strongly implicit procedure to solve the resulting equations, is presented. Modeling of seepage flow through heterogeneous soils, which is common in the field is addressed in the present study. The present model can be applied to both unsaturated and saturated soils and can handle very dry initial condition and steep wetting fronts. The model is validated by comparing experimental results reported in the literature. Newness of this two dimensional model is its application on layered soils with transient seepage face development, which has not been reported in the literature. Application of the two dimensional model for studying unconfined drainage due to sudden drop of water table at seepage face in layered soils is demonstrated. In the present work different sizes of rectangular flow domain with different types of layering are chosen. Sensitivity of seepage height due to problem dimension of layered system is studied. The effect of aspect ratio on seepage face development in case of the flow through layered soil media is demonstrated. The model is also applied to random heterogeneous soils in which the randomness of the model parameters is generated using the turning band technique. The results are discussed in terms of phreatic surface and seepage height development and also flux across the seepage face. Such accurate modeling of seepage face development and quantification of flux moving across the seepage face becomes important while modeling transport problems in variably saturated media.

Dynamics of the Indian-Ocean oxygen minimum zones

In the Indian Ocean, mid-depth oxygen minimum zones (OMZs) occur in the Arabian Sea and the Bay of Bengal. The lower part of the Arabian-Sea OMZ (ASOMZ; below 400 m) intensifies northward across the basin; in contrast, its upper part (above 400 m) is located in the central/eastern basin, well east of the most productive regions along the western boundary. The Bay-of-Bengal OMZ (BBOMZ), although strong, is weaker than the ASOMZ. To investigate the processes that maintain the Indian-Ocean OMZs, we obtain a suite of solutions to a coupled biological/physical model. Its physical component is a variable-density, 6 1/2-layer model, in which each layer corresponds to a distinct dynamical regime or water-mass type. Its biological component has six compartments: nutrients, phytoplankton, zooplankton, two size classes of detritus, and oxygen. Because the model grid is non-eddy resolving (0.5 degrees), the biological model also includes a parameterization of enhanced mixing based on the eddy kinetic energy derived from satellite observations. To explore further the impact of local processes on OMZs, we also obtain analytic solutions to a one-dimensional, simplified version of the biological model. Our control run is able to simulate basic features of the oxygen, nutrient, and phytoplankton fields throughout the Indian Ocean. The model OMZs result from a balance, or lack thereof, between a sink of oxygen by remineralization and subsurface oxygen sources due primarily to northward spreading of oxygenated water from the Southern Hemisphere, with a contribution from Persian-Gulf water in the northern Arabian Sea. The northward intensification of the lower ASOMZ results mostly from horizontal mixing since advection is weak in its depth range. The eastward shift of the upper ASOMZ is due primarily to enhanced advection and vertical eddy mixing in the western Arabian Sea, which spread oxygenated waters both horizontally and vertically. Advection carries small detritus from the western boundary into the central/eastern Arabian Sea, where it provides an additional source of remineralization that drives the ASOMZ to suboxic levels. The model BBOMZ is weaker than the ASOMZ because the Bay lacks a remote source of detritus from the western boundary. Although detritus has a prominent annual cycle, the model OMZs do not because there is not enough time for significant remineralization to occur.

Productividad y concentración de nutrientes del Taiwán Cubano (Pennisetum purpureum X Pennisetum tiphoides), CT 115, en época lluviosa, 2010 en la Fincas Santa Rosa

Con el objetivo de evaluar la Productividad y concentración de nutrientes del pasto cubano (Pennisetum purpureum x Pennisetum typhoides) cv CT – 115, se llevó a cabo la presente investigación en la finca Santa Rosa propiedad de la Universidad Nacional Agraria (UNA) Managua, Nicaragua. Localizada geográficamente a los 12° 08´ 33"de latitud norte y, 86° 10´ 31" de longitud oeste (INETER, 2006). Para ello se utilizó un área total de 33 m2 la cual se subdividió en tres parcelas de 10 m2 cada una donde se realizaron cortes cada 15 días hasta llegar a los 75 días. Se evaluaron las variables : altura, (cm) biomasa fresca (kg ha-1), Materia Seca (%), Biomasa seca (kg ha-1), Nitrógeno (%), Fósforo (%), Potasio (%),Calcio (%), Magnesio (%), Hierro (ppm), Manganeso (ppm), Zinc (ppm). Para la determinación de la producción de biomasa verde y seca (kg ha-1), se utilizó la metodología propuesta por la Red Internacional de evaluación de Pasturas Tropicales (CIAT, 1982). Para determinar los parámetros MS (%) se utilizó la metodología del análisis de Weende o análisis proximal (AOAC. 1990). Para el análisis de correlación se utilizo el programa SAS VERSION 9 para cada uno de las variables evaluadas, para la determinación de las curvas de mejor ajuste se utilizó el programa CVXP32. Durante el ensayo no se aplico ningún tipo de tratamiento agronómico (fertilización, riego). Los resultados obtenidos fueron: Altura 220 cm Biomasa fresca 31,764 kg ha-1 Materia seca 24% Biomasa seca 7529 kg ha-1. Las mayores concentraciones de nutrientes en la planta fueron: N 3.21% P 0.38% K 4.94% Ca 0.24% Mg 0.18% Fe 105.0ppm Mn 53.33ppm Zn 50ppm. Los modelos de mejor ajuste son: Modelo Múltiple Multiplicativo, Modelo Richards, Modelo Harris, Función Rotacional, Tercer Grado Polinomial y Asociación Exponencial.

干旱区绿洲棉田土壤水盐运动数值模拟

利用预报校正法求解土壤水分运动的h型Richards方程.同时利用具有二阶高精度的差分法求解土壤盐分运移的对流扩散方程,此格式可避免利用一般差分格式时出现的数值弥散问题.最后,将上述数学模型对干旱区绿洲棉田,在蒸发条件下土壤水分运动和土壤盐分向上运动积累过程进行数值模拟,计算结果与测量值吻合较好.可为预报土壤盐碱化提供科学依据.

塔里木盆地绿洲农田土壤水盐运动动力学模式研究

采用了改进型Picard迭代差分法求解土壤水分运动的混合型Richards方程,避免了采用Picard迭代法求解h型Richards方程所出现的质量平衡误差问题.同时,采用迎风二次型插值差分法求解土壤盐分运移的对流扩散方程,避免了采用一般差分格式时所出现的数值弥散问题.最后,将上述动力学模型对塔里木盆地绿洲农田在灌既入渗条件下的土壤水盐运动过程进行了数值模拟,模拟结果与测量值吻合较好.

有植被覆盖条件下土壤水盐运动规律研究

利用改进型Picard迭代差分法求解饱和-非饱和土壤中水分运动的混合型Richards方程，避免了传统利用Picard迭代法求解h型Richards方程所出现的质量平衡误差等问题；利用特征 差分法求解盐分运动的对流扩散方程，也可避免利用一般差分格式时出现数值弥散以及数值波动问题。文中将上述数学模型用于模拟土壤在降雨、蒸发以及有植被覆盖条件下土壤水盐运动规律，对若干导致某些土壤盐碱化的原因也作出定量的解释，并给出土壤盐分的动态发展趋势的中长期预测预报。

Robust multiplexed model predictive control

This paper extends the recently developed multiplexed model predictive control (MMPC) concept to ensure satisfaction of hard constraints despite the action of persistent, unknown but bounded disturbances. MMPC uses asynchronous control moves on each input channel instead of synchronised moves on all channels. It offers reduced computation, by dividing the online optimisation into a smaller problem for each channel, and potential performance improvements, as the response to a disturbance is quicker, albeit via only one channel. Robustness to disturbances is introduced using the constraint tightening approach, tailored to suit the asynchronous updates of MMPC and the resulting time-varying optimisations. Numerical results are presented, involving a simple mechanical example and an aircraft control example, showing the potential computational and performance benefits of the new robust MMPC.