## solution first order PDEs with Ics and Bcs ...

Asked by:

Dear all

I need a help how can i solve for example the following PDEs with Initial condition and boundary condition given at x=-1, and x=1.

pde:=diff(u(t,x),t)=diff(u(t,x),x\$2);

ics:=u(0,x)=sin(x);

Bcs:=diff(u(t,-1),x)=0;

Bcs:=diff(u(t,1),x)=0;

Many thanks

## Solving system of Pdes...

Asked by:

Hello,

How can we write a code for solving following system of Pdes ?

```restart;
PDE1:=m1*diff(w1(x,t),t\$2)-S1*diff(w1(x,t),x\$2)+k*(w1(x,t)-w2(x,t))=F1(t)*delta(x-x1);
PDE2:=m2*diff(w2(x,t),t\$2)-S2*diff(w2(x,t),x\$2)+k*(w2(x,t)-w1(x,t))=F2(t)*delta(x-x2);
```

In here,  Ki, mi, xi and k are constants where i=1,2.

Since F1(t) and F2(t) are unspecified (ungiven) functions, solutions u1,u2 which we seek, will be depended on  F1(t) and F2(t).

Thanks for your valuable and praiseworthy suggestions and comments.

## parial derivative PDES...

Asked by:

Dear all;

I have a function  C(S,t)  and I would like to show that this function is the solution of Partial differential equation ( denotes by PDES). I introduce all my function but when I try to simplify my PDES, I can not get zero as result.

Many thanks for your help to simplify the PDES to get PDES=0. So, that C(S,t) is the solution of my PDEs.

simplification.mw

## Solve nonlinear system of equations ...

Asked by:

Hi

I have a nonlinear PDEs, solved using finite difference in the square

I get the following nonlinear system of equation. Is there any idea how correct the code and display the solution.

I will appreciate any help in this question.

restart;
n:=100;
h:=1/(n+1);

# Boundary condition

for j from 0 by 1  to n+1 do
u[0, j] = 0;
u[n+1, j] = 0;
u[j, 0] = 0;
u[j, n+1] = 0 ;
end do;
## Loop for interior point in the square
for i from 1 by 1 to  n do
for j from 1 by 1 to  n do
(u[i+1, j]-u[i, j])*(u[i+1, j]-2*u[i, j]+u[i-1, j])+h*(u[i, j+1]-2*u[i, j]+u[i, j-1]) = 0;
end do;
end do;

How can I solve this system of equations with unknown u[i,j], where i,j=1,..,n

Many thanks for any help

## How do I solve a system with boundary conditions t...

Asked by:

Hello everyone,

I am having trouble trying to solve a system of differential equations. The modeling was made from the equilibrium equations of a pressure vessel. The set of equations is shown below:

As you see it is a set of two second-order partial differential equations. So, we need four boundary conditions. This one is the first. It means that the left end of the pressure vessel is fixed.

This one is the second boundary condition. It means that the right end of the pressure vessel is free.

This one is the third boundary condition. It means that the inner surface of the pressure vessel is subject to an external load:

At last, we have the fourth boundary condition. It means that the outer surface of the pressure vessel is free.

The first test I have been trying to do is the static case. In this case, the time terms (the right side of the two equations shown) is zero.

The maple commands that I am using are the following:

```restart; E := 200*10^9; nu := .33; G := E/(2*(1+nu)); RI := 0.254e-1; RO := 2*RI; p := proc (x) options operator, arrow; 50000000 end proc; sys := [E*(nu*(diff(v(x, r), x))/r+nu*(diff(diff(v(x, r), x), r))+(1-nu)*(diff(diff(u(x, r), x), x)))/(-2*nu^2-nu+1)+G*(diff(diff(u(x, r), r), r)+diff(diff(v(x, r), x), r)+(diff(u(x, r), r))/r+(diff(v(x, r), x))/r) = 0, E*((1-nu)*(diff(diff(v(x, r), r), r))+nu*(diff(diff(u(x, r), x), r))+(1-nu)*(diff(v(x, r), r))/r-(1-nu)*v(x, r)/r^2)/(-2*nu^2-nu+1)+G*(diff(diff(u(x, r), r), x)+diff(diff(v(x, r), x), x)) = 0]; BCs := {E*(nu*v(L, r)/r+nu*(D[2](v))(L, r)+(1-nu)*(D[1](u))(L, r))/(-2*nu^2-nu+1) = 0, E*(nu*v(x, RI)/RI+(1-nu)*(D[2](v))(x, RI)+nu*(D[1](u))(x, RI))/(-2*nu^2-nu+1) = -p(x), E*(nu*v(x, RO)/RO+(1-nu)*(D[2](v))(x, RO)+nu*(D[1](u))(x, RO))/(-2*nu^2-nu+1) = 0, u(0, r) = 0}

sol := pdsolve(sys, BCs, numeric)```

