Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “linear constraints”

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 793 records · Page 44
All sourcesNASA NTRSDOE OSTI
NASA NTRS1985

Theoretical and experimental studies in support of the geophysical fluid flow experiment

Meteorologists and astrophysicists interested in large scale planetary and solar circulations have come to recognize the importance of rotation and stratification in determining the character of these flows. In particular, the effect of latitude-dependent Coriolis force on nonlinear convection is thought to play a crucial role in such phenomena as differential rotation on the Sun, cloud band orientation on Jupiter, and the generation of magnetic fields in thermally driven dynamos. The continuous low-gravity environment of the orbiting space shuttle offers a unique opportunity to make laboratory studies of such large-scale thermally driven flows under the constraint imposed by rotation and sphericity. This is possible because polarization forces in a dielectric liquid, which are linearly dependent on fluid temperature, give rise to an effectively radial buoyancy force when a radial electrostatic field is imposed. The Geophysical Fluid Flow Cell (GFFC) is an implementation of this ideal in which fluid is contained between two rotating hemispheres that are differentially heated and stressed with a large a-c voltage. The experiment, to be flown on Spacelab III (currently set for launch April 29, 1985), will explore non-linear mode selection and high Rayleigh number turbulence in a rotating convecting spherical shell of liquid. Experiments will be carried out in a low driving parameter range where some limited numerical experimentation is currently feasible, as well as in a parameter range significantly beyond numerical computation for many years.

Hart, J.↗
NASA NTRS1989

Stress analysis method for clearance-fit joints with bearing-bypass loads

Within a multi-fastener joint, fastener holes may be subjected to the combined effects of bearing loads and loads that bypass the hole to be reacted elsewhere in the joint. The analysis of a joint subjected to search combined bearing and bypass loads is complicated by the usual clearance between the hole and the fastener. A simple analysis method for such clearance-fit joints subjected to bearing-bypass loading has been developed in the present study. It uses an inverse formulation with a linear elastic finite-element analysis. Conditions along the bolt-hole contact arc are specified by displacement constraint equations. The present method is simple to apply and can be implemented with most general purpose finite-element programs since it does not use complicated iterative-incremental procedures. The method was used to study the effects of bearing-bypass loading on bolt-hole contact angles and local stresses. In this study, a rigid, frictionless bolt was used with a plate having the properties of a quasi-isotropic graphite/epoxy laminate. Results showed that the contact angle as well as the peak stresses around the hole and their locations were strongly influenced by the ratio of bearing and bypass loads. For single contact, tension and compression bearing-bypass loading had opposite effects on the contact angle. For some compressive bearing-bypass loads, the hole tended to close on the fastener leading to dual contact. It was shown that dual contact reduces the stress concentration at the fastener and would, therefore, increase joint strength in compression. The results illustrate the general importance of accounting for bolt-hole clearance and contact to accurately compute local bolt-hole stresses for combined bearings and bypass loading.

Naik, R. A.↗
NASA NTRS1991

Two alternative ways for solving the coordination problem in multilevel optimization

Two techniques for formulating the coupling between levels in multilevel optimization by linear decomposition, proposed as improvements over the original formulation, now several years old, that relied on explicit equality constraints which were shown by application experience as occasionally causing numerical difficulties. The two new techniques represent the coupling without using explicit equality constraints, thus avoiding the above diffuculties and also reducing computational cost of the procedure. The old and new formulations are presented in detail and illustrated by an example of a structural optimization. A generic version of the improved algorithm is also developed for applications to multidisciplinary systems not limited to structures.

Sobieszczanski-Sobieski, Jaroslaw↗
NASA NTRS1996

Efficient Low-Speed Flight in a Wind Field

