Nonlinear analysis of cardiac rhythm fluctuations using DFA method

An approach by which the detrented fluctuation analysis (DFA) method can be used to help diagnose heart failure was demonstrated. DFA was applied to patients suffering from congestive heart failure (CHF) to check correlations between DFA indices and CHF, and determine a correlation between DFA indices and mortality, with a particular attention to the residue parameter, which is a measure of the departure of the DFA from its power law approximation. DFA parameters proved to be useful as a complement to the physiological parameters weber and FE to sort out the patients into three prognostic group.

The stochastic implicit Euler method - A stable coupling scheme for Monte Carlo burnup calculations

Existing Monte Carlo burnup codes use various schemes to solve the coupled criticality and burnup equations. Previous studies have shown that the coupling schemes of the existing Monte Carlo burnup codes can be numerically unstable. Here we develop the Stochastic Implicit Euler method - a stable and efficient new coupling scheme. The implicit solution is obtained by the stochastic approximation at each time step. Our test calculations demonstrate that the Stochastic Implicit Euler method can provide an accurate solution to problems where the methods in the existing Monte Carlo burnup codes fail. © 2013 Elsevier Ltd. All rights reserved.

IGBT converters conducted EMI analysis by controlled multiple-slope switching waveform approximation

IGBTs realise high-performance power converters. Unfortunately, with fast switching of the IGBT-free wheel diode chopper cell, such circuits are intrinsic sources of high-level EMI. Therefore, costly EMI filters or shielding are normally needed on the load and supply side. In order to design these EMI suppression components, designers need to predict the EMI level with reasonable accuracy for a given structure and operating mode. Simplifying the transient IGBT switching current and voltage into a multiple slope switching waveform approximation offers a feasible way to estimate conducted EMI with some accuracy. This method is dependent on the availability of high-fidelity measurements. Also, that multiple slope approximation needs careful and time-costly IGBT parameters optimisation process to approach the real switching waveform. In this paper, Active Voltage Control Gate Drive(AVC GD) is employed to shape IGBT switching into several defined slopes. As a result, Conducted EMI prediction by multiple slope switching approximation could be more accurate, less costly but more friendly for implementation. © 2013 IEEE.

A coupled hydro-mechanical analysis for prediction of hydraulic fracture propagation in saturated porous media using EFG mesh-less method

The details of the Element Free Galerkin (EFG) method are presented with the method being applied to a study on hydraulic fracturing initiation and propagation process in a saturated porous medium using coupled hydro-mechanical numerical modelling. In this EFG method, interpolation (approximation) is based on nodes without using elements and hence an arbitrary discrete fracture path can be modelled.The numerical approach is based upon solving two governing partial differential equations of equilibrium and continuity of pore water simultaneously. Displacement increment and pore water pressure increment are discretized using the same EFG shape functions. An incremental constrained Galerkin weak form is used to create the discrete system of equations and a fully implicit scheme is used for discretization in the time domain. Implementation of essential boundary conditions is based on the penalty method. In order to model discrete fractures, the so-called diffraction method is used.Examples are presented and the results are compared to some closed-form solutions and FEM approximations in order to demonstrate the validity of the developed model and its capabilities. The model is able to take the anisotropy and inhomogeneity of the material into account. The applicability of the model is examined by simulating hydraulic fracture initiation and propagation process from a borehole by injection of fluid. The maximum tensile strength criterion and Mohr-Coulomb shear criterion are used for modelling tensile and shear fracture, respectively. The model successfully simulates the leak-off of fluid from the fracture into the surrounding material. The results indicate the importance of pore fluid pressure in the initiation and propagation pattern of fracture in saturated soils. © 2013 Elsevier Ltd.

SLP: A zero-contact non-invasive method for pulmonary function testing

Structured Light Plethysmography (SLP) is a novel non-invasive method that uses structured light to perform pulmonary function testing that does not require physical contact with a patient. The technique produces an estimate of chest wall volume changes over time. A patient is observed continuously by two cameras and a known pattern of light (i.e. structured light) is projected onto the chest using an off-the-shelf projector. Corner features from the projected light pattern are extracted, tracked and brought into correspondence for both camera views over successive frames. A novel self calibration algorithm recovers the intrinsic and extrinsic camera parameters from these point correspondences. This information is used to reconstruct a surface approximation of the chest wall and several novel ideas for 'cleaning up' the reconstruction are used. The resulting volume and derived statistics (e.g. FVC, FEV) agree very well with data taken with a spirometer. © 2010. The copyright of this document resides with its authors.

