Implicitisation of viscous terms

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

 April 3, 2006, 07:52 Implicitisation of viscous terms #1 Michael Halls-Moore Guest   Posts: n/a Sponsored Links Hi, I'm currently working on programming some computational fluid dynamics software for hypersonic flows (compressible, Mach number ~5). The solver application uses a Finite Volume 2nd Order (time and space) modified Godunov method (the Ben-Artzi/Falcovitz GRP method). This part solves the inviscid part of the system (i.e. the Euler equations) but I am now wanting to introduce viscosity into the simulation. I will be using operator splitting in order to seperate out the inviscid fluxes from the viscous fluxes. Since many of the applications I am studying are steady (i.e. no time-dependent behaviour), I want to use an implicit solver to calculate the viscosity fluxes. One way of doing this is to use the Thin Shear Layer (TSL) approximation, which essentially integrates the boundary-layer equations instead of the full navier-stokes equations. The solver already uses a TSL implicit formulation which produces a tridiagonal matrix, but I want to take it a bit further and calculate Full Navier Stokes. The solver uses a 2D structured orthogonal grid (no stretching, i.e. dx = dy for all cells). I'm currently ploughing through the algebra for adding in viscosity for a horizontal interface on the grid, for the momentum equation. Later on Ill try a vertical interface too, as well as the energy equation terms for both. At the moment, I'm trying to keep it simple! I have a few options as to what to do. My initial guess was to implicitise the u velocity terms (the horizontal velocity component of the field) and explicitise the v (vertical velocity component) velocity terms that arise in the equation. This would lead to a nice tridiagonal matrix solvable by Thomas'/TDMA algorithm. However, I will then have a new set of u values, which presumably I could use when updating the v values on my vertical interface. Question 1 is this: If I update the u values (i.e. do horizontal interfaces first) and use the new u values for updating v, will this give me the same answer as updating v (i.e. doing vertical interfaces first) and then using these new v values for updating u? I can't see any immediate reason why they would be the same, unfortunately. Any thoughts? Question 2: If I implicitise both u AND v for each interface, will this provide any benefits? I know it will create a very ugly solution matrix! Does anyone have any ideas or suggestions as to literature I can peruse to find out more about these issues? I know they're very specific, so it's quite tough to find the details! Thanks a lot, Mike.

 Thread Tools Display Modes Linear Mode

 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 OffTrackbacks are On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post dohnie OpenFOAM Programming & Development 2 March 21, 2012 04:05 nikosb Main CFD Forum 0 January 17, 2010 17:07 titio OpenFOAM Running, Solving & CFD 0 March 19, 2009 17:02 andy Main CFD Forum 7 June 16, 2006 11:38 Brindaban Ghosh Main CFD Forum 2 June 24, 2000 04:22