Optimal control and operation of drum granulation processes

Process optimisation and optimal control of batch and continuous drum granulation processes are studied in this paper. The main focus of the current research has been: (i) construction of optimisation and control relevant, population balance models through the incorporation of moisture content, drum rotation rate and bed depth into the coalescence kernels; (ii) investigation of optimal operational conditions using constrained optimisation techniques; (iii) development of optimal control algorithms based on discretized population balance equations; and (iv) comprehensive simulation studies on optimal control of both batch and continuous granulation processes. The objective of steady state optimisation is to minimise the recycle rate with minimum cost for continuous processes. It has been identified that the drum rotation-rate, bed depth (material charge), and moisture content of solids are practical decision (design) parameters for system optimisation. The objective for the optimal control of batch granulation processes is to maximize the mass of product-sized particles with minimum time and binder consumption. The objective for the optimal control of the continuous process is to drive the process from one steady state to another in a minimum time with minimum binder consumption, which is also known as the state-driving problem. It has been known for some time that the binder spray-rate is the most effective control (manipulative) variable. Although other possible manipulative variables, such as feed flow-rate and additional powder flow-rate have been investigated in the complete research project, only the single input problem with the binder spray rate as the manipulative variable is addressed in the paper to demonstrate the methodology. It can be shown from simulation results that the proposed models are suitable for control and optimisation studies, and the optimisation algorithms connected with either steady state or dynamic models are successful for the determination of optimal operational conditions and dynamic trajectories with good convergence properties. (c) 2005 Elsevier Ltd. All rights reserved.

Optimal control, geometry, and quantum computing

We prove upper and lower bounds relating the quantum gate complexity of a unitary operation, U, to the optimal control cost associated to the synthesis of U. These bounds apply for any optimal control problem, and can be used to show that the quantum gate complexity is essentially equivalent to the optimal control cost for a wide range of problems, including time-optimal control and finding minimal distances on certain Riemannian, sub-Riemannian, and Finslerian manifolds. These results generalize the results of [Nielsen, Dowling, Gu, and Doherty, Science 311, 1133 (2006)], which showed that the gate complexity can be related to distances on a Riemannian manifold.

A new approach to control of MIMO processes with static nonlinearities using an extended IMC framework

In this paper, a new control design method is proposed for stable processes which can be described using Hammerstein-Wiener models. The internal model control (IMC) framework is extended to accommodate multiple IMC controllers, one for each subsystem. The concept of passive systems is used to construct the IMC controllers which approximate the inverses of the subsystems to achieve dynamic control performance. The Passivity Theorem is used to ensure the closed-loop stability. (c) 2005 Elsevier Ltd. All rights reserved.

A D-optimal designed population pharmacokinetic study of itraconazole capsules and solution in adults with cystic fibrosis

