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

Implicit discretization of Jacobian term in FV: Construction of Coefficient matrix

Register Blogs Community New Posts Updated Threads Search

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   January 19, 2017, 07:31
Default Implicit discretization of Jacobian term in FV: Construction of Coefficient matrix
  #1
Senior Member
 
Santiago Lopez Castano
Join Date: Nov 2012
Posts: 354
Rep Power: 15
Santiago is on a distinguished road
I am currently trying to write the linear system coefficients of the following system:

\frac{d\widehat{\mathbf{u}}}{\Delta t} - \underbrace{\frac{1}{|\Omega|}\int_{\partial \Omega}\frac{\partial \, d\widehat{\mathbf{u}}}{\partial \mathbf{n}}\,d\sigma}_{J(x)} = R(\mathbf{u})

But I cannot get my head around on how to come up with the coefficients for the jacobian fluxes. I mean, I see a tensor that needs to be projected along $\mathbf{n}$ and from it you do the discretization. My attempt on this enterprise is (the simplest, I think):

\left(\frac{\partial \, d\widehat{\mathbf{u}}}{\partial \mathbf{n}}\right)_e = \frac{d\widehat{\mathbf{u}}_P-d\widehat{\mathbf{u}}_E}{(\mathbf{x}_P-\mathbf{x}_E)\cdot\mathbf{n}}

What makes me uneasy is the fact that d\widehat{\mathbf{u}} is a vector. Thus, Are the coefficients for each of the vector's components equal to the face-normal projected distance between the neighboring centroids?
Santiago is offline   Reply With Quote

Old   January 19, 2017, 07:45
Default
  #2
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,849
Rep Power: 73
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
I do not understand your symbolism...

1) the first term in the LHS is the time derivative of u??
2) Since the first term is a vector, the second term in the LHS is the integral of a tensor (uu??) projected along the normal unit vector to the surface n.

For example, Int [S] n.uu dSIn your equation the symbol du/dn is still a tensor (why is d(du)?). The Jabian of the flux can be derived from the differential form starting from df/dxi = (df/dui)*(dui/dxi)
FMDenaro is offline   Reply With Quote

Old   January 19, 2017, 08:55
Default
  #3
Senior Member
 
Santiago Lopez Castano
Join Date: Nov 2012
Posts: 354
Rep Power: 15
Santiago is on a distinguished road
Quote:
Originally Posted by FMDenaro View Post
I do not understand your symbolism...

1) the first term in the LHS is the time derivative of u??
2) Since the first term is a vector, the second term in the LHS is the integral of a tensor (uu??) projected along the normal unit vector to the surface n.

For example, Int [S] n.uu dSIn your equation the symbol du/dn is still a tensor (why is d(du)?). The Jabian of the flux can be derived from the differential form starting from df/dxi = (df/dui)*(dui/dxi)
1.) No it's not the time derivative of \hat{u}, d \hat{u} it's just a vector. Besides the equation is already discrete in time.

2.) The the second term is just a diffusion transformed in a flux using the green-gauss theorem.

What I want is to get the coefficients of A, for the system:

[A][d\hat{\mathbf{u}}] = [R(\mathbf{u})]

Part of the diagonal terms are given by the first term of LHS, but my question is how to discretize the second term.

Thanks.
Santiago is offline   Reply With Quote

Old   January 19, 2017, 09:28
Default
  #4
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,849
Rep Power: 73
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
thus, using Gauss you have a surface integral like

Int [S] n. grad u dS

why do you need to introduce a Jacobian? you can directly discretize the integral of the fluxes
FMDenaro is offline   Reply With Quote

Old   January 20, 2017, 05:35
Default
  #5
Senior Member
 
Santiago Lopez Castano
Join Date: Nov 2012
Posts: 354
Rep Power: 15
Santiago is on a distinguished road
I referred the derivative of the vector function as jacobian. My bad.

Sent from my GT-I8190L using CFD Online Forum mobile app
Santiago is offline   Reply With Quote

Reply

Tags
finite volume method, linear system


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 Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Question about heat transfer coefficient setting for CFX Anna Tian CFX 1 June 16, 2013 07:28
ATTENTION! Reliability problems in CFX 5.7 Joseph CFX 14 April 20, 2010 16:45
FVM discretization of diffusion term on crvlnr gr Michail Main CFD Forum 3 March 14, 2008 07:52
Convection Term Discretization Maciej Matyka Main CFD Forum 3 October 7, 2004 05:32
Coefficient matrix Question. Main CFD Forum 1 September 12, 2004 22:19


All times are GMT -4. The time now is 23:50.