Simultaneous synthesis of heat exchanger networks with pressure recovery: optimal integration between heat and work

The optimal integration between heat and work may significantly reduce the energy demand and consequently the process cost. This paper introduces a new mathematical model for the simultaneous synthesis of heat exchanger networks (HENs) in which the pressure levels of the process streams can be adjusted to enhance the heat integration. A superstructure is proposed for the HEN design with pressure recovery, developed via generalized disjunctive programming (GDP) and mixed-integer nonlinear programming (MINLP) formulation. The process conditions (stream temperature and pressure) must be optimized. Furthermore, the approach allows for coupling of the turbines and compressors and selection of the turbines and valves to minimize the total annualized cost, which consists of the operational and capital expenses. The model is tested for its applicability in three case studies, including a cryogenic application. The results indicate that the energy integration reduces the quantity of utilities required, thus decreasing the overall cost.

Retrofit of heat exchanger networks with pressure recovery of process streams at sub-ambient conditions

This paper presents a new mathematical programming model for the retrofit of heat exchanger networks (HENs), wherein the pressure recovery of process streams is conducted to enhance heat integration. Particularly applied to cryogenic processes, HENs retrofit with combined heat and work integration is mainly aimed at reducing the use of expensive cold services. The proposed multi-stage superstructure allows the increment of the existing heat transfer area, as well as the use of new equipment for both heat exchange and pressure manipulation. The pressure recovery of streams is carried out simultaneously with the HEN design, such that the process conditions (streams pressure and temperature) are variables of optimization. The mathematical model is formulated using generalized disjunctive programming (GDP) and is optimized via mixed-integer nonlinear programming (MINLP), through the minimization of the retrofit total annualized cost, considering the turbine and compressor coupling with a helper motor. Three case studies are performed to assess the accuracy of the developed approach, including a real industrial example related to liquefied natural gas (LNG) production. The results show that the pressure recovery of streams is efficient for energy savings and, consequently, for decreasing the HEN retrofit total cost especially in sub-ambient processes.

Anelastic strain recovery and elastic properties of ODP Leg 123 oceanic basalt samples

A knowledge of rock stress is fundamental for improving our understanding of oceanic crustal mechanisms and lithospheric dynamic processes. However, direct measurements of stress in the deep oceans, and in particular stress magnitudes, have proved to be technically difficult. Anelastic strain recovery measurements were conducted on 15 basalt core samples from Sites 765 and 766 during Leg 123. Three sets of experiments were performed: anelastic strain recovery monitoring, dynamic elastic property measurements, and thermal azimuthal anisotropy observations. In addition, a range of other tests and observations were recorded to characterize each of the samples. One common feature of the experimental results and observations is that apparently no consistent orientation trend exists, either between the different measurements on each core sample or between the same sets of measurements on the various core samples. However, some evidence of correspondence between velocity anisotropy and anelastic strain recovery exists, but this is not consistent for all the core samples investigated. Thermal azimuthal anisotropy observations, although showing no conclusive correlations with the other results, were of significant interest in that they clearly exhibited anisotropic behavior. The apparent reproducibility of this behavior may point toward the possibility of rocks that retain a "memory" of their stress history, which could be exploited to derive stress orientations from archived core. Anelastic strain recovery is a relatively new technique. Because use of the method has extended to a wider range of rock types, the literature has begun to include examples of rocks that contracted with time. Strong circumstantial evidence exists to suggest that core-sample contractions result from the slow diffusion of pore fluids from a preexisting microcrack structure that permits the rock to deflate at a greater rate than the expansion caused by anelastic strain recovery. Both expansions and contractions of the Leg 123 cores were observed. The basalt cores have clearly been intersected by an abundance of preexisting fractures, some of which pass right through the samples, but many are intercepted or terminate within the rock matrix. Thus, the behavior of the core samples will be influenced not only by the properties of the rock matrix between the fractures, but also by how these macro- and micro-scale fractures mutually interact. The strain-recovery curves recorded during Leg 123 for each of the 15 basalt core samples may reflect the result of two competing time dependent processes: anelastic strain recovery and pore pressure recovery. Were these the only two processes to influence the gauge responses, then one might expect that given the additional information required, established theoretical models might be used to determine consistent stress orientations and reliable stress magnitudes. However, superimposed upon these competing processes is their respective interaction with the preexisting fractures that intersect each core. Evidence from our experiments and observations suggests that these fractures have a dominating influence on the characteristics of the recovery curves and that their effects are complex.

