ECCOMAS 2012 European Congress on Computational Methods in Applied Sciences and Engineering

Total Page:16

File Type:pdf, Size:1020Kb

ECCOMAS 2012 European Congress on Computational Methods in Applied Sciences and Engineering ECCOMAS 2012 European Congress on Computational Methods in Applied Sciences and Engineering e-book Full Papers September 10-14, 2012 Vienna, Austria ECCOMAS 2012 e-book containing all full papers submitted to the ECCOMAS 2012 Congress (6th European Congress on Computational Methods in Applied Sciences and Engineering), September 10-14, 2012, Vienna, Austria Editors: Josef Eberhardsteiner, Helmut J. Böhm, Franz G. Rammerstorfer © Vienna University of Technology Institute for Mechanics of Materials and Structures Institute of Lightweight Design and Structural Biomechanics © photocredits Gustav Klimt, the Kiss: Belvedere, Vienna (www.belvedere.at); Secession: WienBild ISBN: 978-3-9503537-0-9 ECCOMAS 2012 e-book | I Table of Contents Enhanced formula for a critical velocity of a uniformly moving load 1 Z. Dimitrovová Microscopic model for simulation of traffic flows on multilane highways and crossroads 11 N. Churbanova, A.A. Chechina, I.R. Furmanov, M.A. Trapeznikova Structure-preserving time integrators for thermo-elasticity with heat conduction 18 P.L. Mata, A.J. Lew Analysis of vibro-based isolated building 38 D. Makovicka, D. Makovicka Comparison of some reduction bases approaches for non-linear structural dynamic systems under 47 different excitations F.A. Lülf, D.M. Tran, R. Ohayon Increasing of aerodynamical performances of flying configurations,via global shape optimization 64 A. Nastase Robust vibration control of flexible structures using smart materials with uncertainties 75 K. Zhang, M. Ichchou, G. Scorletti, F. Mieyeville Causality, “superluminality,” and reshaping in undersized waveguides 93 D. Sokolovski , E. Akhmatskaya Simulation of specimen-geometry effect on ceramics strength by randomly distributed damage 104 approach T. Hoshide, H. Sugiyama Analysis of relative slip on press-fitted contact surface between wheel and axle 116 H. Kawashima A finite volume Lagrangian cell centered mimetic approach for computing elasto-plastic deformation of 125 solids in general unstructured grids S.K. Sambasivan, R. Loubere, M. Shashkov Shape optimization of thermoelastic fields for mean compliance minimization 143 E. Katamine, H. Yoshioka Multi-objective optimization of cyclone geometry based on CFD simulations using the desirability 155 function and genetic algorithms K. Elsayed, C. Lacor Study of the effect of mechanical loading on cell cultures in bone tissue engineering 175 M. Cruel, M. Bensidhoum, L. Sudre, G. Puel, T. Hoc On the performance of a parallel direct solver used within a hybrid spectral/finite element solver for 184 incompressible flows X. Dechamps, G. Degrez, M. Rasquin An efficient implementation of simulating annealing in GPUs and its application to calibration of SABR 198 stochastic volatility model A.M. Ferreiro, J.A. García, J.G. López-Salas, C. Vazquez A comprehensive image processing methodology for automatically detecting shield tunnel defects 208 Q. Ai, Y. Yuan Optimization of a Krueger device for a regional aircraft wing configuration 219 D. de Rosa, G. Andreutti, D. Quagliarella Prediction of particle distribution on wall in magnetic separation 230 T. Inaba, T. Sakazume, Y. Yamashita, S. Matsuoka Thermodynamic model for strain-induced crystallization in rubber 237 L. Thien-Nga, J. Guilie, P. Le Tallec Asymptotics of electron flow spin-polarization in quantum waveguides of variable cross-section in 253 presence of magnetic field L.M. Baskin, P. Neittaanmäki, O. Sarafanov Finite element based trajectory planning for distributed-parameter thermal systems 263 T. Meurer ECCOMAS 2012 e-book | II Geometrical and thermal characterization of silica aerogel using classical molecular dynamics 277 J.J. Yeo, T.Y. Ng, Z.S. Liu Description of deformation process due to heat transfer in a two-phase polycrystalline composite 286 structure T. Sadowski, P. Golewski Computational mechanics of knitted structures - critical overview 296 Y.K. Kyosev, W. Renkens Computational modeling of high speed RTM of automotive part involving on-line mixing of highly 304 reactive two-component resin system C. Binetruy, M. Deleglise, S. Comas-Cardona, S.G. Advani Second order linearization effects based on the matrix stiffness 316 A.M. Wahrhaftig, J.A.L. Rocha, R. Brasil Flow over periodic hills - test case for ERCOFTAC knowledge base wiki 325 M. Breuer, C. Rapp, M. Manhart Fluid-structure interaction in turbulent flows: LES predictions and PIV measurements 344 G. De Nayer, M. Breuer, A. Kalmbach Parallel potts model simulation of nanocrystalline grain growth 364 S. Schäfer, D. Zöllner On the motion of continua with non-material boundary conditions - an investigation with one- 375 dimensional examples A. Franze, B.W. Zastrau Numerical modelling of viscoplastic free surface flows in complex 3D geometries 388 K. Nikitin, M. Olshanskii, K. Terehov, Y. Vassilevski Coupling numerical simulations of mold filling processes with numerical simulations on bubble scale for 402 polyurethane foams S.G. Geier, F.W. Bartels, M. Piesche Exact evaluation of singular and near-singular integrals in Galerkin BEM 413 M. Lenoir, N. Salles Nonlinear behavior of thin-layered structures under contact and friction 432 A. Dorninger, E. Karer, J. Gerstmayr Multi-objective optimization for an air pollution problem 448 L.J. Alvarez-Vazquez, N. Garcia-Chan, A. Martinez, M.E. Vazquez-Mendez Comparison of asymptotic and numerical studies of electron flow spin-polarization in quantum 455 waveguide in magnetic field M. Kabardov, P. Neittaanmäki Concurrent multiscale analysis of heterogeneous materials 467 O. Lloberas-Valls, F.P.X. Everdij, D.J. Rixen, A. Simone, L.J. Sluys Remarks on time integration schemes used for high precision control 482 J. Pandtle, R.H. Steinbuch Comparison between CFD and wind tunnel measurements of a roof mounted air ventilated equipment 490 in a high speed train I. Ali, A. Rüter, B. Bahl Numerical investigation of a cavitation phenomenon round a propulsor blade 498 K. Tesch, J.A. Szantyr, P. Flaszynski An approach for modelling the roughness-induced boundary layer transition using transport equations 507 P. Dassler, D. Kožulović, A. Fiala Computational approach to fatigue behaviour of randomly or unidirectional fibre reinforced materials 525 R. Brighenti, A. Carpinteri, D. Scorza Dynamic finite deformation thermo-viscoelasticity using energy-consistent time-integration 536 M. Krüger, M.M. Groß, P. Betsch High-order absorbing boundary conditions and perfectly matched layers for cables and beams laid on 551 elastic supports G. Lancioni Numerical analysis on asymmetrically combined glulam 564 M. Frese, H.J. Blaß ECCOMAS 2012 e-book | III Hydrodynamic applications using ReALE method 574 T. Harribey, J. Breil, P.-H. Maire, M. Shashkov A method for computing waveguide scattering matrices in the presence of point spectrum 588 B.A. Plamenevskii, O. Sarafanov Topology optimization in fluid dynamics using adjoint-based truncated Newton 597 D.I. Papadimitriou, E.M. Papoutsis-Kiachagias, K. Giannakoglou Locally resolving large scale turbulent structures by a hybrid – Lattice Boltzmann and Finite Volume – 611 turbulence model S. Pirker, C. Goniva, C. Kloss, S. Puttinger, P. Seil, S. Schneiderbauer Stochastic modelling and analysis of porous media using first-order perturbation-based homogenization 621 method K.S. Basaruddin, N. Takano Damage-based crack initiation and propagation 631 H. Minnebo, T. Van Hoof, O. Pierard Transport of pollutant solved by radial basis boundary element method 642 K. Kovarik, J. Muzik, D. Sitanyiová Applying CAD/CAE tools in the roll center determination for a double A automotive suspension 653 G.O. Andrade, M.A.A. Nunes, R.D.C. Silva Minimizing the environmental impact of wastewater discharges with SOS 666 A. Martinez, L.J. Alvarez-Vazquez, N. Garcia-Chan, M.E. Vazquez-Mendez High order stochastic Finite Volume Method for the uncertainty quantification in hyperbolic conservation 673 laws with random initial data and flux coefficients S. Tokareva, S. Mishra, C. Schwab Hybrid RANS/LES, PANS and PRNS computations of plane impinging jets 687 S. Kubacki, J. Rokicki, E. Dick Interaction of flexible multibody dynamics with fluids by means of smoothed particle hydrodynamics 703 M. Schörgenhumer, P.G. Gruber, J. Gerstmayr Numerical flow simulations used in industrial problems 723 K. Adamek, J. Kolar Conservative segregated solution method for turbulence model equations in compressible flows 738 C. Morsbach, F. di Mare Approximation of non-convex anisotropic energies via Willmore energy 752 P. Pozzi, P. Reiter Homogenization and modelling of technical textiles 769 S. Fillep, J. Mergheim, P. Steinmann Inversion of seismic reflection data through focusing 781 W.A. Mulder An improved flexibility-based nonlinear frame element endowed with the fiber-free formulation 795 F. Marmo, L. Rosati Modeling of gas separated flows at inlet of suction channels on the basis of stationary discrete vortices 812 O. Averkova, A. Logachev, I. Logachev, K. Logachev A fully unconstrained interior point algorithm for multivariable state and input constrained optimal 832 control problems P. Malisani, F. Chaplais, N. Petit Graph grammar based parallel direct solver for 1D and 2D isogeometric finite element method 849 M. Paszynski, K. Kuznik, V.M. Calo Finite element modeling of timber-concrete composite beams under short-term loadings 859 N. Khorsandnia, H.R. Valipour, K. Crews 2D phase-field analyses of dendrite and axon growth of nerve cell 868 K. Nakagawa, T. Takaki, Y. Morita, E. Nakamachi Simulation of the dynamics of internal combustion engines considering
Recommended publications
  • Program :: Monday 9
    Comandatuba 1 Comandatuba 2 Comandatuba 3 Auditório Transamérica Una Ilhéus São Paulo 3 São Paulo 2 São Paulo 1 Quito Santiago th Room 1 Room 2 Room 3 Room 4 Room 5 Room 6 Room 7 Room 8 Room 9 Room 10 Room 11 8h30-9h15 Opening Ceremony 9h15-10h Plenary Thomas J.R. Hughes 10h-10h30 Coffee Break :: Coffee Break :: Coffee Break :: Coffee Break :: Coffee Break :: Coffee Break :: Coffee Break :: Coffee Break :: Coffee Break :: Coffee Break :: Coffee Break :: Coffee Break :: Coffee Break 10h30-12h30 MS 094 MS 106 MS 099 MS 074 MS 142 MS 133 MS 001 MS 075 MS 085 MS 002 Satellite Symposium 12h30-13h30 Lunch :: Lunch :: Lunch :: Lunch :: Lunch :: Lunch :: Lunch :: Lunch :: Lunch :: Lunch :: Lunch :: Lunch :: Lunch :: Lunch :: Lunch :: Lunch :: Lunch :: Lunch :: Lunch :: Lunch Brasília Semi-Plenary 1a Semi-Plenary 1b Semi-Plenary 2c Wing Semi-Plenary 1d 13h30-14h Pedro M. Reis Kumar Tamma Kam Liu Álvaro Coutinho Semi-Plenary 2a Semi-plenary 2b Semi-Plenary 2c Semi-Plenary 2d 14h-14h30 Eitan Grinspun Ramon Codina Djordje Peric Paulo Lyra 14h30-16h30 MS 094 MS 106 / MS 035 MS 099 MS 074 MS 142 MS 133 MS 001 MS 075 MS 085 MS 002 16h30-17h Coffee Break :: Coffee Break :: Coffee Break :: Coffee Break :: Coffee Break :: Coffee Break :: Coffee Break :: Coffee Break :: Coffee Break :: Coffee Break :: Coffee Break :: Coffee Break :: Coffee Break 17h-19h MS 094 / MS 163 MS 132 MS 099 MS 171 MS 142 / MS 005 MS 116 MS 001 MS 075 MS 085 MS 002 Program :: Monday 9 Brasília 3 Brasília 2 Brasília 1 Buenos Aires Montevideo Room Room Room Room Room Room Room Room
    [Show full text]
  • 14Th U.S. National Congress on Computational Mechanics
    14th U.S. National Congress on Computational Mechanics Montréal • July 17-20, 2017 Congress Program at a Glance Sunday, July 16 Monday, July 17 Tuesday, July 18 Wednesday, July 19 Thursday, July 20 Registration Registration Registration Registration Short Course 7:30 am - 5:30 pm 7:30 am - 5:30 pm 7:30 am - 5:30 pm 7:30 am - 11:30 am Registration 8:00 am - 9:30 am 8:30 am - 9:00 am OPENING PL: Tarek Zohdi PL: Andrew Stuart PL: Mark Ainsworth PL: Anthony Patera 9:00 am - 9.45 am Chair: J.T. Oden Chair: T. Hughes Chair: L. Demkowicz Chair: M. Paraschivoiu Short Courses 9:45 am - 10:15 am Coffee Break Coffee Break Coffee Break Coffee Break 9:00 am - 12:00 pm 10:15 am - 11:55 am Technical Session TS1 Technical Session TS4 Technical Session TS7 Technical Session TS10 Lunch Break 11:55 am - 1:30 pm Lunch Break Lunch Break Lunch Break CLOSING aSPL: Raúl Tempone aSPL: Ron Miller aSPL: Eldad Haber 1:30 pm - 2:15 pm bSPL: Marino Arroyo bSPL: Beth Wingate bSPL: Margot Gerritsen Short Courses 2:15 pm - 2:30 pm Break-out Break-out Break-out 1:00 pm - 4:00 pm 2:30 pm - 4:10 pm Technical Session TS2 Technical Session TS5 Technical Session TS8 4:10 pm - 4:40 pm Coffee Break Coffee Break Coffee Break Congress Registration 2:00 pm - 8:00 pm 4:40 pm - 6:20 pm Technical Session TS3 Poster Session TS6 Technical Session TS9 Reception Opening in 517BC Cocktail Coffee Breaks in 517A 7th floor Terrace Plenary Lectures (PL) in 517BC 7:00 pm - 7:30 pm Cocktail and Banquet 6:00 pm - 8:00 pm Semi-Plenary Lectures (SPL): Banquet in 517BC aSPL in 517D 7:30 pm - 9:30 pm Viewing of Fireworks bSPL in 516BC Fireworks and Closing Reception Poster Session in 517A 10:00 pm - 10:30 pm on 7th floor Terrace On behalf of Polytechnique Montréal, it is my pleasure to welcome, to Montreal, the 14th U.S.
    [Show full text]
  • Analysis-Guided Improvements of the Material Point Method
    ANALYSIS-GUIDED IMPROVEMENTS OF THE MATERIAL POINT METHOD by Michael Dietel Steffen A dissertation submitted to the faculty of The University of Utah in partial fulfillment of the requirements for the degree of Doctor of Philosophy in Computing School of Computing The University of Utah December 2009 Copyright c Michael Dietel Steffen 2009 ° All Rights Reserved THE UNIVERSITY OF UTAH GRADUATE SCHOOL SUPERVISORY COMMITTEE APPROVAL of a dissertation submitted by Michael Dietel Steffen This dissertation has been read by each member of the following supervisory committee and by majority vote has been found to be satisfactory. Chair: Robert M. Kirby Martin Berzins Christopher R. Johnson Steven G. Parker James E. Guilkey THE UNIVERSITY OF UTAH GRADUATE SCHOOL FINAL READING APPROVAL To the Graduate Council of the University of Utah: I have read the dissertation of Michael Dietel Steffen in its final form and have found that (1) its format, citations, and bibliographic style are consistent and acceptable; (2) its illustrative materials including figures, tables, and charts are in place; and (3) the final manuscript is satisfactory to the Supervisory Committee and is ready for submission to The Graduate School. Date Robert M. Kirby Chair, Supervisory Committee Approved for the Major Department Martin Berzins Chair/Dean Approved for the Graduate Council Charles A. Wight Dean of The Graduate School ABSTRACT The Material Point Method (MPM) has shown itself to be a powerful tool in the simulation of large deformation problems, especially those involving complex geometries and contact where typical finite element type methods frequently fail. While these large complex problems lead to some impressive simulations and so- lutions, there has been a lack of basic analysis characterizing the errors present in the method, even on the simplest of problems.
    [Show full text]
  • Boundary Particle Method with High-Order Trefftz Functions
    Copyright © 2010 Tech Science Press CMC, vol.13, no.3, pp.201-217, 2010 Boundary Particle Method with High-Order Trefftz Functions Wen Chen1;2, Zhuo-Jia Fu1;3 and Qing-Hua Qin3 Abstract: This paper presents high-order Trefftz functions for some commonly used differential operators. These Trefftz functions are then used to construct boundary particle method for solving inhomogeneous problems with the boundary discretization only, i.e., no inner nodes and mesh are required in forming the final linear equation system. It should be mentioned that the presented Trefftz functions are nonsingular and avoids the singularity occurred in the fundamental solution and, in particular, have no problem-dependent parameter. Numerical experiments demonstrate the efficiency and accuracy of the present scheme in the solution of inhomogeneous problems. Keywords: High-order Trefftz functions, boundary particle method, inhomoge- neous problems, meshfree 1 Introduction Since the first paper on Trefftz method was presented by Trefftz (1926), its math- ematical theory was extensively studied by Herrera (1980) and many other re- searchers. In 1995 a special issue on Trefftz method, was published in the jour- nal of Advances in Engineering Software for celebrating its 70 years of develop- ment [Kamiya and Kita (1995)]. Qin (2000, 2005) presented an overview of the Trefftz finite element and its application in various engineering problems. The Trefftz method employs T-complete functions, which satisfies the governing dif- ferential operators and is widely applied to potential problems [Cheung, Jin and Zienkiewicz (1989)], two-dimensional elastic problems [Zielinski and Zienkiewicz (1985)], transient heat conduction [Jirousek and Qin (1996)], viscoelasticity prob- 1 Center for Numerical Simulation Software in Engineering and Sciences, Department of Engineer- ing Mechanics, Hohai University, Nanjing, Jiangsu, P.R.China 2 Corresponding author.
    [Show full text]
  • Scalable Large-Scale Fluid-Structure Interaction Solvers in the Uintah Framework Via Hybrid Task-Based Parallelism Algorithms
    1 Scalable Large-scale Fluid-structure Interaction Solvers in the Uintah Framework via Hybrid Task-based Parallelism Algorithms Qingyu Meng, Martin Berzins UUSCI-2012-004 Scientific Computing and Imaging Institute University of Utah Salt Lake City, UT 84112 USA May 7, 2012 Abstract: Uintah is a software framework that provides an environment for solving fluid-structure interaction problems on structured adaptive grids on large-scale science and engineering problems involving the solution of partial differential equations. Uintah uses a combination of fluid-flow solvers and particle-based methods for solids, together with adaptive meshing and a novel asynchronous task-based approach with fully automated load balancing. When applying Uintah to fluid-structure interaction problems with mesh refinement, the combination of adaptive meshing and the move- ment of structures through space present a formidable challenge in terms of achieving scalability on large-scale parallel computers. With core counts per socket continuing to grow along with the prospect of less memory per core, adopting a model that uses MPI to communicate between nodes and a shared memory model on-node is one approach to achieve scalability at full machine capacity on current and emerging large-scale systems. For this approach to be successful, it is necessary to design data-structures that large numbers of cores can simultaneously access without contention. These data structures and algorithms must also be designed to avoid the overhead involved with locks and other synchronization primitives when running on large number of cores per node, as contention for acquiring locks quickly becomes untenable. This scalability challenge is addressed here for Uintah, by the development of new hybrid runtime and scheduling algorithms combined with novel lockfree data structures, making it possible for Uintah to achieve excellent scalability for a challenging fluid-structure problem with mesh refinement on as many as 260K cores.
    [Show full text]
  • The Material-Point Method for Granular Materials
    Comput. Methods Appl. Mech. Engrg. 187 (2000) 529±541 www.elsevier.com/locate/cma The material-point method for granular materials S.G. Bardenhagen a, J.U. Brackbill b,*, D. Sulsky c a ESA-EA, MS P946 Los Alamos National Laboratory, Los Alamos, NM 87545, USA b T-3, MS B216 Los Alamos National Laboratory, Los Alamos, NM 87545, USA c Department of Mathematics and Statistics, University of New Mexico, Albuquerqe, NM 87131, USA Abstract A model for granular materials is presented that describes both the internal deformation of each granule and the interactions between grains. The model, which is based on the FLIP-material point, particle-in-cell method, solves continuum constitutive models for each grain. Interactions between grains are calculated with a contact algorithm that forbids interpenetration, but allows separation and sliding and rolling with friction. The particle-in-cell method eliminates the need for a separate contact detection step. The use of a common rest frame in the contact model yields a linear scaling of the computational cost with the number of grains. The properties of the model are illustrated by numerical solutions of sliding and rolling contacts, and for granular materials by a shear calculation. The results of numerical calculations demonstrate that contacts are modeled accurately for smooth granules whose shape is resolved by the computation mesh. Ó 2000 Elsevier Science S.A. All rights reserved. 1. Introduction Granular materials are large conglomerations of discrete macroscopic particles, which may slide against one another but not penetrate [9]. Like a liquid, granular materials can ¯ow to assume the shape of the container.
    [Show full text]
  • References (In Progress)
    R References (in progress) R–1 Appendix R: REFERENCES (IN PROGRESS) TABLE OF CONTENTS Page §R.1 Foreword ...................... R–3 §R.2 Reference Database .................. R–3 R–2 §R.2 REFERENCE DATABASE §R.1. Foreword Collected references for most Chapters (except those in progress) for books Advanced Finite Element Methods; master-elective and doctoral level, abbrv. AFEM Advanced Variational Methods in Mechanics; master-elective and doctoral level, abbrv. AVMM Fluid Structure Interaction; doctoral level, abbrv. FSI Introduction to Aerospace Structures; junior undergraduate level, abbrv. IAST Matrix Finite Element Methods in Statics; senior-elective and master level, abbrv. MFEMS Matrix Finite Element Methods in Dynamics; senior-elective and master level, abbrv. MFEMD Introduction to Finite Element Methods; senior-elective and master level, abbrv. IFEM Nonlinear Finite Element Methods; master-elective and doctoral level, abbrv. NFEM Margin letters are to facilitate sort; will be removed on completion. Note 1: Many books listed below are out of print. The advent of the Internet has meant that it is easier to surf for used books across the world without moving from your desk. There is a fast search “metaengine” for comparing prices at URL http://www.addall.com: click on the “search for used books” link. Amazon.com has also a search engine, which is poorly organized, confusing and full of unnecessary hype, but does link to online reviews. [Since about 2008, old scanned books posted online on Google are an additional potential source; free of charge if the useful pages happen to be displayed. Such files cannot be downloaded or printed.] Note 2.
    [Show full text]
  • And Smoothed Particle Hydrodynamics (SPH) Computational Techniques
    A strategy to couple the material point method (MPM) and smoothed particle hydrodynamics (SPH) computational techniques The MIT Faculty has made this article openly available. Please share how this access benefits you. Your story matters. Citation Raymond, Samuel J., Bruce Jones, and John R. Williams. “A Strategy to Couple the Material Point Method (MPM) and Smoothed Particle Hydrodynamics (SPH) Computational Techniques.” Computational Particle Mechanics 5, no. 1 (December 10, 2016): 49– 58. As Published http://dx.doi.org/10.1007/s40571-016-0149-9 Publisher Springer International Publishing Version Author's final manuscript Citable link http://hdl.handle.net/1721.1/116233 Terms of Use Creative Commons Attribution-Noncommercial-Share Alike Detailed Terms http://creativecommons.org/licenses/by-nc-sa/4.0/ Computational Particle Mechanics manuscript No. (will be inserted by the editor) A strategy to couple the Material Point Method (MPM) and Smoothed Particle Hydrodynamics (SPH) computational techniques Samuel J. Raymond · Bruce Jones · John R. Williams Received: date / Accepted: date Abstract A strategy is introduced to allow coupling of the Material Point Method (MPM) and Smoothed Particle Hydrodynamics (SPH) for numerical simulations. This new strategy partitions the domain into SPH and MPM regions, particles carry all state variables and as such no special treatment is required for the transition between regions. The aim of this work is to derive and validate the coupling methodology between MPM and SPH. Such coupling allows for general boundary conditions to be used in an SPH simulation without further augmentation. Additionally, as SPH is a purely particle method and MPM is a combination of particles and a mesh, this coupling also permits a smooth transition from particle methods to mesh methods.
    [Show full text]
  • A Self-Singularity-Capturing Scheme for Fractional Differential Equations
    A Self-Singularity-Capturing Scheme for Fractional Differential Equations Jorge L. Suzukia,b and Mohsen Zayernouria,c aDepartment of Mechanical Engineering, Michigan State University, East Lansing, MI, USA bDepartment of Computational Mathematics, Science and Engineering, Michigan State University, East Lansing, MI, USA cDepartment of Probability and Statistics, Michigan State University, East Lansing, MI, USA ARTICLE HISTORY Compiled April 30, 2020 ABSTRACT We develop a two-stage computational framework for robust and accurate time- integration of multi-term linear/nonlinear fractional differential equations. In the first stage, we formulate a self-singularity-capturing scheme, given avail- able/observable data for diminutive time, experimentally obtained or sampled from an approximate numerical solution utilizing a fine grid nearby the initial time. The fractional differential equation provides the necessary knowledge/insight on how the hidden singularity can bridge between the initial and the subsequent short-time so- lution data. In the second stage, we utilize the multi-singular behavior of solution in a variety of numerical methods, without resorting to making any ad-hoc/uneducated guesses for the solution singularities. Particularly, we employed an implicit finite- difference method, where the captured singularities, in the first stage, are taken into account through some Lubich-like correction terms, leading to an accuracy of order O(∆t3−α). We show that this novel framework can control the error even in the presence of strong multi-singularities.
    [Show full text]
  • Singular Boundary Method for Heat Conduction in Layered Materials
    Copyright © 2011 Tech Science Press CMC, vol.24, no.1, pp.1-14, 2011 Singular Boundary Method for Heat Conduction in Layered Materials H. Htike1;2, W. Chen1;2;3 and Y. Gu1;2 Abstract: In this paper, we investigate the application of the singular boundary method (SBM) to two-dimensional problems of steady-state heat conduction in isotropic bimaterials. A domain decomposition technique is employed where the bimaterial is decomposed into two subdomains, and in each subdomain, the solu- tion is approximated separately by an SBM-type expansion. The proposed method is tested and compared on several benchmark test problems, and its relative merits over the other boundary discretization methods, such as the method of fundamental solution (MFS) and the boundary element method (BEM), are also discussed. Keywords: Singular boundary method, origin intensity factor, heat conduction, bimaterials, meshless. 1 Introduction During machining processes, heat is a source that strongly influences the tool per- formance, such as tool wear, tool life, surface quality as well as process quality. The layers of material provide a barrier for the intensive heat flow from the con- tact area into the substrate material. Thus, a clear understanding of the temperature distribution in layered materials is very useful and important [Atluri (2005); Atluri, Liu, and Han (2006);Chang, Liu, and Chang (2010);Chen and Liu (2001)]. The method of fundamental solutions (MFS) [Karageorghis (1992);Fairweather and Karageorghis (1998);Chen, Golberg, and Hon (1998)] is a successful technique for the solution of some elliptic boundary value problems in engineering applica- tions. It is applicable when a fundamental solution of the operator of the governing equation is known.
    [Show full text]
  • Potential Problems by Singular Boundary Method Satisfying Moment Condition
    Copyright © 2009 Tech Science Press CMES, vol.54, no.1, pp.65-85, 2009 Potential Problems by Singular Boundary Method Satisfying Moment Condition Wen Chen1;2, Zhuojia Fu1 and Xing Wei1 Abstract: This study investigates the singular boundary method (SBM), a novel boundary-type meshless method, in the numerical solution of potential problems. Our finding is that the SBM can not obtain the correct solution in some tested cases, in particular, in the cases whose solution includes a constant term. To rem- edy this drawback, this paper presents an improved SBM formulation which is a linear sum of the fundamental solution adding in a constant term. It is stressed that this SBM approximation with the additional constant term has to satisfy the so-called moment condition in order to guarantees the uniqueness of the solution. The efficiency and accuracy of the present SBM scheme are demonstrated through detailed comparisons with the exact solution, the method of fundamental solutions and the regularized meshless method. Keywords: Singular boundary method, fundamental solution, singularity at the origin, moment condition, potential problem 1 Introduction Without a need of harassing mesh generation in the traditional FEM [Zienkiewicz and Taylor (1991)] and BEM [Ong and Lim (2005)], meshless methods [Li and Liu (2002)] have attracted much attention in recent years. This study focuses on the boundary-type meshless numerical techniques, for instance, method of funda- mental solutions (MFS) [Chen, Golberg and Hon (1998); Cheng, Young and Tsai (2000); Fairweather and Karageorghis (1998); Karageorghis (1992); Liu (2008); Marin (2008,2010a,2010b); Poullikkas, Karageorghis and Georgiou (1998); Reut- skiy (2005); Smyrlis and Karageorghis (2001)], boundary knot method (BKM) by Chen and Tanaka (2002), boundary particle method [Chen and Fu (2009): Fu, Chen 1 Center for Numerical Simulation Software in Engineering and Sciences, Department of Engineer- ing Mechanics, College of Civil Engineering, Hohai University, Nanjing, Jiangsu, P.R.China 2 Corresponding author.
    [Show full text]
  • Annual Report 2014
    ANNUAL REPORT 2014 www.svf.stuba.sk Table of Content 2014 ANNUAL REPOrt of the Faculty of Civil Engineering Bratislava © Faculty of Civil Engineering, Slovak University of Technology in Bratislava, 2015 Authorized contributions from the departments Introduction translation by Dagmar Špildová, Slovak University of Technology in Bratislava, Faculty of Civil Engineering, Department of Languages Language consultant: Debra Gambrill JD, Slovak University of Technology in Bratislava, Faculty of Civil Engineering, Department of Languages Editor: Jozef Urbánek Layout: Ivan Pokrývka TABLE OF CONTENT Introduction .....................................................................................................5 Bodies of The Faculty ........................................................................................8 Departments Department of Concrete Structures and Bridges ..................................................14 Department of Transportation Engineering .........................................................22 Department of Theoretical Geodesy ...................................................................26 Department of Surveying ..................................................................................33 Department of Geotechnics ..............................................................................39 Department of Land and Water Resources Management ........................................45 Department of Hydraulic Engineering .................................................................53 Department
    [Show full text]