Search NASASearch

SEARCH · Search NASA

Results for “simulation optimization”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 307 records · Page 17

Optimizing the Resolution of Hydrodynamic Simulations for MCRaT Radiative Transfer Calculations

Despite their discovery about half a century ago, the Gamma-ray burst (GRB) prompt emission mechanism is still not well understood. Theoretical modeling of the prompt emission has advanced considerably due to new computational tools and techniques. One such tool is the PLUTO hydro dynamics code, which is used to numerically simulate GRB outflows. PLUTO uses Adaptive Mesh Refinement to focus computational efforts on the portion of the grid that contains the simulated jet. Another tool is the Monte Carlo Radiation Transfer (MCRaT) code, which predicts electromagnetic signatures of GRBs by conducting photon scatterings within a jet using PLUTO. The effects of the underlying resolution of a PLUTO simulation with respect to MCRaT post-processing radiative transfer results have not yet been quantified. We analyze an analytic spherical outflow and a hydrodynamically simulated GRB jet with MCRaT at varying spatial and temporal resolutions and quantify how decreasing both resolutions affect the resulting mock observations. We find that changing the spatial resolution changes the hydrodynamic properties of the jet, which directly affect the MCRaT mock observable peak energies. We also find that decreasing the temporal resolution artificially decreases the high energy slope of the mock observed spectrum, which increases both the spectral peak energy and the luminosity. We show that the effects are additive when both spatial and temporal resolutions are modified. Our results allow us to understand how decreased hydrodynamic temporal and spatial resolutions affect the results of post-processing radiative transfer calculations, allowing for the optimization of hydrodynamic simulations for radiative transfer codes.

Jose Arita-Escalante

Recursive Branching Simulated Annealing Algorithm

This innovation is a variation of a simulated-annealing optimization algorithm that uses a recursive-branching structure to parallelize the search of a parameter space for the globally optimal solution to an objective. The algorithm has been demonstrated to be more effective at searching a parameter space than traditional simulated-annealing methods for a particular problem of interest, and it can readily be applied to a wide variety of optimization problems, including those with a parameter space having both discrete-value parameters (combinatorial) and continuous-variable parameters. It can take the place of a conventional simulated- annealing, Monte-Carlo, or random- walk algorithm. In a conventional simulated-annealing (SA) algorithm, a starting configuration is randomly selected within the parameter space. The algorithm randomly selects another configuration from the parameter space and evaluates the objective function for that configuration. If the objective function value is better than the previous value, the new configuration is adopted as the new point of interest in the parameter space. If the objective function value is worse than the previous value, the new configuration may be adopted, with a probability determined by a temperature parameter, used in analogy to annealing in metals. As the optimization continues, the region of the parameter space from which new configurations can be selected shrinks, and in conjunction with lowering the annealing temperature (and thus lowering the probability for adopting configurations in parameter space with worse objective functions), the algorithm can converge on the globally optimal configuration. The Recursive Branching Simulated Annealing (RBSA) algorithm shares some features with the SA algorithm, notably including the basic principles that a starting configuration is randomly selected from within the parameter space, the algorithm tests other configurations with the goal of finding the globally optimal solution, and the region from which new configurations can be selected shrinks as the search continues. The key difference between these algorithms is that in the SA algorithm, a single path, or trajectory, is taken in parameter space, from the starting point to the globally optimal solution, while in the RBSA algorithm, many trajectories are taken; by exploring multiple regions of the parameter space simultaneously, the algorithm has been shown to converge on the globally optimal solution about an order of magnitude faster than when using conventional algorithms. Novel features of the RBSA algorithm include: 1. More efficient searching of the parameter space due to the branching structure, in which multiple random configurations are generated and multiple promising regions of the parameter space are explored; 2. The implementation of a trust region for each parameter in the parameter space, which provides a natural way of enforcing upper- and lower-bound constraints on the parameters; and 3. The optional use of a constrained gradient- search optimization, performed on the continuous variables around each branch s configuration in parameter space to improve search efficiency by allowing for fast fine-tuning of the continuous variables within the trust region at that configuration point.

Bolcar, Matthew

From Vehicles to Systems: Understanding Freight Transportation as a Connected Energy, Infrastructure, and Operations System