Calculation of light delay for coupled microrings by FDTD technique and Pade approximation

The Pade approximation with Baker's algorithm is compared with the least-squares Prony method and the generalized pencil-of-functions (GPOF) method for calculating mode frequencies and mode Q factors for coupled optical microdisks by FDTD technique. Comparisons of intensity spectra and the corresponding mode frequencies and Q factors show that the Pade approximation can yield more stable results than the Prony and the GPOF methods, especially the intensity spectrum. The results of the Prony method and the GPOF method are greatly influenced by the selected number of resonant modes, which need to be optimized during the data processing, in addition to the length of the time response signal. Furthermore, the Pade approximation is applied to calculate light delay for embedded microring resonators from complex transmission spectra obtained by the Pade approximation from a FDTD output. The Prony and the GPOF methods cannot be applied to calculate the transmission spectra, because the transmission signal obtained by the FDTD simulation cannot be expressed as a sum of damped complex exponentials. (C) 2009 Optical Society of America

Computation of resonant frequencies and quality factors of cavities by FDTD technique and Pade approximation

The finite-difference time domain (FDTD) technique and the Pade approximation with Baker's algorithm are used to calculate the mode frequencies and quality factors of cavities. Comparing with the fast Fourier transformation/Pade method, we find that the Fade approximation and the Baker's algorithm can obtain exact resonant frequencies and quality factors based on a much shorter time record of the FDTD output.

CALCULATION OF DIELECTRIC-CONSTANTS FOR SEMICONDUCTORS - AN APPLICATION OF ABINITIO PSEUDOPOTENTIAL METHOD

We have used ab initio pseudopotential method to generate basis wavefunctions and eigen energies to carry out first principle calculations of the static macroscopic dielectric constant for GaAs and GaP. The resulted converged random phase approximation (RPA) value is 12.55 and 10.71, in excellent agreement to the experimental value of 12.4 and 10.86, respectively. The inclusion of the exchange correlation contribution makes the calculated result slightly worsen. A convergence test with respect to the number of k points in Brillouin zone (BZ) integration was carried out. Sixty irreducible BZ k points were used to achieve the converged results. Integration with only 10 special k points increased the RPA value by 15%.

Applicatioo of Pade Approximation in Simulating Photonic Crystals

To save finite-difference time-domain(FDTD) computing time, several methods are proposed to convert the time domain FDTD output into frequency domain. The Padé approximation with Baker's algorithm and the program are introduced to simulate photonic crystal structures. For a simple pole system with frequency 160THz and quality factor of 5000,the intensity spectrum obtained by the Padé approximation from a 28-item sequence output is more exact than that obtained by fast Fourier transformation from a 220-item sequence output. The mode frequencies and quality factors are calculated at different wave vectors for the photonic crystal slab from a much shorter FDTD output than that required by the FFT method,and then the band diagrams are obatined. In addition,mode frequencies and Q-factors are calculated for photonic crystal microcavity.

Study on tapered crossed subwavelength gratings by Fourier modal method

Fourier modal method incorporating staircase approximation is used to study tapered crossed subwavelength gratings in this paper. Three intuitive formulations of eigenvalue functions originating from the prototype are presented, and their convergences are compared through numerical calculation. One of them is found to be suitable in modeling the diffraction efficiency of the circular tapered crossed subwavelength gratings without high absorption, and staircase approximation is further proven valid for non-highly-absorption tapered gratings. This approach is used to simulate the "moth-eye" antireflection surface on silicon, and the numerical result agrees well with the experimental one.

An eigenfunction expansion-variational method in prediction of the transverse thermal conductivity of fiber reinforced composites considering interfacial characteristics

