Mathematical programming model for heat exchanger design through optimization of partial objectives

Mathematical programming can be used for the optimal design of shell-and-tube heat exchangers (STHEs). This paper proposes a mixed integer non-linear programming (MINLP) model for the design of STHEs, following rigorously the standards of the Tubular Exchanger Manufacturers Association (TEMA). Bell–Delaware Method is used for the shell-side calculations. This approach produces a large and non-convex model that cannot be solved to global optimality with the current state of the art solvers. Notwithstanding, it is proposed to perform a sequential optimization approach of partial objective targets through the division of the problem into sets of related equations that are easier to solve. For each one of these problems a heuristic objective function is selected based on the physical behavior of the problem. The global optimal solution of the original problem cannot be ensured even in the case in which each of the sub-problems is solved to global optimality, but at least a very good solution is always guaranteed. Three cases extracted from the literature were studied. The results showed that in all cases the values obtained using the proposed MINLP model containing multiple objective functions improved the values presented in the literature.

Digital surface model, hillshade and orthophoto of the world's largest fossil oyster reef, links to GeoTIFFs

The world's largest fossil oyster reef, formed by the giant oyster Crassostrea gryphoides and located in Stetten (north of Vienna, Austria) is studied by Harzhauser et al., 2015, 2016; Djuricic et al., 2016. Digital documentation of the unique geological site is provided by terrestrial laser scanning (TLS) at the millimeter scale. Obtaining meaningful results is not merely a matter of data acquisition with a suitable device; it requires proper planning, data management, and postprocessing. Terrestrial laser scanning technology has a high potential for providing precise 3D mapping that serves as the basis for automatic object detection in different scenarios; however, it faces challenges in the presence of large amounts of data and the irregular geometry of an oyster reef. We provide a detailed description of the techniques and strategy used for data collection and processing in Djuricic et al., 2016. The use of laser scanning provided the ability to measure surface points of 46,840 (estimated) shells. They are up to 60-cm-long oyster specimens, and their surfaces are modeled with a high accuracy of 1 mm. In addition to laser scanning measurements, more than 300 photographs were captured, and an orthophoto mosaic was generated with a ground sampling distance (GSD) of 0.5 mm. This high-resolution 3D information and the photographic texture serve as the basis for ongoing and future geological and paleontological analyses. Moreover, they provide unprecedented documentation for conservation issues at a unique natural heritage site.

High-resolution dynamical downscaling products for the North Slope of Alaska and surrounding areas, links to model result files

Climatic changes are most pronounced in northern high latitude regions. Yet, there is a paucity of observational data, both spatially and temporally, such that regional-scale dynamics are not fully captured, limiting our ability to make reliable projections. In this study, a group of dynamical downscaling products were created for the period 1950 to 2100 to better understand climate change and its impacts on hydrology, permafrost, and ecosystems at a resolution suitable for northern Alaska. An ERA-interim reanalysis dataset and the Community Earth System Model (CESM) served as the forcing mechanisms in this dynamical downscaling framework, and the Weather Research & Forecast (WRF) model, embedded with an optimization for the Arctic (Polar WRF), served as the Regional Climate Model (RCM). This downscaled output consists of multiple climatic variables (precipitation, temperature, wind speed, dew point temperature, and surface air pressure) for a 10 km grid spacing at three-hour intervals. The modeling products were evaluated and calibrated using a bias-correction approach. The ERA-interim forced WRF (ERA-WRF) produced reasonable climatic variables as a result, yielding a more closely correlated temperature field than precipitation field when long-term monthly climatology was compared with its forcing and observational data. A linear scaling method then further corrected the bias, based on ERA-interim monthly climatology, and bias-corrected ERA-WRF fields were applied as a reference for calibration of both the historical and the projected CESM forced WRF (CESM-WRF) products. Biases, such as, a cold temperature bias during summer and a warm temperature bias during winter as well as a wet bias for annual precipitation that CESM holds over northern Alaska persisted in CESM-WRF runs. The linear scaling of CESM-WRF eventually produced high-resolution downscaling products for the Alaskan North Slope for hydrological and ecological research, together with the calibrated ERA-WRF run, and its capability extends far beyond that. Other climatic research has been proposed, including exploration of historical and projected climatic extreme events and their possible connections to low-frequency sea-atmospheric oscillations, as well as near-surface permafrost degradation and ice regime shifts of lakes. These dynamically downscaled, bias corrected climatic datasets provide improved spatial and temporal resolution data necessary for ongoing modeling efforts in northern Alaska focused on reconstructing and projecting hydrologic changes, ecosystem processes and responses, and permafrost thermal regimes. The dynamical downscaling methods presented in this study can also be used to create more suitable model input datasets for other sub-regions of the Arctic.

