Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “optimal start”

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

Integration of Low-Fidelity MDO and CFD-Based Redesign of Low-Boom Supersonic Transports

A mixed-fidelity low-boom multidisciplinary optimization (MDO) problem is formulated for integration of low-fidelity MDO and computational fluid dynamics (CFD) based low-boom inverse design optimization. The mixed-fidelity low-boom MDO problem aims to enforce the weight consistency for CFD-based low-boom inverse design optimization: the optimum low-boom configuration is designed for the weight at the start of cruise of the same configuration for the low-boom overland mission. Moreover, it also seeks the optimum trades among the maximum takeoff gross weight (MTOGW), cruise Mach, and range for the low-boom overland mission while meeting the requirement of flying a transatlantic overwater mission. A block coordinate optimization (BCO) method is developed to find an approximate solution of the mixed-fidelity low-boom MDO problem. The BCO method is successfully applied to generate two CFD-based low-boom configurations that closely match two reversed equivalent area targets with ground noise levels below 70 PLdB, respectively. Moreover, these CFD-based low-boom configurations can be obtained with minor wing modifications from the solutions of the low-fidelity MDO problem. The low-fidelity MDO solutions are Pareto points for constrained multiobjective optimization of MTOGW, low-boom cruise Mach, and low-boom range. The best low-boom concept has a predetermined fuselage shape tailored for passengers and main gear storage, carries 40 passengers at seat pitch of 48 in, flies a low-boom overland mission with cruise Mach of 1.8 and range of 2,950 nm, cruises overwater at Mach 1.8 with range of 3,600 nm, satisfies the specified constraints for landing/cruise/takeoff, has MTOGW of 145,164 lb, trims the low-boom cruise flight with fuel redistributions instead of control surface deflections, and has a reversed equivalent area distribution closely matching a target with ground noise level below 70 PLdB.

low-boom supersonic transports↗

Mesh refinement in finite element analysis by minimization of the stiffness matrix trace

Most finite element packages provide means to generate meshes automatically. However, the user is usually confronted with the problem of not knowing whether the mesh generated is appropriate for the problem at hand. Since the accuracy of the finite element results is mesh dependent, mesh selection forms a very important step in the analysis. Indeed, in accurate analyses, meshes need to be refined or rezoned until the solution converges to a value so that the error is below a predetermined tolerance. A-posteriori methods use error indicators, developed by using the theory of interpolation and approximation theory, for mesh refinements. Some use other criterions, such as strain energy density variation and stress contours for example, to obtain near optimal meshes. Although these methods are adaptive, they are expensive. Alternatively, a priori methods, until now available, use geometrical parameters, for example, element aspect ratio. Therefore, they are not adaptive by nature. An adaptive a-priori method is developed. The criterion is that the minimization of the trace of the stiffness matrix with respect to the nodal coordinates, leads to a minimization of the potential energy, and as a consequence provide a good starting mesh. In a few examples the method is shown to provide the optimal mesh. The method is also shown to be relatively simple and amenable to development of computer algorithms. When the procedure is used in conjunction with a-posteriori methods of grid refinement, it is shown that fewer refinement iterations and fewer degrees of freedom are required for convergence as opposed to when the procedure is not used. The mesh obtained is shown to have uniform distribution of stiffness among the nodes and elements which, as a consequence, leads to uniform error distribution. Thus the mesh obtained meets the optimality criterion of uniform error distribution.

Kittur, Madan G.↗

Machine Learning Techniques in Optimal Design

Many important applications can be formalized as constrained optimization tasks. For example, we are studying the engineering domain of two-dimensional (2-D) structural design. In this task, the goal is to design a structure of minimum weight that bears a set of loads. A solution to a design problem in which there is a single load (L) and two stationary support points (S1 and S2) consists of four members, E1, E2, E3, and E4 that connect the load to the support points is discussed. In principle, optimal solutions to problems of this kind can be found by numerical optimization techniques. However, in practice [Vanderplaats, 1984] these methods are slow and they can produce different local solutions whose quality (ratio to the global optimum) varies with the choice of starting points. Hence, their applicability to real-world problems is severely restricted. To overcome these limitations, we propose to augment numerical optimization by first performing a symbolic compilation stage to produce: (a) objective functions that are faster to evaluate and that depend less on the choice of the starting point and (b) selection rules that associate problem instances to a set of recommended solutions. These goals are accomplished by successive specializations of the problem class and of the associated objective functions. In the end, this process reduces the problem to a collection of independent functions that are fast to evaluate, that can be differentiated symbolically, and that represent smaller regions of the overall search space. However, the specialization process can produce a large number of sub-problems. This is overcome by deriving inductively selection rules which associate problems to small sets of specialized independent sub-problems. Each set of candidate solutions is chosen to minimize a cost function which expresses the tradeoff between the quality of the solution that can be obtained from the sub-problem and the time it takes to produce it. The overall solution to the problem, is then obtained by solving in parallel each of the sub-problems in the set and computing the one with the minimum cost. In addition to speeding up the optimization process, our use of learning methods also relieves the expert from the burden of identifying rules that exactly pinpoint optimal candidate sub-problems. In real engineering tasks it is usually too costly to the engineers to derive such rules. Therefore, this paper also contributes to a further step towards the solution of the knowledge acquisition bottleneck [Feigenbaum, 1977] which has somewhat impaired the construction of rulebased expert systems.

