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 649 records · Page 36

Data-Driven Wind Farm Optimization Incorporating Effects of Turbulence Intensity

We demonstrate assimilation of field data from a nacelle-mounted lidar and meteorological tower into a medium-fidelity Reynolds-Averaged Navier Stokes wind farm flow model to better predict the effects of atmospheric stability. Increased predictability under a variety of atmospheric conditions can lead to more effective control design and optimization of a wind farm. In particular, atmospheric stability affects wind turbine wake propagation and, therefore, aspects of wind farm control and performance, such as active wind farm control, layout optimization, and power output. Accurately modeling wakes in different stability conditions remains a persistent challenge. This paper presents an optimization framework that leverages high-fidelity field or simulation data to correct a lower-fidelity flow model. Optimal model corrections are found by solving a regularized, high-dimensional, gradient-based optimization problem using adjoint flowfield information. We validate the trained model against large eddy simulation results, and perform separate gradient-based layout optimizations of a simulated utility-scale wind farm to maximize power. Using the data-driven model corrections, we find that atmospheric stability significantly impacts layout optimization and power production: the optimal layout for stable conditions produced 9.1% more power than the optimal layout for unstable conditions, and the optimal layout for neutral conditions underperformed by 8.5% in unstable conditions.

17 WIND ENERGY↗

Improving the cross correlation method of indirect flow measurement by water injection and modeling - 245

This paper presents a technique to improve the Cross Correlation Function (CCF) method for indirect flow measurement. Unlike other methods, the CCF method does not need a specific instrument or tool for flow measurements. Instead, the delay time between two sensors can be derived by applying the CCF. A water flow loop was constructed to investigate the CCF method for flow measurement. Because the inherent process temperature variations were small when compared to significant amounts of measurement and process noise, inaccurate flow measurements were obtained. The performance of the CCF method could be improved by increasing the system signal-to-noise ratio, which is the process noise divided by the measurement noise. For this reason, a periodic injection of different temperature water was used to increase the process temperature variation. The results were favorable. However, several design variables needed to be selected. To assist, a CFD simulation was utilized to optimize the distance between sensors and injection point, the rate of injection flow, and other related variables. The results show that CCF optimization of these variables produces better flow measurement performance over a range of flow velocities. (authors)

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Simulating a solar parabolic trough collector plant used for industrial process heat using an optimized operating scheme that utilizes flexible heat integration

This study focuses on the modeling and simulation of a novel design and operation of a solar industrial process heat (SIPH) plant that uses a parabolic trough collector system for generating process heat. The SIPH plant incorporates flexible heat integration (FHI) by having two options for heat sinks as well as a flexible collection temperature for the heat transfer fluid. Leveraging the degrees of freedom created by FHI allows an optimized operating scheme to be created which maximizes the solar heat delivered to the industrial processes. Yearly results for a 27 MW-t plant located in Salt Lake City, UT are simulated for a base case and a case utilizing optimized FHI. For the main process temperature case presented, FHI increases the total solar share by 13.0% and reduces the levelized cost of heat of the SIPH plant by 10.4% relative to the base case. Furthermore, three case study days are presented which highlight how FHI maximizes solar share, specifically during mornings, afternoons, and days with intermittent solar conditions. An environmental analysis is also performed showing that the SIPH plant can reduce harmful emissions by 15.4% for the base case and 17.4% for the FHI case compared to a plant that only utilizes a natural gas furnace for process heat. Overall, the work thoroughly explores an effective operating scheme for a SIPH plant, helping to make it more economical and environmentally beneficial.

14 SOLAR ENERGY↗

Optimized utilization of COMB3 reactive potentials in LAMMPS

An investigation to optimize the application of the third-generation charge optimized many-body (COMB3) interatomic potential and associated input parameters was carried out through the study of solid–liquid interactions in classical molecular dynamics simulations. The rates of these molecular interactions are understood through the wetting rates of water nano-droplets on a bare copper (111) surface. Implementing the Langevin thermostat, the influence of simulation time step, the number of atoms in the system, the frequency at which charge equilibration is performed, and the temperature relaxation rate are all examined. The results indicate that time steps of 0.4 fs are possible when using longer relaxation times for the system temperature, which is almost double the typical time step used for reactive potentials. Additionally, the use of the charge equilibration allows for a fewer atomic layers to be used in the Cu slab. In addition, charge equilibrium schemes do not need to be performed every time step to ensure accurate charge transfer. Interestingly, the rate of wetting for the nanodroplets is dominantly dependent on the temperature relaxation time, which is predicted to significantly change the viscosity of the water droplets. This work provides a pathway for optimizing simulations using the COMB3 reactive interatomic potential.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Multi-scale and multi-physics analysis, design optimization, and experimental validation of heat exchangers utilizing high performance, non-round tubes

