A Finite Element Numerical Algorithm for Modelling and Data Fitting in Complex Systems

Numerical modelling methodologies are important by their application to engineering and scientific problems, because there are processes where analytical mathematical expressions cannot be obtained to model them. When the only available information is a set of experimental values for the variables that determine the state of the system, the modelling problem is equivalent to determining the hyper-surface that best fits the data. This paper presents a methodology based on the Galerkin formulation of the finite elements method to obtain representations of relationships that are defined a priori, between a set of variables: y = z(x1, x2,...., xd). These representations are generated from the values of the variables in the experimental data. The approximation, piecewise, is an element of a Sobolev space and has derivatives defined in a general sense into this space. The using of this approach results in the need of inverting a linear system with a structure that allows a fast solver algorithm. The algorithm can be used in a variety of fields, being a multidisciplinary tool. The validity of the methodology is studied considering two real applications: a problem in hydrodynamics and a problem of engineering related to fluids, heat and transport in an energy generation plant. Also a test of the predictive capacity of the methodology is performed using a cross-validation method.

Majmūʻat al-rasāʼil fī ādāb al-baḥth : manuscript, 1869 or 70

1. Risālah fī ādāb al-muṭālaʻah / Ḥāmid ibn Burhān ibn Abī Dharr al-Ghifārī (ff. 1r-6v) -- 2. Ādāb ʻUthmānīyah (ff. 7r-14r) -- 3. Ādāb Sharīfīyah fī al-munāẓarah / al-Sayyid al-Sharīf al-Jurjānī (ff. 15r-27r) -- 4. Sharḥ ādāb al-baḥth fī ʻilm al-munāẓarah lil-ʻAḍudī (ff. 28r-46v) -- 5. Ḥāshiyah ʻalá Sharḥ Ādāb al-Ḥanafī (ff. 47r-55v).

Oeuvres complètes : mises en ordre, revues et annotées d'après les manuscrits de l'auteur.

t.1. 2. éd., 1862.--t.2. 2. éd., 1862.--t.3. 3. éd., 1864.--t.4. 2. éd., 1863.--t.5. 2. éd., 1863.--t.6. 5. éd., 1864.--t.7. 1864.

Le vignole des ouvriers, ou, Méthode facile pour tracer les cinq ordres d'architecture /

[v. 1]: 3. éd., rev. et augm. de quelques notes.

Parameters affecting the performance of wetting and drying in a two-dimensional finite element long wave hydrodynamic model

The marsh porosity method, a type of thin slot wetting and drying algorithm in a two-dimensional finite element long wave hydrodynamic model, is discussed and analyzed to assess model performance. Tests, including comparisons to simple examples and theoretical calculations, examine the effects of varying the marsh porosity parameters. The findings demonstrate that the wetting and drying concept of marsh porosity, often used in finite element hydrodynamic modeling, can behave in a more complex manner than initially expected.

Suppression of spurious intermediate frequency modes in under-integrated elements by combined stiffness/viscous stabilization

Solutions employing perturbation stiffness or viscous hourglass control with one-point quadrature finite elements often exhibit spurious modes in the intermediate frequency range. These spurious frequencies are demonstrated in several examples and their origin is explained. Then it is shown that by critically damping the hourglass modes, these spurious mid-range frequency modes can be suppressed. Estimates of the hourglass frequency and damping coefficients are provided for the plane 4-node quadrilateral and a 4-node shell element. Results are presented that show almost complete annihilation of spurious intermediate frequency modes for both linear and non-linear problems. Copyright (c) 2005 John Wiley & Sons, Ltd.

Sun and sky:does human vision assume a mixture of point and diffuse illumination when interpreting shape-from-shading?

