Substitution-Newton-Raphson method applied to the modeling of a vapour compression refrigeration system using different representations of the thermodynamic properties of R-134a

This paper gives a detailed presentation of the Substitution-Newton-Raphson method, suitable for large sparse non-linear systems. It combines the Successive Substitution method and the Newton-Raphson method in such way as to take the best advantages of both, keeping the convergence features of the Newton-Raphson with the low requirements of memory and time of the Successive Substitution schemes. The large system is solved employing few effective variables, using the greatest possible part of the model equations in substitution fashion to fix the remaining variables, but maintaining the convergence characteristics of the Newton-Raphson. The methodology is exemplified through a simple algebraic system, and applied to a simple thermodynamic, mechanical and heat transfer modeling of a single-stage vapor compression refrigeration system. Three distinct approaches for reproducing the thermodynamic properties of the refrigerant R-134a are compared: the linear interpolation from tabulated data, the use of polynomial fitted curves and the use of functions derived from the Helmholtz free energy.

Heavy rainfall equations for Santa Catarina, Brazil

Knowledge of intensity-duration-frequency (IDF) relationships of rainfall events is extremely important to determine the dimensions of surface drainage structures and soil erosion control. The purpose of this study was to obtain IDF equations of 13 rain gauge stations in the state of Santa Catarina in Brazil: Chapecó, Urussanga, Campos Novos, Florianópolis, Lages, Caçador, Itajaí, Itá, Ponte Serrada, Porto União, Videira, Laguna and São Joaquim. The daily rainfall data charts of each station were digitized and then the annual maximum rainfall series were determined for durations ranging from 5 to 1440 min. Based on these, with the Gumbel-Chow distribution, the maximum rainfall was estimated for durations ranging from 5 min to 24 h, considering return periods of 2, 5, 10, 20, 25, 50, and 100 years,. Data agreement with the Gumbel-Chow model was verified by the Kolmogorov-Smirnov test, at 5 % significance level. For each rain gauge station, two IDF equations of rainfall events were adjusted, one for durations from 5 to 120 min and the other from 120 to 1440 min. The results show a high variability in maximum intensity of rainfall events among the studied stations. Highest values of coefficients of variation in the annual maximum series of rainfall were observed for durations of over 600 min at the stations of the coastal region of Santa Catarina.

Critical-point yield model to estimate yield damage caused by Cercospora zea-maydis in corn

A model to estimate damage caused by gray leaf spot of corn (Cercospora zea-maydis) was developed from experimental field data gathered during the summer seasons of 2000/01 and during the second crop season [January-seedtime] of 2001, in the southwest of Goiás state. Three corn hybrids were grown over two seasons and on two sites, resulting in 12 experimental plots. A disease intensity gradient (lesions per leaf) was generated through application, three times over the season, of five different doses of the fungicide propiconazol. From tasseling onward, disease intensity on the ear leaf (El), and El - 1, El - 2, El + 1, and El + 2, was evaluated weekly. A manual harvest at the physiological ripening stage was followed by grain drying and cleaning. Finally, grain yield in kg.ha-1 was estimated. Regression analysis, performed between grain yield and all combinations of the number of lesions on each leaf type, generated thirty linear equations representing the damage function. To estimate losses caused by different disease intensities at different corn growth stages, these models should first be validated. Damage coefficients may be used in determining the economic damage threshold.

Critical-point model to estimate yield loss caused by Asian soybean rust

ABSTRACTA model to estimate yield loss caused by Asian soybean rust (ASR) (Phakopsora pachyrhizi) was developed by collecting data from field experiments during the growing seasons 2009/10 and 2010/11, in Passo Fundo, RS. The disease intensity gradient, evaluated in the phenological stages R5.3, R5.4 and R5.5 based on leaflet incidence (LI) and number of uredinium and lesions/cm2, was generated by applying azoxystrobin 60 g a.i/ha + cyproconazole 24 g a.i/ha + 0.5% of the adjuvant Nimbus. The first application occurred when LI = 25% and the remaining ones at 10, 15, 20 and 25-day intervals. Harvest occurred at physiological maturity and was followed by grain drying and cleaning. Regression analysis between the grain yield and the disease intensity assessment criteria generated 56 linear equations of the yield loss function. The greatest loss was observed in the earliest growth stage, and yield loss coefficients ranged from 3.41 to 9.02 kg/ha for each 1% LI for leaflet incidence, from 13.34 to 127.4 kg/ha/1 lesion/cm2 for lesion density and from 5.53 to 110.0 kg/ha/1 uredinium/cm2 for uredinium density.

Turbulent shallow-water model for orographic subgrid-scale perturbations

A parallel pseudo-spectral method for the simulation in distributed memory computers of the shallow-water equations in primitive form was developed and used on the study of turbulent shallow-waters LES models for orographic subgrid-scale perturbations. The main characteristics of the code are: momentum equations integrated in time using an accurate pseudo-spectral technique; Eulerian treatment of advective terms; and parallelization of the code based on a domain decomposition technique. The parallel pseudo-spectral code is efficient on various architectures. It gives high performance onvector computers and good speedup on distributed memory systems. The code is being used for the study of the interaction mechanisms in shallow-water ows with regular as well as random orography with a prescribed spectrum of elevations. Simulations show the evolution of small scale vortical motions from the interaction of the large scale flow and the small-scale orographic perturbations. These interactions transfer energy from the large-scale motions to the small (usually unresolved) scales. The possibility of including the parametrization of this effects in turbulent LES subgrid-stress models for the shallow-water equations is addressed.

A New Model to Determine the Dispersion of Fatigue Damage Evaluations

Reliable predictions of remaining lives of civil or mechanical structures subjected to fatigue damage are very difficult to be made. In general, fatigue damage is extremely sensitive to the random variations of material mechanical properties, environment and loading. These variations may induce large dispersions when the structural fatigue life has to be predicted. Wirsching (1970) mentions dispersions of the order of 30 to 70 % of the mean calculated life. The presented paper introduces a model to estimate the fatigue damage dispersion based on known statistical distributions of the fatigue parameters (material properties and loading). The model is developed by expanding into Taylor series the set of equations that describe fatigue damage for crack initiation.

A simple model for circadian timing by mammals

Circadian timing is structured in such a way as to receive information from the external and internal environments, and its function is the timing organization of the physiological and behavioral processes in a circadian pattern. In mammals, the circadian timing system consists of a group of structures, which includes the suprachiasmatic nucleus (SCN), the intergeniculate leaflet and the pineal gland. Neuron groups working as a biological pacemaker are found in the SCN, forming a biological master clock. We present here a simple model for the circadian timing system of mammals, which is able to reproduce two fundamental characteristics of biological rhythms: the endogenous generation of pulses and synchronization with the light-dark cycle. In this model, the biological pacemaker of the SCN was modeled as a set of 1000 homogeneously distributed coupled oscillators with long-range coupling forming a spherical lattice. The characteristics of the oscillator set were defined taking into account the Kuramoto's oscillator dynamics, but we used a new method for estimating the equilibrium order parameter. Simultaneous activities of the excitatory and inhibitory synapses on the elements of the circadian timing circuit at each instant were modeled by specific equations for synaptic events. All simulation programs were written in Fortran 77, compiled and run on PC DOS computers. Our model exhibited responses in agreement with physiological patterns. The values of output frequency of the oscillator system (maximal value of 3.9 Hz) were of the order of magnitude of the firing frequencies recorded in suprachiasmatic neurons of rodents in vivo and in vitro (from 1.8 to 5.4 Hz).