Modelling of transient heat conduction with diffuse interface methods

We present a survey on different numerical interpolation schemes used for two-phase transient heat conduction problems in the context of interface capturing phase-field methods. Examples are general transport problems in the context of diffuse interface methods with a non-equal heat conductivity in normal and tangential directions to the interface. We extend the tonsorial approach recently published by Nicoli M et al (2011 Phys. Rev. E 84 1-6) to the general three-dimensional (3D) transient evolution equations. Validations for one-dimensional, two-dimensional and 3D transient test cases are provided, and the results are in good agreement with analytical and numerical reference solutions.

Development and Evaluation of Statistical Downscaling Models for Monthly Precipitation

Several statistical downscaling models have been developed in the past couple of decades to assess the hydrologic impacts of climate change by projecting the station-scale hydrological variables from large-scale atmospheric variables simulated by general circulation models (GCMs). This paper presents and compares different statistical downscaling models that use multiple linear regression (MLR), positive coefficient regression (PCR), stepwise regression (SR), and support vector machine (SVM) techniques for estimating monthly rainfall amounts in the state of Florida. Mean sea level pressure, air temperature, geopotential height, specific humidity, U wind, and V wind are used as the explanatory variables/predictors in the downscaling models. Data for these variables are obtained from the National Centers for Environmental Prediction-National Center for Atmospheric Research (NCEP-NCAR) reanalysis dataset and the Canadian Centre for Climate Modelling and Analysis (CCCma) Coupled Global Climate Model, version 3 (CGCM3) GCM simulations. The principal component analysis (PCA) and fuzzy c-means clustering method (FCM) are used as part of downscaling model to reduce the dimensionality of the dataset and identify the clusters in the data, respectively. Evaluation of the performances of the models using different error and statistical measures indicates that the SVM-based model performed better than all the other models in reproducing most monthly rainfall statistics at 18 sites. Output from the third-generation CGCM3 GCM for the A1B scenario was used for future projections. For the projection period 2001-10, MLR was used to relate variables at the GCM and NCEP grid scales. Use of MLR in linking the predictor variables at the GCM and NCEP grid scales yielded better reproduction of monthly rainfall statistics at most of the stations (12 out of 18) compared to those by spatial interpolation technique used in earlier studies.

Evaluating Methods to Predict Streamflow at Ungauged Sites using Regional Flow Duration Curves: A Case Study

Precise information on streamflows is of major importance for planning and monitoring of water resources schemes related to hydro power, water supply, irrigation, flood control, and for maintaining ecosystem. Engineers encounter challenges when streamflow data are either unavailable or inadequate at target locations. To address these challenges, there have been efforts to develop methodologies that facilitate prediction of streamflow at ungauged sites. Conventionally, time intensive and data exhaustive rainfall-runoff models are used to arrive at streamflow at ungauged sites. Most recent studies show improved methods based on regionalization using Flow Duration Curves (FDCs). A FDC is a graphical representation of streamflow variability, which is a plot between streamflow values and their corresponding exceedance probabilities that are determined using a plotting position formula. It provides information on the percentage of time any specified magnitude of streamflow is equaled or exceeded. The present study assesses the effectiveness of two methods to predict streamflow at ungauged sites by application to catchments in Mahanadi river basin, India. The methods considered are (i) Regional flow duration curve method, and (ii) Area Ratio method. The first method involves (a) the development of regression relationships between percentile flows and attributes of catchments in the study area, (b) use of the relationships to construct regional FDC for the ungauged site, and (c) use of a spatial interpolation technique to decode information in FDC to construct streamflow time series for the ungauged site. Area ratio method is conventionally used to transfer streamflow related information from gauged sites to ungauged sites. Attributes that have been considered for the analysis include variables representing hydrology, climatology, topography, land-use/land- cover and soil properties corresponding to catchments in the study area. Effectiveness of the presented methods is assessed using jack knife cross-validation. Conclusions based on the study are presented and discussed. (C) 2015 The Authors. Published by Elsevier B.V.

