Search NASA⌕ Search

SEARCH · Search NASA

Results for “simultaneous 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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Design and Manufacturing of High Performance, Reduced Charge Heat Exchangers

Air-to-fluid heat exchangers (HXs) play a critical role as the main heat transfer component in Heating, Ventilation, Air-Conditioning, and Refrigeration (HVAC&R) systems. However, their airside thermal resistance significantly inhibits their overall performance. Furthermore, these HXs must be continually more compact to meet the latest refrigerant charge limits to reduce emissions. Recent literature suggests that traditional HX geometries (e.g., round or flat tubes with fins), have reached their limits, and more sophisticated shape- and topology-optimized designs are required to achieve the next jump in performance. This research sheds light on the next generation of air-to-refrigerant HXs and aims to address several practical issues to commercialization such as novelty challenges (improved performance for significant charge reduction; modeling expertise & time investment), manufacturing challenges (non-round tube manufacturing; tube-header integration; product qualification, e.g., burst pressure testing, extreme operational environment, etc.), and operational challenges (flow maldistribution, fouling & wetting, noise & vibration). For example, a >20% improvement on one (or more) HX-level performance metrics (e.g., envelope volume, airside pressure drop, face area, capacity, refrigerant charge, weight, cost, etc.) must be achieved before a HX design is considered for commercialization. We present a new, comprehensive and experimentally validated air-to-refrigerant HX optimization framework with simultaneous thermal-hydraulic performance and mechanical strength considerations for novel, non-round, shape- and topology-optimized tubes capable of optimizing single and two-phase HX designs for any refrigerant choice and performance requirement with significant engineering time savings compared to conventional design practices. The framework was exercised for a wide range of applications and refrigerants, resulting in HXs which achieved greater than 20% improved performance, 20% reductions in size, and 25% reductions in refrigerant charge. To enable non-round tube bundle use in next generation HVAC&R equipment, novel manufacturing techniques were investigated, including the development of conventional manufacturing methods for small diameter, non-round tubes and novel tube-header integration strategies. In total, ten HX prototypes were manufactured, nine using conventional methods directly attributed to this project and one using advanced additive manufacturing methods. The five-year manufacturing feasibility of the proposed HXs was found to have a good outlook. The non-round tube HX simulated performance was validated through comprehensive experimental testing, including nine in-house component-level tests, one independent component-level test at an industry partner laboratory, and in-house system-level tests of using a commercially-available, residential packaged A/C unit which was retrofitted with a non-round tube prototype HX. It was found that HX designs proposed by the new framework can successfully predict experimental thermal-hydraulic performance within ±10-20% the first time with no manual design changes, eliminating the need for time-consuming and expensive prototyping efforts. This work will accelerate design and time to market for next generation HXs while simultaneously facilitating industry transition to new refrigerants at lower charge.

42 ENGINEERING↗
NASA NTRS2003

Star Identification Using a Triplet Algorithm

The stars observed by onboard star trackers must be identified in order for the star tracker information to be used for attitude determination. Star identification associates the observed stars in the body frame with catalog stars, the positions of which in an inertial frame are well known. If the attitude is approximately known (from less accurate sensors) identification is relatively simple, but if there is no prior attitude knowledge (the Lost in Space problem) star identification is quite difficult. A new, reliable, and accurate star identification algorithm has been developed for use in a batch, ground-based attitude determination system. The identification algorithm uses input from one or more star trackers, each of which can observe three or more stars simultaneously. If observations are available at different times, rate data is used to relate them. The algorithm includes the following steps that are automatically performed: 1) Observations are examined to determine the optimal set of simultaneous star observations in a single tracker. 2) The selected optimal observations are identified using a triplet match. 3) The identified stars in the selected tracker at the selected time are used to compute an attitude at that time. 4) The rate data is refined by determination of biases that minimize the dispersion of repeated observations of the same stars. 5) The refined rate data is used to propagate the attitude to the times of each star observation. 6) The attitude at each time is used to associate all stars observed at that time, in each tracker, with catalog stars. The overall algorithm is very dependable at producing accurate star identification with no initial attitude input.

Hashmall, Joseph A.↗
DOE OSTI2021

Multiobjective Modular Biorefinery Configuration under Uncertainty

