Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “HOMOGENIZATION METHODS”

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,495 records · Page 83
All sourcesNASA NTRSDOE OSTI
NASA NTRS1956

On the Contribution of Turbulent Boundary Layers to the Noise Inside a Fuselage

The following report deals in preliminary fashion with the transmission through a fuselage of random noise generated on the fuselage skin by a turbulent boundary layer. The concept of attenuation is abandoned and instead the problem is formulated as a sequence of two linear couplings: the turbulent boundary layer fluctuations excite the fuselage skin in lateral vibrations and the skin vibrations induce sound inside the fuselage. The techniques used are those required to determine the response of linear systems to random forcing functions of several variables. A certain degree of idealization has been resorted to. Thus the boundary layer is assumed locally homogeneous, the fuselage skin is assumed flat, unlined and free from axial loads and the 'cabin' air is bounded only by the vibrating plate so that only outgoing waves are considered. Some of the details of the statistical description have been simplified in order to reveal the basic features of the problem. The results, strictly applicable only to the limiting case of thin boundary layers, show that the sound pressure intensity is proportional to the square of the free stream density, the square of cabin air density and inversely proportional to the first power of the damping constant and to the second power of the plate density. The dependence on free stream velocity and boundary layer thickness cannot be given in general without a detailed knowledge of the characteristics of the pressure fluctuations in the boundary layer (in particular the frequency spectrum). For a flat spectrum the noise intensity depends on the fifth power of the velocity and the first power of the boundary layer thickness. This suggests that boundary layer removal is probably not an economical means for decreasing cabin noise. In general, the analysis presented here only reduces the determination of cabin noise intensity to the measurement of the effect of any one of our variables (free stream velocity, boundary layer thickness, plate thickness or the characteristic velocity of propagation in the plate). The plate generates noise by vibrating in resonance over a wide range of frequencies and increasing the damping constant is consequently an effective method of decreasing noise generation. One of the main features of the results is that the relevant quantities upon which noise intensity depends are non-dimensional numbers in which boundary layer and plate properties enter as ratios. This is taken as an indication that in testing models of structures for boundary layer noise it is not sufficient to duplicate in the model the structural characteristics of the fuselage. One must match properly the characteristics of the exicitng pressure fluctuations to that of the structure.

Corcos, G. M.↗
NASA NTRS2016

Development of a Subcell Based Modeling Approach for Modeling the Architecturally Dependent Impact Response of Triaxially Braided Polymer Matrix Composites

Understanding the high velocity impact response of polymer matrix composites with complex architectures is critical to many aerospace applications, including engine fan blade containment systems where the structure must be able to completely contain fan blades in the event of a blade-out. Despite the benefits offered by these materials, the complex nature of textile composites presents a significant challenge for the prediction of deformation and damage under both quasi-static and impact loading conditions. The relatively large mesoscale repeating unit cell (in comparison to the size of structural components) causes the material to behave like a structure rather than a homogeneous material. Impact experiments conducted at NASA Glenn Research Center have shown the damage patterns to be a function of the underlying material architecture. Traditional computational techniques that involve modeling these materials using smeared homogeneous, orthotropic material properties at the macroscale result in simulated damage patterns that are a function of the structural geometry, but not the material architecture. In order to preserve heterogeneity at the highest length scale in a robust yet computationally efficient manner, and capture the architecturally dependent damage patterns, a previously-developed subcell modeling approach where the braided composite unit cell is approximated as a series of four adjacent laminated composites is utilized. This work discusses the implementation of the subcell methodology into the commercial transient dynamic finite element code LS-DYNA (Livermore Software Technology Corp.). Verification and validation studies are also presented, including simulation of the tensile response of straight-sided and notched quasi-static coupons composed of a T700/PR520 triaxially braided [0deg/60deg/-60deg] composite. Based on the results of the verification and validation studies, advantages and limitations of the methodology as well as plans for future work are discussed.

finite element method↗
DOE OSTI2023

Investigating the large-scale environment of wide-angle tailed radio galaxies in the local Universe

