Esempi in Python per simulation

Esempio n. 1

File: test_simulation.py Progetto: pyguerder/hedge-elastodynamics

 def test_simulation_2d_linear_plm(self):
     self.clean_folder()
     simulation(write_output=[],
          allow_features=['mpi'],
          dim=2,
          order=4,
          stfree_tag=TAG_NONE,
          fix_tag=TAG_ALL,
          op_tag=TAG_NONE,
          flux_type="lf",
          max_steps=40,
          output_dir=SimulationTestCase.folder,
          pml=[400, 400, 0, 400, 400, 0],
          sources=numpy.array([200.0,0.0]),
          source_param={'type':'sinus', 'sigma':10, 'fc':7.25, 'td':0.16, 'begin':0, 'end':2},
          final_time=12,
          quiet_output=True,
          nonlinearity_type=None,
          mesh_file = 'Meshes/BiHeterogeneousPeriodicSquare.msh',
          periodicity = [None, None],
          material_files = ['Materials/aluminium.dat', 'Materials/calcite.dat', 'Materials/SiO2.dat'],
          vtu_every=20)
     assert True

Esempio n. 2

File: resonator_3_cells_1d.py Progetto: pyguerder/hedge-elastodynamics

__license__ = "GNU GPLv3 (or more recent equivalent)"


import numpy
from hedge.mesh import TAG_ALL, TAG_NONE

from elastic_wave import main as simulation

simulation(
    write_output=["vtu", "receivers"],
    allow_features=["mpi"],
    dim=1,
    order=4,
    stfree_tag=TAG_ALL,
    fix_tag=TAG_NONE,
    op_tag=TAG_NONE,
    flux_type="lf",
    max_steps=None,
    output_dir="output",
    pml=None,  # [1, 0, 0, 1, 0, 0],
    sources=numpy.array([-7]),
    source_param={"type": "sinus", "sigma": 0.01, "fc": 722.0, "td": 1 / 722.0, "begin": 0, "end": 2 / 722.0},
    final_time=1,
    quiet_output=True,
    nonlinearity_type=None,
    mesh_file="Meshes/Resonator1D.msh",
    periodicity=[None],
    material_files=["Materials/epoxy.dat", "Materials/aluminium.dat", "Materials/aluminium.dat"],
    vtu_every=20,
)

Esempio n. 3

File: mixing_simple.py Progetto: pyguerder/hedge-elastodynamics

__license__ = "GNU GPLv3 (or more recent equivalent)"


from hedge.mesh import TAG_ALL, TAG_NONE

from elastic_wave import main as simulation

import numpy

simulation(write_output=['vtu', 'receivers'],
           allow_features=['mpi', 'cuda'],
           dim=2,
           order=5,
           stfree_tag=TAG_NONE,
           fix_tag=TAG_ALL,
           op_tag=TAG_NONE,
           flux_type="lf",
           max_steps=75e3,
           output_dir='output',
           pml=[1.5, 0, 0, 1.5, 0, 0],
           sources=numpy.array([-20.0, 0.0]),
           source_param={'type': 'SineBurst', 'sigma':1, 'fc':200, 'td':0.0, 'begin':0, 'end':0.1},
           final_time=12,
           quiet_output=True,
           nonlinearity_type='cubic',
           mesh_file = 'Meshes/Mixing2D_simple.msh',
           periodicity = [None, ('minus_y', 'plus_y')],
           material_files = ['Materials/PDMS.dat', 'Materials/PDMS.dat'],
           vtu_every=100)

Esempio n. 4

File: mixing.py Progetto: pyguerder/hedge-elastodynamics

__license__ = "GNU GPLv3 (or more recent equivalent)"


from hedge.mesh import TAG_ALL, TAG_NONE

from elastic_wave import main as simulation

import numpy

