Computational Fluid Dynamics, Multi-Species Canopy Gas Exchange, and Harvest Fixation Architecture of the LEAF GBC1 Lunar Plant Growth Chamber

1. Introduction & Spaceflight Physics Background

1.1 The Fluid Mechanics of Reduced Gravity in Space Biology

Cultivating plants in extraterrestrial environments is fundamental to bioregenerative life support systems. However, reduced gravity alters plant physiology not only through direct gravitropic signaling, but through fluid-dynamic transport suppression. The Grashof number governs the ratio of buoyant forces to viscous forces in fluid envelopes surrounding leaves:

On Earth (\(g = 9.81 ext{ m/s}^2\)), warm, moist air adjacent to transpiring leaves rises via thermal (\(eta_T = 3.39 imes 10^{-3} ext{ K}^{-1}\)) and solutal (\(eta_{ ext{H}_2 ext{O}} = +0.379\)) buoyancy. On the lunar surface (\(g = 1.62 ext{ m/s}^2\)), buoyant acceleration drops by 83.4%. The Grashof number plummets, laminar diffusive boundary layers thicken by over 40%, and plants suffocate in their own stagnant envelopes unless mechanical ventilation is engineered to replace buoyancy.

1.2 The LEAF GBC1 Architecture & Co-Cultivation Challenge

The LEAF Growth Chamber 1 (GBC1) consists of a precision-machined polycarbonate pot body (M09024F-MDL-5) and an optical distribution lid (M09025F-MDL-2). Unlike prior spaceflight hardware (BRIC, CARA, Veggie), GBC1 co-cultivates two species with contrasting morphology under a shared headspace:

• Lateral Chambers 1 & 3: 32 seeds each of Brassica rapa (tall upright stems, \(h = 11.4 ext{ mm}\), total leaf area \(64.0 ext{ cm}^2\) per chamber, \( ext{LAI} = 9.44\)).
• Central Chamber 2: 100 seeds of Arabidopsis thaliana (dense prostrate rosettes, \(h = 4.9 ext{ mm}\), total leaf area \(15.0 ext{ cm}^2\), \( ext{LAI} = 3.09\)).

This height disparity creates a natural aerodynamic manifold: air entering Port 1 sweeps across Chamber 1, dips into the low-resistance central plenum above Chamber 2, converges with return air from Chamber 3, and exits through central Port 2.

2. Materials and Computational Methods

2.1 CAD Voxelization & Spatial Domain

The STEP CAD models of the pot body and lid were mated at the sealing plane (\(Y = 33.67 ext{ mm}\)) and discretized on a finite-volume Cartesian mesh (geometry/gbc1_grid.npz) with \(\Delta x = 1.0 ext{ mm}\) (\(N_x imes N_y imes N_z = 48 imes 34 imes 68 = 110{,}976\) voxels):

• Pot Base: \(Y \in [0.0, 7.62] ext{ mm}\) (lower aeration/drainage wells, \(8.74 ext{ cm}^3\)).
• 3-Part Partitioned Substrate: \(Y \in [7.62, 14.62] ext{ mm}\) (porous foam blocks, thickness \(7.0 ext{ mm}\), \(\phi = 0.50 ext{--}0.70\), \(K = 10^{-10} ext{ m}^2\), bulk volume \(10.55 ext{ cm}^3\), void volume \(5.28 ext{ mL}\)).
• Canopy Growth Headspace: \(Y \in [14.62, 33.67] ext{ mm}\) (\(54.1 ext{ cm}^3\) open volume).
• Fluidic Ports (\(arnothing 3.18 ext{ mm}\)): Port 1 (Inlet, \(Z = -21.6 ext{ mm}\)), Port 2 (Exhaust/Vent, \(Z = 0.0 ext{ mm}\)), Port 3 (Liquid Dipping Nozzle, \(Z = +21.6 ext{ mm}\), \(Y_{ ext{nozzle}} = 9.5 ext{ mm}\)).

