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.

SourceResearch connectionFrom yearThrough year

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

At least 379 records · Page 21

All sources NASA NTRS DOE OSTI

DOE OSTI2022

Data-driven Minimum Entropy Control for Stochastic Nonlinear Systems using the Cumulant-Generating Function

Here, we present a novel minimum entropy control algorithm for a class of stochastic nonlinear systems subjected to non-Gaussian noises. The entropy control can be considered as an optimization problem for the system randomness attenuation, but the mean value has to be considered separately. To overcome this disadvantage, a new representation of the system stochastic properties was given using the cumulant-generating function based on the moment-generating function, in which the mean value and the entropy was reflected by the shape of the cumulant-generating function. Based on the samples of the system output and control input, a time-variant linear model was identified, and the minimum entropy optimization was transformed to system stabilization. Then, an optimal control strategy was developed to achieve the randomness attenuation, and the boundedness of the controlled system output was analyzed. The effectiveness of the presented control algorithm was demonstrated by a numerical example. In this paper, a data-driven minimum entropy design is presented without pre-knowledge of the system model; entropy optimization is achieved by the system stabilization approach in which the stochastic distribution control and minimum entropy are unified using the same identified structure; and a potential framework is obtained since all the existing system stabilization methods can be adopted to achieve the minimum entropy objective.

42 ENGINEERING↗

DOE OSTI2024

Responses of Marginal and Intrinsic Water-Use Efficiency to Changing Aridity Using FLUXNET Observations

According to classic stomatal optimization theory, plant stomata are regulated to maximize carbon assimilation for a given water loss. A key component of stomatal optimization models is marginal water-use efficiency (mWUE), the ratio of the change of transpiration to the change in carbon assimilation. Although the mWUE is often assumed to be constant, variability of mWUE under changing hydrologic conditions has been reported. However, there has yet to be a consensus on the patterns of mWUE variabilities and their relations with atmospheric aridity. We investigate the dynamics of mWUE in response to vapor pressure deficit (VPD) and aridity index using carbon and water fluxes from 115 eddy covariance towers available from the global database FLUXNET. We demonstrate a non-linear mWUE-VPD relationship at a sub-daily scale in general; mWUE varies substantially at both low and high VPD levels. However, mWUE remains relatively constant within the mid-range of VPD. Despite the highly non-linear relationship between mWUE and VPD, the relationship can be informed by the strong linear relationship between ecosystem-level inherent water-use efficiency (IWUE) and mWUE using the slope, m *. We further identify site-specific m * and its variability with changing site-level aridity across six vegetation types. We suggest accurately representing the relationship between IWUE and VPD using Michaelis–Menten or quadratic functions to ensure precise estimation of mWUE variability for individual sites.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS2018

Design of a Modular Monolithic Implicit Solver for Multi-Physics Applications

The design of a modular multi-physics high-order space-time finite-element framework is presented together with its extension to allow monolithic coupling of different physics. One of the main objectives of the framework is to perform efficient high- fidelity simulations of capsule/parachute systems. This problem requires simulating multiple physics including, but not limited to, the compressible Navier-Stokes equations, the dynamics of a moving body with mesh deformations and adaptation, the linear shell equations, non-re effective boundary conditions and wall modeling. The solver is based on high-order space-time - finite element methods. Continuous, discontinuous and C1-discontinuous Galerkin methods are implemented, allowing one to discretize various physical models. Tangent and adjoint sensitivity analysis are also targeted in order to conduct gradient-based optimization, error estimation, mesh adaptation, and flow control, adding another layer of complexity to the framework. The decisions made to tackle these challenges are presented. The discussion focuses first on the "single-physics" solver and later on its extension to the monolithic coupling of different physics. The implementation of different physics modules, relevant to the capsule/parachute system, are also presented. Finally, examples of coupled computations are presented, paving the way to the simulation of the full capsule/parachute system.

Carton De Wiart, Corentin↗

NASA NTRS1976

Optimal control of spin stabilized spacecraft with telescoping appendages

The control of a spin-stabilized spacecraft consisting of a rigid central hub and one or two movable offset telescoping booms (with end masses) is considered. The equations of rotational motion are linearized about either of two desired final states. A control law for the boom and mass position is sought such that a quadratic cost functional involving the weighted components of angular velocity plus the control is minimized when the final time is unspecified and involves the solution of the matrix Riccati algebraic equation. For three-axis control more than one offset boom (orthogonal to each other) is required. For two-axis control with a single boom offset from a symmetrical hub, an analytic solution is obtained; when this system is used for nutation decay the time constant is one order of magnitude smaller than previously achieved using non-optimal control logic. For the general case results are obtained numerically.

