Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Nonlinear programming.”

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 415 records · Page 23

Fork Experiments in the Hot Cell Using Spent Fuel Rods for International Nuclear Safeguards

This work leveraged the rare availability of 25 full-length pressurized water reactor spent fuel rods and 1 irradiated mixed-oxide rod at an Oak Ridge National Laboratory hot cell. This was done to collect measurement data with two Fork detectors to assess the detectors’ capability of verifying operator declaration data and detecting partial defects in spent fuel, which are the two primary goals of international safeguards on spent nuclear fuel. The data can also be used to benchmark the ORIGEN module, which has been adopted in the International Atomic Energy Agency’s (IAEA’s)/European Atomic Energy Community’s (Euratom’s) Integrated Review and Analysis Program to predict the Fork detector count rates in real time. In this project, the authors first calibrated two Fork detectors—a standard one and a modified one—by using known strong neutron and gamma sources. Then, the authors measured all 26 fuel rods at multiple locations along the length. The fuel rods were then assembled into three arrays—2 × 2, 3 × 3, and 5 × 5—by using specially designed support grids to mimic fuel assemblies and measure the arrays with both detectors. For the 5 × 5 array, 4 and 8 fuel rods of the array were replaced in two separate cases with short stainless-steel rods to mimic two partial defect scenarios, and the arrays were measured before and after the replacements. Polyethylene blocks were used in this experiment to mimic water. The results show that the Fork detectors were able to verify operator declarations and detect partial defects in spent fuel, and the authors were the first to demonstrate this through experiments. A discovery was also made that determined the root cause of the nonlinear response to gamma dose in the ion chambers used in IAEA and Euratom’s Fork detectors. After the experiments, both detectors were retrieved from the hot cell for future use. The data collected in this project will be used in a parallel International Nuclear Safeguards Engagement Program (INSEP) project to enhance the safeguards in the Finnish spent fuel encapsulation plant, and the data will be useful to other projects in the future given the increased safeguards needs due to spent fuel transfer and disposal activities worldwide.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Advanced multiphysics coupling for LWR fuel performance analysis

Even the most basic nuclear fuel analysis is a multiphysics undertaking, as a credible simulation must consider at a minimum coupled heat conduction and mechanical deformation. The need for more realistic fuel modeling under a variety of conditions invariably leads to a desire to include coupling between a more complete set of the physical phenomena influencing fuel behavior, including neutronics, thermal hydraulics, and mechanisms occurring at lower length scales. This paper covers current efforts toward coupled multiphysics LWR fuel modeling in three main areas. The first area covered in this paper concerns thermomechanical coupling. The interaction of these two physics, particularly related to the feedback effect associated with heat transfer and mechanical contact at the fuel/clad gap, provides numerous computational challenges. An outline is provided of an effective approach used to manage the nonlinearities associated with an evolving gap in BISON, a nuclear fuel performance application. A second type of multiphysics coupling described here is that of coupling neutronics with thermomechanical LWR fuel performance. DeCART, a high-fidelity core analysis program based on the method of characteristics, has been coupled to BISON. DeCART provides sub-pin level resolution of the multigroup neutron flux, with resonance treatment, during a depletion or a fast transient simulation. Two-way coupling between these codes was achieved by mapping fission rate density and fast neutron flux fields from DeCART to BISON and the temperature field from BISON to DeCART while employing a Picard iterative algorithm. Finally, the need for multiscale coupling is considered. Fission gas production and evolution significantly impact fuel performance by causing swelling, a reduction in the thermal conductivity, and fission gas release. The mechanisms involved occur at the atomistic and grain scale and are therefore not the domain of a fuel performance code. However, it is possible to use lower length scale models such as those used in the mesoscale MARMOT code to compute average properties, e.g. swelling or thermal conductivity. These may then be used by an engineering-scale model. Examples of this type of multiscale, multiphysics modeling are shown.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Coiling free electron matter waves

