Optimal distance from high power laser diode to cylindrical microlens in a coupling system

In a practical coupling system, a cylindrical microlens is used to collimate the emission of a high powerlaser diode (LD) in the dimension perpendicular to the junction plane. Using passive alignment, the LD isplaced in the focus of the cylindrical microlens generally, regardless of the performance of the multimodeoptical fiber and the LD. In this paper, a more complete analysis is arrived at by ray-tracing technique,by which the angle θ of the ray after refraction is computed as a function of the angle θo of the ray beforerefraction. The focus of the cylindrical microlens is not always the optimal position of the LD. In fact, inorder to achieve a higher coupling efficiency, the optimal distance from the LD to the cylindrical microlensis dependent on not only the radius R and the index of refraction n of the cylindrical microlens, but alsothe divergence angle of the LD in the dimension perpendicular to the junction plane and the numericalaperture (NA) of the multimode optical fiber. The results of this discussion are in good agreement withexperimental results.

Computer simulation of the damage evolution of fiber-reinforced polypropylene-matrix composites with matrix defects

The damage evolution of fiber-reinforced polypropylene-matrix composites with matrix defects was studied via a Monte Carlo technique combined with a finite element method. A finite element model was constructed to predict the effects of various matrix defect shapes on the stress distributions. The results indicated that a small matrix defect had almost no effect on fiber stress distributions other than interfacial shear stress distributions. Then, a finite element model with a statistical distribution of the fiber strength was constructed to investigate the influences of the spatial distribution and the volume fraction of matrix defects on composite failure. The results showed that it was accurate to use the shear-lag models and Green's function methods to predict the tensile strength of composites even though the axial stresses in the matrix were neglected.

A biofuel cell with enhanced power output by grape juice

Glucose oxidase and laccase immobilized at multiwalled carbon nanotubes-ionic liquid gel modified electrodes are used as the catalysts of anode and cathode of biofuel cells (BFCs), respectively. The BFC based on glucose and air is proposed. When ferrocene monocarboxylic acid is adopted as the mediator of anode, the power output of the BFC is ca. 4.1 mu W (power density ca. 10.0 mu W cm(-2)), which is higher than the value of 2.7 mu W (power density ca. 6.6 mu W cm(-2)) by taking ferrocene dicarboxylic acid as the mediator. This implies that the mediator with formal potential closing to that of the enzyme does improve the power output. Furthermore, the power output of the BFC is greatly improved by taking grape juice as the fuel of anode rather than glucose. This system also indicates that grape juice as a fuel of the BFC not only is feasible and can also enhances the power output of the BFCs. Besides, it greatly lowers the cost and simplifies the preparation procedure of the BFCs, making the BFC towards "green" bioenergy.

Effective dielectric responses of graded composites of the particles with a general power-law gradation profile

The effective dielectric response of graded spherical composites having general power-law gradient inclusions is investigated under a uniform applied electric field, where the dielectric gradation profile of the spherical inclusions is modeled by the equation epsilon(i) (r) = c(b+r)(k). Analytical solutions of the local electrical potentials are derived in terms of hyper-geometric function and the effective dielectric response of the graded composites is predicted in the dilute limit. From our result, the local potentials of graded spherical composites having both simple power-law dielectric profile epsilon(i)(r) = cr(k) and linear dielectric profile epsilon(i) (r) = c(b+r) are derived exactly by taking the limits b --> 0 and k --> 1, respectively. In the dilute limit, our exact result is used to test the validity of differential effective dipole approximation (DEDA) for estimating the effective response of graded spherical composites, and it is shown that the DEDA is in excellent agreement with exact result. (C) 2005 Elsevier B.V. All rights reserved.

Effective properties of graded elliptical cylindrical composites

The effective property has been investigated theoretically in graded elliptical cylindrical composite's consisting of inhomogeneous graded elliptical cylinders and an isotropic matrix under external uniform electric field. As a theoretical model, the dielectric gradient profile in the elliptical cylinder is modeled by a power-law function of short semi-axis variable parameter (xi(2) - 1) in the elliptical cylindrical coordinates, namely epsilon(i)(xi) = c(k) (xi(2) - 1)(k), where c(k) and k are the parameters, and xi is the long semi-axis space variable in an elliptical cylindrical inclusion region. In the dilute limit, the local analytical potentials in inclusion and matrix regions are derived exactly by means of the hyper-geometric function, and the formulas are given for estimating the effective dielectric responses under the external lfield along (x) over cap- and (y) over cap -directions, respectively. Furthermore, we have demonstrated that our effective response formulas can be reduced to the well-known results of homogeneous isotropic elliptical cylindrical composites if we take the limit k -> 0 in graded elliptical cylindrical composites. (c) 2006 Elsevier B.V. All rights reserved.

