Gaussian extended cubature formulae for polyharmonic functions

The purpose of this paper is to show certain links between univariate interpolation by algebraic polynomials and the representation of polyharmonic functions. This allows us to construct cubature formulae for multivariate functions having highest order of precision with respect to the class of polyharmonic functions. We obtain a Gauss type cubature formula that uses ℳ values of linear functional (integrals over hyperspheres) and is exact for all 2ℳ-harmonic functions, and consequently, for all algebraic polynomials of n variables of degree 4ℳ - 1.

Linear quadratic Gaussian control of discrete-time Markov jump linear systems with horizon defined by stopping times

The linear quadratic Gaussian control of discrete-time Markov jump linear systems is addressed in this paper, first for state feedback, and also for dynamic output feedback using state estimation. in the model studied, the problem horizon is defined by a stopping time τ which represents either, the occurrence of a fix number N of failures or repairs (T N), or the occurrence of a crucial failure event (τ δ), after which the system paralyzed. From the constructive method used here a separation principle holds, and the solutions are given in terms of a Kalman filter and a state feedback sequence of controls. The control gains are obtained by recursions from a set of algebraic Riccati equations for the former case or by a coupled set of algebraic Riccati equation for the latter case. Copyright © 2005 IFAC.

Characteristics of the double-torsion test to determine the R-curve of ceramic materials

Double-torsion tests were carried out on a commercial ceramic floor tile to verify whether this test is suitable for determining the R-curve of ceramics. The instantaneous crack length was obtained by means of compliance calibration, and it was found that the experimental compliance underestimates the real crack length. The load vs. displacement curves were also found to drop after maximum loading, causing the stress intensity factor to decline. The R-curves were calculated by two methods: linear elastic fracture mechanics and the energetic method. It was obtained that the average values of crack resistance, R, and the double of the work of fracture, 2 · γwof, did not depend on notch length, a0, which is a highly relevant finding, indicating that these parameters were less dependent on the test specimen's geometry. The proposal was to use small notches, which produce long stable crack propagation paths that in turn are particularly important in the case of coarse microstructures.

Reduced calibration curve for proton computed tomography

The pCT deals with relatively thick targets like the human head or trunk. Thus, the fidelity of pCT as a tool for proton therapy planning depends on the accuracy of physical formulas used for proton interaction with thick absorbers. Although the actual overall accuracy of the proton stopping power in the Bethe-Bloch domain is about 1%, the analytical calculations and the Monte Carlo simulations with codes like TRIM/SRIM, MCNPX and GEANT4 do not agreed with each other. A tentative to validate the codes against experimental data for thick absorbers bring some difficulties: only a few data is available and the existing data sets have been acquired at different initial proton energies, and for different absorber materials. In this work we compare the results of our Monte Carlo simulations with existing experimental data in terms of reduced calibration curve, i.e. the range - energy dependence normalized on the range scale by the full projected CSDA range for given initial proton energy in a given material, taken from the NIST PSTAR database, and on the final proton energy scale - by the given initial energy of protons. This approach is almost energy and material independent. The results of our analysis are important for pCT development because the contradictions observed at arbitrary low initial proton energies could be easily scaled now to typical pCT energies. © 2010 American Institute of Physics.

Comparison of the GEANT4 releases 8.2 and 9.2 in terms of a pCT reduced calibration curve

The GEANT4 simulations are essential for the development of medical tomography with proton beams pCT. In the case of thin absorbers the latest releases of GEANT4 generate very similar final spectra which agree well with the results of other popular Monte Carlo codes like TRIM/SRIM, or MCNPX. For thick absorbers, however, the disagreements became evident. In a part, these disagreements are due to the known contradictions in the NIST PSTAR and SRIM reference data. Therefore, it is interesting to compare the GEANT4 results with each other, with experiment, and with diverse code results in a reduced form, which is free from this kind of doubts. In this work such comparison is done within the Reduced Calibration Curve concept elaborated for the proton beam tomography. © 2010 IEEE.

Interval multivariate curve resolution in the dereplication of HPLC-DAD data from Jatropha gossypifolia

Introduction Jatropha gossypifolia has been used quite extensively by traditional medicine for the treatment of several diseases in South America and Africa. This medicinal plant has therapeutic potential as a phytomedicine and therefore the establishment of innovative analytical methods to characterise their active components is crucial to the future development of a quality product. Objective To enhance the chromatographic resolution of HPLC-UV-diode-array detector (DAD) experiments applying chemometric tools. Methods Crude leave extracts from J. gossypifolia were analysed by HPLC-DAD. A chromatographic band deconvolution method was designed and applied using interval multivariate curve resolution by alternating least squares (MCR-ALS). Results The MCR-ALS method allowed the deconvolution from up to 117% more bands, compared with the original HPLC-DAD experiments, even in regions where the UV spectra showed high similarity. The method assisted in the dereplication of three C-glycosylflavones isomers: vitexin/isovitexin, orientin/homorientin and schaftoside/isoschaftoside. Conclusion The MCR-ALS method is shown to be a powerful tool to solve problems of chromatographic band overlapping from complex mixtures such as natural crude samples. Copyright © 2013 John Wiley & Sons, Ltd. Extracts from J. gossypifolia were analyzed by HPLC-DAD and, dereplicated applying MCR-ALS. The method assisted in the detection of three C-glycosylflavones isomers: vitexin/isovitexin, orientin/homorientin and schaftoside/isoschaftoside. The application of MCR-ALS allowed solving problems of chromatographic band overlapping from complex mixtures such as natural crude samples. Copyright © 2013 John Wiley & Sons, Ltd.