Here we demonstrate particle beams that spiral in free space devoid of external fields. The beams consist of electrons in two lobes that twist around each other along the optical axis, such that each electron can be described by a two-lobed probability distribution that rotates as it propagates. Furthermore, we demonstrate that this twisting distribution can undergo programmed periods of angular acceleration. These unusual states are produced by preparing each free electron wavefunction in a superposition of non-diffracting Bessel modes that carry orbital angular momentum using nanofabricated diffraction holograms. The holograms can encode nonlinear azimuthal phase so that the resulting electron probability distribution twists in space in a controllable manner, accelerating and decelerating during propagation, without any radiation. This work provides a new platform to explore the dynamics of electrons in magnetic fields, and opens the possibility of producing other types of particle beams with coiling geometries.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

FAST INJECTION SYSTEM R&D FOR THE APS UPGRADE

The MBA upgrade for the APS will operate with bunch swap out and on axis injection. The planned 324 bunch fill pattern places difficult demands on the injection and extraction kickers. The present concept uses dual stripline kickers driven by high Voltage pulsers. Minimizing perturbation on adjacent bunches requires very fast rise and fall times with relatively narrow ~20 nsec, 15 kV pulses. To achieve these requirements we have initiated a multifaceted R&D program. The R&D includes the HV pulser, stripline kicker and HV feedthrough. We have purchased a commercial dual channel HV pulser and are evaluating its performance and reliability. In addition, we are investigating the feasibility of using nonlinear ferrite loaded coaxial cables (shockwave transmission line) to sharpen the leading and trailing edges of high voltage pulses. We are also developing a prototype kicker and high voltage feedthrough. The requirements for injection and extraction, progress on prototype development and results of our HV pulser investigations will be reported.

Lenkszus, F.↗

Final Report - Interaction of radiation and charged particles in miniature plasma structures

The extension of our program to the development of theories and models capable of describing the interactions of intense laser pulses and charged particles in miniature plasma channels is reported. These channels, which have recently been created in the laboratory, have unique dispersion properties that make them interesting for a variety of applications including particle acceleration, high harmonic generation, and THz generation. Our program systematically explored the properties of these channels, including dispersion, losses, and coupling. A particular application that was pursued is the generation of intense pulses of THz radiation by short laser pulses propagating these channels. We also explored the nonlinear dynamics of laser pulses propagating in these channels.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Data-driven modeling and control of dynamical systems using Koopman and Perron-Frobenius operators

This dissertation studies the data-driven modeling and control problem of nonlinear systems by exploiting the linear operator theoretic framework involving Koopman and Perro-Frobenius operator. A systematic linear-operator based controller design procedure has been established, which can be used to solve a variety of nonlinear control problems, including feedback stabilization using control Lyapunov functions, optimal quadratic regulation using Koopman eigenfunctions and convex optimization formulation of optimal control problem using P-F and Koopman operator approximation. As the core of data-driven modeling, we first propose a new algorithm for the finite-dimensional approximation of the linear transfer Koopman and Perron-Frobenius operator from time-series data. We argue that the existing approach for the finite-dimensional approximation of these transfer operators such as Dynamic Mode Decomposition (DMD) and Extended Dynamic Mode Decomposition (EDMD) do not capture two important properties of these operators, namely positivity and Markov property. The algorithm we propose preserves these two properties. We call the proposed algorithm as naturally structured DMD (NSDMD) since it retains the inherent properties of these operators. Naturally structured DMD algorithm leads to a better approximation of the steady-state dynamics of the system regarding computing Koopman and Perron- Frobenius operator eigenfunctions and eigenvalues. However, preserving positivity property is critical for capturing the real transient dynamics of the system. This positivity property of the transfer operators and it's finite-dimensional approximation play an important role for controller and estimator design of nonlinear systems. To solve the feedback stabilization problem for nonlinear control systems, we tried to take advantage of the Koopman operator framework. The Koopman operator approach provides a linear representation for a nonlinear dynamical system and a bilinear representation for a nonlinear control system. The problem of feedback stabilization of a nonlinear control system is then transformed to the stabilization of a bilinear control system. We propose a control Lyapunov function (CLF)-based approach for the design of stabilizing feedback controllers for the bilinear system. The search for finding a CLF for the bilinear control system is formulated as a convex optimization problem. This leads to a schematic procedure for designing CLF-based stabilizing feedback controllers for the bilinear system and hence the original nonlinear system. Another advantage of the proposed controller design approach outlined in this dissertation is that it does not require explicit knowledge of system dynamics. In particular, the bilinear representation of a nonlinear control system in the Koopman eigenfunction space can be obtained from time-series data. Next, we study the optimal quadratic regulation problem for nonlinear systems. The linear operator theoretic framework involving the Koopman operator is used to lift the dynamics of nonlinear control system to an infinite-dimensional bilinear system. The optimal quadratic regulation problem for nonlinear system is formulated in terms of the finite-dimensional approximation of the bilinear system. A convex optimization-based approach is proposed for solving the quadratic regulator problem for bilinear system. We applied a variety of examples and compared the simulation results between our framework and conventional LQR control using linearized model. For more general optimal control problems, we provide a density-function based convex formulation for the optimal control problem of the nonlinear system. The convex formulation relies on the duality result in the stability theory of a dynamical system involving density function and Perron-Frobenius operator. The optimal control problem is formulated as an infinite-dimensional convex optimization program. The finite-dimensional approximation of the optimization problem relies on the recent advances made in the data-driven computation of the Koopman operator, which is dual to the Perron-Frobenius operator. Simulation results are presented to demonstrate the application of the developed framework.