The U.S. freight system may need to handle 50% more cargo by 2050. Upgrading our freight system requires modernizing capital-intensive, long-lived assets including freight trains, ports, and terminal infrastructure. NLR is advancing freight system solutions spanning ALTRIOS, the first digital twin for the full freight rail system; ALTRIOS-LIFTS, which can create digital twins of freight terminals; INFORMES, the first national model of the intermodal freight system; MARINESim, used to simulate and optimize ocean-going vessel operations; and more. These modeling and simulation tools enable data-driven decision-making across freight modes and systems.

24 POWER TRANSMISSION AND DISTRIBUTION

Simulation of Airplane and Rocket Trajectories

Simulation and Optimization of Rocket Trajectories program (SORT) contains comprehensive mathematical models for simulating aircraft dynamics, freely falling objects, and many types of ballistic trajectories. Provides high-fidelity, three-degrees-of-freedom simulation for atmospheric and exoatmospheric flight. It numerically models vehicle subsystems and vehicle environment. Used for wide range of simulations. Written in machine-independent FORTRAN 77.

Wahbah, Magdy M.

Optimization of District Heating Network Parameters in Steady-State Operation

Here we examine the modeling, simulation, and optimization of district heating systems, which are widely used for thermal transport using steam or hot water as a carrier. We propose a generalizable framework to specify network models and scenario parameters, and develop an optimization method for evaluating system states including pressures, fluid flowrates, and temperatures throughout the network. The network modeling includes pipes, thermal plants, pumps, and passive or controllable loads as system components. We propose basic models for thermodynamic fluid transport and enforce the balance of physical quantities in steady-state flow over co-located outgoing and return networks. We formulate an optimization problem with steam and hot water as the outgoing and return carriers, as in legacy twentieth century systems. The physical laws and engineering limitations are specified for each component type, and the thermal network flow optimization problem is formulated and solved for a realistic test network under several scenarios.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI

Simulation and Analysis of Technology and Operational Procedures to Reduce the Combined Effects of Emissions and Contrails

The development and evaluation of concepts and technology to support future air traffic management systems require a hierarchy of models ranging from real-time simulations to extensive field evaluations. Air traffic simulation models such as Airspace Concept Evaluation System, Center Tracon Automation System, Future Air traffic management Concept Evaluation Tool and others are used to design air traffic systems balancing the conflicting objectives of maximizing safety, meeting future demands for airports and airspace and increase efficiency of traffic flows in the presence of uncertain weather. The impact of aviation emissions and contrails on climate imposes another constraint on the design of aircraft and aviation operations. The understanding of the complex interaction between physical climate system, the carbon and other greenhouse gas emissions and aviation activity can be improved by the development of integrated assessment models that include emission and climate models together with air traffic simulations. The modeling of aircraft emissions and their interaction with each other to change the concentration levels of different gasses in the atmosphere and the resulting impact of the radiative forcing on the equilibrium of the Earth's atmosphere is complex and requires the use of coupled atmosphere-ocean general circulation models together with three-dimensional models of carbon cycle and chemistry of other non-CO2 greenhouse gases. These models are computationally intensive and unsuitable for studies involving the generation of multiple scenarios. Simple emission and climate models, based on the input-output relations of linear systems, capture the fundamental emission to climate impact behavior by careful selection of key variables and their dynamics. The impact of various greenhouse gases depends on the total concentration, effect per unit change in atmospheric concentration and the spatial distribution of the gas. All these quantities are influenced by the lifetime of the gas. The impact of a greenhouse gas depends on the interval of assessment, which may vary from a few decades to a few centuries. Climate metrics are aimed at providing a common scale to compare different greenhouse gases. If the metrics are to be used as a tool in developing and evaluating aviation operations, they should be transparent and easy to apply. Global Warming Potential and Aggregate Global Temperature Potential are some of the commonly used metrics. This paper integrates a national-level air traffic simulation and optimization capability with simple climate models and carbon cycle models, and climate metrics to assess the impact of aviation on climate. The capability brings together metrics, which are useful in aviation operations together with metrics used in climate studies. The capability can be used to make trade-offs between extra fuel cost and reduction in climate impact. There is considerable uncertainty in our understanding of the radiative forcing associated with emissions and contrails. The parameters in the simulation can be used to evaluate the effect of various uncertainties in emission models and contrails. It can also be used to evaluate the impact of different decision horizons. Alternatively, the optimization results from the simulation can be used as inputs to other tools that monetize global climate impacts like the FAA's Aviation Environmental Portfolio Management Tool for Impacts.

