Simulation and optimization tools in paper machine concept design

The last decade has shown that the global paper industry needs new processes and products in order to reassert its position in the industry. As the paper markets in Western Europe and North America have stabilized, the competition has tightened. Along with the development of more cost-effective processes and products, new process design methods are also required to break the old molds and create new ideas. This thesis discusses the development of a process design methodology based on simulation and optimization methods. A bi-level optimization problem and a solution procedure for it are formulated and illustrated. Computational models and simulation are used to illustrate the phenomena inside a real process and mathematical optimization is exploited to find out the best process structures and control principles for the process. Dynamic process models are used inside the bi-level optimization problem, which is assumed to be dynamic and multiobjective due to the nature of papermaking processes. The numerical experiments show that the bi-level optimization approach is useful for different kinds of problems related to process design and optimization. Here, the design methodology is applied to a constrained process area of a papermaking line. However, the same methodology is applicable to all types of industrial processes, e.g., the design of biorefiners, because the methodology is totally generalized and can be easily modified.

Object-oriented Programming Applied to the Development of Structural Analysis and Optimization Software

This paper presents the development of a two-dimensional interactive software environment for structural analysis and optimization based on object-oriented programming using the C++ language. The main feature of the software is the effective integration of several computational tools into graphical user interfaces implemented in the Windows-98 and Windows-NT operating systems. The interfaces simplify data specification in the simulation and optimization of two-dimensional linear elastic problems. NURBS have been used in the software modules to represent geometric and graphical data. Extensions to the analysis of three-dimensional problems have been implemented and are also discussed in this paper.

Analysis and Optimization of Machining Process Using Evolutionary Algorithms

To ensure quality of machined products at minimum machining costs and maximum machining effectiveness, it is very important to select optimum parameters when metal cutting machine tools are employed. Traditionally, the experience of the operator plays a major role in the selection of optimum metal cutting conditions. However, attaining optimum values each time by even a skilled operator is difficult. The non-linear nature of the machining process has compelled engineers to search for more effective methods to attain optimization. The design objective preceding most engineering design activities is simply to minimize the cost of production or to maximize the production efficiency. The main aim of research work reported here is to build robust optimization algorithms by exploiting ideas that nature has to offer from its backyard and using it to solve real world optimization problems in manufacturing processes.In this thesis, after conducting an exhaustive literature review, several optimization techniques used in various manufacturing processes have been identified. The selection of optimal cutting parameters, like depth of cut, feed and speed is a very important issue for every machining process. Experiments have been designed using Taguchi technique and dry turning of SS420 has been performed on Kirlosker turn master 35 lathe. Analysis using S/N and ANOVA were performed to find the optimum level and percentage of contribution of each parameter. By using S/N analysis the optimum machining parameters from the experimentation is obtained.Optimization algorithms begin with one or more design solutions supplied by the user and then iteratively check new design solutions, relative search spaces in order to achieve the true optimum solution. A mathematical model has been developed using response surface analysis for surface roughness and the model was validated using published results from literature.Methodologies in optimization such as Simulated annealing (SA), Particle Swarm Optimization (PSO), Conventional Genetic Algorithm (CGA) and Improved Genetic Algorithm (IGA) are applied to optimize machining parameters while dry turning of SS420 material. All the above algorithms were tested for their efficiency, robustness and accuracy and observe how they often outperform conventional optimization method applied to difficult real world problems. The SA, PSO, CGA and IGA codes were developed using MATLAB. For each evolutionary algorithmic method, optimum cutting conditions are provided to achieve better surface finish.The computational results using SA clearly demonstrated that the proposed solution procedure is quite capable in solving such complicated problems effectively and efficiently. Particle Swarm Optimization (PSO) is a relatively recent heuristic search method whose mechanics are inspired by the swarming or collaborative behavior of biological populations. From the results it has been observed that PSO provides better results and also more computationally efficient.Based on the results obtained using CGA and IGA for the optimization of machining process, the proposed IGA provides better results than the conventional GA. The improved genetic algorithm incorporating a stochastic crossover technique and an artificial initial population scheme is developed to provide a faster search mechanism. Finally, a comparison among these algorithms were made for the specific example of dry turning of SS 420 material and arriving at optimum machining parameters of feed, cutting speed, depth of cut and tool nose radius for minimum surface roughness as the criterion. To summarize, the research work fills in conspicuous gaps between research prototypes and industry requirements, by simulating evolutionary procedures seen in nature that optimize its own systems.