Huang, Bowen↗

DFT calculations on spectroscopic and structural properties of a NLO chromophore

The molecular geometry optimization, vibrational frequencies and gauge including atomic orbital (GIAO) {sup 1}H and {sup 13}C NMR chemical shift values of 2-(1′-(4’’’-Methoxyphenyl)-5′-(thien-2″-yl)pyrrol-2′-yl)-1,3-benzothiazole as potential nonlinear optical (NLO) material were calculated using density functional theory (DFT) HSEh1PBE method with 6–311G(d,p) basis set. The best of our knowledge, this study have not been reported to date. Additionally, a detailed vibrational study was performed on the basis of potential energy distribution (PED) using VEDA program. It is noteworthy that NMR chemical shifts are quite useful for understanding the relationship between the molecular structure and electronic properties of molecules. The computed IR and NMR spectra were used to determine the types of the experimental bands observed. Predicted values of structural and spectroscopic parameters of the chromophore were compared with each other so as to display the effects of the different substituents on the spectroscopic and structural properties. Obtained data showed that there is an agreement between the predicted and experimental data.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A fast and accurate domain decomposition nonlinear manifold reduced order model

Here, this paper integrates nonlinear-manifold reduced order models (NM-ROMs) with domain decomposition (DD). NM ROMs approximate the full order model (FOM) state in a nonlinear-manifold by training a shallow, sparse autoencoder using FOM snapshot data. These NM-ROMs can be advantageous over linear-subspace ROMs (LS-ROMs) for problems with slowly decaying Kolmogorov n-width. However, the number of NM-ROM parameters that need to be trained scales with the size of the FOM. Moreover, for “extreme-scale” problems, the storage of high-dimensional FOM snapshots alone can make ROM training expensive. To alleviate the training cost, this paper applies DD to the FOM, computes NM-ROMs on each subdomain, and couples them to obtain a global NM-ROM. This approach has several advantages: Subdomain NM-ROMs can be trained in parallel, involve fewer parameters to be trained than global NM-ROMs, require smaller subdomain FOM dimensional training data, and can be tailored to subdomain specific features of the FOM. The shallow, sparse architecture of the autoencoder used in each subdomain NM-ROM allows application of hyper-reduction (HR), reducing the complexity caused by nonlinearity and yielding computational speedup of the NM-ROM. This paper provides the first application of NM-ROM (with HR) to a DD problem. In particular, this paper details an algebraic DD reformulation of the FOM, training a NM-ROM with HR for each sub domain, and a sequential quadratic programming (SQP) solver to evaluate the coupled global NM-ROM. Theoretical convergence results for the SQP method and a priori and a posteriori error estimates for the DD NM-ROM with HR are provided. The proposed DD NM-ROM with HR approach is numerically compared to a DD LS-ROM with HR on the 2D steady-state Burgers’ equation, showing an order of magnitude improvement in accuracy of the proposed DD NM-ROM over the DD LS-ROM.