Bainum, P. M.↗

DOE OSTI2020

sta-imex (scripts to investigate stability of IMEX methods)

Set of Matlab scripts to investigate properties of IMEX time stepping methods. Scripts include discretization of a system of linearized equations. One script plots numerical dispersion against analytical dispersion of the system. Another set of scripts uses the system to break it into stiff and non stiff part to investigate stability of Implicit-Explicit Runge-Kutta methods. Scripts are based on Thuburn, J. and Woollings, T. J., Vertical Discretizations for Compressible Euler Equation Atmospheric Models Giving Optimal Representation of Normal Modes , 2005 and Lock, S.-J. and Wood, N. and Weller, H., Numerical analyses of Runge-Kutta implicit-explicit schemes for horizontally explicit, vertically implicit solutions of atmospheric models, 2014 expanding ideas from these two papers. Scripts are used to develop IMEX schemes. SAND2020-12594 M Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy's National Nuclear Security Administration under contract DE-NA0003525.

Guba, Oksana↗

NASA NTRS2017

Navigating the Deployment and Downlink Tradespace for Earth Imaging Constellations

Distributed Spacecraft Missions (DSMs) are gaining momentum in their application to Earth Observation (EO) missions owing to their unique ability to increase observation sampling in spatial, spectral, angular and temporal dimensions simultaneously. DSM design includes a much larger number of variables than its monolithic counterpart, therefore, Model-Based Systems Engineering (MBSE) has been often used for preliminary mission concept designs, to understand the trade-offs and interdependencies among the variables. MBSE models are complex because the various objectives a DSM is expected to achieve are almost always conflicting, non-linear and rarely analytical. NASA Goddard Space Flight Center (GSFC) is developing a pre-Phase A tool called Tradespace Analysis Tool for Constellations (TAT-C) to initiate constellation mission design. The tool will allow users to explore the tradespace between various performance, cost and risk metrics (as a function of their science mission) and select Pareto optimal architectures that meet their requirements. This paper will describe the different types of constellations that TAT-C’s Tradespace Search Iterator is capable of enumerating (homogeneous Walker, heterogeneous Walker, precessing type, ad-hoc) and their impact on key performance metrics such as revisit statistics, time to global access and coverage. We will also discuss the ability to simulate phased deployment of the given constellations, as a function of launch availabilities and/or vehicle capability, and show the impact on performance. All performance metrics are calculated by the Data Reduction and Metric Computation module within TAT-C, which issues specific requests and processes results from the Orbit and Coverage module. Our TSI is also capable of generating tradespaces for downlinking imaging data from the constellation, based on permutations of available ground station networks - known (default) or customized (by the user). We will show the impact of changing ground station options for any given constellation, on data latency and required communication bandwidth, which in turn determines the responsiveness of the space system.

Nag, Sreeja↗

DOE OSTI2020

Development and Optimization of a Purification Process to Recover 99 Mo from Low-enriched Uranium

Supported by the DOE-NNSA Office of Material Management & Minimization (M3), a number of domestic entities are pursuing non-highly enriched uranium (non-HEU) production of 99 Mo. As the production technologies of 99Mo pivot toward low-enriched uranium (LEU) or molybdenum targets, new reaction channels and accelerators are being evaluated. Superconducting electron linear accelerators (LINACs) with high-Z converter targets can generate bremsstrahlung photons and neutron fluxes that are capable of inducing photonuclear reactions and LEU fission. A particular advantage of a LINAC is that it does not rely on HEU-fueled reactor cores (which are currently slated for LEU conversion) and can operate on an almost continuous basis. Regarding the chemical purification of 99 Mo from irradiated uranium targets under acidic digestion, there exists a procedure known as Cintichem or modifications thereof with respect to LEU (LEU Modified Cintichem process, LMC). The process relies on a number of selective precipitation steps and column chromatography to purify Mo. It is important to note that LMC prescribes the addition of stable Mo to carry 99 Mo on alpha benzoin oxime, which reduces the specific activity of 99 Mo. This is especially important for processing 99 Mo batches with lower activities (~33 Ci of 99 Mo per batch).

07 ISOTOPE AND RADIATION SOURCES↗

DOE OSTI2023

Development of a continuous synthesis process for carbamazepine using validated in-line Raman spectroscopy and kinetic modelling for disturbance simulation

