Realtime Optimization of MPC Setpoints using Time-Varying Extremum Seeking Control for Vapor Compression Machines

Recently, model predictive control (MPC) has received increased attention in the HVAC community, largely due to its ability to systematically manage constraints while optimally regulating signals of interest to setpoints. For example, in a common formulation of an MPC control problem for variable compressor speed vapor compression machines, the setpoints often include the zone temperature and the evaporator superheat temperature. However, the energy consumption of vapor compression systems has been shown to be sensitive to these setpoints. Further, while superheat temperature is often preferred because it can be easily correlated to heat exchanger efficiency (and therefore cycle efficiency), direct measurement of superheat is not always available. Therefore, identifying alternate signals in the control of vapor compression machines that correlate to efficiency is desired. Conventionally, methods for maximizing the energy efficiency rely on the use of mathematical models of the physics of vapor compression systems. These model-based approaches attempt to describe the influence of commanded inputs on the thermodynamic behavior of the system and the consumed electrical energy, and they are used to predict the combination of inputs that both meet the heat load requirements and minimize energy consumption. However, these models of vapor compression systems rely on simplifying assumptions in order to produce a mathematically tractable representation. Further, they are difficult to derive and calibrate, and often do not describe variations over long time scales, such as those due to refrigerant losses or accumulation of debris on the heat exchangers. In this paper, we consider a model-free extremum seeking algorithm that adjusts setpoints provided to a model predictive controller. While perturbation-based extremum seeking methods have been known for some time, they suffer from slow convergence rates---a problem emphasized by the long time constants associated with thermal systems. Our method uses a new algorithm (time-varying extremum seeking), which has dramatically faster and more reliable convergence properties. In particular, we regulate the compressor discharge temperature using an MPC controller with setpoints selected from a model-free time-varying extremum seeking algorithm. We show that the relationship between compressor discharge temperature and power consumption is convex (a requirement for this class of realtime optimization), and use time-varying extremum seeking to drive these setpoints to values that minimize power. The results are compared to the traditional perturbation-based extremum seeking approach. Further, because the required cooling capacity (and therefore compressor speed) is a function of measured and unmeasured disturbances, the optimal compressor discharge temperature setpoint must vary according to these conditions. We show that the energy optimal discharge temperature is tracked with the time-varying extremum seeking algorithm in the presence of disturbances

Implementation of a Pulsed-Laser Measurement System in the National Transonic Facility

A remotely-adjustable laser transmission and imaging system has been developed for use in a high-pressure, cryogenic wind tunnel. Implementation in the National Transonic Facility has proven the system suitable for velocity and signal lifetime measurements over a range of operating conditions. The measurement system allows for the delivery of high-powered laser pulses through the outer pressure shell and into the test section interior from a mezzanine where the laser is free from environmental disturbances (such as vibrations and excessive condensation) associated with operation of the wind tunnel. Femtosecond laser electronic excitation tagging (FLEET) was utilized to provide freestream velocity measurements, and first results show typical data that may be obtained using the system herein described

Formation Stress Estimation Using Standard Acoustic Logging

In situ formation stress directions and magnitudes are estimated by inverting the borehole flexural and Stoneley dispersions obtained from standard acoustic logging data (dipole and monopole logs). The underlying procedure consists of the following steps: first, we locate stressed zones in the formation by searching for crossovers in flexural dispersions. Second, the fast shear direction is estimated from the cross-dipole waveforms. It corresponds to the direction of the maximum horizontal stress (S[subscript H]). Finally, a multi-frequency inversion of both the Stoneley and flexural dispersions yields the maximum (S[subscript H]) and minimum (S[subscript h]) horizontal stress magnitudes together with the three formation nonlinear elastic constants, C[subscript 111], C[subscript 112] and C[subscript 113], defined about the selected reference (isotropic) state. The inversion method is based on equations that relate S[subscript H] and S[subscript h] with variations in phase velocities of the borehole flexural and Stoneley waves in the stressed state from those in the assumed reference state, the state that is hydrostatically loaded and isotropic. Phase velocities of the borehole flexural and Stoneley modes as a function of frequency can be obtained from processing the cross-dipole and monopole waveforms, respectively. The borehole flexural and Stoneley dispersions in the assumed reference (isotropic) state are obtained from the solution of a standard boundary-value problem. The sensitivity functions for the inversion model are obtained from the eigenfunctions of the boundary-value problem in the reference state. Results for the stress directions and magnitudes obtained from the inversion of the Stoneley and flexural dispersions over a selected bandwidth are consistent with focal mechanism and borehole breakout data present in the world map database (Zoback, 1992).Massachusetts Institute of Technology. Borehole Acoustics and Logging ConsortiumMassachusetts Institute of Technology. Earth Resources Laboratory. Reservoir Delineation Consortiu

