Additional scaled solutions to Richards' equation for infiltration and drainage

Warrick and Hussen developed in the nineties of the last century a method to scale Richards' equation (RE) for similar soils. In this paper, new scaled solutions are added to the method of Warrick and Hussen considering a wider range of soils regardless of their dissimilarity. Gardner-Kozeny hydraulic functions are adopted instead of Brooks-Corey functions used originally by Warrick and Hussen. These functions allow to reduce the dependence of the scaled RE on the soil properties. To evaluate the proposed method (PM), the scaled RE was solved numerically using a finite difference method with a fully implicit scheme. Three cases were considered: constant-head infiltration, constant-flux infiltration, and drainage of an initially uniform wet soil. The results for five texturally different soils ranging from sand to clay (adopted from the literature) showed that the scaled solutions were invariant to a satisfactory degree. However, slight deviations were observed mainly for the sandy soil. Moreover, the scaled solutions deviated when the soil profile was initially wet in the infiltration case or when deeply wet in the drainage condition. Based on the PM, a Philip-type model was also developed to approximate RE solutions for the constant-head infiltration. The model showed a good agreement with the scaled RE for the same range of soils and conditions, however only for Gardner-Kozeny soils. Such a procedure reduces numerical calculations and provides additional opportunities for solving the highly nonlinear RE for unsaturated water flow in soils. (C) 2011 Elsevier B.V. All rights reserved.

Numerical simulation of drop impact and jet buckling problems using the eXtended Pom-Pom model

This work presents numerical simulations of two fluid flow problems involving moving free surfaces: the impacting drop and fluid jet buckling. The viscoelastic model used in these simulations is the eXtended Pom-Pom (XPP) model. To validate the code, numerical predictions of the drop impact problem for Newtonian and Oldroyd-B fluids are presented and compared with other methods. In particular, a benchmark on numerical simulations for a XPP drop impacting on a rigid plate is performed for a wide range of the relevant parameters. Finally, to provide an additional application of free surface flows of XPP fluids, the viscous jet buckling problem is simulated and discussed. (C) 2011 Elsevier B.V. All rights reserved.

Analysis Methodology for Vortex-Induced Motion (VIM) of a Monocolumn Platform Applying the Hilbert-Huang Transform Method

Vortex-induced motion (VIM) is a highly nonlinear dynamic phenomenon. Usual spectral analysis methods, using the Fourier transform, rely on the hypotheses of linear and stationary dynamics. A method to treat nonstationary signals that emerge from nonlinear systems is denoted Hilbert-Huang transform (HHT) method. The development of an analysis methodology to study the VIM of a monocolumn production, storage, and offloading system using HHT is presented. The purposes of the present methodology are to improve the statistics analysis of VIM. The results showed to be comparable to results obtained from a traditional analysis (mean of the 10% highest peaks) particularly for the motions in the transverse direction, although the difference between the results from the traditional analysis for the motions in the in-line direction showed a difference of around 25%. The results from the HHT analysis are more reliable than the traditional ones, owing to the larger number of points to calculate the statistics characteristics. These results may be used to design risers and mooring lines, as well as to obtain VIM parameters to calibrate numerical predictions. [DOI: 10.1115/1.4003493]

0.1T magnetic resonance image in the study of experimental hydrocephalus in rats. Accuracy of the method in the measurements of the ventricular size

PURPOSE: To investigate the accuracy of 1.0T Magnetic Resonance Imaging (MRI) to measure the ventricular size in experimental hydrocephalus in pup rats. METHODS: Wistar rats were subjected to hydrocephalus by intracisternal injection of 20% kaolin (n=13). Ten rats remained uninjected to be used as controls. At the endpoint of experiment animals were submitted to MRI of brain and killed. The ventricular size was assessed using three measures: ventricular ratio (VR), the cortical thickness (Cx) and the ventricles area (VA), performed on photographs of anatomical sections and MRI. RESULTS: The images obtained through MR present enough quality to show the lateral ventricular cavities but not to demonstrate the difference between the cortex and the white matter, as well as the details of the deep structures of the brain. There were no statistically differences between the measures on anatomical sections and MRI of VR and Cx (p=0.9946 and p=0.5992, respectively). There was difference between VA measured on anatomical sections and MRI (p<0.0001). CONCLUSION: The parameters obtained through 1.0T MRI were sufficient in quality to individualize the ventricular cavities and the cerebral cortex, and to calculate the ventricular ratio in hydrocephalus rats when compared to their respective anatomic slice.

