Search NASA⌕ Search

SEARCH · Search NASA

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 1,585 records · Page 88

Optimized Renewable Energy Integration for EV High-Power Dynamic Wireless Charging Systems

Dynamic wireless charging for electric vehicles is an emerging technology to reduce on-board battery size and extend driving range. Due to its unique characteristic of vehicle-speed-related pulse-like load profile, the high-power dynamic wireless charging system (DWCS) introduces high stress to the utility grid. In this paper, an optimization model for renewable energy integration in the DWCS is proposed to mitigate the grid impact and minimize the operation costs of the whole system. As the load profile of DWCS is related to the traffic volume and various approaching vehicle speeds, the annual average daily traffic data and a stochastic model are used to develop 24-hour load profile of DWCS. To find a tradeoff between grid impact mitigation and operation costs minimization, relationships among power demand from power grid, photovoltaic (PV) capacity, wind energy (WE) capacity and energy storage (ES) capacity are analyzed, and the optimization objective and constraints are developed. Numerical simulation results demonstrate that energy storage integration can greatly mitigate the grid impact of DWCS, and optimal ratio of PV and WE can significantly reduce the operation cost of DWCS.

Zeng, Rong↗

Experimental Results for Mars Rotorcraft Airfoils (roamx-0201 and clf5605) at Low Reynolds Number and Compressible Flow in a Mars Wind Tunnel

Experimental results are obtained for a roamx-0201 type airfoil and the clf5605 airfoil at highsubsonic, low Reynolds number conditions using the Tohoku University Mars Wind Tunnel, Japan. The tests are conducted at a Mach number of M = 0.60, and a Reynolds number of Re = 20,000 to reflect representative aerodynamics of a rotor blade for Mars exploration. The angle of attack is varied between α = −2.0 deg and α = 6.0 deg. The roamx-0201 type airfoil is an unconventional airfoil optimized for the chosen tunnel operating conditions using the Evolutionary aLgorithm for Iterative Studies of Aeromechanics (ELISA), developed under the Rotor Optimization for the Advancement of Mars eXploration (ROAMX) project. ELISA is utilized here to optimize aerodynamic airfoil performance using a Genetic Algorithm and two-dimensional high-fidelity CFD simulations, ultimately resulting in a Pareto-optimal airfoil set. The clf5605 airfoil is the outboard airfoil used on the Ingenuity Mars Helicopter and provides a baseline against which the roamx-0201, as well as possible future airfoil profiles for the compressible low Reynolds number regime, can be compared against. Lift and drag data are recorded using a balance, pressure distributions are obtained using Pressure Sensitive Paint (PSP) application, and Schlieren images are obtained to visualize the flowfield. The data is tabulated to aid future research.

Roamx↗

Bayesian optimization for inverse calibration of expensive computer models: A case study for Johnson-Cook model in machining

Inverse model calibration for identifying the constitutive model parameters can be computationally demanding for expensive-to-evaluate simulation models. Here, this paper presents a modified Bayesian optimization (BO) method, denoted as BO-bound, that incorporates theoretical bounds on the quantity of interest. A case study for the inverse calibration of the Johnson Cook (J-C) flow stress model parameters is presented using machining (cutting) force data. The results show fast calibration of the five J-C parameters within 25 simulations. In general, the BO-bound method is applicable for inverse calibration of any expensive simulation models as well as optimization problems with known bounds.

Bayesian optimization↗

Bi-fidelity approximation for uncertainty quantification and sensitivity analysis of irradiated particle-laden turbulence

Particle-laden turbulent flows subject to radiative heating are relevant in many applications, for example concentrated solar power receivers. Efficient and accurate simulations provide valuable insights and enable optimization of such systems. However, as there are many uncertainties inherent in such flows, uncertainty quantification is fundamental to improve the predictive capabilities of the numerical simulations. For large-scale, multi-physics problems exhibiting high-dimensional uncertainty, characterizing the stochastic solution presents a significant computational challenge as most strategies require a large number of high-fidelity solves. This requirement might result in an infeasible number of simulations when a typical converged high-fidelity simulation requires intensive computational resources. To reduce the cost of quantifying high-dimensional uncertainties, we investigate the application of a non-intrusive, bi-fidelity approximation to estimate statistics of quantities of interest associated with an irradiated particle-laden turbulent flow. This method exploits the low-rank structure of the solution to accelerate the stochastic sampling and approximation processes by means of cheaper-to-run, lower fidelity representations. Here, the application of this bi-fidelity approximation results in accurate estimates of the quantities of interest statistics, while requiring a small number of high-fidelity model evaluations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Optimising response to an introduction of African swine fever in wild pigs