Cerbone, Giuseppe↗

Power System Recovery Coordinated with (Non-)Black-Start Generators

Power restoration is an urgent task after a black-out, and recovery efficiency is critical when quantifying system resilience. Multiple elements should be considered to restore the power system quickly and safely. This paper proposes a recovery model to solve a direct-current optimal power flow (DCOPF) based on mixed-integer linear programming (MILP). Since most of the generators cannot start independently, the interaction between black-start (BS) and non-black-start (NBS) generators must be modeled appropriately. The energization status of the NBS is coordinated with the recovery status of transmission lines, and both of them are modeled as binary variables. Also, only after an NBS unit receives the cranking power through connected transmission lines, will it be allowed to participate in the following system dispatch. The amount of cranking power is estimated as a fixed proportion to the maximum generation capacity. The proposed model is validated on several test systems, as well as a 1393-bus representation system of the Puerto Rican electric power grid. Test results demonstrate how the recovery of NBS units and damaged transmission lines can be optimized, resulting in an efficient and well-coordinated recovery procedure.

Zhao, Meng↗

A Control Model: Interpretation of Fitts' Law

The analytical results for several models are given: a first order model where it is assumed that the hand velocity can be directly controlled, and a second order model where it is assumed that the hand acceleration can be directly controlled. Two different types of control-laws are investigated. One is linear function of the hand error and error rate; the other is the time-optimal control law. Results show that the first and second order models with the linear control-law produce a movement time (MT) function with the exact form of the Fitts' Law. The control-law interpretation implies that the effect of target width on MT must be a result of the vertical motion which elevates the hand from the starting point and drops it on the target at the target edge. The time optimal control law did not produce a movement-time formula simular to Fitt's Law.

Connelly, E. M.↗

Quantifying Northern High Latitude Gross Primary Productivity (GPP) Using Carbonyl Sulfide (OCS)

The northern high latitude (NHL, 40°N to 90°N) is where the second peak region of gross primary productivity (GPP) other than the tropics. The summer NHL GPP is about 80% of the tropical peak, but both regions are still highly uncertain. Carbonyl sulfide (OCS) provides an important proxy for photosynthetic carbon uptake. Here we optimize the OCS plant uptake fluxes across the NHL by fitting atmospheric concentration simulation with the GEOS-CHEM global transport model to the aircraft profiles acquired over Alaska during NASA's Carbon in Arctic Reservoirs Vulnerability Experiment (2012–2015). We use the empirical biome-specific linear relationship between OCS plant uptake flux and GPP to derive the six plant uptake OCS fluxes from different GPP data. Such GPP-based fluxes are used to drive the concentration simulations. We evaluate the simulations against the independent observations at two ground sites of Alaska. The optimized OCS fluxes suggest the NHL plant uptake OCS flux of -247 Gg S year -1 , about 25% stronger than the ensemble mean of the six GPP-based OCS fluxes. GPP-based OCS fluxes systematically underestimate the peak growing season across the NHL, while a subset of models predict early start of season in Alaska, consistent with previous studies of net ecosystem exchange. The OCS optimized GPP of 34 PgC yr -1 for NHL is also about 25% more than the ensembles mean from six GPP data. Further work is needed to fully understand the environmental and biotic drivers and quantify their rate of photosynthetic carbon uptake in Arctic ecosystems.

54 ENVIRONMENTAL SCIENCES↗

Developing A Continuous Ozone Record Through the SAGE and Aura Missions With NASA Reanalysis Products

