Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “parameter optimization”

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 703 records · Page 39

Targeted Adaptive Design

Modern advanced manufacturing and advanced materials design often require searches of relatively high-dimensional process control parameter spaces for settings that result in optimal structure, property, and performance parameters. The mapping from the former to the latter must be determined from noisy experiments or from expensive simulations. Here, we abstract this problem to a mathematical framework in which an unknown function from a control space to a design space must be ascertained by means of expensive noisy measurements, which locate control settings generating desired design features within specified tolerances, with quantified uncertainty. We describe targeted adaptive design (TAD), a new algorithm that performs this sampling task efficiently. TAD creates a Gaussian process surrogate model of the unknown mapping at each iterative stage, proposing a new batch of control settings to sample experimentally and optimizing the updated expected log-predictive probability density of the target design. TAD either stops upon locating a solution with uncertainties that fit inside the tolerance box or uses a measure of expected future information to determine that the search space has been exhausted with no solution. TAD thus embodies the exploration-exploitation tension in a manner that recalls, but is essentially different from, Bayesian optimization and optimal experimental design.

97 MATHEMATICS AND COMPUTING↗

Managing Wildfire Risk and Promoting Equity through Optimal Configuration of Networked Microgrids

As climate change increases the risk of large-scale wildfires, wildfire ignitions from electric power lines are a growing concern. To mitigate the wildfire ignition risk, many electric utilities de-energize power lines to prevent electric faults and failures. These preemptive power shutoffs are effective in reducing ignitions, but they could result in wide-scale power outages. Advanced technology, such as networked microgrids, can help reduce the size of the resulting power outages; however, even microgrid technology might not be sufficient to supply power to everyone, thus forcing hard questions about how to prioritize the provision of power among customers. In this paper, we present an optimization problem that configures networked microgrids to manage wildfire risk while maximizing the power served to customers; however, rather than simply maximizing the amount of power served in kilowatts, our formulation also considers the ability of customers to cope with power outages, as measured by social vulnerability, and it discourages the disconnection of particularly vulnerable customer groups. To test our model, we leverage a synthetic but realistic distribution feeder, along with publicly available social vulnerability indices and satellite-based wildfire risk map data, to quantify the parameters in our optimal decision-making model. Our case study results demonstrate the benefits of networked microgrids in limiting load shed and promoting equity during scenarios with high wildfire risk.

distribution systems↗

Design of a High-Pressure Fuel System for Use with Dimethyl Ether

The paper documents the modeling and experimental work on a common rail fuel injection system for Dimethyl Ether, a potential diesel substitute with a low carbon intensity signature. The DME fuel system is deployed on a light duty 2.2L compression ignition engine. The paper describes the injector optimization to shift to higher flows to account for the lower heating value and density of the DME when compared to diesel. The type of the injection system used for the DME application is an advanced rendering of the Common rail noted for a one-piece piston-needle injector construction and a solenoid driven spill valve featuring a pressure balanced poppet. A dedicated high-pressure fuel pump designed to pressurize DME is used. The design results in a fast acting open and close injection event, reduced leakage, with reduced cavitation in the fuel injector volume. Design parameters for system optimization included fill and spill orifices, needle lift, bias spring, and injector hole size. The design model provides good correlation of the instantaneous rates of injection with experiments across a wide range of pressure and injection timings. Proposed performance milestones for the design included similar DME injection duration to the diesel counterpart for same fuel energy injected into the cylinder to retain high engine cycle efficiency. Here, the dedicated DME design provided reduced hydraulic delays of 50%. Tests demonstrated sustained operation at pressures of 1000 bar, with capability to reach 1500bar. Durability tests showed no cavitation-deterioration over a 200-hour test cycle by means of spray imaging and hardware inspection.

