Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Objective function”

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 595 records · Page 33

A novel machine learning-based optimization algorithm (ActivO) for accelerating simulation-driven engine design

A novel design optimization approach (ActivO) that employs an ensemble of machine learning algorithms is presented. The proposed approach is a surrogate-based scheme, where the predictions of a weak leaner and a strong learner are utilized within an active learning loop. The weak learner is used to identify promising regions within the design space to explore, while the strong learner is used to determine the exact location of the optimum within promising regions. For each design iteration, exploration is done by randomly selecting evaluation points within regions where the weak learner-predicted fitness is high. The global optimum obtained by using the strong learner as a surrogate is also evaluated to enable rapid convergence once the most promising region has been identified. First, the performance of ActivO was compared against five other optimizers on a cosine mixture function with 25 local optima and one global optimum. In the second problem, the objective was to minimize indicated specific fuel consumption of a compression-ignition internal combustion (IC) engine while adhering to desired constraints associated with in-cylinder pressure and emissions. In this work, the efficacy of the proposed approach is compared to that of a genetic algorithm, which is widely used within the internal combustion engine community for engine optimization, showing that ActivO reduces the number of function evaluations needed to reach the global optimum, and thereby time-to-design by 80%. Furthermore, the optimization of engine design parameters leads to savings of around 1.9% in energy consumption, while maintaining operability and acceptable pollutant emissions.

97 MATHEMATICS AND COMPUTING↗

The luminosity function of ringed galaxies

ABSTRACT We perform an analysis of the luminosity functions (LFs) of two types of ringed galaxies – polar-ring galaxies and collisional ring galaxies – using data from the Sloan Digital Sky Survey (SDSS). Both classes of galaxies were formed as a result of interaction with their environment and they are very rare objects. We constructed LFs of galaxies by different methods and found their approximations by the Schechter function. The luminosity functions of both types of galaxies show a systematic fall-off at low luminosities. The polar structures around bright (Mr ≤ −20m) and red (g − r > +0.8) galaxies are about twice as common as around blue ones. The LF of collisional rings is shifted towards brighter luminosities compared to polar-ring galaxies. We analysed the published data on the ringed galaxies in several deep fields and confirmed the increase in their volume density with redshift: up to z ∼ 1 their density grows as (1 + z)m, where m ≳ 5.

79 ASTRONOMY AND ASTROPHYSICS↗

Threshold resummation for computing large-x parton distribution through large-momentum effective theory

Parton distribution functions (PDFs) at large x are poorly constrained by high-energy experimental data, but extremely important for probing physics beyond standard model at colliders. We study the calculation of PDFs at large-x through large-momentum P z expansion of the lattice quasi PDFs. Similar to deep-inelastic scattering, there are two distinct perturbative scales in the threshold limit where the matching coefficient can be factorized into a space-like jet function at scale P z |1 – y| and a pair of heavy-light Sudakov form factors at scale P z . The matching formula allows us to derive a full renormalization group resummation of large threshold logarithms, and the result is consistent with the known calculation to the next-to-next to leading order (NNLO). This paves the way for direct large-x PDFs calculations in lattice QCD. As by-products, we find that the space-like jet function is related to a time-like version calculated previously through analytic continuation, and the heavy-light Sudakov form factor, calculated here to NNLO, is a universal object appearing as well in the large momentum expansion of quasi transverse-momentum-dependent PDFs and quasi wave-function amplitudes.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Multi-scale Dynamics of Water Regulation by Bacteria in Synthetic Soil Microsystems (Final Technical Report)

The objective of this project was to develop a multifunctional in situ platform to link gene expression, regulation, and function with the spatiotemporal context of the whole cell environment. In particular, we focus on a synthetic rhizosphere microhabitat and spatial regulation of genes encoding production of exopolysaccharides, or EPS. We have fully met this objective and report results that support the hypothesis: microbial EPS production, acting at pore throats, limits water evaporation and promotes resilience of terrestrial soil systems.

42 ENGINEERING↗

