Computational simulations of latent heat thermal energy storage systems - with innovative and first-principles based simulation for the underlying unsteady melting (and solidification) processes

This thesis develops an effective modeling and simulation procedure for a specific thermal energy storage system commonly used and recommended for various applications (such as an auxiliary energy storage system for solar heating based Rankine cycle power plant). This thermal energy storage system transfers heat from a hot fluid (termed as heat transfer fluid - HTF) flowing in a tube to the surrounding phase change material (PCM). Through unsteady melting or freezing process, the PCM absorbs or releases thermal energy in the form of latent heat. Both scientific and engineering information is obtained by the proposed first-principle based modeling and simulation procedure. On the scientific side, the approach accurately tracks the moving melt-front (modeled as a sharp liquid-solid interface) and provides all necessary information about the time-varying heat-flow rates, temperature profiles, stored thermal energy, etc. On the engineering side, the proposed approach is unique in its ability to accurately solve – both individually and collectively – all the conjugate unsteady heat transfer problems for each of the components of the thermal storage system. This yields critical system level information on the various time-varying effectiveness and efficiency parameters for the thermal storage system.

Proteins in silico-modeling and sampling

Proteins are linear chain molecules made out of amino acids. Only when they fold to their native states, they become functional. This dissertation aims to model the solvent (environment) effect and to develop & implement enhanced sampling methods that enable a reliable study of the protein folding problem in silico. We have developed an enhanced solvation model based on the solution to the Poisson-Boltzmann equation in order to describe the solvent effect. Following the quantum mechanical Polarizable Continuum Model (PCM), we decomposed net solvation free energy into three physical terms– Polarization, Dispersion and Cavitation. All the terms were implemented, analyzed and parametrized individually to obtain a high level of accuracy. In order to describe the thermodynamics of proteins, their conformational space needs to be sampled thoroughly. Simulations of proteins are hampered by slow relaxation due to their rugged free-energy landscape, with the barriers between minima being higher than the thermal energy at physiological temperatures. In order to overcome this problem a number of approaches have been proposed of which replica exchange method (REM) is the most popular. In this dissertation we describe a new variant of canonical replica exchange method in the context of molecular dynamic simulation. The advantage of this new method is the easily tunable high acceptance rate for the replica exchange. We call our method Microcanonical Replica Exchange Molecular Dynamic (MREMD). We have described the theoretical frame work, comment on its actual implementation, and its application to Trp-cage mini-protein in implicit solvent. We have been able to correctly predict the folding thermodynamics of this protein using our approach.

EXPERIMENTAL INVESTIGATION OF CERTAIN INTERNAL CONDENSING AND BOILING FLOWS: THEIR SENSITIVITY TO PRESSURE FLUCTUATIONS AND HEAT TRANSFER ENHANCEMENTS

Space-based (satellite, scientific probe, space station, etc.) and millimeter – to – microscale (such as are used in high power electronics cooling, weapons cooling in aircraft, etc.) condensers and boilers are shear/pressure driven. They are of increasing interest to system engineers for thermal management because flow boilers and flow condensers offer both high fluid flow-rate-specific heat transfer capacity and very low thermal resistance between the fluid and the heat exchange surface, so large amounts of heat may be removed using reasonably-sized devices without the need for excessive temperature differences. However, flow stability issues and degradation of performance of shear/pressure driven condensers and boilers due to non-desirable flow morphology over large portions of their lengths have mostly prevented their use in these applications. This research is part of an ongoing investigation seeking to close the gap between science and engineering by analyzing two key innovations which could help address these problems. First, it is recommended that the condenser and boiler be operated in an innovative flow configuration which provides a non-participating core vapor stream to stabilize the annular flow regime throughout the device length, accomplished in an energy-efficient manner by means of ducted vapor re-circulation. This is demonstrated experimentally. Second, suitable pulsations applied to the vapor entering the condenser or boiler (from the re-circulating vapor stream) greatly reduce the thermal resistance of the already effective annular flow regime. For experiments reported here, application of pulsations increased time-averaged heat-flux up to 900 % at a location within the flow condenser and up to 200 % at a location within the flow boiler, measured at the heat-exchange surface. Traditional fully condensing flows, reported here for comparison purposes, show similar heat-flux enhancements due to imposed pulsations over a range of frequencies. Shear/pressure driven condensing and boiling flow experiments are carried out in horizontal mm-scale channels with heat exchange through the bottom surface. The sides and top of the flow channel are insulated. The fluid is FC-72 from 3M Corporation.