A Family of Domains Associated with mu-Synthesis

We introduce a family of domains-which we call the -quotients-associated with an aspect of -synthesis. We show that the natural association that the symmetrized polydisc has with the corresponding spectral unit ball is also exhibited by the -quotient and its associated unit `` -ball''. Here, is the structured singular value for the case Specifically: we show that, for such an E, the Nevanlinna-Pick interpolation problem with matricial data in a unit `` -ball'', and in general position in a precise sense, is equivalent to a Nevanlinna-Pick interpolation problem for the associated -quotient. Along the way, we present some characterizations for the -quotients.

Timing Recovery Algorithms and Architectures for 2-D Magnetic Recording Systems

We investigate the problem of timing recovery for 2-D magnetic recording (TDMR) channels. We develop a timing error model for TDMR channel considering the phase and frequency offsets with noise. We propose a 2-D data-aided phase-locked loop (PLL) architecture for tracking variations in the position and movement of the read head in the down-track and cross-track directions and analyze the convergence of the algorithm under non-separable timing errors. We further develop a 2-D interpolation-based timing recovery scheme that works in conjunction with the 2-D PLL. We quantify the efficiency of our proposed algorithms by simulations over a 2-D magnetic recording channel with timing errors.

Meshless Local Petrov-Galerkin Method for Rotating Euler-Bernoulli Beam

Free vibration problem of a rotating Euler-Bernoulli beam is solved with a truly meshless local Petrov-Galerkin method. Radial basis function and summation of two radial basis functions are used for interpolation. Radial basis function satisfies the Kronecker delta property and makes it simpler to apply the essential boundary conditions. Interpolation with summation of two radial basis functions increases the node carrying capacity within the sub-domain of the trial function and higher natural frequencies can be computed by selecting the complete domain as a sub-domain of the trial function. The mass and stiffness matrices are derived and numerical results for frequencies are obtained for a fixed-free beam and hinged-free beam simulating hingeless and articulated helicopter blades. Stiffness and mass distribution suitable for wind turbine blades are also considered. Results show an accurate match with existing literature.

On the detection of human influence in extreme precipitation over India

Climate change is expected to influence extreme precipitation which in turn might affect risks of pluvial flooding. Recent studies on extreme rainfall over India vary in their definition of extremes, scales of analyses and conclusions about nature of changes in such extremes. Fingerprint-based detection and attribution (D&A) offer a formal way of investigating the presence of anthropogenic signals in hydroclimatic observations. There have been recent efforts to quantify human effects in the components of the hydrologic cycle at large scales, including precipitation extremes. This study conducts a D&A analysis on precipitation extremes over India, considering both univariate and multivariate fingerprints, using a standardized probability-based index (SPI) from annual maximum one-day (RX1D) and five-day accumulated (RX5D) rainfall. The pattern-correlation based fingerprint method is used for the D&A analysis. Transformation of annual extreme values to SPI and subsequent interpolation to coarser grids are carried out to facilitate comparison between observations and model simulations. Our results show that in spite of employing these methods to address scale and physical processes mismatch between observed and model simulated extremes, attributing changes in regional extreme precipitation to anthropogenic climate change is difficult. At very high (95%) confidence, no signals are detected for RX1D, while for the RX5D and multivariate cases only the anthropogenic (ANT) signal is detected, though the fingerprints are in general found to be noisy. The findings indicate that model simulations may underestimate regional climate system responses to increasing human forcings for extremes, and though anthropogenic factors may have a role to play in causing changes in extreme precipitation, their detection is difficult at regional scales and not statistically significant. (C) 2015 Elsevier B.V. All rights reserved.

Uncertainty Quantification in Ultrasonic Wave Based Detection of Delamination in Composite Structures