We present a statistical analysis of the large-scale (up to two Mpc) environment of a homogeneous and complete sample, both in radio and optical selection, of wide-angle tailed radio galaxies (WATs) in the local Universe (i.e., with redshifts z ≲ 0.15). The analysis is carried out using the parameters obtained from cosmological neighbors within two Mpc of the target source. Results regarding the large-scale environments of the WATs are then compared with those of Fanaroff–Riley type I (FR Is) and type II (FR IIs) radio galaxies listed in two other homogeneous and complete catalogs and selected with the same criterion adopted for the WATs catalog. We obtain indication that at low redshift, WATs inhabit environments with a larger number of galaxies than those of FR Is and FR IIs. In the explored redshift range, the physical size of the galaxy group or cluster in which the WATs reside appears to be almost constant with respect to FR Is and FR IIs, being around one Mpc. From the distribution of the concentration parameter, defined as the ratio between the number of cosmological neighbors lying within 500 kpc and within one Mpc, we conclude that WATs tend to inhabit the central region of the group or cluster in which they reside, in agreement with the general paradigm that WATs are the brightest cluster or group galaxy.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Phonon-Assisted Ballistic Current from First-Principles Calculations

The bulk photovoltaic effect (BPVE) refers to current generation due to illumination by light in a homogeneous bulk material lacking inversion symmetry. In addition to the intensively studied shift current, the ballistic current, which originates from asymmetric carrier generation due to scattering processes, also constitutes an important contribution to the overall kinetic model of the BPVE. In this Letter, we use a perturbative approach to derive a formula for the ballistic current resulting from the intrinsic electron-phonon scattering in a form amenable to first-principles calculation. We then implement the theory and calculate the ballistic current of the prototypical BPVE material BaTiO 3 using quantum-mechanical density functional theory. The magnitude of the ballistic current is comparable to that of the shift current, and the total spectrum (shift plus ballistic) agrees well with the experimentally measured photocurrents. Furthermore, we show that the ballistic current is sensitive to structural change, which could benefit future photovoltaic materials design.

14 SOLAR ENERGY↗
DOE OSTI2022

Report on Depth-Dependent Q from Frequency-Dependent Lg Q

As a phase that travels through the crustal waveguide, Lg samples the velocity and attenuation structure of the crust. We have developed methods (Pasyanos et al., 2009a; Pasyanos et al., 2009b) to estimate Qp and Qs of the crust and upper mantle from the amplitudes of regional phases, and have applied it to a number of regions around the world, and over a broad frequency band (0.5 - 10 Hz). In prior attenuation modeling, Lg phase amplitudes are able to effectively distinguish the high Q (low attenuation) of old stable cratonic crust from the low Q (high attenuation) of more recently active tectonic regions. For a variety of purposes (e.g. waveform modeling, seismic hazard, etc.), it is useful to have more information on the depth-dependence of the Q structure in the crust. One obvious candidate for doing this is using the frequency-dependence information of Lg Q, in the manner of surface waves. Unfortunately, however, Lg samples the whole of the crust without much differentiation with respect to frequency that is characteristic of fundamental mode surface waves. Lg is generally thought to be the sum of supercritically reflected S-waves trapped in the crustal waveguide (e.g. Xie and Lay, 1994). The phase can also be modeled as a superposition of surface waves. At lower frequencies, it is primarily composed of fundamental and low-order surface waves, while at higher frequencies, it is comprised purely of many high-order surface waves. The low-order surface waves that comprise Lg in low frequency bands (< 1 Hz) are preferentially sensitive to shallow earth structure, such as sedimentary basins and upper crust, while high-order surface waves have a more homogenized sensitivity to the whole crust. Furthermore, at the lower frequencies, we might expect to see some differentiation in the sensitivity between early Lg (herein referred to as Lg1) and late Lg (referred to as Lg2) where the latter is comprised of lower-order surface waves and, hence, shallower structure. Therefore, we might expect to see the signature of shallow features, primarily the lower Q of sediments as compared to crystalline crust, at lower frequencies and in later arriving Lg amplitudes.

58 GEOSCIENCES↗
DOE OSTI2020

Steam-Induced Coarsening of Single-Unit-Cell MFI Zeolite Nanosheets and Its Effect on External Surface Brønsted Acid Catalysis