Air-to-refrigerant heat exchangers (HXs) are fundamental components in HVAC&R systems, and considerable research has been dedicated designing continually smaller, lighter, and more efficient HX designs. In recent years, researchers have leveraged advancements in Computational Fluid Dynamics (CFD), Finite Element Analysis (FEA), and optimization algorithms to consider primary tube shape and topology optimization to design highly compact, high performance HXs for a multitude of applications. In this research, we present a computationally efficient, comprehensive, multi-scale, and multi-physics analysis and optimization method for air-to-refrigerant HXs featuring automated CFD and FEA simulations and approximation-assisted optimization. This methodology was utilized to design HXs with shape-optimized, non-round tubes which outperform current state-of-the-art tube-fin HXs without compromising structural integrity. The optimal HXs were shown to deliver similar thermal performance to the baseline HXs while also achieving more than 20% reductions in airside pressure drop and core envelope volume and more than 30% reduction in internal volume. Comprehensive experimental validation of the optimization methodology was conducted through the testing of two prototypes in a standardized wind tunnel facility under multiple operating conditions. For prototype #1 under dry evaporator conditions, the predicted heat load agreed within ± 10% of measured values and the predicted airside pressure drop agreed within ± 30%, while for dehumidifying conditions, the predicted sensible and latent heat loads agreed within ± 10% and ± 20% of the measured values, respectively. For prototype #2, the predicted condenser heat load agreed within ± 3.0% of measured values, and the predicted airside pressure drop agreed within ± 27%. The acceptable agreement between simulation and experimental results for the present application highlights the flexibility of the novel optimization methodology to design next generation HXs with improved performance and reduced volume, weight, and environmental impact.

42 ENGINEERING↗

Optimal control of anthropogenic carbon dioxide emissions through technological options: a modeling study

The anthropogenic emission of carbon dioxide (CO{sub 2}) is the prime culprit for the menace of global warming. To achieve the goal of mitigation of global warming, it is crucial to curb the anthropogenic carbon dioxide emissions. The prime anthropogenic source of CO{sub 2} is fossil fuel burning. In this paper, we propose a nonlinear mathematical model to study the impact of technological options, used for the reduction of CO{sub 2} emissions from fossil fuel burning and industrial processes, on the control of atmospheric CO{sub 2}. In the modeling process, it is assumed that the technological options are implemented to curb the CO{sub 2} emissions from the source at a rate proportional to the anthropogenic CO{sub 2} emissions. Model analysis reveals that the atmospheric level of CO{sub 2} can be effectively reduced by increasing the implementation rate of technological options and their efficiency. The strategies which optimally reduce atmospheric CO{sub 2} levels while minimizing the cost associated with the implementation of technological options are identified using optimal control theory. Numerical simulation has been carried out to illustrate theoretical results.

54 ENVIRONMENTAL SCIENCES↗

Interatomic Potential in the Nonequilibrium Warm Dense Matter Regime

We present a new measurement of lattice disassembly times in femtosecond-laser-heated polycrystalline Au nanofoils. The results are compared with molecular dynamics simulations incorporating a highly optimized, embedded-atom-method interatomic potential. For absorbed energy densities of 0.9 – 4.3 MJ / kg , the agreement between the experiment and simulation reveals a single-crystal-like behavior of homogeneous melting and corroborates the applicability of the interatomic potential in the nonequilibrium warm dense matter regime. For energy densities below 0.9 MJ / kg , the measurement is consistent with nanocrystal behavior where melting is initiated at the grain boundaries.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Active Learning Surrogates for Integrating Electron Microscopy and Computational Insights from Simulations in Autonomous Experiments

