Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “field 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 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

ARC-DPA 1-MeV Neutron Fluence Equivalence Model for GaAs

Here, the legacy data used in the production of the American Society for Testing and Materials (ASTM) E722 standard for gallium arsenide (GaAs) were fit with the athermal recombination corrected-displacements per atom (ARC-DPA) model for the development of a novel 1-MeV neutron equivalent fluence metric. Improving on previous work, the coefficients for the ARC-DPA model were optimized using a conjugate gradient method, and the uncertainties in the observed damage in different benchmark neutron fields were sampled to quantify their contribution to the damage metric. The optimized parameters of −0.474 ± 0.03 for b arcdpa and 0.016 ± 0.002 for c arcdpa provided good agreement and indicated that the solution is not sensitive to the corresponding fluence uncertainties. Compared to the current 1-MeV neutron fluence equivalence standard for displacement damage in GaAs, our results indicate that for a thermal reactor neutron environment, displacement damage has been underestimated by around 25% and that the saturation value for the damage efficiency had been higher than previously modeled.

1-MeV equivalent fluence↗
DOE OSTI2021

Optimal design of acoustic metamaterial cloaks under uncertainty

In this work, we consider the problem of optimal design of an acoustic cloak under uncertainty and develop scalable approximation and optimization methods to solve this problem. The design variable is taken as an infinite-dimensional spatially-varying field that represents the material property, while an additive infinite-dimensional random field represents the variability of the material property or the manufacturing error. Discretization of this optimal design problem results in high-dimensional design variables and uncertain parameters. To solve this problem, we develop a computational approach based on a Taylor approximation and an approximate Newton method for optimization, which is based on a Hessian derived at the mean of the random field. We show our approach is scalable with respect to the dimension of both the design variables and uncertain parameters, in the sense that the necessary number of acoustic wave propagations is essentially independent of these dimensions, for numerical experiments with up to one million design variables and half a million uncertain parameters. Additionally, we demonstrate that, using our computational approach, an optimal design of the acoustic cloak that is robust to material uncertainty is achieved in a tractable manner. The optimal design under uncertainty problem is posed and solved for the classical circular obstacle surrounded by a ring-shaped cloaking region, subjected to both a single-direction single-frequency incident wave and multiple-direction multiple-frequency incident waves. Finally, we apply the method to a deterministic large-scale optimal cloaking problem with complex geometry, to demonstrate that the approximate Newton method’s Hessian computation is viable for large, complex problems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

MyCrunchGPT: A LLM Assisted Framework for Scientific Machine Learning

Scientific machine learning (SciML) has advanced recently across many different areas in computational science and engineering. Here, the objective is to integrate data and physics seamlessly without the need of employing elaborate and computationally taxing data assimilation schemes. However, preprocessing, problem formulation, code generation, postprocessing, and analysis are still time- consuming and may prevent SciML from wide applicability in industrial applications and in digital twin frameworks. Here, we integrate the various stages of SciML under the umbrella of ChatGPT, to formulate MyCrunchGPT, which plays the role of a conductor orchestrating the entire workflow of SciML based on simple prompts by the user. Specifically, we present two examples that demonstrate the potential use of MyCrunchGPT in optimizing airfoils in aerodynamics, and in obtaining flow fields in various geometries in interactive mode, with emphasis on the validation stage. To demonstrate the flow of the MyCrunchGPT, and create an infrastructure that can facilitate a broader vision, we built a web app based guided user interface, that includes options for a comprehensive summary report. The overall objective is to extend MyCrunchGPT to handle diverse problems in computational mechanics, design, optimization and controls, and general scientific computing tasks involved in SciML, hence using it as a research assistant tool but also as an educational tool. While here the examples focus on fluid mechanics, future versions will target solid mechanics and materials science, geophysics, systems biology, and bioinformatics.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

An Adaptive Newton-Based Free-Boundary Grad–Shafranov Solver

