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

Boundary Value Problem for CFD Case

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   July 11, 2017, 07:23
Default Boundary Value Problem for CFD Case
  #1
New Member
 
Chus Mat
Join Date: Feb 2017
Posts: 12
Rep Power: 7
chusma is on a distinguished road
I am currently very frustrated with this problem. Is a Boundary Value Problem for ODEs system as an approach of quasi-steady one dimensional flow in a rocket combustion chamber. The equations have the shape on the picture eq.

The ideal gass state equation closes the system.

The BC's are: At the beginning of the domain as Initial Values: T(0) = value; u(0) = 0, aprox. And at the end of the domain as clossure, the mass flow provides the condition shown in the picture b.c

Must be assumed the values which are not the flow magnitudes (ρ,p,T,u)are known.

I am currently applying a finite differences approach with an Explicit Euler's for the ODEs. I have a lot of problems to choose the best method for this problem. I have tryed with different ideas, the most recent was the Shooting Method with the secant method for the nonliniarrity.

I am very lost with this problem and it is the main part of my Bachelors Degree Thesis, so I would appreciate any help.

Thank you for your attention
Attached Images
File Type: png Eq.PNG (10.9 KB, 13 views)
File Type: png bc.PNG (4.0 KB, 5 views)
chusma is offline   Reply With Quote

Old   July 11, 2017, 09:34
Default
  #2
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,146
Rep Power: 66
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
What's wrong with your solution?
FMDenaro is offline   Reply With Quote

Old   July 11, 2017, 09:38
Default
  #3
New Member
 
Chus Mat
Join Date: Feb 2017
Posts: 12
Rep Power: 7
chusma is on a distinguished road
Quote:
Originally Posted by FMDenaro View Post
What's wrong with your solution?
Thanks for your answer

The problem is that my scheme diverges during the iteration of the shooting method.

I have not so much experience in this kind of programming and I do not know how to avoid or prevent that situation.

If there is another scheme or method that could be more suitable, please tell me.
chusma is offline   Reply With Quote

Old   July 11, 2017, 09:42
Default
  #4
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,146
Rep Power: 66
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
Divergence can be due to a numerical stability issue... check for the step value, does it satisfy the constraint?
FMDenaro is offline   Reply With Quote

Old   July 11, 2017, 13:11
Default
  #5
New Member
 
Chus Mat
Join Date: Feb 2017
Posts: 12
Rep Power: 7
chusma is on a distinguished road
Quote:
Originally Posted by FMDenaro View Post
Divergence can be due to a numerical stability issue... check for the step value, does it satisfy the constraint?
Which constraint are you talking about?

Excuse my ignorance.
chusma is offline   Reply With Quote

Old   July 11, 2017, 13:32
Default
  #6
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,146
Rep Power: 66
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
Quote:
Originally Posted by chusma View Post
Which constraint are you talking about?

Excuse my ignorance.
If you are not familiar with numerical methods for solving ODE, I strongly suggest reading this subject on a good textbook of numerical analysis.
Numerical stability is one of the major issue in explicit methods and you definitely need to take care of that.
FMDenaro is offline   Reply With Quote

Old   July 11, 2017, 13:34
Default
  #7
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,146
Rep Power: 66
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
Note that you are working with a system of first order ODE, so be carefull in you boundary conditions to ensure that the mathematical problem is well posed.
FMDenaro is offline   Reply With Quote

Old   July 11, 2017, 13:38
Default
  #8
New Member
 
Chus Mat
Join Date: Feb 2017
Posts: 12
Rep Power: 7
chusma is on a distinguished road
Quote:
Originally Posted by FMDenaro View Post
Note that you are working with a system of first order ODE, so be carefull in you boundary conditions to ensure that the mathematical problem is well posed.
The problem comes from the shooting method iteration, not from the explicit scheme, however.

The scheme is made to converge but the iteration is not. I would like to fix that
chusma is offline   Reply With Quote

Reply

Tags
boundary value problems, boundary variable, navier-stokes, non linearity

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
Wind turbine simulation Saturn CFX 58 July 3, 2020 01:13
boundary conditions- new case problem Sebaj OpenFOAM 20 November 23, 2017 12:58
Wrong flow in ratating domain problem Sanyo CFX 17 August 15, 2015 06:20
Radiation interface hinca CFX 15 January 26, 2014 17:11
Low Mixing time Problem Mavier CFX 5 April 29, 2013 00:00


All times are GMT -4. The time now is 18:55.