Advanced Topics
- PLL Control
- Harmonic Analysis
- Transformer Coupling
- Energy Efficiency
- Experimental Validation
- Understanding Resonant Action in the Water Fuel Cell
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:
- Phase Detector: Compares phase of two signals, outputs error voltage
- Loop Filter: Averages error signal, sets response speed
- VCO: Voltage-Controlled Oscillator, frequency varies with input voltage
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:
- f < f₀: V leads I (capacitive), phase > 0°
- f = f₀: V and I in phase, phase = 0°
- f > f₀: V lags I (inductive), phase < 0°
Control Law:
- If phase > 0°: Increase frequency (move toward resonance)
- If phase < 0°: Decrease frequency (move toward resonance)
- If phase ≈ 0°: Maintain frequency (at resonance)
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:
- 555 Timer VCO: Simple, wide frequency range, moderate stability
- 74HC4046 PLL IC: Integrated PLL with VCO, easy to use
- DDS (Direct Digital Synthesis): Precise frequency control, programmable
- Microcontroller PWM: Software-adjustable, flexible
VCO Requirements:
- Frequency range covering expected f₀ ± drift range
- Linear frequency vs. voltage response
- Low noise and jitter
- Fast frequency settling
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:
- Initialize VCO near expected f₀
- Enable PLL with wide bandwidth initially
- Wait for lock indication
- Reduce bandwidth for stable operation
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) = (4Vpk/π)[sin(ωt) + (1/3)sin(3ωt) + (1/5)sin(5ωt) + ...]
| Harmonic | Frequency | Relative Amplitude |
|---|---|---|
| 1st (fundamental) | f | 100% |
| 3rd | 3f | 33.3% |
| 5th | 5f | 20% |
| 7th | 7f | 14.3% |
Pulse Train (Variable Duty Cycle)
Pulse train with duty cycle D contains both odd and even harmonics:
an = (2Vpk/nπ) × sin(nπD)
Effect of Duty Cycle:
- D = 50%: Only odd harmonics (even harmonics cancel)
- D = 25%: Strong 2nd harmonic, weak 4th
- D = 33%: No 3rd harmonic (3rd harmonic null)
- Narrow pulse: Wide harmonic spectrum, many significant harmonics
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):
- At 10 kHz (1st): Response = 100%
- At 30 kHz (3rd): Response ≈ 0.6%
- At 50 kHz (5th): Response ≈ 0.2%
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
- Circuit resonance: f₀ = 30 kHz
- Drive frequency: fdrive = 10 kHz
- 3rd harmonic of drive (30 kHz) excites resonance
Advantages:
- Lower switching frequency (easier on semiconductors)
- Different pulse characteristics
- May interact differently with WFC
Disadvantages:
- Harmonic has lower amplitude than fundamental
- Reduced efficiency (energy in unused harmonics)
- More complex analysis
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:
- k = 0: No coupling (independent inductors)
- k = 0.01-0.1: Loose coupling (separate cores, some proximity)
- k = 0.5-0.8: Moderate coupling (shared core, separate windings)
- k = 0.95-0.99: Tight coupling (bifilar, interleaved windings)
- k = 1: Perfect coupling (theoretical ideal transformer)
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₂Cwfc))
With Coupling:
The system has two coupled resonant modes. The frequencies split into:
f₁, f₂ = function of L₁, L₂, C₁, Cwfc, 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:
- Large mode splitting
- Efficient energy transfer between windings
- Built-in inter-winding capacitance
- Lower overall SRF due to capacitance
Measuring Bifilar Coupling:
- Measure Lseries-aid (windings in series, same polarity)
- Measure Lseries-opp (windings in series, opposite polarity)
- Calculate: M = (Lseries-aid - Lseries-opp) / 4
- Calculate: k = M / √(L₁ × L₂)
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:
- Compact design (bifilar combines L1 and L2)
- Intentional transformer action desired
- Specific mode-splitting behavior needed
When to Avoid Coupling:
- Independent tuning of primary and secondary needed
- Simpler analysis desired
- Want predictable single-resonance behavior
Layout Guidelines:
- Toroidal cores have low external field—good for isolation
- Orient coils perpendicular to minimize stray coupling
- Use shielding if isolation is critical
- Measure actual coupling to verify assumptions
Analyzing Coupled VIC Circuits
Coupled Circuit Analysis Steps:
- Measure or estimate coupling coefficient k
- Convert to T-equivalent model
- Analyze as three-inductor circuit
- Or use simulation with mutual inductance
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:
EL = ½LI² (energy in inductor)
EC = ½CV² (energy in capacitor)
At Resonance:
Etotal = EL,max = EC,max = ½CVpeak²
Peak Energy (example):
- C = 10 nF, Vpeak = 1000 V
- E = ½ × 10×10⁻⁹ × 1000² = 5 mJ
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:
ΔEcycle = 2π × Estored / Q
Interpretation:
- Q = 10: Lose 63% of energy per cycle
- Q = 50: Lose 13% of energy per cycle
- Q = 100: Lose 6% of energy per cycle
- Q = 200: Lose 3% of energy per cycle
Energy Retention:
After n cycles: E(n) = E₀ × e^(-2πn/Q)
Power Flow Analysis
Input Power
Pin = Vin × Iin × cos(φ)
For pulsed operation:
Pavg = (1/T) × ∫V(t)I(t)dt
Dissipated Power
Pdiss = Irms² × Rtotal
Where Rtotal = RDCR1 + RDCR2 + Rsol + Rother
Useful Power
Power available for electrochemical work:
Puseful = Pin - Pdiss
Or, for the WFC specifically:
Pwfc = Vwfc × Iwfc × 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:
Pin = PDCR1 + PDCR2 + Psol + Pcore + Pother
Check:
- Sum all loss mechanisms
- Compare to measured input power
- Large discrepancy indicates missing loss or measurement error
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:
- Reduce largest loss first: In example above, Rsol is 40%—optimize water conductivity
- Use larger wire: Each AWG step down reduces DCR by ~25%
- Choose better core: Low-loss ferrite vs. iron powder
- Optimize water conductivity: Not too high (electrolysis), not too low (resistance loss)
- Reduce connection resistance: Good solder joints, clean contacts
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)
- Set LCR meter to inductance mode
- Select test frequency (1 kHz typical)
- Connect inductor, read value
- Repeat at 10 kHz to check for frequency dependence
Method 2: Resonance with Known C
- Connect inductor with known capacitor C
- Drive with function generator, sweep frequency
- Find resonant frequency f₀ (voltage peak)
- Calculate: L = 1/(4π²f₀²C)
Measuring DCR
Four-Wire (Kelvin) Measurement:
For accurate low-resistance measurement, use 4-wire method to eliminate lead resistance:
- Use dedicated low-ohm meter
- Or use LCR meter in R mode
- Allow reading to stabilize (self-heating)
Expected accuracy: ±1-5% compared to calculated value
Measuring WFC Capacitance
- Fill WFC with water at operating temperature
- Measure with LCR meter at 1 kHz and 10 kHz
- Values should be similar (if EDL effects are small)
- Note the ESR reading as well
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:
- Set function generator to low amplitude sine wave
- Start at low frequency (1/10 of expected f₀)
- Slowly increase frequency while watching Ch2 amplitude
- Note frequency of maximum amplitude—this is f₀
- Also note -3dB frequencies (where amplitude = 0.707 × peak)
Calculate Q from Measurement:
Q = f₀ / (fhigh - flow) = f₀ / BW
Phase Measurement Method
- Display both input current and output voltage
- Use X-Y mode or measure phase with oscilloscope
- At resonance, phase difference = 0°
- More accurate than amplitude peak for high-Q circuits
Q Factor Measurement
Method 1: Bandwidth
Measure -3dB bandwidth and calculate:
Q = f₀ / BW
Method 2: Ring-Down
- Excite circuit with single pulse at f₀
- Observe decaying oscillation on oscilloscope
- Count cycles to decay to 1/e (37%)
- Q ≈ π × (number of cycles to 1/e decay)
Alternatively, measure time constant τ:
τ = 2L/R = Q/(πf₀)
Method 3: Voltage Magnification
- Measure input voltage Vin
- Measure output voltage Vout at resonance
- Q ≈ Vout/Vin
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:
- Stray capacitance adding to total C
- Actual L higher than calculated
- Check for loose connections (add L)
Measured f₀ Higher than Calculated:
- Actual L lower (core saturation, wrong μᵣ)
- WFC capacitance overestimated
- Air bubbles reducing effective C
Measured Q Lower than Calculated:
- Additional losses not accounted for
- Core losses at operating frequency
- Poor connections adding resistance
- Radiation losses at high frequency
No Clear Resonance Observed:
- Operating above SRF (choke is capacitive)
- Very low Q (Q < 2) makes resonance hard to see
- Measurement setup loading the circuit
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:
- VIC circuits can develop high voltages at resonance
- Always use proper high-voltage probes
- Keep one hand in pocket when probing live circuits
- Discharge capacitors before handling
⚠️ Gas Production:
- WFC produces hydrogen and oxygen—ensure ventilation
- No open flames or sparks near operating cell
- Use appropriate gas collection if needed
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:
- Voltage ionizes the molecule → creates H+ and OH− carriers
- Conductivity goes up → loss tangent (ε'') rises → Q factor drops
- 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:
- Physical (mechanical) resonance: the water molecule bouncing across the gap at 378 kHz
- 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 |