People readily perceive smooth luminance variations as being due to the shading produced by undulations of a 3-D surface (shape-from-shading). In doing so, the visual system must simultaneously estimate the shape of the surface and the nature of the illumination. Remarkably, shape-from-shading operates even when both these properties are unknown and neither can be estimated directly from the image. In such circumstances humans are thought to adopt a default illumination model. A widely held view is that the default illuminant is a point source located above the observer's head. However, some have argued instead that the default illuminant is a diffuse source. We now present evidence that humans may adopt a flexible illumination model that includes both diffuse and point source elements. Our model estimates a direction for the point source and then weights the contribution of this source according to a bias function. For most people the preferred illuminant direction is overhead with a strong diffuse component.

Optimization of chemical plant simulation using double collocation

A method has been constructed for the solution of a wide range of chemical plant simulation models including differential equations and optimization. Double orthogonal collocation on finite elements is applied to convert the model into an NLP problem that is solved either by the VF 13AD package based on successive quadratic programming, or by the GRG2 package, based on the generalized reduced gradient method. This approach is termed simultaneous optimization and solution strategy. The objective functional can contain integral terms. The state and control variables can have time delays. Equalities and inequalities containing state and control variables can be included into the model as well as algebraic equations and inequalities. The maximum number of independent variables is 2. Problems containing 3 independent variables can be transformed into problems having 2 independent variables using finite differencing. The maximum number of NLP variables and constraints is 1500. The method is also suitable for solving ordinary and partial differential equations. The state functions are approximated by a linear combination of Lagrange interpolation polynomials. The control function can either be approximated by a linear combination of Lagrange interpolation polynomials or by a piecewise constant function over finite elements. The number of internal collocation points can vary by finite elements. The residual error is evaluated at arbitrarily chosen equidistant grid-points, thus enabling the user to check the accuracy of the solution between collocation points, where the solution is exact. The solution functions can be tabulated. There is an option to use control vector parameterization to solve optimization problems containing initial value ordinary differential equations. When there are many differential equations or the upper integration limit should be selected optimally then this approach should be used. The portability of the package has been addressed converting the package from V AX FORTRAN 77 into IBM PC FORTRAN 77 and into SUN SPARC 2000 FORTRAN 77. Computer runs have shown that the method can reproduce optimization problems published in the literature. The GRG2 and the VF I 3AD packages, integrated into the optimization package, proved to be robust and reliable. The package contains an executive module, a module performing control vector parameterization and 2 nonlinear problem solver modules, GRG2 and VF I 3AD. There is a stand-alone module that converts the differential-algebraic optimization problem into a nonlinear programming problem.

State detection of bond wires in IGBT modules using eddy current pulsed thermography

Insulated gate bipolar transistor (IGBT) modules are important safety critical components in electrical power systems. Bond wire lift-off, a plastic deformation between wire bond and adjacent layers of a device caused by repeated power/thermal cycles, is the most common failure mechanism in IGBT modules. For the early detection and characterization of such failures, it is important to constantly detect or monitor the health state of IGBT modules, and the state of bond wires in particular. This paper introduces eddy current pulsed thermography (ECPT), a nondestructive evaluation technique, for the state detection and characterization of bond wire lift-off in IGBT modules. After the introduction of the experimental ECPT system, numerical simulation work is reported. The presented simulations are based on the 3-D electromagnetic-thermal coupling finite-element method and analyze transient temperature distribution within the bond wires. This paper illustrates the thermal patterns of bond wires using inductive heating with different wire statuses (lifted-off or well bonded) under two excitation conditions: nonuniform and uniform magnetic field excitations. Experimental results show that uniform excitation of healthy bonding wires, using a Helmholtz coil, provides the same eddy currents on each, while different eddy currents are seen on faulty wires. Both experimental and numerical results show that ECPT can be used for the detection and characterization of bond wires in power semiconductors through the analysis of the transient heating patterns of the wires. The main impact of this paper is that it is the first time electromagnetic induction thermography, so-called ECPT, has been employed on power/electronic devices. Because of its capability of contactless inspection of multiple wires in a single pass, and as such it opens a wide field of investigation in power/electronic devices for failure detection, performance characterization, and health monitoring.