大别山岩石圈有效弹性厚度初步研究

As powerful tools to study the lithosphere dynamics, the effective elastic thickness (Te) as well as the envelope of yielding stress of lithosphere have been attracted great attention of geoscientists in the past thirty years. The oceanic lithosphere, contrary to the continental lithosphere, has more fruits for its simple structures and evolution process. In continent, the lithosphere commonly is complex and variable in the rheological, thermal structures, and has a complicated history. Therefore, the application of the effective elastic thickness in continent is still a subject to learn in a long time. Te, with the definition of the thickness of an elastic plate in theory flexured by the equal benging of the real stress in the lithosphere plate (Turcotte, 1982), marks the depth of transition between elastic and fluid behaviors of rocks subjected to stress exceeding 100 MPa over the geological timescales (McNutt, 1990). There are three methods often adapted: admittance or isostatic response function, coherence and forwarding. In principle, the models of Te consist of thermal-rheological, non-linear Maxwell, non-linear work hardening and rheological layered models. There is a tentative knowledge of Te that it is affected by the following factors: crustal thickness, crust-mantle decoupling, plate bending, boundary conditions of plate (end forces and bending moments), stress state, sedimentary layer, faulting effect, variation in the mountain belts' strike, foreland basin, inheritance of tectonic evolution, convection of mantle, seismic depth and lithosphere strength. In this thesis, the author introduces the geological sketch of the Dabie collisional orogenic belt and the Hefei Basin. The Dabie Mts. is famous for the ultra-high pressure metamorphism. The crustal materials subducted down to the depth of at least 100 km and exhumed. So that the front subjects arise such as the deeply subduction of continent, and the post-collisional crust-mantle interaction. In a geological journey at June of 1999, the author found the rarely variolitic basaltic andesite in the Dabie Mts. It occurs in Susong Group, near Zhifenghe Countryside, Susong County, Anhui Province. It is just to the south of the boundary between the high-grade Susong melange and the ultra-high grade South Dabie melange. It has a noticeable knobby or pitted appearance in the surface. The size of the varioles is about 1-4 mm. In hand-specimen and under microscope, there are distinct contacts between the varioles and the matrice. The mineralogy of the varioles is primarily radiate plagioclase, with little pyroxene, hornblende and quartz. The pyroxene, hornblende and quartz are in the interstices between plagioclase. The matrix is consisted of glass, and micro-crystals of chlorite, epidote and zoisite. It is clearly subjected and extensive alteration. The andesite has an uncommon chemical composition. The SiO_2 content is about 56.8%, TiO_2 = 0.9%, MgO = 6.4%, (Fe_2O_3)_(Total) = 6.7% ～ 7.6%, 100 Mg/(Mg+Fe) = 64.1 ～ 66.2. Mg# is significantly high. The andesite has higher abundances of large-lithophile trace elements (e.g. K, Ba, Sr, LREE), e.g. La/Nd = 5.56-6.07, low abundances of high-strength-field elements (HFSE, e.g. Ta, Nb, P, Ti), particularly Ta and Nb strongly depleted. These are consistent with the characteristics of subducted-related magmas (Pearcce, 1982; Sun and McDonaugh, 1989). In the spider diagram of trace elements, from Ce to right hand, the abundances of elements decrease quickly, showing a characteristic of the continental margins (Pearce, 1982). There has a strongly enrichment of light-rare-earth elements, with a significant diffraction of REEs (the mean value of (La/Yb)_N is 32.84). No Eu anomaly, but there are anomaly high (La/Yb)_N = 28.63-36.74, (La/Y)_N = 70.33 - 82.84. The elements Y and Yb depleted greatly, Y < 20 ppm, Y_N = 2.74-2.84, Yb_N = 2.18 - 2.35. From the La-(La/Sm) diagram, the andesite is derived from partial melting. But the epsilone value of Nd is -18.7 ～ -19.2, so that the material source may be the mantle materials affected by the crustal materials. The Nd model age is 1.9 Ga indicating that the basaltic andesite was resulted from the post-collisional crust-mantle interaction between the subducted Yangze carton and the mantle of Sino-Korea carton. To obtain the Te of the lithosphere beneath the Dabie Mts. and the Hefei Basin, the author applies the coherence method in this thesis. The author makes two topography-gravity profiles (profiles 7~(th) and 9~(th)) across the Dabie Mts. and the Hefei Basin, and calculates the auto-coherence, across coherence, power spectrum, across power spectrum of the topography and gravity of the two profiles. From the relationships between the coherence and the wave-number of profiles. From the relationships between the coherence and the wave-number of profiles 7~(th) and 9~(th), it is obtained that the characteristic wavelengths respectively are 157 km and 126 km. Consequently the values of effective elastic thickness are 6.5 km and 4.8 km, respectively. However, the Te values merely are the minimum value of the lithosphere because the coherencemethod in a relative small region will generate a systemic underestimation. Why there is a so low Te value? In order to check the strength of the lithosphere beneath the Dabie Mts., the authore tries to outline the yielding-stress envelope of the lithosphere. It is suggested that the elastic layers in the crust and upper mantle are 18 km and 35 km, respectively. Since there exist a low viscosity layer about 3-5 km thickness, so it is reasonable that the decoupling between the crust and mantle occurred. So the effective thickness of the lithosphere can be estimated from the two elastic layers. Te is about 34 km. This is the maximum strength of the lithosphere. We can make an approximately estimation about the strength of the lithosphere beneath the Dabie Mts.: Te is about 20-30 km. The author believes that the following factors should be responsible for the low Te value: (1) the Dabie Mts. has elevated strongly since K_3-J_1. The north part of the Dabie Mts. elevates faster than the south part today; (2) there occur large active striking faults in this area. And in the east, the huge Tan-Lu striking fault anyway tends to decrease the lithosphere strength; (3) the lithosphere beneath the Dabie Mts. is heter-homogeneous in spatio-temporal; (4) the study area just locates in the adjacent region between the eastern China where the lithosphere thickness is significantly reduced and the normal western China. These factors will decrease the lithosphere strength.

