Advanced WFC Concepts

• Foundations of Resonance
• Introduction To Resonance
• LC Circuits
• Q Factor
• Bandwith Ringdown
• Voltage Magnification
• Electric Double Layer
• EDL Introduction
• EDL Capacitance
• Helmholtz Model
• Stern Model
• EDL in WFC
• Electrochemical Impedance
• Impedance Intro
• Randles Circuit
• Cole-Cole Model
• Warburg Impedance
• CPE Elements
• VIC Circuit Theory
• VIC Introduction
• Primary Side
• Secondary Side
• Resonant Charging
• Step Charging
• Choke Design
• Choke Fundamentals
• Core Materials
• Wire Selection
• Bifilar Windings
• Parasitic Effects
• DCR Effects
• Water Fuel Cell Design
• WFC Introduction
• Electrode Geometry
• Water Properties
• Cell Capacitance
• Resonant Matching
• VIC Matrix Calculator
• Calculator Overview
• Component Inputs
• Simulation Tab
• Optimization
• Interpreting Results
• VIC Matrix Calculator Application
• Advanced Topics
• PLL Control
• Harmonic Analysis
• Transformer Coupling
• Energy Efficiency
• Experimental Validation
• Understanding Resonant Action in the Water Fuel Cell
• Appendices
• Complete Formula Reference
• Glossary of Terms
• Wire Gauge Tables
• Core Specifications

Foundations of Resonance

Introduction To Resonance

What is Resonance?

Resonance is a phenomenon that occurs when a system is driven at its natural frequency, causing it to oscillate with maximum amplitude. In electrical circuits, resonance occurs when the inductive and capacitive reactances are equal, creating conditions for energy storage and voltage magnification.

The Physics of Resonance

Every physical system has one or more natural frequencies at which it tends to oscillate. When energy is applied at this frequency, the system absorbs energy most efficiently, leading to large-amplitude oscillations. This principle applies to:

• Mechanical systems: A child on a swing, a vibrating tuning fork
• Acoustic systems: Musical instruments, resonant cavities
• Electrical systems: LC circuits, antennas, oscillators

Electrical Resonance

In electrical circuits containing both inductance (L) and capacitance (C), resonance occurs at a specific frequency where the inductive reactance equals the capacitive reactance:

Resonant Frequency Formula:

f₀ = 1 / (2π√(LC))

Where:

Why Resonance Matters for VIC Circuits

In Stan Meyer's Voltage Intensifier Circuit (VIC), resonance is the key mechanism that enables:

1. Voltage Magnification: At resonance, voltages across reactive components can be many times greater than the input voltage
2. Efficient Energy Transfer: Energy oscillates between the inductor's magnetic field and the capacitor's electric field with minimal loss
3. Impedance Matching: At resonance, the circuit presents a purely resistive impedance to the source

Types of Resonance

Series Resonance

In a series LC circuit, at resonance:

• Impedance is minimum (equals resistance R)
• Current is maximum
• Voltages across L and C can be very high (Q times the source voltage)

Parallel Resonance

In a parallel LC circuit, at resonance:

• Impedance is maximum
• Current from source is minimum
• Circulating current between L and C can be very high

Energy Storage at Resonance

At resonance, energy continuously transfers between the magnetic field of the inductor and the electric field of the capacitor:

Energy in Inductor: E_L = ½LI²

Energy in Capacitor: E_C = ½CV²

At resonance, the total energy remains constant, oscillating between these two forms.

Practical Implications

Understanding resonance is fundamental to designing effective VIC circuits because:

• The primary side (L1-C1) must resonate at the driving frequency
• The secondary side (L2-WFC) should be tuned for optimal energy transfer
• Component values must be carefully calculated to achieve the desired resonant frequency
• The Q factor determines how "sharp" the resonance is and how much voltage magnification occurs

Key Takeaway: Resonance is not just a theoretical concept—it's the working principle behind the VIC's ability to develop high voltages across the water fuel cell while drawing relatively low current from the source.

Next: LC Circuit Fundamentals →

LC Circuits

LC Circuit Fundamentals

An LC circuit consists of an inductor (L) and a capacitor (C) connected together. These circuits form the foundation of resonant systems and are central to understanding how the VIC operates.

Components of an LC Circuit

The Inductor (L)

An inductor stores energy in its magnetic field when current flows through it. Key properties:

• Inductance (L): Measured in Henries (H), represents the inductor's ability to store magnetic energy
• Inductive Reactance: X_L = 2πfL (increases with frequency)
• Current lags voltage by 90° in a pure inductor

The Capacitor (C)

A capacitor stores energy in its electric field between two conductive plates. Key properties:

• Capacitance (C): Measured in Farads (F), represents the capacitor's ability to store electric charge
• Capacitive Reactance: X_C = 1/(2πfC) (decreases with frequency)
• Current leads voltage by 90° in a pure capacitor

Series LC Circuit

Circuit Configuration: L and C connected in series with the source

Total Impedance:

Z = √(R² + (X_L - X_C)²)

At Resonance (X_L = X_C):

Series LC Behavior

Frequency	Condition	Circuit Behavior
f < f₀	XC > XL	Capacitive (current leads voltage)
f = f₀	XC = XL	Resistive (current in phase with voltage)
f > f₀	XL > XC	Inductive (current lags voltage)

Parallel LC Circuit

Circuit Configuration: L and C connected in parallel

At Resonance:

Also called: Tank circuit, because it "tanks" or stores energy

Characteristic Impedance (Z₀)

The characteristic impedance is a fundamental property of any LC circuit:

Z₀ = √(L/C)

This value represents:

Energy Transfer in LC Circuits

In an ideal LC circuit (no resistance), energy oscillates perpetually between the inductor and capacitor:

1. Capacitor fully charged: All energy stored in electric field (E = ½CV²)
2. Current building: Energy transferring to inductor
3. Maximum current: All energy stored in magnetic field (E = ½LI²)
4. Current decreasing: Energy transferring back to capacitor
5. Cycle repeats at the resonant frequency

LC Circuits in the VIC

The VIC uses LC circuits in two critical locations:

Primary Side (L1-C1)

• L1 = Primary choke inductance
• C1 = Tuning capacitor
• Tuned to the driving frequency from the pulse generator
• Develops the initial voltage magnification

Secondary Side (L2-WFC)

• L2 = Secondary choke inductance
• WFC = Water Fuel Cell capacitance
• May be tuned to the same or a harmonic frequency
• Delivers magnified voltage to the water

Design Principle: The relationship between L and C values determines not only the resonant frequency but also the characteristic impedance, which affects how much voltage magnification is achievable.

Practical Considerations

• Component tolerances: Real components have tolerances that affect the actual resonant frequency
• Parasitic elements: Inductors have parasitic capacitance, capacitors have parasitic inductance
• Temperature effects: Component values can drift with temperature
• Losses: Real circuits have resistance that dampens oscillations

Next: Quality Factor (Q) Explained →

Q Factor

Quality Factor (Q) Explained

The Quality Factor, or Q, is one of the most important parameters in resonant circuit design. It quantifies how "sharp" a resonance is and directly determines the voltage magnification achievable in a VIC circuit.

What is Q Factor?

The Q factor is a dimensionless parameter that describes the ratio of energy stored to energy dissipated per cycle in a resonant system. A higher Q means:

• Lower losses relative to stored energy
• Sharper resonance peak
• Higher voltage magnification at resonance
• Narrower bandwidth
• Longer ring-down time when excitation stops

Q Factor Formula

For a series RLC circuit, Q can be calculated several ways:

Primary Definition:

Q = (2π × f₀ × L) / R

Alternative Forms:

Q = X_L / R = (ωL) / R

Q = 1 / (ωCR) = X_C / R

Q = (1/R) × √(L/C) = Z₀ / R

Where:

Physical Meaning of Q

Q can be understood as:

Q = 2π × (Energy Stored / Energy Dissipated per Cycle)

A Q of 100 means the circuit stores 100/(2π) ≈ 16 times more energy than it loses per cycle.

Q Factor and Voltage Magnification

At resonance, the voltage across the inductor (or capacitor) is magnified by the Q factor:

V_L = V_C = Q × V_input

Example: With Q = 50 and V_input = 12V:

V_L = 50 × 12V = 600V across the inductor!

This is why Q factor is so critical in VIC design—it directly determines how much voltage amplification the circuit provides.

Factors Affecting Q

Resistance Sources

Resistance Source	Description	How to Minimize
Wire DCR	DC resistance of the wire	Use larger gauge, shorter length, or copper
Skin Effect	AC resistance increase at high frequency	Use Litz wire or multiple strands
Core Losses	Hysteresis and eddy currents in core	Use appropriate core material for frequency
Capacitor ESR	Equivalent series resistance of capacitor	Use low-ESR capacitors (film, ceramic)
Connection Resistance	Resistance at joints and connections	Use solid connections, avoid corrosion

Wire Material Impact on Q

Different wire materials have vastly different resistivities:

Material	Relative Resistivity	Effect on Q
Copper	1.0× (reference)	Highest Q (best for resonant circuits)
Aluminum	1.6×	Good Q, lighter weight
SS316	~45×	Lower Q, but corrosion resistant
SS430 (Ferritic)	~60×	Much lower Q, magnetic properties
Nichrome	~65×	Very low Q, used for heating elements

Typical Q Values

• Air-core inductors: Q = 50-300 (very low losses)
• Ferrite-core inductors: Q = 20-100 (depends on frequency)
• Iron-powder cores: Q = 50-150
• Practical VIC chokes: Q = 10-50 (with resistance wire, lower)

Q and Bandwidth Relationship

Q is inversely related to bandwidth:

BW = f₀ / Q

Where BW is the -3dB bandwidth (the frequency range where response is within 70.7% of peak).

Example: At f₀ = 10 kHz with Q = 50:

BW = 10,000 / 50 = 200 Hz

Practical Q Measurement

Q can be measured experimentally by:

1. Frequency sweep method: Find f₀ and the -3dB points, then Q = f₀/BW
2. Ring-down method: Count cycles for amplitude to decay to 1/e (37%)
3. LCR meter: Direct measurement at specific frequencies

VIC Design Insight: While higher Q gives more voltage magnification, it also means the circuit is more sensitive to frequency drift and component tolerances. A practical VIC design balances high Q for voltage gain against stability and ease of tuning.

Next: Bandwidth & Ring-Down Decay →

Bandwith Ringdown

Bandwidth & Ring-Down Decay

Understanding bandwidth and ring-down decay is essential for designing VIC circuits that maintain resonance under varying conditions and for predicting how the circuit behaves when the driving signal stops.

Bandwidth Fundamentals

Bandwidth describes the frequency range over which a resonant circuit responds effectively. It's measured as the difference between the upper and lower frequencies where the response drops to 70.7% (-3dB) of the peak value.

Bandwidth Formula:

BW = f₀ / Q

Or equivalently:

BW = R / (2πL)

Where:

Bandwidth and Q Relationship

Q Factor	Bandwidth (at f₀ = 10 kHz)	Frequency Tolerance
Q = 10	1000 Hz	±5% (very forgiving)
Q = 50	200 Hz	±1% (requires tuning)
Q = 100	100 Hz	±0.5% (precise tuning needed)
Q = 200	50 Hz	±0.25% (critical tuning)

Practical Implications of Bandwidth

Narrow Bandwidth (High Q)

• Advantages: Maximum voltage magnification, better selectivity
• Disadvantages: Sensitive to frequency drift, requires precise tuning, may need PLL control

Wide Bandwidth (Low Q)

• Advantages: Easier to tune, more stable, tolerant of component variations
• Disadvantages: Lower voltage magnification, less efficient energy storage

Ring-Down Decay

When the driving signal stops, a resonant circuit doesn't immediately stop oscillating—it "rings down" as stored energy dissipates through resistance. This behavior provides insight into the circuit's Q factor.

Decay Time Constant (τ)

Decay Time Constant:

τ = 2L / R

This is the time for the oscillation amplitude to decay to 1/e (≈37%) of its initial value.

Relationship to Q:

τ = Q / (π × f₀)

Decay Envelope

The amplitude of oscillations during ring-down follows an exponential decay:

A(t) = A₀ × e^-t/τ = A₀ × e^-αt

Where α = R/(2L) is the damping factor.

Damped Oscillation Frequency

During ring-down, the actual oscillation frequency is slightly lower than the natural frequency due to damping:

Damped Frequency:

f_d = √(f₀² - α²/(4π²))

For high-Q circuits (Q > 10), f_d ≈ f₀ (the difference is negligible).

Ring-Down Cycles

A practical measure of how long oscillations persist:

Cycles to 1% Amplitude:

N_1% ≈ Q × 0.733

This is the number of oscillation cycles before amplitude drops to 1% of initial.

Examples:

Ring-Down in VIC Circuits

Understanding ring-down is important for VIC operation because:

Pulsed Operation

• VIC circuits are typically driven by pulsed waveforms
• Between pulses, the circuit rings down
• The ring-down period affects how energy is delivered to the WFC

Step-Charging Considerations

• Each pulse adds energy to the resonant system
• If pulses arrive before ring-down completes, energy accumulates
• This can lead to voltage build-up (step-charging effect)

Measuring Ring-Down

To experimentally determine Q from ring-down:

1. Apply a burst of oscillations at the resonant frequency
2. Stop the driving signal and observe the decay on an oscilloscope
3. Count the number of cycles for amplitude to drop to 37% (1/e)
4. Q ≈ π × (number of cycles to 1/e)

Oscilloscope Tip: Use the "Single" trigger mode to capture the ring-down event. Measure from the point where driving stops to where amplitude reaches ~37% of initial peak.

Summary Table

Parameter	Formula	Depends On
Bandwidth	BW = f₀/Q = R/(2πL)	Resistance, inductance
Decay Time Constant	τ = 2L/R	Inductance, resistance
Damping Factor	α = R/(2L)	Resistance, inductance
Cycles to 1%	N ≈ 0.733 × Q	Q factor only

Design Insight: The VIC Matrix Calculator shows bandwidth and ring-down parameters for your circuit design. Use these to understand how sensitive your circuit will be to frequency variations and how it will behave during pulsed operation.

Next: Voltage Magnification at Resonance →

Voltage Magnification

Voltage Magnification at Resonance

Voltage magnification is the cornerstone of VIC circuit operation. At resonance, the voltage across reactive components (inductors and capacitors) can be many times greater than the input voltage. This is how the VIC develops high voltages across the water fuel cell while drawing modest current from the source.

The Principle of Voltage Magnification

In a series resonant circuit, even though the total impedance is at minimum (just resistance), the individual voltages across L and C can be much larger than the source voltage. This isn't "free energy"—it's the result of energy continuously cycling between the inductor and capacitor.

Key Insight:

At resonance, V_L and V_C are equal in magnitude but opposite in phase. They cancel each other in the circuit loop, but individually each represents a real voltage that can do work.

Voltage Magnification Formula

Q-Based Magnification:

V_output = Q × V_input

Impedance-Based Magnification:

Magnification = Z₀ / R = (1/R) × √(L/C)

Both formulas give the same result since Q = Z₀/R for a series circuit.

Practical Examples

Input Voltage	Q Factor	Output Voltage	Application
12V	10	120V	Low-Q experimental setup
12V	50	600V	Typical VIC circuit
12V	100	1200V	High-Q optimized circuit
24V	50	1200V	Higher input voltage approach

Where the Magnified Voltage Appears

In a Series LC Circuit

• Across the inductor: V_L = Q × V_source (leads current by 90°)
• Across the capacitor: V_C = Q × V_source (lags current by 90°)
• Across resistance: V_R = V_source (in phase with current)

In the VIC Circuit

The water fuel cell acts as the capacitor, so the magnified voltage appears directly across the water:

VIC Voltage Path:

Source → L1 → C1 (series resonance for initial magnification)

Transformed via coupling to → L2 → WFC (secondary resonance)

Result: High voltage across water fuel cell electrodes

Two Approaches to Magnification

Method 1: Maximize Q

Increase Q by reducing resistance:

• Use copper wire instead of resistance wire
• Use larger gauge wire
• Minimize connection resistances
• Use low-ESR capacitors

Method 2: Optimize Z₀/R Ratio

Increase characteristic impedance relative to resistance:

• Increase inductance (more turns, larger core)
• Decrease capacitance (for same resonant frequency, requires more inductance)
• The ratio √(L/C) determines Z₀

Design Trade-off:

For a given resonant frequency f₀ = 1/(2π√LC):

The optimal design balances these factors.

Energy Considerations

Voltage magnification doesn't violate energy conservation:

Power In = Power Dissipated

At steady-state resonance:

The high voltage across L and C represents reactive power—energy that sloshes back and forth but isn't consumed.

Real Power vs. Reactive Power

Type	Symbol	Unit	Description
Real Power	P	Watts (W)	Actually consumed, heats resistors
Reactive Power	Q (or VAR)	Volt-Amperes Reactive	Oscillates, stored in L and C
Apparent Power	S	Volt-Amperes (VA)	Total power flow

Magnification in the VIC Matrix Calculator

The VIC Matrix Calculator displays voltage magnification in several ways:

In Choke Designs

• Q Factor: Calculated from inductance and DCR
• Voltage Magnification: Equals Q for series resonance
• Z₀/R Magnification: Alternative calculation method
• Example Output: Shows actual voltage with 12V input

In Circuit Profiles

• Q_L1C: Q factor of primary side (L1 with C1)
• Q_L2: Q factor of secondary side (L2 with WFC)
• Voltage Magnification: Expected magnification at resonance

Practical Note: Real circuits achieve somewhat less than theoretical magnification due to losses not accounted for in simple models (core losses, radiation, dielectric losses in capacitors, etc.). Expect 70-90% of calculated values in practice.

Safety Warning

⚠️ High Voltage Hazard

Resonant circuits can develop dangerous voltages even from low-voltage sources:

Chapter 1 Complete. Next: The Electric Double Layer (EDL) →

Electric Double Layer

EDL Introduction

What is the Electric Double Layer?

The Electric Double Layer (EDL) is a fundamental electrochemical phenomenon that occurs at the interface between an electrode and an electrolyte solution. Understanding the EDL is crucial for modeling the behavior of water fuel cells in VIC circuits.

The Discovery of the Double Layer

When a metal electrode is immersed in an electrolyte solution, a complex structure spontaneously forms at the interface. This structure, known as the Electric Double Layer, was first described by Hermann von Helmholtz in 1853 and has been refined by many researchers since.

Why Does the Double Layer Form?

Several factors contribute to double layer formation:

1. Charge Separation: The electrode surface may carry an electrical charge (positive or negative)
2. Ion Attraction: Ions of opposite charge in the solution are attracted to the electrode surface
3. Solvent Molecules: Water molecules orient themselves in the electric field near the surface
4. Thermal Motion: The tendency of ions to disperse due to random thermal motion opposes the attraction

Structure of the Double Layer

The EDL consists of several distinct regions:

1. The Electrode Surface

The metal electrode where electronic charge resides.

2. The Inner Helmholtz Plane (IHP)

The plane passing through the centers of specifically adsorbed ions (ions that have lost their solvation shell and are in direct contact with the electrode).

3. The Outer Helmholtz Plane (OHP)

The plane passing through the centers of solvated ions at their closest approach to the electrode.

4. The Diffuse Layer

A region extending into the bulk solution where ion concentration gradually returns to the bulk value.

The Double Layer as a Capacitor

The EDL behaves like a capacitor because:

• Charge is separated across a distance (the Helmholtz layer thickness)
• The layer stores electrical energy in the electric field
• It can be charged and discharged like a conventional capacitor

EDL Capacitance (Simplified Helmholtz Model):

C_dl = ε₀ × ε_r × A / d

Where:

Typical EDL Capacitance Values

Because the separation distance is so small (nanometers), EDL capacitance is remarkably high:

System	Typical Cdl	Notes
Metal in aqueous electrolyte	10-40 µF/cm²	Depends on electrode material and potential
Stainless steel in water	20-30 µF/cm²	Typical for WFC electrodes
Mercury electrode	15-25 µF/cm²	Well-studied reference system

Comparison with Conventional Capacitors

The EDL capacitance is extraordinarily high compared to conventional capacitors:

Example Comparison:

EDL in Water Fuel Cells

In a water fuel cell, the EDL forms at both electrodes:

1. Anode (positive electrode): Attracts negative ions (OH⁻, Cl⁻ if present)
2. Cathode (negative electrode): Attracts positive ions (H⁺, Na⁺ if present)

These two double layers contribute to the total capacitance of the cell and affect how it responds to applied voltages.

Voltage-Dependence of EDL Capacitance

Unlike ideal capacitors, the EDL capacitance varies with applied potential:

• The capacitance reaches a minimum at the potential of zero charge (PZC)
• It increases as the potential deviates from the PZC in either direction
• This non-linear behavior affects VIC circuit operation

Importance for VIC Design

Understanding the EDL is critical because:

• The WFC capacitance determines the resonant frequency with the secondary choke
• The EDL affects how efficiently energy transfers to the water
• The voltage-dependent capacitance can cause resonant frequency shifts
• Proper matching requires accounting for both geometric and EDL capacitance

Key Takeaway: The Electric Double Layer acts as a high-capacitance, nanoscale capacitor at each electrode surface. In a water fuel cell, the total capacitance includes both the geometric (parallel-plate) capacitance of the electrode gap AND the EDL capacitance at each electrode-water interface.

Next: EDL Capacitance in Water →

EDL Capacitance

EDL Capacitance in Water

Calculating the actual capacitance of a water fuel cell requires understanding how the Electric Double Layer contributes to the total capacitance. This page explains how to account for EDL effects in your VIC circuit calculations.

Total WFC Capacitance Model

The total capacitance of a water fuel cell is not simply the geometric parallel-plate capacitance. It includes contributions from multiple components:

Series Combination of Capacitances:

1/C_total = 1/C_geo + 1/C_edl,anode + 1/C_edl,cathode

Where:

Geometric Capacitance

The geometric capacitance depends on electrode geometry and water's dielectric constant:

For Parallel Plate Electrodes:

C_geo = ε₀ × ε_r × A / d

Where ε_r ≈ 80 for water at room temperature

For Concentric Tube Electrodes:

C_geo = (2π × ε₀ × ε_r × L) / ln(r_outer/r_inner)

Where L is the tube length, r is the radius

EDL Capacitance Density

The EDL capacitance is typically specified per unit area:

Electrode Material	Cdl (µF/cm²)	Notes
Stainless Steel 316	20-40	Common WFC electrode
Stainless Steel 304	15-35	Also commonly used
Platinum	25-50	High catalytic activity
Graphite/Carbon	10-20	Lower EDL capacitance
Titanium	30-60	Oxide layer affects value

Calculating Total EDL Capacitance

EDL Capacitance for an Electrode:

C_edl = c_dl × A

Where:

Example Calculation

Given:

Calculate:

Geometric capacitance:

C_geo = (8.854×10⁻¹² × 80 × 0.01) / 0.001 = 7.08 nF

EDL capacitance per electrode:

C_edl = 25 µF/cm² × 100 cm² = 2500 µF = 2.5 mF

Total capacitance:

1/C_total = 1/7.08nF + 1/2.5mF + 1/2.5mF

C_total ≈ 7.08 nF (EDL contribution is negligible when C_edl >> C_geo)

When EDL Matters Most

The EDL capacitance becomes significant when:

Condition	EDL Impact	Reason
Very small electrode gap	Minimal	Cgeo becomes very large
Large electrode gap (>5mm)	Minimal	Cgeo is small, dominates total
Small electrode area	Significant	Cedl becomes comparable to Cgeo
High frequency operation	Significant	EDL may not fully form

Frequency Dependence

The EDL capacitance is not constant with frequency:

• Low frequency (<100 Hz): Full EDL capacitance available
• Medium frequency (100 Hz - 10 kHz): EDL partially developed
• High frequency (>10 kHz): EDL contribution decreases; diffuse layer can't follow

This frequency dependence is modeled using the Cole-Cole relaxation model (covered in Chapter 3).

Effect of Water Purity

The ionic content of water affects both conductivity and EDL behavior:

Water Type	Conductivity	EDL Thickness	Cdl Effect
Deionized	<1 µS/cm	~100 nm	Lower Cdl
Distilled	1-10 µS/cm	~30 nm	Moderate Cdl
Tap water	200-800 µS/cm	~1 nm	Higher Cdl
With electrolyte (NaOH, KOH)	>1000 µS/cm	<1 nm	Highest Cdl

In the VIC Matrix Calculator

The VIC Matrix Calculator's Water Profile settings account for EDL effects:

• Electrode material: Determines specific C_dl
• Water conductivity: Affects EDL thickness and capacitance
• Temperature: Influences dielectric constant and ion mobility
• EDL thickness parameter: Allows fine-tuning based on measurements

Practical Tip: For most VIC calculations using typical electrode gaps (1-3mm), the geometric capacitance dominates. However, for very close electrode spacing or when precise tuning is needed, including EDL effects can improve accuracy.

Next: The Helmholtz Model →

Helmholtz Model

The Helmholtz Model

The Helmholtz model is the simplest description of the Electric Double Layer. While it has limitations, it provides an intuitive understanding of how charge separation occurs at electrode surfaces and remains useful for quick calculations.

Historical Background

In 1853, Hermann von Helmholtz proposed the first model of the electrode-electrolyte interface. He imagined ions arranging themselves in a single, compact layer at the electrode surface—like opposite plates of a capacitor.

The Helmholtz Picture

Key Assumptions:

Visual Representation

    ELECTRODE           SOLUTION

    ┃ + + + + ┃          ⊖   ⊖   ⊖   ⊖
    ┃ + + + + ┃     →    ⊖   ⊖   ⊖   ⊖    (bulk solution)
    ┃ + + + + ┃          ⊖   ⊖   ⊖   ⊖
    ┃ + + + + ┃          ⊖   ⊖   ⊖   ⊖

    |←── d ──→|

    Helmholtz   Inner layer
    layer       of counter-ions

Mathematical Description

The Helmholtz model treats the interface as a simple parallel-plate capacitor:

Helmholtz Capacitance:

C_H = ε₀ε_rA / d

Where:

Note on Dielectric Constant

The relative permittivity (ε_r) in the Helmholtz layer is much lower than bulk water:

Region	εr	Reason
Bulk water	~80	Free rotation of water dipoles
Helmholtz layer	~6-10	Water molecules strongly oriented by electric field
Ice	~3	Fixed molecular orientation

Calculating Helmholtz Capacitance

Example Calculation:

For a typical metal electrode in aqueous solution:

C_H/A = ε₀ε_r/d = (8.854 × 10⁻¹² × 6) / (3 × 10⁻¹⁰)

C_H/A = 0.177 F/m² = 17.7 µF/cm²

Potential Distribution

In the Helmholtz model, the potential drops linearly from the electrode to the ion plane:

φ(x) = φ_electrode - (φ_electrode - φ_solution) × (x/d)

Where x is the distance from the electrode (0 ≤ x ≤ d)

Electric Field in the Layer

The electric field is constant throughout the Helmholtz layer:

E = (φ_electrode - φ_solution) / d = ΔV / d

Example: With ΔV = 1V and d = 0.3 nm:

E = 1V / (3 × 10⁻¹⁰ m) = 3.3 × 10⁹ V/m = 3.3 GV/m

This is an enormous electric field! Such high fields strongly polarize water molecules.

Limitations of the Helmholtz Model

While useful for intuition, the Helmholtz model fails to explain several observations:

Observation	Helmholtz Prediction	Reality
Capacitance vs. concentration	No dependence	Capacitance increases with ion concentration
Capacitance vs. potential	Constant	Varies with applied potential
Temperature dependence	Only through εr	More complex behavior

When to Use the Helmholtz Model

Despite its limitations, the Helmholtz model is appropriate when:

• Quick, order-of-magnitude estimates are needed
• The electrolyte concentration is high (>0.1 M)
• Only the compact layer capacitance is of interest
• Building intuition about EDL behavior

Extension to the VIC Context

In VIC applications, the Helmholtz model helps understand:

1. Maximum possible EDL capacitance: Sets an upper bound on what the interface can contribute
2. Field strength at the electrode: Related to the electrochemical driving force
3. Effect of surface area: Larger electrodes = more capacitance

Key Insight: The Helmholtz model shows that double layer capacitance is fundamentally limited by the minimum approach distance of ions to the electrode (d ≈ 0.3 nm). This explains why EDL capacitance is so high—the "plates" are incredibly close together!

Next: The Stern Layer Model →

Stern Model

The Stern Layer Model

The Stern model combines the best features of the Helmholtz and Gouy-Chapman models, providing a more realistic description of the Electric Double Layer that accounts for both the compact ion layer and the diffuse layer extending into solution.

Why a Better Model Was Needed

The Helmholtz model (single compact layer) and the Gouy-Chapman model (purely diffuse layer) both had shortcomings:

Model	Strength	Weakness
Helmholtz	Predicts correct order of magnitude for C	No concentration or potential dependence
Gouy-Chapman	Explains concentration dependence	Predicts infinite C at high potentials

Otto Stern (1924) resolved these issues by combining both approaches.

The Stern Model Structure

The model divides the double layer into two regions:

1. Stern Layer (Compact Layer)

2. Diffuse Layer (Gouy-Chapman Layer)

Visual Representation

    ELECTRODE    STERN LAYER      DIFFUSE LAYER        BULK

    ┃ + + + ┃   H₂O  ⊖  H₂O     ⊖    ⊖        ⊕  ⊖
    ┃ + + + ┃   H₂O  ⊖  H₂O        ⊖     ⊕      ⊖
    ┃ + + + ┃   H₂O  ⊖  H₂O     ⊖       ⊖   ⊕
    ┃ + + + ┃   H₂O  ⊖  H₂O           ⊖       ⊕  ⊖

              |← IHP  OHP →|←──── λD ────→|

              |←── Stern ──→|←── Diffuse ─→|

    IHP = Inner Helmholtz Plane
    OHP = Outer Helmholtz Plane
    λD = Debye Length

Potential Distribution

The potential varies differently in each region:

In the Stern Layer (0 ≤ x ≤ d):

φ(x) = φ_M - (φ_M - φ_d) × (x/d)

Linear drop from metal potential (φ_M) to diffuse layer potential (φ_d)

In the Diffuse Layer (x > d):

φ(x) = φ_d × exp(-(x-d)/λ_D)

Exponential decay with characteristic length λ_D (Debye length)

The Debye Length

The Debye length (λ_D) characterizes how far the diffuse layer extends:

λ_D = √(ε₀ε_rk_BT / (2n₀e²z²))

For a 1:1 electrolyte in water at 25°C:

λ_D ≈ 0.304 / √c (nm)

Where c is the molar concentration (M).

Debye Length Examples

Concentration	Debye Length	Context
10⁻⁷ M (pure water)	~960 nm	Deionized water
10⁻⁴ M	~30 nm	Distilled water
10⁻³ M	~10 nm	Tap water
10⁻² M	~3 nm	Dilute electrolyte
0.1 M	~1 nm	Concentrated electrolyte

Total Capacitance in Stern Model

The Stern and diffuse layer capacitances are in series:

1/C_total = 1/C_Stern + 1/C_diffuse

Stern Layer Capacitance:

C_Stern = ε₀ε₁A / d

Diffuse Layer Capacitance:

C_diffuse = (ε₀ε_rA / λ_D) × cosh(zeφ_d/2k_BT)

Concentration Effects on Capacitance

The Stern model correctly predicts:

• Low concentration: Diffuse layer is thick (large λ_D), C_diffuse is small, limits total capacitance
• High concentration: Diffuse layer collapses, C_diffuse → ∞, C_total → C_Stern

Practical Implication: In concentrated electrolytes (like tap water with dissolved minerals), the total EDL capacitance approaches the Helmholtz (Stern layer) value. In very pure water, the diffuse layer contribution becomes important.

Temperature Dependence

Temperature affects the Stern model through:

1. Debye length: λ_D ∝ √T (diffuse layer thickens at higher T)
2. Dielectric constant: ε_r decreases with T
3. Thermal voltage: k_BT/e ≈ 26 mV at 25°C

Application to Water Fuel Cells

For VIC circuit design, the Stern model helps predict:

Parameter	Effect on EDL	VIC Design Impact
Adding electrolyte	Compresses diffuse layer	Increases WFC capacitance
Using pure water	Extended diffuse layer	Lower WFC capacitance
Heating water	Thicker diffuse layer	Slightly lower capacitance
Increasing voltage	Higher diffuse layer C	Capacitance increases with V

Key Takeaway: The Stern model shows that EDL capacitance depends on electrolyte concentration. For pure water (low ionic strength), the diffuse layer extends far into solution and reduces total capacitance. Adding small amounts of electrolyte (even tap water impurities) collapses this diffuse layer and increases capacitance toward the Helmholtz limit.

Next: EDL Effects in Water Fuel Cells →

EDL in WFC

EDL Effects in Water Fuel Cells

This page integrates everything we've learned about the Electric Double Layer and applies it specifically to water fuel cell design in VIC circuits. Understanding these effects is crucial for accurate circuit modeling and optimization.

The Complete WFC Electrical Model

A water fuel cell is not a simple capacitor. Its complete electrical model includes:

    ┌────────────────────────────────────────────┐
    │                                            │
    │   ┌─────┐   ┌─────┐   ┌─────┐   ┌─────┐   │
  ──┤   │C_dl1│   │R_ct1│   │R_sol│   │C_dl2│   ├──
    │   │     │   │     │   │     │   │     │   │
    │   └──┬──┘   └──┬──┘   │     │   └──┬──┘   │
    │      │         │      │     │      │      │
    │      └────┬────┘      │     │      └──────┤
    │           │           │     │             │
    │       ┌───┴───┐       │     │      ┌─────┐│
    │       │  W₁   │       │     │      │C_geo││
    │       └───────┘       │     │      └─────┘│
    │                       │     │             │
    │      Anode EDL        │     │  Cathode EDL│
    └────────────────────────────────────────────┘
    

Components:

Frequency-Dependent Behavior

The WFC impedance changes dramatically with frequency:

Frequency Range	Dominant Element	WFC Behavior
Very low (<1 Hz)	Warburg diffusion	Z ~ 1/√f, 45° phase
Low (1-100 Hz)	Charge transfer Rct	Resistive behavior
Medium (100 Hz - 10 kHz)	EDL capacitance Cdl	Capacitive, EDL dominant
High (10 kHz - 1 MHz)	Solution R + geometric C	RC network behavior
Very high (>1 MHz)	Geometric Cgeo	Pure capacitance

EDL Time Constant

The EDL has a characteristic response time:

τ_EDL = R_sol × C_dl

The EDL fully forms in approximately 5×τ_EDL.

Example:

Implication: At frequencies above 1/(2πτ) ≈ 160 Hz, the EDL cannot fully form and its effective capacitance decreases.

Effective WFC Capacitance

At VIC operating frequencies (typically 1-50 kHz), the effective WFC capacitance is:

Simplified Model:

1/C_eff = 1/C_geo + 1/C_dl,eff

Where C_dl,eff is the frequency-reduced EDL capacitance.

Typical VIC Frequency Range:

Non-Linear Capacitance Effects

The EDL capacitance depends on applied voltage:

• Low voltage (<100 mV): Capacitance relatively constant
• Medium voltage (100 mV - 1V): Capacitance increases with voltage
• High voltage (>1V): Electrochemical reactions begin, behavior becomes complex

VIC Implication:

As voltage across the WFC increases during resonant charging, the capacitance changes. This can cause:

Temperature Effects in WFC

Parameter	Temperature Effect	Typical Change
Water εr	Decreases with T	-0.4% per °C
Solution conductivity	Increases with T	+2% per °C
EDL thickness	Increases with T	+0.2% per °C
Reaction rate	Increases with T	~Doubles per 10°C

Practical WFC Design Considerations

Electrode Material Selection

• 316 Stainless Steel: Good corrosion resistance, moderate C_dl
• 304 Stainless Steel: Lower cost, slightly lower performance
• Titanium: Excellent stability, oxide layer affects EDL
• Platinized electrodes: Highest activity, highest C_dl

Electrode Spacing

Trade-offs:

Water Treatment

• Distilled water: Low conductivity, thick diffuse layer, lower total C
• Tap water: Higher conductivity, thinner diffuse layer, higher C
• With electrolyte: Highest conductivity, Helmholtz-dominated C

Measuring WFC Capacitance

To accurately characterize your WFC:

1. Use an LCR meter: Measure at multiple frequencies (100 Hz, 1 kHz, 10 kHz)
2. Perform EIS: Electrochemical Impedance Spectroscopy gives complete picture
3. Measure at operating conditions: Temperature and voltage matter
4. Account for cables: Long leads add inductance and capacitance

Integration with VIC Matrix Calculator

The VIC Matrix Calculator accounts for EDL effects through:

• Water Profile settings: Conductivity, temperature, electrode material
• EDL capacitance model: Calculates C_dl based on electrode area
• Frequency correction: Adjusts effective capacitance for operating frequency
• Cole-Cole parameters: Models frequency dispersion (see Chapter 3)

Design Recommendation: For initial VIC designs, use the geometric capacitance as the primary estimate. Include EDL effects when fine-tuning or when using very close electrode spacing. The Cole-Cole model (next chapter) provides more accurate frequency-dependent behavior.

Chapter 2 Complete. Next: Electrochemical Impedance →

Electrochemical Impedance

Impedance Intro

Introduction to Electrochemical Impedance

Electrochemical Impedance Spectroscopy (EIS) is a powerful technique for characterizing the electrical behavior of electrochemical systems like water fuel cells. Understanding impedance helps us model and predict how the WFC behaves across different frequencies.

What is Impedance?

Impedance (Z) is the AC equivalent of resistance. While resistance applies only to DC circuits, impedance describes how a circuit element opposes current flow at any frequency, including the phase relationship between voltage and current.

Impedance Definition:

Z = V(t) / I(t) = |Z| × e^jθ = Z' + jZ''

Where:

Impedance of Basic Elements

Element	Impedance	Phase	Frequency Dependence
Resistor (R)	Z = R	0°	None
Capacitor (C)	Z = 1/(jωC)	-90°	|Z| decreases with f
Inductor (L)	Z = jωL	+90°	|Z| increases with f

Why Use Impedance for WFC Analysis?

Impedance spectroscopy reveals information that simple DC measurements cannot:

1. Separating processes: Different phenomena occur at different frequencies
2. Non-destructive: Small AC signals don't significantly perturb the system
3. Complete characterization: Maps all electrical behavior across frequency
4. Model fitting: Allows extraction of equivalent circuit parameters

Electrochemical Impedance Spectroscopy (EIS)

EIS measures impedance across a range of frequencies to create a complete picture:

Typical EIS Procedure:

Nyquist Plot