Analysis of Scattered Signal to Estimate Reservoir Fracture Parameters

We detect fracture corridors and determine their orientation and average spacing based on an analysis of seismic coda in the frequency-wave number (f-k ) domain. Fracture corridors have dimensions similar to seismic wavelengths which causes scattering. The distribution of energy in shot records in the f-k domain depends upon the orientation of the records relative to the fracture strike. In the direction normal to fractures, scattered waves propagate with slower apparent velocities than waves propagating along the fracture channels. The associated f-k spectral differences allow the identification of the preferred fracture orientation and spacing. We apply our technique to a fractured reservoir in the Lynx field, in the Canadian foothills. The estimated preferential fracture orientation is about N40 E, which agrees with regional stress measurements. The average fracture spacing is 75 m on the West side of the survey, while fractures are more sparse on the East side. We also apply the Scattering Index methodology (Willis et al., 2006) to the same data, post-stack and pre-stack. This technique has higher resolution to map fracture distribution, intensity and orientation, and therefore complements the spectral method in providing an integrated description of reservoir fractures.United States. Dept. of Energy (award number DE-FC26-06NT42956)Massachusetts Institute of Technology. Earth Resources Laborator

A comparison of LWD and wireline dipole sonic data

Data measured by both wireline and LWD tools in the same borehole are compared. Discrepancies in shear velocities as calculated from the data are on average around 5% and discrepancies between compressional velocities are less than 3%. The consistency of the bias between logs suggest it is related to the calculation of velocity. Comparison of industry and ERL velocity processing show excellent agreement and give an example of possible spread of velocity data due to processing chain. A short section of data in an unconsolidated zone shows velocity differences of just over 10% with an opposite trend to the over all bias. Dispersion analysis of the waveforms show this is consistent with a damaged zone surrounding the borehole wall caused by drilling.Massachusetts Institute of Technology. Earth Resources LaboratoryMassachusetts Institute of Technology. Borehole Acoustics and Logging Consortiu

The Role of Entrepreneur-Venture Fit in Online Home-based Entrepreneurship: A Systematic Literature Review

Home-based businesses and their founders represent an important, but under-researched facet of entrepreneurship. Far from being small, hobby-businesses with little economic impact, home-based business make significant contribution to national economies in terms of both turnover and employment. Online home-based businesses have been recognised as an important and distinct sector of the home-based business domain, offering unique opportunity for innovation and business diversity. The paper presents a systematic literature review of extant research on online home-based entrepreneurs and their businesses. The findings of the review are structured and discussed using the theoretical lens of entrepreneur-venture fit. Use of this lens allows the study to bring coherence to previously fragmented extant studies, providing a basis for future research in this domain. The study also develops a novel model of entrepreneur-venture fit in the specific case of online home-based businesses. This allows us to suggest five positive interactions between entrepreneurial and venture characteristics. It also allows us to suggest a number of previously unidentified negative interactions, which may result in entrepreneurs becoming ‘locked-in’ and suffering multiple sources of stress

Joint flow-seismic inversion for characterizing fractured reservoirs: theoretical approach and numerical modeling

Traditionally, seismic interpretation is performed without any account of the flow behavior. Here, we present a methodology to characterize fractured geologic media by integrating flow and seismic data. The key element of the proposed approach is the identification of the intimate relation between acoustic and flow responses of a fractured reservoir through the fracture compliance. By means of synthetic models, we show that: (1) owing to the strong (but highly uncertain) dependence of fracture permeability on fracture compliance, the modeled flow response in a fractured reservoir is highly sensitive to the geophysical interpretation; and (2) by incorporating flow data (well pressures and production curves) into the inversion workflow, we can simultaneously reduce the error in the seismic interpretation and improve predictions of the reservoir flow dynamics.Eni-MIT Energy Initiative Founding Member Progra

Supersymmetric particle mass measurement with the boost-corrected contransverse mass

