Analysis of dynamic process models for structural insight and model reduction .1. Structural identification measures

An important consideration in the development of mathematical models for dynamic simulation, is the identification of the appropriate mathematical structure. By building models with an efficient structure which is devoid of redundancy, it is possible to create simple, accurate and functional models. This leads not only to efficient simulation, but to a deeper understanding of the important dynamic relationships within the process. In this paper, a method is proposed for systematic model development for startup and shutdown simulation which is based on the identification of the essential process structure. The key tool in this analysis is the method of nonlinear perturbations for structural identification and model reduction. Starting from a detailed mathematical process description both singular and regular structural perturbations are detected. These techniques are then used to give insight into the system structure and where appropriate to eliminate superfluous model equations or reduce them to other forms. This process retains the ability to interpret the reduced order model in terms of the physico-chemical phenomena. Using this model reduction technique it is possible to attribute observable dynamics to particular unit operations within the process. This relationship then highlights the unit operations which must be accurately modelled in order to develop a robust plant model. The technique generates detailed insight into the dynamic structure of the models providing a basis for system re-design and dynamic analysis. The technique is illustrated on the modelling for an evaporator startup. Copyright (C) 1996 Elsevier Science Ltd

Analysis of dynamic models for structural insight and model reduction .2. A multi-stage compressor shutdown case-study

In this second paper, the three structural measures which have been developed are used in the modelling of a three stage centrifugal synthesis gas compressor. The goal of this case study is to determine the essential mathematical structure which must be incorporated into the compressor model to accurately model the shutdown of this system. A simple, accurate and functional model of the system is created via three structural measures. It was found that the model can be correctly reduced into its basic modes and that the order of the differential system can be reduced from 51(st) to 20(th). Of the 31 differential equational 21 reduce to algebraic relations, 8 become constants and 2 can be deleted thereby increasing the algebraic set from 70 to 91 equations. An interpretation is also obtained as to which physical phenomena are dominating the dynamics of the compressor add whether the compressor will enter surge during the shutdown. Comparisons of the reduced model performance against the full model are given, showing the accuracy and applicability of the approach. Copyright (C) 1996 Elsevier Science Ltd

Impact of Stent-Graft Oversizing on the Thoracic Aorta: Experimental Study in a Porcine Model

Purpose: To analyze in an experimental animal model the effect of 4 different levels of stents-graft oversizing on non-atherosclerotic aortas such as those found in young individuals who undergo stent-graft repair for traumatic aortic injuries. Methods: The diameter of the porcine thoracic aorta is similar to the aorta of young adults (18-20 mm), so 25 pigs were randomized into 5 groups: 1 control (without stent-graft) and 4 oversizing groups (A: 10%-19%, B: 20%-29%, C: 30%-39%, and D: >40%). Two types of biomechanical tests were performed on all aortas 4 weeks after endoprosthesis deployment. Results: The results of the detachment test, which analyzed the strength necessary to remove the stent-graft from the aorta, were similar in the 4 groups (A: 42 N, B: 41 N, C: 46 N, and D: 46 N). However, 2 aortas ruptured during the tests (groups C and D). The second test was performed in 3 aortic segments. Maximum shear strength, maximum stress, and maximum tension supported by the aortic wall had a negative and linear correlation with oversizing. There were significant differences in all 4 groups when compared with the control group. Strain, which reflects the elastic properties of the aortic wall, was very similar in all 4 groups, but a great difference was found when compared with the control group (p<0.0001). Conclusion: The study showed an important subacute change in the biomechanical properties of the aortic wall after implantation of an oversized endoprosthesis. This weakness of the aortic wall was confirmed by 2 ruptures during the detachment test. These results partially explain the interaction of stent-grafts with non-atherosclerotic thoracic aortas and may serve as a basis for further studies and the development of specific material to be used in vascular trauma and young patients. J Endovasc Ther. 2011; 18: 576-584

Derivation and external validation of a simple prediction model for the diagnosis of type 2 Diabetes Mellitus in the Brazilian urban population

