A randomized concave programming method for choice network revenue management

Models incorporating more realistic models of customer behavior, as customers choosing froman offer set, have recently become popular in assortment optimization and revenue management.The dynamic program for these models is intractable and approximated by a deterministiclinear program called the CDLP which has an exponential number of columns. However, whenthe segment consideration sets overlap, the CDLP is difficult to solve. Column generationhas been proposed but finding an entering column has been shown to be NP-hard. In thispaper we propose a new approach called SDCP to solving CDLP based on segments and theirconsideration sets. SDCP is a relaxation of CDLP and hence forms a looser upper bound onthe dynamic program but coincides with CDLP for the case of non-overlapping segments. Ifthe number of elements in a consideration set for a segment is not very large (SDCP) can beapplied to any discrete-choice model of consumer behavior. We tighten the SDCP bound by(i) simulations, called the randomized concave programming (RCP) method, and (ii) by addingcuts to a recent compact formulation of the problem for a latent multinomial-choice model ofdemand (SBLP+). This latter approach turns out to be very effective, essentially obtainingCDLP value, and excellent revenue performance in simulations, even for overlapping segments.By formulating the problem as a separation problem, we give insight into why CDLP is easyfor the MNL with non-overlapping considerations sets and why generalizations of MNL posedifficulties. We perform numerical simulations to determine the revenue performance of all themethods on reference data sets in the literature.

Parallel scheduling of multiclass M/M/m queues: Approximate and heavy-traffic optimization of achievable performance

We address the problem of scheduling a multiclass $M/M/m$ queue with Bernoulli feedback on $m$ parallel servers to minimize time-average linear holding costs. We analyze the performance of a heuristic priority-index rule, which extends Klimov's optimal solution to the single-server case: servers select preemptively customers with larger Klimov indices. We present closed-form suboptimality bounds (approximate optimality) for Klimov's rule, which imply that its suboptimality gap is uniformly bounded above with respect to (i) external arrival rates, as long as they stay within system capacity;and (ii) the number of servers. It follows that its relativesuboptimality gap vanishes in a heavy-traffic limit, as external arrival rates approach system capacity (heavy-traffic optimality). We obtain simpler expressions for the special no-feedback case, where the heuristic reduces to the classical $c \mu$ rule. Our analysis is based on comparing the expected cost of Klimov's ruleto the value of a strong linear programming (LP) relaxation of the system's region of achievable performance of mean queue lengths. In order to obtain this relaxation, we derive and exploit a new set ofwork decomposition laws for the parallel-server system. We further report on the results of a computational study on the quality of the $c \mu$ rule for parallel scheduling.

Improving electricity market price scenarios by means of forecasting factor models

In liberalized electricity markets, generation Companies must build an hourly bidthat is sent to the market operator. The price at which the energy will be paid is unknown during the bidding process and has to be forecast. In this work we apply forecasting factor models to this framework and study its suitability.

Ispol'zovanie metodov teorii vosstanovleniia v modeliakh optimal'nogo dobyvaniia pishchi [Restoration theory in models of optimal feeding behavior]

The method of stochastic dynamic programming is widely used in ecology of behavior, but has some imperfections because of use of temporal limits. The authors presented an alternative approach based on the methods of the theory of restoration. Suggested method uses cumulative energy reserves per time unit as a criterium, that leads to stationary cycles in the area of states. This approach allows to study the optimal feeding by analytic methods.

Role of the inflammasome in murine experimental models of arthritis