A Two-level Prediction Model for Deep Reactive Ion Etch (DRIE)

We contribute a quantitative and systematic model to capture etch non-uniformity in deep reactive ion etch of microelectromechanical systems (MEMS) devices. Deep reactive ion etch is commonly used in MEMS fabrication where high-aspect ratio features are to be produced in silicon. It is typical for many supposedly identical devices, perhaps of diameter 10 mm, to be etched simultaneously into one silicon wafer of diameter 150 mm. Etch non-uniformity depends on uneven distributions of ion and neutral species at the wafer level, and on local consumption of those species at the device, or die, level. An ion–neutral synergism model is constructed from data obtained from etching several layouts of differing pattern opening densities. Such a model is used to predict wafer-level variation with an r.m.s. error below 3%. This model is combined with a die-level model, which we have reported previously, on a MEMS layout. The two-level model is shown to enable prediction of both within-die and wafer-scale etch rate variation for arbitrary wafer loadings.

Control and optimization of an SBR for nitrogen removal: from model calibration to plant operation

En aquesta tesis s'ha desenvolupat un sistema de control capaç d'optimitzar el funcionament dels Reactors Discontinus Seqüencials dins el camp de l'eliminació de matèria orgànica i nitrogen de les aigües residuals. El sistema de control permet ajustar en línia la durada de les etapes de reacció a partir de mesures directes o indirectes de sondes. En una primera etapa de la tesis s'ha estudiat la calibració de models matemàtics que permeten realitzar fàcilment provatures de diferents estratègies de control. A partir de l'anàlisis de dades històriques s'han plantejat diferents opcions per controlar l'SBR i les més convenients s'han provat mitjançant simulació. Després d'assegurar l'èxit de l'estratègia de control mitjançant simulacions s'ha implementat en una planta semi-industrial. Finalment es planteja l'estructura d'uns sistema supervisor encarregat de controlar el funcionament de l'SBR no només a nivell de fases sinó també a nivell cicle.

Minimum aberration construction results for nonregular two-level fractional factorial designs

Nonregular two-level fractional factorial designs are designs which cannot be specified in terms of a set of defining contrasts. The aliasing properties of nonregular designs can be compared by using a generalisation of the minimum aberration criterion called minimum G2-aberration.Until now, the only nontrivial designs that are known to have minimum G2-aberration are designs for n runs and m n–5 factors. In this paper, a number of construction results are presented which allow minimum G2-aberration designs to be found for many of the cases with n = 16, 24, 32, 48, 64 and 96 runs and m n/2–2 factors.

Local adaptation in moss life histories: population-level variation and a reciprocal transplant experiment

We examine the extent of population-level differentiation in life history traits of Pogonatum aloides, Polytrichum commune and Polytrichum juniperinum (Polytrichaceae) between upland and lowland localities within Britain. Reciprocal transplant studies are used to estimate the relative importance of genetic versus environmental effects on observed differences. We demonstrate significant life history differentiation between moss populations, and show that at least some of these are genetically determined, although environment and phenotypic plasticity are also significant components of the observed variation. The transplant experiments indicate divergence among populations in plasticity of male reproductive effort and of investment in vegetative shoots by females. Two tradeoffs are identified; one between the number and the size of spores, and the second between reproduction by spores versus vegetative reproduction. The patterns of life history variation observed between populations of Polytrichum juniperinum are consistent with selection along these implied tradeoff curves, and we propose that they reflect selective pressures arising from the spatial and demographic distribution of mortality at upland versus lowland sites. The results underscore the need for more studies of intra-specific life history variation in mosses.

Ocean heat uptake and its consequences for the magnitude of sea level rise and climate change

