A formal representation of assumptions in process modelling

In this work, we present a systematic approach to the representation of modelling assumptions. Modelling assumptions form the fundamental basis for the mathematical description of a process system. These assumptions can be translated into either additional mathematical relationships or constraints between model variables, equations, balance volumes or parameters. In order to analyse the effect of modelling assumptions in a formal, rigorous way, a syntax of modelling assumptions has been defined. The smallest indivisible syntactical element, the so called assumption atom has been identified as a triplet. With this syntax a modelling assumption can be described as an elementary assumption, i.e. an assumption consisting of only an assumption atom or a composite assumption consisting of a conjunction of elementary assumptions. The above syntax of modelling assumptions enables us to represent modelling assumptions as transformations acting on the set of model equations. The notion of syntactical correctness and semantical consistency of sets of modelling assumptions is defined and necessary conditions for checking them are given. These transformations can be used in several ways and their implications can be analysed by formal methods. The modelling assumptions define model hierarchies. That is, a series of model families each belonging to a particular equivalence class. These model equivalence classes can be related to primal assumptions regarding the definition of mass, energy and momentum balance volumes and to secondary and tiertinary assumptions regarding the presence or absence and the form of mechanisms within the system. Within equivalence classes, there are many model members, these being related to algebraic model transformations for the particular model. We show how these model hierarchies are driven by the underlying assumption structure and indicate some implications on system dynamics and complexity issues. (C) 2001 Elsevier Science Ltd. All rights reserved.

Design, modelling and analysis of a six component force balance for hypervelocity wind tunnel testing

A combination of modelling and analysis techniques was used to design a six component force balance. The balance was designed specifically for the measurement of impulsive aerodynamic forces and moments characteristic of hypervelocity shock tunnel testing using the stress wave force measurement technique. Aerodynamic modelling was used to estimate the magnitude and distribution of forces and finite element modelling to determine the mechanical response of proposed balance designs. Simulation of balance performance was based on aerodynamic loads and mechanical responses using convolution techniques. Deconvolution was then used to assess balance performance and to guide further design modifications leading to the final balance design. (C) 2001 Elsevier Science Ltd. All rights reserved.

Modelling cement grinding circuits

Modelling and simulation studies were carried out at 26 cement clinker grinding circuits including tube mills, air separators and high pressure grinding rolls in 8 plants. The results reported earlier have shown that tube mills can be modelled as several mills in series, and the internal partition in tube mills can be modelled as a screen which must retain coarse particles in the first compartment but not impede the flow of drying air. In this work the modelling has been extended to show that the Tromp curve which describes separator (classifier) performance can be modelled in terms of d(50)(corr), by-pass, the fish hook, and the sharpness of the curve. Also the high pressure grinding rolls model developed at the Julius Kruttschnitt Mineral Research Centre gives satisfactory predictions using a breakage function derived from impact and compressed bed tests. Simulation studies of a full plant incorporating a tube mill, HPGR and separators showed that the models could successfully predict the performance of the another mill working under different conditions. The simulation capability can therefore be used for process optimization and design. (C) 2001 Elsevier Science Ltd. All rights reserved.

Historical PCDD Inputs and Their Source Implications from Dated Sediment Cores in Queensland (Australia)

Recent investigations have demonstrated the presence of an unidentified source of polychlorinated dibenzo-p-dioxins (PCDDs) in the coastal zone of Queensland (Australia). The present study provides new information on the possible PCDD sources and their temporal input to this environment. Two estuarine sediment cores were collected in northern Queensland for which radiochemical chronologies were established. Core sections from different depositional ages, up to three centuries, have been analyzed for 2,3,7,8-substituted PCDDs and polychlorinated dibenzofurans (PCDFs). Variations of PCDD concentrations in the sediment cores over several centuries of depositional history were relatively small, and elevated PCDD levels were still present in sediment slices from the early 17th century. PCDD/F isomer patterns and congener profiles in sediments deposited during the last 350 years were almost identical and correlated well to the characteristic profiles observed in surface sediments and soils from the entire Queensland coastline. Profiles were dominated by higher chlorinated PCDDs, in particular octachlorodibenzodioxin (OCDD), whereas PCDF concentrations were below or near the limit of detection. These results indicate the presence of a PCDD source prior to industrialization and production of commercial organochlorine products. Further, the present study demonstrates that PCDD input patterns have been similar along an extensive but localized area over at least several centuries, contributing relatively high concentrations of PCDDs to the coastal system of Queensland.

