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,621 records · Page 90
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

Exploring an Ensemble-Based Approach to Atmospheric Climate Modeling and Testing at Scale

A strict throughput requirement has placed a cap on the degree to which we can depend on the execution of single, long, fine spatial grid simulations to explore global atmospheric climate behavior. Alternatively, running an ensemble of short simulations is computationally more efficient. We test the null hypothesis that the climate statistics of a full-complexity atmospheric model derived from an ensemble of independent short simulation is equivalent to that from an equilibrated long simulation. The climate of short simulation ensembles is statistically distinguishable from that of a long simulation in terms of the distribution of global annual means, largely due to the presence of low-frequency atmospheric intrinsic variability in the long simulation. We also find that model climate statistics of the simulation ensemble are sensitive to the choice of compiler optimizations. While some answer-changing optimization choices do not effect the climate state in terms of mean, variability and extremes, aggressive optimizations can result in significantly different climate states.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2016

TU-AB-BRC-12: Optimized Parallel MonteCarlo Dose Calculations for Secondary MU Checks

Purpose: Secondary MU checks are an important tool used during a physics review of a treatment plan. Commercial software packages offer varying degrees of theoretical dose calculation accuracy, depending on the modality involved. Dose calculations of VMAT plans are especially prone to error due to the large approximations involved. Monte Carlo (MC) methods are not commonly used due to their long run times. We investigated two methods to increase the computational efficiency of MC dose simulations with the BEAMnrc code. Distributed computing resources, along with optimized code compilation, will allow for accurate and efficient VMAT dose calculations. Methods: The BEAMnrc package was installed on a high performance computing cluster accessible to our clinic. MATLAB and PYTHON scripts were developed to convert a clinical VMAT DICOM plan into BEAMnrc input files. The BEAMnrc installation was optimized by running the VMAT simulations through profiling tools which indicated the behavior of the constituent routines in the code, e.g. the bremsstrahlung splitting routine, and the specified random number generator. This information aided in determining the most efficient compiling parallel configuration for the specific CPU’s available on our cluster, resulting in the fastest VMAT simulation times. Our method was evaluated with calculations involving 10{sup 8} – 10{sup 9} particle histories which are sufficient to verify patient dose using VMAT. Results: Parallelization allowed the calculation of patient dose on the order of 10 – 15 hours with 100 parallel jobs. Due to the compiler optimization process, further speed increases of 23% were achieved when compared with the open-source compiler BEAMnrc packages. Conclusion: Analysis of the BEAMnrc code allowed us to optimize the compiler configuration for VMAT dose calculations. In future work, the optimized MC code, in conjunction with the parallel processing capabilities of BEAMnrc, will be applied to provide accurate and efficient secondary MU checks.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2014

Clinical Implementation of Intensity Modulated Proton Therapy for Thoracic Malignancies

Purpose: Intensity modulated proton therapy (IMPT) can improve dose conformality and better spare normal tissue over passive scattering techniques, but range uncertainties complicate its use, particularly for moving targets. We report our early experience with IMPT for thoracic malignancies in terms of motion analysis and management, plan optimization and robustness, and quality assurance. Methods and Materials: Thirty-four consecutive patients with lung/mediastinal cancers received IMPT to a median 66 Gy(relative biological equivalence [RBE]). All patients were able to undergo definitive radiation therapy. IMPT was used when the treating physician judged that IMPT conferred a dosimetric advantage; all patients had minimal tumor motion (<5 mm) and underwent individualized tumor-motion dose-uncertainty analysis and 4-dimensional (4D) computed tomographic (CT)-based treatment simulation and motion analysis. Plan robustness was optimized by using a worst-case scenario method. All patients had 4D CT repeated simulation during treatment. Results: IMPT produced lower mean lung dose (MLD), lung V{sub 5} and V{sub 20}, heart V{sub 40}, and esophageal V{sub 60} than did IMRT (P<.05) and lower MLD, lung V{sub 20}, and esophageal V{sub 60} than did passive scattering proton therapy (PSPT) (P<.05). D{sub 5} to the gross tumor volume and clinical target volume was higher with IMPT than with intensity modulated radiation therapy or PSPT (P<.05). All cases were analyzed for beam-angle-specific motion, water-equivalent thickness, and robustness. Beam angles were chosen to minimize the effect of respiratory motion and avoid previously treated regions, and the maximum deviation from the nominal dose-volume histogram values was kept at <5% for the target dose and met the normal tissue constraints under a worst-case scenario. Patient-specific quality assurance measurements showed that a median 99% (range, 95% to 100%) of the pixels met the 3% dose/3 mm distance criteria for the γ index. Adaptive replanning was used for 9 patients (26.5%). Conclusions: IMPT using 4D CT-based planning, motion management, and optimization was implemented successfully and met our quality assurance parameters for treating challenging thoracic cancers.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2024