An eigenfunction expansion-variational method based on a unit cell is developed to deal with the steady-state heat conduction problem of doubly-periodic fiber reinforced composites with interfacial thermal contact resistance or coating. The numerical results show a rapid convergence of the present method. The present solution provides a unified first-order approximation formula of the effective thermal conductivity for different interfacial characteristics and fiber distributions. A comparison with the present high-order results, available experimental data and micromechanical estimations demonstrates that the first-order approximation formula is a good engineering closed-form formula. An engineering equivalent parameter reflecting the overall influence of the thermal conductivities of the matrix and fibers and the interfacial characteristic on the effective thermal conductivity, is found. The equivalent parameter can greatly simplify the complicated relation of the effective thermal conductivity to the internal structure of a composite. (c) 2010 Elsevier Ltd. All rights reserved.

Analytical method of predicating the instabilities of a micro arch-shaped beam under electrostatic loading

An arch-shaped beam with different configurations under electrostatic loading experiences either the direct pull-in instability or the snap-through first and then the pull-in instability. When the pull-in instability occurs, the system collides with the electrode and adheres to it, which usually causes the system failure. When the snap-through instability occurs, the system experiences a discontinuous displacement to flip over without colliding with the electrode. The snap-through instability is an ideal actuation mechanism because of the following reasons: (1) after snap-through the system regains the stability and capability of withstanding further loading; (2) the system flips back when the loading is reduced, i.e. the system can be used repetitively; and (3) when approaching snap-through instability the system effective stiffness reduces toward zero, which leads to a fast flipping-over response. To differentiate these two types of instability responses for an arch-shaped beam is vital for the actuator design. For an arch-shaped beam under electrostatic loading, the nonlinear terms of the mid-plane stretching and the electrostatic loading make the analytical solution extremely difficult if not impossible and the related numerical solution is rather complex. Using the one mode expansion approximation and the truncation of the higher-order terms of the Taylor series, we present an analytical solution here. However, the one mode approximation and the truncation error of the Taylor series can cause serious error in the solution. Therefore, an error-compensating mechanism is also proposed. The analytical results are compared with both the experimental data and the numerical multi-mode analysis. The analytical method presented here offers a simple yet efficient solution approach by retaining good accuracy to analyze the instability of an arch-shaped beam under electrostatic loading.

Isoscalar Giant Monopole Resonance in Relativistic Continuum Random Phase Approximation

The isoscalar giant monopole resonance (ISGMR) in nuclei is studied in the framework of a fully consistent relativistic continuum random phase approximation (RCRPA). In this method the contribution of the continuum spectrum to nuclear excitations is treated exactly by the single particle Green's function technique. The negative energy states in the Dirac sea are also included in the single particle Green's function in the no-sea approximation. The single particle Green's function is calculated numerically by a proper product of the regular and irregular solutions of the Dirac equation. The strength distributions in the RCRPA calculations, the inverse energy-weighted sum rule m(-1) and the centroid energy of the ISGMR in Sn-120 and Pb-208 are analysed. Numerical results of the RCRPA are checked with the constrained relativistic mean field model and relativistic random phase approximation with a discretized spectrum in the continuum. Good agreement between them is achieved.

Relativistic Continuum Random Phase Approximation and Applications II. Applications

The fully consistent relativistic continuum random phase approximation (RCRPA) has been constructed in the momentum representation in the first part of this paper. In this part we describe the numerical details for solving the Bethe-Salpeter equation. The numerical results are checked by the inverse energy weighted sum rules in the isoscalar giant monopole resonance, which are obtained from the constraint relativistic mean field theory and also calculated with the integration of the RCRPA strengths. Good agreement between the misachieved. We study the effects of the self-consistency violation, particularly the currents and Coulomb interaction to various collective multipole excitations. Using the fully consistent RCRPA method, we investigate the properties of isoscalar and isovector collective multipole excitations for some stable and exotic from light to heavy nuclei. The properties of the resonances, such as the centroid energies and strength distributions are compared with the experimental data as well as with results calculated in other models.

Relativistic Continuum Random Phase Approximation and Applications I. Formalism

A fully consistent relativistic continuum random phase approximation (RCRPA) is constructed in terms of the Green's function technique. In this method the contribution of the continuum spectrum to nuclear excitations is treated exactly by the single particle Green's function, which includes also the negative states in the Dirac sea in the nose aapproximation. The theoretical formalism of RCRPA and numerical details are presented. The single particle Green's function is calculated numerically by a proper product of regular and irregular solutions of the Dirac equation. The numerical details and the formalism of RCRPA in the momentum representation are presented.