Numerical simulation of viscous three-dimensional flow over aerospace vehicles using Euler/boundary-layer zonal approach

This paper presents a viscous three-dimensional simulations coupling Euler and boundary layer codes for calculating flows over arbitrary surfaces. The governing equations are written in a general non orthogonal coordinate system. The Levy-Lees transformation generalized to three-dimensional flows is utilized. The inviscid properties are obtained from the Euler equations using the Beam and Warming implicit approximate factorization scheme. The resulting equations are discretized and approximated by a two-point fmitedifference numerical scheme. The code developed is validated and applied to the simulation of the flowfield over aerospace vehicle configurations. The results present good correlation with the available data.

Optimization procedures in the design of continuous induction hardening and tempering installations for magnetic steel bars

An inverse problem concerning the industrial process of steel bars hardening and tempering is considered. The associated optimization problem is formulated in terms of membership functions and, for the sake of comparison, also in terms of quadratic residuals; both geometric and electromagnetic design variables have been considered. The numerical solution is achieved by coupling a finite difference procedure for the calculation of the electromagnetic and thermal fields to a deterministic strategy of minimization based on modified Flctcher and Reeves method. © 1998 IEEE.

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.

Pressure field in the cathode-anode region of a high-power klystron amplifier

The pressure field of a high-power klystron amplifier in the cathode and anode region was investigated. The investigation was performed using a 1.3 GHz, 100 A and 240 kV high-power klystron with five reentrant coaxial cavities, assembled in cylindrical drift tube 1.2 m long. The diffusion equation in mathematical model was also solved by using a 3-D finite element method code, in order to obtain pressure profile in region of interest. The results show that density profile of molecules between cathode-anode region was determined, where cathode pressure is approximately 10% higher than anode pressure.

A 2-D delaunay refinement algorithm using an initial prerefinement from the boundary mesh

An improvement to the quality bidimensional Delaunay mesh generation algorithm, which combines the mesh refinement algorithms strategy of Ruppert and Shewchuk is proposed in this research. The developed technique uses diametral lenses criterion, introduced by L. P. Chew, with the purpose of eliminating the extremely obtuse triangles in the boundary mesh. This method splits the boundary segment and obtains an initial prerefinement, and thus reducing the number of necessary iterations to generate a high quality sequential triangulation. Moreover, it decreases the intensity of the communication and synchronization between subdomains in parallel mesh refinement. © 2008 IEEE.

Análise da execução de um algoritmo numérico cbs para simulação de escoamento utilizando elementos finitos

This work intends to analyze the application and execution time of a numerical algorithm that simulates incompressible and isothermal flows. It was used the explicit scheme of the Characteristic Based Split (CBS) algorithm and the Artificial Compressibility (AC) scheme for coupling pressure-velocity equations. The discretization was done with the finite element method using a bilinear elements grid. The free software GNU-Octave was used for implementation and execution of routines. The results were analyzed using the classic lid-driven cavity problem. This work shows results for tests with several Reynolds' number. The results for these tests show a good agreement when compared with previous ones obtained from bibliography. The code runtime's analysis shows yet that the matrix's assembly is the part of greater consumption time in the implementation.

Coupling rate measurement of a novel multi-actuated piezoelectric device using optical interferometry

The class of piezoelectric actuators considered in this paper consists of a multi-flexible structure actuated by two or more piezoceramic devices that must generate different output displacements and forces at different specified points of the domain and in different directions. The devices were modeled by finite element using the software ANSYS and the topology optimization method. The following XY actuators were build to achieve maximum displacement in the X and Y directions with a minimum crosstalk between them. The actuator prototypes are composed of an aluminum structure, manufactured by using a wire Electrical Discharge Machining, which are bonded to rectangular PZT5A piezoceramic blocks by using epoxy resin. Multi-actuator piezoelectric device displacements can be measured by using optical interferometry, since it allows dynamic measurements in the kHz range, which is of the order of the first resonance frequency of these piezomechanisms. A Michelson-type interferometer, with a He-Ne laser source, is used to measure the displacement amplitudes in nanometric range. A new optical phase demodulation technique is applied, based on the properties of the triangular waveform drive voltage applied to the XY piezoelectric nanopositioner. This is a low-phase-modulation-depth-like technique that allows the rapid interferometer auto-calibration. The measurements were performed at 100 Hz frequency, and revealed that the device is linear voltage range utilized in this work. The ratio between the generated and coupled output displacements and the drive voltages is equal to 10.97 nm/V and 1.76 nm/V, respectively, which corresponds to a 16% coupling rate. © 2010 IEEE.

Influence of high insertion torque on implant placement - an anisotropic bone stress analysis

The aim of this study was to evaluate the influence of the high values of insertion torques on the stress and strain distribution in cortical and cancellous bones. Based on tomography imaging, a representative mathematical model of a partial maxilla was built using Mimics 11.11 and Solid Works 2010 softwares. Six models were built and each of them received an implant with one of the following insertion torques: 30, 40, 50, 60, 70 or 80 Ncm on the external hexagon. The cortical and cancellous bones were considered anisotropic. The bone/implant interface was considered perfectly bonded. The numerical analysis was carried out using Ansys Workbench 10.0. The convergence of analysis (6%) drove the mesh refinement. Maximum principal stress (σ max) and maximum principal strain (ε max) were obtained for cortical and cancellous bones around to implant. Pearson's correlation test was used to determine the correlation between insertion torque and stress concentration in the periimplant bone tissue, considering the significance level at 5%. The increase in the insertion torque generated an increase in the σ max and ε max values for cortical and cancellous bone. The σmax was smaller for the cancellous bone, with greater stress variation among the insertion torques. The ε max was higher in the cancellous bone in comparison to the cortical bone. According to the methodology used and the limits of this study, it can be concluded that higher insertion torques increased tensile and compressive stress concentrations in the periimplant bone tissue.