During the last quarter of the 20th century, the Stratospheric Aerosol and Gas Experiment (SAGE) missions were crucial in monitoring the loss and subsequent recovery of the stratospheric ozone layer. Due to the employed solar occultation and self-calibration method, the SAGE monitors have produced stable data throughout the lifetime of each instrument. However, over ten years passed between the end of the SAGE II and SAGE III/M3M missions in 2005 and the launch of SAGE III/ISS instrument in 2017, leaving a gap in the data that must be bridged in order to assess trends in the ozone record. The Modern-Era Retrospective analysis for Research and Applications, version 2 (MERRA-2) reanalysis product, with output available starting in 1980, is an attractive candidate for trend analysis due to the statistically optimized combination of multiple observing systems and the regular temporal and spatial coverage. However, changes in the assimilated observation systems can introduce discontinuities within the MERRA-2 ozone record, such as in 2004 when the MERRA-2 system shifted from assimilating ozone retrievals collected by SBUV instruments to those collected by instruments onboard the Aura satellite. In this study, we explore using the SAGE II record as a transfer function to develop a stable reanalysis data product, suitable for trend analysis, from the start of the SAGE II record in 1984 through the present. We follow the procedure outlined by Wargan et al. (2018) to address discontinuities in the MERRA-2 ozone dataset at the 2004 transition and during the Aura record. SAGE II ozone profiles are used to correct discontinuities in upper stratospheric ozone associated with changes in the MERRA-2 meteorological observing system in 1998 and 1995. We will then assess the relative performance of the data from different SAGE sensors using the resulting bias-corrected MERRA-2 ozone fields.

Pamela Wales↗

Tensor decompositions for count data that leverage stochastic and deterministic optimization

There is growing interest to extend low-rank matrix decompositions to multi-way arrays, or tensors. One fundamental low-rank tensor decomposition is the canonical polyadic decomposition (CPD). The challenge of fitting a low-rank, nonnegative CPD model to Poisson-distributed count data is of particular interest. Several popular algorithms use local search methods to approximate the maximum likelihood estimator (MLE) of the Poisson CPD model. Here, this work presents two new algorithms that extend state-of-the-art local methods for Poisson CPD. Hybrid GCP-CPAPR combines Generalized Canonical Decomposition (GCP) with stochastic optimization and CP Alternating Poisson Regression (CPAPR), a deterministic algorithm, to increase the probability of converging to the MLE over either method used alone. Restarted CPAPR with SVDrop uses a heuristic based on the singular values of the CPD model unfoldings to identify convergence toward optimizers that are not the MLE and restarts within the feasible domain of the optimization problem, thus reducing overall computational cost when using a multi-start strategy. We provide empirical evidence that indicates our approaches outperform existing methods with respect to converging to the Poisson CPD MLE.

CPAPR↗

Projected sensitivity of the LUX-ZEPLIN experiment to the two-neutrino and neutrinoless double $\beta$ decays of $^{134}$Xe

The projected sensitivity of the LUX-ZEPLIN (LZ) experiment to two-neutrino and neutrinoless double beta decay of 134 Xe is presented. LZ is a 10-tonne xenon time projection chamber optimized for the detection of dark matter particles, that is expected to start operating in 2021 at Sanford Underground Research Facility, USA. Its large mass of natural xenon provides an exceptional opportunity to search for the double beta decay of 134 Xe, for which xenon detectors enriched in 136 Xe are less effective. For the two-neutrino decay mode, LZ is predicted to exclude values of the half-life up to 1.7 x 10 24 years at 90% confidence level (CL), and has a three-sigma observation potential of 8.7 x 10 23 years, approaching the predictions of nuclear models. For the neutrinoless decay mode LZ, is projected to exclude values of the half-life up to 7.3 x 10 24 years at 90% CL.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

The Heuristic Dynamic Programming Approach in Boost Converters

In this study, a heuristic dynamic programming controller is proposed to control a boost converter. Conventional controllers such as proportional-integral-derivative (PID) or proportional-integral (PI) are designed based on the linearized small-signal model near the operating point. Therefore, the performance of the controller during the start-up, the load change, or the input voltage variation is not optimal since the system model changes by varying the operating point. The heuristic dynamic programming controller optimally controls the boost converter by following the approximate dynamic programming. The advantage of the HDP is that the neural network-based characteristic of the proposed controller enables boost converters to easily cope with large disturbances. An HDP with a well-trained critic and action networks can perform as an optimal controller for the boost converter. To compare the effectiveness of the traditional PI-based and the HDP boost converter, the simulation results are provided.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Deep reinforcement learning with online data augmentation to improve sample efficiency for intelligent HVAC control

