Efficient recovery-based error estimation for the smoothed finite element method for smooth and singular linear elasticity

[EN] An error control technique aimed to assess the quality of smoothed finite element approximations is presented in this paper. Finite element techniques based on strain smoothing appeared in 2007 were shown to provide significant advantages compared to conventional finite element approximations. In particular, a widely cited strength of such methods is improved accuracy for the same computational cost. Yet, few attempts have been made to directly assess the quality of the results obtained during the simulation by evaluating an estimate of the discretization error. Here we propose a recovery type error estimator based on an enhanced recovery technique. The salient features of the recovery are: enforcement of local equilibrium and, for singular problems a ¿smooth + singular¿ decomposition of the recovered stress. We evaluate the proposed estimator on a number of test cases from linear elastic structural mechanics and obtain efficient error estimations whose effectivities, both at local and global levels, are improved compared to recovery procedures not implementing these features.Stephane Bordas would like to thank the partial financial support of the Royal Academy of Engineering and of the Leverhulme Trust for his Senior Research Fellowship Towards the next generation surgical simulators as well as the financial support for Octavio A. Gonzalez-Estrada and Stephane Bordas from the UK Engineering Physical Science Research Council (EPSRC) under grant EP/G042705/1 Increased Reliability for Industrially Relevant Automatic Crack Growth Simulation with the eXtended Finite Element Method. Stephane Bordas also thanks partial financial support of the European Research Council Starting Independent Research Grant (ERC Stg grant agreement No. 279578) and the FP7 Initial Training Network Funding under grant number 289361 "Integrating Numerical Simulation and Geometric Design Technology, INSIST". This work has been carried out within the framework of the research project DPI2010-20542 of the Ministerio de Ciencia e Innovacion (Spain). The financial support from Universitat Politecnica de Valencia, PROMETEO/2012/023 and Generalitat Valenciana are also acknowledged.González Estrada, OA.; Natarajan, S.; J.J. Ródenas; Nguyen-Xuan, H.; Bordas, S. (2013). Efficient recovery-based error estimation for the smoothed finite element method for smooth and singular linear elasticity. Computational Mechanics. 52(1):37-52. https://doi.org/10.1007/s00466-012-0795-6S3752521Liu GR, Dai KY, Nguyen TT (2006) A smoothed finite element method for mechanics problems. Comput Mech 39(6): 859–877. doi: 10.1007/s00466-006-0075-4Liu GR, Nguyen TT, Dai KY, Lam KY (2007) Theoretical aspects of the smoothed finite element method (SFEM). Int J Numer Methods Eng 71(8): 902–930Nguyen-Xuan H, Bordas SPA, Nguyen-Dang H (2008) Smooth finite element methods: convergence, accuracy and properties. Int J Numer Methods Eng 74(2): 175–208. doi: 10.1002/nmeBordas SPA, Natarajan S (2010) On the approximation in the smoothed finite element method (SFEM). Int J Numer Methods Eng 81(5): 660–670. doi: 10.1002/nmeZhang HH, Liu SJ, Li LX (2008) On the smoothed finite element method. Int J Numer Methods Eng 76(8): 1285–1295. doi: 10.1002/nme.2460Nguyen-Thoi T, Liu G, Lam K, Zhang G. (2009) A face-based smoothed finite element method (FS-FEM) for 3D linear and nonlinear solid mechanics using 4-node tetrahedral elements. Int J Numer Methods Eng 78: 324–353Liu G, Nguyen-Thoi T, Lam K (2009) An edge-based smoothed finite element method (ES-FEM) for static, free and forced vibration analyses of solids. J Sound Vib 320: 1100–1130Liu G, Nguyen-Thoi T, Nguyen-Xuan H, Lam K (2009) A node based smoothed finite element method (NS-FEM) for upper bound solution to solid mechanics problems. Comput Struct 87: 14–26Liu G. Smoothed Finite Element Methods. CRC Press, 2010Liu G, Nguyen-Xuan H, Nguyen-Thoi T (2010) A theoretical study on the smoothed FEM (SFEM) models: Properties, accuracy and convergence rates. Int J Numer Methods Biomed Eng 84: 1222–1256Nguyen T, Liu G, Dai K, Lam K (2007) smoothed finite element method. Tsinghua Sci Technol 12: 497–508Hung NX, Bordas S, Hung N (2009) Addressing volumetric locking and instabilities by selective integration in smoothed finite element. Commun Numer Methods Eng 25: 19–34Nguyen-Xuan H, Rabczuk T, Bordas S, Debongnie JF (2008) A smoothed finite element method for plate analysis. Comput Methods Appl Mech Eng 197: 1184–1203Nguyen NT, Rabczuk T, Nguyen-Xuan H, Bordas S (2008) A smoothed finite element method for shell analysis. Comput Methods Appl Mech Eng 198: 165–177Bordas SPA, Rabczuk T, Hung NX, Nguyen VP, Natarajan S, Bog T, óuan DM, Hiep NV (2010) Strain smoothing in FEM and XFEM. Comput Struct 88(23–24): 1419–1443. doi: 10.1016/j.compstruc.2008.07.006Bordas SP, Natarajan S, Kerfriden P, Augarde CE, Mahapatra DR, Rabczuk T, Pont SD (2011) On the performance of strain smoothing for óuadratic and enriched finite element approximations (XFEM/GFEM/PUFEM). Int J Numer Methods Biomed Eng 86: 637–666Liu G, Nguyen-Thoi T, Nguyen-Xuan H, Dai K, Lam K (2009) On the essence and the evaluation of the shape functions for the smoothed finite element method (SFEM). Int J Numer Methods Eng 77: 1863–1869. doi: 10.1002/nme.2587Strouboulis T, Zhang L, Wang D, Babuška I. (2006) A posteriori error estimation for generalized finite element methods. Comput Methods Appl Mech Eng 195(9–12): 852–879Bordas SPA, Duflot M (2007) Derivative recovery and a posteriori error estimate for extended finite elements. Comput Methods Appl Mech Eng 196(35–36): 3381–3399Xiao óZ, Karihaloo BL (2004) Statically admissible stress recovery using the moving least sóuares technique. In: Topping BHV, Soares CAM (eds) Progress in computational structures technology. Saxe-Coburg Publications, Stirling, pp 111–138Ródenas JJ, González-Estrada OA, Tarancón JE, Fuenmayor FJ (2008) A recovery-type error estimator for the extended finite element method based on singular + smooth stress field splitting. Int J Numer Methods Eng 76(4): 545–571. doi: 10.1002/nme.2313Panetier J, Ladevèze P, Chamoin L (2010) Strict and effective bounds in goal-oriented error estimation applied to fracture mechanics problems solved with XFEM. Int J Numer Methods Eng 81(6): 671–700Barros FB, Proenca SPB, de Barcellos CS (2004) On error estimator and p-adaptivity in the generalized finite element method. Int J Numer Methods Eng 60(14):2373–2398. doi: 10.1002/nme.1048Nguyen-Thoi T, Liu G, Nguyen-Xuan H, Nguyen-Tran C (2011) Adaptive analysis using the node-based smoothed finite element method (NS-FEM). Int J Numer Methods Biomed Eng 27(2): 198–218. doi: 10.1002/cnmGonzález-Estrada OA, Ródenas JJ, Bordas SPA, Duflot M, Kerfriden P, Giner E (2012) On the role of enrichment and statical admissibility of recovered fields in a-posteriori error estimation for enriched finite element methods. Eng Comput 29(8)Zienkiewicz OC, Zhu JZ (1987) A simple error estimator and adaptive procedure for practical engineering analysis. Int J Numer Methods Eng 24(2): 337–357Ródenas JJ, González-Estrada OA, Díez P, Fuenmayor FJ (2010) Accurate recovery-based upper error bounds for the extended finite element framework. Comput Methods Appl Mech Eng 199(37–40): 2607–2621Williams ML (1952) Stress singularities resulting from various boundary conditions in angular corners of plate in extension. J Appl Mech 19: 526–534Szabó BA, Babuška I (1991) Finite element analysis. Wiley, New YorkBarber JR. (2010) Elasticity. Series: solid mechanics and its application, 3rd edn. Springer, DordrechtChen JS, Wu CT, Yoon S, You Y (2001) A stabilized conforming nodal integration for Galerki mesh-free methods. Int J Numer Methods Eng 50: 435–466Yoo J, Moran B, Chen J (2004) Stabilized conforming nodal integration in the natural element method. Int J Numer Methods Eng 60: 861–890Zienkiewicz OC, Zhu JZ (1992) The superconvergent patch recovery and a posteriori error estimates. Part 1: The recovery technique. Int J Numer Methods Eng 33(7): 1331–1364Zienkiewicz OC, Zhu JZ (1992) The superconvergent patch recovery and a posteriori error estimates. Part 2: Error estimates and adaptivity. Int J Numer Methods Eng 33(7): 1365–1382Wiberg NE, Abdulwahab F (1993) Patch recovery based on superconvergent derivatives and eóuilibrium. Int J Numer Methods Eng 36(16): 2703–2724. doi: 10.1002/nme.1620361603Blacker T, Belytschko T (1994) Superconvergent patch recovery with eóuilibrium and conjoint interpolant enhancements. Int J Numer Methods Eng 37(3): 517–536Stein E, Ramm E, Rannacher R (2003) Error-controlled adaptive finite elements in solid mechanics. Wiley, ChichesterDuflot M, Bordas SPA (2008) A posteriori error estimation for extended finite elements by an extended global recovery. Int J Numer Methods Eng 76: 1123–1138. doi: 10.1002/nmeBordas SPA, Duflot M, Le P (2008) A simple error estimator for extended finite elements. Commun Numer Methods Eng 24(11): 961–971Ródenas JJ, Tur M, Fuenmayor FJ, Vercher A (2007) Improvement of the superconvergent patch recovery technique by the use of constraint eóuations: the SPR-C technique. Int J Numer Methods Eng 70(6): 705–727. doi: 10.1002/nme.1903Díez P, Ródenas JJ, Zienkiewicz OC (2007) Eóuilibrated patch recovery error estimates: simple and accurate upper bounds of the error. Int J Numer Methods Eng 69(10): 2075–2098. doi: 10.1002/nmeYau J, Wang S, Corten H (1980) A mixed-mode crack analysis of isotropic solids using conservation laws of elasticity. J Appl Mech 47(2): 335–341Ródenas JJ, González-Estrada OA, Fuenmayor FJ, Chinesta F (2010) Upper bounds of the error in X-FEM based on a moving least sóuares (MLS) recovery technique. In: Khalili N, Valliappan S, Li ó, Russell A (eds) 9th World congress on computational mechanics (WCCM9). 4th Asian Pacific Congress on computational methods (APCOM2010). Centre for Infrastructure Engineering and SafetyRódenas JJ, González-Estrada OA, Díez P, Fuenmayor FJ (2007) Upper bounds of the error in the extended finite element method by using an eóuilibrated-stress patch recovery technique. In: International conference on adaptive modeling and simulation (ADMOS 2007). International Center for Numerical Methods in Engineering (CIMNE), pp 210–213Menk A, Bordas S (2010) Numerically determined enrichment function for the extended finite element method and applications to bi-material anisotropic fracture and polycrystals. Int J Numer Methods Eng 83: 805–828Menk A, Bordas S (2011) Crack growth calculations in solder joints based on microstructural phenomena with X-FEM. Comput Mater Sci 3: 1145–1156Ródenas JJ (2001) Error de discretización en el cálculo de sensibilidades mediante el método de los elementos finitos. PhD Thesis, Universidad Politécnica de ValenciaAinsworth M, Oden JT (2000) A posteriori error estimation in finite element analysis. Wiley, Chicheste