A new software tool was needed for flight planning of a high altitude, low speed unmanned aerial vehicle which would be flying in winds close to the actual airspeed of the vehicle. An energy modeled NLP (non-linear programming) formulation was used to obtain results for a variety of missions and wind profiles. The energy constraint derived included terms due to the wind field and the performance index was a weighted combination of the amount of fuel used and the final time. With no emphasis on time and with no winds the vehicle was found to fly at maximum lift to drag velocity, V(sub md). When flying in tail winds the velocity was less than V(sub md), while flying in head winds the velocity was higher than V(sub md). A family of solutions was found with varying times of flight and varying fuel amounts consumed which will aid the operator in choosing a flight plan depending on a desired landing time. At certain parts of the flight, the turning terms in the energy constraint equation were found to be significant. An analysis of a simpler vertical plane cruise optimal control problem was used to explain some of the characteristics of the vertical plane NLP results.

Feldman, Michael A.↗
DOE OSTI2015

A Bound- and Equality-Constrained Quadratic Programming Method for Image Reconstruction

Image deblurring is modeled as an ill-posed, linear inverse problem and can be written in terms of an optimization problem. By adding bound and equality constraints and/or penalty terms to the objective function, the problem is regularized and a viable solution exists; however, there are few methods available for solving such constrained problems. Here we formulate two regularized solutions into constrained optimization problems and adapt an iterative augmented Lagrangian solver. Results are shown on synthetic data in one- and two-dimensions, as well as on a radiograph used to calibrate the transmission curve of a pulsed power X-ray source at a U.S. Department of Energy radiography facility.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2021

Automating Discovery of Physics-Informed Neural State Space Models via Learning and Evolution

Recent works exploring deep learning application to dynamical systems modeling have demonstrated that embedding physical priors into neural networks can yield more effective, physically-realistic, and data-efficient models. However, in the absence of complete prior knowledge of a dynamical system's physical characteristics, determining the optimal structure and optimization strategy for these models can be difficult. In this work, we explore methods for discovering neural state space dynamics models for system identification. Starting with a design space of block-oriented state space models and structured linear maps with strong physical priors, we encode these components into a model genome alongside network structure, penalty constraints, and optimization hyperparameters. Demonstrating the overall utility of the design space, we employ an asynchronous genetic search algorithm that alternates between model selection and optimization and obtains accurate physically consistent models of three physical systems: an aerodynamics body, a continuous stirred tank reactor, and a two tank interacting system.

genetic algorithms, neural architecture search, ne↗
DOE OSTI2014

SU-E-T-234: Modulated Photon Radiotherapy (XMRT):The Impact of Incorporating Energy Modulation Into Intensity Modulated Radiotherapy (IMRT) Optimization

Purpose: To develop a new radiotherapy plan optimization technique that, for a given organ geometry, will find the optimal photon beam energies and fluences to produce a desirable dose distribution. This new modulated (both in energy and fluence) photon radiotherapy (XMRT) was compared with intensity modulated radiotherapy (IMRT) for a simple organ geometry. Methods: The XMRT optimization was formulated using a linear programming approach where the objective function is the mean dose to the healthy organs and dose-point constraints were assigned to each organ of interest. The organ geometry consisted of a target, two organs at risk (OARs), and normal tissue. A seven-equispaced-coplanar beam arrangement was used. For conventional IMRT, only 6 MV beams were available, while XMRT was optimized using 6 and 18 MV beams. A prescribed dose (PD) of 72 GY was assigned to the target, with upper and lower bounds of 110% and 95% of the PD, respectively. Both OARs were assigned a maximum dose of 64 Gy, while the normal tissue was assigned a maximum dose of 66 Gy. A numerical solver, Gurobi, generated solutions for the XMRT and IMRT problems. The dose-volume histograms from IMRT and XMRT solutions were compared. Results: The maximum, minimum, mean, and homogeneity of the dose to the target were comparable between IMRT and XMRT. Though IMRT had improved dose conformity relative to XMRT, XMRT reduced the mean dose to both OARs by more than 1 Gy. For normal tissue, an increase of 5 Gy in mean dose and 27 percent in integral dose was seen for IMRT relative to XMRT. Conclusion: This work demonstrates the benefits of simultaneously modulating photon beam energy and fluence using our XMRT approach in a given phantom geometry. While target coverage was comparable, dose to healthy structures was reduced using XMRT.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2015

