CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Running, Solving & CFD

About sigmaExp of solidEquilibriumDisplacementFoam OpenFOAM 14

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

Reply
 
LinkBack Thread Tools Display Modes
Old   May 24, 2007, 23:04
Default In solidEquilibriumDisplacemen
  #1
New Member
 
Jerry Wei
Join Date: Mar 2009
Posts: 1
Rep Power: 0
jerrywei is on a distinguished road
In solidEquilibriumDisplacementFoam,there are two Sigma tensor:

sigmaExp = (lambda - mu)*gradDcorr + mu*gradDcorr.T()
+ (lambda*I)*tr(gradDcorr);

sigmaD += accFac*(mu*twoSymm(gradDcorr)
+ (lambda*I)*tr(gradDcorr));

and solve equation with sigmaD and sigmaExp:
solve
(
fvm::laplacian(2*mu + lambda, Dcorr, "laplacian(DD,Dcorr)")
+ fvc::div(sigmaExp + sigmaD)
);



But in StressedFoam(OpenFOAM 1.3):

fvVectorMatrix UEqn
(
fvm::d2dt2(U)
==
fvm::laplacian(2*mu + lambda, U, "laplacian(DU,U)")

+ fvc::div
(
mu*gradU.T() + lambda*(I*tr(gradU)) - (mu + lambda)*gradU,
"div(sigma)"
)
);

I can't understand about that the coefficient is not the same between sigmaExp and stressedFOAM:

sigmaExp: (lambda - mu)*gradDcorr

stressedFOAM: - (mu + lambda)*gradU

and is there reference paper about the new solver of solidEquilibriumDisplacementFoam?(The method is cool and good)

thank you very much!!!
Jerry
jerrywei is offline   Reply With Quote

Old   September 15, 2011, 14:16
Default
  #2
Member
 
Santiago
Join Date: Dec 2009
Posts: 85
Rep Power: 7
gascortado is on a distinguished road
I have question related to StressedFoam. The equation to be solved appears to be the one described here (equation #6) : http://powerlab.fsb.hr/ped/kturbo/Op...tressPaper.pdf

However, the solver solves this equation :

fvVectorMatrix UsolidEqn
(
fvm::d2dt2(Usolid)
==
fvm::laplacian(2*mu + lambda, Usolid, "laplacian(DU,U)")
+ fvc::div
(
mu*gradUsolid.T() + lambda*(I*tr(gradUsolid)) - (mu + lambda)*gradUsolid,
"div(sigma)"
)
);


I dont understand where the (mu + lambda)*gradUsolid and fvm::laplacian(2*mu + lambda, Usolid, "laplacian(DU,U)") terms is coming from


Any idea?? I'll appreciate any help. THanks
gascortado is offline   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
OpenFoam vs CFX5 mass balance in OpenFoam tangd OpenFOAM Running, Solving & CFD 33 May 23, 2010 17:36
OpenFOAM with IBM AIX matthias OpenFOAM Installation 20 March 25, 2008 04:36
CheckMesh error using a tutorial from OpenFOAM 114 with openFOAM 13 martapajon OpenFOAM Native Meshers: blockMesh 7 January 21, 2008 13:52
OpenFOAM users in Munich OpenFOAM benutzer in M%c3%bcnchen jaswi OpenFOAM 0 August 3, 2007 14:11
A new Howto on the OpenFOAM Wiki Compiling OpenFOAM under Unix mbeaudoin OpenFOAM Installation 2 April 28, 2006 09:54


All times are GMT -4. The time now is 04:07.