research/cantera
1
0
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity Actions

[1D] Added a minimum grid spacing parameter

Browse source
This commit is contained in:
Ray Speth 2012-12-06 17:31:44 +00:00
parent 6ddb55e020
commit b4f34045a8
4 changed files with 40 additions and 3 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

8
include/cantera/oneD/Sim1D.h
View file

@ -107,6 +107,14 @@ public:
doublereal slope = 0.8, doublereal curve = 0.8, doublereal prune = -0.1);
void setMaxGridPoints(int dom = -1, int npoints = 300);
//! Set the minimum grid spacing in the specified domain(s).
/*!
* @param dom Domain index. If dom == -1, the specified spacing
is applied to all domains.
@param gridmin The minimum allowable grid spacing [m]
*/
void setGridMin(int dom, double gridmin);
void restore(const std::string& fname, const std::string& id);
void getInitialSoln();

13
include/cantera/oneD/refine.h
View file

@ -31,6 +31,18 @@ public:
void setMaxPoints(int npmax) {
m_npmax = npmax;
}
//! Set the minimum allowable spacing between adjacent grid points [m].
void setGridMin(double gridmin) {
m_gridmin = gridmin;
}
//! Returns the the minimum allowable spacing between adjacent
//! grid points [m].
double gridMin() const {
return m_gridmin;
}
int analyze(size_t n, const doublereal* z, const doublereal* x);
int getNewGrid(int n, const doublereal* z, int nn, doublereal* znew);
//int getNewSoln(int n, const doublereal* x, doublereal* xnew);
@ -69,6 +81,7 @@ protected:
Domain1D* m_domain;
size_t m_nv, m_npmax;
doublereal m_thresh;
doublereal m_gridmin; //!< minimum grid spacing [m]
};

14
src/oneD/Sim1D.cpp
View file

@ -624,6 +624,20 @@ void Sim1D::setRefineCriteria(int dom, doublereal ratio,
}
}
void Sim1D::setGridMin(int dom, double gridmin)
{
if (dom >= 0) {
Refiner& r = domain(dom).refiner();
r.setGridMin(gridmin);
} else {
for (size_t n = 0; n < m_nd; n++) {
Refiner& r = domain(n).refiner();
r.setGridMin(gridmin);
}
}
}
void Sim1D::setMaxGridPoints(int dom, int npoints)
{
if (dom >= 0) {

8
src/oneD/refine.cpp
View file

@ -18,7 +18,8 @@ static void r_drawline()
Refiner::Refiner(Domain1D& domain) :
m_ratio(10.0), m_slope(0.8), m_curve(0.8), m_prune(-0.001),
m_min_range(0.01), m_domain(&domain), m_npmax(3000)
m_min_range(0.01), m_domain(&domain), m_npmax(3000),
m_gridmin(5e-6)
{
m_nv = m_domain->nComponents();
m_active.resize(m_nv, true);
@ -107,7 +108,7 @@ int Refiner::analyze(size_t n, const doublereal* z,
dmax = m_slope*(vmax - vmin) + m_thresh;
for (j = 0; j < n-1; j++) {
r = fabs(v[j+1] - v[j])/dmax;
if (r > 1.0) {
if (r > 1.0 && dz[j] >= 2 * m_gridmin) {
m_loc[j] = 1;
m_c[name] = 1;
//if (int(m_loc.size()) + n > m_npmax) goto done;
@ -139,7 +140,8 @@ int Refiner::analyze(size_t n, const doublereal* z,
dmax = m_curve*(smax - smin); // + 0.5*m_curve*(smax + smin);
for (j = 0; j < n-2; j++) {
r = fabs(s[j+1] - s[j]) / (dmax + m_thresh/dz[j]);
if (r > 1.0) {
if (r > 1.0 && dz[j] >= 2 * m_gridmin &&
dz[j+1] >= 2 * m_gridmin) {
m_c[name] = 1;
m_loc[j] = 1;
m_loc[j+1] = 1;