PCDDs in the water/sediment-seagrass-dugong (Dugong dugon) food chain on the Great Barrier Reef (Australia)

Polychlorinated dibenzo-p-dioxin (PCDD) and dibenzofuran (PCDF) concentrations were measured in sediment and seagrass from five locations in or adjacent to the Great Barrier Reef Marine Park. A full spectrum of Cl(5-8)DDs were present in all samples and, in particular, elevated levels of Cl8DD were found. PCDFs could not be quantified in any samples. The PCDD concentrations ranged over two orders of magnitude between sites, and there was a good correlation between sediment and seagrass levels. There were large quantities of sediment present on the seagrass (20-62% on a dry wt. basis), and it was concluded that this was a primary source of the PCDDs in the seagrass samples. The PCDD levels in the seagrass samples were compared with the levels in the tissue of three dugongs stranded in the same region. The relative accumulation of the 2,3,7,8-substituted PCDD congeners in the dugongs decreased by over two orders of magnitude with increasing degree of chlorination. This was attributed to the reduced absorption of the higher chlorinated congeners in the digestive tract, a behaviour that has been observed in other mammals such as domestic cows. (C) 2001 Elsevier Science Ltd. All rights reserved.

Probabilistic modelling of CO2 corrosion laboratory data using neural networks

The present paper addresses two major concerns that were identified when developing neural network based prediction models and which can limit their wider applicability in the industry. The first problem is that it appears neural network models are not readily available to a corrosion engineer. Therefore the first part of this paper describes a neural network model of CO2 corrosion which was created using a standard commercial software package and simple modelling strategies. It was found that such a model was able to capture practically all of the trends noticed in the experimental data with acceptable accuracy. This exercise has proven that a corrosion engineer could readily develop a neural network model such as the one described below for any problem at hand, given that sufficient experimental data exist. This applies even in the cases when the understanding of the underlying processes is poor. The second problem arises from cases when all the required inputs for a model are not known or can be estimated with a limited degree of accuracy. It seems advantageous to have models that can take as input a range rather than a single value. One such model, based on the so-called Monte Carlo approach, is presented. A number of comparisons are shown which have illustrated how a corrosion engineer might use this approach to rapidly test the sensitivity of a model to the uncertainities associated with the input parameters. (C) 2001 Elsevier Science Ltd. All rights reserved.

The modelling of froth zone recovery in batch and continuously operated laboratory flotation cells

This paper proposes an integrated methodology for modelling froth zone performance in batch and continuously operated laboratory flotation cells. The methodology is based on a semi-empirical approach which relates the overall flotation rate constant to the froth depth (FD) in the flotation cell; from this relationship, a froth zone recovery (R,) can be extracted. Froth zone recovery, in turn, may be related to the froth retention time (FRT), defined as the ratio of froth volume to the volumetric flow rate of concentrate from the cell. An expansion of this relationship to account for particles recovered both by true flotation and entrainment provides a simple model that may be used to predict the froth performance in continuous tests from the results of laboratory batch experiments. Crown Copyright (C) 2002 Published by Elsevier Science B.V. All rights reserved.

Combined hydraulic and biological modelling and full-scale validation of SBR process

The biological reactions during the settling and decant periods of Sequencing Batch Reactors (SBRs) are generally ignored as they are not easily measured or described by modelling approaches. However, important processes are taking place, and in particular when the influent is fed into the bottom of the reactor at the same time (one of the main features of the UniFed process), the inclusion of these stages is crucial for accurate process predictions. Due to the vertical stratification of both liquid and solid components, a one-dimensional hydraulic model is combined with a modified ASM2d biological model to allow the prediction of settling velocity, sludge concentration, soluble components and biological processes during the non-mixed periods of the SBR. The model is calibrated on a full-scale UniFed SBR system with tracer breakthrough tests, depth profiles of particulate and soluble compounds and measurements of the key components during the mixed aerobic period. This model is then validated against results from an independent experimental period with considerably different operating parameters. In both cases, the model is able to accurately predict the stratification and most of the biological reactions occurring in the sludge blanket and the supernatant during the non-mixed periods. Together with a correct description of the mixed aerobic period, a good prediction of the overall SBR performance can be achieved.