Säätöventtiilin virtaustekniset ominaisuudet

Virtauslaskennan käyttö jokapäiväisessä insinöörityössä on lisääntynyt viime vuosina nopeaa vauhtia. Virtauslaskennan avulla voidaan tutkia säätöventtiilin virtauskenttää, mikä antaa suunnittelijalle mahdollisuuden korjata virtauskanavan ongelmakohtia jo tuotekehityksen alkuvaiheessa. Tämändiplomityön tavoitteena on määrittää uuden säätöventtiilin mitoituskertoimet jatutkia virtauslaskennan käytettävyyttä säätöventtiilisuunnittelussa. Teoreettisessa tarkastelussa on käsitelty venttiilivirtaukselle ominaisia virtausteknisiä yhtälöitä ja ilmiöitä, säätöventtiilin standardin määräämiä mitoitusyhtälöitä sekä neste- että kaasumelua. Lisäksi kerrotaan yleisimmistä säätöventtiilisovellutuksista ja esitellään suunnitteilla oleva uusi säätöventtiili. Virtauslaskennan avulla tutkittiin venttiilin kapasiteettiaja virtauskenttää. Alustavaa laskentaa tehtiin venttiilin paineenpalautumiskertoimen ja alkavan kavitaation määrittämiseksi. Virtauslaskenta tehtiin Fluent ja Cfdesign -virtauslaskentaohjelmilla. Virtauslaskennan antamia tuloksia verrattiin laboratoriossa saatuihin mittaustuloksiin. Laboratoriokokeiden avulla määritettiin uuden säätöventtiilin mitoituskertoimet. Lisäksi mitattiin säätöventtiilin aiheuttamaa neste- ja kaasumelua.

Computational Analysis and Optimization of Real Gas Flow In Small Centrifugal Compressors

Small centrifugal compressors are more and more widely used in many industrialsystems because of their higher efficiency and better off-design performance comparing to piston and scroll compressors as while as higher work coefficient perstage than in axial compressors. Higher efficiency is always the aim of the designer of compressors. In the present work, the influence of four partsof a small centrifugal compressor that compresses heavy molecular weight real gas has been investigated in order to achieve higher efficiency. Two parts concern the impeller: tip clearance and the circumferential position of the splitter blade. The other two parts concern the diffuser: the pinch shape and vane shape. Computational fluid dynamics is applied in this study. The Reynolds averaged Navier-Stokes flow solver Finflo is used. The quasi-steady approach is utilized. Chien's k-e turbulence model is used to model the turbulence. A new practical real gas model is presented in this study. The real gas model is easily generated, accuracy controllable and fairly fast. The numerical results and measurements show good agreement. The influence of tip clearance on the performance of a small compressor is obvious. The pressure ratio and efficiency are decreased as the size of tip clearance is increased, while the total enthalpy rise keeps almost constant. The decrement of the pressure ratio and efficiency is larger at higher mass flow rates and smaller at lower mass flow rates. The flow angles at the inlet and outlet of the impeller are increased as the size of tip clearance is increased. The results of the detailed flow field show that leakingflow is the main reason for the performance drop. The secondary flow region becomes larger as the size of tip clearance is increased and the area of the main flow is compressed. The flow uniformity is then decreased. A detailed study shows that the leaking flow rate is higher near the exit of the impeller than that near the inlet of the impeller. Based on this phenomenon, a new partiallyshrouded impeller is used. The impeller is shrouded near the exit of the impeller. The results show that the flow field near the exit of the impeller is greatly changed by the partially shrouded impeller, and better performance is achievedthan with the unshrouded impeller. The loading distribution on the impeller blade and the flow fields in the impeller is changed by moving the splitter of the impeller in circumferential direction. Moving the splitter slightly to the suction side of the long blade can improve the performance of the compressor. The total enthalpy rise is reduced if only the leading edge of the splitter ismoved to the suction side of the long blade. The performance of the compressor is decreased if the blade is bended from the radius direction at the leading edge of the splitter. The total pressure rise and the enthalpy rise of thecompressor are increased if pinch is used at the diffuser inlet. Among the fivedifferent pinch shape configurations, at design and lower mass flow rates the efficiency of a straight line pinch is the highest, while at higher mass flow rate, the efficiency of a concave pinch is the highest. The sharp corner of the pinch is the main reason for the decrease of efficiency and should be avoided. The variation of the flow angles entering the diffuser in spanwise direction is decreased if pinch is applied. A three-dimensional low solidity twisted vaned diffuser is designed to match the flow angles entering the diffuser. The numerical results show that the pressure recovery in the twisted diffuser is higher than in a conventional low solidity vaned diffuser, which also leads to higher efficiency of the twisted diffuser. Investigation of the detailed flow fields shows that the separation at lower mass flow rate in the twisted diffuser is later than in the conventional low solidity vaned diffuser, which leads to a possible wider flow range of the twisted diffuser.

