A thermoelastic system of memory type in noncylindrical domains

In this paper, we consider a hyperbolic thermoelastic system of memory type in domains with moving boundary. The problem models vibrations of an elastic bar under thermal effects according to the heat conduction law of Gurtin and Pipkin. Global existence is proved by using the penalty method of Lions. (c) 2007 Elsevier Inc. All rights reserved.

THE ROLE OF THE SUPERIOR COLLICULUS IN PREDATORY HUNTING

Combining the results of behavioral, neuronal immediate early gene activation, lesion and neuroanatomical experiments, we have presently investigated the role of the superior colliculus (SC) in predatory hunting. First, we have shown that insect hunting is associated with a characteristic large increase in Fos expression in the lateral part of the intermediate gray layer of the SC (Wig). Next, we have shown that animals with bilateral NMDA lesions of the lateral parts of the SC presented a significant delay in starting to chase the prey and longer periods engaged in other activities than predatory hunting. They also showed a clear deficit to orient themselves toward the moving prey and lost the stereotyped sequence of actions seen for capturing, holding and killing the prey. Our Phaseolus vulgaris-leucoagglutinin analysis revealed that the lateral SCig, besides providing the well-documented descending crossed pathway to premotor sites in brainstem and spinal cord, projects to a number of midbrain and diencephalic sites likely to influence key functions in the context of the predatory behavior, such as general levels of arousal, motivational level to hunt or forage, behavioral planning, appropriate selection of the basal ganglia motor plan to hunt, and motor output of the primary motor cortex. In contrast to the lateral SC lesions, medial SC lesions produced a small deficit in predatory hunting, and compared to what we have seen for the lateral SCig, the medial SCig has a very limited set of projections to thalamic sites related to the control of motor planning or motor output, and provides conspicuous inputs to brainstem sites involved in organizing a wide range of anti-predatory defensive responses. Overall, the present results served to clarify how the different functional domains in the SC may mediate the decision to pursue and hunt a prey or escape from a predator. (C) 2010 IBRO. Published by Elsevier Ltd. All rights reserved.

INDEFINITE QUASILINEAR ELLIPTIC EQUATIONS IN EXTERIOR DOMAINS WITH EXPONENTIAL CRITICAL GROWTH