Employing perichromism for probing the properties of carboxymethyl cellulose films: an expedient, accurate method for the determination of the degree of substitution of the biopolymer derivative

The properties of films of carboxymethyl cellulose, CMC, of different degree of substitution, DS, have been examined by the use of perichromic indicators (probes). The film properties that have been determined are: empirical polarity, E-T(33); "acidity", alpha; "basicity", beta; and dipolarity/polarizability, pi*. This has been achieved by employing the following perichromic probes: 4-nitroaniline, 4-nitroanisole, 4-nitro-N,N-dimethylaniline, and 2,6-dichloro-4-(2,4,6-triphenyl-pyridinium-1-yl)phenolate, WB. The correlations between both E-T(33)- or pi* and DS were found to be linear; that between beta and DS is a second order polynomial; no obvious correlation was found between alpha and DS. The polarities of CMC films are in the range of those of butyl alcohols. As models for CMC, we have employed cellulose plus CMC of high DS; oxidized cellulose with degree of oxidation = 0.5; sodium glucuronate. The former model behaved akin to CMC, but the plots of the perichromic properties versus DS showed different slopes/intercepts. FTIR data and molecular dynamics simulations on the solvation of WB have shown that this difference can be traced to more efficient hydrogen bonding between the film of the model and the probe. This affects the intra-molecular charge-transfer energy of the latter, leading to different responses to the variation of DS. Based on the excellent linear correlation between E-T(33) and DS, for CMC from different origins, we suggest that perichromism is a simple, accurate, and expedient alternative for the determination of DS of the biopolymer derivative.

A new method to analyze the subjective visual vertical in patients with bilateral vestibular dysfunction

OBJECTIVE: The aim of this study was to assess the subjective visual vertical in patients with bilateral vestibular dysfunction and to propose a new method to analyze subjective visual vertical data in these patients. METHODS: Static subjective visual vertical tests were performed in 40 subjects split into two groups. Group A consisted of 20 healthy volunteers, and Group B consisted of 20 patients with bilateral vestibular dysfunction. Each patient performed six measurements of the subjective visual vertical test, and the mean values were calculated and analyzed. RESULTS: Analyses of the numerical values of subjective visual vertical deviations (the conventional method of analysis) showed that the mean deviation was 0.326 +/- 1.13 degrees in Group A and 0.301 +/- 1.87 degrees in Group B. However, by analyzing the absolute values of the subjective visual vertical (the new method of analysis proposed), the mean deviation became 1.35 +/- 0.48 degrees in Group A and 2.152 +/- 0.93 degrees in Group B. The difference in subjective visual vertical deviations between groups was statistically significant (p < 0.05) only when the absolute values and the range of deviations were considered. CONCLUSION: An analysis of the absolute values of the subjective visual vertical more accurately reflected the visual vertical misperception in patients with bilateral vestibular dysfunction.

The Effect of Contact on the Filter Paper Method for Measuring Soil Suction

The filter paper method is one of the most commonly used and critiqued techniques for measuring soil suction. However, many aspects related to its use still require some clarification. The results of a comprehensive study on the effect of the contact between the soil grains and soil water and the filter paper are presented herein. We investigated the influence of the equilibration time, the texture of the porous material and the degree of contact, or lack thereof, between the soil grains and the filter paper using Miamian #42 and three different types of porous material. To enhance the difference between the total suction and the matrix suction, osmotic suction was induced by saturating the specimens with a sodium chloride solution.

Spectral Induced Polarization (SIP) signatures of clayey soils containing toluene