With increasing interest in using biomass as a raw material for fuel, the development of biorefineries is an active area of research. However, increasing market competition, uncertainty, and environmental concerns are a few of the challenges that need to be addressed. In this work, mathematical optimization formulations are proposed to address these challenges simultaneously. A multiobjective deterministic optimization framework is proposed to address economic and environmental objectives. A two-stage stochastic optimization framework is proposed to account for uncertainties in material flow and product yields. To model the environmental objective, an environmental risk metric is proposed. Finally, as modularization is gaining popularity in the process industry due to its ability to save capital cost because of standardization and reduce time to market, a formulation for simultaneously achieving modular process design and biorefinery configuration is proposed. Furthermore, the results demonstrate the efficacy of the proposed approach.

09 BIOMASS FUELS↗
DOE OSTI2025

Explainable and Differentiable Reinforcement Learning for Multi-objective Optimization in Particle Accelerators

Operating particle accelerators involves optimizing multiple goals simultaneously, which can be challenging due to trade-offs among objectives. While evolutionary algorithms like the genetic algorithm (GA) have been used for various Multi-Objective Optimization (MOO) tasks, they are not inherently suited for complex control problems. This talk highlights two variations of Reinforcement Learning (RL) for concurrently optimizing heat load and trip rates at the Continuous Electron Beam Accelerator Facility (CEBAF). The problem involves strict constraints on individual states, actions, and overall energy requirements of the beam. First, this talk highlights how differentiability can be harnessed through a Deep Differentiable Reinforcement Learning (DDRL) approach to address MOO issues within particle accelerators. We examine the DDRL method alongside Model Free Reinforcement Learning (MFRL), GA, and Bayesian Optimization (BO). The performance of these methods is assessed by generating a Pareto-front for two objectives. Our findings indicate that DDRL excels in handling high-dimensional problems more effectively than MFRL, BO, and GA. Next, we will show integration of explainable physics-based constraints into RL algorithms to enhance trans- parency and trust in decision-making processes by enabling users to verify that agents adhere to established physical principles. This surrogate function can be modeled using neural networks or sparse dictionary mod- els. By examining the mathematical form of the learned constraint function, we are able to confirm the agent has learned to use the established physics of each environment provided but the surrogate model. In addi- tion, we find that the introduction of a mathematical functional dictionary based surrogate model enables our reinforcement learning algorithms to reliably converge for difficult high-dimensional accelerator controls environments.