Equilibria in magnetic confinement devices result from force balancing between the Lorentz force and the plasma pressure gradient. In an axisymmetric configuration like a tokamak, such an equilibrium is described by an elliptic equation for the poloidal magnetic flux, commonly known as the Grad–Shafranov equation. It is challenging to develop a scalable and accurate free-boundary Grad–Shafranov solver, since it is a fully nonlinear optimization problem that simultaneously solves for the magnetic field coil current outside the plasma to control the plasma shape. In this work, we develop a Newton-based free-boundary Grad–Shafranov solver using adaptive finite elements and preconditioning strategies. The free-boundary interaction leads to the evaluation of a domain-dependent nonlinear form of which its contribution to the Jacobian matrix is achieved through shape calculus. The optimization problem aims to minimize the distance between the plasma boundary and specified control points while satisfying two nontrivial constraints, which correspond to the nonlinear finite element discretization of the Grad–Shafranov equation and a constraint on the total plasma current involving a nonlocal coupling term. The linear system is solved by a block factorization, and AMG is called for subblock elliptic operators. The unique contributions of this work include the treatment of a global constraint, preconditioning strategies, nonlocal reformulation, and the implementation of adaptive finite elements. Furthermore, it is found that the resulting Newton solver is robust, successfully reducing the nonlinear residual to 1e-6 and lower in a small handful of iterations while addressing the challenging case to find a Taylor state equilibrium where conventional Picard-based solvers fail to converge.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

An analytical approach to evaluating magnetic-field closure and topological changes in FRC devices

Here, we describe mathematical methods based on optimizing a modified non-linear flux function (MFF) to evaluate whether odd-parity perturbations affect the local closure of magnetic field lines in field-reversed configurations. Using the MFF methodology, quantitative formulas are derived that provide the shift of the field minimum and the threshold for field-line opening, a discontinuous change in field topology.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

High-speed characterization of two-phase flow and bubble dynamics in titanium felt porous media for hydrogen production

Bubble dynamics and two-phase flow phenomena are closely related to the performance of proton exchange membrane electrolyzer cells (PEMECs). This paper reports an in-situ study of the oxygen bubble behavior and associated multiphase evolutions in the anode side of PEMECs with titanium (Ti) felt liquid gas diffusion layers (LGDLs) by a high-speed visualization system. The micro oxygen bubble dynamics was captured and analyzed at different locations and virous operating conditions. The results show that the bubble detachment frequency and detachment diameter greatly increase with the operating current density. Additionally, they are significantly impacted by the local pore structure and morphology of Ti felt LGDLs. In the flow channels, there exist only several discrete micro bubbles at a low current density (0.04 A/cm 2 ) and a large flow velocity (133 mm/s). At a current density (0.2 A/cm 2 ) and a flow velocity (67 mm/s), a number of gas slugs are formed in the follow channels, in addition to discrete micro bubbles. At a high current density (1 A/cm 2 ) and a flow velocity (67 mm/s), more bubbles appear in the channel, and the flow field is dominated by slug or annular flows. These investigations can help to better understand the two-phase flow and bubble detachment mechanism, and provide a foundation for electrochemical reaction, multiphase flow studies and optimize the design of gas diffusion layers and flow fields for PEMECs in the future.

08 HYDROGEN↗
DOE OSTI2020

Three-dimensional traction microscopy accounting for cell-induced matrix degradation

Tractions exerted by cells on the extracellular matrix (ECM) are critical in many important physiological and pathological processes such as embryonic morphogenesis, wound healing, and cancer metastasis. Three-dimensional Traction Microscopy (3DTM) is a tool to quantify cellular tractions by first measuring the displacement field in the ECM in response to these tractions, and then using this measurement to infer tractions. Most applications of 3DTM have assumed that the ECM has spatially-uniform mechanical properties, but cells secrete enzymes that can locally degrade the ECM. In this work, a novel computational method is developed to quantify both cellular tractions and ECM degradation. In particular, the ECM is modeled as a hyperelastic, Neo-Hookean solid, whose material parameters are corrupted by a single degradation parameter. The feasibility of determining both the traction and the degradation parameter is first demonstrated by showing the existence and uniqueness of the solution. An inverse problem is then formulated to determine the nodal values of the traction vector and the degradation parameter, with the objective of minimizing the difference between a predicted and measured displacement field, under the constraint that the predicted displacement field satisfies the equation of equilibrium. The inverse problem is solved by means of a gradient-based optimization approach, and the gradient is computed efficiently using appropriately derived adjoint fields. The computational method is validated in-silico using a geometrically realistic neuronal cell model and synthetic traction and degradation fields. It is found that the method accurately recovers both the traction and degradation fields. Moreover, it is found that neglecting ECM degradation can yield significant errors in traction measurements. Finally, our method can extend the range of context where tractions can be appropriately measured.

