Jimenez, J.

Sanchez, V.

Trost, N.

English

KITopen

1NURESAFE WP1.1 TESTBED FOR INTEGRATED COUPLING ANDUNCERTAINTY QUANTIFICATION METHODSJ. Jimenez, N. Trost, V. SanchezPresented by J. JimenezSUBCHANFLOW sensitivity analysis with URANIEvictor.sanchez@kit.eduor javier.jimenez@kit.eduNURESAFE-SP1 3rd Meeting, January 30-31, UPM, Madrid2Outline    Summary from the last meeting in HZDR• SUBCHANFLOW code• URANIE software• BFBT benchmark used for the analysis• URANIE measurement scripts• Results for BFBT single phase pressure drop analysis• Results for BFBT two phase pressure drop analysis New scripts for transient analysis Application to the O2 -1999 Feedwater transient Conclusion and OutlookNURESAFE-SP1 3rd Meeting, January 30-31, UPM, Madrid3SUBCHANFLOW – a KIT TH sub-channel code Sub-channel fuel assembly simulations Fuel assembly or “channel wise” whole core simulations Steady state and transient solution Available fluids: water, lead, lead-bismuth sodium, helium, air Works with SI units Flexible geometry definition (square, hexagonal) Completely programmed in Fortran 95 Can be used with WINDOWS or LINUX Dynamic memory management Modular structure, keyword and table oriented inputNURESAFE-SP1 3rd Meeting, January 30-31, UPM, Madrid4SUBCHANFLOW FeaturesNational InternationalHZDR Dresden Simulation of single and two phase flow testsUse of SUBCHANFLOW at KIT EU THINS Project (Lead/Bismuth) EU CDT Project (Lead) Multi-physics coupling with • MCNP5, SERPENT• DYN3D-SP3 (transport)• COBAYA3 (Diffusion) International Benchmarks e.g. BFBT, PSBT Code improvement and validation:• (Internships, Bachelor, MSc, PhD)EU JRC IE Petten Gen-IV reactors SUBCHANFLOWKIT Subchannel CodeRWTH Aachen INBK: Simulation of gas cooled reactors for transmutation (AGATE Project)TÜV Nord LWR applications of SUBCHANFLOW (Testing Phase) CIEMAT/ Nuclear Safety Department SCF / MCNPX Coupling DNC / Delft EU HPMC Project SCF coupling with time-dependent Monte Carlo    UPM/ Nuclear Engineering Chair EU ESFR Project Coupling with COBAYA3  NURESAFE-SP1 3rd Meeting, January 30-31, UPM, Madrid5URANIE Software URANIE is a software dedicated to uncertainty and optimization. It allows to perform studies on uncertainty propagation, sensitivity analysis or model calibration in an integrated environment, based on ROOT, a software developed at CERN for particle physics data analysis. As a result, URANIE benefits from the numerous features of ROOT, among which:• a C++ interpreter (CINT)• a Python interface (PyROOT)• access to SQL databases• many advanced data visualization features URANIE training course attended 2-4th April 2013 in Saclay.NURESAFE-SP1 3rd Meeting, January 30-31, UPM, Madrid6BFBT benchmark Phase II, Exercise 0 NUPEC BFBT experimental facilityNURESAFE-SP1 3rd Meeting, January 30-31, UPM, Madrid Steady state pressure drop measurements where used for the analysis7Sensitivity study parameters BFBT single phase P70001 – P70012 cases for pressure drop analyses BFBT two phase P60001 – P60012 cases for pressure drop analyses 100 SUBCHANFLOW simulation runs per test scenario using URANIE Normal and Uniform law for boundary parameter distribution, see Table belowNURESAFE-SP1 3rd Meeting, January 30-31, UPM, MadridNo. Parameter Range Distribution1 Outlet pressure ± 1.0 % Normal2 Mass flow rate ± 1.0 % Normal3 Inlet temperature ± 1.5 K Uniform4 Power (only P600xx) ± 1.5 % Normal8Script for SUBCHANFLOW runvoid BFBT_P6x(Int_t nS = 100) {TDataServer *tds = new TDataServer("tdsSCF", "BFBT_P6x");TNormalDistribution *tnp = new TNormalDistribution("OutletPressure", 7.16e6, 71600.0);  // +‐1%TNormalDistribution *tnf = new TNormalDistribution("MassFlowRate", 5.61, 0.0561);       // +‐1%TUniformDistribution *tut = new TUniformDistribution("InletTemperature", 276.3, 279.3); // +‐1.5KTNormalDistribution *tnl = new TNormalDistribution("Power", 1.951e6, 29265.0);          // +‐1.5%tnp‐>setBounds(7.0884e6, 7.2316e6);tnf‐>setBounds(5.5539, 5.6661);tnl‐>setBounds(1.921735e6, 1.980265e6);tds‐>addAttribute(tnp);tds‐>addAttribute(tnf);tds‐>addAttribute(tut);tds‐>addAttribute(tnl);TString sFileName = TString("input.txt");...Specify the number of simulation runsCreate DataServer to hold all informationsSpecify the distribution for each parameterSet upper and lower bound to avoid non-physical valuesFill up the DataServerSpecify the input file of SUBCHANFLOW to be parsed by URANIE and replaced by the sampled random dataNURESAFE-SP1 3rd Meeting, January 30-31, UPM, Madrid9Script for SUBCHANFLOW runNURESAFE-SP1 3rd Meeting, January 30-31, UPM, Madrid...tds‐>getAttribute("OutletPressure")‐>setFileKey(sFileName, "exit_pressure");tds‐>getAttribute("MassFlowRate")‐>setFileKey(sFileName, "inlet_flow_rate");tds‐>getAttribute("InletTemperature")‐>setFileKey(sFileName, "inlet_temperature");tds‐>getAttribute("Power")‐>setFileKey(sFileName, "total_power");TSampling *sampling = new TSampling(tds, "lhs", nS);sampling‐>generateSample();TOutputFileKey *fout = new TOutputFileKey("output.txt");TAttribute *avgpres = new TAttribute("total_ax_pres_loss_uranie");avgpres‐>setDefaultValue(‐200.0);fout‐>addAttribute(avgpres);TCode *mycode = new TCode(tds, "SUBCHANFLOW");mycode‐>addOutputFile(fout);TLauncher *tlch = new TLauncher(tds, mycode);tlch‐>setSave();tlch‐>setClean();tlch‐>setWorkingDirectory(gSystem‐>Getenv("PWD") + TString("/tmpUranie/SCF"));tlch‐>setVarDraw("MassFlowRate:total_ax_pres_loss_uranie","","");...Give URANIE the strings of SUBCHANFLOW input to be changedGenerate the random dataSpecify the output file to extract the data for post-processingTell URANIE which value has to be extracted from SCF output and initialize it with a non-physical value for easier error detectionInitiate a SUBCHANFLOW simulation run10Script for SUBCHANFLOW post-processingNURESAFE-SP1 3rd Meeting, January 30-31, UPM, MadridTCanvas *Canvas = new TCanvas("c1", "Graph",5,64,1270,667);c1‐>Divide(2, 2);c1‐>cd(1);tlch‐>run();c1‐>cd(3);tds‐>draw("OutletPressure:total_ax_pres_loss_uranie");c1‐>cd(2);tds‐>draw("InletTemperature:total_ax_pres_loss_uranie","","");c1‐>cd(4);tds‐>draw("Power:total_ax_pres_loss_uranie","","");TCanvas *Canvas2 = new TCanvas("c2", "Graph",5,64,1270,667);tds‐>draw("total_ax_pres_loss_uranie");tds‐>exportData("BFBT_P6x_Sampling.dat");}Visualize the output of URANIEPlot a histogram of the output dataExport the output as well as the sampled random numbers into a file for post-processing.11Script for SUBCHANFLOW statisticsNURESAFE-SP1 3rd Meeting, January 30-31, UPM, Madrid{using namespace URANIE::DataServer;using namespace URANIE::Sampler;TDataServer *tds = new TDataServer();tds‐>fileDataRead("BFBT_P6x_Sampling.dat");tds‐>addAttribute("OutletPressure", "OutletPressure");tds‐>addAttribute("MassFlowRate", "MassFlowRate");tds‐>addAttribute("InletTemperature", "InletTemperature");tds‐>addAttribute(„Power", „Power");tds‐>addAttribute("AxialPressureLoss", "total_ax_pres_loss_uranie");tds‐>computeStatistic();tds‐>computeCorrelationMatrix()‐>Print();std::cout << "mean: " << tds‐>getAttribute("AxialPressureLoss")‐>getMean()    << std::endl;std::cout << "min: "  << tds‐>getAttribute("AxialPressureLoss")‐>getMinimum() << std::endl;std::cout << "max: "  << tds‐>getAttribute("AxialPressureLoss")‐>getMaximum() << std::endl;std::cout << "std: "  << tds‐>getAttribute("AxialPressureLoss")‐>getStd()     << std::endl;}Read in an URANIE file that contains sampled random dataFill up the DataServer structureCompute the statistics and print the correlation matrixPrint mean, min, max, and standard deviation12Visualization for total axial pressure lossNURESAFE-SP1 3rd Meeting, January 30-31, UPM, MadridCorrelation of massflow rate, inlet temperature and outlet pressure with the total axial pressure loss.Correlation of massflow rate, inlet temperature and outlet pressure with the total axial pressure loss. Correlation matrix can be also extracted URANIE visualization output for P70001 single phase BFBT benchmark13BFBT Two Phase P60001 – P60012 Sensitivity coefficient for BFBT Two Phase P60001 – P60012 pressure loss analysisNURESAFE-SP1 3rd Meeting, January 30-31, UPM, Madrid-1,0-0,8-0,6-0,4-0,20,00,20,40,60,81,001 02 03 04 05 06 07 08 09 10 11 12Sensitivity coefficientP600xxOutlet Pressure Massflowrate Inlet Temperature Power14NURESAFE-SP1 3rd Meeting, January 30-31, UPM, Madrid0,700,750,800,850,900,951,001,051,101 2 3 4 5 6 7 8 9 10 11 12Relative pressureP600xxPressure minimum maximum Relative pressure for BFBT Two Phase P60001 – P60012 pressure loss analysis to the experimental dataBFBT Two Phase P60001 – P6001215New scripts for transient analysis Up to this point it was reported in the last meeting Study of the pump trip and turbine trip transient scenario of the BFBT benchmark Range of variation ± 1.0 % only for academic purposesNURESAFE-SP1 3rd Meeting, January 30-31, UPM, MadridNo. Parameter Range Distribution1 Outlet pressure ± 1.0 % Normal2 Mass flow rate ± 1.0 % Normal3 Inlet temperature ± 1.0 % Uniform4 Power ± 1.0 % Normal5 Cladding Wall Roughness ± 1.0 % Normal6 Spacer grid pressure drop coefficient ± 1.0 % Normal7 Flow Area ± 1.0 % Normal8 Wetted Perimeter ± 1.0 % Normal16New scripts for transient analysisNURESAFE-SP1 3rd Meeting, January 30-31, UPM, Madrid New scripts for transient analyses and post-processing Around 500 lines of URANIE code More details can be found in:• D11.22 - Report on FLICA UQ results for BWR ATWS analysis (t0+18)DRAFT AVAILABLE17BFBT 4102 - Turbine trip results Comparison of results from several codes with the measured dataNURESAFE-SP1 3rd Meeting, January 30-31, UPM, Madrid18BFBT 4102 - Turbine trip results Sensitivity coefficient of the void fractionNURESAFE-SP1 3rd Meeting, January 30-31, UPM, Madrid19BFBT 4102 - Turbine trip results Sensitivity coefficient of the axial pressure dropNURESAFE-SP1 3rd Meeting, January 30-31, UPM, Madrid20BFBT 4102 - Turbine trip results Void fraction results at the three different elevations.NURESAFE-SP1 3rd Meeting, January 30-31, UPM, Madrid21BFBT 4102 - Turbine trip results Total axial pressure dropNURESAFE-SP1 3rd Meeting, January 30-31, UPM, Madrid22Application to the O2 -1999 FW transient Power oscillation during the event (feedwatertransient) Boundary conditions taken from TRACE/PARCS calculation (KIT model with 444 channels)  Modeling the O2 core with SCF and 444 channelsNURESAFE-SP1 3rd Meeting, January 30-31, UPM, Madrid23Boundary Conditions applied The next boundary conditions were introduced into SCF for the simulation of the oscillations.  They have been extracted from the TRACE5p3 results• Power, inlet temperature, pressure, mass flow rate.NURESAFE-SP1 3rd Meeting, January 30-31, UPM, Madrid24Results from SCF for the full transient Void (bundle average at outlet)NURESAFE-SP1 3rd Meeting, January 30-31, UPM, Madrid25Results from SCF for the full transient Outlet Temperature in Celsius for bundle averageNURESAFE-SP1 3rd Meeting, January 30-31, UPM, Madrid26Results from SCF for the full transient Pressure drop in the core for bundle averageNURESAFE-SP1 3rd Meeting, January 30-31, UPM, Madrid27Sensitivity study Zooming and running of a small portion of the transient (12s), when the oscillations start. Sensitivity analysis with parameters taken from the NURESAFE benchmark specifications (D13.11)NURESAFE-SP1 3rd Meeting, January 30-31, UPM, MadridNo. Parameter Range Distribution1 Outlet pressure ± 0.5 % Uniform2 Mass flow rate ± 0.5 % Uniform3 Inlet temperature ± 2.0 % Normal4 Power ± 0.75 % Normal5 Cladding Wall Roughness ± 30.0 % Normal6 Spacer grid pressure drop coefficient ± 5.0 % Uniform7 Gap Conductance ± 35.0 % Uniform8 Fuel Conductivity ± 10.0 % Uniform9 Cladding Conductivity ± 6.25 % Uniform28Results in the zooming area Sensitivity of the void fractionNURESAFE-SP1 3rd Meeting, January 30-31, UPM, Madrid29Results in the zooming area Sensitivity of the axial pressure dropNURESAFE-SP1 3rd Meeting, January 30-31, UPM, Madrid30Results in the zooming area Void fraction resultsNURESAFE-SP1 3rd Meeting, January 30-31, UPM, Madrid31Results in the zooming area Axial pressure drop resultsNURESAFE-SP1 3rd Meeting, January 30-31, UPM, Madrid32Conclusions and Outlook During the first 12 months of the project, investigations on the use of URANIE platform for sensitivity analyses have been conducted. As a proof of principle, studies have been conducted using the SUBCHANFLOW code (similar to FLICA4). In June 2013 we reported in Dresden the scripts to analyze steady state scenarios.  Now the scripts have been extended also to transient simulations. The scripts can be extrapolated to any code with input text file: CTF, FLICA4, DYN3D, COBAYA3, etc, ...NURESAFE-SP1 3rd Meeting, January 30-31, UPM, Madrid33Conclusions and Outlook Work for the next 6 months: Finalizing D11.22 for (t0+18)• The FLICA4 input deck for O2 has been prepared at KIT for comparison studies• FLICA4 will not be used for the coupled computation analyses within WP1.3, only CTF• Therefore, to us it is a non-sense to use FLICA4 within WP1.1 if it is not going to be applied within WP1.3 Hence we will propose to change the title and content of D11.22 to Report on CTF UQ results for BWR ATWS analysis.DISCUSSIONNURESAFE-SP1 3rd Meeting, January 30-31, UPM, Madrid34THANKS FOR YOUR ATTENTIONNURESAFE-SP1 3rd Meeting, January 30-31, UPM, Madrid

SUBCHANFLOW sensitivity analysis with URANIE

https://publikationen.bibliothek.kit.edu/230095170/3816510

SUBCHANFLOW sensitivity analysis with URANIE

Abstract

Similar works

Full text

Available Versions

KITopen