Mitigation of failure modes in the continuous synthesis (CS) of a drug substance (DS) has the potential to widen the adoption of continuous manufacturing (CM) technologies by the pharmaceutical industry. Here, this work demonstrates the development of a robust continuous process for the synthesis of carbamazepine (CBZ), an essential medicine as per the World Health Organization (WHO), facilitated by kinetic modelling and monitored by in-line Raman spectroscopy. Accurate kinetic modelling and the use of validated process analytical technology (PAT) models for quantitative measurement were found to play an important role in developing CS of drug substances. Kinetic data for the formation of CBZ from iminostilbene (ISB) were collected by batch reaction sampling and high-performance liquid chromatography (HPLC) analysis. A non-linear solver and iterative method was applied to determine two sets of Arrhenius parameters simultaneously for the reaction system by minimizing the standard error of the model fit. The start-up and dynamic equilibrium stages for the CS of CBZ using a continuous stirred tank reactor (CSTR) were modelled based on the batch kinetic data and employed to optimize conversion and simulate process disturbances. An in-line Raman spectroscopy method was successfully developed, validated, and integrated to determine the concentrations of CBZ and ISB within the operating range for the CS. The CS kinetic model was evaluated experimentally from startup to dynamic equilibrium over 10 residence times with monitoring by HPLC and in-line Raman spectroscopy. The developed kinetic model in tandem with in-line Raman spectroscopy successfully predicted disturbances due to changes in process variables and can serve as a useful tool in the future design of advanced process control strategies for the continuous synthesis of CBZ.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Increasing the power of weak lensing data with multipole-based intrinsic alignment estimators

ABSTRACT It has long been known that galaxy shapes align coherently with the large-scale density field. Characterizing this effect is essential to interpreting measurements of weak gravitational lensing, the deflection of light from distant galaxies by matter overdensities along the line of sight, as it also produces coherent galaxy alignments that we wish to interpret in terms of a cosmological model. Existing direct measurements of intrinsic alignments using galaxy samples with high-quality shape and redshift measurements typically use well-understood but sub-optimal projected estimators, which do not make good use of the information in the data when comparing those estimators to theoretical models. We demonstrate a more optimal estimator, based on a multipole expansion of the correlation functions or power spectra, for direct measurements of galaxy intrinsic alignments. We show that even using the lowest order multipole alone increases the significance of inferred model parameters using simulated and real data, without any additional modelling complexity. We apply this estimator to measurements of parameters of the non-linear alignment model using data from the Sloan Digital Sky survey, demonstrating consistent results with a factor of ∼2 greater precision in parameter fits to intrinsic alignments models. This result is functionally equivalent to quadrupling the survey area, but without the attendant costs – thereby demonstrating the value in using this new estimator in current and future intrinsic alignments measurements using spectroscopic galaxy samples.

Astronomy & Astrophysics↗

NASA NTRS2016

In-Flight Performance of the Polarization Modulator in the CLASP Rocket Experiment

We developed a polarization modulation unit (PMU), a motor system to rotate a waveplate continuously. In polarization measurements, the continuous rotating waveplate is an important element as well as a polarization analyzer to record the incident polarization in a time series of camera exposures. The control logic of PMU was originally developed for the next Japanese solar observation satellite SOLAR-C by the SOLAR-C working group. We applied this PMU for the Chromospheric Lyman‐alpha SpectroPolarimeter (CLASP). CLASP is a sounding rocket experiment to observe the linear polarization of the Lyman‐alpha emission (121.6 nm vacuum ultraviolet) from the upper chromosphere and transition region of the Sun with a high polarization sensitivity of 0.1 % for the first time and investigate their vector magnetic field by the Hanle effect. The driver circuit was developed to optimize the rotation for the CLASP waveplate (12.5 rotations per minute). Rotation non‐uniformity of the waveplate causes error in the polarization degree (i.e. scale error) and crosstalk between Stokes components. We confirmed that PMU has superior rotation uniformity in the ground test and the scale error and crosstalk of Stokes Q and U are less than 0.01 %. After PMU was attached to the CLASP instrument, we performed vibration tests and confirmed all PMU functions performance including rotation uniformity did not change. CLASP was successfully launched on September 3, 2015, and PMU functioned well as designed. PMU achieved a good rotation uniformity, and the high precision polarization measurement of CLASP was successfully achieved.

CLASP↗

NASA NTRS2016

Streamlining the Design Tradespace for Earth Imaging Constellations

Science Data Processing↗

NASA NTRS2019

KELT-22Ab: A Massive, Short-Period Hot Jupiter Transiting a Near-solar Twin