African swine fever virus (ASFv) is a virulent pathogen that threatens domestic swine industries globally and persists in wild boar populations in some countries. Persistence in wild boar can challenge elimination and prevent disease-free status, making it necessary to address wild swine in proactive response plans. In the United States, invasive wild pigs are abundant and found across a wide range of ecological conditions that could drive different epidemiological dynamics among populations. Information on the size of the control areas required to rapidly eliminate the ASFv in wild pigs and how this area should change with management constraints and local ecology is needed to optimize response planning. We developed a spatially explicit disease transmission model contrasting wild pig movement and contact ecology in two ecosystems in Southeastern United States. We simulated ASFv spread and determined the optimal response area (reported as the radius of a circle) for eliminating ASFv rapidly over a range of detection times (when ASFv was detected relative to the true date of introduction), culling capacities (proportion of wild pigs in the culling zone removed weekly) and wild pig densities. Large radii for response areas (14 km) were needed under most conditions but could be shortened with early detection (≤ 8 weeks) and high culling capacities (≥ 15% weekly). Under most conditions, the ASFv was eliminated in less than 22 weeks using optimal control radii, although ecological conditions with high rates of wild pig movement required higher culling capacities (≥ 10% weekly) for elimination within 1 year. The results highlight the importance of adjusting response plans based on local ecology and show that wild pig movement is a better predictor of the optimal response area than the number of ASFv cases early in the outbreak trajectory. Furthermore, our framework provides a tool for determining optimal control plans in different areas, guiding expectations of response impacts, and planning resources needed for rapid elimination.

59 BASIC BIOLOGICAL SCIENCES↗

Multi-scale simulations of wind farm performance with complex terrain and weather events

Predicting the response of wind farms to local weather is necessary for optimal design and operation. In this work, simulation and analysis of a frontal passage through a utility scale wind farm is achieved for the first time using a seamless multi-scale modeling approach. A generalized actuator disk (GAD) wind turbine model is used to represent turbine wakes, and results are compared to novel radar observations during the frontal passage. The Weather Research and Forecasting (WRF) model is employed with a nested grid setup that allows for coupling between multi-scale atmospheric conditions and turbine response. Starting with mesoscale forcing, the atmosphere is dynamically downscaled to the region of interest, where the interaction between turbulent flows and individual wind turbines is simulated with the GAD at10 m resolution. Several improvements are made to the GAD model to mimic realistic turbine operation, including a yawing capability and a power output calculation. Ultimately, the model is able to capture both the dynamics of the frontal passage and the turbine response; predictions show good agreement with observed turbine wake structure, velocity deficit, and power output after accounting for a phase shift in the mesoscale forcing. This study demonstrates the utility of the WRF-GAD framework for simulating wind farm performance under complex atmospheric conditions.

17 WIND ENERGY↗

Numerical Investigation and Optimization of a Flushwall Injector for Scramjet Applications at Hypervelocity Flow Conditions

An investigation utilizing Reynolds-averaged simulations (RAS) was performed in order to find optimal designs for an interdigitated flushwall injector suitable for scramjet applications at hypervelocity conditions. The flight Mach number, duct height, spanwise width, and injection angle were the design variables selected to maximize two objective functions: the thrust potential and combustion efficiency. A Latin hypercube sampling design-of-experiments method was used to select design points for RAS. A methodology was developed that automated building geometries and generating grids for each design. The ensuing RAS analysis generated the performance database from which the two objective functions of interest were computed using a one-dimensional performance utility. The data were fitted using four surrogate models: an artificial neural network (ANN) model, a cubic polynomial, a quadratic polynomial, and a Kriging model. Variance-based decomposition showed that both objective functions were primarily driven by changes in the duct height. Multiobjective design optimization was performed for all four surrogate models via a genetic algorithm method. Optimal solutions were obtained at the upper and lower bounds of the flight Mach number range. The Kriging model obtained an optimal solution set that predicted high values for both objective functions. Additionally, three challenge points were selected to assess the designs on the Pareto fronts. Further sampling among the designs of the Pareto fronts are required in order to lower the errors and perform more accurate surrogate-based optimization. sed optimization.

