cfdpre 0.1__py3-none-any.whl

cfdpre/__init__.py

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

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.dist-info/METADATA

@@ -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.dist-info/RECORD

@@ -0,0 +1,6 @@
+cfdpre/__init__.py,sha256=kMXOwXEXmkNP2IXpfR3H9RoTL2hCYhFi-7aZvqR1A9g,41
+cfdpre/main.py,sha256=Ml9AWENuCoNarpx5tdrEARHmWYGICcxX203W_1UbiFE,2758
+cfdpre-0.1.dist-info/METADATA,sha256=sFrHbBFfdXaT5sOd0fP37_49WYVXLpSPHZ4bSQcBn3Q,538
+cfdpre-0.1.dist-info/WHEEL,sha256=GJ7t_kWBFywbagK5eo9IoUwLW6oyOeTKmQ-9iHFVNxQ,92
+cfdpre-0.1.dist-info/top_level.txt,sha256=z61rr-JLldcK6fkxASUG-u6oUgwcuX3rVp8S1U2737s,7
+cfdpre-0.1.dist-info/RECORD,,

cfdpre-0.1.dist-info/WHEEL

@@ -0,0 +1,5 @@
+Wheel-Version: 1.0
+Generator: bdist_wheel (0.43.0)
+Root-Is-Purelib: true
+Tag: py3-none-any
+

cfdpre-0.1.dist-info/top_level.txt

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