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 523 records · Page 29
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Simultaneous shape and topology optimization of inflatable soft robots

Simultaneous shape and topology optimization is used to design pressure-activated inflatable soft robots. The pressure loaded boundary is meshed conformingly and shape optimized, while the morphology of the robot is topology optimized. The design objective is to exert maximum force on an object, i.e. to produce soft “grippers”. The robot’s motion is modeled using nearly incompressible finite deformation hyperelasticity. To ensure stability of the robot, the buckling load factors obtained via linearized buckling analyses are constrained. The finite element method is used to evaluate the optimization cost and constraint functions and the adjoint method is employed to compute their sensitivities. The numerical examples produce pressure-driven soft robots with varying complexity. We also compare our simultaneous optimization results to those obtained via sequential topology and then shape optimization.

42 ENGINEERING↗
DOE OSTI2023

Improving initialization and evolution accuracy of cosmological neutrino simulations

Neutrino mass constraints are a primary focus of current and future large-scale structure (LSS) surveys. Non-linear LSS models rely heavily on cosmological simulations — the impact of massive neutrinos should therefore be included in these simulations in a realistic, computationally tractable, and controlled manner. Here, a recent proposal to reduce the related computational cost employs a symmetric neutrino momentum sampling strategy in the initial conditions. We implement a modified version of this strategy into the Hardware/Hybrid Accelerated Cosmology Code (HACC) and perform convergence tests on its internal parameters. We illustrate that this method can impart Script $\mathcal{O}$(1%) numerical artifacts on the total matter field on small scales, similar to previous findings, and present a method to remove these artifacts using Fourier-space filtering of the neutrino density field. Moreover, we show that the converged neutrino power spectrum does not follow linear theory predictions on relatively large scales at early times at the 15% level, prompting a more careful study of systematics in particle-based neutrino simulations. We also present an improved method for backscaling linear transfer functions for initial conditions in massive neutrino cosmologies that is based on achieving the same relative neutrino growth as computed with Boltzmann solvers. Our self-consistent backscaling method yields sub-percent accuracy in the total matter growth function. Comparisons for the non-linear power spectrum with the Mira-Titan emulator at a neutrino mass of m ν = 0.15 eV are in very good agreement with the expected level of errors in the emulator and in the direct N-body simulation.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Towards Improving luminosity using optics tuning and data-driven methods

The results of Run 24 experiments at Relativistic Heavy Ion Collider (RHIC) for improving luminosity using optics tuning are presented in this study. In the first experiment, MADx matching was used to output magnet strengths corresponding to specific s star movements around Interaction Region 8 (IR8). The corresponding Zero Degree Calorimeter (ZDC) signal was measured in place of luminosity, and Bayesian Optimization aids search of optimal movements. It was found that values retrieved from matching were inaccurate, resulting in negative feedback loops. The second experiment focused on calculating accurate s star movements. The matching method was replaced with a linear sensitivity matrix, directly relating optics to power supply, and its null space was used to fit constraints such as hysteresis effects. At the experiment, beam losses were observed at collimators around boundary of IR8, which were fixed for the third experiment. Dynamic mode decomposition was also introduced to improve quality of turn-by-turn (TBT) data as well as accuracy and consistency of optics measurements at IR8. These improvements will be tested in the experiment of next RHIC run for luminosity optimization.

Accelerator Physics↗
NASA NTRS1972

A VTOL translational rate control system study on a six degrees-of-freedom motion simulator

A linearized translational rate system for near hover flight was optimized on a large motion simulator under the constraints of no disturbances and limited control power. Both lateral and longitudinal modes were considered with the primary variables of investigation being control sensitivity and response stiffness and secondarily system damping. Yaw and height control characteristics were represented by an angular rate and acceleration system, respectively. General regions of desired sensitivity and stiffness for the longitudinal and lateral modes were determined under VFR conditions for both the rapid maneuver task and the station keeping/mild maneuver task.

Corliss, L. D.↗
NASA NTRS1982

Control law design to meet constraints using SYNPAC-synthesis package for active controls