The Nyquist plot shows the imaginary part (-Z'') vs. the real part (Z') of impedance:

    -Z'' (Ohms)
        ↑
    500 │          ○ ○
        │       ○       ○
    400 │     ○           ○
        │    ○              ○  (Semicircle = RC parallel)
    300 │   ○                ○
        │  ○                   ○
    200 │ ○                      ○
        │○                         ○
    100 │                            ○ ○ ○ ○
        │                                    ↘ (Warburg tail)
      0 └─────────────────────────────────────────→ Z' (Ohms)
          0   200   400   600   800   1000   1200

        High freq                          Low freq
        ←─────────────────────────────────────────→

Reading a Nyquist Plot:

• High frequency intercept: Solution resistance (R_s)
• Semicircle diameter: Charge transfer resistance (R_ct)
• Semicircle peak frequency: Related to R_ct × C_dl
• 45° line at low frequency: Warburg diffusion impedance

Bode Plot

The Bode plot shows magnitude and phase vs. frequency on logarithmic scales:

Bode Magnitude Plot:

|Z| (log scale) vs. frequency (log scale)

Bode Phase Plot:

Phase angle θ vs. frequency (log scale)

Frequency Ranges and Processes

Different electrochemical processes dominate at different frequencies:

Frequency	Process	Circuit Element
> 100 kHz	Bulk solution, cables	Rs, parasitic L
1 kHz - 100 kHz	Double layer charging	Cdl
1 Hz - 1 kHz	Charge transfer kinetics	Rct
< 1 Hz	Mass transport (diffusion)	ZW (Warburg)

Why This Matters for VIC

Understanding EIS helps VIC design in several ways:

• Accurate modeling: Know the true WFC impedance at your operating frequency
• Frequency selection: Choose operating frequencies that optimize energy transfer
• Tuning: Understand why resonance may shift during operation
• Diagnostics: Identify problems from impedance changes

Practical EIS for WFC Characterization

Equipment Needed:

Alternative for Hobbyists:

An audio frequency generator + oscilloscope can characterize WFC in the 20 Hz - 20 kHz range relevant to most VIC circuits.

Key Takeaway: Electrochemical impedance reveals that a WFC is far more complex than a simple capacitor. Its impedance varies with frequency, voltage, temperature, and time. The equivalent circuit models in the following pages help capture this complexity for VIC design.

Next: The Randles Equivalent Circuit →

Randles Circuit

The Randles Equivalent Circuit

The Randles circuit is the most widely used equivalent circuit model for electrochemical interfaces. It captures the essential elements of an electrode-electrolyte system and serves as the foundation for more complex models used in WFC analysis.

The Classic Randles Circuit

Proposed by John Randles in 1947, this circuit combines resistive, capacitive, and diffusion elements:

         Rs                 Rct
    ────┬────┬────────────┬────┬────
        │    │            │    │
        │    │            │    │
        │  ──┴──        ──┴──  │
        │  │   │        │    │ │
        │  │Cdl│        │ Zw │ │
        │  │   │        │    │ │
        │  ──┬──        ──┬──  │
        │    │            │    │
        └────┴────────────┴────┘

    Rs  = Solution resistance
    Cdl = Double layer capacitance
    Rct = Charge transfer resistance
    Zw  = Warburg diffusion impedance

Component Meanings

Element	Physical Origin	Typical Value (WFC)
Rs	Ionic resistance of electrolyte solution between electrodes	10 Ω - 10 kΩ (depends on conductivity)
Cdl	Electric double layer capacitance at electrode surface	µF to mF range (depends on area)
Rct	Resistance to electron transfer at electrode (reaction kinetics)	1 Ω - 1 MΩ (depends on overpotential)
ZW	Impedance due to diffusion of reactants/products	Frequency-dependent (see Warburg page)

Total Impedance

The total impedance of the Randles circuit is:

Z_total = R_s + [Z_Cdl || (R_ct + Z_W)]

Expanding:

Z_total = R_s + [(R_ct + Z_W)] / [1 + jωC_dl(R_ct + Z_W)]

Frequency Response

The Randles circuit produces a characteristic Nyquist plot:

    -Z''
      ↑
      │           ○ ○ ○
      │        ○         ○
      │      ○             ○           ← Semicircle from Rct||Cdl
      │     ○               ○
      │    ○                 ○
      │   ○                   ○  ○
      │                            ○ ○
      │                                 ○ ○  ← Warburg 45° line
      │                                     ○ ○
      └──────────────────────────────────────────→ Z'
         ↑                    ↑              ↑
         Rs              Rs + Rct      Low freq limit
    (high freq)     (semicircle end)

Time Constants in the Randles Circuit

Double Layer Time Constant:

τ_dl = R_s × C_dl

Determines how quickly the double layer charges through the solution resistance.

Charge Transfer Time Constant:

τ_ct = R_ct × C_dl

Determines the peak frequency of the semicircle: f_peak = 1/(2πτ_ct)

Simplified Cases

Case 1: Fast Kinetics (R_ct → 0)

When the electrochemical reaction is very fast:

• Semicircle disappears
• Only Warburg tail remains at low frequency
• The system is "diffusion-controlled"

Case 2: Slow Kinetics (R_ct → large)

When the electrochemical reaction is slow:

• Large semicircle dominates
• Warburg region may not be visible
• The system is "kinetically-controlled"

Case 3: No Faradaic Reaction (R_ct → ∞)

When no electrochemical reaction occurs (blocking electrode):

• No semicircle
• Purely capacitive behavior at low frequency
• Nyquist plot is a vertical line

Randles Circuit for WFC

In a water fuel cell, the Randles elements have specific meanings:

Element	WFC Interpretation	Effect on VIC
Rs	Water conductivity, electrode gap	Adds to total circuit resistance, reduces Q
Cdl	EDL at each electrode	Part of total WFC capacitance
Rct	Activation barrier for water splitting	Limits DC current, less relevant at high freq
ZW	Diffusion of H₂/O₂ gases, ions	Important at low frequencies only

Extended Randles Circuit

For more accurate WFC modeling, the Randles circuit can be extended:

                   ┌─────────────────────────┐
    Rs             │   Cathode              │
  ──┬──┬──────────┬┴─────────────────────────┴┬──
    │  │          │                           │
    │ Cgeo        │  Rct,c         Rct,a      │
    │  │        ──┴──            ──┴──        │
    │  │        │    │          │    │        │
    │  │        │Cdl,c│         │Cdl,a│       │
    │  │        │    │          │    │        │
    └──┴────────┬────┬──────────┬────┬────────┘
                │    │          │    │
                │ Zw,c│         │ Zw,a│
                └────┘          └────┘

                   Anode

This model includes separate elements for anode and cathode interfaces plus the geometric capacitance.

Parameter Extraction

From an experimental EIS measurement, Randles parameters can be extracted:

1. R_s: High-frequency real-axis intercept
2. R_ct: Diameter of the semicircle
3. C_dl: From peak frequency: C = 1/(2πf_peakR_ct)
4. Warburg coefficient: From slope of the 45° line

Software Tools: Programs like ZView, EC-Lab, and Nova can automatically fit Randles parameters to EIS data. Open-source options include impedance.py (Python) and EIS Spectrum Analyzer.

VIC Design Application: The Randles circuit shows that at VIC operating frequencies (1-50 kHz), the WFC behaves primarily as C_dl in series with R_s. The charge transfer resistance and Warburg impedance become important only at lower frequencies where actual water splitting occurs.

Next: Cole-Cole Relaxation Model →

Cole-Cole Model

Cole-Cole Relaxation Model

The Cole-Cole model describes how the dielectric properties of materials change with frequency. In WFC applications, it provides a more accurate model of capacitance dispersion than the simple Randles circuit, especially for systems with distributed time constants.

Origin of the Cole-Cole Model

Kenneth and Robert Cole (1941) observed that many dielectric materials don't follow simple Debye relaxation. Instead, the relaxation is "stretched" across a broader frequency range. The Cole-Cole model quantifies this behavior with a single additional parameter.

The Cole-Cole Equation

Complex Permittivity:

ε*(ω) = ε_∞ + (ε_s - ε_∞) / [1 + (jωτ)^(1-α)]

Where:

The α Parameter

The Cole-Cole parameter α describes the "spread" of relaxation times:

α Value	Behavior	Physical Meaning
α = 0	Simple Debye relaxation	Single relaxation time, ideal system
α = 0.1-0.3	Slight distribution	Minor surface heterogeneity
α = 0.3-0.5	Moderate distribution	Typical for WFC electrodes
α = 0.5-0.7	Broad distribution	Rough or porous electrodes
α → 1	Extreme distribution	Highly disordered system

Cole-Cole Plot

Plotting -ε'' vs. ε' creates the characteristic Cole-Cole diagram:

    -ε''
      ↑
      │
      │        Debye (α=0)             Cole-Cole (α>0)
      │          ○ ○ ○                    ○ ○ ○
      │       ○       ○                ○         ○
      │      ○         ○              ○           ○
      │     ○           ○            ○             ○
      │    ○             ○          ○               ○
      │                           ○                   ○
      │                         ○                       ○
      └────────────────────────────────────────────────────→ ε'
          ε∞        ε                ε∞        ε
                    ▲ s                        ▲ s
              Perfect                   Depressed
              semicircle                semicircle

     Center on           Center below
     real axis           real axis

The Cole-Cole model produces a depressed semicircle, with the center located below the real axis.

Depression Angle

The depression angle θ relates to α:

θ = α × (π/2) radians = α × 90°

Example: α = 0.3 gives θ = 27° depression

Physical Origins of Distribution

Why do WFC systems show Cole-Cole behavior?

• Surface roughness: Different local environments at electrode surface
• Porous electrodes: Distribution of pore sizes and depths
• Oxide layers: Non-uniform thickness or composition
• Grain boundaries: In polycrystalline electrodes
• Adsorbed species: Non-uniform coverage of adsorbed ions

Impedance Form of Cole-Cole

For circuit modeling, the Cole-Cole element is expressed as impedance:

Z_CC = R / [1 + (jωτ)^(1-α)]

This can be represented as a resistor in parallel with a Constant Phase Element (CPE).

Cole-Cole in the VIC Matrix Calculator

The VIC Matrix Calculator uses the Cole-Cole model for WFC characterization:

Cole-Cole Parameters in the App:

Frequency-Dependent Capacitance

The Cole-Cole model predicts how capacitance varies with frequency:

Effective Capacitance:

C_eff(ω) = C₀ × [1 + (ωτ)^2(1-α)]^-1/2

At low frequency: C_eff → C₀ (full capacitance)

At high frequency: C_eff → C_∞ < C₀ (reduced capacitance)

Practical Example

WFC with Cole-Cole Parameters:

Effective Capacitance at Different Frequencies:

VIC Design Implications

The Cole-Cole model affects VIC design in several ways:

1. Resonant frequency shift: As frequency changes, C_eff changes, shifting resonance
2. Broader resonance: The distribution of time constants broadens the frequency response
3. Q factor reduction: Losses associated with the relaxation reduce circuit Q
4. Frequency selection: Operating below the characteristic frequency maximizes capacitance

Practical Recommendation: For VIC circuits, choose an operating frequency below the Cole-Cole characteristic frequency (f_c = 1/2πτ) to maximize effective WFC capacitance. The VIC Matrix Calculator can help determine optimal operating frequency based on your WFC's Cole-Cole parameters.

Next: Warburg Diffusion Impedance →

Warburg Impedance

Warburg Diffusion Impedance

The Warburg impedance describes mass transport limitations in electrochemical systems. When reactions are fast but reactants or products can't diffuse quickly enough, the Warburg impedance becomes the dominant factor. Understanding this helps predict WFC behavior at low frequencies.

What is Diffusion?

Diffusion is the spontaneous movement of particles from regions of high concentration to low concentration. In electrochemical cells:

• Reactants must diffuse to the electrode surface
• Products must diffuse away from the electrode
• This mass transport takes time and creates a frequency-dependent impedance

The Warburg Element

Semi-Infinite Warburg Impedance:

Z_W = σ/√ω × (1 - j) = σ/√ω - jσ/√ω

Where:

Magnitude and Phase:

|Z_W| = σ√2/√ω (decreases with frequency)

θ = -45° (constant phase)

Warburg Coefficient

The Warburg coefficient depends on the diffusing species:

σ = (RT)/(n²F²A√2) × [1/(D_O^½C_O) + 1/(D_R^½C_R)]

Where:

Nyquist Plot Appearance

    -Z''
      ↑
      │
      │                              Warburg: 45° line
      │                                    ↗
      │                               ↗
      │      Kinetic                ↗
      │      semicircle           ↗
      │        ○ ○ ○          ↗
      │     ○       ○      ↗
      │    ○         ○   ↗
      │   ○           ○↗
      │  ○             ○
      │ ○               ○
      └──────────────────────────────────→ Z'
         Rs          Rs+Rct
                          (transition to diffusion)

     High ←───────── Frequency ──────────→ Low

Types of Warburg Impedance

1. Semi-Infinite Warburg (W)

The classic form, assumes infinite diffusion layer:

• Appears as 45° line on Nyquist plot
• Valid when diffusion layer << electrode separation
• Most common model for thick electrolyte layers

2. Finite-Length Warburg (Wo)

For thin electrolyte layers or porous electrodes:

Z_o = (σ/√ω) × tanh(√(jωτ_D)) / √(jωτ_D)

Where τ_D = L²/D (diffusion time across layer of thickness L)

3. Short Warburg (Ws)

For convection-limited systems:

Z_s = (σ/√ω) × coth(√(jωτ_D)) / √(jωτ_D)

Frequency Dependence

Frequency	|ZW| Behavior	Physical Meaning
Very low	Large	Plenty of time for diffusion to affect response
Medium	Moderate	Partial diffusion limitation
High	Small	Not enough time for concentration gradients

Warburg in Water Fuel Cells

In a WFC, Warburg impedance arises from:

• H₂ diffusion: Hydrogen gas bubbles and dissolved H₂
• O₂ diffusion: Oxygen gas bubbles and dissolved O₂
• Ion migration: H⁺, OH⁻, and electrolyte ions
• Water replenishment: At high current densities

Typical Values for WFC

Parameter	Typical Range	Notes
Warburg coefficient (σ)	1-100 Ω·s-1/2	Higher in pure water
Characteristic frequency	0.01-10 Hz	Depends on diffusion length
Diffusion length	10-1000 µm	Sets electrode spacing limit

Relevance to VIC Operation

Good News for VIC:

At typical VIC operating frequencies (1-50 kHz), the Warburg impedance is negligibly small because:

When Warburg Matters:

Practical Implications

1. Frequency selection: High-frequency operation minimizes diffusion effects
2. Bubble management: Gas bubbles increase Warburg impedance
3. Electrode design: Porous electrodes have complex diffusion paths
4. Stirring/flow: Can reduce diffusion limitations

Measuring Warburg Parameters

To characterize the Warburg element in your WFC:

1. Perform EIS down to very low frequencies (0.01 Hz)
2. Look for the 45° line region in Nyquist plot
3. Measure the slope to determine σ
4. Note the frequency where Warburg transitions to capacitive/resistive

Key Takeaway: The Warburg impedance is important for understanding electrochemical kinetics but becomes negligible at VIC operating frequencies. Focus on the double layer capacitance and solution resistance for high-frequency VIC design. However, be aware that low-frequency or DC operations will encounter significant diffusion effects.

Next: Constant Phase Elements (CPE) →

CPE Elements

Constant Phase Elements (CPE)

The Constant Phase Element (CPE) is a generalized circuit element that better represents real capacitor behavior in electrochemical systems. It accounts for the non-ideal response of electrode surfaces and is essential for accurate WFC modeling.

Why Ideal Capacitors Don't Work

Real electrochemical interfaces rarely behave as ideal capacitors. EIS measurements typically show:

• Depressed semicircles (not perfect)
• Phase angles between -90° and 0° (not exactly -90°)
• Frequency-dependent capacitance

The CPE was introduced to model this non-ideal behavior with a single additional parameter.

CPE Definition

CPE Impedance:

Z_CPE = 1 / [Q(jω)ⁿ]

Where:

Magnitude and Phase:

|Z_CPE| = 1 / (Qωⁿ)

θ = -n × 90°

Special Cases of CPE

n Value	Phase	Equivalent Element	Physical Meaning
n = 1	-90°	Ideal Capacitor	Perfect dielectric, smooth surface
n = 0.5	-45°	Warburg Element	Semi-infinite diffusion
n = 0	0°	Ideal Resistor	Pure resistance
0.7 < n < 1	-63° to -90°	"Leaky" Capacitor	Typical for rough electrodes

Physical Origins of CPE Behavior

Several factors cause electrodes to exhibit CPE rather than ideal capacitor behavior:

1. Surface Roughness

Real electrode surfaces are not atomically flat. Bumps and valleys create a distribution of local capacitances.

2. Porosity

Porous electrodes have different penetration depths for different frequencies, causing distributed charging.

3. Chemical Heterogeneity

Different chemical composition or oxide thickness across the surface creates varying local properties.

4. Fractal Geometry

Some electrode surfaces have fractal characteristics, leading to CPE exponents related to fractal dimension.

Converting CPE to Effective Capacitance

For circuit analysis, it's often useful to extract an "effective capacitance" from CPE parameters:

Brug Formula (for R-CPE parallel):

C_eff = Q^1/n × R^(1-n)/n

Simplified (when n is close to 1):

C_eff ≈ Q at ω = 1 rad/s

At specific frequency:

C_eff(ω) = Q × ω^(n-1)

CPE in Modified Randles Circuit

A more realistic WFC model replaces the ideal C_dl with a CPE:

         Rs                 Rct
    ────┬────┬────────────┬────┬────
        │    │            │    │
        │    │            │    │
        │  ──┴──        ──┴──  │
        │  │   │        │    │ │
        │  │CPE│        │ Zw │ │  ← CPE replaces Cdl
        │  │Q,n│        │    │ │
        │  ──┬──        ──┬──  │
        │    │            │    │
        └────┴────────────┴────┘

This produces the characteristic depressed semicircle seen in real EIS data.

Typical CPE Values for WFC

Electrode Type	n (typical)	Q (typical)
Polished stainless steel	0.85-0.95	10-50 µF·s(n-1)/cm²
Brushed stainless steel	0.75-0.85	20-100 µF·s(n-1)/cm²
Sandblasted electrode	0.65-0.75	50-200 µF·s(n-1)/cm²
Porous electrode	0.50-0.70	100-1000 µF·s(n-1)/cm²

VIC Design Implications

Why CPE Matters for VIC:

Measuring CPE Parameters

To determine Q and n for your WFC:

1. Perform EIS measurement across relevant frequency range
2. Fit data to modified Randles circuit with CPE
3. Extract Q and n from fitting software
4. Validate by checking phase angle: θ should equal -n × 90°

CPE in VIC Matrix Calculator

The VIC Matrix Calculator can incorporate CPE effects:

• CPE exponent (n): Adjust from the Water Profile or Cole-Cole settings
• Effective capacitance: Calculated at operating frequency
• Loss factor: Related to (1-n), represents energy dissipation

Practical Recommendation: If your WFC electrodes are rough or etched (to increase surface area for gas production), expect significant CPE behavior (n = 0.7-0.85). This will broaden your resonance peak but reduce maximum Q factor. Smooth, polished electrodes (n > 0.9) behave more ideally and allow sharper tuning.

Chapter 3 Complete. Next: VIC Circuit Theory →

VIC Circuit Theory

VIC Introduction

What is a VIC Circuit?

The Voltage Intensifier Circuit (VIC) is a resonant circuit topology designed to develop high voltages across a water fuel cell (WFC) while drawing relatively low current from the source. Originally conceived by Stanley Meyer, the VIC uses the principles of resonance and voltage magnification to create conditions favorable for water dissociation.

The Basic Concept

At its core, the VIC is a series resonant circuit that uses inductors (chokes) and capacitors to magnify voltage. Unlike conventional electrolysis that uses brute-force DC current, the VIC aims to:

• Maximize voltage across the water fuel cell
• Minimize current draw from the power source
• Use resonance to achieve efficient energy transfer
• Exploit the capacitive nature of the water cell

The VIC Block Diagram

    ┌──────────┐     ┌──────┐     ┌──────┐     ┌──────┐     ┌─────────┐
    │  Pulse   │────▶│  L1  │────▶│  C1  │────▶│  L2  │────▶│   WFC   │
    │Generator │     │      │     │      │     │      │     │         │
    └──────────┘     └──────┘     └──────┘     └──────┘     └─────────┘
         ▲             ▲            ▲            ▲              ▲
         │             │            │            │              │
    Frequency     Primary       Tuning      Secondary      Water Fuel
     Control       Choke      Capacitor      Choke           Cell

              PRIMARY SIDE          │         SECONDARY SIDE
              (L1-C1 Tank)          │         (L2-WFC Tank)

Key Components

Component	Symbol	Function
Pulse Generator	—	Provides driving signal at resonant frequency
Primary Choke	L1	Current limiting, energy storage, voltage magnification
Tuning Capacitor	C1	Sets primary resonant frequency with L1
Secondary Choke	L2	Further voltage magnification, resonance with WFC
Water Fuel Cell	WFC	Capacitive load where water dissociation occurs

Operating Principle

Step 1: Pulse Excitation

The pulse generator provides a square wave or pulsed DC signal at or near the resonant frequency of the primary tank circuit (L1-C1).

Step 2: Primary Resonance

The L1-C1 combination resonates, building up voltage across C1 that can be many times the input voltage (determined by Q factor).

Step 3: Energy Transfer

The amplified voltage drives current through L2, which further builds up energy and transfers it to the WFC.

Step 4: Secondary Resonance

If L2 and WFC are tuned together, a second stage of voltage magnification occurs, creating very high voltages across the water.

Step 5: Water Interaction

The high voltage across the WFC creates a strong electric field in the water, affecting the molecular bonds of H₂O.

The "Matrix" Concept

The term "VIC Matrix" refers to the interconnected relationship between all circuit parameters. Everything is connected:

• Changing L1 affects the primary resonant frequency
• The resonant frequency must match the pulse generator
• L2 and WFC capacitance determine secondary resonance
• All inductances and capacitances are linked through the desired frequency
• The Q factors determine voltage magnification at each stage

This is why the VIC Matrix Calculator exists—to help navigate these complex interdependencies.

Circuit Variations

Basic VIC (Two-Choke)

Uses separate L1 and L2 chokes with discrete C1 and WFC capacitance.

Transformer-Coupled VIC

L1 and L2 are wound on the same core, creating transformer action between primary and secondary.

Bifilar VIC

Uses bifilar-wound chokes where two windings are wound together, creating inherent capacitance and magnetic coupling.

Single-Choke VIC

Simplified version where one choke resonates directly with the WFC capacitance.

What Makes VIC Different from Electrolysis?

Parameter	Conventional Electrolysis	VIC Approach
Power Type	DC (constant current)	Pulsed/AC (resonant)
Voltage	1.5-3V (above decomposition)	Hundreds to thousands of volts
Current	High (amps)	Low (milliamps)
Frequency	0 Hz (DC)	kHz to MHz range
WFC View	Resistive load	Capacitive load
Energy Mechanism	Electron transfer	Electric field stress

Goals of VIC Design

1. Maximize Q factor: Higher Q = more voltage magnification
2. Achieve resonance: All components tuned to operating frequency
3. Match impedances: Efficient energy transfer between stages
4. Maintain stability: Prevent frequency drift and oscillation problems
5. Deliver energy to WFC: Create conditions for water molecule stress

Key Insight: The VIC treats water not as a resistive medium to push current through, but as a dielectric capacitor to be charged with high voltage. This fundamental difference drives all aspects of VIC design and is why traditional electrolysis equations don't apply.

Next: Primary Side (L1-C1) Analysis →

Primary Side

Primary Side (L1-C1) Analysis

The primary side of the VIC consists of the first inductor (L1) and tuning capacitor (C1). This stage receives the driving signal and provides the first stage of voltage magnification. Understanding its behavior is crucial for successful VIC design.

Primary Tank Circuit

L1 and C1 form a series resonant tank circuit. At the resonant frequency, this circuit:

• Has minimum impedance (ideally just the DC resistance)
• Draws maximum current from the source
• Develops magnified voltage across L1 and C1

                     R1 (DCR of L1)
                        │
    Pulse      ┌────────┴────────┐
    Generator  │                 │
        ○──────┤      L1         ├────────┬────── To L2
               │                 │        │
               └─────────────────┘       ─┴─
                                         ─┬─ C1
                                          │
                                         ─┴─ GND

    V_in ────▶  [  L1 + R1  ] ────▶ [ C1 ] ────▶ V_out

    At resonance: V_C1 = Q × V_in

Resonant Frequency Calculation

Primary Resonant Frequency:

f₀ = 1 / (2π√(L1 × C1))

Rearranging to Find Components:

L1 = 1 / (4π²f₀²C1)

C1 = 1 / (4π²f₀²L1)

Example Calculations

Target f₀	Given L1	Required C1
10 kHz	1 mH	253 nF
10 kHz	10 mH	25.3 nF
25 kHz	1 mH	40.5 nF
50 kHz	500 µH	20.3 nF

Q Factor of Primary Side

The Q factor determines voltage magnification:

Q Factor:

Q_L1C = (2π × f₀ × L1) / R1 = X_L1 / R1

Voltage Magnification:

V_C1 = Q_L1C × V_in

Example:

Characteristic Impedance

The characteristic impedance of the primary tank affects matching:

Z₀ = √(L1 / C1)

Relationship to Q:

Q = Z₀ / R1

Higher Z₀ (more L, less C) means higher Q for same resistance.

Design Trade-offs

Design Choice	Advantages	Disadvantages
High L1, Low C1	Higher Z₀, potentially higher Q	More wire, higher DCR, harder to wind
Low L1, High C1	Less wire, lower DCR, easier construction	Lower Z₀, may need larger capacitor
High frequency	Smaller components, lower SRF concern	Skin effect losses, harder switching
Low frequency	Lower losses, easier switching	Larger components, SRF may be issue

Current and Power Considerations

At resonance, the circuit draws maximum current:

Resonant Current:

I_res = V_in / R1

Power from Source:

P_in = V_in² / R1 = I_res² × R1

Reactive Power (circulating):

P_reactive = V_C1 × I_res = Q × P_in

Note: The reactive power circulates between L1 and C1 but is not consumed.

Bandwidth and Tuning Sensitivity

The 3dB bandwidth of the primary tank:

BW = f₀ / Q_L1C

Example:

f₀ = 10 kHz, Q = 50 → BW = 200 Hz

The driving frequency must be within ±100 Hz of f₀ for good response.

Practical Implication:

High-Q circuits are sensitive to component tolerances and temperature drift. You may need PLL (Phase-Locked Loop) control to maintain resonance.

Component Selection Guidelines

L1 (Primary Choke)

• Inductance: 100 µH to 100 mH typical
• DCR: As low as practical (determines Q)
• SRF: Should be well above operating frequency (10× minimum)
• Core: Ferrite, iron powder, or air-core depending on frequency
• Wire: Copper preferred; resistance wire reduces Q

C1 (Tuning Capacitor)

• Value: Selected to resonate with L1 at desired frequency
• Voltage rating: Must exceed Q × V_in
• Type: Film (polypropylene, polyester) or ceramic
• ESR: Low ESR for minimal losses
• Temperature stability: NPO/C0G ceramic or film preferred

Practical Assembly Tips

1. Measure L1 accurately: Use an LCR meter at multiple frequencies
2. Start with calculated C1: Then fine-tune for best response
3. Use variable capacitor or parallel caps: For easy tuning
4. Check for SRF: Ensure L1's SRF is well above f₀
5. Monitor temperature: Component values drift with heat

VIC Matrix Calculator: The calculator determines optimal L1 and C1 values based on your target frequency and available components. It also shows the expected Q factor and voltage magnification.

Next: Secondary Side (L2-WFC) Analysis →

Secondary Side

Secondary Side (L2-WFC) Analysis

The secondary side of the VIC consists of the second inductor (L2) and the water fuel cell (WFC) acting as a capacitor. This stage receives the amplified signal from the primary and delivers the final voltage to the water. Proper design of this stage is critical for efficient energy transfer to the WFC.

Secondary Tank Circuit

L2 and the WFC capacitance form the secondary resonant tank:

    From               R2 (DCR of L2)
    Primary      ┌────────┴────────┐
        ○────────┤                 ├────────┬────────○
    (V_C1)       │      L2         │        │       (+)
                 │                 │       ─┴─
                 └─────────────────┘       │ │  WFC
                                           │ │ (C_wfc)
                                           ─┬─
                                            │
        ○───────────────────────────────────┴────────○
                                                    (−)

    V_C1 ────▶ [ L2 + R2 ] ────▶ [ WFC ] ────▶ V_WFC

    At secondary resonance: V_WFC = Q_L2 × V_C1 = Q_L2 × Q_L1C × V_in

The WFC as a Capacitor

The water fuel cell presents a complex impedance, but at VIC frequencies, it behaves predominantly as a capacitor:

WFC Capacitance Components:

At typical VIC frequencies (1-50 kHz), C_wfc is dominated by C_geo.

Secondary Resonant Frequency

Secondary Resonance:

f₀_secondary = 1 / (2π√(L2 × C_wfc))

For Maximum Voltage Transfer:

Ideally, f₀_secondary = f₀_primary

This means: L1 × C1 = L2 × C_wfc

Q Factor of Secondary Side

The secondary Q factor determines the second stage of voltage magnification:

Secondary Q Factor:

Q_L2 = (2π × f₀ × L2) / (R2 + R_wfc)

Where R_wfc is the effective resistance of the WFC (solution resistance + losses).

Total Voltage Magnification:

V_WFC = Q_L1C × Q_L2 × V_in

Example:

Cascaded Resonance Effects

When both stages resonate at the same frequency, the effects multiply:

Configuration	Total Magnification	Notes
Only primary resonance	QL1C	L2-WFC not tuned
Only secondary resonance	QL2	L1-C1 not tuned
Dual resonance	QL1C × QL2	Maximum magnification
Harmonic secondary	Variable	Secondary at 2f₀, 3f₀, etc.

Impedance Matching Considerations

For efficient energy transfer between primary and secondary:

Characteristic Impedance Match:

Z₀_primary = √(L1/C1)

Z₀_secondary = √(L2/C_wfc)

Matching these impedances can improve energy transfer, though it's not always achievable or necessary.

Effect of WFC Properties on Secondary

WFC Parameter	Effect on Secondary	Design Response
Higher Cwfc	Lower f₀, lower Z₀	Increase L2 or reduce C1
Higher Rwfc	Lower QL2	Use purer water or optimize gap
Larger electrode area	Higher Cwfc	Requires larger L2
Narrower gap	Higher Cwfc, lower Rwfc	Trade-off between C and R

Bifilar Choke Considerations

When L2 is bifilar wound (or when L1 and L2 are wound together as a bifilar pair):

• Inherent capacitance: The bifilar winding has capacitance between turns
• Magnetic coupling: Energy transfers inductively between windings
• Lower SRF: The inter-winding capacitance lowers self-resonant frequency
• Complex tuning: The system becomes a coupled resonator

Calculating L2 for Given WFC

Given: Target frequency and WFC capacitance

L2 = 1 / (4π²f₀²C_wfc)

Example:

Sanity check: This is a reasonable inductance, achievable with ~500-1000 turns on a ferrite core.

Power Delivery to WFC

The actual power delivered to the WFC depends on its resistive component:

Power in WFC Resistance:

P_wfc = I²_wfc × R_wfc

Where:

I_wfc = V_WFC / Z_wfc ≈ V_WFC × ω × C_wfc

This power heats the water and drives electrochemical reactions.

Voltage Distribution Across WFC

The high voltage across the WFC creates an electric field:

Electric Field in WFC:

E = V_WFC / d

Where d is the electrode gap.

Example:

This is a substantial electric field that can influence molecular behavior in water.

Design Guidelines for L2

1. Match resonant frequency: L2 should resonate with C_wfc at the same frequency as L1-C1
2. Minimize DCR: R2 directly reduces Q_L2 and thus voltage magnification
3. Consider coupling: If using transformer-coupled design, mutual inductance matters
4. Account for WFC changes: C_wfc varies with temperature, voltage, and bubble formation
5. Leave tuning margin: Design L2 slightly higher, fine-tune with small series capacitor if needed

Key Insight: The secondary side is where VIC theory meets reality. The WFC is not an ideal capacitor—it has losses, frequency-dependent behavior, and changes during operation. Successful VIC design must account for these real-world effects.

Next: Resonant Charging Principle →

Resonant Charging

Resonant Charging Principle

Resonant charging is a technique where energy is transferred to a capacitive load (the WFC) in a controlled, oscillatory manner. Unlike direct DC charging, resonant charging can achieve higher efficiency and allows voltage magnification beyond the source voltage.

Conventional vs. Resonant Charging

Aspect	DC Charging (R-C)	Resonant Charging (L-C)
Final voltage	= Vsource	Can exceed Vsource (up to 2× for half-wave)
Energy efficiency	50% max (half lost in R)	Can approach 100% (minimal loss in L)
Charging curve	Exponential (slow)	Sinusoidal (faster)
Peak current	V/R at start	V/Z₀ (controlled by L)

Basic Resonant Charging Circuit

     Switch (S)
    ────○/○────┬───────────────┬────
               │               │
    V_source   │               │
        +      │    ┌─────┐   ─┴─
               │    │  L  │   ─┬─ C (WFC)
               │    └──┬──┘    │
               │       │       │
    ───────────┴───────┴───────┴────
                              GND

    When S closes:
    1. Current builds in L (energy stored in magnetic field)
    2. Current flows into C, charging it
    3. Voltage on C rises
    4. At peak voltage, current reverses (or S opens)

Half-Cycle Resonant Charging

In half-cycle mode, the switch opens when capacitor voltage reaches maximum:

Ideal Half-Cycle Charging (lossless):

V_C,max = 2 × V_source

Charging Time:

t_charge = π√(LC) = π/ω₀ = 1/(2f₀)

This is exactly half the resonant period.

Why 2× Voltage?

Energy Conservation:

Resonant Charging with Losses

Real circuits have losses that reduce the voltage gain:

With Resistance (damped case):

V_C,max = V_source × (1 + e^{-πR/(2√(L/C))})

V_C,max = V_source × (1 + e^-π/(2Q))

Approximation for high Q:

V_C,max ≈ 2V_source × (1 - π/(4Q))

Voltage Gain vs. Q Factor

Q Factor	VC,max/Vsource	Efficiency
∞ (ideal)	2.00	100%
100	1.98	98.4%
50	1.97	96.9%
20	1.92	92.5%
10	1.85	85.5%
5	1.73	73%

Continuous Resonant Excitation

In the VIC, instead of single pulses, we drive the circuit continuously at the resonant frequency:

Steady-State Resonance:

Energy from the source compensates for losses each cycle, maintaining a steady oscillation amplitude.

Voltage Magnification:

V_C = Q × V_source

This is much greater than the 2× from single-pulse resonant charging when Q > 2.

Resonant Charging in VIC Context

The VIC uses resonant charging principles in several ways:

1. Primary tank: C1 is resonantly charged through L1
2. Secondary transfer: Energy transfers resonantly to WFC through L2
3. Cumulative effect: Multiple stages multiply the magnification

Timing and Switching

For optimal resonant charging:

Critical Timing Points:

Zero-Current Switching (ZCS):

Turning off when current is zero minimizes switching losses and eliminates inductive kick.

Energy Flow Analysis

    Time →

    V_C:    ────╱╲    ╱╲    ╱╲    ╱╲────
               ╱  ╲  ╱  ╲  ╱  ╲  ╱  ╲
              ╱    ╲╱    ╲╱    ╲╱    ╲

    I_L:    ──╱╲    ╱╲    ╱╲    ╱╲────
             ╱  ╲  ╱  ╲  ╱  ╲  ╱  ╲
            ╱    ╲╱    ╲╱    ╲╱    ╲

    Energy in C: High → Low  → High → Low
    Energy in L: Low  → High → Low  → High

    Total energy (minus losses) remains constant in steady state.

Advantages of Resonant Charging for WFC

• High voltage: Achieves voltages beyond source capability
• Low current draw: Source only provides loss compensation
• Controlled energy delivery: Sinusoidal rather than impulsive
• Efficient: Minimal resistive losses when Q is high
• Self-limiting: Voltage limited by Q factor, not infinite

Key Principle: Resonant charging exploits the energy storage capability of inductors and capacitors. By timing the energy injection to match the natural oscillation, we can build up substantial energy in the circuit with modest input power—the same principle used in pushing a swing at just the right moment.

Next: Step-Charging Ladder Effect →

Step Charging

Step-Charging Ladder Effect

Step-charging, also known as the "staircase" or "ladder" effect, refers to the progressive buildup of voltage across a capacitor through successive resonant pulses. This technique can achieve voltage levels far beyond what single-pulse resonant charging provides.

The Concept

Instead of maintaining continuous oscillation, step-charging applies discrete pulses that each add a quantum of energy to the capacitor. The voltage builds up incrementally:

    Voltage
       ↑
       │                                    ┌───
       │                               ┌───┘
       │                          ┌───┘
       │                     ┌───┘
       │                ┌───┘
       │           ┌───┘
       │      ┌───┘
       │ ┌───┘
       │─┘
       └─────────────────────────────────────→ Time
         ↑   ↑   ↑   ↑   ↑   ↑   ↑   ↑   ↑
        Pulse Pulse Pulse ...
        1     2     3

    Each pulse adds approximately 2×V_source to capacitor voltage
    (in ideal lossless case with unidirectional diode)

How Step-Charging Works

Step-by-Step Process:

Circuit for Step-Charging

         Switch
    V_s ──○/○───┬───────────────┬────▶│────┬────
               │               │      D     │
               │    ┌─────┐    │           ─┴─
               │    │  L  │   ─┴─          ─┬─ C (WFC)
               │    └──┬──┘   ─┬─           │
               │       │       │            │
    ───────────┴───────┴───────┴────────────┴────

    D = Diode prevents reverse current
    C charges in discrete steps

Voltage After N Pulses

Ideal Case (no losses):

V_C,N = 2N × V_source

With Losses (exponential decay factor):

V_C,N = 2V_s × Σ(e^-π/(2Q))^k for k=0 to N-1

Converges to Maximum:

V_C,max = 2V_s / (1 - e^-π/(2Q))

For high Q: V_C,max ≈ (4Q/π) × V_source

Maximum Voltage vs. Q Factor

Q Factor	Vmax/Vsource	Pulses to 90%
10	~12.7	~6
20	~25.5	~12
50	~63.7	~30
100	~127	~60

Comparison: Continuous vs. Step Charging

Aspect	Continuous Resonance	Step Charging
Max voltage	Q × Vs (AC peak)	(4Q/π) × Vs (DC)
Waveform	Sinusoidal	Staircase
Power delivery	Constant	Pulsed
Complexity	Simpler	Needs diode/timing

Step-Charging in VIC Systems

Meyer's designs allegedly used step-charging principles:

• Unidirectional charging: Diode prevents energy return to source
• Pulse timing: Gated pulses at resonant frequency
• Voltage accumulation: Progressive buildup across WFC
• Controlled discharge: Occasional reset or bleed-off of accumulated voltage

Pulse Train Design

Optimal Pulse Parameters:

Energy Considerations

Energy Stored After N Pulses:

E_C,N = ½C(V_C,N)² = ½C(2NV_s)² = 2CN²V_s²

Energy Delivered per Pulse:

ΔE = E_C,N - E_C,N-1 = 2CV_s²(2N-1)

Each successive pulse adds more energy because it's working against a higher voltage!

Practical Implementation

Driver Circuit Requirements:

1. High-speed switching: MOSFET or IGBT driver
2. Precise timing: Microcontroller or pulse generator
3. High-voltage diode: Fast recovery, rated for expected voltages
4. Voltage monitoring: Feedback to prevent over-voltage

Safety Considerations:

• Voltages can reach dangerous levels quickly
• Energy stored in capacitor can be lethal
• Include bleed resistor for safe discharge
• Implement hardware over-voltage protection

VIC Matrix Simulation

The VIC Matrix Calculator can simulate step-charging behavior:

• Step-charge simulation: Predicts voltage after N pulses
• Loss modeling: Accounts for resistance and dielectric losses
• Time to saturation: How many pulses to reach maximum voltage
• Energy efficiency: Tracks energy delivered vs. stored

Key Insight: Step-charging combines the voltage doubling of resonant charging with the cumulative effect of multiple pulses. With sufficient Q factor, extremely high voltages can be developed across the WFC—voltages that would be impossible to achieve directly from the source.

Chapter 4 Complete. Next: Choke Design & Construction →

Choke Design

Choke Fundamentals

Inductor/Choke Fundamentals

Inductors, commonly called "chokes" in VIC terminology, are the workhorses of the resonant circuit. They store energy in their magnetic field and, together with capacitors, determine the resonant frequency and voltage magnification capability of the VIC.

What is an Inductor?

An inductor is a passive electrical component that stores energy in a magnetic field when current flows through it. The fundamental properties are:

Inductance (L):

Measured in Henries (H), inductance quantifies the magnetic flux linkage per unit current:

L = NΦ/I = N²μA/l

Where:

Key Inductor Parameters

Parameter	Symbol	Units	Importance
Inductance	L	Henry (H)	Determines resonant frequency with C
DC Resistance	DCR, Rdc	Ohms (Ω)	Limits Q factor and causes losses
Self-Resonant Frequency	SRF	Hz	Must be > operating frequency
Quality Factor	Q	Dimensionless	Ratio of reactance to resistance
Saturation Current	Isat	Amps (A)	Max current before inductance drops

Inductor Construction

A practical inductor consists of:

1. Wire: Conductor wound into coils (turns)
2. Core: Material inside the coil (air, ferrite, iron, etc.)
3. Form: Structure that holds the winding

Types of Cores

Core Type	Permeability	Frequency Range	VIC Application
Air core	1 (reference)	Any (no losses)	High-Q, low inductance
Iron powder	10-100	Up to ~10 MHz	Good for VIC frequencies
Ferrite	100-10000	10 kHz - 100 MHz	Most common for VIC
Laminated iron	1000-10000	50/60 Hz to ~10 kHz	Lower VIC frequencies

Inductance Formulas

Single-Layer Solenoid (air core):

L = (N²μ₀A)/l = (N²r²)/(9r + 10l) µH

Where r and l are in inches (Wheeler's formula)

With Magnetic Core:

L = A_L × N² (nH)

Where A_L is the inductance factor of the core (nH/turn²)

Toroidal Core:

L = (μ₀μ_rN²A) / (2πr_mean)

DC Resistance (DCR)

The DC resistance is determined by the wire properties:

R_dc = ρ × l_wire / A_wire

Where:

Q Factor of Inductors

Inductor Q Factor:

Q = ωL/R = 2πfL/R_total

R_total includes:

Self-Resonant Frequency (SRF)

Every inductor has parasitic capacitance between turns and layers:

SRF = 1 / (2π√(LC_parasitic))

Design Rule:

SRF should be at least 10× the operating frequency.

At frequencies above SRF, the inductor acts like a capacitor!

VIC Choke Design Goals

1. Target inductance: Sets resonant frequency with capacitor
2. Low DCR: Maximizes Q factor
3. High SRF: Ensures proper operation at intended frequency
4. Adequate current rating: Won't saturate or overheat
5. Appropriate core: Low losses at operating frequency

Key Tradeoff: More turns = more inductance, but also more wire = more DCR. The design challenge is achieving the target inductance with minimum resistance, which means selecting appropriate wire gauge, core material, and winding technique.

Next: Core Materials & Properties →

Core Materials

Core Materials & Properties

The core material of an inductor dramatically affects its performance. Choosing the right core is essential for achieving the desired inductance, Q factor, and frequency response in VIC applications.

Why Use a Core?

A magnetic core increases inductance by providing a low-reluctance path for magnetic flux:

L = μ₀μᵣN²A/l

The relative permeability (μᵣ) of the core multiplies the inductance compared to an air core.

Core Material Comparison

Material	μᵣ (typical)	Frequency Range	Saturation	Cost
Air	1	Any	N/A	Free
Iron Powder	10-100	1 kHz - 100 MHz	High (0.5-1.5T)	Low
Ferrite (MnZn)	1000-10000	1 kHz - 1 MHz	Low (0.3-0.5T)	Medium
Ferrite (NiZn)	50-1500	100 kHz - 500 MHz	Low (0.3-0.4T)	Medium
Laminated Silicon Steel	2000-6000	50 Hz - 10 kHz	High (1.5-2.0T)	Low
Amorphous Metal	10000-100000	50 Hz - 100 kHz	High (1.5T)	High
Nanocrystalline	15000-100000	1 kHz - 1 MHz	High (1.2T)	High

Core Losses

All magnetic cores dissipate energy through two mechanisms:

1. Hysteresis Loss

Energy lost each time the core is magnetized and demagnetized.

P_h ∝ f × B_maxⁿ (n ≈ 1.6-2.5)

Proportional to frequency and flux density.

2. Eddy Current Loss

Circulating currents induced in the core material.

P_e ∝ f² × B_max²

Proportional to frequency squared - dominates at high frequencies.

Steinmetz Equation

P_core = k × f^α × B^β × Volume

Where k, α, β are material-specific constants from datasheets.

Ferrite Materials for VIC

Ferrites are the most common choice for VIC frequencies (1-50 kHz):

Material	μᵢ	Optimal Frequency	Application
3C90 (TDK)	2300	25-200 kHz	Power transformers
N87 (EPCOS)	2200	25-500 kHz	General purpose
N97 (EPCOS)	2300	25-150 kHz	Low loss
3F3 (Ferroxcube)	2000	100-500 kHz	Higher frequency
77 Material (Fair-Rite)	2000	Up to 1 MHz	EMI/RFI suppression

Iron Powder Cores

Micrometals and Amidon iron powder cores are popular for their:

• High saturation flux density
• Gradual saturation (soft saturation)
• Good temperature stability
• Self-gapping (distributed gap)

Common Iron Powder Mixes

Mix	μ	Color	Frequency Range
Mix 26	75	Yellow/White	DC - 1 MHz
Mix 52	75	Green/Blue	DC - 3 MHz
Mix 2	10	Red/Clear	1 - 30 MHz
Mix 6	8	Yellow	10 - 50 MHz

Core Shapes

Toroidal

Doughnut shape with closed magnetic path. Excellent flux containment, low EMI. Harder to wind but very efficient.

E-Core / EI-Core

E-shaped halves that mate together. Easy to wind on bobbin. Can add air gap easily.

Pot Core

Cylindrical with center post. Shields winding from external fields. Good for sensitive applications.

Rod Core

Simple cylindrical rod. Open magnetic path, lower inductance per turn but no saturation issues.

Core Saturation

When the magnetic flux density exceeds the saturation limit:

• Permeability drops dramatically
• Inductance decreases
• Current increases rapidly
• Core heating increases

Avoiding Saturation:

B_peak = (L × I_peak) / (N × A_e) < B_sat

Always check that peak flux density stays below saturation limit of your core material.

Recommendations for VIC

Frequency Range	Recommended Core	Notes
1-10 kHz	N97/3C90 ferrite or iron powder	Low loss at these frequencies
10-50 kHz	N87/3F3 ferrite	Good balance of μ and loss
50-200 kHz	3F3/3F4 ferrite or Mix 26 powder	Lower permeability, lower loss
>200 kHz	NiZn ferrite or Mix 2 powder	Designed for high frequency

VIC Matrix Calculator: The Choke Design module includes a core database with A_L values and frequency recommendations. Select your core and it will calculate the required turns for your target inductance.

Next: Wire Gauge & Material Selection →

Wire Selection

Wire Gauge & Material Selection

The wire used to wind an inductor directly affects its DC resistance, current capacity, and Q factor. Proper wire selection is essential for maximizing VIC circuit performance.

Wire Gauge Systems

Wire size is commonly specified using the American Wire Gauge (AWG) system:

AWG	Diameter (mm)	Area (mm²)	Ω/m (Copper)	Max Current (A)
18	1.024	0.823	0.0210	2.3
20	0.812	0.518	0.0333	1.5
22	0.644	0.326	0.0530	0.92
24	0.511	0.205	0.0842	0.58
26	0.405	0.129	0.1339	0.36
28	0.321	0.081	0.2128	0.23
30	0.255	0.051	0.3385	0.14
32	0.202	0.032	0.5383	0.09

Note: AWG follows logarithmic progression. Each 3 AWG steps doubles resistance, halves area.

Wire Materials

Material	Resistivity (×10⁻⁸ Ω·m)	Relative to Copper	Use Case
Copper	1.68	1.0× (reference)	Best for high Q
Aluminum	2.65	1.6×	Lightweight applications
SS304	72	~43×	Corrosion resistance
SS316	74	~44×	Better corrosion resistance
SS430 (Ferritic)	~100	~60×	Magnetic, high resistance
Nichrome (80/20)	108	~64×	Heating elements, damping
Kanthal A1	145	~86×	High-temp resistance wire

Effect of Material on Q Factor

Q Factor Relationship:

Q = 2πfL / R

Since R is proportional to resistivity, using high-resistivity wire dramatically reduces Q:

When to Use Resistance Wire

Despite lower Q, resistance wire has valid uses:

• Current limiting: Built-in current limit without separate resistor
• Damping: Prevents excessive ringing
• Safety: Limits power in fault conditions
• Meyer's designs: Some original VIC designs used stainless steel wire

Warning: Using resistance wire in a resonant circuit dramatically reduces voltage magnification. A Q of 2 means you only get 2× voltage gain instead of 50× or 100× with copper.

Skin Effect

At high frequencies, current flows primarily near the wire surface:

Skin Depth (δ):

δ = √(ρ / (π × f × μ₀ × μᵣ))

For Copper:

δ(mm) ≈ 66 / √f(Hz)

Skin Effect Mitigation

• Litz wire: Multiple thin insulated strands twisted together
• Flat/ribbon wire: More surface area for same cross-section
• Use finer gauge: If wire radius ≈ δ, skin effect is minimal

Magnet Wire Types

Insulation Type	Temp Rating	Voltage Rating	Notes
Polyurethane (solderable)	130°C	~100V/layer	Can solder through coating
Polyester-imide	180°C	~200V/layer	Good general purpose
Polyamide-imide	220°C	~300V/layer	High temp applications
Heavy build (HN)	Various	~500V/layer	Thicker insulation
Triple insulated	Various	~3000V	Safety-rated isolation

Wire Selection Guidelines for VIC

For Maximum Q (recommended):

For Current-Limited Applications:

Calculating Wire Length

Wire Length for N Turns:

l_wire ≈ N × π × d_coil

Where d_coil is the average coil diameter.

Resulting DCR:

R_dc = ρ × l_wire / A_wire

VIC Matrix Calculator: The Choke Design tool automatically calculates DCR based on your wire gauge, material, and number of turns. It shows the resulting Q factor and voltage magnification for your design.

Next: Bifilar Winding Technique →

Bifilar Windings

Bifilar Winding Technique

Bifilar winding is a special technique where two wires are wound together in parallel on a core. This configuration creates unique electromagnetic properties that are particularly relevant to VIC designs, including inherent capacitance between windings and special transformer-like coupling.

What is Bifilar Winding?

In a bifilar winding, two conductors are wound side-by-side along the entire length of the coil:

    Standard Winding:         Bifilar Winding:

       ─────────────            ═══════════════
       │ │ │ │ │ │              ║A║B║A║B║A║B║
       └─┘ └─┘ └─┘              ╚═╝ ╚═╝ ╚═╝

    Single wire wound           Two wires (A & B)
    around core                 wound together

    Cross-section view:

    Standard:    Bifilar:
      ○ ○ ○       ○ ● ○ ●
      ○ ○ ○       ● ○ ● ○

    ○ = Wire A    ● = Wire B

Bifilar Winding Properties

Property	Effect	VIC Relevance
High inter-winding capacitance	Built-in C between A and B	May replace discrete capacitor
Near-unity coupling	k ≈ 1 between windings	Efficient energy transfer
Cancellation modes	Some flux cancellation possible	Affects net inductance
Lower SRF	High Cparasitic reduces SRF	Consider in frequency selection

Connection Configurations

1. Series Aiding (Same Direction):

End of A connects to start of B → Fluxes add

L_total = L_A + L_B + 2M ≈ 4L (for k=1)

2. Series Opposing (Opposite Direction):

End of A connects to end of B → Fluxes subtract

L_total = L_A + L_B - 2M ≈ 0 (for k=1)

3. Parallel Connection:

Starts connected, ends connected → Current splits

L_total = L/2 (for identical windings)

4. Transformer Mode:

A is primary, B is secondary → Voltage transformation

V_B/V_A = N_B/N_A = 1 (for bifilar)

Calculating Bifilar Capacitance

Approximate Inter-Winding Capacitance:

C_winding ≈ ε₀ε_r × (l_wire × d_wire) / s

Where:

Typical Values:

For magnet wire on ferrite: 10-100 pF per meter of winding

Bifilar in VIC Context

Meyer's designs reportedly used bifilar chokes in several ways:

As Primary/Secondary Pair

L1 and L2 wound as bifilar on same core:

• Tight coupling between primary and secondary
• Built-in capacitance may serve as C1
• Simpler construction (single winding operation)

As Choke Sets

Matched pairs for symmetrical circuits:

• Identical L values guaranteed
• Common-mode rejection possible
• Push-pull drive configurations

Winding Techniques

Tips for Bifilar Winding:

Measuring Bifilar Parameters

Measurement	Configuration	What It Tells You
LA alone	Measure A, B open	Inductance of winding A
Lseries-aid	A end to B start, measure	LA + LB + 2M
Lseries-opp	A end to B end, measure	LA + LB - 2M
Cwinding	Measure C between A and B	Inter-winding capacitance

Calculating Coupling Coefficient:

M = (L_series-aid - L_series-opp) / 4

k = M / √(L_A × L_B)

For true bifilar winding: k ≈ 0.95-0.99

Advantages and Disadvantages

Advantages:

Disadvantages:

VIC Matrix Calculator: The Choke Design section includes options for bifilar windings. It can calculate the expected inter-winding capacitance and adjust the SRF estimate accordingly. When designing bifilar chokes, the calculator helps ensure compatibility with your target resonant frequency.

Next: Parasitic Capacitance & SRF →

Parasitic Effects

Parasitic Capacitance & SRF

Real inductors have parasitic capacitance between turns and layers that limits their useful frequency range. Understanding these effects is critical for VIC design, as they determine the maximum operating frequency and affect circuit tuning.

Sources of Parasitic Capacitance

Parasitic capacitance in inductors comes from several sources:

1. Turn-to-Turn Capacitance (C_tt)

Capacitance between adjacent turns in the same layer. Depends on wire spacing and insulation.

2. Layer-to-Layer Capacitance (C_ll)

Capacitance between winding layers. Often the largest contributor in multi-layer coils.

3. Winding-to-Core Capacitance (C_wc)

Capacitance between the winding and the magnetic core (if conductive or grounded).

4. Winding-to-Shield Capacitance

In shielded inductors, capacitance to the external shield.

Self-Resonant Frequency (SRF)

The parasitic capacitance resonates with the inductance at the Self-Resonant Frequency:

SRF = 1 / (2π√(L × C_parasitic))

Behavior at SRF:

Above SRF:

The "inductor" behaves as a capacitor! Impedance decreases with frequency.

Impedance vs. Frequency

    |Z|
     ↑
     │                    ╱╲
     │                   ╱  ╲     ← Peak at SRF
     │                  ╱    ╲
     │                 ╱      ╲
     │               ╱         ╲
     │             ╱            ╲
     │           ╱               ╲
     │         ╱                  ╲
     │       ╱                     ╲
     │     ╱                        ╲
     │   ╱   Inductive region        ╲ Capacitive region
     │ ╱      |Z| = 2πfL              ╲ |Z| = 1/(2πfC)
     └────────────────────────────────────────────→ f
                          SRF

    Phase:  +90° ───────────┬─────────── −90°
                           0° (at SRF)

Operating Frequency Guidelines

fop / SRF	Behavior	Recommendation
< 0.1 (< 10%)	Nearly ideal inductor	Preferred range
0.1 - 0.3 (10-30%)	Slight inductance increase	Acceptable with correction
0.3 - 0.7 (30-70%)	Significant deviation	Caution - Q drops
> 0.7 (> 70%)	Near or past SRF	Do not use

Effective Inductance Near SRF

As frequency approaches SRF, the apparent inductance increases:

L_eff = L_dc / [1 - (f/SRF)²]

Example:

Minimizing Parasitic Capacitance

Winding Techniques:

Design Choices:

Calculating Parasitic Capacitance

Turn-to-Turn Capacitance (Simplified)

C_tt ≈ ε₀ε_r × l_turn × d_wire / s

Where s is the spacing between adjacent turn centers.

Layer-to-Layer Capacitance

C_ll ≈ ε₀ε_r × A_layer / t_insulation

Where A_layer is the overlapping area between layers.

Total Parasitic Capacitance

The total equivalent capacitance is complex because the distributed capacitances see different voltages. For a rough estimate:

C_parasitic ≈ C_ll/3 + C_tt/N

The 1/3 factor accounts for voltage distribution across layers.

Measuring SRF

Method 1: Impedance Analyzer

1. Connect inductor to impedance analyzer
2. Sweep frequency and plot |Z|
3. SRF is where impedance peaks

Method 2: Signal Generator + Oscilloscope

1. Connect inductor in series with known resistor
2. Drive with sine wave, sweep frequency
3. Monitor voltage across inductor
4. SRF is where voltage peaks (current minimum)

Method 3: Resonance with Known Capacitor

1. Measure inductance at low frequency
2. Add known capacitor in parallel
3. Find new resonant frequency
4. Calculate parasitic C from the difference

SRF in VIC Design

Problem	Symptom	Solution
Operating too close to SRF	Resonance frequency higher than calculated	Reduce tuning cap or use different choke
Operating above SRF	No resonance, circuit acts capacitive	Must redesign with fewer turns
Low SRF in bifilar winding	Limited usable frequency range	Accept limitation or use separate chokes

VIC Matrix Calculator: The Choke Design module estimates SRF based on winding geometry and displays a warning if your operating frequency is too close to SRF. It also calculates the effective inductance at your operating frequency.

Next: DC Resistance and Q Factor →

DCR Effects

DC Resistance and Q Factor

The DC resistance (DCR) of an inductor is the primary factor limiting its Q factor and thus the voltage magnification achievable in a VIC circuit. Understanding and minimizing DCR is essential for high-performance designs.

What is DCR?

DCR is simply the resistance of the wire used to wind the inductor, measured with direct current:

R_dc = ρ × l_wire / A_wire

Where:

DCR and Inductor Design

For a given inductance, DCR depends on the design choices:

Design Change	Effect on L	Effect on DCR	Net Q Effect
More turns	L ∝ N²	R ∝ N	Q ∝ N (improves)
Larger wire gauge	No change	R decreases	Q improves
Higher μ core	L increases	Fewer turns needed	Variable*
Larger core	L increases	Longer mean turn	Often improves
Copper vs. SS wire	No change	R × 40-60	Q ÷ 40-60

*Core losses may offset wire resistance reduction at high frequencies

Q Factor Calculation

Q Factor at Operating Frequency:

Q = 2πfL / R_total

Total Resistance includes:

R_total = R_dc + R_skin + R_proximity + R_core

At low frequencies, R_dc dominates. At high frequencies, skin effect and core losses become significant.

Voltage Magnification Impact

Since voltage magnification equals Q at resonance:

Example Comparison:

Measuring DCR

Method 1: Multimeter

• Simple and quick
• Set meter to lowest resistance range
• Subtract lead resistance
• Accuracy: ±1-5%

Method 2: 4-Wire (Kelvin) Measurement

• Eliminates lead resistance error
• Required for low DCR (<1 Ω)
• Uses separate sense and current leads
• Accuracy: ±0.1%

Method 3: LCR Meter

• Measures L and DCR together
• Can measure at different frequencies
• Shows equivalent series resistance (ESR)
• Best for complete characterization

Optimizing DCR

Design Strategies:

Practical Limits:

Temperature Effects

Wire resistance increases with temperature:

R(T) = R_20°C × [1 + α(T - 20)]

Where α ≈ 0.00393 /°C for copper

Example:

At 80°C: R = R_20°C × 1.24 (+24% increase)

This means Q drops by ~20% when the choke heats up!

DCR in the VIC System

The total resistance in a VIC circuit includes:

Source	Typical Range	Mitigation
L1 DCR	1-50 Ω	Optimize winding
L2 DCR	1-50 Ω	Optimize winding
Capacitor ESR	0.01-1 Ω	Use low-ESR caps
WFC solution resistance	10-10000 Ω	Electrode design, electrolyte
Connection resistance	0.01-1 Ω	Solid connections
Driver output resistance	0.1-10 Ω	Low Rds(on) MOSFETs

Practical Example

Target: 10 mH inductor at 10 kHz with Q > 50

Required R_max:

Q = 2πfL/R → R = 2πfL/Q = 2π × 10000 × 0.01 / 50 = 12.6 Ω

Wire selection (100 turns on 25mm toroid):

Mean turn length ≈ 80mm, total wire = 8m

Result: 22-30 AWG copper all meet the requirement. 22 AWG gives highest Q but may be harder to wind.

VIC Matrix Calculator: Enter your wire gauge and material in the Choke Design tool. It calculates DCR automatically and shows how it affects Q factor and voltage magnification. The calculator warns if your DCR is too high for effective resonance.

Chapter 5 Complete. Next: Water Fuel Cell Design →

Water Fuel Cell Design

WFC Introduction

Water Fuel Cell Basics

The Water Fuel Cell (WFC) is the heart of the VIC system—the component where electrical energy interacts with water. Understanding the WFC as an electrical component is essential for successful VIC circuit design.

What is a Water Fuel Cell?

A Water Fuel Cell consists of electrodes immersed in water, forming an electrochemical cell. Unlike conventional electrolysis cells designed for maximum current flow, the WFC in a VIC is treated as a capacitive load designed for maximum voltage development.

Basic WFC Components:

WFC as an Electrical Component

Electrically, the WFC presents a complex impedance with both capacitive and resistive components:

    Simplified WFC Equivalent Circuit:

         ┌────────────────────────────────────┐
         │                                    │
    (+)──┤   ┌─────┐    ┌─────┐    ┌─────┐   ├──(−)
         │   │C_edl│    │R_sol│    │C_edl│   │
         │   │     │    │     │    │     │   │
         │   └──┬──┘    └──┬──┘    └──┬──┘   │
         │      │          │          │      │
         │      └────┬─────┴─────┬────┘      │
         │           │           │           │
         │          ─┴─         ─┴─          │
         │          ─┬─ C_geo   ─┬─ R_leak   │
         │           │           │           │
         └───────────┴───────────┴───────────┘

    C_edl = Electric double layer capacitance (each electrode)
    R_sol = Solution resistance (water conductivity)
    C_geo = Geometric capacitance (parallel plate effect)
    R_leak = Leakage/Faradaic resistance

Capacitive vs. Resistive Behavior

Frequency	Dominant Behavior	Phase Angle	VIC Relevance
DC (0 Hz)	Resistive	0°	Conventional electrolysis
Low (1-100 Hz)	Mixed R-C	-20° to -60°	Transition region
Medium (100 Hz - 50 kHz)	Primarily capacitive	-60° to -85°	VIC operating range
High (>50 kHz)	Capacitive	-85° to -90°	Nearly ideal capacitor

Common WFC Configurations

1. Parallel Plate

Two flat plates facing each other with water between them.

2. Concentric Tubes

Inner and outer cylinders with water in the annular gap.

3. Tube Array

Multiple concentric tube pairs in parallel.

4. Spiral/Wound

Flat electrodes wound in a spiral with separator.

Key WFC Parameters

Parameter	Symbol	Typical Range	Effect
Electrode Area	A	10-1000 cm²	C ∝ A, affects gas production
Electrode Gap	d	0.5-5 mm	C ∝ 1/d, R ∝ d
Capacitance	Cwfc	1-100 nF	Sets resonant frequency with L2
Solution Resistance	Rsol	10 Ω - 10 kΩ	Affects Q factor

Water Properties Matter

The water used in the WFC significantly affects electrical behavior:

Water Type	Conductivity	Rsol	Notes
Deionized	<1 µS/cm	Very high	Nearly pure capacitor
Distilled	1-10 µS/cm	High	Low losses
Tap water	100-800 µS/cm	Medium	Variable by location
With NaOH/KOH	>10000 µS/cm	Low	Traditional electrolyte

VIC vs. Traditional Electrolysis

Traditional Electrolysis:

VIC Approach:

Key Insight: In VIC design, the WFC is treated primarily as a capacitor whose value must be matched to the choke inductance for resonance. The resistive component should be minimized for high Q, but some resistance is always present due to water's ionic conductivity.

Next: Electrode Geometry & Spacing →

Electrode Geometry

Electrode Geometry & Spacing

The physical design of WFC electrodes directly determines its electrical characteristics—capacitance, resistance, and field distribution. Proper geometry is essential for achieving target resonant frequencies and efficient operation.

Parallel Plate Electrodes

The simplest configuration with straightforward calculations:

Capacitance:

C = ε₀ε_rA / d

For Water (ε_r ≈ 80):

C (pF) ≈ 708 × A(cm²) / d(mm)

Example:

Concentric Tube Electrodes

Cylindrical geometry provides more surface area:

Capacitance:

C = 2πε₀ε_rL / ln(r_outer/r_inner)

Simplified (for small gap relative to radius):

C ≈ ε₀ε_r × 2πr_avgL / d

Where d = r_outer - r_inner

Example:

Tube Array Configurations

Multiple tubes in parallel increase total capacitance:

    Top View of 9-Tube Array:

           ┌───┐
         ┌─┤   ├─┐
       ┌─┤ └───┘ ├─┐
     ┌─┤ └───────┘ ├─┐
   ┌─┤ └───────────┘ ├─┐
   │ └───────────────┘ │
   │   Alternating     │
   │   + and − tubes   │
   └───────────────────┘

    Each concentric pair adds to total capacitance.
    C_total = C₁ + C₂ + C₃ + ... (tubes in parallel)

Electrode Spacing Trade-offs

Gap Size	Capacitance	Resistance	Field Strength	Practical Issues
Very small (<0.5 mm)	Very high	Low	Very high	Bubble blocking, arcing risk
Small (0.5-1.5 mm)	High	Medium-low	High	Sweet spot
Medium (1.5-3 mm)	Medium	Medium	Medium	Easy to build
Large (>3 mm)	Low	High	Low	Needs more voltage

Electric Field Calculation

Field Strength (uniform field approximation):

E = V / d

Example:

Note: Water breakdown occurs at ~30-70 MV/m, so typical VIC fields are well below breakdown.

Surface Area Considerations

Larger electrode area provides:

• Higher capacitance (more energy storage)
• Lower current density (longer electrode life)
• More sites for gas evolution
• Better heat dissipation

But requires:

• Larger choke inductance (to maintain resonant frequency)
• More water volume
• Larger enclosure

Dimensional Design Process

Step 1: Determine Target Capacitance

From resonant frequency and available inductance:

C_target = 1 / (4π²f₀²L₂)

Step 2: Choose Geometry Type

Plates, tubes, or array based on available materials and space.

Step 3: Select Gap Distance

Balance capacitance needs with practical concerns (1-2 mm typical).

Step 4: Calculate Required Area

A = C × d / (ε₀ε_r)

Step 5: Dimension the Electrodes

For plates: Choose L × W. For tubes: Choose radius and length.

Practical Design Example

Target: f₀ = 10 kHz, L₂ = 50 mH available

Required capacitance:

C = 1/(4π² × 10000² × 0.05) = 5.07 nF

Using parallel plates with 1.5 mm gap:

A = 5.07 × 10⁻⁹ × 0.0015 / (8.854×10⁻¹² × 80) = 10.7 cm²

Electrode size: ~3.3 cm × 3.3 cm plates (quite small!)

For more practical size, use 1 mm gap:

A = 7.1 cm² → 2.7 × 2.7 cm plates

Note: Very small WFC! May need to increase L₂ for practical electrode sizes.

Edge Effects

Real electrodes have fringing fields at edges that increase effective capacitance:

• For parallel plates, add ~0.9d to each edge dimension
• For tubes, end effects can add 5-10% to capacitance
• Guard rings can reduce edge effects in precision applications

Electrode Alignment

Critical Requirements:

Gas Evolution Considerations

When gas is produced, it affects the electrical characteristics:

• Bubbles displace water, reducing effective capacitance
• Bubble layer increases resistance
• Vertical orientation helps bubbles rise and escape
• Perforated electrodes allow better bubble release

VIC Matrix Calculator: The Water Profile section calculates WFC capacitance from your electrode dimensions. Enter geometry type, dimensions, and spacing to get accurate capacitance values for circuit design.

Next: Water Conductivity & Dielectric Properties →

Water Properties

Water Conductivity & Dielectric Properties

Water's electrical properties—conductivity and dielectric constant—directly affect WFC performance in VIC circuits. Understanding these properties helps predict circuit behavior and optimize design.

Dielectric Constant of Water

Water has an exceptionally high dielectric constant due to its polar molecular structure:

Relative Permittivity (ε_r):

Temperature Dependence:

ε_r(T) ≈ 87.74 - 0.40 × T(°C)

Why Water's ε_r is High

Water molecules are polar (have positive and negative ends). In an electric field, they align with the field, effectively multiplying the field's ability to store charge. This is why water-based capacitors have such high capacitance per unit volume.

Comparison with Other Materials

Material	εr	Relative Capacitance
Vacuum/Air	1	1× (reference)
PTFE (Teflon)	2.1	2.1×
Glass	4-10	4-10×
Ceramic	10-1000	10-1000×
Water	80	80×

Water Conductivity

Conductivity measures how easily current flows through water:

Conductivity (σ) Units:

1 S/m = 10,000 µS/cm = 10 mS/cm

Resistivity (ρ = 1/σ):

ρ (Ω·cm) = 1,000,000 / σ (µS/cm)

Conductivity of Different Waters

Water Type	σ (µS/cm)	ρ (Ω·cm)	Source
Ultra-pure (Type I)	0.055	18,000,000	Lab grade
Deionized	0.1-5	200,000-10,000,000	DI systems
Distilled	1-10	100,000-1,000,000	Distillation
Rain water	5-30	33,000-200,000	Natural
Tap water (typical)	200-800	1,250-5,000	Municipal
Well water	300-1500	670-3,300	Ground water
Sea water	50,000	20	Ocean
0.1M NaOH	~20,000	~50	Electrolyte

Calculating Solution Resistance

For Parallel Plates:

R_sol = ρ × d / A = d / (σ × A)

Example:

Effect on Q Factor

Solution resistance directly impacts circuit Q:

Q_total = 2πfL / (R_choke + R_sol + R_other)

Example Impact:

Insight: Very pure water has high Q losses! For VIC resonance, moderate conductivity may be optimal.

Frequency Dependence

Both ε_r and σ vary with frequency:

Frequency	εr Effect	σ Effect
DC - 1 MHz	Constant (~80)	Constant (DC value)
1 MHz - 1 GHz	Begins to decrease	May increase
>1 GHz	Decreases significantly	High dielectric loss

For VIC frequencies (1-100 kHz), these effects are negligible.

Temperature Effects Summary

• ε_r: Decreases ~0.4% per °C (capacitance drops as water heats)
• σ: Increases ~2% per °C (resistance drops as water heats)
• Net effect: Resonant frequency increases slightly with temperature

Measuring Water Properties

Conductivity Meters:

DIY Measurement:

VIC Matrix Calculator: Enter water conductivity in the Water Profile section. The calculator computes solution resistance and shows its impact on circuit Q. Temperature compensation is also available.

Next: Calculating WFC Capacitance →

Cell Capacitance

Calculating WFC Capacitance

Accurate calculation of WFC capacitance is essential for VIC circuit design. This page provides formulas and methods for determining the effective capacitance of various electrode configurations.

Total WFC Capacitance Model

The WFC has multiple capacitance contributions:

Series Model (simplified):

1/C_total = 1/C_edl,anode + 1/C_geo + 1/C_edl,cathode

For Practical VIC Frequencies:

At kHz frequencies, C_edl >> C_geo, so:

C_total ≈ C_geo

The geometric capacitance dominates for typical electrode gaps (>0.5 mm).

Geometric Capacitance Formulas

Parallel Plates

C = ε₀ε_rA / d

Quick Formula for Water:

C (nF) = 0.0708 × A(cm²) / d(mm)

Example:

Concentric Cylinders

C = 2πε₀ε_rL / ln(r_o/r_i)

Quick Formula for Water:

C (nF) = 4.45 × L(cm) / ln(r_o/r_i)

Thin Gap Approximation (when gap << radius):

C (nF) ≈ 0.0708 × 2πr_avg(cm) × L(cm) / d(mm)

Multiple Tubes (Array)

C_total = n × C_{single tube pair}

Where n is the number of tube pairs in parallel.

Meyer's 9-Tube Array Example:

Capacitance Calculator Table

Area (cm²)	Gap 0.5mm	Gap 1.0mm	Gap 1.5mm	Gap 2.0mm
25	3.54 nF	1.77 nF	1.18 nF	0.89 nF
50	7.08 nF	3.54 nF	2.36 nF	1.77 nF
100	14.2 nF	7.08 nF	4.72 nF	3.54 nF
200	28.3 nF	14.2 nF	9.44 nF	7.08 nF
500	70.8 nF	35.4 nF	23.6 nF	17.7 nF

Including EDL Effects

For more accurate modeling at lower frequencies or smaller gaps:

EDL Capacitance per Electrode:

C_edl = c_dl × A

Where c_dl ≈ 20-40 µF/cm² for stainless steel in water.

Total with EDL:

1/C_total = 1/C_geo + 2/C_edl

(Factor of 2 because both electrodes have EDL)

Example:

Measuring WFC Capacitance

Method 1: LCR Meter

• Most accurate method
• Measure at 1 kHz and 10 kHz (should be similar)
• Provides both C and R (ESR)
• Temperature affects reading

Method 2: RC Time Constant

1. Connect WFC in series with known resistor R
2. Apply step voltage
3. Measure time to reach 63% of final voltage
4. C = τ / R

Method 3: Resonant Frequency

1. Connect WFC with known inductor L
2. Drive with variable frequency
3. Find resonant peak
4. C = 1 / (4π²f₀²L)

Capacitance Variations

WFC capacitance can change during operation:

Factor	Effect on C	Typical Change
Temperature increase	C decreases (εr drops)	-0.4%/°C
Gas bubble formation	C decreases (less water)	-5% to -30%
Water level drop	C decreases	Proportional
Electrode coating	C may decrease	Variable
Applied voltage	Minor change	±5%

Design Workflow

1. Determine Required C

C_wfc = 1 / (4π²f₀²L₂)

2. Choose Electrode Gap

1-2 mm is typical. Smaller = higher C, larger = lower C.

3. Calculate Required Area

A = C × d / (ε₀ε_r) = C(nF) × d(mm) / 0.0708 (cm²)

4. Design Electrodes

Choose plate dimensions or tube sizes to achieve area.

5. Verify by Measurement

Build prototype and measure actual capacitance.

VIC Matrix Calculator: The Water Profile section calculates WFC capacitance automatically. Enter electrode type, dimensions, and gap. The calculator also shows how the capacitance affects resonant frequency and provides warnings if values are outside recommended ranges.

Next: Matching WFC to Circuit →

Resonant Matching

Matching WFC to Circuit

For optimal VIC performance, the WFC must be properly matched to the circuit—its capacitance must resonate with the secondary choke at the desired operating frequency. This page covers the matching process and strategies for achieving good resonance.

The Matching Problem

In a VIC circuit, we have three interdependent parameters:

f₀ = 1 / (2π√(L₂ × C_wfc))

Design Challenge:

Matching Strategies

Strategy 1: Design L₂ for Given WFC

When WFC geometry is fixed (existing cell):

L₂ = 1 / (4π²f₀²C_wfc)

Example:

Strategy 2: Design WFC for Given L₂

When using a pre-wound or available choke:

C_wfc = 1 / (4π²f₀²L₂)

Strategy 3: Tune with Additional Capacitor

When exact match isn't achievable:

If C_wfc is too low:

Add capacitor in parallel with WFC

C_total = C_wfc + C_tune

If C_wfc is too high:

Add capacitor in series with WFC (less common)

1/C_total = 1/C_wfc + 1/C_series

Impedance Matching Considerations

Beyond frequency matching, impedance levels affect energy transfer:

Secondary Characteristic Impedance:

Z₀ = √(L₂/C_wfc)

Example Comparison:

Higher Z₀ = Higher voltage for same energy

Primary-Secondary Matching

For dual-resonant VIC with both L1-C1 and L2-WFC tanks:

Configuration	Condition	Effect
Same frequency	f₀pri = f₀sec	Maximum voltage magnification
Slight offset	f₀sec ≈ 0.95-1.05 × f₀pri	Broader response, easier tuning
Harmonic	f₀sec = 2× or 3× f₀pri	Secondary resonates on harmonic

Finding Resonance

Method 1: Frequency Sweep

1. Connect oscilloscope across WFC
2. Sweep generator frequency slowly
3. Watch for voltage peak
4. Note frequency of maximum amplitude

Method 2: Phase Measurement

1. Monitor current and voltage simultaneously
2. At resonance, current and voltage are in phase (phase = 0°)
3. Below resonance: capacitive (current leads)
4. Above resonance: inductive (current lags)

Method 3: Minimum Current

For a series resonant circuit driven from a voltage source:

• Current is minimum at anti-resonance (parallel resonance)
• May need to reconfigure measurement

Troubleshooting Mismatch

Symptom	Likely Cause	Solution
No clear resonance peak	Very low Q (high losses)	Reduce water conductivity, lower DCR
Resonance far from expected	Wrong L or C values	Measure components, recalculate
Resonance drifts during operation	Temperature change, bubbles	Allow warmup, improve gas venting
Multiple resonance peaks	Coupled modes, parasitics	Check for stray coupling

Fine Tuning Tips

For L₂ Adjustment:

For C_wfc Adjustment:

For Frequency Adjustment:

Complete Matching Checklist

VIC Matrix Calculator: The Simulation tab performs complete matching analysis. Enter your choke and WFC parameters, and it calculates resonant frequency, Q factor, voltage magnification, and shows warnings if components are mismatched.

Chapter 6 Complete. Next: The VIC Matrix Calculator →

VIC Matrix Calculator

Calculator Overview

VIC Matrix Calculator Overview

The VIC Matrix Calculator is a comprehensive design tool that integrates all the concepts covered in this educational series. It allows you to design, simulate, and optimize complete VIC circuits by calculating component values, resonant frequencies, Q factors, and system behavior.

Calculator URL: https://matrix.stanslegacy.com

What the Calculator Does

The calculator brings together multiple design domains:

1. Choke Design Module

Calculate inductance, DCR, parasitic capacitance, and SRF for custom wound chokes.

2. Water Profile Module

Model the WFC as an electrical component with all relevant parameters.

3. Circuit Profile Module

Combine chokes and WFC into complete VIC circuits for analysis.

4. Simulation Module

Visualize circuit behavior and optimize performance.

Design Workflow

The recommended workflow for using the calculator:

Key Features

Feature	Description	Benefit
Real-time Calculations	Results update instantly as you change parameters	Rapid design iteration
Warning System	Alerts for out-of-range values or design issues	Avoid common mistakes
Saved Profiles	Store and recall choke, water, and circuit configurations	Compare designs easily
Interconnected Models	Changes propagate through entire system	See full system impact
Educational Notes	Tooltips and explanations throughout	Learn while designing

Input vs. Output Parameters

You Provide (Inputs):

Calculator Provides (Outputs):

Accuracy and Limitations

Parameter	Typical Accuracy	Notes
Inductance	±10-20%	Core properties vary; always verify
DCR	±5%	Depends on wire tables accuracy
WFC Capacitance	±15%	Fringe effects, water purity affect results
Q Factor	±20-30%	Multiple loss mechanisms; use as estimate
Resonant Frequency	±10-15%	Depends on L and C accuracy

Important: The calculator provides design estimates. Always verify critical parameters with measurements on actual components. Real-world results may vary due to manufacturing tolerances, stray inductance/capacitance, and environmental factors.

Getting Started

To begin using the VIC Matrix Calculator:

1. Navigate to the application dashboard
2. Start with the module that matches your first design decision:
   • If you have specific chokes → Start with Choke Design
   • If you have a specific WFC → Start with Water Profile
   • If you have target frequency → Work backwards from Circuit Profile
3. Follow the guided workflow to complete your design

Tip: The following pages in this chapter provide detailed guidance on each module. Work through them in order for the best understanding of the calculator's capabilities.

Next: Component Input Parameters →

Component Inputs

Component Input Parameters

This page details all input parameters used across the VIC Matrix Calculator modules. Understanding what each parameter means and how to determine its value is essential for accurate calculations.

Choke Design Inputs

Core Parameters

Parameter	Symbol	Units	Description
Core Type	—	—	Toroid, E-core, rod, bobbin, or air-core
Core Material	—	—	Ferrite mix, iron powder, or air
Relative Permeability	μᵣ	—	Material permeability (1 for air, 2000+ for ferrite)
AL Value	Aₗ	nH/turn²	Inductance factor (from core datasheet)
Outer Diameter	OD	mm	Core outer diameter (toroids)
Inner Diameter	ID	mm	Core inner diameter (toroids)
Height	H	mm	Core height/thickness

Finding Core Parameters:

Wire Parameters

Parameter	Symbol	Units	Description
Wire Gauge	AWG	AWG	American Wire Gauge number
Wire Material	—	—	Copper, aluminum, silver
Number of Turns	N	turns	Total turns wound on core
Number of Layers	nlayers	—	Winding layers (affects parasitic C)
Winding Style	—	—	Single, bifilar, or multi-filar

Bifilar-Specific Parameters

Parameter	Description
Choke Role	Primary (L1), Secondary (L2), or Bifilar Set
Coupling Coefficient	k value between bifilar windings (typically 0.95-0.99)
Inter-winding Insulation	Thickness and material of insulation between wires

Water Profile Inputs

Electrode Geometry

Parameter	Symbol	Units	Description
Electrode Type	—	—	Parallel plates, concentric tubes, tube array
Electrode Area	A	cm²	Active electrode surface area
Electrode Gap	d	mm	Distance between electrodes
Inner Radius	ri	mm	Inner tube radius (cylindrical)
Outer Radius	ro	mm	Outer tube radius (cylindrical)
Tube Length	L	cm	Submerged tube length
Number of Tubes	n	—	Tube pairs in array

Water Properties

Parameter	Symbol	Units	Description
Water Conductivity	σ	µS/cm	Electrical conductivity of water
Water Temperature	T	°C	Operating temperature
Dielectric Constant	εr	—	Relative permittivity (~80 for water at 20°C)

Measuring Conductivity:

Circuit Profile Inputs

Component Selection

Parameter	Description
Primary Choke (L1)	Select from saved choke designs or enter values
Secondary Choke (L2)	Select from saved choke designs or enter values
Water Profile (WFC)	Select from saved water profiles or enter values
Primary Capacitor (C1)	Capacitance value for primary resonance
Tuning Capacitor	Optional capacitor in parallel with WFC

Operating Parameters

Parameter	Symbol	Units	Description
Operating Frequency	fop	kHz	Pulse generator frequency
Input Voltage	Vin	V	Peak pulse voltage
Duty Cycle	D	%	Pulse on-time percentage
Source Resistance	Rs	Ω	Driver output impedance

Direct Value Entry

If you have measured values for components (rather than designing from scratch), you can enter them directly:

For Chokes:

For WFC:

Best Practice: When possible, measure actual component values and compare to calculated values. This helps identify measurement errors and improves your understanding of the calculator's accuracy for your specific components.

Next: Simulation Tab Explained →

Simulation Tab

Simulation Tab Explained

The Simulation tab provides visual analysis of your VIC circuit design. It generates frequency response curves, time-domain waveforms, and key performance metrics that help you understand and optimize circuit behavior.

Simulation Overview

The simulation performs several types of analysis:

1. Frequency Domain Analysis

Sweeps through a frequency range to show how the circuit responds at different frequencies.

2. Impedance Analysis

Shows how circuit impedance varies with frequency, identifying resonant points.

3. Time Domain Analysis

Simulates actual voltage and current waveforms during pulse operation.

4. Ring-down Analysis

Shows how oscillations decay after excitation stops.

Frequency Response Display

The frequency response plot shows amplitude vs. frequency:

Amplitude
    ↑
    │
    │              ╱╲
    │             ╱  ╲          ← Secondary resonance
    │            ╱    ╲
    │           ╱      ╲
    │   ╱╲     ╱        ╲
    │  ╱  ╲   ╱          ╲
    │ ╱    ╲ ╱            ╲
    │╱      ╳              ╲
    └─────────────────────────→ Frequency (kHz)
         ↑           ↑
    Primary      Secondary
    resonance    resonance

Key Features in Plot

Feature	What It Means	Ideal Characteristic
Peak Height	Voltage magnification at resonance	Higher = more voltage gain
Peak Sharpness	Q factor (sharp = high Q)	Depends on application
Peak Location	Resonant frequency f₀	Should match design target
-3dB Bandwidth	Frequency range at 70.7% of peak	Narrower = higher Q
Multiple Peaks	Primary and secondary resonances	Aligned for max transfer

Calculated Metrics

The simulation calculates and displays these key values:

Resonance Parameters

Q Factor Metrics

Performance Metrics

Impedance Plot

Shows circuit impedance magnitude and phase vs. frequency:

|Z| (Ω)                          Phase
   ↑                               ↑
   │      ╱╲                       │         ╱────
   │     ╱  ╲    ← Peak at        │        ╱
   │    ╱    ╲     resonance      │       ╱
   │   ╱      ╲                    │──────╳  ← 0° at f₀
   │  ╱        ╲                   │     ╱
   │ ╱          ╲                  │    ╱
   │╱            ╲                 │───╱────
   └──────────────────→ f         └──────────────→ f

Interpreting Impedance

• Peak impedance: Maximum at parallel resonance
• Minimum impedance: At series resonance points
• Phase = 0°: Indicates resonant frequency
• Positive phase: Inductive behavior (current lags)
• Negative phase: Capacitive behavior (current leads)

Time Domain Waveforms

The time-domain view shows actual voltage and current over time:

Waveforms Displayed:

What to Look For:

Ring-Down Display

Shows oscillation decay after excitation stops:

Voltage
   ↑
   │╱╲
   │  ╲╱╲
   │    ╲╱╲
   │      ╲╱╲
   │        ╲╱╲
   │          ╲╱╲
   │            ╲╱─── → Envelope decay
   │              ╲
   └────────────────────→ Time

   ←─── τ ───→
   (63% decay)

Ring-Down Metrics

Metric	Formula	Significance
Time Constant (τ)	τ = 2L/R	Time to decay to 37%
Ring-down Cycles	n ≈ 0.733 × Q	Oscillations before decay
Settling Time	~5τ for 99% decay	Time to reach steady state

Warning Indicators

The simulation flags potential issues:

Warning	Meaning	Action
⚠️ Near SRF	Operating frequency close to choke SRF	Reduce frequency or redesign choke
⚠️ Low Q	Q factor below recommended threshold	Reduce losses (DCR, water R)
⚠️ Frequency Mismatch	Primary and secondary not aligned	Adjust C1 or component values
⚠️ High Voltage	Magnified voltage exceeds safe limits	Verify insulation ratings

Using Simulation Results

Design Iteration Process:

Pro Tip: Save your circuit profile before making changes. This allows you to compare different configurations side-by-side and roll back if needed.

Next: Circuit Optimization Strategies →

Optimization

Circuit Optimization Strategies

This page covers practical strategies for optimizing your VIC circuit design using the calculator. Learn how to achieve specific goals like maximizing Q, hitting a target frequency, or optimizing voltage magnification.

Optimization Goals

Different applications may prioritize different characteristics:

Goal	Optimize For	Trade-offs
Maximum Voltage	High Q, matched resonance	Narrower bandwidth, critical tuning
Stable Operation	Moderate Q, wide bandwidth	Lower peak voltage
Frequency Flexibility	Lower Q, broader response	Reduced magnification
Energy Efficiency	Minimize losses (DCR, Rsol)	May require larger components

Strategy 1: Maximizing Q Factor

Q determines voltage magnification and selectivity. To maximize Q:

Reduce Choke DCR:

Reduce Solution Resistance:

Increase L or Decrease C:

Q Factor Relationships:

Q = 2πf₀L/R = Z₀/R = √(L/C)/R

To double Q: halve R, or quadruple L (while quartering C to maintain f₀)

Strategy 2: Hitting Target Frequency

When you need a specific resonant frequency:

Approach A: Fixed L, Adjust C

Approach B: Fixed C, Adjust L

Approach C: Adjust Both

Fine-Tuning Frequency

Adjustment	Effect on f₀	Typical Range
Add parallel capacitor	Decreases f₀	1-50 nF typical
Adjust core gap (if gapped)	Changes L → changes f₀	±20% L adjustment
Add/remove turns	Changes L significantly	L ∝ N²
Change water level	Changes C → changes f₀	Proportional to area

Strategy 3: Matching Primary to Secondary

For maximum energy transfer, align primary and secondary resonances:

Exact Match (f₀_pri = f₀_sec):

Slight Offset (5-10% difference):

Calculator Approach:

Strategy 4: Optimizing for Available Components

When working with existing components:

Step 1: Characterize What You Have

Step 2: Calculate Natural Resonance

f₀ = 1/(2π√LC)

This is where your circuit wants to resonate.

Step 3: Evaluate Performance

Step 4: Adjust as Needed

Sensitivity Analysis

Understanding how sensitive your design is to variations:

Parameter Change	Effect on f₀	Effect on Q
L +10%	f₀ -5%	Q +5%
C +10%	f₀ -5%	Q -5%
R +10%	No change	Q -10%
Temperature +10°C	f₀ +2% (due to εr drop)	Q +5% (Rsol drops)

Common Optimization Mistakes

❌ Chasing Extreme Q

Very high Q makes the circuit sensitive to drift and hard to tune. Q of 50-100 is often more practical than Q > 200.

❌ Ignoring SRF

A design that works on paper fails if operating frequency is too close to SRF. Always check this!

❌ Forgetting Water Resistance

Solution resistance often dominates losses. Pure distilled water has higher resistance than you might expect.

❌ Not Accounting for Parasitics

Real circuits have stray inductance and capacitance. Leave margin for these effects.

❌ Over-constraining the Design

If you fix too many parameters, you may have no degrees of freedom for optimization.

Optimization Checklist

Remember: Optimization is iterative. The calculator makes it easy to try variations quickly. Don't expect to find the optimal design on the first try—explore the design space!

Next: Interpreting Calculation Results →

Interpreting Results

Interpreting Calculation Results

Understanding what the calculator's output values mean and how to use them for practical circuit construction. This page helps you translate numbers into actionable design decisions.

Understanding Output Values

Inductance Results

Output	Typical Range	What It Means
L (inductance)	1-100 mH	Primary choke property, affects f₀ and Q
DCR	0.1-50 Ω	Wire resistance, major Q limiter
SRF	50 kHz - 1 MHz	Maximum usable frequency
Cparasitic	10-500 pF	Stray capacitance, determines SRF
Wire Length	1-50 m	Total wire needed for winding

Capacitance Results

Output	Typical Range	What It Means
Cwfc	1-100 nF	WFC capacitance, sets resonance with L
Rsolution	0.1-100 Ω	Water resistance, affects Q
Z₀ (characteristic)	100-10,000 Ω	√(L/C), impedance at resonance

Circuit Results

Output	Typical Range	Interpretation
f₀ (resonant freq)	1-100 kHz	Where circuit resonates naturally
Q factor	5-200	Resonance sharpness, voltage gain
Bandwidth	50 Hz - 5 kHz	Usable frequency range around f₀
Vmagnification	5× - 200×	Voltage gain at resonance
Ring-down τ	0.1-10 ms	Decay time constant
Ring-down cycles	3-150	Oscillations during decay

What "Good" Values Look Like

✓ Well-Designed VIC Circuit:

✗ Warning Signs:

Translating Results to Construction

Wire Length and Turns

The calculator provides wire length and turn count. When winding:

Component Selection

Calculated Value	Selection Guidance
C1 = 47.3 nF	Use 47 nF standard value (within 1%)
C1 = 31.2 nF	Use 33 nF or parallel 22+10 nF
L = 15.7 mH	Wind for 16 mH, fine-tune with parallel C

Understanding Accuracy Limits

Know what to expect from calculated vs. measured values:

Parameter	Expected Accuracy	Why Variation Occurs
Inductance	±10-20%	Core μᵣ varies, winding geometry imperfect
DCR	±5%	Wire tables accurate, but length varies
SRF	±30%	Parasitic C is hard to model precisely
Cwfc	±15%	Fringe effects, water purity variation
Rsolution	±20%	Conductivity varies with temperature
f₀ (calculated)	±15%	Depends on L and C accuracy
Q factor	±25%	Multiple loss mechanisms combine

Comparing Calculated vs. Measured

When Measured f₀ is Lower Than Calculated:

When Measured f₀ is Higher Than Calculated:

When Measured Q is Lower Than Calculated:

Using Results for Troubleshooting

Observation	Calculator Check	Likely Issue
No resonance found	Check SRF vs. operating frequency	Operating above SRF
Very weak resonance	Check calculated Q	High losses, low Q
Resonance at wrong frequency	Verify L and C inputs	Input error or mismeasurement
Less voltage gain than expected	Compare Q values	Actual losses higher
Resonance drifts during use	Check temperature effects	Water heating, capacitance changing

Results Summary Checklist

Before building, verify these from your results:

Final Advice: The calculator gives you an excellent starting point. Always plan to measure your actual circuit and iterate. The goal is to get close enough that minor tuning (adjusting C1, trimming frequency) achieves optimal performance.

Chapter 7 Complete. Next: Advanced Topics →

VIC Matrix Calculator Application

The VIC Matrix Calculator (v6) can be found at the following url:

https://matrix.stanslegacy.com

 

Advanced Topics

PLL Control

PLL-Based Frequency Control

Phase-Locked Loop (PLL) circuits can automatically track and maintain resonance in VIC systems, compensating for drift due to temperature changes, water level variations, and other factors. This page covers PLL fundamentals and their application to VIC circuits.

Why PLL Control?

VIC resonant frequency can drift during operation due to:

Factor	Effect on f₀	Typical Drift
Water temperature rise	f₀ increases (εr drops)	+0.2%/°C
Gas bubble formation	f₀ increases (C drops)	+2-10%
Water level change	f₀ changes (C changes)	Variable
Core temperature rise	f₀ may shift (μ changes)	±1%

A PLL can continuously adjust the drive frequency to maintain optimal resonance despite these variations.

PLL Fundamentals

Basic PLL Components:

Reference ──→ [Phase      ] ──→ [Loop    ] ──→ [VCO     ] ──→ Output
Signal        [Detector   ]     [Filter  ]     [        ]     Frequency
                   ↑                                │
                   └────────────────────────────────┘
                            Feedback
    

Components Explained:

PLL for VIC Resonance Tracking

For VIC applications, the PLL tracks the resonant frequency by sensing the phase relationship between drive signal and cell response:

           ┌──────────────────────────────────────┐
           │                                      │
Drive ──→ [VIC Circuit] ──→ Vwfc ──→ [Phase    ] ──→ [Loop   ] ──→ [VCO]
Signal                              [Detector ]     [Filter ]         │
  ↑                                      ↑                           │
  └──────────────────────────────────────┴───────────────────────────┘
                              Feedback Loop

Phase Detection Methods

Method	Description	Pros/Cons
XOR Phase Detector	Digital XOR of drive and response	Simple, but needs square waves
Analog Multiplier	Multiply drive × response	Works with sinusoids, more complex
Zero-Crossing Detector	Compare zero-crossing times	Digital-friendly, noise sensitive
I/Q Demodulation	Quadrature phase detection	Most accurate, most complex

Resonance Tracking Logic

At resonance, the phase relationship between drive current and WFC voltage is 0°:

Phase vs. Frequency:

Control Law:

Loop Filter Design

The loop filter determines how quickly the PLL responds to changes:

Parameter	Fast Response	Slow Response
Tracking speed	Quick adaptation	Slow adaptation
Noise rejection	Poor	Good
Stability	May oscillate	More stable
Best for	Rapid changes	Gradual drift

Design Tip: For VIC applications, a medium-speed loop (bandwidth ~100-500 Hz) usually works well. Fast enough to track bubble-induced changes, slow enough to reject noise.

VCO Implementation

The VCO generates the variable-frequency drive signal:

Common VCO Options:

VCO Requirements:

Complete PLL-VIC System

                    PLL CONTROLLER
     ┌────────────────────────────────────────┐
     │                                        │
     │  [Phase Det] ──→ [Loop Filter] ──→ Vctrl
     │       ↑                           │    │
     │       │                           │    │
     └───────┼───────────────────────────┼────┘
             │                           │
             │                           ↓
     Vsense  │                        [VCO]
       ↑     │                           │
       │     │                           ↓
       │     │                     [Driver Stage]
       │     │                           │
       │     │      ┌────────────────────┘
       │     │      ↓
       │     └── [L1] ──── [C1] ──────────┐
       │                                  │
       │         ┌────────────────────────┘
       │         │
       │         ↓
       └──── [L2] ──── [WFC]
                    ↑
              Resonating
               Circuit

Practical Considerations

Startup Sequence:

Lock Detection:

Monitor loop filter output—stable voltage indicates lock. Large variations indicate searching or loss of lock.

Capture Range:

PLL can only lock if initial frequency is within "capture range." If f₀ drifts too far, may need frequency sweep to re-acquire.

Alternatives to PLL

Method	Description	When to Use
Fixed Frequency	No tracking, fixed drive	Stable systems, low Q
Frequency Sweep	Periodically sweep through range	Testing, characterization
Peak Detector	Track amplitude maximum	Simpler than phase tracking
Self-Oscillation	Circuit sets own frequency	Simple, but less control

VIC Matrix Calculator Note: The VIC5 PLL module provides calculations for PLL component selection, including VCO tuning range, loop filter values, and expected tracking bandwidth. Use these calculations when implementing automatic resonance tracking.

Next: Harmonic Analysis →

Harmonic Analysis

Harmonic Analysis

VIC circuits are typically driven by non-sinusoidal waveforms (pulses, square waves), which contain harmonics. Understanding how these harmonics interact with the resonant circuit is important for predicting actual performance and potential interference effects.

Fourier Analysis Basics

Any periodic waveform can be decomposed into a sum of sinusoids:

Fourier Series:

f(t) = a₀ + Σ[aₙcos(nωt) + bₙsin(nωt)]

Where n = 1, 2, 3... are the harmonic numbers (n=1 is fundamental)

Harmonic Content of Common Waveforms

Square Wave

50% duty cycle square wave contains only odd harmonics:

V(t) = (4V_pk/π)[sin(ωt) + (1/3)sin(3ωt) + (1/5)sin(5ωt) + ...]

Pulse Train (Variable Duty Cycle)

Pulse train with duty cycle D contains both odd and even harmonics:

a_n = (2V_pk/nπ) × sin(nπD)

Effect of Duty Cycle:

Resonant Circuit Response to Harmonics

A resonant circuit acts as a bandpass filter. It responds most strongly to frequencies near f₀:

Response
    │
    │          Fundamental
    │              ↓
    │             ╱╲
    │            ╱  ╲              3rd harmonic
    │           ╱    ╲                 ↓
    │          ╱      ╲              (small response)
    │         ╱        ╲           ┌─┐
    │        ╱          ╲          │ │
    └───────────────────────────────────────→ f
             f₀          3f₀

Response at Harmonic Frequencies:

H(nf) = 1 / √[1 + Q²(n - 1/n)²]

For high Q circuits, harmonics far from f₀ are strongly attenuated.

Example (Q=50, f₀=10 kHz):

Harmonic Resonance

If a harmonic happens to fall near f₀, it can cause problems or opportunities:

Scenario	Effect	Action
Drive at f₀	Fundamental resonates	Normal operation
Drive at f₀/2	2nd harmonic resonates	May be useful or problematic
Drive at f₀/3	3rd harmonic resonates	Subharmonic driving
Harmonic hits SRF	Choke self-resonates	Avoid—causes problems

Sub-Harmonic Driving

It's possible to drive the circuit at a sub-multiple of f₀ and let a harmonic excite resonance:

Example: 3rd Harmonic Drive

Advantages:

Disadvantages:

Pulse Shaping for Harmonic Control

Adjusting pulse shape can control harmonic content:

Technique	Effect
Slower edges (rise/fall time)	Reduces high-order harmonics
Duty cycle = 1/n	Eliminates nth harmonic
Trapezoidal waveform	Controlled harmonic rolloff
Sine wave drive	No harmonics (pure fundamental)

Harmonic Interaction with Multiple Resonances

In dual-resonant VIC (primary + secondary), harmonics may interact with both:

Response
    │
    │     Primary        Secondary
    │     resonance      resonance
    │         ↓              ↓
    │        ╱╲            ╱╲
    │       ╱  ╲          ╱  ╲
    │      ╱    ╲        ╱    ╲
    │     ╱      ╲      ╱      ╲
    │    ╱        ╲    ╱        ╲
    │   ╱          ╲  ╱          ╲
    └──────────────────────────────────→ f
          f₀,pri       f₀,sec

If f₀,sec = 3 × f₀,pri, then:

• Fundamental drives primary resonance
• 3rd harmonic drives secondary resonance
• This is sometimes called "harmonic matching"

Practical Harmonic Considerations

EMI Concerns:

Harmonics can cause electromagnetic interference. Shield appropriately and consider filtering if needed.

Measurement:

Use an oscilloscope with FFT function or spectrum analyzer to view harmonic content of your signals.

Design Rule:

For clean resonance, ensure no significant harmonics fall within the passband (f₀ ± f₀/Q) of unintended resonances.

Harmonic Analysis in VIC Matrix Calculator

Calculator Feature: The simulation can show frequency response across a range that includes harmonics. When analyzing a design, check whether any harmonics of your drive frequency fall near the circuit's resonant points or SRF values.

Next: Transformer Coupling Effects →

Transformer Coupling

Transformer Coupling Effects

In VIC circuits, the primary (L1) and secondary (L2) chokes may be magnetically coupled, either intentionally (bifilar winding) or unintentionally (proximity). This coupling significantly affects circuit behavior and must be understood for accurate analysis.

Magnetic Coupling Fundamentals

When two inductors share magnetic flux, they become coupled:

Mutual Inductance:

M = k × √(L₁ × L₂)

Where k is the coupling coefficient (0 ≤ k ≤ 1)

Coupling Coefficient:

Coupled Inductor Equivalent Circuit

Coupled inductors can be modeled as a transformer with leakage inductances:

    Ideal Coupled Inductors:          Equivalent T-Model:

         L₁          L₂                  L₁(1-k)    L₂(1-k)
    ○────UUUU────●────UUUU────○      ○────UUUU──●──UUUU────○
                 │                              │
              M (mutual)                    k√(L₁L₂)
                                               │
                                              ─┴─

T-Model Components

Component	Formula	Represents
Lleak1	L₁(1-k)	Primary leakage inductance
Lleak2	L₂(1-k)	Secondary leakage inductance
Lm	k√(L₁L₂)	Magnetizing inductance

Effect on VIC Circuit Behavior

Resonant Frequency Shifts

Coupling changes the effective inductances seen by each resonant tank:

Without Coupling (k=0):

f₀,pri = 1/(2π√(L₁C₁))
f₀,sec = 1/(2π√(L₂C_wfc))

With Coupling:

The system has two coupled resonant modes. The frequencies split into:

f₁, f₂ = function of L₁, L₂, C₁, C_wfc, and k

Exact formulas are complex—use simulation for accurate prediction.

Mode Splitting

Coupled resonators exhibit "mode splitting"—two distinct resonant frequencies instead of one:

    Uncoupled (k=0):              Coupled (k>0):

    Response                      Response
        │                             │
        │     ╱╲                      │   ╱╲    ╱╲
        │    ╱  ╲                     │  ╱  ╲  ╱  ╲
        │   ╱    ╲                    │ ╱    ╲╱    ╲
        └────────────→ f              └──────────────→ f
             f₀                          f₁    f₂

    Single resonance            Split into two modes

Mode Splitting (equal resonators):

When f₀,pri = f₀,sec = f₀:

f₁ ≈ f₀ / √(1+k) (lower mode)
f₂ ≈ f₀ / √(1-k) (upper mode)

Separation increases with coupling coefficient k.

Energy Transfer

Coupling provides a path for energy transfer between primary and secondary:

Coupling	Energy Transfer	VIC Behavior
k = 0 (none)	Only through shared current path	Independent resonances
k = 0.1-0.3	Moderate magnetic coupling	Slight interaction
k = 0.5-0.8	Strong coupling	Significant mode splitting
k > 0.9	Very tight coupling	Behaves more like transformer

Bifilar Winding Coupling

Bifilar chokes have inherently high coupling (k ≈ 0.95-0.99):

Effects of Bifilar Coupling:

Measuring Bifilar Coupling:

Stray Coupling

Even separate chokes may have unintended coupling if placed close together:

Configuration	Typical k	Mitigation
Toroids touching	0.01-0.05	Separate by >2× diameter
Air-core coils aligned	0.1-0.3	Orient perpendicular
Coils on same rod	0.5-0.9	Use separate cores

Design Considerations

When to Use Coupling:

When to Avoid Coupling:

Layout Guidelines:

Analyzing Coupled VIC Circuits

Coupled Circuit Analysis Steps:

Simulation Tip: When k > 0.1, coupled effects become significant. Always include coupling in simulation if windings share a core or are in close proximity.

VIC Matrix Calculator: The Choke Design module includes coupling coefficient input for bifilar windings. The simulation accounts for mutual inductance effects when analyzing coupled systems.

Next: Energy Efficiency Analysis →

Energy Efficiency

Energy Efficiency Analysis

Understanding energy flow in VIC circuits helps optimize performance and evaluate system efficiency. This page covers how to analyze energy storage, transfer, and dissipation in resonant VIC systems.

Energy in Resonant Circuits

In an LC resonant circuit, energy oscillates between the inductor and capacitor:

Energy Storage:

E_L = ½LI² (energy in inductor)
E_C = ½CV² (energy in capacitor)

At Resonance:

E_total = E_L,max = E_C,max = ½CV_peak²

Peak Energy (example):

Energy Flow Diagram

                    Input Power
                         │
                         ↓
    ┌─────────────────────────────────────────────┐
    │              VIC CIRCUIT                     │
    │                                              │
    │  ┌──────┐      ┌──────┐      ┌──────┐       │
    │  │ L1   │──────│ L2   │──────│ WFC  │       │
    │  │ DCR  │      │ DCR  │      │ ESR  │       │
    │  └──────┘      └──────┘      └──────┘       │
    │      │             │             │          │
    │      ↓             ↓             ↓          │
    │  Heat Loss    Heat Loss     Heat Loss       │
    │  (copper)     (copper)      (solution)      │
    │                                  │          │
    │                                  ↓          │
    │                           Electrochemical   │
    │                           Work (desired)    │
    └─────────────────────────────────────────────┘

Loss Mechanisms

Loss Type	Formula	How to Minimize
Choke DCR Loss	P = I²RDCR	Use larger wire, copper
Solution Resistance	P = I²Rsol	Optimize water conductivity
Core Loss	P ∝ f^α × B^β	Choose low-loss core material
Skin Effect Loss	Increases R at high f	Use Litz wire at high f
Dielectric Loss	P = ωCV² × tan(δ)	Use low-loss capacitors

Q Factor and Efficiency

Q factor is directly related to energy efficiency per cycle:

Energy Loss Per Cycle:

ΔE_cycle = 2π × E_stored / Q

Interpretation:

Energy Retention:

After n cycles: E(n) = E₀ × e^(-2πn/Q)

Power Flow Analysis

Input Power

P_in = V_in × I_in × cos(φ)

For pulsed operation:

P_avg = (1/T) × ∫V(t)I(t)dt

Dissipated Power

P_diss = I_rms² × R_total

Where R_total = R_DCR1 + R_DCR2 + R_sol + R_other

Useful Power

Power available for electrochemical work:

P_useful = P_in - P_diss

Or, for the WFC specifically:

P_wfc = V_wfc × I_wfc × cos(φ_wfc)

Efficiency Calculations

Efficiency Type	Formula	Typical Values
Resonant Tank η	η = Q/(Q+1) ≈ 1 - 1/Q	90-99% for high Q
Power Transfer η	η = Pwfc/Pin	50-90%
Voltage Multiplication η	Vout/Vin (at resonance)	10-100× typical

Energy Balance Verification

To verify your analysis is correct, energy must balance:

Steady State:

P_in = P_DCR1 + P_DCR2 + P_sol + P_core + P_other

Check:

Loss Breakdown Example

Component	Resistance	Power Loss (at 1A)	% of Total
L1 DCR	2.5 Ω	2.5 W	25%
L2 DCR	3.0 Ω	3.0 W	30%
Rsolution	4.0 Ω	4.0 W	40%
Other (core, leads)	0.5 Ω	0.5 W	5%
Total	10 Ω	10 W	100%

Improving Efficiency

High-Impact Improvements:

Diminishing Returns:

Once a loss mechanism is <10% of total, further improvement has limited benefit. Focus on the dominant losses.

Thermal Considerations

All dissipated power becomes heat:

Component	Heat Concern	Mitigation
Choke windings	Wire insulation damage	Adequate wire size, ventilation
Ferrite core	Curie temp, permeability change	Keep below rated temperature
Water/WFC	Boiling, capacitance drift	Monitor temperature, allow cooling
Capacitors	ESR heating, life reduction	Use low-ESR types, derate

VIC Matrix Calculator: The simulation module calculates expected power dissipation in each component. Use this to identify thermal hotspots and verify your design won't overheat during operation.

Next: Experimental Validation Methods →

Experimental Validation

Experimental Validation Methods

Theoretical calculations and simulations must be validated with actual measurements. This page covers practical techniques for measuring VIC circuit parameters and comparing results to predictions.

Essential Test Equipment

Equipment	Purpose	Key Specifications
Oscilloscope	Waveform viewing, frequency measurement	2+ channels, 100+ MHz bandwidth
Function Generator	Provide test signals	1 Hz - 1 MHz, variable duty cycle
LCR Meter	Measure L, C, R	Multiple test frequencies (1 kHz, 10 kHz)
Multimeter	DC resistance, voltage	True RMS, low-ohm capability
Current Probe	Non-contact current measurement	AC/DC, appropriate bandwidth
High-Voltage Probe	Measure high voltages safely	1000:1 or 100:1, rated voltage

Component Verification

Measuring Inductance

Method 1: LCR Meter (Preferred)

Method 2: Resonance with Known C

Measuring DCR

Four-Wire (Kelvin) Measurement:

For accurate low-resistance measurement, use 4-wire method to eliminate lead resistance:

Expected accuracy: ±1-5% compared to calculated value

Measuring WFC Capacitance

Expected accuracy: ±10-20% compared to calculated value

Resonant Frequency Measurement

Frequency Sweep Method

Setup:

Function ──→ [VIC    ] ──→ Oscilloscope
Generator    [Circuit]     Ch1: Input
                          Ch2: Output (across WFC)
    

Procedure:

Calculate Q from Measurement:

Q = f₀ / (f_high - f_low) = f₀ / BW

Phase Measurement Method

Q Factor Measurement

Method 1: Bandwidth

Measure -3dB bandwidth and calculate:

Q = f₀ / BW

Method 2: Ring-Down

Alternatively, measure time constant τ:

τ = 2L/R = Q/(πf₀)

Method 3: Voltage Magnification

Caution: This assumes lossless input coupling. Actual Q may be higher due to source impedance effects.

Comparing Calculated vs. Measured

Parameter	Acceptable Difference	If Larger Difference
Inductance	±20%	Check core μᵣ, turn count
DCR	±10%	Check wire gauge, connections
WFC Capacitance	±20%	Check geometry, water level
Resonant Frequency	±15%	Check L and C values
Q Factor	±30%	Look for missing losses

Troubleshooting Discrepancies

Measured f₀ Lower than Calculated:

Measured f₀ Higher than Calculated:

Measured Q Lower than Calculated:

No Clear Resonance Observed:

Documentation Template

Record for Each Test:

Date: ___________
Circuit ID: ___________

COMPONENT VALUES (Calculated / Measured):
L1: _______ mH / _______ mH
L2: _______ mH / _______ mH
DCR1: _______ Ω / _______ Ω
DCR2: _______ Ω / _______ Ω
C_wfc: _______ nF / _______ nF
C1: _______ nF / _______ nF

RESONANCE (Calculated / Measured):
f₀_primary: _______ kHz / _______ kHz
f₀_secondary: _______ kHz / _______ kHz

PERFORMANCE (Calculated / Measured):
Q: _______ / _______
Bandwidth: _______ Hz / _______ Hz
V_magnification: _______ / _______

NOTES:
_________________________________
    

Safety Considerations

⚠️ High Voltage Warning:

⚠️ Gas Production:

Best Practice: Always compare measured values to calculator predictions. This builds confidence in both your construction skills and the calculator's accuracy. Document discrepancies—they often reveal important lessons about real-world effects.

Chapter 8 Complete. See Appendices for reference tables and formulas. →

Understanding Resonant Action in the Water Fuel Cell

This article explains the principle of Resonant Action — the mechanism by which Stan Meyer's Water Fuel Cell achieves water dissociation through matched mechanical and electrical resonance, rather than brute-force electrolysis. We walk through the physics, the patent language, and the math to arrive at a complete, actionable design chain.

---

Why Water's Dielectric Properties Matter

The Voltage Intensifier Circuit (VIC) operates in the 1 kHz – 100 kHz range, where both dipolar and ionic mechanisms in water are fully active. At these frequencies, water's dielectric constant remains very high (~78–80), making it an excellent capacitor dielectric inside the gas processor tubes.

The dipolar relaxation cutoff for water doesn't occur until ~17–20 GHz — far above VIC operating range. This means at our target frequencies, water molecules can physically respond to the applied electric field. This is the basis of Stan's Electrical Polarization Process (EPP).

Patents #5,149,407 and WO8912704A1 describe this explicitly:

"Water molecules are broken down into hydrogen and oxygen gas atoms in a capacitive cell by a polarization and resonance process dependent upon the dielectric properties of water."

Complex Permittivity

Water's permittivity has two components that matter for VIC design:

• Real part (ε') — determines the cell's capacitance and therefore your resonant frequency
• Imaginary part (ε'') — the loss tangent, which directly reduces your circuit's Q factor

Because permittivity changes with temperature, conductivity, and frequency, your water "capacitor" is a moving target. This is why VIC tuning can drift during operation, and why water purity matters — too many dissolved ions dump current into conductance instead of polarization.

---

The Ionization-Conductivity Feedback Loop

Applying voltage to water creates a chain reaction:

1. Voltage ionizes the molecule → creates H⁺ and OH⁻ carriers
2. Conductivity goes up → loss tangent (ε'') rises → Q factor drops
3. Resonance degrades

This is precisely why the VIC uses pulsed voltage rather than continuous DC. Hit the molecule hard and fast, then let it rest. The rest period allows electrical polarization to weaken the covalent bond before excessive ionization destroys the resonant condition.

Apply continuous voltage and conductivity keeps climbing — the cell stops acting like a capacitor and starts acting like a resistor. You've built an expensive water heater, not a fuel cell.

Per Patent #4,936,961, the key is that electrical polarization weakens the covalent bond before full ionization occurs. The WFC operates in the narrow window between polarization and brute-force electrolysis.

---

Corrugated Geometry: Momentary Entrapment

Corrugated cell surfaces serve a dual purpose that goes beyond simple surface area increase:

• Peak of corrugation → intense local electric field → strong EPP → bonds weakened at focal points
• Bulk water between peaks → lower average field → lower ionization → conductivity stays manageable

This gives you localized electrical polarization without destroying the Q factor in the bulk medium. You can run higher effective field gradients than smooth tubes at the same voltage, before conductivity kills your resonance.

Patent EP0103656A2 — Resonant Cavity for Hydrogen Generator

Filed December 14, 1982, this is one of Stan's earliest European filings. The patent text on the corrugated exciter (Figure 6) is explicit about why corrugations matter:

"Instead of a forward direct line back-and-forth path of the atom flow, the corrugations of the convex 47 and concave 49 surfaces causes the atoms to move in forward and backward / back-and-forth path."

"The increased surface area provided by the corrugations and creating the resonant cavity, thus enhances the sub-atomic action."

The corrugations aren't just field concentrators — they force molecules into an oscillatory path, increasing residence time in the high-gradient zone. This is Momentary Entrapment to assist Resonant Action: the geometry traps the molecule long enough for multiple resonant cycles to act on it, rather than letting it blow straight through the gap in a single cycle.

A water molecule at room temperature moves at roughly 600 m/s thermally. In a 1 cm gap, it transits in about 16 microseconds — barely one cycle at 60 kHz. The corrugation multiplies the effective interaction time by 5–10x, turning a single glancing pass into meaningful resonant coupling.

---

The Key Insight: Cavity Spacing = Wavelength

The critical passage comes from Patent #4,798,661 (Gas Generator Voltage Control Circuit):

"The phenomena that the spacing between two objects is related to the wavelength of a physical motion between the two objects is utilized herein."

"The pulsing voltage on the plate exciters applying a physical force is matched in repetition rate to the wavelength of the spacing of the plate exciters. The physical motion of the hydrogen and oxygen charged atoms being attracted to the opposite polarity zones will go into resonance. The self sustaining resonant motion of the hydrogen and oxygen atoms of the water molecule greatly enhances their disassociation from the water molecule."

The plate spacing is not arbitrary. It is the wavelength. Charged ions get attracted across the gap, overshoot, get pulled back, overshoot again. When the spacing matches the wavelength of that motion at the pulse frequency, they enter self-sustaining resonance.

The governing relationship:

spacing = drift velocity / pulse frequency

The drift velocity here is not the thermal velocity (~600 m/s) — it's the velocity of charged ions under the applied electric field. This is controllable, and it's how you tune the system.

---

Calculating Resonant Action for a 1/16" Gap

Using F = ma and the cavity spacing relationship, we can calculate the force and frequency needed for Stan's standard 1/16" tube gap:

Parameter	Value
Gap	1/16" = 1.587 mm
λ (spacing)	0.001587 m
f = v / λ	600 / 0.001587 = ~378 kHz
m(H2O)	2.99 × 10−26 kg
Amplitude (gap/2)	0.794 mm
ω = 2πf	2.376 × 106 rad/s
F = m · A · ω²	~1.34 × 10−16 N per molecule
E = F / q	~838 V/m
V = E × d	~1.3 volts to sustain resonance

The sustaining voltage appears tiny — and that's the point. You don't need kilovolts to sustain resonance. You need kilovolts to overcome damping, collisions, and initiate resonance in the first place. Once the molecule is oscillating resonantly, minimal energy maintains it.

---

Dual Resonance: The Unified System

This is the insight that ties everything together. There are two resonances that must be matched:

1. Physical (mechanical) resonance: the water molecule bouncing across the gap at 378 kHz
2. Electrical resonance: the VIC's LC tank circuit ringing at 378 kHz

When both are matched, maximum energy couples into the molecule at peak vulnerability.

Calculating the Choke Inductance

If mechanical resonance = 378 kHz and water cell capacitance ≈ 800 pF (typical for a 3" concentric tube cell), then:

f = 1 / (2π√LC)

Solving for L:
L = 1 / ((2πf)² × C)
L = 1 / ((2π × 378,000)² × 800 × 10−12)
L ≈ 221 μH

This is notably lower than the 500 μH – 2 mH values seen in most replication attempts. The reason: most builders tune to 40–70 kHz without matching the physical gap. Change the gap, you change everything.

---

The Dual Voltage Waveform

Stan's patent language from #4,798,661 describes the waveform strategy:

"The pulsating d.c. voltage and the duty cycle pulses have a maximum amplitude of the level that would cause electron leakage. Varying of the amplitude to an amplitude of maximum level to an amplitude below the maximum level of the pulses, provide an average amplitude below the maximum limit; but with the force of the maximum limit."

This is achieved with two variacs (0–120V each) and a flip-flop switching circuit:

• Peak voltage (Va): Hits the electron leakage threshold — maximum force. This kicks the molecule into oscillation at the resonant frequency. Think of it like striking a tuning fork.
• Low voltage (Vb): The duty cycle sustain level. Keeps the molecule oscillating without crossing into electron leakage territory. Like keeping a pendulum swinging with just enough push.

The flip-flop switches between these two voltage levels at the resonant frequency. You're not pulsing ON/OFF — you're pulsing between two precise voltage levels. The peak delivers maximum force while the duty cycle keeps average energy below the leakage threshold.

Finding Your Electron Leakage Threshold

As you increase the peak variac setting, watch for these indicators:

• Gas production climbs while current stays low — you're in the polarization regime
• Current draw suddenly climbs faster than gas production — you've crossed into electrolysis
• Water temperature begins rising (ohmic heating)
• A sharp "knee" appears on your ammeter curve

Back off just below that knee — that's your Va max. Lock it in, then use the second variac to set the lower sustain level.

---

The Complete Design Chain

Every parameter in the WFC connects to every other parameter. It is one unified system:

Gap spacing (1/16" = 1.587 mm)
  → Molecular resonant frequency (378 kHz)
    → Choke inductance (221 μH for 800 pF cell)
      → Drive frequency matches mechanical + electrical resonance
        → Peak voltage set at electron leakage threshold
          → Dual-variac waveform: peak force + duty cycle sustain
            → Molecular resonance driving (NOT electrolysis)

Most replication attempts treat these as separate problems — picking a gap, picking a frequency, winding a choke to whatever value, and hoping it works. The design chain above shows they are all interdependent. Start with your gap, derive everything else.

---

Volt-Seconds & Transformer Design

When designing the step-up transformer for the VIC, the core saturation limit is governed by volt-seconds:

B_peak = (V_in × t_on) / (N_primary × A_e)
N_min  = (V_in × t_on) / (B_sat × A_e)

A common question is whether turns ratio alone matters. It doesn't — 5:1, 50:10, and 500:100 are not the same design, even though the ratio is identical:

Configuration	Characteristics
5 : 1	Low inductance, requires higher frequency (100 kHz+), tight winding, low copper loss
50 : 10	10× primary inductance, handles lower frequencies, more copper, more inter-winding capacitance
500 : 100	Large core required, parasitic capacitance degrades pulse edges

The key relationships:

• Higher frequency = shorter t_on = fewer volt-seconds per cycle = fewer turns needed
• More turns = less flux per turn = lower frequency operation on the same core
• Optimum = where copper loss and core loss curves intersect

---

Patents Referenced

Patent	Title	Relevance
US #4,936,961	Method for Production of Fuel Gas	Primary VIC patent; EPP mechanism
US #4,798,661	Gas Generator Voltage Control Circuit	Cavity spacing = wavelength; dual voltage waveform
US #5,149,407	Process & Apparatus for Production of Fuel Gas	Polarization dependent on dielectric properties
EP0103656A2	Resonant Cavity for Hydrogen Generator	Corrugated exciter geometry (1982)
WO8912704A1	Process & Apparatus for Production of Fuel Gas	World patent; dielectric-dependent dissociation
Serial 06/367,052	Earlier corrugated surface exciter	Referenced as prior design in EP0103656A2

Appendices

Complete Formula Reference

Complete Formula Reference

This appendix provides a comprehensive reference of all formulas used in VIC circuit design and analysis. Formulas are organized by category for easy lookup.

1. Resonance Formulas

Formula	Equation	Units
Resonant Frequency	f₀ = 1 / (2π√(LC))	Hz
Angular Frequency	ω₀ = 2πf₀ = 1/√(LC)	rad/s
Period	T = 1/f₀ = 2π√(LC)	seconds
Inductance (given f₀, C)	L = 1 / (4π²f₀²C)	Henries
Capacitance (given f₀, L)	C = 1 / (4π²f₀²L)	Farads

2. Q Factor and Magnification

Formula	Equation	Notes
Q Factor (inductive)	Q = 2πfL / R = ωL/R	At frequency f
Q Factor (capacitive)	Q = 1 / (2πfCR) = 1/(ωCR)	At frequency f
Q from Z₀	Q = Z₀/R = (1/R)√(L/C)	Series RLC
Voltage Magnification	Vout = Q × Vin	At resonance
Characteristic Impedance	Z₀ = √(L/C)	Ohms

3. Bandwidth and Damping

Formula	Equation	Notes
Bandwidth (-3dB)	BW = f₀/Q = R/(2πL)	Hz
Decay Time Constant	τ = 2L/R	seconds
Damping Factor	α = R/(2L)	rad/s
Damped Frequency	fd = √(f₀² - α²/(4π²))	Hz
Ringdown Cycles (to 1%)	N ≈ 0.733 × Q	cycles

4. Capacitance Formulas

Formula	Equation	Notes
Parallel Plate	C = ε₀εrA/d	ε₀ = 8.854×10⁻¹² F/m
Concentric Cylinders	C = 2πε₀εrL / ln(ro/ri)	L = length
Capacitors in Series	1/Ctotal = 1/C₁ + 1/C₂ + ...
Capacitors in Parallel	Ctotal = C₁ + C₂ + ...
Energy in Capacitor	E = ½CV²	Joules

5. Inductance Formulas

Formula	Equation	Notes
Solenoid (air core)	L = μ₀N²A/l	μ₀ = 4π×10⁻⁷ H/m
Wheeler's Formula	L(µH) = N²r² / (9r + 10l)	r, l in inches
AL Method	L = AL × N²	AL in nH/turn²
Inductors in Series	Ltotal = L₁ + L₂ (no coupling)
Mutual Inductance	M = k√(L₁L₂)	k = coupling coefficient
Energy in Inductor	E = ½LI²	Joules

6. Resistance and Wire

Formula	Equation	Notes
Wire Resistance	R = ρL/A	ρ = resistivity
Wire Area (AWG)	A = π(d/2)²	d from wire tables
Skin Depth	δ = √(ρ/(πfμ))	meters
Copper Skin Depth	δ(mm) ≈ 66/√f(Hz)	Quick approximation
Power Dissipation	P = I²R = V²/R	Watts

7. Impedance Formulas

Element	Impedance	Phase
Resistor	Z = R	0°
Capacitor	Z = 1/(jωC) = -j/(2πfC)	-90°
Inductor	Z = jωL = j2πfL	+90°
CPE	Z = 1/(Q(jω)n)	-n×90°
Warburg	Z = σ/√ω × (1-j)	-45°

8. Electric Double Layer

Formula	Equation	Notes
Helmholtz Capacitance	CH = ε₀εrA/d	d ≈ 0.3 nm
Debye Length	λD ≈ 0.304/√c (nm)	c in mol/L
Total EDL (series)	1/C = 1/CStern + 1/Cdiff

9. Cole-Cole Model

Complex Permittivity:

ε* = ε_∞ + (ε_s - ε_∞) / [1 + (jωτ)^(1-α)]

Effective Capacitance:

C_eff(ω) = C₀ × [1 + (ωτ)^2(1-α)]^-1/2

10. Step Charging

Formula	Equation	Notes
Ideal N pulses	VC,N = 2N × Vs	Lossless
Maximum voltage	Vmax ≈ (4Q/π) × Vs	With losses
Half-cycle time	t = π√(LC)	For single pulse

Physical Constants

Constant	Symbol	Value
Permittivity of free space	ε₀	8.854 × 10⁻¹² F/m
Permeability of free space	μ₀	4π × 10⁻⁷ H/m
Relative permittivity (water)	εr	~80 at 20°C
Copper resistivity	ρCu	1.68 × 10⁻⁸ Ω·m
Elementary charge	e	1.602 × 10⁻¹⁹ C
Boltzmann constant	kB	1.381 × 10⁻²³ J/K

Reference complete. Use with the VIC Matrix Calculator for automated calculations.

Glossary of Terms

Appendix B: Wire Gauge & Material Tables

Complete reference tables for wire properties used in VIC choke design. All values at 20°C (68°F) unless noted.

AWG Wire Gauge Reference

AWG	Diameter (mm)	Diameter (in)	Area (mm²)	Area (kcmil)	Cu Ω/1000ft	Cu Ω/km
10	2.588	0.1019	5.261	10.38	0.9989	3.277
12	2.053	0.0808	3.309	6.530	1.588	5.211
14	1.628	0.0641	2.081	4.107	2.525	8.286
16	1.291	0.0508	1.309	2.583	4.016	13.17
18	1.024	0.0403	0.823	1.624	6.385	20.95
20	0.812	0.0320	0.518	1.022	10.15	33.31
22	0.644	0.0253	0.326	0.642	16.14	52.96
24	0.511	0.0201	0.205	0.404	25.67	84.22
26	0.405	0.0159	0.129	0.254	40.81	133.9
28	0.321	0.0126	0.081	0.160	64.90	212.9
30	0.255	0.0100	0.051	0.101	103.2	338.6
32	0.202	0.0080	0.032	0.063	164.1	538.3
34	0.160	0.0063	0.020	0.040	260.9	856.0
36	0.127	0.0050	0.013	0.025	414.8	1361
38	0.101	0.0040	0.008	0.016	659.6	2164
40	0.080	0.0031	0.005	0.010	1049	3441

Highlighted rows indicate commonly used gauges for VIC chokes.

Wire Material Resistivity

Material	Resistivity ρ (Ω·m)	Relative to Cu	Temp Coefficient α (/°C)
Silver (Ag)	1.59 × 10⁻⁸	0.95×	0.0038
Copper (Cu)	1.68 × 10⁻⁸	1.00× (reference)	0.00393
Gold (Au)	2.44 × 10⁻⁸	1.45×	0.0034
Aluminum (Al)	2.65 × 10⁻⁸	1.58×	0.00429
Brass	6-9 × 10⁻⁸	4-5×	0.002
Steel	1.0 × 10⁻⁷	6×	0.005
Stainless Steel	6.9 × 10⁻⁷	41×	0.001
Nichrome	1.1 × 10⁻⁶	65×	0.0004

Temperature Correction

Resistance at Temperature T:

R(T) = R₂₀ × [1 + α(T - 20)]

Example (Copper wire):

Magnet Wire Specifications

Magnet wire has enamel insulation. Overall diameter includes insulation:

AWG	Bare Dia. (mm)	Overall Dia. (mm)	Turns/cm	Turns/inch
18	1.024	1.09	9.2	23.3
20	0.812	0.87	11.5	29.2
22	0.644	0.70	14.3	36.3
24	0.511	0.56	17.9	45.4
26	0.405	0.45	22.2	56.4
28	0.321	0.36	27.8	70.6
30	0.255	0.29	34.5	87.6
32	0.202	0.24	41.7	106

Current Capacity Guidelines

For chassis wiring (in open air):

For coils, derate by 50% due to limited cooling. Magnet wire rated for higher temperature can handle more current.

Skin Depth Reference

At high frequencies, current flows near the wire surface. Skin depth δ:

δ = √(ρ / πfμ₀μᵣ)

Skin Depth in Copper:

Use wire diameter ≤ 2×δ for effective use of conductor cross-section. For larger currents at high frequencies, use Litz wire.

Quick Reference: DCR Calculation

For Copper Wire:

DCR (Ω) = Length (m) × Resistance (Ω/km) / 1000

DCR (Ω) = Length (ft) × Resistance (Ω/1000ft) / 1000

For Other Materials:

DCR_material = DCR_Cu × (ρ_material/ρ_Cu)

Wire Gauge Tables

Appendix C: Core Specifications

Reference specifications for magnetic cores commonly used in VIC choke design. Includes ferrite toroids, iron powder cores, and E-cores.

Core Material Overview

Material Type	μᵣ Range	Frequency Range	Best For
MnZn Ferrite	800-10,000	1 kHz - 2 MHz	High L, moderate f
NiZn Ferrite	15-1,500	500 kHz - 100 MHz	High frequency
Iron Powder	8-100	10 kHz - 10 MHz	High current, low cost
MPP (Molypermalloy)	14-550	DC - 1 MHz	Low loss, stable
Kool Mµ	26-125	DC - 500 kHz	High current, moderate loss
Air Core	1	Any	No saturation, linear

Common Ferrite Materials

MnZn Ferrite Materials

Material	μᵢ	Bsat (mT)	Frequency	Notes
Fair-Rite 77	2000	480	<1 MHz	General purpose, high μ
Fair-Rite 78	2300	480	<500 kHz	Very high μ
TDK N87	2200	490	<500 kHz	Popular, low loss
TDK N97	2300	410	<300 kHz	Very low loss
Ferroxcube 3C90	2300	470	<200 kHz	Low loss at high B
Ferroxcube 3F3	2000	440	<500 kHz	Higher frequency

Iron Powder Core Mix Chart

Iron powder cores (Micrometals/Amidon) are identified by color code:

Mix	Color	μᵣ	Frequency Range	Application
-26	Yellow/White	75	DC - 1 MHz	EMI/RFI filters
-2	Red/Clear	10	250 kHz - 10 MHz	RF, resonant circuits
-6	Yellow/Clear	8.5	3 - 40 MHz	Higher frequency
-1	Blue/Clear	20	500 kHz - 5 MHz	Medium frequency
-3	Gray/Clear	35	50 kHz - 500 kHz	Medium μ, low f
-52	Green/Blue	75	DC - 200 kHz	High μ, DC bias

Common Toroid Sizes

FT (Ferrite Toroid) Series

Size	OD (mm)	ID (mm)	H (mm)	Aₗ (77 mat)	Aₗ (43 mat)
FT-37	9.5	4.7	3.2	884	440
FT-50	12.7	7.1	4.8	1140	570
FT-82	21.0	13.0	6.4	2170	557
FT-114	29.0	19.0	7.5	2640	603
FT-140	35.5	23.0	12.7	3170	885
FT-240	61.0	35.5	12.7	4820	1075

Aₗ values in nH/turn². Highlighted sizes are commonly used for VIC chokes.

T (Iron Powder Toroid) Series

Size	OD (mm)	ID (mm)	H (mm)	Aₗ (-2 mix)	Aₗ (-26 mix)
T-37	9.5	4.9	3.2	4.0	27
T-50	12.7	7.7	4.8	4.9	33
T-68	17.5	9.4	4.8	5.7	38
T-80	20.2	12.6	6.4	8.5	55
T-94	24.0	14.5	7.9	8.4	70
T-106	26.9	14.0	11.1	13.5	90
T-130	33.0	19.7	11.1	11.0	96
T-200	50.8	31.8	14.0	12.0	120

Inductance Calculations

Using Aₗ Value:

L (nH) = Aₗ × N²

N = √(L / Aₗ)

Example:

Saturation Considerations

Saturation Flux Density (B_sat):

Calculating Peak Flux:

B = (V × t) / (N × A_e)

Where A_e is effective core area. Keep B < 0.5 × B_sat for linear operation.

Temperature Effects

Material	Curie Temp (°C)	Max Operating (°C)	μ vs. Temp
MnZn Ferrite	200-250	100-120	Peaks near 80°C, then drops
NiZn Ferrite	300-500	150	Relatively stable
Iron Powder	770 (iron)	125 (coating limited)	Stable

Core Selection Guide for VIC

For Primary Choke (L1):

For Secondary Choke (L2):

For High Frequency (>100 kHz):

Quick Reference: Turns Calculation

Approximate values. Verify with actual Aₗ from manufacturer datasheet.

Core Specifications

Glossary of Terms

A comprehensive glossary of technical terms used throughout the VIC Matrix educational content and calculator.

A

A_L (Inductance Factor)

A core specification in nH/turn² that allows quick calculation of inductance: L = A_L × N²

Alpha (α) - Cole-Cole

Distribution parameter (0-1) in the Cole-Cole model. α=0 is ideal Debye relaxation; higher values indicate broader distribution of relaxation times.

Alpha (α) - Damping

Damping factor in an RLC circuit: α = R/(2L). Determines how quickly oscillations decay.

Amplitude

The maximum value of an oscillating quantity, such as voltage or current.

B

Bandwidth (BW)

The frequency range over which a resonant circuit responds effectively. BW = f₀/Q for a series RLC circuit.

Bifilar Winding

A winding technique where two wires are wound together in parallel, creating tight magnetic coupling and significant inter-winding capacitance.

Blocking Electrode

An electrode where no Faradaic (electrochemical) reactions occur, behaving purely as a capacitor.

C

Capacitance (C)

The ability to store electric charge. Measured in Farads (F). C = Q/V where Q is charge and V is voltage.

Characteristic Impedance (Z₀)

The ratio √(L/C) for an LC circuit. Represents the impedance level of the resonant system.

Charge Transfer Resistance (R_ct)

The resistance associated with electron transfer at an electrode surface during electrochemical reactions.

Choke

An inductor used in a circuit to block or impede certain frequencies while allowing others to pass. In VIC context, the resonating inductors.

Cole-Cole Model

A mathematical model describing frequency-dependent dielectric behavior with distributed relaxation times.

Constant Phase Element (CPE)

A circuit element with impedance Z = 1/[Q(jω)ⁿ], used to model non-ideal capacitor behavior in electrochemical systems.

Coupling Coefficient (k)

A measure of magnetic coupling between inductors (0-1). k = M/√(L₁L₂) where M is mutual inductance.

D

DCR (DC Resistance)

The resistance of an inductor measured with direct current. Primary contributor to inductor losses.

Debye Length (λ_D)

The characteristic thickness of the diffuse layer in an electrochemical double layer. Decreases with increasing ion concentration.

Diffuse Layer

The outer region of the electric double layer where ion concentration gradually returns to bulk values.

Dielectric

An insulating material that can be polarized by an electric field. Water is a dielectric with high permittivity (ε_r ≈ 80).

Double Layer

See Electric Double Layer (EDL).

E

EDL (Electric Double Layer)

The structure formed at an electrode-electrolyte interface, consisting of a compact layer of ions and a diffuse layer extending into solution.

EIS (Electrochemical Impedance Spectroscopy)

A technique for characterizing electrochemical systems by measuring impedance across a range of frequencies.

ESR (Equivalent Series Resistance)

The resistive component of a capacitor's impedance, causing power dissipation.

F

Faradaic Reaction

An electrochemical reaction involving electron transfer at an electrode, such as water electrolysis.

Ferrite

A ceramic magnetic material used for inductor cores, suitable for high-frequency applications.

Frequency (f)

The number of complete oscillation cycles per second. Measured in Hertz (Hz).

G-H

Helmholtz Layer

The compact inner layer of the EDL, where ions are closest to the electrode surface.

Hysteresis

Energy loss in magnetic materials due to the lag between applied field and magnetization.

I

Impedance (Z)

The total opposition to alternating current, including both resistance and reactance. Measured in Ohms (Ω).

Inductance (L)

The property of a conductor that opposes changes in current by storing energy in a magnetic field. Measured in Henries (H).

IHP (Inner Helmholtz Plane)

The plane passing through the centers of specifically adsorbed ions in the EDL.

L-M

LC Circuit

A circuit containing an inductor and capacitor, capable of oscillating at a resonant frequency.

Mutual Inductance (M)

The inductance linking two coils, allowing energy transfer between them.

N-O

Nyquist Plot

A plot of imaginary vs. real impedance (-Z'' vs Z') used in EIS analysis.

OHP (Outer Helmholtz Plane)

The plane of closest approach for solvated (hydrated) ions in the EDL.

P

Parasitic Capacitance

Unintended capacitance in an inductor, arising from turn-to-turn and layer-to-layer effects.

Permittivity (ε)

A measure of how much electric field is reduced in a material compared to vacuum. ε = ε₀ε_r.

Permeability (μ)

A measure of how well a material supports magnetic field formation. μ = μ₀μ_r.

PLL (Phase-Locked Loop)

A control system that maintains frequency lock with a reference signal, used to track resonance.

Q

Q Factor (Quality Factor)

A dimensionless parameter indicating the "sharpness" of resonance. Q = ωL/R = Z₀/R. Higher Q means narrower bandwidth and higher voltage magnification.

R

Randles Circuit

An equivalent circuit model for electrochemical cells consisting of R_s, C_dl, R_ct, and Z_W.

Reactance

The imaginary part of impedance. Inductive reactance X_L = ωL; capacitive reactance X_C = 1/(ωC).

Resonance

The condition where inductive and capacitive reactances are equal, resulting in maximum energy storage and voltage magnification.

Ringdown

The decay of oscillations after excitation stops, characterized by the time constant τ = 2L/R.

S

Self-Resonant Frequency (SRF)

The frequency at which an inductor's parasitic capacitance resonates with its inductance. Above SRF, the inductor behaves as a capacitor.

Skin Effect

The tendency of AC current to flow near the surface of a conductor, increasing effective resistance at high frequencies.

Solution Resistance (R_s)

The ionic resistance of the electrolyte between electrodes.

Step Charging

A technique using multiple resonant pulses to progressively build voltage on a capacitor.

Stern Layer

The combined compact and diffuse layer model of the EDL.

T

Tank Circuit

A parallel LC circuit that "tanks" or stores energy, oscillating between magnetic and electric forms.

Tau (τ) - Time Constant

The characteristic time for decay. For an RLC circuit: τ = 2L/R.

Toroidal Core

A doughnut-shaped magnetic core providing a closed magnetic path and good field containment.

V

VIC (Voltage Intensifier Circuit)

A resonant circuit configuration using chokes and capacitors to develop high voltage across a water fuel cell.

Voltage Magnification

The ratio of voltage across a reactive element to the source voltage at resonance. Equals Q for a series RLC circuit.

W

Warburg Impedance (Z_W)

Impedance arising from diffusion of electroactive species, characterized by 45° phase angle and Z ∝ 1/√ω.

WFC (Water Fuel Cell)

An electrochemical cell where water serves as the medium between electrodes, acting as a capacitive-resistive load in VIC circuits.

Z

Z₀ (Characteristic Impedance)

The natural impedance level of an LC circuit: Z₀ = √(L/C). Also Q × R for a series RLC circuit.

Zero-Current Switching (ZCS)

A switching technique where transistors turn off when current is zero, minimizing switching losses.

Glossary compiled for the VIC Matrix educational series.