simulation(write_output=['vtu', 'receivers'],
           allow_features=['mpi', 'cuda'],
           dim=2,
           order=3,
           stfree_tag=TAG_NONE,
           fix_tag=TAG_ALL,
           op_tag=TAG_NONE,
           flux_type="lf",
           max_steps=1e6,
           output_dir='Mixing2D_rand_heterogeneous_5e4_cuda_cq_40Hz',
           pml=[5, 0, 0, 5, 0, 0],
           sources=numpy.array([-19, 0.0]),
           source_param={'type': 'SineBurst', 'sigma':1, 'fc':40, 'td':0.0, 'begin':0, 'end':0.4},
           final_time=1,
           quiet_output=True,
           nonlinearity_type='cubic',
           mesh_file = 'Meshes/Mixing2D_rand.msh',
           periodicity = [None, ('minus_y', 'plus_y')],
           material_files = ['Materials/PDMS.dat', 'Materials/polymer.dat'],
           vtu_every=1000)

Esempio n. 5

File: resonator.py Progetto: pyguerder/hedge-elastodynamics

__copyright__ = "Copyright (C) 2010-2011 the authors"
__license__ = "GNU GPLv3 (or more recent equivalent)"


import numpy
from hedge.mesh import TAG_ALL, TAG_NONE

from elastic_wave import main as simulation

simulation(write_output=['vtu', 'receivers'],
           allow_features=['mpi', 'cuda'],
           dim=2,
           order=4,
           stfree_tag=TAG_NONE,
           fix_tag=TAG_ALL,
           op_tag=TAG_NONE,
           flux_type="lf",
           max_steps=None,
           output_dir='output',
           pml=None,#[0, 0, 0, 0.4, 0, 0],
           sources=numpy.array([-12.0,0.0]),
           source_param={'type': 'sinus', 'sigma':1, 'fc':869.0, 'td':1.1e-3, 'begin':0, 'end':2},
           final_time=1,
           quiet_output=True,
           nonlinearity_type='classical',
           mesh_file = 'Meshes/Resonator.msh',
           periodicity = [None, ('minus_y', 'plus_y')],
           material_files = ['Materials/epoxy.dat', 'Materials/steel.dat', 'Materials/steel.dat'],
           vtu_every=20)

Esempio n. 6

File: simulation_3d.py Progetto: pyguerder/hedge-elastodynamics

__copyright__ = "Copyright (C) 2010-2011 the authors"
__license__ = "GNU GPLv3 (or more recent equivalent)"


import numpy
from hedge.mesh import TAG_ALL, TAG_NONE

from elastic_wave import main as simulation

simulation(write_output=['vtu', 'receivers'],
           allow_features=['cuda', 'mpi'],
           dim=3,
           order=4,
           stfree_tag=TAG_NONE,
           fix_tag=TAG_ALL,
           op_tag=TAG_NONE,
           flux_type="lf",
           max_steps=None,
           output_dir='output',
           pml=[0, 0, 0, 500, 0, 0],
           sources=numpy.array([1000.0,0.0,0.0]),
           source_param={'type': 'Ricker', 'sigma':10, 'fc':7.25, 'td':0.16, 'begin':0, 'end':2},
           final_time=12,
           quiet_output=True,
           nonlinearity_type=None,
           mesh_file = 'Meshes/HeterogeneousPeriodicCube.msh',
           periodicity = [None, None, None],
           material_files = ['Materials/aluminium.dat', 'Materials/calcite.dat'],
           vtu_every=40)

Esempio n. 7

File: simulation_3d_Plate.py Progetto: pyguerder/hedge-elastodynamics

__license__ = "GNU GPLv3 (or more recent equivalent)"


import numpy
from hedge.mesh import TAG_ALL, TAG_NONE

from elastic_wave import main as simulation