Development of a Conceptual Site Model for Evaluating Seawater Intrusion under Sea Level Rise Scenarios Using Analytical Methods

Island County is located in the Puget Sound of Washington State and includes several islands, the largest of which is Whidbey Island. Central Whidbey Island was chosen as the project site, as residents use groundwater for their water supply and seawater intrusion near the coast is known to contaminate this resource. In 1989, Island County adopted a Saltwater Intrusion Policy and used chloride concentrations in existing wells in order to define and map “risk zones.” In 2005, this method of defining vulnerability was updated with the use of water level elevations in conjunction with chloride concentrations. The result of this work was a revised map of seawater intrusion vulnerability that is currently in use by Island County. This groundwater management strategy is defined as trigger-level management and is largely a reactive tool. In order to evaluate trends in the hydrogeologic processes at the site, including seawater intrusion under sea level rise scenarios, this report presents a workflow where groundwater flow and discharge to the sea are quantified using a revised conceptual site model. The revised conceptual site model used several simplifying assumptions that allow for first-order quantitative predictions of seawater intrusion using analytical methods. Data from water well reports included lithologic and well construction information, static water levels, and aquifer tests for specific capacity. Results from specific capacity tests define the relationship between discharge and drawdown and were input for a modified Theis equation to solve for transmissivity (Arihood, 2009). Components of the conceptual site model were created in ArcGIS and included interpolation of water level elevation, creation of groundwater basins, and the calculation of net recharge and groundwater discharge for each basin. The revised conceptual site model was then used to hypothesize regarding hydrogeologic processes based on observed trends in groundwater flow. Hypotheses used to explain a reduction in aquifer thickness and hydraulic gradient were: (1) A large increase in transmissivity occurring near the coast. (2) The reduced aquifer thickness and hydraulic gradient were the result of seawater intrusion. (3) Data used to create the conceptual site model were insufficient to resolve trends in groundwater flow. For Hypothesis 2, analytical solutions for groundwater flow under Dupuit assumptions were applied in order to evaluate seawater intrusion under projected sea level rise scenarios. Results indicated that a rise in sea level has little impact on the position of a saltwater wedge; however, a reduction in recharge has significant consequences. Future work should evaluate groundwater flow using an expanded monitoring well network and aquifer recharge should be promoted by reducing surface water runoff.

A rapidly evolving secretome builds and patterns a sea shell

Background: Instructions to fabricate mineralized structures with distinct nanoscale architectures, such as seashells and coral and vertebrate skeletons, are encoded in the genomes of a wide variety of animals. In mollusks, the mantle is responsible for the extracellular production of the shell, directing the ordered biomineralization of CaCO3 and the deposition of architectural and color patterns. The evolutionary origins of the ability to synthesize calcified structures across various metazoan taxa remain obscure, with only a small number of protein families identified from molluskan shells. The recent sequencing of a wide range of metazoan genomes coupled with the analysis of gene expression in non-model animals has allowed us to investigate the evolution and process of biomineralization in gastropod mollusks. Results: Here we show that over 25% of the genes expressed in the mantle of the vetigastropod Haliotis asinina encode secreted proteins, indicating that hundreds of proteins are likely to be contributing to shell fabrication and patterning. Almost 85% of the secretome encodes novel proteins; remarkably, only 19% of these have identifiable homologues in the full genome of the patellogastropod Lottia scutum. The spatial expression profiles of mantle genes that belong to the secretome is restricted to discrete mantle zones, with each zone responsible for the fabrication of one of the structural layers of the shell. Patterned expression of a subset of genes along the length of the mantle is indicative of roles in shell ornamentation. For example, Has-sometsuke maps precisely to pigmentation patterns in the shell, providing the first case of a gene product to be involved in molluskan shell pigmentation. We also describe the expression of two novel genes involved in nacre (mother of pearl) deposition. Conclusion: The unexpected complexity and evolvability of this secretome and the modular design of the molluskan mantle enables diversification of shell strength and design, and as such must contribute to the variety of adaptive architectures and colors found in mollusk shells. The composition of this novel mantle-specific secretome suggests that there are significant molecular differences in the ways in which gastropods synthesize their shells.