This paper is concerned with the existence of solutions for the quasilinear problem {-div(vertical bar del u vertical bar(N-2) del u) + vertical bar u vertical bar(N-2) u = a(x)g(u) in Omega u = 0 on partial derivative Omega, where Omega subset of R(N) (N >= 2) is an exterior domain; that is, Omega = R(N)\omega, where omega subset of R(N) is a bounded domain, the nonlinearity g(u) has an exponential critical growth at infinity and a(x) is a continuous function and changes sign in Omega. A variational method is applied to establish the existence of a nontrivial solution for the above problem.

Semilinear parabolic problems in thin domains with a highly oscillatory boundary

In this paper, we study the behavior of the solutions of nonlinear parabolic problems posed in a domain that degenerates into a line segment (thin domain) which has an oscillating boundary. We combine methods from linear homogenization theory for reticulated structures and from the theory on nonlinear dynamics of dissipative systems to obtain the limit problem for the elliptic and parabolic problems and analyze the convergence properties of the solutions and attractors of the evolutionary equations. (C) 2011 Elsevier Ltd. All rights reserved.

Dynamics in dumbbell domains II. The limiting problem

In this work we continue the analysis of the asymptotic dynamics of reaction-diffusion problems in a dumbbell domain started in [J.M. Arrieta, AN Carvalho, G. Lozada-Cruz, Dynamics in dumbbell domains I. Continuity of the set of equilibria, J. Differential Equations 231 (2) (2006) 551-597]. Here we study the limiting problem, that is, an evolution problem in a ""domain"" which consists of an open, bounded and smooth set Omega subset of R(N) with a curve R(0) attached to it. The evolution in both parts of the domain is governed by a parabolic equation. In Omega the evolution is independent of the evolution in R(0) whereas in R(0) the evolution depends on the evolution in Omega through the continuity condition of the solution at the junction points. We analyze in detail the linear elliptic and parabolic problem, the generation of linear and nonlinear semigroups, the existence and structure of attractors. (C) 2009 Elsevier Inc. All rights reserved.

Dynamics in dumbbell domains III. Continuity of attractors

In this paper we conclude the analysis started in [J.M. Arrieta, AN Carvalho, G. Lozada-Cruz, Dynamics in dumbbell domains I. Continuity of the set of equilibria, J. Differential Equations 231 (2006) 551-597] and continued in [J.M. Arrieta, AN Carvalho, G. Lozada-Cruz, Dynamics in dumbbell domains II. The limiting problem, J. Differential Equations 247 (1) (2009) 174-202 (this issue)] concerning the behavior of the asymptotic dynamics of a dissipative reaction-diffusion equation in a dumbbell domain as the channel shrinks to a line segment. In [J.M. Arrieta, AN Carvalho. G. Lozada-Cruz, Dynamics in dumbbell domains I. Continuity of the set of equilibria, J. Differential Equations 231 (2006) 551-597], we have established an appropriate functional analytic framework to address this problem and we have shown the continuity of the set of equilibria. In [J.M. Arrieta, AN Carvalho, G. Lozada-Cruz. Dynamics in dumbbell domains II. The limiting problem, J. Differential Equations 247 (1) (2009) 174-202 (this issue)], we have analyzed the behavior of the limiting problem. In this paper we show that the attractors are Upper semicontinuous and, moreover, if all equilibria of the limiting problem are hyperbolic, then they are lower semicontinuous and therefore, continuous. The continuity is obtained in L(p) and H(1) norms. (C) 2008 Elsevier Inc. All rights reserved.

Infeasibility handling in genetic algorithm using nested domains for production planning

In this paper we present a genetic algorithm with new components to tackle capacitated lot sizing and scheduling problems with sequence dependent setups that appear in a wide range of industries, from soft drink bottling to food manufacturing. Finding a feasible solution to highly constrained problems is often a very difficult task. Various strategies have been applied to deal with infeasible solutions throughout the search. We propose a new scheme of classifying individuals based on nested domains to determine the solutions according to the level of infeasibility, which in our case represents bands of additional production hours (overtime). Within each band, individuals are just differentiated by their fitness function. As iterations are conducted, the widths of the bands are dynamically adjusted to improve the convergence of the individuals into the feasible domain. The numerical experiments on highly capacitated instances show the effectiveness of this computational tractable approach to guide the search toward the feasible domain. Our approach outperforms other state-of-the-art approaches and commercial solvers. (C) 2009 Elsevier Ltd. All rights reserved.

Front tracking with moving-least-squares surfaces

The representation of interfaces by means of the algebraic moving-least-squares (AMLS) technique is addressed. This technique, in which the interface is represented by an unconnected set of points, is interesting for evolving fluid interfaces since there is]to surface connectivity. The position of the surface points can thus be updated without concerns about the quality of any surface triangulation. We introduce a novel AMLS technique especially designed for evolving-interfaces applications that we denote RAMLS (for Robust AMLS). The main advantages with respect to previous AMLS techniques are: increased robustness, computational efficiency, and being free of user-tuned parameters. Further, we propose a new front-tracking method based on the Lagrangian advection of the unconnected point set that defines the RAMLS surface. We assume that a background Eulerian grid is defined with some grid spacing h. The advection of the point set makes the surface evolve in time. The point cloud can be regenerated at any time (in particular, we regenerate it each time step) by intersecting the gridlines with the evolved surface, which guarantees that the density of points on the surface is always well balanced. The intersection algorithm is essentially a ray-tracing algorithm, well-studied in computer graphics, in which a line (ray) is traced so as to detect all intersections with a surface. Also, the tracing of each gridline is independent and can thus be performed in parallel. Several tests are reported assessing first the accuracy of the proposed RAMLS technique, and then of the front-tracking method based on it. Comparison with previous Eulerian, Lagrangian and hybrid techniques encourage further development of the proposed method for fluid mechanics applications. (C) 2008 Elsevier Inc. All rights reserved.

An implicit technique for solving 3D low Reynolds number moving free surface flows

This paper describes the development of an implicit finite difference method for solving transient three-dimensional incompressible free surface flows. To reduce the CPU time of explicit low-Reynolds number calculations, we have combined a projection method with an implicit technique for treating the pressure on the free surface. The projection method is employed to uncouple the velocity and the pressure fields, allowing each variable to be solved separately. We employ the normal stress condition on the free surface to derive an implicit technique for calculating the pressure at the free surface. Numerical results demonstrate that this modification is essential for the construction of methods that are more stable than those provided by discretizing the free surface explicitly. In addition, we show that the proposed method can be applied to viscoelastic fluids. Numerical results include the simulation of jet buckling and extrudate swell for Reynolds numbers in the range [0.01, 0.5]. (C) 2008 Elsevier Inc. 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.

