Second law analysis and simulation techniques for the energy optimization of buildings

The research activity described in this thesis is focused mainly on the study of finite-element techniques applied to thermo-fluid dynamic problems of plant components and on the study of dynamic simulation techniques applied to integrated building design in order to enhance the energy performance of the building. The first part of this doctorate thesis is a broad dissertation on second law analysis of thermodynamic processes with the purpose of including the issue of the energy efficiency of buildings within a wider cultural context which is usually not considered by professionals in the energy sector. In particular, the first chapter includes, a rigorous scheme for the deduction of the expressions for molar exergy and molar flow exergy of pure chemical fuels. The study shows that molar exergy and molar flow exergy coincide when the temperature and pressure of the fuel are equal to those of the environment in which the combustion reaction takes place. A simple method to determine the Gibbs free energy for non-standard values of the temperature and pressure of the environment is then clarified. For hydrogen, carbon dioxide, and several hydrocarbons, the dependence of the molar exergy on the temperature and relative humidity of the environment is reported, together with an evaluation of molar exergy and molar flow exergy when the temperature and pressure of the fuel are different from those of the environment. As an application of second law analysis, a comparison of the thermodynamic efficiency of a condensing boiler and of a heat pump is also reported. The second chapter presents a study of borehole heat exchangers, that is, a polyethylene piping network buried in the soil which allows a ground-coupled heat pump to exchange heat with the ground. After a brief overview of low-enthalpy geothermal plants, an apparatus designed and assembled by the author to carry out thermal response tests is presented. Data obtained by means of in situ thermal response tests are reported and evaluated by means of a finite-element simulation method, implemented through the software package COMSOL Multyphysics. The simulation method allows the determination of the precise value of the effective thermal properties of the ground and of the grout, which are essential for the design of borehole heat exchangers. In addition to the study of a single plant component, namely the borehole heat exchanger, in the third chapter is presented a thorough process for the plant design of a zero carbon building complex. The plant is composed of: 1) a ground-coupled heat pump system for space heating and cooling, with electricity supplied by photovoltaic solar collectors; 2) air dehumidifiers; 3) thermal solar collectors to match 70% of domestic hot water energy use, and a wood pellet boiler for the remaining domestic hot water energy use and for exceptional winter peaks. This chapter includes the design methodology adopted: 1) dynamic simulation of the building complex with the software package TRNSYS for evaluating the energy requirements of the building complex; 2) ground-coupled heat pumps modelled by means of TRNSYS; and 3) evaluation of the total length of the borehole heat exchanger by an iterative method developed by the author. An economic feasibility and an exergy analysis of the proposed plant, compared with two other plants, are reported. The exergy analysis was performed by considering the embodied energy of the components of each plant and the exergy loss during the functioning of the plants.

Article Accurate Experimental Values for the Free Energies of Hydration of H+, OH-, and H3O+

Accurate experimental values for the free energies of hydration, or the free energies of solvation, of the H+, OH-, and H3O+ ions are of fundamental importance. By use of the most accurate value for the free energy of solvation of H+, the known value for the free energy of solvation of water, and the known values for the gas phase and aqueous phase deprotonation of water, the corresponding experimental free energy of solvation for OH- is −106.4 ± 0.5 kcal/mol. Similarly, by use of the known values for ΔGf 0 for H3O, H2O+, and OH-, the known values for ΔGs for H+ and OH-, and the known value for the aqueous phase autoionization of water, we obtain an experimental free energy of solvation value for H3O+ of −103.4 ± 0.5 kcal/mol. These values are in excellent agreement with the commonly accepted values and with the value for ΔGs(OH-) obtained from embedding clusters of OH-(H2O)n in a dielectric continuum.

Global Search for Minimum Energy (H2O)n Clusters, n = 3−5