Shell-side transfer in shell-and-tube heat exchangers

Local mass transfer coefficients were determined by using the electrochemical technique. A simple model of a heat exchanger with segmental nickel tube joined to p.v.c. rods replaced the exchanger tubes. Measurements were made for both no-Ieakage, semi-leakage and total leakage configurations. Baffle-spacings of 47.6 mm, 66.6 mm, 97 mm and 149.2 mm wer studied. Also studied were the overall exchanger pressure drops for each configuration. The comparison of the heat transfer data with this work showed good agreement at high flow rates for the no-leakage case, but the agreement became poor for lower flow rates and leakage configurations. This disagreement was explained by non-analogous driving forces existing in the two systems. The no-leakage data showed length-wise variation of transfer coefficients along the exchanger length. The end compartments showing transfer coefficients inferior by up to 26% compared to tbe internal compartments, depending on Reynolds number. With the introduction of leakage streams this variation however became smaller than the experimental accuracy. A model is outlined to show the characteristic behaviour of individual electrode segments within the compartment. This was able to discriminate between cross and window zones for the no- leakage case, but no such distinction could be made for the leakage case. A flow area was found which, when incorporated in the Reynolds number, enabled the correlation of baffle-cut and baffle-spacing parameters for the no-leakage case . This area is the free flow area determined at the baffle edge. Addition of the leakage area to this flow area resulted in correlation of all commercial leakage geometrical parameters. The procedures used to correlate the pressure drop data from a total of eighteen different configurations on a single curve are also outlined.

Local shell-side coefficients in shell and tube exchanges the use of a mass transfer technique

A diffusion-controlled electrochemical mass transfer technique has been employed in making local measurements of shell-side coefficients in segmentally baffled shell and tube heat exchangers. Corresponding heat transfer data are predicted through the Chilton and Colburn heat and mass transfer analogy. Mass transfer coefficients were measured for baffle spacing lengths of individual tubes in an internal baffle compartment. Shell-side pressure measurements were also made. Baffle compartment average coefficients derived from individual tube coefficients are shown to be in good agreement with reported experimental bundle average heat transfer data for a heat exchanger model of similar geometry. Mass transfer coefficients of individual tubes compare favourably with those obtained previously by another mass transfer technique. Experimental data are reported for a variety of segmental baffle configurations over the shell-side Reynolds number range 100 to 42 000. Baffles with zero clearances were studied at three baffle cuts and two baffle spacings. Baffle geometry is shown to have a large effect on the distribution of tube coefficients within the baffle compartment. Fluid "jetting" is identified with some baffle configurations. No simple characteristic velocity is found to correlate zonal or baffle compartment average mass transfer data for the effect of both baffle cut and baffle spacing. Experiments with baffle clearances typical of commercial heat exchangers are also reported. The effect of leakage streams associated with these baffles is identified. Investigations were extended to double segmental baffles for which no data had previously been published. The similarity in the shell-side characteristics of this baffle arrangement and two parallel single segmental baffle arrangements is demonstrated. A general relationship between the shell-side mass transfer performance and pressure drop was indicated by the data for all the baffle configurations examined.

Local heat transfer coefficient in a baffled shell and tube heat exchanger

Local shell side coefficient measurements in the end conpartments of a model shell and tube heat exchanger have been made using an electrochemical technique. Limited data are also reported far the second compartment. The end compartment average coefficients have been found to be smaller than reported data for a corresponding internal conpartment. The second compartment data. have been shown to lie between those for the end compartments and the reported internal compartment data. Experimental data are reported fcr two port types and two baffle orientations. with data for the case of an inlet compartment impingement baffle also being given . Port type is shown to have a small effect on compartment coefficients, these being largely unaffected. Likewise, the outlet compartment average coefficients are slightly snaller than those for the inlet compartment, with the distribution of individual tube coefficients being similar. Baffle orientation has been shown to have no effect on average coefficients, but the distribution of the data is substantially affected. The use of an impingement baffle in the inlet compartment lessens the efect of baffle orientation on distribution . Recommendations are made for future work.

