A simple way to introduce fibers into FEM models

This communication proposes a simple way to introduce fibers into finite element modelling. This is a promising formulation to deal with fiber-reinforced composites by the finite element method (FEM), as it allows the consideration of short or long fibers placed arbitrarily inside a continuum domain (matrix). The most important feature of the formulation is that no additional degree of freedom is introduced into the pre-existent finite element numerical system to consider any distribution of fiber inclusions. In other words, the size of the system of equations used to solve a non-reinforced medium is the same as the one used to solve the reinforced counterpart. Another important characteristic is the reduced work required by the user to introduce fibers, avoiding `rebar` elements, node-by-node geometrical definitions or even complex mesh generation. An additional characteristic of the technique is the possibility of representing unbounded stresses at the end of fibers using a finite number of degrees of freedom. Further studies are required for non-linear applications in which localization may occur. Along the text the linear formulation is presented and the bounded connection between fibers and continuum is considered. Four examples are presented, including non-linear analysis, to validate and show the capabilities of the formulation. Copyright (c) 2007 John Wiley & Sons, Ltd.

Optimal synthesis of anaerobic digester networks

The objective of this paper is to develop a mathematical model for the synthesis of anaerobic digester networks based on the optimization of a superstructure that relies on a non-linear programming formulation. The proposed model contains the kinetic and hydraulic equations developed by Pontes and Pinto [Chemical Engineering journal 122 (2006) 65-80] for two types of digesters, namely UASB (Upflow Anaerobic Sludge Blanket) and EGSB (Expanded Granular Sludge Bed) reactors. The objective function minimizes the overall sum of the reactor volumes. The optimization results show that a recycle stream is only effective in case of a reactor with short-circuit, such as the UASB reactor. Sensitivity analysis was performed in the one and two-digester network superstructures, for the following parameters: UASB reactor short-circuit fraction and the EGSB reactor maximum organic load, and the corresponding results vary considerably in terms of digester volumes. Scenarios for three and four-digester network superstructures were optimized and compared with the results from fewer digesters. (C) 2009 Elsevier B.V. All rights reserved.

Age-related gray matter volume changes in the brain during non-elderly adulthood

Previous magnetic resonance imaging (MRI) studies described consistent age-related gray matter (GM) reductions in the fronto-parietal neocortex, insula and cerebellum in elderly subjects, but not as frequently in limbic/paralimbic structures. However, it is unclear whether such features are already present during earlier stages of adulthood, and if age-related GM changes may follow non-linear patterns at such age range. This voxel-based morphometry study investigated the relationship between GM volumes and age specifically during non-elderly life (18-50 years) in 89 healthy individuals (48 males and 41 females). Voxelwise analyses showed significant (p < 0.05, corrected) negative correlations in the right prefrontal cortex and left cerebellum, and positive correlations (indicating lack of GM loss) in the medial temporal region, cingulate gyrus, insula and temporal neocortex. Analyses using ROI masks showed that age-related dorsolateral prefrontal volume decrements followed non-linear patterns, and were less prominent in females compared to males at this age range. These findings further support for the notion of a heterogeneous and asynchronous pattern of age-related brain morphometric changes, with region-specific non-linear features. (C) 2009 Elsevier Inc. All rights reserved.

Interaction of equatorial waves through resonance with the diurnal cycle of tropical heating

In this work, a new theoretical mechanism is presented in which equatorial Rossby and inertio-gravity wave modes may interact with each other through resonance with the diurnal cycle of tropical deep convection. We have adopted the two-layer incompressible equatorial primitive equations forced by a parametric heating that roughly represents deep convection activity in the tropical atmosphere. The heat source was parametrized in the simplest way according to the hypothesis that it is proportional to the lower-troposphere moisture convergence, with the background moisture state function mimicking the structure of the ITCZ. In this context, we have investigated the possibility of resonant interaction between equatorially trapped Rossby and inertio-gravity modes through the diurnal cycle of the background moisture state function. The reduced dynamics of a single resonant duo shows that when this diurnal variation is considered, a Rossby wave mode can undergo significant amplitude modulations when interacting with an inertio-gravity wave mode, which is not possible in the context of the resonant triad non-linear interaction. Therefore, the results suggest that the diurnal variation of the ITCZ can be a possible dynamical mechanism that leads the Rossby waves to be significantly affected by high frequency modes.

