Optimal placement of active/passive members in truss structures using simulated annealing
Explore the source record for details and available documents.
Chen, Gun-Shing↗
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.
Transverse Mode Dynamics of VCSELs Through Space-Time Domain Simulation
Modeling and simulation are important to understand laser operation and to optimize and design device functions. Numerical simulation of VCSEL (Vertical Cavity Surface Emitting Lasers) has been largely based on solving time-independent Helmholtz equation or time dependent coupled mode equations. There are various advantages for choosing these approaches. However, the disadvantages are also apparent. The former cannot handle dynamical mode competition seen in VCSELs, while the latter assumes a given type and number of modes a priori. Furthermore, the microscopic physics of heterstructures and electron-hole plasma is very often represented by a few parameters such as linear gain coefficients and the linewidth enhancement factor. These are over simplification of space and frequency (wavelength) dependent gain and refractive index functions. When the space-time dynamical operation of VCSELs becomes important, these simple approximations become questionable. In this paper, we apply a recently developed model for edge-emitting lasers to a gain guided VCSEL for space-time domain simulation. This model takes into account the actual nonlinear dependence of gain and refractive index on frequency and carrier density within the frame work of the effective Bloch equations. The corresponding partial differential equations are solved directly by finite difference methods. Laser behavior with increasing pumping current is investigated in detail. Special attention is paid to the dynamical competition of the transverse modes.
Goorjian, Peter M.↗
Range Process Simulation Tool
Range Process Simulation Tool (RPST) is a computer program that assists managers in rapidly predicting and quantitatively assessing the operational effects of proposed technological additions to, and/or upgrades of, complex facilities and engineering systems such as the Eastern Test Range. Originally designed for application to space transportation systems, RPST is also suitable for assessing effects of proposed changes in industrial facilities and large organizations. RPST follows a model-based approach that includes finite-capacity schedule analysis and discrete-event process simulation. A component-based, scalable, open architecture makes RPST easily and rapidly tailorable for diverse applications. Specific RPST functions include: (1) definition of analysis objectives and performance metrics; (2) selection of process templates from a processtemplate library; (3) configuration of process models for detailed simulation and schedule analysis; (4) design of operations- analysis experiments; (5) schedule and simulation-based process analysis; and (6) optimization of performance by use of genetic algorithms and simulated annealing. The main benefits afforded by RPST are provision of information that can be used to reduce costs of operation and maintenance, and the capability for affordable, accurate, and reliable prediction and exploration of the consequences of many alternative proposed decisions.
Phillips, Dave↗
Evaluation of Surface Crack Formation in Photovoltaic Backsheets Using Fragmentation and Finite Element Simulations: Preprint
Backsheet cracking is among the most commonly observed degradation modes of photovoltaic (PV) modules in the field. Cracks can reduce the ability of backsheets to fulfil their functions, for example, protection of the modules from the environment or electrical insulation. This work presents an evaluation of the degradation and cracking propensity of two backsheets during accelerated IEC TS 62788-7-2 A3 (International Electrotechnical Commission) aging with the fragmentation test: a co-extruded polyamide backsheet (AAA), and a laminated multilayer backsheet with a polyethylene terephthalate core and outer layer and ethylene vinyl acetate inner layer (PPE). Results show the surface embrittlement of the AAA outer layer during exposure. A longer time of exposure causes cracks to form at lower strains during stretching and creates deeper cracks. Accordingly, Young's modulus of the outer layer increases, as measured by cross-sectional nanoindentation. PPE exhibits cracking after exposure as well. While no similar increase of modulus or crack depth can be observed, the outer layer of PPE exhibits more obvious signs of erosion during exposure, including progressive morphological changes and thickness losses. A finite element model to simulate surface crack formation is devised, based on initially zero-thickness decohesion elements. Decohesion criteria define a critical stress, at which these elements grow, and cracks can begin to form. These criteria are obtained via parameter optimization by comparison between simulation and experiment. The model is used to interpret crack formation in both backsheets.
backsheet↗
Bayesian optimization scheme for the design of a nanofibrous high power target
High Power Targetry (HPT) R&D is critical in the context of increasing beam intensity and energy for next generation accelerators. Many target concepts and novel materials are being developed and tested for their ability to withstand extreme beam environments; the HPT R&D Group at Fermilab is developing an electrospun nanofiber material for this purpose. The performance of these nanofiber targets is sensitive to their construction parameters, such as the packing density of the fibers. Lowering the density improves the survival of the target, but reduces the secondary particle yield. Optimizing the lifetime and production efficiency of the target poses an interesting design problem, and in this paper we study the applicability of Bayesian optimization to its solution. We first describe how to encode the nanofiber target design problem as the optimization of an objective function, and how to evaluate that function with computer simulations. We then explain the optimization loop setup. Thereafter, we present the optimal design parameters suggested by the algorithm, and close with discussions of limitations and future refinements.
43 PARTICLE ACCELERATORS↗
Bayesian Optimization Scheme for the Design of a Nanofibrous High Power Target
High Power Targetry (HPT) R\&D is critical in the context of increasing beam intensity and energy for next generation accelerators. Many target concepts and novel materials are being developed and tested for their ability to withstand extreme beam environments; the HPT R\&D Group at Fermilab is developing an electrospun nanofiber material for this purpose. The performance of these nanofiber targets is sensitive to their construction parameters, such as the packing density of the fibers. Lowering the density improves the survival of the target, but reduces the secondary particle yield. Optimizing the lifetime and production efficiency of the target poses an interesting design problem, and in this paper we study the applicability of Bayesian optimization to its solution. We first describe how to encode the nanofiber target design problem as the optimization of an objective function, and how to evaluate that function with computer simulations. We then explain the optimization loop setup. Thereafter, we present the optimal design parameters suggested by the algorithm, and close with discussions of limitations and future refinements.
43 PARTICLE ACCELERATORS↗
Development of Multiobjective Optimization Techniques for Sonic Boom Minimization
A discrete, semi-analytical sensitivity analysis procedure has been developed for calculating aerodynamic design sensitivities. The sensitivities of the flow variables and the grid coordinates are numerically calculated using direct differentiation of the respective discretized governing equations. The sensitivity analysis techniques are adapted within a parabolized Navier Stokes equations solver. Aerodynamic design sensitivities for high speed wing-body configurations are calculated using the semi-analytical sensitivity analysis procedures. Representative results obtained compare well with those obtained using the finite difference approach and establish the computational efficiency and accuracy of the semi-analytical procedures. Multidisciplinary design optimization procedures have been developed for aerospace applications namely, gas turbine blades and high speed wing-body configurations. In complex applications, the coupled optimization problems are decomposed into sublevels using multilevel decomposition techniques. In cases with multiple objective functions, formal multiobjective formulation such as the Kreisselmeier-Steinhauser function approach and the modified global criteria approach have been used. Nonlinear programming techniques for continuous design variables and a hybrid optimization technique, based on a simulated annealing algorithm, for discrete design variables have been used for solving the optimization problems. The optimization procedure for gas turbine blades improves the aerodynamic and heat transfer characteristics of the blades. The two-dimensional, blade-to-blade aerodynamic analysis is performed using a panel code. The blade heat transfer analysis is performed using an in-house developed finite element procedure. The optimization procedure yields blade shapes with significantly improved velocity and temperature distributions. The multidisciplinary design optimization procedures for high speed wing-body configurations simultaneously improve the aerodynamic, the sonic boom and the structural characteristics of the aircraft. The flow solution is obtained using a comprehensive parabolized Navier Stokes solver. Sonic boom analysis is performed using an extrapolation procedure. The aircraft wing load carrying member is modeled as either an isotropic or a composite box beam. The isotropic box beam is analyzed using thin wall theory. The composite box beam is analyzed using a finite element procedure. The developed optimization procedures yield significant improvements in all the performance criteria and provide interesting design trade-offs. The semi-analytical sensitivity analysis techniques offer significant computational savings and allow the use of comprehensive analysis procedures within design optimization studies.
Chattopadhyay, Aditi↗
A Method for Combining Experimentation and Molecular Dynamics Simulation to Improve Cohesive Zone Models for Metallic Microstructures
Fracture processes within a material begin at the nanometer length scale at which the formation, propagation, and interaction of fundamental damage mechanisms occur. Physics-based modeling of these atomic processes quickly becomes computationally intractable as the system size increases. Thus, a multiscale modeling method, based on the aggregation of fundamental damage processes occurring at the nanoscale within a cohesive zone model, is under development and will enable computationally feasible and physically meaningful microscale fracture simulation in polycrystalline metals. This method employs atomistic simulation to provide an optimization loop with an initial prediction of a cohesive zone model (CZM). This initial CZM is then applied at the crack front region within a finite element model. The optimization procedure iterates upon the CZM until the finite element model acceptably reproduces the near-crack-front displacement fields obtained from experimental observation. With this approach, a comparison can be made between the original CZM predicted by atomistic simulation and the converged CZM that is based on experimental observation. Comparison of the two CZMs gives insight into how atomistic simulation scales.
Hochhalter, J. D.↗
Optimizing Asynchronous Multi-Level Checkpoint/Restart Configurations with Machine Learning
With the emergence of versatile storage systems, multi-level checkpointing (MLC) has become a common approach to gain efficiency. However, multi-level checkpoint/restart can cause enormous I/O traffic on HPC systems. To use multilevel checkpointing efficiently, it is important to optimize check-point/restart configurations. Current approaches, namely modeling and simulation, are either inaccurate or slow in determining the optimal configuration for a large scale system. In this paper, we show that machine learning models can be used in combination with accurate simulation to determine the optimal checkpoint configurations. We also demonstrate that more advanced techniques such as neural networks can further improve the performance in optimizing checkpoint configurations.
checkpoint-restart↗
Validation and optimization of batch and continuous particle separation processes
Recent advances in simulation capabilities coupled with advanced manufacturing can be leveraged to create highly efficient compact reactors with low emissions, low cost, high reliability, and high flexibility for energy conversion. This type of reactor is especially important for conversion of bio-mass, a storable, renewable and fast-growing energy source, considering the distributed nature of raw material and cost of biomass collection and transportation. To achieve this goal, it is critical to developing new process technologies that are significantly more efficient than traditional processes for which the added manufacture provides a fast way of realization and evaluation. Based upon the flagship MFiX multiphase computational fluid dynamic software suite, NETL has developed an Optimization Toolset to facilitate the CFD simulation based reactor design and optimization, and process intensification. This technical report presents the recent development and application of the Optimization Toolset and the framework for process optimization of particle separation in batch and continuous processes. In these processes, a fluidized bed reactor under proper operating condition was used to separate particles of different physical properties. MFiX-DEM simulations were conducted for a batch system to simulate the process covering a wide range of superficial gas velocities from which the optimal operational parameter can be determined through analyzing the system response. After successful demonstration of the framework, a continuous system for particle separation was designed based on CFD simulations. Similar optimization exercises were conducted for the continuous particle separation considering multiple operating parameters. The design was then fabricated in the laboratory through 3D printing for testing. The particle separation performance was carefully measured through high-speed video for particle tracking. Quantitative comparison between numerical results and experimental measurements were conducted for confirmation and validation.
30 DIRECT ENERGY CONVERSION↗
Mitigating the noise of DESI mocks using analytic control variates
In order to address fundamental questions related to the expansion history of the Universe and its primordial nature with the next generation of galaxy experiments, we need to model reliably large-scale structure observables such as the correlation function and the power spectrum. Cosmological N-body simulations provide a reference through which we can test our models, but their output suffers from sample variance on large scales. Fortunately, this is the regime where accurate analytic approximations exist. To reduce the variance, which is key to making optimal use of these simulations, we can leverage the accuracy and precision of such analytic descriptions using Control Variates (CV). The power of control variates stems from utilizing inexpensive but highly correlated surrogates of the statistics one wishes to measure. The stronger the correlation between the surrogate and the statistic of interest, the larger the variance reduction delivered by the method. We apply two control variate formulations to mock catalogs generated in anticipation of upcoming data from the Dark Energy Spectroscopic Instrument (DESI) to test the robustness of its analysis pipeline. Our CV-reduced measurements offer a factor of 5-10 improvement in the measurement error compared with the raw measurements. We explore the relevant properties of the galaxy samples that dictate this reduction and comment on the improvements we find on some of the derived quantities relevant to Baryon Acoustic Oscillation (BAO) analysis.
79 ASTRONOMY AND ASTROPHYSICS↗
HPC Network Simulation Tuning via Automatic Extraction of Hardware Parameters
Popular HPC network interconnection simulators such as SST/macro provide a variety of configurable parameters to explore the design space of hardware components such as network interface cards (NIC), switches, and links among them. While such knobs provide flexibility to explore design trade-offs for novel hardware, manually configuring simulations for matching configurations of the existing hardware to focus on topology exploration can be cumbersome and error-prone, leading to widely inaccurate simulations. This challenge is compounded when specifications of various (proprietary) technologies are not readily available or intentionally omitted. In this work, we propose a framework to autotune the multiple network models’ simulation configurations within SST/macro using Tree-structured Parzen Estimator-based Bayesian optimization to observe the effect on simulation accuracy across different message regimes. These regimes consist of small to large message sizes and latency to bandwidth-bound messages. We provide a detailed analysis of the simulation error for four representative HPC systems. Our Bayesian optimization based autotuning framework for network models achieves a maximum of 5x improvement in accuracy over best-effort manual configurations based on available hardware specifications.
Simulation, autotuning↗
Optimization of Structural Flood Mitigation Strategies
The dynamics of flooding are primarily influenced by the shape, height, and roughness (friction) of the underlying topography. For this reason, mechanisms to mitigate floods frequently employ structural measures that either modify topographic elevation, for example, through the placement of levees and sandbags, or increase roughness, for example, through revegetation projects. However, the configuration of these measures is typically decided in an ad hoc manner, limiting their overall effectiveness. The advent of high–performance surface–water modeling software and improvements in black–box optimization suggest that a more principled design methodology may be possible. This paper proposes a new computational approach to the problem of designing structural mitigation strategies under physical and budgetary constraints. It presents the development of a problem discretization amenable to simulation–based, derivative–free optimization. However, meta–heuristics alone are found to be insufficient for obtaining quality solutions in a reasonable amount of time. As a result, this paper proposes novel numerical and physics–based procedures to improve convergence to a high–quality mitigation. Furthermore the efficiency of the approach is demonstrated on hypothetical dam break scenarios of varying complexity under various mitigation budget constraints. In particular, experimental results show that, on average, the final proposed algorithm results in a 65% improvement in solution quality compared to a direct implementation.
54 ENVIRONMENTAL SCIENCES↗
Adaptive control of laser-wakefield accelerators driven by mid-IR laser pulses
There has been growing interest both in studying high intensity ultrafast laser plasma interactions with adaptive control systems as well as using long wavelength driver beams. We demonstrate the coherent control of the dynamics of laser-wakefield acceleration driven by ultrashort (~ 100 fs) mid-infrared (~ 3.9 μ m) laser pulses. The critical density at this wavelength is 7.3 × 10 19 cm -3 , which is achievable with an ordinary gas target system. Interactions between mid-infrared laser pulses and such near-critical-density plasma may be beneficial due to much higher absorption of laser energy. In addition, the normalized vector potential of the laser field a 0 increases with longer laser wavelength, lowering the required peak laser intensity to drive non-linear laser-wakefield acceleration. Here, MeV level, collimated electron beams with non-thermal, peaked energy spectra are generated. Optimization of electron beam qualities are realized through adaptive control of the laser wavefront. A genetic algorithm controlling a deformable mirror improves the electron total charge, energy spectra, beam pointing and stability at various plasma density profiles. Particle-in-cell simulations reveal that the optimal wavefront causes an earlier injection on the density up-ramp and thus higher energy gain as well as less filamentation during the interaction, which leads to the improvement in electron beam collimation and energy spectra.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Rapid torque-limited line-of-sight pointing of SCOLE (Spacecraft Control Laboratory Experiment) configuration
The design concept of a control for rapid torque-limited slewing of a rigid-mast version of the NASA SCOLE configuration is presented and demonstrated by means of numerical simulation. The time-optimal control problem for the system is decomposed into separate single-axis problems, expanding analytically the implicit nonlinear transcendental expression for the SCOLE line-of-sight error, and the final Euler attitude angles and slew angles are determined. The simulation results are presented in tables and graphs, and it is found that bang-bang or bang-pause-bang slew maneuvers with control moment applied to the Shuttle and control force applied to the reflector, and with a 5-deg/s slew-rate limit, produce the best pointing accuracy and the shortest slew times, although the specified line-of-sight error of 0.02 deg cannot be achieved using such open-loop single-axis maneuvers.
Lin, J. G.↗
Analytical modeling and sensor monitoring for optimal processing of polymeric composite material systems
Process simulation models and cure monitoring sensors are discussed for use in optimal processing of fiber-reinforced composites. Analytical models relate the specified temperature and pressure cure cycle to the thermal, chemical, and physical processes occurring in the composite during consolidation and cure. Frequency-dependent electromagnetic sensing (FDEMS) is described as an in situ sensor for monitoring the composite curing process and for verification of process simulation models. A model for resin transfer molding of textile composites is used to illustrate the predictive capabilities of a process simulation model. The model is used to calculate the resin infiltration time, fiber volume fraction, resin viscosity, and resin degree of cure. Results of the model are compared with in situ FDEMS measurements.
Loos, Alfred C.↗
Optimization of echo-enabled harmonic generation toward coherent EUV and soft X-ray free-electron laser at NSLS-II
Prebunching via echo-enabled harmonic generation (EEHG) is an efficient way to reduce the radiator length and improve the longitudinal coherence as well as output stability in storage-ring-based free-electron lasers. We propose a conceptual design, which uses two straight sections to seed coherent extreme-ultraviolet (EUV) and soft X-ray emission with nearly MHz repetition rate. To take the large energy spread (10 –3 ) of a storage ring into account and utilize the existing bending magnets between the two straight sections as the first chicane, we implement a special modeling tool, named EEHG optimizer. This tool has been successfully applied to maximize the prebunching with a reasonably low energy modulation, thereby generating intense coherent X-ray pulses within a short undulator length (a few meters) limited by the available space of a storage ring. Numerical simulations confirm that the optimized EEHG parameters can be directly applied to generate a 10 MW scale peak power with fully coherent ultrafast EUV to soft X-ray pulses based on the NSLS-II parameters. This method can be easily extended to other types of diffraction-limited storage rings.
36 MATERIALS SCIENCE↗
Lie-algebraic classical simulations for quantum computing
The classical simulation of quantum dynamics plays an important role in our understanding of quantum complexity and in the development of quantum technologies. Efficient techniques such as those based on the Gottesman-Knill theorem for Clifford circuits, tensor networks for low entanglement-generating circuits, or Wick's theorem for fermionic Gaussian states have become central tools in quantum computing. In this work, we contribute to this body of knowledge by presenting a framework for classical simulations, dubbed “𝔤-sim”, which is based on the underlying Lie algebraic structure of the dynamical process. When the dimension of the algebra grows at most polynomially in the system size, there exist observables for which the simulation is efficient. Indeed, we show that 𝔤-sim enables new regimes for classical simulations, is able to deal with certain forms of noise in the evolution, as well as can be used to tackle several paradigmatic variational and nonvariational quantum computing tasks. For the former, we perform Lie-algebraic simulations to train and optimize parametrized quantum circuits (thus effectively showing that some variational models can be dequantized), design enhanced parameter initialization strategies, solve tasks of quantum circuit synthesis, and train a quantum-phase classifier. For the latter, we report large-scale noiseless and noisy simulations on benchmark problems. By comparing the limitations of 𝔤-sim and certain Wick's theorem-based simulations, we find that the two methods become inefficient for different types of states or observables, hinting at the existence of distinct, nonequivalent resources for classical simulation.
97 MATHEMATICS AND COMPUTING↗