Projected sequential Gaussian processes: flexible interpolation for large data sets

Recently within the machine learning and spatial statistics communities many papers have explored the potential of reduced rank representations of the covariance matrix, often referred to as projected or fixed rank approaches. In such methods the covariance function of the posterior process is represented by a reduced rank approximation which is chosen such that there is minimal information loss. In this paper a sequential framework for inference in such projected processes is presented, where the observations are considered one at a time. We introduce a C++ library for carrying out such projected, sequential estimation which adds several novel features. In particular we have incorporated the ability to use a generic observation operator, or sensor model, to permit data fusion. We can also cope with a range of observation error characteristics, including non-Gaussian observation errors. Inference for the variogram parameters is based on maximum likelihood estimation. We illustrate the projected sequential method in application to synthetic and real data sets. We discuss the software implementation and suggest possible future extensions.

A projected process kriging algorithm for sensor networks with heterogeneous error characteristics

Large monitoring networks are becoming increasingly common and can generate large datasets from thousands to millions of observations in size, often with high temporal resolution. Processing large datasets using traditional geostatistical methods is prohibitively slow and in real world applications different types of sensor can be found across a monitoring network. Heterogeneities in the error characteristics of different sensors, both in terms of distribution and magnitude, presents problems for generating coherent maps. An assumption in traditional geostatistics is that observations are made directly of the underlying process being studied and that the observations are contaminated with Gaussian errors. Under this assumption, sub–optimal predictions will be obtained if the error characteristics of the sensor are effectively non–Gaussian. One method, model based geostatistics, assumes that a Gaussian process prior is imposed over the (latent) process being studied and that the sensor model forms part of the likelihood term. One problem with this type of approach is that the corresponding posterior distribution will be non–Gaussian and computationally demanding as Monte Carlo methods have to be used. An extension of a sequential, approximate Bayesian inference method enables observations with arbitrary likelihoods to be treated, in a projected process kriging framework which is less computationally intensive. The approach is illustrated using a simulated dataset with a range of sensor models and error characteristics.

Studies of the shellside performance of shell-and-tube heat exchangers

Accurate prediction of shellside pressure drop in a baffled shell-and-tube heat exchanger is very difficult because of the complicated shellside geometry. Ideally, all the shellside fluid should be alternately deflected across the tube bundle as it traverses from inlet to outlet. In practice, up to 60% of the shellside fluid may bypass the tube bundle or leak through the baffles. This short-circuiting of the main flow reduces the efficiency of the exchanger. Of the various shellside methods, it is shown that only the multi-stream methods, which attempt to obtain the shellside flow distribution, predict the pressure drop with any degree of accuracy, the various predictions ranging from -30% to +70%, generally overpredicting. It is shown that the inaccuracies are mainly due to the manner in which baffle leakage is modelled. The present multi-stream methods do not allow for interactions of the various flowstreams, and yet it is shown that three main effects are identified, a) there is a strong interaction between the main cross flow and the baffle leakage streams, enhancing the crossflow pressure drop, b) there is a further short-circuit not considered previously i.e. leakage in the window, and c) the crossflow does not penetrate as far, on average, as previously supposed. Models are developed for each of these three effects, along with a new windowflow pressure drop model, and it is shown that the effect of baffle leakage in the window is the most significant. These models developed to allow for various interactions, lead to an improved multi-stream method, named the "STREAM-INTERACTION" method. The overall method is shown to be consistently more accurate than previous methods, with virtually all the available shellside data being predicted to within ±30% and over 60% being within ±20%. The method is, thus, strongly recommended for use as a design method.

Projected sequential Gaussian processes:a C++ tool for interpolation of large datasets with heterogeneous noise