A trust region algorithm for PDE-constrained optimization with bound constraints using reduced order modeling

This paper presents a novel trust region algorithm that relies on proper orthogonal de-composition techniques to construct accurate reduced order models during optimization.The algorithm samples high fidelity snapshots to compute the POD functions that are used to generate reduced order models. The reduced order models are used to replace the computationally intensive high fidelity finite element evaluations during optimization.The proposed algorithm employs a trust region framework to detect loss in predictive ac-curacy in the reduced order model and automatically update the POD functions during optimization. The trust region framework allows the algorithm to use sound mathematical metrics to effectively improve the accuracy and robustness of the reduced order models during optimization. The algorithm also employs a projected gradient algorithm to model bound constraints and compute optimal and feasible controls.This paper also presents an accurate Hessian formulation for topology optimization problems. The proposed trust region framework relies on a quadratic model to update the control. This quadratic model needs reliable second order derivative information to predict the behavior of the objective function within a suitable trust region. If a nonlinear Hessian formulation is used, the computational effort increases due to additional finite element evaluations. The proposed linear Hessian formulation reduces the computational effort and enables the calculation of the second order derivative information without additional finite element model evaluations. Examples in topology optimization are presented to demonstrate the applicability of the proposed algorithm and linear Hessian formulation for large-scale PDE-constrained optimization.

bound constraints↗
DOE OSTI2020

Towards a self-consistent analysis of the anisotropic galaxy two- and three-point correlation functions on large scales: application to mock galaxy catalogues

We establish a practical method for the joint analysis of anisotropic galaxy two- and three-point correlation functions (2PCF and 3PCF, respectively) on the basis of the decomposition formalism of the 3PCF using tripolar spherical harmonics. We perform such an analysis with MultiDark-Patchy mock catalogues to demonstrate and understand the benefit of the anisotropic 3PCF. We focus on scales above $80\, h^{-1}\, {\rm Mpc}$, and use information from the shape and the baryon acoustic oscillation (BAO) signals of the 2PCF and 3PCF. We also apply density field reconstruction to increase the signal-to-noise ratio of BAO in the 2PCF measurement, but not in the 3PCF measurement. In particular, we study in detail the constraints on the angular diameter distance and the Hubble parameter. We build a model of the bispectrum or 3PCF that includes the non-linear damping of the BAO signal in redshift space. We carefully account for various uncertainties in our analysis including theoretical models of the 3PCF, window function corrections, biases in estimated parameters from the fiducial values, the number of mock realizations to estimate the covariance matrix, and bin size. The joint analysis of the 2PCF and 3PCF monopole and quadrupole components shows a $30{{\ \rm per\ cent}}$ and $20{{\ \rm per\ cent}}$ improvement in Hubble parameter constraints before and after reconstruction of the 2PCF measurements, respectively, compared to the 2PCF analysis alone. Furthermore, this study clearly shows that the anisotropic 3PCF increases cosmological information from galaxy surveys and encourages further development of the modelling of the 3PCF on smaller scales than we consider.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1990

Time-optimal aircraft pursuit-evasion with a weapon envelope constraint

The optimal pursuit-evasion problem between two aircraft, including nonlinear point-mass vehicle models and a realistic weapon envelope, is analyzed. Using a linear combination of flight time and the square of the vehicle acceleration as the performance index, a closed-form solution is obtained in nonlinear feedback form. Due to its modest computational requirements, this guidance law can be used for onboard real-time implementation.

Menon, P. K. A.↗
NASA NTRS1989

An algorithm for the solution of dynamic linear programs