We present the discovery of KELT-22Ab, a hot Jupiter from the KELT-South survey. KELT-22Ab transits the moderately bright (V11.1) Sun-like G2V star TYC 7518-468-1. The planet has an orbital period of P = 1.3866529 0.0000027 days, a radius of = -R 1.285+ R P 0.071 J 0.12 , and a relatively large mass of = - M 3.47+ M P 0.14 J 0.15 . The star has = - + R 1.099 0.046 R 0.079 , = - + M 1.092 0.041 M 0.045 , = - T 5767+ eff 49 50 K, = - log g 4.393+0.0600.039 (cgs), and [m/H]=+ -0.259+0.0830.085; thus other than its slightly super-solar metallicity, it appears to be a near-solar twin. Surprisingly, KELT-22A exhibits kinematics and a Galactic orbit that are somewhat atypical for thin-disk stars. Nevertheless, the star is rotating rapidly for its estimated age, and shows evidence of chromospheric activity. Imaging reveals a slightly fainter companion to KELT-22A that is likely bound, with a projected separation of 6″ (∼1400 au). In addition to the orbital motion caused by the transiting planet, we detect a possible linear trend in the radial velocity of KELT-22A, suggesting the presence of another relatively nearby body that is perhaps non-stellar. KELT-22Ab is highly irradiated (as a consequence of the small semimajor axis of a R = 4.97), and is mildly inflated. At such small separations, tidal forces become significant. The configuration of this system is optimal for measuring the rate of tidal dissipation within the host star. Our models predict that, due to tidal forces, the semimajor axis is decreasing rapidly, and KELT-22Ab is predicted to spiral into the star within the next Gyr.

Knicole Colon↗

DOE OSTI2026

Global energy balance simulations of shattered pellet injection

The dynamics of shattered pellet injection (SPI) shutdowns are simulated using a time-dependent global energy balance model, based on a modification of the KPRAD framework. The new SPI particle source in the model calculates the ablation of individual pellet fragments that enter the plasma as a temporally resolved plume, thus capturing the effects of earlier fragments on the ablation of those that follow, which has a significant impact on the overall assimilation. Despite the reduced physics and the global averaging of all quantities, results from a large number of DIII-D, KSTAR, and JET experiments are well reproduced, including the plasma cooling timescales, particle assimilation, and current quench (CQ) rates. Cooling timescales and CQ rates are in good agreement for pellets containing as little as ∼1% neon, while particle assimilations are most accurate for neon fractions above ∼15% by number of atoms. Below this, the assimilation tends to be overestimated due to the lack of radial transport in the particle balance, which becomes important in the low-Z limit. Predictive simulations of mixed-composition dual-SPI shutdowns in ITER are compared against those with the 3D non-linear magnetohydrodynamic code JOREK and are found to reproduce overall trends observed in the higher-fidelity modeling across a range of injection scenarios. The general success of the model points to the critical role of energy balance in determining SPI particle assimilation and the subsequent disruption dynamics and highlights the value of these simulations for experimental interpretation and for optimizing the deployment of computationally expensive, higher-fidelity models.

energy balance model↗

DOE OSTI2020

Reification of latent microstructures: On supervised unsupervised and semi-supervised deep learning applications for microstructures in materials informatics

Machine learning (ML), including deep learning (DL), has become increasingly popular in the last few years due to its continually outstanding performance. In this context, we apply machine learning techniques to "learn" the microstructure using both supervised and unsupervised DL techniques. In particular, we focus (1) on the localization problem bridging (micro)structure (localized) property using supervised DL and (2) on the microstructure reconstruction problem in latent space using unsupervised DL. The goal of supervised and semi-supervised DL is to replace crystal plasticity finite element model (CPFEM) that maps from (micro)structure (localized) property, and implicitly the (micro)structure (homogenized) property relationships, while the goal of unsupervised DL is (1) to represent high-dimensional microstructure images in a non-linear low-dimensional manifold, and (2) to discover a way to interpolate microstructures via latent space associating with latent microstructure variables. At the heart of this report is the applications of several common DL architectures, including convolutional neural networks (CNN), autoencoder (AE), and generative adversarial network (GAN), to multiple microstructure datasets, and the quest of neural architecture search for optimal DL architectures.

36 MATERIALS SCIENCE↗

NASA NTRS2011

Antenna Near-Field Probe Station Scanner

A miniaturized antenna system is characterized non-destructively through the use of a scanner that measures its near-field radiated power performance. When taking measurements, the scanner can be moved linearly along the x, y and z axis, as well as rotationally relative to the antenna. The data obtained from the characterization are processed to determine the far-field properties of the system and to optimize the system. Each antenna is excited using a probe station system while a scanning probe scans the space above the antenna to measure the near field signals. Upon completion of the scan, the near-field patterns are transformed into far-field patterns. Along with taking data, this system also allows for extensive graphing and analysis of both the near-field and far-field data. The details of the probe station as well as the procedures for setting up a test, conducting a test, and analyzing the resulting data are also described.

