Rheolef  7.2
an efficient C++ finite element environment
 
Loading...
Searching...
No Matches
combustion_exact.icc

The combustion problem – its exact solution.

The combustion problem – its exact solution

#include "lambda2alpha.h"
struct u_exact {
Float operator() (const point& x) const {
return 2*log(cosh(a)/cosh(a*(1-2*x[0]))); }
u_exact (Float lambda, bool is_upper)
: a(lambda2alpha(lambda,is_upper)) {}
u_exact (Float a1) : a(a1) {}
};
struct grad_u {
point operator() (const point& x) const {
return point(4*a*tanh(a*(1-2*x[0]))); }
grad_u (Float lambda, bool is_upper)
: a(lambda2alpha(lambda,is_upper)) {}
grad_u (Float a1) : a(a1) {}
};
see the Float page for the full documentation
see the point page for the full documentation
The combustion problem – inversion of the parameter function.
Float lambda2alpha(Float lambda, bool up=false)
point operator()(const point &x) const
point operator()(const point &x) const