Migração 3-D Kirchhoff-Gaussian-Beam (KGB) pré-empilhamento no domínio da profundidade

O Feixe Gaussiano (FG) é uma solução assintótica da equação da elastodinâmica na vizinhança paraxial de um raio central, a qual se aproxima melhor do campo de ondas do que a aproximação de ordem zero da Teoria do Raio. A regularidade do FG na descrição do campo de ondas, assim como a sua elevada precisão em algumas regiões singulares do meio de propagação, proporciona uma forte alternativa na solução de problemas de modelagem e imageamento sísmicos. Nesta Tese, apresenta-se um novo procedimento de migração sísmica pré-empilhamento em profundidade com amplitudes verdadeiras, que combina a flexibilidade da migração tipo Kirchhoff e a robustez da migração baseada na utilização de Feixes Gaussianos para a representação do campo de ondas. O algoritmo de migração proposto é constituído por dois processos de empilhamento: o primeiro é o empilhamento de feixes (“beam stack”) aplicado a subconjuntos de dados sísmicos multiplicados por uma função peso definida de modo que o operador de empilhamento tenha a mesma forma da integral de superposição de Feixes Gaussianos; o segundo empilhamento corresponde à migração Kirchhoff tendo como entrada os dados resultantes do primeiro empilhamento. Pelo exposto justifica-se a denominação migração Kirchhoff-Gaussian-Beam (KGB). As principais características que diferenciam a migração KGB, durante a realização do primeiro empilhamento, de outros métodos de migração que também utilizam a teoria dos Feixes Gaussianos, são o uso da primeira zona de Fresnel projetada para limitar a largura do feixe e a utilização, no empilhamento do feixe, de uma aproximação de segunda ordem do tempo de trânsito de reflexão. Como exemplos são apresentadas aplicações a dados sintéticos para modelos bidimensionais (2-D) e tridimensionais (3-D), correspondentes aos modelos Marmousi e domo de sal da SEG/EAGE, respectivamente.

Análise do efeito da discretização do modelo de velocidades nas migrações Kirchhoff e Kirchhoff-Gaussian- Beam 2D pré-empilhamento em profundidade

O Feixe Gaussiano (FG) é uma solução assintótica da equação da elastodinâmica na vizinhança paraxial de um raio central, a qual se aproxima melhor do campo de ondas do que a aproximação de ordem zero da Teoria do Raio. A regularidade do FG na descrição do campo de ondas, assim como a sua elevada precisão em algumas regiões singulares do meio de propagação, proporciona uma forte alternativa no imageamento sísmicos. Nesta dissertação, apresenta-se um novo procedimento de migração sísmica pré-empilhamento em profundidade com amplitudes verdadeiras, que combina a flexibilidade da migração tipo Kirchhoff e a robustez da migração baseada na utilização de Feixes Gaussianos para a representação do campo de ondas. O algoritmo de migração proposto é constituído por dois processos de empilhamento: o primeiro é o empilhamento de feixes (“beam stack”) aplicado a subconjuntos de dados sísmicos multiplicados por uma função peso definida de modo que o operador de empilhamento tenha a mesma forma da integral de superposição de Feixes Gaussianos; o segundo empilhamento corresponde à migração Kirchhoff tendo como entrada os dados resultantes do primeiro empilhamento. Pelo exposto justifica-se a denominação migração Kirchhoff-Gaussian-Beam (KGB).Afim de comparar os métodos Kirchhoff e KGB com respeito à sensibilidade em relação ao comprimento da discretização, aplicamos no conjunto de dados conhecido como Marmousi 2-D quatro grids de velocidade, ou seja, 60m, 80m 100m e 150m. Como resultado, temos que ambos os métodos apresentam uma imagem muito melhor para o menor intervalo de discretização da malha de velocidade. O espectro de amplitude das seções migradas nos fornece o conteúdo de frequência espacial das seções das imagens obtidas.

Growth curve of female collared peccaries (Pecari tajacu) raised in captivity in the Brazilian Amazon Region

