A Pre-processing Moving Mesh Method for Discontinuous Galerkin Approximations of Advection-Diffusion-Reaction Problems

We propose a pre-processing mesh re-distribution algorithm based upon harmonic maps employed in conjunction with discontinuous Galerkin approximations of advection-diffusion-reaction problems. Extensive two-dimensional numerical experiments with different choices of monitor functions, including monitor functions derived from goal-oriented a posteriori error indicators are presented. The examples presented clearly demonstrate the capabilities and the benefits of combining our pre-processing mesh movement algorithm with both uniform, as well as, adaptive isotropic and anisotropic mesh refinement.

Estudo experimental e numérico de juntas adesivas em T com aderentes de alumínio

Nos dias de hoje, a ligação adesiva de estruturas complexas que não poderiam ou não seriam tão fáceis de ser fabricadas numa só peça é cada vez mais usual. As juntas adesivas têm vindo a substituir muitos outros métodos de ligação, como por exemplo ligações aparafusadas, rebitas ou soldadas, devido às vantagens de facilidade na sua fabricação, resistência superior e capacidade de unir materiais diferentes. Por esta razão as juntas adesivas têm vindo a ser aplicadas cada vez mais em várias industrias como aeroespacial, aeronáutica, automóvel, naval e calçado. O tipo de adesivo a usar em determinada aplicação é principalmente escolhido consoante as suas características mecânicas e o tipo de resposta pretendida às solicitações impostas. Como exemplo de adesivo resistente e frágil existe o Araldite® AV138. Por outro lado, o adesivo Araldite® 2015 é menos resistente, mas apresenta maior ductilidade e flexibilidade. Além dos adesivos Araldite® comerciais, existem adesivos de poliuretano que combinam características de elevada resistência com características de grande ductilidade e flexibilidade, como por exemplo o Sikaforce® 7752. Esta dissertação tem como objetivo estudar experimentalmente e numericamente, através de modelos de dano coesivo (MDC), o comportamento de diferentes configurações de junta em T quando sujeitas a solicitações de arrancamento. Consideram-se os adesivos anteriormente mencionados para testar as juntas sob diferentes tipos de adesivos. A junta em T é constituída por 2 aderentes em L de alumínio e um aderente base também em alumínio, unidos por uma camada de adesivo. Experimentalmente é feito um estudo da resistência da junta com a variação da espessura dos aderentes em L (tP2). Com a análise numérica são estudadas as distribuições de tensões, evolução do dano, modos de rotura e resistência. Além disso, realizou-se um estudo numérico da existência ou não de adesivo de preenchimento na zona da curvatura dos aderentes em L nas tensões e na resistência da junta. Mostrouse que a variação da geometria nos aderentes em L, a presença de adesivo de preenchimento e o tipo de adesivo têm uma influência direta na resistência de junta. Os ensaios experimentais validaram os resultados numéricos e permitiram concluir que os MDC são uma técnica precisa para o estudo das geometrias das juntas em T.

An A Posteriori Error Indicator for Discontinuous Galerkin Approximations of Fourth Order Elliptic Problems

We introduce a residual-based a posteriori error indicator for discontinuous Galerkin discretizations of the biharmonic equation with essential boundary conditions. We show that the indicator is both reliable and efficient with respect to the approximation error measured in terms of a natural energy norm, under minimal regularity assumptions. We validate the performance of the indicator within an adaptive mesh refinement procedure and show its asymptotic exactness for a range of test problems.

Análise hidromecânica do problema de fraturamento hidráulico

Dissertação (mestrado)—Universidade de Brasília, Faculdade de Tecnologia, Departamento de Engenharia Civil e Ambiental, 2016.

Second-order elliptic PDE with discontinuous boundary data

We shall consider the weak formulation of a linear elliptic model problem with discontinuous Dirichlet boundary conditions. Since such problems are typically not well-defined in the standard H^1-H^1 setting, we will introduce a suitable saddle point formulation in terms of weighted Sobolev spaces. Furthermore, we will discuss the numerical solution of such problems. Specifically, we employ an hp-discontinuous Galerkin method and derive an L^2-norm a posteriori error estimate. Numerical experiments demonstrate the effectiveness of the proposed error indicator in both the h- and hp-version setting. Indeed, in the latter case exponential convergence of the error is attained as the mesh is adaptively refined.

Triangular mesh parameterization with trimmed surfaces