Structure-aware methods for expensive derivative-free nonsmooth composite optimization

We present new methods for solving a broad class of bound-constrained nonsmooth composite minimization problems. These methods are specially designed for objectives that are some known mapping of outputs from a computationally expensive function. We provide accompanying implementations of these methods: in particular, a novel manifold sampling algorithm (MS-P) with subproblems that are in a sense primal versions of the dual problems solved by previous manifold sampling methods and a method (GOOMBAH) that employs more difficult optimization subproblems. For these two methods, we provide rigorous convergence analysis and guarantees. We demonstrate extensive testing of these methods. Open-source implementations of the methods developed in this manuscript can be found at https://github.com/POptUS/ IBCDFO/.

97 MATHEMATICS AND COMPUTING↗

or-topas: Operations Research Toolkit for Pyomo Alternative Solutions

SAND2026-16702O OR-TOPAS: Operations Research Toolkit for Pyomo Alternative Solutions is a tool that enhances optimization applications defined by the Pyomo modeling library. It offers functions to generate optimal or near-optimal solutions, operating independently of Pyomo’s solver interface. Users can configure these functions with specific solver names and options, resulting in a custom solution object that returns a list of solutions. The OR-TOPAS library includes methods tailored to the properties of the model, such as binary integer programs versus linear programs, and specific solver interfaces like Gurobi. It does not provide models for specific applications, but it is applicable to a wide range of Pyomo optimization models. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.

