Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Dynamic parameter”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,477 records · Page 82
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

A neural-network-enhanced parameter-varying framework for multi-objective model predictive control applied to buildings

Management of the electrical grid is becoming more complex due to the increased penetration of alternative energy generation technologies and a broadening diversity of electric loads. This complexity creates challenges in balancing demand and generation that can increase the potential for grid instabilities. One effective way to address this issue is to leverage previously unexploited demand flexibility through advanced control strategies. In this work, we propose an advanced control method, called adaptive neural parameter-varying model predictive control (ANPV-MPC), to control the temperature and energy consumption of a building via its Heating, Ventilation, and Air Conditioning system. ANPV-MPC combines key ideas in parameter-varying control, adaptive control, and online learning strategies to bridge the gap between computationally efficient linear model predictive control and more accurate nonlinear model predictive control. The novelty in ANPV-MPC is the use of a physics-inspired Bayesian neural network to estimate the coefficients of the parameter-varying linear control model. The Bayesian neural network additionally provides uncertainty estimates, triggering online training to capture evolving building system conditions. We show that ANPV-MPC can approximate the building system dynamics with a 28.39% higher accuracy than traditional linear model predictive control, resulting in 36.23% better control performance without increasing complexity of the optimal control problem. ANPV-MPC also adapts in real time to previously unseen conditions using online learning, further improving its performance.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2024

On time-dependent backgrounds in 1 + 1 dimensional string theory

In perturbative string theory, one is generally interested in asymptotic observables, such as the S-matrix in flat spacetime, and boundary correlation functions in anti-de Sitter spacetime. However, there are backgrounds in which such observables do not exist. We study examples of such backgrounds in 1 + 1 dimensional string theory. In these examples, the Liouville wall accelerates and can become spacelike in the past and/or future. When that happens, the corresponding null infinity, at which the standard scattering states are defined, is shielded by the Liouville wall. We compute scattering and particle production amplitudes in these backgrounds in the region in parameter space where the wall remains timelike, and discuss the continuation of this picture to the spacelike regime. We also discuss the physics from the point of view of the dynamics of free fermions in backgrounds with a time-dependent Fermi surface.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Linear two-pool models are insufficient to infer soil organic matter decomposition temperature sensitivity from incubations

Abstract Terrestrial carbon (C)-climate feedbacks depend strongly on how soil organic matter (SOM) decomposition responds to temperature. This dependency is often represented in land models by the parameter Q 10 , which quantifies the relative increase of microbial soil respiration per 10 °C temperature increase. Many studies have conducted paired laboratory soil incubations and inferred “active” and “slow” pool Q 10 values by fitting linear two-pool models to measured respiration time series. Using a recently published incubation study (Qin et al. in Sci Adv 5(7):eaau1218, 2019) as an example, here we first show that the very high parametric equifinality of the linear two-pool models may render such incubation-based Q 10 estimates unreliable. In particular, we show that, accompanied by the uncertain initial active pool size, the slow pool Q 10 can span a very wide range, including values as high as 100, although all parameter combinations are producing almost equally good model fit with respect to the observations. This result is robust whether or not interactions between the active and slow pools are considered (typically these interactions are not considered when interpreting incubation data, but are part of the predictive soil carbon models). This very large parametric equifinality in the context of interpreting incubation data is consistent with the poor temporal extrapolation capability of linear multi-pool models identified in recent studies. Next, using a microbe-explicit SOM model (RESOM), we show that the inferred two pools and their associated parameters (e.g., Q 10 ) could be artificial constructs and are therefore unreliable concepts for integration into predictive models. We finally discuss uncertainties in applying linear two-pool (or more generally multiple-pool) models to estimate SOM decomposition parameters such as temperature sensitivities from laboratory incubations. We also propose new observations and model structures that could enable better process understanding and more robust predictive capabilities of soil carbon dynamics.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Calibrating Synchronous-Generator-Interfaced DG Models in Microgrids Using Multiple Event Data

