Model reduction for nonlinear aerodynamics and aeroelasticity using a Discrete Empirical Interpolation Method

A novel surrogate model is proposed in lieu of Computational Fluid Dynamics (CFD) solvers, for fast nonlinear aerodynamic and aeroelastic modeling. A nonlinear function is identified on selected interpolation points by
a discrete empirical interpolation method (DEIM). The flow field is then reconstructed using a least square approximation of the flow modes extracted
by proper orthogonal decomposition (POD). The aeroelastic reduce order
model (ROM) is completed by introducing a nonlinear mapping function
between displacements and the DEIM points. The proposed model is investigated to predict the aerodynamic forces due to forced motions using
a N ACA 0012 airfoil undergoing a prescribed pitching oscillation. To investigate aeroelastic problems at transonic conditions, a pitch/plunge airfoil
and a cropped delta wing aeroelastic models are built using linear structural models. The presence of shock-waves triggers the appearance of limit
cycle oscillations (LCO), which the model is able to predict. For all cases
tested, the new ROM shows the ability to replicate the nonlinear aerodynamic forces, structural displacements and reconstruct the complete flow
field with sufficient accuracy at a fraction of the cost of full order CFD
model.

DYNAMIC MODEL BASED STATE ESTIMATION IN A HEAVY DUTY DIESEL AFTERTREATMENT SYSTEM FOR ONBOARD DIAGNOSTICS AND CONTROLS

Estimating un-measurable states is an important component for onboard diagnostics (OBD) and control strategy development in diesel exhaust aftertreatment systems. This research focuses on the development of an Extended Kalman Filter (EKF) based state estimator for two of the main components in a diesel engine aftertreatment system: the Diesel Oxidation Catalyst (DOC) and the Selective Catalytic Reduction (SCR) catalyst. One of the key areas of interest is the performance of these estimators when the catalyzed particulate filter (CPF) is being actively regenerated. In this study, model reduction techniques were developed and used to develop reduced order models from the 1D models used to simulate the DOC and SCR. As a result of order reduction, the number of states in the estimator is reduced from 12 to 1 per element for the DOC and 12 to 2 per element for the SCR. The reduced order models were simulated on the experimental data and compared to the high fidelity model and the experimental data. The results show that the effect of eliminating the heat transfer and mass transfer coefficients are not significant on the performance of the reduced order models. This is shown by an insignificant change in the kinetic parameters between the reduced order and 1D model for simulating the experimental data. An EKF based estimator to estimate the internal states of the DOC and SCR was developed. The DOC and SCR estimators were simulated on the experimental data to show that the estimator provides improved estimation of states compared to a reduced order model. The results showed that using the temperature measurement at the DOC outlet improved the estimates of the CO , NO , NO2 and HC concentrations from the DOC. The SCR estimator was used to evaluate the effect of NH3 and NOX sensors on state estimation quality. Three sensor combinations of NOX sensor only, NH3 sensor only and both NOX and NH3 sensors were evaluated. The NOX only configuration had the worst performance, the NH3 sensor only configuration was in the middle and both the NOX and NH3 sensor combination provided the best performance.

Dynamic equivalent state estimation for multi-area power systems with synchronized phasor measurement units

This paper proposes a new approach for state estimation of angles and frequencies of equivalent areas in large power systems with synchronized phasor measurement units. Defining coherent generators and their correspondent areas, generators are aggregated and system reduction is performed in each area of inter-connected power systems. The structure of the reduced system is obtained based on the characteristics of the reduced linear model and measurement data to form the non-linear model of the reduced system. Then a Kalman estimator is designed for the reduced system to provide an equivalent dynamic system state estimation using the synchronized phasor measurement data. The method is simulated on two test systems to evaluate the feasibility of the proposed method.

Post-transcriptional regulation in the nucleus and cytoplasm: study of mean time to threshold (MTT) and narrow escape problem