Simultaneous synthesis of work exchange networks with heat integration

The optimal integration of work and its interaction with heat can represent large energy savings in industrial plants. This paper introduces a new optimization model for the simultaneous synthesis of work exchange networks (WENs), with heat integration for the optimal pressure recovery of process gaseous streams. The proposed approach for the WEN synthesis is analogous to the well-known problem of synthesis of heat exchanger networks (HENs). Thus, there is work exchange between high-pressure (HP) and low-pressure (LP) streams, achieved by pressure manipulation equipment running on common axes. The model allows the use of several units of single-shaft-turbine-compressor (SSTC), as well as stand-alone compressors, turbines and valves. Helper motors and generators are used to respond to any demand and excess of energy. Moreover, between the WEN stages the streams are sent to the HEN to promote thermal recovery, aiming to enhance the work integration. A multi-stage superstructure is proposed to represent the process. The WEN superstructure is optimized in a mixed-integer nonlinear programming (MINLP) formulation and solved with the GAMS software, with the goal of minimizing the total annualized cost. Three examples are conducted to verify the accuracy of the proposed method. In all case studies, the heat integration between WEN stages is essential to improve the pressure recovery, and to reduce the total costs involved in the process.

Regional subsidence modelling in Murcia city (SE Spain) using 1-D vertical finite element analysis and 2-D interpolation of ground surface displacements

Subsidence is a hazard that may have natural or anthropogenic origin causing important economic losses. The area of Murcia city (SE Spain) has been affected by subsidence due to groundwater overexploitation since the year 1992. The main observed historical piezometric level declines occurred in the periods 1982–1984, 1992–1995 and 2004–2008 and showed a close correlation with the temporal evolution of ground displacements. Since 2008, the pressure recovery in the aquifer has led to an uplift of the ground surface that has been detected by the extensometers. In the present work an elastic hydro-mechanical finite element code has been used to compute the subsidence time series for 24 geotechnical boreholes, prescribing the measured groundwater table evolution. The achieved results have been compared with the displacements estimated through an advanced DInSAR technique and measured by the extensometers. These spatio-temporal comparisons have showed that, in spite of the limited geomechanical data available, the model has turned out to satisfactorily reproduce the subsidence phenomenon affecting Murcia City. The model will allow the prediction of future induced deformations and the consequences of any piezometric level variation in the study area.

Gas hydrate investigation at the Amsterdam mud volcano

