水环境学——显式差分计算均匀河段各断面BOD变化

最新推荐文章于 2022-02-24 21:10:08 发布

iSuwen

最新推荐文章于 2022-02-24 21:10:08 发布

阅读量1.9k

点赞数 2

分类专栏：爱进步 C++程序文章标签： include c

本文链接：https://blog.csdn.net/Superwen_go/article/details/7639971

版权

爱进步 C++程序专栏收录该内容

27 篇文章 8 订阅

订阅专栏

#include "stdafx.h"
#include <iostream>
#include <fstream>
#include <iomanip>
#include <cmath>
void main()
{
	using namespace std;
	const int Nx = 8,  //等于x/dx
				  Nt = 10;//等于t/dt
	const double dx = 1.0, dt = 0.1, //要保证dt<=dx/u，E*dt<=dx*dx/2
						x = 8.0,//均匀河段长
						t = 1.0,//排污断面持续时间
						u = 5,//流速
						E = 2,//离散系数
						K1 = 0.015,//BOD降解系数
						C0 = 20;//排污断面持续时间内保持的BOD浓度
	double alpha  = E*dt/pow(dx, 2),
				 beta = u*dt/dx - 2*alpha - K1*dt,
				 gamma = 1 - u*dt/dx + alpha, 
				 C[Nx + 2][Nt + 1] = {0};
	ofstream outfile;
	outfile.open("output.txt");
	outfile<<"alpha  = "<<alpha<<endl
			<<"beta = "<<beta<<endl
			<<"gamma = "<<gamma<<endl;
	for(int i = 0; i < 2; i++)
		for(int j = 0; j<Nt + 1; j++)
			C[i][j] = C0;
	for(int i = 2; i < Nx + 2; i++)
		for(int j = 1; j<Nt + 1; j++)
			C[i][j] = alpha*C[i - 2][j - 1] +beta*C[i - 1][j - 1] + gamma*C[i][j - 1];
	for(int j = 0; j<Nt + 1; j++)
	{
		for(int i = 1; i < Nx + 2; i++)
			outfile<<setw(10)<<C[i][j];
		outfile<<endl;
	}
	outfile.close();
}