The algorithm's objective is to efficiently solve Dynamic Linear Programs (DLP) by taking advantage of their special staircase structure. This algorithm constitutes a stepping stone to an improved algorithm for solving Dynamic Quadratic Programs, which, in turn, would make the nonlinear programming method of Successive Quadratic Programs more practical for solving trajectory optimization problems. The ultimate goal is to being trajectory optimization solution speeds into the realm of real-time control. The algorithm exploits the staircase nature of the large constraint matrix of the equality-constrained DLPs encountered when solving inequality-constrained DLPs by an active set approach. A numerically-stable, staircase QL factorization of the staircase constraint matrix is carried out starting from its last rows and columns. The resulting recursion is like the time-varying Riccati equation from multi-stage LQR theory. The resulting factorization increases the efficiency of all of the typical LP solution operations over that of a dense matrix LP code. At the same time numerical stability is ensured. The algorithm also takes advantage of dynamic programming ideas about the cost-to-go by relaxing active pseudo constraints in a backwards sweeping process. This further decreases the cost per update of the LP rank-1 updating procedure, although it may result in more changes of the active set that if pseudo constraints were relaxed in a non-stagewise fashion. The usual stability of closed-loop Linear/Quadratic optimally-controlled systems, if it carries over to strictly linear cost functions, implies that the saving due to reduced factor update effort may outweigh the cost of an increased number of updates. An aerospace example is presented in which a ground-to-ground rocket's distance is maximized. This example demonstrates the applicability of this class of algorithms to aerospace guidance. It also sheds light on the efficacy of the proposed pseudo constraint relaxation scheme.

Psiaki, Mark L.↗
DOE OSTI2015

Internal transport barrier triggered by non-linear lower hybrid wave deposition under condition of beam-driven toroidal rotation

By using gyro-Landau fluid transport model (GLF23), time-dependent integrated modeling is carried out using TRANSP to explore the dynamic process of internal transport barrier (ITB) formation in the neutral beam heating discharges. When the current profile is controlled by LHCD (lower hybrid current drive), with appropriate neutral beam injection, the nonlinear interplay between the transport determined gradients in the plasma temperature (T{sub i,e}) and toroidal velocity (V{sub ϕ}) and the E×B flow shear (including q-profile) produces transport bifurcations, generating spontaneously a stepwise growing ITB. In the discharge, the constraints imposed by the wave propagation condition causes interplay of the LH driven current distribution with the plasma configuration modification, which constitutes non-linearity in the LH wave deposition. The non-linear effects cause bifurcation in LHCD, generating two distinct quasi-stationary reversed magnetic shear configurations. The change of current profile during the transition period between the two quasi-stationary states results in increase of the E×B shearing flow arising from toroidal rotation. The turbulence transport suppression by sheared E×B flow during the ITB development is analysed, and the temporal evolution of some parameters characterized the plasma confinement is examined. Ample evidence shows that onset of the ITB development is correlated with the enhancement of E×B shearing rate caused by the bifurcation in LHCD. It is suggested that the ITB triggering is associated with the non-linear effects of the LH power deposition.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Effective data sampling strategies and boundary condition constraints of physics-informed neural networks for identifying material properties in solid mechanics

Material identification is critical for understanding the relationship between mechanical properties and the associated mechanical functions. However, material identification is a challenging task, especially when the characteristic of the material is highly nonlinear in nature, as is common in biological tissue. In this work, we identify unknown material properties in continuum solid mechanics via physics-informed neural networks (PINNs). To improve the accuracy and efficiency of PINNs, we develop efficient strategies to nonuniformly sample observational data. We also investigate different approaches to enforce Dirichlet-type boundary conditions (BCs) as soft or hard constraints. Finally, we apply the proposed methods to a diverse set of time-dependent and time-independent solid mechanic examples that span linear elastic and hyperelastic material space. The estimated material parameters achieve relative errors of less than 1%. As such, this work is relevant to diverse applications, including optimizing structural integrity and developing novel materials.

97 MATHEMATICS AND COMPUTING↗
NASA NTRSDate not supplied

Multifidelity Optimization with Transonic Flutter Constraints

This work considers static and dynamic aeroelastic optimization of a cantilevered platewing in transonic flow. Low- and high-fidelity aeroelastic predictions are fed into a standard trust region model management scheme in order to efficiently solve this expensive optimization problem with multifidelity methods. Differences in fidelity are entirely driven by different aerodynamic solvers: linear panel methods for the low-fidelity method, and inviscid CFD-based solvers for the high-fidelity response. The optimization problem utilizes shape, sizing, and trim variables to satisfy stress, trim, and flutter constraints, and is solved across a range of subsonic and transonic Mach numbers. The multifidelity method is able to provide a speed-up for all cases considered here, despite sizable inaccuracies in the low-fidelity response for transonic flows.