A risk score model was developed based in a population of 1,224 individuals from the general population without known diabetes aging 35 years or more from an urban Brazilian population sample in order to select individuals who should be screened in subsequent testing and improve the efficacy of public health assurance. External validation was performed in a second, independent, population from a different city ascertained through a similar epidemiological protocol. The risk score was developed by multiple logistic regression and model performance and cutoff values were derived from a receiver operating characteristic curve. Model`s capacity of predicting fasting blood glucose levels was tested analyzing data from a 5-year follow-up protocol conducted in the general population. Items independently and significantly associated with diabetes were age, BMI and known hypertension. Sensitivity, specificity and proportion of further testing necessary for the best cutoff value were 75.9, 66.9 and 37.2%, respectively. External validation confirmed the model`s adequacy (AUC equal to 0.72). Finally, model score was also capable of predicting fasting blood glucose progression in non-diabetic individuals in a 5-year follow-up period. In conclusion, this simple diabetes risk score was able to identify individuals with an increased likelihood of having diabetes and it can be used to stratify subpopulations in which performing of subsequent tests is necessary and probably cost-effective.

Comparison of HRT-3 glaucoma probability score and subjective stereophotograph assessment for prediction of progression in glaucoma

PURPOSE. To assess whether baseline Glaucoma Probability Score (GPS; HRT-3; Heidelberg Engineering, Dossenheim, Germany) results are predictive of progression in patients with suspected glaucoma. The GPS is a new feature of the confocal scanning laser ophthalmoscope that generates an operator-independent, three-dimensional model of the optic nerve head and gives a score for the probability that this model is consistent with glaucomatous damage. METHODS. The study included 223 patients with suspected glaucoma during an average follow-up of 63.3 months. Included subjects had a suspect optic disc appearance and/or elevated intraocular pressure, but normal visual fields. Conversion was defined as development of either repeatable abnormal visual fields or glaucomatous deterioration in the appearance of the optic disc during the study period. The association between baseline GPS and conversion was investigated by Cox regression models. RESULTS. Fifty-four (24.2%) eyes converted. In multivariate models, both higher values of GPS global and subjective stereophotograph assessment ( larger cup-disc ratio and glaucomatous grading) were predictive of conversion: adjusted hazard ratios (95% CI): 1.31 (1.15 - 1.50) per 0.1 higher global GPS, 1.34 (1.12 - 1.62) per 0.1 higher CDR, and 2.34 (1.22 - 4.47) for abnormal grading, respectively. No significant differences ( P > 0.05 for all comparisons) were found between the c-index values ( equivalent to area under ROC curve) for the multivariate models (0.732, 0.705, and 0.699, respectively). CONCLUSIONS. GPS values were predictive of conversion in our population of patients with suspected glaucoma. Further, they performed as well as subjective assessment of the optic disc. These results suggest that GPS could potentially replace stereophotograph as a tool for estimating the likelihood of conversion to glaucoma.

Composite pre-heating: Effects on marginal adaptation, degree of conversion and mechanical properties

Objectives. This study evaluated the effect of composite pre-polymerization temperature and energy density on the marginal adaptation (MA), degree of conversion (DC), flexural strength (FS), and polymer cross-linking (PCL) of a resin composite (Filtek Z350, 3M/ESPE). Methods. For MA, class V cavities (4mmx2mmx2mm) were prepared in 40 bovine incisors. The adhesive system Adper Single Bond 2 (3M/ESPE) was applied. Before being placed in the cavities, the resin composite was either kept at room-temperature (25 degrees C) or previously pre-heated to 68 degrees C in the Calset (TM) device (AdDent Inc., Danbury, CT, USA). The composite was then light polymerized for 20 or 40s at 600mW/cm(2) (12 or 24 J/cm(2), respectively). The percentage of gaps was analyzed by scanning electron microscopy, after sectioning the restorations and preparing epoxy resin replicas. DC (n = 3) was obtained by FT-Raman spectroscopy on irradiated and non-irradiated composite surfaces. FS (n = 10) was measured by the three-point-bending test. KHN (n = 6) was measured after 24h dry storage and again after immersion in 100% ethanol solution for 24 h, to calculate PCL density. Data were analyzed by appropriate statistical analyses. Results. The pre-heated composite showed better MA than the room-temperature groups. A higher number of gaps were observed in the room-temperature groups, irrespective of the energy density, mainly in the axial wall (p < 0.05). Composite pre-heating and energy density did not affect the DC, FS and PCL (p > 0.05). Significance. Pre-heating the composite prior to light polymerization similar in a clinical situation did not alter the mechanical properties and monomer conversion of the composite, but provided enhanced composite adaptation to cavity walls. (C) 2010 Academy of Dental Materials. Published by Elsevier Ltd. All rights reserved.