Detection of pulmonary amylase activity in exhaled breath condensate

Amylase activity in exhaled breath condensate (EBC) is usually interpreted as an indication of oropharyngeal contamination despite the fact that amylase can be found in pulmonary excretions. The aim of this study was to recruit and refine an amylase assay in order to detect amylase activity in any EBC sample and to develop a method to identify EBC samples containing amylase of pulmonary origin. EBC was collected from 40 volunteers with an EcoScreen condenser. Amylase assays and methods to discriminate between oropharyngeal and pulmonary proteins were tested and developed using matched EBC and saliva samples. Our refined 2-chloro-4-nitrophenyl-α-D-maltotriosid (CNP-G3) assay was 40-fold more sensitive than the most sensitive commercial assay and allowed detection of amylase activity in 30 µl of EBC. We developed a dot-blot assay which allowed detection of salivary protein in saliva diluted up to 150 000-fold. By plotting amylase activity against staining intensity we identified a few EBC samples with high amylase activity which were aligned with diluted saliva. We believe that EBC samples aligned with diluted saliva contain amylase activity introduced during EBC collection and that all other EBC samples contain amylase activity of pulmonary origin and are basically free of oropharyngeal protein contamination.

Probabilistic approach to cloud and snow detection on Advanced Very High Resolution Radiometer (AVHRR) imagery

Derivation of probability estimates complementary to geophysical data sets has gained special attention over the last years. Information about a confidence level of provided physical quantities is required to construct an error budget of higher-level products and to correctly interpret final results of a particular analysis. Regarding the generation of products based on satellite data a common input consists of a cloud mask which allows discrimination between surface and cloud signals. Further the surface information is divided between snow and snow-free components. At any step of this discrimination process a misclassification in a cloud/snow mask propagates to higher-level products and may alter their usability. Within this scope a novel probabilistic cloud mask (PCM) algorithm suited for the 1 km × 1 km Advanced Very High Resolution Radiometer (AVHRR) data is proposed which provides three types of probability estimates between: cloudy/clear-sky, cloudy/snow and clear-sky/snow conditions. As opposed to the majority of available techniques which are usually based on the decision-tree approach in the PCM algorithm all spectral, angular and ancillary information is used in a single step to retrieve probability estimates from the precomputed look-up tables (LUTs). Moreover, the issue of derivation of a single threshold value for a spectral test was overcome by the concept of multidimensional information space which is divided into small bins by an extensive set of intervals. The discrimination between snow and ice clouds and detection of broken, thin clouds was enhanced by means of the invariant coordinate system (ICS) transformation. The study area covers a wide range of environmental conditions spanning from Iceland through central Europe to northern parts of Africa which exhibit diverse difficulties for cloud/snow masking algorithms. The retrieved PCM cloud classification was compared to the Polar Platform System (PPS) version 2012 and Moderate Resolution Imaging Spectroradiometer (MODIS) collection 6 cloud masks, SYNOP (surface synoptic observations) weather reports, Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observations (CALIPSO) vertical feature mask version 3 and to MODIS collection 5 snow mask. The outcomes of conducted analyses proved fine detection skills of the PCM method with results comparable to or better than the reference PPS algorithm.

Novel heating/cooling stage designed for fluid inclusion microthermometry of large stalagmite sections

Liquid–vapour homogenisation temperatures of fluid inclusions in stalagmites are used for quantitative temperature reconstructions in paleoclimate research. Specifically for this application, we have developed a novel heating/cooling stage that can be operated with large stalagmite sections of up to 17 × 35 mm2 to simplify and improve the chronological reconstruction of paleotemperature time-series. The stage is designed for use of an oil immersion objective and a high-NA condenser front lens to obtain high-resolution images for bubble radius measurements. The temperature accuracy of the stage is better than ± 0.1 °C with a precision (reproducibility) of ± 0.02 °C.

Slow conduction in mixed cultured strands of primary ventricular cells and stem cell-derived cardiomyocytes