Artificial Intelligence (AI) combined with simulations and experiments has great potential to accelerate scientific discovery across technology and pharmaceuticals. However, the gap between simulations and experiments is challenging due to disparities in time and scale, making it difficult to estimate properties like energy and electronic states from experiments, and to provide feedback based on theoretical insights.Our research addresses the challenge by developing unique deep kernel based surrogate models that learns from microscopic images, mapping structural features to energy differences from defect formation. We start with full-training using simulated images to determine optimal settings, establishing a baseline for active learning. Using these settings from the baseline, active learning is trained, and predicts structures along simulation trajectories based on uncertainty and energetic stability, thus reducing data requirements, simulation time and computational costs. The results demonstrate that the model achieves a low average error margin of approximately 0.03 meV, indicating good performance. To enhance feature extraction and reconstruction capabilities, we developed an autoencoder-decoder as additional surrogate to create latent space to capture essential features, enabling precise comparisons between simulations and experiments. The results from this model achieved a reconstruction loss of around 0.2 and accurately reconstructed molecular structures.Overall, this work advances the steering of experiments through computational simulations by employing a surrogate models that actively predicts the trajectories of structural evolution, achieving time-to-solution comparable to experimental measurements.

Saranathan, Gayathri [Hewlett-Packard]↗

Studies of Front-End Electronics for High-Precision Timing Measurements with LGADs [Poster]

In Time-of-Arrival (ToA) measurements, simulations can be used to optimize precision detection and implementation methods of front-end electronics. The simulation tested in this project takes the response of a filter (in the s-domain) from a noisy input current signal from the detector, an LGAD (Low-Gain Avalanche Diode), and uses a threshold to determine the ToA value for the signal. Measuring the ToA values for a large sample of events and comparing them to a reference statistically produces the timing precision of the simulation. Here, the sub-algorithm to model LGAD events with piecewise linear (PWL) approximations (for more efficient Laplace transforms) was optimized. Using the Fourier Discrete Cosine Transform, 25 points could be selected to describe the primary behavior of the signal. Results showed that timing precision is on par with previous methods. We also found that the new algorithm is faster and equally consistent in producing accurate PWL approximations. Testi ng with larger samples is required to fully verify the new algorithm.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Optimization Through Multi-Fidelity Modeling

We present a novel method for optimizing parameter selection for simulations with an evaluation budget. We start with an existing method for building a multi-fidelity model out of many low-fidelity simulations and few high-fidelity simulations. We propose a novel method to simplify parameter selection without sacrificing performance. We verify these results and compare with existing literature. Next, we propose a novel algorithm which uses this difference model to suggest new points in the parameter design space to simulate. We add each point we simulate to the model to improve its quality for the next iteration. The algorithm trades off reducing the uncertainty of the existing model with optimization of the objective. The first is more useful when a large fraction of the computation budget remains. The second is more useful when a small fraction of the computation budget remains. Our method converges to the optimum by using a high-fidelity evaluation for just 16 of the 427 points. Our method is general enough to work if there is no low-fidelity model. Furthermore, it is agnostic to the underlying physics of the problem. Therefore, both the low-fidelity and high-fidelity models can be generated by any arbitrary function, including simulations and physical experiments.

97 MATHEMATICS AND COMPUTING↗

A Component-Sizing Methodology for a Hybrid Electric Vehicle Using an Optimization Algorithm

Many leading companies in the automotive industry have been putting tremendous effort into developing new powertrains and technologies to make their products more energy efficient. Evaluating the fuel economy benefit of a new technology in specific powertrain systems is straightforward; and, in an early concept phase, obtaining a projection of energy efficiency benefits from new technologies is extremely useful. However, when carmakers consider new technology or powertrain configurations, they must deal with a trade-off problem involving factors such as energy efficiency and performance, because of the complexities of sizing a vehicle’s powertrain components, which directly affect its energy efficiency and dynamic performance. As powertrains of modern vehicles become more complicated, even more effort is required to design the size of each component. This study presents a component-sizing process based on the forward-looking vehicle simulator “Autonomie” and the optimization algorithm “POUNDERS”; the supervisory control strategy based on Pontryagin’s Minimum Principle (PMP) assures sufficient computational system efficiency. We tested the process by applying it to a single power-split hybrid electric vehicle to determine optimal values of gear ratios and each component size, where we defined the optimization problem as minimizing energy consumption when the vehicle’s dynamic performance is given as a performance constraint. The suggested sizing process will be helpful in determining optimal component sizes for vehicle powertrain to maximize fuel efficiency while dynamic performance is satisfied. Indeed, this process does not require the engineer’s intuition or rules based on heuristics required in the rule-based process.

33 ADVANCED PROPULSION SYSTEMS↗

Optimization of Ag Electrocatalyst Performance for CO2 to CO Conversion: Pairing Atomic Simulations with Experiments

