760 likes | 773 Views
Introducing a semi-analytical method utilizing Fourier series for plate and Stokes flow problems with circular boundaries. The research includes direct and indirect boundary integral equation methods, numerical examples, and conclusions for practical engineering applications.
E N D
A semi-analytical approach for Stokes flow and plate problems with circular boundaries 報 告 者:蕭嘉俊 指導教授:陳正宗 博士 呂學育 博士 日 期:2005/6/16 地 點:河工二館307
Outlines • Introduction • Direct boundary integral equation method • Indirect boundary integral equation method • Numerical examples • Conclusions
Outlines • Introduction • Direct boundary integral equation method • Indirect boundary integral equation method • Numerical examples • Conclusions
Engineering problems with arbitrary boundaries Degenerate boundary (Chebyshev polynomials) Elliptic boundary (Mathieu function) Circular boundary (Fourier series) Straight boundary (Legendre polynomials) Degenerate boundary Circular boundary
Motivation BEM/BIEM Improper integral Desingular (Regular) Singular Indirect (Interior) Direct Null-field approach Fictitious boundary method Contour Limiting process : collocation point Fictitious boundary
Degenerate kernel Field point Present approach Motivation BEM/BIEM Improper integral Desingular (Regular) Singular Indirect (Exterior) Direct Null-field Contour Limiting process Fictitious boundary method : collocation point Fictitious boundary ill-posed CPV & HPV
Literature review Torsion bar with circular holes Laplace problems Steady state heat conduction of tube (Shen) Electromagnetic wave Engineering problems Helmholtz problems Membrane vibration Water wave and Acoustic problems (Chen) Plane elasticity:Airy stress function Biharmonic problems Solid mechanics:plate problem Fluid mechanics:Stokes flow (Hsiao)
Literature review • Plane elasticity: Jeffery (1921), Howland and Knight (1939), Green (1940) and Ling (1948) 2.Solid mechanics (Plate problem): Bird and Steele (1991) 3.Viscous flow (Stokes Flow): Kamal (1966), DiPrima and Stuart (1972), Mills (1977) and Ingham and Kelmanson (1984)
Purpose • A semi-analytical approach in conjunction with Fourier series, degenerate kernels and adaptive observer system is extended to biharmonic problems. • Advantages: 1. Mesh free. 2. Accurate. 3. Free of CPV and HPV.
Outlines • Introduction • Direct boundary integral equation method • Indirect boundary integral equation method • Numerical examples • Conclusions • Further research
Problem statement Governing equation: Essential boundary condition: :slope :lateral displacement, Natural boundary condition: : moment, : shear force
Direct boundary integral equations • BIEs are derived from the Rayleigh-Green identity: BIE for the domain point Null-field integral equation Interior problem
Slope Displacement : Poisson ratio Moment Displacement Shear force Displacement Boundary integral equation for the domain point
Slope Displacement : Poisson ratio Moment Displacement Shear force Displacement Null-field integral equation
Relation among the kernels Continuous (Separable form of degenerate kernel) is the fundamental solution, which satisfies
r Degenerate kernels x x s O
Fourier series The boundary densities are expanded in terms of Fourier series: M: truncating terms of Fourier series
Adaptive observer system : Collocation point : Radius of the jth circle : Origin of the jth circle : Boundary of the jth circle
Vector decomposition for normal derivative True normal direction Tangential direction Radial direction : normal derivative : tangential derivative
Linear algebraic system Null-field integral equations for and formulations H: number of circular boundaries Collocation circle index Routing circle index
Fourier series Degenerate kernels Flowchart of the present method Analytical Collocation method Matching B.C. Adaptive observer system Numerical Linear algebraic system Potential Fourier coefficients BIE for domain point
Stokes flow problems (Eccentric case) Governing equation: Essential boundary condition: on on (Stationary) : stream function : normal derivative of stream function
Linear algebraic system Given Unknown Unknown constant
Constraint equation Vorticity: Constraint:
Trapezoid integral Inner circle Outer circle Vector decomposition Numerical Analytical Trapezoid integral Series sum
Outlines • Introduction • Direct boundary integral equation method • Indirect boundary integral equation method • Numerical examples • Conclusions • Further research
Indirect boundary integral equation Indirect boundary integral equation is originated from the physical concept of superposition : Vorticity : single layer fictitious densities : double layer fictitious densities
Stokes flow problems (Eccentric case) Governing equation: Essential boundary condition: on on (Stationary) : stream function : normal derivative of stream function
Linear algebraic system Unknown Given Unknown constant :Collocation point
Constraint equation Vorticity: Constraint:
Trapezoid integral Inner circle Outer circle Vector decomposition Numerical Analytical Trapezoid integral Series sum
Outlines • Introduction • Direct boundary integral equation method • Indirect boundary integral equation method • Numerical examples • Conclusions
Plate problems (Case 1) Geometric data: and Essential boundary conditions: and on and on Exact solution:
Contour plot of displacement (No. of nodes=1,920, No. of elements=3,600) FEM mesh Exact solution FEM (ABAQUS) Present method (M=10)
Boundary densities for outer circle Exact solution: Exact solution:
Boundary densities for inner circle Exact solution: Exact solution:
Plate problems (Case 2) Geometric data: Essential boundary conditions: on and on and on and on and (Bird & Steele, 1991)
Contour plot of displacement Present method (N=41) Present method (N=21) Present method (N=61) Present method (N=81)
Contour plot of displacement Present method (N=101) Bird and Steele (1991) (No. of nodes=3,462, No. of elements=6,606) FEM mesh FEM (ABAQUS)
Stokes flow problems (Case 1) Boundary conditions: Exact solution: (Mills, 1977) where
Contour plot of Streamline Exact solution (Mills, 1977) Present method (N=161) Analytical solution, P=100 (Wu, 2004) (Null-field + collocation) Exact solution: (Mills, 1977) (Closed-form solution) (Null-field equation solution) (Wu, 2004) Analytical solution: (Trefftz solution)
Stokes flow problems (Case 2) Governing equation: Angular velocity: Boundary conditions: on and (Stationary) on and Eccentricity:
Comparison of stream function n: number of boundary nodes N: number of collocation points
BIE (Kelmanson) Present method Analytical solution Comparison for (160) (28) u1 (320) (640) (36) (∞) (44) DOF of BIE (Kelmanson) DOF of present method