Difference between revisions of "Presets/energy.xml"

Revision as of 14:20, 2 April 2015

Introduction

energy.xml is a preset used to compute the fluid energy components. If you wish to print the results in the screen and in an output file include energy_report.xml as well.

This tool is computing the following energy components variation rates for each particle:

${\displaystyle{\frac{dE_{a}^{{kin}}}{dt}}=m_{a}{\mathbf{u}}_{a}\cdot{\frac{d{\mathbf{u}}_{a}}{dt}},}$

${\displaystyle{\frac{dE_{a}^{{pot}}}{dt}}=-m_{a}{\mathbf{g}}\cdot{\mathbf{u}}_{a},}$

${\displaystyle{\frac{dW_{a}}{dt}}=-{\frac{m_{a}}{\rho_{a}}}\nabla p_{a}\cdot{\mathbf{u}}_{a}-{\frac{m_{a}}{\rho_{a}}}p_{a}\nabla\cdot{\mathbf{u}}_{a},}$

${\displaystyle T{\frac{dS_{a}}{dt}}=-\mu{\frac{m_{a}}{\rho_{a}}}{\mathbf{u}}_{a}\cdot\Delta{\mathbf{u}}_{a}-\delta{\frac{p_{a}m_{a}\Delta t}{\rho_{0}\rho_{a}}}\Delta p_{a},}$

where ${\displaystyle E_{a}^{{kin}}}$ , ${\displaystyle E_{a}^{{pot}}}$ are the kinetic and the potential energy components for a generic particle a, ${\displaystyle W_{a}}$ is the work of the external forces, and ${\displaystyle S_{a}}$ is the entropy.

These values per particle are later reduced to the total fluid variation:

${\displaystyle{\frac{dE^{{kin}}}{dt}}=\sum_{{a\in{\mathrm{Fluid}}}}{\frac{dE_{a}^{{kin}}}{dt}},}$

${\displaystyle{\frac{dE^{{pot}}}{dt}}=\sum_{{a\in{\mathrm{Fluid}}}}{\frac{dE_{a}^{{pot}}}{dt}},}$

${\displaystyle{\frac{dW}{dt}}=\sum_{{a\in{\mathrm{Fluid}}}}{\frac{dW_{a}}{dt}},}$

${\displaystyle T{\frac{dS}{dt}}=\sum_{{a\in{\mathrm{Fluid}}}}T{\frac{dS_{a}}{dt}}.}$

Usage

To use it the following tag should be added to include the preset:

<Include file="/usr/share/aquagpusph/resources/Presets/energy.xml" />

Where /usr/share/aquagpusph/ should be conveniently modified.

This preset should be loaded after cfd.xml.

Variables defined

When this preset is loaded some variables are generated (on top of the default variables of AQUAgpusph):

Tools

The following tools are set to be executed each time step, just before the "t = t + dt" tool. (see Tools and cfd.xml to learn more about this):