Microgrids, consisting of distributed generators (DGs), loads, and energy storage systems, play an important role in future smart grids. In order to evaluate how microgrids operate and their dynamic response, it is necessary to develop accurate dynamic models of DGs in microgrids. This paper proposes a systematic method to calibrate synchronous-generator-interfaced DG models. Challenges associated with model validation in microgrids are analyzed. Underlying parameters are categorized into two groups, i.e., the steady-state parameters and time constants, which are estimated in two successive stages using multiple event data. This methodology ensures that the validated model with the estimated parameters are applicable for all events. A modified unscented Kalman filter (UKF) is proposed to deal with the cases in which measurement of power angle is unavailable. The effectiveness of the proposed method is validated using field test data of a 2.6 MVA diesel generator in a real microgrid.

Wang, Zhiwen↗
DOE OSTI2024

Quantification of Hydrogen Isotopes Utilizing Raman Spectroscopy Paired with Chemometric Analysis for Application across Multiple Systems

On-line and real-time analysis of a chemical process is a major analytical challenge that can drastically change the way the chemical industry or chemical research operates. With in situ analyses, new and powerful understanding of chemistry can be gained; however, building robust tools for long-term monitoring faces many challenges that include compensating for instrument drift, instrument replacement, and sensor or probe replacement. Accounting for these changes by recollecting calibration data and rebuilding quantification models can be costly and time consuming. Here, in this study, methods to overcome these challenges are demonstrated with an application of Raman spectroscopy to monitoring hydrogen isotopes with varied speciation within dynamic gas streams. Specifically, chemical data science tools such as chemometric modeling are leveraged along with several examples of calibration transfer approaches. Furthermore, the optimization of instrument and sensors cell parameters for targeted gas phase analyses is discussed. While the particular focus on hydrogen is highly beneficial within the nuclear energy sector, mechanisms built and demonstrated here are widely applicable to optical spectroscopy monitoring in numerous other chemical systems that can be leveraged in other hazardous processes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

How Cooperatively Folding Are Homopolymer Molecular Knots?

Detailed thermodynamic analysis of complex systems with multiple stable configurational states allows for insight into the cooperativity of each individual transition. In this work, we derive a heat capacity decomposition comprising contributions from each individual configurational state, which together sum to a baseline heat capacity, and contributions from each state-to-state transition. We apply this analysis framework to a series of replica exchange molecular dynamics simulations of linear and 1-1 coarse-grained homo-oligomer models, which fold into stable, configurationally well-defined molecular knots, in order to better understand the parameters leading to stable and cooperative folding of these knots. We find that a stiff harmonic backbone bending angle potential is key to achieving knots with specific 3D structures. Tuning the backbone equilibrium angle in small increments yields a variety of knot topologies. Populations of different knotted states as a function of temperature can also be manipulated by tuning the backbone torsion stiffness or by adding side-chain beads. We find that sharp total heat capacity peaks for the homo-oligomer knots are largely due to a coil-to-globule transition rather than a cooperative knotting step. Furthermore, in some cases, the cooperativity of globule-to-knot and coil-to-globule transitions is comparable, suggesting that highly cooperative folding to knotted structures can be achieved by refining the model parameters or adding sequence specificity.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Structural distortion and electron redistribution in dual-emitting gold nanoclusters

Deciphering the complicated excited-state process is critical for the development of luminescent materials with controllable emissions in different applications. Here we report the emergence of a photo-induced structural distortion accompanied by an electron redistribution in a series of gold nanoclusters. Such unexpected slow process of excited-state transformation results in near-infrared dual emission with extended photoluminescent lifetime. We demonstrate that this dual emission exhibits highly sensitive and ratiometric response to solvent polarity, viscosity, temperature and pressure. Thus, a versatile luminescent nano-sensor for multiple environmental parameters is developed based on this strategy. Furthermore, we fully unravel the atomic-scale structural origin of this unexpected excited-state transformation, and demonstrate control over the transition dynamics by tailoring the bi-tetrahedral core structures of gold nanoclusters. Overall, this work provides a substantial advance in the excited-state physical chemistry of luminescent nanoclusters and a general strategy for the rational design of next-generation nano-probes, sensors and switches.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2025

Computational alchemy clarifies origins of alloy strengthening

