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

FVM Source term discretisation

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   November 29, 2015, 07:14
Default FVM Source term discretisation
  #1
Member
 
Victor Koppejan
Join Date: May 2015
Posts: 40
Rep Power: 11
vkoppejan is on a distinguished road
Hi Everyone

I'm looking for some literature references on FVM source term linearisation.

I have Versteeg and Patankar, they give a nice description, but I'm looking for some material which incoorporates langmuir (or orther non-linear type)adsorption istherms in coupled PDE's (coupled in at the PDE level, not block coupled solver).

Thanks in advance for your help.

Victor
vkoppejan is offline   Reply With Quote

Old   November 29, 2015, 08:04
Default
  #2
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,893
Rep Power: 73
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
Quote:
Originally Posted by vkoppejan View Post
Hi Everyone

I'm looking for some literature references on FVM source term linearisation.

I have Versteeg and Patankar, they give a nice description, but I'm looking for some material which incoorporates langmuir (or orther non-linear type)adsorption istherms in coupled PDE's (coupled in at the PDE level, not block coupled solver).

Thanks in advance for your help.

Victor

Lagmuir circulation in ocean is governed by the Craik-Leibovich equations model which, at the best of my knowledge, appear in differential form. I am not aware of numerical solutions based on the integral form but you can derive it by integrating the equations over a Finite Volume and using a proprer discretization for the volume integral
FMDenaro is offline   Reply With Quote

Old   November 29, 2015, 08:08
Default
  #3
Member
 
Victor Koppejan
Join Date: May 2015
Posts: 40
Rep Power: 11
vkoppejan is on a distinguished road
Hi Filipo,

Thanks for your reply. I think there's a misunderstanding here.

I was mentioning Langmuir isotherms for adsorption of a chemical species. They become non-linear at higher species concentrations.

Victor
vkoppejan is offline   Reply With Quote

Old   November 29, 2015, 08:13
Default
  #4
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,893
Rep Power: 73
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
hello, sorry for that...
however the answer is the same, you have first to integrate the term over a Finite Volume. If it is a local production term, it does not appear in divergence form and you cannot use Gauss to produce a surface integral.
That means you have to discretize the volume integral, for example as shown in the book of Peric & Ferziger
FMDenaro is offline   Reply With Quote

Old   November 29, 2015, 08:18
Default
  #5
Member
 
Victor Koppejan
Join Date: May 2015
Posts: 40
Rep Power: 11
vkoppejan is on a distinguished road
I understand that, however there is also the part about implicit/explicit treatment and how that affects the solver speed/stability and possbily the final species concentration (which is bounded).

Patankars book (1980) covers this in chapter 7.

I was wondering if you know of any applications or papers which use this approach for adsorption of a chemical species.
vkoppejan is offline   Reply With Quote

Old   November 29, 2015, 08:29
Default
  #6
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,893
Rep Power: 73
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
Quote:
Originally Posted by vkoppejan View Post
I understand that, however there is also the part about implicit/explicit treatment and how that affects the solver speed/stability and possbily the final species concentration (which is bounded).

Patankars book (1980) covers this in chapter 7.

I was wondering if you know of any applications or papers which use this approach for adsorption of a chemical species.
the time integration has nothing to do with the FV-based spatial discretization... if the source term is stiff, you should consider special techniques.
An old report is here http://ntrs.nasa.gov/archive/nasa/ca...9880008959.pdf

but more recent papers can be found in the international journals
FMDenaro is offline   Reply With Quote

Old   November 29, 2015, 08:30
Default
  #7
Member
 
Victor Koppejan
Join Date: May 2015
Posts: 40
Rep Power: 11
vkoppejan is on a distinguished road
Thanks, I'll look into the report.
vkoppejan is offline   Reply With Quote

Old   November 29, 2015, 10:26
Default
  #8
Senior Member
 
Michael Prinkey
Join Date: Mar 2009
Location: Pittsburgh PA
Posts: 363
Rep Power: 25
mprinkey will become famous soon enough
The linearization of source terms, especially for reacting flows, is a significant problem, at least for segregated solvers. I'm sorry that I can't give you a reference. This is knowledge I gleaned from the School of Hard Knocks while working on fuel cell modeling.

The first problem is a practical one--many (most) linear solvers and preconditioners require the A matrix to be diagonally dominate. So, whatever linearization you apply, it will be different for source and sink terms. I apologize for not being more clear about which is which--I'd need to do a little math to remind myself. But for reacting flow, you will have have both sources and sinks. That is a significant problem for conservation. In particular, species conservation becomes a convergence criteria instead of a guarantee based on the FV scheme. This has huge solution stability implication too, especially for combustion.

The second problem is related. Even if you have linear solvers that can handle A matrices that are not diagonally dominate, a segregated solver will still give you incompatible linearized source terms for different species--you can only linearize a source term with respect to ONE of the involved species at a time because you only solve one linearized species field at a time. Consider, 2 H2 + O2 -> 2 H2O. The source/sink terms for each of those reactions will enter the H2, O2, and H2O transport equations and each source/sink will vary linearly with the mass/mole fraction/density of H2, O2, and H2O respectively. But they will never balance perfectly, so species conservation again becomes a convergence criteria.

Solutions to this. Fractional step treatment of reactions is good. Each FV has a stiff ODE system that can use small/adaptive timesteps and can guarantee conservation by simply overwriting the composition/energy data in each FV. I guess it is possible to use a block solver for all species and energy together and linearize the reaction terms with respect to all of the species variables. Then the terms will be linearized consistently. Of course, that block solver needs to handle matrices that are not diagonally dominate as well.

The best place to start (if you haven't) is to just do the species reaction terms explicitly and dial down the timestep. Maybe there are other smarter ways to accomplish per-iteration species/energy conservation in a segregated implicit framework, but I don't know them. I'd be happy to learn about them though.
mprinkey is offline   Reply With Quote

Reply

Thread Tools Search this Thread
Search this Thread:

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


Similar Threads
Thread Thread Starter Forum Replies Last Post
Source Term due to evaporation in energy transport equation styleworker OpenFOAM Programming & Development 3 September 7, 2022 04:09
[Other] How to use finite area method in official OpenFOAM 2.2.0? Detian Liu OpenFOAM Meshing & Mesh Conversion 4 November 3, 2015 04:04
[swak4Foam] Swak4FOAM 0.2.3 / OF2.2.x installation error FerdiFuchs OpenFOAM Community Contributions 27 April 16, 2014 16:14
[swak4Foam] build problem swak4Foam OF 2.2.0 mcathela OpenFOAM Community Contributions 14 April 23, 2013 14:59
DxFoam reader update hjasak OpenFOAM Post-Processing 69 April 24, 2008 02:24


All times are GMT -4. The time now is 10:15.