A new adaptation of mapping method to study mixing of multiphase flows in mixers with complex geometries

The main objective of the present work is to modify the traditional mapping method for the simulation of distributive mixing of multiphase flows in geometries involving moving parts such as, internal mixers or twin-screw extruders without a limitation on their geometrical periodicity. The periodicity condition, limits the results of traditional mapping method to tracking mapping mesh between specific discrete time intervals or distances for that geometry is repeated, hence, result is only for fixed orientation of rotors. Imaginary domain method is introduced to track mapping mesh from one state to the next free of geometrical periodicity limitations. In this work the method is introduced and its applicability and accuracy is discussed in details. A two-dimensional (2D) simulation of mixing of two Newtonian fluids with different viscosities in an intermeshing internal mixer is used as a test case study. In this example the key issues of ability to predict mixing state in details for all orientations of rotors is presented. To reduce diffusion errors of mapping method in the boundaries of the rotors, mapping mesh refinement technique that relies upon one single reference mesh is also presented.     

Automating application software generation

Many approaches have been proposed to enhance software productivity and reliability. These approaches typically fall into three categories: the engineering approach, the formal approach and the knowledge-based approach. But the optimal gain in software productivity cannot be obtained if one relies on only one of these approaches. This paper describes the work in knowledge-based software engineering conducted by the authors for the past 10 years. The final goal of the research is to develop a paradigm for software engineering which integrates the three approaches mentioned above. A knowledge-based tool which can support the whole process of software development is provided in this paper.     

Non‐reflecting artificial boundaries for transient scalar wave propagation in a two‐dimensional infinite homogeneous layer

This paper presents an exact non‐reflecting boundary condition for dealing with transient scalar wave propagation problems in a two‐dimensional infinite homogeneous layer. In order to model the complicated geometry and material properties in the near field, two vertical artificial boundaries are considered in the infinite layer so as to truncate the infinite domain into a finite domain. This treatment requires the appropriate boundary conditions, which are often referred to as the artificial boundary conditions, to be applied on the truncated boundaries. Since the infinite extension direction is different for these two truncated vertical boundaries, namely one extends toward x →∞ and another extends toward x→‐ ∞, the non‐reflecting boundary condition needs to be derived on these two boundaries. Applying the variable separation method to the wave equation results in a reduction in spatial variables by one. The reduced wave equation, which is a time‐dependent partial differential equation with only one spatial variable, can be further changed into a linear first‐order ordinary differential equation by using both the operator splitting method and the modal radiation function concept simultaneously. As a result, the non‐reflecting artificial boundary condition can be obtained by solving the ordinary differential equation whose stability is ensured. Some numerical examples have demonstrated that the non‐reflecting boundary condition is of high accuracy in dealing with scalar wave propagation problems in infinite and semi‐infinite media. Copyright © 2003 John Wiley & Sons, Ltd.     

动态矩阵控制系统稳定域的根轨迹分析

本文使用根轨迹技术从理论上分析了纯滞后工业过程参数变化对系统稳定性的影响,得到了一个与设计参数有关的稳定域,从而为设计参数自动整定提供了参考。     

THE ANALYSIS OF THE SEPARATED LAYERS ALGORITHM BY DOMAIN DECOMPOSITION METHOD

From the principle of of the Domain Decomposition Method (DDM), we analyse the 2nd-order linear elliptic partial differential problems and link the Separated-Layers Algorithm (SLA) with DDM. The mathematical properties of SLA and numerical example are presented to obtain satisfactory computation results. For general linear differential ones, also are the structure of SLA and its characteristics discussed.     

Automated Support for Building and Extending Expert Models

Building a knowledge-based system is like developing a scientific theory. Although a knowledge base does not constitute a theory of some natural phenomenon, it does represent a theory of how a class of professionals approaches an application task. As when scientists develop a natural theory, builders of expert systems first must formulate a model of the behavior that they wish to understand and then must corroborate and extend that model with the aid of specific examples. Thus there are two interrelated phases of knowledge-base construction: (1) model building and (2) model extension. Computer-based tools can assist developers with both phases of the knowledge-acquisition process. Workers in the area of knowledge acquisition have developed computer-based tools that emphasize either the building of new models or the extension of existing models. The PROTÉGÉ knowledge-acquisition system addresses these two activities individually and facilitates the construction of expert systems when the same general model can be applied to a variety of application tasks.     