Commonly used methods to assess crystallinity, micro-/mesoporosity, Brønsted acid site density and distribution (in micro- vs. mesopores), and catalytic activity suggest nearly invariant structure and function for aluminosilicate zeolite MFI two-dimensional nanosheets before and after superheated steam treatment. Yet, pronounced reaction rate decrease for benzyl alcohol alkylation with mesitylene, a reaction that cannot take place in the zeolite micropores, is observed. Transmission electron microscopy images reveal pronounced changes in nanosheet thickness, aspect ratio and roughness indicating that nanosheet coarsening and the associated changes in the external (mesoporous) surface structure are responsible for the changes in the external surface catalytic activity. Finally, superheated steam treatment of hierarchical zeolites can be used to alter nanosheet morphology and regulate external surface catalytic activity while preserving micro- and mesoporosity, and micropore reaction rates.

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

Atomic-Scale Surface Studies of Bulk Metallic Glasses. Final Report

Bulk metallic glasses (BMGs) are of both scientific and technological interest because the absence of periodic atomic arrangements provides them with unique physical, chemical, and mechanical properties. Their high strength, superior elasticity, and an ability to be easily formed into virtually unlimited shapes with feature sizes from centimeters to Angstroms by blow molding and thermoplastic forming makes them an attractive choice for more and more practical applications and products. Due to their complex internal structure, however, experiments that yield insight into their exact atomic arrangements have been scarce. As a result, glass physics is one of the last remaining unexplored fields of materials science despite the scientific and technological importance of glasses in our daily lives, and the question how to characterize and control matter away from equilibrium, as glasses are, was listed as one of five Grand Challenges in a recent DoE report. The main reason for the slow progress in glass physics is the lack of experimental tools that enable access to atomic-scale structural and behavioral information for disordered materials. Such atomic-scale knowledge is mandatory to establish structure-property relationships that ultimately could allow to custom-design alloys featuring specific desired characteristics. With no such relationships available, theory development in glass remains basic and the few that exist are often untested. The aim of this research was to enable progress in our understanding of BMGs by developing a new approach that will allow a meaningful application of local surface science methods to specially prepared BMG samples to obtain a wealth of quantitative information on their atomic arrangements. Key was the availability of specially prepared samples whose surfaces feature large atomically flat terraces despite being entirely amorphous, which we have produced from a Pt 57.5 Cu 14.7 Ni 5.3 P 22.5 alloy (‘Pt-BMG’) both under ambient conditions as well as in ultrahigh vacuum using a unique setup that has been specially developed within this grant. Our approach starts with the in-situ preparation of oxide crystals that are terminated by large terraces, from which exact mirror images out of BMG will be produced using thermoplastic forming (TPF); for the research within this grant, we have successfully used (001)-oriented SrTiO 3 single crystal surfaces as well as (100)-, (110)-, and (111)-oriented single crystals made from LaAlO 3 . Since the resulting BMG replicas display all features of the original crystal with sub-Angstrom fidelity, thereby mimicking the original crystal’s termination by atomically flat terraces without being crystalline themselves, they are ideally suited for further investigation. The following atomic-scale local studies were then carried within this proposal: (i) high-resolution surface imaging and local spectroscopy using scanning probe microscopy, which showed disordered atom-like features and revealed changes in the gradient of the local surface potential on a 1-2 nm length scale; (ii) characterization of atomic-scale plastic flow with affected volumes as low as 1000 atoms, which showed local hardness near or above the theoretically predicted maximum and, once plastic deformation was initiated, homogeneous flow of the atoms involved; (iii) characterization of surface relaxation processes and the onset of crystallization induced by annealing, which showed that upon heating over the material’s glass transition temperature, the surface rearranges and relaxes towards a more stable, denser packed glass, which increases on-terrace surface roughness, while surface tension smoothens step edges; and (iv) studies that investigate the dependence of the material’s mechanical properties and structure on processing parameters, revealing that relaxed glasses get denser, harder, and more elastic. In combination, this information allows to combine structural models with mechanical properties and preparation history, thereby facilitating the development of preparation-structure-property relationships for metallic glasses. With the availability of such information, bulk metallic glasses can be further optimized to be used in more and more applications in industry.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Fourier analysis of continuous fractional diffusion synthetic acceleration schemes in slab geometry