Under increasing greenhouse gas concentrations, ocean heat uptake moderates the rate of climate change, and thermal expansion makes a substantial contribution to sea level rise. In this paper we quantify the differences in projections among atmosphere-ocean general circulation models of the Coupled Model Intercomparison Project in terms of transient climate response, ocean heat uptake efficiency and expansion efficiency of heat. The CMIP3 and CMIP5 ensembles have statistically indistinguishable distributions in these parameters. The ocean heat uptake efficiency varies by a factor of two across the models, explaining about 50% of the spread in ocean heat uptake in CMIP5 models with CO2 increasing at 1%/year. It correlates with the ocean global-mean vertical profiles both of temperature and of temperature change, and comparison with observations suggests the models may overestimate ocean heat uptake and underestimate surface warming, because their stratification is too weak. The models agree on the location of maxima of shallow ocean heat uptake (above 700 m) in the Southern Ocean and the North Atlantic, and on deep ocean heat uptake (below 2000 m) in areas of the Southern Ocean, in some places amounting to 40% of the top-to-bottom integral in the CMIP3 SRES A1B scenario. The Southern Ocean dominates global ocean heat uptake; consequently the eddy-induced thickness diffusivity parameter, which is particularly influential in the Southern Ocean, correlates with the ocean heat uptake efficiency. The thermal expansion produced by ocean heat uptake is 0.12 m YJ−1, with an uncertainty of about 10% (1 YJ = 1024 J).

Climatic simulations of the eastern Andes low-level jet and its dependency on convective parameterizations

The South American low level jet (SALLJ) of the Eastern Andes is investigated with Regional Climate Model version 3 (RegCM3) simulations during the 2002-2003 austral summer using two convective parameterizations (Grell and Emanuel). The simulated SALLJ is compared with the special observations of SALLJEX (SALLJ Experiment). Both the Grell and Emanuel schemes adequately simulate the low level flow over South America. However, there are some intensity differences. Due to the larger (smaller) convective activity, the Emanuel (Grell) scheme simulates more intense (weaker) low level wind than analysis in the tropics and subtropics. The objectives criteria of Sugahara (SJ) and Bonner (BJ) were used for LLJ identification. When applied to the observations, both criteria suggest a larger frequency of the SALLJ in Santa Cruz, followed by Mariscal, Trinidad and Asuncin. In Mariscal and Asuncin, the diurnal cycle indicates that SJ occurs mainly at 12 UTCs (morning), while the BJ criterion presents the SALLJ as more homogenously distributed. The concentration into two of the four-times-a-day observations does not allow conclusions about the diurnal cycle in Santa Cruz and Trinidad. The simulated wind profiles result in a lower than observed frequency of SALLJ using both the SJ and BJ criteria, with fewer events obtained with the BJ. Due to the stronger simulated winds, the Emanuel scheme produces an equal or greater relative frequency of SALLJ than the Grell scheme. However, the Grell scheme using the SJ criterion simulates the SALLJ diurnal cycle closer to the observed one. Although some discrepancies between observed and simulated mean vertical profiles of the horizontal wind are noted, there is large agreement between the composites of the vertical structure of the SALLJ, especially when the SJ criterion is used with the Grell scheme. On an intraseasonal scale, a larger southward displacement of SALLJ in February and December when compared with January has been noted. The Grell and Emanuel schemes simulated this observed oscillation in the low-level flow. However, the spatial pattern and intensity of rainfall and circulation anomalies simulated by the Grell scheme are closer to the analyses than those obtained with the Emanuel scheme.

Four-level systems and a universal quantum gate

We discuss the possibility of implementing a universal quantum XOR gate by using two coupled quantum dots subject to external magnetic fields that are parallel and slightly different. We consider this system in two different field configurations. In the first case, parallel external fields with the intensity difference at each spin being proportional to the time-dependent interaction between the spins. A general exact solution describing this system is presented and analyzed to adjust field parameters. Then we consider parallel fields with intensity difference at each spin being constant and the interaction between the spins switching on and off adiabatically. In both cases we adjust characteristics of the external fields (their intensities and duration) in order to have the parallel pulse adequate for constructing the XOR gate. In order to provide a complete theoretical description of all the cases, we derive relations between the spin interaction, the inter-dot distance, and the external field. (C) 2008 WILEYNCH Verlag GmbH & Co. KGaA. Weinheim.

Optical Characterization and Optimization of Display Components : Some Applications to Liquid-Crystal-Based and Electrochromics-Based Devices