Model based procedure for scale-up of wet overflow ball mills Part 1: Outline of the methodology

Bond's method for ball mill scale-up only gives the mill power draw for a given duty. This method is incompatible with computer modelling and simulation techniques. It might not be applicable for the design of fine grinding ball mills and ball mills preceded by autogenous and semi-autogenous grinding mills. Model-based ball mill scale-up methods have not been validated using a wide range of full-scale circuit data. Their accuracy is therefore questionable. Some of these methods also need expensive pilot testing. A new ball mill scale-up procedure is developed which does not have these limitations. This procedure uses data from two laboratory tests to determine the parameters of a ball mill model. A set of scale-up criteria then scales-up these parameters. The procedure uses the scaled-up parameters to simulate the steady state performance of full-scale mill circuits. At the end of the simulation, the scale-up procedure gives the size distribution, the volumetric flowrate and the mass flowrate of all the streams in the circuit, and the mill power draw.

Model-based procedure for scale-up of wet, overflow ball mills - Part II: Worked example

A new ball mill scale-up procedure is developed which uses laboratory data to predict the performance of MI-scale ball mill circuits. This procedure contains two laboratory tests. These laboratory tests give the data for the determination of the parameters of a ball mill model. A set of scale-up criteria then scales-up these parameters. The procedure uses the scaled-up parameters to simulate the steady state performance of the full-scale mill circuit. At the end of the simulation, the scale-up procedure gives the size distribution, the volumetric flowrate and the mass flowrate of all the streams in the circuit, and the mill power draw. A worked example shows how the new ball mill scale-up procedure is executed. This worked example uses laboratory data to predict the performance of a full-scale re-grind mill circuit. This circuit consists of a ball mill in closed circuit with hydrocyclones. The MI-scale ball mill has a diameter (inside liners) of 1.85m. The scale-up procedure shows that the full-scale circuit produces a product (hydrocyclone overflow) that has an 80% passing size of 80 mum. The circuit has a recirculating load of 173%. The calculated power draw of the full-scale mill is 92kW (C) 2001 Elsevier Science Ltd. All rights reserved.

A generalised dynamic model for char particle gasification with structure evolution and peripheral fragmentation

A generalised model for the prediction of single char particle gasification dynamics, accounting for multi-component mass transfer with chemical reaction, heat transfer, as well as structure evolution and peripheral fragmentation is developed in this paper. Maxwell-Stefan analysis is uniquely applied to both micro and macropores within the framework of the dusty-gas model to account for the bidisperse nature of the char, which differs significantly from the conventional models that are based on a single pore type. The peripheral fragmentation and random-pore correlation incorporated into the model enable prediction of structure/reactivity relationships. The occurrence of chemical reaction within the boundary layer reported by Biggs and Agarwal (Chem. Eng. Sci. 52 (1997) 941) has been confirmed through an analysis of CO/CO2 product ratio obtained from model simulations. However, it is also quantitatively observed that the significance of boundary layer reaction reduces notably with the reduction of oxygen concentration in the flue gas, operational pressure and film thickness. Computations have also shown that in the presence of diffusional gradients peripheral fragmentation occurs in the early stages on the surface, after which conversion quickens significantly due to small particle size. Results of the early commencement of peripheral fragmentation at relatively low overall conversion obtained from a large number of simulations agree well with experimental observations reported by Feng and Bhatia (Energy & Fuels 14 (2000) 297). Comprehensive analysis of simulation results is carried out based on well accepted physical principles to rationalise model prediction. (C) 2001 Elsevier Science Ltd. AH rights reserved.

