{"id":29032501,"url":"https://github.com/cfdby/laval-nozzle-analysis","last_synced_at":"2025-06-26T11:02:15.241Z","repository":{"id":295144701,"uuid":"989303079","full_name":"CFDBY/Laval-Nozzle-Analysis","owner":"CFDBY","description":"This project presents a 2D CFD simulation of compressible flow through a **De Laval nozzle**, performed in **ANSYS Fluent 2025 R1**. The aim was to observe **choked flow** behavior, check **mass conservation**, and compute **wall total force**. ","archived":false,"fork":false,"pushed_at":"2025-06-19T00:34:58.000Z","size":806,"stargazers_count":0,"open_issues_count":0,"forks_count":0,"subscribers_count":0,"default_branch":"main","last_synced_at":"2025-06-26T11:01:54.447Z","etag":null,"topics":["analysis","cfd","laval-nozzles","matlab"],"latest_commit_sha":null,"homepage":"","language":"MATLAB","has_issues":true,"has_wiki":null,"has_pages":null,"mirror_url":null,"source_name":null,"license":"mit","status":null,"scm":"git","pull_requests_enabled":true,"icon_url":"https://github.com/CFDBY.png","metadata":{"files":{"readme":"README.md","changelog":null,"contributing":null,"funding":null,"license":"LICENSE","code_of_conduct":null,"threat_model":null,"audit":null,"citation":null,"codeowners":null,"security":null,"support":null,"governance":null,"roadmap":null,"authors":null,"dei":null,"publiccode":null,"codemeta":null,"zenodo":null}},"created_at":"2025-05-23T20:52:57.000Z","updated_at":"2025-06-19T00:35:01.000Z","dependencies_parsed_at":"2025-05-26T09:45:43.836Z","dependency_job_id":null,"html_url":"https://github.com/CFDBY/Laval-Nozzle-Analysis","commit_stats":null,"previous_names":["cfdby/laval-nozzle-analysis"],"tags_count":0,"template":false,"template_full_name":null,"purl":"pkg:github/CFDBY/Laval-Nozzle-Analysis","repository_url":"https://repos.ecosyste.ms/api/v1/hosts/GitHub/repositories/CFDBY%2FLaval-Nozzle-Analysis","tags_url":"https://repos.ecosyste.ms/api/v1/hosts/GitHub/repositories/CFDBY%2FLaval-Nozzle-Analysis/tags","releases_url":"https://repos.ecosyste.ms/api/v1/hosts/GitHub/repositories/CFDBY%2FLaval-Nozzle-Analysis/releases","manifests_url":"https://repos.ecosyste.ms/api/v1/hosts/GitHub/repositories/CFDBY%2FLaval-Nozzle-Analysis/manifests","owner_url":"https://repos.ecosyste.ms/api/v1/hosts/GitHub/owners/CFDBY","download_url":"https://codeload.github.com/CFDBY/Laval-Nozzle-Analysis/tar.gz/refs/heads/main","sbom_url":"https://repos.ecosyste.ms/api/v1/hosts/GitHub/repositories/CFDBY%2FLaval-Nozzle-Analysis/sbom","host":{"name":"GitHub","url":"https://github.com","kind":"github","repositories_count":262053537,"owners_count":23251214,"icon_url":"https://github.com/github.png","version":null,"created_at":"2022-05-30T11:31:42.601Z","updated_at":"2022-07-04T15:15:14.044Z","host_url":"https://repos.ecosyste.ms/api/v1/hosts/GitHub","repositories_url":"https://repos.ecosyste.ms/api/v1/hosts/GitHub/repositories","repository_names_url":"https://repos.ecosyste.ms/api/v1/hosts/GitHub/repository_names","owners_url":"https://repos.ecosyste.ms/api/v1/hosts/GitHub/owners"}},"keywords":["analysis","cfd","laval-nozzles","matlab"],"created_at":"2025-06-26T11:00:40.479Z","updated_at":"2025-06-26T11:02:15.225Z","avatar_url":"https://github.com/CFDBY.png","language":"MATLAB","funding_links":[],"categories":[],"sub_categories":[],"readme":"# 🚀 Laval Nozzle CFD Analysis (ANSYS Fluent)\n\n**Project Update Notice ( June 2025 )**\nThis project is currently being revised due todiscrepancies identified in the initial simulation results. Specifically , some of the wall force and mass flow rate values were found to be unrealistic when compared with the expected physical behavior under the given boundary conditions.\n\n🔧 Update Explanation (Added Note):\nThe discrepancy was traced back to the use of a 2D planar model instead of an axisymmetric one. Since the nozzle geometry represents a rotationally symmetric body, solving it in a non-axisymmetric setup led to incorrect evaluation of parameters such as wall force and mass flow rate.\n\nThe updated simulation now uses an axisymmetric configuration, which accurately reflects the physical geometry and ensures more realistic results in accordance with the applied boundary conditions.\n\n\nThis project presents a 2D CFD simulation of compressible flow through a **De Laval nozzle**, performed in **ANSYS Fluent 2025 R1**. The aim was to observe **choked flow** behavior, check **mass conservation**, and compute **wall total force**.\n\n---\n\n## 📘 Description\n\nInitially, a **transient density-based simulation** was attempted using a very small time step (`1e-6`) and 10,000 iterations.\nHowever, due to long computation time and hardware limitations, this approach was aborted after only 3% progress.\n\nThe simulation was later **successfully completed as steady-state**, using a **pressure-based SIMPLE algorithm** with 2000 iterations. Convergence was achieved and post-processing results were collected.\n\n---\n\n## 📀 Geometry \u0026 Boundary Conditions\n\n* **Length**: 100 mm\n* **Height**: 20 mm\n* **Inlet Pressure**: 200,000 Pa\n* **Outlet Pressure**: 101,325 Pa\n\n![Dimensions](./LavalNozzle_CFD_Project/dimensions.png)\n![Geometry](./LavalNozzle_CFD_Project/geometry.png)\n\n---\n\n## 🤩 Mesh Preview\n\nA structured mesh of approximately 900k elements was used. Boundary layer refinement was applied with biasing near walls.\n\n![Mesh](./LavalNozzle_CFD_Project/mesh.png)\n\n---\n\n## ⚙️ Solver Settings\n\n| Parameter       | Value          |\n| --------------- | -------------- |\n| Solver Type     | Pressure-Based |\n| Simulation Type | Steady         |\n| Algorithm       | SIMPLE         |\n| Initialization  | Inlet-based    |\n| Iterations      | 2000           |\n\n![Scaled Residuals - Steady](./LavalNozzle_CFD_Project/scaled_residuals_steady.png)\n\n---\n\n## 📊 Mass Flow Rate Analysis (MATLAB)\n\nTo verify conservation of mass, the following MATLAB code was used:\n\n```matlab\n% Laval Nozzle - Inlet/Outlet Mass Flow Rate Difference\n\nmass_flow_inlet = 0.0964;   % kg/s (from Fluent)\nmass_flow_outlet = 0.0880;  % kg/s (from Fluent)\n\ndiff = abs(mass_flow_inlet - mass_flow_outlet);\npercentage_error = (diff / mass_flow_inlet) * 100;\n\nfprintf('Mass Flow Rate Difference: %.5f kg/s\\n', diff);\nfprintf('Percentage Error: %.5f%%\\n', percentage_error);\n```\n\n### ✅ Result:\n\n* **Mass Flow Rate Difference** = `0.00084 kg/s`\n* **Percentage Error** ≈ `9%`\n\n\u003e 🔍 This minor difference is acceptable in compressible flow simulations and confirms a balanced solution.\n\n\n---\n\n## 🧲 Force \u0026 Flow Visualizations\n\n| Output                | Value        |\n| --------------------- | ------------ |\n| Inlet Mass Flow Rate  | 0.0964 kg/s  |\n| Outlet Mass Flow Rate | 0.0880 kg/s  |\n\n### 🔽 Contours \u0026 Results:\n\n* **Mach Number**\n  ![Mach Number](./LavalNozzle_CFD_Project/mach_number_contour.png)\n\n* **Velocity Field**\n  ![Velocity](./LavalNozzle_CFD_Project/velocity_contour.png)\n\n* **Density**\n  ![Density](./LavalNozzle_CFD_Project/denisty_contour.png)\n\n* **Static Pressure**\n  ![Static Pressure](./LavalNozzle_CFD_Project/static_pressure.png)\n  ![Close-up Static Pressure](./LavalNozzle_CFD_Project/static_pressure_close.png)\n\n* **Static Temperature**\n  ![Static Temperature](./LavalNozzle_CFD_Project/static_temperature.png)\n\n  * **Updated Wall Force**\n ![Updated Wall Force](./LavalNozzle_CFD_Project/updated_wall_force.jpg)\n  * **Updated Inlet Mass Flow**\n  ![Updated Inlet Mass Flow](./LavalNozzle_CFD_Project/updated_inlet.jpg)\n\n  * **Updated Outlet Mass Flow**\n  ![Updated Outlet Mass Flow](./LavalNozzle_CFD_Project/updated_outlet.jpg)\n\n  \n\n\u003e 🧐 Created by Burak Yorukcu\n\u003e 💻 Software: ANSYS Fluent 2025 R1 \u0026 MATLAB\n\n\n---\n","project_url":"https://awesome.ecosyste.ms/api/v1/projects/github.com%2Fcfdby%2Flaval-nozzle-analysis","html_url":"https://awesome.ecosyste.ms/projects/github.com%2Fcfdby%2Flaval-nozzle-analysis","lists_url":"https://awesome.ecosyste.ms/api/v1/projects/github.com%2Fcfdby%2Flaval-nozzle-analysis/lists"}