Back to Research

    Closed-Loop Differential Buoyancy Energy Conversion

    Classical fluid mechanics and electrodynamics — preprint

    Lord Sofiane AouadeneIndependent Researcher and Curious MindDecember 2025Preprint — Not Peer-Reviewed

    Abstract

    A closed-loop energy conversion system based on differential buoyancy in immiscible fluids is analyzed using classical mechanics, fluid dynamics, and electromagnetic induction. A sealed buoyant body circulates through vertical columns of distinct fluid density, performing net work per cycle proportional to the maintained density gradient. Electrical energy extraction occurs exclusively through induction in external coils, without mechanical shafts or fluid exchange. The analysis explicitly incorporates conduit-wall effects on hydrodynamic drag, fluid circulation losses (pumping losses), and the coupling between mechanical and electromagnetic forces. All force balances, terminal velocities, dissipation mechanisms, and energetic limits are explicitly derived. A worked numerical example demonstrates achievable operating regimes and clarifies scaling behavior. The system is shown to obey conservation laws, with the net work traceable to the external energy required to maintain the fluid density stratification.

    1. Introduction

    Buoyancy-induced motion in stratified fluids is a fundamental outcome of Archimedes' principle. Density gradients give rise to mechanically accessible energy when mass redistributes under gravity. The key scientific question addressed here is whether such motion, when constrained within a closed and sealed system, can be converted into electrical power in a manner consistent with conservation laws, including a full accounting of all dissipative forces and the energy required to maintain the driving gradient.

    This paper presents a fully explicit classical analysis of such a system. Emphasis is placed on algebraic transparency, a comprehensive energetic accounting, and the identification of all dissipation channels, including conduit effects and fluid circulation resistance, so that the analysis may be scrutinized, replicated, and experimentally tested.

    2. Assumptions and Notation

    All analysis is conducted under the following clearly stated assumptions:

    • Two immiscible, incompressible fluids of densities ρH\rho_{\mathrm{H}} (high) and ρL\rho_{\mathrm{L}} (low) remain stably stratified, maintained by an external energy input (e.g., heat or chemical separation).
    • A rigid body of fixed volume VV and density ρB\rho_{\mathrm{B}} satisfies ρL<ρB<ρH\rho_{\mathrm{L}} < \rho_{\mathrm{B}} < \rho_{\mathrm{H}}.
    • Gravity is uniform with magnitude gg.
    • Motion is confined to sealed conduits of internal area AconduitA_{\mathrm{conduit}}; no mass leaves or enters the system.
    • Electrical energy is extracted solely via electromagnetic induction (external coils surrounding the conduits).

    3. Work Available from Differential Buoyancy

    3.1 Local Force Balance

    For a body fully immersed in a fluid of density ρ\rho, Archimedes' principle gives the buoyant force ρVg\rho V g, while gravity contributes ρBVg\rho_{\mathrm{B}} V g. The net vertical force (positive upward) is:

    FBuoyancy=(ρρB)VgF_{\mathrm{Buoyancy}} = (\rho - \rho_{\mathrm{B}})\,V\,g

    3.2 Gross Work per Cycle

    Consider a closed trajectory where the body rises through the denser fluid column of height HH and descends through the lighter fluid column of equal height. The mechanical work performed during ascent and descent sums to the gross mechanical work available per cycle:

    Wgross=(ρHρB)VgH+(ρBρL)VgH=(ρHρL)VgHW_{\mathrm{gross}} = (\rho_{\mathrm{H}} - \rho_{\mathrm{B}})\,V\,g\,H + (\rho_{\mathrm{B}} - \rho_{\mathrm{L}})\,V\,g\,H = (\rho_{\mathrm{H}} - \rho_{\mathrm{L}})\,V\,g\,H

    This gross work represents the reduction in gravitational potential energy of the fluid that is cyclically displaced by the body, and is the upper limit for net work extraction. Importantly, WgrossW_{\mathrm{gross}} is not "free energy" — it is supplied by the maintenance of the density gradient (external energy input).

    4. Drag, Dissipation, and Terminal Velocity

    At steady motion, the net buoyant force is balanced by the sum of hydrodynamic drag FDF_{\mathrm{D}} and electromagnetic drag FEMF_{\mathrm{EM}}. We must explicitly include conduit effects that increase drag relative to free-fluid coefficients.

    4.1 Hydrodynamic Drag and Conduit Effects

    Due to motion within a conduit, the drag coefficient must be replaced by an effective drag coefficient CD,effC_{\mathrm{D,eff}} which accounts for the blockage ratio β=A/Aconduit\beta = A / A_{\mathrm{conduit}}. The commonly used quadratic drag model becomes:

    FD=12CD,effρAv2F_{\mathrm{D}} = \tfrac{1}{2}\,C_{\mathrm{D,eff}}\,\rho\,A\,v^{2}

    where vv is the instantaneous relative velocity between the body and surrounding fluid. Note that CD,effC_{\mathrm{D,eff}} increases with blockage ratio and with viscous confinement effects.

    4.2 Fluid Circulation Losses (Pumping Losses)

    Because the body displaces fluid within a closed conduit, the fluid must circulate around the body each cycle. Viscous shear and pressure losses in the return passages produce power loss PPumpingP_{\mathrm{Pumping}}. Over a cycle, these losses integrate to WpumpingW_{\mathrm{pumping}} and must be subtracted from the gross work to determine the net extractable energy.

    4.3 Electromagnetic Drag and Terminal Velocity

    Electrical power PEMP_{\mathrm{EM}} generated by coils surrounding the conduit appears as an effective electromagnetic drag force FEMF_{\mathrm{EM}} on the moving magnetic buoy. The steady-state (terminal) velocity satisfies the force balance:

    FBuoyancy=FD+FEMF_{\mathrm{Buoyancy}} = F_{\mathrm{D}} + F_{\mathrm{EM}}

    If both drags scale approximately like v2v^2, define an effective quadratic coefficient ΓEM\Gamma_{\mathrm{EM}} so that FEM12ΓEMv2F_{\mathrm{EM}} \approx \tfrac{1}{2} \Gamma_{\mathrm{EM}} v^2. Then the terminal speed is:

    vt=2FBuoyancyρACD,eff+ΓEMv_{\mathrm{t}} = \sqrt{\frac{2\,F_{\mathrm{Buoyancy}}}{\rho A C_{\mathrm{D,eff}} + \Gamma_{\mathrm{EM}}}}

    5. Numerical Worked Example

    Parameters used in this worked example (SI units):

    ParameterValue
    Ball radiusr=0.05 mr = 0.05\ \mathrm{m}
    Conduit radiusRconduit=0.06 mR_{\mathrm{conduit}} = 0.06\ \mathrm{m}
    Blockage ratioβ0.69\beta \approx 0.69
    Ball densityρB=850 kgm3\rho_{\mathrm{B}} = 850\ \mathrm{kg\,m^{-3}}
    Dense fluidρH=1000 kgm3\rho_{\mathrm{H}} = 1000\ \mathrm{kg\,m^{-3}}
    Light fluidρL=800 kgm3\rho_{\mathrm{L}} = 800\ \mathrm{kg\,m^{-3}}
    Effective drag coefficientCD,eff2.0C_{\mathrm{D,eff}} \approx 2.0
    Column heightH=3 mH = 3\ \mathrm{m}

    For a sphere, V=43πr3V = \tfrac{4}{3}\pi r^{3} and A=πr2A = \pi r^{2}. With r=0.05 mr=0.05\ \mathrm{m}:

    V5.24×104 m3,A7.85×103 m2V \approx 5.24\times 10^{-4}\ \mathrm{m^3}, \quad A \approx 7.85\times 10^{-3}\ \mathrm{m^2}

    Gross Work per cycle:

    Wgross=(ρHρL)VgH200×5.24×104×9.81×33.08 JW_{\mathrm{gross}} = (\rho_{\mathrm{H}} - \rho_{\mathrm{L}})\,V\,g\,H \approx 200 \times 5.24 \times 10^{-4} \times 9.81 \times 3 \approx 3.08\ \mathrm{J}

    Terminal Velocity (hydrodynamic-only estimate):

    FBuoyancy=150×5.24×104×9.810.77 NF_{\mathrm{Buoyancy}} = 150 \times 5.24 \times 10^{-4} \times 9.81 \approx 0.77\ \mathrm{N}
    v2×0.771000×7.85×103×2.00.35 ms1v_{\uparrow} \approx \sqrt{\frac{2 \times 0.77}{1000 \times 7.85 \times 10^{-3} \times 2.0}} \approx 0.35\ \mathrm{m\,s^{-1}}

    6. Full Energetic Consistency and Power Output

    6.1 Power Budget

    For a steady cycle rate ff (cycles per second) and per-cycle net extractable energy Wnet=WgrossWdissipationWpumpingW_{\mathrm{net}} = W_{\mathrm{gross}} - W_{\mathrm{dissipation}} - W_{\mathrm{pumping}}, the electrical power is PEM=fWnetP_{\mathrm{EM}} = f\,W_{\mathrm{net}}. Equivalently:

    PEM=PgrossPHydrodynamicPPumping=FEMvtP_{\mathrm{EM}} = P_{\mathrm{gross}} - P_{\mathrm{Hydrodynamic}} - P_{\mathrm{Pumping}} = F_{\mathrm{EM}} \cdot v_{\mathrm{t}}

    Operational efficiency depends critically on reducing pumping losses, minimizing hydrodynamic drag via shape optimization, and optimizing electromagnetic coupling so that FEMF_{\mathrm{EM}} removes a useful portion of mechanical energy without stalling the motion.

    6.2 Thermodynamic Consistency

    The system obeys the First Law of Thermodynamics. The energy extracted electrically appears as a reduction of the mechanical energy of the moving body and as heat dissipated by viscous processes. The density gradient is maintained by an external power source PExternalP_{\mathrm{External}}. Balance over time demands:

    dWEMdt+dQDissipateddt=PExternal\frac{dW_{\mathrm{EM}}}{dt} + \frac{dQ_{\mathrm{Dissipated}}}{dt} = P_{\mathrm{External}}

    7. Conclusion

    A closed-loop buoyancy-driven generator operating on differential fluid density has been analyzed with explicit inclusion of critical engineering and physics constraints: conduit-wall effects, fluid circulation losses, and coupled electromagnetic drag. The available work is explicitly linked to the density gradient, and realistic operating speeds and outputs have been quantified, showing a significant reduction in expected velocity due to the confined geometry. The final output is constrained by the efficiency of the hydrodynamic-to-electromagnetic energy transfer and the energy required to maintain the density stratification. The system remains conservative, scalable in principle, and experimentally testable.

    References

    1. L. D. Landau and E. M. Lifshitz, Fluid Mechanics, Pergamon (1987).
    2. G. K. Batchelor, An Introduction to Fluid Dynamics, Cambridge University Press (2000).
    3. J. D. Jackson, Classical Electrodynamics, Wiley (1999).
    4. H. Schlichting & K. Gersten, Boundary-Layer Theory, Springer (2017).
    5. S. F. Hoerner, Fluid-Dynamic Drag, Hoerner Fluid Dynamics (1965).
    6. S. M. V. S. C. Sankararaman, Hydrodynamic Resistance of Bodies in Confined Flow, Journal of Fluids Engineering (2010).
    7. A. Bejan, Convection Heat Transfer, Wiley (2013).

    Compiled December 2025