Black-box decomposition approach for computational hemodynamics: One-dimensional models

Increasing efforts exist in integrating different levels of detail in models of the cardiovascular system. For instance, one-dimensional representations are employed to model the systemic circulation. In this context, effective and black-box-type decomposition strategies for one-dimensional networks are needed, so as to: (i) employ domain decomposition strategies for large systemic models (1D-1D coupling) and (ii) provide the conceptual basis for dimensionally-heterogeneous representations (1D-3D coupling, among various possibilities). The strategy proposed in this article works for both of these two scenarios, though the several applications shown to illustrate its performance focus on the 1D-1D coupling case. A one-dimensional network is decomposed in such a way that each coupling point connects two (and not more) of the sub-networks. At each of the M connection points two unknowns are defined: the flow rate and pressure. These 2M unknowns are determined by 2M equations, since each sub-network provides one (non-linear) equation per coupling point. It is shown how to build the 2M x 2M non-linear system with arbitrary and independent choice of boundary conditions for each of the sub-networks. The idea is then to solve this non-linear system until convergence, which guarantees strong coupling of the complete network. In other words, if the non-linear solver converges at each time step, the solution coincides with what would be obtained by monolithically modeling the whole network. The decomposition thus imposes no stability restriction on the choice of the time step size. Effective iterative strategies for the non-linear system that preserve the black-box character of the decomposition are then explored. Several variants of matrix-free Broyden`s and Newton-GMRES algorithms are assessed as numerical solvers by comparing their performance on sub-critical wave propagation problems which range from academic test cases to realistic cardiovascular applications. A specific variant of Broyden`s algorithm is identified and recommended on the basis of its computer cost and reliability. (C) 2010 Elsevier B.V. All rights reserved.

Sound waves and solitons in hot and dense nuclear matter

Assuming that nuclear matter can be treated as a perfect fluid, we study the propagation of perturbations in the baryon density. The equation of state is derived from a relativistic mean field model, which is a variant of the non-linear Walecka model. The expansion of the Euler and continuity equations of relativistic hydrodynamics around equilibrium configurations leads to differential equations for the density perturbation. We solve them numerically for linear and spherical perturbations and follow the propagation of the initial pulses. For linear perturbations we find single soliton solutions and solutions with one or more solitons followed by ""radiation"". Depending on the equation of state a strong damping may occur. We consider also the evolution of perturbations in a medium without dispersive effects. In this case we observe the formation and breaking of shock waves. We study all these equations also for matter at finite temperature. Our results may be relevant for the analysis of RHIC data. They suggest that the shock waves formed in the quark gluon plasma phase may survive and propagate in the hadronic phase. (C) 2009 Elseiver. B.V. All rights reserved.

Shock wave formation in hot nuclear matter

Assuming that nuclear matter can be treated as a perfect fluid, we study the propagation of perturbations in the baryon density at high temperature. The equation of state is derived from the non-linear Walecka model. The expansion of the Euler and continuity equations of relativistic hydrodynamics around equilibrium configurations lead to the breaking wave equation for the density perturbation. We solve it numerically for this perturbation and follow the propagation of the initial pulses.

Static Hopfions in the extended Skyrme-Faddeev model

We construct static soliton solutions with non-zero Hopf topological charges to a theory which is an extension of the Skyrme-Faddeev model by the addition of a further quartic term in derivatives. We use an axially symmetric ansatz based on toroidal coordinates, and solve the resulting two coupled non-linear partial differential equations in two variables by a successive over-relaxation (SOR) method. We construct numerical solutions with Hopf charge up to four, and calculate their analytical behavior in some limiting cases. The solutions present an interesting behavior under the changes of a special combination of the coupling constants of the quartic terms. Their energies and sizes tend to zero as that combination approaches a particular special value. We calculate the equivalent of the Vakulenko and Kapitanskii energy bound for the theory and find that it vanishes at that same special value of the coupling constants. In addition, the model presents an integrable sector with an in finite number of local conserved currents which apparently are not related to symmetries of the action. In the intersection of those two special sectors the theory possesses exact vortex solutions (static and time dependent) which were constructed in a previous paper by one of the authors. It is believed that such model describes some aspects of the low energy limit of the pure SU(2) Yang-Mills theory, and our results may be important in identifying important structures in that strong coupling regime.