Modern concepts for the treatment of myocardial diseases focus on novel cell therapeutic strategies involving stem cell-derived cardiomyocytes (SCMs). However, functional integration of SCMs requires similar electrophysiological properties as primary cardiomyocytes (PCMs) and the ability to establish intercellular connections with host myocytes in order to contribute to the electrical and mechanical activity of the heart. The aim of this project was to investigate the properties of cardiac conduction in a co-culture approach using SCMs and PCMs in cultured cell strands. Murine embryonic SCMs were pooled with fetal ventricular cells and seeded in predefined proportions on microelectrode arrays to form patterned strands of mixed cells. Conduction velocity (CV) was measured during steady state pacing. SCM excitability was estimated from action potentials measured in single cells using the patch clamp technique. Experiments were complemented with computer simulations of conduction using a detailed model of cellular architecture in mixed cell strands. CV was significantly lower in strands composed purely of SCMs (5.5 ± 1.5 cm/s, n = 11) as compared to PCMs (34.9 ± 2.9 cm/s, n = 21) at similar refractoriness (100% SCMs: 122 ± 25 ms, n = 9; 100% PCMs: 139 ± 67 ms, n = 14). In mixed strands combining both cell types, CV was higher than in pure SCMs strands, but always lower than in 100% PCM strands. Computer simulations demonstrated that both intercellular coupling and electrical excitability limit CV. These data provide evidence that in cultures of murine ventricular cardiomyocytes, SCMs cannot restore CV to control levels resulting in slow conduction, which may lead to reentry circuits and arrhythmias.

Prognostic Implication of Lymphovascular Invasion Detected by Double Immunostaining for D2-40 and MITF1 in Primary Cutaneous Melanoma.

BACKGROUND Lymphovascular invasion (LVI) is associated with adverse outcomes in primary cutaneous melanoma (PCM). Detection of LVI by hematoxylin and eosin staining alone is 0%-6%, but targeting lymphovascular structures increases the detection rate. OBJECTIVE To examine the prognostic significance of LVI detected by immunostaining for D2-40 and microphthalmia-associated transcription factor 1 (MITF1) in PCM. METHODS The authors retrospectively analyzed 120 PCM samples. We compared the LVI detection rates of immunostaining for D2-40 only (22%), double staining for D2-40 and MITF1 (38%), and hematoxylin and eosin, and examined the association of LVI with clinicopathologic variables and clinical outcomes. RESULTS Immunolabeling with both methods significantly increased the LVI detection rate. Double staining for D2-40 and MITF1 as well as D2-40-detected LVI was significantly associated with increased Breslow thickness, number of mitoses, and sentinel lymph node (SLN) metastasis. D2-40-detected LVI was also associated with ulceration. Although the difference was not significant, double staining for D2-40 and MITF1 allowed for easier detection of LVI than D2-40 alone. LIMITATIONS This study was conducted in a tertiary referral institution; therefore, a referral bias cannot be excluded. CONCLUSIONS Immunolabeling increased detection of LVI in PCM. Because LVI is a positive predictive marker for SLN metastasis, the authors propose using anti-D2-40 and anti-MITF1 in the evaluation of LVI in patients with PCM with a certain risk of SLN metastasis.

Game Theoretic Bidding Strategies for Auctions in Green Electricity Markets

Electricity markets in the United States presently employ an auction mechanism to determine the dispatch of power generation units. In this market design, generators submit bid prices to a regulation agency for review, and the regulator conducts an auction selection in such a way that satisfies electricity demand. Most regulators currently use an auction selection method that minimizes total offer costs ["bid cost minimization" (BCM)] to determine electric dispatch. However, recent literature has shown that this method may not minimize consumer payments, and it has been shown that an alternative selection method that directly minimizes total consumer payments ["payment cost minimization" (PCM)] may benefit social welfare in the long term. The objective of this project is to further investigate the long term benefit of PCM implementation and determine whether it can provide lower costs to consumers. The two auction selection methods are expressed as linear constraint programs and are implemented in an optimization software package. Methodology for game theoretic bidding simulation is developed using EMCAS, a real-time market simulator. Results of a 30-day simulation showed that PCM reduced energy costs for consumers by 12%. However, this result will be cross-checked in the future with two other methods of bid simulation as proposed in this paper.

An Experimental Study of Complex-Offer Auctions: Payment Cost Minimization vs. Offer Cost Minimization

