
Utafiti huu unatengeneza framework mpya ya one-dimensional thermal-hydraulics kwa ajili ya kumodel circuits zote za nuclear reactors zinazojumuisha pipes, pumps, junctions, heat sources na hydraulic losses ndani ya mazingira ya OpenFOAM. Badala ya kubadilisha existing single-phase mass, momentum na energy equations za foamForNuclear platform, method inaongeza boundary conditions mpya, multi-branch junction models, pressure-jump calculations, four-region pump model na automatic mesh-generation tools zinazounganisha one-dimensional components huku zikihifadhi mass, momentum na energy. Katika analytical validation problems, flow, pressure na temperature results zimetoa relative errors chini ya %1; katika comparison na TRACE kwa simplified two-loop reactor circuit, steady-state core flow rates zimekubaliana ndani ya takriban %3. Hata hivyo, study inahusu single-phase flows pekee na imecompareiwa na analytical solutions pamoja na another system code badala ya experimental reactor system.
Main objective ya framework mpya ni kujenga common foundation inayoweza kuunganisha low-computational-cost one-dimensional system models na high-fidelity CFD au porous-media models katika same OpenFOAM infrastructure. Developed branchConnector1D boundary condition huhifadhi mass flow rate badala ya kusawazisha velocities wakati connected pipes zina cross-sections tofauti. pressureCoupling1D inaunganisha dynamic pressure, local hydraulic loss, pump pressure na hydrostatic differences ndani ya single pressure-jump equation. Four-region homologous pump model inaweza kuwakilisha flow kukaribia zero, reverse flow na inertia-driven pump rundown wakati wa pump trip.
Kwa mtazamo wa Uturuki: Study ina methodological value kwa nuclear-energy engineering, reactor-safety analysis, research reactors, small modular reactor studies, energy-systems education na open-source engineering software development nchini Uturuki. Ili approach itumike katika research na design projects za Uturuki, target reactor inapaswa kujengwa upya kwa real geometry, material, pump, valve, heat exchanger na control data; ivalidateiwe kwa experimental loops na international benchmark problems; na software quality assurance, version control na uncertainty analyses zikamilishwe. Research hii haiwezi kutafsiriwa kuwa developed framework iko tayari kwa licensing analyses, inaweza kumodel two-phase accidents, inaweza kuchukua nafasi ya TRACE au RELAP5 directly, au imethibitisha safety ya real nuclear facility.
Kwa nini one-dimensional reactor system codes zinahitajika?
Thermal-hydraulic behavior ya nuclear power plants inaamuliwa na kufanya kazi pamoja kwa reactor core, hot na cold legs, pumps, steam generators, pressure vessels, connecting pipes na safety systems. Kumodel entire plant kwa high-resolution three-dimensional CFD kunaweza kuunda very high computational cost kwa routine safety analyses.
One-dimensional system codes huwakilisha pipes na equipment kama control volumes zilizopangwa along axis. Hivyo pressure, flow rate, temperature na energy distribution zinaweza kufuatiliwa katika entire reactor circuit kwa relatively low computational cost. TRACE, RELAP5 na ATHLET ni established examples za approach hii.
Hata hivyo, katika new reactor designs, three-dimensional phenomena kama natural circulation, large pools, thermal stratification na complex mixing zones pia zimekuwa muhimu. Kwa hiyo badala ya single fidelity level, multi-fidelity methods zinatengenezwa ambapo baadhi ya regions zinasolveiwa kwa one-dimensional system model, baadhi kwa porous-media approach, na critical regions kwa CFD.
Ni deficiency gani iliyokuwepo ndani ya foamForNuclear?
foamForNuclear ni OpenFOAM-based multiphysics na multifidelity nuclear-analysis platform. Kupitia development na integration ya GeN-Foam na OFFBEAT solvers, platform inalenga kuendesha different physics kama flow, heat transfer, neutronics na fuel behavior katika common environment.
Kulingana na study, ingawa foamForNuclear iliweza kuendesha CFD na porous-media-based thermal-hydraulic regions ndani ya same simulation, haikuwa na local one-dimensional system infrastructure inayohitajika kujenga entire plant circuit. Components kama pump, valve, elbow, T-junction, diameter change na local loss zilihitaji external system code au OpenModelica coupling.
Watafiti walitengeneza local framework inayofanya kazi ndani ya OpenFOAM ili kuondoa dependency hii. Kauli kwamba work ni “first implementation directly embedded in OpenFOAM” ni priority claim ya researchers wenyewe; study haitoi comprehensive independent software-priority review.
Developed framework inabadilisha nini?
Study haiandiki upya existing single-phase porous-media thermal-hydraulic equations za foamForNuclear. Field equations zinazotatua conservation of mass, momentum na enthalpy zinabaki vilevile. New work inaongeza hasa infrastructure hii:
- Boundary conditions zinazounganisha one-dimensional pipes katika different orientations na cross-sections,
- Velocity coupling inayohifadhi mass flow kwenye area changes,
- Junctions ambapo three au more branches zinakutana,
- Pressure jump inayojumuisha dynamic pressure, local loss, pump na hydrostatic effects,
- Pump characteristics zinazofunika forward na reverse flow,
- Pump-trip na inertia-driven rundown model,
- Mass-flow-weighted enthalpy mixing katika junctions,
- Automatic one-dimensional pipe-network generation tools kupitia Python.
Kwa nini existing one-dimensional approach ya OpenFOAM haikutosha?
OpenFOAM inaweza kutumia empty boundary condition kwenye baadhi ya faces za thin three-dimensional mesh ili kuunda two- au one-dimensional solutions. Lakini kwa kuwa empty condition inaruhusu solution kubadilika only in directions tangent to corresponding face, inasababisha problem katika pipe networks zinazobadilisha direction. Kwa mfano, circuit yenye 90° elbow haiwezi kuwakilishwa kama straight mesh katika one direction.
Kama solution, slip wall inaweza kutumika kwenye outer faces za pipes na cyclic, au cyclicAMI kwa nonconformal meshes, katika junctions. Ingawa method hii inaweza kuunganisha pipes kwa different angles, study imebainisha three main problems:
- Mass conservation: Kutumia same velocity kwenye boundary ya two pipes zenye different cross-sectional areas hakuhakikishi same mass flow inapita.
- Applying same boundary type to all fields: Pressure, velocity na enthalpy zinapaswa coupled; lakini auxiliary field kama hydraulic diameter si lazima iwe continuous across junction.
- Complex mesh generation: Kuweka manually translation vector, rotation angle na rotation center kwa kila connection huongeza risk of error katika large circuits.
branchConnector1D inafanya nini?
Watafiti walitengeneza connection boundary condition inayoitwa branchConnector1D ili different field variables zipate different behavior. Fields zinazohitaji physical coupling kama pressure, velocity, enthalpy na selected turbulence variables zina-matched kwa cyclic coupling, huku fields zinazoweza kubaki different across junction, kama hydraulic diameter, zikitumia zero gradient.
Distinction hii inahakikisha only variables ambazo lazima zihifadhiwe au ziwe continuous zina strong coupling na opposite pipe. Hivyo artificial connection inayolazimisha fields zote kuwa same value inaepukwa.
Mass inahifadhiwaje katika pipes zenye different diameters?
Wakati cross-sections, densities au effective flow areas za two pipes ni tofauti, mass flow rates zinapaswa kusawazishwa badala ya velocities. Fundamental relation iliyotumika katika study ni:
\[ \alpha_P \rho_P A_P U_P = \alpha_N \rho_N A_N U_N \]
Hapa \(\alpha\) inaonyesha fraction inayochukuliwa na phase katika control volume, \(\rho\) density, \(A\) cross-sectional area na \(U\) velocity. \(P\) na \(N\) zinaonyesha control volumes kwenye pande mbili za junction.
Velocity contribution kutoka neighboring side inascaleiwa kwa area, density na phase fraction. Kwa hivyo velocity inaweza kuongezeka kwenye smaller pipe na kupungua kwenye larger pipe huku total mass flow kupitia junction ikibaki same.
Multi-branch junctions zimemodeliwaje?
Standard cyclic OpenFOAM connection imedesigniwa binary relation kati ya two faces. Katika T-junction, main pipe moja inaweza kugawanyika kuwa two branches au two branches kuungana kuwa one pipe. Figure 1 kwenye page 7 ya study inaonyesha arrangement ambapo two separate connection faces zimeundwa mwisho wa main pipe na kuunganishwa na upper na lower branches.
Mass conservation kwenye three-branch junction imeandikwa:
\[ \alpha_a\rho_aA_aU_a+ \alpha_b\rho_bA_bU_b+ \alpha_c\rho_cA_cU_c=0 \]
Velocities zinachukuliwa positive zinapoelekea junction na negative zinapoondoka junction. Effective neighbor velocity ya one boundary face inahesabiwa si kutoka directly matched face pekee, bali pia kutoka mass flows zinazotoka all other branches katika junction.
Kwa kuwa contribution kutoka other branches inaongezwa kama explicit source term, formulation ni semi-implicit badala ya fully implicit. Optional relaxation inaweza kutumika kwa explicit term hii ili kuboresha convergence katika highly unbalanced flow distributions.
Kwa kuwa kuunda multiple connection faces mwisho wa one cell kunaweza kuharibu geometric cell volume, watafiti walirecalculate volumes za cells zilizo adjacent to junction pamoja na face contributions katika pressure–velocity coupling.
Pressure jump inajumuisha effects zipi?
Kuhifadhi mass flow kupitia junction haimaanishi pressures kwenye pande mbili lazima ziwe equal. Area change, elbow, valve, pump na elevation difference zinaweza kuunda physical pressure jump.
Kwa purpose hii boundary condition inayoitwa pressureCoupling1D ilitengenezwa. Katika OpenFOAM solvers zenye buoyancy, transported modified pressure imefafanuliwa kama:
\[ p_{rgh}=p-\rho gz \]
Pressure jump kati ya pande mbili imeelezwa katika study kwa components hizi:
\[ p_{rgh,1}-p_{rgh,2} = \frac{1}{2}\rho\left(U_1^2-U_2^2\right) -\frac{1}{2}K\rho_uU_u^2 +S_p +gz\left(\rho_2-\rho_1\right) \]
First term inawakilisha dynamic-pressure difference, second local hydraulic loss, \(S_p\) pump pressure increase, na last term hydrostatic correction kutokana na density na elevation.
Kwa kuwa local-loss coefficient \(K\) na pump source zinaweza kutolewa kama time-dependent functions, transient events kama valve movement au pump trip zinaweza kumodeliwa. Flow ikireverse, loss term inaelekezwa upya relative to new inlet side.
Ni pump models gani mbili zimetengenezwa?
First model inafafanua pump head kama polynomial ya volumetric flow:
\[ h(Q)=\sum_i c_iQ^i \]
Pump speed ikibadilika, coefficients zinascaleiwa kulingana na affinity laws. Model hii ni simple na useful karibu na design operating point katika forward-flow conditions. Lakini kwa kuwa \(h(Q)\) ni single-valued, haiwezi kuwakilisha kikamilifu zero flow, reverse flow au reverse rotation.
Second model ni homologous au four-region pump approach inayotumika katika system codes. Flow rate, rotational speed na head zinanondimensionalizeiwa relative to nominal values:
\[ q^*=\frac{Q}{Q_R}, \qquad \omega^*=\frac{\omega}{\omega_R}, \qquad h^*=\frac{h}{H_R} \]
Forward na reverse flow pamoja na forward na reverse rotation zinagawanywa katika four operating regions. Katika kila region, pump head na hydraulic torque zinahesabiwa kutoka piecewise-linear tables. Arrangement hii inawakilisha conventional eight-octant pump characteristics katika four regions kwa kuunganisha octants zenye similar normalization.
Pump trip na inertia-driven rundown zinahesabiwaje?
Pump speed inaweza kupewa na user kama time schedule au kuhesabiwa kutoka rotor dynamics baada ya pump trip. Motor torque ikiondolewa, rotational speed hubadilika kwa balance hii:
\[ I\frac{d\omega}{dt} = -T_{\mathrm{hy}}\left(q^*,\omega^*\right) -T_f(\omega) \]
Hapa \(I\) ni rotor moment of inertia, \(T_{\mathrm{hy}}\) hydraulic torque inayotolewa na fluid, na \(T_f\) mechanical friction torque.
Rotor equation ikisolveiwa pamoja na four-region pump characteristics, pump inaweza kupita continuously kutoka normal operating point kwenda low flow, zero flow na reverse-flow region. Feature hii ni muhimu kwa reactor transients zenye loss of pump na flow reversal.
Energy conservation inahakikishwaje kwenye junctions?
Katika point ambapo multiple branches zinaungana, outlet enthalpy haiwezi kuwekwa sawa na enthalpy ya one neighboring pipe. Outlet fluid ni mixture ya all streams zinazoingia junction. Study inahesabu inlet enthalpy kwa mass-flow-weighted average:
\[ h_{\mathrm{mix}} = \frac{\sum_i\dot{m}_ih_i} {\sum_i\dot{m}_i} \]
\[ \dot{m}_i = \max\left(\rho_i\alpha_iA_iU_i,0\right) \]
Summation inafanywa only kwa branches zinazoflow kuelekea junction. Branch ikiflow kutoka junction kwenda nje, enthalpy inasafirishwa kwa zero-gradient behavior. Flow direction ikibadilika wakati wa transient, branches zinazoshiriki katika mixing zinabainishwa dynamically tena.
Python interface inarahisishaje circuit setup?
foamForNuclear Python API imepanuliwa ili one-dimensional pipe segments ziweze kufafanuliwa kwa start position, direction, length, hydraulic diameter, number of cells na elbow radius.
Start ya new pipe inaweza kutolewa kama existing pipe object badala ya coordinate vector. Katika hali hii new pipe inaunganishwa automatically kwenye outlet ya previous pipe; required rotation, translation na boundary conditions zinatengenezwa na software.
Kwa independent pipes inlet na outlet faces, na kwa connected pipes branchConnector1D connections zinagenerateiwa automatically. Hivyo manual calculation ya rotation center na geometric transformation kwa kila junction katika large circuit haihitajiki.
First analytical validation: flow splitting na recombination
Katika first validation problem, water flow yenye inlet flow rate 570 L·s−1 iligawanyika katika two parallel branches zenye different diameters, lengths na friction factors, kisha ikaungana tena. Kuna 9 m elevation difference kati ya two T-junctions.
Results zilizopatikana kwenye 1090-cell mesh ni:
| Variable | Analytical solution | FFN result | Relative error |
|---|---|---|---|
| First-branch flow rate | 121,8 L·s−1 | 122,2 L·s−1 | %0,33 |
| Second-branch flow rate | 448,2 L·s−1 | 447,8 L·s−1 | %0,089 |
| A–B pressure difference | 147 kPa | 147,15 kPa | %0,10 |
Results hizi zinaonyesha mass flow inayogawanyika katika different branches na distributed friction losses zimehesabiwa kwa agreement na analytical solution.
Second validation: pumped na natural-circulation closed loop
Second problem ni closed loop iliyoundwa na four pipes zenye diameters 80, 90, 100 na 110 mm. Kila pipe ina length 40 m. Katika first configuration flow inaendeshwa na pump; katika second configuration na buoyancy inayotokana na heating lower pipe hadi 350 K na cooling upper pipe hadi 300 K.
| Condition | Analytical flow rate | FFN flow rate | Reported error |
|---|---|---|---|
| Pumped loop | 0,0976 m3·s−1 | 0,0976 m3·s−1 | <%0,05 |
| Natural-circulation loop | 0,7345 m3·s−1 | 0,7341 m3·s−1 | <%0,05 |
Source inaonyesha analytical na numerical results ziko very close. Hata hivyo, magnitudes za flow rates zilizotolewa katika m3·s−1 kwenye table zinahitaji unit au scale explanation zikizingatiwa pamoja na stated pipe diameters, lengths na pump curve. Kwa mfano, katika source equation \(H_{\mathrm{pump}}=0{,}15-15Q^2\), ukitumia \(Q=0{,}0976\), pump head ni only approximately 0,0071 m. Kwa hiyo flow-rate unit au decimal notation katika Table 3 inahitaji kuangaliwa na authors.
Third validation: multiple branches na energy mixing
Final analytical problem ni more complex pipe network yenye two successive flow splits, two recombinations na localized 30 MW heat source. Inlet velocity ni 1 m·s−1, elbow loss coefficient 0,9 na T-junction loss coefficient 1,2.
Figure 5 kwenye page 18 inaonyesha temperature katika heated upper branch ikiongezeka na temperature ikisambazwa kwenda other branches baada ya mixing junctions. Figure 6 kwenye page 19 inaonyesha different velocity magnitudes katika parallel branches.
| Variable | Analytical solution | FFN result | Relative error |
|---|---|---|---|
| Lower-branch flow rate \(Q_b\) | 0,4367 m3·s−1 | 0,4369 m3·s−1 | %0,06 |
| Heated-branch flow rate \(Q_h\) | 0,1442 m3·s−1 | 0,1444 m3·s−1 | %0,16 |
| Unheated upper sub-branch flow rate \(Q_l\) | 0,2045 m3·s−1 | 0,2047 m3·s−1 | %0,08 |
| Outlet temperature | 309,11 K | 309,20 K | %0,03 |
Comparison hii inaonyesha si total flow pekee, bali pia distribution katika parallel branches na outlet enthalpy inayotokana na mixing streams za different temperatures zina agreement na analytical result.
Ni reactor circuit gani ilijengwa kwa TRACE comparison?
Katika system-level assessment, simplified two-loop primary circuit inayofanana na pressurized-water reactor ilitumika. Model ina:
- Reactor core,
- Lower na upper plena,
- Two symmetric hot legs,
- Two steam-generator representations,
- Two main coolant pumps,
- Two cold legs,
- Cross-connection kati ya plena
. Figure 7 kwenye page 20 inaonyesha TRACE nodalization pamoja na OpenFOAM/FFN pipe network.
Distributed friction ilihesabiwa kwa Churchill correlation, huku convection katika core na steam generators ikihesabiwa kwa El-Genk correlation. Core iliwakilishwa na lumped fuel-rod model yenye one-dimensional radial heat conduction ndani ya fuel na conjugate heat transfer kwenda coolant.
Nominal conditions za main pumps ni:
| Pump parameter | Value |
|---|---|
| Nominal volumetric flow rate | 4,65 m3·s−1 |
| Nominal angular speed | 124,2 rad·s−1 |
| Rotor moment of inertia | 590 kg·m2 |
Figure 8 kwenye page 21 inaonyesha homologous pump-head na torque curves zilizotumika kwa pamoja katika codes zote mbili.
FFN na TRACE zilikubaliana kiasi gani katika steady state?
Four conditions ambapo pumps zote mbili zilifanya kazi kwa nominal na half speeds zililinganishwa:
| Right pump speed | Left pump speed | TRACE core flow | FFN core flow | Difference |
|---|---|---|---|---|
| %100 | %100 | 4929,2 kg·s−1 | 4817,5 kg·s−1 | −%2,3 |
| %100 | %50 | 3538,5 kg·s−1 | 3459,2 kg·s−1 | −%2,2 |
| %50 | %100 | 3538,5 kg·s−1 | 3459,2 kg·s−1 | −%2,2 |
| %50 | %50 | 2454,5 kg·s−1 | 2399,1 kg·s−1 | −%2,3 |
FFN ilipredict total flow approximately %2,2–2,3 lower katika operating points zote. Watafiti wanaunganisha systematic difference hii si na mass-conservation error katika connection method, bali na differences katika friction na local hydraulic-loss models kati ya codes mbili.
Katika asymmetric cases ambapo one pump ilikuwa half speed, loop-flow ratios pia zilikuwa close. Right pump ikiwa %100 na left pump %50, TRACE ilihesabu right/left flow ratio 3,714 na FFN 3,870.
Zero-power pump-trip transient
Katika first transient, core power haikutumika na only hydraulic behavior ilichunguzwa. Right pump ilipotrip, speed yake ilipungua gradually kutokana na rotor inertia, right-loop flow ikapungua, na pressure difference iliyoundwa na operating left loop kupitia common core ikasababisha flow reversal katika right loop.
Katika Figure 9 kwenye page 23, TRACE solid lines na FFN dashed lines zinaoverlap kwa kiasi kikubwa wakati wa pump rundown, passage through zero flow na reverse flow. Study haijatoa numerical error norm au maximum-difference value katika transient yote; assessment ya “very good agreement” inategemea hasa visual comparison ya curves.
Nominal-power pump-trip transient
Katika second transient, core heat generation pia ilijumuishwa. Right pump ilipotrip, coolant flow ilipungua, core outlet temperature ikaongezeka, na resulting density differences zikaongeza buoyancy. Hivyo natural-circulation effect iliingia katika circuit behavior pamoja na hydraulic resistance na pump inertia.
Figure 10 kwenye page 24 inalinganisha four variables:
- Right-loop mass flow rate,
- Right-pump rotational speed,
- Core outlet temperature,
- Heat-transfer coefficient ya right steam generator.
FFN na TRACE curves zilibaki close wakati wa pump rundown, flow-direction reversal, core-temperature rise na transition to mixed-convection regime. Hata hivyo, study haijaripoti integral error, time lag, peak-temperature difference au uncertainty interval kwa transient hii.
Main conclusions zinazoungwa mkono na study
- One-dimensional pipes zenye different cross-sections na orientations ziliunganishwa ndani ya OpenFOAM huku mass flow ikihifadhiwa.
- Flow split na merging katika junctions zenye three au more branches zilitoa errors chini ya %1 relative to analytical solutions.
- Mass-flow-weighted enthalpy mixing ilireproduce outlet temperature katika multi-branch network kwa %0,03 error.
- Pressure-jump model iliwakilisha dynamic pressure, local loss, pump na hydrostatic effects katika same boundary condition.
- Four-region pump model iliwakilisha transitions kwenda low-flow, zero-flow na reverse-flow regions baada ya pump trip.
- Python API iligenerate connection geometries na boundary conditions za pipe networks automatically.
- Katika simplified reactor circuit, FFN na TRACE steady-state core flow rates zilikubaliana ndani ya %3.
- Katika zero- na nominal-power pump-trip transients, FFN ilifuatilia closely TRACE flow na temperature trends.
Study haionyeshi nini?
- Framework haimodel two-phase flow, boiling, condensation, void fraction au critical heat flux.
- Hakuna pipe break, loss-of-coolant accident au steam-generating reactor transient iliyovalidateiwa.
- Results hazikulinganishwa na experimental thermal-hydraulic facility.
- TRACE comparison si independent physical validation; inaonyesha consistency ya results kutoka two computational codes.
- Only one simplified reactor circuit na two pump-trip scenarios zilichunguzwa.
- Hakuna quantitative transient error norms, uncertainty intervals au sensitivity analysis zilizotolewa.
- Two-way coupling ya one-dimensional na high-fidelity CFD regions katika same problem haikuonyeshwa.
- Software haikuonyeshwa kuwa ready kwa regulatory nuclear-safety analyses.
- Wall time, memory use, parallel scalability na computational cost relative to TRACE hazikuripotiwa numerically.
- Source-code version, commit ID na validation input files zilizotumika katika work hii hazikushirikiwa.
Strengths za study
Main strength ya study ni kwamba haikuishia kutambulisha new solver infrastructure pekee, bali ilijaribu kila fundamental component kwa separate analytical problems. Flow splitting, natural circulation, pump pressure, local losses na temperature mixing zilitathminiwa kwa isolated validation problems.
Kushughulikia junctions kwa separate physical rules kwa mass, momentum na energy ni muhimu methodologically. Hasa dynamic calculation ya enthalpy mixing kulingana na flow direction inaendana na flow reversals katika transients.
Kuwepo kwa four-region homologous model pamoja na polynomial pump model kunaipanua framework beyond normal operating point hadi pump-trip na reverse-flow conditions.
Python API inatoa practical software layer inayoweza kupunguza orientation, connection na transformation errors zinazoweza kutokea wakati complex circuits zinajengwa manually.
Main limitations na points zinazohitaji clarification
Most important scope limitation ni single-phase formulation. Katika nuclear-plant safety analyses, boiling, evaporation, condensation, two-phase pressure loss na interphase transfer ni decisive katika many events. Authors wanasema two-phase extension inaendelezwa, lakini version hii haitoi two-phase equation wala result.
Analytical validations zimetoa low errors; lakini detailed mesh na time-step convergence tables zimeelezwa only to limited extent. Katika first problem, 1090 cells zinasemwa sufficient kwa spatial convergence, lakini results katika different mesh resolutions hazijaonyeshwa.
Katika TRACE comparison, numerical differences zimetolewa kwa steady state, lakini transient assessment inategemea kwa kiasi kikubwa graphs kuoverlap. Metrics kama peak-value difference, root-mean-square error au event-timing error hazijatolewa.
Unit au scale ya closed-loop flow rates zilizoripotiwa katika Table 3 inapaswa kuangaliwa upya pamoja na geometry na pump curve katika same section. Point hii haiondoi closeness ya analytical na FFN values; lakini ikiwa physical unit ya values zote mbili imetolewa incorrectly, validation result inahitaji kutafsiriwa upya.
Ingawa study ni OpenFOAM-based software development, haitoi open-access information kuhusu software version, source-code tag, input files au automated test package. Hii inalimit independent reproduction.
Mbinu na Matokeo ya Utafiti
Scope ya model
| Feature | Implementation katika study |
|---|---|
| Base platform | OpenFOAM-based foamForNuclear |
| Dimension | One-dimensional pipe na system network |
| Flow regime | Single phase |
| Conserved quantities | Mass, momentum na energy/enthalpy |
| Geometric components | Pipes, elbows, junctions, branches na diameter changes |
| Hydraulic components | Local losses, pumps na time-dependent pressure jumps |
| Thermal components | Heat sources, heat sinks na ideal mixing katika junctions |
| Circuit generation | foamForNuclear Python API |
Developed software components
| Component | Task | Main limitation |
|---|---|---|
branchConnector1D | Create cyclic au zero-gradient connection on field basis | Contribution ya other branches katika multi-branch junction inahesabiwa semi-implicitly |
| Mass-conserving velocity coupling | Preserve \(\alpha\rho AU\) across different areas | Applied to current single-phase formulation |
pressureCoupling1D | Apply dynamic-pressure, loss, pump na hydrostatic jump | Accuracy ya loss coefficients inategemea user inputs |
| Polynomial pump model | Create \(h(Q)\) curve at normal operating point | Cannot represent zero na reverse flow |
| Homologous pump model | Represent forward/reverse rotation na forward/reverse flow | Requires correct homologous head na torque tables |
| Inertial pump trip | Calculate rotor speed after motor torque is removed | Inertia na friction model lazima zifafanuliwe na user |
| Enthalpy mixing | Weight streams entering junction by mass flow | Assumes ideal na instantaneous mixing inside junction |
| Python pipe API | Automatically create mesh, orientation na boundary conditions | Full input files na version identity hazikushirikiwa katika study |
Analytical validation matrix
| Validation problem | Feature tested | Mesh au condition | Highest reported error |
|---|---|---|---|
| Flow split na merge | Junction mass conservation na friction loss | 1090 cells, 570 L·s−1 inlet | %0,33 |
| Pumped closed loop | Pump head na distributed losses | 320 cells, four pipes za different diameters | <%0,05 |
| Natural-circulation loop | Temperature-dependent density na buoyancy | 300–350 K temperature difference | <%0,05 |
| Multiple branching na heating | Flow distribution na enthalpy mixing | 30 MW heat source | %0,16 |
Main features za TRACE comparison model
| Subsystem | Modeling approach |
|---|---|
| Reactor circuit | Two-loop simplified PWR-like primary system |
| Core | Lumped fuel rod na one-dimensional radial conduction |
| Heat transfer | Conjugate fuel–coolant transfer na El-Genk correlation |
| Steam generators | Heat sinks with fixed secondary-side temperature |
| Friction | Churchill correlation na component-specific roughness |
| Pumps | Four-region homologous head na torque curves |
| Transient event | Trip ya right main coolant pump |
| Compared scenarios | Zero power na nominal power |
Main quantitative findings
| Finding | Result | Interpretation limit |
|---|---|---|
| Analytical hydraulic validation | Reported errors %0,33 au lower | Selected idealized problems |
| Energy mixing | Outlet-temperature error %0,03 | Constant-property na ideal-mixing assumption |
| TRACE nominal steady state | FFN core flow %2,3 lower | Correlations katika codes mbili si completely identical |
| Asymmetric pump conditions | Total-flow difference below %3 | Only four speed combinations |
| Zero-power pump trip | Flow rundown na reversal visually close | No quantitative transient error metric provided |
| Nominal-power transient | Flow, pump speed, temperature na heat-transfer curves close | Single transient event na single reactor configuration |
Reproducibility na numerical reliability
Study inatoa fundamental equations, boundary conditions, analytical comparison values na main modeling assumptions za reactor circuit kwa detail. Hii inarahisisha methodological evaluation.
Hata hivyo, information ifuatayo haipo au ni limited kwa independent reproduction:
- OpenFOAM na foamForNuclear versions zilizotumika,
- Source-code commit au release tag,
- TRACE input model na version,
- Diameters, lengths, roughness na loss coefficients za all pipes,
- Complete list ya time step na linear-solver tolerances,
- Transient mesh na time-step sensitivity,
- Wall time, CPU use na parallel scalability,
- Quantitative transient error measures kwa TRACE comparison,
- Input files zinazowezesha analytical tests kuendeshwa automatically.
Kwa kuwa hii ni deterministic computational study, experimental replication au classical statistical-significance test haitarajiwi. Badala yake numerical uncertainty inapaswa kutathminiwa kupitia mesh resolution, time step, iteration tolerance, closure correlations na sensitivity kwa input parameters. Analyses hizi zimewasilishwa only to limited extent katika current study.
Dokezo la Chanzo na Mbinu
Full original title ya study: Development of a One-Dimensional Thermal-Hydraulics Framework for Nuclear Reactor Analysis in OpenFOAM
Authors na order: Giovanni Nervi, Alessandro Scolaro, Thomas Guilbaud, Matteo Bettiol, Mathieu Hursin.
Equal contribution au co-first authorship: Hakuna equal-contribution au co-first-authorship statement katika study.
Corresponding au contact author: Uploaded study haina star au email inayobainisha corresponding author. SSRN record page inamtaja Giovanni Nervi kama “Contact Author”.
Institutions: Kwa Giovanni Nervi na Mathieu Hursin: Laboratory for Reactor Physics and Systems Behaviour, École Polytechnique Fédérale de Lausanne, Lausanne, Switzerland; kwa Alessandro Scolaro: Laboratory for Reactor Physics and Thermal-Hydraulics, Paul Scherrer Institute, Villigen, Switzerland; kwa Thomas Guilbaud: Transmutex SA, Geneva, Switzerland; kwa Matteo Bettiol: Laboratory for Simulation and Modelling, Paul Scherrer Institute, Villigen, Switzerland.
DOI: 10.2139/ssrn.7197741. DOI hii ni ya SSRN preprint record na si DOI ya peer-reviewed journal article.
Journal au conference: Hakuna verified peer-reviewed journal au conference publication iliyotajwa kwa reviewed version.
Publication platform: SSRN.
Original publisher: Hakuna verified final journal publisher kwa study. SSRN ndiyo platform ambayo preprint imesambazwa.
Publication year: 2026.
Source type: Research preprint katika computational nuclear engineering, single-phase system thermal hydraulics na scientific-software development.
Peer-review status: Study haijapitia peer review. Kila page ina warnings “This preprint research paper has not been peer reviewed” na “Preprint not peer reviewed”.
Official links:Official SSRN record page na SSRN DOI link.
Funding: Study iliungwa mkono na ESFR-SIMPLE project iliyofadhiliwa na European Union chini ya Grant Agreement 101059543.
Author contributions: Giovanni Nervi: conceptualization, data curation, formal analysis, investigation, methodology, software, validation, visualization na original draft; Alessandro Scolaro: conceptualization, formal analysis, methodology, project administration, software, supervision na review; Thomas Guilbaud: conceptualization, software na review; Matteo Bettiol: conceptualization, software, validation na review; Mathieu Hursin: conceptualization, project administration na supervision.
Data and code access: Reviewed study haitaji open repository ya validation input files, software commit ID, TRACE model au raw numerical results. Example Python API code imetolewa; hata hivyo example hii haitoshi kureproduce entire study independently.
Makala hii ya Kiswahili imeandaliwa kwa kutegemea method descriptions, 27 fundamental equations, four main result tables, T-junction na circuit schematics, temperature na velocity distributions, TRACE nodalization, homologous pump curves na pump-trip transient graphs za uploaded 27-page study. Hakuna new numerical result isiyokuwepo katika study au scientific-performance finding kutoka external source iliyoongezwa. External-source use imewekewa kikomo na bibliographic verification ya SSRN record, DOI, contact author na institutional identity information.
Main limitations za study ni lack of peer review, model kuwa restricted to single-phase flows, absence of physical experimental validation, comparison ya only one simplified reactor circuit na TRACE, absence of quantitative transient error measures, kutoripoti performance na scalability measurements, lack of reproducible code na input files, na hitaji la clarification ya flow-rate unit au scale katika closed-loop validation.
Results zinaonyesha developed framework inatumia conservation equations correctly katika selected single-phase analytical tests na inaweza kutoa trends similar to TRACE katika simplified pump-trip transients. Findings hazimaanishi licensing adequacy kwa real nuclear facility, two-phase accident analysis, experimental validation au regulatory approval for use.

Acha maoni
Anwani yako ya barua pepe haitachapishwa. Sehemu za lazima zimewekewa alama ya *