Zaman, Afroz J.↗

NASA NTRS2016

Nonlinear Transient Growth and Boundary Layer Transition

Parabolized stability equations (PSE) are used in a variational approach to study the optimal, non-modal disturbance growth in a Mach 3 at plate boundary layer and a Mach 6 circular cone boundary layer. As noted in previous works, the optimal initial disturbances correspond to steady counter-rotating streamwise vortices, which subsequently lead to the formation of streamwise-elongated structures, i.e., streaks, via a lift-up effect. The nonlinear evolution of the linearly optimal stationary perturbations is computed using the nonlinear plane-marching PSE for stationary perturbations. A fully implicit marching technique is used to facilitate the computation of nonlinear streaks with large amplitudes. To assess the effect of the finite-amplitude streaks on transition, the linear form of plane- marching PSE is used to investigate the instability of the boundary layer flow modified by spanwise periodic streaks. The onset of bypass transition is estimated by using an N- factor criterion based on the amplification of the streak instabilities. Results show that, for both flow configurations of interest, streaks of sufficiently large amplitude can lead to significantly earlier onset of transition than that in an unperturbed boundary layer without any streaks.

Laminar turbulent transition↗

DOE OSTI2026

Constraining primordial non-Gaussianity from DESI DR1 quasars and Planck PR4 CMB lensing

We present the first measurement of local-type primordial non-Gaussianity from the cross-correlation between 1.2 million spectroscopically confirmed quasars from the first data release (DR1) of the Dark Energy Spectroscopic Instrument (DESI) and the Planck PR4 CMB lensing reconstructions. The analysis is performed in three tomographic redshift bins covering 0.8 < z < 3.5, covering a sky fraction of ∼20%. We adopt a catalog-based pseudo-C ℓ estimator and apply linear imaging weights validated on noiseless mocks. Compared to previous analyses using photometric quasar samples, our results benefit from the high purity of the DESI spectroscopic sample, the reduced noise of PR4 lensing, and the absence of excess large-scale power in the spectroscopic quasar auto-correlation. Fitting simultaneously for the non-Gaussianity parameter f NL and the linear bias amplitude in each redshift bin, we obtain f NL = 2 +28 -34 for a response parameter p = 1.6, and f NL = 6 +20 -24 for p = 1.0. These results improve the constraints on f NL by ∼35% compared to the previous analysis based on the Legacy Imaging Survey DR9. Additionally, we derive an optimal weighting scheme to maximize the constraining power. In this case, and assuming p = 1.6, we obtain f NL = 19 +25 -31 . Our results demonstrate the statistical power of DESI quasars for probing inflationary physics, and highlight the promise of future DESI data releases.

cosmological parameters from CMBR↗

DOE OSTI2021

Topological network analysis of patient similarity for precision management of acute blood pressure in spinal cord injury

Background: Predicting neurological recovery after spinal cord injury (SCI) is challenging. Using topological data analysis, we have previously shown that mean arterial pressure (MAP) during SCI surgery predicts long-term functional recovery in rodent models, motivating the present multicenter study in patients. Methods: Intra-operative monitoring records and neurological outcome data were extracted (n = 118 patients). We built a similarity network of patients from a low-dimensional space embedded using a non-linear algorithm, Isomap, and ensured topological extraction using persistent homology metrics. Confirmatory analysis was conducted through regression methods. Results: Network analysis suggested that time outside of an optimum MAP range (hypotension or hypertension) during surgery was associated with lower likelihood of neurological recovery at hospital discharge. Logistic and LASSO (least absolute shrinkage and selection operator) regression confirmed these findings, revealing an optimal MAP range of 76–[104-117] mmHg associated with neurological recovery. Conclusions: We show that deviation from this optimal MAP range during SCI surgery predicts lower probability of neurological recovery and suggest new targets for therapeutic intervention. Funding: NIH/NINDS: R01NS088475 (ARF); R01NS122888 (ARF); UH3NS106899 (ARF); Department of Veterans Affairs: 1I01RX002245 (ARF), I01RX002787 (ARF); Wings for Life Foundation (ATE, ARF); Craig H. Neilsen Foundation (ARF); and DOD: SC150198 (MSB); SC190233 (MSB); DOE: DE-AC02-05CH11231 (DM).

59 BASIC BIOLOGICAL SCIENCES↗