Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “core model proposal”

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 163 records · Page 9

BATSE Observations of the Large-Scale Isotropy of Gamma-Ray Bursts

We use dipole and quadrupole statistics to test the large-scale isotropy of the first 1005 gamma-ray bursts observed by the Burst and Transient Source Experiment (BATSE). In addition to the entire sample of 1005 gamma-ray bursts, many subsets are examined. We use a variety of dipole and quadrupole statistics to search for Galactic and other predicted anisotropies and for anisotropies in a coordinate-system independent manner. We find the gamma-ray burst locations to be consistent with isotropy, e.g., for the total sample the observed Galactic dipole moment (cos theta) differs from the value predicted for isotropy by 0.9 sigma and the observed Galactic quadrupole moment (sin(exp 2) b - 1/3) by 0.3 sigma. We estimate for various models the anisotropies that could have been detected. If one-half of the locations were within 86 deg of the Galactic center, or within 28 deg of the Galactic plane, the ensuing dipole or quadrupole moment would have typically been detected at the 99% confidence level. We compare the observations with the dipole and quadrupole moments of various Galactic models. Several Galactic gamma-ray bursts models have moments within 2 sigma of the observations; most of the Galactic models proposed to date are no longer in acceptable agreement with the data. Although a spherical dark matter halo distribution could be consistent with the data, the required core radius is larger than the core radius of the dark matter halo used to explain the Galaxy's rotation curve. Gamma-ray bursts are much more isotropic than any observed Galactic population, strongly favoring but not requiring an origin at cosmological distances.

Briggs, Michael S.↗

Nuclear data uncertainty propagation and modeling uncertainty impact evaluation in neutronics core simulation

Uncertainty analysis is a critical requirement in reactor simulation as it is used to quantify the reliability of best-estimate calculation. A comprehensive uncertainty analysis should characterize all sources of uncertainties in a computationally-feasible and scientifically-defendable manner. Here we employ a well-established reduced order modeling (ROM) based uncertainty quantification methodology to propagate uncertainties throughout neutronic calculations. ROM relies on recent advances in randomized data mining techniques applied to large data streams. In our proposed implementation, the nuclear data uncertainties are first propagated from multi-group level through lattice physics calculation to generate few-group parameter uncertainties, described using a vector of mean values and a covariance matrix. Employing an ROM-based compression of the covariance matrix, the few-group uncertainties are then propagated through downstream core simulation in a computationally efficient manner. This straightforward approach, albeit efficient as compared to brute force forward and/or adjoint-based methods, often employs a number of assumptions that have been unquestioned in the literature of neutronic uncertainty analysis. This manuscript argues that these assumptions could introduce another source of uncertainty referred to as modeling uncertainties, whose magnitude needs to be quantified in tandem with nuclear data uncertainties. Thus, our primary goal is to explore the interactions between these two uncertainty sources in order to assess whether modeling uncertainties have an impact on parameter uncertainties. To explore this endeavor, the impact of a number of modeling assumptions on core attributes uncertainties is quantified. The study employs a CANDU reactor model, with Serpent and NEWT as lattice physics solvers and NESTLE-C as core simulator. The modeling assumptions investigated include those related with the uncertainty propagation method employed, e.g., deterministic vs. stochastic, the few-group energy structure employed to represent the cross-sections, the resonance treatment in lattice physics calculation, the reference values for the cross-section, and the number of samples employed to render ROM compression. Results indicate that some of the modeling assumptions could have a non-negligible impact on the core responses propagated uncertainties, highlighting the need for a more comprehensive approach to combine parameter and modeling uncertainties.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

UO 2 -liquid metal suspension fuel concept for enhanced passive safety of LWRs: A heat pipe case study