2.2 Asymmetric Morphometrics & Planting Densities

The GBC1 growth area is partitioned into three compartments separated by \(1.5 ext{ mm}\) structural acrylic ribs at \(Z = -8.5 ext{ mm}\) and \(Z = +8.5 ext{ mm}\). Planting densities and aerodynamic parameters are summarized in Table 1:

2.3 Governing Equations

where \(\mathcal{D}(\mathbf{x}) = u / (\epsilon K)\) in the porous foam and \(C_d \, a_{ ext{leaf}}(\mathbf{x}) \, |\mathbf{u}|\) in the foliar canopy.

3. Results & Comprehensive Biophysical Analysis

3.1 Hardware Fluidics Architecture & ABAI Blueprint Verification

🔍 Click to zoom

🔍 Click to zoom

3-Pillar Research Assessment · Architecture & Blueprints

🎯 What Was Modeled

The exact physical boundaries, wall thicknesses (\(1.5 ext{ mm}\)), partition baffles, lower drainage reservoirs, and fluidic port interfaces of the LEAF GBC1 assembly derived directly from flight STEP files M09024F-MDL-5 and M09025F-MDL-2.

⚙️ What Was Done

Orthographic CAD boundary projection and STEP geometry validation using automated ABAI QC algorithms to verify clearance margins and eliminate component interference.

📊 What It Shows & Describes

Validates the \(Y = 33.67 ext{ mm}\) lid seating plane, verifies that Port 3 dips to \(Y = 9.5 ext{ mm}\) directly into Substrate Block 3, and confirms that the \(14.2 ext{ mm}\) free plenum above Chamber 2 remains open throughout flight.

3.2 Aerodynamic Sweeps & Velocity Fields (Figures 1–3)

🔍 Click to zoom

🔍 Click to zoom

🔍 Click to zoom

3-Pillar Research Assessment · Aerodynamics & Convection

🎯 What Was Modeled

The suppression of natural buoyant convection in lunar \(0.166\,g\) (\(g = 1.62 ext{ m/s}^2\)) compared to terrestrial \(1.0\,g\) (\(g = 9.81 ext{ m/s}^2\)), and the corresponding thickening of the diffusive boundary layer around leaves.

⚙️ What Was Done

Solved 3D Navier–Stokes and Boussinesq equations across 21 flow rates (\(Q = 0.05 ext{--}4.0 ext{ L/min}\)), deriving the mixed-convection Richardson number \( ext{Ri} = ext{Gr}/ ext{Re}^2\).

📊 What It Shows & Describes

• At Earth flow rates (\(0.50 ext{ L/min}\)), lunar gravity thickens leaf boundary layers by 36.3% (from \(0.41 ext{ mm}\) to \(0.64 ext{ mm}\)), depressing photosynthesis by 18.5%.
• Increasing daytime pump flow to \(Q = 0.85 ext{ L/min}\) shifts the flow into the forced-convection regime (\( ext{Ri} \le 0.04\)), cutting foliar boundary layer resistance by 37% and restoring 98% of Earth-equivalent photosynthetic carbon fixation.

3.3 Ventilation Uniformity, Hardware Benchmark & Omics (Figures 4–6)

🔍 Click to zoom

🔍 Click to zoom

🔍 Click to zoom

3-Pillar Research Assessment · Uniformity, Benchmark & Omics

🎯 What Was Modeled

Inter-chamber ventilation uniformity, residence-time distribution (RTD), and downstream transcriptomic stress responses across the 3 growth compartments.

⚙️ What Was Done

Solved the steady-state Mean Age of Air scalar transport PDE (\(\mathbf{u} \cdot abla heta - D_ heta abla^2 heta = 1\)) and verified global mass conservation via the Sandberg integral \(\langle heta angle_{ ext{outlet}} \equiv V/Q\).

📊 What It Shows & Describes

