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.

Thermodynamic evaluation and modeling of proton and water exchange associated with benzamidine and berenil binding to ß-trypsin

Serine-proteases are involved in vital processes in virtually all species. They are important targets for researchers studying the relationships between protein structure and activity, for the rational design of new pharmaceuticals. Trypsin was used as a model to assess a possible differential contribution of hydration water to the binding of two synthetic inhibitors. Thermodynamic parameters for the association of bovine ß-trypsin (homogeneous material, observed 23,294.4 ± 0.2 Da, theoretical 23,292.5 Da) with the inhibitors benzamidine and berenil at pH 8.0, 25ºC and with 25 mM CaCl2, were determined using isothermal titration calorimetry and the osmotic stress method. The association constant for berenil was about 12 times higher compared to the one for benzamidine (binding constants are K = 596,599 ± 25,057 and 49,513 ± 2,732 M-1, respectively; the number of binding sites is the same for both ligands, N = 0.99 ± 0.05). Apparently the driving force responsible for this large difference of affinity is not due to hydrophobic interactions because the variation in heat capacity (DCp), a characteristic signature of these interactions, was similar in both systems tested (-464.7 ± 23.9 and -477.1 ± 86.8 J K-1 mol-1 for berenil and benzamidine, respectively). The results also indicated that the enzyme has a net gain of about 21 water molecules regardless of the inhibitor tested. It was shown that the difference in affinity could be due to a larger number of interactions between berenil and the enzyme based on computational modeling. The data support the view that pharmaceuticals derived from benzamidine that enable hydrogen bond formation outside the catalytic binding pocket of ß-trypsin may result in more effective inhibitors.

Thermodynamic properties of water sorption of jackfruit (Artocarpus heterophyllus Lam.) as a function of moisture content

The Jackfruit tree is one of the most significant trees in tropical home gardens and perhaps the most widespread and useful tree in the important genus Artocarpus. The fruit is susceptible to mechanical and biological damage in the mature state, and some people find the aroma of the fruit objectionable, particularly in confined spaces. The dehydration process could be an alternative for the exploitation of this product, and the relationship between moisture content and water activity provides useful information for its processing and storage. The aim of this study was to determine the thermodynamic properties of the water sorption of jackfruit (Artocarpus heterophyllus Lam.) as a function of moisture content. Desorption isotherms of the different parts of the jackfruit (pulp, peduncle, mesocarp, peel, and seed) were determined at four different temperatures (313.15, 323.15, 333.15, and 343.15 K) in a water activity range of 0.02-0.753 using the static gravimetric method. Theoretical and empirical models were used to model the desorption isotherms. An analytical solution of the Clausius-Clapeyron equation was proposed to calculate the isosteric heat of sorption, the differential entropy, and Gibbs' free energy using the Guggenhein-Anderson-de Boer and Oswin models considering the effect of temperature on the hygroscopic equilibrium.

Kinetic, thermodynamic properties, and optimization of barley hydration

The hydration kinetics of five barley cultivars was studied at six different temperatures ranging from 10 to 35 ºC for 32 hours applying the Peleg model. Response Surface was used to describe dynamic of the process and identify the hydration time for each cultivar. The activation energy (Ea), enthalpy (ΔH*), entropy (ΔS*), and Gibbs free energy (ΔG*) were estimated from the adjusted parameters and Arrhenius equation. Temperature had significant effect on the hydration of the five cultivars. At low temperatures, the stabilization time for hydration was faster. Peleg constants K1 and K2 decreased with increasing temperature. The cultivar BRS BRAU showed the lowest value of initial absorption rate (R0 = 0.149 kg.h-1) at 10 ºC, while the cultivar BRS BOREMA had the highest value of R0 (0.367 kg.h-1 at 35 ºC). The equilibrium moisture content (Me) increased with increasing temperature. The cultivars BRS CAUE and BRS BRAU showed the lowest values of Ea, ΔH*, ΔS* showed negative values, and ΔG* increased with increasing temperature, confirming the effect of temperature on hydration.

Study of thermodynamic water properties and moisture sorption hysteresis of mango skin

The equilibrium moisture content for adsorption and desorption isotherms of mango skin was determined using the static gravimetric method at temperatures of 20, 26, 33, 38 and 44 oC in the 0.056 to 0.873 water activity range. Both sorption curves show a decrease in equilibrium moisture content as the temperature increasing. The hysteresis effect was observed at constant water activity. The Guggenheim, Anderson, and de Boer (GAB) model presented the best fitting accuracy among a group of models and was used to determine the thermodynamic properties of water sorption. Integral enthalpy and integral entropy areas showed inverted values for the adsorption and desorption isotherms over the wide range of water activity studied. These values confirm, in energetic terms, the difference between adsorption and desorption isotherms observed in the hysteresis phenomenon. Finally, the Gibbs free energy revealed that the sorption process was spontaneous for both sorption isotherms.