42 ENGINEERING↗
DOE OSTI2025

Optimizing Material Selection and Operational Conditions for XHV Systems: Lessons from AISI 1020 and 316L Comparative Studies

In this study, AISI 1020 low-carbon steel is investigated as a cost-effective alternative to SS316 stainless steel for reaching extreme high vacuum (XHV) conditions. After being baked at 400°C, a vacuum chamber made of the low-carbon steel material exhibited an outgassing rate approximately 2000 times smaller than a similar chamber made of stainless steel. Its activation energy for hydrogen diffusion (27 kJ/mol) is less than half that of stainless steel (60.3 kJ/mol), indicating more efficient hydrogen removal during bakeout. MolFlow+ simulations supported the experimental data and demonstrated the importance of system geometry optimization and minimizing stainless steel content for achieving optimal vacuum performance. AISI 1020's magnetic properties, typically considered disadvantageous for accelerator applications, could benefit spin-polarized electron sources by shielding photocathodes from stray fields while simultaneously providing improved vacuum through reduced outgassing. To optimize AISI 1020's performance in XHV systems, practical considerations include pre-baking protocols and careful system design to minimize stainless steel components.

Al-Allaq, Aiman H.↗
DOE OSTI2022

Optimization of 3D controlled ELM-free state with recovered global confinement for KSTAR with n = 1 resonant magnetic field perturbation

Mitigation of deleterious heat flux from edge-localized modes (ELMs) on fusion reactors is often attempted with 3D perturbations of the confining magnetic fields. However, the established technique of resonant magnetic perturbations (RMPs) also degrades plasma performance, complicating implementation on future fusion reactors. Here, we introduce an adaptive real-time control scheme on the KSTAR tokamak as a viable approach to achieve an ELM-free state and simultaneously recover high-confinement (βN ~ 1.91, β p ~ 1.53, and H 98 ~ 0.9), demonstrating successful handling of a volatile complex system through adaptive measures. We show that, by exploiting a salient hysteresis process to adaptively minimize the RMP strength, stable ELM suppression can be achieved while actively encouraging confinement recovery. This is made possible by a self-organized transport response in the plasma edge which reinforces the confinement improvement through a widening of the ion temperature pedestal and promotes control stability, in contrast to the deteriorating effect on performance observed in standard RMP experiments. These results establish the real-time approach as an up-and-coming solution toward an optimized ELM-free state, which is an important step for the operation of ITER and reactor-grade tokamak plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

In-Depth Optimization with the OpenACC-to-FPGA Framework on an Arria 10 FPGA

The reconfigurable computing paradigm that uses field programmable gate arrays (FPGAs) has received renewed interest in the high-performance computing field due to FPGAs' unique combination of performance and energy efficiency. However, difficulties in programming and optimizing FPGAs have prevented them from being widely accepted as general-purpose computing devices. In accelerator-based heterogeneous computing, portability across diverse heterogeneous devices is also an important issue, but the unique architectural features in FPGAs make this difficult to achieve. To address these issues, a directive-based, high-level FPGA programming and optimization framework was previously developed. In this work, developed optimizations were combined holistically using the directive-based approach to show that each individual benchmark requires a unique set of optimizations to maximize performance. The relationships between FPGA resource usages and runtime performance were also explored.

Lambert, Jacob↗
DOE OSTI2025

Kβ X-ray Emission Spectra Analysis Using Bayesian Optimization

The Kβ X-ray emission spectrum of 3 d transition metals is rich with electronic and structural information due to strong exchange interactions with the valence shell of the metal, and has become crucial for understanding their spin and oxidation states. The spectrum is commonly treated using crystal-field multiplet theory, a semi-empirical theory that uses tunable parameters to control the strength of the effects present in X-ray emission spectroscopy (XES). However, determining the experimental values of these parameters remains a challenge. We present a methodology that applies Bayesian optimization to crystal-field multiplet theory to determine parameter values. The algorithm is tested on the X-ray emission spectra of a collection of Mn, Co, and Ni oxides. We are able to find optimal values for the four most impactful parameters: Slater−Condon reduction factors F dd , F pd , and G pd , and crystal field splitting 10 Dq . The algorithm produces significantly improved accuracy compared to current analysis methods, and probes interparameter dependencies by modeling the error landscape. This advancement enhances XES analysis by offering an approach of obtaining quantitative electronic structural information on 3 d transition metal valence shells, facilitating applications across various scientific fields.