The future role of government in knowledge-based economies

In this article, we draw together aspects of contemporary theories of knowledge (particularly organisational knowledge) and complexity theory to demonstrate how appropriate conceptual rigor enables both the role of government and the directions of policy development in knowledge-based economies to be identified. Specifically we ask, what is the role of government in helping shape the knowledge society of the future? We argue that knowledge policy regimes must go beyond the modes of policy analysis currently used in innovation, information and technology policy because they are based in an industrial rather than post-industrial analytical framework. We also argue that if we are to develop knowledge-based economies, more encompassing images of the future than currently obtain in policy discourse are required. We therefore seek to stimulate and provoke an array of lines of thought about government and policy for such economies. Our objective is to focus on ideas more than argument and persuasion.

Assessment of the ability to model proteins with leucine-rich repeats in light of the latest structural information

The three-dimensional structures of leucine-rich repeat (LRR) -containing proteins from five different families were previously predicted based on the crystal structure of the ribonuclease inhibitor. using an approach that combined homology-based modeling, structure-based sequence alignment of LRRs, and several rational assumptions. The structural models have been produced based on very limited sequence similarity, which, in general. cannot yield trustworthy predictions. Recently, the protein structures from three of these five families have been determined. In this report we estimate the quality of the modeling approach by comparing the models with the experimentally determined structures. The comparison suggests that the general architecture, curvature, interior/exterior orientations of side chains. and backbone conformation of the LRR structures can be predicted correctly. On the other hand. the analysis revealed that, in some cases. it is difficult to predict correctly the twist of the overall super-helical structure. Taking into consideration the conclusions from these comparisons, we identified a new family of bacterial LRR proteins and present its structural model. The reliability of the LRR protein modeling suggests that it would be informative to apply similar modeling approaches to other classes of solenoid proteins.

The copolymerization of N-vinyl-2-pyrrolidone with 2-hydroxyethyl methacrylate

The bulk free radical copolymerization of 2-hydroxyethyl methacrylate (HEMA) with N-vinyl-2-pyrrolidone (VP) was carried out to low conversions at 50 degreesC, using benzoyl peroxide (BPO) as initiator. The compositions of the copolymers; were determined using C-13 NMR spectroscopy. The conversion of monomers to polymers was studied using FT-NIR spectroscopy in order to predict the extent of conversion of monomer to polymer. From model fits to the composition data, a statistical F-test revealed that die penultimate model describes die copolymerization better than die terminal model. Reactivity ratios were calculated by using a non-linear least squares analysis (NLLS) and r(H) = 8.18 and r(V) = 0.097 were found to be the best fit values of the reactivity ratios for the terminal model and r(HH) = 12.0, r(VH) = 2.20, r(VV) = 0.12 and r(HV) = 0.03 for the penultimate model. Predictions were made for changes in compositions as a function of conversion based upon the terminal and penultimate models.

Entanglement in the steady state of a collective-angular-momentum (Dicke) model

We model the behavior of an ion trap with all ions driven simultaneously and coupled collectively to a heat bath. The equations for this system are similar to the irreversible dynamics of a collective angular momentum system known as the Dicke model. We show how the steady state of the ion trap as a dissipative many-body system driven far from equilibrium can exhibit quantum entanglement. We calculate the entanglement of this steady state for two ions in the trap and in the case of more than two ions we calculate the entanglement between two ions by tracing over all the other ions. The entanglement in the steady state is a maximum for the parameter values corresponding roughly to a bifurcation of a fixed point in the corresponding semiclassical dynamics. We conjecture that this is a general mechanism for entanglement creation in driven dissipative quantum systems.

