mirror of
https://github.com/Cantera/cantera.git
synced 2026-08-17 16:35:17 -05:00
Add author contact information to sample
This commit is contained in:
committed by
Bryan W. Weber
parent
e78975e92d
commit
57b4641d4b
@@ -1,10 +1,10 @@
|
||||
function [F] = PFR_solver(x,soln_vector,gas,mdot,A_in,dAdx,k)
|
||||
|
||||
% This function defines the spatial derivatives for an ideal gas plug-flow
|
||||
% reactor, where the cross-sectional area and pressure are allowed to vary,
|
||||
% axially. The model is set up by the example file 'Plug_Flow_Reactor.m',
|
||||
% which points the integrator to this function. The integrator integrates the
|
||||
% derivatives spatially, to solve the density, temperature, and species mass
|
||||
% This function defines the spatial derivatives for an ideal gas plug-flow
|
||||
% reactor, where the cross-sectional area and pressure are allowed to vary,
|
||||
% axially. The model is set up by the example file 'Plug_Flow_Reactor.m',
|
||||
% which points the integrator to this function. The integrator integrates the
|
||||
% derivatives spatially, to solve the density, temperature, and species mass
|
||||
% fraction profiles as a function of distance x.
|
||||
|
||||
rho = soln_vector(1);
|
||||
|
||||
@@ -12,8 +12,13 @@
|
||||
% *S.R Turns, An Introduction to Combustion - Concepts and Applications,
|
||||
% McGraw Hill Education, India, 2012, 206-210.*
|
||||
%
|
||||
% The current example is written for methane combustion, but can be readily
|
||||
% The current example is written for methane combustion, but can be readily
|
||||
% adapted for other chemistries.
|
||||
%
|
||||
% Developed by Ashwin Kumar/Dr.Joseph Meadows (mgak@vt.edu/jwm84@vt.edu) on 3-June-2020
|
||||
% Research Assistant/Assistant Professor
|
||||
% Advanced Propulsion and Power Laboratory
|
||||
% Virginia Tech
|
||||
|
||||
%% Clear all variables, close all figures, clear the command line:
|
||||
clear all
|
||||
@@ -80,7 +85,7 @@ nsp = nSpecies(gas_calc);
|
||||
|
||||
% Initialize arrays for T, Y, and rho at each location:
|
||||
T_calc = zeros(length(x_calc),1);
|
||||
Y_calc = zeros(length(x_calc), nsp);
|
||||
Y_calc = zeros(length(x_calc),nsp);
|
||||
rho_calc = zeros(length(x_calc),1);
|
||||
|
||||
T_calc(1) = temperature(gas_calc);
|
||||
|
||||
Reference in New Issue
Block a user