I am getting the following error:

Error, (in pdsolve/numeric/process_IBCs) initial/boundary conditions must depend upon exactly one of the independent variables: 0.1459531181e12*v(L, r)/r+0.1459531181e12*(D[2](v))(L, r)+0.2963290579e12*(D[1](u))(L, r) = 0

In this case, my boundary conditions do depend on more than one independent variable. How do I proceed?

Thank you in advance,

Pedro Guaraldi

## Solve PDEs boundary condition ...

Asked by:

Dear all ;

I have a Partial differential equation

restart; with(PDEtools);

pde[2] := (diff(u(x, y), x))*(diff(u(x, y), x, x))+diff(u(x, y), y, y);
where x and y in the square [0,1]

with boundary condition

bc[2] := u(0, y) = 0, u(1, y) = 0, u(x, 0) = 0, u(x, 1) = 0;

Is there a simple code to compute the solution

Many thanks for any help

## How do I obtain a simplified system automatically...

Asked by:

HI everyone,

As can be seen from the attached file, the first three equations of Eq. (5) will render some of the other equations (and other terms) redundant. How can I obtain a simplified system automatically?

Thanks.

Pdesample.mw

## Problem to solve system of linear PDEs with BSc an...

Asked by:

Hello

Hope everything going fine with you. I am facing problem to fine the exact (numerical) solution of the attached system of linear PDEs associated with BSc and ICs. I tried to solve it without BCs and ICs, with BCs and with ICs also all the times I failed. Please solve it either general, with ICs or BCs. You can try to solve it numrically. In attached file H(t) represent the unit step function. I am waiting your positive response.

With my best regards and sincerely.

Muhammad Usman

School of Mathematical Sciences
Peking University, Beijing, China

## New developments on exact solutions for PDEs with...

by: Maple

Hi

New developments (after the release of Maple 2016) happened in the project on exact solutions for "Partial Differential Equations & Boundary Conditions". This is work in collaboration with Katherina von Bulow and the improvements are of wide range, representing a noticeable step forward in the capabilities of the Maple system for this kind of problem. As usual, these improvements can be installed in current Maple 2016 by downloading the updated library from the Maplesoft R&D webpage for Differential Equations and Mathematical functions (the update is distributed merged with the updates of the Physics package)

The improvements cover:

 • PDE&BC in semi-infinite domains for which a bounded solution is sought
 • PDE & BC problems in bounded spatial domains via eigenfunction (Fourier) expansions
 • Implementation of another algebraic method for tackling linear PDE & BC
 • Improvements in solving PDE & BC solutions by first finding the PDE's general solution.
 • Improvements in solving PDE & BC problems by using a Fourier transform.
 • PDE & BC problems that used to require the option HINT = `+` are now solved automatically

What follows is a set of examples solved now with these new developments, organized in sections according to the kind of problem. Where relevant, the sections include a subsection on "How it works step by step".

PDE&BC in semi-infinite domains for which a bounded solution is sought can now also be solved via Laplace transforms

Maple is now able to solve more PDE&BC problems via Laplace transforms.

How it works: Laplace transforms act to change derivatives with respect to one of the independent variables of the domain into multiplication operations in the transformed domain. After applying a Laplace transform to the original problem, we can simplify the problem using the transformed BC, then solve the problem in the transformed domain, and finally apply the inverse Laplace transform to arrive at the final solution. It is important to remember to give pdsolve any necessary restrictions on the variables and constants of the problem, by means of the "assuming" command.

A new feature is that we can now tell pdsolve that the dependent variable is bounded, by means of the optional argument HINT = boundedseries.

 >

