NUTILS 8.1 "Idiyappam"

arguments

L=2.5,2.5

R=0.4,0.7

mu=1

beta=100

gammaskeleton=0.05

gammaghost=0.0005

pbar=1

nelems=12,12

degree=3

maxrefine=2

seed=123456

plotting=matplotlib

start 2023-08-15 02:16:40

pre-processing

NutilsDeprecationWarning: function.norm2 is deprecated; use numpy.linalg.norm instead In /usr/local/lib/python3.11/dist-packages/nutils/function.py:1226

Inclusion at [1.25 1.25] with radius 0.59.

NutilsDeprecationWarning: function.min is deprecated; use numpy.minimum instead In /usr/local/lib/python3.11/dist-packages/nutils/function.py:1226

Inclusion at (0, 0) with radius 0.52.

Inclusion at (0, 2.5) with radius 0.41.

Inclusion at (2.5, 0) with radius 0.69.

Inclusion at (2.5, 2.5) with radius 0.67.

porosity: 0.6541

immersed_domain.png

skeleton_mesh.png

ghost_mesh.png

immersed iga solver

solve

solving 624 dof system to machine precision using arnoldi solver

constructing direct preconditioner

peak memory use 5,034k

residual decreased by 14.0 orders using 0 krylov vectors

residual decreased by 0.2 orders using 1 krylov vectors

residual decreased by 0.1 orders using 2 krylov vectors

solver returned with residual 3e-15

post-processing

velocity.png

pressure.png

(in/out)flow area : 1.5702311197916665 / 1.141560872395833

(in/out)flow pressure: 1.0039625064394604 / 0.0015642550295126747

(in/out)flow flux : 0.017499848135085476 / 0.01749984813499137

(in/out)flow velocity: 0.011144759465350108 / 0.015329754687775728

finish 2023-08-15 02:16:51, elapsed 0:11