Background: Oral itraconazole (ITRA) is used for the treatment of allergic bronchopulmonary aspergillosis in patients with cystic fibrosis (CF) because of its antifungal activity against Aspergillus species. ITRA has an active hydroxy-metabolite (OH-ITRA) which has similar antifungal activity. ITRA is a highly lipophilic drug which is available in two different oral formulations, a capsule and an oral solution. It is reported that the oral solution has a 60% higher relative bioavailability. The influence of altered gastric physiology associated with CF on the pharmacokinetics (PK) of ITRA and its metabolite has not been previously evaluated. Objectives: 1) To estimate the population (pop) PK parameters for ITRA and its active metabolite OH-ITRA including relative bioavailability of the parent after administration of the parent by both capsule and solution and 2) to assess the performance of the optimal design. Methods: The study was a cross-over design in which 30 patients received the capsule on the first occasion and 3 days later the solution formulation. The design was constrained to have a maximum of 4 blood samples per occasion for estimation of the popPK of both ITRA and OH-ITRA. The sampling times for the population model were optimized previously using POPT v.2.0.[1] POPT is a series of applications that run under MATLAB and provide an evaluation of the information matrix for a nonlinear mixed effects model given a particular design. In addition it can be used to optimize the design based on evaluation of the determinant of the information matrix. The model details for the design were based on prior information obtained from the literature, which suggested that ITRA may have either linear or non-linear elimination. The optimal sampling times were evaluated to provide information for both competing models for the parent and metabolite and for both capsule and solution simultaneously. Blood samples were assayed by validated HPLC.[2] PopPK modelling was performed using FOCE with interaction under NONMEM, version 5 (level 1.1; GloboMax LLC, Hanover, MD, USA). The PK of ITRA and OH‑ITRA was modelled simultaneously using ADVAN 5. Subsequently three methods were assessed for modelling concentrations less than the LOD (limit of detection). These methods (corresponding to methods 5, 6 & 4 from Beal[3], respectively) were (a) where all values less than LOD were assigned to half of LOD, (b) where the closest missing value that is less than LOD was assigned to half the LOD and all previous (if during absorption) or subsequent (if during elimination) missing samples were deleted, and (c) where the contribution of the expectation of each missing concentration to the likelihood is estimated. The LOD was 0.04 mg/L. The final model evaluation was performed via bootstrap with re-sampling and a visual predictive check. The optimal design and the sampling windows of the study were evaluated for execution errors and for agreement between the observed and predicted standard errors. Dosing regimens were simulated for the capsules and the oral solution to assess their ability to achieve ITRA target trough concentration (Cmin,ss of 0.5-2 mg/L) or a combined Cmin,ss for ITRA and OH-ITRA above 1.5mg/L. Results and Discussion: A total of 241 blood samples were collected and analysed, 94% of them were taken within the defined optimal sampling windows, of which 31% where taken within 5 min of the exact optimal times. Forty six per cent of the ITRA values and 28% of the OH-ITRA values were below LOD. The entire profile after administration of the capsule for five patients was below LOD and therefore the data from this occasion was omitted from estimation. A 2-compartment model with 1st order absorption and elimination best described ITRA PK, with 1st order metabolism of the parent to OH-ITRA. For ITRA the clearance (ClItra/F) was 31.5 L/h; apparent volumes of central and peripheral compartments were 56.7 L and 2090 L, respectively. Absorption rate constants for capsule (kacap) and solution (kasol) were 0.0315 h-1 and 0.125 h-1, respectively. Comparative bioavailability of the capsule was 0.82. There was no evidence of nonlinearity in the popPK of ITRA. No screened covariate significantly improved the fit to the data. The results of the parameter estimates from the final model were comparable between the different methods for accounting for missing data, (M4,5,6)[3] and provided similar parameter estimates. The prospective application of an optimal design was found to be successful. Due to the sampling windows, most of the samples could be collected within the daily hospital routine, but still at times that were near optimal for estimating the popPK parameters. The final model was one of the potential competing models considered in the original design. The asymptotic standard errors provided by NONMEM for the final model and empirical values from bootstrap were similar in magnitude to those predicted from the Fisher Information matrix associated with the D-optimal design. Simulations from the final model showed that the current dosing regimen of 200 mg twice daily (bd) would provide a target Cmin,ss (0.5-2 mg/L) for only 35% of patients when administered as the solution and 31% when administered as capsules. The optimal dosing schedule was 500mg bd for both formulations. The target success for this dosing regimen was 87% for the solution with an NNT=4 compared to capsules. This means, for every 4 patients treated with the solution one additional patient will achieve a target success compared to capsule but at an additional cost of AUD $220 per day. The therapeutic target however is still doubtful and potential risks of these dosing schedules need to be assessed on an individual basis. Conclusion: A model was developed which described the popPK of ITRA and its main active metabolite OH-ITRA in adult CF after administration of both capsule and solution. The relative bioavailability of ITRA from the capsule was 82% that of the solution, but considerably more variable. To incorporate missing data, using the simple Beal method 5 (using half LOD for all samples below LOD) provided comparable results to the more complex but theoretically better Beal method 4 (integration method). The optimal sparse design performed well for estimation of model parameters and provided a good fit to the data.

Double-strand break repair gene polymorphisms and risk of breast or ovarian cancer

Deficiencies in DNA repair have been hypothesized to increase cancer risk and excess cancer incidence is a feature of inherited diseases caused by defects in DNA damage recognition and repair. We investigated, using a case-control design, whether the double-strand break repair gene polymorphisms RAD51 5' untranslated region -135 G > C, XRCC2 R188H G > A, and XRCC3 T241M C > T were associated with risk of breast or ovarian cancer in Australian women. Sample sets included 1,456 breast cancer cases and 793 age-matched controls ages under 60 years of age, 549 incident ovarian cancer cases, and 335 controls of similar age distribution. For the total sample and the subsample of Caucasian women, there were no significant differences in genotype distribution between breast cancer cases and controls or between ovarian cancer cases and combined control groups. The crude odds ratios (OR) and 95% confidence intervals (95% CI) associated with the RAD51 GC/CC genotype frequency was OR, 1.10; 95% CI, 0.80-1.41 for breast cancer and OR, 1.22; 95% CI, 0.92-1.62 for ovarian cancer. Similarly, there were no increased risks associated with the XRCC2 GA/AA genotype (OR, 0.98; 95% CI, 0.76-1.26 for breast cancer and OR, 0.93; 95% CI, 0.69-1.25 for ovarian cancer) or the XRCC3 CT/TT genotype (OR, 0.92; 95% Cl, 0.77-1.10 for breast cancer and OR, 0.87; 95% CI, 0.71-1.08 for ovarian cancer). Results were little changed after adjustment for age and other measured risk factors. Although there was little statistical power to detect modest increases in risk for the homozygote variant genotypes, particularly for the rare RAD51 and XRCC2 variants, the data suggest that none of these variants play a major role in the etiology of breast or ovarian cancer.