地震波场模拟及叠前深度偏移的李群方法

The processes of seismic wave propagation in phase space and one way wave extrapolation in frequency-space domain, if without dissipation, are essentially transformation under the action of one parameter Lie groups. Consequently, the numerical calculation methods of the propagation ought to be Lie group transformation too, which is known as Lie group method. After a fruitful study on the fast methods in matrix inversion, some of the Lie group methods in seismic numerical modeling and depth migration are presented here. Firstly the Lie group description and method of seismic wave propagation in phase space is proposed, which is, in other words, symplectic group description and method for seismic wave propagation, since symplectic group is a Lie subgroup and symplectic method is a special Lie group method. Under the frame of Hamiltonian, the propagation of seismic wave is a symplectic group transformation with one parameter and consequently, the numerical calculation methods of the propagation ought to be symplectic method. After discrete the wave field in time and phase space, many explicit, implicit and leap-frog symplectic schemes are deduced for numerical modeling. Compared to symplectic schemes, Finite difference (FD) method is an approximate of symplectic method. Consequently, explicit, implicit and leap-frog symplectic schemes and FD method are applied in the same conditions to get a wave field in constant velocity model, a synthetic model and Marmousi model. The result illustrates the potential power of the symplectic methods. As an application, symplectic method is employed to give synthetic seismic record of Qinghai foothills model. Another application is the development of Ray+symplectic reverse-time migration method. To make a reasonable balance between the computational efficiency and accuracy, we combine the multi-valued wave field & Green function algorithm with symplectic reverse time migration and thus develop a new ray+wave equation prestack depth migration method. Marmousi model data and Qinghai foothills model data are processed here. The result shows that our method is a better alternative to ray migration for complex structure imaging. Similarly, the extrapolation of one way wave in frequency-space domain is a Lie group transformation with one parameter Z and consequently, the numerical calculation methods of the extrapolation ought to be Lie group methods. After discrete the wave field in depth and space, the Lie group transformation has the form of matrix exponential and each approximation of it gives a Lie group algorithm. Though Pade symmetrical series approximation of matrix exponential gives a extrapolation method which is traditionally regarded as implicit FD migration, it benefits the theoretic and applying study of seismic imaging for it represent the depth extrapolation and migration method in a entirely different way. While, the technique of coordinates of second kind for the approximation of the matrix exponential begins a new way to develop migration operator. The inversion of matrix plays a vital role in the numerical migration method given by Pade symmetrical series approximation. The matrix has a Toepelitz structure with a helical boundary condition and is easy to inverse with LU decomposition. A efficient LU decomposition method is spectral factorization. That is, after the minimum phase correlative function of each array of matrix had be given by a spectral factorization method, all of the functions are arranged in a position according to its former location to get a lower triangular matrix. The major merit of LU decomposition with spectral factorization (SF Decomposition) is its efficiency in dealing with a large number of matrixes. After the setup of a table of the spectral factorization results of each array of matrix, the SF decomposition can give the lower triangular matrix by reading the table. However, the relationship among arrays is ignored in this method, which brings errors in decomposition method. Especially for numerical calculation in complex model, the errors is fatal. Direct elimination method can give the exact LU decomposition But even it is simplified in our case, the large number of decomposition cost unendurable computer time. A hybrid method is proposed here, which combines spectral factorization with direct elimination. Its decomposition errors is 10 times little than that of spectral factorization, and its decomposition speed is quite faster than that of direct elimination, especially in dealing with a large number of matrix. With the hybrid method, the 3D implicit migration can be expected to apply on real seismic data. Finally, the impulse response of 3D implicit migration operator is presented.

