Hi, everyone. I am starting to solve a chemomechanical problem, but I don't know whether the fully coupled equations can be solved in sfepy or not. Could you have a look at them?
The equation are as followed:
Governing Equations
∇∙σ+ρb=ρu ̈
Where σ is the stress tensor, ρ is the mass density, b is the body force per unit mass, u is the displacement vector.
∂c/∂t+∇∙j=∂r/∂t, where j=-(D*c_0)/RT ∇Ψ
Where c is the concentration, D is the diffusion coefficient, c0 is the reference concentration, …
[View More]R is Boltzmann constant (gas constant), T is the absolute temperature, Ψ is the chemical potential.
Constitutive relations
σ=2Gε+(λe+bc)I
Ψ=be+gc
Where b=-β(2G+3λ), β is the chemical dilation coefficient, e is the volumetric strain and g represents the scalar chemistry modulus.
[View Less]
Hello SfePy folks,
It’s been a long time (e.g., 10+ years) since I tinkered with SFEPY, but I dug out my old calculations and with a bit of tweaking, they still work. Yeah!
Thank you for keeping on for so long with this excellent tool for python!
I’m hoping you can lend me some guidance for a new calculation. I’m interested in using SfePy to answer questions about the behavior of PCB traces. For example, I’d like to be able to extract a trace from a Gerber file and convert it into a mesh …
[View More]that SfePy can grok. I can figure out the file conversion, I think. If nothing else I can fake it using the gmsh “geo” format. However, then I’d like to compute the current flow through the mesh with a fixed voltage difference between two boundaries while maintaining a zero current boundary condition on all the other boundaries. Is that a thing that could be done? Any suggestions on how to go about setting up the boundary conditions?
Thank you,
-Steve
[View Less]