Liquid Crystal Orientation and Shape Optimization for the Active Response of Liquid Crystal Elastomers

Liquid crystal elastomers (LCEs) are responsive materials that can undergo large reversible deformations upon exposure to external stimuli, such as electrical and thermal fields. Controlling the alignment of their liquid crystals mesogens to achieve desired shape changes unlocks a new design paradigm that is unavailable when using traditional materials. While experimental measurements can provide valuable insights into their behavior, computational analysis is essential to exploit their full potential. Accurate simulation is not, however, the end goal; rather, it is the means to achieve their optimal design. Such design optimization problems are best solved with algorithms that require gradients, i.e., sensitivities, of the cost and constraint functions with respect to the design parameters, to efficiently traverse the design space. In this work, a nonlinear LCE model and adjoint sensitivity analysis are implemented in a scalable and flexible finite element-based open source framework and integrated into a gradient-based design optimization tool. To display the versatility of the computational framework, LCE design problems that optimize both the material, i.e., liquid crystal orientation, and structural shape to reach a target actuated shapes or maximize energy absorption are solved. Multiple parameterizations, customized to address fabrication limitations, are investigated in both 2D and 3D. The case studies are followed by a discussion on the simulation and design optimization hurdles, as well as potential avenues for improving the robustness of similar computational frameworks for applications of interest.

42 ENGINEERING↗
DOE OSTI2013

Monte Carlo study for optimal conditions in single-shot imaging with femtosecond x-ray laser pulses

Intense x-ray pulses from x-ray free electron lasers (XFELs) enable the unveiling of atomic structure in material and biological specimens via ultrafast single-shot exposures. As the radiation is intense enough to destroy the sample, a new sample must be provided for each x-ray pulse. These single-particle delivery schemes require careful optimization, though systematic study to find such optimal conditions is still lacking. We have investigated two major single-particle delivery methods: particle injection as flying objects and membrane-mount as fixed targets. The optimal experimental parameters were searched for via Monte Carlo simulations to discover that the maximum single-particle hit rate achievable is close to 40%.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Cloud-resolving model intercomparison of an MC3E squall line case: Part I-Convective updrafts: CRM Intercomparison of a Squall Line

A constrained model intercomparison study of a mid-latitude mesoscale squall line is performed using the Weather Research & Forecasting (WRF) model at 1-km horizontal grid spacing with eight cloud microphysics schemes, to understand specific processes that lead to the large spread of simulated cloud and precipitation at cloud-resolving scales, with a focus of this paper on convective cores. Various observational data are employed to evaluate the baseline simulations. All simulations tend to produce a wider convective area than observed, but a much narrower stratiform area, with most bulk schemes overpredicting radar reflectivity. The magnitudes of the virtual potential temperature drop, pressure rise, and the peak wind speed associated with the passage of the gust front are significantly smaller compared with the observations, suggesting simulated cool pools are weaker. Simulations also overestimate the vertical velocity and Ze in convective cores as compared with observational retrievals. The modeled updraft velocity and precipitation have a significant spread across the eight schemes even in this strongly dynamically-driven system. The spread of updraft velocity is attributed to the combined effects of the low-level perturbation pressure gradient determined by cold pool intensity and buoyancy that is not necessarily well correlated to differences in latent heating among the simulations. Variability of updraft velocity between schemes is also related to differences in ice-related parameterizations, whereas precipitation variability increases in no-ice simulations because of scheme differences in collision-coalescence parameterizations.

Fan, Jiwen↗
DOE OSTI2014

Development of electromagnetic acoustic transducer (EMAT) phased arrays for SFR inspection

A long-standing problem for Sodium cooled Fast Reactor (SFR) instrumentation is the development of efficient under-sodium visualization systems adapted to the hot and opaque sodium environment. Electromagnetic Acoustic Transducers (EMAT) are potential candidates for a new generation of Ultrasonic Testing (UT) probes well-suited for SFR inspection that can overcome drawbacks of classical piezoelectric probes in sodium environment. Based on the use of new CIVA simulation tools, we have designed and optimized an advanced EMAT probe for under-sodium visualization. This has led to the development of a fully functional L-wave EMAT sensing system composed of 8 elements and a casing withstanding 200° C sodium inspection. Laboratory experiments demonstrated the probe's ability to sweep an ultrasonic beam to an angle of 15 degrees. Testing in a specialized sodium facility has shown that it was possible to obtain pulse-echo signals from a target under several different angles from a fixed position.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Adjoint-optimized nanoscale light extractor for nitrogen-vacancy centers in diamond