Shenoy, Rajiv R.↗

Linear Regression Based Multi-fidelity Surrogate for Disturbance Amplification in Multi-phase Explosion

When simulations are very expensive and many are required, as for optimization or uncertainty quantification, a way to reduce cost is using surrogates. With multiple simulations to predict the quantity of interest, some being very expensive and accurate (high-fidelity simulations) and others cheaper but less accurate (low-fidelity simulations), it may be worthwhile to use multifidelity surrogates (MFSs). Moreover, if we can afford just a few high-fidelity simulations or experiments, MFS becomes necessary. Co-Kriging, which is probably the most popular MFS, replaces both low-fidelity and high-fidelity simulations by a single MFS. A recently proposed linear regression–based MFS (LR-MFS) offers the option to correct the LF simulations instead of correcting the LF surrogate in the MFS. When the low-fidelity simulation is cheap enough for use in an application, such as optimization, this may be an attractive option. Here, we explore the performance of LR-MFS using exact and surrogate-replaced low-fidelity simulations. The concern studied is a cylindrical dispersal of 100-μ m-diameter solid particles after detonation and the quantity of interest is a measure of the amplification of the departure from axisymmetry. We find very substantial accuracy improvements for this problem using the LR-MFS with exact low-fidelity simulations. Inspired by these results, we also compare the performance of co-Kriging to the use of Kriging to correct exact low-fidelity simulations and find a similar accuracy improvement when simulations are directly used. For this problem, further improvements in accuracy are achievable by taking advantage of inherent parametric symmetries. These findings may alert users of MFSs to the possible advantages of using exact low-fidelity simulations when this is affordable.

42 ENGINEERING↗

Development of high B s Fe–Ni-based metal amorphous nanocomposite by optimization of glass-forming ability

Thermocalc simulations identified compositions with good glass-forming ability (GFA) by locating minima in liquidus temperature, and solidification range in (Fe 70 Ni 30 ) x (B-Si-Nb) 100-x alloys with x = 82% and x = 85%, increased compared to previously developed x = 80% alloys. 3 compositions in the x = 82% system and 1 in the x = 85% system were successfully cast as amorphous magnetic ribbon. Magnetic properties of as-cast alloys showed improvements in Curie temperature and saturation induction, especially for the 85% alloy. Crystallization behavior in the 85% alloy was determined to occur as a 2-step process. Annealing studies showed (Fe 70 Ni 30 ) 82 Nb 2 Si 0 B 16 and (Fe 70 Ni 30 ) 82 Nb 3 Si 0 B 15 alloys to exhibit low coercivity and high B s with annealing. TEM results show significant crystallization in as-cast ribbon. The (Fe 70 Ni 30 ) 85 Nb 0.5 Si 0 B 14.5 alloy as-cast has small, uniformly distributed grains, allowing a useful nanocrystalline alloy to be produced as cast. Lastly, the parameter ΔT xg , the temperature range in which the material can be thermomechanically processed, was in the range of 12–23 °C.

36 MATERIALS SCIENCE↗

An approach to thermal-vacuum test optimization

An analytical model is presented that simulates failure rate characteristics of spacecraft equipment progressing through component and system level thermal-vacuum test programs and on into orbit. Costs for the test programs are part of the model as well as costs attributable to the mission. A lost value function combines the performance of the equipment with these costs and allows comparisons to be made among various test combinations and the resulting predicted mission performance. A minimum overall cost program is obtained by selecting that test program that yields the lowest lost value figure.

Kruger, R.↗

Implicit neural representations for experimental steering of advanced experiments