Risk Factors of Streptococcus suis Infection in Vietnam. A Case-Control Study

Background: Streptococcus suis infection, an emerging zoonosis, is an increasing public health problem across South East Asia and the most common cause of acute bacterial meningitis in adults in Vietnam. Little is known of the risk factors underlying the disease. Methods and Findings: A case-control study with appropriate hospital and matched community controls for each patient was conducted between May 2006 and June 2009. Potential risk factors were assessed using a standardized questionnaire and investigation of throat and rectal S. suis carriage in cases, controls and their pigs, using real-time PCR and culture of swab samples. We recruited 101 cases of S. suis meningitis, 303 hospital controls and 300 community controls. By multivariate analysis, risk factors identified for S. suis infection as compared to either control group included eating "high risk" dishes, including such dishes as undercooked pig blood and pig intestine (OR1 = 2.22; 95% CI = [1.15-4.28] and OR2 = 4.44; 95% CI = [2.15-9.15]), occupations related to pigs (OR1 = 3.84; 95% CI = [1.32-11.11] and OR2 = 5.52; 95% CI = [1.49-20.39]), and exposures to pigs or pork in the presence of skin injuries (OR1 = 7.48; 95% CI = [1.97-28.44] and OR2 = 15.96; 95% CI = [2.97-85.72]). S. suis specific DNA was detected in rectal and throat swabs of 6 patients and was cultured from 2 rectal samples, but was not detected in such samples of 1522 healthy individuals or patients without S. suis infection. Conclusions: This case control study, the largest prospective epidemiological assessment of this disease, has identified the most important risk factors associated with S. suis bacterial meningitis to be eating 'high risk' dishes popular in parts of Asia, occupational exposure to pigs and pig products, and preparation of pork in the presence of skin lesions. These risk factors can be addressed in public health campaigns aimed at preventing S. suis infectio

