Structural changes in latosols of the cerrado region: II - soil compressive behavior and modeling of additional compaction

Currently in Brazil, as in other parts of the world, the concern is great with the increase of degraded agricultural soil, which is mostly related to the occurrence of soil compaction. Although soil texture is recognized as a very important component in the soil compressive behaviors, there are few studies that quantify its influence on the structural changes of Latosols in the Brazilian Cerrado region. This study aimed to evaluate structural changes and the compressive behavior of Latosols in Rio Verde, Goiás, through the modeling of additional soil compaction. The study was carried out using five Latosols with very different textures, under different soil compaction levels. Water retention and soil compression curves, and bearing capacity models were determined from undisturbed samples collected on the B horizons. Results indicated that clayey and very clayey Latosols were more susceptible to compression than medium-textured soils. Soil compression curves at density values associate with edaphic functions were used to determine the beneficial pressure (σ b) , i.e., pressure with optimal water retention, and critical pressure (σcrMAC), i.e., pressure with macroporosity below critical levels. These pressure values were higher than the preconsolidation pressure (σp), and therefore characterized as additional compaction. Based on the compressive behavior of these Latosols, it can be concluded that the combined preconsolidation pressure, beneficial pressure and critical pressure allow a better understanding of compression processes of Latosols.

Using numerical classification of profiles based on Vis-NIR spectra to distinguish soils from the Piracicaba Region, Brazil

Considering that information from soil reflectance spectra is underutilized in soil classification, this paper aimed to evaluate the relationship of soil physical, chemical properties and their spectra, to identify spectral patterns for soil classes, evaluate the use of numerical classification of profiles combined with spectral data for soil classification. We studied 20 soil profiles from the municipality of Piracicaba, State of São Paulo, Brazil, which were morphologically described and classified up to the 3rd category level of the Brazilian Soil Classification System (SiBCS). Subsequently, soil samples were collected from pedogenetic horizons and subjected to soil particle size and chemical analyses. Their Vis-NIR spectra were measured, followed by principal component analysis. Pearson's linear correlation coefficients were determined among the four principal components and the following soil properties: pH, organic matter, P, K, Ca, Mg, Al, CEC, base saturation, and Al saturation. We also carried out interpretation of the first three principal components and their relationships with soil classes defined by SiBCS. In addition, numerical classification of the profiles based on the OSACA algorithm was performed using spectral data as a basis. We determined the Normalized Mutual Information (NMI) and Uncertainty Coefficient (U). These coefficients represent the similarity between the numerical classification and the soil classes from SiBCS. Pearson's correlation coefficients were significant for the principal components when compared to sand, clay, Al content and soil color. Visual analysis of the principal component scores showed differences in the spectral behavior of the soil classes, mainly among Argissolos and the others soils. The NMI and U similarity coefficients showed values of 0.74 and 0.64, respectively, suggesting good similarity between the numerical and SiBCS classes. For example, numerical classification correctly distinguished Argissolos from Latossolos and Nitossolos. However, this mathematical technique was not able to distinguish Latossolos from Nitossolos Vermelho férricos, but the Cambissolos were well differentiated from other soil classes. The numerical technique proved to be effective and applicable to the soil classification process.

Numerical procedure for estimating temperature in solarized soils

The objective of this work was to develop a simplified numerical procedure for the estimation of accumulated monthly hours of solarized soil temperatures. The proposed model requires monthly means of daily solar radiation and maximum air temperature as input data, and a daily pattern of temperature variation assumed to be sine-shaped. The procedure was verified using observations made during the years 1992 and 1993 in Jaguariúna, SP. The proposed procedure can predict monthly temperature hours at 10 cm depth in the solarized soil, with acceptable accuracy, in the region for which it was developed.

Method of numerical correction of errors occasioned by delay of records during the monitoring of environmental variables of interest for animal production

Thermal and air conditions inside animal facilities change during the day due to the influence of the external environment. For statistical and geostatistical analyses to be representative, a large number of points spatially distributed in the facility area must be monitored. This work suggests that the time variation of environmental variables of interest for animal production, monitored within animal facility, can be modeled accurately from discrete-time records. The aim of this study was to develop a numerical method to correct the temporal variations of these environmental variables, transforming the data so that such observations are independent of the time spent during the measurement. The proposed method approached values recorded with time delays to those expected at the exact moment of interest, if the data were measured simultaneously at the moment at all points distributed spatially. The correction model for numerical environmental variables was validated for environmental air temperature parameter, and the values corrected by the method did not differ by Tukey's test at 5% significance of real values recorded by data loggers.

Numerical evaluation of the modulus of longitudinal elasticity in structural round timber elements of the Eucalyptus genus

