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 1,729 records · Page 96

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↗

Xyce Parallel Electronic Simulator Users' Guide Version 6.6.

This manual describes the use of the Xyce Parallel Electronic Simulator. Xyce has been de- signed as a SPICE-compatible, high-performance analog circuit simulator, and has been written to support the simulation needs of the Sandia National Laboratories electrical designers. This development has focused on improving capability over the current state-of-the-art in the following areas: Capability to solve extremely large circuit problems by supporting large-scale parallel com- puting platforms (up to thousands of processors). This includes support for most popular parallel and serial computers. A differential-algebraic-equation (DAE) formulation, which better isolates the device model package from solver algorithms. This allows one to develop new types of analysis without requiring the implementation of analysis-specific device models. Device models that are specifically tailored to meet Sandia's needs, including some radiation- aware devices (for Sandia users only). Object-oriented code design and implementation using modern coding practices. Xyce is a parallel code in the most general sense of the phrase -- a message passing parallel implementation -- which allows it to run efficiently a wide range of computing platforms. These include serial, shared-memory and distributed-memory parallel platforms. Attention has been paid to the specific nature of circuit-simulation problems to ensure that optimal parallel efficiency is achieved as the number of processors grows. The information herein is subject to change without notice. Copyright c 2002-2016 Sandia Corporation. All rights reserved. Acknowledgements The BSIM Group at the University of California, Berkeley developed the BSIM3, BSIM4, BSIM6, BSIM-CMG and BSIM-SOI models. The BSIM3 is Copyright c 1999, Regents of the University of California. The BSIM4 is Copyright c 2006, Regents of the University of California. The BSIM6 is Copyright c 2015, Regents of the University of California. The BSIM-CMG is Copyright c 2012 and 2016, Regents of the University of California. The BSIM-SOI is Copyright c 1990, Regents of the University of California. All rights reserved. The Mextram model has been developed by NXP Semiconductors until 2007, Delft University of Technology from 2007 to 2014, and Auburn University since April 2015. Copyrights c of Mextram are with Delft University of Technology, NXP Semiconductors and Auburn University. The MIT VS Model Research Group developed the MIT Virtual Source (MVS) model. Copyright c 2013 Massachusetts Institute of Technology (MIT). The EKV3 MOSFET model was developed by the EKV Team of the Electronics Laboratory-TUC of the Technical University of Crete. Trademarks Xyce TM Electronic Simulator and Xyce TM are trademarks of Sandia Corporation. Orcad, Orcad Capture, PSpice and Probe are registered trademarks of Cadence Design Systems, Inc. Microsoft, Windows and Windows 7 are registered trademarks of Microsoft Corporation. Medici, DaVinci and Taurus are registered trademarks of Synopsys Corporation. Amtec and TecPlot are trademarks of Amtec Engineering, Inc. All other trademarks are property of their respective owners. Contacts World Wide Web http://xyce.sandia.gov https://info.sandia.gov/xyce (Sandia only) Email xyce@sandia.gov (outside Sandia) xyce-sandia@sandia.gov (Sandia only) Bug Reports (Sandia only) http://joseki-vm.sandia.gov/bugzilla http://morannon.sandia.gov/bugzilla

47 OTHER INSTRUMENTATION↗

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↗

A GRACE ‐based assessment of interannual groundwater dynamics in the C ommunity L and M odel

Abstract The estimation of groundwater storage variations is important for quantifying available water resources and managing storage surpluses to alleviate storage deficiencies during droughts. This is particularly true in semi‐arid regions, where multiyear droughts can be common. To complement the local information provided by soil moisture and well level measurements, land models such as the Community Land Model (CLM) can be used to simulate regional scale water storage variations. CLM includes a bulk aquifer model to simulate saturated water storage dynamics below the model soil column. Aquifer storage increases when it receives recharge from the overlying soil column, and decreases due to lateral flow (i.e., base flow) and capillary rise. In this study, we examine the response of the CLM aquifer model to transitions between low and high recharge inputs, and show that the model simulates unrealistic long‐period behavior relative to total water storage (TWS) observations from the Gravity Recovery and Climate Experiment (GRACE). We attribute the model's poor response to large wetting events to the lack of a finite lower boundary in the bulk aquifer model. We show that by removing the bulk aquifer model and adding a zero‐flux boundary condition at the base of the soil column, good agreement with GRACE observations can be achieved. In addition, we examine the sensitivity of simulated total water storage to the depth at which the zero‐flux boundary is applied, i.e., the thickness of the soil column. Based on comparisons to GRACE, an optimal soil thickness map is constructed. Simulations using the modified CLM with the derived soil thickness map are shown to perform as well or better than standard CLM simulations. The improvements in simulated, climatically induced, long‐period water storage variability will reduce the uncertainty in GRACE‐based estimates of anthropogenic groundwater depletion.

Swenson, S. C.↗

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.↗

Cogging torque minimization in transverse flux machines

This study presents the design considerations in cogging torque minimization in two types of transverse flux machines. The machines have a double stator-single rotor configuration with flux concentrating ferrite magnets. One of the machines has pole windings across each leg of an E-Core stator. Another machine has quasi-U-shaped stator cores and a ring winding. The flux in the stator back iron is transverse in both machines. Different methods of cogging torque minimization are investigated. Key methods of cogging torque minimization are identified and used as design variables for optimization using a design of experiments (DOE) based on the Taguchi method. A multi-level DOE is proposed as an optimization method to reach an optimum solution with minimum simulations. The case study is analyzed in a two level DOE optimization. Finite element analysis is used to study the different effects. Two prototypes are fabricated for validating the FEA results.

24 POWER TRANSMISSION AND DISTRIBUTION↗

A singular perturbation analysis of minimum time long range intercept

The application of singular perturbation theory for deriving algorithms suitable for on board real time computation of optimal aircraft trajectory control was investigated. Minimum time intercept in three dimensions was selected as a pilot problem formulation, and data for an early version F-4 aircraft was used to represent aircraft aerodynamic and propulsion characteristics. Results include the derivation of nonlinear, near optimal feedback control laws that are readily implementable in an aircraft flight computer, performance results for a variety of initial conditions (including optimal intercept at short ranges), and a point mass three dimensional simulation suitable for representing fighter aircraft under closed loop optimal control.

Calise, A.↗

Improved optical standing-wave beam splitters for dilute Bose–Einstein condensates

Bose–Einstein condensate (BEC)-based atom interferometry exploits low temperatures and long coherence lengths to facilitate high-precision measurements. Progress in atom interferometry promises improvements in navigational devices like gyroscopes and accelerometers, as well as applications in fundamental physics such as accurate determination of physical constants. Previous work demonstrates that beam splitters and mirrors for coherent manipulation of dilute BEC momentum in atom interferometers can be implemented with sequences of non-resonant standing-wave light pulses. While previous work focuses on the optimization of the optical pulses’ amplitude and duration to produce high-order momentum states with high fidelity, we explore how varying the shape of the optical pulses affects optimal beam-splitter performance, as well as the effect of pulse shape on the sensitivity of optimized parameters in achieving high fidelity in high-momentum states. Here, in simulations of two-pulse beam splitters utilizing optimized square, triangle, and sinc-squared pulse shapes applied to dilute BECs, we, in some cases, reduce parameter sensitivity by an order of magnitude while maintaining fidelity.

74 ATOMIC AND MOLECULAR PHYSICS↗