In this paper we consider the problem of guided wave scattering from delamination in laminated composite and further the problem of estimating delamination size and layer-wise location from the guided wave measurement. Damage location and region/size can be estimated from time of flight and wave packet spread, whereas depth information can be obtained from wavenumber modulation in the carrier packet. The key challenge is that these information are highly sensitive to various uncertainties. Variation in reflected and transmitted wave amplitude in a bar due to boundary/interface uncertainty is studied to illustrate such effect. Effect of uncertainty in material parameters on the time of flight are estimated for longitudinal wave propagation. To evaluate the effect of uncertainty in delamination detection, we employ a time domain spectral finite element (tSFEM) scheme where wave propagation is modeled using higher-order interpolation with shape function have spectral convergence properties. A laminated composite beam with layer-wise placement of delamination is considered in the simulation. Scattering due to the presence of delamination is analyzed. For a single delamination, two identical waveforms are created at the two fronts of the delamination, whereas waves in the two sub-laminates create two independent waveforms with different wavelengths. Scattering due to multiple delaminations in composite beam is studied.

A mode-shape classification technique for robust speech rate estimation and syllable nuclei detection

Acoustic feature based speech (syllable) rate estimation and syllable nuclei detection are important problems in automatic speech recognition (ASR), computer assisted language learning (CALL) and fluency analysis. A typical solution for both the problems consists of two stages. The first stage involves computing a short-time feature contour such that most of the peaks of the contour correspond to the syllabic nuclei. In the second stage, the peaks corresponding to the syllable nuclei are detected. In this work, instead of the peak detection, we perform a mode-shape classification, which is formulated as a supervised binary classification problem - mode-shapes representing the syllabic nuclei as one class and remaining as the other. We use the temporal correlation and selected sub-band correlation (TCSSBC) feature contour and the mode-shapes in the TCSSBC feature contour are converted into a set of feature vectors using an interpolation technique. A support vector machine classifier is used for the classification. Experiments are performed separately using Switchboard, TIMIT and CTIMIT corpora in a five-fold cross validation setup. The average correlation coefficients for the syllable rate estimation turn out to be 0.6761, 0.6928 and 0.3604 for three corpora respectively, which outperform those obtained by the best of the existing peak detection techniques. Similarly, the average F-scores (syllable level) for the syllable nuclei detection are 0.8917, 0.8200 and 0.7637 for three corpora respectively. (C) 2016 Elsevier B.V. All rights reserved.

Computation of Viscous Incompressible Flow using Pressure Correction Method on Unstructured Chimera Grid

In this paper, a pressure correction algorithm for computing incompressible flows is modified and implemented on unstructured Chimera grid. Schwarz method is used to couple the solutions of different sub-domains. A new interpolation to ensure consistency between primary variables and auxiliary variables is proposed. Other important issues such as global mass conservation and order of accuracy in the interpolations are also discussed. Two numerical simulations are successfully performed. They include one steady case, the lid-driven cavity and one unsteady case, the flow around a circular cylinder. The results demonstrate a very good performance of the proposed scheme on unstructured Chimera grids. It prevents the decoupling of pressure field in the overlapping region and requires only little modification to the existing unstructured Navier–Stokes (NS) solver. The numerical experiments show the reliability and potential of this method in applying to practical problems.

Transonic Aeroelastic Numerical Simulation in Aeronautical Engineering

A lower-upper symmetric Gauss-Seidel (LU-SGS) subiteration scheme is constructed for time-marching of the fluid equations. The Harten-Lax-van Leer-Einfeldt-Wada (HLLEW) scheme is used for the spatial discretization. The same subiteration formulation is applied directly to the structural equations of motion in generalized coordinates. Through subiteration between the fluid and structural equations, a fully implicit aeroelastic solver is obtained for the numerical simulation of fluid/structure interaction. To improve the ability for application to complex configurations, a multiblock grid is used for the flow field calculation and transfinite interpolation (TFI) is employed for the adaptive moving grid deformation. The infinite plate spline (IPS) and the principal of virtual work are utilized for the data transformation between the fluid and structure. The developed code was first validated through the comparison of experimental and computational results for the AGARD 445.6 standard aeroelastic wing. Then, the flutter character of a tail wing with control surface was analyzed. Finally, flutter boundaries of a complex aircraft configuration were predicted.