Bayesian optimization↗
DOE OSTI2023

Scale-up Testing of Advanced Polaris Membrane in CO2 Capture Technology

This final technical report describes work conducted by Membrane Technology and Research, Inc. (MTR) for the U.S. Department of Energy (DOE), National Energy Technology Lab (NETL) on the scale-up and testing of advanced Polaris™ membrane CO2 capture technology at the Technology Centre Mongstad (TCM) under award number DE-FE0031591. The work was performed from August 1, 2018 through January 31, 2023. The overall goal of this project was to design, build and operate an advanced Polaris membrane CO2 capture system at TCM. MTR was assisted in this project by Trimeric Corporation (Trimeric), an engineering design services company, the Carbon Capture Simulation for Industry Impact (CCSI2), a partnership among national laboratories, industry, and academic institutions, and the Technology Centre Mongstad (TCM), who provided the host site for the slipstream field test. This report details the work conducted to scale-up MTR’s second-generation (Gen-2) Polaris membrane and advanced planar membrane modules to a final form factor optimized for commercial use; validate their performance in an engineering-scale field test at TCM; and to show the potential of the MTR process to meet DOE CO2 capture targets from large source point emitters. Work for this project included membrane optimization and scale-up, advanced planar module design and fabrication, design and fabrication of an engineering-scale field test membrane skid, operation of the field test skid processing Residue Fluid Catalytic Cracker (RFCC) industrial flue gas at TCM, and a detailed techno-economic analysis (TEA) of the MTR membrane post-combustion process for CO2 capture. This project validated recent membrane technology advancements at the engineering-scale, moves the MTR advanced post-combustion capture technology to TRL-6, and mitigates risk in future Large Pilot or Demonstration scale-up activities.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2024

The influence of laser focusing conditions on the direct laser acceleration of electrons

Direct laser acceleration of electrons during a high-energy, picosecond laser interaction with an underdense plasma has been demonstrated to be substantially enhanced by controlling the laser focusing geometry. Experiments using the OMEGA EP facility measured electrons accelerated to maximum energies exceeding 120 times the ponderomotive energy under certain laser focusing, pulse energy, and plasma density conditions. Two-dimensional particle-in-cell simulations show that the laser focusing conditions alter the laser field evolution, channel fields generation, and electron oscillation, all of which contribute to the final electron energies. The optimal laser focusing condition occurs when the transverse oscillation amplitude of the accelerated electron in the channel fields matches the laser beam width, resulting in efficient energy gain. Through this observation, a simple model was developed to calculate the optimal laser focal spot size in more general conditions and is validated by experimental data.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Automatic Qubit Characterization and Gate Optimization with QubiC

As the size and complexity of a quantum computer increases, quantum bit (qubit) characterization and gate optimization become complex and time-consuming tasks. Current calibration techniques require complicated and verbose measurements to tune up qubits and gates, which cannot easily expand to the large-scale quantum systems. We develop a concise and automatic calibration protocol to characterize qubits and optimize gates using QubiC, which is an open source FPGA (field-programmable gate array) based control and measurement system for superconducting quantum information processors. We propose multi-dimensional loss-based optimization of single-qubit gates and full XY-plane measurement method for the two-qubit CNOT gate calibration. We demonstrate the QubiC automatic calibration protocols are capable of delivering high-fidelity gates on the state-of-the-art transmon-type processor operating at the Advanced Quantum Testbed at Lawrence Berkeley National Laboratory. Finally, the single-qubit and two-qubit Clifford gate infidelities measured by randomized benchmarking are of 4.9(1.1) × 10 -4 and 1.4(3) × 10 -2 , respectively.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Co-Simulation of Electric Power Distribution Systems and Buildings including Ultra-Fast HVAC Models and Optimal DER Control

