Prediction of dynamical, nonlinear, and unstable process behavior

Process scheduling techniques consider the current load situation to allocate computing resources. Those techniques make approximations such as the average of communication, processing, and memory access to improve the process scheduling, although processes may present different behaviors during their whole execution. They may start with high communication requirements and later just processing. By discovering how processes behave over time, we believe it is possible to improve the resource allocation. This has motivated this paper which adopts chaos theory concepts and nonlinear prediction techniques in order to model and predict process behavior. Results confirm the radial basis function technique which presents good predictions and also low processing demands show what is essential in a real distributed environment.

An improved finite element space for discontinuous pressures

We consider incompressible Stokes flow with an internal interface at which the pressure is discontinuous, as happens for example in problems involving surface tension. We assume that the mesh does not follow the interface, which makes classical interpolation spaces to yield suboptimal convergence rates (typically, the interpolation error in the L(2)(Omega)-norm is of order h(1/2)). We propose a modification of the P(1)-conforming space that accommodates discontinuities at the interface without introducing additional degrees of freedom or modifying the sparsity pattern of the linear system. The unknowns are the pressure values at the vertices of the mesh and the basis functions are computed locally at each element, so that the implementation of the proposed space into existing codes is straightforward. With this modification, numerical tests show that the interpolation order improves to O(h(3/2)). The new pressure space is implemented for the stable P(1)(+)/P(1) mini-element discretization, and for the stabilized equal-order P(1)/P(1) discretization. Assessment is carried out for Poiseuille flow with a forcing surface and for a static bubble. In all cases the proposed pressure space leads to improved convergence orders and to more accurate results than the standard P(1) space. In addition, two Navier-Stokes simulations with moving interfaces (Rayleigh-Taylor instability and merging bubbles) are reported to show that the proposed space is robust enough to carry out realistic simulations. (c) 2009 Elsevier B.V. All rights reserved.

An unstructured CVFEM and moving interface algorithm for non-Newtonian Hele-Shaw flows in injection molding

Purpose - The purpose of this paper is to develop a novel unstructured simulation approach for injection molding processes described by the Hele-Shaw model. Design/methodology/approach - The scheme involves dual dynamic meshes with active and inactive cells determined from an initial background pointset. The quasi-static pressure solution in each timestep for this evolving unstructured mesh system is approximated using a control volume finite element method formulation coupled to a corresponding modified volume of fluid method. The flow is considered to be isothermal and non-Newtonian. Findings - Supporting numerical tests and performance studies for polystyrene described by Carreau, Cross, Ellis and Power-law fluid models are conducted. Results for the present method are shown to be comparable to those from other methods for both Newtonian fluid and polystyrene fluid injected in different mold geometries. Research limitations/implications - With respect to the methodology, the background pointset infers a mesh that is dynamically reconstructed here, and there are a number of efficiency issues and improvements that would be relevant to industrial applications. For instance, one can use the pointset to construct special bases and invoke a so-called ""meshless"" scheme using the basis. This would require some interesting strategies to deal with the dynamic point enrichment of the moving front that could benefit from the present front treatment strategy. There are also issues related to mass conservation and fill-time errors that might be addressed by introducing suitable projections. The general question of ""rate of convergence"" of these schemes requires analysis. Numerical results here suggest first-order accuracy and are consistent with the approximations made, but theoretical results are not available yet for these methods. Originality/value - This novel unstructured simulation approach involves dual meshes with active and inactive cells determined from an initial background pointset: local active dual patches are constructed ""on-the-fly"" for each ""active point"" to form a dynamic virtual mesh of active elements that evolves with the moving interface.

Integrable maps with non-trivial topology: application to divertor configurations

We explore a method for constructing two-dimensional area-preserving, integrable maps associated with Hamiltonian systems, with a given set of fixed points and given invariant curves. The method is used to find an integrable Poincare map for the field lines in a large aspect ratio tokamak with a poloidal single-null divertor. The divertor field is a superposition of a magnetohydrodynamic equilibrium with an arbitrarily chosen safety factor profile, with a wire carrying an electric current to create an X-point. This integrable map is perturbed by an impulsive perturbation that describes non-axisymmetric magnetic resonances at the plasma edge. The non-integrable perturbed map is applied to study the structure of the open field lines in the scrape-off layer, reproducing the main transport features obtained by integrating numerically the magnetic field line equations, such as the connection lengths and magnetic footprints on the divertor plate.

A Trotter-Kato product formula for a class of non-autonomous evolution equations

