Search NASA⌕ Search

SEARCH · Search NASA

Results for “(non)linear optimization”

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

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

At least 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
NASA NTRS1980

The linear and non-linear aerodynamics of three-surface aircraft concepts

It is noted that most modern fighter aircraft rely on vortex interaction to provide lift enhancement at maneuvering angles of attack. It is shown that the close-coupled horizontal canard in a three-surface configuration provides a control surface which in addition to its other control functions, can be used to optimize this vortex interaction. Attention is given to a study intended to provide a detailed understanding of the aerodynamics of this type of configuration. The discussion examines the results of this investigation and hypotheses are presented to explain the linear and nonlinear aerodynamic phenomena observed.

Agnew, J. W.↗
DOE OSTI2022

SNICAR-ADv4: a physically based radiative transfer model to represent the spectral albedo of glacier ice

Abstract. Accurate modeling of cryospheric surface albedo is essential for our understanding of climate change as snow and ice surfaces regulate the global radiative budget and sea-level through their albedo and mass balance. Although significant progress has been made using physical principles to represent the dynamic albedo of snow, models of glacier ice albedo tend to be heavily parameterized and not explicitly connected with physical properties that govern albedo, such as the number and size of air bubbles, specific surface area (SSA), presence of abiotic and biotic light absorbing constituents (LACs), and characteristics of any overlying snow. Here, we introduce SNICAR-ADv4, an extension of the multi-layer two-stream delta-Eddington radiative transfer model with the adding–doubling solver that has been previously applied to represent snow and sea-ice spectral albedo. SNICAR-ADv4 treats spectrally resolved Fresnel reflectance and transmittance between overlying snow and higher-density glacier ice, scattering by air bubbles of varying sizes, and numerous types of LACs. SNICAR-ADv4 simulates a wide range of clean snow and ice broadband albedo (BBA), ranging from 0.88 for (30 µm) fine-grain snow to 0.03 for bare and bubble-free ice under direct light. Our results indicate that representing ice with a density of 650 kg m−3 as snow with no refractive Fresnel layer, as done previously, generally overestimates the BBA by an average of 0.058. However, because most naturally occurring ice surfaces are roughened “white ice”, we recommend modeling a thin snow layer over bare ice simulations. We find optimal agreement with measurements by representing cryospheric media with densities less than 650 kg m−3 as snow and larger-density media as bubbly ice with a Fresnel layer. SNICAR-ADv4 also simulates the non-linear albedo impacts from LACs with changing ice SSA, with peak impact per unit mass of LACs near SSAs of 0.1–0.01 m2 kg−1. For bare, bubble-free ice, LACs actually increase the albedo. SNICAR-ADv4 represents smooth transitions between snow, firn, and ice surfaces and accurately reproduces measured spectral albedos of a variety of glacier surfaces. This work paves the way for adapting SNICAR-ADv4 to be used in land ice model components of Earth system models.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1992

Generating Far-Infrared Radiation By Two-Wave Mixing

Far-infrared radiation 1 to 6 GHz generated by two-wave mixing in asymmetrically grown GaAs/AlxGa1-xAs multiple-quantum-well devices. Two near-infrared semiconductor diode lasers phase-locked. Outputs amplified, then combined in semiconductor nonlinear multiple-quantum-well planar waveguide. Necessary to optimize design of device with respect to three factors: high degree of confinement of electromagnetic field in nonlinear medium to maximize power density, phase matching to extend length of zone of interaction between laser beams in non-linear medium, and nonlinear susceptibility. Devices used as tunable local oscillators in heterodyne-detection radiometers.

Borenstain, Shmuel↗
NASA NTRS1993

Controls-structures-interaction dynamics during RCS control of the Orbiter/SRMS/SSF configuration

