Numerical approximation of the Ginzburg-Landau equation with memory effects in the dynamics of phase transitions

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

A generalized alternating-direction implicit scheme for incompressible magnetohydrodynamic viscous flows at low magnetic Reynolds number

This paper presents numerical simulations of incompressible fluid flows in the presence of a magnetic field at low magnetic Reynolds number. The equations governing the flow are the Navier-Stokes equations of fluid motion coupled with Maxwell's equations of electromagnetics. The study of fluid flows under the influence of a magnetic field and with no free electric charges or electric fields is known as magnetohydrodynamics. The magnetohydrodynamics approximation is considered for the formulation of the non-dimensional problem and for the characterization of similarity parameters. A finite-difference technique is used to discretize the equations. In particular, an extension of the generalized Peaceman and Rachford alternating-direction implicit (ADI) scheme for simulating two-dimensional fluid flows is presented. The discretized conservation equations are solved in stream function-vorticity formulation. We compare the ADI and generalized ADI schemes, and show that the latter is more efficient in simulating low Reynolds number and magnetic Reynolds number problems. Numerical results demonstrating the applicability of this technique are also presented. The simulation of incompressible magneto hydrodynamic fluid flows is illustrated by numerical solution for two-dimensional cases. (c) 2007 Elsevier B.V. All rights reserved.

Interpolating EFGM for computing continuous and discontinuous electromagnetic fields

Purpose - This paper proposes an interpolating approach of the element-free Galerkin method (EFGM) coupled with a modified truncation scheme for solving Poisson's boundary value problems in domains involving material non-homogeneities. The suitability and efficiency of the proposed implementation are evaluated for a given set of test cases of electrostatic field in domains involving different material interfaces.Design/methodology/approach - the authors combined an interpolating approximation with a modified domain truncation scheme, which avoids additional techniques for enforcing the Dirichlet boundary conditions and for dealing with material interfaces usually employed in meshfree formulations.Findings - the local electric potential and field distributions were correctly described as well as the global quantities like the total potency and resistance. Since, the treatment of the material interfaces becomes practically the same for both the finite element method (FEM) and the proposed EFGM, FEM-oriented programs can, thus, be easily extended to provide EFGM approximations.Research limitations/implications - the robustness of the proposed formulation became evident from the error analyses of the local and global variables, including in the case of high-material discontinuity.Practical implications - the proposed approach has shown to be as robust as linear FEM. Thus, it becomes an attractive alternative, also because it avoids the use of additional techniques to deal with boundary/interface conditions commonly employed in meshfree formulations.Originality/value - This paper reintroduces the domain truncation in the EFGM context, but by using a set of interpolating shape functions the authors avoided the use of Lagrange multipliers as well Mathematics in Engineering high-material discontinuity.

Criteria for mixed grids in computational fluid dynamics

The finite volume method is used as a numerical method for solving the fluid flow equations. This method is appropriate to employ under structured and unstructured meshes. Mixed grids, combining both types of grids, are investigated. The coupling of different grids is done by overlapping strategy. The computational effort for the mixed grid is evaluated by the CPU-time, with different percentage of covering area of the unstructured mesh. The present scheme is tested for the driven cavity problem, where the incompressible fluid is integrated by calculating the velocity fields and computing the pressure field in each time step. Several schemes for unstructured grid are examined, and the compatibility condition is applied to check their consistency. A scheme to verify the compatibility condition for the unstructured grids is presented. (c) 2006 IMACS. Published by Elsevier B.V. All rights reserved.

Stress analysis of an upper central incisor restored with different posts

The effect of different anatomic shapes and materials of posts in the stress distribution on an endodontically treated incisor was evaluated in this work. This study compared three post shapes (tapered, cylindrical and two-stage cylindrical) made of three different materials (stainless steel, titanium and carbon fibre on Bisphenol A-Glycidyl Methacrylate (Bis-GMA) matrix). Two-dimensional stress analysis was performed using the Finite Element Method. A static load of 100N was applied at 45degrees inclination with respect to the incisor's edge. The stress concentrations did not significantly affect the region adjacent to the alveolar bone crest at the palatine portion of the tooth, regardless of the post shape or material. However, stress concentrations on the post/dentin interface on the palatine side of the tooth root presented significant variations for different post shapes and materials. Post shapes had relatively small impact on the stress concentrations while post materials introduced higher variations on them. Stainless steel posts presented the highest level of stress concentration, followed by titanium and carbon/Bis-GMA posts.

Automation and cybernetics: Control of a flexible one-link manipulator