Smart homes and virtual power plant (VPP) controls are growing fields of research with potential for improved electric power grid operation. A novel testbed for the co-simulation of electric power distribution systems and distributed energy resources (DERs) is employed to evaluate VPP scenarios and propose an optimization procedure. DERs of specific interest include behind-the-meter (BTM) solar photovoltaic (PV) systems as well as heating, ventilation, and air-conditioning (HVAC) systems. The simulation of HVAC systems is enabled by a machine learning procedure that produces ultra-fast models for electric power and indoor temperature of associated buildings that are up to 133 times faster than typical white-box implementations. Hundreds of these models, each with different properties, are randomly populated into a modified IEEE 123-bus test system to represent a typical U.S. community. Advanced VPP controls are developed based on the Consumer Technology Association (CTA) 2045 standard to leverage HVAC systems as generalized energy storage (GES) such that BTM solar PV is better utilized locally and occurrences of distribution system power peaks are reduced, while also maintaining occupant thermal comfort. An optimization is performed to determine the best control settings for targeted peak power and total daily energy increase minimization with example peak load reductions of 25+%.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2021

Control Design for Inhomogeneous-Broadening Compensation in Single-Photon Transducers

We report a transducer of single photons between microwave and optical frequencies can be used to realize quantum communication over optical fiber links between distant superconducting quantum computers. A promising scalable approach to constructing such a transducer is to use ensembles of quantum emitters interacting simultaneously with electromagnetic fields at optical and microwave frequencies. However, inhomogeneous broadening in the transition frequencies of the emitters can be detrimental to this collective action. In this paper, we utilize a gradient-based optimization strategy to design the temporal shape of the laser field driving the transduction system to mitigate the effects of inhomogeneous broadening. We study the improvement of transduction efficiencies as a function of inhomogeneous broadening in different single-emitter cooperativity regimes and correlate it with a restoration of super-radiance effects in the emitter ensembles. Furthermore, to assess the optimality of our pulse designs, we provide certifiable bounds on the design problem and compare them to the achieved performance.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

Optimizing the HSX stellarator for microinstability by coil-current adjustments

Abstract The optimization of helically symmetric experiment (HSX) for reduced microinstability has been achieved by examining a large set of configurations within a neighborhood of the standard operating configuration. This entailed generating a database of more than 10 6 magnetic-field configurations for HSX by varying the currents in external coils. Using a set of volume-averaged metrics and gyrokinetic simulations, this database has helped to identify a set of configurations that can be used to regulate trapped-electron-mode stability in HSX. This set of configurations is also found to correlate flux-surface elongation and triangularity with an increase in magnetic-well depth, an increase in rotational transform, and low neoclassical heat-flux relative to the standard quasi-helically-symmetric configuration. These results demonstrate sensitivity of plasma behavior in response to changes in a 3D magnetic field to both neoclassical and gyrokinetic models, and the experimental potential in HSX to explore turbulence optimization. This perturbative optimization approach is not unique to HSX, and can readily be deployed on existing fusion devices to identify novel magnetic-fields to be used in turbulence-optimization experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Algebraic Multigrid with Filtering: An Efficient Preconditioner for Interior Point Methods in Large-Scale Contact Mechanics Optimization

Large-scale contact mechanics simulations are crucial in many engineering fields such as structural design and manufacturing. In the frictionless case, contact can be modeled by minimizing an energy functional; however, these problems are often nonlinear, nonconvex, and increasingly difficult to solve as mesh resolution increases. In this work, we employ a Newton-based interior-point (IP) filter line-search method, an effective approach for large-scale constrained optimization. While this method converges rapidly, each iteration requires solving a large saddle-point linear system that becomes ill-conditioned as the optimization process converges, largely due to IP treatment of the contact constraints. Such ill-conditioning can hinder solver scalability and increase iteration counts with mesh refinement. Here, to address this, we introduce a novel preconditioner, algebraic multigrid with filtering (AMGF), tailored to the Schur complement of the saddle-point system. Building on the classical AMG solver, commonly used for elasticity, we augment it with a specialized subspace correction that filters near null space components introduced by contact interface constraints. Through theoretical analysis and numerical experiments on a range of linear and nonlinear contact problems, we demonstrate that the proposed solver achieves mesh independent convergence and maintains robustness against the ill-conditioning that notoriously plagues IP methods. These results indicate that AMGF makes contact mechanics simulations more tractable and broadens the applicability of Newton-based IP methods in challenging engineering scenarios. More broadly, AMGF is well suited for problems, optimization or otherwise, where solver performance is limited by a low-dimensional subspace, such as those arising from localized constraints, interface conditions, or model heterogeneities. This makes the method widely applicable beyond contact mechanics and constrained optimization.

Mathematics and Computing↗