We performed laboratory experiments to investigate the sensitivity of the Spectral Induced Polarization (SIP) method to toluene contamination in clayey soils. We used mixtures of quartzitic sand and montmorillonite as soil samples, artificially contaminated with varying amounts of toluene. Care was taken to quantify the experimental uncertainty resulting from packing since such effects must be quantified if variations in SIP signatures between samples are to be reliably interpreted in terms of the effects of hydrocarbon concentration. The SIP response of all samples following addition of toluene was monitored for a period of 40 days following sample preparation. Stepwise regression was used to examine the statistical significance of correlations between (i) clay content and (ii) toluene concentration and SIP parameters. Both single-frequency real and imaginary conductivity measurements, along with the integral chargeability, normalized chargeability, DC conductivity and time constant obtained from a Debye decomposition fitting, were examined in this regression analysis. The SIP measurements show a clear time dependence following sample preparation, indicating that samples containing toluene may take significant time to reach an equilibrium electrical response. SIP measurements are significantly related to toluene content shortly after sample preparation, when the expected dependence of SIP on clay concentration is apparently suppressed. However, for the state of electrical equilibrium after 40 days (interpreted to indicate surface chemistry at equilibrium) there is no significant relation between SIP measurements and toluene content; instead SIP measurements are then significantly correlated with clay concentration. The total chargeability, normalized chargeability and relaxation time obtained from the Debye decomposition show no correlation with toluene content, indicating that this procedure, which likely integrates over multiple mechanisms, may not be suitable for understanding relationships between SIP and hydrocarbon contamination. We find only small low-frequency polarization signals observed in relation to toluene concentration (2 mrad at 0.01 Hz), which initially decreases the interfacial polarization. Unlike earlier works, our results do not support the use of the SIP method as a tool for monitoring toluene contamination in clay soils.

Calibration of a Species-Specific Spectral Vegetation Index for Leaf Area Index (LAI) Monitoring: Example with MODIS Reflectance Time-Series on Eucalyptus Plantations

The leaf area index (LAI) is a key characteristic of forest ecosystems. Estimations of LAI from satellite images generally rely on spectral vegetation indices (SVIs) or radiative transfer model (RTM) inversions. We have developed a new and precise method suitable for practical application, consisting of building a species-specific SVI that is best-suited to both sensor and vegetation characteristics. Such an SVI requires calibration on a large number of representative vegetation conditions. We developed a two-step approach: (1) estimation of LAI on a subset of satellite data through RTM inversion; and (2) the calibration of a vegetation index on these estimated LAI. We applied this methodology to Eucalyptus plantations which have highly variable LAI in time and space. Previous results showed that an RTM inversion of Moderate Resolution Imaging Spectroradiometer (MODIS) near-infrared and red reflectance allowed good retrieval performance (R-2 = 0.80, RMSE = 0.41), but was computationally difficult. Here, the RTM results were used to calibrate a dedicated vegetation index (called "EucVI") which gave similar LAI retrieval results but in a simpler way. The R-2 of the regression between measured and EucVI-simulated LAI values on a validation dataset was 0.68, and the RMSE was 0.49. The additional use of stand age and day of year in the SVI equation slightly increased the performance of the index (R-2 = 0.77 and RMSE = 0.41). This simple index opens the way to an easily applicable retrieval of Eucalyptus LAI from MODIS data, which could be used in an operational way.

Augmented Lagrangian method with nonmonotone penalty parameters for constrained optimization

At each outer iteration of standard Augmented Lagrangian methods one tries to solve a box-constrained optimization problem with some prescribed tolerance. In the continuous world, using exact arithmetic, this subproblem is always solvable. Therefore, the possibility of finishing the subproblem resolution without satisfying the theoretical stopping conditions is not contemplated in usual convergence theories. However, in practice, one might not be able to solve the subproblem up to the required precision. This may be due to different reasons. One of them is that the presence of an excessively large penalty parameter could impair the performance of the box-constraint optimization solver. In this paper a practical strategy for decreasing the penalty parameter in situations like the one mentioned above is proposed. More generally, the different decisions that may be taken when, in practice, one is not able to solve the Augmented Lagrangian subproblem will be discussed. As a result, an improved Augmented Lagrangian method is presented, which takes into account numerical difficulties in a satisfactory way, preserving suitable convergence theory. Numerical experiments are presented involving all the CUTEr collection test problems.

Evaluation of Inductively Coupled Plasma Mass Spectrometry for Determining Ca, Cu, Fe, Mg, Mn, Se and Zn in Bovine Semen Samples using a Simple Sample Dilution Method

