Solving mixed integer nonlinear programming problems using spiral dynamics optimization algorithm
NASA Astrophysics Data System (ADS)
Kania, Adhe; Sidarto, Kuntjoro Adji
2016-02-01
Many engineering and practical problem can be modeled by mixed integer nonlinear programming. This paper proposes to solve the problem with modified spiral dynamics inspired optimization method of Tamura and Yasuda. Four test cases have been examined, including problem in engineering and sport. This method succeeds in obtaining the optimal result in all test cases.
Mixed integer evolution strategies for parameter optimization.
Li, Rui; Emmerich, Michael T M; Eggermont, Jeroen; Bäck, Thomas; Schütz, M; Dijkstra, J; Reiber, J H C
2013-01-01
Evolution strategies (ESs) are powerful probabilistic search and optimization algorithms gleaned from biological evolution theory. They have been successfully applied to a wide range of real world applications. The modern ESs are mainly designed for solving continuous parameter optimization problems. Their ability to adapt the parameters of the multivariate normal distribution used for mutation during the optimization run makes them well suited for this domain. In this article we describe and study mixed integer evolution strategies (MIES), which are natural extensions of ES for mixed integer optimization problems. MIES can deal with parameter vectors consisting not only of continuous variables but also with nominal discrete and integer variables. Following the design principles of the canonical evolution strategies, they use specialized mutation operators tailored for the aforementioned mixed parameter classes. For each type of variable, the choice of mutation operators is governed by a natural metric for this variable type, maximal entropy, and symmetry considerations. All distributions used for mutation can be controlled in their shape by means of scaling parameters, allowing self-adaptation to be implemented. After introducing and motivating the conceptual design of the MIES, we study the optimality of the self-adaptation of step sizes and mutation rates on a generalized (weighted) sphere model. Moreover, we prove global convergence of the MIES on a very general class of problems. The remainder of the article is devoted to performance studies on artificial landscapes (barrier functions and mixed integer NK landscapes), and a case study in the optimization of medical image analysis systems. In addition, we show that with proper constraint handling techniques, MIES can also be applied to classical mixed integer nonlinear programming problems. PMID:22122384
Henriques, David; Rocha, Miguel; Saez-Rodriguez, Julio; Banga, Julio R.
2015-01-01
Motivation: Systems biology models can be used to test new hypotheses formulated on the basis of previous knowledge or new experimental data, contradictory with a previously existing model. New hypotheses often come in the shape of a set of possible regulatory mechanisms. This search is usually not limited to finding a single regulation link, but rather a combination of links subject to great uncertainty or no information about the kinetic parameters. Results: In this work, we combine a logic-based formalism, to describe all the possible regulatory structures for a given dynamic model of a pathway, with mixed-integer dynamic optimization (MIDO). This framework aims to simultaneously identify the regulatory structure (represented by binary parameters) and the real-valued parameters that are consistent with the available experimental data, resulting in a logic-based differential equation model. The alternative to this would be to perform real-valued parameter estimation for each possible model structure, which is not tractable for models of the size presented in this work. The performance of the method presented here is illustrated with several case studies: a synthetic pathway problem of signaling regulation, a two-component signal transduction pathway in bacterial homeostasis, and a signaling network in liver cancer cells. Supplementary information: Supplementary data are available at Bioinformatics online. Contact: julio@iim.csic.es or saezrodriguez@ebi.ac.uk PMID:26002881
NASA Astrophysics Data System (ADS)
Wang, Bin; Chiang, Hsiao-Dong
Many applications of smart grid can be formulated as constrained optimization problems. Because of the discrete controls involved in power systems, these problems are essentially mixed-integer nonlinear programs. In this paper, we review the Trust-Tech-based methodology for solving mixed-integer nonlinear optimization. Specifically, we have developed a two-stage Trust-Tech-based methodology to systematically compute all the local optimal solutions for constrained mixed-integer nonlinear programming (MINLP) problems. In the first stage, for a given MINLP problem this methodology starts with the construction of a new, continuous, unconstrained problem through relaxation and the penalty function method. A corresponding dynamical system is then constructed to search for a set of local optimal solutions for the unconstrained problem. In the second stage, a reduced constrained NLP is defined for each local optimal solution by determining and fixing the values of integral variables of the MINLP problem. The Trust-Tech-based method is used to compute a set of local optimal solutions for these reduced NLP problems, from which the optimal solution of the original MINLP problem is determined. A numerical simulation of several testing problems is provided to illustrate the effectiveness of our proposed method.
Goodman, G.V.R.
1987-01-01
The lack of available techniques prompted the development of a mixed integer model to optimize the scheduling of equipment and the distribution of overburden in a typical mountaintop removal operation. Using this format, a (0-1) integer model and transportation model were constructed to determine the optimal equipment schedule and optimal overburden distribution, respectively. To solve this mixed integer program, the model was partitioned into its binary and real-valued components. Each problem was successively solved and their values added to form estimates of the value of the mixed integer program. Optimal convergence was indicated when the difference between two successive estimates satisfied some pre-specific accuracy value. The performance of the mixed integer model was tested against actual field data to determine its practical applications. To provide the necessary input information, production data was obtained from a single seam, mountaintop removal operation located in the Appalachian coal field. As a means of analyzing the resultant equipment schedule, the total idle time was calculated for each machine type and each lift location. Also, the final overburden assignments were analyzed by determining the distribution of spoil material for various overburden removal productivities. Subsequent validation of the mixed integer model was conducted in two distinct areas. The first dealt with changes in algorithmic data and their effects on the optimality of the model. The second area concerned variations in problem structure, specifically those dealing with changes in problem size and other user-inputed values such as equipment productivities or required reclamation.
Enhanced index tracking modeling in portfolio optimization with mixed-integer programming z approach
NASA Astrophysics Data System (ADS)
Siew, Lam Weng; Jaaman, Saiful Hafizah Hj.; Ismail, Hamizun bin
2014-09-01
Enhanced index tracking is a popular form of portfolio management in stock market investment. Enhanced index tracking aims to construct an optimal portfolio to generate excess return over the return achieved by the stock market index without purchasing all of the stocks that make up the index. The objective of this paper is to construct an optimal portfolio using mixed-integer programming model which adopts regression approach in order to generate higher portfolio mean return than stock market index return. In this study, the data consists of 24 component stocks in Malaysia market index which is FTSE Bursa Malaysia Kuala Lumpur Composite Index from January 2010 until December 2012. The results of this study show that the optimal portfolio of mixed-integer programming model is able to generate higher mean return than FTSE Bursa Malaysia Kuala Lumpur Composite Index return with only selecting 30% out of the total stock market index components.
Optimization of a wood dryer kiln using the mixed integer programming technique: A case study
Gustafsson, S.I.
1999-07-01
When wood is to be utilized as a raw material for furniture, buildings, etc., it must be dried from approximately 100% to 6% moisture content. This is achieved at least partly in a drying kiln. Heat for this purpose is provided by electrical means, or by steam from boilers fired with wood chips or oil. By making a close examination of monitored values from an actual drying kiln it has been possible to optimize the use of steam and electricity using the so called mixed integer programming technique. Owing to the operating schedule for the drying kiln it has been necessary to divide the drying process in very short time intervals, i.e., a number of minutes. Since a drying cycle takes about two or three weeks, a considerable mathematical problem is presented and this has to be solved.
A Mixed-Integer Optimization Framework for De Novo Peptide Identification
DiMaggio, Peter A.
2009-01-01
A novel methodology for the de novo identification of peptides by mixed-integer optimization and tandem mass spectrometry is presented in this article. The various features of the mathematical model are presented and examples are used to illustrate the key concepts of the proposed approach. Several problems are examined to illustrate the proposed method's ability to address (1) residue-dependent fragmentation properties and (2) the variability of resolution in different mass analyzers. A preprocessing algorithm is used to identify important m/z values in the tandem mass spectrum. Missing peaks, resulting from residue-dependent fragmentation characteristics, are dealt with using a two-stage algorithmic framework. A cross-correlation approach is used to resolve missing amino acid assignments and to identify the most probable peptide by comparing the theoretical spectra of the candidate sequences that were generated from the MILP sequencing stages with the experimental tandem mass spectrum. PMID:19412358
Passenger-Oriented Optimization for Train Rescheduling on the Basis of Mixed Integer Programming
NASA Astrophysics Data System (ADS)
Chigusa, Kenji; Sato, Keisuke; Koseki, Takafumi
Whenever there is a disruption in areas with a high frequency of train operations, the delays of trains propagate and train passengers have to suffer inconvenience. A rescheduling of the train operations is then conducted by train dispatchers to revert to the original operational plan. The task is quite difficult and a severe burden for the dispatchers because of a variety of options for reverting to the original plan, the necessity for a swift action, and the absence of a proper evaluation criterion. Consequently, a computer-aided rescheduling-support system is required. In this paper, we present a mathematical programming approach for train rescheduling, focusing on minimizing the passengers' arrival delay at their destinations. We simultaneously model the train operations and passenger flow on the basis of mixed integer programming (MIP), and we obtain an optimal rescheduling plan in reasonable time in the case of a small disruption.
NASA Astrophysics Data System (ADS)
Baran, Ismet; Tutum, Cem C.; Hattel, Jesper H.
2013-08-01
In this paper thermo-chemical simulation of the pultrusion process of a composite rod is first used as a validation case to ensure that the utilized numerical scheme is stable and converges to results given in literature. Following this validation case, a cylindrical die block with heaters is added to the pultrusion domain of a composite part and thermal contact resistance (TCR) regions at the die-part interface are defined. Two optimization case studies are performed on this new configuration. In the first one, optimal die radius and TCR values are found by using a hybrid genetic algorithm based on a sequential combination of a genetic algorithm (GA) and a local search technique to fit the centerline temperature of the composite with the one calculated in the validation case. In the second optimization study, the productivity of the process is improved by using a mixed integer genetic algorithm (MIGA) such that the total number of heaters is minimized while satisfying the constraints for the maximum composite temperature, the mean of the cure degree at the die exit and the pulling speed.
Designing cost-effective biopharmaceutical facilities using mixed-integer optimization.
Liu, Songsong; Simaria, Ana S; Farid, Suzanne S; Papageorgiou, Lazaros G
2013-01-01
Chromatography operations are identified as critical steps in a monoclonal antibody (mAb) purification process and can represent a significant proportion of the purification material costs. This becomes even more critical with increasing product titers that result in higher mass loads onto chromatography columns, potentially causing capacity bottlenecks. In this work, a mixed-integer nonlinear programming (MINLP) model was created and applied to an industrially relevant case study to optimize the design of a facility by determining the most cost-effective chromatography equipment sizing strategies for the production of mAbs. Furthermore, the model was extended to evaluate the ability of a fixed facility to cope with higher product titers up to 15 g/L. Examination of the characteristics of the optimal chromatography sizing strategies across different titer values enabled the identification of the maximum titer that the facility could handle using a sequence of single column chromatography steps as well as multi-column steps. The critical titer levels for different ratios of upstream to dowstream trains where multiple parallel columns per step resulted in the removal of facility bottlenecks were identified. Different facility configurations in terms of number of upstream trains were considered and the trade-off between their cost and ability to handle higher titers was analyzed. The case study insights demonstrate that the proposed modeling approach, combining MINLP models with visualization tools, is a valuable decision-support tool for the design of cost-effective facility configurations and to aid facility fit decisions. 2013. PMID:23956206
NASA Astrophysics Data System (ADS)
Tang, Jiafu; Liu, Yang; Fung, Richard; Luo, Xinggang
2008-12-01
Manufacturers have a legal accountability to deal with industrial waste generated from their production processes in order to avoid pollution. Along with advances in waste recovery techniques, manufacturers may adopt various recycling strategies in dealing with industrial waste. With reuse strategies and technologies, byproducts or wastes will be returned to production processes in the iron and steel industry, and some waste can be recycled back to base material for reuse in other industries. This article focuses on a recovery strategies optimization problem for a typical class of industrial waste recycling process in order to maximize profit. There are multiple strategies for waste recycling available to generate multiple byproducts; these byproducts are then further transformed into several types of chemical products via different production patterns. A mixed integer programming model is developed to determine which recycling strategy and which production pattern should be selected with what quantity of chemical products corresponding to this strategy and pattern in order to yield maximum marginal profits. The sales profits of chemical products and the set-up costs of these strategies, patterns and operation costs of production are considered. A simulated annealing (SA) based heuristic algorithm is developed to solve the problem. Finally, an experiment is designed to verify the effectiveness and feasibility of the proposed method. By comparing a single strategy to multiple strategies in an example, it is shown that the total sales profit of chemical products can be increased by around 25% through the simultaneous use of multiple strategies. This illustrates the superiority of combinatorial multiple strategies. Furthermore, the effects of the model parameters on profit are discussed to help manufacturers organize their waste recycling network.
Mixed integer programming model for optimizing the layout of an ICU vehicle
2009-01-01
Background This paper presents a Mixed Integer Programming (MIP) model for designing the layout of the Intensive Care Units' (ICUs) patient care space. In particular, this MIP model was developed for optimizing the layout for materials to be used in interventions. This work was developed within the framework of a joint project between the Madrid Technical Unverstity and the Medical Emergency Services of the Madrid Regional Government (SUMMA 112). Methods The first task was to identify the relevant information to define the characteristics of the new vehicles and, in particular, to obtain a satisfactory interior layout to locate all the necessary materials. This information was gathered from health workers related to ICUs. With that information an optimization model was developed in order to obtain a solution. From the MIP model, a first solution was obtained, consisting of a grid to locate the different materials needed for the ICUs. The outcome from the MIP model was discussed with health workers to tune the solution, and after slightly altering that solution to meet some requirements that had not been included in the mathematical model, the eventual solution was approved by the persons responsible for specifying the characteristics of the new vehicles. According to the opinion stated by the SUMMA 112's medical group responsible for improving the ambulances (the so-called "coaching group"), the outcome was highly satisfactory. Indeed, the final design served as a basis to draw up the requirements of a public tender. Results As a result from solving the Optimization model, a grid was obtained to locate the different necessary materials for the ICUs. This grid had to be slightly altered to meet some requirements that had not been included in the mathematical model. The results were discussed with the persons responsible for specifying the characteristics of the new vehicles. Conclusion The outcome was highly satisfactory. Indeed, the final design served as a basis to draw up the requirements of a public tender. The authors advocate this approach to address similar problems within the field of Health Services to improve the efficiency and the effectiveness of the processes involved. Problems such as those in operation rooms or emergency rooms, where the availability of a large amount of material is critical are eligible to be dealt with in a simmilar manner. PMID:19995438
NASA Astrophysics Data System (ADS)
Skulovich, Olya; Bent, Russell; Judi, David; Perelman, Lina Sela; Ostfeld, Avi
2015-06-01
Despite their potential catastrophic impact, transients are often ignored or presented ad hoc when designing water distribution systems. To address this problem, we introduce a new piece-wise function fitting model that is integrated with mixed integer programming to optimally place and size surge tanks for transient control. The key features of the algorithm are a model-driven discretization of the search space, a linear approximation nonsmooth system response surface to transients, and a mixed integer linear programming optimization. Results indicate that high quality solutions can be obtained within a reasonable number of function evaluations and demonstrate the computational effectiveness of the approach through two case studies. The work investigates one type of surge control devices (closed surge tank) for a specified set of transient events. The performance of the algorithm relies on the assumption that there exists a smooth relationship between the objective function and tank size. Results indicate the potential of the approach for the optimal surge control design in water systems.
Winebrake, James J; Corbett, James J; Wang, Chengfeng; Farrell, Alexander E; Woods, Pippa
2005-04-01
Emissions from passenger ferries operating in urban harbors may contribute significantly to emissions inventories and commuter exposure to air pollution. In particular, ferries are problematic because of high emissions of oxides of nitrogen (NOx) and particulate matter (PM) from primarily unregulated diesel engines. This paper explores technical solutions to reduce pollution from passenger ferries operating in the New York-New Jersey Harbor. The paper discusses and demonstrates a mixed-integer, non-linear programming model used to identify optimal control strategies for meeting NOx and PM reduction targets for 45 privately owned commuter ferries in the harbor. Results from the model can be used by policy-makers to craft programs aimed at achieving least-cost reduction targets. PMID:15887889
Mixed Integer Programming and Heuristic Scheduling for Space Communication Networks
NASA Technical Reports Server (NTRS)
Cheung, Kar-Ming; Lee, Charles H.
2012-01-01
We developed framework and the mathematical formulation for optimizing communication network using mixed integer programming. The design yields a system that is much smaller, in search space size, when compared to the earlier approach. Our constrained network optimization takes into account the dynamics of link performance within the network along with mission and operation requirements. A unique penalty function is introduced to transform the mixed integer programming into the more manageable problem of searching in a continuous space. The constrained optimization problem was proposed to solve in two stages: first using the heuristic Particle Swarming Optimization algorithm to get a good initial starting point, and then feeding the result into the Sequential Quadratic Programming algorithm to achieve the final optimal schedule. We demonstrate the above planning and scheduling methodology with a scenario of 20 spacecraft and 3 ground stations of a Deep Space Network site. Our approach and framework have been simple and flexible so that problems with larger number of constraints and network can be easily adapted and solved.
Ko, Andi Setiady; Chang, Ni-Bin
2008-07-01
Energy supply and use is of fundamental importance to society. Although the interactions between energy and environment were originally local in character, they have now widened to cover regional and global issues, such as acid rain and the greenhouse effect. It is for this reason that there is a need for covering the direct and indirect economic and environmental impacts of energy acquisition, transport, production and use. In this paper, particular attention is directed to ways of resolving conflict between economic and environmental goals by encouraging a power plant to consider co-firing biomass and refuse-derived fuel (RDF) with coal simultaneously. It aims at reducing the emission level of sulfur dioxide (SO(2)) in an uncertain environment, using the power plant in Michigan City, Indiana as an example. To assess the uncertainty by a comparative way both deterministic and grey nonlinear mixed integer programming (MIP) models were developed to minimize the net operating cost with respect to possible fuel combinations. It aims at generating the optimal portfolio of alternative fuels while maintaining the same electricity generation simultaneously. To ease the solution procedure stepwise relaxation algorithm was developed for solving the grey nonlinear MIP model. Breakeven alternative fuel value can be identified in the post-optimization stage for decision-making. Research findings show that the inclusion of RDF does not exhibit comparative advantage in terms of the net cost, albeit relatively lower air pollution impact. Yet it can be sustained by a charge system, subsidy program, or emission credit as the price of coal increases over time. PMID:17395362
Mixed Integer Programming and Heuristic Scheduling for Space Communication
NASA Technical Reports Server (NTRS)
Lee, Charles H.; Cheung, Kar-Ming
2013-01-01
Optimal planning and scheduling for a communication network was created where the nodes within the network are communicating at the highest possible rates while meeting the mission requirements and operational constraints. The planning and scheduling problem was formulated in the framework of Mixed Integer Programming (MIP) to introduce a special penalty function to convert the MIP problem into a continuous optimization problem, and to solve the constrained optimization problem using heuristic optimization. The communication network consists of space and ground assets with the link dynamics between any two assets varying with respect to time, distance, and telecom configurations. One asset could be communicating with another at very high data rates at one time, and at other times, communication is impossible, as the asset could be inaccessible from the network due to planetary occultation. Based on the network's geometric dynamics and link capabilities, the start time, end time, and link configuration of each view period are selected to maximize the communication efficiency within the network. Mathematical formulations for the constrained mixed integer optimization problem were derived, and efficient analytical and numerical techniques were developed to find the optimal solution. By setting up the problem using MIP, the search space for the optimization problem is reduced significantly, thereby speeding up the solution process. The ratio of the dimension of the traditional method over the proposed formulation is approximately an order N (single) to 2*N (arraying), where N is the number of receiving antennas of a node. By introducing a special penalty function, the MIP problem with non-differentiable cost function and nonlinear constraints can be converted into a continuous variable problem, whose solution is possible.
Mixed-Integer Formulations for Constellation Scheduling
NASA Astrophysics Data System (ADS)
Valicka, C.; Hart, W.; Rintoul, M.
Remote sensing systems have expanded the set of capabilities available for and critical to national security. Cooperating, high-fidelity sensing systems and growing mission applications have exponentially increased the set of potential schedules. A definitive lack of advanced tools places an increased burden on operators, as planning and scheduling remain largely manual tasks. This is particularly true in time-critical planning activities where operators aim to accomplish a large number of missions through optimal utilization of single or multiple sensor systems. Automated scheduling through identification and comparison of alternative schedules remains a challenging problem applicable across all remote sensing systems. Previous approaches focused on a subset of sensor missions and do not consider ad-hoc tasking. We have begun development of a robust framework that leverages the Pyomo optimization modeling language for the design of a tool to assist sensor operators planning under the constraints of multiple concurrent missions and uncertainty. Our scheduling models have been formulated to address the stochastic nature of ad-hoc tasks inserted under a variety of scenarios. Operator experience is being leveraged to select appropriate model objectives. Successful development of the framework will include iterative development of high-fidelity mission models that consider and expose various schedule performance metrics. Creating this tool will aid time-critical scheduling by increasing planning efficiency, clarifying the value of alternative modalities uniquely provided by multi-sensor systems, and by presenting both sets of organized information to operators. Such a tool will help operators more quickly and fully utilize sensing systems, a high interest objective within the current remote sensing operations community. Preliminary results for mixed-integer programming formulations of a sensor scheduling problem will be presented. Assumptions regarding sensor geometry and sensing activity time constraints, durations, priorities, etc. will be outlined. Finally, solver speed and stochastic programming details for uncertain activities and scheduling impediments will be discussed.
Smalley, Hannah K; Keskinocak, Pinar; Swann, Julie; Hinman, Alan
2015-11-17
In addition to improved sanitation, hygiene, and better access to safe water, oral cholera vaccines can help to control the spread of cholera in the short term. However, there is currently no systematic method for determining the best allocation of oral cholera vaccines to minimize disease incidence in a population where the disease is endemic and resources are limited. We present a mathematical model for optimally allocating vaccines in a region under varying levels of demographic and incidence data availability. The model addresses the questions of where, when, and how many doses of vaccines to send. Considering vaccine efficacies (which may vary based on age and the number of years since vaccination), we analyze distribution strategies which allocate vaccines over multiple years. Results indicate that, given appropriate surveillance data, targeting age groups and regions with the highest disease incidence should be the first priority, followed by other groups primarily in order of disease incidence, as this approach is the most life-saving and cost-effective. A lack of detailed incidence data results in distribution strategies which are not cost-effective and can lead to thousands more deaths from the disease. The mathematical model allows for what-if analysis for various vaccine distribution strategies by providing the ability to easily vary parameters such as numbers and sizes of regions and age groups, risk levels, vaccine price, vaccine efficacy, production capacity and budget. PMID:26458806
Fast scaffolding with small independent mixed integer programs
Salmela, Leena; Mäkinen, Veli; Välimäki, Niko; Ylinen, Johannes; Ukkonen, Esko
2011-01-01
Motivation: Assembling genomes from short read data has become increasingly popular, but the problem remains computationally challenging especially for larger genomes. We study the scaffolding phase of sequence assembly where preassembled contigs are ordered based on mate pair data. Results: We present MIP Scaffolder that divides the scaffolding problem into smaller subproblems and solves these with mixed integer programming. The scaffolding problem can be represented as a graph and the biconnected components of this graph can be solved independently. We present a technique for restricting the size of these subproblems so that they can be solved accurately with mixed integer programming. We compare MIP Scaffolder to two state of the art methods, SOPRA and SSPACE. MIP Scaffolder is fast and produces better or as good scaffolds as its competitors on large genomes. Availability: The source code of MIP Scaffolder is freely available at http://www.cs.helsinki.fi/u/lmsalmel/mip-scaffolder/. Contact: leena.salmela@cs.helsinki.fi PMID:21998153
Multi-objective mixed integer strategy for the optimisation of biological networks.
Sendín, J O H; Exler, O; Banga, J R
2010-05-01
In this contribution, the authors consider multi-criteria optimisation problems arising from the field of systems biology when both continuous and integer decision variables are involved. Mathematically, they are formulated as mixed-integer non-linear programming problems. The authors present a novel solution strategy based on a global optimisation approach for dealing with this class of problems. Its usefulness and capabilities are illustrated with two metabolic engineering case studies. For these problems, the authors show how the set of optimal solutions (the so-called Pareto front) is successfully and efficiently obtained, providing further insight into the systems under consideration regarding their optimal manipulation. PMID:20500003
PySP : modeling and solving stochastic mixed-integer programs in Python.
Woodruff, David L.; Watson, Jean-Paul
2010-08-01
Although stochastic programming is a powerful tool for modeling decision-making under uncertainty, various impediments have historically prevented its widespread use. One key factor involves the ability of non-specialists to easily express stochastic programming problems as extensions of deterministic models, which are often formulated first. A second key factor relates to the difficulty of solving stochastic programming models, particularly the general mixed-integer, multi-stage case. Intricate, configurable, and parallel decomposition strategies are frequently required to achieve tractable run-times. We simultaneously address both of these factors in our PySP software package, which is part of the COIN-OR Coopr open-source Python project for optimization. To formulate a stochastic program in PySP, the user specifies both the deterministic base model and the scenario tree with associated uncertain parameters in the Pyomo open-source algebraic modeling language. Given these two models, PySP provides two paths for solution of the corresponding stochastic program. The first alternative involves writing the extensive form and invoking a standard deterministic (mixed-integer) solver. For more complex stochastic programs, we provide an implementation of Rockafellar and Wets Progressive Hedging algorithm. Our particular focus is on the use of Progressive Hedging as an effective heuristic for approximating general multi-stage, mixed-integer stochastic programs. By leveraging the combination of a high-level programming language (Python) and the embedding of the base deterministic model in that language (Pyomo), we are able to provide completely generic and highly configurable solver implementations. PySP has been used by a number of research groups, including our own, to rapidly prototype and solve difficult stochastic programming problems.
ERIC Educational Resources Information Center
Han, Kyung T.; Rudner, Lawrence M.
2014-01-01
This study uses mixed integer quadratic programming (MIQP) to construct multiple highly equivalent item pools simultaneously, and compares the results from mixed integer programming (MIP). Three different MIP/MIQP models were implemented and evaluated using real CAT item pool data with 23 different content areas and a goal of equal information…
Linderoth, Jeff T.; Luedtke, James R.
2013-05-30
The mathematical modeling of systems often requires the use of both nonlinear and discrete components. Problems involving both discrete and nonlinear components are known as mixed-integer nonlinear programs (MINLPs) and are among the most challenging computational optimization problems. This research project added to the understanding of this area by making a number of fundamental advances. First, the work demonstrated many novel, strong, tractable relaxations designed to deal with non-convexities arising in mathematical formulation. Second, the research implemented the ideas in software that is available to the public. Finally, the work demonstrated the importance of these ideas on practical applications and disseminated the work through scholarly journals, survey publications, and conference presentations.
Learning oncogenetic networks by reducing to mixed integer linear programming.
Shahrabi Farahani, Hossein; Lagergren, Jens
2013-01-01
Cancer can be a result of accumulation of different types of genetic mutations such as copy number aberrations. The data from tumors are cross-sectional and do not contain the temporal order of the genetic events. Finding the order in which the genetic events have occurred and progression pathways are of vital importance in understanding the disease. In order to model cancer progression, we propose Progression Networks, a special case of Bayesian networks, that are tailored to model disease progression. Progression networks have similarities with Conjunctive Bayesian Networks (CBNs) [1],a variation of Bayesian networks also proposed for modeling disease progression. We also describe a learning algorithm for learning Bayesian networks in general and progression networks in particular. We reduce the hard problem of learning the Bayesian and progression networks to Mixed Integer Linear Programming (MILP). MILP is a Non-deterministic Polynomial-time complete (NP-complete) problem for which very good heuristics exists. We tested our algorithm on synthetic and real cytogenetic data from renal cell carcinoma. We also compared our learned progression networks with the networks proposed in earlier publications. The software is available on the website https://bitbucket.org/farahani/diprog. PMID:23799047
Automatic design of synthetic gene circuits through mixed integer non-linear programming.
Huynh, Linh; Kececioglu, John; Köppe, Matthias; Tagkopoulos, Ilias
2012-01-01
Automatic design of synthetic gene circuits poses a significant challenge to synthetic biology, primarily due to the complexity of biological systems, and the lack of rigorous optimization methods that can cope with the combinatorial explosion as the number of biological parts increases. Current optimization methods for synthetic gene design rely on heuristic algorithms that are usually not deterministic, deliver sub-optimal solutions, and provide no guaranties on convergence or error bounds. Here, we introduce an optimization framework for the problem of part selection in synthetic gene circuits that is based on mixed integer non-linear programming (MINLP), which is a deterministic method that finds the globally optimal solution and guarantees convergence in finite time. Given a synthetic gene circuit, a library of characterized parts, and user-defined constraints, our method can find the optimal selection of parts that satisfy the constraints and best approximates the objective function given by the user. We evaluated the proposed method in the design of three synthetic circuits (a toggle switch, a transcriptional cascade, and a band detector), with both experimentally constructed and synthetic promoter libraries. Scalability and robustness analysis shows that the proposed framework scales well with the library size and the solution space. The work described here is a step towards a unifying, realistic framework for the automated design of biological circuits. PMID:22536398
NASA Technical Reports Server (NTRS)
Laird, Philip
1992-01-01
We distinguish static and dynamic optimization of programs: whereas static optimization modifies a program before runtime and is based only on its syntactical structure, dynamic optimization is based on the statistical properties of the input source and examples of program execution. Explanation-based generalization is a commonly used dynamic optimization method, but its effectiveness as a speedup-learning method is limited, in part because it fails to separate the learning process from the program transformation process. This paper describes a dynamic optimization technique called a learn-optimize cycle that first uses a learning element to uncover predictable patterns in the program execution and then uses an optimization algorithm to map these patterns into beneficial transformations. The technique has been used successfully for dynamic optimization of pure Prolog.
A Mixed Integer Linear Program for Airport Departure Scheduling
NASA Technical Reports Server (NTRS)
Gupta, Gautam; Jung, Yoon Chul
2009-01-01
Aircraft departing from an airport are subject to numerous constraints while scheduling departure times. These constraints include wake-separation constraints for successive departures, miles-in-trail separation for aircraft bound for the same departure fixes, and time-window or prioritization constraints for individual flights. Besides these, emissions as well as increased fuel consumption due to inefficient scheduling need to be included. Addressing all the above constraints in a single framework while allowing for resequencing of the aircraft using runway queues is critical to the implementation of the Next Generation Air Transport System (NextGen) concepts. Prior work on airport departure scheduling has addressed some of the above. However, existing methods use pre-determined runway queues, and schedule aircraft from these departure queues. The source of such pre-determined queues is not explicit, and could potentially be a subjective controller input. Determining runway queues and scheduling within the same framework would potentially result in better scheduling. This paper presents a mixed integer linear program (MILP) for the departure-scheduling problem. The program takes as input the incoming sequence of aircraft for departure from a runway, along with their earliest departure times and an optional prioritization scheme based on time-window of departure for each aircraft. The program then assigns these aircraft to the available departure queues and schedules departure times, explicitly considering wake separation and departure fix restrictions to minimize total delay for all aircraft. The approach is generalized and can be used in a variety of situations, and allows for aircraft prioritization based on operational as well as environmental considerations. We present the MILP in the paper, along with benefits over the first-come-first-serve (FCFS) scheme for numerous randomized problems based on real-world settings. The MILP results in substantially reduced delays as compared to FCFS, and the magnitude of the savings depends on the queue and departure fix structure. The MILP assumes deterministic aircraft arrival times at the runway queues. However, due to taxi time uncertainty, aircraft might arrive either earlier or later than these deterministic times. Thus, to incorporate this uncertainty, we present a method for using the MILP with "overlap discounted rolling planning horizon". The approach is based on valuing near-term decision results more than future ones. We develop a model of taxitime uncertainty based on real-world data, and then compare the baseline FCFS delays with delays using the above MILP in a simple rolling-horizon method and in the overlap discounted scheme.
NASA Astrophysics Data System (ADS)
Onoyama, Takashi; Kubota, Sen; Maekawa, Takuya; Komoda, Norihisa
Adequate response performance is required for the planning of a cooperative logistic network covering multiple enterprises, because this process needs a human expert's evaluation from many aspects. To satisfy this requirement, we propose an accurate model based on mixed integer programming for optimizing cooperative logistics networks where “round transportation” exists together with “depot transportation” including lower limit constraints of loading ratio for round transportation vehicles. Furthermore, to achieve interactive response performance, a dummy load is introduced into the model instead of integer variables. The experimental result shows the proposed method obtains an accurate solution within interactive response time.
A Mixed Integer Linear Program for Solving a Multiple Route Taxi Scheduling Problem
NASA Technical Reports Server (NTRS)
Montoya, Justin Vincent; Wood, Zachary Paul; Rathinam, Sivakumar; Malik, Waqar Ahmad
2010-01-01
Aircraft movements on taxiways at busy airports often create bottlenecks. This paper introduces a mixed integer linear program to solve a Multiple Route Aircraft Taxi Scheduling Problem. The outputs of the model are in the form of optimal taxi schedules, which include routing decisions for taxiing aircraft. The model extends an existing single route formulation to include routing decisions. An efficient comparison framework compares the multi-route formulation and the single route formulation. The multi-route model is exercised for east side airport surface traffic at Dallas/Fort Worth International Airport to determine if any arrival taxi time savings can be achieved by allowing arrivals to have two taxi routes: a route that crosses an active departure runway and a perimeter route that avoids the crossing. Results indicate that the multi-route formulation yields reduced arrival taxi times over the single route formulation only when a perimeter taxiway is used. In conditions where the departure aircraft are given an optimal and fixed takeoff sequence, accumulative arrival taxi time savings in the multi-route formulation can be as high as 3.6 hours more than the single route formulation. If the departure sequence is not optimal, the multi-route formulation results in less taxi time savings made over the single route formulation, but the average arrival taxi time is significantly decreased.
A Two-Stage Stochastic Mixed-Integer Programming Approach to the Smart House Scheduling Problem
NASA Astrophysics Data System (ADS)
Ozoe, Shunsuke; Tanaka, Yoichi; Fukushima, Masao
A “Smart House” is a highly energy-optimized house equipped with photovoltaic systems (PV systems), electric battery systems, fuel cell cogeneration systems (FC systems), electric vehicles (EVs) and so on. Smart houses are attracting much attention recently thanks to their enhanced ability to save energy by making full use of renewable energy and by achieving power grid stability despite an increased power draw for installed PV systems. Yet running a smart house's power system, with its multiple power sources and power storages, is no simple task. In this paper, we consider the problem of power scheduling for a smart house with a PV system, an FC system and an EV. We formulate the problem as a mixed integer programming problem, and then extend it to a stochastic programming problem involving recourse costs to cope with uncertain electricity demand, heat demand and PV power generation. Using our method, we seek to achieve the optimal power schedule running at the minimum expected operation cost. We present some results of numerical experiments with data on real-life demands and PV power generation to show the effectiveness of our method.
Synchronic interval Gaussian mixed-integer programming for air quality management.
Cheng, Guanhui; Huang, Guohe Gordon; Dong, Cong
2015-12-15
To reveal the synchronism of interval uncertainties, the tradeoff between system optimality and security, the discreteness of facility-expansion options, the uncertainty of pollutant dispersion processes, and the seasonality of wind features in air quality management (AQM) systems, a synchronic interval Gaussian mixed-integer programming (SIGMIP) approach is proposed in this study. A robust interval Gaussian dispersion model is developed for approaching the pollutant dispersion process under interval uncertainties and seasonal variations. The reflection of synchronic effects of interval uncertainties in the programming objective is enabled through introducing interval functions. The proposition of constraint violation degrees helps quantify the tradeoff between system optimality and constraint violation under interval uncertainties. The overall optimality of system profits of an SIGMIP model is achieved based on the definition of an integrally optimal solution. Integer variables in the SIGMIP model are resolved by the existing cutting-plane method. Combining these efforts leads to an effective algorithm for the SIGMIP model. An application to an AQM problem in a region in Shandong Province, China, reveals that the proposed SIGMIP model can facilitate identifying the desired scheme for AQM. The enhancement of the robustness of optimization exercises may be helpful for increasing the reliability of suggested schemes for AQM under these complexities. The interrelated tradeoffs among control measures, emission sources, flow processes, receptors, influencing factors, and economic and environmental goals are effectively balanced. Interests of many stakeholders are reasonably coordinated. The harmony between economic development and air quality control is enabled. Results also indicate that the constraint violation degree is effective at reflecting the compromise relationship between constraint-violation risks and system optimality under interval uncertainties. This can help decision makers mitigate potential risks, e.g. insufficiency of pollutant treatment capabilities, exceedance of air quality standards, deficiency of pollution control fund, or imbalance of economic or environmental stress, in the process of guiding AQM. PMID:26367068
Mixed Integer Programming and Heuristic Scheduling for Space Communication Networks
NASA Technical Reports Server (NTRS)
Lee, Charles H.; Cheung, Kar-Ming
2012-01-01
In this paper, we propose to solve the constrained optimization problem in two phases. The first phase uses heuristic methods such as the ant colony method, particle swarming optimization, and genetic algorithm to seek a near optimal solution among a list of feasible initial populations. The final optimal solution can be found by using the solution of the first phase as the initial condition to the SQP algorithm. We demonstrate the above problem formulation and optimization schemes with a large-scale network that includes the DSN ground stations and a number of spacecraft of deep space missions.
Solution of Mixed-Integer Programming Problems on the XT5
Hartman-Baker, Rebecca J; Busch, Ingrid Karin; Hilliard, Michael R; Middleton, Richard S; Schultze, Michael
2009-01-01
In this paper, we describe our experience with solving difficult mixed-integer linear programming problems (MILPs) on the petaflop Cray XT5 system at the National Center for Computational Sciences at Oak Ridge National Laboratory. We describe the algorithmic, software, and hardware needs for solving MILPs and present the results of using PICO, an open-source, parallel, mixed-integer linear programming solver developed at Sandia National Laboratories, to solve canonical MILPs as well as problems of interest arising from the logistics and supply chain management field.
NASA Astrophysics Data System (ADS)
Guo, P.; Huang, G. H.; Li, Y. P.
2010-01-01
In this study, an inexact fuzzy-chance-constrained two-stage mixed-integer linear programming (IFCTIP) approach is developed for flood diversion planning under multiple uncertainties. A concept of the distribution with fuzzy boundary interval probability is defined to address multiple uncertainties expressed as integration of intervals, fuzzy sets and probability distributions. IFCTIP integrates the inexact programming, two-stage stochastic programming, integer programming and fuzzy-stochastic programming within a general optimization framework. IFCTIP incorporates the pre-regulated water-diversion policies directly into its optimization process to analyze various policy scenarios; each scenario has different economic penalty when the promised targets are violated. More importantly, it can facilitate dynamic programming for decisions of capacity-expansion planning under fuzzy-stochastic conditions. IFCTIP is applied to a flood management system. Solutions from IFCTIP provide desired flood diversion plans with a minimized system cost and a maximized safety level. The results indicate that reasonable solutions are generated for objective function values and decision variables, thus a number of decision alternatives can be generated under different levels of flood flows.
Guo, P.; Huang, G.H.
2010-03-15
In this study, an interval-parameter semi-infinite fuzzy-chance-constrained mixed-integer linear programming (ISIFCIP) approach is developed for supporting long-term planning of waste-management systems under multiple uncertainties in the City of Regina, Canada. The method improves upon the existing interval-parameter semi-infinite programming (ISIP) and fuzzy-chance-constrained programming (FCCP) by incorporating uncertainties expressed as dual uncertainties of functional intervals and multiple uncertainties of distributions with fuzzy-interval admissible probability of violating constraint within a general optimization framework. The binary-variable solutions represent the decisions of waste-management-facility expansion, and the continuous ones are related to decisions of waste-flow allocation. The interval solutions can help decision-makers to obtain multiple decision alternatives, as well as provide bases for further analyses of tradeoffs between waste-management cost and system-failure risk. In the application to the City of Regina, Canada, two scenarios are considered. In Scenario 1, the City's waste-management practices would be based on the existing policy over the next 25 years. The total diversion rate for the residential waste would be approximately 14%. Scenario 2 is associated with a policy for waste minimization and diversion, where 35% diversion of residential waste should be achieved within 15 years, and 50% diversion over 25 years. In this scenario, not only landfill would be expanded, but also CF and MRF would be expanded. Through the scenario analyses, useful decision support for the City's solid-waste managers and decision-makers has been generated. Three special characteristics of the proposed method make it unique compared with other optimization techniques that deal with uncertainties. Firstly, it is useful for tackling multiple uncertainties expressed as intervals, functional intervals, probability distributions, fuzzy sets, and their combinations; secondly, it has capability in addressing the temporal variations of the functional intervals; thirdly, it can facilitate dynamic analysis for decisions of facility-expansion planning and waste-flow allocation within a multi-facility, multi-period and multi-option context.
Zhang, Huiling; Huang, Qingsheng; Bei, Zhendong; Wei, Yanjie; Floudas, Christodoulos A
2016-03-01
In this article, we present COMSAT, a hybrid framework for residue contact prediction of transmembrane (TM) proteins, integrating a support vector machine (SVM) method and a mixed integer linear programming (MILP) method. COMSAT consists of two modules: COMSAT_SVM which is trained mainly on position-specific scoring matrix features, and COMSAT_MILP which is an ab initio method based on optimization models. Contacts predicted by the SVM model are ranked by SVM confidence scores, and a threshold is trained to improve the reliability of the predicted contacts. For TM proteins with no contacts above the threshold, COMSAT_MILP is used. The proposed hybrid contact prediction scheme was tested on two independent TM protein sets based on the contact definition of 14 Å between Cα-Cα atoms. First, using a rigorous leave-one-protein-out cross validation on the training set of 90 TM proteins, an accuracy of 66.8%, a coverage of 12.3%, a specificity of 99.3% and a Matthews' correlation coefficient (MCC) of 0.184 were obtained for residue pairs that are at least six amino acids apart. Second, when tested on a test set of 87 TM proteins, the proposed method showed a prediction accuracy of 64.5%, a coverage of 5.3%, a specificity of 99.4% and a MCC of 0.106. COMSAT shows satisfactory results when compared with 12 other state-of-the-art predictors, and is more robust in terms of prediction accuracy as the length and complexity of TM protein increase. COMSAT is freely accessible at http://hpcc.siat.ac.cn/COMSAT/. Proteins 2016; 84:332-348. © 2016 Wiley Periodicals, Inc. PMID:26756402
NASA Astrophysics Data System (ADS)
Irmeilyana, Puspita, Fitri Maya; Indrawati
2016-02-01
The pricing for wireless networks is developed by considering linearity factors, elasticity price and price factors. Mixed Integer Nonlinear Programming of wireless pricing model is proposed as the nonlinear programming problem that can be solved optimally using LINGO 13.0. The solutions are expected to give some information about the connections between the acceptance factor and the price. Previous model worked on the model that focuses on bandwidth as the QoS attribute. The models attempt to maximize the total price for a connection based on QoS parameter. The QoS attributes used will be the bandwidth and the end to end delay that affect the traffic. The maximum goal to maximum price is achieved when the provider determine the requirement for the increment or decrement of price change due to QoS change and amount of QoS value.
A mixed integer bi-level DEA model for bank branch performance evaluation by Stackelberg approach
NASA Astrophysics Data System (ADS)
Shafiee, Morteza; Lotfi, Farhad Hosseinzadeh; Saleh, Hilda; Ghaderi, Mehdi
2015-11-01
One of the most complicated decision making problems for managers is the evaluation of bank performance, which involves various criteria. There are many studies about bank efficiency evaluation by network DEA in the literature review. These studies do not focus on multi-level network. Wu (Eur J Oper Res 207:856-864, 2010) proposed a bi-level structure for cost efficiency at the first time. In this model, multi-level programming and cost efficiency were used. He used a nonlinear programming to solve the model. In this paper, we have focused on multi-level structure and proposed a bi-level DEA model. We then used a liner programming to solve our model. In other hand, we significantly improved the way to achieve the optimum solution in comparison with the work by Wu (2010) by converting the NP-hard nonlinear programing into a mixed integer linear programming. This study uses a bi-level programming data envelopment analysis model that embodies internal structure with Stackelberg-game relationships to evaluate the performance of banking chain. The perspective of decentralized decisions is taken in this paper to cope with complex interactions in banking chain. The results derived from bi-level programming DEA can provide valuable insights and detailed information for managers to help them evaluate the performance of the banking chain as a whole using Stackelberg-game relationships. Finally, this model was applied in the Iranian bank to evaluate cost efficiency.
A mixed integer bi-level DEA model for bank branch performance evaluation by Stackelberg approach
NASA Astrophysics Data System (ADS)
Shafiee, Morteza; Lotfi, Farhad Hosseinzadeh; Saleh, Hilda; Ghaderi, Mehdi
2016-11-01
One of the most complicated decision making problems for managers is the evaluation of bank performance, which involves various criteria. There are many studies about bank efficiency evaluation by network DEA in the literature review. These studies do not focus on multi-level network. Wu (Eur J Oper Res 207:856-864, 2010) proposed a bi-level structure for cost efficiency at the first time. In this model, multi-level programming and cost efficiency were used. He used a nonlinear programming to solve the model. In this paper, we have focused on multi-level structure and proposed a bi-level DEA model. We then used a liner programming to solve our model. In other hand, we significantly improved the way to achieve the optimum solution in comparison with the work by Wu (2010) by converting the NP-hard nonlinear programing into a mixed integer linear programming. This study uses a bi-level programming data envelopment analysis model that embodies internal structure with Stackelberg-game relationships to evaluate the performance of banking chain. The perspective of decentralized decisions is taken in this paper to cope with complex interactions in banking chain. The results derived from bi-level programming DEA can provide valuable insights and detailed information for managers to help them evaluate the performance of the banking chain as a whole using Stackelberg-game relationships. Finally, this model was applied in the Iranian bank to evaluate cost efficiency.
Drawing and Labeling High-Quality Metro Maps by Mixed-Integer Programming.
Nollenburg, M; Wolff, A
2011-05-01
Metro maps are schematic diagrams of public transport networks that serve as visual aids for route planning and navigation tasks. It is a challenging problem in network visualization to automatically draw appealing metro maps. There are two aspects to this problem that depend on each other: the layout problem of finding station and link coordinates and the labeling problem of placing nonoverlapping station labels. In this paper, we present a new integral approach that solves the combined layout and labeling problem (each of which, independently, is known to be NP-hard) using mixed-integer programming (MIP). We identify seven design rules used in most real-world metro maps. We split these rules into hard and soft constraints and translate them into an MIP model. Our MIP formulation finds a metro map that satisfies all hard constraints (if such a drawing exists) and minimizes a weighted sum of costs that correspond to the soft constraints. We have implemented the MIP model and present a case study and the results of an expert assessment to evaluate the performance of our approach in comparison to both manually designed official maps and results of previous layout methods. PMID:20498505
Using a mixed integer linear programming technique to optimize a fracture treatment design
Rueda, J.I.; Rahim, Z.; Holditch, S.A.
1994-12-31
The design of the optimum fracture treatment has three basic steps. The first step involves calculating fracture dimensions and conductivity for various fracture fluid pumping schedules. The second step is determining the oil and gas production rates and recoveries using the values of propped fracture length obtained from the fracture treatment design. The third step requires one to determine the optimum fracture treatment design by maximizing the economic benefit of the treatment. Since a fracture treatment design involves selection of fracturing fluids, additives, proppant materials, injection rate, pump schedule, and fracture dimensions, the determination of the optimum combination of all variables can be quite complicated. In this research, the authors have developed two methods to investigate reasonable combinations of design and treatment parameters to determine the most profitable fracture treatment design. Depending upon the post-fracture treatment production rates, the related expenses, and the economic constraints, the optimum treatment can be easily determined.
Mixed-integer programming methods for transportation and power generation problems
NASA Astrophysics Data System (ADS)
Damci Kurt, Pelin
This dissertation conducts theoretical and computational research to solve challenging problems in application areas such as supply chain and power systems. The first part of the dissertation studies a transportation problem with market choice (TPMC) which is a variant of the classical transportation problem in which suppliers with limited capacities have a choice of which demands (markets) to satisfy. We show that TPMC is strongly NP-complete. We consider a version of the problem with a service level constraint on the maximum number of markets that can be rejected and show that if the original problem is polynomial, its cardinality-constrained version is also polynomial. We propose valid inequalities for mixed-integer cover and knapsack sets with variable upper bound constraints, which appear as substructures of TPMC and use them in a branch-and-cut algorithm to solve this problem. The second part of this dissertation studies a unit commitment (UC) problem in which the goal is to minimize the operational cost of power generators over a time period subject to physical constraints while satisfying demand. We provide several exponential classes of multi-period ramping and multi-period variable upper bound inequalities. We prove the strength of these inequalities and describe polynomial-time separation algorithms. Computational results show the effectiveness of the proposed inequalities when used as cuts in a branch-and-cut algorithm to solve the UC problem. The last part of this dissertation investigates the effects of uncertain wind power on the UC problem. A two-stage robust model and a three-stage stochastic program are compared.
NASA Astrophysics Data System (ADS)
Purnomo, Muhammad Ridwan Andi; Satrio Wiwoho, Yoga
2016-01-01
Facility layout becomes one of production system factor that should be managed well, as it is designated for the location of production. In managing the layout, designing the layout by considering the optimal layout condition that supports the work condition is essential. One of the method for facility layout optimization is Mixed Integer Programming (MIP). In this study, the MIP is solved using Lingo 9.0 software and considering quantitative and qualitative objectives to be achieved simultaneously: minimizing material handling cost, maximizing closeness rating, and minimizing re-layout cost. The research took place in Rekayasa Wangdi as a make to order company, focusing on the making of concrete brick dough stirring machine with 10 departments involved. The result shows an improvement in the new layout for 333,72 points of objective value compared with the initial layout. As the conclusion, the proposed MIP is proven to be used to model facility layout problem under multi objective consideration for a more realistic look.
Zou, Meng; Zhang, Peng-Jun; Wen, Xin-Yu; Chen, Luonan; Tian, Ya-Ping; Wang, Yong
2015-07-15
Multi-biomarker panels can capture the nonlinear synergy among biomarkers and they are important to aid in the early diagnosis and ultimately battle complex diseases. However, identification of these multi-biomarker panels from case and control data is challenging. For example, the exhaustive search method is computationally infeasible when the data dimension is high. Here, we propose a novel method, MILP_k, to identify serum-based multi-biomarker panel to distinguish colorectal cancers (CRC) from benign colorectal tumors. Specifically, the multi-biomarker panel detection problem is modeled by a mixed integer programming to maximize the classification accuracy. Then we measured the serum profiling data for 101 CRC patients and 95 benign patients. The 61 biomarkers were analyzed individually and further their combinations by our method. We discovered 4 biomarkers as the optimal small multi-biomarker panel, including known CRC biomarkers CEA and IL-10 as well as novel biomarkers IMA and NSE. This multi-biomarker panel obtains leave-one-out cross-validation (LOOCV) accuracy to 0.7857 by nearest centroid classifier. An independent test of this panel by support vector machine (SVM) with threefold cross validation gets an AUC 0.8438. This greatly improves the predictive accuracy by 20% over the single best biomarker. Further extension of this 4-biomarker panel to a larger 13-biomarker panel improves the LOOCV to 0.8673 with independent AUC 0.8437. Comparison with the exhaustive search method shows that our method dramatically reduces the searching time by 1000-fold. Experiments on the early cancer stage samples reveal two panel of biomarkers and show promising accuracy. The proposed method allows us to select the subset of biomarkers with best accuracy to distinguish case and control samples given the number of selected biomarkers. Both receiver operating characteristic curve and precision-recall curve show our method's consistent performance gain in accuracy. Our method also shows its advantage in capturing synergy among selected biomarkers. The multi-biomarker panel far outperforms the simple combination of best single features. Close investigation of the multi-biomarker panel illustrates that our method possesses the ability to remove redundancy and reveals complementary biomarker combinations. In addition, our method is efficient and can select multi-biomarker panel with more than 5 biomarkers, for which the exhaustive methods fail. In conclusion, we propose a promising model to improve the clinical data interpretability and to serve as a useful tool for other complex disease studies. Our small multi-biomarker panel, CEA, IL-10, IMA, and NSE, may provide insights on the disease status of colorectal diseases. The implementation of our method in MATLAB is available via the website: http://doc.aporc.org/wiki/MILP_k. PMID:25980368
Final Report---Optimization Under Nonconvexity and Uncertainty: Algorithms and Software
Jeff Linderoth
2011-11-06
the goal of this work was to develop new algorithmic techniques for solving large-scale numerical optimization problems, focusing on problems classes that have proven to be among the most challenging for practitioners: those involving uncertainty and those involving nonconvexity. This research advanced the state-of-the-art in solving mixed integer linear programs containing symmetry, mixed integer nonlinear programs, and stochastic optimization problems. The focus of the work done in the continuation was on Mixed Integer Nonlinear Programs (MINLP)s and Mixed Integer Linear Programs (MILP)s, especially those containing a great deal of symmetry.
Zhang, Xiaodong; Huang, Gordon
2013-02-15
Greenhouse gas (GHG) emissions from municipal solid waste (MSW) management facilities have become a serious environmental issue. In MSW management, not only economic objectives but also environmental objectives should be considered simultaneously. In this study, a dynamic stochastic possibilistic multiobjective programming (DSPMP) model is developed for supporting MSW management and associated GHG emission control. The DSPMP model improves upon the existing waste management optimization methods through incorporation of fuzzy possibilistic programming and chance-constrained programming into a general mixed-integer multiobjective linear programming (MOP) framework where various uncertainties expressed as fuzzy possibility distributions and probability distributions can be effectively reflected. Two conflicting objectives are integrally considered, including minimization of total system cost and minimization of total GHG emissions from waste management facilities. Three planning scenarios are analyzed and compared, representing different preferences of the decision makers for economic development and environmental-impact (i.e. GHG-emission) issues in integrated MSW management. Optimal decision schemes under three scenarios and different p(i) levels (representing the probability that the constraints would be violated) are generated for planning waste flow allocation and facility capacity expansions as well as GHG emission control. The results indicate that economic and environmental tradeoffs can be effectively reflected through the proposed DSPMP model. The generated decision variables can help the decision makers justify and/or adjust their waste management strategies based on their implicit knowledge and preferences. PMID:23313898
NASA Astrophysics Data System (ADS)
Vafaeinezhad, Moghadaseh; Kia, Reza; Shahnazari-Shahrezaei, Parisa
2015-11-01
Cell formation (CF) problem is one of the most important decision problems in designing a cellular manufacturing system includes grouping machines into machine cells and parts into part families. Several factors should be considered in a cell formation problem. In this work, robust optimization of a mathematical model of a dynamic cell formation problem integrating CF, production planning and worker assignment is implemented with uncertain scenario-based data. The robust approach is used to reduce the effects of fluctuations of the uncertain parameters with regards to all possible future scenarios. In this research, miscellaneous cost parameters of the cell formation and demand fluctuations are subject to uncertainty and a mixed-integer nonlinear programming model is developed to formulate the related robust dynamic cell formation problem. The objective function seeks to minimize total costs including machine constant, machine procurement, machine relocation, machine operation, inter-cell and intra-cell movement, overtime, shifting labors between cells and inventory holding. Finally, a case study is carried out to display the robustness and effectiveness of the proposed model. The tradeoff between solution robustness and model robustness is also analyzed in the obtained results.
Optimization using Extremal Dynamics
NASA Astrophysics Data System (ADS)
Boettcher, Stefan
2001-03-01
We explore a new heuristic for finding high-quality solutions to NP-hard optimization problems which we have recently introduced [see ``Nature's Way of Optimizing," Artificial Intelligence 119, 275-286 (2000) and cond-mat/0010337]. The method, called extremal optimization, is inspired by self-organized criticality, a concept introduced to describe emergent complexity in physical systems. Extremal optimization successively replaces extremely undesirable elements of a single sub-optimal solution with new, random ones. Large fluctuations ensue that efficiently explore many local optima. With only one adjustable parameter, its performance has proved competitive with more elaborate methods, especially near phase transitions which are believed to contain the hardest instances. In particular, extremal optimization is superior to simulated annealing in the partitioning of sparse graphs, it finds the overlap of all ground-states at the phase transition of the 3-coloring problem, and it provides independent confirmation for the ground-state energy of spin glasses, previously obtained with elaborate genetic algorithms.
Adaptive critics for dynamic optimization.
Kulkarni, Raghavendra V; Venayagamoorthy, Ganesh Kumar
2010-06-01
A novel action-dependent adaptive critic design (ACD) is developed for dynamic optimization. The proposed combination of a particle swarm optimization-based actor and a neural network critic is demonstrated through dynamic sleep scheduling of wireless sensor motes for wildlife monitoring. The objective of the sleep scheduler is to dynamically adapt the sleep duration to node's battery capacity and movement pattern of animals in its environment in order to obtain snapshots of the animal on its trajectory uniformly. Simulation results show that the sleep time of the node determined by the actor critic yields superior quality of sensory data acquisition and enhanced node longevity. PMID:20223635
Lu, H W; Huang, G H; He, L; Zeng, G M
2009-01-01
Municipal solid waste (MSW) should be properly disposed in order to help protect environmental quality and human health, as well as to preserve natural resources. During MSW disposal processes, a large amount of greenhouse gas (GHG) is emitted, leading to a significant impact on climate change. In this study, an inexact dynamic optimization model (IDOM) is developed for MSW-management systems under uncertainty. It grounds upon conventional mixed-integer linear programming (MILP) approaches, and integrates GHG components into the modeling framework. Compared with the existing models, IDOM can not only deal with the complex tradeoff between system cost minimization and GHG-emission mitigation, but also provide optimal allocation strategies under various emission-control standards. A case study is then provided for demonstrating applicability of the developed model. The results indicate that desired waste-flow patterns with a minimized system cost and GHG-emission amount can be obtained. Of more importance, the IDOM solution is associated with over 5.5 million tonnes of TEC reduction, which is of significant economic implication for real implementations. Therefore, the proposed model could be regarded as a useful tool for realizing comprehensive MSW management with regard to mitigating climate-change impacts. PMID:18096299
Optimal operation of a potable water distribution network.
Biscos, C; Mulholland, M; Le Lann, M V; Brouckaert, C J; Bailey, R; Roustan, M
2002-01-01
This paper presents an approach to an optimal operation of a potable water distribution network. The main control objective defined during the preliminary steps was to maximise the use of low-cost power, maintaining at the same time minimum emergency levels in all reservoirs. The combination of dynamic elements (e.g. reservoirs) and discrete elements (pumps, valves, routing) makes this a challenging predictive control and constrained optimisation problem, which is being solved by MINLP (Mixed Integer Non-linear Programming). Initial experimental results show the performance of this algorithm and its ability to control the water distribution process. PMID:12448464
Optimal dynamic detection of explosives
Moore, David Steven; Mcgrane, Shawn D; Greenfield, Margo T; Scharff, R J; Rabitz, Herschel A; Roslund, J
2009-01-01
The detection of explosives is a notoriously difficult problem, especially at stand-off distances, due to their (generally) low vapor pressure, environmental and matrix interferences, and packaging. We are exploring optimal dynamic detection to exploit the best capabilities of recent advances in laser technology and recent discoveries in optimal shaping of laser pulses for control of molecular processes to significantly enhance the standoff detection of explosives. The core of the ODD-Ex technique is the introduction of optimally shaped laser pulses to simultaneously enhance sensitivity of explosives signatures while reducing the influence of noise and the signals from background interferents in the field (increase selectivity). These goals are being addressed by operating in an optimal nonlinear fashion, typically with a single shaped laser pulse inherently containing within it coherently locked control and probe sub-pulses. With sufficient bandwidth, the technique is capable of intrinsically providing orthogonal broad spectral information for data fusion, all from a single optimal pulse.
New numerical methods for open-loop and feedback solutions to dynamic optimization problems
NASA Astrophysics Data System (ADS)
Ghosh, Pradipto
The topic of the first part of this research is trajectory optimization of dynamical systems via computational swarm intelligence. Particle swarm optimization is a nature-inspired heuristic search method that relies on a group of potential solutions to explore the fitness landscape. Conceptually, each particle in the swarm uses its own memory as well as the knowledge accumulated by the entire swarm to iteratively converge on an optimal or near-optimal solution. It is relatively straightforward to implement and unlike gradient-based solvers, does not require an initial guess or continuity in the problem definition. Although particle swarm optimization has been successfully employed in solving static optimization problems, its application in dynamic optimization, as posed in optimal control theory, is still relatively new. In the first half of this thesis particle swarm optimization is used to generate near-optimal solutions to several nontrivial trajectory optimization problems including thrust programming for minimum fuel, multi-burn spacecraft orbit transfer, and computing minimum-time rest-to-rest trajectories for a robotic manipulator. A distinct feature of the particle swarm optimization implementation in this work is the runtime selection of the optimal solution structure. Optimal trajectories are generated by solving instances of constrained nonlinear mixed-integer programming problems with the swarming technique. For each solved optimal programming problem, the particle swarm optimization result is compared with a nearly exact solution found via a direct method using nonlinear programming. Numerical experiments indicate that swarm search can locate solutions to very great accuracy. The second half of this research develops a new extremal-field approach for synthesizing nearly optimal feedback controllers for optimal control and two-player pursuit-evasion games described by general nonlinear differential equations. A notable revelation from this development is that the resulting control law has an algebraic closed-form structure. The proposed method uses an optimal spatial statistical predictor called universal kriging to construct the surrogate model of a feedback controller, which is capable of quickly predicting an optimal control estimate based on current state (and time) information. With universal kriging, an approximation to the optimal feedback map is computed by conceptualizing a set of state-control samples from pre-computed extremals to be a particular realization of a jointly Gaussian spatial process. Feedback policies are computed for a variety of example dynamic optimization problems in order to evaluate the effectiveness of this methodology. This feedback synthesis approach is found to combine good numerical accuracy with low computational overhead, making it a suitable candidate for real-time applications. Particle swarm and universal kriging are combined for a capstone example, a near optimal, near-admissible, full-state feedback control law is computed and tested for the heat-load-limited atmospheric-turn guidance of an aeroassisted transfer vehicle. The performance of this explicit guidance scheme is found to be very promising; initial errors in atmospheric entry due to simulated thruster misfirings are found to be accurately corrected while closely respecting the algebraic state-inequality constraint.
An optimal spacecraft scheduling model for the NASA deep space network
NASA Technical Reports Server (NTRS)
Webb, W. A.
1985-01-01
A computer model is described which uses mixed-integer linear programming to provide optimal DSN spacecraft schedules given a mission set and specified scheduling requirements. A solution technique is proposed which uses Bender's method and a heuristic starting algorithm.
Final Report-Optimization Under Uncertainty and Nonconvexity: Algorithms and Software
Jeff Linderoth
2008-10-10
The goal of this research was to develop new algorithmic techniques for solving large-scale numerical optimization problems, focusing on problems classes that have proven to be among the most challenging for practitioners: those involving uncertainty and those involving nonconvexity. This research advanced the state-of-the-art in solving mixed integer linear programs containing symmetry, mixed integer nonlinear programs, and stochastic optimization problems.
Multiobjective optimization of dynamic aperture
Yang, L.; Li, Y.; Guo, W.; Krinsky, S.
2011-05-02
Dynamic aperture (DA) is one of the key nonlinear properties for a storage ring. Although there have been both analytical and numerical methods to find the aperture, the reverse problem of how to optimize it is still a challenging problem. A general and flexible way of optimizing the DA is highly demanded in accelerator design and operation. In this paper, we discuss the use of multiobjective optimization for DA. First we consider using objective functions based only on numerical tracking results. Data mining of these results demonstrated a correlation between DA and low-order nonlinear driving terms. Next we considered using objective functions which included both numerical tracking results and analytical estimates of low-order nonlinear driving terms. This resulted in faster convergence. The National Synchrotron Light Source II (NSLS-II) lattice was taken as an example to illustrate this method. This multiobjective approach is not limited by particular linear or nonlinear lattice settings, and can also be applied for optimizing other properties of a storage ring.
Demystifying optimal dynamic treatment regimes.
Moodie, Erica E M; Richardson, Thomas S; Stephens, David A
2007-06-01
A dynamic regime is a function that takes treatment and covariate history and baseline covariates as inputs and returns a decision to be made. Murphy (2003, Journal of the Royal Statistical Society, Series B 65, 331-366) and Robins (2004, Proceedings of the Second Seattle Symposium on Biostatistics, 189-326) have proposed models and developed semiparametric methods for making inference about the optimal regime in a multi-interval trial that provide clear advantages over traditional parametric approaches. We show that Murphy's model is a special case of Robins's and that the methods are closely related but not equivalent. Interesting features of the methods are highlighted using the Multicenter AIDS Cohort Study and through simulation. PMID:17688497
TRACKING CODE DEVELOPMENT FOR BEAM DYNAMICS OPTIMIZATION
Yang, L.
2011-03-28
Dynamic aperture (DA) optimization with direct particle tracking is a straight forward approach when the computing power is permitted. It can have various realistic errors included and is more close than theoretical estimations. In this approach, a fast and parallel tracking code could be very helpful. In this presentation, we describe an implementation of storage ring particle tracking code TESLA for beam dynamics optimization. It supports MPI based parallel computing and is robust as DA calculation engine. This code has been used in the NSLS-II dynamics optimizations and obtained promising performance.
Baliban, Richard C.; DiMaggio, Peter A.; Plazas-Mayorca, Mariana D.; Garcia, Benjamin A.; Floudas, Christodoulos A.
2012-01-01
A novel protein identification framework, PILOT_PROTEIN, has been developed to construct a comprehensive list of all unmodified proteins that are present in a living sample. It uses the peptide identification results from the PILOT_SEQUEL algorithm to initially determine all unmodified proteins within the sample. Using a rigorous biclustering approach that groups incorrect peptide sequences with other homologous sequences, the number of false positives reported is minimized. A sequence tag procedure is then incorporated along with the untargeted PTM identification algorithm, PILOT_PTM, to determine a list of all modification types and sites for each protein. The unmodified protein identification algorithm, PILOT_PROTEIN, is compared to the methods SEQUEST, InsPecT, X!Tandem, VEMS, and ProteinProspector using both prepared protein samples and a more complex chromatin digest. The algorithm demonstrates superior protein identification accuracy with a lower false positive rate. All materials are freely available to the scientific community at http://pumpd.princeton.edu. PMID:22788846
Streak camera dynamic range optimization
Wiedwald, J.D.; Lerche, R.A.
1987-09-01
The LLNL optical streak camera is used by the Laser Fusion Program in a wide range of applications. Many of these applications require a large recorded dynamic range. Recent work has focused on maximizing the dynamic range of the streak camera recording system. For our streak cameras, image intensifier saturation limits the upper end of the dynamic range. We have developed procedures to set the image intensifier gain such that the system dynamic range is maximized. Specifically, the gain is set such that a single streak tube photoelectron is recorded with an exposure of about five times the recording system noise. This ensures detection of single photoelectrons, while not consuming intensifier or recording system dynamic range through excessive intensifier gain. The optimum intensifier gain has been determined for two types of film and for a lens-coupled CCD camera. We have determined that by recording the streak camera image with a CCD camera, the system is shot-noise limited up to the onset of image intensifier nonlinearity. When recording on film, the film determines the noise at high exposure levels. There is discussion of the effects of slit width and image intensifier saturation on dynamic range. 8 refs.
Optimizing Dynamical Network Structure for Pinning Control.
Orouskhani, Yasin; Jalili, Mahdi; Yu, Xinghuo
2016-01-01
Controlling dynamics of a network from any initial state to a final desired state has many applications in different disciplines from engineering to biology and social sciences. In this work, we optimize the network structure for pinning control. The problem is formulated as four optimization tasks: i) optimizing the locations of driver nodes, ii) optimizing the feedback gains, iii) optimizing simultaneously the locations of driver nodes and feedback gains, and iv) optimizing the connection weights. A newly developed population-based optimization technique (cat swarm optimization) is used as the optimization method. In order to verify the methods, we use both real-world networks, and model scale-free and small-world networks. Extensive simulation results show that the optimal placement of driver nodes significantly outperforms heuristic methods including placing drivers based on various centrality measures (degree, betweenness, closeness and clustering coefficient). The pinning controllability is further improved by optimizing the feedback gains. We also show that one can significantly improve the controllability by optimizing the connection weights. PMID:27067020
Optimizing Dynamical Network Structure for Pinning Control
NASA Astrophysics Data System (ADS)
Orouskhani, Yasin; Jalili, Mahdi; Yu, Xinghuo
2016-04-01
Controlling dynamics of a network from any initial state to a final desired state has many applications in different disciplines from engineering to biology and social sciences. In this work, we optimize the network structure for pinning control. The problem is formulated as four optimization tasks: i) optimizing the locations of driver nodes, ii) optimizing the feedback gains, iii) optimizing simultaneously the locations of driver nodes and feedback gains, and iv) optimizing the connection weights. A newly developed population-based optimization technique (cat swarm optimization) is used as the optimization method. In order to verify the methods, we use both real-world networks, and model scale-free and small-world networks. Extensive simulation results show that the optimal placement of driver nodes significantly outperforms heuristic methods including placing drivers based on various centrality measures (degree, betweenness, closeness and clustering coefficient). The pinning controllability is further improved by optimizing the feedback gains. We also show that one can significantly improve the controllability by optimizing the connection weights.
Optimizing Dynamical Network Structure for Pinning Control
Orouskhani, Yasin; Jalili, Mahdi; Yu, Xinghuo
2016-01-01
Controlling dynamics of a network from any initial state to a final desired state has many applications in different disciplines from engineering to biology and social sciences. In this work, we optimize the network structure for pinning control. The problem is formulated as four optimization tasks: i) optimizing the locations of driver nodes, ii) optimizing the feedback gains, iii) optimizing simultaneously the locations of driver nodes and feedback gains, and iv) optimizing the connection weights. A newly developed population-based optimization technique (cat swarm optimization) is used as the optimization method. In order to verify the methods, we use both real-world networks, and model scale-free and small-world networks. Extensive simulation results show that the optimal placement of driver nodes significantly outperforms heuristic methods including placing drivers based on various centrality measures (degree, betweenness, closeness and clustering coefficient). The pinning controllability is further improved by optimizing the feedback gains. We also show that one can significantly improve the controllability by optimizing the connection weights. PMID:27067020
Dynamic optimization and adaptive controller design
NASA Astrophysics Data System (ADS)
Inamdar, S. R.
2010-10-01
In this work I present a new type of controller which is an adaptive tracking controller which employs dynamic optimization for optimizing current value of controller action for the temperature control of nonisothermal continuously stirred tank reactor (CSTR). We begin with a two-state model of nonisothermal CSTR which are mass and heat balance equations and then add cooling system dynamics to eliminate input multiplicity. The initial design value is obtained using local stability of steady states where approach temperature for cooling action is specified as a steady state and a design specification. Later we make a correction in the dynamics where material balance is manipulated to use feed concentration as a system parameter as an adaptive control measure in order to avoid actuator saturation for the main control loop. The analysis leading to design of dynamic optimization based parameter adaptive controller is presented. The important component of this mathematical framework is reference trajectory generation to form an adaptive control measure.
Efficient dynamic optimization of logic programs
NASA Technical Reports Server (NTRS)
Laird, Phil
1992-01-01
A summary is given of the dynamic optimization approach to speed up learning for logic programs. The problem is to restructure a recursive program into an equivalent program whose expected performance is optimal for an unknown but fixed population of problem instances. We define the term 'optimal' relative to the source of input instances and sketch an algorithm that can come within a logarithmic factor of optimal with high probability. Finally, we show that finding high-utility unfolding operations (such as EBG) can be reduced to clause reordering.
Optimizing dissolution dynamic nuclear polarization.
Bornet, Aurélien; Jannin, Sami
2016-03-01
This article is a short review of some of our recent developments in dissolution dynamic nuclear polarization (d-DNP). We present the basic principles of d-DNP, and motivate our choice to step away from conventional approaches. We then introduce a modified d-DNP recipe that can be summed up as follows. PMID:26920826
Optimizing dissolution dynamic nuclear polarization
NASA Astrophysics Data System (ADS)
Bornet, Aurélien; Jannin, Sami
2016-03-01
This article is a short review of some of our recent developments in dissolution dynamic nuclear polarization (d-DNP). We present the basic principles of d-DNP, and motivate our choice to step away from conventional approaches. We then introduce a modified d-DNP recipe that can be summed up as follows:
Dynamic programming in applied optimization problems
NASA Astrophysics Data System (ADS)
Zavalishchin, Dmitry
2015-11-01
Features of the use dynamic programming in applied problems are investigated. In practice such problems as finding the critical paths in network planning and control, finding the optimal supply plan in transportation problem, objects territorial distribution are traditionally solved by special methods of operations research. It should be noted that the dynamic programming is not provided computational advantages, but facilitates changes and modifications of tasks. This follows from the Bellman's optimality principle. The features of the multistage decision processes construction in applied problems are provided.
Optimal BLS: Optimizing transit-signal detection for Keplerian dynamics
NASA Astrophysics Data System (ADS)
Ofir, Aviv
2015-08-01
Transit surveys, both ground- and space-based, have already accumulated a large number of light curves that span several years. We optimize the search for transit signals for both detection and computational efficiencies by assuming that the searched systems can be described by Keplerian, and propagating the effects of different system parameters to the detection parameters. Importnantly, we mainly consider the information content of the transit signal and not any specific algorithm - and use BLS (Kovács, Zucker, & Mazeh 2002) just as a specific example.We show that the frequency information content of the light curve is primarily determined by the duty cycle of the transit signal, and thus the optimal frequency sampling is found to be cubic and not linear. Further optimization is achieved by considering duty-cycle dependent binning of the phased light curve. By using the (standard) BLS, one is either fairly insensitive to long-period planets or less sensitive to short-period planets and computationally slower by a significant factor of ~330 (for a 3 yr long dataset). We also show how the physical system parameters, such as the host star's size and mass, directly affect transit detection. This understanding can then be used to optimize the search for every star individually.By considering Keplerian dynamics explicitly rather than implicitly one can optimally search the transit signal parameter space. The presented Optimal BLS enhances the detectability of both very short and very long period planets, while allowing such searches to be done with much reduced resources and time. The Matlab/Octave source code for Optimal BLS is made available.
Estimation of Optimal Dynamic Treatment Regimes
Laber, Eric B.
2014-01-01
Background Recent advances in medical research suggest that the optimal treatment rules should be adaptive to patients over time. This has led to an increasing interest in studying dynamic treatment regimes (DTRs), a sequence of individualized treatment rules, one per stage of clinical intervention, which map present patient information to a recommended treatment. There has been a recent surge of statistical work for estimating optimal DTRs from randomized and observational studies. The purpose of this paper is to review recent methodological progress and applied issues associated with estimating optimal DTRs. Methods We discuss Sequential Multiple Assignment Randomized Trials (SMARTs), a clinical trial design used to study treatment sequences. We use a common estimator of an optimal DTR that applies to SMART data as a platform to discuss several practical and methodological issues. Results We provide a limited survey of practical issues associated with modeling SMART data. We review some existing estimators of optimal dynamic treatment regimes and discuss practical issues associated with these methods including: model building; missing data; statistical inference; and choosing an outcome when only non-responders are re-randomized. We mainly focus on the estimation and inference of DTRs using SMART data. DTRs can also be constructed from observational data, which may be easier to obtain in practice, however, care must be taken to account for potential confounding. PMID:24872361
Optimizing Motion Planning for Hyper Dynamic Manipulator
NASA Astrophysics Data System (ADS)
Aboura, Souhila; Omari, Abdelhafid; Meguenni, Kadda Zemalache
2012-01-01
This paper investigates the optimal motion planning for an hyper dynamic manipulator. As case study, we consider a golf swing robot which is consisting with two actuated joint and a mechanical stoppers. Genetic Algorithm (GA) technique is proposed to solve the optimal golf swing motion which is generated by Fourier series approximation. The objective function for GA approach is to minimizing the intermediate and final state, minimizing the robot's energy consummation and maximizing the robot's speed. Obtained simulation results show the effectiveness of the proposed scheme.
Pareto optimization in algebraic dynamic programming.
Saule, Cédric; Giegerich, Robert
2015-01-01
Pareto optimization combines independent objectives by computing the Pareto front of its search space, defined as the set of all solutions for which no other candidate solution scores better under all objectives. This gives, in a precise sense, better information than an artificial amalgamation of different scores into a single objective, but is more costly to compute. Pareto optimization naturally occurs with genetic algorithms, albeit in a heuristic fashion. Non-heuristic Pareto optimization so far has been used only with a few applications in bioinformatics. We study exact Pareto optimization for two objectives in a dynamic programming framework. We define a binary Pareto product operator [Formula: see text] on arbitrary scoring schemes. Independent of a particular algorithm, we prove that for two scoring schemes A and B used in dynamic programming, the scoring scheme [Formula: see text] correctly performs Pareto optimization over the same search space. We study different implementations of the Pareto operator with respect to their asymptotic and empirical efficiency. Without artificial amalgamation of objectives, and with no heuristics involved, Pareto optimization is faster than computing the same number of answers separately for each objective. For RNA structure prediction under the minimum free energy versus the maximum expected accuracy model, we show that the empirical size of the Pareto front remains within reasonable bounds. Pareto optimization lends itself to the comparative investigation of the behavior of two alternative scoring schemes for the same purpose. For the above scoring schemes, we observe that the Pareto front can be seen as a composition of a few macrostates, each consisting of several microstates that differ in the same limited way. We also study the relationship between abstract shape analysis and the Pareto front, and find that they extract information of a different nature from the folding space and can be meaningfully combined. PMID:26150892
Application of optimal prediction to molecular dynamics
Barber IV, John Letherman
2004-12-01
Optimal prediction is a general system reduction technique for large sets of differential equations. In this method, which was devised by Chorin, Hald, Kast, Kupferman, and Levy, a projection operator formalism is used to construct a smaller system of equations governing the dynamics of a subset of the original degrees of freedom. This reduced system consists of an effective Hamiltonian dynamics, augmented by an integral memory term and a random noise term. Molecular dynamics is a method for simulating large systems of interacting fluid particles. In this thesis, I construct a formalism for applying optimal prediction to molecular dynamics, producing reduced systems from which the properties of the original system can be recovered. These reduced systems require significantly less computational time than the original system. I initially consider first-order optimal prediction, in which the memory and noise terms are neglected. I construct a pair approximation to the renormalized potential, and ignore three-particle and higher interactions. This produces a reduced system that correctly reproduces static properties of the original system, such as energy and pressure, at low-to-moderate densities. However, it fails to capture dynamical quantities, such as autocorrelation functions. I next derive a short-memory approximation, in which the memory term is represented as a linear frictional force with configuration-dependent coefficients. This allows the use of a Fokker-Planck equation to show that, in this regime, the noise is {delta}-correlated in time. This linear friction model reproduces not only the static properties of the original system, but also the autocorrelation functions of dynamical variables.
Utilizing parallel optimization in computational fluid dynamics
NASA Astrophysics Data System (ADS)
Kokkolaras, Michael
1998-12-01
General problems of interest in computational fluid dynamics are investigated by means of optimization. Specifically, in the first part of the dissertation, a method of optimal incremental function approximation is developed for the adaptive solution of differential equations. Various concepts and ideas utilized by numerical techniques employed in computational mechanics and artificial neural networks (e.g. function approximation and error minimization, variational principles and weighted residuals, and adaptive grid optimization) are combined to formulate the proposed method. The basis functions and associated coefficients of a series expansion, representing the solution, are optimally selected by a parallel direct search technique at each step of the algorithm according to appropriate criteria; the solution is built sequentially. In this manner, the proposed method is adaptive in nature, although a grid is neither built nor adapted in the traditional sense using a-posteriori error estimates. Variational principles are utilized for the definition of the objective function to be extremized in the associated optimization problems, ensuring that the problem is well-posed. Complicated data structures and expensive remeshing algorithms and systems solvers are avoided. Computational efficiency is increased by using low-order basis functions and concurrent computing. Numerical results and convergence rates are reported for a range of steady-state problems, including linear and nonlinear differential equations associated with general boundary conditions, and illustrate the potential of the proposed method. Fluid dynamics applications are emphasized. Conclusions are drawn by discussing the method's limitations, advantages, and possible extensions. The second part of the dissertation is concerned with the optimization of the viscous-inviscid-interaction (VII) mechanism in an airfoil flow analysis code. The VII mechanism is based on the concept of a transpiration velocity boundary condition, whose convergence to steady state is accelerated. The number of variables in the associated optimization problem is reduced by means of function approximation concepts to ensure high number of parallel processors to number of necessary function evaluations ratio. Numerical results are presented for the NACA-0012 and the supercritical RAE-2822 airfoils subject to transonic flow conditions using a parallel direct search technique. They exhibit a satisfactory level of accuracy. Speed-up depends on the number of available computational units and increases for more challenging flow conditions and airfoil geometries. The enhanced code constitutes a useful tool for airfoil flow analysis and design and an acceptable alternative to computationally expensive high fidelity codes.
Optimal bolt preload for dynamic loading
NASA Astrophysics Data System (ADS)
Duffey, T. A.
1992-12-01
A simple spring-mass model is developed for closure bolting systems. The model takes into account the effects of bolt prestress. An analytical solution is developed for the case of an initially peaked, exponentially decaying internal pressure pulse acting on the closure. The dependence of peak bolt stresses and deflections on bolt prestress level is investigated and an optimal prestress that minimizes peak bolt stress is found in certain cases. Vulnerability curves are developed for bolted-closure systems to provide rapid evaluation of the dynamic capacity of designs for a range in bolt prestress.
Dynamic optimization theory with multiple objectives
NASA Technical Reports Server (NTRS)
Jones, John, Jr.
1990-01-01
Let V(t) be a vector-valued function for t belonging to closed interval a,b open interval, a real interval. The main purpose of this paper is to establish the existence of a closed interval alpha,beta contained in closed interval a,b for which there exists a t(sub O) belonging to closed interval alpha,beta contained in closed interval a,b such that V(t(sub O)) = 0, the zero vector. Use of such information in the dynamic optimization theory with multiple objectives present is needed. Examples of such systems will be given.
Robust optimization with transiently chaotic dynamical systems
NASA Astrophysics Data System (ADS)
Sumi, R.; Molnár, B.; Ercsey-Ravasz, M.
2014-05-01
Efficiently solving hard optimization problems has been a strong motivation for progress in analog computing. In a recent study we presented a continuous-time dynamical system for solving the NP-complete Boolean satisfiability (SAT) problem, with a one-to-one correspondence between its stable attractors and the SAT solutions. While physical implementations could offer great efficiency, the transiently chaotic dynamics raises the question of operability in the presence of noise, unavoidable on analog devices. Here we show that the probability of finding solutions is robust to noise intensities well above those present on real hardware. We also developed a cellular neural network model realizable with analog circuits, which tolerates even larger noise intensities. These methods represent an opportunity for robust and efficient physical implementations.
Optimal dynamic remapping of data parallel computations
NASA Technical Reports Server (NTRS)
Nicol, David M.; Reynolds, Paul F., Jr.
1990-01-01
A large class of data parallel computations is characterized by a sequence of phases, with phase changes occurring unpredictably. Dynamic remapping of the workload to processors may be required to maintain good performance. The problem considered, for which the utility of remapping and the future behavior of the workload are uncertain, arises when phases exhibit stable execution requirements during a given phase, but requirements change radically between phases. For these situations, a workload assignment generated for one phase may hinder performance during the next phase. This problem is treated formally for a probabilistic model of computation with at most two phases. The authors address the fundamental problem of balancing the expected remapping performance gain against the delay cost, and they derive the optimal remapping decision policy. The promise of the approach is shown by application to multiprocessor implementations of an adaptive gridding fluid dynamics program and to a battlefield simulation program.
Optimal Empirical Prognostic Models of Climate Dynamics
NASA Astrophysics Data System (ADS)
Loskutov, E. M.; Mukhin, D.; Gavrilov, A.; Feigin, A. M.
2014-12-01
In this report the empirical methodology for prediction of climate dynamics is suggested. We construct the dynamical models of data patterns connected with climate indices, from observed spatially distributed time series. The models are based on artificial neural network (ANN) parameterization and have a form of discrete stochastic evolution operator mapping some sequence of systems state on the next one [1]. Different approaches to reconstruction of empirical basis (phase variables) for system's phase space representation, which is appropriate for forecasting the climate index of interest, are discussed in the report; for this purpose both linear and non-linear data expansions are considered. The most important point of the methodology is finding the optimal structural parameters of the model such as dimension of variable vector, i.e. number of principal components used for modeling, the time lag used for prediction, and number of neurons in ANN determining the quality of approximation. Actually, we need to solve the model selection problem, i.e. we want to obtain a model of optimal complexity in relation to analyzed time series. We use MDL approach [2] for this purpose: the model providing best data compression is chosen. The method is applied to space-distributed time-series of sea surface temperature and sea level pressure taken from IRI datasets [3]: the ability of proposed models to predict different climate indices (incl. Multivariate ENSO index, Pacific Decadal Oscillation index, North-Atlantic Oscillation index) is investigated. References:1. Molkov Ya. I., E. M. Loskutov, D. N. Mukhin, and A. M. Feigin, Random dynamical models from time series. Phys. Rev. E, 85, 036216, 2012.2. Molkov, Ya.I., D.N. Mukhin, E.M. Loskutov, A.M. Feigin, and G.A. Fidelin, Using the minimum description length principle for global reconstruction of dynamic systems from noisy time series. Phys. Rev. E, 80, 046207, 2009.3. IRI/LDEO Climate Data Library (http://iridl.ldeo.columbia.edu/)
Optimally designed fields for controlling molecular dynamics
NASA Astrophysics Data System (ADS)
Rabitz, Herschel
1991-10-01
This research concerns the development of molecular control theory techniques for designing optical fields capable of manipulating molecular dynamic phenomena. Although is has been long recognized that lasers should be capable of manipulating dynamic events, many frustrating years of intuitively driven laboratory studies only serve to illustrate the point that the task is complex and defies intuition. The principal new component in the present research is the recognition that this problem falls into the category of control theory and its inherent complexities require the use of modern control theory tools largely developed in the engineering disciplines. Thus, the research has initiated a transfer of the control theory concepts to the molecular scale. Although much contained effort will be needed to fully develop these concepts, the research in this grant set forth the basic components of the theory and carried out illustrative studies involving the design of optical fields capable of controlling rotational, vibrational and electronic degrees of freedom. Optimal control within the quantum mechanical molecular realm represents a frontier area with many possible ultimate applications. At this stage, the theoretical tools need to be joined with merging laboratory optical pulse shaping capabilities to illustrate the power of the concepts.
An Optimal Dynamic Threat Evaluation and Weapon Scheduling Technique
NASA Astrophysics Data System (ADS)
Naeem, H.; Masood, A.
Real time scheduling problems demand high level of flexibility and robustness under complex dynamic scenarios. Threat Evaluation (TE) and Weapon Assignment (WA), together TEWA is one such complex dynamic system having optimal or near optimal utilization of scarce defensive resources of supreme priority. Several static solutions of TEWA have been proposed. This paper discusses an optimal dynamic multi-air threat evaluation and weapon allocation algorithm using a variant of Stable Marriage Algorithm (SMA). WA uses a new dynamic weapon scheduling algorithm, allowing multiple engagements using shoot-look-shoot strategy, to compute near-optimal solution. Testing part of this paper shows feasibility of this approach for a range of scenarios.
Dynamic optimization of a copolymerization reactor using tabu search.
Anand, P; Rao, M Bhagvanth; Venkateswarlu, Ch
2015-03-01
A novel multistage dynamic optimization strategy based on meta-heuristic tabu search (TS) is proposed and evaluated through sequential and simultaneous implementation procedures by applying it to a semi-batch styrene-acrylonitrile (SAN) copolymerization reactor. The adaptive memory and responsive exploration features of TS are exploited to design the dynamic optimization strategy and compute the optimal control policies for temperature and monomer addition rate so as to achieve the desired product quality parameters expressed in terms of single and multiple objectives. The dynamic optimization results of TS sequential and TS simultaneous implementation strategies are analyzed and compared with those of a conventional optimization technique based on iterative dynamic programming (IDP). The simulation results demonstrate the usefulness of TS for optimal control of transient dynamic systems. PMID:25466914
Chaotic dynamics in optimal monetary policy
NASA Astrophysics Data System (ADS)
Gomes, O.; Mendes, V. M.; Mendes, D. A.; Sousa Ramos, J.
2007-05-01
There is by now a large consensus in modern monetary policy. This consensus has been built upon a dynamic general equilibrium model of optimal monetary policy as developed by, e.g., Goodfriend and King [ NBER Macroeconomics Annual 1997 edited by B. Bernanke and J. Rotemberg (Cambridge, Mass.: MIT Press, 1997), pp. 231 282], Clarida et al. [J. Econ. Lit. 37, 1661 (1999)], Svensson [J. Mon. Econ. 43, 607 (1999)] and Woodford [ Interest and Prices: Foundations of a Theory of Monetary Policy (Princeton, New Jersey, Princeton University Press, 2003)]. In this paper we extend the standard optimal monetary policy model by introducing nonlinearity into the Phillips curve. Under the specific form of nonlinearity proposed in our paper (which allows for convexity and concavity and secures closed form solutions), we show that the introduction of a nonlinear Phillips curve into the structure of the standard model in a discrete time and deterministic framework produces radical changes to the major conclusions regarding stability and the efficiency of monetary policy. We emphasize the following main results: (i) instead of a unique fixed point we end up with multiple equilibria; (ii) instead of saddle-path stability, for different sets of parameter values we may have saddle stability, totally unstable equilibria and chaotic attractors; (iii) for certain degrees of convexity and/or concavity of the Phillips curve, where endogenous fluctuations arise, one is able to encounter various results that seem intuitively correct. Firstly, when the Central Bank pays attention essentially to inflation targeting, the inflation rate has a lower mean and is less volatile; secondly, when the degree of price stickiness is high, the inflation rate displays a larger mean and higher volatility (but this is sensitive to the values given to the parameters of the model); and thirdly, the higher the target value of the output gap chosen by the Central Bank, the higher is the inflation rate and its volatility.
Optimal control of HIV/AIDS dynamic: Education and treatment
NASA Astrophysics Data System (ADS)
Sule, Amiru; Abdullah, Farah Aini
2014-07-01
A mathematical model which describes the transmission dynamics of HIV/AIDS is developed. The optimal control representing education and treatment for this model is explored. The existence of optimal Control is established analytically by the use of optimal control theory. Numerical simulations suggest that education and treatment for the infected has a positive impact on HIV/AIDS control.
Multi-objective optimization for deepwater dynamic umbilical installation analysis
NASA Astrophysics Data System (ADS)
Yang, HeZhen; Wang, AiJun; Li, HuaJun
2012-08-01
We suggest a method of multi-objective optimization based on approximation model for dynamic umbilical installation. The optimization aims to find out the most cost effective size, quantity and location of buoyancy modules for umbilical installation while maintaining structural safety. The approximation model is constructed by the design of experiment (DOE) sampling and is utilized to solve the problem of time-consuming analyses. The non-linear dynamic analyses considering environmental loadings are executed on these sample points from DOE. Non-dominated Sorting Genetic Algorithm (NSGA-II) is employed to obtain the Pareto solution set through an evolutionary optimization process. Intuitionist fuzzy set theory is applied for selecting the best compromise solution from Pareto set. The optimization results indicate this optimization strategy with approximation model and multiple attribute decision-making method is valid, and provide the optimal deployment method for deepwater dynamic umbilical buoyancy modules.
Structural optimization with constraints from dynamics in Lagrange
NASA Technical Reports Server (NTRS)
Pfeiffer, F.; Kneppe, G.; Ross, C.
1990-01-01
Structural optimization problems are mostly solved under constraints from statics, such as stresses, strains, or displacements under static loads. But in the design process, dynamic quantities like eigenfrequencies or accelerations under dynamic loads become more and more important. Therefore, it is obvious that constraints from dynamics must be considered in structural optimization packages. This paper addresses the dynamics branch in MBB-LAGRANGE. It will concentrate on two topics, namely on the different formulations for eigenfrequency constraints and on frequency response constraints. For the latter the necessity of a system reduction is emphasized. The methods implemented in LAGRANGE are presented and examples are given.
Dynamic systems of regional economy management optimization
NASA Astrophysics Data System (ADS)
Trofimov, S.; Kudzh, S.
One of the most actual problems of the Russian economic life is a regional economic systems formation. The hierarchy of economic and branch priorities should follow from the general idea of an industrial policy. The matter is that the concept of an industrial policy is defined by the system of priorities mainly incorporated in it. The problem of priorities is not solved yet neither on federal, nor at a regional level. It is necessary to recognize, that a substantiation of this or that variant of priorities - objectively a challenge. Such substantiation can be received with the help of dynamic structural modeling and management technology. At formation of the regional industrial policy program the special attention is given to creation of modern type commercial structures. In regions there are headquarters and branches of many largest corporations, holdings and banks. Besides it, many regional enterprises already became inter-regional or even the transnational companies. In this connection an assistance of transformation of the industrial enterprises and their groups in vertically integrated companies and modern type holdings can become a prominent aspect of an industrial policy. Regional economic structures should be reconstructed gradually on the general model of the world class competitive companies. Assistance to creation of new corporational control systems, the organization of headquarters and the central services work - all this can be included into the sphere of regional administration industrial policy. The special attention should be turned on necessity of development of own system of the corporate structures, capable to provide to the region an independent participation in use of the natural resources and industrial-technological potential, at the stage of a regional industrial policy program formation. Transformation of the industrial enterprises and their groups into modern type vertically-integrated companies and holdings can become one of the major directions of an industrial policy of region. The situational-analytical centers (SAC) of regional administration The major component of SAC is dynamic modeling, analysis, forecasting and optimization systems, based on modern intellectual information technologies. Spheres of SAC are not only financial streams management and investments optimization, but also strategic forecasting functions, which provide an optimum choice, "aiming", search of optimum ways of regional development and corresponding investments. It is expedient to consider an opportunity of formation of the uniform organizational-methodical center of an industrial policy of region. This organization can be directly connected to the scheduled-analytical services of the largest economic structures, local authorities, the ministries and departments. Such "direct communication" is capable to provide an effective regional development strategic management. Anyway, the output on foreign markets demands concentration of resources and support of authorities. Offered measures are capable to provide a necessary coordination of efforts of a various level economic structures. For maintenance of a regional industrial policy an attraction of all newest methods of strategic planning and management is necessary. Their activity should be constructed on the basis of modern approaches of economic systems management, cause the essence of an industrial policy is finally reduced to an effective regional and corporate economic activities control centers formation. Opportunities of optimum regional economy planning and management as uniform system Approaches to planning regional economic systems can be different. We will consider some most effective methods of planning and control over a regional facilities condition. All of them are compact and evident, that allows to put them into the group of average complexity technologies. At the decision of problems of a regional resource management is rather perspective the so-called "topographical" approach, which is used by intellectual information technology "Dynamics of systems". According to it the realistic plan of regional economic system is created in the virtual space -directly on a computer desktop. And economic objects are displayed on evident schemes according to their real "geographical" structure. Each enterprise, bank, business-unit or the detached division of the company receives its own "Module" (area in working space of a spreadsheet). In result the general plan of regional economic system appears at planners. A whole real picture of all economic system functioning is recreated by such way. The idea of a method is obvious: the operator sees actual functioning of regional economy. This is promoted by "the friendly interface", allowing to display real objects as a clear symbols. The regional economy can be considered as a set of the separate enterprises connected by various economic communications. Constant monitoring of an infrastructure development, tracking of a cargoes transportation condition, supervision over following the ecological specifications by the regional enterprises, growth of housing and industrial building level, condition of communications, etc. is necessary for carrying out with the help of modern technologies of space shooting and satellite navigating systems. It will allow to obtain the data in an operative mode and will also help to the quickly modeling of a situation development variants, and to accept anticipatory administrative decisions. Other sources of the information are statistical directories and reports on a social condition in region: about a migration level and the population incomes, consumer's basket structure, demographic parameters - age of the capable population, a sexual and national attributes, etc. It is possible to attribute financial parameters to the third group: the regional budget condition data, a gain of investments into the regional economy, a growth of incomes in the regional budget from the enterprises taxes, etc. What purposes are the specified methods applied to? The regional companies management and their reforming. One of the major directions - transformation of the enterprises into the viable industrial companies, an effective economic structures formation. In this connection, the technology of financial streams management allows to solve another important problem. The created (or reorganized) structure can be checked up in various economic modes even before practical realization of the project. The received data can be applied to an indicative plans substantiation and control over their performance. They can be used at negotiations with the enterprises heads, which are spent for interests of a regional industrial policy. Such approach can become base for concrete (as figures, plans and schemes) coordination of long-term investment priorities of the largest firms and companies, their production programs coordination, the mutual markets formation and the cooperation communications stabilization, payment chains, the social and economic development programs creation. It opens opportunities for the regional branches, enterprises, the major economic centers coordinated management. An effective planning of price pools, investment programs, the control over tax revenues and financial streams becomes possible. All this are necessary elements of a regional authorities economic policy. The offered approach considerably surpasses by it's efficiency auditor, "paper" or traditional computer technology. Bulky computer centers are not capable to cope with the given class problems. At the same time the small mobile commissions of experts "armed" with modern intellectual technologies, are capable to solve rather scale analytical problems. On the basis of the offered methods can be generated means of automation of the regional industrial complex development management. They are focused on concrete objects of various scale: territorial complexes, branches, groups of the companies, the separate enterprises. Thus, the program of a regional industrial policy can be put on an effective scheduled-analytical basis.
An Optimization Framework for Dynamic, Distributed Real-Time Systems
NASA Technical Reports Server (NTRS)
Eckert, Klaus; Juedes, David; Welch, Lonnie; Chelberg, David; Bruggerman, Carl; Drews, Frank; Fleeman, David; Parrott, David; Pfarr, Barbara
2003-01-01
Abstract. This paper presents a model that is useful for developing resource allocation algorithms for distributed real-time systems .that operate in dynamic environments. Interesting aspects of the model include dynamic environments, utility and service levels, which provide a means for graceful degradation in resource-constrained situations and support optimization of the allocation of resources. The paper also provides an allocation algorithm that illustrates how to use the model for producing feasible, optimal resource allocations.
An Optimization Framework for Dynamic Hybrid Energy Systems
Wenbo Du; Humberto E Garcia; Christiaan J.J. Paredis
2014-03-01
A computational framework for the efficient analysis and optimization of dynamic hybrid energy systems (HES) is developed. A microgrid system with multiple inputs and multiple outputs (MIMO) is modeled using the Modelica language in the Dymola environment. The optimization loop is implemented in MATLAB, with the FMI Toolbox serving as the interface between the computational platforms. Two characteristic optimization problems are selected to demonstrate the methodology and gain insight into the system performance. The first is an unconstrained optimization problem that optimizes the dynamic properties of the battery, reactor and generator to minimize variability in the HES. The second problem takes operating and capital costs into consideration by imposing linear and nonlinear constraints on the design variables. The preliminary optimization results obtained in this study provide an essential step towards the development of a comprehensive framework for designing HES.
Dynamics systems vs. optimal control--a unifying view.
Schaal, Stefan; Mohajerian, Peyman; Ijspeert, Auke
2007-01-01
In the past, computational motor control has been approached from at least two major frameworks: the dynamic systems approach and the viewpoint of optimal control. The dynamic system approach emphasizes motor control as a process of self-organization between an animal and its environment. Nonlinear differential equations that can model entrainment and synchronization behavior are among the most favorable tools of dynamic systems modelers. In contrast, optimal control approaches view motor control as the evolutionary or development result of a nervous system that tries to optimize rather general organizational principles, e.g., energy consumption or accurate task achievement. Optimal control theory is usually employed to develop appropriate theories. Interestingly, there is rather little interaction between dynamic systems and optimal control modelers as the two approaches follow rather different philosophies and are often viewed as diametrically opposing. In this paper, we develop a computational approach to motor control that offers a unifying modeling framework for both dynamic systems and optimal control approaches. In discussions of several behavioral experiments and some theoretical and robotics studies, we demonstrate how our computational ideas allow both the representation of self-organizing processes and the optimization of movement based on reward criteria. Our modeling framework is rather simple and general, and opens opportunities to revisit many previous modeling results from this novel unifying view. PMID:17925262
Parallel dynamic programming for on-line flight path optimization
NASA Technical Reports Server (NTRS)
Slater, G. L.; Hu, K.
1989-01-01
Parallel systolic algorithms for dynamic programming(DP) and their respective hardware implementations are presented for a problem in on-line trajectory optimization. The method is applied to a model for helicopter flight path optimization through a complex constraint region. This problem has application to an air traffic control problem and also to a terrain following/threat avoidance problem.
Review of dynamic optimization methods in renewable natural resource management
Williams, B.K.
1989-01-01
In recent years, the applications of dynamic optimization procedures in natural resource management have proliferated. A systematic review of these applications is given in terms of a number of optimization methodologies and natural resource systems. The applicability of the methods to renewable natural resource systems are compared in terms of system complexity, system size, and precision of the optimal solutions. Recommendations are made concerning the appropriate methods for certain kinds of biological resource problems.
First principles molecular dynamics without self-consistent field optimization
Souvatzis, Petros; Niklasson, Anders M. N.
2014-01-28
We present a first principles molecular dynamics approach that is based on time-reversible extended Lagrangian Born-Oppenheimer molecular dynamics [A. M. N. Niklasson, Phys. Rev. Lett. 100, 123004 (2008)] in the limit of vanishing self-consistent field optimization. The optimization-free dynamics keeps the computational cost to a minimum and typically provides molecular trajectories that closely follow the exact Born-Oppenheimer potential energy surface. Only one single diagonalization and Hamiltonian (or Fockian) construction are required in each integration time step. The proposed dynamics is derived for a general free-energy potential surface valid at finite electronic temperatures within hybrid density functional theory. Even in the event of irregular functional behavior that may cause a dynamical instability, the optimization-free limit represents a natural starting guess for force calculations that may require a more elaborate iterative electronic ground state optimization. Our optimization-free dynamics thus represents a flexible theoretical framework for a broad and general class of ab initio molecular dynamics simulations.
Dynamic positioning configuration and its first-order optimization
NASA Astrophysics Data System (ADS)
Xue, Shuqiang; Yang, Yuanxi; Dang, Yamin; Chen, Wu
2014-02-01
Traditional geodetic network optimization deals with static and discrete control points. The modern space geodetic network is, on the other hand, composed of moving control points in space (satellites) and on the Earth (ground stations). The network configuration composed of these facilities is essentially dynamic and continuous. Moreover, besides the position parameter which needs to be estimated, other geophysical information or signals can also be extracted from the continuous observations. The dynamic (continuous) configuration of the space network determines whether a particular frequency of signals can be identified by this system. In this paper, we employ the functional analysis and graph theory to study the dynamic configuration of space geodetic networks, and mainly focus on the optimal estimation of the position and clock-offset parameters. The principle of the D-optimization is introduced in the Hilbert space after the concept of the traditional discrete configuration is generalized from the finite space to the infinite space. It shows that the D-optimization developed in the discrete optimization is still valid in the dynamic configuration optimization, and this is attributed to the natural generalization of least squares from the Euclidean space to the Hilbert space. Then, we introduce the principle of D-optimality invariance under the combination operation and rotation operation, and propose some D-optimal simplex dynamic configurations: (1) (Semi) circular configuration in 2-dimensional space; (2) the D-optimal cone configuration and D-optimal helical configuration which is close to the GPS constellation in 3-dimensional space. The initial design of GPS constellation can be approximately treated as a combination of 24 D-optimal helixes by properly adjusting the ascending node of different satellites to realize a so-called Walker constellation. In the case of estimating the receiver clock-offset parameter, we show that the circular configuration, the symmetrical cone configuration and helical curve configuration are still D-optimal. It shows that the given total observation time determines the optimal frequency (repeatability) of moving known points and vice versa, and one way to improve the repeatability is to increase the rotational speed. Under the Newton's law of motion, the frequency of satellite motion determines the orbital altitude. Furthermore, we study three kinds of complex dynamic configurations, one of which is the combination of D-optimal cone configurations and a so-called Walker constellation composed of D-optimal helical configuration, the other is the nested cone configuration composed of n cones, and the last is the nested helical configuration composed of n orbital planes. It shows that an effective way to achieve high coverage is to employ the configuration composed of a certain number of moving known points instead of the simplex configuration (such as D-optimal helical configuration), and one can use the D-optimal simplex solutions or D-optimal complex configurations in any combination to achieve powerful configurations with flexile coverage and flexile repeatability. Alternately, how to optimally generate and assess the discrete configurations sampled from the continuous one is discussed. The proposed configuration optimization framework has taken the well-known regular polygons (such as equilateral triangle and quadrangular) in two-dimensional space and regular polyhedrons (regular tetrahedron, cube, regular octahedron, regular icosahedron, or regular dodecahedron) into account. It shows that the conclusions made by the proposed technique are more general and no longer limited by different sampling schemes. By the conditional equation of D-optimal nested helical configuration, the relevance issues of GNSS constellation optimization are solved and some examples are performed by GPS constellation to verify the validation of the newly proposed optimization technique. The proposed technique is potentially helpful in maintenance and quadratic optimization of single GNSS of which the orbital inclination and the orbital altitude change under the precession, as well as in optimally nesting GNSSs to perform global homogeneous coverage of the Earth.
Computer aided analysis and optimization of mechanical system dynamics
NASA Technical Reports Server (NTRS)
Haug, E. J.
1984-01-01
The purpose is to outline a computational approach to spatial dynamics of mechanical systems that substantially enlarges the scope of consideration to include flexible bodies, feedback control, hydraulics, and related interdisciplinary effects. Design sensitivity analysis and optimization is the ultimate goal. The approach to computer generation and solution of the system dynamic equations and graphical methods for creating animations as output is outlined.
Integrated Network Decompositions and Dynamic Programming for Graph Optimization (INDDGO)
2012-05-31
The INDDGO software package offers a set of tools for finding exact solutions to graph optimization problems via tree decompositions and dynamic programming algorithms. Currently the framework offers serial and parallel (distributed memory) algorithms for finding tree decompositions and solving the maximum weighted independent set problem. The parallel dynamic programming algorithm is implemented on top of the MADNESS task-based runtime.
Bridging Developmental Systems Theory and Evolutionary Psychology Using Dynamic Optimization
ERIC Educational Resources Information Center
Frankenhuis, Willem E.; Panchanathan, Karthik; Clark Barrett, H.
2013-01-01
Interactions between evolutionary psychologists and developmental systems theorists have been largely antagonistic. This is unfortunate because potential synergies between the two approaches remain unexplored. This article presents a method that may help to bridge the divide, and that has proven fruitful in biology: dynamic optimization. Dynamic
Neural dynamic optimization for control systems. I. Background.
Seong, C Y; Widrow, B
2001-01-01
The paper presents neural dynamic optimization (NDO) as a method of optimal feedback control for nonlinear multi-input-multi-output (MIMO) systems. The main feature of NDO is that it enables neural networks to approximate the optimal feedback solution whose existence dynamic programming (DP) justifies, thereby reducing the complexities of computation and storage problems of the classical methods such as DP. This paper mainly describes the background and motivations for the development of NDO, while the two other subsequent papers of this topic present the theory of NDO and demonstrate the method with several applications including control of autonomous vehicles and of a robot arm, respectively. PMID:18244815
Neural dynamic optimization for control systems.II. Theory.
Seong, C Y; Widrow, B
2001-01-01
The paper presents neural dynamic optimization (NDO) as a method of optimal feedback control for nonlinear multi-input-multi-output (MIMO) systems. The main feature of NDO is that it enables neural networks to approximate the optimal feedback solution whose existence dynamic programming (DP) justifies, thereby reducing the complexities of computation and storage problems of the classical methods such as DP. This paper mainly describes the theory of NDO, while the two other companion papers of this topic explain the background for the development of NDO and demonstrate the method with several applications including control of autonomous vehicles and of a robot arm, respectively. PMID:18244816
A generalized gradient algorithm for dynamic optimization
NASA Technical Reports Server (NTRS)
Zhao, Yiyuan; Bryson, A. E.; Slattery, R.
1989-01-01
A gradient algorithm is developed that determines optimal trajectories with path equality constraints and terminal constraints. A generalized gradient is formed which improves both the performance index and the path equality constraints simultaneously. The algorithm is extended to treat terminal constraints by using Bryson's impulse response technique. The main features of this algorithm are its numerical stability and smooth convergence near the optimum.
Optimal entangling capacity of dynamical processes
Campbell, Earl T.
2010-10-15
We investigate the entangling capacity of dynamical operations when provided with local ancilla. A comparison is made between the entangling capacity with and without the assistance of prior entanglement. An analytic solution is found for the log-negativity entangling capacity of two-qubit gates, which equals the entanglement of the Choi matrix isomorphic to the unitary operator. Surprisingly, the availability of prior entanglement does not affect this result, a property we call resource independence of the entangling capacity. We prove several useful upper bounds on the entangling capacity that hold for general qudit dynamical operations and for a whole family of entanglement monotones including log negativity and log robustness. The log-robustness entangling capacity is shown to be resource independent for general dynamics. We provide numerical results supporting a conjecture that the log-negativity entangling capacity is resource independent for all two-qudit unitary operators.
Structural optimization of rotor blades with integrated dynamics and aerodynamics
NASA Technical Reports Server (NTRS)
Chattopadhyay, Aditi; Walsh, Joanne L.
1988-01-01
The problem of structural optimization of helicopter rotor blades with integrated dynamic and aerodynamic design considerations is addressed. Results of recent optimization work on rotor blades for minimum weight with constraints on multiple coupled natural flap-lag frequencies, blade autorotational inertia and centrifugal stress has been reviewed. A strategy has been defined for the ongoing activities in the integrated dynamic/aerodynamic optimization of rotor blades. As a first step, the integrated dynamic/airload optimization problem has been formulated. To calculate system sensitivity derivatives necessary for the optimization recently developed, Global Sensitivity Equations (GSE) are being investigated. A need for multiple objective functions for the integrated optimization problem has been demonstrated and various techniques for solving the multiple objective function optimization are being investigated. The method called the Global Criteria Approach has been applied to a test problem with the blade in vacuum and the blade weight and the centrifugal stress as the multiple objectives. The results indicate that the method is quite effective in solving optimization problems with conflicting objective functions.
Practical synchronization on complex dynamical networks via optimal pinning control.
Li, Kezan; Sun, Weigang; Small, Michael; Fu, Xinchu
2015-07-01
We consider practical synchronization on complex dynamical networks under linear feedback control designed by optimal control theory. The control goal is to minimize global synchronization error and control strength over a given finite time interval, and synchronization error at terminal time. By utilizing the Pontryagin's minimum principle, and based on a general complex dynamical network, we obtain an optimal system to achieve the control goal. The result is verified by performing some numerical simulations on Star networks, Watts-Strogatz networks, and Barabási-Albert networks. Moreover, by combining optimal control and traditional pinning control, we propose an optimal pinning control strategy which depends on the network's topological structure. Obtained results show that optimal pinning control is very effective for synchronization control in real applications. PMID:26274112
Practical synchronization on complex dynamical networks via optimal pinning control
NASA Astrophysics Data System (ADS)
Li, Kezan; Sun, Weigang; Small, Michael; Fu, Xinchu
2015-07-01
We consider practical synchronization on complex dynamical networks under linear feedback control designed by optimal control theory. The control goal is to minimize global synchronization error and control strength over a given finite time interval, and synchronization error at terminal time. By utilizing the Pontryagin's minimum principle, and based on a general complex dynamical network, we obtain an optimal system to achieve the control goal. The result is verified by performing some numerical simulations on Star networks, Watts-Strogatz networks, and Barabási-Albert networks. Moreover, by combining optimal control and traditional pinning control, we propose an optimal pinning control strategy which depends on the network's topological structure. Obtained results show that optimal pinning control is very effective for synchronization control in real applications.
Dynamic optimization problems with bounded terminal conditions
NASA Technical Reports Server (NTRS)
Lee, A. Y.
1987-01-01
Bounded terminal conditions of nonlinear optimization problems are converted to equality terminal conditions via Valentine's device. In so doing, additional unknown parameters are introduced into the problem. The transformed problems can still be easily solved using the sequential gradient-restoration algorithm (SGRA) via a simple augmentation of the unknown parameter vector pi. Three example problems with bounded terminal conditions are solved to verify this technique.
Energy-optimal path planning by stochastic dynamically orthogonal level-set optimization
NASA Astrophysics Data System (ADS)
Subramani, Deepak N.; Lermusiaux, Pierre F. J.
2016-04-01
A stochastic optimization methodology is formulated for computing energy-optimal paths from among time-optimal paths of autonomous vehicles navigating in a dynamic flow field. Based on partial differential equations, the methodology rigorously leverages the level-set equation that governs time-optimal reachability fronts for a given relative vehicle-speed function. To set up the energy optimization, the relative vehicle-speed and headings are considered to be stochastic and new stochastic Dynamically Orthogonal (DO) level-set equations are derived. Their solution provides the distribution of time-optimal reachability fronts and corresponding distribution of time-optimal paths. An optimization is then performed on the vehicle's energy-time joint distribution to select the energy-optimal paths for each arrival time, among all stochastic time-optimal paths for that arrival time. Numerical schemes to solve the reduced stochastic DO level-set equations are obtained, and accuracy and efficiency considerations are discussed. These reduced equations are first shown to be efficient at solving the governing stochastic level-sets, in part by comparisons with direct Monte Carlo simulations. To validate the methodology and illustrate its accuracy, comparisons with semi-analytical energy-optimal path solutions are then completed. In particular, we consider the energy-optimal crossing of a canonical steady front and set up its semi-analytical solution using a energy-time nested nonlinear double-optimization scheme. We then showcase the inner workings and nuances of the energy-optimal path planning, considering different mission scenarios. Finally, we study and discuss results of energy-optimal missions in a wind-driven barotropic quasi-geostrophic double-gyre ocean circulation.
Combining optimal control theory and molecular dynamics for protein folding.
Arkun, Yaman; Gur, Mert
2012-01-01
A new method to develop low-energy folding routes for proteins is presented. The novel aspect of the proposed approach is the synergistic use of optimal control theory with Molecular Dynamics (MD). In the first step of the method, optimal control theory is employed to compute the force field and the optimal folding trajectory for the Cα atoms of a Coarse-Grained (CG) protein model. The solution of this CG optimization provides an harmonic approximation of the true potential energy surface around the native state. In the next step CG optimization guides the MD simulation by specifying the optimal target positions for the Cα atoms. In turn, MD simulation provides an all-atom conformation whose Cα positions match closely the reference target positions determined by CG optimization. This is accomplished by Targeted Molecular Dynamics (TMD) which uses a bias potential or harmonic restraint in addition to the usual MD potential. Folding is a dynamical process and as such residues make different contacts during the course of folding. Therefore CG optimization has to be reinitialized and repeated over time to accomodate these important changes. At each sampled folding time, the active contacts among the residues are recalculated based on the all-atom conformation obtained from MD. Using the new set of contacts, the CG potential is updated and the CG optimal trajectory for the Cα atoms is recomputed. This is followed by MD. Implementation of this repetitive CG optimization-MD simulation cycle generates the folding trajectory. Simulations on a model protein Villin demonstrate the utility of the method. Since the method is founded on the general tools of optimal control theory and MD without any restrictions, it is widely applicable to other systems. It can be easily implemented with available MD software packages. PMID:22238629
Combining Optimal Control Theory and Molecular Dynamics for Protein Folding
Arkun, Yaman; Gur, Mert
2012-01-01
A new method to develop low-energy folding routes for proteins is presented. The novel aspect of the proposed approach is the synergistic use of optimal control theory with Molecular Dynamics (MD). In the first step of the method, optimal control theory is employed to compute the force field and the optimal folding trajectory for the atoms of a Coarse-Grained (CG) protein model. The solution of this CG optimization provides an harmonic approximation of the true potential energy surface around the native state. In the next step CG optimization guides the MD simulation by specifying the optimal target positions for the atoms. In turn, MD simulation provides an all-atom conformation whose positions match closely the reference target positions determined by CG optimization. This is accomplished by Targeted Molecular Dynamics (TMD) which uses a bias potential or harmonic restraint in addition to the usual MD potential. Folding is a dynamical process and as such residues make different contacts during the course of folding. Therefore CG optimization has to be reinitialized and repeated over time to accomodate these important changes. At each sampled folding time, the active contacts among the residues are recalculated based on the all-atom conformation obtained from MD. Using the new set of contacts, the CG potential is updated and the CG optimal trajectory for the atoms is recomputed. This is followed by MD. Implementation of this repetitive CG optimization - MD simulation cycle generates the folding trajectory. Simulations on a model protein Villin demonstrate the utility of the method. Since the method is founded on the general tools of optimal control theory and MD without any restrictions, it is widely applicable to other systems. It can be easily implemented with available MD software packages. PMID:22238629
Dynamical optimization theory of a diversified portfolio
NASA Astrophysics Data System (ADS)
Marsili, Matteo; Maslov, Sergei; Zhang, Yi-Cheng
We propose and study a simple model of dynamical redistribution of capital in a diversified portfolio. We consider a hypothetical situation of a portfolio composed of N uncorrelated stocks. Each stock price follows a multiplicative random walk with identical drift and dispersion. The rules of our model naturally give rise to power law tails in the distribution of capital fractions invested in different stocks. The exponent of this scale free distribution is calculated in both discrete and continuous time formalism. It is demonstrated that the dynamical redistribution strategy results in a larger typical growth rate of the capital than a static “buy-and-hold” strategy. In the large N limit the typical growth rate is shown to asymptotically approach that of the expectation value of the stock price. The finite dimensional variant of the model is shown to describe the partition function of directed polymers in random media.
Optimal dynamic remapping of parallel computations
NASA Technical Reports Server (NTRS)
Nicol, David M.; Reynolds, Paul F., Jr.
1987-01-01
A large class of computations are characterized by a sequence of phases, with phase changes occurring unpredictably. The decision problem was considered regarding the remapping of workload to processors in a parallel computation when the utility of remapping and the future behavior of the workload is uncertain, and phases exhibit stable execution requirements during a given phase, but requirements may change radically between phases. For these problems a workload assignment generated for one phase may hinder performance during the next phase. This problem is treated formally for a probabilistic model of computation with at most two phases. The fundamental problem of balancing the expected remapping performance gain against the delay cost was addressed. Stochastic dynamic programming is used to show that the remapping decision policy minimizing the expected running time of the computation has an extremely simple structure. Because the gain may not be predictable, the performance of a heuristic policy that does not require estimnation of the gain is examined. The heuristic method's feasibility is demonstrated by its use on an adaptive fluid dynamics code on a multiprocessor. The results suggest that except in extreme cases, the remapping decision problem is essentially that of dynamically determining whether gain can be achieved by remapping after a phase change. The results also suggest that this heuristic is applicable to computations with more than two phases.
Dynamic optimization of district energy grid
NASA Astrophysics Data System (ADS)
Salsbery, Scott
The University of Iowa Power Plant operates utility generation and distribution for campus facilities, including electricity, steam, and chilled water. It is desirable to evaluate the optimal load combination of boilers, engines and chillers to meet the demand at minimal cost, particularly for future demand scenarios. An algorithm has been developed which takes into account the performance of individual units as part of the mix which ultimately supplies the campus and determine the degree that each should be operating to most efficiently meet demand. The algorithm is part of an integrated simulation tool which is specifically designed to apply traditional optimization techniques for a given (both current and possible) circumstance. The second component is to couple the algorithm with accurate estimates and historical data through which expected demand could be predicted. The simulation tool can account for any theoretical circumstance, which will be highly beneficial for strategic planning. As part of the process it is also necessary to determine the unique operating characteristics of the system components. The algorithms rely upon performance curves of individual system components (boiler, chiller, etc.) and those must be developed and refined when possible from experimental testing and commissioning or manufacturer supplied data.
Optimal motor control may mask sensory dynamics
Kiemel, Tim; Cowan, Noah J.; Jeka, John J.
2009-01-01
Properties of neural controllers for closed-loop sensorimotor behavior can be inferred with system identification. Under the standard paradigm, the closed-loop system is perturbed (input), measurements are taken (output), and the relationship between input and output reveals features of the system under study. Here we show that under common assumptions made about such systems (e.g. the system implements optimal control with a penalty on mechanical, but not sensory, states) important aspects of the neural controller (its zeros mask the modes of the sensors) remain hidden from standard system identification techniques. Only by perturbing or measuring the closed-loop system “between” the sensor and the control can these features be exposed with closed-loop system identification methods; while uncommon, there exist noninvasive techniques such as galvanic vestibular stimulation that perturb between sensor and controller in this way. PMID:19408009
Synthesis of standing-up trajectories using dynamic optimization.
Kuzelicki, Jernej; Zefran, Milos; Burger, Helena; Bajd, Tadej
2005-01-01
Dynamic optimization as a tool to compute standing-up trajectories was investigated. Sit-to-stand manoeuvres in five intact persons and five trans-femoral amputees were measured. Movements and ground reaction forces acting on the body were recorded. A five-segment 3D dynamic model of standing-up was developed. In each particular subject, the optimization criterion which yielded trajectories that best resemble the measured standing-up movement was determined. Since the intact persons used considerably different criteria in choosing the standing-up trajectories than the amputees, the optimal trajectories were computed by minimizing cost functionals (CF) with distinctive structures for each group of individuals. In intact persons, a unique cost functional was found which yielded realistic standing-up manoeuvres. In amputees, subject-specific sets of parameters indicating slightly different preferences in optimizing the effort of particular muscle groups were used. PMID:15536029
Global dynamic optimization approach to predict activation in metabolic pathways
2014-01-01
Background During the last decade, a number of authors have shown that the genetic regulation of metabolic networks may follow optimality principles. Optimal control theory has been succesfully used to compute optimal enzyme profiles considering simple metabolic pathways. However, applying this optimal control framework to more general networks (e.g. branched networks, or networks incorporating enzyme production dynamics) yields problems that are analytically intractable and/or numerically very challenging. Further, these previous studies have only considered a single-objective framework. Results In this work we consider a more general multi-objective formulation and we present solutions based on recent developments in global dynamic optimization techniques. We illustrate the performance and capabilities of these techniques considering two sets of problems. First, we consider a set of single-objective examples of increasing complexity taken from the recent literature. We analyze the multimodal character of the associated non linear optimization problems, and we also evaluate different global optimization approaches in terms of numerical robustness, efficiency and scalability. Second, we consider generalized multi-objective formulations for several examples, and we show how this framework results in more biologically meaningful results. Conclusions The proposed strategy was used to solve a set of single-objective case studies related to unbranched and branched metabolic networks of different levels of complexity. All problems were successfully solved in reasonable computation times with our global dynamic optimization approach, reaching solutions which were comparable or better than those reported in previous literature. Further, we considered, for the first time, multi-objective formulations, illustrating how activation in metabolic pathways can be explained in terms of the best trade-offs between conflicting objectives. This new methodology can be applied to metabolic networks with arbitrary topologies, non-linear dynamics and constraints. PMID:24393148
Experimental Testing of Dynamically Optimized Photoelectron Beams
NASA Astrophysics Data System (ADS)
Rosenzweig, J. B.; Cook, A. M.; Dunning, M.; England, R. J.; Musumeci, P.; Bellaveglia, M.; Boscolo, M.; Catani, L.; Cianchi, A.; Di Pirro, G.; Ferrario, M.; Fillipetto, D.; Gatti, G.; Palumbo, L.; Serafini, L.; Vicario, C.; Jones, S.
2006-11-01
We discuss the design of and initial results from an experiment in space-charge dominated beam dynamics which explores a new regime of high-brightness electron beam generation at the SPARC photoinjector. The scheme under study employs the tendency of intense electron beams to rearrange to produce uniform density, giving a nearly ideal beam from the viewpoint of space charge-induced emittance. The experiments are aimed at testing the marriage of this idea with a related concept, emittance compensation. We show that this new regime of operating photoinjector may be the preferred method of obtaining highest brightness beams with lower energy spread. We discuss the design of the experiment, including developing of a novel time-dependent, aerogel-based imaging system. This system has been installed at SPARC, and first evidence for nearly uniformly filled ellipsoidal charge distributions recorded.
Optimal dynamic bandwidth allocation for complex networks
NASA Astrophysics Data System (ADS)
Jiang, Zhong-Yuan; Liang, Man-Gui; Li, Qian; Guo, Dong-Chao
2013-03-01
Traffic capacity of one network strongly depends on the link’s bandwidth allocation strategy. In previous bandwidth allocation mechanisms, once one link’s bandwidth is allocated, it will be fixed throughout the overall traffic transmission process. However, the traffic load of every link changes from time to time. In this paper, with finite total bandwidth resource of the network, we propose to dynamically allocate the total bandwidth resource in which each link’s bandwidth is proportional to the queue length of the output buffer of the link per time step. With plenty of data packets in the network, the traffic handling ability of all links of the network achieves full utilization. The theoretical analysis and the extensive simulation results on complex networks are consistent. This work is valuable for network service providers to improve network performance or to do reasonable network design efficiently.
Experimental Testing of Dynamically Optimized Photoelectron Beams
Rosenzweig, J. B.; Cook, A. M.; Dunning, M.; England, R. J.; Musumeci, P.; Bellaveglia, M.; Boscolo, M.; Catani, L.; Cianchi, A.; Di Pirro, G.; Ferrario, M.; Fillipetto, D.; Gatti, G.; Palumbo, L.; Vicario, C.; Serafini, L.; Jones, S.
2006-11-27
We discuss the design of and initial results from an experiment in space-charge dominated beam dynamics which explores a new regime of high-brightness electron beam generation at the SPARC photoinjector. The scheme under study employs the tendency of intense electron beams to rearrange to produce uniform density, giving a nearly ideal beam from the viewpoint of space charge-induced emittance. The experiments are aimed at testing the marriage of this idea with a related concept, emittance compensation. We show that this new regime of operating photoinjector may be the preferred method of obtaining highest brightness beams with lower energy spread. We discuss the design of the experiment, including developing of a novel time-dependent, aerogel-based imaging system. This system has been installed at SPARC, and first evidence for nearly uniformly filled ellipsoidal charge distributions recorded.
Aerospace applications on integer and combinatorial optimization
NASA Technical Reports Server (NTRS)
Padula, S. L.; Kincaid, R. K.
1995-01-01
Research supported by NASA Langley Research Center includes many applications of aerospace design optimization and is conducted by teams of applied mathematicians and aerospace engineers. This paper investigates the benefits from this combined expertise in formulating and solving integer and combinatorial optimization problems. Applications range from the design of large space antennas to interior noise control. A typical problem. for example, seeks the optimal locations for vibration-damping devices on an orbiting platform and is expressed as a mixed/integer linear programming problem with more than 1500 design variables.
Aerospace applications of integer and combinatorial optimization
NASA Technical Reports Server (NTRS)
Padula, S. L.; Kincaid, R. K.
1995-01-01
Research supported by NASA Langley Research Center includes many applications of aerospace design optimization and is conducted by teams of applied mathematicians and aerospace engineers. This paper investigates the benefits from this combined expertise in solving combinatorial optimization problems. Applications range from the design of large space antennas to interior noise control. A typical problem, for example, seeks the optimal locations for vibration-damping devices on a large space structure and is expressed as a mixed/integer linear programming problem with more than 1500 design variables.
Aerospace Applications of Integer and Combinatorial Optimization
NASA Technical Reports Server (NTRS)
Padula, S. L.; Kincaid, R. K.
1995-01-01
Research supported by NASA Langley Research Center includes many applications of aerospace design optimization and is conducted by teams of applied mathematicians and aerospace engineers. This paper investigates the benefits from this combined expertise in formulating and solving integer and combinatorial optimization problems. Applications range from the design of large space antennas to interior noise control. A typical problem, for example, seeks the optimal locations for vibration-damping devices on an orbiting platform and is expressed as a mixed/integer linear programming problem with more than 1500 design variables.
MULTIOBJECTIVE DYNAMIC APERTURE OPTIMIZATION AT NSLS-II
Yang, L.; Li, Y.; Guo, W.; Krinsky, S.
2011-03-28
In this paper we present a multiobjective approach to the dynamic aperture (DA) optimization. Taking the NSLS-II lattice as an example, we have used both sextupoles and quadrupoles as tuning variables to optimize both on-momentum and off-momentum DA. The geometric and chromatic sextupoles are used for nonlinear properties while the tunes are independently varied by quadrupoles. The dispersion and emittance are fixed during tunes variation. The algorithms, procedures, performances and results of our optimization of DA will be discussed and they are found to be robust, general and easy to apply to similar problems.
Dynamics and optimal number of replicas in parallel tempering simulations
NASA Astrophysics Data System (ADS)
Nadler, Walter; Hansmann, Ulrich H. E.
2007-12-01
We study the dynamics of parallel tempering simulations, also known as the replica exchange technique, which has become the method of choice for simulation of proteins and other complex systems. Recent results for the optimal choice of the control parameter discretization allow a treatment independent of the system in question. By analyzing mean first passage times across the control parameter space, we find an expression for the optimal number of replicas in simulations covering a given temperature range. Our results suggest a particular protocol to optimize the number of replicas in actual simulations.
Dynamically optimal strategies for managing resistance to genetically modified crops.
Qiao, Fangbin; Wilen, James; Rozelle, Scott
2008-06-01
This paper develops a dynamic model of the evolution of pest a population and pest resistance to characterize the socially optimal refuge strategy for managing a pest's resistance to genetically modified crops. Previous theoretical economic analyses of this problem focus on steady states; we also address refuge policies along the optimal path to the final equilibrium. To elaborate on our theoretical analysis of the resistance problem, we develop a simulation model calibrated to cotton (Gossypium spp.) production in China. Our results show the importance of fitness cost as a determinant of the qualitative nature of optimal refuge policies. PMID:18613595
Developing learning algorithms via optimized discretization of continuous dynamical systems.
Tao, Qing; Sun, Zhengya; Kong, Kang
2012-02-01
Most of the existing numerical optimization methods are based upon a discretization of some ordinary differential equations. In order to solve some convex and smooth optimization problems coming from machine learning, in this paper, we develop efficient batch and online algorithms based on a new principle, i.e., the optimized discretization of continuous dynamical systems (ODCDSs). First, a batch learning projected gradient dynamical system with Lyapunov's stability and monotonic property is introduced, and its dynamical behavior guarantees the accuracy of discretization-based optimizer and applicability of line search strategy. Furthermore, under fair assumptions, a new online learning algorithm achieving regret O(?T) or O(logT) is obtained. By using the line search strategy, the proposed batch learning ODCDS exhibits insensitivity to the step sizes and faster decrease. With only a small number of line search steps, the proposed stochastic algorithm shows sufficient stability and approximate optimality. Experimental results demonstrate the correctness of our theoretical analysis and efficiency of our algorithms. PMID:21880573
Bridging Developmental Systems Theory and Evolutionary Psychology Using Dynamic Optimization
ERIC Educational Resources Information Center
Frankenhuis, Willem E.; Panchanathan, Karthik; Clark Barrett, H.
2013-01-01
Interactions between evolutionary psychologists and developmental systems theorists have been largely antagonistic. This is unfortunate because potential synergies between the two approaches remain unexplored. This article presents a method that may help to bridge the divide, and that has proven fruitful in biology: dynamic optimization. Dynamic…
Aerodynamic design optimization with sensitivity analysis and computational fluid dynamics
NASA Technical Reports Server (NTRS)
Baysal, Oktay
1995-01-01
An investigation was conducted from October 1, 1990 to May 31, 1994 on the development of methodologies to improve the designs (more specifically, the shape) of aerodynamic surfaces of coupling optimization algorithms (OA) with Computational Fluid Dynamics (CFD) algorithms via sensitivity analyses (SA). The study produced several promising methodologies and their proof-of-concept cases, which have been reported in the open literature.
Multiobjective Optimization of Low-Energy Trajectories Using Optimal Control on Dynamical Channels
NASA Technical Reports Server (NTRS)
Coffee, Thomas M.; Anderson, Rodney L.; Lo, Martin W.
2011-01-01
We introduce a computational method to design efficient low-energy trajectories by extracting initial solutions from dynamical channels formed by invariant manifolds, and improving these solutions through variational optimal control. We consider trajectories connecting two unstable periodic orbits in the circular restricted 3-body problem (CR3BP). Our method leverages dynamical channels to generate a range of solutions, and approximates the areto front for impulse and time of flight through a multiobjective optimization of these solutions based on primer vector theory. We demonstrate the application of our method to a libration orbit transfer in the Earth-Moon system.
A dynamic optimization model for solid waste recycling.
Anghinolfi, Davide; Paolucci, Massimo; Robba, Michela; Taramasso, Angela Celeste
2013-02-01
Recycling is an important part of waste management (that includes different kinds of issues: environmental, technological, economic, legislative, social, etc.). Differently from many works in literature, this paper is focused on recycling management and on the dynamic optimization of materials collection. The developed dynamic decision model is characterized by state variables, corresponding to the quantity of waste in each bin per each day, and control variables determining the quantity of material that is collected in the area each day and the routes for collecting vehicles. The objective function minimizes the sum of costs minus benefits. The developed decision model is integrated in a GIS-based Decision Support System (DSS). A case study related to the Cogoleto municipality is presented to show the effectiveness of the proposed model. From optimal results, it has been found that the net benefits of the optimized collection are about 2.5 times greater than the estimated current policy. PMID:23158873
Neural dynamic optimization for control systems.III. Applications.
Seong, C Y; Widrow, B
2001-01-01
For pt.II. see ibid., p. 490-501. The paper presents neural dynamic optimization (NDO) as a method of optimal feedback control for nonlinear multi-input-multi-output (MIMO) systems. The main feature of NDO is that it enables neural networks to approximate the optimal feedback solution whose existence dynamic programming (DP) justifies, thereby reducing the complexities of computation and storage problems of the classical methods such as DP. This paper demonstrates NDO with several applications including control of autonomous vehicles and of a robot-arm, while the two other companion papers of this topic describes the background for the development of NDO and present the theory of the method, respectively. PMID:18244817
Analysis and Optimization of Pulse Dynamics for Magnetic Stimulation
Goetz, Stefan M.; Truong, Cong Nam; Gerhofer, Manuel G.; Peterchev, Angel V.; Herzog, Hans-Georg; Weyh, Thomas
2013-01-01
Magnetic stimulation is a standard tool in brain research and has found important clinical applications in neurology, psychiatry, and rehabilitation. Whereas coil designs and the spatial field properties have been intensively studied in the literature, the temporal dynamics of the field has received less attention. Typically, the magnetic field waveform is determined by available device circuit topologies rather than by consideration of what is optimal for neural stimulation. This paper analyzes and optimizes the waveform dynamics using a nonlinear model of a mammalian axon. The optimization objective was to minimize the pulse energy loss. The energy loss drives power consumption and heating, which are the dominating limitations of magnetic stimulation. The optimization approach is based on a hybrid global-local method. Different coordinate systems for describing the continuous waveforms in a limited parameter space are defined for numerical stability. The optimization results suggest that there are waveforms with substantially higher efficiency than that of traditional pulse shapes. One class of optimal pulses is analyzed further. Although the coil voltage profile of these waveforms is almost rectangular, the corresponding current shape presents distinctive characteristics, such as a slow low-amplitude first phase which precedes the main pulse and reduces the losses. Representatives of this class of waveforms corresponding to different maximum voltages are linked by a nonlinear transformation. The main phase, however, scales with time only. As with conventional magnetic stimulation pulses, briefer pulses result in lower energy loss but require higher coil voltage than longer pulses. PMID:23469168
Dynamic characterization of bolted joints using FRF decoupling and optimization
NASA Astrophysics Data System (ADS)
Tol, ?erife; Ozguven, H. Nevzat
2015-03-01
Mechanical connections play a significant role in predicting dynamic characteristics of assembled structures. Therefore, equivalent dynamic models for joints are needed. Due to the complexity of joints, it is difficult to describe joint dynamics with analytical models. Reliable models are generally obtained using experimental measurements. In this paper an experimental identification method based on FRF decoupling and optimization algorithm is proposed for modeling joints. In the method the FRFs of two substructures connected with a joint are measured, while the FRFs of the substructures are obtained numerically or experimentally. Then the joint properties are calculated in terms of translational, rotational and cross-coupling stiffness and damping values by using FRF decoupling. In order to eliminate the numerical errors associated with matrix inversion an optimization algorithm is used to update the joint values obtained from FRF decoupling. The validity of the proposed method is demonstrated with experimental studies with bolted joints.
Optimal experimental dynamical decoupling of both longitudinal and transverse relaxations
NASA Astrophysics Data System (ADS)
Zhen, Xing-Long; Zhang, Fei-Hao; Feng, Guanru; Li, Hang; Long, Gui-Lu
2016-02-01
Both longitudinal and transverse relaxations exist in the practical environment. Their simultaneous eliminations are extremely demanding in real applications. Previous experimental work has focused mainly on the suppression of transverse relaxation. In this paper we investigate the performance of three important dynamical decoupling schemes—quadratic dynamical decoupling, periodic dynamical decoupling, and concatenated dynamical decoupling—in an environment with hybrid errors. We propose a method to engineer arbitrary environment by modulating the control field. The technique developed here is universal and can be applied to other quantum information processing systems. Three-dimensional filter functions technique is utilized to analyze the fidelity decay of a one-qubit state protected by dynamical decoupling sequences. This enables us to quantitatively compare the performance of different dynamical decoupling sequences and demonstrate the superiority of quadratic dynamical decoupling in experiments for the first time. Our work reveals that quadratic dynamical decoupling is optimal conditioned on the appropriate noise properties. The difference of constructing dynamical decoupling sequences with various Pauli pulses is also investigated.
Optimizing legacy molecular dynamics software with directive-based offload
NASA Astrophysics Data System (ADS)
Michael Brown, W.; Carrillo, Jan-Michael Y.; Gavhane, Nitin; Thakkar, Foram M.; Plimpton, Steven J.
2015-10-01
Directive-based programming models are one solution for exploiting many-core coprocessors to increase simulation rates in molecular dynamics. They offer the potential to reduce code complexity with offload models that can selectively target computations to run on the CPU, the coprocessor, or both. In this paper, we describe modifications to the LAMMPS molecular dynamics code to enable concurrent calculations on a CPU and coprocessor. We demonstrate that standard molecular dynamics algorithms can run efficiently on both the CPU and an x86-based coprocessor using the same subroutines. As a consequence, we demonstrate that code optimizations for the coprocessor also result in speedups on the CPU; in extreme cases up to 4.7X. We provide results for LAMMPS benchmarks and for production molecular dynamics simulations using the Stampede hybrid supercomputer with both Intel® Xeon Phi™ coprocessors and NVIDIA GPUs. The optimizations presented have increased simulation rates by over 2X for organic molecules and over 7X for liquid crystals on Stampede. The optimizations are available as part of the "Intel package" supplied with LAMMPS.
Optimal scales to observe habitat dynamics: a coral reef example.
Habeeb, Rebecca L; Johnson, Craig R; Wotherspoon, Simon; Mumby, Peter J
2007-04-01
A new technique to estimate the characteristic length scales (CLSs) of real ecological systems provides an objective means to identify the optimal scale(s) of observation to best detect underlying dynamical trends. Application of the technique to natural systems has focused on identifying appropriate scales to measure the dynamics of species as descriptors of community and ecosystem dynamics. However, ecosystem monitoring is often based not on assessing single species, but on species assemblages, functional groups, or habitat types. We asked whether the concept of CLSs based on dynamic interactions among species could be extended to examine interactions among habitat types and thus to identify optimal scales for observing habitat dynamics. A time series of three spatial maps of benthic habitats on a Caribbean coral reef was constructed from aerial photographs, Compact Airborne Spectrographic Imager (CASI) images, and IKONOS satellite images, providing the short time sequence required for this technique. We estimated the CLS based on the dynamics of three distinct habitat types: dense stands of seagrass, sparse stands of seagrass, and Montastrea patch reefs. Despite notable differences in the areal extent of and relative change in these habitats over the 21-year observation period, analyses based on each habitat type indicated a similar CLS of -300 m. We interpret the consistency of CLSs among habitats to indicate that the dynamics of the three habitat types are linked. The results are encouraging, and they indicate that CLS techniques can be used to identify the appropriate scale at which to monitor ecosystem trends on the basis of the dynamics of only one of a disparate suite of habitat types. PMID:17494384
Set-valued dynamic treatment regimes for competing outcomes
Laber, Eric B.; Lizotte, Daniel J.; Ferguson, Bradley
2014-01-01
Summary Dynamic treatment regimes operationalize the clinical decision process as a sequence of functions, one for each clinical decision, where each function maps up-to-date patient information to a single recommended treatment. Current methods for estimating optimal dynamic treatment regimes, for example Q-learning, require the specification of a single outcome by which the ‘goodness’ of competing dynamic treatment regimes is measured. However, this is an over-simplification of the goal of clinical decision making, which aims to balance several potentially competing outcomes, e.g., symptom relief and side-effect burden. When there are competing outcomes and patients do not know or cannot communicate their preferences, formation of a single composite outcome that correctly balances the competing outcomes is not possible. This problem also occurs when patient preferences evolve over time. We propose a method for constructing dynamic treatment regimes that accommodates competing outcomes by recommending sets of treatments at each decision point. Formally, we construct a sequence of set-valued functions that take as input up-to-date patient information and give as output a recommended subset of the possible treatments. For a given patient history, the recommended set of treatments contains all treatments that produce non-inferior outcome vectors. Constructing these set-valued functions requires solving a non-trivial enumeration problem. We offer an exact enumeration algorithm by recasting the problem as a linear mixed integer program. The proposed methods are illustrated using data from the CATIE schizophrenia study. PMID:24400912
Power distribution system planning with reliability modeling and optimization
Tang, Y.
1996-02-01
A new approach for the systemized optimization of power distribution systems is presented in this paper. Distribution system reliability is modeled in the optimization objective function via outage costs and costs of switching devices, along with the nonlinear costs of investment, maintenance and energy losses of both the substations and the feeders. The optimization model established is multi-stage, mixed-integer and nonlinear, which is solved by a network-flow programming algorithm. A multi-stage interlacing strategy and a nonlinearity iteration method are also designed. Supported by an extensive database, the planning software tool has been applied to optimize the power distribution system of a developing city.
Optimization Research of Generation Investment Based on Linear Programming Model
NASA Astrophysics Data System (ADS)
Wu, Juan; Ge, Xueqian
Linear programming is an important branch of operational research and it is a mathematical method to assist the people to carry out scientific management. GAMS is an advanced simulation and optimization modeling language and it will combine a large number of complex mathematical programming, such as linear programming LP, nonlinear programming NLP, MIP and other mixed-integer programming with the system simulation. In this paper, based on the linear programming model, the optimized investment decision-making of generation is simulated and analyzed. At last, the optimal installed capacity of power plants and the final total cost are got, which provides the rational decision-making basis for optimized investments.
Artificial bee colony algorithm for constrained possibilistic portfolio optimization problem
NASA Astrophysics Data System (ADS)
Chen, Wei
2015-07-01
In this paper, we discuss the portfolio optimization problem with real-world constraints under the assumption that the returns of risky assets are fuzzy numbers. A new possibilistic mean-semiabsolute deviation model is proposed, in which transaction costs, cardinality and quantity constraints are considered. Due to such constraints the proposed model becomes a mixed integer nonlinear programming problem and traditional optimization methods fail to find the optimal solution efficiently. Thus, a modified artificial bee colony (MABC) algorithm is developed to solve the corresponding optimization problem. Finally, a numerical example is given to illustrate the effectiveness of the proposed model and the corresponding algorithm.
Sequential activation of metabolic pathways: a dynamic optimization approach.
Oyarzún, Diego A; Ingalls, Brian P; Middleton, Richard H; Kalamatianos, Dimitrios
2009-11-01
The regulation of cellular metabolism facilitates robust cellular operation in the face of changing external conditions. The cellular response to this varying environment may include the activation or inactivation of appropriate metabolic pathways. Experimental and numerical observations of sequential timing in pathway activation have been reported in the literature. It has been argued that such patterns can be rationalized by means of an underlying optimal metabolic design. In this paper we pose a dynamic optimization problem that accounts for time-resource minimization in pathway activation under constrained total enzyme abundance. The optimized variables are time-dependent enzyme concentrations that drive the pathway to a steady state characterized by a prescribed metabolic flux. The problem formulation addresses unbranched pathways with irreversible kinetics. Neither specific reaction kinetics nor fixed pathway length are assumed.In the optimal solution, each enzyme follows a switching profile between zero and maximum concentration, following a temporal sequence that matches the pathway topology. This result provides an analytic justification of the sequential activation previously described in the literature. In contrast with the existent numerical approaches, the activation sequence is proven to be optimal for a generic class of monomolecular kinetics. This class includes, but is not limited to, Mass Action, Michaelis-Menten, Hill, and some Power-law models. This suggests that sequential enzyme expression may be a common feature of metabolic regulation, as it is a robust property of optimal pathway activation. PMID:19412635
Shape Optimization of Vehicle Radiator Using Computational Fluid Dynamics (cfd)
NASA Astrophysics Data System (ADS)
Maddipatla, Sridhar; Guessous, Laila
2002-11-01
Automotive manufacturers need to improve the efficiency and lifetime of all engine components. In the case of radiators, performance depends significantly on coolant flow homogeneity across the tubes and overall pressure drop between the inlet and outlet. Design improvements are especially needed in tube-flow uniformity to prevent premature fouling and failure of heat exchangers. Rather than relying on ad-hoc geometry changes, the current study combines Computational Fluid Dynamics with shape optimization methods to improve radiator performance. The goal is to develop an automated suite of virtual tools to assist in radiator design. Two objective functions are considered: a flow non-uniformity coefficient,Cf, and the overall pressure drop, dP*. The methodology used to automate the CFD and shape optimization procedures is discussed. In the first phase, single and multi-variable optimization methods, coupled with CFD, are applied to simplified 2-D radiator models to investigate effects of inlet and outlet positions on the above functions. The second phase concentrates on CFD simulations of a simplified 3-D radiator model. The results, which show possible improvements in both pressure and flow uniformity, validate the optimization criteria that were developed, as well as the potential of shape optimization methods with CFD to improve heat exchanger design. * Improving Radiator Design Through Shape Optimization, L. Guessous and S. Maddipatla, Paper # IMECE2002-33888, Proceedings of the 2002 ASME International Mechanical Engineering Congress and Exposition, November 2002
Tensor-optimized antisymmetrized molecular dynamics in nuclear physics
NASA Astrophysics Data System (ADS)
Myo, Takayuki; Toki, Hiroshi; Ikeda, Kiyomi; Horiuchi, Hisashi; Suhara, Tadahiro
2015-07-01
We develop a new formalism to treat nuclear many-body systems using the bare nucleon-nucleon interaction. It has become evident that the tensor interaction plays an important role in nuclear many-body systems due to the role of the pion in strongly interacting systems. We take the antisymmetrized molecular dynamics (AMD) as a basic framework and add a tensor correlation operator acting on the AMD wave function using the concept of the tensor-optimized shell model. We demonstrate a systematical and straightforward formulation utilizing the Gaussian integration and differentiation method and the antisymmetrization technique to calculate all the matrix elements of the many-body Hamiltonian. We can include the three-body interaction naturally and calculate the matrix elements systematically in the progressive order of the tensor correlation operator. We call the new formalism "tensor-optimized antisymmetrized molecular dynamics".
Optimizing airport runway improvement program - A dynamic programming approach
NASA Technical Reports Server (NTRS)
Yu, J. C.; Gibson, D. R.
1975-01-01
In order to reduce the air traffic delay in the terminal area, an immediate remedy is to increase airport capacity by an expansion of the existing runway system. The runway expansion program is often limited by budgetary constraints; the expensive facilities for a long-term improvement cannot be built at once. When a runway improvement strategy is being considered for a longer planning horizon, the investiment decision depends upon the interrelations of its composite periods. The problem, therefore, is to determine how time factor and investment decisions interact to yield an optimal improvement scheme that meets demand at a minimum cost. With this objective in mind, a dynamic programming methodology is employed to determine the optimal planning scheme. Also, an example runway improvement problem is tested to illustrate how a dynamic programming model is practical in actual application.
A Dynamically Optimized SSVEP Brain-Computer Interface (BCI) Speller.
Yin, Erwei; Zhou, Zongtan; Jiang, Jun; Yu, Yang; Hu, Dewen
2015-06-01
The aim of this study was to design a dynamically optimized steady-state visually evoked potential (SSVEP) brain-computer interface (BCI) system with enhanced performance relative to previous SSVEP BCIs in terms of the number of items selectable on the interface, accuracy, and speed. In this approach, the row/column (RC) paradigm was employed in a SSVEP speller to increase the number of items. The target is detected by subsequently determining the row and column coordinates. To improve spelling accuracy, we added a posterior processing after the canonical correlation analysis (CCA) approach to reduce the interfrequency variation between different subjects and named the new signal processing method CCA-RV, and designed a real-time biofeedback mechanism to increase attention on the visual stimuli. To achieve reasonable online spelling speed, both fixed and dynamic approaches for setting the optimal stimulus duration were implemented and compared. Experimental results for 11 subjects suggest that the CCA-RV method and the real-time biofeedback effectively increased accuracy compared with CCA and the absence of real-time feedback, respectively. In addition, both optimization approaches for setting stimulus duration achieved reasonable online spelling performance. However, the dynamic optimization approach yielded a higher practical information transfer rate (PITR) than the fixed optimization approach. The average online PITR achieved by the proposed adaptive SSVEP speller, including the time required for breaks between selections and error correction, was 41.08 bit/min. These results indicate that our BCI speller is promising for use in SSVEP-based BCI applications. PMID:24801483
Rethinking design parameters in the search for optimal dynamic seating.
Pynt, Jennifer
2015-04-01
Dynamic seating design purports to lessen damage incurred during sedentary occupations by increasing sitter movement while modifying muscle activity. Dynamic sitting is currently defined by O'Sullivan et al. ( 2013a) as relating to 'the increased motion in sitting which is facilitated by the use of specific chairs or equipment' (p. 628). Yet the evidence is conflicting that dynamic seating creates variation in the sitter's lumbar posture or muscle activity with the overall consensus being that current dynamic seating design fails to fulfill its goals. Research is needed to determine if a new generation of chairs requiring active sitter involvement fulfills the goals of dynamic seating and aids cardio/metabolic health. This paper summarises the pursuit of knowledge regarding optimal seated spinal posture and seating design. Four new forms of dynamic seating encouraging active sitting are discussed. These are 1) The Core-flex with a split seatpan to facilitate a walking action while seated 2) the Duo balans requiring body action to create rocking 3) the Back App and 4) Locus pedestal stools both using the sitter's legs to drive movement. Unsubstantiated claims made by the designers of these new forms of dynamic seating are outlined. Avenues of research are suggested to validate designer claims and investigate whether these designs fulfill the goals of dynamic seating and assist cardio/metabolic health. Should these claims be efficacious then a new definition of dynamic sitting is suggested; 'Sitting in which the action is provided by the sitter, while the dynamic mechanism of the chair accommodates that action'. PMID:25892386
Optimally combining dynamical decoupling and quantum error correction.
Paz-Silva, Gerardo A; Lidar, D A
2013-01-01
Quantum control and fault-tolerant quantum computing (FTQC) are two of the cornerstones on which the hope of realizing a large-scale quantum computer is pinned, yet only preliminary steps have been taken towards formalizing the interplay between them. Here we explore this interplay using the powerful strategy of dynamical decoupling (DD), and show how it can be seamlessly and optimally integrated with FTQC. To this end we show how to find the optimal decoupling generator set (DGS) for various subspaces relevant to FTQC, and how to simultaneously decouple them. We focus on stabilizer codes, which represent the largest contribution to the size of the DGS, showing that the intuitive choice comprising the stabilizers and logical operators of the code is in fact optimal, i.e., minimizes a natural cost function associated with the length of DD sequences. Our work brings hybrid DD-FTQC schemes, and their potentially considerable advantages, closer to realization. PMID:23559088
Optimally combining dynamical decoupling and quantum error correction
Paz-Silva, Gerardo A.; Lidar, D. A.
2013-01-01
Quantum control and fault-tolerant quantum computing (FTQC) are two of the cornerstones on which the hope of realizing a large-scale quantum computer is pinned, yet only preliminary steps have been taken towards formalizing the interplay between them. Here we explore this interplay using the powerful strategy of dynamical decoupling (DD), and show how it can be seamlessly and optimally integrated with FTQC. To this end we show how to find the optimal decoupling generator set (DGS) for various subspaces relevant to FTQC, and how to simultaneously decouple them. We focus on stabilizer codes, which represent the largest contribution to the size of the DGS, showing that the intuitive choice comprising the stabilizers and logical operators of the code is in fact optimal, i.e., minimizes a natural cost function associated with the length of DD sequences. Our work brings hybrid DD-FTQC schemes, and their potentially considerable advantages, closer to realization. PMID:23559088
Confronting dynamics and uncertainty in optimal decision making for conservation
Williams, Byron K.; Johnson, Fred A.
2013-01-01
The effectiveness of conservation efforts ultimately depends on the recognition that decision making, and the systems that it is designed to affect, are inherently dynamic and characterized by multiple sources of uncertainty. To cope with these challenges, conservation planners are increasingly turning to the tools of decision analysis, especially dynamic optimization methods. Here we provide a general framework for optimal, dynamic conservation and then explore its capacity for coping with various sources and degrees of uncertainty. In broadest terms, the dynamic optimization problem in conservation is choosing among a set of decision options at periodic intervals so as to maximize some conservation objective over the planning horizon. Planners must account for immediate objective returns, as well as the effect of current decisions on future resource conditions and, thus, on future decisions. Undermining the effectiveness of such a planning process are uncertainties concerning extant resource conditions (partial observability), the immediate consequences of decision choices (partial controllability), the outcomes of uncontrolled, environmental drivers (environmental variation), and the processes structuring resource dynamics (structural uncertainty). Where outcomes from these sources of uncertainty can be described in terms of probability distributions, a focus on maximizing the expected objective return, while taking state-specific actions, is an effective mechanism for coping with uncertainty. When such probability distributions are unavailable or deemed unreliable, a focus on maximizing robustness is likely to be the preferred approach. Here the idea is to choose an action (or state-dependent policy) that achieves at least some minimum level of performance regardless of the (uncertain) outcomes. We provide some examples of how the dynamic optimization problem can be framed for problems involving management of habitat for an imperiled species, conservation of a critically endangered population through captive breeding, control of invasive species, construction of biodiversity reserves, design of landscapes to increase habitat connectivity, and resource exploitation. Although these decision making problems and their solutions present significant challenges, we suggest that a systematic and effective approach to dynamic decision making in conservation need not be an onerous undertaking. The requirements are shared with any systematic approach to decision making--a careful consideration of values, actions, and outcomes.
Confronting dynamics and uncertainty in optimal decision making for conservation
NASA Astrophysics Data System (ADS)
Williams, Byron K.; Johnson, Fred A.
2013-06-01
The effectiveness of conservation efforts ultimately depends on the recognition that decision making, and the systems that it is designed to affect, are inherently dynamic and characterized by multiple sources of uncertainty. To cope with these challenges, conservation planners are increasingly turning to the tools of decision analysis, especially dynamic optimization methods. Here we provide a general framework for optimal, dynamic conservation and then explore its capacity for coping with various sources and degrees of uncertainty. In broadest terms, the dynamic optimization problem in conservation is choosing among a set of decision options at periodic intervals so as to maximize some conservation objective over the planning horizon. Planners must account for immediate objective returns, as well as the effect of current decisions on future resource conditions and, thus, on future decisions. Undermining the effectiveness of such a planning process are uncertainties concerning extant resource conditions (partial observability), the immediate consequences of decision choices (partial controllability), the outcomes of uncontrolled, environmental drivers (environmental variation), and the processes structuring resource dynamics (structural uncertainty). Where outcomes from these sources of uncertainty can be described in terms of probability distributions, a focus on maximizing the expected objective return, while taking state-specific actions, is an effective mechanism for coping with uncertainty. When such probability distributions are unavailable or deemed unreliable, a focus on maximizing robustness is likely to be the preferred approach. Here the idea is to choose an action (or state-dependent policy) that achieves at least some minimum level of performance regardless of the (uncertain) outcomes. We provide some examples of how the dynamic optimization problem can be framed for problems involving management of habitat for an imperiled species, conservation of a critically endangered population through captive breeding, control of invasive species, construction of biodiversity reserves, design of landscapes to increase habitat connectivity, and resource exploitation. Although these decision making problems and their solutions present significant challenges, we suggest that a systematic and effective approach to dynamic decision making in conservation need not be an onerous undertaking. The requirements are shared with any systematic approach to decision making—a careful consideration of values, actions, and outcomes.
Optimal forwarding ratio on dynamical networks with heterogeneous mobility
NASA Astrophysics Data System (ADS)
Gan, Yu; Tang, Ming; Yang, Hanxin
2013-05-01
Since the discovery of non-Poisson statistics of human mobility trajectories, more attention has been paid to understand the role of these patterns in different dynamics. In this study, we first introduce the heterogeneous mobility of mobile agents into dynamical networks, and then investigate packet forwarding strategy on the heterogeneous dynamical networks. We find that the faster speed and the higher proportion of high-speed agents can enhance the network throughput and reduce the mean traveling time in random forwarding. A hierarchical structure in the dependence of high-speed is observed: the network throughput remains unchanged at small and large high-speed value. It is also interesting to find that a slightly preferential forwarding to high-speed agents can maximize the network capacity. Through theoretical analysis and numerical simulations, we show that the optimal forwarding ratio stems from the local structural heterogeneity of low-speed agents.
Human opinion dynamics: An inspiration to solve complex optimization problems
Kaur, Rishemjit; Kumar, Ritesh; Bhondekar, Amol P.; Kapur, Pawan
2013-01-01
Human interactions give rise to the formation of different kinds of opinions in a society. The study of formations and dynamics of opinions has been one of the most important areas in social physics. The opinion dynamics and associated social structure leads to decision making or so called opinion consensus. Opinion formation is a process of collective intelligence evolving from the integrative tendencies of social influence with the disintegrative effects of individualisation, and therefore could be exploited for developing search strategies. Here, we demonstrate that human opinion dynamics can be utilised to solve complex mathematical optimization problems. The results have been compared with a standard algorithm inspired from bird flocking behaviour and the comparison proves the efficacy of the proposed approach in general. Our investigation may open new avenues towards understanding the collective decision making. PMID:24141795
Optimal flow rates and well locations for soil vapor extraction design
NASA Astrophysics Data System (ADS)
Sawyer, Charles S.; Kamakoti, Madhavi
1998-07-01
A mixed-integer programming model to determine the optimum number of wells, their locations and pumping rates for soil vapor extraction (SVE) is developed by coupling an air flow simulation model (AIR3D) to the GAMS optimization software. The model was tested for sensitivity of the vertical discretization of the domain, the number of potential well locations, the number of constraints, and the screen length of the wells. It was shown that these variables affected the optimal solution. It was also shown that the installation costs of the wells in the model influenced the optimal design. This was demonstrated by comparing the results of the mixed-integer programming model to a linear programming model in which the installation costs of the wells were neglected. The mixed-integer programming model could be useful in the design process in cases of short remediation times when the installation costs of wells could be significant. Numerous test cases with results are presented to demonstrate the applicability and usefulness of the model.
Optimized Uncertainty Quantification Algorithm Within a Dynamic Event Tree Framework
J. W. Nielsen; Akira Tokuhiro; Robert Hiromoto
2014-06-01
Methods for developing Phenomenological Identification and Ranking Tables (PIRT) for nuclear power plants have been a useful tool in providing insight into modelling aspects that are important to safety. These methods have involved expert knowledge with regards to reactor plant transients and thermal-hydraulic codes to identify are of highest importance. Quantified PIRT provides for rigorous method for quantifying the phenomena that can have the greatest impact. The transients that are evaluated and the timing of those events are typically developed in collaboration with the Probabilistic Risk Analysis. Though quite effective in evaluating risk, traditional PRA methods lack the capability to evaluate complex dynamic systems where end states may vary as a function of transition time from physical state to physical state . Dynamic PRA (DPRA) methods provide a more rigorous analysis of complex dynamic systems. A limitation of DPRA is its potential for state or combinatorial explosion that grows as a function of the number of components; as well as, the sampling of transition times from state-to-state of the entire system. This paper presents a method for performing QPIRT within a dynamic event tree framework such that timing events which result in the highest probabilities of failure are captured and a QPIRT is performed simultaneously while performing a discrete dynamic event tree evaluation. The resulting simulation results in a formal QPIRT for each end state. The use of dynamic event trees results in state explosion as the number of possible component states increases. This paper utilizes a branch and bound algorithm to optimize the solution of the dynamic event trees. The paper summarizes the methods used to implement the branch-and-bound algorithm in solving the discrete dynamic event trees.
Optimizing spread dynamics on graphs by message passing
NASA Astrophysics Data System (ADS)
Altarelli, F.; Braunstein, A.; Dall'Asta, L.; Zecchina, R.
2013-09-01
Cascade processes are responsible for many important phenomena in natural and social sciences. Simple models of irreversible dynamics on graphs, in which nodes activate depending on the state of their neighbors, have been successfully applied to describe cascades in a large variety of contexts. Over the past decades, much effort has been devoted to understanding the typical behavior of the cascades arising from initial conditions extracted at random from some given ensemble. However, the problem of optimizing the trajectory of the system, i.e. of identifying appropriate initial conditions to maximize (or minimize) the final number of active nodes, is still considered to be practically intractable, with the only exception being models that satisfy a sort of diminishing returns property called submodularity. Submodular models can be approximately solved by means of greedy strategies, but by definition they lack cooperative characteristics which are fundamental in many real systems. Here we introduce an efficient algorithm based on statistical physics for the optimization of trajectories in cascade processes on graphs. We show that for a wide class of irreversible dynamics, even in the absence of submodularity, the spread optimization problem can be solved efficiently on large networks. Analytic and algorithmic results on random graphs are complemented by the solution of the spread maximization problem on a real-world network (the Epinions consumer reviews network).
Dynamic optimization of bioprocesses: efficient and robust numerical strategies.
Banga, Julio R; Balsa-Canto, Eva; Moles, Carmen G; Alonso, Antonio A
2005-06-29
The dynamic optimization (open loop optimal control) of non-linear bioprocesses is considered in this contribution. These processes can be described by sets of non-linear differential and algebraic equations (DAEs), usually subject to constraints in the state and control variables. A review of the available solution techniques for this class of problems is presented, highlighting the numerical difficulties arising from the non-linear, constrained and often discontinuous nature of these systems. In order to surmount these difficulties, we present several alternative stochastic and hybrid techniques based on the control vector parameterization (CVP) approach. The CVP approach is a direct method which transforms the original problem into a non-linear programming (NLP) problem, which must be solved by a suitable (efficient and robust) solver. In particular, a hybrid technique uses a first global optimization phase followed by a fast second phase based on a local deterministic method, so it can handle the nonconvexity of many of these NLPs. The efficiency and robustness of these techniques is illustrated by solving several challenging case studies regarding the optimal control of fed-batch bioreactors and other bioprocesses. In order to fairly evaluate their advantages, a careful and critical comparison with several other direct approaches is provided. The results indicate that the two-phase hybrid approach presents the best compromise between robustness and efficiency. PMID:15888349
Active flutter control using discrete optimal constrained dynamic compensators
NASA Technical Reports Server (NTRS)
Broussard, J. R.; Halyo, N.
1983-01-01
A method for synthesizing digital active flutter suppression controllers using the concept of optimal output feedback is presented. A recently developd convergent algorithm is employed to determine constrained control law parameters that minimize an infinite-time discrete quadratic performance index. Low-order compensator dynamics are included in the control law and the compensator parameters are computed along with the output feedback gain as part of the optimization process. An input noise adjustment procedure is used to improve the stability margins of the digital active flutter controller. Results from investigations into sample rate variation, prefilter pole variation, and effects of varying flight condtions are discussed. The study indicates that a digital control law which accommodates computation delay can stabilize the wing with reasonable rms performance and adequate stability margins.
Optimal Dynamic Regimes: Presenting a Case for Predictive Inference
Arjas, Elja; Saarela, Olli
2010-01-01
Dynamic treatment regime is a decision rule in which the choice of the treatment of an individual at any given time can depend on the known past history of that individual, including baseline covariates, earlier treatments, and their measured responses. In this paper we argue that finding an optimal regime can, at least in moderately simple cases, be accomplished by a straightforward application of nonparametric Bayesian modeling and predictive inference. As an illustration we consider an inference problem in a subset of the Multicenter AIDS Cohort Study (MACS) data set, studying the effect of AZT initiation on future CD4-cell counts during a 12-month follow-up. PMID:20648215
Khawaja, Sajid Gul; Mushtaq, Mian Hamza; Khan, Shoab A; Akram, M Usman; Jamal, Habib Ullah
2015-01-01
With the increase of transistors' density, popularity of System on Chip (SoC) has increased exponentially. As a communication module for SoC, Network on Chip (NoC) framework has been adapted as its backbone. In this paper, we propose a methodology for designing area-optimized application specific NoC while providing hard Quality of Service (QoS) guarantees for real time flows. The novelty of the proposed system lies in derivation of a Mixed Integer Linear Programming model which is then used to generate a resource optimal Network on Chip (NoC) topology and architecture while considering traffic and QoS requirements. We also present the micro-architectural design features used for enabling traffic and latency guarantees and discuss how the solution adapts for dynamic variations in the application traffic. The paper highlights the effectiveness of proposed method by generating resource efficient NoC solutions for both industrial and benchmark applications. The area-optimized results are generated in few seconds by proposed technique, without resorting to heuristics, even for an application with 48 traffic flows. PMID:25898016
Khawaja, Sajid Gul; Mushtaq, Mian Hamza; Khan, Shoab A.; Akram, M. Usman; Jamal, Habib ullah
2015-01-01
With the increase of transistors' density, popularity of System on Chip (SoC) has increased exponentially. As a communication module for SoC, Network on Chip (NoC) framework has been adapted as its backbone. In this paper, we propose a methodology for designing area-optimized application specific NoC while providing hard Quality of Service (QoS) guarantees for real time flows. The novelty of the proposed system lies in derivation of a Mixed Integer Linear Programming model which is then used to generate a resource optimal Network on Chip (NoC) topology and architecture while considering traffic and QoS requirements. We also present the micro-architectural design features used for enabling traffic and latency guarantees and discuss how the solution adapts for dynamic variations in the application traffic. The paper highlights the effectiveness of proposed method by generating resource efficient NoC solutions for both industrial and benchmark applications. The area-optimized results are generated in few seconds by proposed technique, without resorting to heuristics, even for an application with 48 traffic flows. PMID:25898016
Optimization of phonon dynamics protocols in ion traps
NASA Astrophysics Data System (ADS)
Dutta, T.; Mukherjee, M.; Sengupta, K.
2016-03-01
We develop a theory to address the non equilibrium dynamics of phonons in a one-dimensional finite size trapped ion system for non linear ramp and periodic protocols. Our analysis, which is based on our earlier proposal of dynamics-induced cooling and entanglement generation between phonons in these systems when subjected to a linear ramp protocol [1], identifies the optimal protocol within the above-mentioned classes, which minimizes both the cooling and entanglement generation time. We also introduce single-/two-site addressing to achieve cooling/entanglement, which is expected to lead to simpler implementation of these protocols. Finally, we discuss the effect of noise due to the fluctuation of the intensity of the laser used to generate the trap on entanglement generation. We also discuss realistic experimental setups that may serve as test beds for our theory.
Dynamic Simulation and Optimization of Nuclear Hydrogen Production Systems
Paul I. Barton; Mujid S. Kaximi; Georgios Bollas; Patricio Ramirez Munoz
2009-07-31
This project is part of a research effort to design a hydrogen plant and its interface with a nuclear reactor. This project developed a dynamic modeling, simulation and optimization environment for nuclear hydrogen production systems. A hybrid discrete/continuous model captures both the continuous dynamics of the nuclear plant, the hydrogen plant, and their interface, along with discrete events such as major upsets. This hybrid model makes us of accurate thermodynamic sub-models for the description of phase and reaction equilibria in the thermochemical reactor. Use of the detailed thermodynamic models will allow researchers to examine the process in detail and have confidence in the accurary of the property package they use.
Optimized dynamical decoupling for power-law noise spectra
Pasini, S.; Uhrig, G. S.
2010-01-15
We analyze the suppression of decoherence by means of dynamical decoupling in the pure-dephasing spin-boson model for baths with power law spectra. The sequence of ideal pi pulses is optimized according to the power of the bath. We expand the decoherence function and separate the canceling divergences from the relevant terms. The proposed sequence is chosen to be the one minimizing the decoherence function. By construction, it provides the best performance. We analytically derive the conditions that must be satisfied. The resulting equations are solved numerically. The solutions are very close to the Carr-Purcell-Meiboom-Gill sequence for a soft cutoff of the bath while they approach the Uhrig dynamical-decoupling sequence as the cutoff becomes harder.
Clustering molecular dynamics trajectories for optimizing docking experiments.
De Paris, Renata; Quevedo, Christian V; Ruiz, Duncan D; Norberto de Souza, Osmar; Barros, Rodrigo C
2015-01-01
Molecular dynamics simulations of protein receptors have become an attractive tool for rational drug discovery. However, the high computational cost of employing molecular dynamics trajectories in virtual screening of large repositories threats the feasibility of this task. Computational intelligence techniques have been applied in this context, with the ultimate goal of reducing the overall computational cost so the task can become feasible. Particularly, clustering algorithms have been widely used as a means to reduce the dimensionality of molecular dynamics trajectories. In this paper, we develop a novel methodology for clustering entire trajectories using structural features from the substrate-binding cavity of the receptor in order to optimize docking experiments on a cloud-based environment. The resulting partition was selected based on three clustering validity criteria, and it was further validated by analyzing the interactions between 20 ligands and a fully flexible receptor (FFR) model containing a 20 ns molecular dynamics simulation trajectory. Our proposed methodology shows that taking into account features of the substrate-binding cavity as input for the k-means algorithm is a promising technique for accurately selecting ensembles of representative structures tailored to a specific ligand. PMID:25873944
Clustering Molecular Dynamics Trajectories for Optimizing Docking Experiments
De Paris, Renata; Quevedo, Christian V.; Ruiz, Duncan D.; Norberto de Souza, Osmar; Barros, Rodrigo C.
2015-01-01
Molecular dynamics simulations of protein receptors have become an attractive tool for rational drug discovery. However, the high computational cost of employing molecular dynamics trajectories in virtual screening of large repositories threats the feasibility of this task. Computational intelligence techniques have been applied in this context, with the ultimate goal of reducing the overall computational cost so the task can become feasible. Particularly, clustering algorithms have been widely used as a means to reduce the dimensionality of molecular dynamics trajectories. In this paper, we develop a novel methodology for clustering entire trajectories using structural features from the substrate-binding cavity of the receptor in order to optimize docking experiments on a cloud-based environment. The resulting partition was selected based on three clustering validity criteria, and it was further validated by analyzing the interactions between 20 ligands and a fully flexible receptor (FFR) model containing a 20 ns molecular dynamics simulation trajectory. Our proposed methodology shows that taking into account features of the substrate-binding cavity as input for the k-means algorithm is a promising technique for accurately selecting ensembles of representative structures tailored to a specific ligand. PMID:25873944
Aircraft path planning for optimal imaging using dynamic cost functions
NASA Astrophysics Data System (ADS)
Christie, Gordon; Chaudhry, Haseeb; Kochersberger, Kevin
2015-05-01
Unmanned aircraft development has accelerated with recent technological improvements in sensing and communications, which has resulted in an "applications lag" for how these aircraft can best be utilized. The aircraft are becoming smaller, more maneuverable and have longer endurance to perform sensing and sampling missions, but operating them aggressively to exploit these capabilities has not been a primary focus in unmanned systems development. This paper addresses a means of aerial vehicle path planning to provide a realistic optimal path in acquiring imagery for structure from motion (SfM) reconstructions and performing radiation surveys. This method will allow SfM reconstructions to occur accurately and with minimal flight time so that the reconstructions can be executed efficiently. An assumption is made that we have 3D point cloud data available prior to the flight. A discrete set of scan lines are proposed for the given area that are scored based on visibility of the scene. Our approach finds a time-efficient path and calculates trajectories between scan lines and over obstacles encountered along those scan lines. Aircraft dynamics are incorporated into the path planning algorithm as dynamic cost functions to create optimal imaging paths in minimum time. Simulations of the path planning algorithm are shown for an urban environment. We also present our approach for image-based terrain mapping, which is able to efficiently perform a 3D reconstruction of a large area without the use of GPS data.
Data-driven optimization of dynamic reconfigurable systems of systems.
Tucker, Conrad S.; Eddy, John P.
2010-11-01
This report documents the results of a Strategic Partnership (aka University Collaboration) LDRD program between Sandia National Laboratories and the University of Illinois at Urbana-Champagne. The project is titled 'Data-Driven Optimization of Dynamic Reconfigurable Systems of Systems' and was conducted during FY 2009 and FY 2010. The purpose of this study was to determine and implement ways to incorporate real-time data mining and information discovery into existing Systems of Systems (SoS) modeling capabilities. Current SoS modeling is typically conducted in an iterative manner in which replications are carried out in order to quantify variation in the simulation results. The expense of many replications for large simulations, especially when considering the need for optimization, sensitivity analysis, and uncertainty quantification, can be prohibitive. In addition, extracting useful information from the resulting large datasets is a challenging task. This work demonstrates methods of identifying trends and other forms of information in datasets that can be used on a wide range of applications such as quantifying the strength of various inputs on outputs, identifying the sources of variation in the simulation, and potentially steering an optimization process for improved efficiency.
Optimizing Groundwater Remediation Designs Under Uncertainty Using Dynamic Surrogate Models
NASA Astrophysics Data System (ADS)
Yan, S.; Minsker, B.
2005-12-01
Computational cost is a critical issue for large-scale water resource optimization problems that often involve time-consuming simulation models. This issue is compounded when optimizing under uncertainty, since Monte Carlo simulations are often required to evaluate objective function values over multiple parameter realizations. In order to improve computational efficiency, we propose a soft computing approach, in which the time-consuming numerical models are approximated and replaced by dynamic surrogates embedded within a noisy genetic algorithm (GA) optimization framework. The surrogates are trained to predict the distribution of the objectives online, using Monte Carlo simulation results created during the GA run. The surrogates are then adaptively updated to improve their prediction performance and correct the GA's convergence as the search progresses. Latin Hypercube sampling method is used to efficiently sample parameters for the Monte Carlo simulation and the sampling results are archived so that the estimation of the objective function distributions is progressively improved. The GA is modified to incorporate hypothesis tests to produce reliable solutions. The method is applied to two groundwater remediation design case studies, where the primary source of uncertainty stems from hydraulic conductivity values in the aquifers. Our preliminary results show that the technique can lead to reliable and cost-effective solutions with significantly less computational effort.
Assessment of optimally filtered recent geodetic mean dynamic topographies
NASA Astrophysics Data System (ADS)
Siegismund, F.
2013-01-01