Rajput, Kishansingh [Thomas Jefferson National Acc↗
DOE OSTI2026

Bayesian Optimization of Non-Invariant Systems with Constraints Developed for Application to the ECR Ion Source VENUS

In this work, we consider the optimization of non-invariant systems with both safety and control constraints. We present a new approach based on Bayesian optimization for the dynamic, safe and controlled optimization of such systems. Although there are other possible use cases, we focus on the application to the electron cyclotron resonance ion source VENUS. From experimental data, we have observed that VENUS behaves to first order as a non-invariant dynamic system with moving areas of instability. Our novel approach aims at providing a tool that can maintain system optimization in a safe way. This is accomplished by making sure the objective function, the beam current in the case of VENUS, does not fall under an operational minimum, while simultaneously requiring the optimization to avoid areas where VENUS is unstable. We compare the result of our approach on synthetic data modeled to mimic the behavior of VENUS with two methods from the literature, a standard Bayesian optimizer and a safe Bayesian optimizer, both adapted to deal with dynamic systems. A cross Student T-test is conducted to show the significance of the improvement given by the new method we introduce here, regarding the two preexisting methods we compared to. The results of the tests conducted on synthetic data show that the proposed method succeeds at maintaining the system optimized and obeys the predefined constraints better than the literature methods explored.

Bayesian optimization↗
NASA NTRS2015

Coupled Low-thrust Trajectory and System Optimization via Multi-Objective Hybrid Optimal Control

The optimization of low-thrust trajectories is tightly coupled with the spacecraft hardware. Trading trajectory characteristics with system parameters ton identify viable solutions and determine mission sensitivities across discrete hardware configurations is labor intensive. Local independent optimization runs can sample the design space, but a global exploration that resolves the relationships between the system variables across multiple objectives enables a full mapping of the optimal solution space. A multi-objective, hybrid optimal control algorithm is formulated using a multi-objective genetic algorithm as an outer loop systems optimizer around a global trajectory optimizer. The coupled problem is solved simultaneously to generate Pareto-optimal solutions in a single execution. The automated approach is demonstrated on two boulder return missions.

Spacecraft Systems↗
DOE OSTI2024

Manufacturability-based optical design optimization for advanced Kirkpatrick–Baez X-ray focusing mirrors

The advanced Kirkpatrick–Baez (AKB) mirror setup is an effective and compelling solution to provide stable X-ray nano-focusing for synchrotron radiation or free-electron laser beamlines. We propose an AKB mirror design optimization approach to mitigate the difficulties associated with mirror fabrication by minimizing the total slope ranges of the four curved mirrors while achieving the expected focusing performance. In the optimization, we have considered geometry constraints to ensure the beam acceptance with the required clear aperture, the diffraction-limited focal size with the adequate numerical aperture, and the desired mirror gaps for adjustment and the necessary working distance for the sample stage. Additionally, practical constraints linked to mirror metrology and fabrication, such as mirror length limits and curvature uncertainty in measurement, are taken into account. Furthermore, progressive objective optimization eliminates the need for any initial guess, fully automating the AKB optimization process. This approach facilitates the development of an elegant Wolter-I or Wolter-III type AKB design solution that satisfies these multiple constraints. In cases where constraints cannot be simultaneously satisfied, the optimization results provide valuable insights into areas where trade-offs need to be considered. Simulations with ray tracing and wavefront propagation validate the optimized AKB design showing high tolerance to the beam incident angle.

36 MATERIALS SCIENCE↗
NASA NTRSDate not supplied

End-to-End Performance Optimization of a Crewed Lunar Landing Mission Staged from a Near Rectilinear Halo Orbit

A crewed lunar landing mission staged from a Near Rectilinear Halo Orbit (NRHO) presents a unique set of challenges not considered before for the Apollo-type direct missions. The NRHO does not have a fixed orbital plane and its period also changes with time, due to which the mission performance is very sensitive to the timing and location of the NRHO departure and arrival burns for a near-polar surface landing site. For such type of mission staged from the NRHO, an end-to-end optimization approach for the entire mission, from the NRHO departure to surface landing and from surface lift-off to the NRHO insertion, is well-suited for extracting the maximum performance out of this design. In this paper, an integrated approach for performing end-to-end optimization of a NRHO-staged crewed lunar landing mission is proposed that optimizes the in-space and powered flight segments simultaneously using Copernicus trajectory optimization tool. Copernicus is augmented with analytical approximations for the in-space, powered descent, and powered ascent segments in this work to improve the convergence of its multiple-shooting based direct optimization method. New analytical solutions are derived for the different phases of the powered flight segments by solving associated guidance problems. These analytical solutions act as initial guesses for starting the numerical optimization of the integrated mission in Copernicus. An example lunar landing mission is simulated to validate the proposed end-to-end optimization approach and the results are compared with the conventional patched trajectory approach.

Bharat Mahajan↗
DOE OSTI2023

Microgrid Tiered Circuits Effects for a Planned Housing Community in Puerto Rico

Puerto Rico faced a double strike from hurricanes Irma and Maria in 2017. The resulting damage required a comprehensive rebuild of electric infrastructure. There are plans and pilot projects to rebuild with microgrids to increase resilience. This paper provides a techno-economic analysis technique and case study of a potential future community in Puerto Rico that combines probabilistic microgrid design analysis with tiered circuits in building energy modeling. Tiered circuits in buildings allow electric load reduction via remote disconnection of non-critical circuits during an emergency. When coupled to a microgrid, tiered circuitry can reduce the chances of a microgrid's storage and generation resources being depleted. The analysis technique is applied to show 1) Approximate cost savings due to a tiered circuit structure and 2) Approximate cost savings gained by simultaneously considering resilience and sustainability constraints in the microgrid optimization. The analysis technique uses a resistive capacitive thermal model with load profiles for four tiers (tier 1-3 and non-critical loads). Three analyses were conducted using: 1) open-source software called Tiered Energy in Buildings and 2) the Microgrid Design Toolkit. For a fossil fuel based microgrid 30% of the total microgrid costs of 1.18 million USD were calculated where the non-tiered case keeps all loads 99.9% available and the tiered case keeps tier 1 at 99.9%, tier 2 at 95%, tier 3 at 80% availability, with no requirement on non-critical loads. The same comparison for a sustainable microgrid showed 8% cost savings on a 5.10 million USD microgrid due to tiered circuits. The results also showed 6-7% cost savings when our analysis technique optimizes sustainability and resilience simultaneously in comparison to doing microgrid resilience analysis and renewables net present value analysis independently. Though highly specific to our case study, similar assessments using our analysis technique can elucidate value of tiered circuits and simultaneous consideration of sustainability and resilience in other locations.

building energy modeling↗
NASA NTRSDate not supplied

Simultaneous Aerosol and Ocean Polarimeter Products Using Coupled Atmosphere-Ocean Vector Radiative Transfer and Neural Networks: The PACE-MAPP Algorithm

We describe the PACE-MAPP algorithm that simultaneously retrieves aerosol and ocean optical parameters using multiangle and multi-channel polarimeter measurements from the SPEXone, Hyper-Angular Rainbow Polarimeter 2 (HARP2), and Ocean Color Instrument (OCI) instruments onboard the NASA Plankton, Aerosol, Cloud, ocean Ecosystem (PACE) observing system PACE-MAPP is adapted from the Research Scanning Polarimeter (RSP) Microphysical Aerosol Properties from Polarimetry (RSP-MAPP) algorithm. A key feature of the MAPP family of algorithms is the use of a coupled vector radiative transfer model such that the atmosphere and ocean are always considered together as one system. Consequently, conservation of energy ensures that negative water-leaving radiances do not occur. PACE-MAPP uses optimal estimation to simultaneously characterize the optical and microphysical properties of aerosol and ocean constituents, find the optimal solution, and reliably account for the uncertainties of each parameter. This coupled approach, together with multiangle, multi-channel polarimeter measurements, will enable retrievals of aerosol and water properties across the Earth’s oceans. The PACE-MAPP algorithm provides aerosol and ocean products for both the open ocean and coastal areas and is designed to be accurate, modular, and efficient by using fast neural networks that replace the time-consuming vector radiative transfer calculations. We provide an overview of the PACE-MAPP framework and also describe its modular components including its aerosol and hydrosol models, ocean bio-optical models, and thin cirrus model.

Snorre Stamnes↗
NASA NTRS2019

Sizing and Topology Design of an Aeroelastic Wingbox Under Uncertainty

The goals of this work are to use a nested optimizer to conduct simultaneous sizing (inner level) and topology (outer level) design of a wingbox, considering uncertainties in the safety factors used to define the aeroelastic constraints. These uncertainties, propagated via sampling-driven polynomial chaos, are explicitly introduced at the inner level of the method, during gradient-based sizing optimization, resulting in a stochastic optimal sizing distribution. Measures of robustness in the total structural mass are then passed to the outer level, where a global optimizer evolves the topology parameters. The results demonstrate design choices needed to improve robustness in the face of uncertain safety factors, and the various physical mechanisms driving this process.

Stanford, Bret K.↗
DOE OSTI2025

Marine energy supported multi-energy system planning and operation optimization for sustainable coastal community

The growing need for sustainable energy solutions in coastal areas necessitates the development of integrated systems that leverage abundant marine resources. In this study, a standalone Marine Energy Supported Multi-Energy System (MRE-MES) is designed for sustainable coastal community development, utilizing renewable marine resources, including offshore wind, wave, and solar energy, to address the energy needs of electricity, heat, freshwater, and hydrogen. The proposed MRE-MES incorporates a co-optimization model that simultaneously balances capacity planning and operational efficiency to minimize costs and environmental impacts. The system is tested under different renewable energy penetration levels and demand uncertainties, using a two-stage stochastic programming to account for variability in renewable resources and consumption needs. The experimental results indicate that in the optimal system capacity configuration, the percentage of total renewable energy generation is around 80 %, with or without capacity limitation constraints on PV, water tank, and hydrogen storage. Compared to the worst-case scenario in Monte Carlo experiments, two-stage stochastic optimization results in a more robust decision that effectively mitigates the risks posed by future uncertain demand conditions. In conclusion, the findings highlight the viability of marine energy for providing a resilient, comprehensive energy solution to coastal communities.

Capacity planning↗
DOE OSTI2020

Dynamic Modeling of Sequential Service Restoration in Islanded Single Master Microgrids

This paper presents a methodology for generating restoration sequences in distribution systems and microgrids with emphasis on deriving optimal switching times for remote control switches (RCS). Microgrid architecture considered are assumed to be operating in islanded single master mode. Control methods for handling the unbalance inherent in distribution systems using single-phase controlled PQ inverters were developed and incorporated in the proposed method. The problem of considering simultaneously both energy optimization and dynamic stability of the system, which is lacking in existing methods, was addressed by studying the transient behavior of microgrids based on extended EMTP simulations and then approximate and incorporate these behaviors into the optimization formulation for restoration. The approach is based on multi-objective optimization of final restored load and total restoration time. The methodology was studied on a modified IEEE 123 node test feeder. An EMTP simulation in PSCAD of the formed microgrids was used to validate the accuracy of the dynamic and steady-state characteristics of the resulting systems. By considering the transient behavior of DGs during restoration, a better restoration sequence that does not trigger protection relays could be realized.

42 ENGINEERING↗
NASA NTRS1985

Single step optimization of feedback-decoupled spacecraft attitude maneuvers

A next-state optimization technique for monitoring spacecraft maneuvers is examined. The optimization technique permits simultaneous trajectory planning and tracking, incorporates actuator saturation and attitude constraints, and prevents overshoot and undershoot. The global linearization of the dynamics and kinematics for a spacecraft driven by reaction wheels is considered. The equations utilized to derive the spacecraft trajectories are presented. The applications of the next-state optimization procedure to a spacecraft maneuvered by momentum transfer devices and the tracking of the nominal trajectory of a critically damped harmonic oscillator response are described.

Dwyer, T. A. W., III↗
DOE OSTI2025

Safe Reinforcement Learning-Based Transient Stability Control for Islanded Microgrids With Topology Reconfiguration

This paper proposes a safe reinforcement learning (RL)-based transient stability emergency control (TSEC) method for islanded microgrids. RL requires extensive interaction with the environment to learn control strategies, hence, a data-driven approach is used as a substitute for time-consuming time-domain simulation calculations. Deep sigma point processes (DSPP), which is a Gaussian process model, is utilized to predict the normal distribution of transient stability of microgrids and to construct a transient stability chance constraint. Reward-constrained policy optimization (RCPO) can simultaneously achieve objective prediction, policy learning, and constraint cost coefficient update across multiple timescales. RCPO interacts with the DSPP-based microgrid environment through a multi-process parallel manner, greatly increasing the training speed. Case studies on a real islanded microgrid demonstrate that the proposed method can efficiently and quickly obtain the optimal emergency control strategy while adhering to all hard constraints.

14 SOLAR ENERGY↗
DOE OSTI2022

Integrated Optimization and Control of a Hybrid Gas Turbine/sCO 2 Power System

During phase-I, the project team led by Echogen Power Systems (EPS) had two primary objectives based on investigating the application of gas turbines with supercritical carbon dioxide (sCO 2 ) power cycles. The first objective was to improve the overall efficiency and performance of a hybrid gas turbine/sCO 2 power system through a joint optimization of the two subsystems (gas turbine and sCO 2 power cycle) using non-linear optimization techniques that simultaneously evaluate thermal performance of the combined cycle. The hybrid power system included several points of interaction, including (but not limited to) gas turbine exhaust, fuel heating, inlet chilling and turbine cooling. The second objective was to establish a baseline transient response model of the hybrid power system and a notional microgrid and begin steps to integrate the control systems of the three major elements (gas turbine, sCO 2 cycle and grid controller). The project team established a baseline performance for a combined cycle power plant using a production gas turbine and scaled sCO 2 power cycle only utilizing exhaust heat recovery. Echogen’s non-linear techno-economic optimization code was extended by adding gas turbine component models derived from a in-house developed gas turbine design code. With the two cycles coupled by the gas turbine exhaust, design parameters of both cycles were allowed to vary simultaneously to determine performance opportunity versus isolated designs. Returning to the baseline gas turbine/sCO 2 power cycle transient models: Echogen had in-house developed sCO 2 cycle transient model in GT-Suite system simulation software, and had partnered with Siemens Finspång for gas turbine transient model, and Siemens PTI group to provide micro-grid load profile as well as hybrid power cycle generated load (power and frequency) analysis. The transient model for the SGT-750 Siemens gas turbine was a “black-box” functional mock-up interface (FMI) model developed by Siemens Industrial Turbomachinery in Finspång, Sweden. The SGT-750 is a twin-shaft gas turbine that produces 40 MW electricity with an efficiency of about 40% at ISO conditions. At 100% gas turbine throttle (load), the SGT-750 has average exhaust conditions of 114.6 kg/s and 469.8°C. The transient model for sCO 2 power cycle was developed by Echogen in GT-SUITE 1D system simulation software platform. The basic CO 2 flow circuit has single-shaft turbomachinery with net 11.5 MW electrical power output at design conditions. The power turbine has a double-ended shaft with one end connected to synchronous generator through a fixed-ratio gearbox. The other end of power turbine is connected to the compressor through a continuously variable transmission. The major components of the sCO 2 power cycle modeled include air cooled condenser/cooler, CO 2 compressor, recuperator, two waste heat exchanger coils, power turbine, continuous variable transmission, gearbox and generator. Integration of SGT-750 transient model and sCO 2 power cycle transient model was done in Matlab Simulink. In the integrated model, the gas turbine and sCO 2 power cycle interacted at two points, first one being the gas turbine exhaust gas flow rate and temperature, which were inputs to sCO 2 power cycle model. The second point was the distribution of micro-grid load demand signal between the SGT-750 generator and sCO 2 cycle generator. For a given combined-cycle load demand, the gas turbine load demand was equal to the total demand minus the sCO 2 cycle power generated. In the present study the integrated model was simulated for two cases of grid load demand: (i) for a step change, both positive-step and negative-step, in grid load demand (ii) for a micro-grid load demand curve provided by Siemens PTI group. Finally, the time series plots representing load demand versus integrated system response were presented including the sCO 2 power cycle control system performance plots. The actual generated power and frequency of both the generators, gas turbine and sCO 2 power cycle, was supplied to Siemens PTI group for dynamic grid assessment, results of which are provided in appendices.

03 NATURAL GAS↗
DOE OSTI2025

An Adaptive Newton-Based Free-Boundary Grad–Shafranov Solver

Equilibria in magnetic confinement devices result from force balancing between the Lorentz force and the plasma pressure gradient. In an axisymmetric configuration like a tokamak, such an equilibrium is described by an elliptic equation for the poloidal magnetic flux, commonly known as the Grad–Shafranov equation. It is challenging to develop a scalable and accurate free-boundary Grad–Shafranov solver, since it is a fully nonlinear optimization problem that simultaneously solves for the magnetic field coil current outside the plasma to control the plasma shape. In this work, we develop a Newton-based free-boundary Grad–Shafranov solver using adaptive finite elements and preconditioning strategies. The free-boundary interaction leads to the evaluation of a domain-dependent nonlinear form of which its contribution to the Jacobian matrix is achieved through shape calculus. The optimization problem aims to minimize the distance between the plasma boundary and specified control points while satisfying two nontrivial constraints, which correspond to the nonlinear finite element discretization of the Grad–Shafranov equation and a constraint on the total plasma current involving a nonlocal coupling term. The linear system is solved by a block factorization, and AMG is called for subblock elliptic operators. The unique contributions of this work include the treatment of a global constraint, preconditioning strategies, nonlocal reformulation, and the implementation of adaptive finite elements. Furthermore, it is found that the resulting Newton solver is robust, successfully reducing the nonlinear residual to 1e-6 and lower in a small handful of iterations while addressing the challenging case to find a Taylor state equilibrium where conventional Picard-based solvers fail to converge.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Implicit shock tracking for unsteady flows by the method of lines

Here, a recently developed high-order implicit shock tracking (HOIST) framework for resolving discontinuous solutions of inviscid, steady conservation laws is extended to the unsteady case. Central to the framework is an optimization problem which simultaneously computes a discontinuity-aligned mesh and the corresponding high-order approximation to the flow, which provides nonlinear stabilization and a high-order approximation to the solution. This work extends the implicit shock tracking framework to the case of unsteady conservation laws using a method of lines discretization via a diagonally implicit Runge-Kutta method by “solving a steady problem at each timestep”. We formulate and solve an optimization problem that produces a feature-aligned mesh and solution at each Runge-Kutta stage of each timestep, and advance this solution in time by standard Runge-Kutta update formulas. A Rankine-Hugoniot based prediction of the shock location together with a high-order, untangling mesh smoothing procedure provides a high-quality initial guess for the optimization problem at each time, which results in rapid convergence of the sequential quadratic programing (SQP) optimization solver. This method is shown to deliver highly accurate solutions on coarse, high-order discretizations without nonlinear stabilization and recover the design accuracy of the Runge-Kutta scheme. We demonstrate this framework on a series of inviscid, unsteady conservation laws in both one- and two- dimensions. We also verify that our method is able to recover the design order of accuracy of our time integrator in the presence of a strong discontinuity.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