grad(3rheolef) rheolef-6.1 grad(3rheolef)NAME
grad -- gradient operator
SYNOPSIS
form(const space V, const space& M, "grad");
DESCRIPTION
Assembly the form associated to the gradient operator on a finite element space V:
/
|
b(u, q) = | grad(u).q dx
|
/ Omega
The V space may be a either P1 or P2 finite element space, while the M space may be P0 or P1d respectively. See also form(2) and
space(2).
EXAMPLE
The following piece of code build the divergence form associated to the P1 approximation:
geo omega("square");
space V(omega, "P1");
space M(omega, "P0", "vector");
form b(V, M, "grad");
SEE ALSO form(2), space(2)rheolef-6.1 rheolef-6.1 grad(3rheolef)
Check Out this Related Man Page
s_grad_grad(3rheolef) rheolef-6.1 s_grad_grad(3rheolef)NAME
s_grad_grad -- grad_grad-like operator for the Stokes stream function computation
SYNOPSIS
form(const space V, const space& V, "s_grad_grad");
DESCRIPTION
Assembly the form associated to the -div(grad) variant operator on a finite element space V. The V space may be a either P1 or P2 finite
element space. See also form(2) and space(2). On cartesian coordinate systems, the form coincide with the "grad_grad" one (see
grad_grad(3)):
/
|
a(u,v) = | grad(u).grad(v) dx
|
/ Omega
The stream function on tri-dimensionnal cartesian coordinate systems is such that
u = curl psi
div psi = 0
where u is the velocity field. Taking the curl of the first relation, using the identity:
curl(curl(psi)) = -div(grad(psi)) + grad(div(psi))
and using the div(psi)=0 relation leads to:
-div(grad(psi)) = curl(u)
This relation leads to a variational formulation involving the the "grad_grad" and the "curl" forms (see grad_grad(3), curl(3)).
In the axisymmetric case, the stream function psi is scalar ans is defined from the velocity field u=(ur,uz) by (see Batchelor, 6th ed.,
1967, p 543):
d psi d psi
uz = (1/r) ----- and ur = - (1/r) -----
d r d r
See also http://en.wikipedia.org/wiki/Stokes_stream_function . Multiplying by rot(xi)=(d xi/dr, -d xi/dz), and integrating with r dr dz,
we get a well-posed variationnal problem:
a(psi,xi) = b(xi,u)
with
/
| (d psi d xi d psi d xi)
a(psi,xi) = | (----- ---- + ----- ----) dr dz
| ( d r d r d z d z )
/ Omega
and
/
| (d xi d xi )
b(xi,u) = | (---- ur - ---- uz) r dr dz
| (d z d r )
/ Omega
Notice that a is symmetric definite positive, but without the 'r' weight as is is usual for axisymmetric standard forms. The b form is
named "s_curl", for the Stokes curl variant of the "curl" operator (see s_curl(3)) as it is closely related to the "curl" operator, but
differs by the r and 1/r factors, as:
( d (r xi) d xi )
curl(xi) = ( (1/r) -------- ; - -----)
( d r d z )
while
( d xi d xi )
s_curl(xi) = ( ---- ; - ---- )
( d r d z )
EXAMPLE
The following piece of code build the form associated to the P1 approximation:
geo g("square");
space V(g, "P1");
form a(V, V, "s_grad_grad");
SEE ALSO form(2), space(2), grad_grad(3), grad_grad(3), curl(3), s_curl(3)rheolef-6.1 rheolef-6.1 s_grad_grad(3rheolef)
Man Page
12 More Discussions You Might Find Interesting
1. Post Here to Contact Site Administrators and Moderators
ive been using Ubuntu 7.10 for about 3 weeks now and i love it more than any other linux operating system. ive been watching the countdown to the upgrade of ubuntu 8.04. has any one tried out the 8.04 beta at all. if so, did u think it is worth upgrading my box to 8.04 or should i keep it at... (13 Replies)
I need to write a script that will take the input from a file and convert the number from centigrade to fahrenheit and vice versa.
This is what I have but it doesn't seem to be correct.
Also the data file has 11 numbers inside of it and the output needs to be listed as so:
Fahrenheit Temperature... (18 Replies)
I am writing a software product and hope that it will work on a variety of Linux distributions. At the moment, I am trying to create some kind of Linux version of patches/upgrades of installed software. Gathering information on available updates isn't hard, nor is installation of updates, but I... (27 Replies)
Hi All,
We have patched our solaris 10 u4 server with lustarter patch in order to get it prepared for solaris live upgrade.The patch installation went pretty fine and we did reboot the server after patching.But after patching we have noticed O.S is not able to detect tape library.Not sure what... (22 Replies)
Dear friends,
i am having some Problems with a Sco Openserver in a Box (normally on VMWare, importing it on a Virtualbox does the same).
When i boot up the system on xx/09/2015 i got fork failed - too many processes during startup.
The programs which should run on the server, do not run.
I... (67 Replies)
Hello everybody,
I upgraded from Ubuntu 12.04 to 14.04, so some of my scripts in bash does NOT work, for example:
#!/bin/bash
# Some codes goes here
for ((offset=$thumbStart; offset<=$thumbEnd; offset+=10)) # LINE 459
do
printf -v outfile... (16 Replies)
Hello,
Long story short, we built an alternate boot env. back in July and applied the latest CPU to about 15 servers. Of the 15 servers, 7 servers failed to start the zones after the luactivate / reboot. The zones failed to rename from <zone>-<boot_env> back to <zone>. This is fixed in patch... (16 Replies)
Hey,
Upgrade (step 1) the posts, by putting a "user info" icon in the top right and making it so it toggles the user info.
The user info icon has a badge which shows the number of posts. I will have to dig around in the code more to get the total posts thanks and other badges working; but... (22 Replies)
Three days ago we received an expected notice from our long time data center that they were going dark on Sept 12th.
About one and a half hours ago, after three days of marathon work, I just cut over the unix.com to a new data center with a completely new OS and Ubuntu distribution. (22 Replies)