Messenger RNAs (mRNAs) can be repressed and degraded by small non-coding RNA molecules. In this paper, we formulate a coarsegrained Markov-chain description of the post-transcriptional regulation of mRNAs by either small interfering RNAs (siRNAs) or microRNAs (miRNAs). We calculate the probability of an mRNA escaping from its domain before it is repressed by siRNAs/miRNAs via cal- culation of the mean time to threshold: when the number of bound siRNAs/miRNAs exceeds a certain threshold value, the mRNA is irreversibly repressed. In some cases,the analysis can be reduced to counting certain paths in a reduced Markov model. We obtain explicit expressions when the small RNA bind irreversibly to the mRNA and we also discuss the reversible binding case. We apply our models to the study of RNA interference in the nucleus, examining the probability of mRNAs escaping via small nuclear pores before being degraded by siRNAs. Using the same modelling framework, we further investigate the effect of small, decoy RNAs (decoys) on the process of post-transcriptional regulation, by studying regulation of the tumor suppressor gene, PTEN : decoys are able to block binding sites on PTEN mRNAs, thereby educing the number of sites available to siRNAs/miRNAs and helping to protect it from repression. We calculate the probability of a cytoplasmic PTEN mRNA translocating to the endoplasmic reticulum before being repressed by miRNAs. We support our results with stochastic simulations

An analysis of assessment outcomes from eight years' operation of the Australian border weed risk assessment system

The majority of Australian weeds are exotic plant species that were intentionally introduced for a variety of horticultural and agricultural purposes. A border weed risk assessment system (WRA) was implemented in 1997 in order to reduce the high economic costs and massive environmental damage associated with introducing serious weeds. We review the behaviour of this system with regard to eight years of data collected from the assessment of species proposed for importation or held within genetic resource centres in Australia. From a taxonomic perspective, species from the Chenopodiaceae and Poaceae were most likely to be rejected and those from the Arecaceae and Flacourtiaceae were most likely to be accepted. Dendrogram analysis and classification and regression tree (TREE) models were also used to analyse the data. The latter revealed that a small subset of the 35 variables assessed was highly associated with the outcome of the original assessment. The TREE model examining all of the data contained just five variables: unintentional human dispersal, congeneric weed, weed elsewhere, tolerates or benefits from mutilation, cultivation or fire, and reproduction by vegetative propagation. It gave the same outcome as the full WRA model for 71% of species. Weed elsewhere was not the first splitting variable in this model, indicating that the WRA has a capacity for capturing species that have no history of weediness. A reduced TREE model (in which human-mediated variables had been removed) contained four variables: broad climate suitability, reproduction in less or than equal to 1 year, self-fertilisation, and tolerates and benefits from mutilation, cultivation or fire. It yielded the same outcome as the full WRA model for 65% of species. Data inconsistencies and the relative importance of questions are discussed, with some recommendations made for improving the use of the system.

Studies on characteristics of lightning generated currents in an interconnected lightning protection system

An in-depth knowledge about the characteristics of lightning generated currents will facilitate evaluation of the interception efficacy of lightning protection systems. In addition, it would aid in extraction of valuable statistics (from measured current data) on local lightning parameters. Incidentally, present day knowledge on characteristics of lightning induced current in typical lightning protection systems is rather limited. This is particularly true with closely interconnected protection systems, like the one employed in Indian Satellite Launch Pad-II. This system is taken as a specific example in the present study. Various aspects suggest that theoretical modelling would be the best possible approach for the intended work. From the survey of pertinent literature, it is concluded that electromagnetic modelling of lightning return-stroke with current source at the channel base is best suited for this study. Numerical electromagnetic code was used for the required electromagnetic field solution and Fourier transform techniques were employed for computing time-domain results. A validation for the numerical modelling is provided by laborator experiments on a reduced scale model of the system. Apart from ascertaining the influence of various parameters, salient characteristics of tower base currents for different kinds of events are deduced. This knowledge can be used in identifying the type of event, as well as its approximate location. A method for estimation of injected stroke current has also been proposed.

Investigations on voltages and currents in the lightning protection system of the Indian satellite launch pad-I during a stroke interception

A reliable protection against direct lightning hit is very essential for satellite launch pads. In view of this, suitable protection systems are generally employed. The evaluation of efficacy of the lightning protection schemes among others requires an accurate knowledge of the consequential potential rise at the struck point and the current injected into soil at the earth termination. The present work has made a detailed effort to deduce these quantities for the lightning protection scheme of the Indian satellite launch pad-I. A reduced scale model of the system with a frequency domain approach is employed for the experimental study. For further validation of the experimental approach, numerical simulations using numerical electromagnetic code-2 are also carried out on schemes involving single tower. The study results on the protection system show that the present design is quite safe with regard to top potential rise. It is shown that by connecting ground wires to the tower, its base current and, hence, the soil potential rise can be reduced. An evaluation of an alternate design philosophy involving insulated mast scheme is also made. The potential rise in that design is quantified and the possibility of a flashover to supporting tower is briefly looked into. The supporting tower is shown to have significant induced currents.

