From d7d3a9495ec762f35dda5430e2d3a65223947b9d Mon Sep 17 00:00:00 2001 From: Hursanay Date: Sun, 9 Aug 2026 11:45:58 +0200 Subject: [PATCH 1/2] fvPatchFields: Added turbulentSyntheticEddyInlet condition Velocity inlet generating synthetic turbulent fluctuations with a divergence-free synthetic eddy method, reproducing a target Reynolds stress tensor and integral length scale field. Target profiles default to bundled DNS channel flow data (Re_tau = 180, 395 and 590) mapped in wall units, or may be supplied as registered fields. Optional swirl, mass flow or mean velocity forcing as Function1, and self-correction of the inlet Reynolds stresses. Reference: Fyhn, E. H., & Fyhn, H. (2026). turbulentSyntheticEddyInlet: A self-correcting, divergence-free synthetic-eddy boundary condition for OpenFOAM. SoftwareX, 35, 102907. --- src/finiteVolume/Make/files | 2 + .../dns/dnsCatalog.H | 3722 +++++++++++++++++ .../eddy/syntheticEddy.C | 86 + .../eddy/syntheticEddy.H | 159 + .../eddy/syntheticEddyI.H | 66 + ...lentSyntheticEddyInletFvPatchVectorField.C | 2478 +++++++++++ ...lentSyntheticEddyInletFvPatchVectorField.H | 468 +++ 7 files changed, 6981 insertions(+) create mode 100644 src/finiteVolume/fields/fvPatchFields/derived/turbulentSyntheticEddyInlet/dns/dnsCatalog.H create mode 100644 src/finiteVolume/fields/fvPatchFields/derived/turbulentSyntheticEddyInlet/eddy/syntheticEddy.C create mode 100644 src/finiteVolume/fields/fvPatchFields/derived/turbulentSyntheticEddyInlet/eddy/syntheticEddy.H create mode 100644 src/finiteVolume/fields/fvPatchFields/derived/turbulentSyntheticEddyInlet/eddy/syntheticEddyI.H create mode 100644 src/finiteVolume/fields/fvPatchFields/derived/turbulentSyntheticEddyInlet/turbulentSyntheticEddyInletFvPatchVectorField.C create mode 100644 src/finiteVolume/fields/fvPatchFields/derived/turbulentSyntheticEddyInlet/turbulentSyntheticEddyInletFvPatchVectorField.H diff --git a/src/finiteVolume/Make/files b/src/finiteVolume/Make/files index 00c90aecdf..d1fa43591c 100644 --- a/src/finiteVolume/Make/files +++ b/src/finiteVolume/Make/files @@ -320,6 +320,8 @@ $(derivedFvPatchFields)/totalTemperature/totalTemperatureFvPatchScalarField.C $(derivedFvPatchFields)/translatingWallVelocity/translatingWallVelocityFvPatchVectorField.C $(derivedFvPatchFields)/turbulentInlet/turbulentInletFvPatchFields.C $(derivedFvPatchFields)/turbulentKineticEnergy/turbulentKineticEnergyFvPatchScalarField.C +$(derivedFvPatchFields)/turbulentSyntheticEddyInlet/turbulentSyntheticEddyInletFvPatchVectorField.C +$(derivedFvPatchFields)/turbulentSyntheticEddyInlet/eddy/syntheticEddy.C $(derivedFvPatchFields)/uniformDensityHydrostaticPressure/uniformDensityHydrostaticPressureFvPatchScalarField.C $(derivedFvPatchFields)/uniformFixedGradient/uniformFixedGradientFvPatchFields.C $(derivedFvPatchFields)/uniformFixedValue/uniformFixedValueFvPatchFields.C diff --git a/src/finiteVolume/fields/fvPatchFields/derived/turbulentSyntheticEddyInlet/dns/dnsCatalog.H b/src/finiteVolume/fields/fvPatchFields/derived/turbulentSyntheticEddyInlet/dns/dnsCatalog.H new file mode 100644 index 0000000000..ae080439d7 --- /dev/null +++ b/src/finiteVolume/fields/fvPatchFields/derived/turbulentSyntheticEddyInlet/dns/dnsCatalog.H @@ -0,0 +1,3722 @@ +/*---------------------------------------------------------------------------*\ + ========= | + \\ / F ield | OpenFOAM: The Open Source CFD Toolbox + \\ / O peration | Website: https://openfoam.org + \\ / A nd | Copyright (C) 2026 OpenFOAM Foundation + \\/ M anipulation | +------------------------------------------------------------------------------- +License + This file is part of OpenFOAM. + + OpenFOAM is free software: you can redistribute it and/or modify it + under the terms of the GNU General Public License as published by + the Free Software Foundation, either version 3 of the License, or + (at your option) any later version. + + OpenFOAM is distributed in the hope that it will be useful, but WITHOUT + ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or + FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License + for more details. + + You should have received a copy of the GNU General Public License + along with OpenFOAM. If not, see . + +Description + DNS channel flow profiles used by + turbulentSyntheticEddyInletFvPatchVectorField. + + AUTO-GENERATED from preprocessed DNS profiles. DO NOT EDIT BY HAND. + + One DNS profile per Re_tau in {180, 395, 590}, preprocessed to + dimensionless form: r=y/h, yplus, U=U/Ub, L=L/h, R=R/Utau^2. + +\*---------------------------------------------------------------------------*/ + +#ifndef dnsCatalog_H +#define dnsCatalog_H + +namespace Foam { + +struct DNSDataset +{ + int Re_tau; + int N; + const scalar* r; + const scalar* yplus; + const scalar* U; + const scalar* L; + const scalar* Lkeps; + const scalar* Ruu; + const scalar* Rvv; + const scalar* Rww; + const scalar* Ruv; + const scalar* Ruw; + const scalar* Rvw; +}; + +namespace dnsCatalogDetail { + +static const scalar Re180_r[65] = { + 0.000000000000000e+00, + 3.011800000000000e-04, + 1.204500000000000e-03, + 2.709500000000000e-03, + 4.815300000000000e-03, + 7.520500000000000e-03, + 1.082300000000000e-02, + 1.472200000000000e-02, + 1.921500000000000e-02, + 2.429800000000000e-02, + 2.996900000000000e-02, + 3.622400000000000e-02, + 4.306000000000000e-02, + 5.047200000000000e-02, + 5.845600000000000e-02, + 6.700700000000000e-02, + 7.611999999999999e-02, + 8.579000000000001e-02, + 9.601100000000000e-02, + 1.067800000000000e-01, + 1.180800000000000e-01, + 1.299100000000000e-01, + 1.422700000000000e-01, + 1.551500000000000e-01, + 1.685300000000000e-01, + 1.824200000000000e-01, + 1.967900000000000e-01, + 2.116500000000000e-01, + 2.269900000000000e-01, + 2.427900000000000e-01, + 2.590500000000000e-01, + 2.757500000000000e-01, + 2.928900000000000e-01, + 3.104600000000000e-01, + 3.284400000000000e-01, + 3.468300000000000e-01, + 3.656100000000000e-01, + 3.847700000000000e-01, + 4.043000000000000e-01, + 4.241900000000000e-01, + 4.444300000000000e-01, + 4.650000000000000e-01, + 4.859000000000000e-01, + 5.071000000000000e-01, + 5.286000000000000e-01, + 5.503900000000000e-01, + 5.724399999999999e-01, + 5.947600000000000e-01, + 6.173200000000000e-01, + 6.401000000000000e-01, + 6.631100000000000e-01, + 6.863200000000000e-01, + 7.097200000000000e-01, + 7.332900000000000e-01, + 7.570200000000000e-01, + 7.809000000000000e-01, + 8.049100000000000e-01, + 8.290400000000000e-01, + 8.532700000000000e-01, + 8.775900000000000e-01, + 9.019800000000000e-01, + 9.264400000000000e-01, + 9.509300000000001e-01, + 9.754600000000000e-01, + 1.000000000000000e+00 +}; + +static const scalar Re180_yplus[65] = { + 0.000000000000000e+00, + 5.364800000000000e-02, + 2.145600000000000e-01, + 4.826300000000000e-01, + 8.577100000000000e-01, + 1.339600000000000e+00, + 1.927900000000000e+00, + 2.622400000000000e+00, + 3.422600000000000e+00, + 4.328000000000000e+00, + 5.338100000000000e+00, + 6.452300000000000e+00, + 7.670000000000000e+00, + 8.990200000000000e+00, + 1.041200000000000e+01, + 1.193600000000000e+01, + 1.355900000000000e+01, + 1.528100000000000e+01, + 1.710200000000000e+01, + 1.901900000000000e+01, + 2.103300000000000e+01, + 2.314100000000000e+01, + 2.534200000000000e+01, + 2.763500000000000e+01, + 3.001900000000000e+01, + 3.249200000000000e+01, + 3.505300000000000e+01, + 3.770100000000000e+01, + 4.043200000000000e+01, + 4.324700000000000e+01, + 4.614300000000000e+01, + 4.911800000000000e+01, + 5.217100000000000e+01, + 5.530000000000000e+01, + 5.850300000000000e+01, + 6.177800000000000e+01, + 6.512300000000000e+01, + 6.853600000000000e+01, + 7.201600000000001e+01, + 7.555900000000000e+01, + 7.916400000000000e+01, + 8.282800000000000e+01, + 8.655000000000000e+01, + 9.032700000000000e+01, + 9.415700000000000e+01, + 9.803700000000001e+01, + 1.019700000000000e+02, + 1.059400000000000e+02, + 1.099600000000000e+02, + 1.140200000000000e+02, + 1.181200000000000e+02, + 1.222500000000000e+02, + 1.264200000000000e+02, + 1.306200000000000e+02, + 1.348400000000000e+02, + 1.391000000000000e+02, + 1.433700000000000e+02, + 1.476700000000000e+02, + 1.519900000000000e+02, + 1.563200000000000e+02, + 1.606600000000000e+02, + 1.650200000000000e+02, + 1.693800000000000e+02, + 1.737500000000000e+02, + 1.781200000000000e+02 +}; + +static const scalar Re180_U[65] = { + 0.000000000000000e+00, + 3.420790404637327e-03, + 1.367512605504171e-02, + 3.073730590285152e-02, + 5.456211395975811e-02, + 8.506855684754994e-02, + 1.221150555901406e-01, + 1.654241919235773e-01, + 2.144857713394407e-01, + 2.684579728988390e-01, + 3.260971974875006e-01, + 3.857899549725495e-01, + 4.457314322804528e-01, + 5.041487034944536e-01, + 5.595622045401364e-01, + 6.108941495184652e-01, + 6.575769470388240e-01, + 6.994128967292002e-01, + 7.365933215302509e-01, + 7.695008673232907e-01, + 7.984544056760815e-01, + 8.240279054105949e-01, + 8.466677867216978e-01, + 8.666929211771524e-01, + 8.845497289718250e-01, + 9.005570816734781e-01, + 9.150976251634260e-01, + 9.283626823823258e-01, + 9.406711248979395e-01, + 9.522142756509243e-01, + 9.630559089548328e-01, + 9.734511220638744e-01, + 9.833361406644968e-01, + 9.929022876973573e-01, + 1.002149563162456e+00, + 1.011141741373344e+00, + 1.019942596643576e+00, + 1.028552128973150e+00, + 1.036970338362068e+00, + 1.045197224810327e+00, + 1.053232788317930e+00, + 1.061013254571323e+00, + 1.068602397884059e+00, + 1.076000218256138e+00, + 1.083206715687559e+00, + 1.090158115864771e+00, + 1.096918193101326e+00, + 1.103423173083671e+00, + 1.109673055811806e+00, + 1.115667841285732e+00, + 1.121471303819001e+00, + 1.127019669098060e+00, + 1.132376711436462e+00, + 1.137414882207102e+00, + 1.142261730037084e+00, + 1.146725931985753e+00, + 1.150807488053106e+00, + 1.154506398239146e+00, + 1.157758888230318e+00, + 1.160564958026624e+00, + 1.162924607628063e+00, + 1.164774062721082e+00, + 1.166113323305683e+00, + 1.166878615068312e+00, + 1.167133712322521e+00 +}; + +static const scalar Re180_L[65] = { + 1.158991000000000e+00, + 1.158991000000000e+00, + 1.158991000000000e+00, + 1.158991000000000e+00, + 1.158991000000000e+00, + 1.158991000000000e+00, + 1.158991000000000e+00, + 1.158991000000000e+00, + 1.158991000000000e+00, + 1.158991000000000e+00, + 1.158991000000000e+00, + 1.166956559842209e+00, + 1.175676925049310e+00, + 1.185131327731755e+00, + 1.195299120092915e+00, + 1.195390630894309e+00, + 1.195488086295006e+00, + 1.195591486295006e+00, + 1.195700830894309e+00, + 1.195815939953542e+00, + 1.188013246000000e+00, + 1.179842254727273e+00, + 1.171310778545455e+00, + 1.162422693636364e+00, + 1.153181876181818e+00, + 1.138905702593660e+00, + 1.124119556292027e+00, + 1.108831108261287e+00, + 1.093065894299945e+00, + 1.080249621472053e+00, + 1.067064567791920e+00, + 1.053519838959601e+00, + 1.039619987825124e+00, + 1.025374120088544e+00, + 1.010799763213940e+00, + 1.005082938431752e+00, + 9.992439219748306e-01, + 9.932862050338819e-01, + 9.872115332042594e-01, + 9.810268888673765e-01, + 9.747463308474577e-01, + 9.796293871186440e-01, + 9.845897406779661e-01, + 9.896233933898305e-01, + 9.947276798305085e-01, + 9.998986018644067e-01, + 9.937723741883117e-01, + 9.875797542207793e-01, + 9.813091415584415e-01, + 9.749761347402597e-01, + 9.685807337662338e-01, + 9.621385371753246e-01, + 9.556339464285715e-01, + 9.490825600649351e-01, + 9.424999766233766e-01, + 9.358549990259740e-01, + 9.291944228896104e-01, + 9.224870511363636e-01, + 9.157484823051948e-01, + 9.089943149350649e-01, + 9.022245490259740e-01, + 8.954235860389610e-01, + 8.886226230519481e-01, + 8.818060615259740e-01, + 8.749895000000000e-01 +}; + +static const scalar Re180_Lkeps[65] = { + 1.506088265314988e-43, + 1.257506035466936e-07, + 8.076065230013989e-06, + 9.223566937681549e-05, + 5.176752453310932e-04, + 1.961290785886347e-03, + 5.783847104483866e-03, + 1.432945140783534e-02, + 3.114523669195107e-02, + 6.051439765446127e-02, + 1.051458436340214e-01, + 1.623680856789475e-01, + 2.234332433427494e-01, + 2.788669566376978e-01, + 3.238760368502057e-01, + 3.582793689550575e-01, + 3.838424042346824e-01, + 4.024576412687668e-01, + 4.154074318480789e-01, + 4.239347080633978e-01, + 4.295059405768852e-01, + 4.333956543310591e-01, + 4.372099195296222e-01, + 4.420850331271139e-01, + 4.485513506374851e-01, + 4.562116416286846e-01, + 4.652244063008562e-01, + 4.756960254998844e-01, + 4.860363267844045e-01, + 4.968506340056358e-01, + 5.084110068690831e-01, + 5.194214895772298e-01, + 5.305694634881002e-01, + 5.418390505068157e-01, + 5.543388640363788e-01, + 5.691131558333951e-01, + 5.828337935384047e-01, + 5.962948043849389e-01, + 6.110952488536966e-01, + 6.249750100075371e-01, + 6.377294111168368e-01, + 6.481030243412845e-01, + 6.587459263283010e-01, + 6.689905512346738e-01, + 6.763292058443233e-01, + 6.831240011412941e-01, + 6.912658185080417e-01, + 6.980691582508718e-01, + 7.010460442163724e-01, + 7.010215180795789e-01, + 7.001387348690200e-01, + 6.978691499051558e-01, + 6.924895964642808e-01, + 6.842102922485065e-01, + 6.747707730396310e-01, + 6.660075485866894e-01, + 6.568895153235088e-01, + 6.463741958262547e-01, + 6.352164894094358e-01, + 6.258428802319814e-01, + 6.189596352113642e-01, + 6.122914701138181e-01, + 6.069126156817437e-01, + 6.041080561502340e-01, + 6.029857482329249e-01 +}; + +static const scalar Re180_Ruu[65] = { + 2.535200000000000e-28, + 3.804600000000000e-04, + 6.073000000000000e-03, + 3.064100000000000e-02, + 9.644000000000000e-02, + 2.341700000000000e-01, + 4.812800000000000e-01, + 8.769000000000000e-01, + 1.450800000000000e+00, + 2.206500000000000e+00, + 3.107400000000000e+00, + 4.076200000000000e+00, + 5.013600000000000e+00, + 5.826100000000000e+00, + 6.449500000000000e+00, + 6.856600000000000e+00, + 7.053100000000000e+00, + 7.065500000000000e+00, + 6.930300000000000e+00, + 6.686400000000000e+00, + 6.369800000000000e+00, + 6.011500000000000e+00, + 5.637300000000000e+00, + 5.266000000000000e+00, + 4.909300000000000e+00, + 4.573000000000000e+00, + 4.261500000000000e+00, + 3.976400000000000e+00, + 3.715700000000000e+00, + 3.478900000000000e+00, + 3.265800000000000e+00, + 3.073800000000000e+00, + 2.899400000000000e+00, + 2.740900000000000e+00, + 2.597800000000000e+00, + 2.468700000000000e+00, + 2.349500000000000e+00, + 2.239000000000000e+00, + 2.137600000000000e+00, + 2.043800000000000e+00, + 1.956600000000000e+00, + 1.875100000000000e+00, + 1.799300000000000e+00, + 1.727800000000000e+00, + 1.657800000000000e+00, + 1.588900000000000e+00, + 1.520900000000000e+00, + 1.452400000000000e+00, + 1.383400000000000e+00, + 1.315200000000000e+00, + 1.249000000000000e+00, + 1.185200000000000e+00, + 1.122700000000000e+00, + 1.061600000000000e+00, + 1.002600000000000e+00, + 9.462500000000000e-01, + 8.931600000000000e-01, + 8.434800000000000e-01, + 7.978700000000000e-01, + 7.576100000000000e-01, + 7.236200000000000e-01, + 6.965400000000000e-01, + 6.772600000000000e-01, + 6.661500000000000e-01, + 6.625600000000000e-01 +}; + +static const scalar Re180_Rvv[65] = { + 3.349100000000000e-29, + 6.469500000000000e-10, + 1.582100000000000e-07, + 3.759800000000000e-06, + 3.383700000000000e-05, + 1.768300000000000e-04, + 6.500800000000000e-04, + 1.865000000000000e-03, + 4.448800000000000e-03, + 9.204700000000000e-03, + 1.702300000000000e-02, + 2.877700000000000e-02, + 4.522800000000000e-02, + 6.695200000000000e-02, + 9.428100000000000e-02, + 1.272400000000000e-01, + 1.655100000000000e-01, + 2.084200000000000e-01, + 2.549800000000000e-01, + 3.039600000000000e-01, + 3.539800000000000e-01, + 4.036200000000000e-01, + 4.516300000000000e-01, + 4.969800000000000e-01, + 5.387999999999999e-01, + 5.763900000000000e-01, + 6.091900000000000e-01, + 6.368600000000000e-01, + 6.593000000000000e-01, + 6.765200000000000e-01, + 6.886700000000000e-01, + 6.961300000000000e-01, + 6.992800000000000e-01, + 6.985400000000000e-01, + 6.944399999999999e-01, + 6.874400000000001e-01, + 6.780200000000000e-01, + 6.667500000000000e-01, + 6.542900000000000e-01, + 6.413100000000000e-01, + 6.281700000000000e-01, + 6.148700000000000e-01, + 6.012200000000000e-01, + 5.870300000000001e-01, + 5.722100000000000e-01, + 5.567800000000001e-01, + 5.409000000000000e-01, + 5.249300000000000e-01, + 5.091700000000000e-01, + 4.937100000000000e-01, + 4.786700000000000e-01, + 4.642100000000000e-01, + 4.505000000000000e-01, + 4.375900000000000e-01, + 4.255000000000000e-01, + 4.143600000000000e-01, + 4.044400000000000e-01, + 3.959500000000000e-01, + 3.890000000000000e-01, + 3.836000000000000e-01, + 3.796600000000000e-01, + 3.770200000000000e-01, + 3.754200000000000e-01, + 3.746000000000000e-01, + 3.743600000000000e-01 +}; + +static const scalar Re180_Rww[65] = { + 2.700300000000000e-28, + 1.110600000000000e-04, + 1.710000000000000e-03, + 8.129200000000000e-03, + 2.357500000000000e-02, + 5.169600000000000e-02, + 9.441500000000000e-02, + 1.513900000000000e-01, + 2.202900000000000e-01, + 2.976300000000000e-01, + 3.799900000000000e-01, + 4.646800000000000e-01, + 5.499000000000001e-01, + 6.342900000000000e-01, + 7.164700000000001e-01, + 7.949400000000000e-01, + 8.682900000000000e-01, + 9.351600000000000e-01, + 9.941700000000000e-01, + 1.044300000000000e+00, + 1.085300000000000e+00, + 1.117600000000000e+00, + 1.142100000000000e+00, + 1.159700000000000e+00, + 1.171900000000000e+00, + 1.179500000000000e+00, + 1.182200000000000e+00, + 1.179700000000000e+00, + 1.172000000000000e+00, + 1.159300000000000e+00, + 1.142000000000000e+00, + 1.121200000000000e+00, + 1.098800000000000e+00, + 1.075600000000000e+00, + 1.051800000000000e+00, + 1.028300000000000e+00, + 1.005400000000000e+00, + 9.821600000000000e-01, + 9.571400000000000e-01, + 9.296200000000000e-01, + 9.006000000000000e-01, + 8.713700000000000e-01, + 8.414300000000000e-01, + 8.097200000000000e-01, + 7.762700000000000e-01, + 7.425000000000000e-01, + 7.096700000000000e-01, + 6.783200000000000e-01, + 6.478900000000000e-01, + 6.178000000000000e-01, + 5.880100000000000e-01, + 5.584200000000000e-01, + 5.292100000000000e-01, + 5.008000000000000e-01, + 4.739200000000000e-01, + 4.493900000000000e-01, + 4.274600000000000e-01, + 4.082000000000000e-01, + 3.918600000000000e-01, + 3.785100000000000e-01, + 3.677200000000000e-01, + 3.592100000000000e-01, + 3.528800000000000e-01, + 3.487600000000000e-01, + 3.472700000000000e-01 +}; + +static const scalar Re180_Ruv[65] = { + -5.800900000000000e-31, + -1.112000000000000e-07, + -7.190300000000000e-06, + -8.323700000000000e-05, + -4.769600000000000e-04, + -1.853800000000000e-03, + -5.597400000000000e-03, + -1.405100000000000e-02, + -3.043700000000000e-02, + -5.818700000000000e-02, + -9.974600000000000e-02, + -1.553700000000000e-01, + -2.226500000000000e-01, + -2.970200000000000e-01, + -3.730400000000000e-01, + -4.456600000000000e-01, + -5.111500000000000e-01, + -5.673800000000000e-01, + -6.135900000000000e-01, + -6.500800000000000e-01, + -6.777800000000000e-01, + -6.978600000000000e-01, + -7.114800000000000e-01, + -7.196100000000000e-01, + -7.230799999999999e-01, + -7.226100000000000e-01, + -7.189400000000000e-01, + -7.127200000000000e-01, + -7.043700000000001e-01, + -6.943400000000000e-01, + -6.830100000000000e-01, + -6.704700000000000e-01, + -6.567700000000000e-01, + -6.419800000000000e-01, + -6.262799999999999e-01, + -6.098100000000000e-01, + -5.925200000000000e-01, + -5.745700000000000e-01, + -5.562000000000000e-01, + -5.374700000000000e-01, + -5.185400000000000e-01, + -4.994700000000000e-01, + -4.802000000000000e-01, + -4.606200000000000e-01, + -4.406000000000000e-01, + -4.201500000000000e-01, + -3.993400000000000e-01, + -3.781800000000000e-01, + -3.566400000000000e-01, + -3.348300000000000e-01, + -3.128700000000000e-01, + -2.906600000000000e-01, + -2.681700000000000e-01, + -2.455700000000000e-01, + -2.229900000000000e-01, + -2.006600000000000e-01, + -1.786500000000000e-01, + -1.567600000000000e-01, + -1.347000000000000e-01, + -1.124100000000000e-01, + -8.999300000000000e-02, + -6.758200000000000e-02, + -4.513100000000000e-02, + -2.261200000000000e-02, + 0.000000000000000e+00 +}; + +static const scalar Re180_Ruw[65] = { + -3.817200000000000e-31, + -9.711800000000000e-07, + -1.533400000000000e-05, + -7.582900000000000e-05, + -2.312300000000000e-04, + -5.361000000000000e-04, + -1.034500000000000e-03, + -1.742500000000000e-03, + -2.633900000000000e-03, + -3.622100000000000e-03, + -4.531300000000000e-03, + -5.128600000000000e-03, + -5.236100000000000e-03, + -4.769700000000000e-03, + -3.666900000000000e-03, + -1.877100000000000e-03, + 6.070400000000000e-04, + 3.701300000000000e-03, + 6.995300000000000e-03, + 9.703300000000000e-03, + 1.113600000000000e-02, + 1.110300000000000e-02, + 1.033200000000000e-02, + 1.006400000000000e-02, + 1.062500000000000e-02, + 1.101200000000000e-02, + 9.882600000000000e-03, + 6.930900000000000e-03, + 3.259500000000000e-03, + -2.774600000000000e-04, + -3.145500000000000e-03, + -4.543400000000000e-03, + -4.096800000000000e-03, + -1.612700000000000e-03, + 2.639400000000000e-03, + 7.584700000000000e-03, + 1.163500000000000e-02, + 1.411300000000000e-02, + 1.510600000000000e-02, + 1.511800000000000e-02, + 1.544000000000000e-02, + 1.709000000000000e-02, + 2.001200000000000e-02, + 2.305600000000000e-02, + 2.495600000000000e-02, + 2.541400000000000e-02, + 2.478100000000000e-02, + 2.361600000000000e-02, + 2.251800000000000e-02, + 2.163600000000000e-02, + 2.069200000000000e-02, + 1.989900000000000e-02, + 1.970600000000000e-02, + 1.988400000000000e-02, + 2.009900000000000e-02, + 2.031700000000000e-02, + 2.050800000000000e-02, + 2.054200000000000e-02, + 2.036800000000000e-02, + 2.028300000000000e-02, + 2.023900000000000e-02, + 2.025200000000000e-02, + 2.041400000000000e-02, + 2.065300000000000e-02, + 2.078100000000000e-02 +}; + +static const scalar Re180_Rvw[65] = { + 8.847700000000000e-32, + -3.356300000000000e-10, + -2.148200000000000e-08, + -2.461900000000000e-07, + -1.400400000000000e-06, + -5.398200000000000e-06, + -1.595300000000000e-05, + -3.793400000000000e-05, + -7.352000000000000e-05, + -1.157400000000000e-04, + -1.463500000000000e-04, + -1.441700000000000e-04, + -9.835300000000000e-05, + -1.613400000000000e-05, + 8.220000000000001e-05, + 1.787700000000000e-04, + 2.636300000000000e-04, + 3.356000000000000e-04, + 4.038100000000000e-04, + 4.748300000000000e-04, + 5.603400000000000e-04, + 6.914400000000000e-04, + 8.623000000000000e-04, + 1.022700000000000e-03, + 1.143100000000000e-03, + 1.234200000000000e-03, + 1.271300000000000e-03, + 1.189800000000000e-03, + 1.014000000000000e-03, + 8.525000000000000e-04, + 7.024300000000000e-04, + 5.223300000000000e-04, + 3.818100000000000e-04, + 4.527500000000000e-04, + 7.665699999999999e-04, + 1.281700000000000e-03, + 2.032100000000000e-03, + 2.922700000000000e-03, + 3.823600000000000e-03, + 4.549200000000000e-03, + 4.970600000000000e-03, + 5.015000000000000e-03, + 4.794200000000000e-03, + 4.422100000000000e-03, + 3.891900000000000e-03, + 3.314100000000000e-03, + 2.694100000000000e-03, + 1.979800000000000e-03, + 1.206000000000000e-03, + 3.854400000000000e-04, + -3.453300000000000e-04, + -9.595100000000000e-04, + -1.557900000000000e-03, + -2.141600000000000e-03, + -2.487100000000000e-03, + -2.556600000000000e-03, + -2.476400000000000e-03, + -2.489200000000000e-03, + -2.695300000000000e-03, + -3.022900000000000e-03, + -3.191600000000000e-03, + -2.976800000000000e-03, + -2.327100000000000e-03, + -1.282600000000000e-03, + 0.000000000000000e+00 +}; + +static const scalar Re395_r[129] = { + 0.000000000000000e+00, + 7.529800000000000e-05, + 3.011800000000000e-04, + 6.776200000000000e-04, + 1.204500000000000e-03, + 1.881900000000000e-03, + 2.709500000000000e-03, + 3.687400000000000e-03, + 4.815300000000000e-03, + 6.093000000000000e-03, + 7.520500000000000e-03, + 9.097400000000000e-03, + 1.082300000000000e-02, + 1.269900000000000e-02, + 1.472200000000000e-02, + 1.689500000000000e-02, + 1.921500000000000e-02, + 2.168300000000000e-02, + 2.429800000000000e-02, + 2.706000000000000e-02, + 2.996900000000000e-02, + 3.302400000000000e-02, + 3.622400000000000e-02, + 3.956900000000000e-02, + 4.306000000000000e-02, + 4.669400000000000e-02, + 5.047200000000000e-02, + 5.439300000000000e-02, + 5.845600000000000e-02, + 6.266099999999999e-02, + 6.700700000000000e-02, + 7.149400000000000e-02, + 7.611999999999999e-02, + 8.088600000000000e-02, + 8.579000000000001e-02, + 9.083200000000000e-02, + 9.601100000000000e-02, + 1.013300000000000e-01, + 1.067800000000000e-01, + 1.123600000000000e-01, + 1.180800000000000e-01, + 1.239300000000000e-01, + 1.299100000000000e-01, + 1.360300000000000e-01, + 1.422700000000000e-01, + 1.486400000000000e-01, + 1.551500000000000e-01, + 1.617800000000000e-01, + 1.685300000000000e-01, + 1.754100000000000e-01, + 1.824200000000000e-01, + 1.895400000000000e-01, + 1.967900000000000e-01, + 2.041600000000000e-01, + 2.116500000000000e-01, + 2.192600000000000e-01, + 2.269900000000000e-01, + 2.348300000000000e-01, + 2.427900000000000e-01, + 2.508600000000000e-01, + 2.590500000000000e-01, + 2.673500000000000e-01, + 2.757500000000000e-01, + 2.842700000000000e-01, + 2.928900000000000e-01, + 3.016200000000000e-01, + 3.104600000000000e-01, + 3.194000000000000e-01, + 3.284400000000000e-01, + 3.375800000000000e-01, + 3.468300000000000e-01, + 3.561700000000000e-01, + 3.656100000000000e-01, + 3.751400000000000e-01, + 3.847700000000000e-01, + 3.944900000000000e-01, + 4.043000000000000e-01, + 4.142000000000000e-01, + 4.241900000000000e-01, + 4.342700000000000e-01, + 4.444300000000000e-01, + 4.546700000000000e-01, + 4.650000000000000e-01, + 4.754100000000000e-01, + 4.859000000000000e-01, + 4.964600000000000e-01, + 5.071000000000000e-01, + 5.178199999999999e-01, + 5.286000000000000e-01, + 5.394600000000001e-01, + 5.503900000000000e-01, + 5.613800000000000e-01, + 5.724399999999999e-01, + 5.835700000000000e-01, + 5.947600000000000e-01, + 6.060100000000000e-01, + 6.173200000000000e-01, + 6.286800000000000e-01, + 6.401000000000000e-01, + 6.515800000000000e-01, + 6.631100000000000e-01, + 6.746900000000000e-01, + 6.863200000000000e-01, + 6.979900000000000e-01, + 7.097200000000000e-01, + 7.214800000000000e-01, + 7.332900000000000e-01, + 7.451300000000000e-01, + 7.570200000000000e-01, + 7.689400000000000e-01, + 7.809000000000000e-01, + 7.928900000000000e-01, + 8.049100000000000e-01, + 8.169600000000000e-01, + 8.290400000000000e-01, + 8.411400000000000e-01, + 8.532700000000000e-01, + 8.654200000000000e-01, + 8.775900000000000e-01, + 8.897800000000000e-01, + 9.019800000000000e-01, + 9.142000000000000e-01, + 9.264400000000000e-01, + 9.386800000000000e-01, + 9.509300000000001e-01, + 9.631900000000000e-01, + 9.754600000000000e-01, + 9.877300000000000e-01, + 1.000000000000000e+00 +}; + +static const scalar Re395_yplus[129] = { + 0.000000000000000e+00, + 2.953500000000000e-02, + 1.181400000000000e-01, + 2.657900000000000e-01, + 4.724800000000000e-01, + 7.381600000000000e-01, + 1.062800000000000e+00, + 1.446400000000000e+00, + 1.888800000000000e+00, + 2.390000000000000e+00, + 2.949900000000000e+00, + 3.568400000000000e+00, + 4.245500000000000e+00, + 4.981000000000000e+00, + 5.774800000000000e+00, + 6.626800000000000e+00, + 7.536900000000000e+00, + 8.504899999999999e+00, + 9.530700000000000e+00, + 1.061400000000000e+01, + 1.175500000000000e+01, + 1.295300000000000e+01, + 1.420900000000000e+01, + 1.552100000000000e+01, + 1.689000000000000e+01, + 1.831500000000000e+01, + 1.979700000000000e+01, + 2.133500000000000e+01, + 2.292900000000000e+01, + 2.457800000000000e+01, + 2.628300000000000e+01, + 2.804300000000000e+01, + 2.985800000000000e+01, + 3.172700000000000e+01, + 3.365100000000000e+01, + 3.562800000000000e+01, + 3.766000000000000e+01, + 3.974400000000000e+01, + 4.188200000000000e+01, + 4.407300000000000e+01, + 4.631600000000000e+01, + 4.861100000000000e+01, + 5.095800000000000e+01, + 5.335600000000000e+01, + 5.580500000000000e+01, + 5.830500000000000e+01, + 6.085500000000000e+01, + 6.345600000000000e+01, + 6.610500000000000e+01, + 6.880400000000000e+01, + 7.155100000000000e+01, + 7.434699999999999e+01, + 7.719100000000000e+01, + 8.008199999999999e+01, + 8.302000000000000e+01, + 8.600500000000000e+01, + 8.903500000000000e+01, + 9.211199999999999e+01, + 9.523399999999999e+01, + 9.840000000000001e+01, + 1.016100000000000e+02, + 1.048600000000000e+02, + 1.081600000000000e+02, + 1.115000000000000e+02, + 1.148900000000000e+02, + 1.183100000000000e+02, + 1.217800000000000e+02, + 1.252800000000000e+02, + 1.288300000000000e+02, + 1.324200000000000e+02, + 1.360400000000000e+02, + 1.397100000000000e+02, + 1.434100000000000e+02, + 1.471500000000000e+02, + 1.509200000000000e+02, + 1.547400000000000e+02, + 1.585800000000000e+02, + 1.624700000000000e+02, + 1.663900000000000e+02, + 1.703400000000000e+02, + 1.743300000000000e+02, + 1.783400000000000e+02, + 1.823900000000000e+02, + 1.864800000000000e+02, + 1.905900000000000e+02, + 1.947300000000000e+02, + 1.989100000000000e+02, + 2.031100000000000e+02, + 2.073400000000000e+02, + 2.116000000000000e+02, + 2.158900000000000e+02, + 2.202000000000000e+02, + 2.245400000000000e+02, + 2.289000000000000e+02, + 2.332900000000000e+02, + 2.377000000000000e+02, + 2.421400000000000e+02, + 2.466000000000000e+02, + 2.510800000000000e+02, + 2.555800000000000e+02, + 2.601000000000000e+02, + 2.646400000000000e+02, + 2.692000000000000e+02, + 2.737800000000000e+02, + 2.783800000000000e+02, + 2.830000000000000e+02, + 2.876300000000000e+02, + 2.922700000000000e+02, + 2.969400000000000e+02, + 3.016100000000000e+02, + 3.063000000000000e+02, + 3.110100000000000e+02, + 3.157200000000000e+02, + 3.204500000000000e+02, + 3.251900000000000e+02, + 3.299300000000000e+02, + 3.346900000000000e+02, + 3.394600000000000e+02, + 3.442300000000000e+02, + 3.490100000000000e+02, + 3.538000000000000e+02, + 3.585900000000000e+02, + 3.633900000000000e+02, + 3.681900000000000e+02, + 3.730000000000000e+02, + 3.778100000000000e+02, + 3.826200000000000e+02, + 3.874300000000000e+02, + 3.922400000000000e+02 +}; + +static const scalar Re395_U[129] = { + 0.000000000000000e+00, + 1.683537709989837e-03, + 6.731757022374556e-03, + 1.513918635410320e-02, + 2.690081008737908e-02, + 4.200636900398164e-02, + 6.044389401598693e-02, + 8.217063838080939e-02, + 1.071461351821347e-01, + 1.352449939750451e-01, + 1.662962277891985e-01, + 2.000148583406911e-01, + 2.360133151634138e-01, + 2.737615386493055e-01, + 3.126439757051341e-01, + 3.519709788838526e-01, + 3.910529007384139e-01, + 4.292285916501612e-01, + 4.659622924453555e-01, + 5.008208361324631e-01, + 5.335192444275801e-01, + 5.638979294917208e-01, + 5.918998956681040e-01, + 6.176049368762229e-01, + 6.410871474698924e-01, + 6.625175144194551e-01, + 6.820100290384720e-01, + 6.997926739540665e-01, + 7.160364361365806e-01, + 7.307983112427952e-01, + 7.443632775566141e-01, + 7.567883307348179e-01, + 7.682444577477489e-01, + 7.787886542521879e-01, + 7.886489028752580e-01, + 7.978252036169590e-01, + 8.063745521340716e-01, + 8.144679353969383e-01, + 8.221623490623399e-01, + 8.294007974734953e-01, + 8.362972719439662e-01, + 8.429657637873140e-01, + 8.493492773467582e-01, + 8.555618039358601e-01, + 8.616033435546199e-01, + 8.675308918598180e-01, + 8.732874531946738e-01, + 8.789300232159682e-01, + 8.844586019237012e-01, + 8.899301849746534e-01, + 8.952307810552632e-01, + 9.004743814790923e-01, + 9.056609862461408e-01, + 9.107335996996275e-01, + 9.157492174963338e-01, + 9.206508439794783e-01, + 9.255524704626230e-01, + 9.303971012889870e-01, + 9.352417321153507e-01, + 9.400293672849338e-01, + 9.447600067977363e-01, + 9.494906463105386e-01, + 9.541642901665602e-01, + 9.588379340225819e-01, + 9.633975865650419e-01, + 9.680142347642829e-01, + 9.725168916499621e-01, + 9.769625528788607e-01, + 9.813512184509786e-01, + 9.856828883663158e-01, + 9.899575626248721e-01, + 9.941752412266478e-01, + 9.983929198284234e-01, + 1.002553602773418e+00, + 1.006771281375194e+00, + 1.010988959976969e+00, + 1.015206638578745e+00, + 1.019424317180521e+00, + 1.023641995782296e+00, + 1.027859674384072e+00, + 1.032077352985848e+00, + 1.036295031587623e+00, + 1.040398718875837e+00, + 1.044502406164051e+00, + 1.048492102138704e+00, + 1.052481798113357e+00, + 1.056414498431228e+00, + 1.060290203092320e+00, + 1.064051916439849e+00, + 1.067756634130598e+00, + 1.071347360507785e+00, + 1.074881091228192e+00, + 1.078300830635037e+00, + 1.081663574385101e+00, + 1.084912326821604e+00, + 1.088161079258107e+00, + 1.091352836037829e+00, + 1.094430601503990e+00, + 1.097508366970150e+00, + 1.100472141122749e+00, + 1.103378919618568e+00, + 1.106171706800825e+00, + 1.108907498326300e+00, + 1.111586294194996e+00, + 1.114151098750130e+00, + 1.116715903305263e+00, + 1.119166716546836e+00, + 1.121560534131627e+00, + 1.123897356059638e+00, + 1.126120186674087e+00, + 1.128286021631756e+00, + 1.130337865275863e+00, + 1.132218721949628e+00, + 1.134042582966612e+00, + 1.135752452670034e+00, + 1.137348331059895e+00, + 1.138830218136195e+00, + 1.140141118242152e+00, + 1.141395022691329e+00, + 1.142534935826944e+00, + 1.143560857648997e+00, + 1.144472788157489e+00, + 1.145270727352420e+00, + 1.145954675233789e+00, + 1.146524631801596e+00, + 1.146923601399062e+00, + 1.147265575339746e+00, + 1.147436562310089e+00, + 1.147493557966869e+00 +}; + +static const scalar Re395_L[129] = { + 5.343185000000000e-01, + 5.343185000000000e-01, + 5.343185000000000e-01, + 5.343185000000000e-01, + 5.343185000000000e-01, + 5.343185000000000e-01, + 5.343185000000000e-01, + 5.343185000000000e-01, + 5.343185000000000e-01, + 5.343185000000000e-01, + 5.343185000000000e-01, + 5.343185000000000e-01, + 5.343185000000000e-01, + 5.343258566593407e-01, + 5.401655728439561e-01, + 5.464334466021978e-01, + 5.531287422681319e-01, + 5.602499885098902e-01, + 5.677938490769231e-01, + 5.717387276923077e-01, + 5.758937230769231e-01, + 5.802562861538462e-01, + 5.848300584615385e-01, + 5.896077569230769e-01, + 5.945930230769231e-01, + 5.997822153846153e-01, + 6.051789753846153e-01, + 6.094132404075546e-01, + 6.137989111630219e-01, + 6.183359068190855e-01, + 6.230269787375745e-01, + 6.278693755566600e-01, + 6.328630972763419e-01, + 6.387067543522267e-01, + 6.447231490890688e-01, + 6.509052761133604e-01, + 6.572593894736842e-01, + 6.637741823295454e-01, + 6.694304374715909e-01, + 6.752269084801137e-01, + 6.811609497727272e-01, + 6.872325613494318e-01, + 6.934417432102272e-01, + 6.997858497727273e-01, + 7.062648810369317e-01, + 7.128788370028408e-01, + 7.196250720880681e-01, + 7.263489329864724e-01, + 7.331965704214359e-01, + 7.401734573361082e-01, + 7.472744237513007e-01, + 7.545020546566077e-01, + 7.618537650624350e-01, + 7.693241732314410e-01, + 7.728104502183406e-01, + 7.763524981441048e-01, + 7.799479437772926e-01, + 7.835991603493450e-01, + 7.873037746288210e-01, + 7.910606000000000e-01, + 7.918954444782940e-01, + 7.927406920030464e-01, + 7.935989433358720e-01, + 7.944675977151562e-01, + 7.953492559025134e-01, + 7.962387163747144e-01, + 7.971411806549886e-01, + 7.980514472201066e-01, + 7.989747175932977e-01, + 7.999083910129474e-01, + 8.008498667174410e-01, + 8.018043462300076e-01, + 8.027666280274182e-01, + 8.037393128712871e-01, + 8.047198000000000e-01, + 8.098343334444676e-01, + 8.149756445509482e-01, + 8.201838998124610e-01, + 8.254323215669931e-01, + 8.307209098145447e-01, + 8.360630533861221e-01, + 8.414319746197124e-01, + 8.468544511773286e-01, + 8.523304830589706e-01, + 8.578332926026255e-01, + 8.633762686392998e-01, + 8.689728000000000e-01, + 8.691830656699252e-01, + 8.693948332374928e-01, + 8.696081027027027e-01, + 8.698228740655549e-01, + 8.700386466935021e-01, + 8.702559212190915e-01, + 8.704741970097758e-01, + 8.706939746981024e-01, + 8.709147536515239e-01, + 8.711370345025877e-01, + 8.713603166187465e-01, + 8.715846000000000e-01, + 8.629285521076588e-01, + 8.542340328913517e-01, + 8.455010423510787e-01, + 8.367295804868395e-01, + 8.279196472986345e-01, + 8.190712427864636e-01, + 8.101843669503266e-01, + 8.012782554522067e-01, + 7.923529082921038e-01, + 7.833698541460519e-01, + 7.743868000000000e-01, + 7.598716407039612e-01, + 7.452945830740374e-01, + 7.307175254441134e-01, + 7.160785694803047e-01, + 7.014086643495532e-01, + 6.867387592188018e-01, + 6.720069557541654e-01, + 6.572442031225865e-01, + 6.424814504910074e-01, + 6.276877486924860e-01, + 6.128630977270220e-01, + 5.980384467615582e-01, + 5.831828466291515e-01, + 5.683272464967450e-01, + 5.534406971973961e-01, + 5.385541478980470e-01, + 5.236675985986980e-01, + 5.087810492993490e-01, + 4.938945000000000e-01 +}; + +static const scalar Re395_Lkeps[129] = { + 1.664289358938353e-43, + 1.079207111946708e-08, + 6.917646913520677e-07, + 7.908359493903670e-06, + 4.457241969515345e-05, + 1.702303697729449e-04, + 5.073887856147327e-04, + 1.273178400079191e-03, + 2.815811120930599e-03, + 5.654802237749245e-03, + 1.050958553686484e-02, + 1.826059914399258e-02, + 2.974688742461301e-02, + 4.535327514167637e-02, + 6.453606255349646e-02, + 8.576111871085369e-02, + 1.070554545743771e-01, + 1.268251507443336e-01, + 1.443411137909917e-01, + 1.595013232996176e-01, + 1.725929155193667e-01, + 1.839473054080221e-01, + 1.938582581014135e-01, + 2.025427487164213e-01, + 2.101329090311211e-01, + 2.167672792710871e-01, + 2.226851829515993e-01, + 2.282037880705980e-01, + 2.336173611392537e-01, + 2.390778909178276e-01, + 2.446290872157351e-01, + 2.503922888236600e-01, + 2.563282600177446e-01, + 2.623700541806437e-01, + 2.687069242713673e-01, + 2.755870431522100e-01, + 2.831416145077278e-01, + 2.911589619570361e-01, + 2.995878349054479e-01, + 3.083265372548887e-01, + 3.172630295002534e-01, + 3.263486503189463e-01, + 3.354496195886431e-01, + 3.446920403253757e-01, + 3.537261660519386e-01, + 3.627232236378541e-01, + 3.723251787850627e-01, + 3.823383293109757e-01, + 3.924705802639217e-01, + 4.027010690107726e-01, + 4.128827545497974e-01, + 4.234264113943166e-01, + 4.340748218199580e-01, + 4.452924588662474e-01, + 4.568728895152261e-01, + 4.681922538193385e-01, + 4.793516421833623e-01, + 4.898598346604647e-01, + 5.001652593040559e-01, + 5.104300119900885e-01, + 5.206986983509322e-01, + 5.311361728500357e-01, + 5.408730041453185e-01, + 5.496856759048514e-01, + 5.576137391393682e-01, + 5.650229293959703e-01, + 5.732961244135784e-01, + 5.821891989591726e-01, + 5.911926787183273e-01, + 6.007120907840836e-01, + 6.100941035104538e-01, + 6.177821884753104e-01, + 6.244200363605632e-01, + 6.312414195975462e-01, + 6.383998599932088e-01, + 6.449359079832541e-01, + 6.512600681310129e-01, + 6.568871324829316e-01, + 6.620676770034978e-01, + 6.666357952600277e-01, + 6.713630340054566e-01, + 6.770059631597830e-01, + 6.820498425208822e-01, + 6.853109336087452e-01, + 6.886122910621065e-01, + 6.930784203322555e-01, + 6.988976207153511e-01, + 7.059190005229651e-01, + 7.122354391583325e-01, + 7.173533697592511e-01, + 7.223710058003835e-01, + 7.263834494776019e-01, + 7.298420005543480e-01, + 7.319232210321809e-01, + 7.322815511105965e-01, + 7.314357163140471e-01, + 7.299247998200248e-01, + 7.278620585377512e-01, + 7.251538547502624e-01, + 7.232352298788507e-01, + 7.211565886072751e-01, + 7.185127537378586e-01, + 7.162759804641324e-01, + 7.144541612653507e-01, + 7.116155587431247e-01, + 7.071175712181420e-01, + 7.018782054845693e-01, + 6.971310215642395e-01, + 6.926447443432519e-01, + 6.888649356014379e-01, + 6.856210638607619e-01, + 6.818460646161992e-01, + 6.770932597151812e-01, + 6.719347082875970e-01, + 6.680174666526330e-01, + 6.642689977904167e-01, + 6.585379901977008e-01, + 6.515872438367379e-01, + 6.436311683332884e-01, + 6.360458192130479e-01, + 6.297166760075837e-01, + 6.236128368572218e-01, + 6.171128788599741e-01, + 6.102622105242579e-01, + 6.029272432322437e-01, + 5.963797132441645e-01, + 5.909875544297788e-01, + 5.882037311171485e-01, + 5.874336317676262e-01 +}; + +static const scalar Re395_Ruu[129] = { + 5.337600000000000e-28, + 1.368400000000000e-04, + 2.182800000000000e-03, + 1.102000000000000e-02, + 3.471800000000000e-02, + 8.445400000000000e-02, + 1.743800000000000e-01, + 3.212300000000000e-01, + 5.433600000000000e-01, + 8.585100000000000e-01, + 1.280200000000000e+00, + 1.813200000000000e+00, + 2.449000000000000e+00, + 3.164400000000000e+00, + 3.922900000000000e+00, + 4.681100000000000e+00, + 5.395300000000000e+00, + 6.028600000000000e+00, + 6.555800000000000e+00, + 6.964200000000000e+00, + 7.252400000000000e+00, + 7.427400000000000e+00, + 7.501600000000000e+00, + 7.489700000000000e+00, + 7.406900000000000e+00, + 7.267900000000000e+00, + 7.086500000000000e+00, + 6.875000000000000e+00, + 6.644100000000000e+00, + 6.402500000000000e+00, + 6.157300000000000e+00, + 5.914100000000000e+00, + 5.676700000000000e+00, + 5.447700000000000e+00, + 5.229100000000000e+00, + 5.022300000000000e+00, + 4.827900000000000e+00, + 4.646200000000000e+00, + 4.476800000000000e+00, + 4.319400000000000e+00, + 4.173000000000000e+00, + 4.037000000000000e+00, + 3.910600000000000e+00, + 3.793000000000000e+00, + 3.683400000000000e+00, + 3.582100000000000e+00, + 3.489600000000000e+00, + 3.405300000000000e+00, + 3.327800000000000e+00, + 3.255700000000000e+00, + 3.188600000000000e+00, + 3.126900000000000e+00, + 3.070200000000000e+00, + 3.017800000000000e+00, + 2.967800000000000e+00, + 2.919100000000000e+00, + 2.871200000000000e+00, + 2.823900000000000e+00, + 2.777600000000000e+00, + 2.732400000000000e+00, + 2.688300000000000e+00, + 2.645300000000000e+00, + 2.603100000000000e+00, + 2.561200000000000e+00, + 2.519600000000000e+00, + 2.479300000000000e+00, + 2.441100000000000e+00, + 2.405500000000000e+00, + 2.371400000000000e+00, + 2.338000000000000e+00, + 2.304600000000000e+00, + 2.270800000000000e+00, + 2.236500000000000e+00, + 2.202300000000000e+00, + 2.168200000000000e+00, + 2.133500000000000e+00, + 2.098300000000000e+00, + 2.063100000000000e+00, + 2.027700000000000e+00, + 1.992000000000000e+00, + 1.956800000000000e+00, + 1.922900000000000e+00, + 1.889600000000000e+00, + 1.856200000000000e+00, + 1.823200000000000e+00, + 1.790800000000000e+00, + 1.759100000000000e+00, + 1.727400000000000e+00, + 1.694800000000000e+00, + 1.661000000000000e+00, + 1.626400000000000e+00, + 1.591000000000000e+00, + 1.555000000000000e+00, + 1.518400000000000e+00, + 1.480800000000000e+00, + 1.442500000000000e+00, + 1.403700000000000e+00, + 1.364700000000000e+00, + 1.326000000000000e+00, + 1.288100000000000e+00, + 1.251300000000000e+00, + 1.216200000000000e+00, + 1.182700000000000e+00, + 1.150000000000000e+00, + 1.117500000000000e+00, + 1.084300000000000e+00, + 1.050700000000000e+00, + 1.017300000000000e+00, + 9.847800000000000e-01, + 9.537700000000000e-01, + 9.245400000000000e-01, + 8.970600000000000e-01, + 8.711300000000000e-01, + 8.464500000000000e-01, + 8.229600000000000e-01, + 8.006200000000000e-01, + 7.792900000000000e-01, + 7.588900000000000e-01, + 7.393200000000000e-01, + 7.205100000000000e-01, + 7.026100000000000e-01, + 6.857600000000000e-01, + 6.702300000000000e-01, + 6.562000000000000e-01, + 6.438900000000000e-01, + 6.337500000000000e-01, + 6.262600000000000e-01, + 6.217400000000000e-01, + 6.202400000000000e-01 +}; + +static const scalar Re395_Rvv[129] = { + 4.901800000000000e-29, + 9.473400000000000e-11, + 2.360300000000000e-08, + 5.787800000000001e-07, + 5.440000000000000e-06, + 3.001400000000000e-05, + 1.175600000000000e-04, + 3.620400000000000e-04, + 9.320300000000000e-04, + 2.087900000000000e-03, + 4.185200000000000e-03, + 7.661000000000000e-03, + 1.300500000000000e-02, + 2.072400000000000e-02, + 3.130000000000000e-02, + 4.515700000000000e-02, + 6.263600000000000e-02, + 8.397100000000000e-02, + 1.092800000000000e-01, + 1.385400000000000e-01, + 1.716200000000000e-01, + 2.082400000000000e-01, + 2.480300000000000e-01, + 2.905000000000000e-01, + 3.351000000000000e-01, + 3.812100000000000e-01, + 4.281800000000000e-01, + 4.753500000000000e-01, + 5.221100000000000e-01, + 5.679000000000000e-01, + 6.122900000000000e-01, + 6.549300000000000e-01, + 6.955400000000000e-01, + 7.338800000000000e-01, + 7.697100000000000e-01, + 8.028500000000000e-01, + 8.331600000000000e-01, + 8.606100000000000e-01, + 8.852000000000000e-01, + 9.070000000000000e-01, + 9.260699999999999e-01, + 9.424700000000000e-01, + 9.562600000000000e-01, + 9.675400000000000e-01, + 9.765100000000000e-01, + 9.834100000000000e-01, + 9.885000000000000e-01, + 9.919700000000000e-01, + 9.939600000000000e-01, + 9.945500000000000e-01, + 9.937900000000000e-01, + 9.918300000000000e-01, + 9.888700000000000e-01, + 9.851200000000000e-01, + 9.807500000000000e-01, + 9.758300000000000e-01, + 9.704400000000000e-01, + 9.646500000000000e-01, + 9.584600000000000e-01, + 9.517900000000000e-01, + 9.445900000000000e-01, + 9.367900000000000e-01, + 9.284500000000000e-01, + 9.197600000000000e-01, + 9.108700000000000e-01, + 9.018500000000000e-01, + 8.926800000000000e-01, + 8.833500000000000e-01, + 8.738500000000000e-01, + 8.641799999999999e-01, + 8.543600000000000e-01, + 8.443900000000000e-01, + 8.342100000000000e-01, + 8.236900000000000e-01, + 8.128500000000000e-01, + 8.018100000000000e-01, + 7.907400000000000e-01, + 7.796800000000000e-01, + 7.685600000000000e-01, + 7.573100000000000e-01, + 7.458800000000000e-01, + 7.342500000000000e-01, + 7.224900000000000e-01, + 7.106400000000000e-01, + 6.986200000000000e-01, + 6.864200000000000e-01, + 6.741700000000000e-01, + 6.620100000000000e-01, + 6.500700000000000e-01, + 6.384300000000001e-01, + 6.271200000000000e-01, + 6.161300000000000e-01, + 6.054500000000000e-01, + 5.951200000000000e-01, + 5.852300000000000e-01, + 5.758900000000000e-01, + 5.671500000000000e-01, + 5.589600000000000e-01, + 5.512000000000000e-01, + 5.437200000000000e-01, + 5.364200000000000e-01, + 5.291900000000000e-01, + 5.219700000000000e-01, + 5.147300000000000e-01, + 5.075000000000000e-01, + 5.002900000000000e-01, + 4.931200000000000e-01, + 4.860500000000000e-01, + 4.792200000000000e-01, + 4.726800000000000e-01, + 4.664400000000000e-01, + 4.604700000000000e-01, + 4.547800000000000e-01, + 4.493200000000000e-01, + 4.441000000000000e-01, + 4.391500000000000e-01, + 4.345800000000000e-01, + 4.304400000000000e-01, + 4.267100000000000e-01, + 4.233900000000000e-01, + 4.204600000000000e-01, + 4.179600000000000e-01, + 4.159100000000000e-01, + 4.143100000000000e-01, + 4.131300000000000e-01, + 4.123400000000000e-01, + 4.118800000000000e-01, + 4.116600000000000e-01, + 4.116000000000000e-01 +}; + +static const scalar Re395_Rww[129] = { + 5.957500000000000e-28, + 5.294300000000000e-05, + 8.271500000000000e-04, + 4.039200000000000e-03, + 1.214700000000000e-02, + 2.785200000000000e-02, + 5.356400000000000e-02, + 9.094700000000000e-02, + 1.406200000000000e-01, + 2.020900000000000e-01, + 2.738400000000000e-01, + 3.536900000000000e-01, + 4.391900000000000e-01, + 5.280400000000000e-01, + 6.183400000000000e-01, + 7.087100000000000e-01, + 7.981900000000000e-01, + 8.861000000000000e-01, + 9.717800000000000e-01, + 1.054500000000000e+00, + 1.133300000000000e+00, + 1.207500000000000e+00, + 1.276200000000000e+00, + 1.338700000000000e+00, + 1.394800000000000e+00, + 1.444400000000000e+00, + 1.487800000000000e+00, + 1.525100000000000e+00, + 1.556800000000000e+00, + 1.583300000000000e+00, + 1.604700000000000e+00, + 1.621600000000000e+00, + 1.634400000000000e+00, + 1.643900000000000e+00, + 1.651200000000000e+00, + 1.656700000000000e+00, + 1.660600000000000e+00, + 1.662300000000000e+00, + 1.661800000000000e+00, + 1.659400000000000e+00, + 1.655300000000000e+00, + 1.650000000000000e+00, + 1.643900000000000e+00, + 1.637300000000000e+00, + 1.629900000000000e+00, + 1.620900000000000e+00, + 1.610500000000000e+00, + 1.599100000000000e+00, + 1.587200000000000e+00, + 1.575400000000000e+00, + 1.563700000000000e+00, + 1.551000000000000e+00, + 1.536700000000000e+00, + 1.521200000000000e+00, + 1.505200000000000e+00, + 1.488900000000000e+00, + 1.472100000000000e+00, + 1.455100000000000e+00, + 1.437900000000000e+00, + 1.420600000000000e+00, + 1.403500000000000e+00, + 1.386700000000000e+00, + 1.369600000000000e+00, + 1.351500000000000e+00, + 1.332200000000000e+00, + 1.311800000000000e+00, + 1.291000000000000e+00, + 1.270000000000000e+00, + 1.248800000000000e+00, + 1.227600000000000e+00, + 1.206600000000000e+00, + 1.185500000000000e+00, + 1.164100000000000e+00, + 1.142300000000000e+00, + 1.120200000000000e+00, + 1.098300000000000e+00, + 1.076700000000000e+00, + 1.055100000000000e+00, + 1.033900000000000e+00, + 1.013100000000000e+00, + 9.927900000000000e-01, + 9.725400000000000e-01, + 9.518700000000000e-01, + 9.305700000000000e-01, + 9.093599999999999e-01, + 8.892600000000001e-01, + 8.704200000000000e-01, + 8.524300000000000e-01, + 8.346900000000000e-01, + 8.169000000000000e-01, + 7.993900000000000e-01, + 7.824100000000000e-01, + 7.658700000000001e-01, + 7.493100000000000e-01, + 7.325800000000000e-01, + 7.157200000000000e-01, + 6.986800000000000e-01, + 6.815900000000000e-01, + 6.649200000000000e-01, + 6.491600000000000e-01, + 6.342400000000000e-01, + 6.200100000000000e-01, + 6.066000000000000e-01, + 5.937700000000000e-01, + 5.810100000000000e-01, + 5.681500000000000e-01, + 5.554100000000000e-01, + 5.430400000000000e-01, + 5.311300000000000e-01, + 5.197100000000000e-01, + 5.088200000000001e-01, + 4.985300000000000e-01, + 4.888300000000000e-01, + 4.798200000000000e-01, + 4.716400000000000e-01, + 4.639200000000000e-01, + 4.561800000000000e-01, + 4.484200000000000e-01, + 4.407100000000000e-01, + 4.334500000000000e-01, + 4.270800000000000e-01, + 4.216900000000000e-01, + 4.170900000000000e-01, + 4.131900000000000e-01, + 4.098400000000000e-01, + 4.068900000000000e-01, + 4.043100000000000e-01, + 4.024700000000000e-01, + 4.017900000000000e-01 +}; + +static const scalar Re395_Ruv[129] = { + 1.120600000000000e-30, + -2.441400000000000e-08, + -1.572800000000000e-06, + -1.811300000000000e-05, + -1.032400000000000e-04, + -4.002700000000000e-04, + -1.214300000000000e-03, + -3.100900000000000e-03, + -6.948600000000000e-03, + -1.401500000000000e-02, + -2.585800000000000e-02, + -4.413500000000000e-02, + -7.027600000000001e-02, + -1.051300000000000e-01, + -1.486900000000000e-01, + -1.999900000000000e-01, + -2.572400000000000e-01, + -3.180800000000000e-01, + -3.800100000000000e-01, + -4.407300000000000e-01, + -4.983700000000000e-01, + -5.516200000000000e-01, + -5.996899999999999e-01, + -6.422400000000000e-01, + -6.792400000000000e-01, + -7.109000000000000e-01, + -7.375900000000000e-01, + -7.598000000000000e-01, + -7.781000000000000e-01, + -7.930199999999999e-01, + -8.050300000000000e-01, + -8.145200000000000e-01, + -8.218400000000000e-01, + -8.273300000000000e-01, + -8.313300000000000e-01, + -8.341200000000000e-01, + -8.359000000000000e-01, + -8.368300000000000e-01, + -8.370000000000000e-01, + -8.364900000000000e-01, + -8.352300000000000e-01, + -8.332100000000000e-01, + -8.304400000000000e-01, + -8.268700000000000e-01, + -8.225200000000000e-01, + -8.176200000000000e-01, + -8.124900000000000e-01, + -8.072600000000000e-01, + -8.018000000000000e-01, + -7.959200000000000e-01, + -7.896100000000000e-01, + -7.829600000000000e-01, + -7.761300000000000e-01, + -7.692900000000000e-01, + -7.624300000000001e-01, + -7.555200000000000e-01, + -7.485100000000000e-01, + -7.414600000000000e-01, + -7.344500000000000e-01, + -7.273700000000000e-01, + -7.201200000000000e-01, + -7.127100000000000e-01, + -7.052100000000000e-01, + -6.976100000000000e-01, + -6.898500000000000e-01, + -6.819200000000000e-01, + -6.739600000000000e-01, + -6.659800000000000e-01, + -6.579400000000000e-01, + -6.497700000000000e-01, + -6.414800000000001e-01, + -6.329800000000000e-01, + -6.241500000000000e-01, + -6.150000000000000e-01, + -6.055100000000000e-01, + -5.956100000000000e-01, + -5.853000000000000e-01, + -5.746000000000000e-01, + -5.635000000000000e-01, + -5.520400000000000e-01, + -5.404000000000000e-01, + -5.288300000000000e-01, + -5.174700000000000e-01, + -5.062600000000000e-01, + -4.950700000000000e-01, + -4.838600000000000e-01, + -4.726900000000000e-01, + -4.616400000000000e-01, + -4.507200000000000e-01, + -4.399900000000000e-01, + -4.294900000000000e-01, + -4.191000000000000e-01, + -4.087600000000000e-01, + -3.985300000000000e-01, + -3.884300000000000e-01, + -3.783500000000000e-01, + -3.681000000000000e-01, + -3.576100000000000e-01, + -3.469500000000000e-01, + -3.362100000000000e-01, + -3.254400000000000e-01, + -3.146800000000000e-01, + -3.039300000000000e-01, + -2.931500000000000e-01, + -2.823100000000000e-01, + -2.713300000000000e-01, + -2.601900000000000e-01, + -2.489300000000000e-01, + -2.375300000000000e-01, + -2.260500000000000e-01, + -2.146200000000000e-01, + -2.032600000000000e-01, + -1.918800000000000e-01, + -1.804300000000000e-01, + -1.689800000000000e-01, + -1.575200000000000e-01, + -1.460100000000000e-01, + -1.344200000000000e-01, + -1.227200000000000e-01, + -1.109400000000000e-01, + -9.912000000000000e-02, + -8.730100000000000e-02, + -7.544099999999999e-02, + -6.345199999999999e-02, + -5.125200000000000e-02, + -3.877200000000000e-02, + -2.600100000000000e-02, + -1.303100000000000e-02, + 0.000000000000000e+00 +}; + +static const scalar Re395_Ruw[129] = { + -6.128200000000000e-30, + 5.946400000000000e-08, + 8.322800000000000e-07, + 3.361600000000000e-06, + 1.073200000000000e-05, + 2.829700000000000e-05, + 6.686500000000000e-05, + 1.450600000000000e-04, + 2.909700000000000e-04, + 5.472300000000000e-04, + 9.737300000000000e-04, + 1.652900000000000e-03, + 2.684700000000000e-03, + 4.172100000000000e-03, + 6.188900000000000e-03, + 8.756200000000000e-03, + 1.183400000000000e-02, + 1.538700000000000e-02, + 1.944900000000000e-02, + 2.413300000000000e-02, + 2.950100000000000e-02, + 3.541800000000000e-02, + 4.151100000000000e-02, + 4.733000000000000e-02, + 5.251200000000000e-02, + 5.692600000000000e-02, + 6.065900000000000e-02, + 6.383500000000000e-02, + 6.641800000000000e-02, + 6.820700000000000e-02, + 6.900100000000001e-02, + 6.889600000000000e-02, + 6.828500000000000e-02, + 6.748100000000000e-02, + 6.645000000000000e-02, + 6.493699999999999e-02, + 6.287500000000000e-02, + 6.071000000000000e-02, + 5.918200000000000e-02, + 5.852700000000000e-02, + 5.843700000000000e-02, + 5.851200000000000e-02, + 5.837600000000000e-02, + 5.808400000000000e-02, + 5.812500000000000e-02, + 5.853700000000000e-02, + 5.880900000000000e-02, + 5.849400000000000e-02, + 5.787600000000000e-02, + 5.703400000000000e-02, + 5.574600000000000e-02, + 5.425900000000000e-02, + 5.295700000000000e-02, + 5.181900000000000e-02, + 5.087200000000000e-02, + 5.051400000000000e-02, + 5.056400000000000e-02, + 5.022900000000000e-02, + 4.917100000000000e-02, + 4.756500000000000e-02, + 4.571600000000000e-02, + 4.391600000000000e-02, + 4.213400000000000e-02, + 4.007700000000000e-02, + 3.752000000000000e-02, + 3.435800000000000e-02, + 3.085800000000000e-02, + 2.709400000000000e-02, + 2.337100000000000e-02, + 2.050000000000000e-02, + 1.879500000000000e-02, + 1.789100000000000e-02, + 1.757200000000000e-02, + 1.751700000000000e-02, + 1.709500000000000e-02, + 1.591300000000000e-02, + 1.430600000000000e-02, + 1.281200000000000e-02, + 1.190700000000000e-02, + 1.164000000000000e-02, + 1.150300000000000e-02, + 1.109700000000000e-02, + 1.021800000000000e-02, + 9.000400000000000e-03, + 7.875300000000000e-03, + 7.161700000000000e-03, + 7.012100000000000e-03, + 7.380200000000000e-03, + 8.229900000000000e-03, + 9.395400000000000e-03, + 1.058500000000000e-02, + 1.183400000000000e-02, + 1.372800000000000e-02, + 1.656300000000000e-02, + 1.983500000000000e-02, + 2.255500000000000e-02, + 2.396500000000000e-02, + 2.409900000000000e-02, + 2.337800000000000e-02, + 2.223600000000000e-02, + 2.087100000000000e-02, + 1.922400000000000e-02, + 1.727300000000000e-02, + 1.511600000000000e-02, + 1.305500000000000e-02, + 1.136000000000000e-02, + 9.978700000000000e-03, + 8.767700000000000e-03, + 7.886799999999999e-03, + 7.615600000000000e-03, + 7.942400000000001e-03, + 8.637900000000000e-03, + 9.584000000000001e-03, + 1.054000000000000e-02, + 1.115900000000000e-02, + 1.120900000000000e-02, + 1.061700000000000e-02, + 9.526000000000000e-03, + 8.299800000000000e-03, + 7.228700000000000e-03, + 6.450200000000000e-03, + 6.169600000000000e-03, + 6.427700000000000e-03, + 7.023800000000000e-03, + 7.781800000000000e-03, + 8.748900000000000e-03, + 9.862700000000000e-03, + 1.073100000000000e-02, + 1.103200000000000e-02 +}; + +static const scalar Re395_Rvw[129] = { + 6.279300000000000e-32, + -2.852500000000000e-11, + -1.714100000000000e-09, + -1.905100000000000e-08, + -1.041100000000000e-07, + -3.827200000000000e-07, + -1.096900000000000e-06, + -2.628800000000000e-06, + -5.459600000000000e-06, + -1.009100000000000e-05, + -1.711900000000000e-05, + -2.776400000000000e-05, + -4.458700000000000e-05, + -7.104400000000001e-05, + -1.089300000000000e-04, + -1.546000000000000e-04, + -1.955200000000000e-04, + -2.102900000000000e-04, + -1.726200000000000e-04, + -6.049000000000000e-05, + 1.326800000000000e-04, + 3.908600000000000e-04, + 6.781800000000000e-04, + 9.551100000000000e-04, + 1.198900000000000e-03, + 1.401600000000000e-03, + 1.553200000000000e-03, + 1.636000000000000e-03, + 1.646500000000000e-03, + 1.602600000000000e-03, + 1.507300000000000e-03, + 1.326800000000000e-03, + 1.030300000000000e-03, + 6.411900000000000e-04, + 2.364200000000000e-04, + -1.361200000000000e-04, + -4.756900000000000e-04, + -7.841200000000000e-04, + -1.053900000000000e-03, + -1.265500000000000e-03, + -1.412600000000000e-03, + -1.493000000000000e-03, + -1.511800000000000e-03, + -1.519600000000000e-03, + -1.622800000000000e-03, + -1.829500000000000e-03, + -2.014500000000000e-03, + -2.082600000000000e-03, + -2.110100000000000e-03, + -2.162800000000000e-03, + -2.244100000000000e-03, + -2.322000000000000e-03, + -2.351900000000000e-03, + -2.301800000000000e-03, + -2.085100000000000e-03, + -1.620900000000000e-03, + -8.580200000000000e-04, + 1.178300000000000e-04, + 1.033700000000000e-03, + 1.644000000000000e-03, + 1.985100000000000e-03, + 2.229100000000000e-03, + 2.414300000000000e-03, + 2.512800000000000e-03, + 2.500600000000000e-03, + 2.372600000000000e-03, + 2.195200000000000e-03, + 2.051700000000000e-03, + 1.999900000000000e-03, + 1.970000000000000e-03, + 1.818900000000000e-03, + 1.433600000000000e-03, + 7.537500000000000e-04, + -7.928999999999999e-05, + -8.372000000000000e-04, + -1.341300000000000e-03, + -1.577100000000000e-03, + -1.517900000000000e-03, + -1.176200000000000e-03, + -7.049800000000000e-04, + -2.068900000000000e-04, + 2.052900000000000e-04, + 5.097100000000000e-04, + 7.621100000000000e-04, + 9.229200000000000e-04, + 9.301700000000000e-04, + 7.281400000000000e-04, + 2.735200000000000e-04, + -2.930300000000000e-04, + -7.456500000000000e-04, + -9.538400000000000e-04, + -9.538900000000000e-04, + -1.007500000000000e-03, + -1.297000000000000e-03, + -1.762000000000000e-03, + -2.230500000000000e-03, + -2.563700000000000e-03, + -2.718300000000000e-03, + -2.709500000000000e-03, + -2.572200000000000e-03, + -2.464600000000000e-03, + -2.547200000000000e-03, + -2.787100000000000e-03, + -3.009400000000000e-03, + -3.120200000000000e-03, + -3.177200000000000e-03, + -3.241000000000000e-03, + -3.317100000000000e-03, + -3.373800000000000e-03, + -3.413500000000000e-03, + -3.489700000000000e-03, + -3.672000000000000e-03, + -3.988200000000000e-03, + -4.389600000000000e-03, + -4.768500000000000e-03, + -5.041700000000000e-03, + -5.168700000000000e-03, + -5.104000000000000e-03, + -4.828700000000000e-03, + -4.405600000000000e-03, + -3.898400000000000e-03, + -3.346200000000000e-03, + -2.792700000000000e-03, + -2.237600000000000e-03, + -1.685500000000000e-03, + -1.185800000000000e-03, + -7.655400000000000e-04, + -3.883500000000000e-04, + 0.000000000000000e+00 +}; + +static const scalar Re590_r[129] = { + 0.000000000000000e+00, + 7.529800000000000e-05, + 3.011800000000000e-04, + 6.776200000000000e-04, + 1.204500000000000e-03, + 1.881900000000000e-03, + 2.709500000000000e-03, + 3.687400000000000e-03, + 4.815300000000000e-03, + 6.093000000000000e-03, + 7.520500000000000e-03, + 9.097400000000000e-03, + 1.082300000000000e-02, + 1.269900000000000e-02, + 1.472200000000000e-02, + 1.689500000000000e-02, + 1.921500000000000e-02, + 2.168300000000000e-02, + 2.429800000000000e-02, + 2.706000000000000e-02, + 2.996900000000000e-02, + 3.302400000000000e-02, + 3.622400000000000e-02, + 3.956900000000000e-02, + 4.306000000000000e-02, + 4.669400000000000e-02, + 5.047200000000000e-02, + 5.439300000000000e-02, + 5.845600000000000e-02, + 6.266099999999999e-02, + 6.700700000000000e-02, + 7.149400000000000e-02, + 7.611999999999999e-02, + 8.088600000000000e-02, + 8.579000000000001e-02, + 9.083200000000000e-02, + 9.601100000000000e-02, + 1.013300000000000e-01, + 1.067800000000000e-01, + 1.123600000000000e-01, + 1.180800000000000e-01, + 1.239300000000000e-01, + 1.299100000000000e-01, + 1.360300000000000e-01, + 1.422700000000000e-01, + 1.486400000000000e-01, + 1.551500000000000e-01, + 1.617800000000000e-01, + 1.685300000000000e-01, + 1.754100000000000e-01, + 1.824200000000000e-01, + 1.895400000000000e-01, + 1.967900000000000e-01, + 2.041600000000000e-01, + 2.116500000000000e-01, + 2.192600000000000e-01, + 2.269900000000000e-01, + 2.348300000000000e-01, + 2.427900000000000e-01, + 2.508600000000000e-01, + 2.590500000000000e-01, + 2.673500000000000e-01, + 2.757500000000000e-01, + 2.842700000000000e-01, + 2.928900000000000e-01, + 3.016200000000000e-01, + 3.104600000000000e-01, + 3.194000000000000e-01, + 3.284400000000000e-01, + 3.375800000000000e-01, + 3.468300000000000e-01, + 3.561700000000000e-01, + 3.656100000000000e-01, + 3.751400000000000e-01, + 3.847700000000000e-01, + 3.944900000000000e-01, + 4.043000000000000e-01, + 4.142000000000000e-01, + 4.241900000000000e-01, + 4.342700000000000e-01, + 4.444300000000000e-01, + 4.546700000000000e-01, + 4.650000000000000e-01, + 4.754100000000000e-01, + 4.859000000000000e-01, + 4.964600000000000e-01, + 5.071000000000000e-01, + 5.178199999999999e-01, + 5.286000000000000e-01, + 5.394600000000001e-01, + 5.503900000000000e-01, + 5.613800000000000e-01, + 5.724399999999999e-01, + 5.835700000000000e-01, + 5.947600000000000e-01, + 6.060100000000000e-01, + 6.173200000000000e-01, + 6.286800000000000e-01, + 6.401000000000000e-01, + 6.515800000000000e-01, + 6.631100000000000e-01, + 6.746900000000000e-01, + 6.863200000000000e-01, + 6.979900000000000e-01, + 7.097200000000000e-01, + 7.214800000000000e-01, + 7.332900000000000e-01, + 7.451300000000000e-01, + 7.570200000000000e-01, + 7.689400000000000e-01, + 7.809000000000000e-01, + 7.928900000000000e-01, + 8.049100000000000e-01, + 8.169600000000000e-01, + 8.290400000000000e-01, + 8.411400000000000e-01, + 8.532700000000000e-01, + 8.654200000000000e-01, + 8.775900000000000e-01, + 8.897800000000000e-01, + 9.019800000000000e-01, + 9.142000000000000e-01, + 9.264400000000000e-01, + 9.386800000000000e-01, + 9.509300000000001e-01, + 9.631900000000000e-01, + 9.754600000000000e-01, + 9.877300000000000e-01, + 1.000000000000000e+00 +}; + +static const scalar Re590_yplus[129] = { + 0.000000000000000e+00, + 4.421500000000000e-02, + 1.768500000000000e-01, + 3.978900000000000e-01, + 7.073000000000000e-01, + 1.105000000000000e+00, + 1.591000000000000e+00, + 2.165200000000000e+00, + 2.827500000000000e+00, + 3.577800000000000e+00, + 4.416000000000000e+00, + 5.341900000000000e+00, + 6.355500000000000e+00, + 7.456500000000000e+00, + 8.644900000000000e+00, + 9.920400000000001e+00, + 1.128300000000000e+01, + 1.273200000000000e+01, + 1.426800000000000e+01, + 1.589000000000000e+01, + 1.759700000000000e+01, + 1.939100000000000e+01, + 2.127000000000000e+01, + 2.323500000000000e+01, + 2.528400000000000e+01, + 2.741800000000000e+01, + 2.963700000000000e+01, + 3.193900000000000e+01, + 3.432500000000000e+01, + 3.679400000000000e+01, + 3.934600000000000e+01, + 4.198100000000000e+01, + 4.469800000000000e+01, + 4.749600000000000e+01, + 5.037600000000000e+01, + 5.333600000000000e+01, + 5.637700000000000e+01, + 5.949800000000000e+01, + 6.269800000000000e+01, + 6.597700000000000e+01, + 6.933499999999999e+01, + 7.277100000000000e+01, + 7.628400000000001e+01, + 7.987400000000000e+01, + 8.354100000000000e+01, + 8.728300000000000e+01, + 9.110100000000000e+01, + 9.499400000000000e+01, + 9.895999999999999e+01, + 1.030000000000000e+02, + 1.071100000000000e+02, + 1.113000000000000e+02, + 1.155600000000000e+02, + 1.198800000000000e+02, + 1.242800000000000e+02, + 1.287500000000000e+02, + 1.332900000000000e+02, + 1.378900000000000e+02, + 1.425700000000000e+02, + 1.473100000000000e+02, + 1.521100000000000e+02, + 1.569800000000000e+02, + 1.619200000000000e+02, + 1.669200000000000e+02, + 1.719900000000000e+02, + 1.771100000000000e+02, + 1.823000000000000e+02, + 1.875500000000000e+02, + 1.928600000000000e+02, + 1.982300000000000e+02, + 2.036500000000000e+02, + 2.091400000000000e+02, + 2.146800000000000e+02, + 2.202800000000000e+02, + 2.259300000000000e+02, + 2.316400000000000e+02, + 2.374000000000000e+02, + 2.432200000000000e+02, + 2.490800000000000e+02, + 2.550000000000000e+02, + 2.609700000000000e+02, + 2.669800000000000e+02, + 2.730500000000000e+02, + 2.791600000000000e+02, + 2.853200000000000e+02, + 2.915200000000000e+02, + 2.977700000000000e+02, + 3.040600000000000e+02, + 3.103900000000000e+02, + 3.167700000000000e+02, + 3.231900000000000e+02, + 3.296400000000000e+02, + 3.361400000000000e+02, + 3.426700000000000e+02, + 3.492400000000000e+02, + 3.558400000000000e+02, + 3.624800000000000e+02, + 3.691600000000000e+02, + 3.758700000000000e+02, + 3.826000000000000e+02, + 3.893700000000000e+02, + 3.961700000000000e+02, + 4.030000000000000e+02, + 4.098600000000000e+02, + 4.167400000000000e+02, + 4.236500000000000e+02, + 4.305800000000000e+02, + 4.375400000000000e+02, + 4.445200000000000e+02, + 4.515200000000000e+02, + 4.585400000000000e+02, + 4.655800000000000e+02, + 4.726400000000000e+02, + 4.797100000000000e+02, + 4.868100000000000e+02, + 4.939100000000000e+02, + 5.010400000000000e+02, + 5.081700000000000e+02, + 5.153200000000001e+02, + 5.224700000000000e+02, + 5.296400000000000e+02, + 5.368099999999999e+02, + 5.440000000000000e+02, + 5.511900000000001e+02, + 5.583800000000000e+02, + 5.655800000000000e+02, + 5.727800000000000e+02, + 5.799900000000000e+02, + 5.871900000000001e+02 +}; + +static const scalar Re590_U[129] = { + 0.000000000000000e+00, + 2.371075499776431e-03, + 9.487840037653017e-03, + 2.134402163620501e-02, + 3.792670109407034e-02, + 5.921389968694006e-02, + 8.515951570041008e-02, + 1.156456145162374e-01, + 1.504256055690046e-01, + 1.890813570303956e-01, + 2.309267038488144e-01, + 2.750396117361790e-01, + 3.202889804902494e-01, + 3.654579392773358e-01, + 4.094207485596639e-01, + 4.512446527202203e-01, + 4.902595687008060e-01, + 5.261063319822262e-01, + 5.586348439594444e-01, + 5.879576785862380e-01, + 6.142249344676437e-01, + 6.376510381822850e-01, + 6.586648428874096e-01, + 6.773735618723290e-01, + 6.942060482942909e-01, + 7.093767287319191e-01, + 7.231000297638371e-01, + 7.355903779686688e-01, + 7.470085932803818e-01, + 7.575154956329441e-01, + 7.672719049603233e-01, + 7.763314279071754e-01, + 7.848012777628124e-01, + 7.927886678165459e-01, + 8.004008113576878e-01, + 8.076377083862385e-01, + 8.146065721915093e-01, + 8.213074027735006e-01, + 8.277938067768681e-01, + 8.341193908462677e-01, + 8.402305483370437e-01, + 8.462344925385079e-01, + 8.521312234506602e-01, + 8.578671344288447e-01, + 8.636030454070291e-01, + 8.691781364512458e-01, + 8.746996208508067e-01, + 8.801138919610555e-01, + 8.854745564266485e-01, + 8.907280076029296e-01, + 8.959278521345548e-01, + 9.010204833768681e-01, + 9.060595079745256e-01, + 9.110449259275271e-01, + 9.159767372358726e-01, + 9.209085485442180e-01, + 9.257331465632517e-01, + 9.305041379376295e-01, + 9.352215226673515e-01, + 9.398853007524173e-01, + 9.445490788374832e-01, + 9.491056436332371e-01, + 9.536086017843354e-01, + 9.581115599354334e-01, + 9.625609114418755e-01, + 9.669566563036618e-01, + 9.713524011654481e-01, + 9.757481460272344e-01, + 9.800366775997087e-01, + 9.843252091721829e-01, + 9.885601341000014e-01, + 9.927950590278199e-01, + 9.969763773109824e-01, + 1.001157695594145e+00, + 1.005339013877308e+00, + 1.009466725515814e+00, + 1.013540830509665e+00, + 1.017561328858859e+00, + 1.021528220563398e+00, + 1.025387898978625e+00, + 1.029193970749196e+00, + 1.032946435875111e+00, + 1.036645294356371e+00, + 1.040290546192974e+00, + 1.043882191384921e+00, + 1.047420229932213e+00, + 1.050904661834848e+00, + 1.054335487092827e+00, + 1.057712705706151e+00, + 1.061089924319474e+00, + 1.064359929643486e+00, + 1.067629934967498e+00, + 1.070846333646854e+00, + 1.074009125681554e+00, + 1.077118311071598e+00, + 1.080120283172330e+00, + 1.083068648628406e+00, + 1.085963407439826e+00, + 1.088804559606590e+00, + 1.091645711773355e+00, + 1.094379650650807e+00, + 1.097006376238947e+00, + 1.099633101827088e+00, + 1.102206220770573e+00, + 1.104725733069401e+00, + 1.107138032078918e+00, + 1.109550331088435e+00, + 1.111855416808640e+00, + 1.114160502528845e+00, + 1.116411981604394e+00, + 1.118556247390631e+00, + 1.120646906532213e+00, + 1.122683959029138e+00, + 1.124613798236751e+00, + 1.126436424155053e+00, + 1.128205443428699e+00, + 1.129813642768376e+00, + 1.131368235463398e+00, + 1.132762008224453e+00, + 1.134102174340851e+00, + 1.135281520523281e+00, + 1.136353653416400e+00, + 1.137264966375551e+00, + 1.138069066045390e+00, + 1.138712345781261e+00, + 1.139194805583164e+00, + 1.139570052095756e+00, + 1.139784478674380e+00, + 1.139838085319036e+00 +}; + +static const scalar Re590_L[129] = { + 4.197654000000000e-01, + 4.197654000000000e-01, + 4.197654000000000e-01, + 4.197654000000000e-01, + 4.197654000000000e-01, + 4.197654000000000e-01, + 4.197654000000000e-01, + 4.197654000000000e-01, + 4.197654000000000e-01, + 4.197654000000000e-01, + 4.197654000000000e-01, + 4.197790882139738e-01, + 4.270813901528384e-01, + 4.350133499344978e-01, + 4.435749675589519e-01, + 4.527632896261879e-01, + 4.598816012354804e-01, + 4.674512721013727e-01, + 4.754754366631467e-01, + 4.839488708553326e-01, + 4.928663506124604e-01, + 5.022374808097561e-01, + 5.104690223414634e-01, + 5.190773135121951e-01, + 5.280535927024390e-01, + 5.374022407219512e-01, + 5.471232575707317e-01, + 5.559513010916580e-01, + 5.650997804325437e-01, + 5.745665010916581e-01, + 5.843514630690010e-01, + 5.931462315583126e-01, + 6.022126380744417e-01, + 6.115493349578164e-01, + 6.211596591265508e-01, + 6.310369367444169e-01, + 6.411845047295285e-01, + 6.515990261637716e-01, + 6.606481274226804e-01, + 6.699195844329897e-01, + 6.794144159793813e-01, + 6.891297945360825e-01, + 6.990628925773196e-01, + 7.092115724253536e-01, + 7.176203933682557e-01, + 7.262011972865374e-01, + 7.349562772865375e-01, + 7.438833402619173e-01, + 7.529778000000000e-01, + 7.580348690498590e-01, + 7.631795605832550e-01, + 7.684243920978363e-01, + 7.737568460959549e-01, + 7.791644050799623e-01, + 7.846721040451552e-01, + 7.902674254938852e-01, + 7.959503694261524e-01, + 8.017084183443085e-01, + 8.075666072436500e-01, + 8.134999011288805e-01, + 8.195083000000000e-01, + 8.214133126409366e-01, + 8.233457073937553e-01, + 8.253015725281873e-01, + 8.272848197745013e-01, + 8.292876256721595e-01, + 8.313178136816999e-01, + 8.333714720728534e-01, + 8.354486008456201e-01, + 8.375492000000000e-01, + 8.363682262340216e-01, + 8.351720000393313e-01, + 8.339648792527040e-01, + 8.327446849557523e-01, + 8.315135960668634e-01, + 8.302694336676499e-01, + 8.290143766764996e-01, + 8.277462461750246e-01, + 8.264694000000000e-01, + 8.230434842062025e-01, + 8.195886333127952e-01, + 8.161106343397001e-01, + 8.125979132470733e-01, + 8.090620440747587e-01, + 8.054972398028343e-01, + 8.019092874512220e-01, + 7.982924000000000e-01, + 7.932744660384690e-01, + 7.882246215271032e-01, + 7.831348888284437e-01, + 7.780132455799494e-01, + 7.728676694190790e-01, + 7.676822050709151e-01, + 7.624728078103750e-01, + 7.572314999999999e-01, + 7.501552098214286e-01, + 7.430360330357142e-01, + 7.358739696428571e-01, + 7.286797412946429e-01, + 7.214640696428571e-01, + 7.142055113839286e-01, + 7.069147881696429e-01, + 6.995919000000000e-01, + 6.947305324402309e-01, + 6.898549917971971e-01, + 6.849581915292663e-01, + 6.800472181780709e-01, + 6.751149852019785e-01, + 6.701685791426216e-01, + 6.652080000000000e-01, + 6.588169364297254e-01, + 6.524076647011309e-01, + 6.459801848142165e-01, + 6.395436008481422e-01, + 6.330797046445881e-01, + 6.266158084410339e-01, + 6.201246000000000e-01, + 6.149744971329077e-01, + 6.098099479744632e-01, + 6.046453988160185e-01, + 5.994664033662217e-01, + 5.942874079164249e-01, + 5.890939661752758e-01, + 5.839005244341265e-01, + 5.787070826929774e-01, + 5.735064178061521e-01, + 5.683057529193268e-01, + 5.630978648868253e-01, + 5.578972000000000e-01 +}; + +static const scalar Re590_Lkeps[129] = { + 1.793956702328780e-43, + 2.500670261071508e-08, + 1.606266380362454e-06, + 1.839039797513212e-05, + 1.035877293089363e-04, + 3.943728324386658e-04, + 1.169683018957784e-03, + 2.919131141896699e-03, + 6.414886315180999e-03, + 1.271444696028502e-02, + 2.288710240337253e-02, + 3.726547911631772e-02, + 5.468006876559973e-02, + 7.275358992816129e-02, + 8.934493506513441e-02, + 1.035078665136140e-01, + 1.153006304070393e-01, + 1.251828292948646e-01, + 1.336199344509259e-01, + 1.409816205473265e-01, + 1.475783629303827e-01, + 1.536809063755710e-01, + 1.595118703687390e-01, + 1.651990446299342e-01, + 1.708890436568559e-01, + 1.768861818054518e-01, + 1.833402068873848e-01, + 1.902467235181357e-01, + 1.975368191896557e-01, + 2.052130152228798e-01, + 2.133550300796443e-01, + 2.218696368266096e-01, + 2.307744015448879e-01, + 2.399211338480710e-01, + 2.492801123858805e-01, + 2.588846790953635e-01, + 2.687477515618852e-01, + 2.788081324311310e-01, + 2.891159141019418e-01, + 2.996796062356756e-01, + 3.102839132436510e-01, + 3.208024687307798e-01, + 3.313186532805980e-01, + 3.420203242389459e-01, + 3.530584767274762e-01, + 3.644854615538338e-01, + 3.761403020737028e-01, + 3.877578524412519e-01, + 3.994858705252315e-01, + 4.112941916808132e-01, + 4.230942793387439e-01, + 4.352830676354775e-01, + 4.476380911171339e-01, + 4.597791442777099e-01, + 4.718933012132248e-01, + 4.842806424783027e-01, + 4.968945629406032e-01, + 5.092695327286920e-01, + 5.211746674274160e-01, + 5.326826018446754e-01, + 5.441531126518806e-01, + 5.556227697991355e-01, + 5.672658820668975e-01, + 5.786007197598678e-01, + 5.892343647520798e-01, + 5.994306317748098e-01, + 6.098393056870663e-01, + 6.207536870923682e-01, + 6.312914480611881e-01, + 6.408101982175313e-01, + 6.500340699664456e-01, + 6.592621479870368e-01, + 6.678919267974072e-01, + 6.753040824391408e-01, + 6.820962979883348e-01, + 6.882701629027054e-01, + 6.944693628055562e-01, + 7.007027155097920e-01, + 7.061609036642372e-01, + 7.114985003422737e-01, + 7.172871976931973e-01, + 7.228413564506918e-01, + 7.279215438598963e-01, + 7.319483098924342e-01, + 7.348634096697936e-01, + 7.368635564239492e-01, + 7.380371107724003e-01, + 7.391948619128029e-01, + 7.405817638407975e-01, + 7.413347805532612e-01, + 7.409722308358556e-01, + 7.406699152827051e-01, + 7.408914127052421e-01, + 7.408540020232459e-01, + 7.402277269916649e-01, + 7.399175678056427e-01, + 7.388470242387847e-01, + 7.360736519705360e-01, + 7.319533743651223e-01, + 7.279350855552498e-01, + 7.244735165294008e-01, + 7.208881689858692e-01, + 7.173948230069089e-01, + 7.135145593133115e-01, + 7.091451533995295e-01, + 7.050708618788807e-01, + 7.014229992966290e-01, + 6.973589964146462e-01, + 6.930076016077142e-01, + 6.885049028253482e-01, + 6.832776022228736e-01, + 6.780928262716091e-01, + 6.726114618257262e-01, + 6.668552812718236e-01, + 6.608584119370413e-01, + 6.539869415053823e-01, + 6.462275059225702e-01, + 6.382272891230178e-01, + 6.310676902078552e-01, + 6.249359524016720e-01, + 6.180955623890663e-01, + 6.106270645802421e-01, + 6.034085346643908e-01, + 5.973389413372151e-01, + 5.931827622130115e-01, + 5.903035904643219e-01, + 5.879367636852111e-01, + 5.867196750936549e-01, + 5.866317370877272e-01 +}; + +static const scalar Re590_Ruu[129] = { + 7.800300000000000e-28, + 3.206100000000000e-04, + 5.113600000000000e-03, + 2.578500000000000e-02, + 8.109600000000000e-02, + 1.968300000000000e-01, + 4.048200000000000e-01, + 7.397899999999999e-01, + 1.232100000000000e+00, + 1.895700000000000e+00, + 2.715700000000000e+00, + 3.642800000000000e+00, + 4.601400000000000e+00, + 5.507900000000000e+00, + 6.291500000000000e+00, + 6.907300000000000e+00, + 7.339000000000000e+00, + 7.594000000000000e+00, + 7.694300000000000e+00, + 7.668400000000000e+00, + 7.545600000000000e+00, + 7.352900000000000e+00, + 7.113200000000000e+00, + 6.844800000000000e+00, + 6.561900000000000e+00, + 6.276300000000000e+00, + 5.996600000000000e+00, + 5.728200000000000e+00, + 5.474600000000000e+00, + 5.237600000000000e+00, + 5.018300000000000e+00, + 4.816600000000000e+00, + 4.631800000000000e+00, + 4.462500000000000e+00, + 4.307600000000000e+00, + 4.166200000000000e+00, + 4.037100000000000e+00, + 3.919400000000000e+00, + 3.812500000000000e+00, + 3.715800000000000e+00, + 3.628300000000000e+00, + 3.548600000000000e+00, + 3.475300000000000e+00, + 3.407100000000000e+00, + 3.343600000000000e+00, + 3.284500000000000e+00, + 3.229100000000000e+00, + 3.177300000000000e+00, + 3.128700000000000e+00, + 3.083200000000000e+00, + 3.040700000000000e+00, + 3.000800000000000e+00, + 2.962600000000000e+00, + 2.925500000000000e+00, + 2.889400000000000e+00, + 2.855000000000000e+00, + 2.821600000000000e+00, + 2.788500000000000e+00, + 2.755300000000000e+00, + 2.721800000000000e+00, + 2.688500000000000e+00, + 2.655400000000000e+00, + 2.622300000000000e+00, + 2.589300000000000e+00, + 2.556800000000000e+00, + 2.524200000000000e+00, + 2.491100000000000e+00, + 2.457700000000000e+00, + 2.424100000000000e+00, + 2.390200000000000e+00, + 2.356200000000000e+00, + 2.322200000000000e+00, + 2.287700000000000e+00, + 2.251900000000000e+00, + 2.214800000000000e+00, + 2.177100000000000e+00, + 2.139300000000000e+00, + 2.101700000000000e+00, + 2.063700000000000e+00, + 2.025200000000000e+00, + 1.986100000000000e+00, + 1.946300000000000e+00, + 1.906200000000000e+00, + 1.865800000000000e+00, + 1.825200000000000e+00, + 1.784800000000000e+00, + 1.744800000000000e+00, + 1.705300000000000e+00, + 1.666100000000000e+00, + 1.626900000000000e+00, + 1.587400000000000e+00, + 1.548000000000000e+00, + 1.509400000000000e+00, + 1.471700000000000e+00, + 1.434600000000000e+00, + 1.397900000000000e+00, + 1.361700000000000e+00, + 1.326000000000000e+00, + 1.290800000000000e+00, + 1.256100000000000e+00, + 1.222400000000000e+00, + 1.189500000000000e+00, + 1.157300000000000e+00, + 1.126200000000000e+00, + 1.096000000000000e+00, + 1.066500000000000e+00, + 1.037100000000000e+00, + 1.007800000000000e+00, + 9.786000000000000e-01, + 9.496300000000000e-01, + 9.209200000000000e-01, + 8.926800000000000e-01, + 8.652100000000000e-01, + 8.388900000000000e-01, + 8.135400000000000e-01, + 7.885600000000000e-01, + 7.640200000000000e-01, + 7.405900000000000e-01, + 7.188000000000000e-01, + 6.986900000000000e-01, + 6.799400000000000e-01, + 6.625200000000000e-01, + 6.467300000000000e-01, + 6.327800000000000e-01, + 6.209900000000000e-01, + 6.116500000000000e-01, + 6.049900000000000e-01, + 6.011200000000000e-01, + 5.999000000000000e-01 +}; + +static const scalar Re590_Rvv[129] = { + 5.417900000000001e-29, + 5.247400000000000e-10, + 1.289400000000000e-07, + 3.091300000000000e-06, + 2.816800000000000e-05, + 1.494500000000000e-04, + 5.590300000000000e-04, + 1.634200000000000e-03, + 3.975000000000000e-03, + 8.387100000000000e-03, + 1.580800000000000e-02, + 2.720100000000000e-02, + 4.344800000000000e-02, + 6.526200000000000e-02, + 9.312700000000000e-02, + 1.272500000000000e-01, + 1.675400000000000e-01, + 2.135700000000000e-01, + 2.646300000000000e-01, + 3.197400000000000e-01, + 3.777700000000000e-01, + 4.375200000000000e-01, + 4.977800000000000e-01, + 5.574700000000000e-01, + 6.155600000000000e-01, + 6.712200000000000e-01, + 7.237900000000000e-01, + 7.728100000000000e-01, + 8.179900000000000e-01, + 8.591200000000000e-01, + 8.961200000000000e-01, + 9.290000000000000e-01, + 9.579100000000000e-01, + 9.830300000000000e-01, + 1.004600000000000e+00, + 1.022800000000000e+00, + 1.038000000000000e+00, + 1.050300000000000e+00, + 1.059900000000000e+00, + 1.067100000000000e+00, + 1.072200000000000e+00, + 1.075500000000000e+00, + 1.077300000000000e+00, + 1.077800000000000e+00, + 1.076900000000000e+00, + 1.074700000000000e+00, + 1.071500000000000e+00, + 1.067300000000000e+00, + 1.062500000000000e+00, + 1.056900000000000e+00, + 1.050900000000000e+00, + 1.044500000000000e+00, + 1.037800000000000e+00, + 1.031100000000000e+00, + 1.024300000000000e+00, + 1.017200000000000e+00, + 1.009800000000000e+00, + 1.002000000000000e+00, + 9.938000000000000e-01, + 9.854100000000000e-01, + 9.768700000000000e-01, + 9.681800000000000e-01, + 9.593600000000000e-01, + 9.504500000000000e-01, + 9.415400000000000e-01, + 9.326000000000000e-01, + 9.236100000000000e-01, + 9.145500000000000e-01, + 9.054200000000000e-01, + 8.962599999999999e-01, + 8.870600000000000e-01, + 8.778200000000000e-01, + 8.685600000000000e-01, + 8.592900000000000e-01, + 8.498599999999999e-01, + 8.401000000000000e-01, + 8.298700000000000e-01, + 8.192199999999999e-01, + 8.083500000000000e-01, + 7.973000000000000e-01, + 7.860400000000000e-01, + 7.746100000000000e-01, + 7.630900000000000e-01, + 7.514600000000000e-01, + 7.396700000000000e-01, + 7.278000000000000e-01, + 7.158500000000000e-01, + 7.038200000000000e-01, + 6.917400000000000e-01, + 6.797299999999999e-01, + 6.678200000000000e-01, + 6.560100000000000e-01, + 6.443400000000000e-01, + 6.328700000000000e-01, + 6.216100000000000e-01, + 6.105800000000000e-01, + 5.998200000000000e-01, + 5.893900000000000e-01, + 5.792900000000000e-01, + 5.694500000000000e-01, + 5.597800000000001e-01, + 5.502600000000000e-01, + 5.408900000000000e-01, + 5.316500000000000e-01, + 5.224700000000000e-01, + 5.133400000000000e-01, + 5.042800000000000e-01, + 4.953200000000000e-01, + 4.864700000000000e-01, + 4.777500000000000e-01, + 4.692000000000000e-01, + 4.608800000000000e-01, + 4.528400000000000e-01, + 4.451200000000000e-01, + 4.377500000000000e-01, + 4.307700000000000e-01, + 4.242600000000000e-01, + 4.182400000000000e-01, + 4.127200000000000e-01, + 4.077800000000000e-01, + 4.035100000000000e-01, + 3.999500000000000e-01, + 3.970400000000000e-01, + 3.947200000000000e-01, + 3.929800000000000e-01, + 3.917900000000000e-01, + 3.910700000000000e-01, + 3.906900000000000e-01, + 3.905600000000000e-01 +}; + +static const scalar Re590_Rww[129] = { + 8.529200000000000e-28, + 1.312600000000000e-04, + 2.029700000000000e-03, + 9.738100000000000e-03, + 2.859100000000000e-02, + 6.362700000000000e-02, + 1.181600000000000e-01, + 1.929200000000000e-01, + 2.859100000000000e-01, + 3.931300000000000e-01, + 5.096700000000000e-01, + 6.309300000000000e-01, + 7.533300000000001e-01, + 8.744900000000000e-01, + 9.929000000000000e-01, + 1.107400000000000e+00, + 1.216600000000000e+00, + 1.319300000000000e+00, + 1.414000000000000e+00, + 1.499700000000000e+00, + 1.575800000000000e+00, + 1.642100000000000e+00, + 1.698600000000000e+00, + 1.746000000000000e+00, + 1.785000000000000e+00, + 1.816700000000000e+00, + 1.841700000000000e+00, + 1.860900000000000e+00, + 1.874600000000000e+00, + 1.883300000000000e+00, + 1.887400000000000e+00, + 1.887300000000000e+00, + 1.883700000000000e+00, + 1.877000000000000e+00, + 1.867700000000000e+00, + 1.856100000000000e+00, + 1.843100000000000e+00, + 1.829000000000000e+00, + 1.814000000000000e+00, + 1.797900000000000e+00, + 1.780300000000000e+00, + 1.761400000000000e+00, + 1.741800000000000e+00, + 1.722400000000000e+00, + 1.703300000000000e+00, + 1.684600000000000e+00, + 1.666000000000000e+00, + 1.647000000000000e+00, + 1.628000000000000e+00, + 1.609000000000000e+00, + 1.589600000000000e+00, + 1.570300000000000e+00, + 1.551900000000000e+00, + 1.534200000000000e+00, + 1.517000000000000e+00, + 1.500100000000000e+00, + 1.483300000000000e+00, + 1.466300000000000e+00, + 1.448600000000000e+00, + 1.430100000000000e+00, + 1.411700000000000e+00, + 1.393800000000000e+00, + 1.376600000000000e+00, + 1.359500000000000e+00, + 1.341700000000000e+00, + 1.323400000000000e+00, + 1.305200000000000e+00, + 1.287600000000000e+00, + 1.270400000000000e+00, + 1.253200000000000e+00, + 1.235900000000000e+00, + 1.218400000000000e+00, + 1.200500000000000e+00, + 1.181900000000000e+00, + 1.163100000000000e+00, + 1.144500000000000e+00, + 1.126400000000000e+00, + 1.108600000000000e+00, + 1.091200000000000e+00, + 1.074200000000000e+00, + 1.057400000000000e+00, + 1.040400000000000e+00, + 1.023100000000000e+00, + 1.005400000000000e+00, + 9.874300000000000e-01, + 9.694199999999999e-01, + 9.512300000000000e-01, + 9.327299999999999e-01, + 9.141800000000000e-01, + 8.957900000000000e-01, + 8.778800000000000e-01, + 8.601700000000000e-01, + 8.423000000000000e-01, + 8.245200000000000e-01, + 8.072000000000000e-01, + 7.903300000000000e-01, + 7.730000000000000e-01, + 7.548100000000000e-01, + 7.360900000000000e-01, + 7.172800000000000e-01, + 6.985300000000000e-01, + 6.798900000000000e-01, + 6.615000000000000e-01, + 6.431500000000000e-01, + 6.248200000000000e-01, + 6.070000000000000e-01, + 5.901800000000000e-01, + 5.741600000000000e-01, + 5.585000000000000e-01, + 5.434099999999999e-01, + 5.291500000000000e-01, + 5.157600000000000e-01, + 5.030100000000000e-01, + 4.907400000000000e-01, + 4.789300000000000e-01, + 4.676500000000000e-01, + 4.568100000000000e-01, + 4.465900000000000e-01, + 4.373600000000000e-01, + 4.290900000000000e-01, + 4.214400000000000e-01, + 4.144400000000000e-01, + 4.083600000000000e-01, + 4.035100000000000e-01, + 3.999800000000000e-01, + 3.975100000000000e-01, + 3.957700000000000e-01, + 3.948300000000000e-01, + 3.946000000000000e-01 +}; + +static const scalar Re590_Ruv[129] = { + -3.411000000000000e-32, + -8.364300000000000e-08, + -5.416500000000000e-06, + -6.283800000000000e-05, + -3.610000000000000e-04, + -1.407900000000000e-03, + -4.272200000000000e-03, + -1.080600000000000e-02, + -2.366400000000000e-02, + -4.592300000000000e-02, + -8.024600000000000e-02, + -1.279100000000000e-01, + -1.881200000000000e-01, + -2.580500000000000e-01, + -3.335000000000000e-01, + -4.099100000000000e-01, + -4.832500000000000e-01, + -5.506600000000000e-01, + -6.104600000000000e-01, + -6.620600000000000e-01, + -7.056100000000000e-01, + -7.417200000000000e-01, + -7.712100000000000e-01, + -7.949500000000000e-01, + -8.138100000000000e-01, + -8.286300000000000e-01, + -8.401300000000000e-01, + -8.489000000000000e-01, + -8.553700000000000e-01, + -8.599000000000000e-01, + -8.628000000000000e-01, + -8.643400000000000e-01, + -8.646700000000000e-01, + -8.638800000000000e-01, + -8.621500000000000e-01, + -8.596700000000000e-01, + -8.566300000000000e-01, + -8.532300000000000e-01, + -8.496300000000000e-01, + -8.459100000000001e-01, + -8.419500000000000e-01, + -8.376000000000000e-01, + -8.328300000000000e-01, + -8.276600000000000e-01, + -8.221900000000000e-01, + -8.164600000000000e-01, + -8.104300000000000e-01, + -8.040700000000000e-01, + -7.974500000000000e-01, + -7.906700000000000e-01, + -7.838500000000000e-01, + -7.770500000000000e-01, + -7.701500000000000e-01, + -7.629700000000000e-01, + -7.555400000000000e-01, + -7.480400000000000e-01, + -7.405700000000000e-01, + -7.331100000000000e-01, + -7.255700000000000e-01, + -7.179600000000000e-01, + -7.102600000000000e-01, + -7.024200000000000e-01, + -6.943400000000000e-01, + -6.860000000000001e-01, + -6.775200000000000e-01, + -6.689700000000000e-01, + -6.603400000000000e-01, + -6.516400000000000e-01, + -6.429200000000000e-01, + -6.342100000000001e-01, + -6.254300000000000e-01, + -6.165100000000000e-01, + -6.075000000000000e-01, + -5.984900000000000e-01, + -5.894100000000000e-01, + -5.802100000000000e-01, + -5.709500000000000e-01, + -5.617000000000000e-01, + -5.523900000000000e-01, + -5.429100000000000e-01, + -5.332500000000000e-01, + -5.235100000000000e-01, + -5.136800000000000e-01, + -5.037199999999999e-01, + -4.936900000000000e-01, + -4.836000000000000e-01, + -4.733200000000000e-01, + -4.628600000000000e-01, + -4.523200000000000e-01, + -4.417600000000000e-01, + -4.311200000000000e-01, + -4.203700000000000e-01, + -4.095400000000000e-01, + -3.986200000000000e-01, + -3.876600000000000e-01, + -3.768000000000000e-01, + -3.660600000000000e-01, + -3.552700000000000e-01, + -3.443600000000000e-01, + -3.333400000000000e-01, + -3.222400000000000e-01, + -3.109800000000000e-01, + -2.995300000000000e-01, + -2.879700000000000e-01, + -2.764200000000000e-01, + -2.649800000000000e-01, + -2.536800000000000e-01, + -2.423800000000000e-01, + -2.310200000000000e-01, + -2.196300000000000e-01, + -2.081900000000000e-01, + -1.966700000000000e-01, + -1.851200000000000e-01, + -1.735700000000000e-01, + -1.620100000000000e-01, + -1.504000000000000e-01, + -1.387300000000000e-01, + -1.270700000000000e-01, + -1.155100000000000e-01, + -1.040400000000000e-01, + -9.262200000000000e-02, + -8.116900000000001e-02, + -6.964700000000000e-02, + -5.806200000000000e-02, + -4.642200000000000e-02, + -3.476800000000000e-02, + -2.313800000000000e-02, + -1.154300000000000e-02, + 0.000000000000000e+00 +}; + +static const scalar Re590_Ruw[129] = { + -2.061400000000000e-30, + -3.664900000000000e-07, + -6.064600000000000e-06, + -3.085600000000000e-05, + -9.620100000000000e-05, + -2.303400000000000e-04, + -4.601200000000000e-04, + -8.066800000000000e-04, + -1.281000000000000e-03, + -1.896900000000000e-03, + -2.672800000000000e-03, + -3.613500000000000e-03, + -4.643000000000000e-03, + -5.583700000000000e-03, + -6.194100000000000e-03, + -6.249900000000000e-03, + -5.538400000000000e-03, + -3.915100000000000e-03, + -1.406100000000000e-03, + 1.726700000000000e-03, + 5.152000000000000e-03, + 8.662400000000001e-03, + 1.222000000000000e-02, + 1.572200000000000e-02, + 1.890300000000000e-02, + 2.141000000000000e-02, + 2.309200000000000e-02, + 2.393800000000000e-02, + 2.411400000000000e-02, + 2.385900000000000e-02, + 2.340900000000000e-02, + 2.294700000000000e-02, + 2.251100000000000e-02, + 2.198900000000000e-02, + 2.113900000000000e-02, + 1.968300000000000e-02, + 1.770200000000000e-02, + 1.544600000000000e-02, + 1.317900000000000e-02, + 1.123600000000000e-02, + 9.750699999999999e-03, + 8.284700000000001e-03, + 6.802000000000000e-03, + 5.899700000000000e-03, + 5.736000000000000e-03, + 5.947400000000000e-03, + 6.180600000000000e-03, + 6.473300000000000e-03, + 6.582100000000000e-03, + 6.246600000000000e-03, + 5.810000000000000e-03, + 5.399200000000000e-03, + 4.811600000000000e-03, + 4.235900000000000e-03, + 3.821300000000000e-03, + 3.503500000000000e-03, + 3.652500000000000e-03, + 4.460600000000000e-03, + 5.719000000000000e-03, + 6.914300000000000e-03, + 7.767100000000000e-03, + 8.382199999999999e-03, + 8.889200000000000e-03, + 9.284300000000001e-03, + 9.585500000000000e-03, + 9.706299999999999e-03, + 9.650300000000001e-03, + 9.582400000000000e-03, + 9.735600000000001e-03, + 9.658200000000000e-03, + 8.910000000000000e-03, + 7.502400000000000e-03, + 5.474300000000000e-03, + 3.320800000000000e-03, + 1.854900000000000e-03, + 1.386100000000000e-03, + 1.473000000000000e-03, + 1.838300000000000e-03, + 2.630700000000000e-03, + 4.104400000000000e-03, + 5.972600000000000e-03, + 7.790400000000000e-03, + 9.528800000000000e-03, + 1.114600000000000e-02, + 1.263400000000000e-02, + 1.401900000000000e-02, + 1.504400000000000e-02, + 1.544100000000000e-02, + 1.534100000000000e-02, + 1.494200000000000e-02, + 1.459200000000000e-02, + 1.424600000000000e-02, + 1.366100000000000e-02, + 1.303700000000000e-02, + 1.270700000000000e-02, + 1.284800000000000e-02, + 1.346600000000000e-02, + 1.419500000000000e-02, + 1.489100000000000e-02, + 1.544100000000000e-02, + 1.558100000000000e-02, + 1.556400000000000e-02, + 1.577500000000000e-02, + 1.631200000000000e-02, + 1.705700000000000e-02, + 1.789300000000000e-02, + 1.859500000000000e-02, + 1.902000000000000e-02, + 1.899400000000000e-02, + 1.836600000000000e-02, + 1.717300000000000e-02, + 1.565400000000000e-02, + 1.421600000000000e-02, + 1.298700000000000e-02, + 1.180500000000000e-02, + 1.064400000000000e-02, + 9.686699999999999e-03, + 8.933200000000001e-03, + 8.270100000000001e-03, + 7.662600000000000e-03, + 7.219900000000000e-03, + 6.947900000000000e-03, + 6.800000000000000e-03, + 6.515700000000000e-03, + 6.077000000000000e-03, + 5.638800000000000e-03, + 5.149900000000000e-03, + 4.663800000000000e-03, + 4.469300000000000e-03 +}; + +static const scalar Re590_Rvw[129] = { + -1.252100000000000e-31, + -1.275800000000000e-10, + -8.253900000000000e-09, + -9.286200000000000e-08, + -4.967400000000000e-07, + -1.726800000000000e-06, + -4.417100000000000e-06, + -8.739600000000000e-06, + -1.334900000000000e-05, + -1.416400000000000e-05, + -3.728300000000000e-06, + 2.765800000000000e-05, + 8.819599999999999e-05, + 1.798200000000000e-04, + 2.954500000000000e-04, + 4.222300000000000e-04, + 5.485800000000000e-04, + 6.687600000000000e-04, + 7.783500000000000e-04, + 8.742300000000000e-04, + 9.527200000000000e-04, + 1.005200000000000e-03, + 1.016800000000000e-03, + 9.789200000000001e-04, + 9.080600000000000e-04, + 8.405300000000000e-04, + 7.960400000000000e-04, + 7.760400000000000e-04, + 7.617300000000000e-04, + 7.539300000000000e-04, + 7.730200000000000e-04, + 8.343400000000000e-04, + 9.082000000000000e-04, + 9.583100000000000e-04, + 9.882199999999999e-04, + 1.020800000000000e-03, + 1.007400000000000e-03, + 9.114800000000000e-04, + 7.460400000000000e-04, + 5.596899999999999e-04, + 3.807300000000000e-04, + 2.123600000000000e-04, + 4.315900000000000e-05, + -2.032200000000000e-04, + -5.130999999999999e-04, + -7.063300000000000e-04, + -6.862200000000000e-04, + -5.752600000000000e-04, + -3.995100000000000e-04, + -1.218400000000000e-04, + 1.987400000000000e-04, + 5.752000000000000e-04, + 1.026000000000000e-03, + 1.418400000000000e-03, + 1.744900000000000e-03, + 2.127000000000000e-03, + 2.545300000000000e-03, + 2.852400000000000e-03, + 3.049800000000000e-03, + 3.265900000000000e-03, + 3.492000000000000e-03, + 3.603300000000000e-03, + 3.583200000000000e-03, + 3.487700000000000e-03, + 3.411700000000000e-03, + 3.425600000000000e-03, + 3.548700000000000e-03, + 3.775700000000000e-03, + 4.061200000000000e-03, + 4.352800000000000e-03, + 4.655100000000000e-03, + 4.987000000000000e-03, + 5.301400000000000e-03, + 5.495300000000000e-03, + 5.523600000000000e-03, + 5.458100000000000e-03, + 5.472400000000000e-03, + 5.680500000000000e-03, + 5.976000000000000e-03, + 6.199000000000000e-03, + 6.343900000000000e-03, + 6.471200000000000e-03, + 6.581600000000000e-03, + 6.719600000000000e-03, + 6.859200000000000e-03, + 6.907400000000000e-03, + 6.875000000000000e-03, + 6.813000000000000e-03, + 6.664200000000000e-03, + 6.421900000000000e-03, + 6.184800000000000e-03, + 6.057400000000000e-03, + 6.108600000000000e-03, + 6.265700000000000e-03, + 6.391100000000000e-03, + 6.361800000000000e-03, + 6.125000000000000e-03, + 5.772700000000000e-03, + 5.323400000000000e-03, + 4.797200000000000e-03, + 4.348600000000000e-03, + 4.041400000000000e-03, + 3.781700000000000e-03, + 3.461900000000000e-03, + 3.052300000000000e-03, + 2.670600000000000e-03, + 2.423800000000000e-03, + 2.275200000000000e-03, + 2.184300000000000e-03, + 2.143700000000000e-03, + 2.095500000000000e-03, + 1.993900000000000e-03, + 1.858000000000000e-03, + 1.685200000000000e-03, + 1.504600000000000e-03, + 1.355200000000000e-03, + 1.241100000000000e-03, + 1.181000000000000e-03, + 1.181100000000000e-03, + 1.196800000000000e-03, + 1.174600000000000e-03, + 1.122400000000000e-03, + 1.063400000000000e-03, + 9.533800000000000e-04, + 7.563000000000000e-04, + 5.133100000000000e-04, + 2.866300000000000e-04, + 1.174400000000000e-04, + 0.000000000000000e+00 +}; + +} // namespace dnsCatalogDetail + +static const DNSDataset dnsDatasets[] = { + { 180, 65, dnsCatalogDetail::Re180_r, dnsCatalogDetail::Re180_yplus, dnsCatalogDetail::Re180_U, dnsCatalogDetail::Re180_L, dnsCatalogDetail::Re180_Lkeps, dnsCatalogDetail::Re180_Ruu, dnsCatalogDetail::Re180_Rvv, dnsCatalogDetail::Re180_Rww, dnsCatalogDetail::Re180_Ruv, dnsCatalogDetail::Re180_Ruw, dnsCatalogDetail::Re180_Rvw }, + { 395, 129, dnsCatalogDetail::Re395_r, dnsCatalogDetail::Re395_yplus, dnsCatalogDetail::Re395_U, dnsCatalogDetail::Re395_L, dnsCatalogDetail::Re395_Lkeps, dnsCatalogDetail::Re395_Ruu, dnsCatalogDetail::Re395_Rvv, dnsCatalogDetail::Re395_Rww, dnsCatalogDetail::Re395_Ruv, dnsCatalogDetail::Re395_Ruw, dnsCatalogDetail::Re395_Rvw }, + { 590, 129, dnsCatalogDetail::Re590_r, dnsCatalogDetail::Re590_yplus, dnsCatalogDetail::Re590_U, dnsCatalogDetail::Re590_L, dnsCatalogDetail::Re590_Lkeps, dnsCatalogDetail::Re590_Ruu, dnsCatalogDetail::Re590_Rvv, dnsCatalogDetail::Re590_Rww, dnsCatalogDetail::Re590_Ruv, dnsCatalogDetail::Re590_Ruw, dnsCatalogDetail::Re590_Rvw }, +}; +static const int dnsDatasetsN = sizeof(dnsDatasets)/sizeof(dnsDatasets[0]); + +} // namespace Foam + +#endif diff --git a/src/finiteVolume/fields/fvPatchFields/derived/turbulentSyntheticEddyInlet/eddy/syntheticEddy.C b/src/finiteVolume/fields/fvPatchFields/derived/turbulentSyntheticEddyInlet/eddy/syntheticEddy.C new file mode 100644 index 0000000000..916c4363ca --- /dev/null +++ b/src/finiteVolume/fields/fvPatchFields/derived/turbulentSyntheticEddyInlet/eddy/syntheticEddy.C @@ -0,0 +1,86 @@ +/*---------------------------------------------------------------------------*\ + ========= | + \\ / F ield | OpenFOAM: The Open Source CFD Toolbox + \\ / O peration | Website: https://openfoam.org + \\ / A nd | Copyright (C) 2026 OpenFOAM Foundation + \\/ M anipulation | +------------------------------------------------------------------------------- +License + This file is part of OpenFOAM. + + OpenFOAM is free software: you can redistribute it and/or modify it + under the terms of the GNU General Public License as published by + the Free Software Foundation, either version 3 of the License, or + (at your option) any later version. + + OpenFOAM is distributed in the hope that it will be useful, but WITHOUT + ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or + FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License + for more details. + + You should have received a copy of the GNU General Public License + along with OpenFOAM. If not, see . + +\*---------------------------------------------------------------------------*/ + +#include "syntheticEddy.H" + +// * * * * * * * * * * * * * * * * Constructors * * * * * * * * * * * * * * // + +Foam::syntheticEddy::syntheticEddy() +: + position0_(Zero), + x_(0), + sigma_(Zero), + alpha_(Zero), + Rpg_(tensor::I), + c1_(-1), + Uconv_(0) +{} + + +Foam::syntheticEddy::syntheticEddy +( + const point& position0, + scalar x, + const vector& sigma, + const vector& alpha, + const tensor& Rpg, + scalar c1, + scalar Uconv +) +: + position0_(position0), + x_(x), + sigma_(sigma), + alpha_(alpha), + Rpg_(Rpg), + c1_(c1), + Uconv_(Uconv) +{} + + +// * * * * * * * * * * * * * * * Member Functions * * * * * * * * * * * * * // + +Foam::vector Foam::syntheticEddy::uPrime(const point& xp, const vector& n) const +{ + if (sigma_[2] <= 0) + { + return vector::zero; + } + + vector u(vector::zero); + + const vector rp(cmptDivide(Rpg_.T() & (xp - position(n)), sigma_)); + if (mag(rp) < scalar(1)) + { + const vector q(sigma_*(1 - magSqr(rp))); + const vector uPrimep(cmptMultiply(q, rp^alpha_)); + u += c1_*(Rpg_ & uPrimep); + } + + return u; +} + + +// ************************************************************************* // diff --git a/src/finiteVolume/fields/fvPatchFields/derived/turbulentSyntheticEddyInlet/eddy/syntheticEddy.H b/src/finiteVolume/fields/fvPatchFields/derived/turbulentSyntheticEddyInlet/eddy/syntheticEddy.H new file mode 100644 index 0000000000..cffde29a04 --- /dev/null +++ b/src/finiteVolume/fields/fvPatchFields/derived/turbulentSyntheticEddyInlet/eddy/syntheticEddy.H @@ -0,0 +1,159 @@ +/*---------------------------------------------------------------------------*\ + ========= | + \\ / F ield | OpenFOAM: The Open Source CFD Toolbox + \\ / O peration | Website: https://openfoam.org + \\ / A nd | Copyright (C) 2026 OpenFOAM Foundation + \\/ M anipulation | +------------------------------------------------------------------------------- +License + This file is part of OpenFOAM. + + OpenFOAM is free software: you can redistribute it and/or modify it + under the terms of the GNU General Public License as published by + the Free Software Foundation, either version 3 of the License, or + (at your option) any later version. + + OpenFOAM is distributed in the hope that it will be useful, but WITHOUT + ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or + FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License + for more details. + + You should have received a copy of the GNU General Public License + along with OpenFOAM. If not, see . + +Class + Foam::syntheticEddy + +Description + Single eddy used by turbulentSyntheticEddyInletFvPatchVectorField. + +SourceFiles + syntheticEddyI.H + syntheticEddy.C + +\*---------------------------------------------------------------------------*/ + +#ifndef syntheticEddy_H +#define syntheticEddy_H + +#include "vector.H" +#include "point.H" +#include "tensor.H" + +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +namespace Foam +{ + +class syntheticEddy; + +bool operator==(const syntheticEddy& a, const syntheticEddy& b); +bool operator!=(const syntheticEddy& a, const syntheticEddy& b); + +/*---------------------------------------------------------------------------*\ + Class syntheticEddy Declaration +\*---------------------------------------------------------------------------*/ + +class syntheticEddy +{ + // Private Data + + //- Reference position on the patch + point position0_; + + //- Streamwise distance from the reference position + scalar x_; + + //- Length scales in the principal frame + vector sigma_; + + //- Intensities in the principal frame + vector alpha_; + + //- Principal-to-global rotation tensor + tensor Rpg_; + + //- Kernel normalisation factor + scalar c1_; + + //- Per-eddy streamwise convection velocity. Set at spawn time. The + // BC's convectEddies() advances x_ by deltaT*Uconv_ each step. + scalar Uconv_; + + +public: + + // Constructors + + //- Construct null + syntheticEddy(); + + //- Construct from components + syntheticEddy + ( + const point& position0, + scalar x, + const vector& sigma, + const vector& alpha, + const tensor& Rpg, + scalar c1, + scalar Uconv + ); + + + // Member Functions + + //- Return the reference position + inline const point& position0() const noexcept; + + //- Return the streamwise distance + inline scalar x() const noexcept; + + //- Return the convection velocity + inline scalar Uconv() const noexcept; + + //- Return the current position + inline point position(const vector& n) const; + + //- Return the eddy volume + inline scalar volume() const; + + //- Advance the eddy by dx + inline void move(scalar dx); + + //- Return the velocity fluctuation at a point + vector uPrime(const point& xp, const vector& n) const; + + + // Friend Operators + + friend bool operator==(const syntheticEddy& a, const syntheticEddy& b) + { + return + a.position0_ == b.position0_ + && a.x_ == b.x_ + && a.sigma_ == b.sigma_ + && a.alpha_ == b.alpha_ + && a.Rpg_ == b.Rpg_ + && a.c1_ == b.c1_ + && a.Uconv_ == b.Uconv_; + } + + friend bool operator!=(const syntheticEddy& a, const syntheticEddy& b) + { + return !(a == b); + } +}; + + +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +} // End namespace Foam + +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +#include "syntheticEddyI.H" + +#endif + +// ************************************************************************* // diff --git a/src/finiteVolume/fields/fvPatchFields/derived/turbulentSyntheticEddyInlet/eddy/syntheticEddyI.H b/src/finiteVolume/fields/fvPatchFields/derived/turbulentSyntheticEddyInlet/eddy/syntheticEddyI.H new file mode 100644 index 0000000000..f4779740de --- /dev/null +++ b/src/finiteVolume/fields/fvPatchFields/derived/turbulentSyntheticEddyInlet/eddy/syntheticEddyI.H @@ -0,0 +1,66 @@ +/*---------------------------------------------------------------------------*\ + ========= | + \\ / F ield | OpenFOAM: The Open Source CFD Toolbox + \\ / O peration | Website: https://openfoam.org + \\ / A nd | Copyright (C) 2026 OpenFOAM Foundation + \\/ M anipulation | +------------------------------------------------------------------------------- +License + This file is part of OpenFOAM. + + OpenFOAM is free software: you can redistribute it and/or modify it + under the terms of the GNU General Public License as published by + the Free Software Foundation, either version 3 of the License, or + (at your option) any later version. + + OpenFOAM is distributed in the hope that it will be useful, but WITHOUT + ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or + FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License + for more details. + + You should have received a copy of the GNU General Public License + along with OpenFOAM. If not, see . + +\*---------------------------------------------------------------------------*/ + +#include "mathematicalConstants.H" + +// * * * * * * * * * * * * * * * Member Functions * * * * * * * * * * * * * // + +inline const Foam::point& Foam::syntheticEddy::position0() const noexcept +{ + return position0_; +} + + +inline Foam::scalar Foam::syntheticEddy::x() const noexcept +{ + return x_; +} + + +inline Foam::scalar Foam::syntheticEddy::Uconv() const noexcept +{ + return Uconv_; +} + + +inline Foam::point Foam::syntheticEddy::position(const vector& n) const +{ + return position0_ + n*x_; +} + + +inline Foam::scalar Foam::syntheticEddy::volume() const +{ + return Foam::constant::mathematical::pi*4.0/3.0*cmptProduct(sigma_); +} + + +inline void Foam::syntheticEddy::move(Foam::scalar dx) +{ + x_ += dx; +} + + +// ************************************************************************* // diff --git a/src/finiteVolume/fields/fvPatchFields/derived/turbulentSyntheticEddyInlet/turbulentSyntheticEddyInletFvPatchVectorField.C b/src/finiteVolume/fields/fvPatchFields/derived/turbulentSyntheticEddyInlet/turbulentSyntheticEddyInletFvPatchVectorField.C new file mode 100644 index 0000000000..b1d08eb7af --- /dev/null +++ b/src/finiteVolume/fields/fvPatchFields/derived/turbulentSyntheticEddyInlet/turbulentSyntheticEddyInletFvPatchVectorField.C @@ -0,0 +1,2478 @@ +/*---------------------------------------------------------------------------*\ + ========= | + \\ / F ield | OpenFOAM: The Open Source CFD Toolbox + \\ / O peration | Website: https://openfoam.org + \\ / A nd | Copyright (C) 2026 OpenFOAM Foundation + \\/ M anipulation | +------------------------------------------------------------------------------- +License + This file is part of OpenFOAM. + + OpenFOAM is free software: you can redistribute it and/or modify it + under the terms of the GNU General Public License as published by + the Free Software Foundation, either version 3 of the License, or + (at your option) any later version. + + OpenFOAM is distributed in the hope that it will be useful, but WITHOUT + ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or + FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License + for more details. + + You should have received a copy of the GNU General Public License + along with OpenFOAM. If not, see . + +\*---------------------------------------------------------------------------*/ + +#include "turbulentSyntheticEddyInletFvPatchVectorField.H" +#include "volFields.H" +#include "surfaceFields.H" +#include "fvcMeshPhi.H" +#include "addToRunTimeSelectionTable.H" +#include "mathematicalConstants.H" +#include "wallPolyPatch.H" +#include "cyclicPolyPatch.H" +#include "Switch.H" +#include "Pstream.H" +#include "scalarMatrices.H" +#include "meshSearch.H" +#include "clock.H" +#include "clockTime.H" +#include "ListOps.H" + +// * * * * * * * * * * * * * * * * Constructors * * * * * * * * * * * * * * // + +Foam::turbulentSyntheticEddyInletFvPatchVectorField:: +turbulentSyntheticEddyInletFvPatchVectorField +( + const fvPatch& p, + const DimensionedField& iF, + const dictionary& dict +) +: + fixedValueFvPatchVectorField(p, iF, dict, false), + Utau_(dict.lookupOrDefault("Utau", 0.0)), + UtauIsInDict_(dict.found("Utau", false, false)), + yPlusScale_(0), + avNu_(0), + U_(size(), Zero), + massFlowMode_(dict.found("massFlow")), + eddyConvectionFactor_ + ( + dict.lookupOrDefault("eddyConvectionFactor", 0.5) + ), + inFlow_ + ( + Function1::New + ( + massFlowMode_ ? "massFlow" : "Umean", + this->db().time().userUnits(), + iF.dimensions(), + dict + ) + ), + Umean_(0), + rhoInf_(dict.lookupOrDefault("rhoInf", 1)), + centerCoord_(Zero), + patchNormal_(Zero), + delta_(1), + wallPatchNames_(dict.lookupOrDefault("wallPatches", wordList())), + wallPatchesIsInDict_(dict.found("wallPatches")), + maxVertexWallDist_(0), + swirlNumber_(dict.lookupOrDefault("swirlNumber", 0)), + swirlAxis_(dict.lookupOrDefault("swirlAxis", Zero)), + swirlCenter_(dict.lookupOrDefault("swirlCenter", Zero)), + swirlAxisIsInDict_(dict.found("swirlAxis")), + swirlCenterIsInDict_(dict.found("swirlCenter")), + useUserR_(false), + useUserL_(false), + customRName_(dict.lookupOrDefault("customRName", word())), + customLName_(dict.lookupOrDefault("customLName", word())), + customUmeanName_(dict.lookupOrDefault("customUmeanName", word())), + inletRcorrection_(dict.lookupOrDefault("inletRcorrection", false)), + mcOverlaps_(dict.lookupOrDefault("inletRcorrectionOverlaps", 500)), + nEddy_(0), + eddyDensity_(dict.lookupOrDefault("eddyDensity", 5)), + Lscale_(dict.lookupOrDefault("Lscale", 1)), + Rscale_(dict.lookupOrDefault("Rscale", 1)), + nCellPerEddy_(dict.lookupOrDefault