ENO adaptive method for solving one-dimensional conservation laws

In this work an efficient third order non-linear finite difference scheme for solving adaptively hyperbolic systems of one-dimensional conservation laws is developed. The method is based oil applying to the solution of the differential equation an interpolating wavelet transform at each time step, generating a multilevel representation for the solution, which is thresholded and a sparse point representation is generated. The numerical fluxes obtained by a Lax-Friedrichs flux splitting are evaluated oil the sparse grid by an essentially non-oscillatory (ENO) approximation, which chooses the locally smoothest stencil among all the possibilities for each point of the sparse grid. The time evolution of the differential operator is done on this sparse representation by a total variation diminishing (TVD) Runge-Kutta method. Four classical examples of initial value problems for the Euler equations of gas dynamics are accurately solved and their sparse solutions are analyzed with respect to the threshold parameters, confirming the efficiency of the wavelet transform as an adaptive grid generation technique. (C) 2008 IMACS. Published by Elsevier B.V. All rights reserved.

Probing the dependence of the properties of cellulose acetates and their films on the degree of biopolymer substitution: use of solvatochromic indicators and thermal analysis

Although cellulose acetates, CAs, are extensively employed there is scant information about the systematic dependence of their properties on their degree of substitution, DS; this is the subject of the present work. Nine CAs samples, DS from 0.83 to 3.0 were synthesized; their films were prepared. The following solvatochromic probes have been employed in order to determine the empirical polarity, E (T)(33); ""acidity, alpha""; ""basicity, beta"", and ""dipolarity/polarizability, pi*"" of the casted films: 2,6-dichloro-4-(2,4,6-triphenyl-pyridinium-1-yl) phenolate, WB; 4-nitroaniline; 4-nitroanisole; 4-nitro-N,N-dimethylaniline; 2,6-diphenyl-4-(2,4,6-triphenyl-pyridinium-1-yl)phenolate, RB. Additionally, two systems, ethanol plus ethyl acetate (EtOH-EtAc), and cellulose plus cellulose triacetate, CTA, were employed as models for CAs of different DS. Regarding the model systems, the following was observed: (i) For EtOH-EtAc, the dependence of all solvatochromic parameters on the ""equivalent-DS"" of the binary mixture was non-linear because of preferential solvation; (ii) The dependence of E (T)(33) on equivalent DS of the cellulose-CTA films is linear, but the slope is smaller than that of the corresponding plot for CAs. This is attributed to the more efficient hydrogen bonding in the model system, a conclusion corroborated by IR measurements. The dependence of solvatochromic parameters of CAs on their DS is described by the simple equations; a consequence of the substitution of the OH by the ester group. The thermal properties of bulk CAs samples were investigated by DSC and TGA; their dependence on DS is described by simple equations. The relevance of these data to the processing and applications of CAs is briefly discussed.

Teoria para o trabalho de expansão aplicada às brisas do nordeste Brasileiro

Neste artigo investiga-se numericamente a interação entre brisas marítima-terrestre e de vale-montanha, que ocorrem no nordeste brasileiro devido à presença do Planalto da Borborema, aplicando o trabalho de expansão associado ao ramo inferior das circulações. Resultados da teoria da máquina térmica são comparados aos resultados 3D, obtidos com a versão brasileira do modelo RAMS. Os resultados indicaram que o efeito do contraste de temperatura no trabalho ligado às circulações (Wa), isoladamente, contribui para a formação de brisas marítimas mais intensas e de brisas terrestres menos intensas. Na realidade, o que se observa são brisas terrestres com intensidades iguais ou até maiores que das brisas marítimas. Assim, a contribuição da montanha para a intensidade das circulações de brisa no período noturno mostra-se extremamente não-linear. O trabalho de expansão que realmente está ligado às circulações, contribui em apenas 7% para o trabalho total. Dessa forma, a maior parte do trabalho total está associada ao trabalho de compressão que a atmosfera realiza para compensar a perda por resfriamento radiativo, e muito da energia disponibilizada para as circulações é gasta para vencer os processos dissipativos.