The Gaussian-3 (G3) model chemistry method has been used to calculate the relative ΔG° values for all possible conformers of neutral clusters of water, (H2O)n, where n = 3−5. A complete 12-fold conformational search around each hydrogen bond produced 144, 1728, and 20 736 initial starting structures of the water trimer, tetramer, and pentamer. These structures were optimized with PM3, followed by HF/6-31G* optimization, and then with the G3 model chemistry. Only two trimers are present on the G3 potential energy hypersurface. We identified 5 tetramers and 10 pentamers on the potential energy and free-energy hypersurfaces at 298 K. None of these 17 structures were linear; all linear starting models folded into cyclic or three-dimensional structures. The cyclic pentamer is the most stable isomer at 298 K. On the basis of this and previous studies, we expect the cyclic tetramers and pentamers to be the most significant cyclic water clusters in the atmosphere.

Use of intrinsic binding energy for catalysis by an RNA enzyme

The contribution of several individual ribozyme⋅substrate base pairs to binding and catalysis has been investigated using hammerhead ribozyme substrates that were truncated at their 3′ or 5′ ends. The base pairs at positions 1.1–2.1 and 15.2–16.2, which flank the conserved core, each contribute 104-fold in the chemical step, without affecting substrate binding. In contrast, base pairs distal to the core contribute to substrate binding but have no effect on the chemical step. These results suggest a “fraying model” in which each ribozyme⋅substrate helix can exist in either an unpaired (“open”) state or a helical (“closed”) state, with the closed state required for catalysis. The base pairs directly adjacent to the conserved core contribute to catalysis by allowing the closed state to form. Once the number of base pairs is sufficient to ensure that the closed helical state predominates, additional residues provide stabilization of the helix, and therefore increase binding, but have no further effect on the chemical step. Remarkably, the >5 kcal/mol free energy contribution to catalysis from each of the internal base pairs is considerably greater than the free energy expected for formation of a base pair. It is suggested that this unusually large energetic contribution arises because free energy that is typically lost in constraining residues within a base pair is expressed in the transition state, where it is used for positioning. This extends the concept of “intrinsic binding energy” from protein to RNA enzymes, suggesting that intrinsic binding energy is a fundamental feature of biological catalysis.

Equilibrium constants and free energies in unfolding of proteins in urea solutions

A novel thermodynamic approach to the reversible unfolding of proteins in aqueous urea solutions has been developed based on the premise that urea ligands are bound cooperatively to the macromolecule. When successive stoichiometric binding constants have values larger than expected from statistical effects, an equation for moles of bound urea can be derived that contains imaginary terms. For a very steep unfolding curve, one can then show that the fraction of protein unfolded, B̄, depends on the square of the urea concentration, U, and is given by \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} \begin{equation*}\bar {B}=\frac{{\mathit{A}}^{{\mathit{2}}}_{{\mathit{1}}}{\mathit{e}}^{{\mathrm{{\lambda}}}n\bar {B}}{\mathit{U}}^{{\mathit{2}}}}{{\mathrm{1\hspace{.167em}+\hspace{.167em}}}{\mathit{A}}^{{\mathrm{2}}}_{{\mathrm{1}}}{\mathit{e}}^{{\mathrm{{\lambda}}}\bar {B}}{\mathit{U}}^{{\mathrm{2}}}}{\mathrm{.}}\end{equation*}\end{document} A12 is the binding constant as B̄→ 0, and λ is a parameter that reflects the augmentation in affinities of protein for urea as the moles bound increases to the saturation number, n. This equation provides an analytic expression that reproduces the unfolding curve with good precision, suggests a simple linear graphical procedure for evaluating A12 and λ, and leads to the appropriate standard free energy changes. The calculated ΔG° values reflect the coupling of urea binding with unfolding of the protein. Some possible implications of this analysis to protein folding in vivo are described.

Recent improvements in prediction of protein structure by global optimization of a potential energy function

