Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Simulation 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 955 records · Page 53

Initial conceptual demonstration of control co-design for WEC optimization

Abstract While some engineering fields have benefited from systematic design optimization studies, wave energy converters have yet to successfully incorporate such analyses into practical engineering workflows. The current iterative approach to wave energy converter design leads to sub-optimal solutions. This short paper presents an open-source MATLAB toolbox for performing design optimization studies on wave energy converters where power take-off behavior and realistic constraints can be easily included. This tool incorporates an adaptable control co-design approach, in that a constrained optimal controller is used to simulate device dynamics and populate an arbitrary objective function of the user’s choosing. A brief explanation of the tool’s structure and underlying theory is presented. To demonstrate the capabilities of the tool, verify its functionality, and begin to explore some basic wave energy converter design relationships, three conceptual case studies are presented. In particular, the importance of considering (and constraining) the magnitudes of device motion and forces in design optimization is shown.

16 TIDAL AND WAVE POWER↗

Reduced energetic particle transport models enable comprehensive time-dependent tokamak simulations

Time-dependent integrated simulations through codes such as TRANSP are becoming an indispensable tool for the interpretation of existing experiments and predictions of optimized scenarios. For many practical cases, quantitative simulations need to include the effect of plasma instabilities on the evolution of a tokamak discharge. An example is the degradation in energetic particle (EP) confinement induced by instabilities, which in turn affects important source terms for heating, non-inductive current, and momentum in a simulation. The reduced-physics "kick model" provides phase-space resolved transport probability matrices to TRANSP that are used to account for enhanced EP transport by instabilities in addition to neoclassical transport. The model has recovered the measured Alfvén eigenmode (AE) spectrum on NSTX, NSTX-U and DIII-D, and has reproduced details of phase-space resolved fast ion diagnostic data measured on DIII-D for EP-driven modes and tearing modes. In general, the kick model has proven the potential of phase-space resolved EP simulations to unravel details of EP transport for detailed theory/experiment comparison and for scenario planning based on optimization of Neutral Beam (NB) injection parameters. In this work, the extension of the kick model to low-frequency instabilities such as Tearing Modes and fishbones, in addition to AEs, is assessed. The goal is to enable TRANSP simulations that retain the main effects of multiple types of instabilities through a common framework. Results from the NSTX/NSTX-U and DIII-D tokamaks show that the extension to "multi-mode" scenarios can expand the range of applicability of the model for more reliable, quantitative integrated simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

AC Resistance Reduction Using Orthogonal Air Gaps in High-Frequency Inductors: Preprint

This paper presents a relatively simple technique to reduce winding losses due to fringing fields in high-frequency inductors. In high-frequency power electronics, ac inductor winding losses are affected by skin and proximity effects, including uneven current distribution due to fringing magnetic fields around airgaps. It is well known how fringing effects can be mitigated using distributed airgaps, at the expense of non-standard core or winding geometry. In planar magnetics, the orthogonal-airgap approach proposed in this paper combines airgaps in core segments parallel with the windings with conventionally placed airgaps in segments perpendicular to the windings. The approach is developed using a 1D analytical framework and validated by 2D finite-element simulations. Analytical guidelines are presented to optimize the airgaps to achieve minimum ac resistance. As a case study, a planar inductor is designed for an 8~kW SiC-based buck converter operating at 250~kHz. It is shown how the orthogonal airgaps result in more than 45\% reduction in ac resistance and substantially reduced inductor losses compared to the design using standard airgaps. The results are verified by loss measurements on an experimental converter prototype.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Numerical investigation on a new designed hollow-core photonic crystal fiber with large modal separation

We investigate the characteristics of optical modes in a new designed hollow-core photonic crystal fiber (PCF) with large modal separation. The large effective refractive index difference between the adjacent second-order vector modes can be maintained higher than 0.001 within the wavelength range from 1 to 2 μm. The maximum value can arrive 0.01199 between HE{sub 21} and TM{sub 01} mode at 2 μm. The fundamental mode of the proposed PCF has flat dispersion in a wide spectral range, which can be applied on wideband supercontinuum generation. The influence on the characteristics by the PCF structure has also been investigated. The increase of the size of the central air hole and the duty cycle affects the flatness of the dispersion curve to some extent. Our simulated results could provide an optimized fiber structure for the experimental study on stable transmission of vector modes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Socorro Electronic Structure Software