用宽频带远震体波波形研究中国及邻近地区地壳构造及地质意义初探

The theory and approach of the broadband teleseismic body waveform inversion are expatiated in this paper, and the defining the crust structure's methods are developed. Based on the teleseismic P-wave data, the theoretic image of the P-wave radical component is calculated via the convolution of the teleseismic P-wave vertical component and the transform function, and thereby a P-wavefrom inversion method is built. The applied results show the approach effective, stable and its resolution high. The exact and reliable teleseismic P waveforms recorded by CDSN and IRIS and its geodynamics are utilized to obtain China and its vicinage lithospheric transfer functions, this region ithospheric structure is inverted through the inversion of reliable transfer functions, the new knowledge about the deep structure of China and its vicinage is obtained, and the reliable seismological evidence is provided to reveal the geodynamic evolution processes and set up the continental collisional theory. The major studies are as follows: Two important methods to study crustal and upper mantle structure -- body wave travel-time inversion and waveform modeling are reviewed systematically. Based on ray theory, travel-time inversion is characterized by simplicity, crustal and upper mantle velocity model can be obtained by using 1-D travel-time inversion preliminary, which introduces the reference model for studying focal location, focal mechanism, and fine structure of crustal and upper mantle. The large-scale lateral inhomogeneity of crustal and upper mantle can be obtained by three-dimensional t ravel-time seismic tomography. Based on elastic dynamics, through the fitting between theoretical seismogram and observed seismogram, waveform modeling can interpret the detail waveform and further uncover one-dimensional fine structure and lateral variation of crustal and upper mantle, especially the media characteristics of singular zones of ray. Whatever travel-time inversion and waveform modeling is supposed under certain approximate conditions, with respective advantages and disadvantages, and provide convincing structure information for elucidating physical and chemical features and geodynamic processes of crustal and upper mantle. Because the direct wave, surface wave, and refraction wave have lower resolution in investigating seismic velocity transitional zone, which is inadequate to study seismic discontinuities. On the contrary, both the converse and reflected wave, which sample the discontinuities directly, must be carefully picked up from seismogram to constrain the velocity transitional zones. Not only can the converse wave and reflected wave study the crustal structure, but also investigate the upper mantle discontinuities. There are a number of global and regional seismic discontinuities in the crustal and upper mantle, which plays a significant role in understanding physical and chemical properties and geodynamic processes of crustal and upper mantle. The broadband teleseismic P waveform inversion is studied particularly. The teleseismic P waveforms contain a lot of information related to source time function, near-source structure, propagation effect through the mantle, receiver structure, and instrument response, receiver function is isolated form teleseismic P waveform through the vector rotation of horizontal components into ray direction and the deconvolution of vertical component from the radial and tangential components of ground motion, the resulting time series is dominated by local receiver structure effect, and is hardly irrelevant to source and deep mantle effects. Receiver function is horizontal response, which eliminate multiple P wave reflection and retain direct wave and P-S converted waves, and is sensitive to the vertical variation of S wave velocity. Velocity structure beneath a seismic station has different response to radial and vertical component of an accident teleseismic P wave. To avoid the limits caused by a simplified assumption on the vertical response, the receiver function method is mended. In the frequency domain, the transfer function is showed by the ratio of radical response and vertical response of the media to P wave. In the time domain, the radial synthetic waveform can be obtained by the convolution of the transfer function with the vertical wave. In order to overcome the numerical instability, generalized reflection and transmission coefficient matrix method is applied to calculate the synthetic waveform so that all multi-reflection and phase conversion response can be included. A new inversion method, VFSA-LM method, is used in this study, which successfully combines very fast simulated annealing method (VFSA) with damped least square inversion method (LM). Synthetic waveform inversion test confirms its effectiveness and efficiency. Broadband teleseismic P waveform inversion is applied in lithospheric velocity study of China and its vicinage. According to the data of high quality CDSN and IRIS, we obtained an outline map showing the distribution of Asian continental crustal thickness. Based on these results gained, the features of distribution of the crustal thickness and outline of crustal structure under the Asian continent have been analyzed and studied. Finally, this paper advances the principal characteristics of the Asian continental crust. There exist four vast areas of relatively minor variations in the crustal thickness, namely, northern, eastern southern and central areas of Asian crust. As a byproduct, the earthquake location is discussed, Which is a basic issue in seismology. Because of the strong trade-off between the assumed initial time and focal depth and the nonlinear of the inversion problems, this issue is not settled at all. Aimed at the problem, a new earthquake location method named SAMS method is presented, In which, the objective function is the absolute value of the remnants of travel times together with the arrival times and use the Fast Simulated Annealing method is used to inverse. Applied in the Chi-Chi event relocation of Taiwan occurred on Sep 21, 2000, the results show that the SAMS method not only can reduce the effects of the trade-off between the initial time and focal depth, but can get better stability and resolving power. At the end of the paper, the inverse Q filtering method for compensating attenuation and frequency dispersion used in the seismic section of depth domain is discussed. According to the forward and inverse results of synthesized seismic records, our Q filtrating operator of the depth domain is consistent with the seismic laws in the absorbing media, which not only considers the effect of the media absorbing of the waves, but also fits the deformation laws, namely the frequency dispersion of the body wave. Two post stacked profiles about 60KM, a neritic area of China processed, the result shows that after the forward Q filtering of the depth domain, the wide of the wavelet of the middle and deep layers is compressed, the resolution and signal noise ratio are enhanced, and the primary sharp and energy distribution of the profile are retained.

