# Linearized NS euqations: how to solve them?(problem with Matrix operations..)

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

 November 18, 2009, 07:58 Linearized NS euqations: how to solve them?(problem with Matrix operations..) #1 Member   matteo lombardi Join Date: Apr 2009 Posts: 67 Rep Power: 10 Sponsored Links Hello, I am trying to solve the steady linearized navier-Stokes equation (which may arise from Newton Formulation of the NS or adjoint problems or other..). Let's be more general and let's say that i I have the following system: 1)Au+Bp=F 2)Cu=0 where u is the velocity vector, p is the pressure, A,B and C are the corrispective Matrix and F is the know-term of the momentum equations (sources+explicit terms+...) Possible classical startegies: a)Can I solve the problem with the classical Pressure Matrix formulation? that would mean getting u from equation 1) thus u=A^-1F - A^-1Bp thus i get the eqaution for p is: CA^-1Bp=CA^-1F Once p known, I can get U. (Of course I need good preconditioner and so on but let's remain simple at the moment..) To implement this method, I would need to be able to create those matrices, and be able to do matrices operations and solve linear problems passing to OF directly the matrices... Is it possible in openFoam? From what i know OF can solve only segregated problems (u and p segregated, and even u is solved segregated between components, i.e. solve first Ux, then Uy, then Uz). Am I right? Thus probably solvng A^1, where A is a vector matrix, is not possible, right? Am I correct about my assumptions about OF capabilities? b)if I had a stokes problem, i.e. no convective terms, I could apply the divergence operator to the momentum equations, get a laplacian equation for the pressure, and then correct u. Unfortunately I have some convective terms, thus some terms on the momentum equations don't go away.. thus i cannot "segregate" u and p with this approach.. Am i correct? any other ideas? c)I could try to adapt a strategy like the simple method to my problems, using some sorts of predictor/corrector strategy. Thus i solve the momentum euqations with the old U, then get the new P, correct U, get new P, correct U and so on until a certain convergence.. Does this make any sense?(I have tried to implement thsi strategy but so far quite unsuccessfully..) Could anyone point me to any article,book, notes about those problem is usually solved with Finite Volumes and segregated approach? (Mainly I want to solve the Oseen equation with finite Volume in a segregated way.) thank you very much, regards, Matteo

 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 lakeat OpenFOAM Running, Solving & CFD 42 August 26, 2009 21:47