impact of aviation on the environment

Optimal washout for control of a moving base simulator

A general form of an optimal washout filter is derived using state-space linear optimal control theory, and this is applied to the design of washout filters of various types of moving base motion simulators, including the NASA's vertical motion simulator. Attention is given to the linear elements of a washout filter. One of the nonlinearities considered is braking which may be required near the end of the simulator excursion to prevent a crash. Although the general form of the optimal washout filter is applicable to time-variant system, the applications analyzed in the study are restricted to time-invariant cases.

Kurosaki, M.

Design and Analysis of Short Period 2G-HTS Undulators

Recent advancements in second-generation high-temperature superconducting (2G-HTS) materials have opened new possibilities for superconducting undulator (SCU) technology, particularly in achieving shorter periods and higher magnetic fields. Here, this paper presents the design and analysis of a short-period, as small as 10 mm, 2G-HTS undulator, focusing on magnetic performance, mechanical force minimization, and quench protection. A force-balanced coil pack design is introduced to significantly reduce mechanical stresses on the conductors, eliminating the need for heavy support structures. Rapid prototyping techniques such as3D printing is employed to refine the design for scalability and manufacturability. Correction schemes for end-field effects are developed using a combination of 2G-HTS tapes and NbTi conductors to optimize field integrals. Simulation results demonstrated that optimized magnetic performance and enhanced mechanical stability, positioning 2G-HTS undulators as promising candidates for storage ring and free-electron based light sources.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND

OTIS Modeling of IBEX Flight

When a launch vehicle carries a payload into space, it must follow a very particular trajectory to achieve its desired orbit. Modern rockets have very tight delta-v budgets and many forces work to rob the rocket of energy, placing limits on the performance of the launch vehicle. Trajectory optimization evolved out of this need to economize on every bit of energy of a launch vehicle. Today NASA uses Optimal Trajectory by Implicit Simulation (OTIS) to optimize trajectories based on given parameters. The Interstellar Boundary EXplorer (IBEX) mission profile is unique and required special alterations to the standard OTIS data output.

deZarraga, Gaston

Fuzzy logic mode switching in helicopters

The application of fuzzy logic to a wide range of control problems has been gaining momentum internationally, fueled by a concentrated Japanese effort. Advanced Research & Development within the Engineering Department at Sikorsky Aircraft undertook a fuzzy logic research effort designed to evaluate how effective fuzzy logic control might be in relation to helicopter operations. The mode switching module in the advanced flight control portion of Sikorsky's motion based simulator was identified as a good candidate problem because it was simple to understand and contained imprecise (fuzzy) decision criteria. The purpose of the switching module is to aid a helicopter pilot in entering and leaving coordinated turns while in flight. The criteria that determine the transitions between modes are imprecise and depend on the varied ranges of three flight conditions (i.e., simulated parameters): Commanded Rate, Duration, and Roll Attitude. The parameters were given fuzzy ranges and used as input variables to a fuzzy rulebase containing the knowledge of mode switching. The fuzzy control program was integrated into a real time interactive helicopter simulation tool. Optimization of the heading hold and turn coordination was accomplished by interactive pilot simulation testing of the handling quality performance of the helicopter dynamic model. The fuzzy logic code satisfied all the requirements of this candidate control problem.

Sherman, Porter D.

A sensitivity equation approach to shape optimization in fluid flows

A sensitivity equation method to shape optimization problems is applied. An algorithm is developed and tested on a problem of designing optimal forebody simulators for a 2D, inviscid supersonic flow. The algorithm uses a BFGS/Trust Region optimization scheme with sensitivities computed by numerically approximating the linear partial differential equations that determine the flow sensitivities. Numerical examples are presented to illustrate the method.

Borggaard, Jeff

Discrete-State Simulated Annealing For Traveling-Wave Tube Slow-Wave Circuit Optimization

