Hello Fenics community, I am solving Stokes equations in steady state and I want to extract the components of the velocity field for post processing, I have been looking and trying solutions for hours now but to no avail. The doc in this regard is very sparse and sometimes outdated. Here is what I …

There are several ways of doing this. First, you can split the fields, Lets say you have the combined vector up, containing the Mixed space for velocity and pressure. Then the first component can be obtained by: u0=up.sub(0).sub(0,deepcopy=True) Note that if you just save the velocity field to fi…

Thanks Doken, I already splitted the mixed space and have access to u and p and I saved them as .pvd and I am able to play with them on paraview. However, I would like to do post-processing in python once the solution is computed and access the x and y components of the velocity field.

[image] dokken: u0=up.sub(0).sub(0,deepcopy=True) Then do as sketched above or ux=u.sub(0,deepcopy=True) uy=u.sub(1,deepcopy=True) if u is the velocity field resulting from the split.

Thanks, I will try your solution dokken. In the mean time, I found an (ugly) work around where I save the velocity field as pvd then use meshio to reload the file and access the different components of the field from the x,y positions of the nodes to the ux,uy components of the velocity.

Hello dokken! If I have two scalar functions ux and uy, one vector function u, how to replace the components of u with ux and uy? I have tried temp = u.sub(0) temp = ux It seems u is not modified. Can deepcopy influence function other than mixed element function space?

Access to the velocity field components

I/O

dokken November 18, 2020, 4:18pm 8

Use a FunctionAssigner, see Interpolation of Indexed Functions of Mixed Elements

Topic		Replies	Views
Save results in text file I/O	2	703	April 26, 2022
How to extract/operate on individual components of a VectorFunctionSpace Uncategorized	2	901	September 17, 2021
Writing vector filed components in XDMF I/O	5	1409	June 4, 2021
Convert FEniCS solution and coordinates to (numpy) 2D arrays I/O	9	4510	February 16, 2023
How to get the solution obtained by fenics as a vector? Uncategorized	13	4372	June 15, 2022

Access to the velocity field components

Related topics