# buoyantPimpleFoam General Assumptions (thermophysical model, energy equation, etc)

 User Name Remember Me Password
 Register Blogs Members List Search Today's Posts Mark Forums Read

 October 7, 2017, 03:37 buoyantPimpleFoam General Assumptions (thermophysical model, energy equation, etc) #1 New Member   Join Date: Oct 2017 Posts: 12 Rep Power: 7 I am trying to use buoyantPimpleFoam to simulate Rayleigh-Taylor instability of a gas (eventually a reacting simulation). I was looking through the source files and I had some general questions on the assumptions that the solver takes, and hope that someone can shed some light on it for me. 1. The energy equation is states as Code: fvScalarMatrix EEqn ( fvm::ddt(rho, he) + fvm::div(phi, he) + fvc::ddt(rho, K) + fvc::div(phi, K) + ( he.name() == "e" ? fvc::div ( fvc::absolute(phi/fvc::interpolate(rho), U), p, "div(phiv,p)" ) : -dpdt ) - fvm::laplacian(turbulence->alphaEff(), he) == rho*(U&g) + radiation->Sh(thermo, he) + fvOptions(rho, he) ); All the terms look correct (though it is missing mechanical stress term), but the thermal diffusion term seems a bit odd to me. If we use Fourier's law of heat diffusion, we would have: But in the solver energy equation, this is simply states as . Obviously this does not matter for a thermally perfect gas as the specific heat is constant, but for gases with varying specific heats, there should be a gradient there. I have never seen the form used in buoyantPimpleFoam. Is this additional gradient term simply just ignored because it is negligible? 2. The divergence of gradients in OpenFoam (e.g. or ) are simply written as laplacians, i.e. laplacian(mu,D(U)). Again, this is only true if the transport properties do not vary. This is definitely not the case for reacting flows. Is this simply ignored in OpenFoam calculations or is this handled behind the scenes? 3. The thermophysical properties for gases seem to me to be overdefined. For example, a sample thermophysical model could be: Code: thermoType { type heRhoThermo; //Thermophysical Model Based on h/e and rho mixture pureMixture; // Passive Gas Mixtures transport const; // Constant transport properties (mu, Pr, etc.) thermo hConst; // Constant Cp and heat of fusion equationOfState perfectGas; // Ideal Gas Law specie specie; energy sensibleEnthalpy; // Solve Energy Eqn. In Terms of Enthalpy } From my understanding, heRhoThermo calculates the enthalpy and density, but the density and the enthalpy are already defined by the perfect gas and 'hConst' assumptions. What additional things is the thermophysical type doing to calculate these? Is the type even necessary at this point? 4. Why is read from the initial conditions and then re-defined, i.e. in createFields.H Code: #include "readGravitationalAcceleration.H" #include "readhRef.H" #include "gh.H" Info<< "Reading field p_rgh\n" << endl; volScalarField p_rgh ( IOobject ( "p_rgh", runTime.timeName(), mesh, IOobject::MUST_READ, IOobject::AUTO_WRITE ), mesh ); // Force p_rgh to be consistent with p p_rgh = p - rho*gh; If you've correctly defined Prgh, then it should be consistent already. I'm not quite sure what they are doing here. Thank you in advance.

 October 23, 2018, 05:30 #2 Senior Member   Jianrui Zeng Join Date: May 2018 Location: China Posts: 157 Rep Power: 7 Hello,I am also confused about these questions above. Have you found the answer of it? Any idea will be appreciated.

 Tags assumptions, buoyantpimplefoam, consistency, laplacian operator

 Thread Tools Search this Thread Search this Thread: Advanced Search 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 Off Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post Joseph CFX 14 April 20, 2010 15:45 Rich Main CFD Forum 0 December 16, 2009 14:01 Beginner Main CFD Forum 1 January 7, 2009 05:36 Fabio Main CFD Forum 0 June 1, 2007 06:06 Miguel Baritto CFX 4 August 31, 2006 12:02

All times are GMT -4. The time now is 01:43.

 Contact Us - CFD Online - Privacy Statement - Top