Darcy`s law and the differential equation of motion

This note addresses the relation between the differential equation of motion and Darcy`s law. It is shown that, in different flow conditions, three versions of Darcy`s law can be rigorously derived from the equation of motion.

Scaling of structures subject to impact loads when using a power law constitutive equation

It is well known that structures subjected to dynamic loads do not follow the usual similarity laws when the material is strain rate sensitive. As a consequence, it is not possible to use a scaled model to predict the prototype behaviour. In the present study, this problem is overcome by changing the impact velocity so that the model behaves exactly as the prototype. This exact solution is generated thanks to the use of an exponential constitutive law to infer the dynamic flow stress. Furthermore, it is shown that the adopted procedure does not rely on any previous knowledge of the structure response. Three analytical models are used to analyze the performance of the technique. It is shown that perfect similarity is achieved, regardless of the magnitude of the scaling factor. For the class of material used, the solution outlined has long been sought, inasmuch as it allows perfect similarity for strain rate sensitive structures subject to impact loads. (C) 2009 Elsevier Ltd. All rights reserved.

Modelling grain boundary migration during geometric dynamic recrystallization

High-angle grain boundary migration is predicted during geometric dynamic recrystallization (GDRX) by two types of mathematical models. Both models consider the driving pressure due to curvature and a sinusoidal driving pressure owing to subgrain walls connected to the grain boundary. One model is based on the finite difference solution of a kinetic equation, and the other, on a numerical technique in which the boundary is subdivided into linear segments. The models show that an initially flat boundary becomes serrated, with the peak and valley migrating into both adjacent grains, as observed during GDRX. When the sinusoidal driving pressure amplitude is smaller than 2 pi, the boundary stops migrating, reaching an equilibrium shape. Otherwise, when the amplitude is larger than 2 pi, equilibrium is never reached and the boundary migrates indefinitely, which would cause the protrusions of two serrated parallel boundaries to impinge on each other, creating smaller equiaxed grains.

AN ASSESSMENT ON THE USE OF THE DEBYE-HUCKEL EQUATION FOR THE THERMODYNAMIC MODELING OF AQUEOUS SYSTEMS CONTAINING POLYMERS AND SALTS

In this work, a study on the role of the long-range term of excess Gibbs energy models in the modeling of aqueous systems containing polymers and salts is presented. Four different approaches on how to account for the presence of polymer in the long-range term were considered, and simulations were conducted considering aqueous solutions of three different salts. The analysis of water activity curves showed that, in all cases, a liquid-phase separation may be introduced by the sole presence of the polymer in the long-range term, regardless of how it is taken into account. The results lead to the conclusion that there is no single exact solution for this problem, and that any kind of approach may introduce inconsistencies.

An extension of the Pitzer equation for the excess Gibbs energy of aqueous electrolyte systems to aqueous polyelectrolyte solutions

Pitzer`s equation for the excess Gibbs energy of aqueous solutions of low-molecular electrolytes is extended to aqueous solutions of polyelectrolytes. The model retains the original form of Pitzer`s model (combining a long-range term, based on the Debye-Huckel equation, with a short-range term similar to the virial equation where the second osmotic virial coefficient depends on the ionic strength). The extension consists of two parts: at first, it is assumed that a constant fraction of the monomer units of the polyelectrolyte is dissociated, i.e., that fraction does not depend on the concentration of the polyelectrolyte, and at second, a modified expression for the ionic strength (wherein each charged monomer group is taken into account individually) is introduced. This modification is to account for the presence of charged polyelectrolyte chains, which cannot be regarded as punctual charges. The resulting equation was used to correlate osmotic coefficient data of aqueous solutions of a single polyelectrolyte as well as of binary mixtures of a single polyelectrolyte and a salt with low-molecular weight. It was additionally applied to correlate liquid-liquid equilibrium data of some aqueous two-phase systems that might form when a polyelectrolyte and another hydrophilic but neutral polymer are simultaneously dissolved in water. A good agreement between the experimental data and the correlation result is observed for all investigated systems. (c) 2008 Elsevier B.V. All rights reserved.

Coarse-grid higher order finite-difference time-domain algorithm with low dispersion errors

Higher order (2,4) FDTD schemes used for numerical solutions of Maxwell`s equations are focused on diminishing the truncation errors caused by the Taylor series expansion of the spatial derivatives. These schemes use a larger computational stencil, which generally makes use of the two constant coefficients, C-1 and C-2, for the four-point central-difference operators. In this paper we propose a novel way to diminish these truncation errors, in order to obtain more accurate numerical solutions of Maxwell`s equations. For such purpose, we present a method to individually optimize the pair of coefficients, C-1 and C-2, based on any desired grid size resolution and size of time step. Particularly, we are interested in using coarser grid discretizations to be able to simulate electrically large domains. The results of our optimization algorithm show a significant reduction in dispersion error and numerical anisotropy for all modeled grid size resolutions. Numerical simulations of free-space propagation verifies the very promising theoretical results. The model is also shown to perform well in more complex, realistic scenarios.

Extension of a Recent Method for Obtaining Exact Solutions of the Bruce and Klute Equation

A method based on a specific power-law relationship between the hydraulic head and the Boltzmann variable was recently presented. We generalized this relationship to a range of powers and extended the solution to include the saturated zone. As a result, the new solution satisfies the Bruce and Klute equation exactly.

Time stability of soil water storage measured by neutron probe and the effects of calibration procedures in a small watershed

The knowledge of soil water storage (SWS) of soil profiles is crucial for the adoption of vegetation restoration practices. With the aim of identifying representative sites to obtain the mean SWS of a watershed, a time stability analysis of neutron probe evaluations of SWS was performed by the means of relative differences and Spearman rank correlation coefficients. At the same time, the effects of different neutron probe calibration procedures were explored on time stability analysis. mean SWS estimation. and preservation of the spatial variability of SWS. The selected watershed, with deep gullies and undulating slopes which cover an area of 20 ha, is characterized by an Ust-Sandiic Entisol and an Aeolian sandy soil. The dominant vegetation species are bunge needlegrass (Stipa bungeana Trim) and korshinsk peashrub (Carugano Korshinskii kom.). From June 11, 2007 to July 23,2008, SWS of the top1 m soil layer was evaluated for 20 dates, based on neutron probe data of 12 sampling sites. Three calibration procedures were employed: type 1, most complete, with each site having its own linear calibration equation (TrE); type II. with TrE equations extended over the whole field: and type III, with one single linear calibration curve for the whole field (UnE) and also correcting its intercept based on site specific relative difference analysis (RdE) and on linear fitting of data (RcE), both maintaining the same slope. A strong time stability of SWS estimated by TrE equations was identified. Soil particle size and soil organic matter content were recognized as the influencing factors for spatial variability of SWS. Land use influenced neither the spatial variability nor the time stability of SWS. Time stability analysis identified one site to represent the mean SWS of the whole watershed with mean absolute percentage errors of less than 10%, therefore. this site can be used as a predictor for the mean SWS of the watershed. Some equations of type II were found to be unsatisfactory to yield reliable mean SWS values or in preserving the associated soil spatial variability. Hence, it is recommended to be cautious in extending calibration equations to other sites since they might not consider the field variability. For the equations with corrected intercept (type III), which consider the spatial variability of calibration in a different way in relation to TrE, it was found that they can yield satisfactory means and standard deviation of SWS, except for the RdE equations, which largely leveled off the SWS values in the watershed. Correlation analysis showed that the neutron probe calibration was linked to soil bulk density and to organic matter content. Therefore, spatial variability of soil properties should be taken into account during the process of neutron probe calibration. This study provides useful information on the mean SWS observation with a time stable site and on distinct neutron probe calibration procedures, and it should be extended to soil water management studies with neutron probes, e.g., the process of vegetation restoration in wider area and soil types of the Loess Plateau in China. (C) 2009 Elsevier B.V. All rights reserved.

Estimating hourly net radiation for leaf wetness duration using the Penman-Monteith equation

Leaf wetness duration (LWD) is related to plant disease occurrence and is therefore a key parameter in agrometeorology. As LWD is seldom measured at standard weather stations, it must be estimated in order to ensure the effectiveness of warning systems and the scheduling of chemical disease control. Among the models used to estimate LWD, those that use physical principles of dew formation and dew and/or rain evaporation have shown good portability and sufficiently accurate results for operational use. However, the requirement of net radiation (Rn) is a disadvantage foroperational physical models, since this variable is usually not measured over crops or even at standard weather stations. With the objective of proposing a solution for this problem, this study has evaluated the ability of four models to estimate hourly Rn and their impact on LWD estimates using a Penman-Monteith approach. A field experiment was carried out in Elora, Ontario, Canada, with measurements of LWD, Rn and other meteorological variables over mowed turfgrass for a 58 day period during the growing season of 2003. Four models for estimating hourly Rn based on different combinations of incoming solar radiation (Rg), airtemperature (T), relative humidity (RH), cloud cover (CC) and cloud height (CH), were evaluated. Measured and estimated hourly Rn values were applied in a Penman-Monteith model to estimate LWD. Correlating measured and estimated Rn, we observed that all models performed well in terms of estimating hourly Rn. However, when cloud data were used the models overestimated positive Rn and underestimated negative Rn. When only Rg and T were used to estimate hourly Rn, the model underestimated positive Rn and no tendency was observed for negative Rn. The best performance was obtained with Model I, which presented, in general, the smallest mean absolute error (MAE) and the highest C-index. When measured LWD was compared to the Penman-Monteith LWD, calculated with measured and estimated Rn, few differences were observed. Both precision and accuracy were high, with the slopes of the relationships ranging from 0.96 to 1.02 and R-2 from 0.85 to 0.92, resulting in C-indices between 0.87 and 0.93. The LWD mean absolute errors associated with Rn estimates were between 1.0 and 1.5h, which is sufficient for use in plant disease management schemes.

Validation of a root water uptake model to estimate transpiration constraints

Experimental results obtained from a greenhouse trial with common bean (Phaseolus vulgaris L) plants performed to test model hypotheses regarding the onset of limiting hydraulic conditions and the shape of the transpiration reduction curve in the falling rate phase are presented. According to these hypotheses based on simulations with an upscaled single-root model, the matric flux potential at the onset of limiting hydraulic conditions is as a function of root length density and potential transpiration rate, while the relative transpiration in the falling rate phase equals the relative matric flux potential. Transpiration of bean plants in water stressed pots with four different soils was determined daily by weighing and compared to values obtained from non-stressed pots. This procedure allowed determining the onset of the falling rate phase and corresponding soil hydraulic conditions. At the onset of the falling rate phase, the value of matric flux potential M(I) showed to differ in order of magnitude from the model predicted value for three out of four soils. This difference between model and experiment can be explained by the heterogeneity of the root distribution which is not considered by the model. An empirical factor to deal with this heterogeneity should be included in the model to improve predictions. Comparing the predictions of relative transpiration in the falling rate phase using a linear shape with water content, pressure head or matric flux potential, the matric flux potential based reduction function, in agreement with the hypothesis, showed the best performance, while the pressure head based equation resulted in the highest deviations between observed and predicted values of relative transpiration rates. (C) 2010 Elsevier B.V. All rights reserved.

Modeling of transpiration reduction in van Genuchten-Mualem type soils

We derive an analytic expression for the matric flux potential (M) for van Genuchten-Mualem (VGM) type soils which can also be written in terms of a converging infinite series. Considering the first four terms of this series, the accuracy of the approximation was verified by comparing it to values of M estimated by numerical finite difference integration. Using values of the parameters for three soils from different texture classes, the proposed four-term approximation showed an almost perfect match with the numerical solution, except for effective saturations higher than 0.9. Including more terms reduced the discrepancy but also increased the complexity of the equation. The four-term equation can be used for most applications. Cases with special interest in nearly saturated soils should include more terms from the infinite series. A transpiration reduction function for use with the VGM equations is derived by combining the derived expression for M with a root water extraction model. The shape of the resulting reduction function and its dependency on the derivative of the soil hydraulic diffusivity D with respect to the soil water content theta is discussed. Positive and negative values of dD/d theta yield concave and convex or S-shaped reduction functions, respectively. On the basis of three data sets, the hydraulic properties of virtually all soils yield concave reduction curves. Such curves based solely on soil hydraulic properties do not account for the complex interactions between shoot growth, root growth, and water availability.

Comparative analysis of the gas-phase reactions of cylindrospermopsin and the difference in the alkali metal cation mobility

Cylindrospermopsin (CYN) belongs to a group of toxins produced by several strains of freshwater cyanobacteria. It is a compact zwitterionic molecule composed of a uracil section and a tricyclic guanidinium portion with a primarily hepatotoxic effect. Using low multi-stage and high-resolution mass spectrometry, the gas-phase reactions of this toxin have been investigated. Our data show that collision-induced dissociation (CID) spectra of CYN are dominated by neutral losses, and three major initial fragmentation pathways are clearly distinguishable. Interestingly, comparative analysis of protonated and cationizated molecules showed a significant difference in the balance of the SO(3) and terminal ring elimination. These data indicate that the differential ion mobility of H(+), Li(+), Na(+) and K(+) leads to different fragmentation pathways, giving rise to mass spectra with different profiles. Copyright (C) 2008 John Wiley & Sons, Ltd.

Existence results for a damped second order abstract functional differential equation with impulses

In this work we study the existence and regularity of mild solutions for a damped second order abstract functional differential equation with impulses. The results are obtained using the cosine function theory and fixed point criterions. (C) 2009 Elsevier Ltd. All rights reserved.

New Specific Equation to Estimate Resting Energy Expenditure in Severely Obese Patients

Calculating the estimated resting energy expenditure (REE) in severely obese patients is useful, but there is controversy concerning the effectiveness of available prediction equations (PE) using body weight (BW). We evaluated the efficacy of REE equations against indirect calorimetry (IC) in severely obese subjects and aimed to develop a new equation based on body composition compartments. One hundred and twenty severely obese patients had their REE measured (MREE) by IC and compared to the most commonly used PE (Harris-Benedict (HB), Ireton-Jones, Owen, and Mifflin St. Jeor). In a random sample (n = 60), a new REE equation based on fat-free mass (FFM) was developed and validated. All PE studied failed to estimate REE in severe obesity (low concordance correlation coefficient (CCC) and limits of agreement of nearly 50% of the sample +/- 10% of MREE). The HB equation using actual BW exhibited good results for all samples when compared to IC (2,117 +/- 518 kcal/day by HB vs. 2,139 +/- 423 kcal/day by MREE, P > 0.01); these results were blunted when patients were separated by gender (2,771 vs. 2,586 kcal/day, P < 0.001 in males and 1,825 vs. 1,939 kcal/day, P < 0.001 in females). A new resting energy expenditure equation prediction was developed using FFM, Horie-Waitzberg, & Gonzalez, expressed as 560.43 + (5.39 x BW) + (14.14 x FFM). The new resting energy expenditure equation prediction, which uses FFM and BW, demonstrates higher accuracy, precision, CCC, and limits of agreement than the standard PE in patients when compared to MREE (2,129 +/- 45 kcal/day vs. 2,139 +/- 423 kcal/day, respectively, P = 0.1). The new equation developed to estimate REE, which takes into account both FFM and BW, provides better results than currently available equations.

The use of sodium-chloride difference and chloride-sodium ratio as strong ion difference surrogates in the evaluation of metabolic acidosis in critically ill patients

Purpose: Inorganic apparent strong ion difference (SIDai) improves chloride-associated acidosis recognition in dysnatremic patients. We investigated whether the difference between sodium and chloride (Na+-C1-) or the ratio between chloride and sodium (Cl-/Na+) could be used as SIDai surrogates in mixed and dysnatremic patients. Patients and Methods: Two arterial blood samples were collected from 128 patients. Physicochemical analytical approach was used. Correlation, agreement, accuracy, sensitivity, and specificity were measured to examine whether Na(+)-C1(-) and CI(-)/Na(+) could be used instead of SIDai in the diagnosis of acidosis. Results: Na(+)-C1(-) and CF/Na+ were well correlated with SIDai (R = 0.987, P < 0.001 and R = 0.959, P < 0.001, respectively). Bias between Na(+)-C1(-) and SIDai was high (6.384 with a limit of agreement of 4.4638.305 mEq/L). Accuracy values for the identification of SIDai acidosis (<38.9 mEq/L) were 0.989 (95% confidence interval [CI], 0.980-0.998) for Na+-C1- and 0.974 (95% CI, 0.959-0.989) for Cr/Na+. Receiver operator characteristic curve showed that values revealing SIDai acidosis were less than 32.5 mEq/L for Nata- and more than 0.764 for C17Na+ with sensitivities of 94.0% and 92.0% and specificities of 97.0% and 90.0%, respectively. Nata- was a reliable S IDai surrogate in dysnatremic patients. Conclusions: Nata- and CI-/Na+ are good tools to disclose S IDai acidosis. In patients with dysnatremia, Nata- is an accurate tool to diagnose SIDai acidosis. (C) 2010 Elsevier Inc. All rights reserved.