In this article dedicated to Professor V. Lakshmikantham on the occasion of the celebration of his 84th birthday, we announce new results concerning the existence and various properties of an evolution system UA+B(t, s)(0 <= s <= t <= T) generated by the sum -(A(t)+B(t)) of two linear, time-dependent and generally unbounded operators defined on time-dependent domains in a complex and separable Banach space B. In particular, writing G(B) for the algebra of all linear bounded operators on B, we can express UA+B(t, s)(0 <= s <= t <= T) as the strong limit in L(B) of a product of the holomorphic contraction semigroups generated by -A(t) and -B(t), thereby getting a product formula of the Trotter-Kato type under very general conditions which allow the domain D(A(t)+B(t)) to evolve with time provided there exists a fixed set D subset of boolean AND D-t epsilon[0,D-T](A(t)+B(t)) everywhere dense in B. We then mention several possible applications of our product formula to various classes of non-autonomous parabolic initial-boundary value problems, as well as to evolution problems of Schrodinger type related to the theory of time-dependent singular perturbations of self-adjoint operators in quantum mechanics. We defer all the proofs and all the details of the applications to a separate publication. (C) 2008 Elsevier Ltd. All rights reserved.

Time-resolved study electronic and thermal contributions to the nonlinear refractive index of Nd3+: SBN laser crystals

The propagation of an optical beam through dielectric media induces changes in the refractive index, An, which causes self-focusing or self-defocusing. In the particular case of ion-doped solids, there are thermal and non-thermal lens effects, where the latter is due to the polarizability difference, Delta alpha, between the excited and ground states, the so-called population lens (PL) effect. PL is a pure electronic contribution to the nonlinearity, while the thermal lens (TL) effect is caused by the conversion of part of the absorbed energy into heat. In time-resolved measurements such as Z-scan and TL transient experiments, it is not easy to separate these two contributions to nonlinear refractive index because they usually have similar response times. In this work, we performed time-resolved measurements using both Z-scan and mode mismatched TL in order to discriminate thermal and electronic contributions to the laser-induced refractive index change of the Nd3+-doped Strontium Barium Niobate (SrxBa1-xNb2O6) laser crystal. Combining numerical simulations with experimental results we could successfully distinguish between the two contributions to An. (C) 2007 Elsevier B.V. All rights reserved.

Asymptotic Skewness in Exponential Family Nonlinear Models

In this article, we give an asymptotic formula of order n(-1/2), where n is the sample size, for the skewness of the distributions of the maximum likelihood estimates of the parameters in exponencial family nonlinear models. We generalize the result by Cordeiro and Cordeiro ( 2001). The formula is given in matrix notation and is very suitable for computer implementation and to obtain closed form expressions for a great variety of models. Some special cases and two applications are discussed.

A criterion for unitary similarity of upper triangular matrices in general position

Each square complex matrix is unitarily similar to an upper triangular matrix with diagonal entries in any prescribed order. Let A = [a(ij)] and B = [b(ij)] be upper triangular n x n matrices that are not similar to direct sums of square matrices of smaller sizes, or are in general position and have the same main diagonal. We prove that A and B are unitarily similar if and only if parallel to h(A(k))parallel to = parallel to h(B(k))parallel to for all h is an element of C vertical bar x vertical bar and k = 1, ..., n, where A(k) := [a(ij)](i.j=1)(k) and B(k) := [b(ij)](i.j=1)(k) are the leading principal k x k submatrices of A and B, and parallel to . parallel to is the Frobenius norm. (C) 2011 Elsevier Inc. All rights reserved.

Movements and mechanical stresses in gas-filled flat plate solar collectors

Sealed gas filled flat plate solar collectors will have stresses in the material since volume and pressure varies in the gas when the temperature changes. Several geometries were analyzed and it could be seen that it is possible reducing the stresses and improve the safety factor of the weakest point in the construction by using larger area and/or reducing the distance between glass and absorber and/or change width and height relationship so the tubes are getting longer. Further it could be shown that the safety factor won't always get improved with reinforcements. It is so because when an already strong part of the collector gets reinforced it will expose weaker parts for higher stresses. The finite element method was used for finding out the stresses.

Testing common nonlinear features in vector nonlinear autoregressive models