Solid solution strengthening (SSS) is widely used to enhance mechanical properties of metals. Originally developed for dilute alloys, classical SSS theories are presently challenged by the rise of complex concentrated alloys (CCA) with nearly equiatomic compositions. Here, we propose and develop a method of “computational alchemy” in which interatomic interactions are modified to systematically vary two key physical parameters defining SSS - atomic size misfit and elastic stiffness misfit - over a maximally wide range of two misfits. The resulting alchemical alloys are subjected to massive (~10 8 atoms) molecular dynamics (MD) simulations reproducing full complexity of plastic strength response. At variance with prevailing views, stiffness misfit is observed to contribute to SSS on par if not more than size misfit. Furthermore, depending on exactly how two misfits are combined, they result in synergistic (amplification) or antagonistic (compensation) effect on alloy strengthening. Unlike real CCAs in which each component element comes with its own specific size and stiffness, our alchemical model alloys span the space of two misfits continuously revealing trends in alloy strengthening unrecognized so far. Our study demonstrates unique value of intentionally unrealistic models for gaining deep physical insights into material behaviors that are difficult to reveal otherwise.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Liquid state theory of the structure of model polymerized ionic liquids

We employ polymer integral equation theory to study a simplified model of semiflexible polymerized ionic liquids (PolyILs) that interact via hard core repulsions and short range screened Coulomb interactions. The multi-scale structure in real and Fourier space of PolyILs (ions chosen to mimic Li, Na, K, Br, PF 6 , and TFSI) are determined as a function of melt density, Coulomb interaction strength, and ion size. Comparisons with a homopolymer melt, a neutral polymer–solvent-like athermal mixture, and an atomic ionic liquid are carried out to elucidate the distinct manner that ions mediate changes of polymer packing, the role of excluded volume effects, and the influence of chain connectivity, respectively. The effect of Coulomb strength depends in a rich manner on ion size and density, reflecting the interplay of steric packing, ion adsorption, and charge layering. Ion-mediated bridging of monomers is found, which intensifies for larger ions. Intermediate range charge layering correlations are characterized by a many-body screening length that grows with PolyIL density, cooling, and Coulomb strength, in disagreement with Debye–Hückel theory, but in accord with experiments. Qualitative differences in the collective structure, including an ion-size-dependent bifurcation of the polymer structure factor peak and pair correlation function, are predicted. The monomer cage order parameter increases significantly, but its collective ion counterpart decreases, as ions become smaller. Such behaviors allow one to categorize PolyILs into two broad classes of small and large ions. Furthermore, dynamical implications of the predicted structural results are qualitatively discussed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

The influence of laser energy on deuterium emission characteristics from a Zircaloy-4 plasma

Laser-produced plasma coupled with optical emission spectroscopy (OES) is a promising technique for detecting certain isotopes, with unique capabilities such as standoff and rapid detection and minimal to no sample preparation requirements. The key figure-of-merit for isotopic analysis using optical spectroscopy tools is the linewidth relative to the isotope shift. Although the isotopes of hydrogen (1H, 2H, and 3H) possess large isotopic shifts (1H–2H ≈ 180 pm, 1H–3H ≈ 240 pm), being a light element, the H transitions are susceptible to various broadening mechanisms in the plasma environment. One of the critical parameters that influence the linewidth of a transition in an LPP is the incident laser energy. In the present study, we evaluated the role of laser energy on plume expansion dynamics, deuterium emission intensity, and linewidth in a nanosecond laser-produced Zircaloy-4 plasma. The changes in 2Hα emission intensity and linewidth were investigated for varying laser fluence and time after plasma onset. Spatially resolved and spatially integrated OES were performed and compared to investigate the emission spectral features and linewidth of 2Hα. Monochromatic two-dimensional time-resolved imaging was also performed to understand the morphology of the deuterium and protium emission relative to all species in the plume. Our results showed that 1Hα and 2Hα emissions predominantly occur closer to the target. Measurements of 2Hα linewidth approached similar values at later times of plasma evolution regardless of the laser energy. The linewidths of the 2Hα transition showed insignificant differences between spatially resolved and spatially integrated measurements.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Deriving effective electrode–ion interactions from free-energy profiles at electrochemical interfaces