Recent improvements of a hierarchical ab initio or de novo approach for predicting both α and β structures of proteins are described. The united-residue energy function used in this procedure includes multibody interactions from a cumulant expansion of the free energy of polypeptide chains, with their relative weights determined by Z-score optimization. The critical initial stage of the hierarchical procedure involves a search of conformational space by the conformational space annealing (CSA) method, followed by optimization of an all-atom model. The procedure was assessed in a recent blind test of protein structure prediction (CASP4). The resulting lowest-energy structures of the target proteins (ranging in size from 70 to 244 residues) agreed with the experimental structures in many respects. The entire experimental structure of a cyclic α-helical protein of 70 residues was predicted to within 4.3 Å α-carbon (Cα) rms deviation (rmsd) whereas, for other α-helical proteins, fragments of roughly 60 residues were predicted to within 6.0 Å Cα rmsd. Whereas β structures can now be predicted with the new procedure, the success rate for α/β- and β-proteins is lower than that for α-proteins at present. For the β portions of α/β structures, the Cα rmsd's are less than 6.0 Å for contiguous fragments of 30–40 residues; for one target, three fragments (of length 10, 23, and 28 residues, respectively) formed a compact part of the tertiary structure with a Cα rmsd less than 6.0 Å. Overall, these results constitute an important step toward the ab initio prediction of protein structure solely from the amino acid sequence.

DNA-strand exchange promoted by RecA protein in the absence of ATP: implications for the mechanism of energy transduction in protein-promoted nucleic acid transactions.

DNA-strand exchange promoted by Escherichia coli RecA protein normally requires the presence of ATP and is accompanied by ATP hydrolysis, thereby implying a need for ATP hydrolysis. Previously, ATP hydrolysis was shown not to be required; here we demonstrate furthermore that a nucleoside triphosphate cofactor is not required for DNA-strand exchange. A gratuitous allosteric effector consisting of the noncovalent complex of ADP and aluminum fluoride, ADP.AIF4-, can both induce the high-affinity DNA-binding state of RecA protein and support the homologous pairing and exchange of up to 800-900 bp of DNA. These results demonstrate that induction of the functionally active, high-affinity DNA-binding state of RecA protein is needed for RecA protein-promoted DNA-strand exchange and that there is no requirement for a high-energy nucleotide cofactor for the exchange of DNA strands. Consequently, the free energy needed to activate the DNA substrates for DNA-strand exchange is not derived from ATP hydrolysis. Instead, the needed free energy is derived from ligand binding and is transduced to the DNA via the associated ligand-induced structural transitions of the RecA protein-DNA complex; ATP hydrolysis simply destroys the effector ligand. This concept has general applicability to the mechanism of energy transduction by proteins.

Tuning the metal-to-metal charge transfer energy of cyano-bridged dinuclear complexes

The metal-to-metal charge transfer (MMCT) transitions of a series of Class II mixed valence dinuclear complexes bearing cyano bridging ligands may be varied systematically by variations to either the hexacyanometallate(II) donor or Co-III acceptor moieties. Specifically, the new dinuclear species trans-[(LCoNCFe)-Co-14S(CN)(5)](-) (L-14S = 6-methyl-1,11-diaza-4,8-dithia- cyclotetradecane-6-amine) and trans-[(LCoNCRu)-Co-14(CN)(5)]-(L-14 = 6-methyl-1,4,8,11-tetraazacyclotetradecane-6-amine) have been prepared and their spectroscopic and electrochemical properties are compared with the relative trans-[(LCoNCFe)-Co-14(CN)(5)](-). The crystal structures of Na{trans-[(LCoNCFe)-Co-14S(CN)(5)]}.51/2H(2)O.1/2EtOH, Na{trans-[(LCoNCRu)-Co-14(CN)(5)]}.3H(2)O and Na{trans-[(LCoNCRu)-Co-14(CN)(5)]}.8H(2)O are also reported. The ensuing changes to the MMCT energy have been examined within the framework of Hush theory, and it was found that the free energy change between the redox isomers was the dominant effect in altering the energy of the MMCT transition.

Can surface energy measurements predict the impact of catalyst hydrophobicity upon fatty acid esterification over sulfonic acid functionalised periodic mesoporous organosilicas?

