Design a laboratory experiment to validate the analytical predictions of flexural mode frequencies and coupling efficiencies for a fluid-filled viscoelastic shell, using a tissue-mimicking phantom.
The strongest reviewer criticism is the absence of experimental validation. Even a simple phantom experiment comparing analytical predictions with measured frequencies/amplitudes would dramatically strengthen the paper.
Construction:
- Cast a hollow oblate spheroidal silicone shell using 3D-printed molds
- Inner mold: oblate spheroid (a=18cm, c=12cm)
- Outer mold: slightly larger oblate spheroid (a=19cm, c=13cm) → h=10mm
- Material: Ecoflex 00-30 or Dragon Skin 10 (Shore 00-30 to 10A)
- Fill with degassed water through a sealed port
Material Calibration:
- Measure E of the silicone with a tensile test (Instron or similar)
- Expected range: 0.05-0.5 MPa (Ecoflex) or 0.1-2 MPa (Dragon Skin)
- Measure loss tangent with DMA (Dynamic Mechanical Analysis)
- Target: E ≈ 0.1 MPa to match relaxed abdominal wall
Advantages:
- E can be tuned by silicone formulation
- Geometry is well-controlled and measurable
- Water-filled → known fluid properties
- Can inject air bubbles for gas pocket experiments
- Can be constrained (partially clamped) to test BC effects
Disadvantages:
- Isotropic (real tissue is anisotropic)
- Uniform wall (real wall is multi-layer)
- No organs or mesentery
Construction:
- Double balloon technique: inflate balloon to spheroid shape in gelatin mold
- Multiple dipping to build up wall thickness
- Control E through gelatin concentration (10-20% w/v)
- Fill with water
Material properties:
- E = 0.01-0.1 MPa (tunable with concentration)
- More biologically realistic viscoelastic behavior
- Limited shelf life (days), temperature-sensitive
Construction:
- Large latex balloon (~15-20cm diameter when inflated)
- Fill with water to oblate spheroid shape (support in a bowl)
- Very thin wall (h ≈ 0.3-1mm)
Advantages: Fast, cheap, immediately available Disadvantages: Geometry not well controlled, very thin wall
-
Laser Doppler Vibrometer (LDV) — Gold standard
- Polytec PSV-500 or similar
- Scan the shell surface for mode shapes
- Frequency resolution: 0.01 Hz
- Non-contact → no mass loading
-
Accelerometers — Simpler
- Lightweight MEMS accelerometers (< 1g mass)
- Mount to shell surface with wax
- Multiple points for mode identification
- Concern: mass loading on soft shell
-
Hydrophone — Internal pressure measurement
- B&K 8103 or similar miniature hydrophone
- Insert through sealed port
- Measures internal pressure oscillation
Mechanical (WBV simulation):
- Electrodynamic shaker (Brüel & Kjær 4809 or similar)
- Phantom sits on rigid plate mounted to shaker
- Sine sweep 1-20 Hz, acceleration 0.1-2.0 m/s²
- Control with accelerometer on plate
Airborne acoustic:
- Large subwoofer (18" or 21") in sealed enclosure
- Or: dedicated infrasound source (Rotary Woofer, pneumatic)
- Calibration microphone (G.R.A.S. 40AZ or similar)
- SPL range: 90-130 dB at 2-20 Hz
- Test in anechoic chamber or large room (λ ≈ 50m at 7 Hz!)
Impact (for modal identification):
- Instrumented impact hammer with soft tip
- Impulse response → FFT → natural frequencies
- Quick screening method before detailed sweeps
- Fill phantom, seal, equilibrate temperature (22°C)
- Mount phantom on vibration-isolated stand
- Impact with instrumented hammer at 10 locations
- Record response at 3-5 locations (accelerometers or LDV)
- Compute FRFs → identify natural frequencies and mode shapes
- Compare with analytical predictions
- Mount phantom on shaker platform
- Sine sweep: 1-20 Hz, 0.1 m/s² to 2.0 m/s²
- Record:
- Base acceleration (control accelerometer)
- Shell surface displacement (LDV)
- Internal pressure (hydrophone)
- Compute transfer functions
- Compare with model predictions:
- Resonant frequency
- Transmissibility T(f)
- Internal pressure amplitude
- Place phantom in anechoic environment (or large room)
- Expose to infrasound: 2-20 Hz, 90-130 dB
- Record:
- Incident pressure (reference microphone)
- Shell surface displacement (LDV)
- Internal pressure (hydrophone)
- Measure coupling coefficient: ξ_measured / ξ_predicted
- KEY TEST: Verify (ka)^n coupling penalty
- Inject known volumes of air (1, 5, 10, 20, 50 mL)
- Repeat Phase 3 at each air volume
- Measure:
- Change in natural frequencies
- Change in airborne coupling
- Local displacement near gas pocket (LDV pointed at bubble)
- Compare with gas pocket resonance model predictions
- Partially clamp phantom (rigid supports simulating spine/pelvis)
- Repeat Phase 1 with different constraint configurations
- Measure frequency shift vs. constraint geometry
- Compare with BC multiplier estimates
- n=2 flexural mode at 4-10 Hz (depending on silicone E)
- Breathing mode at much higher frequency (>100 Hz for silicone)
- Mechanical coupling >> airborne coupling (by 10³ factor)
- Gas injection increases airborne coupling
- Non-linear effects at large amplitude (WBV)
- Mode coupling between flexural and sloshing modes
- The breathing mode might be measurable as internal pressure oscillation
- Gas pocket resonance might be at unexpected frequency (not Minnaert)
| Item | Cost Range |
|---|---|
| Silicone (Ecoflex 00-30, 2 gallons) | $80-120 |
| 3D-printed molds (resin, ~500g) | $50-100 |
| Accelerometers (3× MEMS) | $50-200 |
| Hydrophone (miniature) | $200-500 |
| Impact hammer (PCB 086C03) | ~$1500 (may borrow) |
| Shaker + amplifier (existing lab) | $0 (shared facility) |
| Subwoofer (18") | $200-500 |
| Calibration microphone | $500-1000 |
| Data acquisition system | $0 (existing lab) |
| Total (excluding shared equipment) | $600-2500 |
- Phantom fabrication: 1-2 weeks
- Material characterization: 1 week
- Modal testing: 1 week
- Mechanical excitation: 1 week
- Acoustic excitation: 1-2 weeks
- Gas pocket experiments: 1 week
- Data analysis and comparison: 2 weeks
- Total: 8-10 weeks
The minimum set needed for a compelling validation section:
- Measured f₂ vs. predicted f₂ — the core frequency validation
- Transmissibility T(f) curve — compare with model and ISO 2631
- Airborne displacement vs. SPL — confirm weak coupling
- Mechanical displacement vs. acceleration — confirm strong coupling
- Coupling ratio: ξ_mech / ξ_air — the key novel result
If time permits: 6. Mode shapes from scanning LDV 7. Gas pocket coupling enhancement 8. BC sensitivity