During the assembly flights of the Space Station Freedom (SSF), the Orbiter will either dock with the SSF and retract to the final berthed position, or will grapple the SSF using the Shuttle Remote Manipulator System (SRMS) and maneuver the SRMS coupled vehicles to their final berthed position. The SRMS method is expected to take approximately one to one and a half hours to complete and require periodic attitude corrections by either the Orbiter or the SSF reaction control system (RCS) or continuous control by a control moment gyro (CMG) system with RCS desaturation as required. Free drift of the attached vehicles is not currently thought to be acceptable because the desired system attitude will quickly deteriorate due to unbalanced gravity gradient and aerodynamic torques resulting in power generation problems, thermodynamic control problems, and communications problems. This paper deals with the simulation and control of the SRMS during trunnion/latch interaction dynamics and during RCS maneuvers. The SRMS servo drive joints have highly non-linear elastic characteristics which tend to degrade sensitive control strategies. In addition the system natural frequencies are extremely low and depend on the drive joint deflections and SRMS geometric position. The lowest mean period of oscillation for the Orbiter/SRMS/SSF(MB6) system in brakes hold mode positioned near the final berthed position is approximately 120 seconds. A detailed finite element model of the SRMS has been developed and used in a newly developed SRMS systems dynamics simulation to investigate the non-linear transient response dynamics of the Orbiter/SRMS/SSF systems. The present SRMS control strategy of brakes only recommended by the Charles Draper Labs is contrasted with a robust controller developed by the authors. The robust controller uses an optimal inear quadratic regulator (LQR) to optimally place the closed-loop poles of a multivariable continuous-time system within the common region of an open sector with the sector angle plus or minus 45 degrees from the negative real axis, and the left-hand side of a parallel to the imaginary axis in the complex s-plane. This guarantees that the critical damping ratio for the desired control modes is equal to or in excess of 0.707. The matrix sign function is used for solving the Riccati equations which appear in the controller design procedure. Fast and stable algorithms have recently been developed for the computation of the matrix sign function. Simulation results are given which demonstrate the potential CSI involvement for the current SRMS control system and the proposed control system.

Schliesing, J. A.↗
DOE OSTI2021

Physics-guided logistic classification for tool life modeling and process parameter optimization in machining

This paper describes a physics-guided logistic classification method for tool life modeling and process parameter optimization in machining. Tool life is modeled using a classification method since the exact tool life cannot be measured in a typical production environment where tool wear can only be directly measured when the tool is replaced. Here, in this study, laboratory tool wear experiments are used to simulate tool wear data normally collected during part production. Two states are defined: tool not worn (class 0) and tool worn (class 1). The non-linear reduction in tool life with cutting speed is modeled by applying a logarithmic transformation to the inputs for the logistic classification model. A method for interpretability of the logistic model coefficients is provided by comparison with the empirical Taylor tool life model. The method is validated using tool wear experiments for milling. Results show that the physics-guided logistic classification method can predict tool life using limited datasets. A method for pre-process optimization of machining parameters using a probabilistic machining cost model is presented. The proposed method offers a robust and practical approach to tool life modeling and process parameter optimization in a production environment.

Machine learning↗
DOE OSTI2025

Uniform secondary current distribution at disk electrodes under Tafel kinetics enabled by concentric current-shielding rings

Non-uniform secondary current distribution at rotating disk electrodes (RDE) is a common problem when using resistive electrolyte media or large applied currents. In a recent publication, we have shown that auxiliary electrodes such as the ring of a rotating ring disk electrode (RRDE) can help suppress current non-uniformities at the disk enabling reliable electroanalytical measurements. However, this previous work considered linear kinetics where current distribution non-uniformities were moderate. In the present contribution, we consider current distribution non-uniformities encountered under Tafel kinetics. We show, for the case of 2,5-dihydroxy-1,4-benzoquinone (DHBQ) reduction, that optimally chosen ring conditions serve to provide effective shielding at the disk edge rendering the overall disk current distribution to be uniform. Numerical modeling and scaling analysis (using the Wagner number) are presented to aid a user in determining the optimal ring current density for achieving uniform disk current distribution under Tafel kinetics. This approach is especially useful when studying soluble-soluble redox transitions for which, unlike deposit distribution in electrodeposition, the current distribution non-uniformity is not visually apparent.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Amplitude- and frequency-dependent activation of layer II/III neurons by intracortical microstimulation

Intracortical microstimulation (ICMS) has been used for the development of brain machine interfaces. However, further understanding about the spatiotemporal responses of neurons to different electrical stimulation parameters is necessary to inform the design of optimal therapies. In this study, we employed in vivo electrophysiological recording, two-photon calcium imaging, and electric field simulation to evaluate the acute effect of ICMS on layer II/III neurons. Our results show that stimulation frequency non-linearly modulates neuronal responses, whereas the magnitude of responses is linearly correlated to the electric field strength and stimulation amplitude before reaching a steady state. Temporal dynamics of neurons’ responses depends more on stimulation frequency and their distance to the stimulation electrode. In addition, amplitude-dependent post-stimulation suppression was observed within ~500 μm of the stimulation electrode, as evidenced by both calcium imaging and local field potentials. These findings provide insights for selecting stimulation parameters to achieve desirable spatiotemporal specificity of ICMS.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Structure optimization with stochastic density functional theory