This work proposes a novel fuel concept for use in light-water type reactors (LWRs): UO 2 particles suspended in a low-melting-point liquid metal (LM). This new fuel form offers lower energy density and higher thermal conductivity, thus enhancing the fuel’s thermal performance compared to sintered UO 2 pellets. Two candidate alloys for the LM are proposed based on previous studies into accident-tolerant fuels: Bi-Pb-Sn and Pb-Sn. Ten research questions are provided to guide future development of this fuel. A simple heat pipe model was constructed by coupling a reactor physics code to a heat transport code and other physics modules. Scoping calculations were performed on this model core as a preliminary investigation into five of the research questions. A fuel with 13 wt% UO 2 was shown to maintain a steady-state peak coolant channel wall temperature below 400 °C, and maintain criticality via self-regulating reactivity feedback. Shortcomings of the present model are discussed in this work, along with proposed future model improvements and fuel development. The appendix provides an example of a reactor concept that could be developed to utilize this fuel, but more modeling and experimental work are needed before a full reactor concept can be developed.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Proposal and application of ROM-Lasso method for sensitivity coefficient evaluation

We propose a novel method for evaluating sensitivity coefficients of neutronics parameters to cross sections, so-called the reduced-order modeling technique ROM-Lasso. In this method, cross sections of interest are randomly sampled, and corresponding perturbed core analyses are performed. Then, the sensitivity coefficient vector of the higher-level model is expanded via the active subspace bases obtained with the lower-level model whose dimensional complexity is smaller than that of the higher-level model, and the expansion coefficients are estimated by the Lasso regression. A unique feature of the ROM-Lasso method allows the use of different bases optimized for each neutronics parameter. We conducted a verification calculation for an accelerator-driven system and demonstrated that the ROM-Lasso method can reproduce the sensitivity coefficients with a much smaller number of forward calculations than the direct method. The proposed method can be used to practically evaluate sensitivity coefficients. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Tungsten and molybdenum isotopic evidence for an impact origin of pallasites

The origin of pallasites—stony-iron meteorites mainly composed of olivine and Fe-Ni metal—is debated and proposed formation scenarios broadly range from models that explain pallasite formation by internal processes in the mantle of a differentiated planetesimal to those that involve impact–induced mixing of core and mantle materials. Here, the origin of pallasites is examined by studying the nebular source regions of their precursor material using Mo isotopes and their history of metal-silicate segregation using Hf-W chronometry. We report new Mo and W isotopic data for a large suite of pallasite metal samples, alongside Pt isotope data to quantify superimposed cosmic ray exposure effects. Most main-group pallasites exhibit uniform pre-exposure 182 W and Mo isotopic compositions that bear an excellent similarity to those of IIIAB iron meteorites. Four main-group pallasites and the IIIAB iron Thunda have more radiogenic pre-exposure 182 W compositions, but display the same Mo isotopic composition as other main-group pallasites and IIIAB irons. This strong chronological and genetic link strongly suggests that main-group pallasite metal originated in the IIIAB parent body core. This, combined with prior Pd-Ag chronometric evidence for an early collisional disruption of the IIIAB parent body, implies that main-group pallasites formed by impact–induced mixing of metal and silicates rather than by an internal process on the IIIAB parent body. This mixing led to elevated 182 W compositions in some pallasites, which are best accounted for by partial re-equilibration of IIIAB metal with radiogenic 182 W from the colliding body. Altogether, our results support models that explain main-group pallasite formation by injection of pallasite metal into the mantle of another differentiated body, implying that pallasite silicates did not primarily derive from the IIIAB mantle, but instead from that of the colliding body.

58 GEOSCIENCES↗

Dynamical instability of collapsed dark matter halos

A self-interacting dark matter halo can experience gravothermal collapse, resulting in a central core with an ultrahigh density. It can further contract and collapse into a black hole, a mechanism proposed to explain the origin of supermassive black holes. We study dynamical instability of the core in general relativity. Here we use a truncated Maxwell-Boltzmann distribution to model the dark matter distribution and solve the Tolman-Oppenheimer-Volkoff equation. For given model parameters, we obtain a series of equilibrium configurations and examine their dynamical instability based on considerations of total energy, binding energy, fractional binding energy, and adiabatic index. Our numerical results indicate that the core can collapse into a black hole when the fractional binding energy reaches 0.035 with a central gravitational redshift of 0.5. We further show for the instability to occur in the classical regime, the boundary temperature of the core should be at least 10% of the mass of dark matter particles; for a 10 9 M ⊙ seed black hole, the particle mass needs to be larger than a few keV. These results can be used to constrain different collapse models, in particular, those with dissipative dark matter interactions.