States that control is of the essence in cybernetics. Summarizes the dynamic equations for a flexible one-link manipulator moving in the horizontal plane. Employs the finite element method, based on elementary beam theory, during the process of formulation. Develops and instruments a one-link flexible manipulator in order to control its vibration modes. Uses a simple second-order vibration model which permits vibrations on the rod to be estimated using the hub angle. The validation of the dynamic model and the structural analysis of the flexible manipulator is reached using proper infrared cameras and active light sources for determining actual positions of objects in space. Shows that the performance of the control is satisfactory, even under perturbation action.

Dynamic modeling and simulation of a flexible robotic manipulator

This work focuses on the dynamic modeling of a flexible robotic manipulator with two flexible links and two revolute joints, which rotates in the horizontal plane. The dynamic equations are derived using the Newton-Euler formulation and the finite element method, based on elementary beam theory. Computer simulation results are presented to illustrate this study. The dynamic model becomes necessary for use in future design and control applications.

Nutrient digestibility of degermed, dehulled corn, citrus pulp, and soy protein concentrate by barrows

This study was carried out to determine apparent ileal digestibility (AID) and apparent total tract digestibility (ATTD) of DM, CP, GE, and their respective digestible content of degermed dehulled corn (Zea mays), citrus pulp, and soy (Glycine max) protein concentrate by pigs using the difference method. Thirty-two barrows (28.1 +/- 1.6 kg of BW) were fed a corn-soybean meal basal diet or 1 of 3 diets formulated by replacing 30% of the basal diet with 30% of 1 of the test feedstuffs for 11 d. Chromic oxide (0.3%) was included in the diets. Feces were collected from days 7 to 11 by grab sampling and ileal digesta were collected after pigs were slaughtered on day 12. The AID of DM and AID and ATTD of GE of degermed corn (77.4, 88.7, and 77.7%) were greater (P < 0.05) than those observed in citrus pulp (50.3, 86.5, and 55.8%) and in soy protein concentrate (63.5, 85.1, and 59.4%), which did not differ (P > 0.05). The ATTD of CP, total digestible CP, and total DE of soy protein concentrate (87.5%, 500 g/kg, and 3739 kcal/kg) were higher (P < 0.05) than the values in degermed corn (81.7%, 57.5 g/kg, and 3330 kcal/kg), which were greater (P < 0.05) than those in citrus pulp (60.5%, 39.5 g/kg, and 3223 kcal/kg). Total and ileal digestible DM, AID of CP, and ileal DE of degermed corn (782 g/kg, 673 g/kg, 70.7%, and 2913 kcal/kg) and soy protein concentrate (778 g/kg, 570 g/kg, 78.7%, and 2878 kcal/kg) were similar (P > 0.05) and greater (P < 0.05) than those in citrus pulp (737 g/kg, 436 g/kg, 50.6%, and 2081 kcal/kg). Ileal digestible CP of degermed corn (49.8 g/kg) and citrus pulp (33.0 g/kg) did not differ (P > 0.05) but were smaller (P < 0.05) than the value found in soy protein concentrate (434 g/kg). The DM and energy from degermed corn are more efficiently digested by the pig than those from soy protein concentrate and citrus pulp. Soy protein concentrate was the best protein source evaluated in this study.

Application of a model updating technique for health condition monitoring and damage detection of flexible structures

This paper discusses the application of a damage detection methodology to monitor the location and extent of partial structural damage. The methodology combines, in an iterative way, the model updating technique based on frequency response functions (FRF) with monitoring data aiming at identifying the damage area of the structure. After the updating procedure reaches a good correlation between the models, it compares the parameters of the damage structure with those of the undamaged one to find the deteriorated area. The influence of the FEM mesh size on the evaluation of the extent of the damage has also been discussed. The methodology is applied using real experimental data from a spatial frame structure.

PCM thermal insulation in buildings

This paper presents the results of a numerical and experimental study of phase change material (PCM) filled walls and roofs under real operational conditions to achieve passive thermal comfort. The numerical part of the study was based on a one-dimensional model for the phase change problem controlled by pure conduction. Real radiation data was used to determine the external face temperature. The numerical treatment was based upon using finite difference approximations and the ADI scheme. The results obtained were compared with field measurements. The experimental set-up consisted of a small room with movable roof and side wall. The roof was constructed in the traditional way but with the phase change material enclosed. Thermocouples were distributed across the cross section of the roof. Another roof, identical but without the PCM, was also used during comparative tests. The movable wall was also constructed as is done traditionally but with the PCM enclosed. Again, thermocouples were distributed across the wall thickness to enable measurement of the local temperatures. Another wall, identical but without the PCM, was also used during comparative tests. The PCM used in the numerical and experimental tests was composed of a mixture of two commercial grades of glycol in order to obtain the required fusion temperature range. Comparison between the simulation results and the experiments indicated good agreement. Field tests also indicated that the PCM used was adequate and that the concept was effective in maintaining the indoor temperature very close to the established comfort limits. Further economical analysis indicated that the concept could effectively help in reducing the electric energy consumption and improving the energy demand pattern. © 1997 by John Wiley & Sons, Ltd.