Understanding ion adsorption at electrified metal–electrolyte interfaces is essential for accurate modeling of electrochemical systems. Here, in this study, we systematically investigate the free energy profiles of Na + , Cl − , and F − ions at the Au(111)–water interface using enhanced sampling molecular dynamics with both classical force fields and machine-learned interatomic potentials (MLIPs). Our classical metadynamics results reveal a strong dependence of predicted ion adsorption on the Lennard-Jones parameters, highlighting that—without due care—standard mixing rules can lead to qualitatively incorrect descriptions of ion–metal interactions. We present a systematic methodology for tuning the cross term LJ parameters to control adsorption energetics in agreement with more accurate models. As a surrogate for an ab initio model, we employed the recently released Universal Models for Atoms MLIP, which validates classical trends and displays strong specific adsorption for chloride, weak adsorption for fluoride, and no specific adsorption for sodium, in agreement with experimental and theoretical expectations. By integrating molecular-level adsorption free energies into continuum models of the electric double layer, we show that specific ion adsorption substantially alters the interfacial ion population, the potential of zero charge, and the differential capacitance of the system. Our results underscore the critical importance of force field parameterization and advanced interatomic potentials for the predictive modeling of ion-specific effects at electrified interfaces and provide a robust framework for bridging molecular simulations and continuum electrochemical models.