97 MATHEMATICS AND COMPUTING↗

THERMAL DIAGNOSTICS WITH THE ATMOSPHERIC IMAGING ASSEMBLY ON BOARD THE SOLAR DYNAMICS OBSERVATORY: A VALIDATED METHOD FOR DIFFERENTIAL EMISSION MEASURE INVERSIONS

We present a new method for performing differential emission measure (DEM) inversions on narrow-band EUV images from the Atmospheric Imaging Assembly (AIA) on board the Solar Dynamics Observatory. The method yields positive definite DEM solutions by solving a linear program. This method has been validated against a diverse set of thermal models of varying complexity and realism. These include (1) idealized Gaussian DEM distributions, (2) 3D models of NOAA Active Region 11158 comprising quasi-steady loop atmospheres in a nonlinear force-free field, and (3) thermodynamic models from a fully compressible, 3D MHD simulation of active region (AR) corona formation following magnetic flux emergence. We then present results from the application of the method to AIA observations of Active Region 11158, comparing the region's thermal structure on two successive solar rotations. Additionally, we show how the DEM inversion method can be adapted to simultaneously invert AIA and Hinode X-ray Telescope data, and how supplementing AIA data with the latter improves the inversion result. The speed of the method allows for routine production of DEM maps, thus facilitating science studies that require tracking of the thermal structure of the solar corona in time and space.

79 ASTRONOMY AND ASTROPHYSICS↗

FY2014 FES (Fusion Energy Sciences) Theory & Simulation Performance Target, Final Report

We report here the work done for the FY14 OFES Theory Performance Target as given below: "Understanding alpha particle confinement in ITER, the world's first burning plasma experiment, is a key priority for the fusion program. In FY 2014, determine linear instability trends and thresholds of energetic particle-driven shear Alfven eigenmodes in ITER for a range of parameters and profiles using a set of complementary simulation models (gyrokinetic, hybrid, and gyrofluid). Carry out initial nonlinear simulations to assess the effects of the unstable modes on energetic particle transport". In the past year (FY14), a systematic study of the alpha-driven Alfven modes in ITER has been carried out jointly by researchers from six institutions involving seven codes including the transport simulation code TRANSP (R. Budny and F. Poli, PPPL), three gyrokinetic codes: GEM (Y. Chen, Univ. of Colorado), GTC (J. McClenaghan, Z. Lin, UCI), and GYRO (E. Bass, R. Waltz, UCSD/GA), the hybrid code M3D-K (G.Y. Fu, PPPL), the gyro-fluid code TAEFL (D. Spong, ORNL), and the linear kinetic stability code NOVA-K (N. Gorelenkov, PPPL). A range of ITER parameters and profiles are specified by TRANSP simulation of a hybrid scenario case and a steady-state scenario case. Based on the specified ITER equilibria linear stability calculations are done to determine the stability boundary of alpha-driven high-n TAEs using the five initial value codes (GEM, GTC, GYRO, M3D-K, and TAEFL) and the kinetic stability code (NOVA-K). Both the effects of alpha particles and beam ions have been considered. Finally, the effects of the unstable modes on energetic particle transport have been explored using GEM and M3D-K.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Highlights of the 2016 Summer Laboratory Student Intern Program at Sandia National Labs [Slides]

Combined picosecond (10 -12 s) and femtosecond (10 -15 s) laser pulses can give sensitive, low-noise measurements of important quantities in reacting flows, such as species concentrations and temperature. Emmanuel’s work focused on the development of an instrument for tailoring the time profile of picosecond laser pulses for use in nonlinear optical spectroscopic methods created from broad bandwidth femtosecond pulses. In addition to constructing the device, Emmanuel produced a LabView-based automation code, building off skills he developed in a previous CCI internship at Sandia.

