With COMSOL, That Time is Now Presented by: James D. Freels, - - PowerPoint PPT Presentation

▶
with comsol that time is now
SMART_READER_LITE
LIVE PREVIEW

With COMSOL, That Time is Now Presented by: James D. Freels, - - PowerPoint PPT Presentation

May 25, 2023 10 likes •222 views

With COMSOL, That Time is Now Presented by: James D. Freels, Ph.D. Senior Research Staff Oak Ridge National Laboratory Oak Ridge, TN 37831-6392 freelsjd@ornl.gov Presented to: COMSOL Conference 2014 Boston Marriott Newton Thursday,

slide-1
SLIDE 1

ORNL is managed by UT-Battelle for the US Department of Energy

With COMSOL, That Time is Now

Presented by:

James D. Freels, Ph.D.

Senior Research Staff Oak Ridge National Laboratory Oak Ridge, TN 37831-6392 freelsjd@ornl.gov

Presented to:

COMSOL Conference 2014

Boston Marriott Newton Thursday, October 9, 2014 10:30am Keynote Presentation

slide-2
SLIDE 2

2 James D. Freels, Keynote Presentation, COMSOL Conference 2014 Boston, October 9, 2014

Outline

  • Early simulation tools
  • Career inspiration
  • Engineering choices
  • Projects prior to COMSOL
  • Current COMSOL-

involved projects

  • Inspiring others to use

COMSOL

slide-3
SLIDE 3

3 James D. Freels, Keynote Presentation, COMSOL Conference 2014 Boston, October 9, 2014

1972 Freshmen Engineering Students Learned How To Use a Slide Rule

  • Last year of slide rule course at VPI&SU.
  • The instructor used “teaching” slide rules.
  • Calculators became common the next year.
  • COMSOL could not be envisioned yet.

TI SR-10 Sears ESR

HP-35 HP-45

slide-4
SLIDE 4

4 James D. Freels, Keynote Presentation, COMSOL Conference 2014 Boston, October 9, 2014

1972 Freshmen Engineering Students Learned How To Program in FORTRAN

  • FORTRAN is still in use today.
  • Many legacy FORTRAN codes.
  • Now is a good time to to

replace legacy code with COMSOL models.

300 baud acoustic coupler DEC VT-100 terminal FORTRAN 80-column punch card Keypunch machine

slide-5
SLIDE 5

5 James D. Freels, Keynote Presentation, COMSOL Conference 2014 Boston, October 9, 2014

1977: My Initial Inspiration for a Career in Simulation

  • Senior Lab UTK-NE: ORNL Health-Physics Research Reactor (HPRR)
  • Compared pulse experiments to FORTRAN code
  • Stack of FORTRAN cards about 4-inches deep, tweaked coefficients to match data

These experiments could not be repeated today, but this would have been perfect for the COMSOL application builder !

slide-6
SLIDE 6

6 James D. Freels, Keynote Presentation, COMSOL Conference 2014 Boston, October 9, 2014

Model noding diagram blowdown: results vs test data

1984 : RELAP5 Nuclear Pow er Plant Simulation and Testing of Gamma Thermometer Level Monitoring

  • In response to the TMI accident, need improved level monitoring.
  • An NRC code RELAP5 was used to simulate tests for the new instruments.
  • RELAP5 and similar codes are still used today.
  • COMSOL Pipe Flow Module has demonstrated similar models (Wang et al.)

test and simulation overlay TC response as level uncovers

slide-7
SLIDE 7

7 James D. Freels, Keynote Presentation, COMSOL Conference 2014 Boston, October 9, 2014

Engineering Choices: Technology that Works Requires High Quality

  • 1972 : Nuclear Engineering (NE)
  • 1977 : computer simulation
  • 1979-1982 : utilized codes that needed

a lot of improvements

  • received excellent mentoring
  • 1983-1992 : finite element (FE) methods,