Major features of SYNPAC (Synthesis Package for Active Controls) are described. SYNPAC employs constrained optimization techniques which allow explicit inclusion of design criteria (constraints) in the control law design process. Interrelationships are indicated between this constrained optimization approach, classical and linear quadratic Gaussian design techniques. Results are presented that were obtained by applying SYNPAC to the design of a combined stability augmentation/gust load alleviation control law for the DAST ARW-2.

Adams, W. M., Jr.↗
NASA NTRS1981

Entry trajectory shaping for Shuttle-deployed experiments

The availability of the Space Transportation System has provided a unique opportunity for deployment of atmospheric entry experiments. These experiments require a trajectory design that provides a well-defined and controlled flight environment. Items addressed in this paper include control of thermal and aerodynamic loading consistent with the test objectives, maintenance of acceptable velocity-altitude conditions, and accommodation of tracking and recovery requirements. Various strategies, including constant-drag flight, equilibrium glide, and controlled linear velocity-altitude flight segments are employed to effect desired trajectory conditions. Algorithms to match thermal constraints to entry conditions and vehicle aerodynamic configuration are also discussed.

Letts, W. R., Jr.↗
NASA NTRS1983

Extrapolation of optimum design based on sensitivity derivatives

In many applications, it is useful to have information regarding the effect of the perturbations of the design problem constants (parameters) on the optimum solution of the problem. Sobieszczanski-Sobieski et al. (1981) have, therefore, proposed to utilize an optimum sensitivity analysis as a structural design tool to assess the parameter perturbation effects by extrapolations based on optimum sensitivity derivatives. The present investigation is concerned with an example of optimum sensitivity analysis involving a composite panel. The case was selected in connection with the nonlinearity of the constraints and, the tendency for a set of active constraints to change when the parameters are disturbed. Under these adverse conditions, sensitivity derivatives are used in linear extrapolations to extend optimization information over a limited range of variation for a typical parameter.

Barthelemy, J.-F. M.↗
NASA NTRS1989

Use of the quasilinearization algorithm for the simulation of LSS slewing

The use of the Maximum Principle for the large angle slewing of large space structures (LSS) usually results in the so-called two-point boundary-value problem, in which many requirements (e.g., minimum time, small amplitude, and limited control power, etc.) must be satisfied simultaneously. The successful solution of this problem depends largely on the use of an efficient numerical algorithm. There are many candidate algorithms available for this problem (e.g., quasilinearization, gradient, etc.). Here researchers discuss only the quasilinearization method which has been used for several cases of large angle slewing of LSS. The basic idea of this algorithm is to make a series of successive approximations of the solution from a particular solvable case (linear or nonlinear) to a more general practical case. For the rigid spacecraft slewing problem with no constraints on the controls, the solution procedure can be found in the literature. This procedure needs to be modified if a minimum time for the slewing problem is desired with control limits given. Recently, an indirect method for finding the minimum time was developed to meet all these requirements. For the general mixed (including both rigid and flexible parts) problem, an additional constraint of small vibrational amplitude on the flexible parts is imposed.

Li, Feiyue↗
NASA NTRS1990

Inverse methods for assessing ship-of-opportunity networks and estimating circulation and winds from tropical expendable bathythermograph data

Inverse methods for estimating the surface ciculation of the equatorial Pacific by combining a linear reduced-gravity shallow-water model with the Tropical Ocean-Global Atmosphere ship-of-opportunity expendable bathythermograph (TOGA SOP XBT) observing program are examined. It is demonstrated that a simple linear model of the upper circulation of the equatorial Pacific can be successfully used as a weak constraint when smoothing the TOGA SOP XBT data. A circulation is sought as the weighted least squares fit to the dynamics and the data. The solution method is an expansion in representer functions, and the generalized inverse problem is thereby reduced from a functional problem to an algebraic problem for the coefficients of the representer. A specific inverse calculation using synthetic forcing and data is presented.

Bennett, Andrew F.↗
NASA NTRS2005

Aeroelastic Sizing for High-Speed Research (HSR) Longitudinal Control Alternatives Project (LCAP)