Numerical and experimental analysis of natural convection in a cavity with flush mounted heat sources on a side wall

This work considers a problem of interest in several technological applications such as the thermal control of electronic equipment. It is also important to study the heat transfer performance of these components under off-normal conditions, such as during failure of cooling fans. The effect of natural convection on the flow and heat transfer in a cavity with two flush mounted heat sources on the left vertical wall, simulating electronic components, is studied numerically and experimentally. The influence of the power distribution, spacing between the heat sources and cavity aspect ratio have been investigated. An analysis of the average Nusselt number of the two heat sources was performed to investigate the behavior of the heat transfer coefficients. The results obtained numerically and experimentally, after an error analysis, showed a good agreement.

Simulação numérica de trocadores de calor multitubulares equipados com fitas helicoidais externamente aos tubos

Classical shell-and-tube heat exchangers are usually equipped with segmental baffles. These baffles serve two basic functions: (a) they provide tube supports, thereby preventing or reducing mechanical problems, such as sagging or vibration; (b) they direct the fluid flow over the tubes so as to introduce a cross-flow component, thereby increasing the heat transfer. Segmented baffles have several sources of performance loss, some due to various leakage flows and others caused by stagnation zones. A new concept of longitudinal flow heat exchanger - based on placing twisted tapes along the tube bundle subchannels - was developed to mitigate drawbacks of other types of tubular heat exchangers. In this paper, a numerical model has been implemented in order to simulate the thermal-hydraulic feature of tubular heat exchangers equipped either with segmental baffles or with subchannel twisted tapes. The tube bundle has been described by means of an equivalent porous medium type model, allowing a macroscopic description of the shell-side flow. The basic equations - continuity, momentum and energy - have been solved by using the finite volume method. Typical numerical results have been compared with experimental data, reaching a very good agreement. A comparative analysis of different types of heat exchangers has been carried out, revealing the satisfactory thermal-hydraulic efficiency level of the twisted tapes heat exchangers.

Study on transient ice formation of laminar flow inside externally supercooled rectangular duct

The transient process of solidification of laminar liquid flow (water) submitted to super-cooling was investigated both theoretically and experimentally. In this study an alternative analytical formulation and numerical approach were adopted resulting in the unsteady model with temperature dependent thermophysical properties in the solid region. The proposed model is based upon the fundamental equations of energy balance in the solid and liquid regions as well as across the solidification front. The basic equations and the associated boundary and initial conditions were made dimensionless by using the Landau transformation to immobilize the moving front and render the problem to a fixed plane type problem. A laminar velocity profile is admitted in the liquid domain and the resulting equations were discretized using the finite difference approach. The numerical predictions obtained were compared with the available results based on other models and concepts such as Neumann analytical model, the apparent thermal capacity model due to Bonacina and the conventional fixed grid energy model due to Goodrich. To obtain further comparisons and more validation of the model and the numerical solution, an experimental rig was constructed and instrumented permitting very well controlled experimental measurements. The numerical predictions were compared with the experimental results and the agreement was found satisfactory.

Optimal design of smart structures using bonded piezoelectrics for vibration control

Smart material technology has become an area of increasing interest for the development of lighter and stronger structures which are able to incorporate actuator and sensor capabilities for collocated control. In the design of actively controlled structures, the determination of the actuator locations and the controller gains, is a very important issue. For that purpose, smart material modelling, modal analysis methods, control and optimization techniques are the most important ingredients to be taken into account. The optimization problem to be solved in this context presents two interdependent aspects. The first one is related to the discrete optimal actuator location selection problem, which is solved in this paper using genetic algorithms. The second is represented by a continuous variable optimization problem, through which the control gains are determined using classical techniques. A cantilever Euler-Bernoulli beam is used to illustrate the presented methodology.

A FEM approach to the equations of magneto-aerodynamics

In this work, a Finite Element Method treatment is outlined for the equations of Magnetoaerodynamics. In order to provide a good basis for numerical treatment of Magneto-aerodynamics, a full version of the complete equations is presented and FEM contribution matrices are deduced, as well as further terms of stabilization for the compressible flow case.