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 523 records · Page 29
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Microhole drilling technology utilizing a golden section search algorithm

A fundamental issue in microhole drilling is that delivering high weight-on-bit (WOB), high torque rotational horsepower to a conventional drill bit does not scale down to the hole sizes necessary to realize the envisioned cost savings An optimization algorithm called a golden section search (GSS) was used to systematically identify the preferred WOB for a given set of conditions. This research focused on implementing and evaluating two low WOB drilling technologies for microhole drilling: - Laser-assisted mechanical drill, which was tested in the laboratory - Lightly modified commercial off-the-shelf (COTS) percussive hammer, which was tested in a limited field test. Data were collected for microhole GSS using WOB optimization via simulation as well as at the Blue Canyon Dome Site in Socorro, NM. Information on the attached files and folders are as follows: - the .tdms files are LabView data files, which can be opened within Excel using a .tdms add-in or using a Matlab .tdms converter - the .tdms_index files are part of the .tdms file structure - sampling rate, column headers, and length data within the .tdms files follow SOP when utilizing Excel and/or Matlab as described above

15 GEOTHERMAL ENERGY↗
DOE OSTI2025

Optimizing temperature distributions for training neural quantum states using parallel tempering

Parametrized artificial neural networks (ANNs) can be very expressive ansatzes for variational algorithms, reaching state-of-the-art energies on many quantum many-body Hamiltonians. Nevertheless, the training of the ANN can be slow and stymied by the presence of local minima in the parameter landscape. One approach to mitigate this issue is to use parallel tempering methods, and in this work, we focus on the role played by the temperature distribution of the parallel tempering replicas. Using an adaptive method that adjusts the temperatures in order to equate the exchange probability between neighboring replicas, we show that this temperature optimization can significantly increase the success rate of the variational algorithm with negligible computational cost by eliminating bottlenecks in the replicas' random walk. Furthermore, we demonstrate this using two different neural networks, a restricted Boltzmann machine and a feedforward network, which we use to study a toy problem based on a permutation invariant Hamiltonian with a pernicious local minimum and the 𝐽 1 −𝐽 2 model on a rectangular lattice.

Neural network simulations↗
NASA NTRS2015

Space Synthetic Biology Project