79 ASTRONOMY AND ASTROPHYSICS↗

Design, Manufacture, and Testing of an Open-Source Benchmark Composite Hydrokinetic Turbine Blade: Preprint

In a trend toward clean energy alternatives, recent years have seen great strides in the marine energy space. Consequently, there is a pressing need for the design, development, and validation of novel energy harvesting technologies such as hydrokinetic devices, which capture kinetic energy from waves, tides, and currents. However, these devices span numerous concepts and designs that often lack solid benchmark research that can be freely referenced throughout their development. This work focuses on the design process of an open-source composite hydrokinetic turbine blade for a three-bladed marine turbine rotor assembly with a diameter of 2.5 m. The proposed blade consists of two structural composite skins that are bonded with an adhesive and filled with a foam core. This study also explores and contrasts the efficiency and resolution of low-fidelity rapid design methodologies and comprehensive high-fidelity approaches in the context of blade design, modeling, and analysis efforts, a key objective in this research. Blade hydrodynamic loads were modeled and applied to finite-element blade models to study deformations and potential failure. Ongoing and upcoming efforts will result in blade manufacture and structural testing at the National Renewable Energy Laboratory. In future work, multiple blades will be deployed at the Living Bridge site at the University of New Hampshire and will be compared to rigid aluminum blades of the same geometry, developed by Sandia National Laboratories. Ultimately, this research will lay foundational groundwork for researchers and manufacturers, establishing a baseline composite blade design that will serve as a benchmark in the development of future hydrokinetic turbine blades.

blade design↗

COMSOL Results for the Nominal Steady-State Operation of the Proposed 95-MW LEU Silicide Core for HFIR Conversion

Engineering design studies are being performed to determine the feasibility of converting the High Flux Isotope Reactor (HFIR) from highly enriched uranium (HEU) to low-enriched uranium (LEU) fuel at Oak Ridge National Laboratory. This activity is sponsored by the Office of Reactor Conversion and Uranium Supply (ORCUS) under the auspices of the US Department of Energy National Nuclear Security Administration’s Office of Material Management and Minimization. HFIR is a very high flux, pressurized, light water–cooled and moderated, flux trap–type research reactor with a core made of involute shaped U 3 O 8 /Al cermet fuel plates and coolant channels. HFIR currently operates at a thermal power of 85 MW and supports key national and international missions in neutron scattering, isotope production, materials/fuels irradiation, neutron activation analysis, gamma irradiation, and neutrino research. Advanced multiphysics computational fluid dynamics models have been developed in the COMSOL Multiphysics software to simulate the steady-state operating conditions for the proposed low-and high-density LEU U 3 Si 2 -Al (uranium silicide dispersion) fuel designs. The COMSOL models for HFIR inner and outer fuel element models incorporate various essential inputs and physics such as spatially dependent nuclear heat deposition, multilayer heat conduction, conjugate heat transfer, turbulent flows (using Reynolds-averaged Navier Stokes turbulence models), structural mechanics (thermal–structural interactions and fuel swelling), and oxide layer build-up. This report presents the best-estimate thermal hydraulics results for the low- and high-density optimized silicide LEU core designs at 95 MW steady-state nominal operation.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

A Parallel Saturation Algorithm on Shared Memory Architectures

Symbolic state-space generators are notoriously hard to parallelize. However, the Saturation algorithm implemented in the SMART verification tool differs from other sequential symbolic state-space generators in that it exploits the locality of ring events in asynchronous system models. This paper explores whether event locality can be utilized to efficiently parallelize Saturation on shared-memory architectures. Conceptually, we propose to parallelize the ring of events within a decision diagram node, which is technically realized via a thread pool. We discuss the challenges involved in our parallel design and conduct experimental studies on its prototypical implementation. On a dual-processor dual core PC, our studies show speed-ups for several example models, e.g., of up to 50% for a Kanban model, when compared to running our algorithm only on a single core.

Ezekiel, Jonathan↗

An improved lunar moment of inertia determination - A proposed strategy