Roncoroni, Fabrice [Lawrence Berkeley National Lab↗
DOE OSTI2023

Chaotic saddles and interior crises in a dissipative nontwist system

Here, we consider a dissipative version of the standard nontwist map. Nontwist systems present a robust transport barrier, called the shearless curve, that becomes the shearless attractor when dissipation is introduced. This attractor can be regular or chaotic depending on the control parameters. Chaotic attractors can undergo sudden and qualitative changes as a parameter is varied. These changes are called crises, and at an interior crisis the attractor suddenly expands. Chaotic saddles are nonattracting chaotic sets that play a fundamental role in the dynamics of nonlinear systems; they are responsible for chaotic transients, fractal basin boundaries, and chaotic scattering, and they mediate interior crises. In this work we discuss the creation of chaotic saddles in a dissipative nontwist system and the interior crises they generate. We show how the presence of two saddles increases the transient times and we analyze the phenomenon of crisis induced intermittency.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Carbon-Based Multi-Phase Rotating Detonation Engine

Abstract This article focuses on extending an H2/air rotating detonation engine's detonability limits by introducing solid carbon particles into the combustor. Carbon black particles consisting of 1% volatility and a carbon concentration of 99% were used as a solid-phase mixing agent for enhanced reaction wave dynamics. Carbon black was found to sustain detonations over multiple operational regimes formerly unattainable without carbon particles. The experiments confirm that detonations were attainable over a wide range of operational parameters, including the total mass flux flowing through the annulus (≅120–270 kg/s m2), the hydrogen/air equivalence ratios (0.65–1.0), and carbon additions (0–20 g). Chemiluminescence imaging was used to visualize the detonation wave within the annulus, quantify detonation wave velocities, and define a detonability map. The detonability map demonstrates the advantage of carbon addition, shows that detonation-based combustion can be sustained at leaner equivalence ratios, reduces hydrogen consumption dependency. The detonation wave velocities decreased as the H2/air equivalence ratio was reduced, where, in general, the detonation wave velocities decreased with respect to the Chapman–Jouguet velocity, suggesting a decrease in the detonation waves efficiency with reduced H2 concentrations. However, an extraordinary phenomenon was witnessed at very lean H2/air equivalence ratios and low mass flux conditions, where the detonation wave velocity increased upward of 100 m/s. This variation is a direct effect of the carbon particles, which drive the detonation wave. Thus, the results demonstrate that carbon particles’ addition provides an economically feasible solution to sustain high-efficiency energy production.

Energy & Fuels↗
DOE OSTI2022

Assessing the Effectiveness of Generalized Likelihood Ratio Test Detector Schemes in Seismic Event Detection and the Avoidance of Nontarget Signals

Cross-correlation techniques have played a long-standing and pivotal role in seismic event monitoring. However, the performance of correlation-based detectors is challenged by nuisance seismicity, or nontarget signals. Such detections are a problem when the mission is to automatically map events to the correct source region. Using aftershocks of the 2014 $M_w$ 6.0 South Napa, California, earthquake, we demonstrate the effectiveness of utilizing a dynamic correlation processor framework in a generalized likelihood ratio test (GLRT) detector configuration to minimize nontarget detections. A GLRT maximizes a detection statistic with respect to one or more unknown parameters. In this case, the detection statistic is a template signal match against the waveform in a window sliding over a data stream, and the unknown parameter is an index variable indicating group membership of the template event or events. Detected events are assigned to the event group that yields the largest detection statistic. In this work, our results show that a GLRT detector will outperform a suite of independently operating correlation and subspace detectors in terms of having a lower nontarget detection rate at a given missed detection rate. We also show that a GLRT detector composed of a few high-rank subspace detectors has a slightly higher nontarget detection rate, but a significantly lower missed detection rate, than a GLRT detector composed of many low-rank subspace detectors. The high-rank GLRT configuration produced impressive results even with marginal data (single channel, single station, and very low time bandwidth product), which bodes well for the utility of building efficient aftershock classification systems and global monitoring systems at larger scales. However, future work is required to assess performance at the regional scale and to assess the performance of the system at detecting target events not used in the detector template creation.

58 GEOSCIENCES↗
DOE OSTI2021

Open Source Fault-tolerant Grid Frequency Measurement for Solar Inverters

The Discrete Fourier transform (DFT) based measurement algorithms are one of the most common measurement algorithms for grid parameter estimation such as rms, phase angle, frequency. Over the past few years, many DFT based algorithms have been developed to enhance its measurement accuracy under steady-state and/or dynamic grid conditions. For example, an adaptive band-pass filter utilizing exponential modulation filter has been proposed to reduce measurement errors at the presence of large frequency deviations. Measurement accuracy of different algorithms including FIR filter, extended Kalman filtering (EKF), and enhanced DFT method have been compared in detail under different grid conditions. Two artificial signals that have 90-degree phase difference were constructed by the Clarke transformation to address the frequency spectrum leakage of DFT. A multi-module approach was developed to enhance both steady-state and dynamic measurement accuracies, in which each module was developed to eliminate some specific errors. Besides DFT-based measurement algorithms, some signal model-based algorithms have been developed to further improve the accuracy under dynamic conditions. However, a key drawback of the state-of-the-art algorithms is that they cannot perform measurements accurately during system transient faults. In the Blue Cut Fire event, there was a phase angle jump of about 26 degrees in the voltage waveform during the transient fault. The phase angle jump fault will cause waveform discontinuity, and these algorithms will fail to provide reliable measurements during this period because they typically assume the waveform to be measured is continuous, no matter what method (DFT, PLL, EKF, FIR, or Taylor WLS) is used for estimation. In fact, the measurement errors during the system transient faults like phase-jump is not required in the IEEE Standard. As a result, although a measurement instrument can pass the strict IEEE Standard, it could still be the source of the problem in the future if we have similar system transient faults, which could happen again. Therefore, developing the fault-tolerant measurement technology is the key to solve the problem.

14 SOLAR ENERGY↗
DOE OSTI2021

Recurring Planetary Debris Transits and Circumstellar Gas around White Dwarf ZTF J0328–1219

We present follow-up photometry and spectroscopy of ZTF J0328-1219, strengthening its status as a white dwarf exhibiting transiting planetary debris. Using TESS and Zwicky Transient Facility photometry, along with follow-up high-speed photometry from various observatories, we find evidence for two significant periods of variability at 9.937 and 11.2 hr. We interpret these as most likely the orbital periods of different debris clumps. Changes in the detailed dip structures within the light curves are observed on nightly, weekly, and monthly timescales, reminiscent of the dynamic behavior observed in the first white dwarf discovered to harbor a disintegrating asteroid, WD 1145+017. We fit previously published spectroscopy along with broadband photometry to obtain new atmospheric parameters for the white dwarf, with M * = 0.731 ± 0.023 M⊙, T eff = 7630 ± 140 K, and [Ca/He] = - 9.55 ± 0.12. With new high-resolution spectroscopy, we detect prominent and narrow Na D absorption features likely of circumstellar origin, with velocities 21.4 ± 1.0 km s -1 blueshifted relative to atmospheric lines. We attribute the periodically modulated photometric signal to dusty effluents from small orbiting bodies such as asteroids or comets, but we are unable to identify the most likely material that is being sublimated, or otherwise ejected, as the environmental temperatures range from roughly 400 to 700 K.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Turbulent Regimes in Collisions of 3D Alfvén-wave Packets

Using three-dimensional gyrofluid simulations, we revisit the problem of Alfvén-wave (AW) collisions as building blocks of the Alfvénic turbulent cascade and their interplay with magnetic reconnection at magnetohydrodynamic (MHD) scales. Depending on the large-scale value of the nonlinearity parameter χ 0 (the ratio between the AW linear propagation time and nonlinear turnover time), different regimes are observed. For strong nonlinearities (χ 0 ~ 1), turbulence is consistent with a dynamically aligned, critically balanced cascade—fluctuations exhibit a scale-dependent alignment $\sin {\theta }_{{k}_{\perp }}\propto {k}_{\perp }^{-1/4}$, resulting in a ${k}_{\perp }^{-3/2}$ spectrum and ${k}_{\parallel }\propto {k}_{\perp }^{1/2}$ spectral anisotropy. At weaker nonlinearities (small χ 0 ), a spectral break marking the transition between a large-scale weak regime and a small-scale ${k}_{\perp }^{-11/5}$ tearing-mediated range emerges, implying that dynamic alignment occurs also for weak nonlinearities. At χ 0 < 1 the alignment angle ${\theta }_{{k}_{\perp }}$ shows a stronger scale dependence than in the χ 0 ~ 1 regime, namely $\sin {\theta }_{{k}_{\perp }}\propto {k}_{\perp }^{-1/2}$ at χ 0 ~ 0.5, and $\sin {\theta }_{{k}_{\perp }}\propto {k}_{\perp }^{-1}$ at χ 0 ~ 0.1. Dynamic alignment in the weak regime also modifies the large-scale spectrum, scaling approximately as ${k}_{\perp }^{-3/2}$ for χ 0 ~ 0.5 and as ${k}_{\perp }^{-1}$ for χ 0 ~ 0.1. A phenomenological theory of dynamically aligned turbulence at weak nonlinearities that can explain these spectra and the transition to the tearing-mediated regime is provided; at small χ 0 , the strong scale dependence of the alignment angle combines with the increased lifetime of turbulent eddies to allow tearing to onset and mediate the cascade at scales that can be larger than those predicted for a critically balanced cascade by several orders of magnitude. Such a transition to tearing-mediated turbulence may even supplant the usual weak-to-strong transition.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Applicability study of Bayesian optimization in core neutronic design using a toy model

At the Japan Atomic Energy Agency (JAEA), an innovative design approach named ARKADIA (Advanced Reactor Knowledge- and AI-aided Design Integration Approach through the whole plant life cycle) for advanced nuclear reactors is currently under development. One task in ARKADIA is to build a system that automatically optimizes core and fuel designs by conducting core neutronic and thermal-hydraulic calculations, fuel integrity evaluations, and plant dynamic analyses. This system will be implemented to automatically find an optimal design that minimizes (or maximizes) objective function defined by core performance while varying the core and fuel design parameters such as fuel pin diameter, core height and diameter. In this study, as the first step of system development, we focused only on core neutronic design and conducted a study of automatic optimization. As the optimization algorithm, Bayesian optimization (BO), an effective method for optimization problems with expensive computational cost of objective function, was utilized. The applicability of BO was studied based on single- and two-objective optimization examples of core neutronic design in a toy model. As a result, in the former, it was shown that BO can give the optimal solution, which matches the reference solution calculated by a brute force calculation well, with a small number of required calculations. Usability on core neutronic designs, where the computational cost per case is high, was confirmed. In the latter, it was found that BO can give a Pareto solutions-set that shows good agreement with the reference solution. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