Consider the problem of a falling cable lying on a table that is suddenly removed (cf. David J. Logan's Applied Partial Differential Equations p.115).

 >

If we ask pdsolve to solve this problem without the condition of boundedness of the solution, we obtain:

 >
 (1.1)

New: If we now ask for a bounded solution, by means of the option HINT = boundedseries, pdsolve simplifies the problem accordingly.

 >
 (1.2)

And we can check this answer against the original problem, if desired:

 >
 (1.3)

How it works, step by step

Let us see the process this problem undergoes to be solved by pdsolve, step by step.

First, the Laplace transform is applied to the PDE:

 >
 >
 (1.1.1)

and the result is simplified using the initial conditions:

 >
 (1.1.2)

Next, we call the function "laplace(u(x,t),t,s)" by the new name U:

 >
 (1.1.3)

And this equation, which is really an ODE, is solved:

 >
 (1.1.4)

Now, since we want a BOUNDED solution, the term with the positive exponential must be zero, and we are left with:

 >
 (1.1.5)

Now, the initial solution must also be satisfied. Here it is, in the transformed domain:

 >
 (1.1.6)

Or, in the new variable U,

 >
 (1.1.7)

And by applying it to bounded_solution_U, we find the relationship

 >
 (1.1.8)
 >
 (1.1.9)

so that our solution now becomes

 >
 (1.1.10)

to which we now apply the inverse Laplace transform to obtain the solution to the problem:

 >
 (1.1.11)

Four other related examples

A few other examples:

 >
 >
 (1.2.1)
 >
 (1.2.2)
 >
 >
 (1.2.3)
 >
 (1.2.4)
 >
 >
 (1.2.5)
 >
 (1.2.6)

The following is an example from page 76 in Logan's book:

 >
 >
 (1.2.7)

More PDE&BC problems in bounded spatial domains can now be solved via eigenfunction (Fourier) expansions

The code for solving PDE&BC problems in bounded spatial domains has been expanded. The method works by separating the variables by product, so that the problem is transformed into an ODE system (with initial and/or boundary conditions) problem, one of which is a Sturm-Liouville problem (a type of eigenvalue problem) which has infinitely many solutions - hence the infinite series representation of the solutions.

 >

Here is a simple example for the heat equation:

 >
 >
 (2.1)
 >
 (2.2)

Now, consider the displacements of a string governed by the wave equation, where c is a constant (cf. Logan p.28).

 >
 >
 (2.3)
 >
 (2.4)

Another wave equation problem (cf. Logan p.130):

 >
 >
 (2.5)
 >
 (2.6)

Here is a problem with periodic boundary conditions (cf. Logan p.131). The function  stands for the concentration of a chemical dissolved in water within a tubular ring of circumference . The initial concentration is given by , and the variable  is the arc-length parameter that varies from 0 to .

 >
 >
 (2.7)
 >
 (2.8)

The following problem is for heat flow with both boundaries insulated (cf. Logan p.166, 3rd edition)

 >
 >
 (2.9)
 >
 (2.10)

This is a problem in a bounded domain with the presence of a source. A source term represents an outside influence in the system and leads to an inhomogeneous PDE (cf. Logan p.149):

 >
 >
 (2.11)

Current pdetest is unable to verify that this solution cancells the  mainly because it currently fails in identifying that there is a fourier expansion in it, but its subroutines for testing the boundary conditions work well with this problem

 >
 >
 (2.12)

Consider a heat absorption-radiation problem in the bounded domain :

 >
 >
 (2.13)
 >
 (2.14)

Consider the nonhomogeneous wave equation problem (cf. Logan p.213, 3rd edition):

 >
 >
 (2.15)
 >
 (2.16)

Consider the following Schrödinger equation with zero potential energy (cf. Logan p.30):

 >
 >
 (2.17)
 >
 (2.18)

Another method has been implemented for linear PDE&BC

This method is for problems of the form

or

where M is an arbitrary linear differential operator of any order which only depends on the spatial variables .

Here are some examples:

 >
 >
 (3.1)
 >
 (3.2)

Here are two examples for which the derivative with respect to t is of the second order, and two initial conditions are given:

 >
 >
 (3.3)
 >
 (3.4)
 >
 >
 (3.5)
 >
 (3.6)

More PDE&BC problems are now solved via first finding the PDE's general solution.

The following are examples of PDE&BC problems for which pdsolve is successful in first calculating the PDE's general solution, and then fitting the initial or boundary condition to it.

 >

If we ask pdsolve to solve the problem, we get:

 >
 (4.1)

and we can check this answer by using pdetest:

 >
 (4.2)

How it works, step by step:

The general solution for just the PDE is:

 >
 (4.1.1)

Substituting in the condition , we get:

 (4.1.2)
 >
 (4.1.3)

We then isolate one of the functions above (we can choose either one, in this case), convert it into a function operator, and then apply it to gensol

 >
 (4.1.4)
 >
 (4.1.5)

Three other related examples

 >
 >
 (4.2.1)
 >
 (4.2.2)
 >
 >
 (4.2.3)
 >
 (4.2.4)
 >
 >
 (4.2.5)
 >
 (4.2.6)

More PDE&BC problems are now solved by using a Fourier transform.

 >

Consider the following problem with an initial condition:

 >

pdsolve can solve this problem directly:

 >
 (5.1)

And we can check this answer against the original problem, if desired:

 >
 (5.2)

How it works, step by step

Similarly to the Laplace transform method, we start the solution process by first applying the Fourier transform to the PDE:

 >
 >
 (5.1.1)

Next, we call the function "fourier(u(x,t),x,s1)" by the new name U:

 >
 (5.1.2)

And this equation, which is really an ODE, is solved:

 >
 (5.1.3)

Now, we apply the Fourier transform to the initial condition :

 (5.1.4)
 >
 (5.1.5)

Or, in the new variable U,

 >
 (5.1.6)

Now, we evaluate solution_U at t = 0:

 >
 (5.1.7)

and substitute the transformed initial condition into it:

 >
 (5.1.8)

Putting this into our solution_U, we get

 >
 (5.1.9)

Finally, we apply the inverse Fourier transformation to this,

 >
 (5.1.10)

PDE&BC problems that used to require the option HINT = `+` to be solved are now solved automatically

The following two PDE&BC problems used to require the option HINT = `+` in order to be solved. This is now done automatically within pdsolve.

 >
 >
 (6.1)
 >
 (6.2)
 >
 >
 (6.3)
 >
 (6.4)
 >

Download PDE_and_BC_update.mw

Edgardo S. Cheb-Terrab
Physics, Differential Equations and Mathematical Functions, Maplesoft

## How do I solve following PDE system?...

Asked by:

I have the following PDE system steaming from Flash Photolysis:

when I use pdsolve(pdesys,[J,C]) I get:

The solution appears to be either C(x,t) = 0 or J(x,t) = 0. This are the obvious solutions (0 = 0). I have the analytical solution to this PDE system where neither C(x,t) nor J(x,t) are 0.

How to solve this system in maple? Thanks.

## Numerical solution of PDE system ...

Asked by:

I have a system of PDEs with 3 dependent variables (U, V, W) and 3 independent variables (x, y, t).

I need an explicit numerical solution, i.e., a table with the values of U, V and W, corresponding to several values of x, y and t.

Can anybody help me?

## Two independent variables...

Asked by:

HI, I am trying to solve two PDEs but in boundry conditions there is arising an error plz help.
Nazi.mw

## Solution to any PDE system with initial or boundar...

Asked by:

Hello All,

I looked through the Maple help on PDE systems and pdsolve and the physics problems that appear there. THere are a number of single-PDE cases with initial / boundary conditions; but I couldn't find PDE systems with ics/bcs.

Would you have a (simple) example of a PDE system with its initial / boundary conditions? I am attempting to build understanding of the syntax and different options of "pdsolve". Examples seem to be a great way to learn how to solve PDE systems. One can then pdsolve the PDE system without ics; then add them, try different options etc.

Thank you!

## Problem with PDEs ...

Asked by:

I'm trying to build a Maple procedure that will generate vector fields on a metric with certain properties. Working with metric g over the coordinates {u,v,w}, call the field X = (a(u,v,w), b(u,v,w), c(u,v,w)). The field should satisfy <X, X> = 0 and have the directional covariant derivative of X in the direction of each coordinate vector field = 0 (with resepct to the Levi-Civita conenction).

Basically, these conditions yield a system of 3 PDEs and an algebraic expressionin terms of a,b,c. I've been trying to solve them using pdsolve, but I'm getting the error message:

>Error, (in pdsolve/sys) the input system cannot contain equations in the arbitrary parameters alone; found equation depending only on _F1(u,v,w): _F1(u,v,w)

I've attached my worksheet. Can anyone help me out?

Thanks! ppwaves.mw

## How can PDE be solved in terms of Bessel function?...

Asked by:

Good day everyone,

please how can one solve this pde in terms of Bessel function or any other analytic solution with the plot.

See the file ID.mw

Thanks.

 1 2 Page 1 of 2
﻿