Siirola, John [Sandia National Lab. (SNL-CA), Live↗

Heterologous expression of Arabidopsis laccase2, laccase4 and peroxidase52 driven under developing xylem specific promoter DX15 improves saccharification in populus

Secondary cell wall holds considerable potential as it has gained immense momentum to replace the lignocellulosic feedstock into fuels. Lignin one of the components of secondary cell wall tightly holds the polysaccharides thereby enhancing the recalcitrance and complexity in the biomass. Laccases (LAC) and peroxidases (PRX) are the major phenyl-oxidases playing key functions during the polymerization of monolignols into lignin. Yet, the functions of laccase and peroxidases gene families remained largely unknown. Hence, the objective of this conducted study is to understand the role of specific LAC and PRX in Populus wood formation and to further investigate how the altered Lac and Prx expression affects biomass recalcitrance and plant growth. This study of heterologous expression of Arabidopsis Lac and Prx genes was conducted in poplar to avoid any otherwise occurring co-suppression mechanism during the homologous overexpression of highly expressed native genes. In the pursuit of optimizing lignocellulosic biomass for biofuel production, the present study focuses on harnessing the enzymatic potential of Arabidopsis thaliana Laccase2, Laccase4, and Peroxidase52 through heterologous expression. We overexpressed selected Arabidopsis laccase2 (AtLac2), laccase4 (AtLac4), and peroxidase52 (AtPrx52) genes, based on their high transcript expression respective to the differentiating xylem tissues in the stem, in hybrid poplar (cv. 717) expressed under the developing xylem tissue-specific promoter, DX15 characterized the transgenic populus for the investigation of growth phenotypes and recalcitrance efficiency. Bioinformatics analyses conducted on AtLac2 and AtLac4 and AtPrx52, revealed the evolutionary relationship between the laccase gene and peroxidase gene homologs, respectively. Transgenic poplar plant lines overexpressing the AtLac2 gene (AtLac2-OE) showed an increase in plant height without a change in biomass yield as compared to the controls; whereas AtLac4-OE and AtPrx52-OE transgenic lines did not show any such observable growth phenotypes compared to their respective controls. The changes in the levels of lignin content and S/G ratios in the transgenic poplar resulted in a significant increase in the saccharification efficiency as compared to the control plants. Overall, saccharification efficiency was increased by 35–50%, 21–42%, and 8–39% in AtLac2-OE, AtLac4-OE, and AtPrx52-OE transgenic poplar lines, respectively, as compared to their controls. Moreover, the bioengineered plants maintained normal growth and development, underscoring the feasibility of this approach for biomass improvement without compromising overall plant fitness. This study also sheds light on the potential of exploiting regulatory elements of DX15 to drive targeted expression of lignin-modifying enzymes, thereby providing a promising avenue for tailoring biomass for improved biofuel production. These findings contribute to the growing body of knowledge in synthetic biology and plant biotechnology, offering a sustainable solution to address the challenges associated with lignocellulosic biomass recalcitrance.

09 BIOMASS FUELS↗

Rupture Dynamics Simulations of Shallow Crustal Earthquakes on Reverse Slip Faults

The main goal of this study was to perform dynamic rupture simulations in order constrain the characterization of kinematic rupture models implemented in techniques for simulating strong ground motion for crustal earthquakes. First, we investigated the rupture process of the 2008 Iwate-Miyagi Nairiku earthquake by developing a dynamic rupture model on a reverse fault using a trial-and-error technique that produced a slip distribution and near-fault ground motion that matched the recorded ones. The simulations were performed in the frequency range 0-2 Hz, using a 3D staggered grid finite-difference method and a linear slip weakening friction law. Constrained by the observed slip distribution and consistent with dynamic rupture models, the derived kinematic rupture model of the Iwate-Miyagi Nairiku earthquake contains areas with relatively large slip rate, representing strong motion generation areas, set against lower amplitude heterogeneous background slip, and a relatively low slip rate in the weak zone of the top-most crust (upper 3km). Second, we performed rupture dynamics modeling to constrain shallow slip characterization in rupture models for strike-slip crustal earthquakes. The objective was to establish general rules about the characterization of slip rate function and slip at shallow depths and in the strong motion generation areas (SMGAs). The simulations of spontaneous rupture were performed in the frequency range 0-2.0Hz, using a 3D staggered-grid finite-difference method and a layer over half-space 1D crustal velocity model with a minimum shear-wave velocity of 2.8 km/s. In order to account for changes in material ductility and reduction of stress drop, observed in the shallow crust (upper 3-5 km), and the transition from ductile state to brittle state in the upper seismogenic zone, in our stress models we included a shallow weak zone (<4km). In this zone the stress drop was set to zero at the free surface and gradually increased with depth, while the slip weakening distance was set to 75 cm at the free surface and decreased to 50 cm at the base of the weak zone. As in the case of reverse faulting during the Iwate-Miyagi Nairiku earthquake, from these computations we found a systematic change in the shape of the slip-rate function, from Kostrov-type in the deeper part of the fault to a more symmetric cosine-type in the upper few kilometers, near the free-surface. Moreover, the average slip duration in the weak zone, with respect to slip duration in the deeper parts of the fault, increases by at most a factor of 1.5. We found a systematic and gradual change in the shape of the slip-rate function from Kostrov-type in the asperity areas (SMGAs), to more symmetric cosine-type in the upper few km near the free-surface, on the long period motion generation area (LMGA). The effective rise time in the LMGA, with respect to that in the SMGAs, increases by at most a factor of 2. Effective rise time is the time difference between the time at which the slip rate drops to a level that is equal to 25% of its peak and the onset time of the signal. In addition, the slip in the LMGA located above the SMGAs is almost the same as the one in the SMGAs, and about 1.5 times larger than the average slip.

58 GEOSCIENCES↗

First-Principles Studies of Tritium Species Diffusivity Across the Interface of Nickel-Plated Zircaloy-4

In this study, we will calculate the 3H2 and 3H2O diffusion mechanism in the Ni coated layer and the dominant 3H species dissolution into the Zircaloy-4 getter. We will optimize the Ni/Zr interface and introduce Sn and O impurities. To achieve the objectives, we will perform the investigation by using spin-polarized density functional theory (DFT) methods to capture the effects of Ni magnetic properties on 3H diffusion. We will introduce NiO/Ni(OH) on the Ni side closer to the interface. Then, the 3H diffusion pathways and barriers from Ni side to the zircalory-4 will be explored. We will investigate how NiO/Ni(OH) affects the 3H diffusion across the interface.

Paudel, Hari P.↗

Energy Clusters Offshore: A Technology Feasibility Review

Globally, governments, companies, and other organizations have committed to achieving net-zero emissions targets in the coming decades. To achieve decarbonization at the scale and pace required to meet these targets, future energy systems will need renewable energy to serve 100% of the existing direct electricity demand, support additional electrification, and decarbonize the wider economy. An energy cluster offshore (ECO) is a concept that seeks to meet this challenge by integrating and optimizing large-scale renewable electricity generation, storage, and fuel production technologies and pairing them with other complementary uses, such as carbon capture or water desalination. This research project explores the techno-economic feasibility of ECO concepts by taking a holistic view of complex, multidisciplinary hybrid plant designs while considering different configurations and objectives. We will outline the most promising technology combinations and configurations, their functional requirements, and opportunities for optimization.

14 SOLAR ENERGY↗

Modeling the Effects of Global Change on Ecosystem Processes in a Tropical Rainforest

Ongoing land-use change and climate change in wet tropical forests can potentially drive shifts in tree species composition, representing a change in individual species within a functional group, tropical evergreen trees. The impacts on the global carbon cycle are potentially large, but unclear. We explored the differential effects of species within this functional group, in comparison with the effects of climate change, using the Century model as a research tool. Simulating effects of individual tree species on biome-level biogeochemical cycles constituted a novel application for Century. Background and Objectives: A unique, long-term, replicated field experiment containing five evergreen tree species in monodominant stands under similar environmental conditions in a Costa Rican wet forest provided data for model evaluation. Our objectives were to gain insights about this forest’s biogeochemical cycles and effects of tree species within this functional group, in comparison with climate change. Materials and Methods: We calibrated Century, using long-term meteorological, soil, and plant data from the field-based experiment. In modeling experiments, we evaluated effects on forest biogeochemistry of eight plant traits that were both observed and modeled. Climate-change simulation experiments represented two climate-change aspects observed in this region. Results: Model calibration revealed that unmodeled soil processes would be required to sustain observed P budgets. In species-traits experiments, three separate plant traits (leaf death rate, leaf C:N, and allocation to fine roots) resulted in modeled biomass C stock changes of >50%, compared with a maximum 21% change in the climate-change experiments. Conclusions: Modeled ecosystem properties and processes in Century were sensitive to changes in plant traits and nutrient limitations to productivity. Realistic model output was attainable for some species, but unusual plant traits thwarted predictions for one species. Including more plant traits and soil processes could increase realism, but less-complex models provide an accessible means for exploring plant-soil-atmosphere interactions.

54 ENVIRONMENTAL SCIENCES↗

MOOSE Reactor Module Meshing Enhancements to Support Reactor Physics Analysis

The U.S. Department of Energy Office of Nuclear Energy Advanced Modeling and Simulation (NEAMS) program develops an integrated suite of advanced reactor physics tools built upon the Multiphysics Object-Oriented Simulation Environment (MOOSE) framework. Each code generally requires an input finite element mesh on which the physics solution is calculated, reported, and transferred to other physics codes. The meshing process is often burdensome for the complex geometries present in reactors due to lack of easy-to-use, open-source meshing tools. To address the bottleneck associated with meshing complex geometries found in nuclear reactors, NEAMS initiated the development of the MOOSE Reactor Module starting in FY21. The Reactor Module builds off the existing MOOSE Mesh System to include targeted meshing capabilities such as the ability to generate hexagonal pin cells, assemblies with ducts, rotating control drums, cores, peripheral zones around a core, as well as the automatic labeling (“reporting IDs”) of pin, assembly, and planar zones to simplify post-processing of results. As a Physics Module in MOOSE, the Reactor Module is open-source, available with any MOOSE installation, directly compatible with MOOSE-based tools, and can be invoked from MOOSEbased applications to generate meshes. Functionality from the Reactor Module has been applied to several advanced reactor concepts to demonstrate user workflow improvements and accuracy. The primary objective of the Reactor Module is to improve useability of MOOSEbased tools by streamlining mesh generation and output inspection processes. During FY22, the functionality of the Reactor Module (and accompanying Mesh System) has been expanded based on user needs. First, the Reactor Geometry Mesh Builder capability developed primarily in FY21 has been refactored and merged to the public MOOSE repository. This capability wraps underlying Reactor Module mesh generators into a “Pin – Assembly – Core” workflow appropriate for conventional Cartesian and hexagonal geometries, and notably assigns material IDs during mesh generation stage and generates only the minimal number of blocks needed in order to reduce computational burden. Biasing and boundary layer options have been added to the base mesh generators as required by thermal hydraulics solvers. The reporting ID functionality has been expanded to differentiate ring-wise and azimuthal sectors within a pin for use with depletion algorithms, and VectorPostProcessor and Reporter objects are now available to integrate solution variables across zones based on ID combinations. Functionality to trim hexagonal meshes along the center or periphery has been developed so users may leverage symmetry and reflective boundary conditions to reduce the mesh size. A flexible and powerful tool to fill the space between two sidesets has been introduced to the framework and can be used for transition layers such as stitching two assemblies together with different numbers of pins, or for complex geometries which do not follow conventional Cartesian/hexagonal patterns. Finally, additional verification problems were performed with NEAMS physics tools in complement with existing NEAMS work.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Application of Cable Condition Monitoring Technologies to Assess Age-Related Degradation of Industrial Cables Installed in Harsh Environments

The aging of electrical cables has been the subject of substantial research and development (R&D) projects performed by national and international laboratories, universities, and private organizations for many years. This R&D was conducted to develop guidance, equipment, and techniques to support aging management of in-service cables in industrial facilities such as nuclear power plants, research reactors, waste facilities, and fuel fabrication plants. Through these research efforts, condition monitoring technologies have been developed that can determine the severity of age-related degradation that occurs in industrial cables and insulation polymers during service. This paper summarizes the results of aging assessments that were performed for cables installed in two U.S. nuclear power plants, one a pressurized water reactor and one a boiling water reactor. These cables had been in service for over 40 years and during plant operation were exposed to harsh environmental conditions including elevated temperatures and radiation. For these assessments, a comprehensive series of measurements was performed to assess the aged condition of the cables. These cables came from different manufacturers, were manufactured in different years, and were constructed with a variety of jacket and insulation polymers including chloro-sulfonated polyethylene (CSPE), cross-linked polyethylene (XLPE)/cross-linked polyolefin (XLPO), neoprene, and ethylene propylene rubber (EPR). The goal of these assessments was to determine the current aged condition of the cable polymers and provide an estimate of how long the cable insulation materials could remain exposed to their in-service environmental conditions before reaching their end-of-life condition. Both nuclear power plants have received license renewals to extend their operation from 40 to 60 years, and the utilities need objective evidence to show that critical components such as cables will be able to function safely and reliably during the extended operating period. Furthermore, the results of these assessments showed that the cables exhibited different aged conditions depending on the type of polymers they were constructed with and the environment they were exposed to during service. Some of the cables and insulation polymers showed signs of significant age-related degradation and were estimated to have approximately 5 years of remaining service life. Other cables exhibited no signs of significant age-related degradation and were estimated to have 50 years or more of remaining service life. Using the results of these cable aging assessments, plant personnel were able to (1) determine the overall aged condition of cables and insulation polymers using objective test results, (2) identify aged or degraded cables before they caused operability issues, and (3) avoid unnecessary and costly replacement of cables that can continue to operate safely and reliably.

36 MATERIALS SCIENCE↗

Assembly and Functional Test of NRAD Heated Instrumentation Rig

This is a technical report (TEV-4296) which documents the configuration and materials of construction for a heated sensor test conducted in the NRAD reactor. It is being reviewed for external release so that operators of university reactors can evaluate whether this level of documentation will be acceptable for similar tests conducted in their reactors. The Objective section of this TEV is reproduced below. The purpose of this TEV is to document assembly and functional testing of a device which is designed to allow advanced sensors to be temporarily placed at the periphery of the NRAD reactor. While in this location, the operational characteristics of the sensors can be evaluated while they are exposed to a neutron flux. This device (or “rig”) consists of a cartridge heater surrounded by six guide tubes, which are in turn surrounded by insulation and an aluminum dry tube. The heater allows for the sensors to be tested while at elevated temperature.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Acceleration of Pd-V intermetallic diffusion by hydrogen

Vanadium-based membranes have great potential for hydrogen purification due to their perfect selectivity, high permeability, and relatively low cost. With appropriate surface cleaning, V efficiently permeates hydrogen at elevated temperature, but performance declines due to its affinity to absorb impurities. Here, the application of palladium thin films maintains a clean surface that catalyzes hydrogen dissociation and recombination. Hydrogen permeation in Pd-V-Pd membranes initially reach theoretical permeability, but declines due to Pd-V interdiffusion. The objective of this work was to quantify the intermetallic diffusion process as a function of temperature and ambient. Pd-V composites were subjected to various annealing treatments and characterized using Auger electron spectroscopy, X-ray diffraction, and energy dispersive X-ray spectroscopy, as well as correlated to measurements of membrane permeability. In an inert environment Pd-V interdiffusion was observable as low as 300 °C, and the diffusion coefficient had an activation energy of 44 kJ/mol. Furthermore, the presence of hydrogen at partial pressures > 10 kPa accelerated interdiffusion six-fold at T = 400 °C. Membrane performance degraded with an activation energy 75 kJ/mol, suggesting that intermetallic diffusion leads to both a loss of catalytic activity and as well as degradation of bulk permeability. These findings provide a baseline for evaluating hydrogen permeable interdiffusion barriers to overcome these challenges.

08 HYDROGEN↗

Finding the Pareto front for high-entropy-alloy catalysts

Finding catalysts that have both high activity and high stability presents a long-standing challenge. Since optimizing activity and stability are conflicting objectives, the best one can do is find the Pareto front that yields optimal tradeoffs between these features. On the Pareto front, there is a trade-off where a portion of catalytic activity must be sacrificed to gain further stability and vice versa . Here, we provide a method to optimize the front by designing a multi-objective genetic algorithm that combines machine learning, graph neural network calculations, and density functional calculations. The application considered is the oxygen evolution reaction catalyzed by high-entropy alloys. We find that the Pareto front generally contains alloys with diverse elements, but that enhancing stability inevitably inflicts a toll on activity. We compare the general conclusions of our work to a survey of 545 experiments.

Zhang, Chengyi [Univ. of Auckland (New Zealand)]↗

Adaptive 3D convolutional neural network-based reconstruction method for 3D coherent diffraction imaging

We present a novel adaptive machine-learning based approach for reconstructing three-dimensional (3D) crystals from coherent diffraction imaging. We represent the crystals using spherical harmonics (SH) and generate the corresponding synthetic diffraction patterns. We utilize 3D convolutional neural networks (CNNs) to learn a mapping between 3D diffraction volumes and the SH, which describe the boundary of the physical volumes from which they were generated. We use the 3D CNN-predicted SH coefficients as the initial guesses, which are then fine-tuned using adaptive model-independent feedback for improved accuracy. We also adaptively tune the locations, intensities, and decay rates of collections of radial basis functions in order to reproduce the non-uniform internal structure of 3D objects and demonstrate the method for a synthetic volume that has an internal void and a density ramp.

36 MATERIALS SCIENCE↗

Real Vector Framework

SAND2025-11463O The Real Vector Framework (RVF) is a modern and flexible C++ vector math library for developing scientific computing software that involves vector computations. RVF allows an opt-in approach to functionality that parallels the familiar base-class and override structures of object-oriented programming. Users can reuse and customize the code without inheritance entanglements and dynamic dispatch, while enabling seamless interoperability between diverse container types. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.

von Winckel, Gregory [Sandia National Lab. (SNL-CA↗