Skip to content

Simulations & Solvers

FDTD, Low-Frequency, Neuron, Mode-Matching, Flow, Acoustics, etc...

244 Topics 734 Posts

Subcategories


  • FDTD, Low-Frequency, Neuron, Mode-Matching, Flow, Acoustics, etc...

    11 28
    11 Topics
    28 Posts
    AntoninoMCA
    @KK The titration procedure provides a scaling factor that must be applied to the stimulation source (incident E-field, input current, or voltage applied at electrode pairs) in order to initiate an action potential in a fiber or neuron. Here are some examples depending on the stimulation source: Source: Incident E-field (e.g., induced by TMS) Assume the incident electric field is 1 V/m, and the titration procedure yields a titration factor, tf. The threshold E-field is therefore tf [V/m]. If this 1 V/m field results from a coil current rate of change dI/dt = A [A/s], then the threshold dI/dt required is A × tf [A/s]. Source: Applied voltage Let the applied voltage across an electrode pair be V = V₀. This value is also used as the Dirichlet boundary condition in low-frequency (LF) simulations. If the titration factor is tf, the threshold voltage becomes V₀ × tf [V]. Source: Applied current Let the applied current to an electrode pair be I = I₀, typically derived from LF simulations via a current flux integrator. With a titration factor tf, the threshold current is I₀ × tf [A]. The titration factor tf is a dimensionless number that gains physical meaning when applied to the relevant stimulation quantity (electric field, voltage, or current). @LJ was suggesting the correct course of action. I hope this clarification is helpful!
  • Orientation of antennas affect simulation results

    4
    0 Votes
    4 Posts
    729 Views
    V
    Thanks Ofi. How can I fix this or is this something I need to account for?
  • tutorial not running

    1
    0 Votes
    1 Posts
    253 Views
    No one has replied
  • FDFD solver tutorials

    1
    0 Votes
    1 Posts
    253 Views
    No one has replied
  • Thermal external input file

    1
    0 Votes
    1 Posts
    260 Views
    No one has replied
  • Modeling Ferrite in an FDTD Simulation

    2
    0 Votes
    2 Posts
    419 Views
    ofliO
    Hi @arc Please check "Console" window (open it via menu VIEW | Console if it is not open) to see what kind of error you are getting. The reason for failed simulation might be related to something else. You can easily test with an existing tutorial example. For example, open Dipole Antenna Tutorial, model a rectangular brick next to the dipole, clone one of the existing simulation settings and drag&drop this new brick into the simulation set its material such that rel/ Permeability o 1e4 generate grid, voxel and run to see if it also fails I hope it helps.
  • Python API for Debye dispersive model

    2
    0 Votes
    2 Posts
    470 Views
    C
    Yes, You need to import XMaterials and specify the material model to be linear dispersive, as in the following example (which should be easy to extend), where parameters for one Debye pole are assigned. import XMaterials as xm # Adding a new MaterialSettings material_settings = emfdtd.MaterialSettings() material_settings.ElectricProps.MaterialModel = material_settings.ElectricProps.MaterialModel.enum.LinearDispersive # Specify settings for dispersive poles material_settings.raw.ElectricDispersiveSettings.StartFrequency = 10.e6 material_settings.raw.ElectricDispersiveSettings.EndFrequency = 100000.e6 material_settings.raw.ElectricDispersiveSettings.Conductivity = 0. material_settings.raw.ElectricDispersiveSettings.Permittivity = 11.098 debye_pole_01 = xm.LinearDispersionPole() debye_pole_01.Active = True debye_pole_01.Type = xm.LinearDispersionPole.ePoleType.kDebye debye_pole_01[xm.LinearDispersionPole.ePoleProperty.kDebyeAmplitude] = 1.0 debye_pole_01[xm.LinearDispersionPole.ePoleProperty.kDebyeStaticPermittivity] = 22.67 debye_pole_01[xm.LinearDispersionPole.ePoleProperty.kDebyeInfinityPermittivity] = 11.098 debye_pole_01[xm.LinearDispersionPole.ePoleProperty.kDebyeDamping] = 3.234e-11 # Add dispersive poles to the Linear Electric Dispersion Viewer material_settings.raw.ElectricDispersiveSettings.Poles = [debye_pole_01]
  • Ohmic QS simulation with anisotropic material fails

    2
    1 Votes
    2 Posts
    446 Views
    B
    I suspect that this is because the grid from your LF simulation does not match the grid that you use in your isotropic and anisotropic simulations. The isotropic simulation runs fine because it does not use your cache file. The anisotropic simulation has stored those conductivity values expecting a specific grid. Try right clicking on the grid settings folder in the simulation from which you are creating your anisotropy tensor, select "Copy Grid Configuration", then "Paste Grid Configuration" on the anisotropic simulation. If you want to include features in your anisotropic simulation that you don't don't want to have simulated in your initial LF simulation, you can still include include objects in your initial simulation so that they are considered for the gridding without assigning them material properties, using them as a boundary condition or voxeling them. Just drag the object directly into the grid settings folder. That way your grid can always match.
  • FDTD simulation failure in sensor post-processing stage

    2
    0 Votes
    2 Posts
    391 Views
    ofliO
    Hi @yiyang did it happen with one particular simulation project or does it always happen? Please (if you can) share the project with Sim4Life Application team via sending an email to s4-support@zmt.swiss to check. Thanks
  • 0 Votes
    1 Posts
    328 Views
    No one has replied
  • This topic is deleted!

    1
    0 Votes
    1 Posts
    6 Views
    No one has replied
  • Acoustic simulation

    2
    0 Votes
    2 Posts
    411 Views
    M
    We use standard units in the acoustic solver, which means we're solving the wave equation in pressure which has units of Pa. Note though that the acoustic equation is Linear (careful, this doesn't apply to the nonlinear solver), which means that you can arbitrarily scale the input signal by a scalar and the output (in pressure) will be exactly the same and scaled by that same constant. (You need to be careful when you consider energy related quantitites which are something something pressure squared) Essentially, find the scaling factor between pressure and voltage (assuming linearity), then run your simulation with an arbitritrary amplitude and then scale the pressure output.
  • This topic is deleted!

    1
    0 Votes
    1 Posts
    1 Views
    No one has replied
  • Import Voxels - How to Assign Materials?

    Solved
    6
    0 Votes
    6 Posts
    1k Views
    D
    @dbsim4 Answering my own question: In the Modelling tab (not simulation), each imported label/voxel type can be assigned a material in the Controller window. Auto Assign seems to work decently well if the materials in the file are well named. Probably should be done before dragging the model into the Simulation.
  • 0 Votes
    7 Posts
    1k Views
    J
    @montanaro Thank you! The acoustic reflection simulation is now behaving and giving results which map to what you'd expect. Thank you so much for helping! Jean
  • PdhAddCounter for processor total failed

    1
    0 Votes
    1 Posts
    286 Views
    No one has replied
  • Sensors placement

    1
    0 Votes
    1 Posts
    245 Views
    No one has replied
  • Running the simulation causes the software to hang

    1
    0 Votes
    1 Posts
    263 Views
    No one has replied
  • Unable to see 1g and 10g SAR Values

    2
    0 Votes
    2 Posts
    426 Views
    ofliO
    Dear Shreya, I assume you meant simulated 1g and 10g SAR values. If you are referring to dipole antenna tutorial, it is not possible to extract SAR values since the simulation project does not contain any lossy object - it is just a dipole antenna in free-space. Please check "SAR in a Flat Phantom" tutorial and section 3.1.4 of Sim4Life tutorial document for details related to SAR simulation and extractions.
  • This topic is deleted!

    1
    0 Votes
    1 Posts
    1 Views
    No one has replied
  • Issue with simulation of Temporal Interferece on a head model

    3
    0 Votes
    3 Posts
    671 Views
    Paria M.P
    Dear AntoninoMC, Thank you for your comprehensive solution!!! Following your advice, I assigned distinct priorities to each tissue, and there is no more need for boolean operations (subtraction). This approach proved to be completely successful.