A strategy for determining an improved lunar moment of inertia is proposed. An improved uncertainty in the lunar inhomogeneity parameter could reduce the core density error from 4.20 gr/cu cm to 0.1 gm/cu cm for the case of a lunar density model having a 300 km core radius. The current error of 0.0025 for the lunar inhomogeneity parameter is mostly due to the uncertainties in the C20 and C22 gravity harmonics. An optimum reduction strategy for obtaining an order of magnitude improvement in the gravity estimates is based on covariance analyses of Doppler data. The long-arc reduction method is shown to provide better results than the short-arc technique, and the use of an orbit with a 7000 km semimajor axis along with certain terms of up to degree and order eight is recommended.

Ananda, M. P.↗

Supersonic shock wave/vortex interaction

Although shock wave/vortex interaction is a basic and important fluid dynamics problem, very little research has been conducted on this topic. Therefore, a detailed experimental study of the interaction between a supersonic streamwise turbulent vortex and a shock wave was carried out at the Penn State Gas Dynamics Laboratory. A vortex is produced by replaceable swirl vanes located upstream of the throat of various converging-diverging nozzles. The supersonic vortex is then injected into either a coflowing supersonic stream or ambient air. The structure of the isolated vortex is investigated in a supersonic wind tunnel using miniature, fast-response, five-hole and total temperature probes and in a free jet using laser Doppler velocimetry. The cases tested have unit Reynolds numbers in excess of 25 million per meter, axial Mach numbers ranging from 2.5 to 4.0, and peak tangential Mach numbers from 0 (i.e., a pure jet) to about 0.7. The results show that the typical supersonic wake-like vortex consists of a non-isentropic, rotational core, where the reduced circulation distribution is self similar, and an outer isentropic, irrotational region. The vortex core is also a region of significant turbulent fluctuations. Radial profiles of turbulent kinetic energy and axial-tangential Reynolds stress are presented. The interactions between the vortex and both oblique and normal shock waves are investigated using nonintrusive optical diagnostics (i.e. schlieren, planar laser scattering, and laser Doppler velocimetry). Of the various types, two Mach 2.5 overexpanded-nozzle Mach disc interactions are examined in detail. Below a certain vortex strength, a 'weak' interaction exists in which the normal shock is perturbed locally into an unsteady 'bubble' shock near the vortex axis, but vortex breakdown (i.e., a stagnation point) does not occur. For stronger vortices, a random unsteady 'strong' interaction results that causes vortex breakdown. The vortex core reforms downstream of the rear stagnation point, and the reduced circulation distribution once again becomes self-similar in this region. A-new model of this interaction is proposed. Finally, a curve defining the approximate limits of supersonic vortex breakdown is presented.

Settles, G. S.↗

Jet Noise Physics and Modeling Using First-principles Simulations

An extensive analysis of our jet DNS database has provided for the first time the complex correlations that are the core of many statistical jet noise models, including MGBK. We have also for the first time explicitly computed the noise from different components of a commonly used noise source as proposed in many modeling approaches. Key findings are: (1) While two-point (space and time) velocity statistics are well-fitted by decaying exponentials, even for our low-Reynolds-number jet, spatially integrated fourth-order space/retarded-time correlations, which constitute the noise "source" in MGBK, are instead well-fitted by Gaussians. The width of these Gaussians depends (by a factor of 2) on which components are considered. This is counter to current modeling practice, (2) A standard decomposition of the Lighthill source is shown by direct evaluation to be somewhat artificial since the noise from these nominally separate components is in fact highly correlated. We anticipate that the same will be the case for the Lilley source, and (3) The far-field sound is computed in a way that explicitly includes all quadrupole cancellations, yet evaluating the Lighthill integral for only a small part of the jet yields a far-field noise far louder than that from the whole jet due to missing nonquadrupole cancellations. Details of this study are discussed in a draft of a paper included as appendix A.

Freund, Jonathan B.↗

White Paper: Research & Development for the Time at Temperature Approach