We investigated gas hydrate in situ inventories as well as the composition and principal transport mechanisms of fluids expelled at the Amsterdam mud volcano (AMV; 2,025 m water depth) in the Eastern Mediterranean Sea. Pressure coring (the only technique preventing hydrates from decomposition during recovery) was used for the quantification of light hydrocarbons in near-surface deposits. The cores (up to 2.5 m in length) were retrieved with an autoclave piston corer, and served for analyses of gas quantities and compositions, and pore-water chemistry. For comparison, gravity cores from sites at the summit and beyond the AMV were analyzed. A prevalence of thermogenic light hydrocarbons was inferred from average C1/C2+ ratios <35 and d13C-CH4 values of -50.6 per mil. Gas venting from the seafloor indicated methane oversaturation, and volumetric gas-sediment ratios of up to 17.0 in pressure cores taken from the center demonstrated hydrate presence at the time of sampling. Relative enrichments in ethane, propane, and iso-butane in gas released from pressure cores, and from an intact hydrate piece compared to venting gas suggest incipient crystallization of hydrate structure II (sII). Nonetheless, the co-existence of sI hydrate can not be excluded from our dataset. Hydrates fill up to 16.7% of pore volume within the sediment interval between the base of the sulfate zone and the maximum sampling depth at the summit. The concave-down shapes of pore-water concentration profiles recorded in the center indicate the influence of upward-directed advection of low-salinity fluids/fluidized mud. Furthermore, the SO42- and Ba2+ pore-water profiles in the central part of the AMV demonstrate that sulfate reduction driven by the anaerobic oxidation of methane is complete at depths between 30 cm and 70 cm below seafloor. Our results indicate that methane oversaturation, high hydrostatic pressure, and elevated pore-water activity caused by low salinity promote fixing of considerable proportions of light hydrocarbons in shallow hydrates even at the summit of the AMV, and possibly also of other MVs in the region. Depending on their crystallographic structure, however, hydrates will already decompose and release hydrocarbon masses if sediment temperatures exceed ca. 19.3°C and 21.0°C, respectively. Based on observations from other mud volcanoes, the common occurrence of such temperatures induced by heat flux from below into the immediate subsurface appears likely for the AMV.

Large eddy simulation of shock boundary layer interaction control using micro-vortex generators

