simtk.org

Code: Select all

params_db['...']

Code: Select all

build_folder/SimVascular-build/sv --python -- your_script.py

Code: Select all

import sv

# create 1d simulation
oned_simulation = sv.simulation.OneDimensional()

# Create 1D simulation parameters
params = sv.simulation.OneDimensionalParameters()

# Mesh parameters
mesh_params = params.MeshParameters()

# Model parameters
model_params = params.ModelParameters()
model_params.name = geo
model_params.inlet_face_names = ['inflow']
model_params.outlet_face_names = params_db['outlets']
model_params.centerlines_file_name = params_db['centerlines_input_file']

# Fluid properties
fluid_props = params.FluidProperties()

# Set wall properties
material = params.WallProperties.OlufsenMaterial()

# Set boundary conditions
bcs = params.BoundaryConditions()

# add inlet
bcs.add_velocities(face_name='inlet', file_name=params_db['inflow_input_file'])

# add outlets
for cp in params_db['outlets']:
    t = params_db['bc_def']['bc_type'][cp]
    bc = params_db['bc_def']['bc'][cp]
    if 'Po' in bc and bc['Po'] != 0.0:
        raise ValueError('RCR reference pressure not implemented')
    if t == 'rcr':
        bcs.add_rcr(face_name=cp, Rp=bc['Rp'], C=bc['C'], Rd=bc['Rd'])
    elif t == 'resistance':
        bcs.add_resistance(face_name=cp, resistance=bc['R'])
    elif t == 'coronary':
        raise ValueError('Coronary BC not implemented')

# Set solution parameters
solution_params = params.Solution()
solution_params.time_step = params_db['dt']
solution_params.num_time_steps = params_db['num_dts']
solution_params.save_frequency = params_db['save_data_freq']

# Write a 1D solver input file
oned_simulation.write_input_file(model=model_params, mesh=mesh_params, fluid=fluid_props,
                                 material=material, boundary_conditions=bcs, solution=solution_params,
                                 directory=params_db['fpath_1d'])