Deep Reinforcement Learning (DRL) has started showing success in real-world applications such as building energy optimization. Much of the research in this space utilized simulated environments to train RL-agent in an offline mode. Very few research have used DRL-based control in real-world systems due to two main reasons: 1) sample efficiency challenge---DRL approaches need to perform a lot of interactions with the environment to collect sufficient experiences to learn from, which is difficult in real systems, and 2) comfort or safety related constraints---user's comfort must never or at least rarely be violated. In this work, we propose a novel deep Reinforcement Learning framework with online Data Augmentation (RLDA) to address the sample efficiency challenge of real-world RL. We used a time series Generative Adversarial Network (TimeGAN) architecture as a data generator. We further evaluated the proposed RLDA framework using a case study of an intelligent HVAC control. With a ≈28% improvement in the sample efficiency, RLDA framework lays the way towards increased adoption of DRL-based intelligent control in real-world building energy management systems.

Kurte, Kuldeep↗

Recovery Simulator and Analysis Formulation: Mathematical Framework for Enhanced Resilience and Resource Allocation

This report introduces recovery simulator and analysis (RSA), a framework aimed at enhancing the resilience of electrical grids post-disruption. The RSA model leverages an optimization problem formulation that focuses on maximizing the load served (or optionally customers served) through a coordinated and cooptimized recovery of non-black start generation, transmission lines, feeders and substations subject to labor budget constraints. By integrating advanced linear programming techniques, the simulator selects efficient reocovery pathways, optimizing both short-term and long-term grid recovery strategies. The mathematical framework guides decision-making through a comprehensive evaluation of potential recovery actions, factoring in the trade-offs between labor constraints and load (or optionally customer) restoration efficacy. This enables grid operators to simulate diverse outage scenarios and delineate optimal recovery pathways, thereby prioritizing critical repair tasks and ensuring resource allocation is both economical and effective. The intended use case of RSA is to allow planners to explore many recovery scenarios quickly and determine assets most critical across a wide range of scenarios, and therefore strong candidates for hardening or additional investment. RSA might also be used in an operational setting, following a single event, for exploring efficient recovery pathways.

24 POWER TRANSMISSION AND DISTRIBUTION↗

The Newton Modified Barrier Method for QP Problems

The Modified Barrier Functions (MBF) have elements of both Classical Lagrangians (CL) and Classical Barrier Functions (CBF). The MBF methods find an unconstrained minimizer of some smooth barrier function in primal space and then update the Lagrange multipliers, while the barrier parameter either remains fixed or can be updated at each step. The numerical realization of the MBF method leads to the Newton MBF method, where the primal minimizer is found by using Newton's method. This minimizer is then used to update the Lagrange multipliers. In this paper, we examine the Newton MBF method for the Quadratic Programming (QP) problem. It will be shown that under standard second-order optimality conditions, there is a ball around the primal solution and a cut cone in the dual space such that for a set of Lagrange multipliers in this cut cone, the method converges quadratically to the primal minimizer from any point in the aforementioned ball, and continues to do so after each Lagrange multiplier update. The Lagrange multipliers remain within the cut cone and converge linearly to their optimal values. Any point in this ball will be called a "hot start". Starting at such a "hot start", at most Omicron(1n 1n epsilon(exp -1)) Newton steps are sufficient to perform the primal minimization which is necessary for the Lagrange multiplier update. Here, epsilon > 0 is the desired accuracy. Because of the linear convergence of the Lagrange multipliers, this means that only Omicron(1n epsilon(exp -1))omicron(ln 1n epsilon(exp-1)) Newton steps are required to reach an epsilon-approximation to the solution from any "hot start". In order to reach the "hot start", one has to perform Omicron(square root(m) 1n C) Newton steps, where m characterizes the size of the problem and C > 0 is the condition number of the QP problem. This condition number will be characterized explicitly in terms of key parameters of the QP problem, which in turn depend on the input data and the size of the problem.

Melman, A.↗

A review of some exact solutions to the planar equations of motion of a thrusting spacecraft

With the complexities in computing optimal low thrust trajectories, easily-computed, good sub-optimal trajectories provide both a practical alternative for mission designers and a starting point for optimisation. The present paper collects in one place for easy reference and comparison several exact solutions that have been obtained in the literature over the last few decades: the logarithmic spiral, Pinkham's variant thereof, Forbes spiral, the exponential sinusoid, the case of constant radial thrust, Markopoulos's Keplerian thrust arcs, Lawden's spiral, and the analogous Bishop and Azimov spiral.

low-thrust trajectories electric propulsion solar ↗

Making High-Tensile-Strength Amalgam Components

