Hi Sylvain, thanks again for your detailed answer. I also talked with Arjama; I guess the best solution is to import both Huygens sources (Even if this means loosing the information about the input power) and to simulate the exposure field not with the BC coil, but with the two imported sources. Then, the other simulations (as in the youtube tutorial) can be linked to the field sensor used when simulating the two Huygens-imported sources.

@bryn Thank you very much for your reply. Both the visualization of the electric field in the nerve region and the extraction of the field data as a .mat file have been successfully resolved. I really appreciate your kind support and clear explanations.

This is a simplified example (partly generated using the "To Python" action in the context menu). The example does the following creates an EM simulation as heat source creates a thermal simulation creates a "continue" simulation using the thermal distribution from the first simulation extracts the temperature field from pathlib import Path import numpy as np import s4l_v1.document as document import s4l_v1.materials.database as database import s4l_v1.model as model import s4l_v1.simulation.emfdtd as emfdtd import s4l_v1.simulation.thermal as thermal import s4l_v1.units as units from s4l_v1 import Unit def create_emfdtd_simulation(): simulation = emfdtd.Simulation() entity__lines1 = model.AllEntities()["Edge Source"] entity__sphere1 = model.AllEntities()["Sphere 1"] material_settings = simulation.AddMaterialSettings([entity__sphere1]) simulation.LinkMaterialWithDatabase(material_settings, database["IT'IS LF 5.0"]["Liver"]) edge_source_settings = emfdtd.EdgeSourceSettings() simulation.Add(edge_source_settings, [entity__lines1]) edge_sensor_settings = emfdtd.EdgeSensorSettings() simulation.Add(edge_sensor_settings, [entity__lines1]) automatic_grid_settings = [x for x in simulation.AllSettings if isinstance(x, emfdtd.AutomaticGridSettings) and x.Name == "Automatic"][0] simulation.Add(automatic_grid_settings, [entity__lines1, entity__sphere1]) automatic_voxeler_settings = [x for x in simulation.AllSettings if isinstance(x, emfdtd.AutomaticVoxelerSettings) and x.Name == "Automatic Voxeler Settings"][0] simulation.Add(automatic_voxeler_settings, [entity__lines1, entity__sphere1]) simulation.UpdateAllMaterials() simulation.UpdateGrid() return simulation def create_thermal_transient(em_overall_field): simulation = thermal.TransientSimulation() entity__sphere1 = model.AllEntities()["Sphere 1"] material_settings = simulation.AddMaterialSettings([entity__sphere1]) simulation.LinkMaterialWithDatabase(material_settings, database["IT'IS LF 5.0"]["Liver"]) initial_condition_settings = thermal.InitialConditionSettings() initial_condition_settings.InitialTemperature = 37.0, Unit("C") simulation.Add(initial_condition_settings, [entity__sphere1]) transient_heat_source_settings = simulation.AddHeatSourceSettings([em_overall_field]) automatic_grid_settings = [x for x in simulation.AllSettings if isinstance(x, thermal.AutomaticGridSettings) and x.Name == "Automatic"][0] simulation.Add(automatic_grid_settings, [entity__sphere1]) automatic_voxeler_settings = [x for x in simulation.AllSettings if isinstance(x, thermal.AutomaticVoxelerSettings) and x.Name == "Automatic Voxeler Settings"][0] simulation.Add(automatic_voxeler_settings, [entity__sphere1]) simulation.UpdateAllMaterials() simulation.UpdateGrid() return simulation if __name__ == "__main__": # create dummy model line = model.CreatePolyLine([model.Vec3(0,0,20), model.Vec3(5,0,20)]) line.Name = "Edge Source" sphere = model.CreateSolidSphere(model.Vec3(0,0,0), radius=15) sphere.Name = "Sphere 1" sphere.MaterialName = "Liver" # create heat source em_sim = create_emfdtd_simulation() em_sim.Name = "EM Heat Source" document.AllSimulations.Add( em_sim ) # setup initial thermal simulation th_sim = create_thermal_transient(em_sim.AllComponents["Overall Field"]) th_sim.Name = "T Initial Simulation" document.AllSimulations.Add( th_sim ) # setup simulation using result from previous simulation th_continue = th_sim.CloneWithDiscretization() th_continue.Name = "T Continue Simulation" init_cond_settings = th_continue.GlobalInitialConditionSettings init_cond_settings.InitializationOptions = init_cond_settings.InitializationOptions.enum.ContinueSimulation init_cond_settings.InputSensor = "Overall Field" init_cond_settings.InputSimulation = th_sim.Name document.AllSimulations.Add( th_continue ) # run simulations project_file_path = Path.home() / "t_continue.smash" em_sim.CreateVoxels(project_file_path) em_sim.RunSimulation() th_sim.CreateVoxels(project_file_path) th_sim.RunSimulation() th_continue.CreateVoxels(project_file_path) th_continue.RunSimulation() # extract some results th_continue_extractor = th_continue.Results() th_sensor_extractor = th_continue_extractor["Overall Field"] th_sensor_extractor.Outputs["T(x,y,z,t)"].Update() field = th_sensor_extractor.Outputs["T(x,y,z,t)"].Data print(field)

