Edgar J. Kaiser
Certified Consultant
Please login with a confirmed email address before reporting spam
Posted:
7 years ago
1 juin 2017, 09:23 UTC−4
Updated:
7 years ago
1 juin 2017, 09:27 UTC−4
Hi Madlen,
just set up an integration on the seabed domain. If rho is the specific heat capacity in J/m3 and T0 the start ambient temperature:
Wth is the stored heat and V is the seabed domain. You can set up the integration in the definitions section and use the integration operator to set up a global variable, or you can do it in post processing in case your model is already solved.
Cheers
Edgar
--
Edgar J. Kaiser
emPhys Physical Technology
www.emphys.com
Hi Madlen,
just set up an integration on the seabed domain. If rho is the specific heat capacity in J/m3 and T0 the start ambient temperature:
[math]W_{th} = \int_{V} \rho (T - T_0)dV}[/math]
Wth is the stored heat and V is the seabed domain. You can set up the integration in the definitions section and use the integration operator to set up a global variable, or you can do it in post processing in case your model is already solved.
Cheers
Edgar
--
Edgar J. Kaiser
emPhys Physical Technology
http://www.emphys.com
Please login with a confirmed email address before reporting spam
Posted:
7 years ago
4 juil. 2017, 06:58 UTC−4
Hi Edgar,
sorry for the late replay and thanks for the answer. I have now tried it but didnt succeed yet.
I added an integration node (named intop1) to my component definitions and selected the domain for it. Where in Comsol can I add the formula you gave earlier?
I also tried in the post processing but am not sure if I should take a surface or volume integration. My model is in 2D, but in total I would like to obtain a value for the heat content of the sediment in J along the whole cable which is some kilometers long in reality.
And one question about the formula: For me the heat capacity is Cp in J/kgK. Shouldn't I add also the densitiy (rho) which is in kg/m³ to the integration formula? Then after the integration the unit that remains would be only J.
Thank you for answering these beginner questions,
Best Madlen
Hi Edgar,
sorry for the late replay and thanks for the answer. I have now tried it but didnt succeed yet.
I added an integration node (named intop1) to my component definitions and selected the domain for it. Where in Comsol can I add the formula you gave earlier?
I also tried in the post processing but am not sure if I should take a surface or volume integration. My model is in 2D, but in total I would like to obtain a value for the heat content of the sediment in J along the whole cable which is some kilometers long in reality.
And one question about the formula: For me the heat capacity is Cp in J/kgK. Shouldn't I add also the densitiy (rho) which is in kg/m³ to the integration formula? Then after the integration the unit that remains would be only J.
Thank you for answering these beginner questions,
Best Madlen
Edgar J. Kaiser
Certified Consultant
Please login with a confirmed email address before reporting spam
Posted:
7 years ago
4 juil. 2017, 16:42 UTC−4
Updated:
7 years ago
4 juil. 2017, 16:43 UTC−4
Hi Madlen,
You can do it in postprocessing of course. The volume integral in 2D is just the surface integral multiplied by the model length, that is typically set in the physics top node.
And you are right, the
in my formula has to be the gravimetric density multiplied by the heat capacity. Both can be a function of the spatial variables.
You can also define a variable:
heat = intop1(Cp*rho*(T-T0))
To get the heat in the volume you can again multiply by the model length.
Cheers
Edgar
--
Edgar J. Kaiser
emPhys Physical Technology
www.emphys.com
Hi Madlen,
You can do it in postprocessing of course. The volume integral in 2D is just the surface integral multiplied by the model length, that is typically set in the physics top node.
And you are right, the [math]\rho[/math] in my formula has to be the gravimetric density multiplied by the heat capacity. Both can be a function of the spatial variables.
You can also define a variable:
heat = intop1(Cp*rho*(T-T0))
To get the heat in the volume you can again multiply by the model length.
Cheers
Edgar
--
Edgar J. Kaiser
emPhys Physical Technology
http://www.emphys.com