The Longitudinal Control Alternatives Project (LCAP) compared three high-speed civil transport configurations to determine potential advantages of the three associated longitudinal control concepts. The three aircraft configurations included a conventional configuration with a layout having a horizontal aft tail, a configuration with a forward canard in addition to a horizontal aft tail, and a configuration with only a forward canard. The three configurations were aeroelastically sized and were compared on the basis of operational empty weight (OEW) and longitudinal control characteristics. The sized structure consisted of composite honeycomb sandwich panels on both the wing and the fuselage. Design variables were the core depth of the sandwich and the thicknesses of the composite material which made up the face sheets of the sandwich. Each configuration was sized for minimum structural weight under linear and nonlinear aeroelastic loads subject to strain, buckling, ply-mixture, and subsonic and supersonic flutter constraints. This report describes the methods that were used and the results that were generated for the aeroelastic sizing of the three configurations.

Walsh, Joanne L.↗
NASA NTRS2008

5-kWe Free-piston Stirling Engine Convertor

The high reliability, long life, and efficient operation of Free-Piston Stirling Engines (FPSEs) make them an attractive power system to meet future space power requirements with less mass, better efficiency, and less total heat exchanger area than other power convertor options. FPSEs are also flexible in configuration as they can be coupled with many potential heat sources and various heat input systems, heat rejection systems, and power management and distribution systems. Development of a 5-kWe Stirling Convertor Assembly (SCA) is underway to demonstrate the viability of an FPSE for space power. The design is a scaled-down version of the successful 12.5-kWe Component Test Power Converter (CTPC) developed under NAS3-25463. The ultimate efficiency target is 25% overall convertor efficiency (electrical power out over heat in). For the single cylinder prototype now in development, cost and time constraints required use of economical and readily available materials (steel versus beryllium) and components (a commercially available linear alternator) and thus lower efficiency. The working gas is helium at 150 bar mean pressure. The design consists of a displacer suspended on internally pumped gas bearings and a power piston/alternator supported on flexures. Non-contacting clearance seals are used between internal volumes. Heat to and from the prototype convertor is done via pumped liquid loops passing through shell and tube heat exchangers. The preliminary and detail designs of the convertor, controller, and support systems (heating loop, cooling loop, and helium supply system) are complete and all hardware is on order. Assembly and test of the prototype at Foster- Miller is planned for early 2008, when work will focus on characterizing convertor dynamics and steady-state operation to determine maximum power output and system efficiency. The device will then be delivered to Auburn University where assessments will include start-up and shutdown characterization and transient response to temperature and load variations. Future activities may include testing at NASA GRC.

Chapman, Peter A.↗
NASA NTRS2009

Alternative Representations of Convective Processes in the NASA GEOS-5 AGCM

The gap in explicit resolution of phenomena between global climate models and cloud resolving models is shrinking at a steady pace with global integrations of several km in resolution now practical for at least some time scales. In moving toward finer resolution the long standing problem of convective parameterization is being examined along with the assumption of convective quasi-equilibrium and how this can be reconciled with the stochastic and intermittent nature of convection. In this context we examine the nature of parameterized convection in the NASA Goddard Earth Observing System (GEOS-5) Atmospheric General Circulation Model. Our analysis uses both coarse (2.5 degree) and fine scale 0.25 degree spatial resolution integrations. Two basic formulations are compared: the default option is the Relaxed Arakawa-Schubert (RAS) scheme which invokes a sequence of linearly entraining plumes and quasi-equilibrium closure. An optional modification of this method (the "Stochastic Tokioka" constraint) places a random lower limit on plume entrainment. An alternative representation is the Kain-Fritsch parameterization which was originally developed for mesoscale numerical modeling strategies. Here entrainment is determined by a crude buoyancy sorting approach that allows the plume spectrum to be more responsive to ambient vertical stratification of moisture. Diagnostics of the model behavior are referenced to recent observational evidence of continuous phase transition behavior. In particular we examine the relationship between column water vapor and probablility of convective presence and intensity. Sensitivity of the statistics of convective behavior to parcel mixing/entrainment formulations and parcel initial thermodynamics are considered. Observational statistics from A-Train and TRMM sensors provide validation of the model integrations.

Robertson, Franklin↗
NASA NTRS2024

Global Optimization via Quadratic Disjunctive Programming for Water Networks Design with Energy Recovery