simulation(write_output=['vtu', 'receivers'],
           allow_features=['cuda','mpi'],
           dim=3,
           order=4,
           stfree_tag=TAG_ALL,
           fix_tag=TAG_NONE,
           op_tag=TAG_NONE,
           flux_type="lf",
           max_steps=None,
           output_dir='output',
           #pml=[0, 0, 0, 500, 0, 0],
           pml=None,
           sources=numpy.array([-0.4,-0.3,0.003]),
           source_param={'type': 'Ricker', 'sigma':5e-3, 'fc':20e3, 'td':58e-6, 'amp':1e-2, 'begin':0, 'end':10e-3},
           final_time=40e-3,
           quiet_output=True,
           nonlinearity_type=None,
           mesh_file = 'Meshes/3DPlate.msh',
           periodicity = [None, None, None],
           material_files = ['Materials/aluminium.dat', 'Materials/aluminium.dat'],
           vtu_every=100)

Esempio n. 8

File: simulation_1d.py Progetto: pyguerder/hedge-elastodynamics

__authors__ = ["Olivier Bou Matar <*****@*****.**>",
               "Pierre-Yves Guerder <*****@*****.**>"]
__copyright__ = "Copyright (C) 2010-2011 the authors"
__license__ = "GNU GPLv3 (or more recent equivalent)"


from hedge.mesh import TAG_ALL, TAG_NONE

from elastic_wave import main as simulation

simulation(write_output=['vtu', 'receivers'],
           allow_features='mpi',
           dim=1,
           order=4,
           stfree_tag=TAG_NONE,
           fix_tag=TAG_ALL,
           op_tag=TAG_NONE,
           flux_type="lf",
           max_steps=None,
           output_dir='output',
           pml=None,
           sources=[0.0],
           source_param={'type': 'Ricker', 'sigma':10, 'fc':7.25, 'td':0.16, 'begin':0, 'end':2},
           final_time=12,
           quiet_output=True,
           nonlinearity_type=None,
           mesh_file = '',
           periodicity = [None],
           material_files = ['Materials/aluminium.dat'],
           vtu_every=20)

Esempio n. 9

File: wave_guide.py Progetto: pyguerder/hedge-elastodynamics

__license__ = "GNU GPLv3 (or more recent equivalent)"


from hedge.mesh import TAG_ALL, TAG_NONE

from elastic_wave import main as simulation

import numpy

simulation(write_output=['vtu', 'receivers'],
           allow_features=['mpi'],
           dim=2,
           order=2,
           stfree_tag=TAG_NONE,
           fix_tag=TAG_ALL,
           op_tag=TAG_NONE,
           flux_type="lf",
           max_steps=1e4,
           output_dir='wave_guide',
           pml=[1, 0, 0, 1, 0, 0],
           sources=numpy.array([-5, 0.0]),
           source_param={'type': 'SineBurst', 'sigma':1, 'fc':1100, 'td':0.0, 'begin':0, 'end':0.4},
           final_time=1,
           quiet_output=True,
           nonlinearity_type=None,
           mesh_file = 'Meshes/wave_guide.msh',
           periodicity = [None, ('minus_y', 'plus_y')],
           material_files = ['Materials/PDMS.dat', 'Materials/polymer.dat'],
           vtu_every=10)

Esempio n. 10

File: simulation_2d.py Progetto: pyguerder/hedge-elastodynamics

__copyright__ = "Copyright (C) 2010-2011 the authors"
__license__ = "GNU GPLv3 (or more recent equivalent)"


from hedge.mesh import TAG_ALL, TAG_NONE

from elastic_wave import main as simulation


simulation(write_output=['vtu','receiver'],
           allow_features=['mpi', 'cuda'],
           dim=2,
           order=5,
           stfree_tag=TAG_NONE,
           fix_tag=TAG_ALL,
           op_tag=TAG_NONE,
           flux_type="lf",
           max_steps=None,
           output_dir='output',
           pml=[400, 400, 0, 400, 400, 0],
           sources=None,
           source_param={'type': 'Ricker', 'sigma':10, 'fc':7.25, 'td':0.16, 'begin':0, 'end':2},
           final_time=12,
           quiet_output=True,
           nonlinearity_type='classical',
           mesh_file = 'Meshes/BiHeterogeneousPeriodicSquare.msh',
           periodicity = [None, None],
           material_files = ['Materials/epoxy.dat', 'Materials/aluminium.dat', 'Materials/steel.dat'],
           vtu_every=100)