cfdpre 0.1__tar.gz

cfdpre-0.1/PKG-INFO

@@ -0,0 +1,28 @@
+Metadata-Version: 2.1
+Name: cfdpre
+Version: 0.1
+Summary: CFD PreProcessing Library
+Author: Pushkar Sheth
+Author-email: siglyserdev@gmail.com
+Description-Content-Type: text/markdown
+Requires-Dist: numpy
+Requires-Dist: CoolProp
+
+# CFD-Pre
+CFD PreProcessing Library. 
+
+Provides the following functionality:
+1. Calculate Boundary Layer Mesh dimensions.
+
+
+
+
+## Project Log
+January 2025:
+1. Created Library
+
+
+## Project Road Map:
+
+1. Documentation for existing functionality.
+2. Include example data within library.

cfdpre-0.1/README.md

@@ -0,0 +1,18 @@
+# CFD-Pre
+CFD PreProcessing Library. 
+
+Provides the following functionality:
+1. Calculate Boundary Layer Mesh dimensions.
+
+
+
+
+## Project Log
+January 2025:
+1. Created Library
+
+
+## Project Road Map:
+
+1. Documentation for existing functionality.
+2. Include example data within library.

cfdpre-0.1/cfdpre/__init__.py

@@ -0,0 +1,3 @@
+# __init__.py
+
+from .main import yhgrcalc

cfdpre-0.1/cfdpre/main.py

@@ -0,0 +1,99 @@
+# main.py
+
+
+import numpy as np
+import CoolProp.CoolProp as CP
+
+
+
+def calculate_f(r, N, delta99, yH):
+
+    """Calculates the function f(r) for root finding."""
+
+    return r**N - r * (delta99 / yH) + (delta99 / yH - 1)
+
+
+
+def calculate_df_dr(r, N, delta99, yH):
+
+    """Calculates the derivative of f(r) with respect to r."""
+
+    return N * r**(N-1) - (delta99 / yH)
+
+
+
+def newton_raphson(N, delta99, yH, initial_guess=5, tolerance=1e-12, max_iterations=1000):
+
+    """Finds the growth ratio using the Newton-Raphson method."""
+
+    r = initial_guess
+
+    for _ in range(max_iterations):
+
+        f_r = calculate_f(r, N, delta99, yH)
+
+        df_dr = calculate_df_dr(r, N, delta99, yH)
+
+        r_new = r - f_r / df_dr
+       
+
+        if np.abs(r_new - r) < tolerance:
+
+            return r_new
+
+        r = r_new
+
+    return r
+
+
+
+
+def yhgrcalc(fluid, temperature_c, pressure_bar, massflow_kgpersec, hydraulicdia_mm, target_yplus, num_layers):
+    
+    pressure_pa = pressure_bar * 1e5
+    hydraulicdia_m = hydraulicdia_mm/1000
+    temperature_k = temperature_c + 273.15
+    
+    dynvisc_pas = CP.PropsSI('V', 'T', temperature_k, 'P', pressure_pa, fluid)
+    dynvisc_nsm2 = dynvisc_pas
+    
+    thermal_conductivity_wpermk = CP.PropsSI('L', 'T', temperature_k, 'P', pressure_pa, fluid)
+    specific_heat_cp_jperkgk = CP.PropsSI('C', 'T', temperature_k, 'P', pressure_pa, fluid)
+    density_kgperm3 =  CP.PropsSI('D', 'T', temperature_k, 'P', pressure_pa, fluid)
+    kinevisc_m2s = dynvisc_nsm2 / density_kgperm3
+    volflowrate_m3persec = massflow_kgpersec / density_kgperm3
+    flowvelocity_mpersec = volflowrate_m3persec / ((np.pi/4) * np.power(hydraulicdia_m, 2))
+    reynolds = flowvelocity_mpersec * hydraulicdia_m / kinevisc_m2s
+    prandtl = specific_heat_cp_jperkgk * dynvisc_nsm2 / thermal_conductivity_wpermk
+    cf = np.power(((2 * np.log10(reynolds)- 0.65)), (-2.3)) # skin friction coefficient
+    tau_wall = 0.5 * density_kgperm3 * np.square(flowvelocity_mpersec) * cf # wall shear stress
+    u_tau = np.sqrt(tau_wall/density_kgperm3) #friction velocity
+    yp_m = (target_yplus * dynvisc_nsm2)  /(u_tau * density_kgperm3)
+    yh_m = yp_m * 2
+
+    
+    
+
+    if reynolds < 5e5:
+
+        delta99 = 4.91 * hydraulicdia_m / np.sqrt(reynolds)
+
+    else:
+
+        delta99 = 0.38 * hydraulicdia_m * reynolds**(-1/5)
+
+    
+    growth_ratio = newton_raphson(num_layers, delta99, yh_m)
+    final_layer_thickness_m = yh_m * growth_ratio**(num_layers - 1)
+    propdict = [density_kgperm3, specific_heat_cp_jperkgk, thermal_conductivity_wpermk, dynvisc_nsm2]
+    
+    return yh_m, growth_ratio, final_layer_thickness_m, reynolds, propdict
+
+
+
+    
+
+
+
+
+    

cfdpre-0.1/cfdpre.egg-info/PKG-INFO

@@ -0,0 +1,28 @@
+Metadata-Version: 2.1
+Name: cfdpre
+Version: 0.1
+Summary: CFD PreProcessing Library
+Author: Pushkar Sheth
+Author-email: siglyserdev@gmail.com
+Description-Content-Type: text/markdown
+Requires-Dist: numpy
+Requires-Dist: CoolProp
+
+# CFD-Pre
+CFD PreProcessing Library. 
+
+Provides the following functionality:
+1. Calculate Boundary Layer Mesh dimensions.
+
+
+
+
+## Project Log
+January 2025:
+1. Created Library
+
+
+## Project Road Map:
+
+1. Documentation for existing functionality.
+2. Include example data within library.

cfdpre-0.1/cfdpre.egg-info/SOURCES.txt

@@ -0,0 +1,9 @@
+README.md
+setup.py
+cfdpre/__init__.py
+cfdpre/main.py
+cfdpre.egg-info/PKG-INFO
+cfdpre.egg-info/SOURCES.txt
+cfdpre.egg-info/dependency_links.txt
+cfdpre.egg-info/requires.txt
+cfdpre.egg-info/top_level.txt

cfdpre-0.1/cfdpre.egg-info/dependency_links.txt

@@ -0,0 +1 @@
+

cfdpre-0.1/cfdpre.egg-info/requires.txt

@@ -0,0 +1,2 @@
+numpy
+CoolProp

cfdpre-0.1/cfdpre.egg-info/top_level.txt

@@ -0,0 +1 @@
+cfdpre

cfdpre-0.1/setup.cfg

@@ -0,0 +1,4 @@
+[egg_info]
+tag_build = 
+tag_date = 0
+

cfdpre-0.1/setup.py

@@ -0,0 +1,25 @@
+# setup.py
+
+from setuptools import setup, find_packages
+
+with open("README.md", "r") as f:
+    description = f.read()
+
+
+setup(
+    name='cfdpre',
+    version='0.1',
+    description='CFD PreProcessing Library',
+    author='Pushkar Sheth',
+    author_email='siglyserdev@gmail.com',
+    packages=find_packages(),
+    install_requires=[
+        'numpy',
+        'CoolProp'
+
+        
+    ],
+
+    long_description=description,
+    long_description_content_type="text/markdown",
+)