Results-driven researcher with strong problem-solving abilities and 16 peer-reviewed publications. Demonstrated innovation through multiple patents and the development of novel models and methods. Years of recorded experience in both industry R&D and academic research with strong collaborative skills across academic and industrial teams. A highly motivated self-starter who thrives in challenging multidisciplinary projects, combining technical excellence with effective communication.
I work on the numerical side of multiphase flow — building the classical computational fluid dynamics and lattice Boltzmann models that let a computer resolve what an interface does when two fluids meet, deform, and change phase. That focus has carried me from 16 peer-reviewed publications and several patents to several years of moving between academic research and industrial R&D, always on the same question: how to translate the physics into computational models that engineers can use to make better decisions.
That question is the spine of my PhD at NTNU (2021–2025) — Pore Scale Simulations for Wettability Description, with Carl Fredrik Berg supervising and Eirik Grude Flekkøy co-supervising — where I characterize wettability directly from 3D pore-scale images and reproduce it with multiphase simulations. It traces back to my MSc at the University of Tehran, where I started working on computational physics.
Four threads run through my work: rigorous multiphase fluid flow models, efficient numerical solvers, the physics of porous media and reservoir engineering, and translating all of it into engineering practice.
Phase-field modeling such as the Conservative Allen–Cahn and Cahn–Hilliard formulations, ternary fluid systems, and interface dynamics at high density and viscosity, solid-fluid interactions, and phase change.
Development of lattice Boltzmann and hybrid LB–finite difference solvers for different physics scenarios equipped with parallel computing such as OpenMP and MPI. Examples are droplet evaporation in binary and ternary systems, pool boiling, and condensation on liquid-impregnated surfaces.
Pore-scale simulation for wettability description, spatial characterization of wetting from 3D images, and reservoir-scale workflows with OPM Flow, Eclipse 100, ResInsight, and LBPM — bridging digital rock physics and field-scale modeling.
Applied simulation for industry: turbulent mixing tanks, gas turbine combustion, recuperator heat exchangers, ingot cooling, and mechanical/HVAC system design under three levels of professional engineering permits.
Twelve representative CFD and design studies.
A two-dimensional transient study of incompressible water spilling over an air-inflated rubber dam, aimed at optimizing the geometry and installation angle of a downstream flow deflector — maximizing the free-jet throw length while minimizing flow-induced structural vibration and wall erosion. Geometry is meshed in Gambit and solved in ANSYS Fluent; the air–water interface is captured with the Volume of Fluid (VOF) method with geometric reconstruction, turbulence is closed with standard k-ε and RNG models, and surface tension enters through the Continuum Surface Force formulation. The inflated dam's own shape is validated against Anwar's hydrostatic equation, and the predicted jet throw length is checked against both Chanson's analytical correlation and Sorouri-Nezhad et al.'s experimental correlation. A dimensionless cavitation parameter (σ) monitors structural risk.
A transient, three-dimensional study of the isothermal turbulent flow inside a cylindrical mixing tank (T = 300 mm, four wall baffles) driven by a standard six-blade Rushton turbine at 250 rpm. Geometry is built in SolidWorks, meshed with Gambit, and solved in ANSYS Fluent with the impeller rotation captured by the Sliding Mesh technique — a fully transient method that resolves the blade-baffle interaction directly, rather than approximating it the way steady Multiple-Reference-Frame models do. The single-phase field is governed by the incompressible Navier–Stokes equations closed with the standard k-ε turbulence model; the water–oil configuration uses an Eulerian–Eulerian two-fluid framework with inter-phase momentum exchange dominated by drag, closed via the Schiller–Naumann model. Before the parametric study begins, the CFD model is checked for grid independence and validated on three independent fronts against published experimental and numerical data.
A comparative runner-design study for two classical reaction turbines. The Kaplan runner is sized from head, flow rate, and specific-speed similarity (de Leva–de Siervo correlations), with blade angles refined via Carter's deviation rule and NACA compressor-cascade correlations across five streamlines, then meshed as a single-blade periodic domain and solved in ANSYS CFX (k-ε) to check the resulting pressure and blade-loading distributions. The Francis runner blade surface is derived two independent ways for cross-validation: a graphical hub-to-shroud streamline construction (inscribed circles, conformally mapped into 3D across roughly 60 angular stations per blade passage), and a numerical solution of the meridional stream-function equation — including the guide-vane exit swirl — on the meridional plane (finite-difference, SOR).
A well-placement optimization study for CO2 storage in a saline aquifer, using a synthetic case built on the Norne field's geology (Norwegian Sea). OPM Flow's dedicated CO2STORE module computes CO2 PVT properties — density, viscosity, enthalpy — from analytic correlations and converts them to black-oil equivalents internally, giving compositional-simulator accuracy at black-oil performance; results are visualized and quantified in ResInsight. Injection layouts with 1 to 5 wells at different grid positions are screened against a fixed injection schedule and a 350 bar per-well pressure limit, tracking how much CO2 stays trapped in the field versus migrates into the neighboring aquifers.
A numerical-methods study of steady two-dimensional heat conduction on a unit square with a sinusoidal Dirichlet condition on one edge, solved on both structured finite-volume grids and Gambit-generated unstructured triangular grids to compare mesh types under an identical discretization. Two explicit relaxation schemes (Gauss-Seidel and SOR) and two implicit Krylov-subspace solvers (GMRES and BiCGSTAB) are implemented in Fortran 90, with the implicit Jacobian built matrix-free via finite-difference perturbation and stored in Compressed Row Storage; Reverse Cuthill-McKee reordering shrinks the sparse bandwidth on the unstructured mesh, and SOR/ILU(0) preconditioning is benchmarked against the unpreconditioned baseline.
A first-principles design of a single-stage, single-suction centrifugal pump — impeller blades and volute casing — sized from head, flow rate, and speed using the correlations in Gülich's Centrifugal Pumps and cross-checked against Lobanov/Stepanov design charts. The meridional half-section and blade geometry are built up analytically (impeller eye and outlet diameters via two independent methods, blade angles from inlet/outlet velocity triangles, blade surface unrolled by Kaplan's conformal-mapping construction) and modeled in CATIA and EES, then the same input parameters are fed into CF Turbo to regenerate the impeller and volute independently — its meridional section, velocity distributions, and 3D geometry closely match the hand-derived design, validating the analytical procedure.
An exergy and thermoeconomic study of a three-evaporator Variable Refrigerant Flow (VRF) cooling cycle (R-410a; one compressor, one air-cooled condenser, three evaporators at 4°C/0°C/−4°C, five expansion valves, and two refrigerant separators), modeled state-by-state in EES and linked to a genetic algorithm in MATLAB. The GA searches the condenser temperature and evaporator cooling-capacity design space for the combination that jointly maximizes COP and minimizes total refrigerant mass flow; independently, an exergy balance quantifies irreversibility at every component, and an economic model — capital-cost correlations plus electricity cost — locates the condenser temperature that minimizes total annualized cost.
A two-dimensional, axisymmetric CFD study of the supersonic flow inside a converging-diverging nozzle, solved in ANSYS Fluent (density-based, k-ω SST) across the compressible Navier–Stokes, energy, and ideal-gas equations. Because a full CFD sweep of every candidate geometry is too slow to embed inside an optimizer, a two-layer neural network is trained on 125 Fluent cases to predict the thrust-to-weight ratio directly from three geometric parameters — inlet, throat, and outlet height — and a genetic algorithm then searches that trained surrogate for the geometry that maximizes thrust-to-weight, with the result verified by re-running it through Fluent.
A full computational fluid dynamics study of a Rayleigh step hydrodynamic bearing — solving the coupled Navier–Stokes and energy equations directly in ANSYS Fluent, rather than the thin-film Reynolds-equation approximation classical lubrication theory relies on. Oil viscosity is coupled to temperature through a custom UDF, and the solution is validated against two independent published results before sweeping the bearing's step position, film-height ratio, runner speed, and minimum film thickness to see how each reshapes the pressure and temperature fields — and ultimately, load capacity and friction.
A from-scratch Fortran solver for the classical 1D shock-tube (Riemann) problem, comparing three second-order finite-volume flux schemes for the compressible Euler equations — Roe's approximate Riemann solver, Steger-Warming flux-vector splitting, and AUSM — against the closed-form exact Riemann solution. Three canonical initial conditions probe different flow physics: a stationary contact discontinuity, a weak acoustic wave, and a strong normal shock, with pressure, temperature, density, velocity, and entropy tracked along the tube for each scheme.
A custom solver for the classical water-hammer problem: an upstream reservoir feeding two long pipelines through a surge tank and a control valve into a downstream reservoir. The 1D unsteady continuity and momentum equations for compressible pipe flow are converted into compatibility equations along their characteristic lines and marched forward on a fixed x-t grid using the Method of Characteristics — the standard technique for hydraulic transients in pipelines, turbines, pumps, and surge tanks. Four linear valve-closing durations (5, 10, 15, and 20 s) are compared to see how closing speed governs the resulting pressure surge and flow oscillation.
Design and construction-supervision work covering the full mechanical (MEP) scope of residential and office buildings totaling roughly 70,000 m², carried out under a licensed Design & Supervision Engineer permit. The scope spans fire protection (wet-pipe sprinkler and standpipe systems), sanitary sewage and roof stormwater drainage, domestic hot and cold water distribution, ducted HVAC, hydronic underfloor radiant heating, water heating radiators, and other systems — each system sized, drawn, and detailed in AutoCAD from first-principles heat-load, flow, and pressure calculations, then carried through construction-phase site supervision.
Sixteen peer-reviewed articles in journals including Journal of Computational Physics, Physical Review E, Physics of Fluids, and Transport in Porous Media.
Two patents granted in Iran and one in Norway for novel engineering inventions.
Bronze Medal for an invention at the 4th International Invention and Innovation Competition for IFIA INV members.
Best design out of 40 submissions at the 2nd Festival of Innovations in the Steel Industry (2019), Iranian Iron and Steel Association.
Ranked first in the Mechanical Engineering BSc program at Shahrood University of Technology, winning the top prize three consecutive years.
"A note on the summation relation in phase-field equation" selected as a Featured Article by Physics of Fluids (2023).
Member of Iran's National Elites Foundation, IFIA, and SPE; reviewer for Physics of Fluids and Journal of Computational Physics.
I'm open to research collaboration, consulting on multiphase-flow and CFD problems, and opportunities in computational physics and scientific software development. Send me an email or connect with me via the contact form.
Research profiles