Synthetic biology is an effort to make genetic engineering more useful by standardizing sections of genetic code. By standardizing genetic components, biological engineering will become much more similar to traditional fields of engineering, in which well-defined components and subsystems are readily available in markets. Specifications of the behavior of those components and subsystems can be used to model a system which incorporates them. Then, the behavior of the novel system can be simulated and optimized. Finally, the components and subsystems can be purchased and assembled to create the optimized system, which most often will exhibit behavior similar to that indicated by the model. The Space Synthetic Biology project began in 2012 as a multi-Center effort. The purpose of this project was to harness Synthetic Biology principals to enable NASA's missions. A central target for application was to Environmental Control & Life Support (ECLS). Engineers from NASA Marshall Space Flight Center's (MSFC's) ECLS Systems Development Branch (ES62) were brought into the project to contribute expertise in operational ECLS systems. Project lead scientists chose to pursue the development of bioelectrochemical technologies to spacecraft life support. Therefore, the ECLS element of the project became essentially an effort to develop a bioelectrochemical ECLS subsystem. Bioelectrochemical systems exploit the ability of many microorganisms to drive their metabolisms by direct or indirect utilization of electrical potential gradients. Whereas many microorganisms are capable of deriving the energy required for the processes of interest (such as carbon dioxide (CO2) fixation) from sunlight, it is believed that subsystems utilizing electrotrophs will exhibit smaller mass, volume, and power requirements than those that derive their energy from sunlight. In the first 2 years of the project, MSFC personnel conducted modeling, simulation, and conceptual design efforts to assist the project in selecting the best approaches to the application of bioelectrochemical technologies to ECLS. Figure 1 shows results of simulation of charge transport in an experimental system. Figure 2 shows one of five conceptual designs for ECLS subsystems based on bioelectrochemical reactors. Also during the first 2 years, some work was undertaken to gather fundamental data (conductivities, overpotentials) relevant to the modeling efforts.

Howard, David↗
DOE OSTI2021

Subsurface Characterization of Hydraulic Fracture Test Site-2 (HFTS-2), Delaware Basin

Hydraulic Fracturing Test Site-2 (HFTS-2) is a field-based research experiment performed in the Wolfcamp Formation of the Permian (Delaware) Basin. This paper focuses on integration, advanced geological characterization, and 3D subsurface modeling of the comprehensive HFTS-2 dataset. The study showcases a multidisciplinary reservoir characterization approach that incorporates geology, petrophysics, geochemistry, geomechanics, microseismic, and subsurface engineering analysis. Subsurface characterization of organic-rich mudstone formations requires understanding complex hydraulic fracture network growth in relation to inherent lithology, geomechanical properties, and interaction with pre-existing natural fractures. This paper presents a characterization workflow incorporating pre- and post-stimulation subsurface data, unique to the HFTS-2 dataset. The study integrated: (1) rock properties from logs, cores, and thin sections; (2) natural and hydraulic fracture descriptions from cores and image logs; (3) local and regional stresses; (4) geomechanics; (5) microseismic; (6) fiber optic (FO) and bottomhole pressure gauge (BHPG) response; and (7) produced fluids analysis. During a stimulation treatment, creation of the stimulated rock volume (SRV) is influenced by several subsurface factors. Key contributing factors include structural context, stress conditions, lithology, facies architecture, pre-existing natural fractures, and geomechanical properties. The HFTS-2 subsurface data integration indicates that the SRV is comprised of a complex juxtaposition of hydraulic fracture swarms, as evidenced by image logs analysis, core description, and microseismic monitoring. The HFTS-2 microseismic event density was used to generate 3D heat maps that serve as a representative SRV footprint, corroborated by secondary datasets. These maps were further integrated with petrophysical and geomechanical characteristics, as well as responses from FO and BHPG, to estimate the lateral and vertical dimensions of the effective fractures. The geological characterization for the HFTS-2 dataset combined with 3D modeling for petrophysical and geomechanical properties provides a strong foundation for subsurface simulation and optimization studies. Downloaded from http://onepetro.org/URTECONF/proceedings-pdf/21URTC/1-21URTC/D011S005R001/2477501/urtec-2021-5243-ms.pdf/1 by Carol Worster on 28 February 2022 URTeC 5243 The workflow improved our understanding of HFTS-2 hydraulic fracture propagation and characteristics in relation to offset pressure depletion and interaction with pre-existing natural fractures. Analysis showed that fracture geometry varies by stage and by well, and a complex fracture network is generated with varying fracture density. The multidisciplinary workflow presented herein for integration and characterization serves as a foundation to evaluate completion efficiency and estimate areal and vertical stimulation and depletion extent for the project. Furthermore, the workflow and learnings can also be transferred to other unconventional plays.

58 GEOSCIENCES↗
DOE OSTI2019

Ten questions on urban building energy modeling

Buildings in cities consume up to 70% of all primary energy. To achieve cities’ energy and climate goals, it is necessary to reduce energy use and associated greenhouse gas emissions in buildings through energy conservation and efficiency improvements. Computational tools empowered with rich urban datasets can model performance of buildings at the urban scale to provide quantitative insights for stakeholders and inform their decision making on urban energy planning, as well as building energy retrofits at scale, to achieve efficiency, sustainability, and resilience of urban buildings. Designing and operating urban buildings as a group (from a city block to a district to an entire city) rather than as single individuals requires simulation and optimization to account for interactions among buildings and between buildings and their surrounding urban environment, and for district energy systems serving multiple buildings with diverse thermal loads across space and time. When hundreds or more buildings are involved in typical urban building energy modeling (UBEM) to estimate annual energy demand, evaluate design or retrofit options, and quantify impacts of extreme weather events or climate change, it is crucial to integrate urban datasets and UBEM tools in a seamless automatic workflow with cloud or high-performance computing for users including urban planners, designers and researchers. This paper presents ten questions that highlight significant UBEM research and applications. The proposed answers in this work aim to stimulate discussion and provide insights into the current and future research on UBEM, and more importantly, to inspire new and important questions from young researchers in the field.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
NASA NTRS1990

Optimal fabrication processes for unidirectional metal-matrix composites: A computational simulation

A method is proposed for optimizing the fabrication process of unidirectional metal matrix composites. The temperature and pressure histories are optimized such that the residual microstresses of the composite at the end of the fabrication process are minimized and the material integrity throughout the process is ensured. The response of the composite during the fabrication is simulated based on a nonlinear micromechanics theory. The optimal fabrication problem is formulated and solved with non-linear programming. Application cases regarding the optimization of the fabrication cool-down phases of unidirectional ultra-high modulus graphite/copper and silicon carbide/titanium composites are presented.

Saravanos, D. A.↗
NASA NTRS1990

Optimal fabrication processes for unidirectional metal-matrix composites - A computational simulation

A method is proposed for optimizing the fabrication process of unidirectional metal matrix composites. The temperature and pressure histories are optimized such that the residual microstresses of the composite at the end of the fabrication process are minimized and the material integrity throughout the process is ensured. The response of the composite during the fabrication is simulated based on a nonlinear micromechanics theory. The optimal fabrication problem is formulated and solved with nonlinear programming. Application cases regarding the optimization of the fabrication cool-down phases of unidirectional ultra-high modulus graphite/copper and silicon carbide/titanium composites are presented.

Saravanos, D. A.↗
DOE OSTI2022

Comparison of Some RANS Solvers

We will take a look at solving the Reynolds-averaged Navier-Stokes (RANS) equations that are encountered in the context of wind farm performance simulations and optimizations. We will compare some of the more popular ways to solve these equations with a focus on using iterative solvers for the linear solve. We will compare their performance and reliability to a direct solve as we scale the problem both by adding more and parallel resources and by increasing the size of the domain, both in two and three dimensions. There are many strategies that can be applied to solving the RANS equations, some are very efficient, while others are very insensitive to, for example, the Reynolds number. The first contender we will consider is the Pressure-Convection-Diffusion (PCD) preconditioner. Early results suggest that PCD is indeed a very efficient solver, in particular in two dimensions, as long as the Reynold's number remains small. Next we will try to reorder our degrees of freedom such that we can use GMRES with ILU for our linear solve. Another popular choice we will consider for solving the RANS equations is SIMPLE (and its derivatives). For all of our implementations we make use of either FEniCS or Firedrake, basing our work on both existing implementations of some of these solvers while also writing new extensions for others.

CFD↗
DOE OSTI2019

A Parallel-in-Time Multigrid Approach to Constrained Optimization

Valuable computer simulations such as weather prediction, traffic flow, and resource allocation are time dependent problems. One such example is the viscous Burgers’ equation which can model various systems such as fluids and traffic flow. Further, many of these models require optimization routines to meet the needs of many businesses. Unfortunately, traditional sequential-in-time algorithms have become bottlenecks in solution model generation. This is because computer architectures have begun to rely more and more on concurrency (parallel processing) as opposed to processor power to provide greater performance. To avoid this bottleneck, many researchers use parallelization methods to take advantage of the increasing number of processors in computers. This approach allows for variables in a problem to be solved for simultaneously. As such, parallelization is a crucial requirement to achieve high-performance computing.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Impact of aquifer properties on the extent and timeline of CO 2 trapping

Abstract Geologic CO 2 sequestration in porous saline aquifers is a promising approach to reducing atmospheric concentrations of CO 2 . Reactive transport simulations provide the opportunity to analyze which factors influence geochemical reactivity in the reservoir, understand those most important for promoting CO 2 trapping, and assess individual sites. Field‐scale aquifer characterization is time and resource intensive such that here, reactive transport simulations are leveraged to enhance understanding of selected aquifer properties including porosity, permeability, depth of storage, and carbonate mineralogy on the overall CO 2 trapping potential to better select sites promoting geochemical reactivity for CO 2 trapping. There are different mechanisms for sequestrating CO 2 . Once injected, CO 2 will dissolve into the brine to create an acidic environment, resulting in the dissolution of pre‐injection formation minerals. Released ions can reprecipitate as secondary minerals. The dissolved CO 2 and mineralized CO 2 are considered as a more secure form of CO 2 trapping in this study compared to the free supercritical CO 2 . Here, a framework leveraging a controlled set of field scale simulations is developed to facilitate rapid, optimized site selection. Simulations vary aquifer properties to understand the impact of each unique property on CO 2 trapping, tracking, and comparing the amount of supercritical, aqueous, and mineralized CO 2 . The rate at which the CO 2 injected into the aquifer is converted to aqueous or mineralized CO 2 is newly defined here as the sequestration efficiency and used to compare simulation results. The reservoir depth and fraction of carbonate minerals in the formation are shown to be more important factors than reservoir porosity and permeability in affecting CO 2 trapping. However, the impact of aquifer properties on the evolution of injected CO 2 depends on the stage of the sequestration project. © 2023 Society of Chemical Industry and John Wiley & Sons, Ltd.

Energy & Fuels↗
DOE OSTI2023

Novel gas exposure system for the controlled exposure of plants to gaseous hydrogen fluoride

Abstract Plants can serve as sensitive bioindicators of the presence of contaminant vapors in the atmosphere. This work describes a novel laboratory-based gas exposure system capable of calibrating plants as bioindicators for the detection and delineation of the atmospheric contaminant hydrogen fluoride (HF) as a preparatory step for monitoring release emissions. To evaluate changes in plant phenotype and stress-induced physiological effects attributed to HF alone, the gas exposure chamber must have additional controls to simulate otherwise optimal plant growth conditions including variables such as light intensity, photoperiod, temperature, and irrigation. The exposure system was designed to maintain constant growth conditions during a series of independent experiments that varied between optimal (control) and stressful (HF exposure) conditions. The system was also designed to ensure the safe handling and application of HF. An initial system calibration introduced HF gas into the exposure chamber and monitored HF concentrations by cavity ring-down spectroscopy for a 48-h period. Stable concentrations inside the exposure chamber were observed after approximately 15 h, and losses of HF to the system ranged from 88 to 91%. A model plant species ( Festuca arundinacea ) was then exposed to HF for 48 h. Visual phenotype stress-induced responses aligned with symptoms reported in the literature for fluoride exposure (tip dieback and discoloration along the dieback transition margin). Fluoride concentrations in exposed tissues compared to control tissues confirmed enhanced fluoride uptake due to HF exposure. The system described herein can be applied to other reactive atmospheric pollutants of interest in support of bioindicator research.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Machine learning application to single channel design of molten salt reactor

This study proposes a robust approach to quickly design a nuclear reactor core and explores the best performing machine learning (ML) technique for predicting feature parameters of the core. Here we implemented the approach into a hypothetical channel of molten salt reactors to demonstrate the applicability of the method. We prepared a Python tool, named Plankton, which couples to a reactor physics code and an optimization tool, and imports ML methods. The tool performs three consecutive phases: reactor database generation, machine learning application, and design optimization. We identified the extra trees method as the best performing estimator. With the estimator, we found nine optimum designs in total, one for each fuel-salt pair, and estimated all the performance metrics of the designs with a <5% prediction error compared to their actual values. U-Pu-NaCl fuel-salt gave promising results with the highest conversion ratio, the most negative feedback coefficient, and the lowest fast flux.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2022

Kinetic model development and Bayesian uncertainty quantification for the complete reduction of Fe-based oxygen carriers with CH 4 , CO, and H 2 for chemical looping combustion

In this work, three kinetic models are developed and calibrated for the complete multi-step reduction of an Fe-based oxygen carrier (OC) particle with CH 4 , CO, and H 2 , using data from thermogravimetric analysis. The complete reduction rate profiles exhibit complex dynamics whose trajectory is significantly different depending on the reducing gas. A Bayesian model building and parameter estimation framework is applied for simultaneous parameter and model structure uncertainty quantification. The final models show excellent agreement between model predictions and calibration data, as well as new data not used for calibration (for the reduction of the OC with CH 4 ). Parameter uncertainty is quantified by determining joint posterior distribution, and model structure uncertainty is addressed by incorporating Gaussian process stochastic functions (represented by Bayesian smoothing splines) into the kinetic models. The final kinetic models with discrepancy functions are readily employable in equation-oriented simulation and optimization platforms.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Machine learning modeling and model predictive control of a closed-circuit reverse osmosis system

Closed-circuit reverse osmosis (CCRO) offers a flexible and energy-efficient alternative to conventional reverse osmosis by operating in a semi-batch mode that recycles brine, enabling higher recovery rates and reduced specific energy consumption (SEC). However, developing accurate, system-level dynamic models for CCRO remains challenging due to its nonlinear, multi-phase operation and sensitivity to variable feed water conditions. Traditional modeling approaches, such as NARMAX (nonlinear autoregressive moving average with exogenous inputs), often struggle to generalize across varying inlet feed concentrations, necessitating frequent parameter re-estimation and limiting their utility for real-time control applications. To address these limitations, we developed a long short-term memory (LSTM) neural network model trained on an extensive experimental data set from a CCRO pilot plant. The model accepts three inputs, feed flow rate, recirculation flow rate, and initial feed conductivity, and predicts three key outputs: reject conductivity, feed pump power draw, and recirculation pump power draw. We validated the LSTM model against experimental data, demonstrating its ability to distinguish between different feed conductivities and adapt to variable flow rates. Subsequently, we incorporated the LSTM model within a nonlinear model predictive control (MPC) scheme and conducted closed-loop simulations to optimize the integrated SEC (iSEC). In conclusion, the results project up to a 6% reduction in iSEC by using MPC to optimize performance over the entire experiment duration, without requiring any random excitation for data collection or parameter re-estimation.

Desalination↗
DOE OSTI2020

Minimizing total energy cost and tardiness penalty for a scheduling-layout problem in a flexible job shop system: A comparison of four metaheuristic algorithms

Job scheduling and machine layout are interrelated in improving energy consumption (EC) and productivity measures such as tardiness and represent two important decisions that must be made by manufacturers. This interdependency can be explained by transportation time, which connects scheduling and layout. Scheduling and layout, however, have not been thoroughly studied in conjunction using an integrated model in the context of sustainable manufacturing. Hence, we propose an energy-aware optimization model in which scheduling is integrated with layout in a single-level framework. More specifically, a single objective function is defined to minimize the facility energy cost and the job tardiness penalty, which control EC and tardiness respectively in a flexible job shop system. In order to model machine EC more accurately, we also consider three different machine states: a processing state and two idle states. Our case studies show that the integrated model exhibits better performance in controlling manufacturing EC and job tardiness than a non-integrated model in which machine locations are uncontrollable and transportation times between machines are unchangeable. To deal with large-sized problems, we also introduce four new metaheuristics. The performances of these new metaheuristics are compared in terms of objective function values and CPU times using various case studies. The results indicate that a hybrid ant colony optimization and simulated annealing (ACO-SA) algorithm provides better performance than the other algorithms. Specifically, our case studies show that the integrated model using ACO-SA can improve the objective function value by around 5% when compared to the non-integrated model.

42 ENGINEERING↗
DOE OSTI2023

Equation-based and data-driven modeling: Open-source software current state and future directions

Here, a review of current trends in scientific computing reveals a broad shift to open-source and higher-level programming languages such as Python and growing career opportunities over the next decade. Open-source modeling tools accelerate innovation in equation-based and data-driven applications. Significant resources have been deployed to develop data-driven tools (PyTorch, TensorFlow, Scikit-learn) from tech companies that rely on machine learning services to meet business needs while keeping the foundational tools open. Open-source equation-based tools such as Pyomo, CasADi, Gekko, and JuMP are also gaining momentum according to user community and development pace metrics. Integration of data-driven and principles-based tools is emerging. New compute hardware, productivity software, and training resources have the potential to radically accelerate progress. However, long-term support mechanisms are still necessary to sustain the momentum and maintenance of critical foundational packages.

97 MATHEMATICS AND COMPUTING↗