Stability analysis and antiwindup design of uncertain discrete-time switched linear systems subject to actuator saturation

This paper investigates the stability analysis and antiwindup design problem for a class of discrete-time switched linear systems with time-varying norm-bounded uncertainties and saturating actuators by using the switched Lyapunov function approach.Supposing that a set of linear dynamic output controllers have been designed to stabilize the switched system without considering its input saturation,we design antiwindup compensation gains in order to enlarge the domain of attraction of the closed-loop system in the presence of saturation.Then,in terms of a sector condition,the antiwindup compensation gains which aim to maximize the estimation of domain of attraction of the closed-loop system are presented by solving a convex optimization problem with linear matrix inequality(LMI)constraints.A numerical example is given to demonstrate the effectiveness of the proposed design method.     

Dynamics of fluid circulation in coupled free and heterogeneous porous domains

Fluid flow in coupled free and porous domain, particularly when the porous medium is heterogeneous, is encountered in many hydro-environmental conditions, e.g., leakage from underground pipe, combined groundwater lake-subsurface interactions. One of the most difficult problems in the study of coupled flow behaviour has been the development of a universally applicable modelling scheme for combining the flow regimes. This is because the free/porous interfacial properties (e.g., shear-stress; velocity slip) that govern the coupled flow behaviour are difficult to determine experimentally under hydro-environmental conditions. On the other hand, the implications of various forms of heterogeneity in the porous media properties can be very different on the fluid-flow behaviour. Difficulties may also arise in direct coupling of the model equations that govern the fluid flow in the individual regions (e.g., Navier-Stokes for free-flow region and the Darcy's equation for the porous flow region). Consequently, models of coupled free and porous flow for hydro-environmental conditions are not very well developed at the moment. While there are some indications that fluids in coupled free and porous domains may circulate (i.e., development of flow cells), there is a lack of appropriate 3D analysis on how heterogeneities in porous media may affect such flow patterns. In this paper, we aim to analyse how porous media heterogeneity affects the dynamics of flow circulation in the porous side of a coupled free and porous domain. For this purpose, we analyse flow patterns in several model domains made up of two porous layers with differing permeabilities. The governing model equations are discretised and solved using the standard finite volume method on a staggered cell-centred mesh. The temporal discretisation is done using the explicit method. An in-house graphical user interface (GUI) has been created specifically to aid in the visualisation of otherwise complex flow patterns. The GUI contains many post-processing options and provides a comprehensive tool for the analysis of hydrodynamics and contaminant motion (not discussed in this paper) in coupled free and porous flow domains. This GUI is described in this paper briefly. The effects of altering the aspect ratio (i.e., multi-scale) of the domain on the coupled flow pattern have also been discussed.     

Specificity of the pheromone system ofAdoxophyes orana andClepsis spectrana

The pheromone system of the summer fruit tortrix mothAdoxophyes orana (Fischer von Röslerstamm) consists of a mixture ofcis-9- andcis-11-tetradecen-1-ol acetate (cis-9- andcis-11-TDA, respectively) in a ratio of 9:1. Substitution of one or both of these compounds by related unsaturated acetates reduced the attractancy. Onlycis-11-TDA could be replaced bycis-11-tridecen-1-ol acetate (cis-11-TriDA), although the 9:1 mixture ofcis-9-TDA andcis-11-TriDA was less attractive in the field than the pheromone system itself. The major component of the pheromone system of the leaf rollerClepsis spectrana (Treitschke)(cis-9-TDA:cis-11-TDA=1:9) could also be replaced bycis-11-TriDA, again with some reduction in attractancy.     

Total curvature variation fairing for medial axis regularization

We present a new fairing method for planar curves, which is particularly well suited for the regularization of the medial axis of a planar domain. It is based on the concept of total variation regularization. The original boundary (given as a closed B-spline curve or several such curves for multiply connected domains) is approximated by another curve that possesses a smaller number of curvature extrema. Consequently, the modified curve leads to a smaller number of branches of the medial axis. In order to compute the medial axis, we use the state-of-the-art algorithm from [1] which is based on arc spline approximation and a domain decomposition approach. We improve this algorithm by using a different decomposition strategy that allows to reduce the number of base cases from 13 to only 5. Moreover, the algorithm reduces the number of conic arcs in the output by approx. 50%.