This dissertation is focused on theoretical and experimental studies of optical properties of materials and multilayer structures composing liquid crystal displays (LCDs) and electrochromic (EC) devices. By applying spectroscopic ellipsometry, we have determined the optical constants of thin films of electrochromic tungsten oxide (WOx) and nickel oxide (NiOy), the films’ thickness and roughness. These films, which were obtained at spattering conditions possess high transmittance that is important for achieving good visibility and high contrast in an EC device. Another application of the general spectroscopic ellipsometry relates to the study of a photo-alignment layer of a mixture of azo-dyes SD-1 and SDA-2. We have found the optical constants of this mixture before and after illuminating it by polarized UV light. The results obtained confirm the diffusion model to explain the formation of the photo-induced order in azo-dye films. We have developed new techniques for fast characterization of twisted nematic LC cells in transmissive and reflective modes. Our techniques are based on the characteristics functions that we have introduced for determination of parameters of non-uniform birefringent media. These characteristic functions are found by simple procedures and can be utilised for simultaneous determination of retardation, its wavelength dispersion, and twist angle, as well as for solving associated optimization problems. Cholesteric LCD that possesses some unique properties, such as bistability and good selective scattering, however, has a disadvantage – relatively high driving voltage (tens of volts). The way we propose to reduce the driving voltage consists of applying a stack of thin (~1µm) LC layers. We have studied the ability of a layer of a surface stabilized ferroelectric liquid crystal coupled with several retardation plates for birefringent color generation. We have demonstrated that in order to accomplish good color characteristics and high brightness of the display, one or two retardation plates are sufficient.

HISTOCHEMICAL RESPONSES of RATS EXERCISED IN TWO WEEKLY FREQUENCIES and INGESTING STANDARD OR HYPERCALORIC DIET

This study investigated if overfed rats present morphological and histochemical muscle adaptation similar to normally fed, both submitted to two different weekly frequencies of training. Thirty male Wistar rats were fed either with standard chow (SCO) or with hypercaloric diet (HCO). They were subdivided into six subgroups: sedentary (SCO and HCO), trained twice/week (SC2 and HC2) and trained five times/week (SC5 and HC5). The trained groups swam 60 min/day, during 10 weeks. Twenty four hours after the last training, samples of Gastrocnemius were excised and stained with HE, NADH-TR and m-ATPase, and the capillary density was calculated. Total heart mass (HM) and the mass of atrium (AM), left (LV) and right (RV) ventricles were excised and weighted. The comparisons were made by ANOVA and by Covariance analysis, adjusting the variables by body weight. The results showed that the HCO achieved higher BM, however, absolute HM did not differ post training. Irrespective of the diet, rats that were trained twice a week presented significantly greater increase in the AM. In general, the SC5 and HC5 groups showed higher HM, LV, RV, proportion of oxidative fibres and capillary density, compared to the sedentary and twice week trained groups. A higher proportion of injuries (splitting) was noted in the HC2 and HC5 compared to SC2 and SC5. These results indicate that the frequency of training influenced the skeletal and heart adaptation and larger changes were observed in the 5x/week group, which ingested the standard diet. The 5x/week training groups also presented large amount of muscle fibres damage.

Deviation from the sech2 superradiant emission law in a two-level atomic system

The atomic superradiant emission is treated in the single-particle mean-field approximation. A single-particle Hamiltonian, which represents a dressed two-level atom in a radiation field, can be obtained and it is verified that it describes the transient regime of the emission process. While the line-shape emission for a bare atom follows the sech2 law, for the dressed atom the line shape deviates appreciably from this law and it is verified that the deviation depends crucially on the ratio of the dynamic frequency shift to the transition frequency. This kind of deviation is observed in experimental results. © 1990 The American Physical Society.

Electronic correlations in a two-level dimer

The exact solution for the full electronic Hamiltonian for a two-level dimer is obtained. The parameter constellation (roughly 20) is reparametrized via orthogonal Gaussian atomic orbitals, yielding a five-parameter model. With the dimer embedded in a thermal bath, the specific heat and several temperature-dependent dynamical susceptibilities are computed. © 2001 Elsevier Science B.V. All rights reserved.

Two-level network design with intermediate facilities: An application to electrical distribution systems

We consider the two-level network design problem with intermediate facilities. This problem consists of designing a minimum cost network respecting some requirements, usually described in terms of the network topology or in terms of a desired flow of commodities between source and destination vertices. Each selected link must receive one of two types of edge facilities and the connection of different edge facilities requires a costly and capacitated vertex facility. We propose a hybrid decomposition approach which heuristically obtains tentative solutions for the vertex facilities number and location and use these solutions to limit the computational burden of a branch-and-cut algorithm. We test our method on instances of the power system secondary distribution network design problem. The results show that the method is efficient both in terms of solution quality and computational times. © 2010 Elsevier Ltd.