Sulfonic acid functionalised periodic mesoporous organosilicas (PrSO3 H-PMOs) with tunable hydrophobicity were synthesised via a surfactant-templating route, and characterised by porosimetry, TEM, XRD, XPS, inverse gas chromatography (IGC) and ammonia pulse chemisorption. IGC reveals that incorporation of ethyl or benzyl moieties into a mesoporous SBA-15 silica framework significantly increases the non-specific dispersive surface energy of adsorption for alkane adsorption, while decreasing the free energy of adsorption of methanol, reflecting increased surface hydrophobicity. The non-specific dispersive surface energy of adsorption of PMO-SO3H materials is strongly correlated with their activity towards palmitic acid esterification with methanol, demonstrating the power of IGC as an analytical tool for identifying promising solid acid catalysts for the esterification of free fatty acids. A new parameter [-ΔGCNP-P], defined as the per carbon difference in Gibbs free energy of adsorption between alkane and polar probe molecules, provides a simple predictor of surface hydrophobicity and corresponding catalyst activity in fatty acid esterification. © 2014 Elsevier B.V.

Interactions of phenyldithioesters with gold nanoparticles (AuNPs): Implications for AuNP functionalization and molecular barcoding of AuNP assemblies

The interactions of phenyldithioesters with gold nanoparticles (AuNPs) have been studied by monitoring changes in the surface plasmon resonance (SPR), depolarised light scattering, and surface enhanced Raman spectroscopy (SERS). Changes in the SPR indicated that an AuNP-phenyldithioester charge transfer complex forms in equilibrium with free AuNPs and phenyldithioester. Analysis of the Langmuir binding isotherms indicated that the equilibrium adsorption constant, Kads, was 2.3 ± 0.1 × 106 M−1, which corresponded to a free energy of adsorption of 36 ± 1 kJ mol−1. These values are comparable to those reported for interactions of aryl thiols with gold and are of a similar order of magnitude to moderate hydrogen bonding interactions. This has significant implications in the application of phenyldithioesters for the functionalization of AuNPs. The SERS results indicated that the phenyldithioesters interact with AuNPs through the C═S bond, and the molecules do not disassociate upon adsorption to the AuNPs. The SERS spectra are dominated by the portions of the molecule that dominate the charge transfer complex with the AuNPs. The significance of this in relation to the use of phenyldithioesters for molecular barcoding of nanoparticle assemblies is discussed.

Dynamics of transformation and fragmentation of composite liquid nano-particles

