|
[Sponsors] |
July 8, 2015, 23:16 |
Convection operator (Theoretical question)
|
#1 |
New Member
Join Date: Jul 2015
Posts: 7
Rep Power: 10 |
Hello,
I am trying to write convection term from Navier-Stokes equation in vector-matrix form. Let V,W be two velocity solutions (2-dimensional in my case). From my calculation I got that (V convection operator W) = Hessian(W)*V. Can somebody confirm/disprove whether above is true? Also want to make sure that convection operator is linear in both arguments. The last seems obvious but need to know for sure. Thank you. |
|
July 9, 2015, 03:27 |
|
#2 | |
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,769
Rep Power: 71 |
Quote:
I am not sure about the sense of your question ... independently from the number of dimensions, the convective operator is defined by the nabla operator along with a velocity vector field: Div (v * ) if you consider the momentum equation, the term is Div (v v) and is non linear. By exploiting the continuity equation you can recast it as v *Grad v. |
||
July 9, 2015, 20:30 |
|
#3 | |
New Member
Join Date: Jul 2015
Posts: 7
Rep Power: 10 |
Quote:
For example in weak formulation I have the term convection operator((3/2*Vel_{n+1}-1/2*Vel_n)*Vel_{n+1}, VelocityTestFunction) Vel_n denotes velocity solution at time step n. After I substitute Vel_n = \sum_i{c_i*VelocityBasisFunction_i} and VelocityTestFunction=VelocityBasisFunction_j into the above term, can I take sum outside convection operator outside convection operator for both 1st and 2nd arguments? Thats what I mean by 'linear in both operands'. If yes, than it becomes sum_i( sum_j(convection operator(VelocityBasisFunction_i, VelocityBasisFunction_j) )). Then to compute convection operator(VelocityBasisFunction_i, VelocityBasisFunction_j) I need to write it in the code in matrix-vector form using grad(VelocityBasisFunction_i), VelocityBasisFunction_j values etc. components (everything happens in cartesian coordinate system). |
||
July 10, 2015, 04:30 |
|
#4 | |
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,769
Rep Power: 71 |
Quote:
sorry, but I still don't get the focus point of the problem.... weak formulation means you have an integral form, therefore the convective operator is actually Int [S] n.v (*) dS S being the surface of the volume and (*) is the general function (scalar or vector) on which you apply the operator. You can use CN time integration on the weak form but you have to consider the discretization of the integral of the flux. Then you have to linearize the term to get a linear algebric system. However, very often the NS code use the convective term in the explicit Adam-Bashforth integration to avoid get into a non-linear system. The diffusion is discretized by the CN. |
||
July 11, 2015, 23:47 |
|
#5 | |
New Member
Join Date: Jul 2015
Posts: 7
Rep Power: 10 |
Quote:
Where can I read about this? |
||
July 12, 2015, 04:21 |
|
#6 | |
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,769
Rep Power: 71 |
Quote:
is a hystorical methods .... http://ntrs.nasa.gov/archive/nasa/ca...9840014260.pdf |
||
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Unanswered question | niklas | OpenFOAM | 2 | July 31, 2013 16:03 |
convection, diffusion, advection, dispersion | Lam | Main CFD Forum | 7 | May 18, 2007 12:04 |
heat convection in a forced fluid | Daniele | COMSOL | 0 | December 20, 2005 08:09 |
natural convection in a sealed enclosure | James | Main CFD Forum | 4 | April 2, 2001 15:48 |
question | K.L.Huang | Siemens | 1 | March 29, 2000 04:57 |