Advanced WFC Concepts 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: f₀ = resonant frequency (Hz) L = inductance (Henries) C = capacitance (Farads) Why Resonance Matters for VIC Circuits In Stan Meyer's Voltage Intensifier Circuit (VIC), resonance is the key mechanism that enables: Voltage Magnification: At resonance, voltages across reactive components can be many times greater than the input voltage Efficient Energy Transfer: Energy oscillates between the inductor's magnetic field and the capacitor's electric field with minimal loss 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 ): Z = R (minimum impedance) Current = V/R (maximum current) Voltage across L = Voltage across C = Q × V source Series LC Behavior Frequency Condition Circuit Behavior f < f₀ X C > X L Capacitive (current leads voltage) f = f₀ X C = X L Resistive (current in phase with voltage) f > f₀ X L > X C Inductive (current lags voltage) Parallel LC Circuit Circuit Configuration: L and C connected in parallel At Resonance: Impedance approaches infinity (in ideal case) Current from source is minimum Large circulating current flows between L and C 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: The impedance at resonance for a parallel LC circuit The ratio of voltage to current in a traveling wave A design parameter for matching circuits Energy Transfer in LC Circuits In an ideal LC circuit (no resistance), energy oscillates perpetually between the inductor and capacitor: Capacitor fully charged: All energy stored in electric field (E = ½CV²) Current building: Energy transferring to inductor Maximum current: All energy stored in magnetic field (E = ½LI²) Current decreasing: Energy transferring back to capacitor 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: f₀ = resonant frequency (Hz) L = inductance (Henries) R = total series resistance (Ohms) C = capacitance (Farads) ω = 2πf₀ (angular frequency) Z₀ = √(L/C) (characteristic impedance) 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: Frequency sweep method: Find f₀ and the -3dB points, then Q = f₀/BW Ring-down method: Count cycles for amplitude to decay to 1/e (37%) 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: BW = bandwidth in Hz f₀ = resonant frequency Q = quality factor R = total series resistance L = inductance 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: Q = 10: ≈7.3 cycles to 1% Q = 50: ≈36.7 cycles to 1% Q = 100: ≈73.3 cycles to 1% 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: Apply a burst of oscillations at the resonant frequency Stop the driving signal and observe the decay on an oscilloscope Count the number of cycles for amplitude to drop to 37% (1/e) 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): Higher L with lower C → Higher Z₀ → Higher magnification (but more wire, more DCR) Lower L with higher C → Lower Z₀ → Lower magnification (but less wire, less DCR) The optimal design balances these factors. Energy Considerations Voltage magnification doesn't violate energy conservation: Power In = Power Dissipated At steady-state resonance: Current through circuit: I = V source /R Power from source: P = V source × I = V source ²/R Power dissipated in R: P = I²R = V source ²/R (same!) 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: A 12V source with Q=50 produces 600V peaks These voltages can cause electric shock or burns Energy stored in capacitors remains after power is removed Always discharge capacitors before handling circuits Use appropriate insulation and safety equipment 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: Charge Separation: The electrode surface may carry an electrical charge (positive or negative) Ion Attraction: Ions of opposite charge in the solution are attracted to the electrode surface Solvent Molecules: Water molecules orient themselves in the electric field near the surface 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: ε₀ = permittivity of free space (8.854 × 10⁻¹² F/m) ε r = relative permittivity of the layer (~6-10 for water near electrode) A = electrode area d = thickness of the double layer (~0.3-0.5 nm) Typical EDL Capacitance Values Because the separation distance is so small (nanometers), EDL capacitance is remarkably high: System Typical C dl 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: Parallel plate capacitor (1mm gap, air): ~0.0088 µF/cm² Electric Double Layer (~0.3nm gap, water): ~20 µF/cm² EDL is about 2,000× higher capacitance per unit area! EDL in Water Fuel Cells In a water fuel cell, the EDL forms at both electrodes: Anode (positive electrode): Attracts negative ions (OH⁻, Cl⁻ if present) 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: C geo = geometric (parallel-plate) capacitance C edl,anode = double layer capacitance at anode C edl,cathode = double layer capacitance at cathode 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 C dl (µ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: c dl = specific EDL capacitance (µF/cm²) A = electrode surface area (cm²) Example Calculation Given: Electrode area: 100 cm² Electrode gap: 1 mm c dl : 25 µF/cm² (for stainless steel) 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 C geo becomes very large Large electrode gap (>5mm) Minimal C geo is small, dominates total Small electrode area Significant C edl becomes comparable to C geo 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 C dl Effect Deionized <1 µS/cm ~100 nm Lower C dl Distilled 1-10 µS/cm ~30 nm Moderate C dl Tap water 200-800 µS/cm ~1 nm Higher C dl With electrolyte (NaOH, KOH) >1000 µS/cm <1 nm Highest C dl 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: The electrode surface carries a uniform charge Counter-ions in solution form a single plane at a fixed distance from the electrode No ions exist between the electrode and this plane The potential drops linearly between the electrode and ion plane 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 = ε₀ε r A / d Where: ε₀ = 8.854 × 10⁻¹² F/m (permittivity of free space) ε r = relative permittivity of the inner layer (~6-10) A = electrode surface area d = distance from electrode to ion centers (~0.3-0.5 nm) 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: ε r = 6 (strongly oriented water) d = 0.3 nm = 3 × 10⁻¹⁰ m 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: Maximum possible EDL capacitance: Sets an upper bound on what the interface can contribute Field strength at the electrode: Related to the electrochemical driving force 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) A layer of specifically adsorbed ions and solvent molecules Extends from electrode surface to the Outer Helmholtz Plane (OHP) No free charges within this region Potential drops linearly (like Helmholtz) 2. Diffuse Layer (Gouy-Chapman Layer) Begins at the OHP and extends into solution Ion concentration follows Boltzmann distribution Potential decays exponentially Thickness characterized by the Debye length 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 = √(ε₀ε r k B T / (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 = ε₀ε 1 A / d Diffuse Layer Capacitance: C diffuse = (ε₀ε r A / λ D ) × cosh(zeφ d /2k B T) 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: Debye length: λ D ∝ √T (diffuse layer thickens at higher T) Dielectric constant: ε r decreases with T Thermal voltage: k B T/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: C dl1 , C dl2 : Double layer capacitances at each electrode R ct1 , R ct2 : Charge transfer resistances (reaction kinetics) W₁, W₂: Warburg impedances (diffusion) R sol : Solution resistance C geo : Geometric capacitance 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 R ct Resistive behavior Medium (100 Hz - 10 kHz) EDL capacitance C dl Capacitive, EDL dominant High (10 kHz - 1 MHz) Solution R + geometric C RC network behavior Very high (>1 MHz) Geometric C geo 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: R sol = 100 Ω (tap water, small cell) C dl = 10 µF τ EDL = 100 × 10×10⁻⁶ = 1 ms Full formation time ≈ 5 ms 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: At 1 kHz: C dl,eff ≈ 0.3-0.7 × C dl (DC) At 10 kHz: C dl,eff ≈ 0.1-0.3 × C dl (DC) At 50 kHz: C dl,eff ≈ 0.05-0.15 × C dl (DC) 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: Resonant frequency shift during operation Detuning from optimal operating point Need for adaptive frequency control (PLL) 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: Narrow gap (0.5-1mm): Higher C geo , but higher R sol , risk of bridging Wide gap (3-5mm): Lower C geo , lower R sol , easier construction Optimal (1-2mm): Balances capacitance, resistance, and practicality 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: Use an LCR meter: Measure at multiple frequencies (100 Hz, 1 kHz, 10 kHz) Perform EIS: Electrochemical Impedance Spectroscopy gives complete picture Measure at operating conditions: Temperature and voltage matter 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: |Z| = impedance magnitude (Ohms) θ = phase angle between voltage and current Z' = real part (resistance-like) Z'' = imaginary part (reactance-like) j = √(-1) (imaginary unit) 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: Separating processes: Different phenomena occur at different frequencies Non-destructive: Small AC signals don't significantly perturb the system Complete characterization: Maps all electrical behavior across frequency 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: Apply small AC voltage (5-50 mV) superimposed on DC bias Sweep frequency from high to low (e.g., 1 MHz to 0.01 Hz) Measure current response at each frequency Calculate impedance Z = V/I at each frequency Plot results as Nyquist or Bode diagrams 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) Flat regions indicate resistive behavior Slope of -1 indicates capacitive behavior Slope of +1 indicates inductive behavior Bode Phase Plot: Phase angle θ vs. frequency (log scale) θ = 0° indicates resistive θ = -90° indicates capacitive θ = +90° indicates inductive Frequency Ranges and Processes Different electrochemical processes dominate at different frequencies: Frequency Process Circuit Element > 100 kHz Bulk solution, cables R s , parasitic L 1 kHz - 100 kHz Double layer charging C dl 1 Hz - 1 kHz Charge transfer kinetics R ct < 1 Hz Mass transport (diffusion) Z W (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: Potentiostat with EIS capability (or dedicated EIS analyzer) Three-electrode setup (working, counter, reference) Shielded cables to minimize noise Faraday cage for low-frequency measurements 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) R s Ionic resistance of electrolyte solution between electrodes 10 Ω - 10 kΩ (depends on conductivity) C dl Electric double layer capacitance at electrode surface µF to mF range (depends on area) R ct Resistance to electron transfer at electrode (reaction kinetics) 1 Ω - 1 MΩ (depends on overpotential) Z W 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 R s Water conductivity, electrode gap Adds to total circuit resistance, reduces Q C dl EDL at each electrode Part of total WFC capacitance R ct Activation barrier for water splitting Limits DC current, less relevant at high freq Z W 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: R s : High-frequency real-axis intercept R ct : Diameter of the semicircle C dl : From peak frequency: C = 1/(2πf peak R ct ) 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: ε ∞ = high-frequency (optical) permittivity ε s = static (DC) permittivity τ = characteristic relaxation time α = Cole-Cole parameter (0 ≤ α < 1) ω = angular frequency (2πf) 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: alpha (α) Distribution parameter (0-1) tau (τ) Characteristic time constant (seconds) epsilon_s Static permittivity epsilon_inf High-frequency permittivity Frequency-Dependent Capacitance The Cole-Cole model predicts how capacitance varies with frequency: Effective Capacitance: C eff (ω) = C 0 × [1 + (ωτ) 2(1-α) ] -1/2 At low frequency: C eff → C 0 (full capacitance) At high frequency: C eff → C ∞ < C 0 (reduced capacitance) Practical Example WFC with Cole-Cole Parameters: τ = 10 µs (characteristic frequency ~16 kHz) α = 0.4 (moderate distribution) C 0 = 10 nF (DC capacitance) Effective Capacitance at Different Frequencies: Frequency ωτ C eff 100 Hz 0.006 ~10 nF (98%) 1 kHz 0.063 ~9.5 nF (95%) 10 kHz 0.63 ~7.5 nF (75%) 50 kHz 3.14 ~4 nF (40%) VIC Design Implications The Cole-Cole model affects VIC design in several ways: Resonant frequency shift: As frequency changes, C eff changes, shifting resonance Broader resonance: The distribution of time constants broadens the frequency response Q factor reduction: Losses associated with the relaxation reduce circuit Q 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: σ = Warburg coefficient (Ω·s -1/2 ) ω = angular frequency (rad/s) j = imaginary unit 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: R = gas constant (8.314 J/mol·K) T = temperature (K) n = number of electrons transferred F = Faraday constant (96485 C/mol) A = electrode area D O , D R = diffusion coefficients of oxidized/reduced species C O , C R = bulk concentrations 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 |Z W | 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: |Z W | ∝ 1/√f decreases rapidly with frequency At 10 kHz: |Z W | is ~100× smaller than at 1 Hz Diffusion processes can't keep up with rapid voltage changes When Warburg Matters: Very low frequency operation (<10 Hz) Step-charging with long dwell times DC bias measurements Diagnosing electrode fouling or gas buildup Practical Implications Frequency selection: High-frequency operation minimizes diffusion effects Bubble management: Gas bubbles increase Warburg impedance Electrode design: Porous electrodes have complex diffusion paths Stirring/flow: Can reduce diffusion limitations Measuring Warburg Parameters To characterize the Warburg element in your WFC: Perform EIS down to very low frequencies (0.01 Hz) Look for the 45° line region in Nyquist plot Measure the slope to determine σ 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ω) n ] Where: Q = CPE coefficient (units: S·s n or F·s (n-1) ) n = CPE exponent (0 ≤ n ≤ 1) ω = angular frequency (rad/s) Magnitude and Phase: |Z CPE | = 1 / (Qω n ) θ = -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: Frequency-dependent capacitance: C eff = Qω (n-1) means capacitance varies with operating frequency Resonant frequency prediction: Must account for CPE when calculating f₀ Q factor effects: The lossy nature of CPE (when n < 1) reduces circuit Q Surface treatment: Smoother electrodes (higher n) behave more like ideal capacitors Measuring CPE Parameters To determine Q and n for your WFC: Perform EIS measurement across relevant frequency range Fit data to modified Randles circuit with CPE Extract Q and n from fitting software 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 Maximize Q factor: Higher Q = more voltage magnification Achieve resonance: All components tuned to operating frequency Match impedances: Efficient energy transfer between stages Maintain stability: Prevent frequency drift and oscillation problems 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: f₀ = 10 kHz, L1 = 10 mH, R1 = 10 Ω X L1 = 2π × 10,000 × 0.01 = 628 Ω Q = 628 / 10 = 62.8 With 12V input: V C1 = 62.8 × 12 = 754V 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 Measure L1 accurately: Use an LCR meter at multiple frequencies Start with calculated C1: Then fine-tune for best response Use variable capacitor or parallel caps: For easy tuning Check for SRF: Ensure L1's SRF is well above f₀ 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: Geometric capacitance: C geo = ε₀ε r A/d EDL capacitance: C edl (in series, at each electrode) Effective capacitance: C wfc = f(C geo , C edl , frequency) 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: Q L1C = 30, Q L2 = 20, V in = 12V V WFC = 30 × 20 × 12 = 7,200V theoretical Cascaded Resonance Effects When both stages resonate at the same frequency, the effects multiply: Configuration Total Magnification Notes Only primary resonance Q L1C L2-WFC not tuned Only secondary resonance Q L2 L1-C1 not tuned Dual resonance Q L1C × Q L2 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 C wfc Lower f₀, lower Z₀ Increase L2 or reduce C1 Higher R wfc Lower Q L2 Use purer water or optimize gap Larger electrode area Higher C wfc Requires larger L2 Narrower gap Higher C wfc , lower R wfc 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: f₀ = 10 kHz C wfc = 5 nF (typical small WFC) L2 = 1 / (4π² × 10⁴² × 5×10⁻⁹) = 50.7 mH 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: V WFC = 1000V, d = 1mm E = 1000V / 0.001m = 1 MV/m = 10 kV/cm This is a substantial electric field that can influence molecular behavior in water. Design Guidelines for L2 Match resonant frequency: L2 should resonate with C wfc at the same frequency as L1-C1 Minimize DCR: R2 directly reduces Q L2 and thus voltage magnification Consider coupling: If using transformer-coupled design, mutual inductance matters Account for WFC changes: C wfc varies with temperature, voltage, and bubble formation 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 = V source Can exceed V source (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: Initially: All energy in source (voltage V s ) Quarter cycle: Energy split between L (current max) and C (V = V s ) Half cycle: All energy in C, current = 0 For energy to be conserved: ½CV c ² = C×V s ² (accounting for work done by source) This gives V c = 2V s 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 V C,max /V source 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: Primary tank: C1 is resonantly charged through L1 Secondary transfer: Energy transfers resonantly to WFC through L2 Cumulative effect: Multiple stages multiply the magnification Timing and Switching For optimal resonant charging: Critical Timing Points: Turn-on: When capacitor voltage is minimum (or at desired starting point) Turn-off: When current through inductor reaches zero (zero-current switching) Period: Should match or be a harmonic of the resonant frequency 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: Pulse 1: Capacitor charges from 0 to 2V s (resonant half-cycle) Hold: Diode prevents discharge back through inductor Pulse 2: Starting from 2V s , capacitor charges to ~4V s Hold: Energy stored, waiting for next pulse Continue: Each pulse adds ~2V s (minus losses) 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 V max /V source 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 × V s (AC peak) (4Q/π) × V s (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: Pulse duration: π√(LC) = half resonant period Pulse frequency: f pulse < f resonant /2 Duty cycle: Typically 10-50% Gap between pulses: Allow ring-down and settling 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: High-speed switching: MOSFET or IGBT driver Precise timing: Microcontroller or pulse generator High-voltage diode: Fast recovery, rated for expected voltages 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: N = number of turns Φ = magnetic flux I = current μ = permeability of core material A = cross-sectional area of core l = magnetic path length Key Inductor Parameters Parameter Symbol Units Importance Inductance L Henry (H) Determines resonant frequency with C DC Resistance DCR, R dc 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 I sat Amps (A) Max current before inductance drops Inductor Construction A practical inductor consists of: Wire: Conductor wound into coils (turns) Core: Material inside the coil (air, ferrite, iron, etc.) 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 = (μ₀μ r N²A) / (2πr mean ) DC Resistance (DCR) The DC resistance is determined by the wire properties: R dc = ρ × l wire / A wire Where: ρ = resistivity of wire material (Ω·m) l wire = total wire length ≈ N × π × d coil A wire = wire cross-sectional area Q Factor of Inductors Inductor Q Factor: Q = ωL/R = 2πfL/R total R total includes: DC resistance of wire Skin effect losses (increases with frequency) Proximity effect losses Core losses (hysteresis + eddy currents) 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 Target inductance: Sets resonant frequency with capacitor Low DCR: Maximizes Q factor High SRF: Ensures proper operation at intended frequency Adequate current rating: Won't saturate or overheat 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 (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: Copper wire Q = 100 → SS316 wire Q ≈ 2.3 Copper wire Q = 50 → Nichrome wire Q ≈ 0.8 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) 1 kHz δ ≈ 2.1 mm 10 kHz δ ≈ 0.66 mm 100 kHz δ ≈ 0.21 mm 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): Use copper magnet wire Choose gauge based on skin depth at operating frequency Use largest gauge that fits the core/bobbin Consider Litz wire for frequencies >50 kHz For Current-Limited Applications: Use stainless steel or nichrome Calculate required resistance: R = V max /I limit Accept reduced Q factor as tradeoff 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 C parasitic 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: l wire = length of each wire d wire = wire diameter s = spacing between wires (≈ insulation thickness × 2) ε r = dielectric constant of insulation 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: Keep wires parallel: Twist them together before winding or use a jig Maintain tension: Even tension prevents gaps and loose spots Mark the wires: Use different colors or tag ends carefully Wind in layers: Complete one layer before starting next Insulate between layers: Add tape for voltage isolation Measuring Bifilar Parameters Measurement Configuration What It Tells You L A alone Measure A, B open Inductance of winding A L series-aid A end to B start, measure L A + L B + 2M L series-opp A end to B end, measure L A + L B - 2M C winding 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: Built-in capacitance may simplify circuit Excellent magnetic coupling Matched characteristics between windings Compact construction Disadvantages: Lower SRF due to high parasitic capacitance Difficult to adjust windings independently Insulation must handle full voltage difference More complex to wind correctly 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: Impedance is maximum (parallel resonance) Inductor is neither inductive nor capacitive Phase angle crosses through 0° 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 f op / 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: L dc = 10 mH, SRF = 100 kHz At 30 kHz: L eff = 10 / [1 - 0.09] = 11.0 mH (+10%) At 50 kHz: L eff = 10 / [1 - 0.25] = 13.3 mH (+33%) At 70 kHz: L eff = 10 / [1 - 0.49] = 19.6 mH (+96%) Minimizing Parasitic Capacitance Winding Techniques: Single-layer winding: Eliminates layer-to-layer capacitance Space-wound turns: Increases turn-to-turn distance Honeycomb/basket winding: Crosses turns to reduce adjacent voltage Bank winding: Winds in sections to reduce voltage across layers Progressive winding: Keeps voltage gradient low between adjacent turns Design Choices: Use fewer turns (requires higher permeability core) Use thinner insulation (but watch voltage ratings) Use air-core (eliminates winding-to-core capacitance) Choose toroidal cores (natural progressive winding) 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 Connect inductor to impedance analyzer Sweep frequency and plot |Z| SRF is where impedance peaks Method 2: Signal Generator + Oscilloscope Connect inductor in series with known resistor Drive with sine wave, sweep frequency Monitor voltage across inductor SRF is where voltage peaks (current minimum) Method 3: Resonance with Known Capacitor Measure inductance at low frequency Add known capacitor in parallel Find new resonant frequency 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: ρ = resistivity of wire material (Ω·m) l wire = total wire length (m) A wire = wire cross-sectional area (m²) 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: Scenario L DCR Q @ 10kHz V out (12V in) 22 AWG Copper 10 mH 5 Ω 126 1,508 V 26 AWG Copper 10 mH 13 Ω 48 580 V 22 AWG SS316 10 mH 220 Ω 2.9 34 V 22 AWG Nichrome 10 mH 320 Ω 2.0 24 V 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: Use the largest wire that fits: Fill the available winding area Choose copper: Unless current limiting is specifically needed Use higher permeability core: Fewer turns needed for same L Optimize core size: Larger cores have more room for thicker wire Consider parallel windings: Two parallel wires = half the DCR Practical Limits: Wire must fit on the core with proper insulation Multiple layers increase parasitic capacitance Very thick wire is hard to wind neatly Cost and availability of materials 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 R ds(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 22 AWG copper: 8m × 0.053 Ω/m = 0.42 Ω ✓ 26 AWG copper: 8m × 0.134 Ω/m = 1.07 Ω ✓ 30 AWG copper: 8m × 0.339 Ω/m = 2.71 Ω ✓ 22 AWG SS316: 8m × 2.3 Ω/m = 18.4 Ω ✗ (Q = 34) 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: Electrodes: Conductive plates or tubes (typically stainless steel) Electrolyte: Water (pure, tap, or with additives) Container: Housing to hold electrodes and water Connections: Electrical leads to the VIC circuit 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. Advantages: Simple to build, easy to calculate Disadvantages: Limited surface area, edge effects Typical spacing: 1-5 mm 2. Concentric Tubes Inner and outer cylinders with water in the annular gap. Advantages: Larger surface area, uniform field Disadvantages: Harder to machine precisely Typical gap: 0.5-3 mm 3. Tube Array Multiple concentric tube pairs in parallel. Advantages: Maximum surface area, scalable Disadvantages: Complex construction, uniform spacing critical Stanley Meyer's design: Used 9 tube pairs 4. Spiral/Wound Flat electrodes wound in a spiral with separator. Advantages: Very large surface area in compact volume Disadvantages: Complex to build, water flow issues 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 C wfc 1-100 nF Sets resonant frequency with L2 Solution Resistance R sol 10 Ω - 10 kΩ Affects Q factor Water Properties Matter The water used in the WFC significantly affects electrical behavior: Water Type Conductivity R sol 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: DC voltage applied Current flows continuously Higher conductivity = more efficient Faraday's law determines gas production VIC Approach: High-frequency pulsed/AC voltage Capacitive charging dominates Lower conductivity may be preferred Electric field stress is the focus 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 = ε₀ε r A / d For Water (ε r ≈ 80): C (pF) ≈ 708 × A(cm²) / d(mm) Example: 10 cm × 10 cm plates = 100 cm² 2 mm gap C = 708 × 100 / 2 = 35,400 pF = 35.4 nF Concentric Tube Electrodes Cylindrical geometry provides more surface area: Capacitance: C = 2πε₀ε r L / ln(r outer /r inner ) Simplified (for small gap relative to radius): C ≈ ε₀ε r × 2πr avg L / d Where d = r outer - r inner Example: Inner tube: 20 mm OD Outer tube: 22 mm ID Length: 100 mm Gap: 1 mm C ≈ 708 × π × 2.1 × 10 / 1 = 46.7 nF 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: V = 1000 V (from VIC magnification) d = 1 mm = 0.001 m E = 1000 / 0.001 = 1,000,000 V/m = 1 MV/m 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: Parallelism: Plates must be parallel for uniform field Concentricity: Tubes must be truly concentric Uniform gap: Variations cause hot spots and non-uniform current Insulating spacers: Use non-conductive materials (PTFE, ceramic) 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 ): Pure water at 20°C: ε r ≈ 80 Pure water at 25°C: ε r ≈ 78.5 Pure water at 100°C: ε r ≈ 55 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: Siemens per meter (S/m) Microsiemens per centimeter (µS/cm) - most common Millisiemens per centimeter (mS/cm) 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: Tap water: σ = 500 µS/cm = 0.05 S/m Electrode area: 100 cm² = 0.01 m² Gap: 2 mm = 0.002 m R sol = 0.002 / (0.05 × 0.01) = 4 Ω Effect on Q Factor Solution resistance directly impacts circuit Q: Q total = 2πfL / (R choke + R sol + R other ) Example Impact: Water Type R sol Q (if R choke =5Ω) Distilled (σ=5 µS/cm) ~400 Ω Q ≈ 1.5 Tap (σ=500 µS/cm) ~4 Ω Q ≈ 70 Electrolyte (σ=20000 µS/cm) ~0.1 Ω Q ≈ 125 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: TDS meters (approximate, assume NaCl) True conductivity meters (more accurate) Laboratory grade (calibrated, temperature compensated) DIY Measurement: Use known electrode geometry cell Measure AC resistance at 1 kHz (to avoid polarization) Calculate σ from geometry and resistance 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 = ε₀ε r A / d Quick Formula for Water: C (nF) = 0.0708 × A(cm²) / d(mm) Example: A = 50 cm², d = 1 mm C = 0.0708 × 50 / 1 = 3.54 nF Concentric Cylinders C = 2πε₀ε r L / 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: 9 concentric tube pairs Each pair: C ≈ 5 nF Total: C = 9 × 5 = 45 nF 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: A = 100 cm², d = 1 mm, c dl = 25 µF/cm² C geo = 7.08 nF C edl = 25 µF/cm² × 100 cm² = 2500 µF = 2.5 mF 1/C = 1/7.08nF + 2/2.5mF ≈ 1/7.08nF C total ≈ 7.08 nF (EDL negligible) 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 Connect WFC in series with known resistor R Apply step voltage Measure time to reach 63% of final voltage C = τ / R Method 3: Resonant Frequency Connect WFC with known inductor L Drive with variable frequency Find resonant peak 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: f₀ is set by the pulse generator (typically 1-50 kHz) C wfc is constrained by electrode geometry and water properties L₂ must be designed to complete the resonant match Matching Strategies Strategy 1: Design L₂ for Given WFC When WFC geometry is fixed (existing cell): Measure C wfc with LCR meter Choose target frequency f₀ Calculate required L₂: L₂ = 1 / (4π²f₀²C wfc ) Example: C wfc = 10 nF (measured) f₀ = 10 kHz (desired) L₂ = 1 / (4π² × 10⁴² × 10⁻⁸) = 25.3 mH Strategy 2: Design WFC for Given L₂ When using a pre-wound or available choke: Measure L₂ with LCR meter Choose target frequency f₀ Calculate required C wfc : C wfc = 1 / (4π²f₀²L₂) Design electrodes to achieve that capacitance 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: L₂ C wfc f₀ Z₀ 10 mH 25 nF 10 kHz 632 Ω 50 mH 5 nF 10 kHz 3162 Ω 100 mH 2.5 nF 10 kHz 6325 Ω 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 Connect oscilloscope across WFC Sweep generator frequency slowly Watch for voltage peak Note frequency of maximum amplitude Method 2: Phase Measurement Monitor current and voltage simultaneously At resonance, current and voltage are in phase (phase = 0°) Below resonance: capacitive (current leads) 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: Add/remove turns (large adjustment) Adjust core gap if gapped (medium) Use adjustable ferrite slug (fine) For C wfc Adjustment: Add parallel capacitor (increases C) Change water level (changes effective area) Adjust electrode spacing (if possible) For Frequency Adjustment: PLL feedback to track resonance Variable frequency oscillator Multiple operating modes Complete Matching Checklist ☐ Measure or calculate C wfc ☐ Measure or calculate L₂ ☐ Calculate expected f₀ = 1/(2π√(L₂C)) ☐ Verify f₀ is within driver frequency range ☐ Calculate Z₀ = √(L₂/C) ☐ Estimate R total (DCR + solution R) ☐ Calculate Q = Z₀/R ☐ Build circuit and measure actual resonance ☐ Fine-tune as needed ☐ Verify Q meets design goals 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. Core selection (ferrite, iron powder, air core) Wire gauge and material selection Bifilar winding support Multi-layer winding calculations 2. Water Profile Module Model the WFC as an electrical component with all relevant parameters. Electrode geometry (plates, tubes, arrays) Water conductivity effects Temperature compensation EDL and solution resistance 3. Circuit Profile Module Combine chokes and WFC into complete VIC circuits for analysis. Primary and secondary resonance Q factor and bandwidth Voltage magnification Ring-down characteristics 4. Simulation Module Visualize circuit behavior and optimize performance. Frequency response plots Time-domain waveforms Impedance analysis Sensitivity analysis Design Workflow The recommended workflow for using the calculator: Define Requirements: Target frequency, available components, constraints Design/Select Chokes: Use Choke Design module or enter measured values Configure Water Profile: Enter WFC geometry and water properties Create Circuit Profile: Combine components and select topology Run Simulation: Analyze resonance, Q, and system behavior Optimize: Adjust parameters to improve performance Build & Verify: Construct circuit and compare to predictions 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): Core dimensions and material properties Wire gauge, material, and turn count Electrode geometry and spacing Water conductivity and temperature Operating frequency or frequency range Calculator Provides (Outputs): Inductance (L), DCR, parasitic capacitance Self-resonant frequency (SRF) WFC capacitance and ESR Resonant frequency (f₀) Q factor, bandwidth, ring-down time Voltage magnification ratio Impedance characteristics Frequency response curves 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: Navigate to the application dashboard 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 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: Check manufacturer datasheet for Aₗ and μᵣ Measure physical dimensions with calipers For unknown cores, estimate μᵣ from material type 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 n layers — 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 r i mm Inner tube radius (cylindrical) Outer Radius r o 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: Use a TDS or conductivity meter Distilled water: 1-10 µS/cm Tap water: 200-800 µS/cm If unknown, 500 µS/cm is a reasonable tap water estimate 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 f op kHz Pulse generator frequency Input Voltage V in V Peak pulse voltage Duty Cycle D % Pulse on-time percentage Source Resistance R s Ω 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: Inductance (measured at low frequency) DC Resistance (measured with ohmmeter) Self-Resonant Frequency (if known) For WFC: Capacitance (measured with LCR meter) ESR or solution resistance 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 Primary f₀: Resonant frequency of L1-C1 tank Secondary f₀: Resonant frequency of L2-C wfc tank Match Status: How well primary and secondary are tuned Q Factor Metrics Primary Q: Q factor of primary circuit Secondary Q: Q factor of secondary circuit System Q: Effective Q of coupled system Performance Metrics Voltage Magnification: V out /V in at resonance Bandwidth: -3dB frequency range Ring-down Time: Time constant τ = 2L/R Ring-down Cycles: Oscillation cycles during decay 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: Input Voltage: The driving pulse waveform Primary Current: Current through L1 WFC Voltage: Voltage across the water cell WFC Current: Current through the cell What to Look For: Voltage build-up during resonance Ring-down oscillations after pulse ends Phase relationship between V and I Settling time and stability 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: Run initial simulation with your component values Check if resonant frequency matches your target Evaluate Q factor—is it sufficient for your goals? Look for warnings and address them Adjust parameters and re-simulate Compare before/after to verify improvements 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, R sol ) May require larger components Strategy 1: Maximizing Q Factor Q determines voltage magnification and selectivity. To maximize Q: Reduce Choke DCR: Use larger wire gauge (lower AWG number) Use copper instead of aluminum Minimize wire length (fewer turns with higher-μ core) Consider Litz wire for high frequencies Reduce Solution Resistance: Increase water conductivity slightly (add small amount of electrolyte) Increase electrode area Decrease electrode gap (but watch capacitance change) Ensure good electrode contact Increase L or Decrease C: Higher L/C ratio raises Z₀ = √(L/C) Q = Z₀/R, so higher Z₀ means higher Q Must maintain same f₀ = 1/(2π√LC) 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 Design or select choke for desired L Calculate required C: C = 1/(4π²f₀²L) If C wfc ≠ required C: Add parallel capacitor if C wfc is too low Modify electrode geometry if adjustment is large Approach B: Fixed C, Adjust L Measure or calculate WFC capacitance Calculate required L: L = 1/(4π²f₀²C) Design choke for that inductance Approach C: Adjust Both Start with practical component ranges Use calculator to explore L/C combinations Choose combination that also optimizes Q 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 ): Maximum voltage transfer at resonance Narrow combined response Requires precise tuning Slight Offset (5-10% difference): Broader frequency response More tolerant of drift Slightly reduced peak transfer Calculator Approach: Design secondary (L2 + WFC) first—this is usually more constrained Calculate secondary f₀ Select C1 to tune primary to match: C1 = 1/(4π²f₀²L1) Verify with simulation Strategy 4: Optimizing for Available Components When working with existing components: Step 1: Characterize What You Have Measure L of available chokes Measure C of your WFC Note DCR values Step 2: Calculate Natural Resonance f₀ = 1/(2π√LC) This is where your circuit wants to resonate. Step 3: Evaluate Performance Is f₀ in your driver's range? Is Q acceptable at this frequency? Are there SRF issues? Step 4: Adjust as Needed Add tuning capacitor if f₀ is too high Consider different choke if f₀ is way off Accept the natural f₀ if performance is good 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% (R sol 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 ☐ Define your primary optimization goal ☐ Identify fixed constraints (available components, frequency range) ☐ Calculate baseline performance ☐ Identify largest loss contributor (DCR vs R sol ) ☐ Make targeted improvements to dominant loss ☐ Verify SRF is >3× operating frequency ☐ Check that primary/secondary are reasonably matched ☐ Run simulation to verify improvements ☐ Consider sensitivity to variations ☐ Document final design parameters 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 C parasitic 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 C wfc 1-100 nF WFC capacitance, sets resonance with L R solution 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₀ V magnification 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: Q factor: 30-100 (good balance of gain vs. stability) f₀: Within your driver's frequency range Operating frequency: < 30% of SRF (preferably < 10%) Primary/Secondary f₀ match: Within 5-10% Bandwidth: Wide enough to accommodate drift Voltage magnification: As needed for your application ✗ Warning Signs: Q < 10: Very low—circuit barely resonates Q > 300: Extremely sharp—hard to tune, sensitive to drift f op > 0.5 × SRF: Operating too close to SRF DCR > Z₀/10: Resistance dominates, poor Q Primary/Secondary mismatch > 20%: Poor energy transfer Translating Results to Construction Wire Length and Turns The calculator provides wire length and turn count. When winding: Add 10-20% to wire length for lead connections and margins Count turns carefully —L varies as N², so turn count is critical Verify L after winding —actual may differ from calculated 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 C wfc ±15% Fringe effects, water purity variation R solution ±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: Actual L is higher than calculated Stray capacitance adding to C total WFC capacitance underestimated When Measured f₀ is Higher Than Calculated: Actual L is lower than calculated Core saturation reducing effective L WFC capacitance overestimated When Measured Q is Lower Than Calculated: Additional losses not accounted for (core loss, skin effect) Poor connections adding resistance Water conductivity different than assumed 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: ☐ f₀ is within driver frequency range ☐ f₀ is < 30% of SRF (ideally < 10%) ☐ Q is in acceptable range (typically 20-150) ☐ Voltage magnification won't exceed component ratings ☐ Wire gauge handles expected current ☐ Primary and secondary frequencies are matched ☐ No warning indicators are present ☐ Results are saved for reference 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: 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] ──→ V wfc ──→ [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] ──→ V ctrl │ ↑ │ │ │ │ │ │ └───────┼───────────────────────────┼────┘ │ │ │ ↓ V sense │ [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) = (4V pk /π)[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: a n = (2V pk /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: f drive = 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 L leak1 L₁(1-k) Primary leakage inductance L leak2 L₂(1-k) Secondary leakage inductance L m 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: Large mode splitting Efficient energy transfer between windings Built-in inter-winding capacitance Lower overall SRF due to capacitance Measuring Bifilar Coupling: Measure L series-aid (windings in series, same polarity) Measure L series-opp (windings in series, opposite polarity) Calculate: M = (L series-aid - L series-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: 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): C = 10 nF, V peak = 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²R DCR Use larger wire, copper Solution Resistance P = I²R sol 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: 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 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 η η = P wfc /P in 50-90% Voltage Multiplication η V out /V in (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: 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% R solution 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, R sol 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₀ / (f high - f low ) = 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 V in Measure output voltage V out at resonance Q ≈ V out /V in 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(H 2 O) 2.99 × 10 −26 kg Amplitude (gap/2) 0.794 mm ω = 2πf 2.376 × 10 6 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 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 V out = Q × V in 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 f d = √(f₀² - α²/(4π²)) Hz Ringdown Cycles (to 1%) N ≈ 0.733 × Q cycles 4. Capacitance Formulas Formula Equation Notes Parallel Plate C = ε₀ε r A/d ε₀ = 8.854×10⁻¹² F/m Concentric Cylinders C = 2πε₀ε r L / ln(r o /r i ) L = length Capacitors in Series 1/C total = 1/C₁ + 1/C₂ + ... Capacitors in Parallel C total = 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 A L Method L = A L × N² A L in nH/turn² Inductors in Series L total = 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 C H = ε₀ε r A/d d ≈ 0.3 nm Debye Length λ D ≈ 0.304/√c (nm) c in mol/L Total EDL (series) 1/C = 1/C Stern + 1/C diff 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 V C,N = 2N × V s Lossless Maximum voltage V max ≈ (4Q/π) × V s 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 k B 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): R₂₀ = 10 Ω at 20°C At 50°C: R = 10 × [1 + 0.00393(50-20)] = 10 × 1.118 = 11.18 Ω At 80°C: R = 10 × [1 + 0.00393(80-20)] = 10 × 1.236 = 12.36 Ω 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): AWG Max Current (A) AWG Max Current (A) 10 15 24 1.4 12 9.3 26 0.9 14 5.9 28 0.55 16 3.7 30 0.35 18 2.3 32 0.22 20 1.8 34 0.14 22 2.1 36 0.09 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: Frequency Skin Depth (mm) Max Useful Wire Dia. 1 kHz 2.1 mm ~4 mm (AWG 6) 10 kHz 0.66 mm ~1.3 mm (AWG 16) 50 kHz 0.30 mm ~0.6 mm (AWG 22) 100 kHz 0.21 mm ~0.4 mm (AWG 26) 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 μᵢ B sat (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: Want L = 10 mH = 10,000,000 nH Using FT-240-77 (Aₗ = 4820 nH/turn²) N = √(10,000,000 / 4820) = 45.6 turns Use 46 turns for L ≈ 10.2 mH Saturation Considerations Saturation Flux Density (B sat ): Material Type B sat (mT) MnZn Ferrite 400-500 NiZn Ferrite 250-350 Iron Powder 800-1000 MPP 750 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): Moderate L (1-50 mH typical) Moderate current handling Consider: FT-82-77, FT-114-77, T-106-26 For Secondary Choke (L2): May need higher L (10-100 mH) for high Q Lower current typically Consider: FT-140-77, FT-240-77 For High Frequency (>100 kHz): Use lower-μ materials to maintain SRF margin Consider: Iron powder -2 or -6 mix, NiZn ferrite Quick Reference: Turns Calculation Desired L FT-82-77 FT-240-77 T-106-26 1 mH 21 turns 14 turns 105 turns 5 mH 48 turns 32 turns 236 turns 10 mH 68 turns 46 turns 333 turns 25 mH 107 turns 72 turns 527 turns 50 mH 152 turns 102 turns 745 turns 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ω) n ], 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.