A simple and fast method for the determination of Ca, Cu, Fe, Mg, Mn, Se and Zn in bovine semen by quadrupole inductively coupled plasma spectrometry (q-ICP-MS) is described. Prior to analysis, samples (200 mu L) were diluted 1:50 in a solution containing 0.01% v/v Triton (R) X-100 and 0.5% v/v nitric acid and directly analyzed by ICP-MS. The limits of detection of the method are 0.3, 0.03, 0.2, 0.04, 0.04, 0.03 and 0.03 mu g L-1 for Ca-44, Cu-63, Fe-57, Mg-24, Zn-64, Se-82 and Mn-55, respectively. For purposes of comparison and method validation, four ordinary bovine semen samples were directly analyzed by ICP-MS and by flame atomic absorption spectrometry (FAAS) or graphite furnace atomic absorption spectrometry (GF AAS), with no statistical difference between the techniques at the 95% level when applying the t-test. Then, the proposed method was applied in the determinations of Ca, Cu, Fe, Mg, Mn, Se and Zn in collected samples of bovine semen from different breeds, which are used in reproduction programs and artificial insemination.

The boundedness of penalty parameters in an augmented Lagrangian method with constrained subproblems

Augmented Lagrangian methods are effective tools for solving large-scale nonlinear programming problems. At each outer iteration, a minimization subproblem with simple constraints, whose objective function depends on updated Lagrange multipliers and penalty parameters, is approximately solved. When the penalty parameter becomes very large, solving the subproblem becomes difficult; therefore, the effectiveness of this approach is associated with the boundedness of the penalty parameters. In this paper, it is proved that under more natural assumptions than the ones employed until now, penalty parameters are bounded. For proving the new boundedness result, the original algorithm has been slightly modified. Numerical consequences of the modifications are discussed and computational experiments are presented.

Processing of high resolution magic angle spinning spectra of breast cancer cells by the filter diagonalization method

Proton nuclear magnetic resonance (H-1 NMR) spectroscopy for detection of biochemical changes in biological samples is a successful technique. However, the achieved NMR resolution is not sufficiently high when the analysis is performed with intact cells. To improve spectral resolution, high resolution magic angle spinning (HR-MAS) is used and the broad signals are separated by a T-2 filter based on the CPMG pulse sequence. Additionally, HR-MAS experiments with a T-2 filter are preceded by a water suppression procedure. The goal of this work is to demonstrate that the experimental procedures of water suppression and T-2 or diffusing filters are unnecessary steps when the filter diagonalization method (FDM) is used to process the time domain HR-MAS signals. Manipulation of the FDM results, represented as a tabular list of peak positions, widths, amplitudes and phases, allows the removal of water signals without the disturbing overlapping or nearby signals. Additionally, the FDM can also be used for phase correction and noise suppression, and to discriminate between sharp and broad lines. Results demonstrate the applicability of the FDM post-acquisition processing to obtain high quality HR-MAS spectra of heterogeneous biological materials.

A simple method to measure critical angles for high-sensitivity differential refractometry

A total internal reflection-based differencial refractometer, capable of measuring the real and imaginary parts of the complex refractive index in real time, is presented. The device takes advantage of the phase difference acquired by s- and p-polarized light to generate an easily detectable minimum at the reflected profile. The method allows to sensitively measuring transparent and turbid liquid samples. (C)2012 Optical Society of America

Assessment of Pharmaceutical Equivalence: Difference Test or Equivalence Test?

Pharmaceutical equivalence is an important step towards the confirmation of similarity and Interchangeability among pharmaceutical products, particularly regarding those that win not be tested for bioequivalence. The aim of this paper is to compare traditional difference testing to two one-side equivalence tests in the assessment of pharmaceutical equivalence, by means of equivalence studies between similar, generic and reference products of acyclovir cream, atropine sulfate injection, meropenem for injection, and metronidazole injection. All tests were performed in accordance with the Brazilian Pharmacopeia or the United States Pharmacopeia. All four possible combinations of results arise in these comparisons of difference testing and equivalence testing. Most of the former did not show significant difference, whereas the latter presented similarity. We concluded that equivalence testing is more appropriate than difference testing, what can make it a useful tool to assess pharmaceutical equivalence in products that will not be tested for bioequivalence.