PhD research

  • continue to apply this FE-based simulation in NE
  • used advanced software tools such as LaTeX and Linux
  • 1994-2004 : sought a code that met established goals
  • 2004-2014 : the code became COMSOL
  • 2014-future : expand COMSOL usage to new problems
slide-8
SLIDE 8

8 James D. Freels, Keynote Presentation, COMSOL Conference 2014 Boston, October 9, 2014

1992: PhD Graduation, the Mentoring Finalized

  • Many long discussions together.
  • We envisioned an ideal computer

software environment.

  • Numerical Analysis Digest

announced release 3.0 of COMSOL.

  • After graduation, I asked AJ if he

ever heard of COMSOL.

  • His response: “Sure! My students

do their projects using COMSOL.”

  • I guessed that the number of

incomplete grades (I) for the semester has been reduced !

  • I have often desired to go back

and repeat my own “wine-glass problem” using COMSOL.

  • So, that is how I started with

COMSOL.

  • A. J. Baker

Professor Emeritus UT-Knoxville

Calculated Transverse Momentum for the Ideal-Gas Viscous Solution of the GAMM Double-Throat Nozzle Problem (Fig. 8.84 of Ref #3).

slide-9
SLIDE 9

9 James D. Freels, Keynote Presentation, COMSOL Conference 2014 Boston, October 9, 2014

A Brief Description of the ORNL High Flux Isotope Reactor

  • HFIR is a DOE research

reactor

  • Does not produce power
  • Main product is neutrons
  • HFIR core (shown at the right)

consists of 2 fuel elements, 9.4 kg of 93% enriched U235

  • HFIR has 4 primary missions:
  • 1. neutron scattering,
  • 2. Irradiation materials

testing

  • 3. Isotope production
  • 4. Neutron activation

analysis

slide-10
SLIDE 10

10 James D. Freels, Keynote Presentation, COMSOL Conference 2014 Boston, October 9, 2014

2004: The First COMSOL Project for HFIR at ORNL

  • used COMSOL to design the pressurizer component
  • connects ambient to 20 K supercritical para-hydrogen loop
  • pressure and temperature dependent material properties
  • weak statement provides a thin-shell for the structure
  • COMSOL simulations match the instrument outputs
  • has been operating successfully for over 7 years
  • Could revisit with present versions and computing capability

Installed HFIR cold source pressurizer

COMSOL simulated surface temperature distribution COMSOL simulated internal natural convection flow streamlines

slide-11
SLIDE 11

11 James D. Freels, Keynote Presentation, COMSOL Conference 2014 Boston, October 9, 2014

2011: HFIR Fuel Plate Thermal Deflection Validation

  • Validation case for the HFIR LEU fuel conversion project
  • Demonstrates COMSOL can simulate thermal-structural interaction sufficiently for this project

Jain, et al. ORNL/TM-2012/138 , June 2012

Cheverton, Kelley, ORNL-TM-2325, 1968.

Top-edge Mid-edge Bottom-edge

slide-12
SLIDE 12

12 James D. Freels, Keynote Presentation, COMSOL Conference 2014 Boston, October 9, 2014

  • equation-based modeling, space-time kinetics in 2D axisymmetric geometry.
  • demonstrates COMSOL could be a major tool for modeling reactor physics.

COMSOL Applied to Nuclear Reactor Kinetics, David Chandler et al.

Free mesh + boundary layer mesh Solved for power excursion transient with control element feedback

Fast Flux Epithermal Flux Thermal Flux

neutron flux distribution

slide-13
SLIDE 13

13 James D. Freels, Keynote Presentation, COMSOL Conference 2014 Boston, October 9, 2014

single bare pellet, 2nd irradiation cycle, COMSOL 4.2a, 3D, ¼ pie slice reduced-length bare pellet, 2nd irradiation cycle, COMSOL 4.2a, 3D, ¼ pie slice partially-loaded (8 pellets) prototype production target, 2 irradiation cycles, COMSOL 4.3 fully-loaded prototype production target (52 pellets), COMSOL 4.3, 2D axisymmetric

individual pellet at maximum temperature in stack: stress contour with 10000x deformation

COMSOL Applied to Pu-238 Production for NASA COMSOL Freels, Jain, Hurt, et al.

slide-14
SLIDE 14

14 James D. Freels, Keynote Presentation, COMSOL Conference 2014 Boston, October 9, 2014

2009-Present: COMSOL is providing fundamental support for the conversion of HFIR to LEU fuel

  • Provide best-estimate and safety-basis accurate 3D simulations of:
  • fluid flow into, through, and out of the core flow channels,
  • fuel-plate and coolant heat transfer,
  • structural response due to thermal expansion (TSI),
  • structural response due to fluid interaction (FSI),
  • oxidation of the fuel-plate surface and subsequent geometry

changes,

  • fuel shrink and swell due to radiation effects,
  • hot spot and hot-streak effects due to manufacturing defects,
  • additional physics as necessary.
  • The approach has been a multi-step process that includes:
  • develop physics test bed in 2D to feed results to the 3D models,
  • develop 3D geometry, mesh, and property inputs,
  • perform separate-effects simulations in 3D to V&V physics

goals (heat transfer, FSI, TSI).

  • provide perturbed estimates and/or separate safety-basis models

from the best-estimate models to answer safety-analysis questions (hot spots, buckling potential, reduced flow, decay heat, etc.).

  • combine separate effects models into the all-physics model
  • provide input data for other parts of the HFIR safety analysis

(transients using RELAP5, etc.).

slide-15
SLIDE 15

15 James D. Freels, Keynote Presentation, COMSOL Conference 2014 Boston, October 9, 2014

Typical LEU Fuel COMSOL Safety-Basis Case (from Prashant K. Jain, see poster) Inner Fuel Element under 100 MW Nominal Conditions

volumetric heat source (kW/cm3) clad surface heat flux (W/cm2) clad surface temperature (̊C) plate deflection (mils) coolant pressure (psia)

slide-16
SLIDE 16

16 James D. Freels, Keynote Presentation, COMSOL Conference 2014 Boston, October 9, 2014

Results for Run 1 from Kennedy’s report with a plate thickness of 40 mils. No comb on leading or trailing edge (data courtesy of John Kennedy, University of Missouri, Columbia, Mechanical Engineering Department).

Typical LEU Fuel COMSOL Separate Effects Case (from Franklin G. Curtis) Fluid-Structure Interaction Flat Plate Analysis Results

slide-17
SLIDE 17

17 James D. Freels, Keynote Presentation, COMSOL Conference 2014 Boston, October 9, 2014

Typical LEU Fuel COMSOL Best Estimate Case HFIR Inner Fuel Element Plate – Old Reference Design

wetted surface temperature Tmax cut plane temperature Tmax cut plane height- enabled velocity Tmax axial cut line mesh convergence history of temperature

Tmax normal cut line mesh convergence history of temperature (left) and velocity (right)

SST turbulence model, temperature-dependent properties

slide-18
SLIDE 18

18 James D. Freels, Keynote Presentation, COMSOL Conference 2014 Boston, October 9, 2014

COMSOL has enabled research and employment

  • pportunities for graduate students and

professionals through our group on ORNL projects

  • Kirk T. Lowe , PhD ME dissertation, Bettis Atomic Power.
  • Prashant K. Jain, PhD NE post-doc, full-time ORNL staff.
  • David Chandler, PhD NE dissertation, full-time ORNL staff.
  • Vaibhav Khane, PhD ChE, summer intern at ORNL, Intel.
  • Isaac T. Bodey, PhD ME dissertation, Savannah River

Remediation.

  • Adam R. Travis, MS ME thesis, Alcoa in Pittsburgh.
  • R. Burns Cunningham, MS NE thesis, Enercon Services.
  • Franklin G. Curtis, PhD candidate ME, finishing up soon.
  • Christopher J. Hurt, PhD candidate NE, currently supported.
  • Several intern students have used COMSOL in their research

while supported by ORNL projects.

  • We engage through user groups and workshops to encourage

COMSOL usage at ORNL.

slide-19
SLIDE 19

19 James D. Freels, Keynote Presentation, COMSOL Conference 2014 Boston, October 9, 2014

In the 1980’s, w hen computers became more common in the work place, my colleagues and I would often marvel and say “someday w e w ill be able to solve time-dependent, 3D, Navier- Stokes equations right on our desktops.”

  • w ith COMSOL, that time is now.

Today, I think that, someday, w e w ill be able to solve most any physics w e w ant in full 3D to the resolution that w e desire. Perhaps then w e w ill reexamine the fundamental equations of physics that w e currently hold true.

  • w ith COMSOL, w e w ill be able to test these

new theories and equations.

slide-20
SLIDE 20

20 James D. Freels, Keynote Presentation, COMSOL Conference 2014 Boston, October 9, 2014

Oak Ridge National Laboratory:

Meeting the challenges of the 21st century

w w w.ornl.gov

Thank you for your attention!

slide-21
SLIDE 21

21 James D. Freels, Keynote Presentation, COMSOL Conference 2014 Boston, October 9, 2014

List of References

1.

  • M. W. Rosenthal, An Account of Oak Ridge National Laboratory’s Thirteen Nuclear Reactors,

ORNL/TM-2009/181, March, 2010. 2.

  • R. A. Hedrick, J. D. Freels, et al., Final Report -- Phase I of Arkansas Power and Light Company

Inadequate Core Cooling Monitoring System Program In-core Test Series, Technology for Energy Corporation, R-84-014, August,1984. 3.

  • J. D. Freels, A Taylor Weak Statement Finite Element Algorithm for Real-Gas Compressible Navier-

Stokes Simulation, May, 1992. 4.

  • F. P. Griffin et al., HFIR HB-4 Cold Source Test Data, ORNL/TM-2008/043, April, 2008.

5.

  • P. K. Jain et al., 3D COMSOL Simulations for Thermal Deflection of HFIR Fuel Plate in the

"Cheverton-Kelley" Experiments, ORNL/TM-2012/138, August, 2012. 6.

  • R. D. Cheverton and W. H. Kelley, Experimental Investigation of HFIR Fuel Plate Deflections Induced

by Temperature and Pressure Differentials, ORNL-TM-2325, 1968. 7.

  • D. Chandler et al., COMSOL-based Nuclear Reactor Kinetics Studies at the HFIR, COMSOL

Conference Boston 2011. 8.

  • J. D. Freels et al., Design and Nuclear-Safety Related Simulations of Bare-Pellet Test Irradiations for

the Production of Pu-238 in the High Flux Isotope Reactor using COMSOL, COMSOL Conference Boston 2012. 9.

  • F. G. Curtis et al., Fluid-Structure Interaction Modeling of High-Aspect Ratio Nuclear Fuel Plates

Using COMSOL, COMSOL Conference Boston 2013.

  • 10. D. Wang, Application of COMSOL Pipe Flow Module to Develop a High Flux Isotope Reactor (HFIR)

System Loop Model, COMSOL Conference Boston 2013.

  • 11. C. J. Hurt, Plutonium-238 Production Target Design Studies, Institute of Nuclear Materials

Management 55th Annual Meeting, July 2014.

  • 12. I. T. Bodey, Thermal Hydraulic Characteristics of Fuel Defects in Plate Type Nuclear Research

Reactors, May 2014.

  • 13. A. R. Travis, Simulating High Flux Isotope Reactor Core Thermal-Hydraulics via Interdimensional

Model Coupling, May 2014.

  • 14. P. K. Jain, 3D Multi-physics Analyses to Support Low Enriched Uranium (LEU) Conversion of HFIR,

COMSOL Conference Boston 2014.