Desenvolvimento de estruturas integradas de filtros e antenas para aplicações em sistemas de comunicações sem fio

This paper presents a study of the integration of filters and microstrip antennas, yielding devices named as filtennas for applications in wireless communications systems. The design of these structures is given from the observation of filtennas based integration between horn antennas and frequency selective surfaces (FSS), used in the band X. The choice of microstrip line structures for the development of a new configuration filtennas justifies the wide application of these transmission lines, in recent decades, always resulting in the production of circuit structures with planar light-weight, compact size, low cost, easy to construct and particularly easy to integrate with other microwave circuits. In addition, the antenna structure considered for the composition of filtennas consists of a planar monopole microstrip to microstrip filters integrated in the feed line of the antenna. In particular, are considered elliptical monopole microstrip (operating in UWB UWB) microstrip filters and (in structures with associated sections in series and / or coupled). In addition, the monopole microstrip has a proper bandwidth and omnidirectional radiation pattern, such that its integration with microstrip filters results in decreased bandwidth, but with slight changes in the radiation pattern. The methods used in the analysis of monopoles, and filters were filtennas finite elements and moments by using commercial software Ansoft Designer and HFSS Ansoft, respectively. Specifically, we analyze the main characteristics of filtennas, such as radiation pattern, gain and bandwidth. Were designed, constructed and measures, several structures filtennas, for validation of the simulated results. Were also used computational tools (CAD) in the process of building prototypes of planar monopoles, filters and filtennas. The prototypes were constructed on substrates of glass-fiber (FR4). Measurements were performed at the Laboratory for Telecommunications UFRN. Comparisons were made between simulated and measured, and found good agreement in the cases considered

Pore water and solid-phase measurements on sediment cores from the Black Sea

The surface sediments in the Black Sea are underlain by extensive deposits of iron (Fe) oxide-rich lake sediments that were deposited prior to the inflow of marine Mediterranean Sea waters ca. 9000 years ago. The subsequent downward diffusion of marine sulfate into the methane-bearing lake sediments has led to a multitude of diagenetic reactions in the sulfate-methane transition zone (SMTZ), including anaerobic oxidation of methane (AOM) with sulfate. While the sedimentary cycles of sulfur (S), methane and Fe in the SMTZ have been extensively studied, relatively little is known about the diagenetic alterations of the sediment record occurring below the SMTZ. Here we combine detailed geochemical analyses of the sediment and pore water with multicomponent diagenetic modeling to study the diagenetic alterations below the SMTZ at two sites in the western Black Sea. We focus on the dynamics of Fe, S and phosphorus (P) and demonstrate that diagenesis has strongly overprinted the sedimentary burial records of these elements. Our results show that sulfate-mediated AOM substantially enhances the downward diffusive flux of sulfide into the deep limnic deposits. During this downward sulfidization, Fe oxides, Fe carbonates and Fe phosphates (e.g. vivianite) are converted to sulfide phases, leading to an enrichment in solid phase S and the release of phosphate to the pore water. Below the sulfidization front, high concentrations of dissolved ferrous Fe (Fe2+) lead to sequestration of downward diffusing phosphate as authigenic vivianite, resulting in a transient accumulation of total P directly below the sulfidization front. Our model results further demonstrate that downward migrating sulfide becomes partly re-oxidized to sulfate due to reactions with oxidized Fe minerals, fueling a cryptic S cycle and thus stimulating slow rates of sulfate-driven AOM (~ 1-100 pmol/cm**3/d) in the sulfate-depleted limnic deposits. However, this process is unlikely to explain the observed release of dissolved Fe2+ below the SMTZ. Instead, we suggest that besides organoclastic Fe oxide reduction, AOM coupled to the reduction of Fe oxides may also provide a possible mechanism for the high concentrations of Fe2+ in the pore water at depth. Our results reveal that methane plays a key role in the diagenetic alterations of Fe, S and P records in Black Sea sediments. The downward sulfidization into the limnic deposits is enhanced through sulfate-driven AOM with sulfate and AOM with Fe oxides may provide a deep source of dissolved Fe2+ that drives the sequestration of P in vivianite below the sulfidization front.