42 ENGINEERING↗

Static actuator-sharing algorithm for concurrent control of multiple plasma properties

Simultaneous regulation of multiple properties in next-generation tokamaks like ITER and fusion pilot plant may require the integration of different plasma control algorithms. Such integration requires the conversion of individual controller commands into physical actuator requests while accounting for the coupling between different plasma properties. This work proposes a tokamak and scenario-agnostic actuator-sharing algorithm (ASA) to perform the above-mentioned command-request conversion and, hence, integrate multiple plasma controllers. The proposed algorithm implicitly solves a quadratic programming (QP) problem formulated to account for the saturation limits and the relation between the controller commands and physical actuator requests. Since the constraints arising in the QP program are linear, the proposed ASA is highly computationally efficient and can be implemented in the tokamak plasma control system in real time. Furthermore, the proposed algorithm is designed to handle real-time changes in the control objectives and actuators’ availability. Nonlinear simulations carried out using the Control Oriented Transport SIMulator illustrate the effectiveness of the proposed algorithm in achieving multiple control objectives simultaneously.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Theoretical and Experimental Studies in Accelerator Physics

This report describes research supported by the US Dept. of Energy Office of High Energy Physics (OHEP), performed by the UCLA Particle Beam Physics Laboratory (PBPL). The UCLA PBPL has, over the last two decades-plus, played a critical role in the development of advanced accelerators, fundamental beam physics, and new applications enabled by these thrusts, such as new types of accelerator-based light sources. As the PBPL mission is broad it is natural that it has been grown within the context of the accelerator science and technology stewardship of the OHEP. Indeed, steady OHEP support for the program has always been central to the success of the PBPL; it has provided stability, and above all has set the over-arching themes for our research directions, which have producing over 500 publications (>120 in high level journals). While other agency support has grown notably in recent years, permitting more vigorous pursuit of the program, it is transient by comparison. Beyond permitting program growth in a time of flat OHEP budgets, the influence of other agency missions is found in push to adapt advanced accelerator methods to applications, in light of the success the field has had in proof-of-principle experiments supported first by the DoE OHEP. This three-pronged PBPL program — advanced accelerators, fundamental beam physics and technology, and revolutionary applications — has produced a generation of students that have had a profound affect on the US accelerator physics community. PBPL graduates, numbering 28 in total, form a significant population group in the accelerator community, playing key roles as university faculty, scientific leaders in national labs (two have been named Panofsky Fellows at SLAC), and vigorous proponents of industrial application of accelerators. Indeed, the development of advanced RF, optical and magnet technology at the PBPL has led directly to the spin-off company, RadiaBeam Technologies, now a leading industrial accelerator firm. We note also that PBPL graduates remain as close elaborators for the program after leaving UCLA. The UCLA PBPL program is a foremost developer of on-campus facilities, such as the Neptune and Pegasus Laboratories, providing a uniquely strong environment for student-based research. In addition, the PBPL is a strong user of off-campus national lab facilities, such as SLAC FACET and NLCTA, and the BNL ATF. UCLA has also vigorously participated in the development of these facilities. The dual emphases on off- and on-campus opportunities permit the PBPL to address in an agile way a wide selection of cutting-edge research topics. The topics embraced by this proposal illustrate this program aspect well. These include: GV/m dielectric wakefield acceleration/coherent Cerenkov radiation experiments at FACET (E-201) and the ATF; synergistic laser-excited dielectric accelerator and light source development; plasma wakefield (PWFA) experiments on “Trojan horse” ionization injection (FACET E-210), quasi-nonlinear PWFA at BNL and the production at Neptune high transformer ratio plasma wakes; the inauguration of a new type of RF photoinjector termed “hybrid” at UCLA, and application to PWFA; space-charge dominated beam and cathode/near cathode physics; the study of advanced IFEL systems, for very high energy gain and utilization of novel OAM modes; the physcis of inverse Compton scattering (ICS), with applications to e+ production and γγ colliders; electron diffraction; and advanced beam diagnostics using coherent imaging techniques. These subjects are addressed under the leadership of PBPL director Prof. James Rosenzweig in Task A, and Prof. Pietro Musumeci in Task J, which was initiated following his OHEP Outstanding Junior Investigator award.