We propose two fractional extensions of continuous diffusion synthetic acceleration (DSA) with fractional derivative order α varying over the interval 2 ≥ α ≥ 1 . We investigate the spectral properties of the corresponding continuous families of fractional preconditioners by performing Fourier analysis for a model infinite homogeneous medium problem in slab geometry. The first family results in a fractional acceleration scheme, FrDSAo, that reduces to traditional DSA for .α = 2 and scattering ratio c limiting to a unit value (c → 1) but is otherwise optimized via the Fourier analysis, to obtain the smallest possible spectral radius, for c < 1 and 2 ≥ α ≥ 1. The second family corresponds to a fractional acceleration scheme, FrDSAs, that reduces to traditional DSA for α = 2 for all values of c. The latter scheme is not optimized but has the advantage of lending itself to a more straightforward implementation. For high values of c, the results of the Fourier analysis point to the existence of an interval 2 > α > ∼1.8 where both FrDSAo and FrDSAs can achieve a lower spectral radius than DSA. For example, DSA has a spectral radius of ∼0.2246 for c = 0.9999 while FrDSAo produces a value of ∼0.1616 at α = 1.92 and FrDSAs results in ∼0.2116 at α =1.93. (author)

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Design of Efficient Molecular Electrocatalysts for Water and Carbon Dioxide Reduction Using Predictive Models of Thermodynamic Properties

Objectives: 1) To design and synthesize water soluble transition metal complexes that function as electrocatalysts for the reduction of water to hydrogen or carbon dioxide to formate 2) To experimentally measure thermochemical properties, such as hydride donor ability, and use these quantities to optimize catalyst performance in solutions of varying pH 3) To incorporate results from detailed mechanistic and kinetic studies to direct improvements in catalyst design 4) To use thermochemical data from objective 2 to generate predictive models for the properties of new complexes to apply to the development of catalysts for other reductive reactions, with an emphasis on carbon dioxide reduction to methanol Description: Wide-spread implementation of renewable but intermittent energy sources, such as solar, requires the development of efficient methods for energy storage. The high energy density of chemical bonds makes chemical fuels an ideal solution for energy storage. Hydrogen and reduced carbon compounds have been proposed as ideal candidates for chemical fuels. However, the generation of chemical fuels from electricity requires competent electrocatalysts. The proposed research focuses on developing electrocatalysts for the reduction of water to hydrogen, and carbon dioxide to formate. Both products can be used directly as an energy carrier in fuel cells, or as a reductant for more saturated chemical fuels. Hydrogen is the most common fuel used in current commercial fuel cells. Additionally, there is interest in formic acid as a liquid carrier for hydrogen because of its increased storage density and ease of dehydrogenation. Formate is also an intermediate in the sequential reduction of CO 2 to methanol, another potential chemical fuel with high energy storage density. Despite the utility of formate as a chemical fuel or precursor, there are very few examples of electrocatalysts for its production from CO 2 , and even fewer demonstrate high product selectivity. In heterogeneous catalysis, the Sabatier principle is used to generate volcano curves that describe optimal thermochemical properties for key intermediates that result in peak catalytic activity. In most cases, such as hydrogen production and oxidation, the most favorable metals (Pt, Re, Rh, and Ir) are rare and expensive. Molecular inorganic complexes provide an opportunity to use electronic and steric ligand effects to tune the critical thermodynamic parameters of abundant metals to values comparable to key surface intermediates on precious metals. This principle will be applied to the design of aqueous homogeneous catalysts for the reduction of H 2 O and CO 2 optimized to function at specific pH ranges. The critical intermediate in both of these reactions is a metal hydride. The strength of this bond, or hydricity (ΔGH - ) dictates the overall thermodynamics of the reduction of H + to H 2 and the sequential reduction of C1 substrates, such as CO 2 , CO, and H 2 CO (the latter to CH 3 OH). ΔGH - will be systematically measured for a series of first row metal complexes to form predictive models for metal and ligand electronic effects. Benefits and Outcomes: This approach addresses thermodynamic requirements essential for energy-efficient catalyst design and transition state barrier considerations for faster catalysis. The advancements made from this research will include the discovery of efficient, stable, and fast catalysts for hydrogen and formate generation from electricity. The ability to store energy efficiently is critical to the widespread use of renewable energy schemes. This provides a viable solution to produce high energy density fuels for stationary and mobile energy applications. The fundamental research conducted is also directly applicable to the design of catalysts for the generation of further reduced products such as methanol.

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