Linear-scaling techniques for Kohn–Sham density functional theory are essential to describe the ground state properties of extended systems. Still, these techniques often rely on the localization of the density matrix or accurate embedding approaches, limiting their applicability. In contrast, stochastic density functional theory (sDFT) achieves linear- and sub-linear scaling by statistically sampling the ground state density without relying on embedding or imposing localization. In return, ground state observables, such as the forces on the nuclei, fluctuate in sDFT, making optimizing the nuclear structure a highly non-trivial problem. In this work, we combine the most recent noise-reduction schemes for sDFT with stochastic optimization algorithms to perform structure optimization within sDFT. We compare the performance of the stochastic gradient descent approach and its variations (stochastic gradient descent with momentum) with stochastic optimization techniques that rely on the Hessian, such as the stochastic Broyden–Fletcher–Goldfarb–Shanno algorithm. In conclusion, we further provide a detailed assessment of the computational efficiency and its dependence on the optimization parameters of each method for determining the ground state structure of bulk silicon with varying supercell dimensions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Non-linear MHD modelling of edge localized modes suppression by resonant magnetic perturbations in ITER

Edge localized modes (ELMs) suppression by resonant magnetic perturbations (RMPs) was studied with the non-linear magneto-hydro-dynamic (MHD) code JOREK for the ITER H-mode scenarios at 15 MA, 12.5 MA, 10 MA/5.3 T. In this study, the main aim was to demonstrate that ELMs can be suppressed by RMPs while the divertor 3D footprints of heat and particle fluxes remain within divertor material limits. The unstable peeling–ballooning modes responsible for ELMs without RMPs were modelled first for each scenario using numerically accessible parameters for ITER. Then the stabilization of ELMs by RMPs was modelled with the same parameters. RMP spectra, optimized by the linear MHD MARS-F code, with main toroidal harmonics N = 2, N = 3, N = 4 have been used as boundary conditions of the computational domain of JOREK, including realistic RMP coils, main plasma, scrape off layer (SOL) divertor and realistic first wall. The model includes all relevant plasma flows: toroidal rotation, two fluid diamagnetic effects and neoclassical poloidal friction. With RMPs, the main toroidal harmonic and the non-linearly coupled harmonics remain dominant at the plasma edge, producing saturated modes and a continuous MHD turbulent transport thereby avoiding ELM crashes in all scenarios considered here. The threshold for ELM suppression was found at a maximum RMP coils current of 45 kAt–60 kAt compared to the coils maximum capability of 90 kAt. In the high beta poloidal steady-state 10 MA/5.3 T scenario, a rotating QH-mode without ELMs was observed even without RMPs. In this scenario with RMPs N = 3, N = 4 at 20 kAt maximum current in RMP coils, similar QH-mode behaviour was observed however with dominant edge harmonic corresponding to the main toroidal number of RMPs. The present MHD modelling was limited in time by few tens of ms after RMPs were switched on until the magnetic energy of the modes saturates. As a consequence the thermal energy was still evolving on this time scale, far from the ITER confinement time scale and hence only the form of 3D footprints on the divertor targets can be indicated within this set-up. Also note, that the divertor physics was missing in this model, so realistic values of fluxes are out of reach in this modelling. However the stationary 3D divertor and particle fluxes could be simply extrapolated from these results to the stationary situation considering that a large power fraction should be radiated in the core and SOL and only about 50 MW power is going to the divertor, which is an arbitrary, but reasonable number used here. The 3D footprints with RMPs show the characteristic splitting with the main RMP toroidal symmetry. The maximum radial extension of the footprints typically was ~20 cm in inner divertor and ~40 cm in outer divertor with stationary heat fluxes decreasing further out from the initial strike point from ~5 MW m ₋2 to ~1 MW m ₋2 assuming a total power in the divertor and walls is 50 MW. The heat fluxes remain within the divertor target and baffle areas, however with rather small margin in the outer divertor which could be an issue for the first wall especially in transient regimes when part of the plasma thermal energy is released due to switching on the RMP coils. This fact should be considered when RMPs are applied with a more favorable application before or soon after the L–H transition, although optimization is required to avoid increasing the L–H power threshold with RMPs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1972

On the application of deterministic optimization methods to stochastic control problems.

