A matrix-calculation-based algorithm for numerical change propagation analysis

Engineering changes (ECs) are raised throughout the lifecycle of engineering products. A single change to one component produces knock-on effects on others necessitating additional changes. This change propagation significantly affects the development time and cost and determines the product's success. Predicting and managing such ECs is, thus, essential to companies. Some prediction tools model change propagation by algorithms, whereof a subgroup is numerical. Current numerical change propagation algorithms either do not account for the exclusion of cyclic propagation paths or are based on exhaustive searching methods. This paper presents a new matrix-calculation-based algorithm which can be applied directly to a numerical product model to analyze change propagation and support change prediction. The algorithm applies matrix multiplications on mutations of a given design structure matrix accounting for the exclusion of self-dependences and cyclic propagation paths and delivers the same results as the exhaustive search-based Trail Counting algorithm. Despite its factorial time complexity, the algorithm proves advantageous because of its straightforward matrix-based calculations which avoid exhaustive searching. Thereby, the algorithm can be implemented in established numerical programs such as Microsoft Excel which promise a wider application of the tools within and across companies along with better familiarity, usability, practicality, security, and robustness. © 1988-2012 IEEE.

Effects of turbulent reynolds number on the performance of algebraic flame surface density models for large eddy simulation in the thin reaction zones regime: A direct numerical simulation analysis

A direct numerical simulation (DNS) database of freely propagating statistically planar turbulent premixed flames with a range of different turbulent Reynolds numbers has been used to assess the performance of algebraic flame surface density (FSD) models based on a fractal representation of the flame wrinkling factor. The turbulent Reynolds number Ret has been varied by modifying the Karlovitz number Ka and the Damköhler number Da independently of each other in such a way that the flames remain within the thin reaction zones regime. It has been found that the turbulent Reynolds number and the Karlovitz number both have a significant influence on the fractal dimension, which is found to increase with increasing Ret and Ka before reaching an asymptotic value for large values of Ret and Ka. A parameterisation of the fractal dimension is presented in which the effects of the Reynolds and the Karlovitz numbers are explicitly taken into account. By contrast, the inner cut-off scale normalised by the Zel'dovich flame thickness ηi/δz does not exhibit any significant dependence on Ret for the cases considered here. The performance of several algebraic FSD models has been assessed based on various criteria. Most of the algebraic models show a deterioration in performance with increasing the LES filter width. © 2012 Mohit Katragadda et al.

Local strain rate and curvature dependences of scalar dissipation rate transport in turbulent premixed flames: A direct numerical simulation analysis

The statistical behaviours of the instantaneous scalar dissipation rate Nc of reaction progress variable c in turbulent premixed flames have been analysed based on three-dimensional direct numerical simulation data of freely propagating statistically planar flame and V-flame configurations with different turbulent Reynolds number Ret. The statistical behaviours of N c and different terms of its transport equation for planar and V-flames are found to be qualitatively similar. The mean contribution of the density-variation term T1 is positive, whereas the molecular dissipation term (-D2) acts as a leading order sink. The mean contribution of the strain rate term T2 is predominantly negative for the cases considered here. The mean reaction rate contribution T3 is positive (negative) towards the unburned (burned) gas side of the flame, whereas the mean contribution of the diffusivity gradient term (D) assumes negative (positive) values towards the unburned (burned) gas side. The local statistical behaviours of Nc, T1, T2, T 3, (-D2), and f(D) have been analysed in terms of their marginal probability density functions (pdfs) and their joint pdfs with local tangential strain rate aT and curvature km. Detailed physical explanations have been provided for the observed behaviour. © 2014 Y. Gao et al.

Numerical and theoretical analysis of the crosstalk in linear optical amplifiers

The dynamic characteristics, including the crosstalk and relaxation oscillation, of linear optical amplifiers (LOAs) are investigated by small-signal analysis under an averaging carrier density approximation and compared with the results of numerical simulation. The good agreement between the numerical simulation and the small-signal analysis indicated the averaging carrier density is an appropriate approximation for analyzing LOAs. Theoretical analyzes also show that the dynamic properties of the vertical laser fields dominate the dynamic performance of LOAs. Based on the small-signal analysis, a concise equation for the crosstalk under high bit rate was derived, which can be applied to measure the differential gain of LOAs.

Symmetry analysis and numerical simulation of mode characteristics for equilateral-polygonal optical microresonators

Mode characteristics for two-dimensional equilateral-polygonal microresonators are investigated based on symmetry analysis and finite-difference time-domain numerical simulation. The symmetries of the resonators can be described by the point group C-Nv, accordingly, the confined modes in these resonators can be classified into irreducible representations of the point group C-Nv. Compared with circular resonators, the modes in equilateral-polygonal resonators have different characteristics due to the break of symmetries, such as the split of double-degenerate modes, high field intensity in the center region, and anomalous traveling-wave modes, which should be considered in the designs of the polygonal resonator microlasers or optical add-drop filters.

Analysis and numerical simulation of eigenmode characteristics for semiconductor lasers with an equilateral triangle micro-resonator

The eigenmodes confined in the equilateral triangle resonator (ETR) are analyzed by deriving the eigenvalues and the mode field distributions and by the finite difference time domain (FDTD) technique. The analytical results show that the one-period-length for the mode light rays inside the ETR is the perimeter of the ETR, and the number of transverse modes is limited by the condition of total internal reflection. In addition, the sum of the longitudinal mode index and the transverse mode index should be an even number, which limits the number of confined modes again. Based on the FDTD technique and the Pade approximation, we calculate the mode resonant frequencies and the quality factors from the local maximum and the width of the spectral distribution of the intensity The numerical results of mode frequencies agree very well with the analytical results, and the quality factor of the fundamental mode is usually higher than that of the higher order transverse modes. The results show that the ETR is suitable to realize single-made operation as semiconductor microcavity lasers.

