|
[Sponsors] |
January 7, 2012, 11:13 |
1D fluid test case
|
#1 |
Senior Member
Join Date: Jun 2010
Posts: 111
Rep Power: 15 |
Hi,
I am trying to think of a 1D fluid test case, where the flow is laminar and the fluid is near-compressible and the process is isothermal. I was thinking of just a moving fluid grid with a pressure applied at an inlet, not sure what BC to use at the other end ... Please let me know if you can think of any examples. I just need to test out a code. Any hints are really appreciated. Thanks, Hooman |
|
January 7, 2012, 13:15 |
|
#2 |
Senior Member
cfdnewbie
Join Date: Mar 2010
Posts: 557
Rep Power: 20 |
Hi Hooman,
If you just want to test your code, you could use the manufactured solution method and generate a solution for your case. It should be very simple for 1d NS, and you can test the order of accuracy of your code! Hope this helps! |
|
January 7, 2012, 13:32 |
|
#3 |
Senior Member
Join Date: Jun 2010
Posts: 111
Rep Power: 15 |
Thanks very much,
I'm not sure what you mean by 'the manufactured solution method', possibly you're referring to available CFD codes. I've written my own code and need to test it for a 1D test case before moving on to 2D. I've already tested a simple 1D case where velocity and density are constants just to test the discretization methods. I now need a test case where density, velocity and pressure all vary. The only 1D example of this sort that I've seen is a fluid piston connected to a solid spring. I'm wondering if something similar can be done with just the fluid and replacing velocity of the fluid where it's connected to the solid with a known function. Also these test cases usually use gas but I'm using a near compressible flow so I assume would need a pressure or velocity inlet. Would this be reasonable? I would really appreciate your help. |
|
January 9, 2012, 12:53 |
|
#4 |
Senior Member
Join Date: Oct 2011
Posts: 239
Rep Power: 16 |
Hi !
The method of manufactured solution consists in choosing an arbitrary solution (or vector solution) and modifying accordingly the governing equations by adding source terms so that the equations are satisfied. Here are two references: Roache, P.J., Verification and validation in computational science and engineering, Albuquerque, NM: Hermosa Publishers, 1998 Roy, C. J., Review of code and solution verification procedures for computational simulation, Journal of Computational Physics, 205, pp. 131-156, 2005 So basically, you define your density, velocity and pressure to some infinitely differentiable function (usually trigonometric functions). You inject these solutions inside the PDE and extract the additional terms that you have to add to your discrete equations. In this way you can test the order of accuracy of your code (except bc which need some special treatments if not set to dirichlet) |
|
January 10, 2012, 03:45 |
|
#5 |
Senior Member
duri
Join Date: May 2010
Posts: 245
Rep Power: 16 |
Test a simple converging or diverging duct, use area at node same as actual 3d area.
|
|
Thread Tools | Search this Thread |
Display Modes | |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Pressure instabilities with interDyMFoam for the floatingObject case | nbadano | OpenFOAM Running, Solving & CFD | 15 | October 15, 2021 07:35 |
How to find standard store-seperation test case? | peter.zhao | Main CFD Forum | 6 | March 6, 2015 12:31 |
rhoSimpleFoam - test case | atareen64 | OpenFOAM Running, Solving & CFD | 5 | May 1, 2011 18:17 |
Help Required With Simple Test Case | steph79 | OpenFOAM Pre-Processing | 4 | August 3, 2010 08:45 |
Test case data | Snegan | Main CFD Forum | 3 | June 14, 2006 09:56 |