A Payment Cost Minimization (PCM) auction has been proposed as an alternative to the Offer Cost Minimization (OCM) auction to be used in wholesale electric power markets with the intention to lower the procurement cost of electricity. Efficiency concerns about this proposal have relied on the assumption of true production cost revelation. Using an experimental approach, I compare the two auctions, strictly controlling for the level of unilateral market power. A specific feature of these complex-offer auctions is that the sellers submit not only the quantities and the minimum prices at which they are willing to sell, but also the start-up fees that are designed to reimburse the fixed start-up costs of the generation plants. I find that both auctions result in start-up fees that are significantly higher than the start-up costs. Overall, the two auctions perform similarly in terms of procurement cost and efficiency. Surprisingly, I do not find a substantial difference between less market power and more market power designs. Both designs result in similar inefficiencies and equally higher procurement costs over the competitive prediction. The PCM auction tends to have lower price volatility than the OCM auction when the market power is minimal but this property vanishes in the designs with market power. These findings lead me to conclude that both the PCM and the OCM auctions do not belong to the class of truth revealing mechanisms and do not easily elicit competitive behavior.

Age determination of sediment cores off western Spitsbergen

The Late Weichselian glacial history of the continental shelf off western Spitsbergen is discussed, based on acoustic sub-bottom records and sediment cores. The outer part of Isfjorden and the inner shelf to the west of this fjord are characterized by a thin veneer (10-20 m) of glacigenic sediments and absence of ice-marginal features. Towards the outer shelf the sediment thickness increases significantly, and exceeds 500 m at the shelf edge. Possible moraine complexes were identified in this outer part. Sediment cores from the inner shelf sampled a firm diamicton, interpreted as till, beneath soft glaciomarine sediments. Radiocarbon dates on shells from the clay resting directly on the till, suggest an age of around 12,500 yrs B.P. for the base of the marine sequence. We argue that grounded ice covered the sites shortly before. In contrast to suggestions that the fjords and coast were partly ice free during the Late Weischselian, we conclude that the ice must have reached out onto the continental shelf.

Sm-Nd isotopic investigation for rocks of Prince Charles Mountains, East Antarctica

The Prince Charles Mountains have been subject to extensive geological and geophysical investigations by former Soviet, Russian and Australian scientists from the early 1970s. In this paper we summarise, and review available geological and isotopic data, and report results of new isotopic studies (Sm-Nd, Pb-Pb, and U-Pb SHRIMP analyses); field geological data obtained during the PCMEGA 2002/2003 are utilised. The structure of the region is described in terms of four tectonic terranes. Those include Archaean Ruker, Palaeoproterozoic Lambert, Mesoproterozoic Fisher, and Meso- to Neoproterozoic Beaver Terranes. Pan-African activities (granite emplacement and probably tectonics) in the Lambert Terrane are reported. We present a summary of the composition of these terranes, discuss their origin and relationships. We also outline the most striking geological features, and problems, and try to draw attention to those rocks and regional geological features which are important in understanding the composition and evolution of the PCM and might suggest targets for further investigations.

Age determination of sediment cores from the western Svalbard margin

The deglaciation of the continental shelf to the west of Spitsbergen and the main fjord, Isfjorden, is discussed based on sub-bottom seismic records and sediment cores. The sea floor on the shelf to the west of Isfjorden is underlain by less than 2 m of glaciomarine sediments over a firm diamicton interpreted as till. In central Isfjorden up to 10 m of deglaciation sediments were recorded, whereas in cores from the innermost tributary, Billefjorden, less than a meter of ice proximal sediments was recognized between the till and the 'normal' Holocene marine sediments. We conclude that the Barents Sea Ice Sheet terminated along the shelf break during the Late Weichselian glacial maximum. Radiocarbon dates from the glaciomarine sediments above the till indicate a stepwise deglaciation. Apparently the ice front retreated from the outermost shelf around 14.8 ka. A dramatic increase in the flux of line-grained glaciomarine sediments around 13 ka is assumed to reflect increased melting and/or current activity due to a climatic warming. This second stage of deglaciation was interrupted by a glacial readvance culminating on the mid-shelf area shortly after 12.4 ka. The glacial readvance, which is correlated with a simultaneous readvance of the Fennoscandian ice sheet along the western coast of Norway, is attributed to the so-called 'Older Dryas' cooling event in the North Atlantic region. Following this glacial readvance the outer part of Isfjorden became rapidly deglaciated around 12.3 ka. During the Younger Dryas the inner fjord branches were occupied by large outlet glaciers and possibly the ice front terminated far out in the main fjord. The remnants of the Barents Sea Ice Sheet melted quickly away as a response to the Holocene warming around 10 ka.