Structural components made of amalgams can be made to have tensile strengths much greater than previously known to be possible. Amalgams, perhaps best known for their use in dental fillings, have several useful attributes, including room-temperature fabrication, corrosion resistance, dimensional stability, and high compressive strength. However, the range of applications of amalgams has been limited by their very small tensile strengths. Now, it has been discovered that the tensile strength of an amalgam depends critically on the sizes and shapes of the particles from which it is made and, consequently, the tensile strength can be greatly increased through suitable choice of the particles. Heretofore, the powder particles used to make amalgams have been, variously, in the form of micron-sized spheroids or flakes. The tensile reinforcement contributed by the spheroids and flakes is minimal because fracture paths simply go around these particles. However, if spheroids or flakes are replaced by strands having greater lengths, then tensile reinforcement can be increased significantly. The feasibility of this concept was shown in an experiment in which electrical copper wires, serving as demonstration substitutes for copper powder particles, were triturated with gallium by use of a mortar and pestle and the resulting amalgam was compressed into a mold. The tensile strength of the amalgam specimen was then measured and found to be greater than 10(exp 4) psi (greater than about 69 MPa). Much remains to be done to optimize the properties of amalgams for various applications through suitable choice of starting constituents and modification of the trituration and molding processes. The choice of wire size and composition are expected to be especially important. Perusal of phase diagrams of metal mixtures could give insight that would enable choices of solid and liquid metal constituents. Finally, whereas heretofore, only binary alloys have been considered for amalgams, ternary additions to liquid or solid components should be considered as means to impart desired properties to amalgams.

Grugel, Richard↗

New Horizons in Cosmology with Spectral Distortions of the Cosmic Microwave Background

Following the pioneering observations with COBE in the early 1990s, studies of the cosmic mi- crowave background (CMB) have primarily focused on temperature and polarization anisotropies. CMB spectral distortions – tiny departures of the CMB energy spectrum from that of a perfect blackbody – provide a second, independent probe of fundamental physics, with a reach deep into the primordial Universe. The theoretical foundation of spectral distortions has seen major advances in recent years, highlighting the immense potential of this emerging field. Spectral distortions probe a fundamental property of the Universe – its thermal history – thereby providing additional insight into processes within the cosmological standard model(I) (CSM) as well as new physics beyond. Spectral distortions are an important tool for understanding inflation and the nature of dark matter. They shed new light on the physics of recombination and reionization, both prominent stages in the evolution of our Universe, and furnish critical information on baryonic feedback processes, in addition to probing primordial correlation functions at scales inaccessible to other tracers. In principle the range of signals is vast: many orders of magnitude of discovery space can be explored by detailed observations of the CMB energy spectrum. Several CSM signals are predicted and provide clear experimental targets that are observable with present-day technology. Confirmation of these signals would extend the reach of the CSM by orders of magnitude in physical scale as the Universe evolves from the initial stages to its present form. Their absence would pose a huge theoretical challenge, immediately pointing to new physics. Here, we advocate for a dedicated effort to measure CMB spectral distortions at the largest angular scales (greater than approximately 1°) within the ESA Voyage 2050 Program. We argue that an L-class mission with a pathfinder would allow a precise measurement of all the expected CSM distortions. With an M-class mission, the primordial distortions (created at z >~ 10(exp 3)) would still be detected at modest significance, while the late-time distortions will continue to be measured to high accuracy. Building on the heritage of COBE/FIRAS, a spectrometer that consists of multiple, cooled (approximately equal to 0.1 K), absolutely-calibrated Fourier Transform Spectrometers (FTS) with wide frequency coverage (ν approximately equal to 10 GHz to a few x THz) and all-sky spectral sensitivity at the level of 0.1 0.5 Jy/sr would be the starting point for the M-class option. A scaled and further optimized version of this concept is being envisioned as the L-class option. Such measurements can only be done from space and would deliver hundreds of absolutely-calibrated maps of the Universe at large scales, opening numerous science opportunities for cosmology and astrophysics. This will provide independent probes of inflation, dark matter and particle physics, recombination and the energy output of our Universe from at late times, turning the long-standing spectral distortion limits of COBE/FIRAS into clear detections.

Chluba, J.↗

Cold Weather Impacts on Electric School Buses

School bus fleets considering electric school buses (ESBs) can start here to learn how cold weather impacts ESBs and considerations for optimizing bus performance when temperatures drop.

ADVANCED PROPULSION SYSTEMS↗

Intelligent driving passive pedal control

A method for assignment of vehicle control includes receiving route data indicating a route between a starting location of a vehicle and a destination location, and determining an optimal vehicle configuration for the route based on a target vehicle speed and a hybrid torque split. The method further includes receiving a driver requested torque value and determining a passive pedal torque value based on the route data and vehicle powertrain data. The method further includes selectively assigning control of the vehicle to a vehicle system or to a driver of the vehicle based on the driver requested torque value and the passive pedal torque value.

Aggoune, Karim↗