"Sub Tiberio quies"? Judaea’s Political Situation under the Prefects (6-41 c.e.), between Reality and Ideology

Deeply conflicting views on the political situation of Judaea under the Roman prefects (6-41 c.e.) have been offered. According to some scholars, this was a period of persistent political unrest and agitation, whilst according to a widespread view it was a quiescent period of political calm (reflected in Tacitus’ phrase sub Tiberio quies). The present article critically examines again the main available sources –particularly Josephus, the canonical Gospels and Tacitus– in order to offer a more reliable historical reconstruction. The conclusions drawn by this survey calls into question some widespread and insufficiently nuanced views on the period. This, in turn, allows a reflection on the non-epistemic factors which might contribute to explain the origin of such views.

A Finite Element to Perform 3D Shakedown Analysis for Limited Kinematic Hardening Materials

The use of the Design by Analysis (DBA) route is a modern trend in pressure vessel and piping international codes in mechanical engineering. However, to apply the DBA to structures under variable mechanical and thermal loads, it is necessary to assure that the plastic collapse modes, alternate plasticity and incremental collapse (with instantaneous plastic collapse as a particular case), be precluded. The tool available to achieve this target is the shakedown theory. Unfortunately, the practical numerical applications of the shakedown theory result in very large nonlinear optimization problems with nonlinear constraints. Precise, robust and efficient algorithms and finite elements to solve this problem in finite dimension has been a more recent achievements. However, to solve real problems in an industrial level, it is necessary also to consider more realistic material properties as well as to accomplish 3D analysis. Limited kinematic hardening, is a typical property of the usual steels and it should be considered in realistic applications. In this paper, a new finite element with internal thermodynamical variables to model kinematic hardening materials is developed and tested. This element is a mixed ten nodes tetrahedron and through an appropriate change of variables is possible to embed it in a shakedown analysis software developed by Zouain and co-workers for elastic ideally-plastic materials, and then use it to perform 3D shakedown analysis in cases with limited kinematic hardening materials

Easy preoperative planning of deeply located brain lesions using external skin reference and 3-Dimensional surface MRI

Open skull surgery of deeply located intracerebral lesions requires precise determination of the treatment area in 3-dimensional (3-D) space. 3-D MRI can give important additional information in presurgical determination of the surgical approach to the target, taking into account highly functional brain areas and important vascular structures. The day before surgery, a grid composed of 9 tubings intersecting at 90° at 1 cm intervals and filled with a Q1SO4 solution is firmly attached to the skin of the patient’s head in the presumed region of the craniotomy. A 3-D turbo-FLASH sequence is then performed in the sagittal plane after intravenous Gd-DOTA injection on a IT Magnetom. 3-D surface reconstruction of the cortical gyri and sulci is performed. Once the gyri are identified, the 3-D program is then implemented in order to perform a color display of the cortical veins and of the tumor boundaries. The surgical access is then chosen by the surgeon, taking into account highly functional areas. Finally, the boundaries of the tumor are projected on the cortex reconstruction and on the external reference placed on the skin. The entry place for surgery as well as the size of craniotomy are drawn on the skin and the tubed grid is removed. The accuracy of this method tested in 9 patients with deeply located brain tumors or arteriovenous malformations was very satisfactory. In daily practice, this method is a valuable technique providing important clinical information in determining the shortest and safest way through the brain tissue, decreasing possible functional deficit and reducing craniotomy size in cases of difficult to access deep brain areas. Our method does not require a stereotactic frame permanently fixed to the head of the patient during surgery. © 1994 S. Karger AG, Basel.

A finite volume approach to dynamic fluid structure interaction

Abstract not available