Socorro can be used to compute the ground-state electron density for a periodically repeated simulation cell in which the external potential is obtained from norm-conserving pseudopotentials or projector-augmented-wave functions. The ground-state electron density is computed by summation over one-electron orbitals which are obtained using the Kohn-Sham formulation of density-functional theory and which are expanded in a plane-wave basis. Various quantities can be computed from the ground-state solution, including atomic forces which can be used to optimize atom positions in the simulation cell and to perform molecular dynamics simulations.

Foiles, StephenM↗

Socorro Electronic Structure Software

Socorro can be used to compute the ground-state electron density for a periodically repeated simulation cell in which the external potential is obtained from norm-conserving pseudopotentials or projector-augmented-wave functions. The ground-state electron density is computed by summation over one-electron orbitals which are obtained using the Kohn-Sham formulation of density-functional theory and which are expanded in a plane-wave basis. Various quantities can be computed from the ground-state solution, including atomic forces which can be used to optimize atom positions in the simulation cell and to perform molecular dynamics simulations.

Foiles, Stephen M↗

Design Paradigm for Modular Multilevel Converter-Based Generator Rectifier Systems

Modular Multilevel Converters (MMC) are being widely considered for medium to high voltage DC applications. Designing such converters through multi-objective optimization is of interest because such an approach allows the trade-off between competing objectives (for example mass and loss) to be explicitly and quantitatively identified. In this work, an optimization based design paradigm for MMC based generator rectifier systems is proposed. Such development typically requires detailed component design and simulation models for the electric machine and converter which are computationally expensive. As an alternative, the proposed work utilizes an electric machine metamodel, inductor metamodel, and high-speed steady-state simulation model for the MMC to facilitate multi-objective optimization minimizing system metrics of interest while satisfying system constraints. A case study is undertaken to demonstrate the functionality of the proposed design paradigm.

42 ENGINEERING↗

High Power Solid Target for Positron Source at CEBAF

The progress in the development of a polarized positron injector for the Continuous Electron Beam Accelerator Facility (CEBAF) at Jefferson Laboratory is presented. The polarized positron beam is generated by a high current polarized electron beam (>1 mA @ 120 MeV) via bremsstrahlung radiation and e+e- pair production in the tungsten target. The simulations show that using an optimized target and positron beamline, the positron injector can provide a cw positron beam with a current larger than 50 nA and a polarization as large as 60%. Injected into the North Linac of CEBAF at an energy of 123 MeV, the positron beam can reach a maximum energy of 12 GeV to perform a rich experimental program. The results of the thermal and structural FEA analysis of the heat load in the target are presented, as well as the simulation results of radiation damage in the target. The performed and planned target material fatigue and radiation damage tests are discussed. This work is supported by the U.S. DOE, Office of Science, Office of Nuclear Physics, Contract No. DE-AC05-06OR23177 and the European Union?s Horizon 2020 research and innovation program under grant agreement No. 824093.

Ushakov, Andriy↗

Mapping the performance of a versatile water-based condensation particle counter (vWCPC) with numerical simulation and experimental study

