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

Question about the DG implementation

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

Reply
 
LinkBack Thread Tools Display Modes
Old   July 4, 2012, 10:20
Default Question about the DG implementation
  #1
New Member
 
Elmiloud Chaabelasri
Join Date: Oct 2010
Posts: 25
Rep Power: 6
phy4me is on a distinguished road
Hello
I am asking your help because I'm stuck and I can't find help.
In fact, I try to develop a numerical code for solving shallow water equations (initially 1D) with the discontinuous Galerkin method. The first step that i can start it is construction of the discretized scheme by following the steps found everywhere in the articles (ideally Cockborne and all.). The physical states are approached by the following formula:

Uh = \sum U_h,i Ph_i

I understand at the moment that it need to inject this formula into the discretized form found above.In the case of i = 3 we have three equations to be solved to derive U_h,i.

My question: how i can go back to calculate the U_h
__________________
----
Chaabelasri
CFD Research staff
http://chaabelasri.blogspot.com/
phy4me is offline   Reply With Quote

Old   July 4, 2012, 13:12
Default
  #2
Senior Member
 
cfdnewbie
Join Date: Mar 2010
Posts: 551
Rep Power: 11
cfdnewbie is on a distinguished road
here's how you do it:
1) multiply your equation (the PDE) by the test function
2) integrate over the domain/cell
3) use integration by parts to separate volume and boundary terms
4) insert your expansion into the equation
5) that will give you a system to solve for your degrees of freedom u_h,i
6) evaluate the sum (the expansion) per element to get U_h
cfdnewbie is offline   Reply With Quote

Old   July 4, 2012, 13:53
Default
  #3
New Member
 
Elmiloud Chaabelasri
Join Date: Oct 2010
Posts: 25
Rep Power: 6
phy4me is on a distinguished road
Thank you for your response.
Indeed, this is the sixth point that interests me, i must evaluate the degree of freedom for use it in the sum. And for each item of the sum, i must raplace the function test by his value. It's that required ?
__________________
----
Chaabelasri
CFD Research staff
http://chaabelasri.blogspot.com/
phy4me is offline   Reply With Quote

Old   July 4, 2012, 13:59
Default
  #4
Senior Member
 
cfdnewbie
Join Date: Mar 2010
Posts: 551
Rep Power: 11
cfdnewbie is on a distinguished road
Quote:
Originally Posted by phy4me View Post
Thank you for your response.
Indeed, this is the sixth point that interests me, i must evaluate the degree of freedom for use it in the sum. And for each item of the sum, i must raplace the function test by his value. It's that required ?
I am not sure I understand what you are asking.... You are solving for the degrees of freedom u, the basis functions are something you know... So you multily your basis by your dof, add them up and evaluate your polynomial where you need it...

Hope that helps!
cfdnewbie is offline   Reply With Quote

Old   July 4, 2012, 14:35
Default
  #5
New Member
 
Elmiloud Chaabelasri
Join Date: Oct 2010
Posts: 25
Rep Power: 6
phy4me is on a distinguished road
In fact, what I meant is that after calculating the degrees of freedom we must return to the above formula of the sum to have the value of U in each element. I try to validate this point, because in some code I see U_h = U_h,n[n = max (i)] the last degrees of freedom
__________________
----
Chaabelasri
CFD Research staff
http://chaabelasri.blogspot.com/
phy4me is offline   Reply With Quote

Old   July 4, 2012, 15:01
Default
  #6
Senior Member
 
cfdnewbie
Join Date: Mar 2010
Posts: 551
Rep Power: 11
cfdnewbie is on a distinguished road
Quote:
Originally Posted by phy4me View Post
In fact, what I meant is that after calculating the degrees of freedom we must return to the above formula of the sum to have the value of U in each element. I try to validate this point, because in some code I see U_h = U_h,n[n = max (i)] the last degrees of freedom
Sorry, i don't understand the equation you posted...
To get U in your cell, just evaluate the (full) basis with the dof you computed....btw, are you referring to a modal or nodal scheme?
cfdnewbie 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
Question Re Engineering Data Source imnull ANSYS 0 March 5, 2012 14:51
internal field question - PitzDaily Case atareen64 OpenFOAM Running, Solving & CFD 2 January 26, 2011 16:26
Question about Table applicaiton. universez OpenFOAM Running, Solving & CFD 0 January 12, 2010 21:31
Please explain the implementation of species transport Eqn in reactingFoam kallipygian OpenFOAM Running, Solving & CFD 0 October 13, 2008 07:29
Poisson Solver question Suresh Main CFD Forum 3 August 12, 2005 04:37


All times are GMT -4. The time now is 07:22.