Trade unions and the economic performance of brazilian establishments

Este artigo analisa, pela primeira vez na literatura, o impacto dos sindicatos de trabalhadores em vários indicadores de desempenho econômico de firmas industriais brasileiras. Realizou-se uma pesquisa retrospectiva sobre a densidade sindical de 1000 estabelecimentos industriais brasileiros e seus resultados foram combinados aos indicadores de desempenho econômico da Pesquisa Industrial Anual (PIA) de 1990 a 2000. Os resultados indicam que a relação entre a densidade sindical na firma e seus salários, emprego e produtividade, é não-linear, ou seja, um aumento no grau de sindicalização leva a um melhor desempenho, porém a taxas decrescentes. Observou-se, também, uma relação negativa entre sindicalização e rentabilidade. Finalmente, estabelecimentos que introduziram mecanismos de 'participação nos lucros' aumentaram sua produtividade e rentabilidade no período e pagaram maiores salários nas firmas onde o grau de sindicalização era maior.

Numerical investigation of the three-dimensional secondary instabilities in the time-developing compressible mixing layer

Mixing layers are present in very different types of physical situations such as atmospheric flows, aerodynamics and combustion. It is, therefore, a well researched subject, but there are aspects that require further studies. Here the instability of two-and three-dimensional perturbations in the compressible mixing layer was investigated by numerical simulations. In the numerical code, the derivatives were discretized using high-order compact finite-difference schemes. A stretching in the normal direction was implemented with both the objective of reducing the sound waves generated by the shear region and improving the resolution near the center. The compact schemes were modified to work with non-uniform grids. Numerical tests started with an analysis of the growth rate in the linear regime to verify the code implementation. Tests were also performed in the non-linear regime and it was possible to reproduce the vortex roll-up and pairing, both in two-and three-dimensional situations. Amplification rate analysis was also performed for the secondary instability of this flow. It was found that, for essentially incompressible flow, maximum growth rates occurred for a spanwise wavelength of approximately 2/3 of the streamwise spacing of the vortices. The result demonstrated the applicability of the theory developed by Pierrehumbet and Widnall. Compressibility effects were then considered and the maximum growth rates obtained for relatively high Mach numbers (typically under 0.8) were also presented.

Effect of therapeutic dose X rays on mechanical and chemical properties of esthetic dental materials

The aim of this study was to investigate the influence of therapeutic dose X rays on the microhardness (MH) and degree of conversion (DC) of two different esthetic restorative dental materials. The materials were photo-activated with a LED light-curing unit using three cure-times: 5, 20 and 40 seconds. The photo-activation was carried out in two distinct periods: before and after irradiation with doses of 5, 35 and 70 Gy, from a 6 MV X rays beam. In accordance with the methodology used, it was conclude that a therapeutic dose does not have a detrimental effect on the photoinitiator molecules, because the photo-activation occurred after they were irradiated. When the irradiation was applied before photo-activation, the materials showed MH improvement, but when photo-activation was performed after irradiation, there was less improvement. However, there was no correlation between MH and DC. Thus, a therapeutic dose applied to cured material can promote linking and breaking of chain bonds in a non-linear way.

Comparison of some theoretical models for fittings of the temperature dependence of the fundamental energy gap in GaAs

In this work we report on a comparison of some theoretical models usually used to fit the dependence on temperature of the fundamental energy gap of semiconductor materials. We used in our investigations the theoretical models of Viña, Pässler-p and Pässler-ρ to fit several sets of experimental data, available in the literature for the energy gap of GaAs in the temperature range from 12 to 974 K. Performing several fittings for different values of the upper limit of the analyzed temperature range (Tmax), we were able to follow in a systematic way the evolution of the fitting parameters up to the limit of high temperatures and make a comparison between the zero-point values obtained from the different models by extrapolating the linear dependence of the gaps at high T to T = 0 K and that determined by the dependence of the gap on isotope mass. Using experimental data measured by absorption spectroscopy, we observed the non-linear behavior of Eg(T) of GaAs for T > ΘD.