An invariant subspace-based approach to the random eigenvalue problem of systems with clustered spectrum

The repeated or closely spaced eigenvalues and corresponding eigenvectors of a matrix are usually very sensitive to a perturbation of the matrix, which makes capturing the behavior of these eigenpairs very difficult. Similar difficulty is encountered in solving the random eigenvalue problem when a matrix with random elements has a set of clustered eigenvalues in its mean. In addition, the methods to solve the random eigenvalue problem often differ in characterizing the problem, which leads to different interpretations of the solution. Thus, the solutions obtained from different methods become mathematically incomparable. These two issues, the difficulty of solving and the non-unique characterization, are addressed here. A different approach is used where instead of tracking a few individual eigenpairs, the corresponding invariant subspace is tracked. The spectral stochastic finite element method is used for analysis, where the polynomial chaos expansion is used to represent the random eigenvalues and eigenvectors. However, the main concept of tracking the invariant subspace remains mostly independent of any such representation. The approach is successfully implemented in response prediction of a system with repeated natural frequencies. It is found that tracking only an invariant subspace could be sufficient to build a modal-based reduced-order model of the system. Copyright (C) 2012 John Wiley & Sons, Ltd.

Hydrodynamic instability and shear layer effect on the response of an acoustically excited laminar premixed flame

Combustion instabilities can cause serious problems which limit the operating envelope of low-emission lean premixed combustion systems. Predicting the onset of combustion instability requires a description of the unsteady heat release driving the instability, i.e., the heat release response transfer function of the system. This study focuses on the analysis of fully coupled two-way interactions between a disturbance field and a laminar premixed flame that incorporates gas expansion effects by solving the conservation equations of a compressible fluid. Results of the minimum and maximum flame front deflections are presented to underline the impact of the hydrodynamic instability on the flame and the shear layer effect on the initial flame front wrinkling which is increased at decreasing gas expansion. These phenomena influence the magnitude of the burning area and burning area rate response of the flame at lower frequency excitation more drastically than reduced-order model (ROM) predictions even for low temperature ratios. It is shown that the general trend of the flame response magnitudes can be well captured at higher frequency excitation, where stretch effects are dominant. The phase response is influenced by the DL mechanism, which cannot be captured by the ROM, and by the resulting discrepancy in the flame pocket formation and annihilation process at the flame tip. (C) 2014 The Combustion Institute. Published by Elsevier Inc. All rights reserved,

Dynamic Stall on Vertical Axis Wind Turbines

In this study the dynamics of flow over the blades of vertical axis wind turbines was investigated using a simplified periodic motion to uncover the fundamental flow physics and provide insight into the design of more efficient turbines. Time-resolved, two-dimensional velocity measurements were made with particle image velocimetry on a wing undergoing pitching and surging motion to mimic the flow on a turbine blade in a non-rotating frame. Dynamic stall prior to maximum angle of attack and a leading edge vortex development were identified in the phase-averaged flow field and captured by a simple model with five modes, including the first two harmonics of the pitch/surge frequency identified using the dynamic mode decomposition. Analysis of these modes identified vortical structures corresponding to both frequencies that led the separation and reattachment processes, while their phase relationship determined the evolution of the flow.

Detailed analysis of the leading edge vortex found multiple regimes of vortex development coupled to the time-varying flow field on the airfoil. The vortex was shown to grow on the airfoil for four convection times, before shedding and causing dynamic stall in agreement with 'optimal' vortex formation theory. Vortex shedding from the trailing edge was identified from instantaneous velocity fields prior to separation. This shedding was found to be in agreement with classical Strouhal frequency scaling and was removed by phase averaging, which indicates that it is not exactly coupled to the phase of the airfoil motion.