On the force exerted on a non-spherical asymmetric dust grain from homogeneous, stationary, isotropic, non-magnetized plasma

Whereas the conventional wisdom suggests that the force between non-magnetized homogeneous, stationary, isotropic plasma, and the dust grain is only possible for the case of relative plasma–grain velocity, it is shown that stationary non-spherical asymmetric dust grain immersed in stationary, non-magnetized, isotropic plasma can experience a force caused by the grain–plasma interactions. The component of the force due to scattering of plasma particles in the limit of infinite Debye length is considered analytically. Both the particle scattering and absorption force components are modeled numerically in the limits of infinite and finite Debye length using a newly developed 2D3V Aspherical Particle-in-Cell code. The code simulates interactions of dust grain of selected non-spherical asymmetric shape with plasmas using dust shape conforming coordinates. Here, the simulations confirm the existence of the force on non-spherical asymmetric grain in stationary non-magnetized plasma and show that the plasma screening effects can lead to reversal of the force direction.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Feasibility Study for a Subcritical Assembly at the Oak Ridge National Laboratory

A conceptual design has been completed for a new subcritical assembly at Oak Ridge National Laboratory (ORNL). This subcritical assembly will support the Nuclear Criticality Safety Program (NCSP) training and education program to augment current subcritical assembly capabilities at Sandia National Laboratories (SNL) and at the National Criticality Experiments Research Center (NCERC). The proposed subcritical assembly uses legacy AGN-201M research reactor fuel plates that are available at the Y-12 National Security Complex. This subcritical assembly contains approximately 620 grams of 235 U as UO 2 particulates distributed homogeneously in radiation-stabilized polyethylene. The fuel plates will have a graphite neutron reflector to obtain a core multiplication, $M$, from 10 to 50, corresponding to a k eff of 0.9 to 0.95, respectively. The subcritical assembly will be able to support at least four experiments for the training courses: the addition of fissile material to increase neutron multiplication (mass), a core separation experiment (interaction), the effect of adding moderators to the core (moderation), and the effect of adding neutron absorbers to the core (poison/absorption). The proposed ORNL subcritical assembly will provide additional student bandwidth to support NCSP training courses and to provide for some backup capacity in the event that nuclear facility operations are disrupted at SNL or NCERC due to some unforeseen issue. The subcritical assembly will be designed to be an inherently safe subcritical assembly—subcritical under all normal and abnormal conditions—and it will provide the capability to conduct hands-on training to support NCSP and general nuclear criticality safety (NCS) staff training and qualification goals.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

A Subcritical Assembly for Training and Education Use at the Oak Ridge National Laboratory

A near-final design has been completed for a new subcritical assembly at Oak Ridge National Laboratory (ORNL). This subcritical assembly will support the Nuclear Criticality Safety Program (NCSP) training and education program to augment current subcritical assembly capabilities at Sandia National Laboratories (SNL) and at the National Criticality Experiments Research Center (NCERC). Other collaborations have been proposed with the Sunrise Consortium and others. This proposed subcritical assembly uses legacy AGN-201M research reactor fuel plates that are available from the Y-12 National Security Complex. This subcritical assembly contains approximately 620 grams of 235 U as UO 2 particulates distributed homogeneously in radiation-stabilized polyethylene. The fuel plates will have a graphite neutron reflector to obtain a core multiplication, M, from 10 to 50, corresponding to a k eff of 0.9 to 0.95, respectively. The subcritical assembly will be able to support at least four experiments for the training courses: (1) the addition of fissile material to increase neutron multiplication (mass), (2) a core separation experiment (interaction), (3) the effect of adding moderators to the core (moderation), and (4) the effect of adding neutron absorbers to the core (poison/absorption). The proposed ORNL subcritical assembly will provide additional student bandwidth to support NCSP training courses and to provide for some backup capacity if nuclear facility operations are disrupted at SNL or NCERC due to some unforeseen issue. The subcritical assembly will be designed to be an inherently safe subcritical assembly—subcritical under all normal and abnormal conditions—and it will provide the capability to conduct hands-on training to support NCSP and general nuclear criticality safety (NCS) staff training and qualification goals.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2026

