1 / 79

Null-field approach for Laplace problems with circular boundaries using degenerate kernels

This research study explores the use of the null-field approach and degenerate kernels to solve Laplace problems with circular boundaries. The mathematical formulation, expansions of fundamental solution and boundary density, adaptive observer system, vector decomposition technique, linear algebraic equation, and numerical examples are discussed. The study concludes with insights on the degenerate scale and overall conclusions.

marialawson
Download Presentation

Null-field approach for Laplace problems with circular boundaries using degenerate kernels

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript


  1. Null-field approach for Laplace problems with circular boundaries using degenerate kernels 研究生:沈文成 指導教授:陳正宗 教授 時間:10:30 ~ 12:00 地點:河工二館 307室 碩士論文口試

  2. Outlines • Motivation and literature review • Mathematical formulation • Expansions of fundamental solution and boundary density • Adaptive observer system • Vector decomposition technique • Linear algebraic equation • Numerical examples • Degenerate scale • Conclusions

  3. Outlines • Motivation and literature review • Mathematical formulation • Expansions of fundamental solution and boundary density • Adaptive observer system • Vector decomposition technique • Linear algebraic equation • Numerical examples • Degenerate scale • Conclusions

  4. Motivation and literature review BEM/BIEM Improper integral Singular and hypersingular Regular Fictitious BEM Bump contour Limit process Fictitious boundary Null-field approach CPV and HPV Collocation point Ill-posed

  5. Present approach Degenerate kernel Fundamental solution No principal value CPV and HPV Advantages of degenerate kernel 1. No principal value 2. Well-posed

  6. Engineering problem with arbitrary geometries Straight boundary Degenerate boundary (Chebyshev polynomial) (Legendre polynomial) Circular boundary (Fourier series) (Mathieu function) Elliptic boundary

  7. Motivation and literature review Analytical methods for solving Laplace problems with circular holes Special solution Bipolar coordinate Conformal mapping Chen and Weng, 2001, “Torsion of a circular compound bar with imperfect interface”, ASME Journal of Applied Mechanics Honein, Honein and Hermann, 1992, “On two circular inclusions in harmonic problem”, Quarterly of Applied Mathematics Lebedev, Skalskaya and Uyand, 1979, “Work problem in applied mathematics”, Dover Publications Limited to doubly connected domain

  8. Fourier series approximation • Ling (1943) - torsion of a circular tube • Caulk et al. (1983) - steady heat conduction with circular holes • Bird and Steele (1992) - harmonic and biharmonic problems with circular holes • Mogilevskaya et al. (2002) - elasticity problems with circular boundaries

  9. Contribution and goal • However, they didn’t employ the null-field integral equationand degenerate kernels to fully capture the circular boundary, although they all employed Fourier series expansion. • To develop a systematic approach for solving Laplace problems with multiple holes is our goal.

  10. Outlines • Motivation and literature review • Mathematical formulation • Expansions of fundamental solution and boundary density • Adaptive observer system • Vector decomposition technique • Linear algebraic equation • Numerical examples • Degenerate scale • Conclusions

  11. Boundary integral equation and null-field integral equation Interior case Exterior case Null-field integral equation

  12. Outlines • Motivation and literature review • Mathematical formulation • Expansions of fundamental solution and boundary density • Adaptive observer system • Vector decomposition technique • Linear algebraic equation • Numerical examples • Degenerate scale • Conclusions

  13. Expansions of fundamental solution and boundary density • Degenerate kernel - fundamental solution • Fourier series expansions - boundary density

  14. Separable form of fundamental solution (1D) Separable property continuous discontinuous

  15. Separable form of fundamental solution (2D)

  16. Boundary density discretization Fourier series Ex . constant element Present method Conventional BEM

  17. Outlines • Motivation and literature review • Mathematical formulation • Expansions of fundamental solution and boundary density • Adaptive observer system • Vector decomposition technique • Linear algebraic equation • Numerical examples • Degenerate scale • Conclusions

  18. collocation point Adaptive observer system

  19. Outlines • Motivation and literature review • Mathematical formulation • Expansions of fundamental solution and boundary density • Adaptive observer system • Vector decomposition technique • Linear algebraic equation • Numerical examples • Degenerate scale • Conclusions

  20. Vector decomposition technique for potential gradient True normal direction Non-concentric case: Special case (concentric case) :

  21. Outlines • Motivation and literature review • Mathematical formulation • Expansions of fundamental solution and boundary density • Adaptive observer system • Vector decomposition technique • Linear algebraic equation • Numerical examples • Degenerate scale • Conclusions

  22. Linear algebraic equation where Index of collocation circle Index of routing circle Column vector of Fourier coefficients (Nth routing circle)

  23. Explicit form of each submatrix [Upk] and vector {tk} Truncated terms of Fourier series Number of collocation points Fourier coefficients

  24. Flowchart of present method Potential gradient Vector decomposition Degenerate kernel Fourier series Adaptive observer system Potential of domain point Collocation point and matching B.C. Analytical Fourier coefficients Linear algebraic equation Numerical

  25. Comparisons of conventional BEM and present method

  26. Outlines • Motivation and literature review • Mathematical formulation • Expansions of fundamental solution and boundary density • Adaptive observer system • Vector decomposition technique • Linear algebraic equation • Numerical examples • Degenerate scale • Conclusions

  27. Numerical examples • Steady state heat conduction problems • Electrostatic potential of wires • Flow of an ideal fluid pass cylinders • A circular bar under torque • An infinite medium under antiplane shear • Half-plane problems

  28. Numerical examples • Steady state heat conduction problems • Electrostatic potential of wires • Flow of an ideal fluid pass cylinders • A circular bar under torque • An infinite medium under antiplane shear • Half-plane problems

  29. Steady state heat conduction problems Case 1 Case 2

  30. Steady state heat conduction problems Case 3 Case 4

  31. Case 1: Isothermal line Exact solution (Carrier and Pearson) FEM-ABAQUS (1854 elements) Present method (M=10) BEM-BEPO2D (N=21)

  32. Relative error of flux on the small circle

  33. Convergence test - Parseval’s sum for Fourier coefficients Parseval’s sum

  34. Case 2: Isothermal line Present method (M=10) Caulk’s data (1983) IMA Journal of Applied Mathematics FEM-ABAQUS (6502 elements)

  35. Case 3: Isothermal line Present method (M=10) Caulk’s data (1983) IMA Journal of Applied Mathematics FEM-ABAQUS (8050 elements)

  36. Case 4: Isothermal line Present method (M=10) Caulk’s data (1983) IMA Journal of Applied Mathematics FEM-ABAQUS (8050 elements)

  37. Numerical examples • Steady state heat conduction problems • Electrostatic potential of wires • Flow of an ideal fluid pass cylinders • A circular bar under torque • An infinite medium under antiplane shear • Half-plane problems

  38. Electrostatic potential of wires Two parallel cylinders held positive and negative potentials Hexagonal electrostatic potential

  39. Contour plot of potential Exact solution (Lebedev et al.) Present method (M=10)

  40. Contour plot of potential Onishi’s data (1991) Present method (M=10)

  41. Numerical examples • Steady state heat conduction problems • Electrostatic potential of wires • Flow of an ideal fluid pass cylinders • A circular bar under torque • An infinite medium under antiplane shear • Half-plane problems

  42. Flow of an ideal fluid pass two parallel cylinders is the velocity of flow far from the cylinders is the incident angle

  43. Velocity field in different incident angle Present method (M=10) Present method (M=10)

  44. Numerical examples • Steady state heat conduction problems • Electrostatic potential of wires • Flow of an ideal fluid pass cylinders • A circular bar under torque • An infinite medium under antiplane shear • Half-plane problems

  45. Torsion bar with circular holes removed The warping function Boundary condition where Torque on

  46. Axial displacement with two circular holes Dashed line: exact solution Solid line: first-order solution Caulk’s data (1983) ASME Journal of Applied Mechanics Present method (M=10)

  47. Axial displacement with three circular holes Dashed line: exact solution Solid line: first-order solution Caulk’s data (1983) ASME Journal of Applied Mechanics Present method (M=10)

  48. Axial displacement with four circular holes Dashed line: exact solution Solid line: first-order solution Caulk’s data (1983) ASME Journal of Applied Mechanics Present method (M=10)

  49. Numerical examples • Steady state heat conduction problems • Electrostatic potential of wires • Flow of an ideal fluid pass cylinders • A circular bar under torque • An infinite medium under antiplane shear • Half-plane problems

  50. Infinite medium under antiplane shear The displacement Boundary condition Total displacement on

More Related