A cassette ligation strategy with thioether replacement of three Gly-Gly peptide bonds: Total chemical synthesis of the 101 residue protein early pregnancy factor [psi(CH2S)(28-29,56-57,76-77)]

The 101 residue protein early pregnancy factor (EPF), also known as human chaperonin 10, was synthesized from four functionalized, but unprotected, peptide segments by a sequential thioether ligation strategy. The approach exploits the differential reactivity of a peptide-NHCH2CH2SH thiolate with XCH2CO-peptides, where X = Cl or I/Br. Initial model studies with short functionalized (but unprotected) peptides showed a significantly faster reaction of a peptide-NHCH2CH2SH thiolate with a BrCH2CO-peptide than with a CICH2CO-peptide, where thiolate displacement of the halide leads to chemoselective formation of a thioether surrogate for the Gly-Gly peptide bond. This rate difference was used as the basis of a novel sequential ligation approach to the synthesis of large polypeptide chains. Thus, ligation of a model bifunctional N-alpha-chloroacetyl, C-terminal thiolated peptide with a second N-alpha-bromoacetyl peptide demonstrated chemoselective bromide displacement by the thiol group. Further investigations showed that the relatively unreactive N-alpha-chloroacetyl peptides could be activated by halide exchange using saturated KI solutions to yield the highly reactive No-iodoacetyl peptides. These findings were used to formulate a sequential thioether ligation strategy for the synthesis of EPF, a 101 amino acid protein containing three Gly-Gly sites approximately equidistantly spaced within the peptide chain. Four peptide segments or cassettes comprising the EPF protein sequence (BrAc-[EPF 78-101] 12, ClAc-[EPF 58-75]-[NHCH2CH2SH] 13, ClAc-[EPF 30-55]-[NHCH2CH2SH] 14, and Ac-[EPF 1-27]-[NHCH2CH2SH] 15) of EPF were synthesized in high yield and purity using Boc SPPS chemistry. In the stepwise sequential ligation strategy, reaction of peptides 12 and 13 was followed by conversion of the N-terminal chloroacetyl functional group to an iodoacetyl, thus activating the product peptide for further ligation with peptide 14. The process of ligation followed by iodoacetyl activation was repeated to yield an analogue of EPF (EPF psi(CH2S)(28-29,56-57,76-77)) 19 in 19% overall yield.

Application of the Williams-Landel-Ferry model to the viscosity-temperature relationship of Australian honeys

The rheological behaviour of nine unprocessed Australian honeys was investigated for the applicability of the Williams-Landel-Ferry (WLF) model. The viscosity of the honeys was obtained over a range of shear rates (0.01-40 s(-1)) from 2degrees to 40 degreesC, and all the honeys exhibited Newtonian behaviour with viscosity reducing as the temperature was increased. The honeys with high moisture were of lower viscosity, The glass transition temperatures of the honeys, as measured with a differential scanning calorimeter (DSC), ranged from -40degrees to -46 degreesC, and four models (WLF. Arrhenius, Vogel-Tammann-Fulcher (VTF), and power-law) were investigated to describe the temperature dependence of the viscosity. The WLF was the most suitable and the correlation coefficient averaged 0.999 +/- 0.0013 as against 0.996 +/- 0.0042 for the Arrhenius model while the mean relative deviation modulus was 0-12% for the WLF model and 10-40% for the Arrhenius one. With the universal values for the WLF constants, the temperature dependence of the viscosity was badly predicted. From non-linear regression analysis, the constants of the WLF models for the honeys were obtained (C-1 = 13.7-21.1: C-2 = 55.9-118.7) and are different from the universal values. These WLF constants will be valuable for adequate modeling of the rheology of the honeys, and they can be used to assess the temperature sensitivity of the honeys. (C) 2002 Elsevier Science Ltd. All rights reserved.