Tensor density matrix theory for determining population and alignment of symmetric top molecule using rotationally unresolved fluorescence

General expressions used for transforming raw laser-induced fluorescence (LIF) intensity into the population and alignment parameters of a symmetric top molecule are derived by employing the density matrix approach. The molecular population and alignment are described by molecular state multipoles. The results are presented for a general excitation-detection geometry and then applied to some special geometries. In general cases, the LIF intensity is a complex function of the initial molecular state multipoles, the dynamic factors and the excitation-detection geometrical factors. It contains a population and 14 alignment multipoles. How to extract all initial state multipoles from the rotationally unresolved emission LIF intensity is discussed in detail.

Separation and detection of compounds in Honeysuckle by integration of ion-exchange chromatography fractionation with reversed-phase liquid chromatography-atmospheric pressure chemical ionization mass spectrometer and matrix-assisted laser desorption/ionization time-of-flight mass spectrometry analysis

A hyphenated method for the isolation and identification of components in a traditional Chinese medicine of Honeysuckle was developed. Ion-exchange chromatography (IEC) was chosen for the fractionation of Honeysuckle extract, and then followed by concentration of all the fractions with rotary vacuum evaporator. Each of the enriched fractions was then further analyzed by reversed-phase liquid chromatography-atmospheric pressure chemical ionization mass spectrometer (RPLC-APCI/MS) and matrix-assisted laser desorption/ionization time-of-flight mass spectrometry (MALDI-TOF/MS) with matrix of oxidized carbon nanotubes, respectively. It can be noted totally more than 117 components were detected by UV detector, APCI/MS and MALDI-TOF/MS in Honeysuckle extract except the, 145 components identified by MALDI-TOF/MS alone with this integrated approach, and 7 of them were preliminary identified according to their UV spectra and mass spectra performed by APCI/MS and MALDI-TOF/MS, respectively. The obtained analytical results not only indicated the approach of integration IEC fractionation with RPLC-APCI/MS and MALDI-TOF/MS is capable of analyzing complex samples, but also exhibited the potential power of the mass spectrometer in detection of low-mass compounds, such as traditional Chinese medicines (TCMs) and complex biological samples. (c) 2005 Elsevier B.V. All rights reserved.