Modelling of Lamb waves in composite laminated plates excited by interdigital transducers

The technique of permanently attaching interdigital transducers (IDT) to either flat or curved structural surfaces to excite single Lamb wave mode has demonstrated great potential for quantitative non-destructive evaluation and smart materials design, In this paper, the acoustic wave field in a composite laminated plate excited by an IDT is investigated. On the basis of discrete layer theory and a multiple integral transform method, an analytical-numerical approach is developed to evaluate the surface velocity response of the plate due to the IDTs excitation. In this approach, the frequency spectrum and wave number spectrum of the output of IDT are obtained directly. The corresponding time domain results are calculated by applying a standard inverse fast Fourier transformation technique. Numerical examples are presented to validate the developed method and show the ability of mode selection and isolation. A new effective way of transfer function estimation and interpretation is presented by considering the input wave number spectrum in addition to the commonly used input frequency spectrum. The new approach enables the simple physical evaluation of the influences of IDT geometrical features such as electrode finger widths and overall dimension and excitation signal properties on the input-output characteristics of IDT. Finally, considering the convenience of Mindlin plate wave theory in numerical computations as well as theoretical analysis, the validity is examined of using this approximate theory to design IDT for the excitation of the first and second anti-symmetric Lamb modes. (C) 2002 Elsevier Science Ltd. All rights reserved.

Shear stress and sediment transport calculations for sheet flow under waves

A simple method is provided for calculating transport rates of not too fine (d(50) greater than or equal to 0.20 mm) sand under sheet flow conditions. The method consists of a Meyer-Peter-type transport formula operating on a time-varying Shields parameter, which accounts for both acceleration-asymmetry and boundary layer streaming. While velocity moment formulae, e.g.., = Constant x calibrated against U-tube measurements, fail spectacularly under some real waves (Ribberink, J.S., Dohmen-Janssen, C.M., Hanes, D.M., McLean, S.R., Vincent, C., 2000. Near-bed sand transport mechanisms under waves. Proc. 27th Int. Conf. Coastal Engineering, Sydney, ASCE, New York, pp. 3263-3276, Fig. 12), the new method predicts the real wave observations equally well. The reason that the velocity moment formulae fail under these waves is partly the presence of boundary layer streaming and partly the saw-tooth asymmetry, i.e., the front of the waves being steeper than the back. Waves with saw-tooth asymmetry may generate a net landward sediment transport even if = 0, because of the more abrupt acceleration under the steep front. More abrupt accelerations are associated with thinner boundary layers and greater pressure gradients for a given velocity magnitude. The two real wave effects are incorporated in a model of the form Q(s)(t) = Q(s)[theta(t)] rather than Q(S)(t) = Q(S)[u(infinity)(t)], i.e., by expressing the transport rate in terms of an instantaneous Shields parameter rather than in terms of the free stream velocity, and accounting for both streaming and accelerations in the 0(t) calculations. The instantaneous friction velocities u(*)(t) and subsequently theta(t) are calculated as follows. Firstly, a linear filter incorporating the grain roughness friction factor f(2.5) and a phase angle phi(tau) is applied to u(infinity)(t). This delivers u(*)(t) which is used to calculate an instantaneous grain roughness Shields parameter theta(2.5)(t). Secondly, a constant bed shear stress is added which corresponds to the streaming related bed shear stress -rho ($) over bar((u) over tilde(w) over tilde)(infinity) . The method can be applied to any u(infinity)(t) time series, but further experimental validation is recommended before application to conditions that differ strongly from the ones considered below. The method is not recommended for rippled beds or for sheet flow with typical prototype wave periods and d(50) < 0.20 turn. In such scenarios, time lags related to vertical sediment movement become important, and these are not considered by the present model. (C) 2002 Elsevier Science B.V. All rights reserved.