Scattering measurements using electrons, neutrons, or photons are essential for obtaining microscopic insights into materials. However, limited facility availability and high-dimensional scattering data necessitate more efficient experimental steering techniques. Here, we report a machine learning method that guides scattering data collection and facilitates real-time estimation of model parameters, given a reliable forward model to simulate experimental signals. We employ implicit neural representations as efficient surrogates that link model parameters with simulated spectroscopies. This enables a Bayesian optimal experimental design framework to estimate the probability distributions of parameters from high-dimensional scattering data. We demonstrate the proposed method using inelastic neutron scattering with simulated and real experimental data, highlighting the method’s ability to provide real-time parameter estimation with quantified uncertainties and to deliver informed experimental guidance that reduces experimental time while maximizing scientific output. This approach paves the way for accelerated discoveries in condensed matter through scattering measurements.

36 MATERIALS SCIENCE↗

Ultrafast Dynamics of Photoexcited Charge Carriers in In{sub 0.53}Ga{sub 0.47}As/In{sub 0.52}Al{sub 0.48}As Superlattices under Femtosecond Laser Excitation

The results of experimental studies of the time dynamics of photoexcited charge carriers in In{sub 0.53}Ga{sub 0.47}As/In{sub 0.52}Al{sub 0.48}As superlattices grown by molecular-beam epitaxy on a GaAs substrate with a metamorphic buffer are reported. On the basis of the results of the numerical simulation of band diagrams, the optimal thickness of the In{sub 0.52}Al{sub 0.48}As barrier layer (4 nm) is chosen. At this thickness, the electron wave functions in In{sub 0.53}Ga{sub 0.47}As substantially overlap the In{sub 0.52}Al{sub 0.48}As barriers. This makes it possible to attain a short lifetime of photoexcited charge carriers (τ ~ 3.4 ps) at the wavelength λ = 800 nm and the pumping power 50 mW without doping of the In{sub 0.53}Ga{sub 0.47}As layer with beryllium. It is shown that an increase in the wavelength to λ = 930 nm (at the same pumping power) yields a decrease in the lifetime of photoexcited charge carriers to τ ~ 2 ps. This effect is attributed to an increase in the capture cross section of trapping states for electrons with lower energies and to a decrease in the occupancy of traps at lower excitation densities.

36 MATERIALS SCIENCE↗

Examining Passenger Flow Choke Points at Airports Using Discrete Event Simulation

The movement of passengers through an airport quickly, safely, and efficiently is the main function of the various checkpoints (check-in, security. etc) found in airports. Human error combined with other breakdowns in the complex system of the airport can disrupt passenger flow through the airport leading to lengthy waiting times, missing luggage and missed flights. In this paper we present a model of passenger flow through an airport using discrete event simulation that will provide a closer look into the possible reasons for breakdowns and their implications for passenger flow. The simulation is based on data collected at Norfolk International Airport (ORF). The primary goal of this simulation is to present ways to optimize the work force to keep passenger flow smooth even during peak travel times and for emergency preparedness at ORF in case of adverse events. In this simulation we ran three different scenarios: real world, increased check-in stations, and multiple waiting lines. Increased check-in stations increased waiting time and instantaneous utilization. while the multiple waiting lines decreased both the waiting time and instantaneous utilization. This simulation was able to show how different changes affected the passenger flow through the airport.

Brown, Jeremy R.↗

Open-Loop Pitch Table Optimization for the Maximum Dynamic Pressure Orion Abort Flight Test

NASA has scheduled the retirement of the space shuttle orbiter fleet at the end of 2010. The Constellation program was created to develop the next generation of human spaceflight vehicles and launch vehicles, known as Orion and Ares respectively. The Orion vehicle is a return to the capsule configuration that was used in the Mercury, Gemini, and Apollo programs. This configuration allows for the inclusion of an abort system that safely removes the capsule from the booster in the event of a failure on launch. The Flight Test Office at NASA's Dryden Flight Research Center has been tasked with the flight testing of the abort system to ensure proper functionality and safety. The abort system will be tested in various scenarios to approximate the conditions encountered during an actual Orion launch. Every abort will have a closed-loop controller with an open-loop backup that will direct the vehicle during the abort. In order to provide the best fit for the desired total angle of attack profile with the open-loop pitch table, the table is tuned using simulated abort trajectories. A pitch table optimization program was created to tune the trajectories in an automated fashion. The program development was divided into three phases. Phase 1 used only the simulated nominal run to tune the open-loop pitch table. Phase 2 used the simulated nominal and three simulated off nominal runs to tune the open-loop pitch table. Phase 3 used the simulated nominal and sixteen simulated off nominal runs to tune the open-loop pitch table. The optimization program allowed for a quicker and more accurate fit to the desired profile as well as allowing for expanded resolution of the pitch table.