This paper studies a special class of vector smooth-transition autoregressive (VSTAR) models that contains common nonlinear features (CNFs), for which we proposed a triangular representation and developed a procedure of testing CNFs in a VSTAR model. We first test a unit root against a stable STAR process for each individual time series and then examine whether CNFs exist in the system by Lagrange Multiplier (LM) test if unit root is rejected in the first step. The LM test has standard Chi-squared asymptotic distribution. The critical values of our unit root tests and small-sample properties of the F form of our LM test are studied by Monte Carlo simulations. We illustrate how to test and model CNFs using the monthly growth of consumption and income data of United States (1985:1 to 2011:11).

Comportamento de fundações superficiais apoiadas em aterros estabilizados com resíduos industriais

Neste trabalho foi estudado o comportamento de fundações superficiais apoiadas em sistema de dupla camada, quando a superior é cimentada. O estudo consistiu-se de três etapas, chamadas de Etapa de Laboratório, Etapa Numérica e Etapa de Campo. Na Etapa de Laboratório foi verificada a viabilidade técnica de utilizar os resíduos industriais cinza pesada e cal de carbureto na estabilização de um solo residual de arenito botucatu. Estudou-se a reatividade da cinza pesada com a cal de carbureto, a influência da temperatura e do tempo de cura no desenvolvimento das reações pozolânicas, a influência de diferentes teores de resíduos na resistência à compressão simples, compressão diametral e durabilidade, objetivando definir uma mistura ótima e, ainda, o impacto ambiental da utilização da mistura ótima, através de ensaios de lixiviação e solubilização. Na Etapa Numérica foi estudado, através do Método dos Elementos Finitos, o comportamento de fundações superficiais apoiadas em dupla camada. O modelo utilizado para representar o comportamento do material cimentado e não-cimentado foi o elástico-plástico com critério de ruptura de Drucker-Prager e fluxo não-associado. Verificou-se, através de análise paramétrica, a influência da espessura da camada cimentada e do diâmetro da fundação, bem como a influência dos parâmetros dos materiais cimentado e não-cimentado na resposta carga x recalque de fundações superficiais. Na Etapa de Campo foram construídos aterros experimentais utilizando a mistura ótima determinada na Etapa de Laboratório e, sobre estes aterros, foram executados provas de carga de placas. A análise dos resultados obtidos nas três etapas levou às seguintes conclusões: é possível utilizar cinza pesada e cal de carbureto para estabilizar o solo residual de botucatu; o comportamento de fundações superficiais sobre solos cimentados é controlado pela relação espessura da camada cimentada diâmetro da fundação; os parâmetros ângulo de atrito e módulo de elasticidade da camada cimentada não influenciam os resultados de prova de carga; a ruptura da fundação é função de dois mecanismos progressivos, os quais são função das tensões de tração geradas na parte inferior da camada cimentada e das tensões cisalhantes existentes logo abaixo das bordas da fundação.

Low cost BIST techniques for linear and non-linear analog circuits

With the ever increasing demands for high complexity consumer electronic products, market pressures demand faster product development and lower cost. SoCbased design can provide the required design flexibility and speed by allowing the use of IP cores. However, testing costs in the SoC environment can reach a substantial percent of the total production cost. Analog testing costs may dominate the total test cost, as testing of analog circuits usually require functional verification of the circuit and special testing procedures. For RF analog circuits commonly used in wireless applications, testing is further complicated because of the high frequencies involved. In summary, reducing analog test cost is of major importance in the electronic industry today. BIST techniques for analog circuits, though potentially able to solve the analog test cost problem, have some limitations. Some techniques are circuit dependent, requiring reconfiguration of the circuit being tested, and are generally not usable in RF circuits. In the SoC environment, as processing and memory resources are available, they could be used in the test. However, the overhead for adding additional AD and DA converters may be too costly for most systems, and analog routing of signals may not be feasible and may introduce signal distortion. In this work a simple and low cost digitizer is used instead of an ADC in order to enable analog testing strategies to be implemented in a SoC environment. Thanks to the low analog area overhead of the converter, multiple analog test points can be observed and specific analog test strategies can be enabled. As the digitizer is always connected to the analog test point, it is not necessary to include muxes and switches that would degrade the signal path. For RF analog circuits, this is specially useful, as the circuit impedance is fixed and the influence of the digitizer can be accounted for in the design phase. Thanks to the simplicity of the converter, it is able to reach higher frequencies, and enables the implementation of low cost RF test strategies. The digitizer has been applied successfully in the testing of both low frequency and RF analog circuits. Also, as testing is based on frequency-domain characteristics, nonlinear characteristics like intermodulation products can also be evaluated. Specifically, practical results were obtained for prototyped base band filters and a 100MHz mixer. The application of the converter for noise figure evaluation was also addressed, and experimental results for low frequency amplifiers using conventional opamps were obtained. The proposed method is able to enhance the testability of current mixed-signal designs, being suitable for the SoC environment used in many industrial products nowadays.

