Abstract — COMSOL Coronary Artery Flow
Comsol Coronary Artery Flow: Biomedical Engineering, Medical Electronics & Bio-Signal Modeling Integrated Clinical Investigation: Bio-Sensors · Analog Front-Ends (AFE) · Diagnostic Waveforms · DSP Algorithms · ~13512 Words Chapter 1: Electrophysiological Foundations, Action Potentials & Clinical Significance Cardiovascular diseases are the leading cause of death worldwide. In the United States, coronary artery disease accounts for 1 in 7 cardiovascular deaths. Cells in the blood and blood vessels are known to respond to alterations of the mechanical forces of their environment, such as blood flow-induced shear stress and tensile strain caused by the dilation and contraction of the blood vessel itself. Computational models of the coronary arteries are used to estimate fluid shear stress and wall tensile strain that are implicated in the development of atherosclerosis. Here we describe the development of two Fluid-Structure Interaction studies of the left anterior descending portion of the coronary artery in COMSOL®. The first is a meso-scale FSI model of a 19mm section of the left coronary artery. The purpose of this study was to couple the macroscopic regime of blood flow in the coronary arteries with the microscopic regime of cellular interactions, in particular endothelial cell-platelet interactions. The 3D geometry of the coronary artery was established using a patient’s CTA data; and a fibrous plaque was inserted into the artery, creating a 70% reduction in vessel diameter, which is clinically classified as a “severe” stenosis. This FSI model was able to provide detailed information in flow patterns, wall shear stress, and wall circumferential strain near or within the plaque with a spatial resolution of 21µm and a temporal resolution of 0.05 s. The average element size in this model was 50µm, which was comparable to vascular endothelial cells (20-100µm or bigger) and circulating platelets (2-4µm), enabling us to observe their interactions driven by blood flow. The second model was a macro-scale coronary artery model with a length of 7.56 cm. This model was also established using a patient-specific geometry (from CTA data). It would be used to assess the spatio-temporal distribution of fluid shear stress and wall tensile strain along the vascular wall and identify mechanical conditions that would promote atherosclerotic plaque growth. The significance of this model include: the patient-specific geometry, large volume (much bigger than most of the reported coronary artery numerical models), transient blood flow, non-linear elastic wall property, and transient squeezing induced by the myocardial contraction.
Sensor Interfacing & Analog Front-End Conditioning
Both of these studies required the use of the CFD and Solid Mechanics Modules and the platelet interaction studies would be performed with the Particle Tracing Module. Since the artery was treated as a hyper-elastic material, the Non-Linear Structural Materials module was employed. The fluid-structure interaction problem will be fully coupled at the interface between the fluid domain, representing the blood volume, and the solid domain, representing the artery wall or plaque material as appropriate.
Bio-Signal Waveform Analysis & Simulation Findings
Both studies provided new perspectives on the combined effects of blood flow and vessel strain on the development of cardiovascular disease states. Recently, I’ve been working through a computational fluid dynamics study examining blood velocity behavior in the right coronary artery under different geometric conditions, including healthy flow, stenosis, stent placement, and combined disease–stent cases. This project required translating physiological concepts into simplified, solvable models. I constructed custom RCA geometries, defined governing equations and boundary conditions, as well as making necessary assumptions to simplify the Navier-Stokes equation to be solvable. Using COMSOL, I analyzed velocity fields and cross-sectional profiles to understand how localized geometric changes influence flow acceleration, symmetry, and downstream recovery. While this work is still being fine-tuned, the process itself was incredibly valuable. It allowed me to bring together skill sets I’ve developed both as a biomedical engineer and as a cardiac CT analyst: interpreting coronary anatomy, thinking critically about stenosis severity and location, and connecting geometry to fluid mechanics and transport phenomena. Through this comparative analysis, the study revealed clear differences in peak velocity magnitude, profile shape, and flow recovery across each case.
Scope & Clinical Applications
- Real-time multi-parameter vital signs telemetry (ECG, Heart Rate, SpO2, Temperature, Respiration)
- Low-noise Analog Front-End (AFE) conditioning with active Driven Right Leg (DRL) 50 Hz cancellation
- Digital signal processing in MATLAB: Pan-Tompkins QRS detection, FIR/IIR notch, and Continuous Wavelet Transforms
- Finite element tissue bio-electromagnetics and orthopedic mechanical simulations in COMSOL & ANSYS
- Complete Arduino firmware .INO sketches, circuit schematics, and university dissertation reports
Bio-Potential Sensors
ECG, EEG, EMG & PPG optical probes
Low-Noise AFE
INA128 / AD8232 instrumentation amps
Arduino Firmware
Real-time ADC interrupt sampling & BLE
MATLAB Bio-DSP
Pan-Tompkins QRS & wavelet transforms
COMSOL Multiphysics
Tissue electro-thermal & SAR simulation
IEC 60601-1 Safety
Galvanic optical & magnetic isolation
Biomedical Sensors, Electronics & Simulation Suites
Engineered using clinical-grade bio-transducers, low-noise analog amplifiers, and multi-physics simulation platforms.
| Subsystem / Software | Diagnostic & Engineering Role | Instrumentation Tier |
|---|---|---|
| Arduino / ESP32 MCU | High-speed SPI/I2C sensor polling, 500 Hz timer interrupt ADC sampling, BLE telemetry streaming | Embedded Hardware |
| MATLAB R2024b | Pan-Tompkins QRS detection, 50 Hz Butterworth notch filtering, wavelet decomposition, SpO2 ratio | Signal Analytics |
| COMSOL Multiphysics | Electromagnetic tissue absorption (SAR), bio-heat thermal transfer, acoustic transducer wave propagation | Tissue Physics |
| ANSYS Mechanical | Orthopedic bone-implant contact mechanics, prosthetic socket stress distribution, dental load FEA | Biomechanics FEA |
| AD8232 / ADS1298 AFE | Medical-grade 3-lead ECG front-end with built-in instrumentation amp, DRL feedback, and high CMRR | Clinical AFE IC |
Biomedical Acquisition & Multi-Physics Simulation Workflow
Comprehensive signal flow from human bio-potential electrodes and optical sensor probes to embedded Arduino firmware, MATLAB DSP, and COMSOL tissue modeling.
(ECG/EEG/PPG/Temp)
(INA128 / DRL Loop)
500Hz ADC Sampling
& QRS Detection
Tissue FEA Models
& EHR Dashboard
Verified Medical Waveforms, Diagnostic Images & Sensor Circuits
Advantages of the COMSOL Coronary Artery Flow Project
Why this biomedical engineering project delivers superior academic grades and provides high technical evaluation standing.
Clinical Diagnostic Precision
Combines high-CMRR analog front-ends with precision digital filtering to extract microvolt bio-potentials with high SNR (> 35 dB).
Rigorous IEC 60601-1 Patient Safety
Incorporates dual galvanic isolation barriers (optical/magnetic) ensuring patient auxiliary leakage currents remain strictly under 10 μA.
Multi-Sensor Vital Signs Fusion
Simultaneously captures ECG, PPG, SpO2, Core Temperature, and Respiration for comprehensive physiological health monitoring.
Real-Time Arduino Embedded Firmware
Features optimized C++ interrupt-driven ADC sampling routines ready for immediate deployment on Arduino, ESP32, and STM32.
Validated Pan-Tompkins QRS Algorithm
Executes derivative, squaring, moving-window integration, and dual-threshold peak detection with > 99% R-peak accuracy.
COMSOL & ANSYS Multi-Physics Coupling
Models complex anatomical tissue structures, bio-electromagnetic Specific Absorption Rates (SAR), and prosthetic stress fields.
Driven Right Leg (DRL) Noise Suppression
Actively senses and neutralizes 50/60 Hz powerline hum and common-mode displacement currents without loading patient electrodes.
DICOM Medical Imaging Pipeline
Processes clinical MRI, CT, and Ultrasound scans using MATLAB contrast enhancement and morphological lesion boundary segmentation.
Direct Telehealth & Hospital IoT Readiness
Streams validated vital signs to mobile dashboards and electronic health records over encrypted BLE and MQTT protocols.
Complete Turnkey Project Package
Delivered with full Arduino .INO sketches, circuit schematics, MATLAB .m algorithms, university-format report, and viva defense slides.
Sample Arduino Firmware, ECG Waves & Sensor Telemetry
Representative excerpts illustrating the embedded sensor acquisition loop, live physiological waveform traces, and clinical parameter schedules.
[*] Initializing Analog Front-End (AFE) Signal Acquisition... [✓] Driven Right Leg (DRL) Active Feedback: 50 Hz Hum Attenuated by 42 dB ──────────────────────────────────────────────────────────────────────── [10:14:02.102] TRACE ECG Lead-II: ---__/\_/\_/\___ (Baseline: 0.02 mV | Normal ST Segment) [10:14:02.115] QRS R-Peak Registered (Fiducial Trigger: 1.24 mV | R-R Interval: 832 ms) [10:14:02.128] VITALS Heart Rate: 72.1 BPM | SpO2: 98.5% | PR Interval: 156 ms [10:14:02.140] TRACE ECG Lead-II: ---__/\_/\_/\___ (Baseline Stable | SNR = 36.8 dB) ──────────────────────────────────────────────────────────────────────── [✓] Diagnostic Evaluation: Sinus Rhythm Verified. Zero Ectopic Ventricular Beats.
[*] Synchronizing vital signs telemetry stream across bio-sensor channels...
| Physiological Parameter | Primary Transducer | Measured Amplitude | Sampling Bandwidth | Clinical Status |
|---|---|---|---|---|
| Cardiac Electrocardiogram (ECG) | Ag/AgCl + AD8232 AFE | 1.24 mV (R-Peak) | 0.05 Hz - 150 Hz | NORMAL SINUS |
| Pulse Oximetry (SpO₂) | MAX30102 Dual Red/IR | 98.5% Saturation | 100 Hz Optical | OPTIMAL |
| Heart Rate (HR) | R-R Peak Derivative | 72.1 BPM | 30 - 220 BPM range | EUKINETIC |
| Core Body Temperature | NTC Medical Thermistor | 36.8 °C | 0.1 Hz Thermal | AFEBRILE |
[✓] Diagnostic Evaluation: Patient physiological parameters within verified clinical tolerance band.
Frequently Asked Questions
Project Deliverables & Support
Everything needed from project inception to final viva defense — delivered with expert biomedical support.
.INO Sketches
Schematics
Algorithms
Report
Slides
Models
Datasets
Q&A
Support
Extension