Accurate airborne aerosol instrumentation is required to determine the spatial distribution of ambient aerosol particles, particularly when dealing with the complex vertical profiles and horizontal variations of atmospheric aerosols. A versatile water-based condensation particle counter (vWCPC) has been developed to provide aerosol concentration measurements under various environments with the advantage of reducing the health and safety concerns associated with using butanol or other chemicals as the working fluid. However, the airborne deployment of vWCPCs is relatively limited due to the lack of characterization of vWCPC performance at reduced pressures. Given the complex combinations of operating parameters in vWCPCs, modeling studies have advantages in mapping vWCPC performance. In this work, we thoroughly investigated the performance of a laminar-flow vWCPC using COMSOL Multiphysics® simulation coupled with MATLAB™. We compared it against a modified vWCPC (vWCPC model 3789, TSI, Shoreview, MN, USA). Our simulation determined the performance of particle activation and droplet growth in the vWCPC growth tube, including the supersaturation, $D$ p,kel,0 (smallest size of particle that can be activated), $D$ p,kel,50 (particle size activated with 50 % efficiency) profile, and final growth particle size D d under wide operating temperatures, inlet pressures $P$ (30–101 kPa), and growth tube geometry (diameter $D$ and initiator length $L$ ini ). The effect of inlet pressure and conditioner temperature on vWCPC 3789 performance was also examined and compared with laboratory experiments. The COMSOL simulation result showed that increasing the temperature difference (Δ$T$) between conditioner temperature $T$ con and initiator $T$ ini will reduce $D$ p,kel,0 and the cut-off size $D$ p,kel,50 of the vWCPC. In addition, lowering the temperature midpoint ($T$ mid = $\frac{T_{con}+T_{ini}}{2}$ increases the supersaturation and slightly decreases the $D$ p,kel . The droplet size at the end of the growth tube is not significantly dependent on raising or lowering the temperature midpoint but significantly decreases at reduced inlet pressure, which indirectly alters the vWCPC empirical cut-off size. Our study shows that the current simulated growth tube geometry ($D$=6.3 mm and $L$ ini =30 mm) is an optimized choice for current vWCPC flow and temperature settings. The current simulation can more realistically represent the $D$ p,kel for 7 nm vWCPC and also achieved good agreement with the 2 nm setting. Using the new simulation approach, we provide an optimized operation setting for the 7 nm setting. This study will guide further vWCPC performance optimization for applications requiring precise particle detection and atmospheric aerosol monitoring.

47 OTHER INSTRUMENTATION↗

A novel high-temperature MAS probe with optimized temperature gradient across sample rotor for in-situ monitoring of high-temperature high-pressure chemical reactions

We present a novel nuclear magnetic resonance (NMR) probe design focused on optimizing the temperature gradient across the sample for high temperature magic angle spinning (MAS) experiments using standard rotors. Computational flow dynamics (CFD) simulations were used to assess and optimize the temperature gradient across the sample under MAS conditions. The chemical shift and linewidth of 207 Pb direct polarization in lead nitrate were used to calibrate the sample temperature and temperature gradient, respectively. A temperature gradient of less than 3 °C across the sample was obtained by heating bearing gas flows and adjusting its temperature and flow rate during variable temperature (VT) experiments. A maximum temperature of 350 °C was achieved in this probe using a Varian 5 mm MAS rotor with standard Vespel drive tips and end caps. Time-resolved 13 C and 1 H MAS NMR experiments were performed at 325 °C and 60 bar to monitor an in-situ mixed phase reverse water gas shift reaction, industrial synthesis of CH 3 OH from a mixture of CO 2 and H 2 with a Cu/ZnO/Al 2 O 3 catalyst, demonstrating the first in-situ NMR monitoring of a chemical system at temperatures higher than 250 °C in a pressurized environment. Lastly, the combination of this high-temperature probe and high-pressure rotors will allow for in-situ NMR studies of a great variety of chemical reactions that are inaccessible to conventional NMR setup.

47 OTHER INSTRUMENTATION↗

Modelica-based modeling and simulation of district cooling systems: A case study

While equation-based object-oriented modeling language Modelica can evaluate practical energy improvements for district cooling systems, few have adopted Modelica for this type of large-scale thermo-fluid system. Further, to our best knowledge, district cooling modeling studies have yet to include hydraulics in piping networks alongside plant models featuring realistic mechanical systems and controls. These are critical details to include when looking to make energy and control improvements in many physical system installations. To fill these gaps, this study released new open-source district cooling models at the Modelica Buildings Library and leveraged these models for a real-world case study at the University of Colorado Boulder. Here, the site includes six buildings connected to a central chiller plant featuring a waterside economizer. Several energy saving strategies are pursued based on the validated model, including control setpoint optimization, equipment modification, and pump setpoint adjustments. Results indicate that a combination of the studied measures can save the campus annually 84.6 MWh of energy, 8.9% of electricity costs, 58.0 metric tons of carbon dioxide emissions, while the waterside economizer cuts down chillers’ run times by 201 days/year, reducing maintenance costs and extending chiller life.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

TRANSP-based optimization towards tokamak scenario development

An optimization approach that incorporates the predictive transport code TRANSP is proposed for tokamak scenario development. Optimization methods are often employed to develop open-loop control strategies to aid access to high performance tokamak scenarios. Generally, the optimization approaches use control-oriented models, i.e. models that are reduced in complexity and prediction accuracy as compared to physics-oriented transport codes such as TRANSP. In the presented approach, an optimization procedure using the TRANSP code to simulate the tokamak plasma is considered for improved predictive capabilities. As a test case, the neutral beam injection (NBI) power is optimized to develop a control strategy that maximizes the noninductive current fraction during the ramp-up phase for NSTX-U. Simulation studies towards the achievement of noninductive ramp-up in NSTX-U have already been carried out with the TRANSP code. The optimization-based approach presented in this work is used to maximize the noninductive current fraction during ramp-up in NSTX-U, demonstrating that the scenario development task can be automated. An additional test case considers optimization of the current ramp rate in DIII-D for obtaining a stationary plasma characterized by a flat loop voltage profile in the flattop phase.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Simulation of Helium Flow Visualization Apparatus for Studies of Blanket Cooling in Fusion Reactors

Flow visualization is essential to understanding helium cooling performance. This article investigates helium flow visualization in an apparatus to support the design of the blanket first wall for a fusion reactor. Helium’s safety advantages make it an attractive coolant, but effective cooling and flow visualization remain challenging. Sophisticated simulations are conducted in different test sections to address three key areas in the design of a test apparatus and cooling channel enhancements: the influence of viewing glass windows on flow patterns, the impact of increased heating, and the effectiveness of cooling structures like baffles. Viewing glass windows introduce flow recirculation region bifurcation and flow asymmetry, affecting flow patterns and necessitating careful analysis in future experiments. Increased heating results in helium flow detachment from the heated surface, leading to hot spots. Cooling structures, particularly baffles, prove effective in maintaining consistent attachment to the heated surface, improving the heat transfer performance. Further, this study also examines a tunable parameter in the turbulence model, highlighting the importance of accurate model tuning for future fusion reactor cooling designs. With an imminent helium flow visualization facility, these simulations will be used to optimize cooling structures on the heated wall to improve flow attachment and heat transfer efficiency. This work serves as a first investigation of the helium flow visualization apparatus for blanket cooling enhancement in fusion reactor design.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Parallelized hybrid optimization methods for nonsmooth problems using NOMAD and linesearch

Two parallelized hybrid methods are presented for single-function optimization problems with side constraints. The optimization problems are difficult not only due to possible existence of local minima and nonsmoothness of functions, but also due to the fact that objective function and constraint values for a solution vector can only be obtained by querying a black box whose execution requires considerable computational effort. Examples are optimization problems in Engineering where objective function and constraint values are computed via complex simulation programs, and where local minima exist and smoothness of functions is not assured. The hybrid methods consist of the well-known method NOMAD and two new methods called DENCON and DENPAR that are based on the linesearch scheme CS-DFN. The hybrid methods compute for each query a set of solution vectors that are evaluated in parallel. The hybrid methods have been tested on a set of difficult optimization problems produced by a certain seeding scheme for multiobjective optimization. We compare computational results with solution by NOMAD, DENCON, and DENPAR as stand-alone methods. It turns out that among the stand-alone methods, NOMAD is significantly better than DENCON and DENPAR. However, the hybrid methods are definitely better than NOMAD.

97 MATHEMATICS AND COMPUTING↗

TRANSP-based Optimization Towards Tokamak Scenario Development

An optimization approach that incorporates the predictive transport code TRANSP is proposed for tokamak scenario development. Optimization methods are often employed to develop open-loop control strategies to aid access to high performance tokamak scenarios. In general, the optimization approaches use control-oriented models, i.e. models that are reduced in complexity and prediction accuracy as compared to physics-oriented transport codes such as TRANSP. In the presented approach, an optimization procedure using the TRANSP code to simulate the tokamak plasma is considered for improved predictive capabilities. As a test case, the neutral beam injection (NBI) power is optimized to develop a control strategy that maximizes the non-inductive current fraction during the ramp-up phase for NSTX-U. Simulation studies towards the achievement of non-inductive ramp up in NSTX-U have already been carried out with the TRANSP code. The optimization-based approach proposed in this work is used to maximize the non-inductive current fraction during ramp-up in NSTX-U, demonstrating that the scenario development task can be automated. An additional test case considers optimization of the current ramp rate in DIII-D for obtaining a stationary plasma characterized bya flat loop voltage profile in the flattop phase.

DIII-D↗

Enabling Efficient Surveillance, Control, and Automation of Geothermal Operations with Advanced Predictive Analytics

Automation and control of geothermal energy production and operations require reliable and efficient predictive tools. While physics-based simulation offers a comprehensive tool for predicting energy production performance in geothermal systems, predicting the behavior of geothermal reservoirs involves complex multi-physics processes with coupling effects, highly uncertain input parameters and subsurface descriptions. Moreover, building, running, and integrating simulation models into standard model calibration and optimization workflows entail significant technical and computational efforts. An emerging alternative to physics-based simulation is data-driven predictive analytics models that have gained popularity in energy industry. In this report, we develop novel predictive models for integration into real-time fault diagnosis and model predictive control algorithms to improve the efficiency of energy production operations in geothermal reservoirs. The report includes two major research Thrust Areas, that is, the surface power plant and the subsurface reservoir.

15 GEOTHERMAL ENERGY↗

Enabling kilometer-scale E3SM land model simulation over North America: A new integrated framework solution

This study introduces a novel framework designed to enhance the performance, scalability, and portability of the kilometer-scale E3SM Land Model (km-ELM) within the E3SM modeling infrastructure. By seamlessly integrating cutting-edge data tools, we address existing challenges such as slow performance, limited scalability, and difficulties in software integration in current data-driven ELM simulation over large geographic areas. Our innovative approach leverages the KiloCraft data toolkit to generate unified inputs for simulations ranging from a single-cite case, to a 72,083-cell regional case to a continental configuration encompassing 21.6 million land grid cells at a 1 km × 1 km resolution. We conduct extensive strong- and weak-scaling experiments on three state-of-the-art supercomputers, utilizing up to 100,800 CPU cores across 2400 compute nodes to evaluate end-to-end metrics including wall-clock time, simulation-years-per-day (SYPD), initialization costs, and I/O throughput. Our results reveal the land (LND) component’s efficient scaling, demonstrating near-ideal weak scaling and strong-scaling parallel efficiencies reaching up to 87% at 50,400 cores. We confirm portability and reproducibility through bitwise-equivalent outputs across different machines using identical inputs over supported machines. Notably, at extreme scales, we identify I/O as a critical bottleneck and that leads to effective solution with the SCORPIO/ADIOS stack. Collectively, these findings validate the deployment of km-ELM at a continental scale with high parallel efficiency and provide essential guidance on configuration, decomposition, and I/O settings for optimized kilometer-scale land simulations in E3SM. This work emphasizes the innovative design and practical solutions that enhance the operational capabilities of km-ELM, focusing on software performance and scalability while leaving detailed scientific evaluations of simulated land processes for future investigations.

E3SM land model (ELM), km-ELM, scalability, perfor↗

Real-time hybrid controls of energy storage and load shedding for integrated power and energy systems of ships

This paper presents an original energy management methodology to enhance the resilience of ship power systems. The integration of various energy storage systems (ESS), including battery energy storage systems (BESS) and super-capacitor energy storage systems (SCESS), in modern ship power systems poses challenges in designing an efficient energy management system (EMS). The EMS proposed in this paper aims to achieve multiple objectives. The primary objective is to minimize shed loads, while the secondary objective is to effectively manage different types of ESS. Considering the diverse ramp-rate characteristics of generators, SCESS, and BESS, the proposed EMS exploits these differences to determine an optimal long-term schedule for minimizing shed loads. Furthermore, the proposed EMS balances the state-of-charge (SoC) of ESS and prioritizes the SCESS’s SoC levels to ensure the efficient operation of BESS and SCESS. For better computational efficiency, we introduce the receding horizon optimization method, enabling real-time EMS implementation. Further, a comparison with the fixed horizon optimization (FHO) validates its effectiveness. Simulation studies and results demonstrate that the proposed EMS efficiently manages generators, BESS, and SCESS, ensuring system resilience under generation shortages. Additionally, the proposed methodology significantly reduces the computational burden compared to the FHO technique while maintaining acceptable resilience performance.

24 POWER TRANSMISSION AND DISTRIBUTION↗