De Ojeda, William [WM International Engineering, L↗

Shifting agricultural practices to produce sustainable, low carbon intensity feedstocks for biofuel production

The carbon intensity (CI) of biofuel's well-to-pump life cycle is calculated by life cycle analysis (LCA) to account for the energy/material inputs of the feedstock production and fuel conversion stages and the associated greenhouse gas (GHG) emissions during these stages. The LCA is used by the California Air Resources Board's Low Carbon Fuel Standard (LCFS) program to calculate CI and monetary credits are issued based on the difference between a given fuel's CI and a reference fuel's CI. Through the Tier 2 certification program under which individual fuel production facilities can submit their own CIs with their facility input data, the LCFS has driven innovative technologies to biofuel conversion facilities, resulting in substantial reductions in GHG emissions as compared to the baseline gasoline or diesel. A similar approach can be taken to allow feedstock petition in the LCFS so that lower-CI feedstock can be rewarded. Here we examined the potential for various agronomic practices to improve the GHG profiles of corn ethanol by performing feedstock-level CI analysis for the Midwestern United States. Our system boundary covers GHG emissions from the cradle-to-farm-gate activities (i.e. farm input manufacturing and feedstock production), along with the potential impacts of soil organic carbon change during feedstock production. We conducted scenario-based CI analysis of ethanol, coupled with regionalized inventory data, for various farming practices to manage corn fields, and identified key parameters affecting cradle-to-farm-gate GHG emissions. The results demonstrate large spatial variations in CI of ethanol due to farm input use and land management practices. In particular, adopting conservation tillage, reducing nitrogen fertilizer use, and implementing cover crops has the potential to reduce GHG emissions per unit corn produced when compared to a baseline scenario of corn–soybean rotation. This work shows a large potential emission offset opportunity by allowing feedstock producers a path to Tier 2 petitions that reward low-CI feedstocks and further reduce biofuels' CI. The prevalence of significant acreage that has not been optimized for CI suggests that policy changes that incentivize optimization of this parameter could provide significant additionality over current trends in farm efficiency and adoption of conservation practice.

09 BIOMASS FUELS↗

GP Cosmology Surrogate v1.0

GP Cosmology Surrogate is a Python library for building and training a generalized multi-output Gaussian process (GP) framework of @takhtaganov2021cosmic. In this approach, the surrogate is constructed sequentially, guided by a Bayesian optimization acquisition function that targets reduction of emulation error in the regions most consistent with the observational data. This adaptive design concentrates computational resources where they have the greatest impact on inference accuracy. The library supports efficient training for separable GP kernels, which allows the use of Kronecker algebra to handle high-dimensional input spaces and large numbers of correlated outputs. This makes it well suited for applications such as modeling cosmological power spectra, large-scale physical simulations, and multi-output hyperparameter tuning. By combining scalable multi-output GP modeling with data-driven adaptive sampling, GPsurrogate enables parameter inference and optimization with substantially fewer simulations than conventional space-filling designs.

Lukic, Zarija [Lawrence Berkeley National Laborato↗

Shortest path network interdiction with asymmetric uncertainty

Abstract This paper considers an extension of the shortest path network interdiction problem that incorporates robustness to account for parameter uncertainty. The shortest path interdiction problem is a game of two players with conflicting agendas and capabilities: an evader, who traverses the arcs of a network from a source node to a sink node using a path of shortest length, and an interdictor, who maximizes the length of the evader's shortest path by interdicting arcs on the network. It is usually assumed that the parameters defining the network are known exactly by both players. We consider the situation where the evader assumes the nominal parameter values while the interdictor uses robust optimization techniques to account for parameter uncertainty or sensor degradation. We formulate this problem as a nonlinear mixed‐integer semi‐infinite bilevel program and show that it can be converted into a mixed‐integer linear program with a second order cone constraint. We use random geometric networks and transportation networks to perform computational studies and demonstrate the unique decision strategies that our variant produces. Solving the shortest path interdiction problem with asymmetric uncertainty protects the interdictor from investing in a strategy that hinges on key interdictions performing as promised. It also provides an alternate strategy that mitigates the risk of these worst‐case possibilities.

Punla‐Green, She'ifa Z.↗

Energy efficiency in industrial drying: A hybrid ultrasonic system with a novel dynamic optimization framework

Drying processes are among the most energy-consuming operations in industrial and manufacturing settings, demanding strategic selection, design, and control for enhanced efficiency. Advancing drying technologies is critical for improving sustainability, lowering energy use, reducing carbon emissions, and minimizing waste. This study explores two innovative strategies aimed at transforming drying processes into sustainable, low-carbon systems by reducing energy consumption, minimizing waste, and maintaining a strong emphasis on preserving product quality. The first strategy showcases a sub-pilot scale hybrid ultrasonic-convective dryer for agrifood products. This technology, powered by electricity (process electrification), integrates non-thermal ultrasonic dehydration with convective heating and is presented as a sustainable and energy-efficient solution that enhances eco-friendly practices. The second strategy involves introducing and implementing a novel, multiobjective, mixed integer dynamic optimization technique to determine the optimal time-dependent process parameter values for the drying operation. This optimization technique yields operating conditions that are piecewise constant in time aiming to maximize the energy efficiency of the hybrid ultrasonic-convective dryer while ensuring strict adherence to product quality constraints. By adopting the hybrid ultrasonic-convective dryer, a notable 35% improvement in energy efficiency was achieved compared to conventional hot-air drying systems for drying apple slices. The proposed optimization framework further enhanced energy efficiency by nearly 14% over the most efficient process on the identical testbed, under static operating conditions. The reported enhancements have been experimentally validated. Regarding drying time (thereby improving production yield), the developed hybrid ultrasonic-convective dryer demonstrates as much as a 41% reduction in total processing time, which is further optimized by an additional 10% using our proposed optimization framework. The research outcomes have profound implications for the design and operation of drying systems, encompassing crucial aspects such as process electrification, cost-effectiveness, energy savings, time efficiency, product yield, product quality, and process automation.

Dynamic optimization↗

Improving the W Coating Uniformity by a COMSOL Model-Based CVD Parameter Study for Denser W f /W Composites

Tungsten (W) has the unique combination of excellent thermal properties, low sputter yield, low hydrogen retention, and acceptable activation. Therefore, W is presently the main candidate for the first wall and armor material for future fusion devices. However, its intrinsic brittleness and its embrittlement during operation bears the risk of a sudden and catastrophic component failure. As a countermeasure, tungsten fiber-reinforced tungsten (W f /W) composites exhibiting extrinsic toughening are being developed. A possible W f /W production route is chemical vapor deposition (CVD) by reducing WF 6 with H 2 on heated W fabrics. The challenge here is that the growing CVD-W can seal gaseous domains leading to strength reducing pores. In previous work, CVD models for W f /W synthesis were developed with COMSOL Multiphysics and validated experimentally. In the present article, these models were applied to conduct a parameter study to optimize the coating uniformity, the relative density, the WF 6 demand, and the process time. A low temperature and a low total pressure increase the process time, but in return lead to very uniform W layers at the micro and macro scales and thus to an optimized relative density of the W f /W composite. High H 2 and low WF 6 gas flow rates lead to a slightly shorter process time and an improved coating uniformity as long as WF 6 is not depleted, which can be avoided by applying the presented reactor model.

36 MATERIALS SCIENCE↗

Separable physics-informed DeepONet: Breaking the curse of dimensionality in physics-informed machine learning

The deep operator network (DeepONet) has shown remarkable potential in solving partial differential equations (PDEs) by mapping between infinite-dimensional function spaces using labeled datasets. However, in scenarios lacking labeled data, the physics-informed DeepONet (PI-DeepONet) approach, which utilizes the residual loss of the governing PDE to optimize the network parameters, faces significant computational challenges, particularly due to the curse of dimensionality. This limitation has hindered its application to high-dimensional problems, making even standard 3D spatial with 1D temporal problems computationally prohibitive. Additionally, the computational requirement increases exponentially with the discretization density of the domain. Here, to address these challenges and enhance scalability for high-dimensional PDEs, we introduce the Separable physics-informed DeepONet (Sep-PI-DeepONet). This framework employs a factorization technique, utilizing sub-networks for individual one-dimensional coordinates, thereby reducing the number of forward passes and the size of the Jacobian matrix required for gradient computations. By incorporating forward-mode automatic differentiation (AD), we further optimize computational efficiency, achieving linear scaling of computational cost with discretization density and dimensionality, making our approach highly suitable for high-dimensional PDEs. We demonstrate the effectiveness of Sep-PI-DeepONet through three benchmark PDE models: the viscous Burgers’ equation, Biot’s consolidation theory, and a parameterized heat equation. Our framework maintains accuracy comparable to the conventional PI-DeepONet while reducing training time by two orders of magnitude. Notably, for the heat equation solved as a 4D problem, the conventional PI-DeepONet was computationally infeasible (estimated 289.35 h), while the Sep-PI-DeepONet completed training in just 2.5 h. These results underscore the potential of Sep-PI-DeepONet in efficiently solving complex, high-dimensional PDEs, marking a significant advancement in physics-informed machine learning.

Neural operator↗

Experimental study of ECH pre-ionization on J-TEXT

An experimental study on electron cyclotron heating (ECH) pre-ionization has been conducted on J-TEXT in support of the joint experiment research for ITER plasma initiation. In this experiment, ECH power was injected to the vessel before the application of loop voltage ionize the neutral gas and form the initial plasma or so-called pre-plasma. The impact of several significant factors, such as magnetic field configuration, pre-fill gas pressure, ECH toroidal injection angle and ECH power on the evolution of pre-plasma are systematically studied, aiming to identify shared features, clarify their potential relationship and optimize the discharge parameters to generate a rather high pre-plasma density. By separating ECH power from the inductive start-up, the effect of pre-plasma on tokamak start-up can be observed. A dynamic magnetic configuration facilitates the transition of pre-plasma to tokamak plasma. To assess the influence of pre-plasma density on tokamak start-up, two kinds of magnetic field configurations are examined. While the effect of different pre-plasma densities on tokamak start-up is negligible, a significant difference is observed between pure ohmic start-up and start-up with pre-ionization. The studies presented here show evolutionary trends and threshold values needed to optimize ECH pre-ionization and also a feasible way to improve pre-plasma density and a configuration to stabilize pre-plasma for transition. Eventually, these results may contribute to multi-machine research and physics analysis that can assist ITER with its optimal preparation for first plasma operation.

ECH↗

Enabling Hyper-Differential Sensitivity Analysis for Ill-Posed Inverse Problems

Inverse problems constrained by partial differential equations (PDEs) play a critical role in model development and calibration. In many applications, there are multiple uncertain parameters in a model that must be estimated. However, high dimensionality of the parameters and computational complexity of the PDE solves make such problems challenging. A common approach is to reduce the dimension by fixing some parameters (which we will call auxiliary parameters) to a best estimate and use techniques from PDE-constrained optimization to estimate the other parameters. In this article, hyper-differential sensitivity analysis (HDSA) is used to assess the sensitivity of the solution of the PDE-constrained optimization problem to changes in the auxiliary parameters. Foundational assumptions for HDSA require satisfaction of the optimality conditions which are not always practically feasible as a result of ill-posedness in the inverse problem. Here we introduce novel theoretical and computational approaches to justify and enable HDSA for ill-posed inverse problems by projecting the sensitivities on likelihood informed subspaces and defining a posteriori updates. Our proposed framework is demonstrated on a nonlinear multiphysics inverse problem motivated by estimation of spatially heterogeneous material properties in the presence of spatially distributed parametric modeling uncertainties.

97 MATHEMATICS AND COMPUTING↗

Impact of n -Butylammonium Bromide on the Chemical and Electronic Structure of Double-Cation Perovskite Thin Films

2D/3D perovskite heterostructures have emerged as a promising material composition to reduce nonradiative recombination in perovskite-based LEDs and solar cells. Such heterostructures can be created by a surface treatment with large organic cations, for example, n-butylammonium bromide (BABr). To understand the impact of the BABr surface treatment on the double-cation (Cs 0.17 FA 0.83 Pb(I 0.6 Br 0.4 ) 3 ) (FA = formamidinium) perovskite thin film and further optimize the corresponding structures, an in-depth understanding of the chemical and electronic properties of the involved surfaces, interfaces, and bulk is required. We study the impact of the BABr treatment with a combination of surface-sensitive X-ray photoelectron spectroscopy and bulk-sensitive resonant inelastic soft X-ray scattering (RIXS). A quantitative analysis of the BABr-treated perovskite thin film shows a modified chemical perovskite surface environment of carbon, nitrogen, bromine, iodine, and lead, indicating that the treatment leads to a perovskite surface with a modified composition and bonding structure. With K-edge RIXS, the local environment at the nitrogen and carbon atoms is probed, allowing us to identify the presence of BABr in the perovskite bulk albeit with a modified bonding environment. This, in turn, identifies a "hidden parameter" for the optimization of the BABr treatment and overall performance of 2D/3D perovskite solar cell absorbers.

2D/3D perovskite heterostructure↗

Stable magnesium zinc oxide by reactive Co-Sputtering for CdTe-based solar cells

Magnesium zinc oxide (MZO) is a promising front contact material for CdTe solar cells. Due to its higher band gap than traditional CdS, MZO can reduce parasitic absorption to significantly increase short-circuit current density while also providing a benefit of conduction band offset tuning through Mg:Zn ratio optimization. MZO has been successfully implemented into CdTe devices, however its stability has been of concern. The MZO stability issue has been attributed to the presence of oxygen in the CdTe device processing ambient, leading to double-diode behavior (S-kink) in the current density-voltage curves. Here we report on MZO thin films deposited by reactive co-sputtering. The reactively co-sputtered MZO thin films have encouraging stability, show no significant variation in work function of the surface over a period of 6 months, as measured by Kelvin probe. Energy conversion efficiencies of around 16% have been achieved both with and without presence of oxygen in device processing ambients across multiple research facilities. Finally, these efficiencies should be possible to increase further by tuning of the thin film deposition and device processing parameters, especially through optimization of the back contact.

14 SOLAR ENERGY↗

Optimal Sizing of an Electric Vehicle Charging Station with Integration of PV and Energy Storage

This paper proposes an optimization model for the optimal configuration of an grid-connected electric vehicle (EV) extreme fast charging station considering integration of photovoltaic (PV) and energy storage. The proposed model minimizes the annualized net cost (i.e., maximizes the annualized net profit) of the extreme fast charging station, including investment and maintenance cost of charging ports, PV and energy storage, net cost of purchasing energy from utility and selling energy to EV customers, degradation cost of energy storage and demand charge. The decision variables are number of charging ports, capacity of invested PV and the power and energy ratings of invested energy storage. The Erlang-loss system is adopted to model the EV mobility. Results of numerical simulations indicate that investment of PV and energy storage could increase the annualized profit of the extreme fast charging station. In addition, the impacts of various parameters on the optimal solution are investigated by sensitivity analysis.

Liu, Guodong↗

Beyond Single-Reference Fixed-Node Approximation in Ab Initio Diffusion Monte Carlo Using Antisymmetrized Geminal Power Applied to Systems with Hundreds of Electrons

Diffusion Monte Carlo (DMC) is an exact technique to project out the ground state (GS) of a Hamiltonian. Since the GS is always bosonic, in Fermionic systems, the projection needs to be carried out while imposing antisymmetric constraints, which is a nondeterministic polynomial hard problem. In practice, therefore, the application of DMC on electronic structure problems is made by employing the fixed-node (FN) approximation, consisting of performing DMC with the constraint of having a fixed, predefined nodal surface. How do we get the nodal surface? The typical approach, applied in systems having up to hundreds or even thousands of electrons, is to obtain the nodal surface from a preliminary mean-field approach (typically, a density functional theory calculation) used to obtain a single Slater determinant. This is known as single reference. In this paper, we propose a new approach, applicable to systems as large as the C 60 fullerene, which improves the nodes by going beyond the single reference. In practice, we employ an implicitly multireference ansatz (antisymmetrized geminal power wave function constraint with molecular orbitals), initialized on the preliminary mean-field approach, which is relaxed by optimizing a few parameters of the wave function determining the nodal surface by minimizing the FN-DMC energy. We highlight the improvements of the proposed approach over the standard single-reference method on several examples and, where feasible, the computational gain over the standard multireference ansatz, which makes the methods applicable to large systems. We also show that physical properties relying on relative energies, such as binding energies, are affordable and reliable within the proposed scheme.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

From atomistic models to machine learning: Predictive design of nanocarbons under extreme conditions

The formation of technologically valuable nanocarbon structures under extreme conditions, such as those produced during high-explosive detonations, remains poorly understood but holds significant potential for the development of controlled synthesis pathways. While detonation shockwaves provide the high-pressure, high-temperature environment required for nanodiamond formation, subsequent cooling and decompression dictate whether the diamond phase is preserved or transformed into other nanocarbon structures. Here, in this study, we employ GPU-accelerated reactive molecular dynamics (ReaxFF) simulations to investigate the graphitization and structural remodeling of detonation nanodiamond under nonlinear quench and pressure-release trajectories. We further investigate how the initial nanodiamond morphology; cuboctahedral, octahedral, or hexagonal prism influences the resulting transformation products. Evolution of nanostructure, allotrope (via simulated x-ray diffraction), carbon hybridization, and ring statistics are tracked during a two-stage quench from 5000 K to 60 GPa. Rapid cooling combined with slow decompression optimizes cubic diamond retention, whereas slow cooling with rapid pressure release promotes surface-to-core graphitization, producing concentric sp 2 -hybridized layers and hollowed inner shells. Octahedral nanodiamonds evolve into carbon nano-onions, initially forming bucky diamonds that progressively transform into fully sp 2 -hybridized structures, while hexagonal prisms preferentially form parallel-stacked graphite layers resembling carbon dots. Transient hexagonal diamond (lonsdaleite) emerges as an interfacial phase, suggesting potential reversibility in the shock-induced graphite-to-diamond transformation pathway transformation route. To extend predictive capabilities, we trained machine learning (ML) regressors on over 10 5 node-hours of molecular dynamics (MD) trajectories. A multilayer perceptron (MLP) model reliably predicts the number of graphitized layers from temperature–pressure trajectories with a coefficient of determination (R 2 ) exceeding 0.90. This high predictive fidelity enables efficient, high-throughput mapping of the synthesis parameter space for optimized graphitization outcomes. Collectively, morphological control combined with optimized quench–decompression conditions promote the selective synthesis of nanocarbon allotropes. This work establishes a data-driven framework for the rational, a priori design of carbon nanomaterials for applications in energy storage, sensing, and biomedicine.

Detonation nanodiamond remodeling↗

Analytical characterization of laser induced plasmas towards uranium isotopic analysis in gaseous uranium hexafluoride

To perform direct enrichment assay on gaseous uranium hexafluoride (UF 6 ) with laser induced breakdown spectroscopy (LIBS), the dominant spectral-line features, evolution of the signal and background of the U II 424.437 nm line, and its Stark width and shift, were studied as a function of UF 6 gas pressure and pulse energy of a nanosecond Nd:YAG laser. In this work, vapor pressure of UF 6 was found to be the most important parameter for LIBS analysis of gaseous UF 6 . Spectral congestion with numerous U lines of high excitation potential was observed and signal-to-background ratio (SBR) was low for measurements with 80 Torr UF 6 . Only when both UF 6 vapor pressure and laser pulse energy were low, for example, less than 20 Torr pressure and 30 mJ pulse energy, the resultant LIBS spectra from gaseous UF 6 resembled those obtained from solid U samples. The experimental data also suggest that U and F atoms recombine back to UF 6 after the laser pulse. The U emission was found to decay fast with a persistent background signal, degrading SBR with delay time. Systematic positive biases were found for UF 6 enrichment assays performed with the 235 U– 238 U line pair at 424.412–424.437 nm, which was confirmed to be caused by self-absorption. Even with optimization of experimental parameters and incorporation of a self-absorption term into the spectral-fitting algorithm, to reduce and compensate for self-absorption, self-absorption is still a main factor limiting accurate UF 6 enrichment assay. The use of another spectral window which contains no resonance lines is a prospective solution for the self-absorption issue.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