The effects of ageing and exercise on skeletal muscle structure and function

Musculoskeletal ageing is associated with profound morphological and functional changes that increase fall risk and disease incidence and is characterised by age-related reductions in motor unit number and atrophy of muscle fibres, particularly type II fibres. Decrements in functional strength and power are relatively modest until the 6th decade, after which the rate of loss exponentially accelerates, particularly beyond the 8th decade of life. Physical activity is a therapeutic modality that can significantly attenuate age-related decline. The underlying signature of ageing, as manifested by perturbed redox homeostasis, leads to a blunting of acute and chronic redox regulated exercise adaptations. Impaired redox regulated exercise adaptations are mechanistically related to altered exercise-induced reactive oxygen and nitrogen species generation and a resultant failure to properly activate redox regulated signaling cascades. Despite the aforementioned specific impairment in redox signaling, exercise induces a plethora of beneficial effects, irrespective of age. There is, therefore, strong evidence for promoting regular physical exercise, especially progressive resistance training as a lifelong habitual practice.

TCP Optimization through FEC, ARQ and Transmission Power Tradeoffs

TCP performance degrades when end-to-end connections extend over wireless connections-links which are characterized by high bit error rate and intermittent connectivity. Such link characteristics can significantly degrade TCP performance as the TCP sender assumes wireless losses to be congestion losses resulting in unnecessary congestion control actions. Link errors can be reduced by increasing transmission power, code redundancy (FEC) or number of retransmissions (ARQ). But increasing power costs resources, increasing code redundancy reduces available channel bandwidth and increasing persistency increases end-to-end delay. The paper proposes a TCP optimization through proper tuning of power management, FEC and ARQ in wireless environments (WLAN and WWAN). In particular, we conduct analytical and numerical analysis taking into "wireless-aware" TCP) performance under different settings. Our results show that increasing power, redundancy and/or retransmission levels always improves TCP performance by reducing link-layer losses. However, such improvements are often associated with cost and arbitrary improvement cannot be realized without paying a lot in return. It is therefore important to consider some kind of net utility function that should be optimized, thus maximizing throughput at the least possible cost.

On the Origin of Power Laws in Internet Topologies

Recent empirical studies have shown that Internet topologies exhibit power laws of the form for the following relationships: (P1) outdegree of node (domain or router) versus rank; (P2) number of nodes versus outdegree; (P3) number of node pairs y = x^α within a neighborhood versus neighborhood size (in hops); and (P4) eigenvalues of the adjacency matrix versus rank. However, causes for the appearance of such power laws have not been convincingly given. In this paper, we examine four factors in the formation of Internet topologies. These factors are (F1) preferential connectivity of a new node to existing nodes; (F2) incremental growth of the network; (F3) distribution of nodes in space; and (F4) locality of edge connections. In synthetically generated network topologies, we study the relevance of each factor in causing the aforementioned power laws as well as other properties, namely diameter, average path length and clustering coefficient. Different kinds of network topologies are generated: (T1) topologies generated using our parametrized generator, we call BRITE; (T2) random topologies generated using the well-known Waxman model; (T3) Transit-Stub topologies generated using GT-ITM tool; and (T4) regular grid topologies. We observe that some generated topologies may not obey power laws P1 and P2. Thus, the existence of these power laws can be used to validate the accuracy of a given tool in generating representative Internet topologies. Power laws P3 and P4 were observed in nearly all considered topologies, but different topologies showed different values of the power exponent α. Thus, while the presence of power laws P3 and P4 do not give strong evidence for the representativeness of a generated topology, the value of α in P3 and P4 can be used as a litmus test for the representativeness of a generated topology. We also find that factors F1 and F2 are the key contributors in our study which provide the resemblance of our generated topologies to that of the Internet.

Systematic delay-driven power optimisation and power-driven delay optimisation of combinational circuits