Degrees of freedom and motor planning in purposive movement

This paper presents empirical evidence suggesting that healthy humans can perform a two degree of freedom visuo-motor pursuit tracking task with the same response time delay as a one degree of freedom task. In contrast, the time delay of the response is influenced markedly by the nature of the motor synergy required to produce it. We suggest a conceptual account of this evidence based on adaptive model theory, which combines theories of intermittency from psychology and adaptive optimal control from engineering. The intermittent response planning stage has a fixed period. It possesses multiple optimal trajectory generators such that multiple degrees of freedom can be planned concurrently, without requiring an increase in the planning period. In tasks which require unfamiliar motor synergies, or are deemed to be incompatible, internal adaptive models representing movement dynamics are inaccurate. This means that the actual response which is produced will deviate from the one which is planned. For a given target-response discrepancy, corrective response trajectories of longer duration are planned, consistent with the principle of speed-accuracy trade-off. Compared to familiar or compatible tasks, this results in a longer response time delay and reduced accuracy. From the standpoint of the intermittency approach, the findings of this study help make possible a more integral and predictive account of purposive action. (c) 2005 Elsevier B.V. All rights reserved.

Risks for ross river virus disease in tropical Australia

Background There are no analytical studies of individual risks for Ross River virus (RRV) disease. Therefore, we set out to determine individual risk and protective factors for RRV disease in a high incidence area and to assess the utility of the case-control design applied for this purpose to an arbovirus disease. Methods We used a prospective matched case-control study of new community cases of RRV disease in the local government areas of Cairns, Mareeba, Douglas, and Atherton, in tropical Queensland, from January I to May 31, 1998. Results Protective measures against mosquitoes reduced the risk for disease. Mosquito coils, repellents, and citronella candles each decreased risk by at least 2-fold, with a dose-response for the number of protective measures used. Light-coloured clothing decreased risk 3-fold. Camping increased the risk 8-fold. Conclusions These risks were substantial and statistically significant, and provide a basis for educational programs on individual protection against RRV disease in Australia. Our study demonstrates the utility of the case-control method for investigating arbovirus risks. Such a risk analysis has not been done before for RRV infection, and is infrequently reported for other arbovirus infections.

Effect of knee joint angle on motor unit synchronization

Activity of the vasti has been argued to vary through knee range of movement due to changes in passive support of the patellofemoral joint and the relative contribution of these muscles to knee extension. Efficient function of the knee is dependent on optimal control of the patellofemoral joint, largely through coordinated activity of the medial and lateral quadriceps. Motor unit synchronization may provide a mechanism to coordinate the activity of vastus medialis (VMO) and vastus lateralis (VL), and may be more critical in positions of reduced passive support for the patellofemoral joint (i.e., full extension). Therefore, the aim of this study was to determine whether the degree of motor unit synchronization between the vasti muscles is dependent on joint angle. Electromyographic (EMG) recordings of single motor unit action potentials (MUAPs) were made from VMO and multiunit recordings from VL during isometric contractions of the quadriceps at 0 degrees, 30 degrees, and 60 degrees of knee flexion. The degree of synchronization between motor unit firing was evaluated by identification of peaks in the rectified EMG averages of VL, triggered from MUA-Ps in VMO. The proportion of cases in which there was a significant peak in the triggered averages was calculated. There was no significant difference in the degree of synchronization between the vasti at different knee angles (p = 0.57). These data suggest that this basic coordinative mechanism between the vasti muscles is controlled consistently throughout knee range of motion, and is not augmented at specific angles where the requirement for dynamic control of stability is increased. (D 2006 Orthopaedic Research Society. Published by Wiley Periodicals, Inc.