Given a 2manifold triangular mesh \(M \subset {\mathbb {R}}^3\), with border, a parameterization of \(M\) is a FACE or trimmed surface \(F=\{S,L_0,\ldots, L_m\}\) -- \(F\) is a connected subset or region of a parametric surface \(S\), bounded by a set of LOOPs \(L_0,\ldots ,L_m\) such that each \(L_i \subset S\) is a closed 1manifold having no intersection with the other \(L_j\) LOOPs -- The parametric surface \(S\) is a statistical fit of the mesh \(M\) -- \(L_0\) is the outermost LOOP bounding \(F\) and \(L_i\) is the LOOP of the ith hole in \(F\) (if any) -- The problem of parameterizing triangular meshes is relevant for reverse engineering, tool path planning, feature detection, redesign, etc -- Stateofart mesh procedures parameterize a rectangular mesh \(M\) -- To improve such procedures, we report here the implementation of an algorithm which parameterizes meshes \(M\) presenting holes and concavities -- We synthesize a parametric surface \(S \subset {\mathbb {R}}^3\) which approximates a superset of the mesh \(M\) -- Then, we compute a set of LOOPs trimming \(S\), and therefore completing the FACE \(F=\ {S,L_0,\ldots ,L_m\}\) -- Our algorithm gives satisfactory results for \(M\) having low Gaussian curvature (i.e., \(M\) being quasi-developable or developable) -- This assumption is a reasonable one, since \(M\) is the product of manifold segmentation preprocessing -- Our algorithm computes: (1) a manifold learning mapping \(\phi : M \rightarrow U \subset {\mathbb {R}}^2\), (2) an inverse mapping \(S: W \subset {\mathbb {R}}^2 \rightarrow {\mathbb {R}}^3\), with \ (W\) being a rectangular grid containing and surpassing \(U\) -- To compute \(\phi\) we test IsoMap, Laplacian Eigenmaps and Hessian local linear embedding (best results with HLLE) -- For the back mapping (NURBS) \(S\) the crucial step is to find a control polyhedron \(P\), which is an extrapolation of \(M\) -- We calculate \(P\) by extrapolating radial basis functions that interpolate points inside \(\phi (M)\) -- We successfully test our implementation with several datasets presenting concavities, holes, and are extremely nondevelopable -- Ongoing work is being devoted to manifold segmentation which facilitates mesh parameterization

Influência das tensões tectônicas na permeabilidade da formação Açu - Bacia Potiguar

The behavior of the fluid flux in oil fields is influenced by different factors and it has a big impact on the recovery of hydrocarbons. There is a need of evaluating and adapting the actual technology to the worldwide reservoirs reality, not only on the exploration (reservoir discovers) but also on the development of those that were already discovered, however not yet produced. The in situ combustion (ISC) is a suitable technique for these recovery of hydrocarbons, although it remains complex to be implemented. The main objective of this research was to study the application of the ISC as an advanced oil recovery technique through a parametric analysis of the process using vertical wells within a semi synthetic reservoir that had the characteristics from the brazilian northwest, in order to determine which of those parameters could influence the process, verifying the technical and economical viability of the method on the oil industry. For that analysis, a commercial reservoir simulation program for thermal processes was used, called steam thermal and advanced processes reservoir simulator (STARS) from the computer modeling group (CMG). This study aims, through the numerical analysis, find results that help improve mainly the interpretation and comprehension of the main problems related to the ISC method, which are not yet dominated. From the results obtained, it was proved that the mediation promoted by the thermal process ISC over the oil recovery is very important, with rates and cumulated production positively influenced by the method application. It was seen that the application of the method improves the oil mobility as a function of the heating when the combustion front forms inside the reservoir. Among all the analyzed parameters, the activation energy presented the bigger influence, it means, the lower the activation energy the bigger the fraction of recovered oil, as a function of the chemical reactions speed rise. It was also verified that the higher the enthalpy of the reaction, the bigger the fraction of recovered oil, due to a bigger amount of released energy inside the system, helping the ISC. The reservoir parameters: porosity and permeability showed to have lower influence on the ISC. Among the operational parameters that were analyzed, the injection rate was the one that showed a stronger influence on the ISC method, because, the higher the value of the injection rate, the higher was the result obtained, mainly due to maintaining the combustion front. In connection with the oxygen concentration, an increase of the percentage of this parameter translates into a higher fraction of recovered oil, because the quantity of fuel, helping the advance and the maintenance of the combustion front for a longer period of time. About the economic analysis, the ISC method showed to be economically feasible when evaluated through the net present value (NPV), considering the injection rates: the higher the injection rate, the higher the financial incomes of the final project