A modification to the contransverse mass (MCT) technique for measuring the masses of pair-produced semi-invisibly decaying heavy particles is proposed in which MCT is corrected for non-zero boosts of the centre-of-momentum (CoM) frame of the heavy states in the laboratory transverse plane. Lack of knowledge of the mass of the CoM frame prevents exact correction for this boost, however it is shown that a conservative correction can nevertheless be derived which always generates an MCT value which is less than or equal to the true value of MCT in the CoM frame. The new technique is demonstrated with case studies of mass measurement with fully leptonic ttbar events and with SUSY events possessing a similar final state.Comment: 33 pages, 33 .eps figures, JHEP3 styl

Dimension- and Context-Specific Expression of Preschoolers’ Disruptive Behaviors Associated with Prenatal Tobacco Exposure

Objective—Precise phenotypic characterization of prenatal tobacco exposure (PTE) −related disruptive behavior (DB) that integrates nuanced measures of both exposures and outcomes is optimal for elucidating underlying mechanisms. Using this approach, our goals were to identify dimensions of DB most sensitive to PTE prior to school entry and assess contextual variation in these dimensions. Methods—A community obstetric sample of N=369 women (79.2% lifetime smokers; 70.2% pregnancy smokers) from two Midwestern cities were assessed for PTE using cotinine-calibrated interview-based reports at 16, 28, and 40 weeks of gestation. A subset of n=244 who completed observational assessments with their 5-year-old children in a subsequent preschool follow-up study constitute the analytic sample. Using two developmentally-meaningful dimensions previously associated with emergent clinical risk for DB—irritability and noncompliance—we assessed children with 2 parent-report scales: the Multidimensional Assessment Profile of Disruptive Behavior (MAP-DB) and the Early Childhood Inventory (ECI). We also assessed children by direct observation across 3 interactional contexts with the Disruptive Behavior Diagnostic Observation Schedule (DB-DOS). We used generalized linear models to examine between-child variability across behavioral dimensions, and mixed effects models to examine directly observed within-child variability by interactional context. Results—Increasing PTE predicted increasing impairment in preschoolers’ modulation of negative affect (irritability), but not negative behavior (noncompliance) across reported (MAP-DB) and observed (DB-DOS) dimensional measures. Moreover, children’s PTE-related irritability was more pronounced when observed with parents than with the examiner. The ECI did not detect PTE-related irritability nor noncompliance. Conclusions—Nuanced, dimension- and context-specific characterization of PTE-related DB described can optimize early identification of at-risk children

Managing understory light conditions in boreal mixedwoods through variation in the intensity and spatial pattern of harvest: A modelling approach

In the context of partial harvesting, adequately managing post-harvest light conditions are essential to obtain a desired composition of tree species regeneration. The objective of this study was to determine how varying the intensity and spatial pattern of harvest would affect understory light conditions in boreal mixedwood stands of northwestern Quebec using the spatially explicit SORTIE-ND light model. The model was evaluated based on comparisons of observed and predicted light levels in both mapped and un-mapped plots. In mapped plots, reasonably accurate predictions of the overall variation in light levels were obtained, but predictions tended to lack spatial precision. In un-mapped plots, SORTIE-ND accurately predicted stand-level mean GLI (Gap Light Index) under a range of harvest intensities. The model was then used to simulate nine silvicultural treatments based on combinations of three intensities of overstory removal (30%, 45% and 60% of basal area) and three harvest patterns (uniform, narrow strips, large gaps). Simulations showed that increasing overstory removal had less impact on light conditions with uniform harvests, and a more marked effect with more aggregated harvest patterns. Whatever the harvest intensity, uniform cuts almost never created high light conditions (GLI > 50%). Gap cuts, on the other hand, resulted in up to 40% of microsites receiving GLI > 50%. Our results suggest that either a 30% strip or gap cut or a 45–60% uniform partial harvest could be used to accelerate the transition from an aspen dominated composition to a mixedwood stand because both types of cut generate the greatest proportion of moderately low light levels (e.g., 15–40% GLI). These light levels tend to favour an accelerated growth response among shade-tolerant conifers, while preventing excessive recruitment of shade-intolerant species. A better understanding of how spatial patterns of harvest interact with tree removal intensity to affect understory light conditions can provide opportunities for designing silvicultural prescriptions that are tailored to species’ traits and better suited to meet a variety of management objectives