Substitution-newton-Raphson method for the solution of electric network equations

The power flow problem, in transmission networks, has been well solved, for most cases, using Newton-Raphson method (NR) and its decoupled versions. Generally speaking, the solution of a non-linear system of equations refers to two methods: NR and Successive Substitution. The proposal of this paper is to evaluate the potential of the Substitution-Newton-Raphson Method (SNR), which combines both methods, on the solution of the power flow problem. Simulations were performed using a two-bus test network in order to observe the characteristics of these methods. It was verified that the NR is faster than SNR, in terms of convergence, considering non-stressed scenarios. For those cases where the power flow in the network is closed to the limits (stressed system), the SNR converges faster. This paper presents the power flow formulation of the SNR and describes its potential for its application in special cases such as stressed scenarios. © 2006 IEEE.

Assessing the effectiveness of 30% sodium chloride aqueous solution for the preservation of fixed anatomical specimens: a 5-year follow-up study

Anatomical specimens used in human or veterinary anatomy laboratories are usually prepared with formaldehyde (a cancerous and teratogenic substance), glycerin (an expensive and viscous fluid), or ethanol (which is flammable). This research aimed to verify the viability of an aqueous 30% sodium chloride solution for preservation of anatomical specimens previously fixed with formaldehyde. Anatomical specimens of ruminant, carnivorous, equine, swine and birds were used. All were previously fixed with an aqueous 20% formaldehyde solution and held for 7days in a 10% aqueous solution of the same active ingredient. During the first phase of the experiment, small specimens of animal tissue previously fixed in formaldehyde were distributed in vials with different concentrations of formaldehyde, with or without 30% sodium chloride solution, a group containing only 30% sodium chloride, and a control group containing only water. During this phase, no contamination was observed in any specimen containing 30% sodium chloride solution, whether alone or in combination with different concentrations of formaldehyde. In the second phase of the experiment, the 30% sodium chloride solution, found to be optimal in the first phase of the experiment, was tested for its long-term preservation properties. For a period of 5years, the preserved specimens were evaluated three times a week for visual contamination, odors, and changes in color and texture. There was no visual contamination or decay found in any specimen. Furthermore, no strange odors, or changes in color or softness were noted. The 30% sodium chloride solution was determined to be effective in the preservation of anatomic specimens previously fixed in formaldehyde.

Study Of The Binding Of Substrate By The Phe557Val Mutant Soybean Lipoxygenase-1

Lipoxygenases are a class of enzymes which consist of non-heme iron dioxygenases that are produced by fungi, plants, and mammals and catalyze the oxygenation of polyunsaturated fatty acid substrates to unsaturated fatty acid hydroperoxide products. The unsaturated fatty acid hydroperoxide products are stereo- and regiospecific. One such lipoxygenase, soybean lipoxygenase-1 (SBLO-1), catalyzes the conversion of linoleate to 13-hydroperoxy-9(Z),11(E)-octadecadienoate (13-HPOD) and a small amount of 9-hydroperoxy-10(E),12(Z)-octadecadienoate (9-HPOD). Although the structure of SBLO-1 is known and it is the most widely studied lipoxygenase, how it binds to substrate is still poorly understood. Two competing binding hypotheses that have been used to understand and explain the binding are the head first binding model and the tail first binding model. The head first binding model predicts linoleate binds with its polar carboxylate group in the binding pocket and the methyl terminus at the surface of the binding pocket. The tail first binding model predicts that linoleate binds with its methyl terminus end in the binding pocket and the polar carboxylate group at the surface of the binding pocket. Both binding models have been used in the explanation of previous work. In previous work the replacement of phenylalanine with valine has been performed to produce the phe557val mutant SBLO-1. The mutant SBLO-1 was then used in the enzymatic oxygenation of linoleate. With this mutant, the amount of 9-HPOD that is formed increases. This result has been interpreted using the head-first binding model in which the smaller valine residue allows linoleate to bind with the polar carboxylate group of linoleate interacting with arginine-707. The work presented in this thesis confirms the regiochemical results of the previous work and further tests the head-first binding model. If head-first binding occurs, the 9-HPOD is expected to have primarily S configuration. Utilizing chiral-phase HPLC, it was found that the 9-HPOD produced by the phe557val mutant SBLO-1 is primarily S, consistent with head-first binding. The head-first binding model was also tested using linoleyl dimethylamine (LDMA), which has been shown to be a good substrate for SBLO-1 at pH 7.0, where LDMA is thought to be positively charged. This model predicts that less of the 9-peroxide should be produced with this substrate. Through the use of gas chromatography/mass spectrometry, it was found that the conversion of LDMA by the phe557val mutant SBLO-1 resulted in the formation of a 46:54 mixture of the 13-peroxide:9-peroxide. The higher amount of 9-peroxide is the opposite of what is expected for the currently proposed model suggesting that the proposed model may not be entirely correct. The results thus far have been consistent with reverse binding but not with the proposed interaction of the polar end of the substrate with arginine-707.