Phonon spectra and thermodynamic properties of rare gas solids based on empirical and semi-empirical (ab initio) two-body potentials : a comparative study /

We study the phonon dispersion, cohesive and thermal properties of raxe gas solids Ne, Ar, Kr, and Xe, using a variety of potentials obtained from different approaches; such as, fitting to crystal properties, purely ab initio calculations for molecules and dimers or ab initio calculations for solid crystalline phase, a combination of ab initio calculations and fitting to either gas phase data or sohd state properties. We explore whether potentials derived with a certain approaxih have any obvious benefit over the others in reproducing the solid state properties. In particular, we study phonon dispersion, isothermal ajid adiabatic bulk moduli, thermal expansion, and elastic (shear) constants as a function of temperatiue. Anharmonic effects on thermal expansion, specific heat, and bulk moduli have been studied using A^ perturbation theory in the high temperature limit using the neaxest-neighbor central force (nncf) model as developed by Shukla and MacDonald [4]. In our study, we find that potentials based on fitting to the crystal properties have some advantage, particularly for Kr and Xe, in terms of reproducing the thermodynamic properties over an extended range of temperatiures, but agreement with the phonon frequencies with the measured values is not guaranteed. For the lighter element Ne, the LJ potential which is based on fitting to the gas phase data produces best results for the thermodynamic properties; however, the Eggenberger potential for Ne, where the potential is based on combining ab initio quantum chemical calculations and molecular dynamics simulations, produces results that have better agreement with the measured dispersion, and elastic (shear) values. For At, the Morse-type potential, which is based on M0ller-Plesset perturbation theory to fourth order (MP4) ab initio calculations, yields the best results for the thermodynamic properties, elastic (shear) constants, and the phonon dispersion curves.

Thermodynamic properties and Debye-Waller factor of fee materials