Heterogeneous datasets arise naturally in most applications due to the use of a variety of sensors and measuring platforms. Such datasets can be heterogeneous in terms of the error characteristics and sensor models. Treating such data is most naturally accomplished using a Bayesian or model-based geostatistical approach; however, such methods generally scale rather badly with the size of dataset, and require computationally expensive Monte Carlo based inference. Recently within the machine learning and spatial statistics communities many papers have explored the potential of reduced rank representations of the covariance matrix, often referred to as projected or fixed rank approaches. In such methods the covariance function of the posterior process is represented by a reduced rank approximation which is chosen such that there is minimal information loss. In this paper a sequential Bayesian framework for inference in such projected processes is presented. The observations are considered one at a time which avoids the need for high dimensional integrals typically required in a Bayesian approach. A C++ library, gptk, which is part of the INTAMAP web service, is introduced which implements projected, sequential estimation and adds several novel features. In particular the library includes the ability to use a generic observation operator, or sensor model, to permit data fusion. It is also possible to cope with a range of observation error characteristics, including non-Gaussian observation errors. Inference for the covariance parameters is explored, including the impact of the projected process approximation on likelihood profiles. We illustrate the projected sequential method in application to synthetic and real datasets. Limitations and extensions are discussed. © 2010 Elsevier Ltd.

Projected Reorganization of Florida Bay Seagrass Communities in Response to the Increased Freshwater Inflow of Everglades Restoration

Historic changes in water-use management in the Florida Everglades have caused the quantity of freshwater inflow to Florida Bay to decline by approximately 60% while altering its timing and spatial distribution. Two consequences have been (1) increased salinity throughout the bay, including occurrences of hypersalinity, coupled with a decrease in salinity variability, and (2) change in benthic habitat structure. Restoration goals have been proposed to return the salinity climates (salinity and its variability) of Florida Bay to more estuarine conditions through changes in upstream water management, thereby returning seagrass species cover to a more historic state. To assess the potential for meeting those goals, we used two modeling approaches and long-term monitoring data. First, we applied the hydrological mass balance model FATHOM to predict salinity climate changes in sub-basins throughout the bay in response to a broad range of freshwater inflow from the Everglades. Second, because seagrass species exhibit different sensitivities to salinity climates, we used the FATHOM-modeled salinity climates as input to a statistical discriminant function model that associates eight seagrass community types with water quality variables including salinity, salinity variability, total organic carbon, total phosphorus, nitrate, and ammonium, as well as sediment depth and light reaching the benthos. Salinity climates in the western sub-basins bordering the Gulf of Mexico were insensitive to even the largest (5-fold) modeled increases in freshwater inflow. However, the north, northeastern, and eastern sub-basins were highly sensitive to freshwater inflow and responded to comparatively small increases with decreased salinity and increased salinity variability. The discriminant function model predicted increased occurrences ofHalodule wrightii communities and decreased occurrences of Thalassia testudinum communities in response to the more estuarine salinity climates. The shift in community composition represents a return to the historically observed state and suggests that restoration goals for Florida Bay can be achieved through restoration of freshwater inflow from the Everglades.

Data compilation on the biological response to ocean acidification: environmental and experimental context of data sets and related literature

The exponential growth of studies on the biological response to ocean acidification over the last few decades has generated a large amount of data. To facilitate data comparison, a data compilation hosted at the data publisher PANGAEA was initiated in 2008 and is updated on a regular basis (doi:10.1594/PANGAEA.149999). By January 2015, a total of 581 data sets (over 4 000 000 data points) from 539 papers had been archived. Here we present the developments of this data compilation five years since its first description by Nisumaa et al. (2010). Most of study sites from which data archived are still in the Northern Hemisphere and the number of archived data from studies from the Southern Hemisphere and polar oceans are still relatively low. Data from 60 studies that investigated the response of a mix of organisms or natural communities were all added after 2010, indicating a welcomed shift from the study of individual organisms to communities and ecosystems. The initial imbalance of considerably more data archived on calcification and primary production than on other processes has improved. There is also a clear tendency towards more data archived from multifactorial studies after 2010. For easier and more effective access to ocean acidification data, the ocean acidification community is strongly encouraged to contribute to the data archiving effort, and help develop standard vocabularies describing the variables and define best practices for archiving ocean acidification data.