The performance of supersonic engine inlets and external aerodynamic surfaces can be critically affected by shock wave / boundary layer interactions (SBLIs), whose severe adverse pressure gradients can cause boundary layer separation. Currently such problems are avoided primarily through the use of boundary layer bleed/suction which can be a source of significant performance degradation. This study investigates a novel type of flow control device called micro-vortex generators (µVGs) which may offer similar control benefits without the bleed penalties. µVGs have the ability to alter the near-wall structure of compressible turbulent boundary layers to provide increased mixing of high speed fluid which improves the boundary layer health when subjected to flow disturbance. Due to their small size,µVGs are embedded in the boundary layer which provide reduced drag compared to the traditional vortex generators while they are cost-effective, physically robust and do not require a power source. To examine the potential of µVGs, a detailed experimental and computational study of micro-ramps in a supersonic boundary layer at Mach 3 subjected to an oblique shock was undertaken. The experiments employed a flat plate boundary layer with an impinging oblique shock with downstream total pressure measurements. The moderate Reynolds number of 3,800 based on displacement thickness allowed the computations to use Large Eddy Simulations without the subgrid stress model (LES-nSGS). The LES predictions indicated that the shock changes the structure of the turbulent eddies and the primary vortices generated from the micro-ramp. Furthermore, they generally reproduced the experimentally obtained mean velocity profiles, unlike similarly-resolved RANS computations. The experiments and the LES results indicate that the micro-ramps, whose height is h≈0.5δ, can significantly reduce boundary layer thickness and improve downstream boundary layer health as measured by the incompressible shape factor, H. Regions directly behind the ramp centerline tended to have increased boundary layer thickness indicating the significant three-dimensionality of the flow field. Compared to baseline sizes, smaller micro-ramps yielded improved total pressure recovery. Moving the smaller ramps closer to the shock interaction also reduced the displacement thickness and the separated area. This effect is attributed to decreased wave drag and the closer proximity of the vortex pairs to the wall. In the second part of the study, various types of µVGs are investigated including micro-ramps and micro-vanes. The results showed that vortices generated from µVGs can partially eliminate shock induced flow separation and can continue to entrain high momentum flux for boundary layer recovery downstream. The micro-ramps resulted in thinner downstream displacement thickness in comparison to the micro-vanes. However, the strength of the streamwise vorticity for the micro-ramps decayed faster due to dissipation especially after the shock interaction. In addition, the close spanwise distance between each vortex for the ramp geometry causes the vortex cores to move upwards from the wall due to induced upwash effects. Micro-vanes, on the other hand, yielded an increased spanwise spacing of the streamwise vortices at the point of formation. This resulted in streamwise vortices staying closer to the wall with less circulation decay, and the reduction in overall flow separation is attributed to these effects. Two hybrid concepts, named “thick-vane” and “split-ramp”, were also studied where the former is a vane with side supports and the latter has a uniform spacing along the centerline of the baseline ramp. These geometries behaved similar to the micro-vanes in terms of the streamwise vorticity and the ability to reduce flow separation, but are more physically robust than the thin vanes. Next, Mach number effect on flow past the micro-ramps (h~0.5δ) are examined in a supersonic boundary layer at M=1.4, 2.2 and 3.0, but with no shock waves present. The LES results indicate that micro-ramps have a greater impact at lower Mach number near the device but its influence decays faster than that for the higher Mach number cases. This may be due to the additional dissipation caused by the primary vortices with smaller effective diameter at the lower Mach number such that their coherency is easily lost causing the streamwise vorticity and the turbulent kinetic energy to decay quickly. The normal distance between the vortex core and the wall had similar growth indicating weak correlation with the Mach number; however, the spanwise distance between the two counter-rotating cores further increases with lower Mach number. Finally, various µVGs which include micro-ramp, split-ramp and a new hybrid concept “ramped-vane” are investigated under normal shock conditions at Mach number of 1.3. In particular, the ramped-vane was studied extensively by varying its size, interior spacing of the device and streamwise position respect to the shock. The ramped-vane provided increased vorticity compared to the micro-ramp and the split-ramp. This significantly reduced the separation length downstream of the device centerline where a larger ramped-vane with increased trailing edge gap yielded a fully attached flow at the centerline of separation region. The results from coarse-resolution LES studies show that the larger ramped-vane provided the most reductions in the turbulent kinetic energy and pressure fluctuation compared to other devices downstream of the shock. Additional benefits include negligible drag while the reductions in displacement thickness and shape factor were seen compared to other devices. Increased wall shear stress and pressure recovery were found with the larger ramped-vane in the baseline resolution LES studies which also gave decreased amplitudes of the pressure fluctuations downstream of the shock.

Enhanced diastolic reflections on arterial pressure pulse during exercise recovery.

During recovery from a maximal or submaximal aerobic exercise, augmentation of central (aortic) systolic pressure by reflected pressure waves is blunted in healthy humans. However, the extent to which reflected pressure waves modify the central pulse in diastole in these conditions remains unknown. We evaluated systolic and diastolic central reflected waves in 11 endurance-trained athletes on recovery from a maximal running test on a treadmill (treadmill-max) and a 4000 m run in field conditions. On both occasions in each subject, the radial pulse was recorded with applanation tonometry in the resting preexercise state and then 5, 15, 25, 35, and 45 min after exercise termination. From the central waveform, as reconstructed by application of a generalized transfer function, we computed a systolic (AIx) and a diastolic index (AId) of pressure augmentation by reflections. At 5 min, both indices were below preexercise. At further time-points, AIx remained low, while AId progressively increased, to overshoot above preexercise at 45 min. The same behavior was observed with both exercise types. Beyond the first few minutes of recovery following either maximal or submaximal aerobic exercise, reflected waves selectively augment the central pressure pulse in diastole, at least in endurance-trained athletes.

Enhanced diastolic reflections on arterial pressure pulse during exercise recovery