Electrical transport in disordered and ordered magnetic domains under pressures and magnetic fields

Magnetic M( T, H, P) and electrical transport.( T, H, P) measurements in a strong spin-lattice-charge coupled La(0.7)Ca(0.3)MnO(3) system have been conducted. The application of H and P leads to the formation of different magnetic domain structures in the vicinity and below the polaronic-to-ferromagnetic transition temperature. The charge mobility is more sensitive to the variation of the spatial wave function overlap between Mn(3+) eg and O(2-) 2p orbitals due to the applied compacting pressure rather than the relative spin orientation between neighbouring Mn ions when the magnetic field is applied. In spite of the presence of different magnetic domain structures due to the sample history, the effect of magnetic field and pressure is less pronounced at lower temperatures on electrical transport properties.

Vesicles with charged domains

This work summarizes results obtained on membranes composed of the ternary mixture dioleoylphosphatidylglycerol (DOPG), egg sphingomyelin (eSM) and cholesterol (Chol). The membrane phase state as a function of composition is characterized from data collected with fluorescence microscopy on giant unilamellar vesicles. The results suggest that the presence of the charged DOPG significantly decreases the composition region of coexistence of liquid ordered and liquid disordered phases as compared to that in the ternary mixture of dioleoylphosphatidycholine, sphingomyelin and cholesterol. The addition of calcium chloride to DOPG:eSM:Chol vesicles, and to a lesser extent the addition of sodium chloride, leads to the stabilization of the two-phase coexistence region, which is expressed in an increase in the miscibility temperature. On the other hand, addition of the chelating agent EDTA has the opposite effect, suggesting that impurities of divalent cations in preparations of giant vesicles contribute to the stabilization of charged domains. We also explore the behavior of these membranes in the presence of extruded unilamellar vesicles made of the positively charged lipid dioleoyltrimethylammoniumpropane (DOTAP). The latter can induce domain formation in DOPG:eSM:Chol vesicles with initial composition in the one-phase region. (C) 2010 Elsevier B.V. All rights reserved.

Reaction microtextures of monazite: Correlation between chemical and age domains in the Nazare Paulista migmatite, SE Brazil

Back-scattered imaging, X-ray element mapping and electron microprobe analyzer (EMPA) chemical dating reveal complex compositional and age zoning in monazite crystals from different layers and textural positions in a garnet-bearing migmatite in SE Brazil. Y-rich (variable Y(2)O(3), averaging 2.5 wt.%) relict cores are preserved in mesosome and melanosome monazite, and correspond to 793 +/- 6 Ma inherited crystals possibly generated in a previous metamorphic event. These cores are overgrown and widely replaced by two generations of monazite, which are present in all migmatite layers. The first, also Y-rich (average 2.5 wt.% Y(2)O(3)), was produced at similar to 635 Ma during prograde metamorphism under subsolidus conditions, while the second has an Y-poor (<1.5 wt.% Y(2)O(3)), low Th/U signature, and precipitated from low Y and HREE anatectic melts produced by reactions in which garnet was inert. Quartz-rich trondhjemitic leucosome represents lower temperature melt (bearing some subsolidus quartz and garnet with included monazite) formed at temperatures below muscovite breakdown; its Y-poor monazite indicates an age of 617 +/- 6 Ma. Granitic leucosomes formed close to peak metamorphic conditions (T>750 degrees C) above muscovite breakdown have their slightly younger character confirmed by a 609 +/- 7 Ma low-Y monazite age. A similar 606 +/- 5 Ma age was obtained for low-Y monazite rims and domains in mesosome and melanosome, and reflects the time of monazite saturation in interstitial granitic melt that was trapped in these layers. Our results confirm that inherited monazite crystals can be preserved during partial melting at temperatures above muscovite breakdown. Moreover, careful textural control aided by X-ray chemical mapping may allow monazite generated at different stages in a similar to 25 Myr prograde metamorphic path to be identified and dated using an electron microprobe. (C) 2008 Elsevier B.V. All rights reserved.

U-Pb-Hf-trace element systematics and geochronology of zircon from a granulite-facies metamorphosed mafic-ultramafic layered complex in Central Brazil