Recent advancements in nuclear power research are greatly improving reactor safety and performance through the development of Accident Tolerant Fuel (ATF) and Low-Enriched Uranium Plus (LEU+). These innovations can address Departure from Nucleate Boiling (DNB) margins, which are vital for reactor safety. DNB happens when the coolant switches to film boiling, significantly decreasing heat transfer and posing a risk of fuel cladding failure. The U.S. Nuclear Regulatory Commission (NRC) employs conservative DNB criteria, which can potentially restrict the operational flexibility and efficiency of reactors. The Time at Temperature (TaT) approach could provide a more detailed and adaptable operational guideline by establishing acceptable time-temperature limits, accounting for the duration a material can withstand elevated temperatures without losing its integrity. This method allows reactors to operate more efficiently and safely, offering additional operational margins, faster power adjustments, and improved fuel cycle economics. TaT criteria allow for higher power levels and more flexible responses to operational transients, particularly applicable for anticipated operational occurrences (AOOs) that result in short durations of post-DNB conditions. It enhances plant operational flexibility, allows faster startup times, and enables quicker power level adjustments, optimizing fuel loading patterns and improving fuel cycle economics. Implementing TaT limits reduces core design constraints, lowers fuel usage, and reduces costs, essential for the long-term sustainability of Light Water Reactors (LWRs). TaT maximizes the use of advanced fuel technologies like ATF and LEU+, further enhancing their economic and environmental benefits. To apply the TaT approach in existing LWRs, collaborative research activities among various DOE-sponsored programs are essential. These efforts should incorporate fuel experiments, physics-based high-fidelity modeling, ML-based surrogate modeling, and optimization techniques. This whitepaper proposes four research and development areas: 1) Investigation of the feasibility of new operations of LWR with updated safety limits; 2) Assessment of reactor operation limits through uncertainty reduction; 3) Evaluation of power uprate in virtual environment; and 4) Lattice and reactor core design for power uprate. Each area includes why this research is in need and a suggested scope of work. These comprehensive research areas ensure practical and beneficial advancements for existing reactors, translating innovations in nuclear fuel and cladding technology into improved reactor performance and safety.

42 - ENGINEERING↗

Planning Amidst Uncertainty: Identifying Core CCS Infrastructure Robust to Storage Uncertainty

Carbon Capture and Storage (CCS) is a critical technology for reducing anthropogenic CO2 emissions, but its large-scale deployment is complicated by uncertainties in geological storage performance. These uncertainties pose significant financial and operational risks, as underperforming storage sites can lead to costly infrastructure modifications, inefficient pipeline routing, and economic shortfalls. To address this challenge, we propose a novel optimization workflow that is based on mixed-integer linear programming and explicitly integrates probabilistic modeling of storage uncertainty into CCS infrastructure design. This workflow generates multiple infrastructure scenarios by sampling storage capacity distributions, optimally solving each scenario using a mixed-integer linear programming model, and aggregating results into a heatmap to identify core infrastructure components that have a low likelihood of underperforming. A risk index parameter is introduced to balance trade-offs between cost, CO2 processing capacity, and risk of underperformance, allowing stakeholders to quantify and mitigate uncertainty in CCS planning. Applying this workflow to a CCS dataset from the US Department of Energy’s Carbon Utilization and Storage Partnership project reveals key insights into infrastructure resilience. Reducing the risk index from 15% to 0% is observed to lead to an 83.7% reduction in CO2 processing capacity and a 77.1% decrease in project profit, quantifying the trade-off between risk tolerance and project performance. Furthermore, our results highlight critical breakpoints, where small adjustments in the risk index produce disproportionate shifts in infrastructure performance, providing actionable guidance for decision-makers. Unlike prior approaches that aimed to cheaply repair underperforming infrastructure, our workflow constructs robust CCS networks from the ground up, ensuring cost-effective infrastructure under storage uncertainty. These findings demonstrate the practical relevance of incorporating uncertainty-aware optimization into CCS planning, equipping decision-makers with a tool to make informed project planning decisions.

Olson, Daniel↗

Comparing simulated demand flexibility against actual performance in commercial office buildings