Molecular Characterization of HIV-1 CRF01_AE in Mekong Delta, Vietnam, and Impact of T-Cell Epitope Mutations on HLA Recognition (ANRS 12159)

To date, 11 HIV-1 subtypes and 48 circulating recombinant forms have been described worldwide. The underlying reason why their distribution is so heterogeneous is not clear. Host genetic factors could partly explain this distribution. The aim of this study was to describe HIV-1 strains circulating in an unexplored area of Mekong Delta, Vietnam, and to assess the impact of optimal epitope mutations on HLA binding.We recruited 125 chronically antiretroviral-naive HIV-1-infected subjects from five cities in the Mekong Delta. We performed high-resolution DNA typing of HLA class I alleles, sequencing of Gag and RT-Prot genes and phylogenetic analysis of the strains. Epitope mutations were analyzed in patients bearing the HLA allele restricting the studied epitope. Optimal wild-type epitopes from the Los Alamos database were used as reference. T-cell epitope recognition was predicted using the immune epitope database tool according to three different scores involved in antigen processing (TAP and proteasome scores) and HLA binding (MHC score). with a Vietnamese specificity held by two different haplotypes. The percentage of homology between Mekong and B consensus HIV-1 sequences was above 85%. Divergent epitopes had TAP and proteasome scores comparable with wild-type epitopes. MHC scores were significantly lower in divergent epitopes with a mean of 2.4 (±0.9) versus 2 (±0.7) in non-divergent ones (p<0.0001).Our study confirms the wide predominance of CRF01_AE in the Mekong Delta where patients harbor a specific HLA pattern. Moreover, it demonstrates the lower MHC binding affinity among divergent epitopes. This weak immune pressure combined with a narrow genetic diversity favors immune escape and could explain why CRF01_AE is still predominant in Vietnam, particularly in the Mekong area