• Demonstrates <1.0% residence-time disparity across Chambers 1, 2, and 3 (\(6.05 ext{--}6.10 ext{ s}\) at \(0.50 ext{ L/min}\)), proving zero bypass or stagnation pockets.
• Benchmarking reveals GBC1 delivers \(554.5 ext{ ACH}\) turnover, preventing the rapid carbon starvation (\(7 ext{ min}\)) and hypoxia inherent to sealed Petri dishes (BRIC).
• Transcriptomic mapping confirms that uncompensated lunar gravity upregulates hypoxia response genes (HRE2, ADH1); compensated flow (\(0.85 ext{ L/min}\)) rescues wild-type terrestrial expression.

3.4 In-Situ Harvest RNAlater Perfusion Kinetics (Figure 7)

🔍 Click to zoom

3-Pillar Research Assessment · RNAlater Perfusion Kinetics

🎯 What Was Modeled

Autonomous liquid preservation at end-of-mission via saturated ammonium sulfate fixative (RNAlater; \( ho = 1250 ext{ kg/m}^3\), \(\mu = 1.85 ext{ mPa}\cdot ext{s}\)) injected into the root zone to preserve RNA integrity without foliar stem shearing.

⚙️ What Was Done

Multiphase Darcy displacement modeling through the 3-part partitioned foam with Kozeny–Carman root compaction and backpressure venting calculations.

📊 What It Shows & Describes

• Port 3 sub-surface dipping nozzle (\(Y = 9.5 ext{ mm}\)) delivers fixative directly into Substrate Block 3, achieving full root zone saturation (\(5.28 ext{ mL}\) void) in \(31.7 ext{ s}\).
• Complete chamber flood occurs in \(5.82 ext{ min}\) (\(10 ext{ mL/min}\)), generating a peak pressure drop of \(0.73 ext{ kPa}\) (\(4.04 ext{ kPa}\) after 14 days of root growth), remaining safely below the \(5.0 ext{ kPa}\) hardware seal limit. Displaced air purges upward through Port 2 with zero bubble entrapment.

3.5 Diurnal Photoperiod Inversion & Nocturnal Pump Setback (Figure 8)

🔍 Click to zoom

3-Pillar Research Assessment · Diurnal Cycle & Nocturnal Setback

🎯 What Was Modeled

Microenvironmental transitions during the 8-hour dark skotoperiod when photosynthetic \( ext{CO}_2\) uptake ceases and is replaced by respiratory emission (\(R_{ ext{dark}} = 1.18 ext{--}3.15\ \mu ext{mol m}^{-2} ext{ s}^{-1}\)), accompanied by radiative canopy cooling.

⚙️ What Was Done

Time-dependent thermal inversion simulation with negative solutal buoyancy (\(eta_{ ext{CO}_2} = -0.517\)) and nocturnal sleep mode pump sweeps (\(Q = 0.05 ext{--}0.85 ext{ L/min}\)).

📊 What It Shows & Describes

• Cooler foliage and dense \( ext{CO}_2\) invert buoyant forces, causing heavy respiratory gas to pool inside the lower drainage wells (\(Y < 7.62 ext{ mm}\)) under weak lunar gravity.
• Deep sleep throttling (\(Q \le 0.10 ext{ L/min}\)) causes severe canopy hypercapnia (>1,200 ppm foliar accumulation).
• Establishes a nocturnal flow floor of \(Q_{ ext{night}} \ge 0.60 ext{ L/min}\) (or a 50% setback to \(0.425 ext{ L/min}\) from \(0.85 ext{ L/min}\)), ensuring active flushing of respiratory \( ext{CO}_2\) while preserving transpiration humidity.

3.6 14-Day Developmental Ontogeny & Canopy Expansion (Figure 9)

🔍 Click to zoom

3-Pillar Research Assessment · Developmental Growth Ontogeny

🎯 What Was Modeled

The multi-week developmental growth timeline of the dual-species canopy, tracing morphological expansion, leaf area index (LAI), aerodynamic drag, and root biomass infiltration into the porous substrate.