Recent research on particle size distributions and particle concentrations near a busy road cannot be explained by the conventional mechanisms for particle evolution of combustion aerosols. Specifically they appear to be inadequate to explain the experimental observations of particle transformation and the evolution of the total number concentration. This resulted in the development of a new mechanism based on their thermal fragmentation, for the evolution of combustion aerosol nano-particles. A complex and comprehensive pattern of evolution of combustion aerosols, involving particle fragmentation, was then proposed and justified. In that model it was suggested that thermal fragmentation occurs in aggregates of primary particles each of which contains a solid graphite/carbon core surrounded by volatile molecules bonded to the core by strong covalent bonds. Due to the presence of strong covalent bonds between the core and the volatile (frill) molecules, such primary composite particles can be regarded as solid, despite the presence of significant (possibly, dominant) volatile component. Fragmentation occurs when weak van der Waals forces between such primary particles are overcome by their thermal (Brownian) motion. In this work, the accepted concept of thermal fragmentation is advanced to determine whether fragmentation is likely in liquid composite nano-particles. It has been demonstrated that at least at some stages of evolution, combustion aerosols contain a large number of composite liquid particles containing presumably several components such as water, oil, volatile compounds, and minerals. It is possible that such composite liquid particles may also experience thermal fragmentation and thus contribute to, for example, the evolution of the total number concentration as a function of distance from the source. Therefore, the aim of this project is to examine theoretically the possibility of thermal fragmentation of composite liquid nano-particles consisting of immiscible liquid v components. The specific focus is on ternary systems which include two immiscible liquid droplets surrounded by another medium (e.g., air). The analysis shows that three different structures are possible, the complete encapsulation of one liquid by the other, partial encapsulation of the two liquids in a composite particle, and the two droplets separated from each other. The probability of thermal fragmentation of two coagulated liquid droplets is discussed and examined for different volumes of the immiscible fluids in a composite liquid particle and their surface and interfacial tensions through the determination of the Gibbs free energy difference between the coagulated and fragmented states, and comparison of this energy difference with the typical thermal energy kT. The analysis reveals that fragmentation was found to be much more likely for a partially encapsulated particle than a completely encapsulated particle. In particular, it was found that thermal fragmentation was much more likely when the volume ratio of the two liquid droplets that constitute the composite particle are very different. Conversely, when the two liquid droplets are of similar volumes, the probability of thermal fragmentation is small. It is also demonstrated that the Gibbs free energy difference between the coagulated and fragmented states is not the only important factor determining the probability of thermal fragmentation of composite liquid particles. The second essential factor is the actual structure of the composite particle. It is shown that the probability of thermal fragmentation is also strongly dependent on the distance that each of the liquid droplets should travel to reach the fragmented state. In particular, if this distance is larger than the mean free path for the considered droplets in the air, the probability of thermal fragmentation should be negligible. In particular, it follows form here that fragmentation of the composite particle in the state with complete encapsulation is highly unlikely because of the larger distance that the two droplets must travel in order to separate. The analysis of composite liquid particles with the interfacial parameters that are expected in combustion aerosols demonstrates that thermal fragmentation of these vi particles may occur, and this mechanism may play a role in the evolution of combustion aerosols. Conditions for thermal fragmentation to play a significant role (for aerosol particles other than those from motor vehicle exhaust) are determined and examined theoretically. Conditions for spontaneous transformation between the states of composite particles with complete and partial encapsulation are also examined, demonstrating the possibility of such transformation in combustion aerosols. Indeed it was shown that for some typical components found in aerosols that transformation could take place on time scales less than 20 s. The analysis showed that factors that influenced surface and interfacial tension played an important role in this transformation process. It is suggested that such transformation may, for example, result in a delayed evaporation of composite particles with significant water component, leading to observable effects in evolution of combustion aerosols (including possible local humidity maximums near a source, such as a busy road). The obtained results will be important for further development and understanding of aerosol physics and technologies, including combustion aerosols and their evolution near a source.

Effect of Fe3O4 addition on removal of ammonium by zeolite NaA

Magnetic zeolite NaA with different Fe3O4 loadings was prepared by hydrothermal synthesis based on metakaolin and Fe3O4. The effect of added Fe3O4 on the removal of ammonium by zeolite NaA was investigated by varying the Fe3O4 loading, pH, adsorption temperature, initial concentration, adsorption time. Langmuir, Freundlich, and pseudo-second-order modeling were used to describe the nature and mechanism of ammonium ion exchange using both zeolite and magnetic zeolite. Thermodynamic parameters such as change in Gibbs free energy, enthalpy and entropy were calculated. The results show that all the selected factors affect the ammonium ion exchange by zeolite and magnetic zeolite, however, the added Fe3O4 apparently does not affect the ion exchange performance of zeolite to the ammonium ion. Freundlich model provides a better description of the adsorption process than Langmuir model. Moreover, kinetic analysis indicates the exchange of ammonium on the two materials follows a pseudosecond-order model. Thermodynamic analysis makes it clear that the adsorption process of ammonium is spontaneous and exothermic. Regardless of kinetic or thermodynamic analysis, all the results suggest that no considerable effect on the adsorption of the ammonium ion by zeolite is found after the addition of Fe3O4. According to the results, magnetic zeolite NaA can be used for the removal of ammonium due to the good adsorption performance and easy separation method from aqueous solution.

In situ imaging of catalytic etching on silver during methanol oxidation conditions by environmental scanning electron microscopy