Simulación de excavaciones urbanas mediante modelos de elementos finitos en depósitos sedimentarios de origen aluvial en Sabaneta

En este trabajo se realizan simulaciones de excavaciones profundas en suelos de origen aluvial en la ciudad de Sabaneta, mediante el empleo de modelos en elementos finitos integrados por el software PLAXIS® -- Los desplazamientos horizontales son comparados con mediciones de inclinómetros instalados en el trasdós del muro diafragma anclado del proyecto Centro Comercial Mayorca Fase III, localizado en el municipio de Sabaneta, Antioquia -- Finalmente, se concluye acerca de la sensibilidad de los parámetros más relevantes según el modelo constitutivo empleado y la viabilidad en su aplicación para la solución del problema evaluado

Simulação do escoamento monofásico em um estágio de uma bomba centrífuga utilizando técnicas de fluidodinâmica computacional

Oil production and exploration techniques have evolved in the last decades in order to increase fluid flows and optimize how the required equipment are used. The base functioning of Electric Submersible Pumping (ESP) lift method is the use of an electric downhole motor to move a centrifugal pump and transport the fluids to the surface. The Electric Submersible Pumping is an option that has been gaining ground among the methods of Artificial Lift due to the ability to handle a large flow of liquid in onshore and offshore environments. The performance of a well equipped with ESP systems is intrinsically related to the centrifugal pump operation. It is the pump that has the function to turn the motor power into Head. In this present work, a computer model to analyze the three-dimensional flow in a centrifugal pump used in Electric Submersible Pumping has been developed. Through the commercial program, ANSYS® CFX®, initially using water as fluid flow, the geometry and simulation parameters have been defined in order to obtain an approximation of what occurs inside the channels of the impeller and diffuser pump in terms of flow. Three different geometry conditions were initially tested to determine which is most suitable to solving the problem. After choosing the most appropriate geometry, three mesh conditions were analyzed and the obtained values were compared to the experimental characteristic curve of Head provided by the manufacturer. The results have approached the experimental curve, the simulation time and the model convergence were satisfactory if it is considered that the studied problem involves numerical analysis. After the tests with water, oil was used in the simulations. The results were compared to a methodology used in the petroleum industry to correct viscosity. In general, for models with water and oil, the results with single-phase fluids were coherent with the experimental curves and, through three-dimensional computer models, they are a preliminary evaluation for the analysis of the two-phase flow inside the channels of centrifugal pump used in ESP systems

Comparative study of human age estimation based on hand-crafted and deep face features

In the past few years, human facial age estimation has drawn a lot of attention in the computer vision and pattern recognition communities because of its important applications in age-based image retrieval, security control and surveillance, biomet- rics, human-computer interaction (HCI) and social robotics. In connection with these investigations, estimating the age of a person from the numerical analysis of his/her face image is a relatively new topic. Also, in problems such as Image Classification the Deep Neural Networks have given the best results in some areas including age estimation. In this work we use three hand-crafted features as well as five deep features that can be obtained from pre-trained deep convolutional neural networks. We do a comparative study of the obtained age estimation results with these features.

An analytical model for web-warping in variable width flexible roll forming

The development of ultra/advanced high strength steels (U/AHSS) has challenged traditional forming methods due to their higher strength and reduced formability. An alternative method is flexible roll forming, which allows the manufacture of sheet metal of high strength and limited ductility into complex and weight-optimized components. However, one major problem in flexible roll forming is the web-warping defect, which is the deviation in height of the web over the length of the profile. The authors’ previous work developed an analytical model to predict the magnitude of web-warping. That model was purely geometric and neglected the effect of material properties. This work develops an analytical solution for the prediction of web-warping that considers both geometric and material parameters. The model results were validated by comparison with numerical and experimental results. The impact of this new model will be the ability to provide a rapid initial design assessment before an intensive numerical analysis of flexible roll forming is conducted.

Separation of longitudinal and flexural wave in a cylindrical structure based on sensor arrangement for non-destructive evaluation

Low strain integrity testing is commonly used to assess the in situ condition of the poles or piles. For poles, it is important to calculate the embedment length and location of damage which is highly influenced by the accurate determination of the wave velocity. In general, depending on impact location and orientation, both longitudinal and bending waves may generate inside the pole, and these two waves have very distinct characteristics and wave velocity. These differences are even more prominent in the low frequency which is usually induced in the low strain non-destructive testing. Consequently, it will be useful if these two waves can be separated for the condition assessment of the poles. In this paper, a numerical analysis is performed on a pole considering that both waves are generated, and a method is proposed to differentiate these two waves based on an appropriate sensor arrangement that includes the location and the orientation of the sensors. Continuous wavelet transform is applied on the numerical signal to calculate the phase velocity of the waves and compared with analytical phase velocity curves. From the results, it can be seen that appropriate location and orientation of the sensors can separate the longitudinal and flexural waves as they match significantly well with the corresponding analytical phase velocity curves of these two waves.