Generalized disjunctive programming (GDP) models with bilinear and concave constraints, often seen in water network design, are challenging optimization problems. This work proposes quadratic and piecewise linear approximations for nonlinear terms to reformulate GDP models into quadratic GDP (QGDP) models that suitable solvers may solve more efficiently. We illustrate the benefits of the quadratic reformulation with a water treatment network design problem in which nonconvexities arise from bilinear terms in the mixers’ mass balances and concave investment cost functions of treatment units. Given the similarities with water network design problems, we suggest quadratic approximation for the GDP model for the optimal design of a large-scale reverse electrodialysis (RED) process. This power technology can recover energy from salinity differences between by-product streams of the water sector, such as desalination brine mixed with regenerated wastewater effluents. The solver Gurobi excels in handling QGDP problems, but weighing the problem’s precision and tractability balance is crucial. The piecewise linear approximation yields more accurate, yet larger QGDP models that may require longer optimization times in large-scale process synthesis problems.

Water Networks↗
DOE OSTI2020

Examining the Net Revenue and Downstream Flow Impact Trade-Offs for a Network of Cascading, Small-Scale Hydropower Facilities: Preprint

In this work, we used a price-taker model to investigate the trade-offs between net revenue and downstream flow impacts for a network of small, cascading hydro facilities. The network consisted of 36 facilities, each with a small amount of local storage (between 2 and 45 minutes). Generator sizes ranged between 0.5 and 1 MW, nominal, and the total capacity of the network was 33.5 MW. We used a multi-integer linear programing model to maximize the net revenue of the combined network subject to operating and environmental constraints. Net revenue optimizations relied on historic price data, and dry, typical, and wet years were studied to help ensure robustness. Energy and ancillary service sales were included in the net revenue calculations, and both unit commitment and dispatch simulations were performed. We found that limiting the downstream flows to ±50% of the river’s natural flows had a negligible impact on net revenues (<1% reduction), irrespective of hydrologic conditions, and even when downstream flows were limited to ±5%, net revenues were only impacted by 4%. These outcomes are significant because they demonstrate how an array of small-scale hydropower facilities can be operated to have minimal impact on natural stream flows—addressing a critical environmental concern.

downstream flow↗
NASA NTRS1999

Engineering Overview of a Multidisciplinary HSCT Design Framework Using Medium-Fidelity Analysis Codes

An objective of the HPCC Program at NASA Langley has been to promote the use of advanced computing techniques to more rapidly solve the problem of multidisciplinary optimization of a supersonic transport configuration. As a result, a software system has been designed and is being implemented to integrate a set of existing discipline analysis codes, some of them CPU-intensive, into a distributed computational framework for the design of a High Speed Civil Transport (HSCT) configuration. The proposed paper will describe the engineering aspects of integrating these analysis codes and additional interface codes into an automated design system. The objective of the design problem is to optimize the aircraft weight for given mission conditions, range, and payload requirements, subject to aerodynamic, structural, and performance constraints. The design variables include both thicknesses of structural elements and geometric parameters that define the external aircraft shape. An optimization model has been adopted that uses the multidisciplinary analysis results and the derivatives of the solution with respect to the design variables to formulate a linearized model that provides input to the CONMIN optimization code, which outputs new values for the design variables. The analysis process begins by deriving the updated geometries and grids from the baseline geometries and grids using the new values for the design variables. This free-form deformation approach provides internal FEM (finite element method) grids that are consistent with aerodynamic surface grids. The next step involves using the derived FEM and section properties in a weights process to calculate detailed weights and the center of gravity location for specified flight conditions. The weights process computes the as-built weight, weight distribution, and weight sensitivities for given aircraft configurations at various mass cases. Currently, two mass cases are considered: cruise and gross take-off weight (GTOW). Weights information is obtained from correlations of data from three sources: 1) as-built initial structural and non-structural weights from an existing database, 2) theoretical FEM structural weights and sensitivities from Genesis, and 3) empirical as-built weight increments, non-structural weights, and weight sensitivities from FLOPS. For the aeroelastic analysis, a variable-fidelity aerodynamic analysis has been adopted. This approach uses infrequent CPU-intensive non-linear CFD to calculate a non-linear correction relative to a linear aero calculation for the same aerodynamic surface at an angle of attack that results in the same configuration lift. For efficiency, this nonlinear correction is applied after each subsequent linear aero solution during the iterations between the aerodynamic and structural analyses. Convergence is achieved when the vehicle shape being used for the aerodynamic calculations is consistent with the structural deformations caused by the aerodynamic loads. To make the structural analyses more efficient, a linearized structural deformation model has been adopted, in which a single stiffness matrix can be used to solve for the deformations under all the load conditions. Using the converged aerodynamic loads, a final set of structural analyses are performed to determine the stress distributions and the buckling conditions for constraint calculation. Performance constraints are obtained by running FLOPS using drag polars that are computed using results from non-linear corrections to the linear aero code plus several codes to provide drag increments due to skin friction, wave drag, and other miscellaneous drag contributions. The status of the integration effort will be presented in the proposed paper, and results will be provided that illustrate the degree of accuracy in the linearizations that have been employed.