We have calculated the equation of state and the various thermodynamic properties of monatomic fcc crystals by minimizing the Helmholtz free energy derived in the high temperature limit for the quasiharmonic theory, QH, and the lowest-order (cubic and quartic), 'A2, anharmonic terms of the perturbation theory, PT. The total energy in each case is obtained by adding the static energy. The calculation of the thermal properties was carried out for a nearest-neighbour central-force model of the fcc lattice by means of the appropriate thermodynamic relations. We have calculated the lattice constant, the thermal expansion, the coefficient of volume expansion, the specific heat at constant volume and at constant pressure, the isothermal and adiabatic bulk moduli, and the Griineisen parameter, for the rare-gas solids Kr and Xe, and gold. Morse potential and modified Morse potential were each used to represent the atomic interaction for the three fcc materials. For most of the calculated thermodynamic properties from the QH theory, the results for Kr and Xe with the modified Morse potential show an improvement over the results for the Morse potential when compared with the experimental data. However, the results of the 'A 2 equation of state with the modified Morse potential are in good agreement with experiment only in the case of the specific heat at constant volume and at constant pressure. For Au we have calculated the lattice contribution from the QH and 'A 2 PT and the electronic contribution to the thermal properties. The electronic contribution was taken into account by using the free electron model. The results of the thermodynamic properties calculated with the modified Morse potential were similar to those obtained with the Morse potential. U sing the minimized equation of state we also calculated the Mossbauer recoilless fraction for Kr and Xe and the Debye-Waller factor (DWF) for Pb, AI, eu, Ag, and Au. The Mossbauer recoilless fraction was obtained for the above two potentials and Lennard-Jones potential. The L-J potential gives the best agreement with experiment for Kr. No experimental data exists for Xe. At low temperature the calculated DWF results for Pb, AI, and eu show a good agreement with experimental values, but at high temperature the experimental DWF results increase very rapidly. For Ag the computed values were below the expected results at all temperatures. The DWF results of the modified Morse potential for Pb, AI, eu and Ag were slightly better than those of the Morse potential. In the case of Au the calculated values were in poor agreement with experimental results. We have calculated the quasiharmonic phonon dispersion curves for Kr, Xe, eu, Ag, and Au. The calculated and experimental results of the frequencies agree quite well for all the materials except for Au where the longitudinal modes show serious discrepancies with the experimental results. In addition, the two lowest-order anharmonic contributions to the phonon frequency were derived using the Green's function method. The A 2 phonon dispersion curves have been calculated only for eu, and the results were similar to those of the QH dispersion curves. Finally, an expression for the Griineisen parameter "( has been derived from the anharmonic frequencies, and calculated for these materials. The "( results are comparable with those obtained from the thermodynamic definition.

Phonon spectra and temperature variation of thermodynamic properties of fcc metals via Finnis-Sinclair type many body potentials: Sutton-Chen and improved Sutton-Chen models

Volume(density)-independent pair-potentials cannot describe metallic cohesion adequately as the presence of the free electron gas renders the total energy strongly dependent on the electron density. The embedded atom method (EAM) addresses this issue by replacing part of the total energy with an explicitly density-dependent term called the embedding function. Finnis and Sinclair proposed a model where the embedding function is taken to be proportional to the square root of the electron density. Models of this type are known as Finnis-Sinclair many body potentials. In this work we study a particular parametrization of the Finnis-Sinclair type potential, called the "Sutton-Chen" model, and a later version, called the "Quantum Sutton-Chen" model, to study the phonon spectra and the temperature variation thermodynamic properties of fcc metals. Both models give poor results for thermal expansion, which can be traced to rapid softening of transverse phonon frequencies with increasing lattice parameter. We identify the power law decay of the electron density with distance assumed by the model as the main cause of this behaviour and show that an exponentially decaying form of charge density improves the results significantly. Results for Sutton-Chen and our improved version of Sutton-Chen models are compared for four fcc metals: Cu, Ag, Au and Pt. The calculated properties are the phonon spectra, thermal expansion coefficient, isobaric heat capacity, adiabatic and isothermal bulk moduli, atomic root-mean-square displacement and Gr\"{u}neisen parameter. For the sake of comparison we have also considered two other models where the distance-dependence of the charge density is an exponential multiplied by polynomials. None of these models exhibits the instability against thermal expansion (premature melting) as shown by the Sutton-Chen model. We also present results obtained via pure pair potential models, in order to identify advantages and disadvantages of methods used to obtain the parameters of these potentials.

Dynamic and Thermodynamic Structure of the Monsoon Boundary Layer over India

The objective of the study is to examine the dynamic and thermodynamic structure and the variations that occur in the surface layer during the pre-monsoon, onset and post-monsoon periods over the Indian region. The variations caused during the occurrence of micro and mesoscale systems, structure and variation in the marine boundary layer over the Indian region is also investigated. The drag coefficient computed indirectly also shows variation during various seasons. The thermodynamic structure of the atmosphere shows variation during the various seasons. The onset monsoon causes lowering of the Lifting Condensation Levels. The outcome of the study is expected to provide a better understanding of the structure and variations in the boundary layer over India, which is useful for many applications especially for numerical modeling studies.

Thermodynamic Parameters of Transfer of 2-Nitrobenzoic Acid & 3-Nitrobenzoic Acid from Water to Aqueous Salt Solutions

The Setschenow parameters of solubility in salt solutions and the thermodynamic parameters (25·C) of transfer from aqueous solution to aqueous salt solutions for 2-nitrobenzoic acid and 3-nitrobenzoic acid have been reported. The data have been rationalized on the basis of the localized hydrolysis model and the structure breaking action of ions of the electrolytes.

Precursor nanoscale modulations in ferromagnets: Modeling and thermodynamic characterization

Using a Ginzburg-Landau model for the magnetic degrees of freedom with coupling to disorder, we demonstrate through simulations the existence of stripelike magnetic precursors recently observed in Co-Ni-Al alloys above the Curie temperature. We characterize these magnetic modulations by means of the temperature dependence of local magnetization distribution, magnetized volume fraction, and magnetic susceptibility. We also obtain a temperature-disorder strength phase diagram in which a magnetic tweed phase exists in a small region between the paramagnetic and dipolar phases.

EVALUATION OF THERMODYNAMIC PARAMETERS OF THE ATMOSPHERE BY A FORTRAN PROGRAM

Thermodynamic parameters of the atmosphere form part of the input to numerical forecasting models. Usually these parameters are evaluated from a thermodynamic diagram. Here, a technique is developed to evaluate these parameters quickly and accurately using a Fortran program. This technique is tested with four sets of randomly selected data and the results are in agreement with the results from the conventional method. This technique is superior to the conventional method in three respects: more accuracy, less computation time, and evaluation of additional parameters. The computation time for all the parameters on a PC AT 286 machine is II sec. This software, with appropriate modifications, can be used, for verifying various lines on a thermodynamic diagram

Thermodynamic functions of element 105 in neutral and ionized states

The basic thermodynamic functions, the entropy, free energy, and enthalpy, for element 105 (hahnium) in electronic configurations d^3 s^2, d^3 sp, and d^4s^1 and for its +5 ionized state (5f^14) have been calculated as a function of temperature. The data are based on the results of the calculations of the corresponding electronic states of element 105 using the multiconfiguration Dirac-Fock method.

Prediction of some thermodynamic properties of selected compounds of element 104

A set of parametrized equations has been published by Bratsch and Lagowski for calculating thermodynamic properties of the lanthanides, actinides, element 104, and certainrelated elements. Since these equations were applied to element 104, new values for the first four ionization energies and radii of the ions of charge +1, +2, +3, and +4 have been calculated for this element. The parametrized equations are used here with these new values to calculate some thermodynamic properties of element 104.