Dependência estrutural da capacidade de rotação plástic em vigas de concreto armado

Summary In this work the structural dependence of plastic rotation capacity in RC beams is evaluated using the Finite Element Method. The objective is to achieve a better understanding of the non-linear behavior of reinforced concrete members and perform extensive parameter studies, using a rational model developed by Bigaj [1] to analyze the phenomenon of plastic rotation capacity in reinforced concrete members. It is assumed that only bending failure is relevant due to sufficient member resistance against shear and torsion. The paper begins with the physical and theoretical background of the phenomenon of plastic hinge development in RC structures. Special emphasis is laid on the issue of structural dependence of deformation capacity of plastic hinges in RC members. Member size dependence and influence of properties of construction materials were emphasized as well. The essential components of the Bigajs model for calculating the plastic rotation capacity are discussed. The behaviour of the plastic hinge is analysed taking into account the strain localisation in the damage zones of the hinge region. The Fictitious Crack Model (FCM) and the Compressive Damage Zone Model (CDZ) are adopted in a Fracture Mechanics approach to model the behaviour of concrete in tension and compression, respectively. The approach is implemented in FEMOOP, a FEM in-house solver under development, and applied to evaluate ductility in 2D beams. The models were generated with GiD, a pre-processor and post-processor developed by CIMNE, and analyzed with the capabilities implemented in FEMOOP. © Universitat Politècnica de Catalunya, Barcelona, España 2010.

Integrating open source codes for solid modeling and automatic mesh generation in FEM applications

The applications of the Finite Element Method (FEM) for three-dimensional domains are already well documented in the framework of Computational Electromagnetics. However, despite the power and reliability of this technique for solving partial differential equations, there are only a few examples of open source codes available and dedicated to the solid modeling and automatic constrained tetrahedralization, which are the most time consuming steps in a typical three-dimensional FEM simulation. Besides, these open source codes are usually developed separately by distinct software teams, and even under conflicting specifications. In this paper, we describe an experiment of open source code integration for solid modeling and automatic mesh generation. The integration strategy and techniques are discussed, and examples and performance results are given, specially for complicated and irregular volumes which are not simply connected. © 2011 IEEE.

Fractal dimension and Shannon's entropy analyses of the architectural complexity caused by the inflammatory reactions induced by highly crystalline polyvinyl alcohol) microspheres implanted in subcutaneous tissues of the Wistar rats

The results of the histopathological analyses after the implantation of highly crystalline PVA microspheres in subcutaneous tissues of Wistar rats are here in reported. Three different groups of PVA microparticles were systematically studied: highly crystalline, amorphous, and commercial ones. In addition to these experiments, complementary analyses of architectural complexity were performed using fractal dimension (FD), and Shannon's entropy (SE) concepts. The highly crystalline microspheres induced inflammatory reactions similar to the ones observed for the commercial ones, while the inflammatory reactions caused by the amorphous ones were less intense. Statistical analyses of the subcutaneous tissues of Wistar rats implanted with the highly crystalline microspheres resulted in FD and SE values significantly higher than the statistical parameters observed for the amorphous ones. The FD and SE parameters obtained for the subcutaneous tissues of Wistar rats implanted with crystalline and commercial microparticles were statistically similar. Briefly, the results indicated that the new highly crystalline microspheres had biocompatible behavior comparable to the commercial ones. In addition, statistical tools such as FD and SE analyses when combined with histopathological analyses can be useful tools to investigate the architectural complexity tissues caused by complex inflammatory reactions. © 2012 WILEY PERIODICALS, INC.

Influence of convergence angle of tooth preparation on the fracture resistance of Y-TZP-based all-ceramic restorations

Objective: To investigate the influence of the convergence angle of tooth preparation on the fracture load of Y-TZP-based ceramic (YZ-Vita YZ) substructure (SB) veneered with a feldspathic porcelain (VM9-Vita VM9). Methods: Finite element stress analysis (FEA) was performed to examine the stress distribution of the system. Eighty YZ SB were fabricated using a CAD-CAM system and divided into four groups (n = 20), according to the total occlusal convergence (TOC) angle: G6-6° TOC; G12-12° TOC; G20-20° TOC; and G20MOD-20° TOC with modified SB. All SB were veneered with VM9, cemented in a fiber reinforced epoxy resin die, and loaded to failure. Half of the specimens from each group (n = 10) were cyclic fatigued (106 cycles) before testing. Failure analysis was performed to determine the fracture origin. Data were statistically analyzed using Anova and Tukey's tests (α = 0.05). Results: The greatest mean load to fracture value was found for the G20MOD, which was predicted by the FEA. Cyclic fatigue did not significantly affect the load of fracture. Catastrophic failure originating from the internal occlusal surface of the SB was the predominant failure mode, except for G20MOD. Significance: The YZ-VM9 restorations resisted greater compression load than the usual physiological occlusal load, regardless of the TOC angle of preparations. Yet, the G20MOD design produced the best performance among the experimental conditions evaluated. © 2013 Academy of Dental Materials. Published by Elsevier Ltd. All rights reserved.

Three-dimensional modeling and highly refined mesh generation of the aorta artery and its tunics

This paper describes strategies and techniques to perform modeling and automatic mesh generation of the aorta artery and its tunics (adventitia, media and intima walls), using open source codes. The models were constructed in the Blender package and Python scripts were used to export the data necessary for the mesh generation in TetGen. The strategies proposed are able to provide meshes of complicated and irregular volumes, with a large number of mesh elements involved (12,000,000 tetrahedrons approximately). These meshes can be used to perform computational simulations by Finite Element Method (FEM). © Published under licence by IOP Publishing Ltd.