Abstract We designed a nanoscale light extractor (NLE) for the efficient outcoupling and beaming of broadband light emitted by shallow, negatively charged nitrogen-vacancy (NV) centers in bulk diamond. The NLE consists of a patterned silicon layer on diamond and requires no etching of the diamond surface. Our design process is based on adjoint optimization using broadband time-domain simulations and yields structures that are inherently robust to positioning and fabrication errors. Our NLE functions like a transmission antenna for the NV center, enhancing the optical power extracted from an NV center positioned 10 nm below the diamond surface by a factor of more than 35, and beaming the light into a ±30° cone in the far field. This approach to light extraction can be readily adapted to other solid-state color centers.

42 ENGINEERING↗
DOE OSTI2019

Designing Critical Experiments Using Gaussian Process Optimization

New nuclear applications may lack similarity between the application and existing ICSBEP benchmarks. Without accurate benchmark data, predictive simulations will have high uncertainty in the estimate of a parameter of interest (such as k eff ). The design of critical experiments tuned for specific applications of interest is one solution to reducing uncertainty in the necessarily nuclear data. We propose an optimization technique to design critical experiments that match the cross-section sensitivities of an application model while minimizing the number of simulation runs in the optimization. The optimization of a critical experiment is slow because the potential design space covers many qualitative variables (types of material and geometries) and corresponding quantitative variables (e.g., material dimensions) and the evaluation of any specific design requires running an expensive code like MCNP1 and Whisper. The optimization of expensive black-box computer codes with quantitative variables has been extensively explored in the literature of design of computer experiments. We propose a modified version of standard Gaussian process (GP) optimization that optimizes over multiple discrete configurations of qualitative variables in parallel which reduces the computational cost of the entire optimization. The goal of this work is to minimize the number of MCNP evaluations while still exploring many different design configurations in parallel.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Incorporating Valve Switching Losses Into a Static Optimal Control Algorithm for the Hybrid Hydraulic-Electric Architecture (HHEA)

The Hybrid Hydraulic-Electric Architecture (HHEA) has, in recent years, been proposed as an energy efficient alternative to conventional load-sensing architectures in mobile machines such as excavators and wheel-loaders. HHEA leverages the advantages of hydraulic power and electric power to eliminate throttling valves while also improving the energy and control performance of the system. The architecture utilizes a set of common pressure rails to provide a majority of power and and a small electric motor driven pump to modulate this power to meet the exact demand. Previous work has developed a computationally efficient Lagrange Multiplier approach for determining the optimal pressure rail selections that minimizes the energy losses in the system. The static model used considers only the energy use for each pressure rail selection but not the losses associated with the valves during the transition. This paper presents an approach to include the switching losses in the model and in the optimization procedure. To capture the switching losses, switching events between different rails and at various input and output flow rates were simulated with consideration of valve spool dynamics. A parameterized model that summarizes the losses is then obtained, allowing switching losses to be added to the previous energy analysis. The performance of the switching loss model was compared with reference data obtained from a high-fidelity simulation model. To incorporate the switching losses into optimal control algorithm, an efficient dynamic programming approach that prevents frequent switching is adopted in place of the Lagrange multiplier approach. The overall effect of switching losses on the energy consumption and optimal control decisions is presented. In general, switching losses contribute to about 9–10% of input energy.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2021

Topology Optimization with a Manufacturability Objective

Part distortion and residual stress are critical factors for metal additive manufacturing (AM) because they can lead to high failure rates during both manufacturing and service. We present a topology optimization approach that incorporates a fast AM process simulation at each design iteration to provide predictions of manufacturing outcomes (i.e., residual stress, distortion, residual elastic energy) that can be optimized or constrained. The details of the approach and implementation are discussed, and an example design is presented that illustrates the efficacy of the method.

42 ENGINEERING↗
DOE OSTI2019

Mixing state evolution of agglomerating particles in an aerosol chamber: Comparison of measurements and particle-resolved simulations

This article introduces a validation study of the stochastic particle-resolved aerosol model PartMC with experimental data from an aerosol chamber experiment. For the experiment, a scanning mobility particle sizer and a single-particle soot photometer were used to monitor the aerosol mixing state evolution of two initially externally-mixed aerosol populations of ammonium sulfate and black carbon particles undergoing agglomeration. We applied an efficient optimization algorithm (ProSRS) to determine several unconstrained simulation parameters and were able to successfully reproduce number concentrations and size distributions of mixed particles that formed by agglomeration. The PartMC modeling method in conjunction with the optimization procedure provides a tool for detailed comparisons of chamber experiments and modeling, where aerosol mixing state is the focus of investigation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2016

Analysis of the effects of mismatched errors on coherent beam combining based on a self-imaging waveguide

