Optimal Shape Determination of A Oscillating Body by ALE Finite Element Method

(2006

)

ALE

Optimal Shape Determination of A Oscillating Body by ALE Finite Element Method

41

^!"#$%

Hiroki YOSHIDA

1

&'(

)*+,-./

,

^{0123456,78/)*+9}

:;<=>?

,

^{@ABCD.EF9GH6/.<}

.

^IJ

,

)*+90123456@H:;<5KL4CDM?+

KNO,PQR/-B

,

^{5KLSTUFV<NO4WX}

?5KL,-;<YZ4[F\]^<

.

^_`4Tab

cde4fg1

,

^{hijk,-;<lmnopq4fg}

,

MH,?TUrst4fg,78/

,

^5KLTUrNO

?hiu,WXV<EF9vwFA8/.<

.

^R@Rx

yu,

,

^+KNO4z{?

,

^{5K,78/+K?|}S~}

EMA.F.€X4F‚ƒ6/DJ

.

^„…†4‡u

?

,

^TaANOS

,

^{5K,78/+K9|}V<F.ˆ}

‰4Š

,

^Tabcde,‹ŒD

,

^hiu,WXV<EF]

^<

.

5KLTUNOWXZŽ4B‘?

, Pironneau

1),2)

9

Stokes

^{_’4“”V<•–]—4jt4s˜F}

™šS›R

,

^{—4œhi5KL4fg,7B569}

bžŸ jk]D<EF,A8JEF,7B„¡rRJ

.

„…†]?

,

^{¢£¤h,5KLS¥¦§.<EF]‡u}

¤hFR

,

^¨©iSz{ªh

X i

^,ǤV<

.

^IJ

,

^Ta

bcde,‹Œ.J

Lagrange

^¬htS§.

,

^­®¢£¤

hSX¯V<

.

^°±²³´,?

, ALE

^{µ¶,7<·¸¹}

rM6Jº»¼

Navier-Stokes

^²³´F

,

^½¾xXˆ‰

4´S§.J

.

^{NO¿À,§.H6<Á±,?}

,

^°±²³

´SÂÃÄÅ_’V<Æ,TaNO4X´rS‚.

,

^z

{ªh

X i

^{4¥¦uA¤hFV<EFA {ÇS‚8J}

4

(

^ÈÉt

)

^S§./.<

.

^{TUr,?Á±tSa§R}

,

NO4¿ÀS‚8J

.

2

ÊËÌÍÎ

2.1 Navier-Stokes

ÏÐÑ

ÒÓ

Γ

^SÂV<

2

^¸¹2ÔÕÖ

Ω ∈ R ²

^SX¯R

,

^Õ

Ö

Ω

^{?º»¼×™5K]ØJM6/.<4FV<}

.

^°

±²³´,?

,

^ي4

ALE

^{µ¶,7<·¸¹rM6Jº}

»¼

Navier-Stokes

^²³´S§.<

.

˙

u i + b j u i,j + p ,i − ν (u i,j + u j,i ) ,j = 0 in Ω (1) u i,i = 0 in Ω (2)

EE]

, u i

^F

p

^{?5ÚF»L]^B}

, ν

^?}יÛhS›

V

.

^IJ

ν

^?

Reynolds

^h4Üh]µM6<

. b j

^?

ALE

¨©ÝÞßà,áV<šâ4ãáڞ]^<

.

2.2

äåæçèé

+K4NO?{Çê³,-./

,

^xëR/½¾xX4

ˆ‰SØìR

,

^{TarM6A;6íAHA.}

.

^EE]î

ï<+K4½¾9xX]^<F.EF?

,

^jkÕÖ,

-;<ÂÃÄÅðñÝo4Vï/4ÄÅ4½¾4òó9

,

xX]^<EFFô£]^<

.

^{½¾xXˆ‰4´?ي}

47,µM6<

.

m

X

e=1

a e (x i ) ^(l) − A 0 = 0 in Ω (3)

EE]

, a e

^{?õÄÅ4½¾SµR}

, A 0

^{?ö÷ÕÖ4½¾}

, m

^?òÄÅh

, l

^{?Tar,ÉøBùRúhSµV}

.

3

Îû

3.1

üýþÿ

„…†]?

,

^5KLbcZS/.<

.

^—E]

,

^b

cS‚5KLS¥¦¢£¤h,§.<

3) .

^¢£¤h

J

?5KL4¬ó,7B

,

^{ي47,X¯M6<}

.

J = 1 2

Z t _f

t 0

q 1 (F 1 − F ˆ 1 ) ² dt + 1 2

Z t _f

t 0

q 2 (F 2 − F ˆ 2 ) ² dt (4)

EE]

, q 1

^F

q 2

^?L

,

^{L,Û<‘¤hSR/}

.<

.

^MH,

, F 1

^F

F 2

^?L

,

^L4{ÇiS

, ˆ F 1

^F

F ˆ 2

^?L

,

^L4‡uiSR

,

^{ي4´,78/µ}

M6<

. F i = −

Z

Γ B

− pδ ij + ν (u i,j + u j,i )}n j

dΓ (5)

3.2

üýþÿ

Tabcde]?

,

^{X¯M6J¢£¤h9TUFA8J}

O9

,

^{TaAO]^<FX¯M6<EE]¢£¤h}

?

,

^°±²³´

(1),(2)

^F½¾b

(3)

^SØJRA9HT

UrM6A;6íAHA.E6?VAƒ

,

^ˆ‰

D4TUrZFH6<

.

^„…†]?

,

^m

¬htS§.

,

^{¢£¤hS­®V<EF]}

,

^ˆ‰AR

4TUrZ,D<EFFV<

.

^´

(1),(2),(3)

^S

IF/

, R k (U i , X j ) = 0

^Fµ¶R

,

^mo¬h

S

λ

^{FV<F­®¢£¤h?}

,

J ^∗ (U i , X j ) = J(U i , X j ) − λ k R k (U i , X j ) (6)

FA<

.

^{E4­®¢£¤h4xªSFB}

,

^ˆ‰S

<EF]Èɲ³´FNO¿À,§.<Á±S

V<

.

3.3

èé

¢£¤hSTUrV<J4ˆ‰?

,

^­®¢£¤

h

J ^∗

^4xªS

(δJ ^∗ = 0)

^{,V<EF]^<}

.

δJ ^∗ = ( ∂J

∂U i

− λ k

∂R k

∂U i

)δU i

+ ( ∂J

∂X j

− λ k

∂R k

∂X j

)δX j + R k δλ k (7)

ˆ‰7B´

(7)

^4õSF<F

,

^Èɲ³´

FÁ±?ي47,H6<

.

− u ˙ i ∗ − (b j u i ∗ ) ,j + p i ∗

− ν(u ^∗ _i,j + u ^∗ _j,i ) ,j = 0 in Ω (8)

u ^∗ _i,i = 0 in Ω (9)

u ^∗ _i = 0 on Γ U (10) s i = 0 on Γ D (11) s 1 = 0, u ^∗ ₂ = 0 on Γ S (12) u ^∗ _i − q i F i = 0 on Γ B (13) u ^∗ _i (t f ) = 0 in Ω (14)

∂J ^∗

∂X j

= grad(J ^∗ ) i = −s j u j,i (15)

EE]

, s i

^{?ي4´]µM6<}

. s i = {b j u ^∗ _i − p ^∗ δ ij + ν(u ^∗ _i,j + u ^∗ _j,i )}n j (16)

3.4

NO¿À,§.H6<Á±?

,

^郻h@HH6<

Ým

s j

^F5Ú42Ô²4i

u j,i

^,78

/{ÇM6<

.

^„

,

^Òӈ‰SV<F

,

^5K,

@6J+K4µ½5Ú?FAB

,

^—45Ú4

iFA8/RI

.

^—E]

, u j,i

^,áR/N

x¸_’Sa§R

,

^{+K4â,-;<5Ú42Ô²}

4iS

,

^{—4â,¤ƒ<ÄÅÔ]4iFR}

/<EFFV<

.

4

!

5K,@6J+K?"KF‘AR

,

^ي4|}²

³´,#

,

^ª$

Y

^“,7<

1

^%&žSÂR/.<F

€XV<

.

m Y ¨ + c Y ˙ + kY = F 2 (17)

EE]

, m, c, k

^?—6'6š(Ûh

,

^)*Ûh

,

^+™Ûh

SµR/.<

.

5

,û-

5.1

.

„…†]?

,

^{TUrstFR/ö÷i,PQRA.Á±}

tSa§RJ

.

^„st]?

,

^­®¢£¤h

J ^∗

^,áR//

01ß2pS3V<EF,7B

,

^45¢£¤h

K

⁹

@6<

. K ^(l) = J ^{∗ (l)}

+ 1 2

X _i ^(l+1) − X _i ^(l) T

W

X _i ^(l+1) − X _i ^(l)

(18)

EE]

, W

^{?‘Ûh‚6SµR/.<}

.

^45¢£¤h

9TU,A<FD

,

^{3RJ/01ß2p9,A}

<J

,

^{45¢£¤h9TU,A<EFF}

,

^­®¢£¤h

9TU,A<EF?ô£]^<EF,A<

. X _i ^(l)

^STa

jFV<F

, δK ^(l) = 0

^9“B”\J

,

^NO¿À47

ðpq?ي4´]µM6<

. X _i ^(l+1) = X _i ^(l) − 1

W grad(J ^∗ ) ^(l) (19)

6

89û

„…†]?

,

^°±²³´

·

^{Èɲ³´F,:;4<=}

rS‚

.

^{2Ô²4<=rFR/?}

,

^>Xr3?¤hÂ

ÃÄÅtSa§R

,

^{5Ú,¤R/?3?¤hÄÅ}

,

^»L,

¤R/?x¸ÄÅS§./.<

.

@Ô²4<=r,?

,

^>X™,A6

,

^@ÔBSCD

 6<

Crank-Nicolson

^t

(

^Cujt

)

^Sa§V<

.

7

DEFGHIJ

NO¿À,§.H6<Á±4K?

,

^«§V<ÂÃÄ

ÅðñÝo,PQV<EF9@8/.<

.

^L

.1

^47

,+K_M9NxAÂÃÄÅðñÝo

(Delaunay

^ðñ

Ýo

)

^S§.<F

,

^Á±4K?L

.2

^47,NOPA

4FA8/RI.

,

^{NO¿ÀS‚.P;<FNO9QR}

,

^{ðñÝo4ST15KL4BCAU4VW9~D/R}

I

.

^{—E]„…†]?L}

.3

^,›V7,

,

^+K_M,

X45YðñÝoS±RJÂÃÄÅðñÝoS«§V

<EFS<

.

^{E4ðñÝoS«§V<EF]Á±4}

K?L

.4

^{47,ZH@A49H6}

,

^ðñÝo4

ST1

,

^{5KL4BCF.8JZSjWRJ}

.

[

.1 Delaunay

^{\]^_[}

.2

^`a