Density-functional theory- based calculations that complement our series of experimental efforts on optimizing Ag electrocatalyst performance for CO2 to CO conversion were presented. Three key findings are drawn out from the combined UHV surface science, STM and electrochemical measurements: (1) Using a series of Ag nanoparticles with 2-6 nm average diameters, CO2 reduction reaction (CO2RR) activity increases, with particle between 2 nm and ∼4 nm demonstrating the highest combination of activity and selectivity; (2) Electronic metal−support interactions (EMSIs) between Ag and C dramatically improve CO2RR performance as evidenced by a scaling relationship between particle size and the relative Ag−C EMSIs strength, which improves the CO2-to-CO Faradaic efficiency of sub-2 nm Ag particles from 2 to ∼100% and increases the CO turnover frequency ∼15-fold compared with similarly sized bare Ag particles; (3) The performance Ag electrocatalysts is improved when supported on S-doped C materials. Computational modeling of 1−10 nm Ag particles predicts a nearly identical size-dependent trend with maximum CO2RR activity predicted for 3.7nm particles. The calculations support the promotional effect of C materials, showing a large charge transfer of 1.02 e from Ag clusters to defective C and a 0.41 eV less endergonic step of forming COOH intermediate on Ag/defective-C compared to the Ag/C system. Calculations indicate a more favorable energetic pathway of CO2-to-CO at the C-S-Ag interface, consistent with experiments.

CO2 utilization↗

A phase transition for finding needles in nonlinear haystacks with LASSO artificial neural networks

To fit sparse linear associations, a LASSO sparsity inducing penalty with a single hyperparameter provably allows to recover the important features (needles) with high probability in certain regimes even if the sample size is smaller than the dimension of the input vector (haystack). More recently learners known as artificial neural networks (ANN) have shown great successes in many machine learning tasks, in particular fitting nonlinear associations. Small learning rate, stochastic gradient descent algorithm and large training set help to cope with the explosion in the number of parameters present in deep neural networks. Yet few ANN learners have been developed and studied to find needles in nonlinear haystacks. Driven by a single hyperparameter, our ANN learner, like for sparse linear associations, exhibits a phase transition in the probability of retrieving the needles, which we do not observe with other ANN learners. To select our penalty parameter, we generalize the universal threshold of Donoho and Johnstone (Biometrika 81(3):425–455, 1994) which is a better rule than the conservative (too many false detections) and expensive cross-validation. In the spirit of simulated annealing, we propose a warm-start sparsity inducing algorithm to solve the high-dimensional, non-convex and non-differentiable optimization problem. We perform simulated and real data Monte Carlo experiments to quantify the effectiveness of our approach.

97 MATHEMATICS AND COMPUTING↗

Performance assessment of 3D printed multi-material energy absorber for automotive bumper: pedestrian lower extremity protection

Designing an energy absorber for automotive bumpers involves balancing low-speed and high-speed impacts to ensure safety, reduce repair costs, and meet regulatory standards. Here, this study explores a novel design using multi-material 3D printing and structural optimization to fabricate a lightweight and cost-efficient energy absorber. The design effectively dissipates energy in low-speed collisions and minimizes force transmission in high-speed pedestrain impacts, helping to meet both safety and performance requirements. The energy absorber design combines 20% carbon fiber-reinforced acrylonitrile butadiene styrene (CF-ABS) and thermoplastic polyurethane (TPU) for optimal stiffness and flexibility. It uses 3D-printed lattice structures optimized through finite element simulations to help meet both low-speed and high-speed impact requirements. Full-scale energy absorbers were 3D-printed using optimized CF-ABS/TPU blends and tested under high-speed impact using the Flexible Pedestrian Legform Impactor (Flex-PLI). For fair comparison, a baseline bumper with a traditional triangular lattice structure, also 3D-printed from the same CF-ABS/TPU materials, was similarly tested. Interestingly, both the optimized and baseline 3D-printed energy absorbers showed nearly identical performance, successfully meeting injury limits. Their performances were also benchmarked against an injection-molded energy absorber. While both 3D-printed and injection-molded designs met injury limits, the 3D-printed absorber exhibited a higher tibia bending moment, indicating an opportunity for further optimization. A Techno-Economic Analysis compared the costs of producing energy absorbers using traditional manufacturing and 3D printing. The analysis highlighted that 3D printing offers cost benefits for low to medium production volumes, with the total cost per energy absorber at ∼ $\$$74, compared to traditional methods that become economical beyond 2000 units.

Additive manufacturing↗

Adaptive local basis set for Kohn–Sham density functional theory in a discontinuous Galerkin framework II: Force, vibration, and molecular dynamics calculations