Moral hazard and nonlinear pricing in a general equilibrium model

The paper analyzes a two period general equilibrium model with individual risk and moral hazard. Each household faces two individual states of nature in the second period. These states solely differ in the household's vector of initial endowments, which is strictly larger in the first state (good state) than in the second state (bad state). In the first period households choose a non-observable action. Higher leveis of action give higher probability of the good state of nature to occur, but lower leveIs of utility. Households have access to an insurance market that allows transfer of income across states of oature. I consider two models of financiaI markets, the price-taking behavior model and the nonlínear pricing modelo In the price-taking behavior model suppliers of insurance have a belief about each household's actíon and take asset prices as given. A variation of standard arguments shows the existence of a rational expectations equilibrium. For a generic set of economies every equilibrium is constraíned sub-optímal: there are commodity prices and a reallocation of financiaI assets satisfying the first period budget constraint such that, at each household's optimal choice given those prices and asset reallocation, markets clear and every household's welfare improves. In the nonlinear pricing model suppliers of insurance behave strategically offering nonlinear pricing contracts to the households. I provide sufficient conditions for the existence of equilibrium and investigate the optimality properties of the modeI. If there is a single commodity then every equilibrium is constrained optimaI. Ir there is more than one commodity, then for a generic set of economies every equilibrium is constrained sub-optimaI.

Inorganic pyrophosphate-phosphohydrolytic activity associated with rat osseous plate alkaline phosphatase

Purified membrane-bound alkaline phosphatase from rat osseous plate hydrolyzed pyrophosphate in the presence of magnesium ions, with a specific activity of 92.7 U/mg. Optimal apparent pH for pyrophosphatase activity was 8.0 and it remained unchanged on increasing the pyrophosphate concentration. In the absence of magnesium ions the enzyme had a K-m = 88 mu M and V = 36.7 U/mg for pyrophosphate and no inhibition by excess substrate was observed. Pyrophosphatase activity was rapidly destroyed at temperatures above 40 degrees C, but magnesium ions apparently protected the enzyme against danaturation. Sodium metavanadate (Ki = 1.0 mM) was a competitive inhibitor of pyrophosphatase activity, while levamisole (Ki = 8.2 mM) and theophylline (Ki = 7.4 mM) were uncompetitive inhibitors. Magnesium ions (K-0.5 = 1.7 mu M) stimulated pyrophosphatase activity, while cobalt (Ki = 48.5 mu M) and zinc (Ki = 22.0 mu M) ions were non-competitive inhibitors. Manganese and calcium ions had no effect on pyrophosphatase activity. The M-w of the pyrophosphatase: protein was 130 kDa by gel filtration, but a value of 65 kDa was obtained by dissociative gel electrophoresis, suggesting that it was a dimer of apparently identical subunits. These results suggested that pyrophosphatase activity stems from the membrane-bound osseous plate alkaline phosphatase and not from a different protein.

Comparison of non-linear growth models to describe the growth curve in West African Dwarf sheep

The objectives of this study were to compare the goodness of fit of four non-linear growth models, i.e. Brody, Gompertz, Logistic and Von Bertalanffy, in West African Dwarf (WAD) sheep. A total of 5274 monthly weight records from birth up to 180 days of age from 889 lambs, collected during 2001 to 2004 in Betecoucou breeding farm in Benin were used. In the preliminary analysis, the General Linear Model Procedure of the Statistical Analysis Systems Institute was applied to the dataset to identify the significant effects of the sex of lamb (male and female), type of birth (single and twin), season of birth (rainy season and dry season), parity of dam (1, 2 and 3) and year of birth (2001, 2002, 2003 and 2004) on the observed birth weight and monthly weight up to 6 months of age. The models parameters (A, B and k), coefficient of determination (112), mean square error (MSE) were calculated using language of technical computing package Matlab(R), 2006. The mean values of A, B and k were substituted into each model to calculate the corresponding Akaike's Information Criterion (AIC). Among the four growth functions, the Brody model has been selected for its accuracy of fit according to the higher R(2), lower MSE and A/C Finally, the parameters A, B and k were adjusted in Matlab(R) 2006 for the sex of lamb, year of birth, season of birth, birth type and the parity of ewe, providing a specific slope of the Brody growth curve. The results of this study suggest that Brody model can be useful for WAD sheep breeding in Betecoucou farm conditions through growth monitoring.