Systematic meta-analyses and field synopsis of genetic and epigenetic studies in paediatric inflammatory bowel disease

We provide a comprehensive field synopsis of genetic and epigenetic associations for paediatric Inflammatory Bowel Disease (IBD). A systematic review was performed and included 84 genetic association studies reporting data for 183 polymorphisms in 71 genes. Meta-analyses were conducted for 20 SNPs in 10 genes of paediatric Crohn’s disease (CD) and for 8 SNPs in 5 genes of paediatric ulcerative colitis (UC). Five epigenetic studies were also included, but formal meta-analysis was not possible. Venice criteria and Bayesian false discovery probability test were applied to assess the credibility of associations. Nine SNPs in 4 genes were considered to have highly credible associations with paediatric CD, of which four variants (rs2066847, rs12521868, rs26313667, rs1800629) were not previously identified in paediatric GWAS. Differential DNA methylation in NOD2 and TNF-α, dysregulated expression in let-7 and miR-124 were associated with paediatric IBD, but not as yet replicated. Highly credible SNPs associated with paediatric IBD have also been implicated in adult IBD, with similar magnitudes of associations. Early onset and distinct phenotypic features of paediatric IBD might be due to distinct epigenetic changes, but these findings need to be replicated. Further progress identifying genetic and epigenetic susceptibility of paediatric IBD will require international collaboration, population diversity and harmonization of protocols

Clinical and virological features of acute HBV-related hepatitis in southern vietnam

Background: Despite the availability of effective vaccines, hepatitis B virus (HBV) infection is still frequent worldwide, and accounts for significant morbidity and mortality. However, data of acute HBeAg negative hepatitis still remain limited. Aims and Methods: To understand clinical pictures of acute HBV hepatitis and its natural evolution, a prospective study was conducted in adult patients. Results: Ninety patients were enrolled between March 2004 and April 2005 at Hospital for Tropical Diseases in Ho Chi Minh city. The prevalence of HBeAg negative was 53%. No significant difference was found in clinical characteristics and laboratory findings between HBeAg positive and negative patients. HBV-DNA was detected in 75% and 88% HBe negative and positive patients, respectively, where the frequency of ALT below 400 U/L was significantly higher in HBeAg negative cases (p = 0.01). Six month follow-up was available in 47 patients. HBsAg positivity was found in 16% of HBeAg negative subjects but only in 4.5% of HBeAg positive cases. Thirty two patients had neither HBsAg nor anti-HBs. Conclusions: The clinical and laboratory feature and the outcome at 6 months of HBV acute hepatitis in Vietnam is similar in HBeAg positive and negative patients

Angiostrongylus cantonensis is an Important Cause of Eosinophilic Meningitis in southern Vietnam

We utilised PCR to demonstrate that Angiostrongylus cantonensis was responsible for 67.3% of 55 cases of eosinophilic meningitis from a cohort of 1690 adult patients with CNS infection at a tertiary hospital in southern Vietnam. Longer duration of illness, depressed consciousness and peripheral blood eosinophilia were associated with PCR positivity