Bret K Stanford↗
NASA NTRS1984

The influence of leading-edge load alleviation on supersonic wing design

A theoretical and experimental program to assess the effect of leading-edge load constraints on wing design and performance was conducted. For a planform characterized by a highly swept leading edge on the inboard region, linear theory was used to design camber surfaces which produced minimum drag-due-to-lift at the design lift coefficient of 0.08 and a design Mach number of 2.4. In an effort to delay the formation of leading edge vortices which often occur on highly swept wings, two approaches were used in the design criteria to limit the loadings on the leading edge. One wing was constrained to have the normal Mach number less than one everywhere along the leading edge and the second wing was constrained to have a pressure coefficient of zero on the leading edge. Force tests were run on the two constrained wings, on a flat reference wing and on an optimized wing with no leading edge constraints. All wings had identical planforms and thicknesses and were tested over a range of Mach numbers from 1.8 to 2.8 and a range in angles of attack from -5 deg to 8 deg. A comparison of the experimental performance of these four models is shown. Correlations of these results with theoretical predictions and flow visualization photographs are also included.

Darden, C. M.↗
DOE OSTI2019

A portalino to the dark sector

“Portal” models that connect the Standard Model to a Dark Sector allow for a wide variety of scenarios beyond the simplest WIMP models. Kinetic mixing of gauge fields in particular has allowed a broad range of new ideas. However, the models that evade CMB constraints are often non-generic, with new mass scales and operators to split states and suppress indirect detection signals. Models with a “portalino”, a neutral fermion that marries a linear combination of a standard model neutrino and dark sector fermion and carries a conserved quantum number, can be simpler. This is especially interesting for interacting dark sectors; then the unmarried linear combination which we identify as the standard model neutrino inherits these interactions too, and provides a new, effective interaction between the dark sector and the standard model. These interactions can be simple Z' type interactions or lepton-flavor changing. Dark matter freezes out into neutrinos, thereby evading CMB constraints, and conventional direct detection signals are largely absent. The model offers different signals, however. The “portalino” mechanism itself predicts small corrections to the standard model neutrino couplings as well as the possibility of discovering the portalino particle in collider experiments. Possible cosmological and astroparticle signatures include monochromatic neutrino signals from annihilation, spectral features in high energy CR neutrinos as well as conventional signals of additional light species and dark matter interactions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2016

MO-FG-BRC-02: Low-Z Switching Linear Accelerator Targets: New Options for Image Guidance and Dose Enhancement in Radiotherapy

Experimental research in medical physics has expanded the limits of our knowledge and provided novel imaging and therapy technologies for patients around the world. However, experimental efforts are challenging due to constraints in funding, space, time and other forms of institutional support. In this joint ESTRO-AAPM symposium, four exciting experimental projects from four different countries are highlighted. Each project is focused on a different aspect of radiation therapy. From the USA, we will hear about a new linear accelerator concept for more compact and efficient therapy devices. From Canada, we will learn about novel linear accelerator target design and the implications for imaging and therapy. From France, we will discover a mature translational effort to incorporate theranostic nanoparticles in MR-guided radiation therapy. From Germany, we will find out about a novel in-treatment imaging modality for particle therapy. These examples of high impact, experimental medical physics research are representative of the diversity of such efforts that are on-going around the globe. J. Robar, Research is supported through collaboration with Varian Medical Systems and Brainlab AGD. Westerly, This work is supported by the Department of Radiation Oncology at the University of Colorado School of Medicine. COI: NONEK. Parodi, Part of the presented work is supported by the DFG (German Research Foundation) Cluster of Excellence MAP (Munich-Centre for Advanced Photonics) and has been carried out in collaboration with IBA.

60 APPLIED LIFE SCIENCES↗