飞行器跨声速气动弹性数值分析

将流体和结构运动方程分别构造为含子迭代的计算格式，发展了一种紧耦合气动弹性分析方法．其中流体计算的空间离散采用改进的HLLEW（Harten-Lax-van Leer-Einfeldt-Wada）格式．TFI（transfinite interpolation）方法用于生成随结构变形的自适应多块动网格．利用所发展的方法，对一翼-身-尾气动外形，数值预测了马赫数在0．3-1．3范围内的气动颤振边界．并详细研究了时间步长、子迭代步数、初始流场、耦合方法、疏密网格对颤振计算结果的影响．

Computer generated hologram from point cloud using graphics processor

Computer generated holography is an extremely demanding and complex task when it comes to providing realistic reconstructions with full parallax, occlusion, and shadowing. We present an algorithm designed for data-parallel computing on modern graphics processing units to alleviate the computational burden. We apply Gaussian interpolation to create a continuous surface representation from discrete input object points. The algorithm maintains a potential occluder list for each individual hologram plane sample to keep the number of visibility tests to a minimum.We experimented with two approximations that simplify and accelerate occlusion computation. It is observed that letting several neighboring hologramplane samples share visibility information on object points leads to significantly faster computation without causing noticeable artifacts in the reconstructed images. Computing a reduced sample set via nonuniform sampling is also found to be an effective acceleration technique. © 2009 Optical Society of America.

Automatic transcription of conversational telephone speech

This paper discusses the Cambridge University HTK (CU-HTK) system for the automatic transcription of conversational telephone speech. A detailed discussion of the most important techniques in front-end processing, acoustic modeling and model training, language and pronunciation modeling are presented. These include the use of conversation side based cepstral normalization, vocal tract length normalization, heteroscedastic linear discriminant analysis for feature projection, minimum phone error training and speaker adaptive training, lattice-based model adaptation, confusion network based decoding and confidence score estimation, pronunciation selection, language model interpolation, and class based language models. The transcription system developed for participation in the 2002 NIST Rich Transcription evaluations of English conversational telephone speech data is presented in detail. In this evaluation the CU-HTK system gave an overall word error rate of 23.9%, which was the best performance by a statistically significant margin. Further details on the derivation of faster systems with moderate performance degradation are discussed in the context of the 2002 CU-HTK 10 × RT conversational speech transcription system. © 2005 IEEE.

A Multimoment Finite-Volume Shallow-Water Model On The Yin-Yang Overset Spherical Grid

A numerical model for shallow-water equations has been built and tested on the Yin-Yang overset spherical grid. A high-order multimoment finite-volume method is used for the spatial discretization in which two kinds of so-called moments of the physical field [i.e., the volume integrated average ( VIA) and the point value (PV)] are treated as the model variables and updated separately in time. In the present model, the PV is computed by the semi-implicit semi-Lagrangian formulation, whereas the VIA is predicted in time via a flux-based finite-volume method and is numerically conserved on each component grid. The concept of including an extra moment (i.e., the volume-integrated value) to enforce the numerical conservativeness provides a general methodology and applies to the existing semi-implicit semi-Lagrangian formulations. Based on both VIA and PV, the high-order interpolation reconstruction can only be done over a single grid cell, which then minimizes the overlapping zone between the Yin and Yang components and effectively reduces the numerical errors introduced in the interpolation required to communicate the data between the two components. The present model completely gets around the singularity and grid convergence in the polar regions of the conventional longitude-latitude grid. Being an issue demanding further investigation, the high-order interpolation across the overlapping region of the Yin-Yang grid in the current model does not rigorously guarantee the numerical conservativeness. Nevertheless, these numerical tests show that the global conservation error in the present model is negligibly small. The model has competitive accuracy and efficiency.