With the proliferation of mobile wireless communication and embedded systems, the energy efficiency becomes a major design constraint. The dissipated energy is often referred as the product of power dissipation and the input-output delay. Most of electronic design automation techniques focus on optimising only one of these parameters either power or delay. Industry standard design flows integrate systematic methods of optimising either area or timing while for power consumption optimisation one often employs heuristics which are characteristic to a specific design. In this work we answer three questions in our quest to provide a systematic approach to joint power and delay Optimisation. The first question of our research is: How to build a design flow which incorporates academic and industry standard design flows for power optimisation? To address this question, we use a reference design flow provided by Synopsys and integrate in this flow academic tools and methodologies. The proposed design flow is used as a platform for analysing some novel algorithms and methodologies for optimisation in the context of digital circuits. The second question we answer is: Is possible to apply a systematic approach for power optimisation in the context of combinational digital circuits? The starting point is a selection of a suitable data structure which can easily incorporate information about delay, power, area and which then allows optimisation algorithms to be applied. In particular we address the implications of a systematic power optimisation methodologies and the potential degradation of other (often conflicting) parameters such as area or the delay of implementation. Finally, the third question which this thesis attempts to answer is: Is there a systematic approach for multi-objective optimisation of delay and power? A delay-driven power and power-driven delay optimisation is proposed in order to have balanced delay and power values. This implies that each power optimisation step is not only constrained by the decrease in power but also the increase in delay. Similarly, each delay optimisation step is not only governed with the decrease in delay but also the increase in power. The goal is to obtain multi-objective optimisation of digital circuits where the two conflicting objectives are power and delay. The logic synthesis and optimisation methodology is based on AND-Inverter Graphs (AIGs) which represent the functionality of the circuit. The switching activities and arrival times of circuit nodes are annotated onto an AND-Inverter Graph under the zero and a non-zero-delay model. We introduce then several reordering rules which are applied on the AIG nodes to minimise switching power or longest path delay of the circuit at the pre-technology mapping level. The academic Electronic Design Automation (EDA) tool ABC is used for the manipulation of AND-Inverter Graphs. We have implemented various combinatorial optimisation algorithms often used in Electronic Design Automation such as Simulated Annealing and Uniform Cost Search Algorithm. Simulated Annealing (SMA) is a probabilistic meta heuristic for the global optimization problem of locating a good approximation to the global optimum of a given function in a large search space. We used SMA to probabilistically decide between moving from one optimised solution to another such that the dynamic power is optimised under given delay constraints and the delay is optimised under given power constraints. A good approximation to the global optimum solution of energy constraint is obtained. Uniform Cost Search (UCS) is a tree search algorithm used for traversing or searching a weighted tree, tree structure, or graph. We have used Uniform Cost Search Algorithm to search within the AIG network, a specific AIG node order for the reordering rules application. After the reordering rules application, the AIG network is mapped to an AIG netlist using specific library cells. Our approach combines network re-structuring, AIG nodes reordering, dynamic power and longest path delay estimation and optimisation and finally technology mapping to an AIG netlist. A set of MCNC Benchmark circuits and large combinational circuits up to 100,000 gates have been used to validate our methodology. Comparisons for power and delay optimisation are made with the best synthesis scripts used in ABC. Reduction of 23% in power and 15% in delay with minimal overhead is achieved, compared to the best known ABC results. Also, our approach is also implemented on a number of processors with combinational and sequential components and significant savings are achieved.

Power system parameters forecasting using Hilbert-Huang transform and machine learning

A novel hybrid data-driven approach is developed for forecasting power system parameters with the goal of increasing the efficiency of short-term forecasting studies for non-stationary time-series. The proposed approach is based on mode decomposition and a feature analysis of initial retrospective data using the Hilbert-Huang transform and machine learning algorithms. The random forests and gradient boosting trees learning techniques were examined. The decision tree techniques were used to rank the importance of variables employed in the forecasting models. The Mean Decrease Gini index is employed as an impurity function. The resulting hybrid forecasting models employ the radial basis function neural network and support vector regression. A part from introduction and references the paper is organized as follows. The second section presents the background and the review of several approaches for short-term forecasting of power system parameters. In the third section a hybrid machine learningbased algorithm using Hilbert-Huang transform is developed for short-term forecasting of power system parameters. Fourth section describes the decision tree learning algorithms used for the issue of variables importance. Finally in section six the experimental results in the following electric power problems are presented: active power flow forecasting, electricity price forecasting and for the wind speed and direction forecasting.