Polycrystalline silver is used to catalytically oxidise methanol to formaldehyde. This paper reports the results of extensive investigations involving the use of environmental scanning electron microscopy (ESEM) to monitor structural changes in silver during simulated industrial reaction conditions. The interaction of oxygen, nitrogen, and water, either singly or in combination, with a silver catalyst at temperatures up to 973 K resulted in the appearance of a reconstructed silver surface. More spectacular was the effect an oxygen/methanol mixture had on the silver morphology. At a temperature of ca. 713 K pinholes were created in the vicinity of defects as a consequence of subsurface explosions. These holes gradually increased in size and large platelet features were created. Elevation of the catalyst temperature to 843 K facilitated the wholescale oxygen induced restructuring of the entire silver surface. Methanol reacted with subsurface oxygen to produce subsurface hydroxyl species which ultimately formed water in the subsurface layers of silver. The resultant hydrostatic pressure forced the silver surface to adopt a "hill and valley" conformation in order to minimise the surface free energy. Upon approaching typical industrial operating conditions widespread explosions occurred on the catalyst and it was also apparent that the silver surface was extremely mobile under the applied conditions. The interaction of methanol alone with silver resulted in the initial formation of pinholes primarily in the vicinity of defects, due to reaction with oxygen species incorporated in the catalyst during electrochemical synthesis. However, dramatic reduction in the hole concentration with time occurred as all the available oxygen became consumed. A remarkable correlation between formaldehyde production and hole concentration was found.

Development of a Solar Assisted Heat Pump Desalination System

In view of the growing global demand for energy and concern expressed for environmental degradation, a clean and "free" energy source, such as solar energy, has been receiving greater attention in recent years for various applications using different techniques. The Direct Expansion Solar Assisted Heat Pump (DX-SAHP) principle is one of the most promising techniques as it makes use of both solar and ambient energy. As the system has capability to function at low temperatures, it has the potential to operate at night in the tropics. The system utilizes multi-effect distillation (MED) principle for the conversion of seawater to fresh water. An experimental setup of the DX-SAHP desalination system has been built at the Department of Mechanical Engineering, National University of Singapore (NUS). This system uses two types of flat-plate solar collectors. One is called evaporator-collector, where no glazing is used, and the efficiency varies between 80 and 90%. The other type of collector is single-glazed, where the maximum efficiency is about 60%, and it is used for feed water heating. For the heat pump cycle, refrigerant R134a is used. The present study provides a comprehensive analyses and performance evaluation of this system under different operating and meteorological conditions of Singapore. The Coefficient of Performance (COP) of the heat pump system reached a maximum value of 10. For a single effect of desalination, the system shows a Performance Ratio (PR) of around 1.3.

Adsorption thermodynamics and kinetics study for the self-assembly of 1,8,15,22-tetraaminophthalocyanatocobalt(II) on glassy carbon surface

Spontaneous adsorption of 1,8,15,22-tetraaminophthalocyanatocobalt(II) (4α-CoIITAPc) on glassy carbon (GC) electrode leads to the formation of a stable self-assembled monolayer (SAM). Since the SAM of 4α-CoIITAPc is redox active, its adsorption on GC electrode was followed by cyclic voltammetry. SAM of 4α-CoIITAPc on GC electrode shows two pairs of well-defined redox peaks corresponding to CoIII/CoII and CoIIIPc−1/CoIIIPc−2. The surface coverage (Γ) value, calculated by integrating the charge under CoII oxidation, was used to study the adsorption thermodynamics and kinetics of 4α-CoIITAPc on GC surface. Cyclic voltammetric studies show that the adsorption of 4α-CoIITAPc on GC electrode has reached the saturation coverage (Γs) within 3 h. The Γs value for the SAM of 4α-CoIITAPc on GC electrode was found to be 2.37 × 10−10 mol cm−2. Gibbs free energy (ΔGads) and adsorption rate constant (kad) for the adsorption of 4α-CoIITAPc on GC surface were found to be −16.76 kJ mol−1 and 7.1 M−1 s−1, respectively. The possible mechanism for the self-assembly of 4α-CoIITAPc on GC surface is through the addition of nucleophilic amines to the olefinic bond on the GC surface in addition to a meager contribution from π stacking. The contribution of π stacking was confirmed from the adsorption of unsubstituted phthalocyanatocobalt(II) (CoPc) on GC electrode. Raman spectra for the SAM of 4α-CoIITAPc on carbon surface shows strong stretching and breathing bands of Pc macrocycle, pyrrole ring and isoindole ring. Raman and CV studies suggest that 4α-CoIITAPc is adopting nearly a flat orientation or little bit tilted orientation.