43 PARTICLE ACCELERATORS↗

Usability and Optimization Improvements in MOOSE

The Multiphysics Object-Oriented Simulation Environment (MOOSE) framework is a foundational capability used by the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program to create over 15 different simulation tools for advanced nuclear reactors. Due to MOOSE’s broad use, improvements to the framework in support of modeling and simulation goals are critical to the program. Such improvements can take many forms, including optimization, improved user experience, streamlined application programming interfaces (APIs), parallelism, and new capabilities. The work described in this report was conducted in direct support of the simulation tools and has already been deployed. The capabilities were implemented in the same order as they are covered in this report: multiple nonlinear systems in the same input file, implementation of generic field transfers and other transfer system enhancements, support for stateful material property redistribution for adaptivity and distributed meshes, and dynamic linking and loading of individually compiled applications. These four additions are fundamental capabilities that will be leveraged by many NEAMS applications.

97 MATHEMATICS AND COMPUTING↗

Feedback Control Approaches for Restoration of Power Grids from Blackouts

The automated restoration of power systems with variable energy resources is a timely problem to tackle. Automated restoration advice can support operators in deciding on strategic actions to restore power grids from a blackout with a mix of conventional and renewable generation resources. To this end, this paper frames the restoration process of power grids with solar resources as a nonlinear dynamic model with algebraic constraints in discrete time which is steered by feedback control loops. We discuss two feedback-control strategies based on greedy and reinforcement learning algorithms, and contrast their performance with restoration plans generated by a mixed-integer linear program. We found that the reinforcement learning algorithm infers restoration actions faster than the greedy one. However, the tuning process of the reinforcement learning parameters is slower than for the greedy one.

machine learning↗

Pyomo.GDP: an ecosystem for logic based modeling and optimization development

We present three core principles for engineering-oriented integrated modeling and optimization tool sets—intuitive modeling contexts, systematic computer-aided reformulations, and flexible solution strategies—and describe how new developments in Pyomo.GDP for Generalized Disjunctive Programming (GDP) advance this vision. We describe a new logical expression system implementation for Pyomo.GDP allowing for a more intuitive description of logical propositions. The logical expression system supports automated reformulation of these logical constraints to linear constraints. We also describe two new logic-based global optimization solver implementations built on Pyomo.GDP that exploit logical structure to avoid “zero-flow” numerical difficulties that arise in nonlinear network design problems when nodes or streams disappear. These new solvers also demonstrate the capability to link to external libraries for expanded functionality within an integrated implementation. We present these new solvers in the context of a flexible array of solution paths available to GDP models. Finally, we present results on a new library of GDP models demonstrating the value of multiple solution approaches.

42 ENGINEERING↗

Energy Impacts of Nonlinear Behavior of PCM When Applied into Building Envelope

Presented at the ASME 2012 6th International Conference on Energy Sustainability & 10th Fuel Cell Science, Engineering and Technology Conference on July 23-26, 2012, this study analyzes the effects a nonlinear enthalpy profile has on thermal performance and expected energy benefits for PCM-enhanced insulation.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

The TRANFT User Manual version 2.0

The Fortran program TRANFT has undergone significant revision. While many of the original features remain others have been modified and new capabilities have been added. A second trans verse direction is now tracked so that the tune and force distributions associated with the beam beam force and nonlinear space charge are correct. Also, general mutibunch instabilites are now treated. For symmetric coupled bunch modes one can track one bunch for purely transverse modes or a few contiguous bunches to include longitudinal modes as well. Another option which tracks all the bunchs in the ring with provision for a clearing or abort gap is also available. In addition to this a new algorithm to model the transverse voltage on crab cavities is introduced which includes the effect of various feedbacks.

43 PARTICLE ACCELERATORS↗