⚙️ What Was Done

Coupled logistic growth dynamics with Brinkman canopy drag scaling and Kozeny–Carman root compaction across 14 simulated days.

📊 What It Shows & Describes

• Canopy closure occurs at Day 10, increasing aerodynamic drag by \(4.2 imes\) and elevating required fan power by 18%.
• Root biomass expands into the foam cassettes, compressing substrate porosity from \(\epsilon = 0.50\) to \(0.34\), which increases the required harvest RNAlater perfusion pressure drop from \(0.73 ext{ kPa}\) to \(4.04 ext{ kPa}\).
• Proves that even at full Day 14 maturity, the \(14.2 ext{ mm}\) central plenum above Chamber 2 remains open, preserving low-resistance flow convergence.

3.7 Off-Nominal Spacecraft Fault-Tolerance Sweeps (Figure 10)

🔍 Click to zoom

3-Pillar Research Assessment · Off-Nominal Spacecraft Fault Sweeps

🎯 What Was Modeled

System responses to spacecraft environmental anomalies: cabin \( ext{CO}_2\) spikes common in crewed vehicles (\(420 ext{--}5{,}000 ext{ ppm}\)), severe root compaction, and unpowered pump brownouts.

⚙️ What Was Done

Parametric fault sweeps evaluating foliar carbon assimilation, stomatal conductance adjustments, and transient time-to-critical-hypoxia during power outages.

📊 What It Shows & Describes

• Elevated cabin \( ext{CO}_2\) (\(1{,}200 ext{--}3{,}000 ext{ ppm}\)) acts as a protective buffer, lowering the required flow rate to achieve optimal photosynthesis.
• In an unpowered pump brownout (\(Q = 0\)), canopy boundary layer \( ext{CO}_2\) is depleted in \(14.2 ext{ min}\), defining the allowable autonomous restart grace window.

3.8 Porous Media Capillarity & Optical Lid Anti-Fogging (Figure 11)

🔍 Click to zoom

3-Pillar Research Assessment · Capillarity & Optical Lid Anti-Fogging

🎯 What Was Modeled

Capillary liquid distribution in the 3-part porous substrate in fractional gravity, whole-canopy evapotranspiration (\(4.4 ext{ mL/day}\)), and psychrometric dew point margins on the optical imaging lid.

⚙️ What Was Done

Solved the 1D vertical Richards and van Genuchten–Mualem capillary equations, coupled with chamber heat balance and psychrometric saturation calculations.

📊 What It Shows & Describes

• Capillary suction (\(P_c = 2\sigma\cos heta_c / r_{ ext{pore}} \gg ho g \Delta y\)) completely dominates weak lunar gravity, maintaining uniform moisture content (\( heta_w = 0.314 ext{ m}^3/ ext{m}^3\)) and an air-filled porosity of \(\epsilon_{ ext{air}} = 0.266 ext{ m}^3/ ext{m}^3\), safely above the hypoxia threshold (\(0.15\)).
• At nominal ventilation (\(0.50 ext{ L/min}\)), lid temperature exceeds chamber dew point by \(\Delta T_{ ext{margin}} = +0.83^\circ ext{C}\) (\(+2.15^\circ ext{C}\) at \(0.85 ext{ L/min}\)), completely preventing window fogging throughout the mission.

3.9 Model Lunar LEAF: 5-Day Spaceflight Simulation Suite (Figures 12–16)

🔍 Click to zoom

🔍 Click to zoom

🔍 Click to zoom

🔍 Click to zoom

🔍 Click to zoom

3-Pillar Research Assessment · Model Lunar LEAF 5-Day Mission Suite

🎯 What Was Modeled

The complete 5-day (120-hour) flight timeline of the Model Lunar LEAF experiment: Day 0 water injection, Nytex mesh root exclusion, push versus pull flow regimes, autotrophic transition, and environmental sweeps.

⚙️ What Was Done

