|
[Sponsors] |
Temperature-dependent properties and mass conservation in porousExplicitSourceReactin |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
October 26, 2010, 06:46 |
Temperature-dependent properties and mass conservation in porousExplicitSourceReactin
|
#1 |
Member
Robertas N.
Join Date: Mar 2009
Location: Kaunas, Lithuania
Posts: 53
Rep Power: 17 |
Hi *,
I modified the 'porousExplicitSourceReactingParcelFoam' to include heat sources in the energy equation and set up a very simple simulation: a box filled with pure N2, no mass flows (U=uniform(0,0,0)) and a uniform field of the heat source, in order to raise the temperature; that's just for testing purposes -- the setup will get more complex in the future. The solver uses the polynomial coefficients defined in constant/thermo.incompressiblePoly for temperature-dependent parameters: density, heat capacity @p=const and some other. As the temperature rises, the density gets lower, therefore, the total mass in the box (calculated as m=sum_over_cells(rho(cell)*V(cell))) diminishes, but the pressure stays constant. Of course, this setup is incorrect -- it must use heat capacity @ constant volume, not @ constant pressure, and assume the compressible fluid (gas in this case). I have a few questions: 1) Is it possible to obtain somewhere the polynomial coefficients for heat capacity at constant volume? Where are the data in thermo.incompressiblePoly taken from? I couldn't find any references -- ChemKin, maybe? 2) Why does the pressure stay constant in the current simulation? Isn't 'porousExplicitSourceReactingParcelFoam' meant to be compressible solver? The modified energy equation looks like this: solve ( fvm::ddt(rho, h) + mvConvection->fvmDiv(phi, h) - fvm::laplacian(turbulence->alphaEff(), h) == pWork() + parcels.Sh() + radiation->Sh(thermo) + Q_heat ); All the terms on the RHS are zero, except for Q_heat, which is defined as a volScalarField with constant power density (in W/m^3). Did I overlook something? Thank you for any pointers, hints, keywords for searching... |
|
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
mass flow in is not equal to mass flow out | saii | CFX | 12 | March 19, 2018 05:21 |
Simulation of a single bubble with a VOF-method | Suzzn | CFX | 21 | January 29, 2018 00:58 |
Constant velocity of the material | Sas | CFX | 15 | July 13, 2010 08:56 |
Mass Properties of Axial Fan Bladegen | Katamish | CFX | 2 | May 12, 2009 05:39 |
Two-Phase Buoyant Flow Issue | Miguel Baritto | CFX | 4 | August 31, 2006 12:02 |