CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Main CFD Forum

finite element method for the Shallow Water Equations help

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Display Modes
Old   December 4, 2011, 22:28
Default finite element method for the Shallow Water Equations help
  #1
New Member
 
John Smith
Join Date: Dec 2011
Posts: 4
Rep Power: 5
alvesker is on a distinguished road
Hello,

I try to solve the Shallow Water Equations using finite element methods. Can anyone explain me how to treat nonlinear terms?

The problem is that we can write linear terms as matrix-vector multiplication where matrix is independent on time. For example for the height field h we have

h=\sum_{i=1,N}c_i^hB_i

after getting the Galerkin equation by multiplying on test function and integrating over domain we will get the term like

\sum_{i=1,N}c_i^h\int B_iB_jdx

or in the matrix form

(\int B_iB_jdx)_{i,j=1,N}(c_i^h)_{i=1,N}

Hence we compute the matrix only once and then multiply by updated coefficient vector

But how to deal with the terms like $uh$? where u is the velocity

uh=\sum_{i,k=1,N}c_i^uc_k^hB_iB_k

after getting the Galerkin equation this will be like

\sum_{i,k=1,N}c_i^uc_k^h\int B_iB_jB_kdx

How these can be represented in the matrix form in order not to compute all these integrals on each timestep?

Last edited by alvesker; December 4, 2011 at 22:54. Reason: typo
alvesker is offline   Reply With Quote

Old   December 5, 2011, 03:31
Default
  #2
Member
 
seyedashraf's Avatar
 
Omid Seyedashraf
Join Date: May 2010
Location: Iran
Posts: 41
Rep Power: 7
seyedashraf is on a distinguished road
Send a message via AIM to seyedashraf Send a message via Yahoo to seyedashraf
Dear alvesker,

you don't have to write the FEM approximated form of each parameter
is more accurate, though, I would suggest writing the the mean value of those parameters to simplify the weak form.

about the integration on each time step
it depends on the grid form you are using to split the domain.

Last edited by seyedashraf; December 5, 2011 at 05:52.
seyedashraf is offline   Reply With Quote

Old   December 5, 2011, 05:48
Default
  #3
New Member
 
John Smith
Join Date: Dec 2011
Posts: 4
Rep Power: 5
alvesker is on a distinguished road
Dear seyedashraf,

Thanks for the quick reply.

Can you explain what do you understand by the mean value?
is it like the following?

hu = \sum_{i=1,N}c_i^hc_i^u B_i

and how write it to the terms like

u\nabla h

where u, h are the velocity component and height field

u = \sum_{i=1,N}c_i^uB_i,\ h=\sum_{i=1,N}c_i^hB_i

So if apply the FEM we have

u\nabla h=\sum_{i=1,N}c_i^uB_i\sum_{i=1,N}c_i^h\nabla B_i=\sum_{i,j=1,N}c_i^uc_j^hB_i\nabla B_j

what is the average in the case when in the multiplication one of unknows one contain gradient (we have these in the advection formulation)

Last edited by alvesker; December 6, 2011 at 00:20. Reason: add
alvesker is offline   Reply With Quote

Old   December 5, 2011, 06:01
Default
  #4
New Member
 
John Smith
Join Date: Dec 2011
Posts: 4
Rep Power: 5
alvesker is on a distinguished road
the grid form - I use triangular grid
alvesker is offline   Reply With Quote

Reply

Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Lattice Boltzmann method vs Finite Element Method and Finite Volume Method solemnpriest Main CFD Forum 3 August 12, 2013 11:00
OpenFOAM-extend, Finite Element Method (FEM) Martin80 OpenFOAM Programming & Development 1 December 3, 2011 07:28
finite element method harten Main CFD Forum 0 November 1, 2007 06:37
BCs of hyperbolic Shallow Water Equation Kang, S. K. Main CFD Forum 4 November 21, 2000 22:20
What is C.V. based finite element method C-H Kuo Main CFD Forum 3 November 5, 1998 10:07


All times are GMT -4. The time now is 16:56.