Geochemistry and age determination of zircons in Cretaceous to Quaternary sediments

In central Antarctica, drainage today and earlier back to the Paleozoic radiates from the Gamburtsev Subglacial Mountains (GSM). Proximal to the GSM past the Permian-Triassic fluvial sandstones in the Prince Charles Mountains (PCM) are Cretaceous, Eocene, and Pleistocene sediment in Prydz Bay (ODP741, 1166, and 1167) and pre-Holocene sediment in AM04 beneath the Amery Ice Shelf. We analysed detrital zircons for U-Pb ages, Hf-isotope compositions, and trace elements to determine the age, rock type, source of the host magma, and "crustal" model age (T(C)DM). These samples, together with others downslope from the GSM and the Vostok Subglacial Highlands (VSH), define major clusters of detrital zircons interpreted as coming from (1) 700 to 460 Ma mafic granitoids and alkaline rock, epsilon-Hf 9 to -28, signifying derivation 2.5 to 1.3 Ga from fertile and recycled crust, and (2) 1200-900 Ma mafic granitoids and alkaline rock, epsilon-Hf 11 to -28, signifying derivation 1.8 to 1.3 Ga from fertile and recycled crust. Minor clusters extend to 3350 Ma. Similar detrital zircons in Permian-Triassic, Ordovician, Cambrian, and Neoproterozoic sandstones located along the PaleoPacific margin of East Antarctica and southeast Australia further downslope from central Antarctica reflect the upslope GSM-VSH nucleus of the central Antarctic provenance as a complex of 1200-900 Ma (Grenville) mafic granitoids and alkaline rocks and older rocks embedded in 700-460 Ma (Pan-Gondwanaland) fold belts. The wider central Antarctic provenance (CAP) is tentatively divided into a central sector with negative ?Hf in its 1200-900 Ma rocks bounded on either side by positive epsilon-Hf. The high ground of the GSM-VSH in the Permian and later to the present day is attributed to crustal shortening by far-field stress during the 320 Ma mid-Carboniferous collision of Gondwanaland and Laurussia. Earlier uplifts in the ~500 Ma Cambrian possibly followed the 700-500 Ma assembly of Gondwanaland, and in the Neoproterozoic the 1000-900 Ma collisional events in the Eastern Ghats-Rayner Province at the end of the 1300-1000 Ma assembly of Rodinia.

Highly stereoselective decarboxylation of (+)-1-Bromo-1-chloro-2,2,2-trifluoropropanoic acid gives (+)-1-Bromo-1-chloro-2,2,2-trifluoroethane ((+)-Halothane) with retention of configuration

The absolute configuration of the title acid (2) has been determined to be S by X-ray crystallography. Thus, decarboxylation of 2 produces (S)-(+)-halothane with 99% retention of configuration. This behavior is compared to other stereoselective decarboxylation reactions of ?-haloacids from the literature that also give high degrees of retention of configuration when in the form of their quaternary ammonium salts, which contain one proton. The proton of the ammonium salt is necessary to protonate the anionic intermediate formed from decarboxylation. In the absence of this relatively acidic proton, we had previously found that using triethylene glycol (TEG) as both solvent and proton source for the decarboxylation reaction of acid 2 caused poor stereoselectivity. This was in contrast to 1,2,2,2-tetrafluoro-1-methoxypropionic acid (6), which showed a high degree of retention of configuration in TEG. To rationalize this differing behavior we report DFT studies at PCM-B3LYP/6-31++G** level of theory (the results were additionally confirmed with 6-311++G** and aug-cc-pVDZ basis sets). The energy barrier to inversion of configuration of the anionic reaction intermediate of acid 2 (11) is 10.23 kcal/mol. However, we find that the anionic intermediate from acid 6 (10) would rather undergo ?-elimination instead of inversion of configuration. Thus the planar transition state required for inversion of configuration is never reached, regardless of the rate of proton transfer to the anion.