Zircon recrystallization is a common process during high-grade metamorphism and promotes partial or complete resetting of the original isotopic and chemical characteristics of the mineral and thus complicates U-Pb geochronological interpretation. In Central Brazil, this complexity may be illustrated by three composite mafic-ultramafic intrusions metamorphosed under amphibolite-to-granulite conditions. Their ages of emplacement and metamorphic ages have been a matter of controversy for the last thirty years. The Serra da Malacacheta and Barro Alto complexes make up the southernmost of these layered bodies and four samples from distinct rock types were investigated in order to verify the consequences of metamorphic alteration of zircon for U-Pb dating. Cathodoluminescent imaging reveals internal features which are typical of concomitant dissolution-reprecipitation processes, such as convolute zoning and inward-moving recrystallization fronts, even in samples in which partially preserved igneous textures are observed. Due to this extensive alteration, LA-ICPMS U-Pb isotopic analysis yielded inconclusive data. However, in situ Hf isotopic and trace-element analyses help to clarify the real meaning of the geochronological data. Low Lu/Hf (<0.004) and homogeneous (176)Hf/(177)Hf(t) values imply that the zircon populations within individual samples have crystallized in a single episode, despite the observed variations in age values. Trace element signatures of zircon grains from garnet-bearing samples reveal that they were unreactive to the development of the peak metamorphism mineral assemblage and, thus, the main chemical feature in such grains is attributed to a coupled dissolution-reprecipitation process. However, in the Cafelandia amphibolite an additional alteration process is identified, probably related to the influx of late-stage fluids. Combined isotopic and geochemical investigation on zircon grains allowed the distinction of two magmatic events. The first corresponds to the crystallization of the Serra da Malacacheta Complex and characterizes a juvenile magmatism at similar to 1.3 Ga. The younger episode, recognized in the Barro Alto Complex, is dated at ca. 800 Ma and is represented by mafic and ultramafic rocks showing intense contamination with continental crust, implying that the emplacement took place, most likely, in a continental back-arc setting. Altered zircon domains as well as titanite grains date the metamorphic event at ca. 760-750 Ma. (C) 2011 Elsevier B.V. All rights reserved.

Dating coeval mafic magmatism and ultrahigh temperature metamorphism in the Anapolis-Itaucu Complex, Central Brazil

Dating granulites has always been of great interest because they represent one of the most extreme settings of an orogen. Owing to the resilience of zircon, even in such severe environments, the link between P-T conditions and geological time is possible. However, a challenge to geochronologists is to define whether the growth of new zircon is related to pre- or post-P-T peak conditions and which processes might affect the (re) crystallization. In this context, the Anapolis-Itaucu Complex, a high-grade complex in central Brazil with ultrahigh temperature (UHT) granulites, may provide valuable information within this topic. The Anapolis-Itaucu Complex (AIC) includes ortho- and paragranulites, locally presenting UHT mineral assemblages, with igneous zircon ages varying between 760 and 650 Ma and metamorphic overgrowths dated at around 650-640 Ma. Also common in the Anapolis-Itaucu Complex are layered mafic-ultramafic complexes metamorphosed under high-grade conditions. This article presents the first geological and geochronological constraints of three of these layered complexes within the AIC, the Damolandia, Taquaral and Goianira-Trindade complexes. U-Pb (LA-MC-ICPMS, SHRIMP and ID-TIMS) zircon analyses reveal a spread of concordant ages spanning within an age interval of similar to 80 Ma with an ""upper"" intercept age of similar to 670 Ma. Under cathodoluminescence imaging, these crystals show partially preserved primary sector zoning, as well as internal textures typical of alteration during high-grade metamorphism, such as inward-moving boundaries. Zircon grains reveal homogeneous initial (176)Hf/(177)Hf values in distinct crystal-scale domains in all samples. Moreover. Hf isotopic ratios show correlation neither with U-Pb ages nor with Th/U ratios, suggesting that zircon grains crystallized during a single growth event. It is suggested, therefore, that the observed spread of concordant U-Pb ages may be related to a memory effect due to coupled dissolution-reprecipitation process during high grade metamorphism. Therefore, understanding the emplacement and metamorphism of this voluminous mafic magmatism is crucial as it may represent an additional heat source for the development of the ultrahigh temperature paragenesis recorded in the paragranulites. (C) 2010 Elsevier B.V. All rights reserved.