Phase field dislocation dynamics formulation coupled with Fourier based micromechanics solver and its application to grain boundary–dislocation interactions

A new phase field dislocation dynamics (PFDD) formulation for homogeneous and heterogeneous materials is presented, which couples micromechanical solvers and the time-dependent Ginzburg–Landau equation. The strain fields are obtained from the micromechanical solver by solving the Lippmann–Schwinger equation and then used to define energy terms to model the evolution of the dislocations. Grain boundary (GB)–dislocation interactions are studied using the coupled PFDD formulation and by describing GBs as inclusions. GB energy and stiffness tensors are computed from molecular statics simulations, and a newly proposed lattice energy term that is dependent on the GB energy is considered in the calculations. Interaction of a screw dislocation with minimum energy and metastable states of low and high angle ⟨110⟩ symmetric tilt grain boundaries are studied. We show good agreement between predictions from our PFDD formulation and molecular dynamics simulations of grain boundary–dislocation interactions.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Can Cooling and Heating Functions Be Modeled with Homogeneous Radiation Fields?

Abstract Cooling and heating functions describe how radiative processes impact the thermal state of a gas as a function of its temperature and other physical properties. In a most general case the functions depend on the detailed distributions of ionic species and on the radiation spectrum. Hence, these functions may vary on a very wide range of spatial and temporal scales. In this paper, we explore cooling and heating functions between 5 ≤ z ≤ 10 in simulated galaxies from the Cosmic Reionization On Computers project. We compare three functions: (1) the actual cooling and heating rates of hydrodynamic cells as a function of cell temperature, (2) the median cooling and heating functions computed using median interstellar medium (ISM) properties (median ISM), and (3) the median of the cooling and heating functions of all gas cells (instantaneous). We find that the median ISM and instantaneous approaches to finding a median cooling and heating function give identical results within the spread due to cell-to-cell variation. However, the actual cooling (heating) rates experienced by the gas at different temperatures in the simulations do not correspond to either summarized cooling (heating) functions. In other words, the thermodynamics of the gas in the simulations cannot be described by a single set of a cooling plus a heating function with a spatially constant radiation field that could be computed with common tools, such as CLOUDY.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Evolution of the superconducting properties from binary to ternary APC-Nb 3 Sn wires

We present a study conducted on binary Tube Type and ternary powder-in-tube Nb 3 Sn wires manufactured using the artificial pinning centres-internal oxidation method. All the specimens are doped with Zr: oxide nano-particles of this element are responsible for the pinning improvement, both by refining the A-15 grain-size and their own point-pinning contribution. Low-field J c magnetometry confirms that the hadron-hadron Future Circular Collider (FCC-hh) specifications are met by one ternary doped-sample. The differences in microstructure were assessed by scanning electron microscopy/transmission electron microscopy to clarify the reasons for the pinning improvement between the two generations. The deviations from the Dew Hughes model are also discussed, underlying some non-linear addition due to competition between the two pinning mechanisms. Finally, we show how the introduction of Ta as a ternary addition influences the A-15 phase by focusing on the radial inhomogeneities, evaluating the T c distribution and Sn composition gradients. The latter are used to model the currents, enabling us to evaluate the individual weights of the pinning mechanisms and their absolute contributions at the High Luminosity-Large Hadron Collider and FCC-hh dipoles operational points.

43 PARTICLE ACCELERATORS↗
DOE OSTI2025

Benders Decomposition Using Graph Modeling and Multi-Parametric Programming