to get the current density in Python you could use a script like import s4l_v1.document as document try: # add a SimulationExtractor for the simulation call "LF" simulation = document.AllSimulations["LF"] simulation_extractor = simulation.Results() # create an EmSensorExtractor em_sensor_extractor = simulation_extractor["Overall Field"] em_sensor_extractor.FrequencySettings.ExtractedFrequency = u"All" document.AllAlgorithms.Add(em_sensor_extractor) # update the pipeline, make sure current density is extracted em_sensor_extractor.Outputs["J(x,y,z,f0)"].Update() # get data current_density_field = em_sensor_extractor.Outputs["J(x,y,z,f0)"].Data except Exception as exc: import traceback traceback.print_exc()

Dear @lorenero_99, Thank you for contacting us! Please have a read to my recent comment to this post, where I tried to be very detail in explaining how the current flux normalization works and what are the limitations. https://forum.zmt.swiss/topic/733/normalization-for-precise-current-control-via-jupyter/5 For what concerns the impedance, if you use the Ohmic-Current Dominated solver, the Ohmic laws apply and you can extract the resistance (impedance) knowing the applied voltages at the electrodes and the current through them. If you need further explanations, please do not hesitate to contact us immediately! All the best, Antonino

Hi @AntoninoMC, I really appreciate this answer as it helped clear some things up for me. I understand the current extractor much better than I did before and I've realized why I may have been running into some issues. The method of using the analysis as a source also seems to be preferable. I am using the EMLF Electro Quasi-Static model for modelling transcutaneous spinal cord stim, for which precise control of the current is of absolute importance. The multiplier method seems to be the way to go for feeding into the NEURON simulation, as it means I don't have to rerun my EMLF sim. One thing I'm wondering, when I select the analysis/cache as a source, I lose the ability to use the contact impedance model in the NEURON setup. If I use a contact impedance model producer to modulate the resulting EM field from the multiplier, I lose the ability to set pulse parameters, which then seems to prevent NEURON from running a simulation. Do you know if there's a way to use the analysis connection with a multiplier AND a contact impedance model simultaneously for driving a NEURON simulation?

If I understand correctly, you would like to have one simulation with a total current of 2.47mA and another one with 1.98mA. Sim4Life EQS solvers use Dirichlet boundary conditions for the electrodes (fixed voltage) as they solve for electric potential and electrode surfaces are assumed to be equipotential, hence the flux should be able to vary to ensure an equipotential surface. Moreover, when you check “Treat as port” for the LF solvers, you should be able to see in the log that they run 1 simulation per port, by setting this port to 1V and all the rest to 0V. With only 2 electrodes, and assuming you have one anode and one cathode, I would suggest the following: Assign Dirichlet boundary conditions to the electrodes (one setting per electrode, for instance 0.5V and -0.5V). Run the simulation Extract the total current of the simulation (select the "Overall Field" and use the "Current Extractor" from the ribbon). Normalize your fields of interest by the scaling factor: desired current/total current. Visualize the output of your normalized field. I hope this helps!