Algorithms based on the global optimization technique of simulated annealing (SA) have proven useful in designing traveling-wave tube (TWT) slow-wave circuits for high RF power efficiency. The characteristic of SA that enables it to determine a globally optimized solution is its ability to accept non-improving moves in a controlled manner. In the initial stages of the optimization, the algorithm moves freely through configuration space, accepting most of the proposed designs. This freedom of movement allows non-intuitive designs to be explored rather than restricting the optimization to local improvement upon the initial configuration. As the optimization proceeds, the rate of acceptance of non-improving moves is gradually reduced until the algorithm converges to the optimized solution. The rate at which the freedom of movement is decreased is known as the annealing or cooling schedule of the SA algorithm. The main disadvantage of SA is that there is not a rigorous theoretical foundation for determining the parameters of the cooling schedule. The choice of these parameters is highly problem dependent and the designer needs to experiment in order to determine values that will provide a good optimization in a reasonable amount of computational time. This experimentation can absorb a large amount of time especially when the algorithm is being applied to a new type of design. In order to eliminate this disadvantage, a variation of SA known as discrete-state simulated annealing (DSSA), was recently developed. DSSA provides the theoretical foundation for a generic cooling schedule which is problem independent, Results of similar quality to SA can be obtained, but without the extra computational time required to tune the cooling parameters. Two algorithm variations based on DSSA were developed and programmed into a Microsoft Excel spreadsheet graphical user interface (GUI) to the two-dimensional nonlinear multisignal helix traveling-wave amplifier analysis program TWA3. The algorithms were used to optimize the computed RF efficiency of a TWT by determining the phase velocity profile of the slow-wave circuit. The mathematical theory and computational details of the DSSA algorithms will be presented and results will be compared to those obtained with a SA algorithm.

Wilson, Jeffrey D.

A practical six-degree of freedom solar sail dynamics model for optimizing solar sail trajectories with torque constraints

Controlled flight of a solar sail-propelled spacecraft ('sailcraft') is a six-degree-of-freedom dynamics problem. Current state-of-the-art tools that simulate and optimize the trajectories flown by sailcraft do not treat the full kinetic (i.e. force and torque-constrained) motion, instead treating a discrete history of commanded sail attitudes, and either neglecting the sail attitude motion over an integration timestep, or treating the attitude evolution kinematically with a spline or similar treatment. The present paper discusses an aspect of developing a next generation sailcraf trajectory designing optimization tool JPL, for NASA's Solar Sail Spaceflight Simulation Software (SS). The aspect discussed in an experimental approach to modeling full six-degree-of-freedom kinetic motion of a solar sail in a trajectory propagator. Early results from implementing this approach in a new trajectory propagation tool are given.

six-degree-of-freedom

L-Band Microwave Satellite Data and Model Simulations Over the Dry Chaco to Estimate Soil Moisture, Soil Temperature, Vegetation and Soil Salinity

The Dry Chaco in South America is a semi-arid ecoregion prone to dryland salinization. In this region, we investigated coarse-scale surface soil moisture (SM), soil temperature, soil salinity and vegetation, using L-band microwave brightness temperature (TB) observations and retrievals from the Soil Moisture Ocean Salinity (SMOS) and Soil Moisture Active Passive (SMAP) satellite missions, Catchment Land Surface Model (CLSM) simulations, and in situ measurements within 26 sampled satellite pixels. Across these 26 sampled pixels, the satellite-based SM outperformed CLSM SM compared to field data, and forward L-band TB simulations derived from in situ SM and temperature performed better than those derived from CLSM estimates relative to SMOS TB observations. The surface salinity for the sampled pixels was on average only 4 mg/g and only locally influenced the TB simulations, when including salinity in the dielectric mixing model of the forward radiative transfer model (RTM) simulations. To explore the potential of retrieving salinity together with other RTM parameters to optimize TB simulations over the entire Dry Chaco, the RTM was inverted using 10 years of multi-angular SMOS TB data and constraints of CLSM SM and temperature. However, the latter modeled SM was not sufficiently accurate and factors such as open surface water were missing in the background constraints, so that the salinity retrievals effectively represented a bulk correction of the dielectric constant, rather than salinity per se. However, the retrieval of vegetation, scattering albedo and surface roughness resulted in realistic values.

L-Band