[

.3

^{\]^_[}

.4

^`a

8

hijkFR/

,

^1h

250

^4º»¼

Navier-

Stokes

^{56,-;<”+K45KLTUNOWXZ}

S

.

^EE]”+K?L,7B

,

^U|}S~E

V4FV<

.

^„hijk]?

,

^NOSªNM<EF

]LFLSbcV<EF9‡u]^<J

,

^‘¤

h

q 1 , q 2

^SF,

1.0

^FzXRJ

.

^‘á‚6?

,

^NO

ªN,Éâ4}(9

,

^TC]

0.02

^S<EF

4A.7,<

.

^E6?

,

^ANO4},7<

,

jk,§.H6<ÂÃÄÅðñÝo4STSEFF

,

,NOSªNM<EF,78/

,

^,TaNO,

_Œ;<F.L,7<4]^<

.

^ö÷4NO,-

;<jkÕÖFÒӈ‰

,

^{ÂÃÄÅðñÝoSL}

.5, 6

^,

›V

.

(15.0, 5.0)

(15.0, -5.0) (-15.0, -5.0)

4.5D

U D

14.5D

u = 0 v = 0

u = 0, t = 0

u = 0, t = 0 14.5D

(-15.0, 5.0)

u = 1 1 t = 0

2 1

2 1

i

[

.5

^!"#

[

.6

^{$%&'\]^_}

8.1

()*+,-

jk,§.H6<ðñÝo?

,

^"K4|},#8/}

M6A;6íAHA.

.

^—E]

,

^´

(17)

^,78/H6

<.¥²

Y

^4ª$S

Laplace

^{²³´,/0V<EF}

]

,

^{â4}(S{ÇV<EFFV<}

.

∇ · φ = 0 in Ω (20)

φ = Y on Γ B (21)

φ = 0 on Γ U (22)

φ = 0 on Γ S (23)

φ = 0 on Γ D (24)

R

,

^{+Kµ½4â4}@H}

,

^ÕÖ34õâ4}

SÇV<EFFV<

.

^{+K4},É8/4â9}

}S‚4]

,

^{°±²³´S{ÇV<˜jk,-./?}

,

âFÄÅ4òh?ªrRA.

.

8.2

NO¿À,§.H6<Á±?

,

^郻h@HH6<

Ým

s j

^F5Ú42Ô²4i

u j,i

^,78

/{ÇM6<

.

^„

,

^Òӈ‰SV<F

,

^5K,

@6J+K4µ½5Ú?FAB

,

^—45Ú4

iFA8/RI

.

^—E]

, u j,i

^,áR/N

x¸_’Sي47,a§V<

. u j,i = ∂u j

∂X i

= Φ 1,i u j1 + Φ 2,i u j2 + Φ 3,i u j3 (25)

´

(47)

^,7B

,

^{+K4â,-;<5Ú42Ô²4}

iS

,

^{—4â,¤ƒ<ÄÅÔ]4iFR/}

<EFFV<

.

9

56

hijk78SL

.7

^@HL

.14

^,›V

.

^„…†]H6

JT9NO4­CLSL

.7,

^{NOTar,É¢£¤h4}

:O;SL

.8,

^{ö÷NOFTaNO,-;<LFL}

4@<=SL

.9,10,

^.¥²

Y

^4|}4@<=SL

.11

,—6'6›V

.

^L

.12

^{?ö÷NOF|}SRA.}

>,H6JNO

,

^{|}SRJ>,H6JNO}

S?@RJL]^<

.

^L

.9,10,11

^,-./

,

^âA?ö÷N

O,-;<õiS

,

^{A?T9NO,-;<õiSµR}

/.<

.

^MH,

,

^{õNO]4·¸¹@Ô}

100

^,-;<»

LKLSL

.13,14

^,›V

.

^{¢£¤h4:O;@H}

,

H6JNO9L

·

^{LSF,TUFR/.<EF9}

ƒ@B

,

^+K4|}

,

^{+KBB,fCV<»LDSEH}

6/.<EF9ƒ@<

.

^{MH,|}SRA.>,}

?ï

,

^{+KƲ9D,ªNV<F.FG9HH6}

<

.

^·¸¹@Ô

50

^@H

100

^{4•–]4ö÷NO@,-;}

<L

,

^L

,

^{+K4|}4npi,?ï}

,

^H6JN

O]4õi?

66.2%, 92.8%, 84.8%

^—6'6E<EF

9J

.

^xyu,

,

^{5K,^<+K4TaNO?5A}

I]^<EF9JH6/.<9

,

^H6JT9NO?

,

⁺

KƲ7Bœ²4K9L‘SM8J4FA8J

.

[

.11

^NOPQRS[

0.4 0.5 0.6 0.7 0.8 0.9 1

0 100 200 300 400 500

Iteration cycle peformance function

[

.7

^Q

0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1

0 20 40 60 80 100

Force

Non-dimensional time final shape initial shape

[

.8

-0.6 -0.4 -0.2 0 0.2 0.4 0.6

0 20 40 60 80 100

Force

Non-dimensional time final shape initial shape

[

.9

-0.04 -0.02 0 0.02 0.04

0 20 40 60 80 100

Displacement

Non-dimensional time final shape initial shape

[

.10

-1 -0.5 0 0.5 1

-1.5 -1 -0.5 0 0.5 1 1.5

Y-coordinate

X-coordinate

final shape (oscillating body) final shape (non-oscillating body) initial shape

[

.12

^PQNOPQ[

[

.13

^PQ[

(T = 100)

[

.14

^NOPQ[

(T = 100)

10

!"#$

,

^%&'

Navier-Stokes

^()*+,-

, ALE

^.

/0123456)789:;*<=->?@ABCD

EFG3HI7

.

^{JKLMNO>?PQ*RS-TU$}

,

VWXY*<Z

Lagrange

^[Y23AR=-\O#

,

^{]^}

VWXY3_`Z

,

^{aEbcdeEfghijkl-\}

Omno7

.

^{\)k$@A>?3l-7p*.qrs2}