Summary : The purpose of this study was to investigate the role of the inflammasome in human and experimental murine models (such as ΑΙΑ and K/BxN) of rheumatoid arthritis (RA)RA, affecting 1% of the population is the most frequent inflammatory disease characterized by synovial hyperplasia and cartilage and bone erosion, leading to joint destruction. In general, women are 3 times more affected by RA suggesting a role of estrogen in this disease. The inflammasome is a multiproteic complex triggering the activation of caspase-1 leading to the activation of IL-1 β, an important pro-inflammatory cytokine implicated in arthritis. The inflammasome has been implicated in several inflammatory diseases and particularly in gout. To highlight a possible role of the inflammasome in murine arthritis, we obtained ASC, caspase-1 and NALP3 +/+ and -/- littermate mice to perform ΑΙΑ and K/BxN arthritis. NALP3 -/- and caspase-1 -/- mice were as arthritic as wild type littermate mice in both ΑΙΑ and K/BxN models implicating that the NALP3 inflammasome is not involved in experimental arthritis. By contrast, ΑΙΑ severity was significantly diminished in ASC- deficient male and female mice, and in the K/BxN model, in ASC-deficient female mice. These results were supported by histological scoring and acute phase protein serum amyloid A (SAA) levels that were equivalent between NALP+/+ and NALP3-/- mice and diminished in ASC -/- mice. In ΑΙΑ and K/BxN murine experimental models, we observed a sexdependent phenotype. We studied the role of estradiol in both the ALA and the K/BxN models. Castrated female or male ASC -/- mice that received estradiol had a decreased arthritis severity. This implies a protective role of estrogen in the absence of ASC. In the ΑΙΑ model, proliferation assay were performed using splenocytes from mBSA- immunized ASC +/+ and -/- mice. The mBSA-induced proliferation was significantly lower in ASC-/- splenocytes. Moreover the CD3-specific proliferation of purified splenic Τ cells was significantly lower in ASC-/- cells. Finally, Τ cells from ASC-/- mice produced significantly decreased levels of IFN-gamma associated with increased levels of IL-10. These results imply a possible role of ASC in the TCR-signaling pathway and Τ cell cytokine production. In parallel the expression of the different inflammasome components were analyzed in biopsies from rheumatoid arthritis (RA) and osteoarthritis (OA) patiens. The expression of the 14 different NALPs, their effector protein ASC, and caspase-1 and -5 was readily measurable by RT-PCR in a similar proportion in RA and OA synovial samples, with the exception of NALP-5 and NALP-13, which weren't found in samples from either disease. The corresponding NALP1, -3, -12 and ASC proteins were expressed at similar levels in both OA and RA biopsies, as determined by immunohistochemistry and Western-blot analysis. By contrast, caspase-1 levels were significantly enhanced in RA synovial tissues compared to those from OA patients. NALP-1, -2, -3, -10, -12 and -14, as well as ASC, caspase-1, and -5 were detected in RNA from unstimulated and stimulated RA synoviocytes. In FLS, only ASC and caspase-1 were expressed at the protein level. NALP1, 3 and 12 were not detected. However, upon stimulation, no secreted IL-Ιβ was detectable in either RA or in OA synoviocytes culture medium. Résumé : Le but de ce projet était d'étudier le rôle de l'inflammasome dans des modèles expérimentaux d'arthrite tels que les modèles ΑΙΑ et K/BxN ainsi que dans la polyarthrite humaine (RA). La polyarthrite est une maladie inflammatoire très fréquente avec 1 % de la population affectée et touche 3 fois plus les femmes que les hommes, suggérant un rôle des hormones sexuelles dans cette pathologie. L'inflammasome est un complexe multiprotéique qui permet l'activation de la caspase-1, une cystéine protéase qui va ensuite cliver et activer rinterleukine-ΐβ (IL-Ιβ). L'inflammasome a été impliqué ces dernières années dans de nombreuses maladies inflammatoires notamment dans la goutte. Pour mettre en évidence un éventuel rôle de l'inflammasome dans l'arthrite expérimentale nous avons obtenu des souris déficientes pour certains des composants de l'inflammasome tels que ASC, NALP3 et caspase-1. Les souris NALP3 déficientes et caspase-1 déficientes sont aussi arthritiques que les souris wild type correspondantes que ce soit dans le modèle ΑΙΑ ou K/BxN. Par contre les souris mâles et femelles ASC-déficientes sont moins arthritiques que les souris +/+ correspondantes dans le modèle ΑΙΑ. Dans le modèle KRN, le même phénotype (diminution de la sévérité de l'arthrite) est observé uniquement chez les femelles ASC-/- Ce phénotype est corrélé avec l'histologie ainsi qu'avec le dosage du serum amyloid A (SAA) qui reflète l'inflammation systémique et qui est diminué chez les souris ASC-déficientes. Nous avons ensuite étudié le rôle de Γ estradiol (une des formes active des estrogènes) dans les modèles K/BxN et ΑΙΑ. Les souris castrées maies ou femelles déficientes pour ASC ayant reçu de l'estradiol ont une arthrite moins sévère ce qui implique que les estradiol ont un effet protecteur en l'absence de ASC. Dans le modèle ΑΙΑ, nous nous sommes aussi intéressés à la réponse immune. Des tests de prolifération ont été effectués sur des splénocytes en présence de mBSA (qui est l'antigène utilisé dans le modèle ΑΙΑ). Les splénocytes ASC -/- ont une proliferation qui est diminuée en présence de l'antigène. De plus la proliferation de cellules Τ spléniques purifiées en présence d'anti-CD3 est diminuée chez les cellules Τ ASC-/-. Ces résultats nous indiquent une éventuelle implication de ASC dans la signalisation par le récépteur des cellules T. En parallèle l'expression des différents composants de l'inflammasome a été analysée dans des biopsies de patients atteints de polyarthrite rhumatoide (RA) et d'arthrose (OA). L'expression des 14 différents NALPs, de l'adaptateur ASC, ainsi que des caspase-1 et -5 était similaires dans les échantillons RA et OA, à l'exception de NALP5 et 13 qui n'étaient pas détéctables. L'expression protéique de NALP1, 3, 12 et ASC effectuée par Western blot et immunohistochimie était similaire dans les biopsies RA et OA. Par contre la quantité de la caspase-1 mesurée par ELISA était augmentée de façon significative dans les extraits protéiques de biopsies RA. NALP-1, -2. -3, -10, -12, and -14 ainsi que ASC, caspase-1 et -5 étaient exprimés de façon similaire par les synoviocytes RA non stimulés et stimulés. Dans les synoviocytes seuls ASC et caspase-1 étaient détéctable au niveau protéique. NALP-1, -3 et -12 n'était pas détéctables. Cependant après stimulation il n'y avait d'IL-Ιβ sécrété que ce soit dans les surnageants de cultures de synoviocytes RA ou OA.

Parallel evolution of facial stripe patterns in the Neolamprologus brichardi/pulcher species complex endemic to Lake Tanganyika.

Colour pattern diversity can be due to random processes or to natural or sexual selection. Consequently, similarities in colour patterns are not always correlated with common ancestry, but may result from convergent evolution under shared selection pressures or drift. Neolamprologus brichardi and Neolamprologus pulcher have been described as two distinct species based on differences in the arrangement of two dark bars on the operculum. Our study uses DNA sequences of the mitochondrial control region to show that relatedness of haplotypes disagrees with species assignment based on head colour pattern. This suggests repeated parallel evolution of particular stripe patterns. The complete lack of shared haplotypes between populations of the same or different phenotypes reflects strong philopatric behaviour, possibly induced by the cooperative breeding mode in which offspring remain in their natal territory and serve as helpers until they disperse to nearby territories or take over a breeding position. Concordant phylogeographic patterns between N. brichardi/N. pulcher populations and other rock-dwelling cichlids suggest that the same colonization routes have been taken by sympatric species and that these routes were affected by lake level fluctuations in the past.

Three mathematical models for bucking-to-order

Abstract

The Caribbean-North America-Cocos Triple Junction and the dynamics of the Polochic-Motagua fault systems: Pull-up and zipper models

The Polochic-Motagua fault systems (PMFS) are part of the sinistral transform boundary between the North American and Caribbean plates. To the west, these systems interact with the subduction zone of the Cocos plate, forming a subduction-subduction-transform triple junction. The North American plate moves westward relative to the Caribbean plate. This movement does not affect the geometry of the subducted Cocos plate, which implies that deformation is accommodated entirely in the two overriding plates. Structural data, fault kinematic analysis, and geomorphic observations provide new elements that help to understand the late Cenozoic evolution of this triple junction. In the Miocene, extension and shortening occurred south and north of the Motagua fault, respectively. This strain regime migrated northward to the Polochic fault after the late Miocene. This shift is interpreted as a ``pull-up'' of North American blocks into the Caribbean realm. To the west, the PMFS interact with a trench-parallel fault zone that links the Tonala fault to the Jalpatagua fault. These faults bound a fore-arc sliver that is shared by the two overriding plates. We propose that the dextral Jalpatagua fault merges with the sinistral PMFS, leaving behind a suturing structure, the Tonala fault. This tectonic ``zipper'' allows the migration of the triple junction. As a result, the fore-arc sliver comes into contact with the North American plate and helps to maintain a linear subduction zone along the trailing edge of the Caribbean plate. All these processes currently make the triple junction increasingly diffuse as it propagates eastward and inland within both overriding plates.

Towards Machine Consciousness: Grounding Abstract Models as pi-Processes

We present a two-level model of concurrent communicating systems (CCS) to serve as a basis formachine consciousness. A language implementing threads within logic programming is ¯rstintroduced. This high-level framework allows for the de¯nition of abstract processes that can beexecuted on a virtual machine. We then look for a possible grounding of these processes into thebrain. Towards this end, we map abstract de¯nitions (including logical expressions representingcompiled knowledge) into a variant of the pi-calculus. We illustrate this approach through aseries of examples extending from a purely reactive behavior to patterns of consciousness.

Field data and prediction models of pesticide spray drift on coffee crop

The objective of this work was to generate drift curves from pesticide applications on coffee plants and to compare them with two European drift-prediction models. The used methodology is based on the ISO 22866 standard. The experimental design was a randomized complete block with ten replicates in a 2x20 split-plot arrangement. The evaluated factors were: two types of nozzles (hollow cone with and without air induction) and 20 parallel distances to the crop line outside of the target area, spaced at 2.5 m. Blotting papers were used as a target and placed in each of the evaluated distances. The spray solution was composed of water+rhodamine B fluorescent tracer at a concentration of 100 mg L-1, for detection by fluorimetry. A spray volume of 400 L ha-1 was applied using a hydropneumatic sprayer. The air-induction nozzle reduces the drift up to 20 m from the treated area. The application with the hollow cone nozzle results in 6.68% maximum drift in the nearest collector of the treated area. The German and Dutch models overestimate the drift at distances closest to the crop, although the Dutch model more closely approximates the drift curves generated by both spray nozzles.

Optimisation of data analysis method for a fluidized bed combustion test rig

Nowadays the used fuel variety in power boilers is widening and new boiler constructions and running models have to be developed. This research and development is done in small pilot plants where more faster analyse about the boiler mass and heat balance is needed to be able to find and do the right decisions already during the test run. The barrier on determining boiler balance during test runs is the long process of chemical analyses of collected input and outputmatter samples. The present work is concentrating on finding a way to determinethe boiler balance without chemical analyses and optimise the test rig to get the best possible accuracy for heat and mass balance of the boiler. The purpose of this work was to create an automatic boiler balance calculation method for 4 MW CFB/BFB pilot boiler of Kvaerner Pulping Oy located in Messukylä in Tampere. The calculation was created in the data management computer of pilot plants automation system. The calculation is made in Microsoft Excel environment, which gives a good base and functions for handling large databases and calculations without any delicate programming. The automation system in pilot plant was reconstructed und updated by Metso Automation Oy during year 2001 and the new system MetsoDNA has good data management properties, which is necessary for big calculations as boiler balance calculation. Two possible methods for calculating boiler balance during test run were found. Either the fuel flow is determined, which is usedto calculate the boiler's mass balance, or the unburned carbon loss is estimated and the mass balance of the boiler is calculated on the basis of boiler's heat balance. Both of the methods have their own weaknesses, so they were constructed parallel in the calculation and the decision of the used method was left to user. User also needs to define the used fuels and some solid mass flowsthat aren't measured automatically by the automation system. With sensitivity analysis was found that the most essential values for accurate boiler balance determination are flue gas oxygen content, the boiler's measured heat output and lower heating value of the fuel. The theoretical part of this work concentrates in the error management of these measurements and analyses and on measurement accuracy and boiler balance calculation in theory. The empirical part of this work concentrates on the creation of the balance calculation for the boiler in issue and on describing the work environment.

Siltanosturin teräsrakenteen suunnittelun avuksi luodun ohjelman jatkokehitys

Diplomityössä tehdään jatkokehitystä KCI Konecranes yrityksen siltanosturin laskentaohjelmaan. Ohjelman tärkeimmät jatkokehityskohteet kartoitettiin käyttäjäkyselyn avulla ja niistä valittiin toivotuimmat, sekä diplomityön lujuusopilliseen aihepiiriin parhaiten soveltuvat. Työhön valitut kaksi aihetta ovat koteloprofiilin kaksiosaisen uuman lujuuslaskennan selvittäminen ja siltanosturin kahdeksanpyöräisenpäätykannattajan elementtimallin suunnittelu. Diplomityössä selvitetään jatkokehityskohteisiin liittyvä teoria, mutta varsinainen ohjelmointi jätetään työn ulkopuolelle. Kaksiosaisella uumalla varustetussa koteloprofiilissa nostovaunun kulkukiskon alla olevan uuman yläosa tehdään paksummaksi, jotta uuma kestäisi nostovaunun pyöräkuormasta aiheutuvan paikallisen jännityksen, eliniin sanotun rusennusjännityksen. Rusennusjännityksen määrittäminen uumalevyissä on kaksiosaisen uuman lujuuslaskennan tärkein tehtävä. Rusennuksen aiheuttamankalvojännityksen ja jännityskeskittymien määrittämiseen erilaisissa konstruktioissa etsittiin sopivimmat menetelmät kirjallisuudesta ja standardeista. Kalvojännitys voidaan määrittää luotettavasti käyttäen joko 45 asteen sääntöä tai standardin mukaista menetelmää ja jännityskonsentraatioiden suuruus saadaan kertomallakalvojännitys jännityskonsentraatiokertoimilla. Menetelmien toimivuus verifioitiin tekemällä kymmeniä uuman elementtimalleja erilaisin dimensioin ja reunaehdoin ja vertaamalla elementtimallien tuloksia käsin laskettuihin. Käsin lasketut jännitykset saatiin vastaamaan tarkasti elementtimallien tuloksia. Kaksiosaisen uuman lommahdus- ja väsymislaskentaa tutkittiin alustavasti. Kahdeksanpyöräisiä päätykannattajia käytetään suurissa siltanostureissa pienentämään pyöräkuormia ja radan rusennusjännityksiä. Kahdeksanpyöräiselle siltanosturin päätykannattajalle suunniteltiin elementtimallit molempiin rakenteesta käytettyihin konstruktioihin: nivelöityyn ja jäykkäkehäiseen malliin. Elementtimallien rakentamisessa hyödynnettiin jo olemassa olevia malleja, jolloin niiden lisääminen ohjelmakoodiin nopeutuu ja ne ovat varmasti yhteensopivia muiden laskentamoduuleiden kanssa. Elementtimallien värähtelyanalyysin reunaehtoja tarkasteltiin. Värähtelyanalyysin reunaehtoihin ei tutkimuksen perusteella tarvitse tehdä muutoksia, mutta staattisen analyysin reunaehdot kaipaavat vielä lisätutkimusta.

Introducing massively parallel computers into operational weather forecasting

Numerical weather prediction and climate simulation have been among the computationally most demanding applications of high performance computing eversince they were started in the 1950's. Since the 1980's, the most powerful computers have featured an ever larger number of processors. By the early 2000's, this number is often several thousand. An operational weather model must use all these processors in a highly coordinated fashion. The critical resource in running such models is not computation, but the amount of necessary communication between the processors. The communication capacity of parallel computers often fallsfar short of their computational power. The articles in this thesis cover fourteen years of research into how to harness thousands of processors on a single weather forecast or climate simulation, so that the application can benefit as much as possible from the power of parallel high performance computers. The resultsattained in these articles have already been widely applied, so that currently most of the organizations that carry out global weather forecasting or climate simulation anywhere in the world use methods introduced in them. Some further studies extend parallelization opportunities into other parts of the weather forecasting environment, in particular to data assimilation of satellite observations.

Combining Generic Programming and Service-Oriented Architectures for the Effective and Timely Development of Complex e-Learning Systems

Peer-reviewed

Parallel-Operating Three-Phase Voltage Source Inverters – Circulating CurrentModeling, Analysis and Mitigation

The maximum realizable power throughput of power electronic converters may be limited or constrained by technical or economical considerations. One solution to this problemis to connect several power converter units in parallel. The parallel connection can be used to increase the current carrying capacity of the overall system beyond the ratings of individual power converter units. Thus, it is possible to use several lower-power converter units, produced in large quantities, as building blocks to construct high-power converters in a modular manner. High-power converters realized by using parallel connection are needed for example in multimegawatt wind power generation systems. Parallel connection of power converter units is also required in emerging applications such as photovoltaic and fuel cell power conversion. The parallel operation of power converter units is not, however, problem free. This is because parallel-operating units are subject to overcurrent stresses, which are caused by unequal load current sharing or currents that flow between the units. Commonly, the term ’circulatingcurrent’ is used to describe both the unequal load current sharing and the currents flowing between the units. Circulating currents, again, are caused by component tolerances and asynchronous operation of the parallel units. Parallel-operating units are also subject to stresses caused by unequal thermal stress distribution. Both of these problemscan, nevertheless, be handled with a proper circulating current control. To design an effective circulating current control system, we need information about circulating current dynamics. The dynamics of the circulating currents can be investigated by developing appropriate mathematical models. In this dissertation, circulating current models aredeveloped for two different types of parallel two-level three-phase inverter configurations. Themodels, which are developed for an arbitrary number of parallel units, provide a framework for analyzing circulating current generation mechanisms and developing circulating current control systems. In addition to developing circulating current models, modulation of parallel inverters is considered. It is illustrated that depending on the parallel inverter configuration and the modulation method applied, common-mode circulating currents may be excited as a consequence of the differential-mode circulating current control. To prevent the common-mode circulating currents that are caused by the modulation, a dual modulator method is introduced. The dual modulator basically consists of two independently operating modulators, the outputs of which eventually constitute the switching commands of the inverter. The two independently operating modulators are referred to as primary and secondary modulators. In its intended usage, the same voltage vector is fed to the primary modulators of each parallel unit, and the inputs of the secondary modulators are obtained from the circulating current controllers. To ensure that voltage commands obtained from the circulating current controllers are realizable, it must be guaranteed that the inverter is not driven into saturation by the primary modulator. The inverter saturation can be prevented by limiting the inputs of the primary and secondary modulators. Because of this, also a limitation algorithm is proposed. The operation of both the proposed dual modulator and the limitation algorithm is verified experimentally.