Recently, we have proposed the adaptive local basis set for electronic structure calculations based on Kohn–Sham density functional theory in a pseudopotential framework. The adaptive local basis set is efficient and systematically improvable for total energy calculations. In this paper, we present the calculation of atomic forces, which can be used for a range of applications such as geometry optimization and molecular dynamics simulation. We demonstrate that, under mild assumptions, the computation of atomic forces can scale nearly linearly with the number of atoms in the system using the adaptive local basis set. We quantify the accuracy of the Hellmann–Feynman forces for a range of physical systems, benchmarked against converged planewave calculations, and find that the adaptive local basis set is efficient for both force and energy calculations, requiring at most a few tens of basis functions per atom to attain accuracies required in practice. Since the adaptive local basis set has implicit dependence on atomic positions, Pulay forces are in general nonzero. However, we find that the Pulay force is numerically small and systematically decreasing with increasing basis completeness, so that the Hellmann–Feynman force is sufficient for basis sizes of a few tens of basis functions per atom. We verify the accuracy of the computed forces in static calculations of quasi-1D and 3D disordered Si systems, vibration calculation of a quasi-1D Si system, and molecular dynamics calculations of H{sub 2} and liquid Al–Si alloy systems, where we show systematic convergence to benchmark planewave results and results from the literature.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Offline Arterial Signal Timing Optimization based on Virtual Phase Link Model - A Real-world Case Study

Conventional signal timing for arterial usually takes bottom-up approaches. Engineers optimize each individual intersection first and then coordinate them by adjusting the offsets. This work is based on the Virtual Phase-Link (VPL) model, a street traffic model designed for online traffic model predictive control, to obtain a top-down offline arterial signal timing. We Studied the Shallowford Rd. in Chattanooga, TN and found that the inconsistency in intersection capacities along the arterial could lead to some intersections becoming bottlenecks. Signal timing is a significant factor that affect the intersection capacities. We realized that the VPL-based model can guarantee the consistency in intersections along an arterial. We therefore adopted the VPL-based model and developed an offline signal timing optimization approach. The proposed timing derived from the VPL-based offline signal timing optimization showed very good results in simulation. The Chattanooga Department of Transportation adopted the optimized timing obtained from the proposed approach and gave positive feedbacks to the research team. We also collected field experiment data, which demonstrated overall energy reductions and speed improvements on some sections of the Shallowford Rd. arterial. We will continue the experiment when the COVID-19 pandemic impact subsides to have a more robust quantitative evaluation.

Wang, Qichao↗

Short-depth QAOA circuits and quantum annealing on higher-order ising models

Abstract We present a direct comparison between QAOA (Quantum Alternating Operator Ansatz), and QA (Quantum Annealing) on 127 qubit problem instances. QAOA with p = 1, 2 rounds is executed on the 127 qubit heavy-hex graph gate-model quantum computer ibm_washington, using on-device grid-searches for angle finding, and QA is executed on two Pegasus-chip D-Wave quantum annealers. The problems are random Ising models whose connectivity matches heavy-hex graphs and the Pegasus graph connectivity, and optionally include hardware-compatible cubic terms ( Z Z Z terms). The QAOA circuits are heavily optimized and of extremely short depth, with a CNOT depth of 6 per round, which allows whole chip usage of the heavy-hex lattice. QAOA and QA are both compared against simulated annealing and the optimal solutions are computed exactly using CPLEX. The noiseless mean QAOA expectation values for p = 1, 2 are computed using classical light-cone based simulations. We find QA outperforms QAOA on the evaluated devices.

127 qubits↗

Stochastic Optimization to Find Optimum Beginning-of-Life Core Configuration of Stable Salt Reactor with Online Refueling

A stochastic optimization method has been developed to find an optimum equilibrium cycle core configuration of the waste-burning stable salt reactor, which is a fast-spectrum molten salt reactor with frequent online refueling. An optimum core configuration was determined with the goal of minimizing radial power peaking. Because of the vast number of potential candidate core configurations, stochastic optimization was applied based on simulated annealing and an additional acceleration method, which screened out unpromising core configurations. It has been demonstrated that the developed stochastic optimization method successfully finds the optimal core configuration regardless of the initial guess and outperforms the gradient descent approach. In addition, it has been observed that the use of a so-called out-in core configuration as the initial guess speeds up convergence of the iterative solution more than five times. Based on the searched optimum equilibrium cycle core configuration, new beginning-of-life (BOL) core configurations have been developed. In conclusion, the new BOL core configurations will be used in developing optimum refueling strategies.

Moltex static salt reactor↗