The flow field over an airfoil undergoing solely pitch motion was shown to develop similarly to the pitch/surge motion; however, flow separation took place earlier, corresponding to the earlier formation of the leading edge vortex. A similar reduced-order model to the pitch/surge case was developed, with similar vortical structures leading separation and reattachment; however, the relative phase lead of the separation mode, corresponding to earlier separation, necessitated that a third frequency to be incorporated into the reattachment mode to provide a relative lag in reattachment.

Finally, the results are returned to the rotating frame and the effects of each flow phenomena on the turbine are estimated, suggesting kinematic criteria for the design of improved turbines.

Compliance verification methodology for renewable generation integration. Application to island power grids

261 p.

Raman scattering characterization of Mn composition and strain in Ga1-xMnxSb/GaSb epitaxial layers

Raman scattering (RS) experiments have been performed for simultaneous determination of Mn composition and strain in Ga1-xMnxSb thin films grown on GaSb substrate by liquid phase epitaxy technique. The Raman spectra obtained from various Ga1-xMnxSb samples show only GaSb-like phonon modes whose frequency positions are found to have Mn compositional dependence. With the combination of epilayer strain model, RS and energy dispersive x-ray (EDX) experiments, the compositional dependence of GaSb-like LO phonon frequency is proposed both in strained and unstrained conditions. The proposed relationships are used to evaluate Mn composition and strain from the Ga1-xMnxSb samples. The results obtained from the RS data are found to be in good agreement with those determined independently by the EDX analysis. Furthermore, the frequency positions of MnSb-like phonon modes are suggested by reduced-mass model. (C) 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim

软件项目资源优化调度研究

信息技术的不断进步使得软件产品的应用领域不断扩大，同时软件产品的规模也在迅速膨胀。软件产品的开发模式已由最初的手工作坊式开发逐渐转变为大规模的工程化软件开发。这使得资源调度问题逐渐成为软件项目管理的核心研究内容之一。 软件项目与传统工业项目相比具有如下两个显著特点：其一是软件项目对人力资源能力的依赖性非常高；其二是软件项目在开发过程中具有的不确定性因素较多，也就是风险较高。这两个特点决定了传统工业调度方法不能很好地适用于软件项目管理，同时也为软件项目资源优化调度问题研究提出了新的挑战。需要根据软件项目的特点研究适合的资源调度方法为软件项目管理工作提供支持。 本文的研究工作旨在通过对软件项目的结构进行分析和描述，在建立软件项目核心要素模型的基础上，考虑软件项目高人力资源能力依赖性和高风险性两大特征，由人力资源能力和风险作为驱动因素，对软件项目中人力资源和项目缓冲两大核心资源进行优化分配和调度，以提高软件项目的资源利用效率和软件项目执行的稳定性。本文的主要贡献有： （1）建立了软件项目资源优化调度研究框架QMMT和项目核心要素模型PTHR。 QMMT研究框架由问题驱动(Question Driven)，模型描述(Model Description)，方法研究(Method Research)和工具验证(Tool Validation)四个模块构成。四个模块之间既存在顺序关系也存在信息反馈机制，框架具有良好的适应性和可扩展性。实践表明，QMMT研究框架对研究软件项目资源优化调度问题具有良好的指导作用。本文中涉及软件项目资源优化调度的多个研究问题均遵循QMMT研究框架。 通过对软件项目所包含的各个要素及要素之间的关系进行定义和描述，我们建立了软件项目核心要素模型PTHR。模型对软件项目的四个要素：项目（Project）、任务（Task）、人力资源（Human Resource）、风险（Risk）以及四个要素之间的关系进行了形式化定义和描述。PTHR模型涵盖了软件项目的核心要素并具有良好的可扩展性，可以为资源优化调度中具体问题的分析、算法的设计、流程的安排以及工具开发提供底层支持。PTHR模型是本文后续方法中相关系列子模型的基础模型。 （2） 提出了软件项目中任务-人员匹配的三维匹配模型3D-THM和基于3D-THM模型的任务人员优化分配方法。 任务人员匹配是人力资源调度的基础。3D-THM（3 Dimensional model for Task Human Matching）模型通过对人力资源的技术能力、性格能力和职业规划进行描述，以及对任务的技术能力需求、性格能力需求和职业规划需求进行描述，设定相应的多因素匹配算法，为任务-人员的全面优化匹配提供支持。实验表明，3D-THM模型较好的描绘了软件项目中任务-人员优化匹配问题，能够体现软件项目的高人力资源能力依赖性。模型实例化后所得到的匹配方法和相应的原型工具可为软件项目资源优化调度以及软件过程建模提供人员优化匹配支持，能够提高项目管理人员的工作效率，提升项目人员对任务分配的满意度。 （3） 提出了基于人员可用性的人力资源调度方法。 在对任务人员进行优化匹配的基础上，通过综合考虑人力资源能力和工作时间实现了基于人员可用性的人力资源调度方法。方法结合软件项目的结构特征，建立了任务人员可用性约束模型THACM（Task Human resources Availability Constraints Model）。基于THACM模型实现了在给定资源集合、任务集合下的人力资源自动分配和项目进度的自动安排。方法可有效解决采用矩阵组织结构的企业所面临的低资源可见性问题，协助其提高人力资源的利用效率。 （4）提出了基于任务优先级的抢占式人力资源调度方法PP-HAS。 在对人力资源可用性进行考虑的基础上，为了解决多项目环境下常见的资源冲突问题，我们提出了基于任务优先级的抢占式人力资源调度方法PP-HAS（Task Priority Based Preemptive Human Resource Scheduling Method）。方法首先建立了综合考虑进度、成本、质量三方面因素的基于价值的任务优先级模型VBTPM（Value Based Task Priority Model），将该任务优先级模型与过程Agent技术结合，通过设计支持抢占的人力资源调度流程，实现了多过程Agent协商下的人力资源优化调度。方法通过抢占和再计划实现了人力资源的动态高效利用，能够为资源冲突的解决以及项目的再计划工作提供决策支持。 （5）提出了风险驱动的软件项目缓冲资源分配方法。 项目缓冲的合理分配是降低风险对项目进度造成影响的重要手段。我们在软件项目资源调度方法中加入对风险因素的考量，基于软件项目中风险的特征，建立了简化的风险模型RRM（Reduced Risk Model）。基于RRM模型提出了风险驱动的项目缓冲分配方法，旨在软件项目的执行效率和稳定性二者之间进行权衡。模拟实验的结果表明，相对于传统关键链项目管理理论中尾部集中的项目缓冲分配方法，风险驱动的项目缓冲分配方法能够在确保对项目平均执行工期产生较小影响的同时，显著降低项目执行时计划变更的发生频率。该缓冲分配方法与项目模拟工具可以帮助项目经理确定合适的项目缓冲时间长度以及缓冲分配方案，进而提高软件项目计划的可信性和执行的稳定性。