A technique is presented by which one can apply the Minimum Principle of Pontryagin to stochastic optimal control problems formulated around linear systems with Gaussian noises and general cost criteria. Using this technique, the stochastic nature of the problem is suppressed but for two expectation operations, the optimization being essentially deterministic. The technique is applied to systems with quadratic and non-quadratic costs to illustrate its use.

Kramer, L. C.↗
NASA NTRS1991

Neural networks: What non-linearity to choose

Neural networks are now one of the most successful learning formalisms. Neurons transform inputs (x(sub 1),...,x(sub n)) into an output f(w(sub 1)x(sub 1) + ... + w(sub n)x(sub n)), where f is a non-linear function and w, are adjustable weights. What f to choose? Usually the logistic function is chosen, but sometimes the use of different functions improves the practical efficiency of the network. The problem of choosing f as a mathematical optimization problem is formulated and solved under different optimality criteria. As a result, a list of functions f that are optimal under these criteria are determined. This list includes both the functions that were empirically proved to be the best for some problems, and some new functions that may be worth trying.

Kreinovich, Vladik YA.↗
NASA NTRSDate not supplied

Pterodactyl: Guidance and Control of a Symmetric Deployable Entry Vehicle using an Aerodynamic Control System

The NASA-funded Pterodactyl project seeks to advance the state-of-the-art for varying entry vehicle types by developing unconventional guidance and control technologies for Deployable Entry Vehicles (DEVs) that can be applied to different entry vehicle configurations. Prior work by the authors [1–5] involved developing both traditional and novel integrated guidance and control solutions for a Pterodactyl Baseline Vehicle (PBV), a variant of an asymmetric DEV called the Lifting Nano ADEPT (LNA) [6]. In the prior studies, two different guidance schemes were designed and implemented for the PBV: (i) traditional bank angle guidance developed using the Fully Numerical Predictor-Corrector Entry Guidance (FNPEG) and (ii) novel angle of attack and sideslip (α - β) guidance developed using FNPEG with Uncoupled Range Control [4]. Using Linear Quadratic Regulator (LQR) optimal control methods with state-feedback integral control designs, these guidance trajectories were designed to be tracked using (i) a conventional propulsive entry vehicle control hardware architecture - reaction control systems (RCS) and (ii) novel non-propulsive entry vehicle control systems - aerodynamic flap control system (FCS) and moving mass control system (MMCS) [1]. The novel FCS and MMCS architectures were designed to track α - β guidance commands while the RCS was designed to track bank angle commands. It was discovered that the asymmetric DEV, the PBV, experienced a non-zero induced roll moment due to sideslip that the FCS and MMCS architectures had limited capability to trim out. These two architectures were designed to provide independent angle of attack and sideslip commands with limited consideration for roll moment generation to trim. As a result, for the PBV, the FCS and MMCS configurations as designed, were limited in providing the control authority needed to track an α - β guidance trajectory [1]. These results form the motivation for the work presented in this paper - utilizing an aerodynamic control system to track α - β guidance commands for a symmetric DEV with the expectation that a symmetric entry vehicle will have zero or significantly reduced roll moment due to sideslip that the FCS can handle when tracking an α - β guidance trajectory. To demonstrate the feasibility of a novel guidance and control architecture on a DEV, we utilize a symmetric DEV, the PBV-II, for (i) the novel α - β guidance development using FNPEG with Uncoupled Range Control and (ii) LQR control design using eight aerodynamic control surfaces. This paper demonstrates that the novel uncoupled α - β guidance tracking can be achieved using aerodynamic control surfaces on a symmetric deployable entry vehicle configuration.

Wendy A Okolo↗
NASA NTRS2020

Pterodactyl: Control Architectures Development for Integrated Control Design of a Mechanically Deployed Entry Vehicle