Currently, the standards that deal with the determination of the properties of rigidity and strength for structural round timber elements do not take in consideration in their calculations and mathematical models the influence of the existing irregularities in the geometry of these elements. This study has as objective to determine the effective value of the modulus of longitudinal elasticity for structural round timber pieces of the Eucalyptus citriodora genus by a technique of optimization allied to the Inverse Analysis Method, to the Finite Element Method and the Least Square Method.

Localization of metastasis within the sentinel lymph node biopsies: a predictor of additional axillary spread of breast cancer?

PURPOSE: To explore the relationship between morphological characteristics and histologic localization of metastasis within sentinel lymph nodes (SLN) and axillary spread in women with breast cancer. METHODS: We selected 119 patients with positive SLN submitted to complete axillary lymph node dissection from July 2002 to March 2007. We retrieved the age of patients and the primary tumor size. In the primary tumor, we evaluated histologic and nuclear grade, and peritumoral vascular invasion (PVI). In SLNs we evaluated the size of metastasis, their localization in the lymph node, number of foci, number of involved lymph nodes, and extranodal extension. RESULTS: Fifty-one (42.8%) patients had confirmed additional metastasis in non-sentinel lymph nodes (NLSN). High histologic grade, PVI, intraparenchymatous metastasis, extranodal neoplastic extension and size of metastasis were associated with positive NLSN. SLN metastasis affecting the capsule were associated to low risk incidence of additional metastasis. After multivariate analysis, PVI and metastasis size in the SLN remained as the most important risk factors for additional metastasis. CONCLUSIONS:The risk of additional involvement of NSLN is higher in patients with PVI and it increases progressively according the histologic localization in the lymph node, from capsule, where the afferent lymphatic channel arrives, to the opposite side of capsule promoting the extranodal extension. Size of metastasis greater than 6.0 mm presents higher risk of additional lymph node metastasis.

Additional information about a mange outbreak by Allopsoroptoides galli (Acari: Psoroptoididae) in commercial laying hens in the state of São Paulo, Brazil

This paper reports additional information about a mange outbreak by the mite Allopsoroptoides galli in a commercial egg-laying hen facility in the state of São Paulo, Brazil. About half of the 76,000 multi-age birds of the flock were affected. Experimental infestations carried out on naive hens resulted in clinical signs similar to those diagnosed in naturally infested hens, such as generalized scaly dermatitis, presence of mucus-like material and yellowish crusts on the skin and around the calami, feather loss and strong unpleasant odor. About 30% drop of egg production was estimated. The possible source of infestation were wild birds identified on the ground and roofs of the sheds.

Dynamics of coherent vortices in mixing layers using direct numerical and large-eddy simulations

Coherent vortices in turbulent mixing layers are investigated by means of Direct Numerical Simulation (DNS) and Large-Eddy Simulation (LES). Subgrid-scale models defined in spectral and physical spaces are reviewed. The new "spectral-dynamic viscosity model", that allows to account for non-developed turbulence in the subgrid-scales, is discussed. Pseudo-spectral methods, combined with sixth-order compact finite differences schemes (when periodic boundary conditions cannot be established), are used to solve the Navier- Stokes equations. Simulations in temporal and spatial mixing layers show two types of pairing of primary Kelvin-Helmholtz (KH) vortices depending on initial conditions (or upstream conditions): quasi-2D and helical pairings. In both cases, secondary streamwise vortices are stretched in between the KH vortices at an angle of 45° with the horizontal plane. These streamwise vortices are not only identified in the early transitional stage of the mixing layer but also in self-similar turbulence conditions. The Re dependence of the "diameter" of these vortices is analyzed. Results obtained in spatial growing mixing layers show some evidences of pairing of secondary vortices; after a pairing of the primary Kelvin-Helmholtz (KH) vortices, the streamwise vortices are less numerous and their diameter has increased than before the pairing of KH vortices.

An experimental/numerical study of the turbulent boundary layer development along a surface with a sudden change in roughness

A theory for the description of turbulent boundary layer flows over surfaces with a sudden change in roughness is considered. The theory resorts to the concept of displacement in origin to specify a wall function boundary condition for a kappa-epsilon model. An approximate algebraic expression for the displacement in origin is obtained from the experimental data by using the chart method of Perry and Joubert(J.F.M., vol. 17, pp. 193-122, 1963). This expression is subsequently included in the near wall logarithmic velocity profile, which is then adopted as a boundary condition for a kappa-epsilon modelling of the external flow. The results are compared with the lower atmospheric observations made by Bradley(Q. J. Roy. Meteo. Soc., vol. 94, pp. 361-379, 1968) as well as with velocity profiles extracted from a set of wind tunnel experiments carried out by Avelino et al.( 7th ENCIT, 1998). The measurements are found to be in good agreement with the theoretical computations. The skin-friction coefficient was calculated according to the chart method of Perry and Joubert(J.F.M., vol. 17, pp. 193-122, 1963) and to a balance of the integral momentum equation. In particular, the growth of the internal boundary layer thickness obtained from the numerical simulation is compared with predictions of the experimental data calculated by two methods, the "knee" point method and the "merge" point method.

