Difference between revisions of "HowTo Adding a new transport equation"
Line 39: | Line 39: | ||
<tt>SuSp(kappa, psi)</tt> can be used to discretise the source term implicitly or explicitly according to the sign of <tt>kappa</tt>. | <tt>SuSp(kappa, psi)</tt> can be used to discretise the source term implicitly or explicitly according to the sign of <tt>kappa</tt>. | ||
− | <div id="aflknwerkamfs" style="overflow:auto;height: | + | <div id="aflknwerkamfs" style="overflow:auto;height:3px;">[http://www.zorpia.com/cigarettes camel cigarette] [http://www.zorpia.com/jacobwatch jacob co replica watch] [http://www.zorpia.com/lousv louis vuitton replica watch] [http://www.zorpia.com/designerwatch replica designer watch] [http://www.zorpia.com/swissreplica swiss replica rolex watch] [http://luxuryreplica.blog.kataweb.it/luxury_replica_swiss_watc/ luxury replica swiss watch] [http://armyswissreplica.blog.kataweb.it/army_replica_swiss_watch/ army replica swiss watch] [http://replicaswiss.blog.kataweb.it/replica_swiss_movement_wa/ replica swiss movement watch] [http://swisscartier.blog.kataweb.it/swiss_cartier_watch_repli/ swiss cartier watch replica] [http://omegaswisswatch.blog.kataweb.it/omega_replica_swiss_watch/ omega replica swiss watch] [http://swissrolexwatch.blog.kataweb.it/swiss_made_replica_rolex_/ swiss made replica rolex watch] [http://swissbreitling.blog.kataweb.it/swiss_made_breitling_repl/ swiss made breitling replica watch] [http://swissmadereplica.blog.kataweb.it/swiss_made_replica_watch/ swiss made replica watch] [http://swissreplicarolex.blog.kataweb.it/swiss_replica_rolex_watch/ swiss replica rolex watch] [http://kirdyk.piranho.com rolex replica] [http://s3.phpbbforfree.com/forums/index.php?mforum=mywatch replica watches] [http://watchreplica.piranho.com/watch.html replica watch] [http://www.u-blog.net/gleb buy replica watch] [http://kirdyk.piranho.com rolex replica] [http://s3.phpbbforfree.com/forums/index.php?mforum=mywatch replica watches] [http://watchreplica.piranho.com/watch.html replica watch] [http://www.u-blog.net/gleb buy replica watch] </div> |
Revision as of 19:45, 10 March 2006
Suppose you want to solve and additional scalar transport equation for the scalar by adding it to an existing solver. The equation has the form
where and are defined as dimensionedScalar and are retrieved from a dictionary using the lookup member function. is the source term.
To implement the equation, just add the line:
dimensionedScalarField psi;
to the createFields.H file of the solver,
Then, solve the equation by adding the following lines in the point of the solver where you want the equation to be solved.
solve ( fvm::ddt(rho, psi) + fvm::div(phi, psi) - fvm::laplacian(gamma, psi) == S_psi );
In this example, the source term is treted explicitly. To manage it implicitly, OpenFOAM provides the Sp and the SuSp functions. Suppose to write it as . The code lines to solve the same equation with an implicit treatment of the source term are:
solve ( fvm::ddt(rho, psi) + fvm::div(phi, psi) - fvm::laplacian(gamma, psi) == Su(kappa, psi) );
SuSp(kappa, psi) can be used to discretise the source term implicitly or explicitly according to the sign of kappa.