The need to return high mass payloads is driving the development of a new class of vehicles, Deployable Entry Vehicles (DEV) for which feasible and optimized control architectures have not been developed. The Pterodactyl project, seeks to advance the current state-of-the-art for entry vehicles by developing a design, test, and build capability for DEVs that can be applied to various entry vehicle configurations. This paper details the efforts on the NASA-funded Pterodactyl project to investigate multiple control techniques for the Lifting Nano-ADEPT (LNA) DEV. We design and implement multiple control architectures on the LNA and evaluate their performance in achieving varying guidance commands during entry.First we present an overview of DEVs and the Lifting Nano-ADEPT (LNA), along with the physical LNA configuration that influences the different control designs. Existing state-of-the-art for entry vehicle control is primarily propulsive as reaction control systems (RCS) are widely employed. In this work, we analyze the feasibility of using both propulsive control systems such as RCS to generate moments, and non-propulsive control systems such as aerodynamic control surfaces and internal moving mass actuations to shift the LNA center of gravity and generate moments. For these diverse control systems, we design different multi-input multi-output (MIMO) state-feedback integral controllers based on linear quadratic regulator (LQR) optimal control methods. The control variables calculated by the controllers vary, depending on the control system being utilized and the outputs to track for the controller are either the (i) bank angle or the (ii) angle of attack and sideslip angle as determined by the desired guidance trajectory. The LQR control design technique allows the relative allocation of the control variables through the choice of the weighting matrices in the cost index. Thus, it is easy to (i) specify which and how much of a control variable to use, and (ii) utilize one control design for different control architectures by simply modifying the choice of the weighting matrices.By providing a comparative analysis of multiple control systems, configurations, and performance, this paper and the Pterodactyl project as a whole will help entry vehicle system designers and control systems engineers determine suitable control architectures for integration with DEVs and other entry vehicle types.

Okolo, Wendy A.↗
DOE OSTI2025

Harnessing Machine Learning to Predict MoS 2 Solid Lubricant Performance

Physical vapor deposited (PVD) molybdenum disulfide (MoS 2 ) solid lubricant coatings are an exemplar material system for machine learning methods due to small changes in process variables often causing large variations in microstructure and mechanical/tribological properties. Here, in this work, a gradient boosted regression tree machine learning method is applied to an existing experimental data set containing process, microstructure, and property information to create deeper insights into the process-structure–property relationships for molybdenum disulfide (MoS 2 ) solid lubricant coatings. The optimized and cross-validated models show good predictive capabilities for density, reduced modulus, hardness, wear rate, and initial coefficients of friction. The contribution of individual deposition variables (i.e., argon pressure, deposition power, target conditioning) on coating properties is highlighted through feature importance. The process-property relationships established herein show linear and non-linear relationships and highlight the influence of uncontrolled deposition variables (i.e., target conditioning) on the tribological performance.

MoS2↗
DOE OSTI2021

Sequence of polyhedral relaxations for nonlinear univariate functions

Here, given a nonlinear, univariate, bounded, and differentiable function f(x), this article develops a sequence of Mixed Integer Linear Programming (MILP) and Linear Programming (LP) relaxations that converge to the graph of f(x) and its convex hull, respectively. Theoretical convergence of the sequence of relaxations to the graph of the function and its convex hull is established. For nonlinear non-convex optimization problems, the relaxations presented in this article can be used to construct tight MILP and LP relaxations. These MILP and the LP relaxations can also be used with MILP-based and spatial branch-and-bound based global optimization algorithms, respectively.

42 ENGINEERING↗
NASA NTRS1988

An integrated approach to the synthesis of geometrically non-linear structures

An integrated approach to the minimum weight design of geometrically nonlinear three-dimensional truss structures with geometric imperfections, subject to inequality constraints on static displacements, stresses, local buckling and cross sectional areas, is investigated. The integrated structural synthesis problem involves design and response quantities as independent variables and equilibrium equations, describing the finite element model, as equality constraints. The nonlinear structural analysis and the optimization are thus merged together into a single process. A computer program developed to compute the constraint values and analytical gradients is coupled with a generalized reduced gradient algorithm to solve the integrated problem. Numerical results for a geometrically nonlinear shallow dome example problem are presented for various types of imperfections. Furthermore, it is found that the algorithm is capable of detecting and guarding against system as well as element elastic instability using equilibrium information only, that is, without imposing system and local buckling inequality constraints.

Smaoui, H.↗
DOE OSTI2023

On unifying randomized methods for inverse problems

This work unifies the analysis of various randomized methods for solving linear and nonlinear inverse problems with Gaussian priors by framing the problem in a stochastic optimization setting. By doing so, we show that many randomized methods are variants of a sample average approximation (SAA). More importantly, we are able to prove a single theoretical result that guarantees the asymptotic convergence for a variety of randomized methods. Additionally, viewing randomized methods as an SAA enables us to prove, for the first time, a single non-asymptotic error result that holds for randomized methods under consideration. Another important consequence of our unified framework is that it allows us to discover new randomization methods. Here, we present various numerical results for linear, nonlinear, algebraic, and PDE-constrained inverse problems that verify the theoretical convergence results and provide a discussion on the apparently different convergence rates and the behavior for various randomized methods.

42 ENGINEERING↗