On the numerical simulation of machining processes

Numerical simulation of machining processes can be traced back to the early seventies when finite element models for continuous chip formation were proposed. The advent of fast computers and development of new techniques to model large plastic deformations have favoured machining simulation. Relevant aspects of finite element simulation of machining processes are discussed in this paper, such as solution methods, material models, thermo-mechanical coupling, friction models, chip separation and breakage strategies and meshing/re-meshing strategies.

Numerical study of Petrov-Galerkin formulations for the shallow water wave equations

The behavior of Petrov-Galerkin formulations for shallow water wave equations is evaluated numerically considering typical one-dimensional propagation problems. The formulations considered here use stabilizing operators to improve classical Galerkin approaches. Their advantages and disadvantages are pointed out according to the intrinsic time scale (free parameter) which has a particular importance in this kind of problem. The influence of the Courant number and the performance of the formulation in dealing with spurious oscillations are adressed.

Numerical study on mixed convection in a horizontal flow past a square porous cavity using UNIFAES scheme

Mixed convection on the flow past a heated length and past a porous cavity located in a horizontal wall bounding a saturated porous medium is numerically simulated. The cavity is heated from below. The steady-state regime is studied for several intensities of the buoyancy effects due to temperature variations. The influences of Péclet and Rayleigh numbers on the flow pattern and the temperature distributions are examined. Local and global Nusselt numbers are reported for the heated surface. The convective-diffusive fluxes at the volume boundaries are represented using the UNIFAES, Unified Finite Approach Exponential-type Scheme, with the Power-Law approximation to reduce the computing time. The conditions established by Rivas for the quadratic order of accuracy of the central differencing to be maintained in irregular grids are shown to be extensible to other quadratic schemes, including UNIFAES, so that accuracy estimates could be obtained.

On the numerical integration of rigid body nonlinear dynamics in presence of parameters singularities

One of the main complexities in the simulation of the nonlinear dynamics of rigid bodies consists in describing properly the finite rotations that they may undergo. It is well known that, to avoid singularities in the representation of the SO(3) rotation group, at least four parameters must be used. However, it is computationally expensive to use a four-parameters representation since, as only three of the parameters are independent, one needs to introduce constraint equations in the model, leading to differential-algebraic equations instead of ordinary differential ones. Three-parameter representations are numerically more efficient. Therefore, the objective of this paper is to evaluate numerically the influence of the parametrization and its singularities on the simulation of the dynamics of a rigid body. This is done through the analysis of a heavy top with a fixed point, using two three-parameter systems, Euler's angles and rotation vector. Theoretical results were used to guide the numerical simulation and to assure that all possible cases were analyzed. The two parametrizations were compared using several integrators. The results show that Euler's angles lead to faster integration compared to the rotation vector. An Euler's angles singular case, where representation approaches a theoretical singular point, was analyzed in detail. It is shown that on the contrary of what may be expected, 1) the numerical integration is very efficient, even more than for any other case, and 2) in spite of the uncertainty on the Euler's angles themselves, the body motion is well represented.

Numerical and experimental investigation of thermal louvers for space applications

Thermal louvers, using movable or rotating shutters over a radiating surface, have gained a wide acceptance as highly efficient devices for controlling the temperature of a spacecraft. This paper presents a detailed analysis of the performance of a rectangular thermal louver with movable blades. The radiative capacity of the louver, determined by its effective emittance, is calculated for different values of the blades opening angle. Experimental results obtained with a prototype of a spacecraft thermal louver show good agreement with the theoretical values.

Solution of aerospace problems using structured and unstructured strategies

Products developed at industries, institutes and research centers are expected to have high level of quality and performance, having a minimum waste, which require efficient and robust tools to numerically simulate stringent project conditions with great reliability. In this context, Computational Fluid Dynamics (CFD) plays an important role and the present work shows two numerical algorithms that are used in the CFD community to solve the Euler and Navier-Stokes equations applied to typical aerospace and aeronautical problems. Particularly, unstructured discretization of the spatial domain has gained special attention by the international community due to its ease in discretizing complex spatial domains. This work has the main objective of illustrating some advantages and disadvantages of numerical algorithms using structured and unstructured spatial discretization of the flow governing equations. Numerical methods include a finite volume formulation and the Euler and Navier-Stokes equations are applied to solve a transonic nozzle problem, a low supersonic airfoil problem and a hypersonic inlet problem. In a structured context, these problems are solved using MacCormack’s implicit algorithm with Steger and Warming’s flux vector splitting technique, while, in an unstructured context, Jameson and Mavriplis’ explicit algorithm is used. Convergence acceleration is obtained using a spatially variable time stepping procedure.