Dean Huang
146 lines
3.3 KiB
Python
146 lines
3.3 KiB
Python
# Copyright 2013 DEVSIM LLC
|
|
#
|
|
# SPDX-License-Identifier: Apache-2.0
|
|
|
|
# cap device to test displacement current at contact
|
|
import devsim
|
|
import test_common
|
|
|
|
device = "MyDevice"
|
|
region = "MyRegion"
|
|
|
|
test_common.CreateSimpleMesh(device, region)
|
|
|
|
###
|
|
### Set parameters on the region
|
|
###
|
|
devsim.set_parameter(
|
|
device=device, region=region, name="Permittivity", value=3.9 * 8.85e-14
|
|
)
|
|
|
|
###
|
|
### Create the Potential solution variable
|
|
###
|
|
devsim.node_solution(device=device, region=region, name="Potential")
|
|
|
|
###
|
|
### Creates the Potential@n0 and Potential@n1 edge model
|
|
###
|
|
devsim.edge_from_node_model(device=device, region=region, node_model="Potential")
|
|
|
|
###
|
|
### Electric field on each edge, as well as its derivatives with respect to
|
|
### the potential at each node
|
|
###
|
|
devsim.edge_model(
|
|
device=device,
|
|
region=region,
|
|
name="ElectricField",
|
|
equation="(Potential@n0 - Potential@n1)*EdgeInverseLength",
|
|
)
|
|
|
|
devsim.edge_model(
|
|
device=device,
|
|
region=region,
|
|
name="ElectricField:Potential@n0",
|
|
equation="EdgeInverseLength",
|
|
)
|
|
|
|
devsim.edge_model(
|
|
device=device,
|
|
region=region,
|
|
name="ElectricField:Potential@n1",
|
|
equation="-EdgeInverseLength",
|
|
)
|
|
|
|
###
|
|
### Model the D Field
|
|
###
|
|
devsim.edge_model(
|
|
device=device, region=region, name="DField", equation="Permittivity*ElectricField"
|
|
)
|
|
|
|
devsim.edge_model(
|
|
device=device,
|
|
region=region,
|
|
name="DField:Potential@n0",
|
|
equation="diff(Permittivity*ElectricField, Potential@n0)",
|
|
)
|
|
|
|
devsim.edge_model(
|
|
device=device,
|
|
region=region,
|
|
name="DField:Potential@n1",
|
|
equation="-DField:Potential@n0",
|
|
)
|
|
|
|
###
|
|
### Create the bulk equation
|
|
###
|
|
devsim.equation(
|
|
device=device,
|
|
region=region,
|
|
name="PotentialEquation",
|
|
variable_name="Potential",
|
|
edge_model="DField",
|
|
variable_update="default",
|
|
)
|
|
|
|
# the topbias is a circuit node, and we want to prevent it from being overridden by a parameter
|
|
devsim.set_parameter(device=device, region=region, name="botbias", value=0.0)
|
|
|
|
for name, equation in (
|
|
("topnode_model", "Potential - topbias"),
|
|
("topnode_model:Potential", "1"),
|
|
("topnode_model:topbias", "-1"),
|
|
("botnode_model", "Potential - botbias"),
|
|
("botnode_model:Potential", "1"),
|
|
):
|
|
devsim.node_model(device=device, region=region, name=name, equation=equation)
|
|
|
|
# attached to circuit node
|
|
devsim.contact_equation(
|
|
device=device,
|
|
contact="top",
|
|
name="PotentialEquation",
|
|
node_model="topnode_model",
|
|
edge_charge_model="DField",
|
|
circuit_node="topbias",
|
|
)
|
|
# attached to ground
|
|
devsim.contact_equation(
|
|
device=device,
|
|
contact="bot",
|
|
name="PotentialEquation",
|
|
node_model="botnode_model",
|
|
edge_charge_model="DField",
|
|
)
|
|
|
|
#
|
|
# Voltage source
|
|
#
|
|
devsim.circuit_element(name="V1", n1=1, n2=0, value=1.0, acreal=1.0)
|
|
devsim.circuit_element(name="R1", n1="topbias", n2=1, value=1e3)
|
|
|
|
#
|
|
devsim.solve(type="dc", absolute_error=1.0, relative_error=1e-10, maximum_iterations=30)
|
|
|
|
test_common.print_circuit_solution()
|
|
#
|
|
print(
|
|
devsim.get_contact_charge(
|
|
device=device, contact="top", equation="PotentialEquation"
|
|
)
|
|
)
|
|
print(
|
|
devsim.get_contact_charge(
|
|
device=device, contact="bot", equation="PotentialEquation"
|
|
)
|
|
)
|
|
#
|
|
devsim.solve(type="ac", frequency=1e10)
|
|
test_common.print_ac_circuit_solution()
|
|
|
|
devsim.solve(type="ac", frequency=1e15)
|
|
test_common.print_ac_circuit_solution()
|