An analytic model for footprint dispersions and its application to mission design

This is the final report on our recent research activities that are complementary to those conducted by our colleagues, Professor Farrokh Mistree and students, in the context of the Taguchi method. We have studied the mathematical model that forms the basis of the Simulation and Optimization of Rocket Trajectories (SORT) program and developed an analytic method for determining mission reliability with a reduced number of flight simulations. This method can be incorporated in a design algorithm to mathematically optimize different performance measures of a mission, thus leading to a robust and easy-to-use methodology for mission planning and design.

Rao, J. R. Jagannatha

Inference of three-dimensional hot-spot and shell morphology in inertial confinement fusion experiments using a convolutional neural network

The performance of inertial confinement fusion (ICF) implosions is sensitive to the three-dimensional (3D) morphology of the hot-spot and shell configurations. The ability to infer shell-mass uniformity and reconstruct 3D hot spots is crucial for quantifying the degradation of ignition criteria and improving symmetry in ICF implosion experiments. In this work, we present a deep-learning convolutional neural network (CNN) for reconstructing 3D hot-spot and shell structures for ICF capsules. The 3D geometry of the hot spot is reconstructed from x-ray images measured from multiple lines of sight on OMEGA. The shell configuration is inferred indirectly through machine learning using a convolutional neural network extensively trained on a dec3d simulation database. This simulation-dependent approach yields consistent agreement between reconstructed 3D shell densities and machine-learning optimized dec3d simulation results. This work demonstrates a CNN framework that successfully reconstructs 3D capsule structures from two-dimensional images in ICF implosions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY

Multiabsorber Transition-Edge Sensors for X-Ray Astronomy

We are developing arrays of position-sensitive microcalorimeters for future x-ray astronomy applications. These position-sensitive devices commonly referred to as hydras consist of multiple x-ray absorbers, each with a different thermal coupling to a single-transition-edge sensor microcalorimeter. Their development is motivated by a desire to achieve very large pixel arrays with some modest compromise in performance. We report on the design, optimization, and first results from devices with small pitch pixels (<75 μm) being developed for a high-angular and energy resolution imaging spectrometer for Lynx. The Lynx x-ray space telescope is a flagship mission concept under study for the National Academy of Science 2020 decadal survey. Broadband full-width-half-maximum (FWHM) resolution measurements on a 9-pixel hydra have demonstrated ΔE(FWHM) = 2.23 ± 0.14 eV at Al-Kα, ΔE(FWHM) = 2.44 ± 0.29 eV at Mn-Kα, and ΔE(FWHM) = 3.39 ± 0.23 eV at Cu-Kα. Position discrimination is demonstrated to energies below <1 keV and the device performance is well-described by a finite-element model. Results from a prototype 20-pixel hydra with absorbers on a 50-μm pitch have shown ΔE(FWHM) = 3.38 ± 0.20 eV at Cr-Kα1. We are now optimizing designs specifically for Lynx and extending the number of absorbers up to 25/hydra. Numerical simulation suggests optimized designs could achieve ∼3 eV while being compatible with the bandwidth requirements of the state-of-the art multiplexed readout schemes, thus making a 100,000 pixel microcalorimeter instrument a realistic goal.

Smith, Stephen J.

A Review of High-Order and Optimized Finite-Difference Methods for Simulating Linear Wave Phenomena

This paper presents a review of high-order and optimized finite-difference methods for numerically simulating the propagation and scattering of linear waves, such as electromagnetic, acoustic, or elastic waves. The spatial operators reviewed include compact schemes, non-compact schemes, schemes on staggered grids, and schemes which are optimized to produce specific characteristics. The time-marching methods discussed include Runge-Kutta methods, Adams-Bashforth methods, and the leapfrog method. In addition, the following fourth-order fully-discrete finite-difference methods are considered: a one-step implicit scheme with a three-point spatial stencil, a one-step explicit scheme with a five-point spatial stencil, and a two-step explicit scheme with a five-point spatial stencil. For each method studied, the number of grid points per wavelength required for accurate simulation of wave propagation over large distances is presented. Recommendations are made with respect to the suitability of the methods for specific problems and practical aspects of their use, such as appropriate Courant numbers and grid densities. Avenues for future research are suggested.

Zingg, David W.