Hello to everybody, I want to know if it is possible to use an output of a simulation inside another simulation. My problem is this: I’m trying to compute the symmetric part of the gradient of the velocity field u, coming from solving the stokes equation inside an obstacle with certain boundary co…

Thanks for the help. I’m trying to see that question but it’s different from me. What I want is to save a function (f_32_0) called u ( which is the velocity field given by stokes system ) inside the code. But I want to save exactly with the same format (f_32_0), not in xdmf or txt or something else,…

Could you provide a minimal example, as I do not understand what you are trying to do, or where you are obtaining an error.

from dolfin import * import numpy as np mesh = UnitCubeMesh(10, 10, 10) P1 = FiniteElement("Lagrange", mesh.ufl_cell(), 1) B = FiniteElement("Bubble", mesh.ufl_cell(), mesh.topology().dim() + 1) V = VectorElement(NodalEnrichedElement(P1, B)) Q = P1 W = FunctionSpace(mesh, V * Q) class Everywher…

Please make your code a minimal example, using built in meshes. Remove all unnecessary code to reproduce the error you are getting. As you have not supplied how you are currently trying to use the output w (or u, p) as input to your problem, it is not clear what error you are getting. As you have …

I’ve updated the question with MWE

Consider the minimal working code, where you do your operation for two different inputs, and compute the combined assembly: from dolfin import * import numpy as np class Everywhere(SubDomain): def inside(self, x, on_boundary): return True def simulation(W: FunctionSpace, boundary_con…

sorry but I’m not understanding. you write u1, p1 = simulation(W, as_vector((x[0], 2*x[1], 3*x[2]**2))) but I don’t really know why you should choose (x[0], 2x[1], 3x[2]**2). My u_1 must be equal to the u field solution obtained in the previously simulation, not be the same as any generic vector

The example was to illustrate how you should structure your code. You can easily replace it with: u1, p1 = simulation(W, as_vector((x[0], 2*x[1], 3*x[2]**2))) u2, p2 = simulation(W, u1)

I tried this approach : # Write `u` to a file: u1 = HDF5File(MPI.comm_world,"u_1.h5","w") u1.write(u,"u_1") u1.close() # Read the contents of the file back into a new function, `f1`: f1 = Function(W.sub(0).collapse()) u1 = HDF5File(MPI.comm_world,"u_1.h5","r") u1.read(f1, "u_1") u1.close() # Save…

Using output from other simulation

I/O

dokken January 27, 2022, 8:17pm 2

See: Loading xdmf data back in - #4 by dokken

Topic		Replies	Views
Saving an xdmf file - dolfinx v 0.7.0 dolfinx io	6	1680	February 26, 2024
Saving multiple functions in one file to visualize it I/O	15	3404	November 17, 2022
How do I save a custom grid？ General	3	257	September 6, 2023
Save solution and load it in the next iteration I/O	2	3486	April 13, 2019
Issue reading/writing functions from separate runs on the same mesh I/O dolfinx , mpi	8	89	November 25, 2025

Using output from other simulation

Related topics