How to give initial guess and initial condition for a nonlinear - time dependent PDE problem (transient non-linear heat equation) in pde tool box?
3 views (last 30 days)
Show older comments
Sravan Kumar Putta
on 4 Mar 2021
Answered: Ravi Kumar
on 4 Mar 2021
I have (transient non-linear heat equation) i.e , k = 0.7 + 0.003*T, tlist = linespace(0,0.5,20); therefore my pde coefficients wills be d = 1, m = a = 0, c = k, f = 1; since c coefficent is non linear (dependent on T, temperature), My Intial Boundary condition is 3 and initial guess is 5. But i see command IC and initial guess is same. i.e setInitialConditions(model,u0); So how do i differentiate them?
This is how my code goes
clear,clc,close all;
model = createpde();
%% Geometry
R1 = [3;4;-1;1;1;-1;-1;-1;1;1];
g = decsg(R1);
heatmodel_geom = geometryFromEdges(model,g);
%% Mesh
msh = generateMesh(model,'GeometricOrder','linear');
[P,E,T] = meshToPet( msh );
%% Specify Coefficients
c = @(~,state) 0.7+0.003*state.u; % nonlinear coefficient
specifyCoefficients(model,'m',0,'d',1,'c',c,'a',0,'f',1);
%% Spacial BC
applyBoundaryCondition(model,'dirichlet','edge',[1 2 3 4],'u',0);
%% Initial guess and initial condition
setInitialConditions(model,3); % Initial condition
setInitialConditions(model,5); % Initial guess to solve non linear pde.... BUT DO I DIFFERENTIATE BOTH THIS COMMANDS TO ITS RESPECTIVE JOBS
0 Comments
Accepted Answer
Ravi Kumar
on 4 Mar 2021
For nonlinear elliptic problem, inputs in setInitialCcnditions is the initial guess for the nonlinear solver. For time-dependent problems initial condiions and initial gues are one and the same.
Regards,
Ravi
0 Comments
More Answers (0)
See Also
Categories
Find more on Geometry and Mesh in Help Center and File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!