Descriptive Simulation: Combining Symbolic and Numerical Methods in the Analysis of Chemical Reaction Mechanisms

The Kineticist's Workbench is a computer program currently under development whose purpose is to help chemists understand, analyze, and simplify complex chemical reaction mechanisms. This paper discusses one module of the program that numerically simulates mechanisms and constructs qualitative descriptions of the simulation results. These descriptions are given in terms that are meaningful to the working chemist (e.g., steady states, stable oscillations, and so on); and the descriptions (as well as the data structures used to construct them) are accessible as input to other programs.

Numerical optimisation techniques applied to evacuation analysis

A common problem faced by fire safety engineers in the field of evacuation analysis concerns the optimal design of an arbitrarily complex structure in order to minimise evacuation times. How does the engineer determine the best solution? In this study we introduce the concept of numerical optimisation techniques to address this problem. The study makes user of the buildingEXODUS evacuation model coupled with classical optimisation theory including Design of Experiments (DoE) and Response Surface Models (RSM). We demonstrate the technique using a relatively simple problem of determining the optimal location for a single exit in a square room.

Numerical modelling of laboratory tests for TBM cutter indentation process analysis

Considering that TBMs are nowadays used for long Trans-Alpine tunnels, the
understanding of rock breaking and chipping due to TBM cutter disks mechanism, for deep tunnelling operations, becomes very interesting. In this paper, the results from carried out laboratory tests that simulate the disk cutter action at the rock tunnel face by means of an indentation tool, acting on a rock
specimen with proper size, and the related three-dimensional and two-dimensional numerical modelling are proposed. The developed numerical models simulate the different test conditions (applied load, boundary conditions) allowing the analysis of the stresses distributions along possible breaking planes.
The influence of a confinement-free area on one side of the specimen, simulating the formation of a groove near the tool, is pointed out.
The obtained results from numerical modelling put in evidence a satisfactory agreement with the experimental observations.

Development and Validation of a Procedure for Numerical Vibration Analysis of an Oscillating Wave Surge Converter

During extreme sea states so called impact events can be observed on the wave energy converter Oyster. In small scale experimental tests these impact events cause high frequency signals in the measured load which decrease confidence in the data obtained. These loads depend on the structural dynamics of the model. Amplification of the loads can occur and is transferred through the structure from the point of impact to the load cell located in the foundation. Since the determination of design data and load cases for Wave Energy Converters originate from scale experiments, this lack of confidence has a direct effect on the development.

Numerical vibration analysis is a valuable tool in the research of the structural load response of Oyster to impact events, but must take into account the effect of the surrounding water. This can be done efficiently by adding an added mass distribution, computed with a linearised potential boundary element method. This paper presents the development and validation of a numerical procedure, which couples the OpenSource boundary element code NEMOH with the Finite Element Analysis tool CodeAster. Numerical results of the natural frequencies and mode shapes of the structure under the influence of added mass due to specific structural modes are compared with experimental results.

Numerical tools for isogeometric analysis in the nonlinear regime

The present work deals with the development of robust numerical tools for Isogeometric Analysis suitable for problems of solid mechanics in the nonlinear regime. To that end, a new solid-shell element, based on the Assumed Natural Strain method, is proposed for the analysis of thin shell-like structures. The formulation is extensively validated using a set of well-known benchmark problems available in the literature, in both linear and nonlinear (geometric and material) regimes. It is also proposed an alternative formulation which is focused on the alleviation of the volumetric locking pathology in linear elastic problems. In addition, an introductory study in the field of contact mechanics, in the context of Isogeometric Analysis, is also presented, with special focus on the implementation of a the Point-to-Segment algorithm. All the methodologies presented in the current work were implemented in a in-house code, together with several pre- and post-processing tools. In addition, user subroutines for the commercial software Abaqus were also implemented.

Identification of descendants of an extinct bovine population from the Algarve region of Portugal using numerical taxonomy analysis of morphological traits

The morphology of a sample of four bulls and 43 cows, presumed to be descendants of the extinct cattle breed ‘Algarvia’ (AG), was used to assign their relationship with animals from other Portuguese autochthonous breeds – Arouquesa (AR), Barrosa˜ (BA), Cachena (CA), Marinhoa (MA), Maronesa (MO), Minhota (MN), Mirandesa (MI), (only bulls), Alentejana (AL), Garvonesa (GA), Mertolenga (ME) and Preta (PR). Standard numerical taxonomic methods were applied to a set of 183 (cows) and 170 (bulls) traits, to derive average pairwise taxonomic distances among the sample of 257 cows and 76 bulls. Distance coefficients (morphological index of distance) ranged from 0.22 to 2.62 (cows) and from 0.49 to 2.13 (bulls). Unweighted pair group method using arithmetic averages (UPGMA)-based phenograms and a principal coordinate analysis showed that bulls were highly clustered and cows showed a tendency to cluster according to their geographical and breed origin. The AG population grouped together with GA, AL, ME and MN breeds in the Red Convex group. The average taxonomic distance among breeds was 1.02, the highest being 1.39 (ME versus BA) and the lowest being 0.64 (MA versus AR). The approach allowed for the identification of a phenotypically differentiated set of animals, comprising 19 cows and four bulls representative of the AG breed, and which can be targeted in further studies aiming at the recovery of this extinct breed.