悬索桥隧道式锚碇与围岩系统的破坏模式研究

The anchorages are unparalleled structures only in a suspension bridge, and as main bearing facilities, play an important role in connecting the superstructures and the ground. The tunnel anchorage, as one alternative type of the anchorages, has more advantages over its counterpart, the gravity anchorage. With the tunnel anchorages adopted, not only can surface excavation be reduced to protect the environment, and natural condition of the rock be utilized and potential bearing capacity of surrounding rock be mobilized to save engineering cost, but also the technological predominance of auxiliary engineering measures, such as prestressed concrete, anchoring piles, rock anchors and collar beam between the two separated anchorages, can be easily cooperated to work together harmoniously under the circumstances of poor rock quality. There are plentiful high mountains and deep canyons in west part of China, and long-span bridge construction is inevitably encountered in order to realize leapfrogging development of the transportation infrastructure. Western mountainous areas usually possess the conditions for constructing tunnel anchorages, and therefore, the tunnel anchorages, which are conformed to the conception of resource conservative and sustainable society, extremely have application and popularization value in western underdeveloped region. The scientific and technological problem about the design, construction and operation of tunnel anchorages should be further investigated. Combining the engineering of western tunnel anchorages for the Balinghe Suspension Bridge, this paper probed into the survey method and in-situ test method for tunnel anchorages, scientific rock quality evaluation of surrounding rock to provide reasonable physical and mechanical parameters for design, construction and operation of tunnel anchorages, bearing capacity estimation for tunnel anchorage, deformation prediction of the anchorage-rockmass system, tunnel-anchorage slope stability analysis and the evaluation of excavation stability and degree of safety of the anchorage tunnel. The following outcomes were obtained: 1. Materials of tunnel anchorages of suspension bridge built (and in progress) at home and abroad were systematically sorted out, with the engineering geological condition and geomechanical property of surrounding rock around the anchorage tunnel, the design size of anchorages and the construction method of anchorage tunnel paid more emphasis on, to unveil the internal relationship between the engineering geological conditions of surrounding rock and the design size and axis angle of anchorages and provide references for future design, construction and study of tunnel anchorages. 2. Physical and mechanical parameters were recommended based on three domestic and foreign methods of rock quality evaluation. 3. In-situ tests, adopting the back-thrust method, of two kinds of reduced scale model, 1/30 and 1/20, for the tunnel anchorages were conducted in the declining exploration drift with rock mass at the test depth being the same as surrounding rock around real anchorages, and reliable field rockmass displacement data were acquired. Attenuation relation between the increment of distance from the anchorage and the decrement of rockmass displacement under maximum test load, and influential scope suffered by anchorage load were obtained. 4. Using similarity theory, the magnitude of real anchorage and rockmass displacement under design load and degree of safety of the anchorage system were deduced. Furthermore, inversion analysis to deformation modulus of slightly weathered dolomite rock, the surrounding rock of anchorage tunnel, was performed by the means of numerical simulation. 5. The influential law of the geometrical size to the limit bearing capacity of tunnel anchorage was studied. 6. Based on engineering geological survey data, accounting for the combination of strata layer and adverse discontinuities, the failure patterns of tunnel anchorage slope were divided into three modes: sliding of splay saddle pier slope, superficial-layer slippage, and deep-layer slippage. Using virtual work principle and taking anchorage load in account, the stability of the three kinds of failure patterns were analyzed in detail. 7. The step-by-step excavation of anchorage tunnel, the numerical overload and the staged decrement of rock strength parameters were numerically simulated to evaluate the excavation stability of surrounding rock around anchorage tunnel, the overload performance of tunnel anchorage, and the safety margin of strength parameters of the surrounding rock.

