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

Spectral Method/B Spline basis function

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

Reply
 
LinkBack Thread Tools Display Modes
Old   October 29, 2004, 11:20
Default Spectral Method/B Spline basis function
  #1
DJ
Guest
 
Posts: n/a
HI ALL, I AM USING SPECTRAL METHOD WITH B SPLINE AS THE BASIS FUNCTION FOR SOLVING THE N-S EQUATION(STEADY STATE LID DRIVEN CAVITY FLOW PROBLEM IN 2D). I AM USING NEWTON-RAPHSON METHOD FOR SOLVING SYSTEM OF NON LINEAR EQUATIONS. THE PROBLEM IS THAT JACOBAIN MATRIX OBTAINED AT EACH STEP OF THE N-R ITERATION IS VERY ILL CONDITIONED. FOR THE VERY FIRST ITERATION, IF I REPLACE ONLY ONE OF THE CONTINUITY EQUATION BY RANDOM EQUATIONS THEN CONDITION NUMBER IMPROVES CONSIDERABLY..IT MEANS THAT ONE OF THE CONTINUITY EQUATION IS LINEAR COMBINATION OF SOME OTHER EQUATION(??). HOWEVER..THE STORY IS ENTIRELY DIFFEENT WITH THE BIGGER SIZE MATRIX(REFINED MESH). AS I MAKE THE MESH MORE AND MORE REFINED, THE CONDITION NUMBER OF THE JACOBIAN BECOMES WORSE AND WORSE. I AM USING STAGGERED KNOT SEQUENCE FOR THE CALCULATION OF U AND P B -SPLINES. AND EVEYTHINGS SEEMS TO BE FINE IN MY CODE SO FAR(TYPO ERRORS ETC ETC..) ANY SUGGESTION..OR COMMENTS ABOUT THE ORIGIN OF THE (ALMOST) SINGULARITY OF THE JACOBIAN MATRIX.. AND REMEDIES TO REMOVE THE SAME. THANX IN ADVANCE,
  Reply With Quote

Reply

Thread Tools
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 On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Compile problem ivanyao OpenFOAM Running, Solving & CFD 1 October 12, 2012 09:31
ParaView for OF-1.6-ext Chrisi1984 OpenFOAM Installation 0 December 31, 2010 07:42
latest OpenFOAM-1.6.x from git failed to compile phsieh2005 OpenFOAM Bugs 25 February 9, 2010 05:37
Error with Wmake skabilan OpenFOAM Installation 3 July 28, 2009 00:35
Axisymmetrical mesh Rasmus Gjesing (Gjesing) OpenFOAM Native Meshers: blockMesh 10 April 2, 2007 14:00


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