A theoretical model of coherent beam combining (CBC) based on a self-imaging waveguide (SIW) is built and the effects of mismatched errors on SIW-based CBC are simulated and analysed numerically. With the combination of the theoretical model and the finite difference beam propagation method, two main categories of errors, assembly and nonassembly errors, are numerically studied to investigate their effect on the beam quality by using the M{sup 2} factor. The optimisation of the SIW and error control principle of the system is briefly discussed. The generalised methodology offers a good reference for investigating waveguide-based high-power coherent combining of fibre lasers in a comprehensive way. (lasers and laser beams)

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Coordinated Optimization of Multiservice Dispatch for Energy Storage Systems with Degradation Model for Utility Applications

The rapid development of energy storage technologies is motivating researchers to understand the possible value of energy storage. This paper evaluates the multiservice dispatch of distribution-sited energy storage systems with realistic degradation models by integrating a multiservice optimization model with a distribution system simulator. This work attempts to define the theoretical value energy storage can bring when multiple value streams are considered in the modeling environment. Use cases presented in this paper deliver multiple services to both distribution and transmission stakeholders that can provide positive net value to the electricity system under prevailing energy storage cost structures. Furthermore, this work was performed in collaboration with a distribution utility to help utility personnel understand the value distributed energy storage can provide by examining different utility use cases without considering implementation costs.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2020

An open source fast fluid dynamics model for data center thermal management

Although computational fluid dynamics (CFD) has been widely adopted to improve data center thermal management, the high computational demand limits its applications, such as multivariate optimal design and operation. Fast fluid dynamics (FFD), which has been applied for fast airflow simulation, shows great potential. However, few research applied FFD for optimal design and operation of data center thermal management. This research improves the FFD model for data centers and conducts a comprehensive evaluation and demonstration. First, the FFD model is improved by solving the advection and diffusion equations together using an upwind scheme instead of a semi-Lagrangian advection solver in the conventional FFD model. Second, new features for data centers are added, such as a pressure correction method to simulate plenum airflow and dynamic boundary conditions for IT racks. The new FFD model is first validated with two indoor environment cases and the results show that the new FFD model has slightly better overall prediction accuracy and faster speed compared to the conventional FFD model. It is also observed that both FFD models achieve acceptable accuracy, except for a few localized disparities with experimental data, which might be due to simplified handling of turbulence viscosity near the boundaries. Furthermore, validation with a real data center shows that the FFD model achieves a similar level of accuracy as CFD when compared to the experimental measurements with some level of uncertainties. It is then demonstrated for data center optimal design and operation, which saves 53.4–58.8% of annual energy while still meeting the thermal requirements. In conclusion, with a much faster speed and comparable accuracy compared to CFD, the FFD model parallelized on a graphics processing unit is promising for practical model-based data center early design and operation.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2019

Designing Critical Experiments using Gaussian Process Optimization [Slides]

When designing new nuclear applications, there may be a lack of similarity between the application and existing ICSBEP benchmarks. Without appropriate benchmark data, predictive simulations will have high uncertainty in parameters of interest, such as K eff . The goal of the project is to create an algorithm to optimally design critical experiments that match cross-section sensitivities of an application model while minimizing the number of simulation runs in the optimization.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Minimizing thickness variation in monolithic U-10Mo fuel foil and Zr interlayer during hot rolling: A microstructure-based finite element method analysis

Low-enriched uranium alloyed with 10 wt. % molybdenum (U-10Mo) has been identified as a promising alternative to highly enriched uranium fuel for the United States’ high performance research reactors. The monolithic U-10Mo fuel plate consists of a metallic U-10Mo fuel foil with a 25 µm Zr interlayer and a relatively thick cladding of aluminum alloy 6061. The Zr interlayer is typically applied during the hot co-rolling process, and this process dictates the uniformity of the Zr interlayer. Thickness variation observed in the U-10Mo and Zr interlayer has been attributed to several sources: the initial grain size of the U-10Mo castings, can materials, rolling temperature, inhomogeneous molybdenum content, and porosity in the cast U-10Mo. This thickness variation limits the ability to meet the dimensional specification; thus, a better understanding of the factors causing the nonuniform thickness is needed. In this work, we used a novel, microstructure-based finite element method to model the hot rolling process to address these concerns. Grain microstructures in U-10Mo were tessellated and explicitly considered in the finite element model. Each grain was assigned a random material property to mimic the grain strength variations induced by different grain orientations. Simulations were performed using six steel can thicknesses, four grain sizes, and with or without a Zr interlayer to investigate the influences of those variables on the thickness nonuniformity. The simulation results showed that a thinner steel can and finer U-10Mo grain size reduce thickness variations in both the U-10Mo fuel foil and Zr interlayer. The direct findings from the simulations and analysis can be used to optimize the hot rolling schedule, reduce fabrication defects, and meet the dimensional specifications. The proposed microstructure-based finite element model can be also coupled with experimental microstructure characterization data, images, and models to simulate multi-pass hot rolling.

36 MATERIALS SCIENCE↗