udf.data['object_type']['spherical_particle']['particle_spec'][0][
'janus_slip_vel'] = 0.01
udf.data['object_type']['spherical_particle']['particle_spec'][0][
'janus_slip_mode'] = 50.0
udf.data['gravity']['g'] = 0.06
udf.data['switch']['init_distribution']['user_specify']['particles'][0]['q'][
'q0'] = q.q0
udf.data['switch']['init_distribution']['user_specify']['particles'][0]['q'][
'q1'] = q.q1
udf.data['switch']['init_distribution']['user_specify']['particles'][0]['q'][
'q2'] = q.q2
udf.data['switch']['init_distribution']['user_specify']['particles'][0]['q'][
'q3'] = q.q3
udf.data['output']['num_snap'] = 800
gammadots = np.arange(0.005, 0.055, 0.005)
if not Path(f'udf').exists():
os.mkdir(f'udf')
for gammadot in gammadots:
gammadot_str = f'{gammadot * 1000:0=3.0f}'
udf.data['constitutive_eq']['shear_navier_stokes']['external_field']['dc'][
'shear_rate'] = gammadot
os.mkdir(f'udf/{gammadot_str}')
udf.to_udf(f'udf/{gammadot_str}/input.udf')
shutil.copy('/Users/taiga/Projects/lab/simulation/define.udf',
f'udf/{gammadot_str}')
import shutil
import numpy as np
from pathlib import Path
from pyudf import Udf
from pyudf.rotation import z_rotation
udf = Udf()
udf.load_jsonconfig()
_degrees = [d * 2 * np.pi / 32 for d in range(1, 33)]
_degrees_str = [f'{i:0=2}' for i in range(1, 33)]
if not Path('udf').exists():
os.mkdir('udf')
for d, d_str in zip(_degrees, _degrees_str):
q = z_rotation(d)
udf.data['switch']['init_distribution']['user_specify']['particles'][0][
'q']['q0'] = q.q0
udf.data['switch']['init_distribution']['user_specify']['particles'][0][
'q']['q1'] = q.q1
udf.data['switch']['init_distribution']['user_specify']['particles'][0][
'q']['q2'] = q.q2
udf.data['switch']['init_distribution']['user_specify']['particles'][0][
'q']['q3'] = q.q3
os.mkdir(f'udf/{d_str}')
udf.to_udf(f'udf/{d_str}/input.udf')
shutil.copy('/Users/taiga/Projects/lab/simulation/define.udf',
f'udf/{d_str}/')