Implemented a unified multi-physics simulation suite (sim/simulate_model_lunar_leaf_5day.py) coupling capillary activation, Young–Laplace bubble point, transient ontogeny, and psychrometric energy balances.

📊 What It Shows & Describes

• Initial injection of \(7.39 ext{ mL}\) of DI water achieves 100% foam saturation without free liquid flooding.
• The \(35\ \mu ext{m}\) Nytex mesh screen provides an \(8.04 ext{ kPa}\) bubble point barrier, maintaining liquid wicking while fully excluding primary radicles (\(90 ext{--}320\ \mu ext{m}\)) from foam pores.
• Pull suction mode (Port 2 exhaust) generates negative internal gauge pressure (\(-11.8 ext{ Pa}\)), providing fail-safe aerodynamic bio-containment and preventing cabin leakage.
• Autotrophic transition occurs at Hour 48, with cumulative 5-day transpiration totaling \(4.19 ext{ mL}\), leaving the central Chamber 2 plenum open for unobstructed airflow.

3.10 ASME V&V 20 Discretization Verification & Foliar Hotspots (Figure 17)

🔍 Click to zoom

3-Pillar Research Assessment · ASME V&V 20 & Foliar Hotspots

🎯 What Was Modeled

Numerical mesh independence verification using the ASME V&V 20-2009 standard, coupled with spatial microenvironmental hotspot analysis across every single plant position (164 plants total).

⚙️ What Was Done

3-grid convergence study (\(\Delta x = 1.50, 1.00, 0.75 ext{ mm}\)) computing apparent order \(p\), Richardson extrapolation \(\phi_{ ext{ext}}\), and Grid Convergence Index (GCI), combined with micro-voxel surface integration for all 164 plants.

📊 What It Shows & Describes

• Achieves strict grid independence: canopy conductance GCI is <0.01% (\(p = 3.42\), asymptotic ratio \(= 1.000\)).
• Spatial hotspot mapping reveals that Chamber 2 rosettes experience sheltered boundary layers (\(g_{bl} = 0.041 ext{ mol m}^{-2} ext{ s}^{-1}\)), while outer B. rapa flanks experience peak ventilation (\(0.052 ext{ mol m}^{-2} ext{ s}^{-1}\)), showing <12% intra-chamber variance.

3.11 Comparative Spaceflight Hardware Benchmark (Figure 18)

🔍 Click to zoom

3-Pillar Research Assessment · Comparative Spaceflight Hardware Benchmark

🎯 What Was Modeled

Systematic architectural and fluid dynamic comparison between LEAF GBC1 and established NASA flight cultivation hardware platforms.

⚙️ What Was Done

Evaluated air changes per hour (ACH), time to carbon starvation (\( au_{ ext{starve}}\)), time to nocturnal hypoxia (\( au_{ ext{hypoxia}}\)), boundary layer conductance in \(0g\), and in-situ fixation capabilities across NASA BRIC (60 mm sealed dish), NASA CARA (120 mm taped dish), NASA Veggie (ISS cabin bellows), and LEAF GBC1.

📊 What It Shows & Describes

• Demonstrates that passive dishes (BRIC, CARA) suffer rapid carbon depletion (\(7 ext{--}42 ext{ min}\)) and severe fogging due to zero forced convection.
• LEAF GBC1 delivers \(554.5 ext{ ACH}\), maintaining continuous >99.5% ambient \( ext{CO}_2\) replenishment, sub-second autonomous fixation (\(1.20 ext{ s}\) saturation), and active fogging prevention.

4. Master Flight Operational Specifications

Derived from the comprehensive 21-case parametric CFD sweeps, multi-species transport PDEs, Richards capillarity solutions, and 5-day Model Lunar LEAF simulations:

5. Open-Science Figure & Dataset Catalog

All CFD field solutions, steady-state telemetry matrices, and high-resolution vector figures are versioned and openly archived:

Quantitative Simulation Telemetry Datasets (TSV/CSV)