Commercial building energy benchmarking has been used as a mechanism to evaluate energy use of a single building over time, relative to other similar buildings, or to simulations of a reference building conforming to various energy standards. Lack of empirical demand flexibility data and consistent flexibility metrics has limited the ability to compare demand flexibility performance with estimated demand flexibility in buildings. In this study, we collected demand response performance data for a total of 831 demand response events from 192 sites as a first step to build such a building demand flexibility dataset, and propose a standard core data schema to consolidate field data from different sources. We also performed parametric simulations of a control strategy called “global temperature adjustment” using commercial office prototype building models. We then compared the simulated demand flexibility performance against the actual data for offices with global temperature adjustment strategy implemented. During demand response events with an average outside air temperature of 34 °C (range 23 °C–42 °C), the measured demand decrease intensity of the demand flexibility metrics were 6.1 watts per square meter (W/m 2 ), 10.0 W/m 2 , 11.1 W/m 2 , 7.1 W/m 2 , and 4.7 W/m 2 for small, small–medium, medium, medium–large, and large office buildings, respectively. Compared to the measured data in medium- and large-size buildings, the simulated demand decrease intensity was 0.7 W/m 2 (17%) lower on average. The discrepancy between simulated and measured peak demand intensities fell within one standard deviation of the mean measured data. Here, the comparison results validate the credibility of simulations in capturing real building data for assessing the technical potential of building demand flexibility.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Martini 3 Coarse-Grained Force Field for Cholesterol

Cholesterol plays a crucial role in biomembranes by regulating various properties, such as fluidity, rigidity, permeability, and organization of lipid bilayers. The latest version of the Martini model, Martini 3, offers significant improvements in interaction balance, molecular packing, and inclusion of new bead types and sizes. However, the release of the new model resulted in the need to reparameterize many core molecules, including cholesterol. Here, we describe the development and validation of a Martini 3 cholesterol model, addressing issues related to its bonded setup, shape, volume, and hydrophobicity. In conclusion, the proposed model mitigates some limitations of its Martini 2 predecessor while maintaining or improving the overall behavior.

59 BASIC BIOLOGICAL SCIENCES↗

On discriminating tokamak sawtooth crash models via localized density and temperature measurements

The core electron temperature drops rapidly during the sawtooth crash in tokamak plasmas, which causes heat loss and may lead to fast particle losses or even a disruption. Several models have been proposed for the periodic crash, including the Kadomtsev model with magnetic reconnection and the quasi-interchange model with the growth of higher-mode-number pressure-driven instabilities. 3D MHD simulations were performed for these two models with a goal to develop intuition and to predict qualitatively how different types of sawtooth will appear in various diagnostics. The structures of electron density ne and electron temperature Te show a dominant (1, 1) mode for the Kadomtsev case and a dominant (4, 4) mode for the quasi-interchange case. The oscillations of ne and Te have a positive correlation near the inversion layer for both cases, while their frequencies and amplitudes are different depending on the dominant modes. Particularly, for the Kadomtsev case, we find a relation between the amount of flux reconnected during a sawtooth event and ne or Te oscillations. Therefore, we connect recently developed measurement capabilities for ne and Te to the internal sawtooth behavior. We propose that this method of analysis can help in identifying the type of sawtooth in future experiments augmented by simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Prediction of cccDNA dynamics in hepatitis B patients by a combination of serum surrogate markers

Quantification of intrahepatic covalently closed circular DNA (cccDNA) is a key for evaluating an elimination of hepatitis B virus (HBV) in infected patients. However, quantifying cccDNA requires invasive methods such as a liver biopsy, which makes it impractical to access the dynamics of cccDNA in patients. Although HBV RNA and HBV core-related antigens (HBcrAg) have been proposed as surrogate markers for evaluating cccDNA activity, they do not necessarily estimate the amount of cccDNA. Here, we employed a recently developed multiscale mathematical model describing intra- and intercellular viral propagation and applied it in HBV-infected patients under treatment. We developed a model that can predict intracellular HBV dynamics by use of extracellular viral markers, including HBsAg, HBV DNA, and HBcrAg in peripheral blood. Importantly, the model prediction of the amount of cccDNA in patients over time was confirmed to be well correlated with the data for quantified cccDNA by paired liver biopsy. Thus, our method combining classic and emerging surrogate markers enables us to predict the decay dynamics of cccDNA in patients undergoing treatment.

60 APPLIED LIFE SCIENCES↗