#tI7OSu-

.

^{vw*!"#lk)7xy3z{}

=-

.

•

^|}~€Y

250

^{E‚rƒ„jr()*+S}

,

⁽

;†3‡u-@A>?$(ˆ‰#$rŠ

,

^‹L*Œ

3ŽI7>?m@ArEOr-

.

•

^{893‘ZrS’“O‘Z7’“#$}

,

^Ӡ*

•–r—m˜™

,

^{lk)->?*aEš›m•–}

Šœ)-\OmjI7

.

žŸ ¡

1. O. Pironneau, On Optimum Profiles in Stokes Flow. J. Fluid Mech, 59 , No.1, 117–128, 1973.

2. J. Matsumoto, T. Umetsu and M. Kawahara, In- compressible Viscous Flow Analysis and Adaptive Finite Element Method Using Linear Bubble Func- tion, J. Appl. Mech, 2 , 223-232, 1999.

3. A. Maruoka and M. Kawahara, Optimal Control in Navier-Stokes Equation, Int. J. Comp. Fluid Dyn, 9 , 313-322, 1998.

4. E. Kuhl, S. Hulshoff and R. de Borst, An Arbitrary Lagrangian Eulerian Finite-Element Approach for Fluid-Structure Interaction Phenomena, Int. J.

Numer. Meth. Eng 2003, 57 :117-142