Hi @Seifeldin_E I assume you’re referring to the second row, where the field is shown on the surface of the spinal cord and peripheral nerves. Here’s the simplest way to achieve this: In the Field Sensor, select Current Density. Add a Masking Filter: Set selection to None (deselect everything). Type Nerve in the search filter and activate all nerve structures. Search for the spinal cord and activate it as well. Add a Surface Viewer — this will extract the surface at the masked regions. How it works: The masking filter replaces all unselected field values with NaN (not-a-number). The surface viewer detects NaN values and extracts the surface surrounding the masked voxels. Alternative method (less robust): Drag a TriangleMesh entity (e.g., Spinal_cord) to the Analysis tab. Select Current Density and the spinal cord mesh (now in Analysis, apply the Model to Grid filter). Add the Interpolation filter. Add a Surface Viewer Be aware: if the field changes abruptly near the surface, this method may interpolate the field on the "wrong" side of the surface. That’s why the masking approach is usually more reliable and easier to use.

The RestoreCamera function has a second argument animate=True. Setting this to False should fix your issue. If you are using s4l_v1.renderer.SaveScreenCapture to save an image of the scene in your script, you might need to give the GUI a chance to refresh during script execution. On Windows, you can do this with win32gui.PumpWaitingMessages() I typically do something like this def refresh_gui() if sys.platform == "win32": import win32gui win32gui.PumpWaitingMessages() for config in all_configurations: # change the model, run a simulation, change camera settings, etc. do_something(config) refresh_gui() s4l_v1.renderer.SaveScreenCapture( width=1024, height=1024, multi_sample=True, transparent_background=True, output_folder="C\temp\screenshots" output_prefix=f"subcase_{config}" )

Hi there, I hope you can help, and that my questions fits within this thread. I have simulated an overall field and I would like to export it to matlab. At first I can live with the UI way i.e. by the import/export menu, but later I'd like to script it too. Now selecting e.g. the E-field and clicking the Imp/Export menu gives me no options but Import. Can anyone show the workflow? And/Or can anyone post a snippet of python code that does the export? Does exporting require a certain license?

Hi, sorry to bother you. Do you happen to know how to calculate the current absorbed by the electrodes and the impedance?

I ran the simulation again. This time it ran completely. Physical results the same but without any errors, i.e. the analysis section could indeed connect ports. So is there a random seed thing and a voxeling issue after all?

Please select the overall field and select the E field, then look under "Field Data Tools". Please "Go to the HTML manual" from "Help" and search for NaN filter. I would recommend using the latest version of Sim4Life v9.0. [image: 1753361215835-bc62f2dd-580a-4b6b-b11d-728454d14416-sim4life_xhh8cgbkwm.png]

Hi, Sorry to bother You. IN ATTACHED VIDEO SHOWS WHAT HAPPENS AFTER CLICKING THE OPTION "CREATE PLOT" AND THE SYSTEM CRASHES. The workstation and memory have been tested, all without errors. Can one of the service engineers help? Kind regards, Thank you, Andrei Andrei Churakov MD,PhD, Associate Professor of the Department of Electronic Engineering and Technology BSUIR E-mail: anchurakov@bsuir.by

This question is quite old. Maybe you found a solution? Otherwise I guess you could use one of these options: you can crop the field, e.g. via the bounding box (or the extent if you have an E-field on a rectilinear grid) you could use the interpolation filter: define your little box grids, and interpolate from the field sensor field to your box grids you could define small field sensors in the simulation setup (not sure it is scalable to hundreds of sensors)

Hi @parsley, Thank you for your feedback and for reporting the issue you experienced. The MATCH tool has already been successfully used/validated by many users and offers several useful features. For example, it allows users to add loading or matching circuits the S-matrix of a simulation and obtain updated scattering parameters without re‑running the full simulation. Have you experienced any problems with this particular workflow? From the information you provided, the issue appears to be related to only the Initial Matching function—the option that generates a matching circuit for a given S-matrix of a simulation at a chosen target frequency. In your case, you reported that it worked with the dipole antenna example but not with your specific antenna model. We would be happy to discuss the details so we can reproduce the problem and resolve it quickly. We will contact you to find out a convenient time to follow up. Thank you again for bringing this to our attention.