Weston, R. P.↗
DOE OSTI2020

Time-dependent relaxed magnetohydrodynamics: Inclusion of cross helicity constraint using phase-space action

In this work, a phase-space version of the ideal magnetohydrodynamic (MHD) Lagrangian is derived from first principles and shown to give a relabeling transformation when a cross-helicity constraint is added in Hamilton's Action Principle. A new formulation of time-dependent “relaxed” magnetohydrodynamics is derived using microscopic conservation of mass and macroscopic constraints on total magnetic helicity, cross helicity, and entropy under variations of density, pressure, fluid velocity, and magnetic vector potential. This gives Euler–Lagrange equations consistent with previous work on both ideal and relaxed MHD equilibria with flow, but generalizes the relaxation concept from statics to dynamics. The application of the new dynamical formalism is illustrated for short-wavelength linear waves, and the interface connection conditions for Multiregion Relaxed MHD (MRxMHD) are derived. The issue of whether E+u×B=0 should be a constraint is discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Fuzzy dark matter and the Dark Energy Survey Year 1 data

ABSTRACT Gravitational weak lensing by dark matter haloes leads to a measurable imprint in the shear correlation function of galaxies. Fuzzy dark matter (FDM), composed of ultralight axion-like particles of mass m ∼ 10−22 eV, suppresses the matter power spectrum and shear correlation with respect to standard cold dark matter. We model the effect of FDM on cosmic shear using the optimized halo model HMCode, accounting for additional suppression of the mass function and halo concentration in FDM as observed in N-body simulations. We combine Dark Energy Survey Year 1 (DES-Y1) data with the Planck cosmic microwave background anisotropies to search for shear correlation suppression caused by FDM. We find no evidence of suppression compared to the preferred cold dark matter model, and thus set a new lower limit to the FDM particle mass. Using a log-flat prior and marginalizing over uncertainties related to the non-linear model of FDM, we find a new, independent 95 per cent C.L. lower limit log10m > −23 combining Planck and DES-Y1 shear, an improvement of almost two orders of magnitude on the mass bound relative to CMB-only constraints. Our analysis is largely independent of baryonic modelling, and of previous limits to FDM covering this mass range. Our analysis highlights the most important aspects of the FDM non-linear model for future investigation. The limit to FDM from weak lensing could be improved by up to three orders of magnitude with $\mathcal {O}(0.1)$ arcmin cosmic shear angular resolution, if FDM and baryonic feedback can be simultaneously modelled to high precision in the halo model.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2005

Small Body GN&C Research Report: A Guidance and Control Technique for Small-Body Proximity Operations with Guaranteed Guidance Resolvability and Required Thruster Silent Time

The guidance and control (G&C) algorithms for enabling small-body proximity operations are developed by using a model predictive control approach along with a convexification of the governing dynamics, control constraints, and trajectory/state constraints. The open-loop guidance is solved ahead of time or in a resolvable, real-time manner through the use of PWG (Pseudo Way-point Generation), a technique developed in this research. The PWG scheme ensures required thruster silent times during trajectory maneuvers. The feedback control is implemented to track the PWG trajectories in a manner that guarantees the resolvability for the open-loop problem, enabling the ability to update the G&C in a model-predictive manner. The schemes incorporate gravity models and thruster ring times into discrete dynamics that are solved as a optimal control problem to minimize fuel consumption or thruster energy expenditure. The optimal control problem is cast as an LMI (Linear Matrix Inequality) and then solved through Semi-Definite Programming techniques in a computationally efficient manner that provides convergence and constraint guarantees.

guidance and control (G&C)↗