Com o objetivo de ajustar modelos não-lineares, foram utilizados registros mensais do peso de 10 fêmeas de cateto (Pecari tajacu) coletados durante dois anos, no criatório do campo experimental Álvaro Adolfo da Embrapa Amazônia Oriental, Belém, PA. Utilizaram-se os modelos de Von Bertalanffy, Brody, Gompertz e Logístico. Os parâmetros foram estimados usando o procedimento NLIN do aplicativo SAS. Os critérios utilizados para verificar o ajuste dos modelos foram: desvio padrão assintótico (ASD); coeficiente de determinação (R²); desvio médio absoluto dos resíduos (ARD) e o índice assintótico (AR). Os modelos Brody e Logístico estimaram, respectivamente, o maior (19,44kg) e o menor (19,18kg) peso assintótico (A), caracterizando a menor (0,0064kg/dia) e a maior (0,0113kg/dia) taxa de maturação (K), haja vista a natureza antagônica entre estes parâmetros, comprovada pela correlação fenotípica variando entre -0,75 à -0,47. O modelo Brody estimou o menor valor para o ARD, fator limitante para caracterizar o menor valor para o AR por este modelo. Considerando o AR, o modelo Brody apresentou o melhor ajuste, contudo, pelos valores encontrados, os demais modelos também apresentaram ajuste adequando aos dados ponderais da referida espécie/sexo. Com base no AR adotado neste trabalho, recomenda-se o modelo Brody para ajustar a curva de crescimento de fêmeas de cateto (Pecari tajacu). Em razão dos valores estimados, sobretudo, para a K, essa característica pode ser incluída em um índice de seleção. Contudo, estudos com grupos mais representativos e criados em outras condições se faz oportuno.

Multivariate curve resolution analysis applied to time-resolved synchrotron X-ray Absorption Spectroscopy monitoring of the activation of copper alumina catalyst

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Influence of salts on the coexistence curve and protein partitioning in nonionic aqueous two-phase micellar systems

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Estudo de temperatura de plasma de álcool através de espectroscopia de emissão natural radical C2*

Non-intrusive methods of diagnosis, such as spectral analysis of the radiation emitted by the system, have been used as a viable alternative for determining the temperature of combustion systems. Among them, the determination of temperature by natural emission spectroscopy has the advantage of requiring relatively simple experimental devices. Once Chemiluminescent species are formed directly in the excited state, the collection and recording of radiation emission spectrum is enough to determine the temperature (CARINHANA, 2008). In this study we used the process of making direct comparisons between the experimental spectra obtained in the laboratory from the plasma of alcohol, and the theoretical spectra plotted from a computer program developed at the IEAv. The objective was to establish a fast and reliable method to measure the rotational temperature of the radical C2*. The results showed that the temperature of the plasma, which in turn can be taken as the rotational temperature of the system, is proportional to the pressure. The temperature values ranged from ca. 2300 ~ 2500 K at a pressure of 19 mmHg to 3100 ~ 3500 K for the pressure of 46 mmHg. The temperature values are somewhat smaller when we consider the theoretical spectrum as a Lorentzian curve. The overlap of the spectra was better when using the profile curve, but still were not exactly superimposed. The solution to improve the overlap of the theoretical with the experimental spectra is the use of a curve that has the convolution of two profiles analyzed: Lorentzian and Gaussian. This curve is called the Voigt profile, which will also be implemented by programmers and studied in a next work

Comparison of growth curve parameters of organs and body components in meat- (Coturnix coturnix coturnix) and laying-type (Coturnix coturnix japonica) quail show interactions between gender and genotype

The objective of this study was to estimate growth parameters of carcass components (wing, thighs and drumsticks, back and breast) and organs (heart, liver, gizzard and gut) in males and females of one meat-type quail strain (Coturnix coturnix coturnix) and two laying strains (Coturnix coturnix japonica) designated either yellow or red.A total of 1350 quail from 1 to 42d old were distributed in a completely randomised design, with 5 replicates of each strain. The carcass component weights and body organs were analysed weekly and evaluated using the Gompertz function; growth rates were evaluated through derivative equations.The meat-type strain presented the highest growth rates in carcass components and organs. Across strains, females showed the highest weight of internal organs at maturity compared to males.Females had greater growth potential in breast, wings and back than males for both yellow and red laying quail.

The Performance of Elliptic Curve Based Group Diffie-Hellman Protocols for Secure Group Communication over Ad Hoc Networks

The security of the two party Diffie-Hellman key exchange protocol is currently based on the discrete logarithm problem (DLP). However, it can also be built upon the elliptic curve discrete logarithm problem (ECDLP). Most proposed secure group communication schemes employ the DLP-based Diffie-Hellman protocol. This paper proposes the ECDLP-based Diffie-Hellman protocols for secure group communication and evaluates their performance on wireless ad hoc networks. The proposed schemes are compared at the same security level with DLP-based group protocols under different channel conditions. Our experiments and analysis show that the Tree-based Group Elliptic Curve Diffie-Hellman (TGECDH) protocol is the best in overall performance for secure group communication among the four schemes discussed in the paper. Low communication overhead, relatively low computation load and short packets are the main reasons for the good performance of the TGECDH protocol.