Stillwater, Ryan A.↗

Active Learning for Metamaterial Optimization on HPC and QC Integrated Systems

Active learning algorithms, integrating machine learning, quantum computing and optics simulation in an iterative loop, offer a promising approach to optimizing metamaterials. However, these algorithms can face difficulties in optimizing highly complex structures due to computational limitations. High-performance computing (HPC) and quantum computing (QC) integrated systems can address these issues by enabling parallel computing. In this study, we develop an active learning algorithm working on HPC-QC integrated systems. We evaluate the performance of optimization processes within active learning (i.e., training a machine learning model, problem-solving with quantum computing, and evaluating optical properties through wave-optics simulation) for highly complex metamaterial cases. Our results showcase that utilizing multiple cores on the integrated system can significantly reduce computational time, thereby enhancing the efficiency of optimization processes. Therefore, we expect that leveraging HPC-QC integrated systems helps effectively tackle large-scale optimization challenges in general.

Kim, Seongmin↗

Analysis of Vector Particle-In-Cell (VPIC) memory usage optimizations on cutting-edge computer architectures

Vector Particle-In-Cell (VPIC) is one of the fastest plasma simulation codes in the world, with particle numbers ranging from one trillion on the first petascale system, Roadrunner, to ten trillion particles on the more recent Blue Waters supercomputer. As supercomputers continue to grow rapidly in size, so too does the gap between computing capability and memory capability. Current memory systems limit VPIC simulations greatly as the maximum number of particles that can be simulated directly depends on the available memory. In this study, we present a suite of VPIC memory optimizations (i.e., particle weight, half-precision, and fixed-point optimizations) that enable a significant increase in the number of particles in VPIC simulations. Here, we assess the optimizations’ impact on memory and runtime performance for a suite of cutting-edge computer architectures such has the NVIDIA V100 GPU, the IBM Power9, and the Fujitsu A64FX architectures. Our optimizations enable a 31.25% reduction in memory usage and up to 40% increase in the number of particles. This paper extends our work on developing particle storage format optimizations Tan et al.

97 MATHEMATICS AND COMPUTING↗

Spatiotemporal quenches for efficient critical ground state preparation in the two-dimensional transverse field Ising model

Quantum simulators have the potential to shed light on the study of quantum many-body systems and materials, offering unique insights into various quantum phenomena. Although adiabatic evolution has been conventionally employed for state preparation, it faces challenges when the system evolves too quickly or the coherence time is limited. In such cases, shortcuts to adiabaticity, such as spatiotemporal quenches, provide a promising alternative. This paper numerically investigates the application of spatiotemporal quenches in the two-dimensional transverse field Ising model with ferromagnetic interactions, focusing on the emergence of the ground state and its correlation properties at criticality when the gap vanishes. We demonstrate the effectiveness of these quenches in rapidly preparing ground states in critical systems. Our simulations reveal the existence of an optimal quench front velocity at the emergent speed of light, leading to minimal excitation energy density and correlation lengths of the order of finite system sizes we can simulate. These findings emphasize the potential of spatiotemporal quenches for efficient ground state preparation in quantum systems, with implications for the exploration of strongly correlated phases and programmable quantum computing.

2-dimensional systems↗

Stratospheric limb radiance simulations for the HCl channel of CLAES

The earthlimb radiance of the HCl detector array of the CLAES instrument is simulated in order to specify and optimize the free spectral range for the etalon used in this spectral region and to select sampling positions which facilitate the recovery of HCl profiles. Considerations entering into the limb radiance computation of each contributing species are discussed. Simulations are then presented for tangent heights of 15 to 5 km for daytime and nighttime conditions, and their significance for HCl retrievals is addressed. The simulation results indicate that the HCl (P2) feature should be observable in the earth limb emission spectrum, and that an etalon with an FSR of 12.6/cm optimally avoids interference from OH emission lines blended with the HCl feature by adjacent transmission peaks.

Mergenthaler, J. L.↗