Appropriate complexity for the prediction of coastal and estuarine geomorphic behaviour at decadal to centennial scales

Coastal and estuarine landforms provide a physical template that not only accommodates diverse ecosystem functions and human activities, but also mediates flood and erosion risks that are expected to increase with climate change. In this paper, we explore some of the issues associated with the conceptualisation and modelling of coastal morphological change at time and space scales relevant to managers and policy makers. Firstly, we revisit the question of how to define the most appropriate scales at which to seek quantitative predictions of landform change within an age defined by human interference with natural sediment systems and by the prospect of significant changes in climate and ocean forcing. Secondly, we consider the theoretical bases and conceptual frameworks for determining which processes are most important at a given scale of interest and the related problem of how to translate this understanding into models that are computationally feasible, retain a sound physical basis and demonstrate useful predictive skill. In particular, we explore the limitations of a primary scale approach and the extent to which these can be resolved with reference to the concept of the coastal tract and application of systems theory. Thirdly, we consider the importance of different styles of landform change and the need to resolve not only incremental evolution of morphology but also changes in the qualitative dynamics of a system and/or its gross morphological configuration. The extreme complexity and spatially distributed nature of landform systems means that quantitative prediction of future changes must necessarily be approached through mechanistic modelling of some form or another. Geomorphology has increasingly embraced so-called ‘reduced complexity’ models as a means of moving from an essentially reductionist focus on the mechanics of sediment transport towards a more synthesist view of landform evolution. However, there is little consensus on exactly what constitutes a reduced complexity model and the term itself is both misleading and, arguably, unhelpful. Accordingly, we synthesise a set of requirements for what might be termed ‘appropriate complexity modelling’ of quantitative coastal morphological change at scales commensurate with contemporary management and policy-making requirements: 1) The system being studied must be bounded with reference to the time and space scales at which behaviours of interest emerge and/or scientific or management problems arise; 2) model complexity and comprehensiveness must be appropriate to the problem at hand; 3) modellers should seek a priori insights into what kind of behaviours are likely to be evident at the scale of interest and the extent to which the behavioural validity of a model may be constrained by its underlying assumptions and its comprehensiveness; 4) informed by qualitative insights into likely dynamic behaviour, models should then be formulated with a view to resolving critical state changes; and 5) meso-scale modelling of coastal morphological change should reflect critically on the role of modelling and its relation to the observable world.