RESUMEDurant la phase de récupération d'un exercice de course à pied d'intensité maximale ou submaximale, une augmentation de la pression artérielle systolique centrale (aortique) résultant de la réflexion des ondes de pouls sur l'arbre vasculaire est constatée chez l'individu en bonne santé. En diastole cependant, l'impact de la réflexion de ces ondes de pouls sur la pression centrale demeure inconnu durant la récupération d'un exercice.Nous avons évalué les ondes de pouls centrales systolique et diastolique chez onze athlètes d'endurance durant la phase de récupération d'un exercice de course à pied dans des conditions d'effort maximal (sur tapis de course) et lors d'un effort submaximal lors d'une course à pied de 4000 mètres en plein air sur terrain mixte.Pour chaque sujet et lors des deux exercices, l'onde de pouls a été mesurée au niveau radial par tonométrie d'aplanation durant une phase de repos précédant l'exercice, puis à 5, 15, 25, 35 et 45 minutes après la fin de l'exercice. En utilisant une fonction mathématique de transfert, l'onde de pouls centrale a été extrapolée à partir de l'onde de pouls radiale. En compilant la forme de l'onde de pouls centrale avec une mesure simultanée de la pression artérielle brachiale, un index d'augmentation de l'onde de pouls en systole (Alx) et en diastole (Als) peut être calculé, reflétant l'augmentation des pressions résultant de la réflexion des ondes sur l'arbre vasculaire périphérique.A 5 minutes de la fin de l'exercice, les deux index ont été mesurés moindres que ceux mesurés lors de la phase précédant celui-ci. Lors des mesures suivantes, Alx est resté bas, alors que Aid a progressivement augmenté pour finalement dépasser la valeur de repos après 45 minutes de récupération. Le même phénomène a été constaté pour les deux modalités d'exercice (maximal ou submaximal). Ainsi, au-delà de quelques minutes de récupération après un exercice de course d'intensité maximale ou submaximale, nous avons montré par ces investigations que les ondes de pouls réfléchies en périphérie augmentent de façon sélective la pression centrale en diastole chez l'athlète d'endurance.ABSTRACTDuring recovery from a maximal or submaximal aerobic exercise, augmentation of central (aortic) systolic pressure by reflected pressure waves is blunted in healthy humans. However, the extent to which reflected pressure waves modify the central pulse in diastole in these conditions remains unknown. We evaluated systolic and diastolic central reflected waves in 11 endurance-trained athletes on recovery from a maximal running test on a treadmill (treadmill-max) and a 4000m run in field conditions. On both occasions in each subject, the radial pulse was recorded with applanation tonometry in the resting preexercise state and then 5, 15, 25, 35, and 45 minutes after exercise termination. From the central waveform, as reconstructed by application of a generalized transfer function, we computed a systolic (Alx) and a diastolic index (Aid) of pressure augmentation by reflections. At 5 minutes, both indices were below preexercise. At further time-points, Alx remained low, while Aid progressively increased, to overshoot above preexercise at 45 minutes. The same behavior was observed with both exercise types. Beyond the first few minutes of recovery following either maximal or submaximal aerobic exercise, reflected waves selectively augment the central pressure pulse in diastole, at least in endurance- trained athletes.

Enhancing the recovery of tiger nut (Cyperus esculentus) oil by mechanical pressing: moisture content, particle size, high pressure and enzymatic pre-treatment effects

Tiger nut (Cyperus esculentus) tuber contains oil that is high in monounsaturated fatty acids, and this oil makes up about 23% of the tuber. The study aimed at evaluating the impact of several factors and enzymatic pre-treatment on the recovery of pressed tiger nut oil. Smaller particles were more favourable for pressing. High pressure pre-treatment did not increase oil recovery but enzymatic treatment did. The highest yield obtained by enzymatic treatment prior to mechanical extraction was 33 % on a dry defatted basis, which represents a recovery of 90 % of the oil. Tiger nut oil consists mainly of oleic acid; its acid and peroxide values reflect the high stability of the oil.

Recovery of high blood pressure after chronic lesions of the commissural NTS in SHR

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)