The Math Forum

Search All of the Math Forum:

Views expressed in these public forums are not endorsed by NCTM or The Math Forum.

Math Forum » Discussions » Software » comp.soft-sys.matlab

Notice: We are no longer accepting new posts, but the forums will continue to be readable.

Topic: Linear system of Hydrodynamic solving by means of PDE toolbox
Replies: 13   Last Post: Jan 16, 2014 10:05 PM

Advanced Search

Back to Topic List Back to Topic List Jump to Tree View Jump to Tree View   Messages: [ Previous | Next ]

Posts: 1,717
Registered: 11/8/10
Re: Linear system of Hydrodynamic solving by means of PDE toolbox
Posted: Jan 8, 2014 9:07 AM
  Click to see the message monospaced in plain text Plain Text   Click to reply to this topic Reply

"Bill Greene" wrote in message <lajkp6$mce$>...
> The answer to your second question is to call the pdegrad function with
> the solution vector.
> Regarding your first question, where are you stuck? If you post the code you
> have so far and describe the point of your confusion, we might be able to help.
> Bill
> "Sergey" wrote in message <lain4d$9uh$>...

> > Alright, If there is no idea how the problem can be solved directly (I've stuck as well) Maybe someone could tell me how to get the Phi_z (= dPhi / dz) from the solution of the Laplace equation of PDE toolbox? It gives result on the unstructured grid, that I can't directly differentiate by the usual scheme like (Phi(2,1) - Phi(1,1))/dz.
> >
> > "Sergey" wrote in message <la2vkf$rrb$>...

> > > Dear Matlab users and developers,
> > >
> > > I need to solve the following PDE system in 2D space:
> > > (1) div(grad( Phi )) = 0 |any x, any z;
> > > (2) Eta_t = Phi_z | any x, z = 0;
> > > (3) Phi_t = -g*Eta | any x, z = 0.
> > > where "?_t" - time derivative, "?_z" - z - axis derivative, "div" - divergence operator, "grad" - gradient operator.
> > > (1) - simple Laplace equation for interior flow potential,
> > > (2) - surface elevation (Eta) and the potential (Phi) connection,
> > > (3) - linearized Bernoulli equation.
> > > It is the linearized equations of hydrodynamics. Of course I could solve it using just a FD scheme, but if there is a possibility to utilize some internal Matlab tools like PDE toolsbox, it is preferable to use, since it is already developed and verified.
> > > I know, that the PDE toolbox provide a possibility to solve the 2D Laplace equation. What about the rest of the system? Can I somehow couple all the equations by means of the PDE toolbox functionality?
> > >
> > > Sincerely,
> > > -Sergey

The OP needs the gradient at the boundary surface z=0.
This is also possible with pdegrad ? What about the accuracy ?

Best wishes

Point your RSS reader here for a feed of the latest messages in this topic.

[Privacy Policy] [Terms of Use]

© The Math Forum at NCTM 1994-2018. All Rights Reserved.