Benders decomposition is a widely used method for solving large and structured optimization problems, but its performance is affected by the repeated solution of subproblems. We propose a flexible and modular algorithmic framework for accelerating Benders decomposition. Specifically, we express the problem structure by using a graph-theoretic modeling abstraction in which nodes represent optimization subproblems and edges represent connectivity between subproblems. A key innovation of our approach is that we embed multiparametric programming (mp) surrogates for node subproblems, which maps the exact analytical map of the subproblem solution space. The use of mp surrogates allows us to replace subproblem solves with fast look-ups and function evaluations for primal and dual variables during the iterative Benders process. We formally show the equivalence between classical Benders cuts and those derived from the mp solution. We implement our framework in the open-source PlasmoBenders.jl software package. To demonstrate the capabilities of the proposed framework, we apply it to a two-stage stochastic programming problem, which aims to make optimal capacity expansion decisions under market uncertainty. We evaluate both single-cut and multicut variants of Benders decomposition and show that the use of mp surrogates achieves substantial speedups in subproblem solve time, while preserving the convergence guarantees of Benders decomposition. We highlight advantages in solution analysis and interpretability that is enabled by mp critical region tracking; specifically, we show that these reveal how decisions evolve geometrically across the Benders search. Our results aim to demonstrate that combining surrogate modeling with graph modeling offers a promising and extensible foundation for structure-exploiting decomposition. In addition, by decomposing the problem into more tractable subproblems, the proposed approach also aims to overcome scalability issues of mp. Finally, the use of mp surrogates provides a unifying and modular optimization framework that enables the representation of heterogeneous node subproblems as modeling objects with a homogeneous structure.

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

Towards a self-consistent model of the convective core boundary in upper main sequence stars: I. 2.5D and 3D simulations

There is strong observational evidence that the convective cores of intermediate-mass and massive main sequence stars are substantially larger than those predicted by standard stellar-evolution models. However, it is unclear what physical processes cause this phenomenon or how to predict the extent and stratification of stellar convective boundary layers. Convective penetration is a thermal-timescale process that is likely to be particularly relevant during the slow evolution on the main sequence. We use our low-Mach-number S EVEN -L EAGUE H YDRO code to study this process in 2.5D and 3D geometries. Starting with a chemically homogeneous model of a 15 M ⊙ zero-age main sequence star, we construct a series of simulations with the luminosity increased and opacity decreased by the same factor, ranging from 10 3 to 10 6 . After reaching thermal equilibrium, all of our models show a clear penetration layer; its thickness becomes statistically constant in time and it is shown to converge upon grid refinement. The penetration layer becomes nearly adiabatic with a steep transition to a radiative stratification in simulations at the lower end of our luminosity range. This structure corresponds to the adiabatic ‘step overshoot’ model often employed in stellar-evolution calculations. The simulations with the highest and lowest luminosity differ by less than a factor of two in the penetration distance. The high computational cost of 3D simulations makes our current 3D data set rather sparse. Depending on how we extrapolate the 3D data to the actual luminosity of the initial stellar model, we obtain penetration distances ranging from 0.09 to 0.44 pressure scale heights, which is broadly compatible with observations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

PIONEER, a high-resolution single-crystal polarized neutron diffractometer

PIONEER is a high Q-resolution, single-crystal, polarized neutron diffractometer at the Second Target Station (STS), Oak Ridge National Laboratory. It will provide the unprecedented capability of measuring tiny crystals (0.001 mm 3 , i.e., x-ray diffraction size), ultra-thin films (10 nm thickness), and weak structural and magnetic transitions. PIONEER benefits from the increased peak brightness of STS cold-neutron sources and uses advanced Montel mirrors that are able to deliver a focused beam with a high brilliance transfer, a homogeneous profile, and a low background. Monte Carlo simulations suggest that the optimized instrument has a high theoretical peak brilliance of 2.9 × 10 12 n cm -2 sr -1 Å -1 s -1 at 2.5 Å at the sample position, within a 5 × 5 mm 2 region and a ±0.3° divergence range. The moderator-to-sample distance is 60 m, providing a nominal wavelength band of 4.3 Å with a wavelength resolution better than 0.2% in the wavelength range of 1.0–6.0 Å. PIONEER is capable of characterizing large-scale periodic structures up to 200 Å. With a sample-to-detector distance of 0.8 m, PIONEER accommodates various sample environments, including low/high temperature, high pressure, and high magnetic/electric field. Additionally, a large cylindrical detector array (4.0 sr) with a radial collimator is planned to suppress the background scattering from sample environments. Bottom detector banks provide an additional 0.4 sr coverage or can be removed if needed to accommodate special sample environments. We present virtual experimental results to demonstrate the scientific performance of PIONEER in measuring tiny samples.

47 OTHER INSTRUMENTATION↗