Directional moisture transfer through a wild silkworm cocoon wall

A silkworm cocoon is a porous biological structure with multiple protective functions. In the current work, the authors have used both experimental and numerical methods to reveal the unique moisture transfer characteristics through a wild Antheraea pernyi silkworm cocoon wall, in comparison with the long-domesticated Bombyx mori silkworm cocoon walls. The water vapor transmission and water vapor permeability (WVP) properties show that the A. pernyi cocoons exhibit directional moisture transfer behavior, with easier moisture transfer from inside out than outside in [e.g., the average WVP is 0.057 g/(h m bar) from inside out and is 0.034 g/(h m bar) from outside in]. Numerical analysis shows that the cubic mineral crystals in the outer section of the A. pernyi cocoon wall create a rough surface that facilitates air turbulence and promotes disturbance amplitude of the flow field, leading to lengthened water vapor transfer path and increased tortuosity of the moist air. It also indicates the vortex of water vapor can be generated in the outer section of cocoon wall, which increases the diffusion distance of water vapor and enhances the turbulence kinetic energy and turbulence eddy dissipation, signifying higher moisture resistance in the outer section. The difference in moisture resistance of the multiple A. pernyi cocoon layers is largely responsible for the unique directional moisture transfer behavior of this wild silkworm cocoon. These findings may inspire a biomimicry approach to develop novel lightweight moisture management materials and structures.

Constitutive modelling of high strength titanium alloy Ti-6Al-4 V for sheet forming applications at room temperature

To enable the design and optimisation of forming processes at room temperature the material behaviour of Ti-6Al-4 V needs to be accurately represented in numerical analysis and this requires an advanced material model. In particular, an accurate representation of the shape and size of the yield locus as well as its evolution during forming is important. In this study a rigorous set of experiments on the quasi-static deformation behaviour of a Ti-6Al-4 V alloy sheet sample at room temperature was conducted for various loading conditions and a constitutive material model developed. To quantify the anisotropy and asymmetry properties, tensile and compression tests were carried out for different specimen orientations. To examine the Bauschinger effect and the transient hardening behaviour in - plane tensile - compression and compression - tensile tests were performed. Balanced biaxial and plane strain tension tests were conducted to construct and validate the yield surface of the Ti-6Al-4 V alloy sheet sample at room temperature. A recently proposed anisotropic elastic-plastic constitutive material model, so-called HAH, was employed to describe the behaviour, in particular for load reversals. The HAH yield surface is composed of a stable component, which includes plastic anisotropy and is distorted by a fluctuating component. The key of the formulation is the use of a suitable yield function that reproduces the experimental observations well for the stable component. Meanwhile, the rapid evolution of the material structure must be captured at the macro - scale level by the fluctuating component embedded in the HAH model. Compared to conventional hardening equations, the proposed model leads to higher accuracy in predicting the Bauschinger effect and the transient hardening behaviour for the Ti-6Al-4 V sheet sample tested at room temperature.

Evaluating temporal changes in hydraulic conductivities near karst-terrain dams: Dokan Dam (Kurdistan-Iraq)

Dam sites provide an outstanding opportunity to explore dynamic changes in the groundwater flow regime because of the high hydraulic gradient rapidly induced in their surroundings. This paper investigates the temporal changes of the hydraulic conductivities of the rocks and engineered structures via a thorough analysis of hydrological data collected at the Dokam Dam, Iraq, and a numerical model that simulates the Darcian component of the seepage. Analysis of the data indicates increased seepage with time and suggests that the hydraulic conductivity of the rocks increased as the conductivity of the grout curtain decreased. Conductivity changes on the order of 10-8m/s, in a 20-yr period were quantified using the numerical analysis. It is postulated that the changes in hydraulic properties in the vicinity of Dokan Dam are due to suspension of fine materials, interbedded in small fissures in the rocks, and re-settlement of these materials along the curtain. Consequently, the importance of the grout curtain to minimize the downstream seepage, not only as a result of the conductivity contrast with the rocks, but also as a barrier to suspended clay sediments, is demonstrated. The numerical analysis also helped us to estimate the proportion of the disconnected karstic conduit flow to the overall flow.