The Effect of the Cu:Fe''' Ratio Upon the Current Efficiency in the Electrolysis of a Copper Sulfate Solution Containing Iron Sulfate

It is a well-known fact that, in the electrolysis of a CuSO4 solution containing iron sulfate, using insoluble anodes, with the depletion of copper, the point is finally reached where the current efficiency becomes zero. This decrease in current efficiency is due to the oxidation of the ferrous sulfate to the ferric condition at the anode, by the oxygen liberated. The resulting ferric sulfate diffuses over to the cathode and there dissolves copper from the cathode according to the chemical equation Cu + Fe2 (SO4)3 = CuSO4 + 2FeSO4. This copper, which has been deposited at the cathode by the electric current, is thus redissolved by the Fe2(SO4)3. The solution of the copper causes at the same time a formation of FeSO4 which in turn diffuses over to the anode and is there oxidized to Fe2(SO4)3; and so the cycle continues, using electric current without rendering useful work. E. H. Larison has noted that a definite amount of ferric salts must be reduced to the ferrous condition before all the copper will remain on the cathode; he does not state, however, just what this point is. L. Addicks has plotted the relation between current efficiency and ferric sulphate content. The existence of the results scattered the points more or less, although the decrease in current efficiency with increased ferric sulphate content is clearly indicated. E. T.Kern has likewise noted that the smaller the amount of copper in the solution, the greater is the reduction of current efficiency. In this work, therefore, it was desired to determine what amount of ferric iron was permissible in a copper sulfate solution of definite concentration before the current efficiency would drop to zero, and what, if any, was the effect of definite Cu:Fe’’’ratio upon the current efficiency of the electrolysis.

Matrix-free time-stepping methods for the solution of TriGlobal instability problems

La inmensa mayoría de los flujos de relevancia ingenieril permanecen sin estudiar en el marco de la teoría de estabilidad global. Esto es debido a dos razones fundamentalmente, las dificultades asociadas con el análisis de los flujos turbulentos y los inmensos recursos computacionales requeridos para obtener la solución del problema de autovalores asociado al análisis de inestabilidad de flujos tridimensionales, también conocido como problema TriGlobal. En esta tesis se aborda el problema asociado con la tridimensionalidad. Se ha desarrollado una metodología general para obtener soluciones de problemas de análisis modal de las inestabilidades lineales globales mediante el acoplamiento de métodos de evolución temporal, desarrollados en este trabajo, con códigos de mecánica de fluidos computacional de segundo orden, utilizados de forma general en la industria. Esta metodología consiste en la resolución del problema de autovalores asociado al análisis de inestabilidad mediante métodos de proyección en subespacios de Krylov, con la particularidad de que dichos subespacios son generados por medio de la integración temporal de un vector inicial usando cualquier código de mecánica de fluidos computacional. Se han elegido tres problemas desafiantes en función de la exigencia de recursos computacionales necesarios y de la complejidad física para la demostración de la presente metodología: (i) el flujo en el interior de una cavidad tridimensional impulsada por una de sus tapas, (ii) el flujo alrededor de un cilindro equipado con aletas helicoidales a lo largo su envergadura y (iii) el flujo a través de una cavidad abierta tridimensinal en ausencia de homogeneidades espaciales. Para la validación de la tecnología se ha obtenido la solución del problema TriGlobal asociado al flujo en la cavidad tridimensional, utilizando el método de evolución temporal desarrollado acoplado con los operadores numéricos de flujo incompresible del código CFD OpenFOAM (código libre). Los resultados obtenidos coinciden plentamente con la literatura. La aplicación de esta metodología al estudio de inestabilidades globales de flujos abiertos tridimensionales ha proporcionado por primera vez, información sobre la transición tridimensional de estos flujos. Además, la metodología ha sido adaptada para resolver problemas adjuntos TriGlobales, permitiendo el control de flujo basado en modificaciones de las inestabilidades globales. Finalmente, se ha demostrado que la cantidad moderada de los recursos computacionales requeridos para la solución del problema de valor propio TriGlobal usando este método numérico, junto a su versatilidad al poder acoplarse a cualquier código aerodinámico, permite la realización de análisis de inestabilidad global y control de flujos complejos de relevancia industrial. Abstract Most flows of engineering relevance still remain unexplored in a global instability theory context for two reasons. First, because of the difficulties associated with the analysis of turbulent flows and, second, for the formidable computational resources required for the solution of the eigenvalue problem associated with the instability analysis of three-dimensional base flows, also known as TriGlobal problem. In this thesis, the problem associated with the three-dimensionality is addressed by means of the development of a general approach to the solution of large-scale global linear instability analysis by coupling a time-stepping approach with second order aerodynamic codes employed in industry. Three challenging flows in the terms of required computational resources and physical complexity have been chosen for demonstration of the present methodology; (i) the flow inside a wall-bounded three-dimensional lid-driven cavity, (ii) the flow past a cylinder fitted with helical strakes and (iii) the flow over a inhomogeneous three-dimensional open cavity. Results in excellent agreement with the literature have been obtained for the three-dimensional lid-driven cavity by using this methodology coupled with the incompressible solver of the open-source toolbox OpenFOAM®, which has served as validation. Moreover, significant physical insight of the instability of three-dimensional open flows has been gained through the application of the present time-stepping methodology to the other two cases. In addition, modifications to the present approach have been proposed in order to perform adjoint instability analysis of three-dimensional base flows and flow control; validation and TriGlobal examples are presented. Finally, it has been demonstrated that the moderate amount of computational resources required for the solution of the TriGlobal eigenvalue problem using this method enables the performance of instability analysis and control of flows of industrial relevance.

The effect of boundary conditions in the numerical solution of 3-D thermoelastic problems

After the extensive research on the capabilities of the Boundary Integral Equation Method produced during the past years the versatility of its applications has been well founded. Maybe the years to come will see the in-depth analysis of several conflictive points, for example, adaptive integration, solution of the system of equations, etc. This line is clear in academic research. In this paper we comment on the incidence of the manner of imposing the boundary conditions in 3-D coupled problems. Here the effects are particularly magnified: in the first place by the simple model used (constant elements) and secondly by the process of solution, i.e. first a potential problem is solved and then the results are used as data for an elasticity problem. The errors add to both processes and small disturbances, unimportant in separated problems, can produce serious errors in the final results. The specific problem we have chosen is especially interesting. Although more general cases (i.e. transient)can be treated, here the domain integrals can be converted into boundary ones and the influence of the manner in which boundary conditions are applied will reflect the whole importance of the problem.

A dynamic programming approach to the formulation and solution of finite element equations

A method for formulating and algorithmically solving the equations of finite element problems is presented. The method starts with a parametric partition of the domain in juxtaposed strips that permits sweeping the whole region by a sequential addition (or removal) of adjacent strips. The solution of the difference equations constructed over that grid proceeds along with the addition removal of strips in a manner resembling the transfer matrix approach, except that different rules of composition that lead to numerically stable algorithms are used for the stiffness matrices of the strips. Dynamic programming and invariant imbedding ideas underlie the construction of such rules of composition. Among other features of interest, the present methodology provides to some extent the analyst's control over the type and quantity of data to be computed. In particular, the one-sweep method presented in Section 9, with no apparent counterpart in standard methods, appears to be very efficient insofar as time and storage is concerned. The paper ends with the presentation of a numerical example

High Resolution Detection of Mechanical Forces Exerted by Locomoting Fibroblasts on the Substrate

We have developed a new approach to detect mechanical forces exerted by locomoting fibroblasts on the substrate. Cells were cultured on elastic, collagen-coated polyacrylamide sheets embedded with 0.2-μm fluorescent beads. Forces exerted by the cell cause deformation of the substrate and displacement of the beads. By recording the position of beads during cell locomotion and after cell removal, we discovered that most forces were radially distributed, switching direction in the anterior region. Deformations near the leading edge were strong, transient, and variable in magnitude, consistent with active local contractions, whereas those in the posterior region were weaker, more stable, and more uniform, consistent with passive resistance. Treatment of cells with cytochalasin D or myosin II inhibitors caused relaxation of the forces, suggesting that they are generated primarily via actin–myosin II interactions; treatment with nocodazole caused no immediate effect on forces. Immunofluorescence indicated that the frontal region of strong deformation contained many vinculin plaques but no apparent concentration of actin or myosin II filaments. Strong mechanical forces in the anterior region, generated by locally activated myosin II and transmitted through vinculin-rich structures, likely play a major role in cell locomotion and in mechanical signaling with the surrounding environment.

Using interior point algorithms for the solution of linear programs with special structural features

Linear Programming (LP) is a powerful decision making tool extensively used in various economic and engineering activities. In the early stages the success of LP was mainly due to the efficiency of the simplex method. After the appearance of Karmarkar's paper, the focus of most research was shifted to the field of interior point methods. The present work is concerned with investigating and efficiently implementing the latest techniques in this field taking sparsity into account. The performance of these implementations on different classes of LP problems is reported here. The preconditional conjugate gradient method is one of the most powerful tools for the solution of the least square problem, present in every iteration of all interior point methods. The effect of using different preconditioners on a range of problems with various condition numbers is presented. Decomposition algorithms has been one of the main fields of research in linear programming over the last few years. After reviewing the latest decomposition techniques, three promising methods were chosen the implemented. Sparsity is again a consideration and suggestions have been included to allow improvements when solving problems with these methods. Finally, experimental results on randomly generated data are reported and compared with an interior point method. The efficient implementation of the decomposition methods considered in this study requires the solution of quadratic subproblems. A review of recent work on algorithms for convex quadratic was performed. The most promising algorithms are discussed and implemented taking sparsity into account. The related performance of these algorithms on randomly generated separable and non-separable problems is also reported.

On the Operational Solution of a System of Fractional Differential Equations

MSC 2010: 26A33, 44A45, 44A40, 65J10

A simplified implementation of the least squares solution for pairwise comparisons matrices

This is a follow up to "Solution of the least squares method problem of pairwise comparisons matrix" by Bozóki published by this journal in 2008. Familiarity with this paper is essential and assumed. For lower inconsistency and decreased accuracy, our proposed solutions run in seconds instead of days. As such, they may be useful for researchers willing to use the least squares method (LSM) instead of the geometric means (GM) method.

The parallel solution of early-exercise Asian options with stochastic volatility

This paper describes an parallel semi-Lagrangian finite difference approach to the pricing of early exercise Asian Options on assets with a stochastic volatility. A multigrid procedure is described for the fast iterative solution of the discrete linear complementarity problems that result. The accuracy and performance of this approach is improved considerably by a strike-price related analytic transformation of asset prices. Asian options are contingent claims with payoffs that depend on the average price of an asset over some time interval. The payoff may depend on this average and a fixed strike price (Fixed Strike Asians) or it may depend on the average and the asset price (Floating Strike Asians). The option may also permit early exercise (American contract) or confine the holder to a fixed exercise date (European contract). The Fixed Strike Asian with early exercise is considered here where continuous arithmetic averaging has been used. Pricing such an option where the asset price has a stochastic volatility leads to the requirement to solve a tri-variate partial differential inequation in the three state variables of asset price, average price and volatility (or equivalently, variance). The similarity transformations [6] used with Floating Strike Asian options to reduce the dimensionality of the problem are not applicable to Fixed Strikes and so the numerical solution of a tri-variate problem is necessary. The computational challenge is to provide accurate solutions sufficiently quickly to support realtime trading activities at a reasonable cost in terms of hardware requirements.

Analysis Of The Fluorescence Of Body Fluids On Different Surfaces And Times.

The use of screening techniques, such as an alternative light source (ALS), is important for finding biological evidence at a crime scene. The objective of this study was to evaluate whether biological fluid (blood, semen, saliva, and urine) deposited on different surfaces changes as a function of the age of the sample. Stains were illuminated with a Megamaxx™ ALS System and photographed with a Canon EOS Utility™ camera. Adobe Photoshop™ was utilized to prepare photographs for analysis, and then ImageJ™ was used to record the brightness values of pixels in the images. Data were submitted to analysis of variance using a generalized linear mixed model with two fixed effects (surface and fluid). Time was treated as a random effect (through repeated measures) with a first-order autoregressive covariance structure. Means of significant effects were compared by the Tukey test. The fluorescence of the analyzed biological material varied depending on the age of the sample. Fluorescence was lower when the samples were moist. Fluorescence remained constant when the sample was dry, up to the maximum period analyzed (60 days), independent of the substrate on which the fluid was deposited, showing the novelty of this study. Therefore, the forensic expert can detect biological fluids at the crime scene using an ALS even several days after a crime has occurred.

Cyclic variability of the circumstellar disk of the Be star zeta Tauri II. Testing the 2D global disk oscillation model

Context. About 2/3 of the Be stars present the so-called V/R variations, a phenomenon characterized by the quasi-cyclic variation in the ratio between the violet and red emission peaks of the HI emission lines. These variations are generally explained by global oscillations in the circumstellar disk forming a one-armed spiral density pattern that precesses around the star with a period of a few years. Aims. This paper presents self-consistent models of polarimetric, photometric, spectrophotometric, and interferometric observations of the classical Be star zeta Tauri. The primary goal is to conduct a critical quantitative test of the global oscillation scenario. Methods. Detailed three-dimensional, NLTE radiative transfer calculations were carried out using the radiative transfer code HDUST. The most up-to-date research on Be stars was used as input for the code in order to include a physically realistic description for the central star and the circumstellar disk. The model adopts a rotationally deformed, gravity darkened central star, surrounded by a disk whose unperturbed state is given by a steady-state viscous decretion disk model. It is further assumed that this disk is in vertical hydrostatic equilibrium. Results. By adopting a viscous decretion disk model for zeta Tauri and a rigorous solution of the radiative transfer, a very good fit of the time-average properties of the disk was obtained. This provides strong theoretical evidence that the viscous decretion disk model is the mechanism responsible for disk formation. The global oscillation model successfully fitted spatially resolved VLTI/AMBER observations and the temporal V/R variations in the H alpha and Br gamma lines. This result convincingly demonstrates that the oscillation pattern in the disk is a one-armed spiral. Possible model shortcomings, as well as suggestions for future improvements, are also discussed.

A note on the likelihood and moments of the skew-normal distribution

In this paper an alternative approach to the one in Henze (1986) is proposed for deriving the odd moments of the skew-normal distribution considered in Azzalini (1985). The approach is based on a Pascal type triangle, which seems to greatly simplify moments computation. Moreover, it is shown that the likelihood equation for estimating the asymmetry parameter in such model is generated as orthogonal functions to the sample vector. As a consequence, conditions for a unique solution of the likelihood equation are established, which seem to hold in more general setting.