Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “parameter 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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
NASA NTRS1984

Optimal design of GaAs-based concentrator space solar cells for 100 AMO, 80 deg C operation

Using a detailed computer code and reasonable values of electrical and optical material parameters from current published literature, parameter optimization studies were performed on three configurations of GaAs-based concentrator solar cells for 100 AMO, 80 C operation. These studies show the possibility of designing GaAs-based solar cells with efficiencies exceeding 22% at 100 AMO 80 C and probable efficiency degradation of less than 15% after a 70% reduction in diffusion length in each cell region.

Goradia, C.↗
NASA NTRS2009

Determining Spacecraft Reaction Wheel Friction Parameters

Software was developed to characterize the drag in each of the Cassini spacecraft's Reaction Wheel Assemblies (RWAs) to determine the RWA friction parameters. This tool measures the drag torque of RWAs for not only the high spin rates (greater than 250 RPM), but also the low spin rates (less than 250 RPM) where there is a lack of an elastohydrodynamic boundary layer in the bearings. RWA rate and drag torque profiles as functions of time are collected via telemetry once every 4 seconds and once every 8 seconds, respectively. Intermediate processing steps single-out the coast-down regions. A nonlinear model for the drag torque as a function of RWA spin rate is incorporated in order to characterize the low spin rate regime. The tool then uses a nonlinear parameter optimization algorithm based on the Nelder-Mead simplex method to determine the viscous coefficient, the Dahl friction, and the two parameters that account for the low spin-rate behavior.

Sarani, Siamak↗
DOE OSTI2024

Measurement of the Pion Exclusive Electro-Production Cross-Section in the E12-19-006 Experiment in Hall-C at Jefferson Lab

One of the most effective methods for exploring the transition from hadronic degrees of freedom to quark-gluon degrees of freedom in Quantum Chromodynamics (QCD) is through the investigation of \exclusive" pion and kaon electro-production reactions at various Q2 and ?t values. The E12-19-006 experiment is conducted within the confi?nes of experimental Hall C at the Thomas Jefferson National Accelerator Facility, USA, for such studies. The primary aim of the experiment is to ?first enhance our comprehension of the pion electro-production cross-section and its form factor at Q2 = 0.38 and 0.42 GeV2. This is the fi?rst run period of the E12-19-006 experiment which ran in summer 2019. A more profound understanding of the pion electro-production reaction, 1H(e,e'?+)n, at low Q2 is deemed essential to employ this electro-production reaction (an indirect technique) for the high Q2 studies, thereby delving deeper into the realm of QCD. Consequently, this dissertation presents a thorough analysis of the experimental data acquired in the ?first run period of the E12-19-006 experiment. In pursuit of precision, a series of systematic studies (target boiling correction study, the elastic reaction cross-section measurements, study for determining vari- ous kinematics o?sets, etc.) are conducted to discern the accuracy of the analyzed data, a prerequisite for the use of Rosenbluth separation technique to separate the pion electro-production cross-section terms in t bins. The separated pion electro-production cross-section through the Rosenbluth separation technique is then used to extract the pion electromagnetic form factor. In this dissertation, the pion electro-production cross-section is carefully dissected into its four constituent components: longitudinal (?L), transverse (?T ), longitudinal-transverse (?LT ), and transverse-transverse (?TT ), using the full version of Rosenbluth separation technique for the Q2 = 0.38 GeV2. The technique is simultaneously fi?tted to the unseparated pion electro-production cross-sections at the three values of polarization of the virtual photon (?), i.e., ? = 0.286, 0.629 and 0.781. An iterative process is applied to re?ne the parameters of the model cross-sections until the yield ratio of experimental and Monte Carlo simulation converges. In this study, 21 iterations are conducted to re?ne the model cross-section parameters. The fi?nal pion electro-production cross-section terms are then determined for 7 t bins using the optimized parameters of the model cross-sections.

Kumar, Vijay↗
DOE OSTI2021

Engineering a Non‐Natural Photoenzyme for Improved Photon Efficiency**

Abstract Photoenzymes are biological catalysts that use light to convert starting materials into products. These catalysts require photon absorption for each turnover, making quantum efficiency an important optimization parameter. Flavin‐dependent “ene”‐reductases (EREDs) display latent photoenzymatic activity for synthetically valuable hydroalkylations; however, protein engineering has not been used to optimize this non‐natural function. We describe a protein engineering platform for the high throughput optimization of photoenzymes. A single round of engineering results in improved catalytic function toward the synthesis of γ, δ, ϵ‐lactams, and acyclic amides. Mechanistic studies show that key mutations can alter the enzyme's excited state dynamics, enhance its photon efficiency, and ultimately increase catalyst performance. Transient absorption spectroscopy reveals that engineered variants display dramatically decreased radical lifetimes, indicating an evolution toward a concerted mechanism.

Nicholls, Bryce T.↗
DOE OSTI2020

Modeling of fluid injection-induced fault reactivation using coupled fluid flow and mechanical interface model

The present study is aimed at developing a numerical model to reproduce coupled hydro-mechanical processes associated with fault reactivation by fluid injection in low permeability rock, as part of the DECOVALEX-2019 project Task B. We proposed a modeling approach for simulating the processes using the TOUGH-FLAC simulator, and modeled a fault reactivation experiment conducted at Mont Terri Rock Laboratory in Switzerland. The first step of the study involved benchmark calculations considering a simplified fault plane and geometry. Fluid flow along a fault was modeled using elements of aperture-sized thickness on the basis of Darcy's law and the cubic law in TOUGH2, whereas the mechanical behavior of a single fault was represented by zero-thickness interface elements in FLAC3D upon which a slip and/or separation is allowed. A methodology to connect a TOUGH2 volume element to a FLAC3D interface element was developed for handling the hydro-mechanical interactions on the fault during fluid injection. Two different fault models for describing the evolutions of hydraulic aperture by elastic fracture opening and failure-induced aperture increase were considered in the benchmark calculations. In the coupling process, the changes in geometrical features and hydrological properties induced by mechanical deformation were continuously updated. The transient responses of the fault and host rock to stepwise pressurization were examined during the simulation. The hydro-mechanical behavior, including the injection flow rate, pressure distribution around the borehole, stress conditions, and displacements in normal and shear directions were monitored in the surrounding rock and along the fault. The results of benchmark calculations suggest that the developed model reasonably represents the hydro-mechanical behavior of a fault and the surrounding rock. This modeling approach was applied to the fault reactivation experiment of the Mont Terri Rock Laboratory. In this interpretive modeling, a parametric study was conducted to examine the effects of input parameters regarding in situ stress and fault properties on the hydro-mechanical responses of the fault to water injection. Then, an optimal parameter set to reproduce the field experiment results was chosen by trial-and-error. The injection flow rate and pressure response during fault reactivation closely matched those obtained at the site, which indicates the capability of the model to appropriately capture the progressive pathway evolution during fault reactivation tests at the site. Finally, the anchor displacements were overestimated by the model, but a fair agreement was obtained in terms of the order of magnitude and the variation tendency.

58 GEOSCIENCES↗
DOE OSTI2024

An extensible perturbed parameter ensemble for the Community Atmosphere Model version 6

This paper documents the methodology and preliminary results from a perturbed parameter ensemble (PPE) technique, where multiple parameters are varied simultaneously and the parameter values are determined with Latin hypercube sampling. This is done with the Community Atmosphere Model version 6 (CAM6), the atmospheric component of the Community Earth System Model version 2 (CESM2). We apply the PPE method to CESM2–CAM6 to understand climate sensitivity to atmospheric physics parameters. The initial simulations vary 45 parameters in the microphysics, convection, turbulence and aerosol schemes with 263 ensemble members. These atmospheric parameters are typically the most uncertain in many climate models. Control simulations and targeted simulations to understand climate forcing due to aerosols and fast climate feedbacks are analyzed. The use of various emulators is explored in the multi-dimensional space mapping input parameters to output metrics. Parameter impacts on various model outputs, such as radiation, cloud and aerosol properties, are evaluated. Machine learning is also used to probe optimal parameter values against observations. Our findings show that PPE is a valuable tool for climate uncertainty analysis. Furthermore, by varying many parameters simultaneously, we find that many different combinations of parameter values can produce results consistent with observations, and thus careful analysis of tuning is important. The CESM2–CAM6 PPE is publicly available and extensible to other configurations to address questions of other model processes in the atmosphere and other model components (e.g., coupling to the land surface).

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2024

An Extensible Perturbed Parameter Ensemble for the Community Atmosphere Model Version 6

This paper documents the methodology and preliminary results from a Perturbed Parameter Ensemble (PPE) technique, where multiple parameters are varied simultaneously and the parameter values are determined with Latin hypercube sampling. This is done with the Community Atmosphere Model version 6 (CAM6), the atmospheric component of the Community Earth System Model version 2 (CESM2). We apply the PPE method to CESM2-CAM6 to understand climate sensitivity to atmospheric physics parameters. The initial simulations vary 45 parameters in the microphysics, convection, turbulence and aerosol schemes with 263 ensemble members. These atmospheric parameters are typically the most uncertain in many climate models. Control simulations are analyzed and targeted simulations to understand climate forcing due to aerosols and fast climate feedbacks. The use of various emulators is explored in the multi- dimensional space mapping input parameters to output metrics. Parameter impacts on various model outputs, such as radiation, cloud and aerosol properties are evaluated. Machine learning is also used to probe optimal parameter values against observations. Our findings show that using PPE is a valuable tool for climate uncertainty analysis. Furthermore, by varying many parameters simultaneously, we find that many different combinations of parameter values can produce results consistent with observations, and thus careful analysis of tuning is important. The CESM2-CAM6 PPE is publicly available, and extensible to other configurations to address questions of other model processes in the atmosphere and other model components (e.g. coupling to the land surface).

Machine learning↗
DOE OSTI2022

Reactive burn model calibration using high-throughput initiation experiments at sub-millimeter length scales

We report a first-of-its-kind model calibration was performed using Sandia National Laboratories’ high-throughput initiation (HTI) experiment for two types of vapor-deposited explosive films consisting of hexanitrostilbene (HNS) or pentaerythritol tetranitrate (PETN). These films exhibit prompt initiation, and they reach steady detonation at sub-millimeter length scales. Following prior work on HNS, we test the hypothesis of approximating these explosive films as fine-grained homogeneous solids with simple Arrhenius kinetics burn models. The model calibration process is described herein using a single-step as well as a two-step Arrhenius rate law, and it consists of systematic parameter sampling leading to a reduction in the model degrees of freedom. Multiple local minima are observed; results are given for seven different optimized parameter sets. Each model set is further evaluated in a two-dimensional simulation of the critical failure thickness for a sustained detonation. Overall, the two-step Arrhenius kinetics model captures the observed behavior for HNS; however, neither model produces a good fit to the PETN data. We hypothesize that the HTI results for PETN correspond to a heterogeneous response, owing to the smaller reaction zone of PETN compared to HNS (i.e., it does not homogenize the fine-grained hot spots as well). Future work should consider using the ignition and growth model for PETN, as well as other reactive burn models such as xHVRB, AWSD, PiSURF, and CREST.

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

A machine learning approach to emulation and biophysical parameter estimation with the Community Land Model, version 5

Abstract. Land models are essential tools for understanding and predicting terrestrial processes and climate–carbon feedbacks in the Earth system, but uncertainties in their future projections are poorly understood. Improvements in physical process realism and the representation of human influence arguably make models more comparable to reality but also increase the degrees of freedom in model configuration, leading to increased parametric uncertainty in projections. In this work we design and implement a machine learning approach to globally calibrate a subset of the parameters of the Community Land Model, version 5 (CLM5) to observations of carbon and water fluxes. We focus on parameters controlling biophysical features such as surface energy balance, hydrology, and carbon uptake. We first use parameter sensitivity simulations and a combination of objective metrics including ranked global mean sensitivity to multiple output variables and non-overlapping spatial pattern responses between parameters to narrow the parameter space and determine a subset of important CLM5 biophysical parameters for further analysis. Using a perturbed parameter ensemble, we then train a series of artificial feed-forward neural networks to emulate CLM5 output given parameter values as input. We use annual mean globally aggregated spatial variability in carbon and water fluxes as our emulation and calibration targets. Validation and out-of-sample tests are used to assess the predictive skill of the networks, and we utilize permutation feature importance and partial dependence methods to better interpret the results. The trained networks are then used to estimate global optimal parameter values with greater computational efficiency than achieved by hand tuning efforts and increased spatial scale relative to previous studies optimizing at a single site. By developing this methodology, our framework can help quantify the contribution of parameter uncertainty to overall uncertainty in land model projections.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

A Behavioral Robotics Approach to Radiation Mapping Using Adaptive Sampling

Radiation mapping is a desirable task to automate because of the inherent risks involved and its tedious nature. A novel system was designed to address this by combining various existing technologies, utilizing behavior-based robotics and Bayesian optimization. The system uses a quadruped robot equipped with a manipulator and gamma detector to take measurements at locations that are selected based on the uncertainty of a surrogate model used to estimate the true radiation field. The robot uses input from the world with depth cameras to avoid collisions with the robot’s body, and unreachable points for the end effector are addressed by both allowing for a soft collision with the environment to occur, prompting the system to abandon that point, and varying the exploration tendency of the optimization based on consecutive collisions. This approach provides unique traversability and adaptability over other strategies in the literature. Experiments were performed by placing a Cesium-137 source on the ground and varying geometric setups and an optimization parameter demonstrating the adaptability to diverse environments and the increased robustness resulting from the designed behavior. The results additionally demonstrate that dynamically adjusting the optimization algorithm’s exploration tendency based on the arm’s collision history improves the system’s ability to navigate cluttered environments and construct accurate radiation maps without getting stuck in unreachable areas.

Adams, Joel↗
NASA NTRS1988

Attitude control system testing on SCOLE

This paper presents implementation of two control policies on SCOLE (Space Control Laboratory Experiment), a laboratory apparatus representing an offset-feed antenna attached to the Space Shuttle by a flexible mast. In the first case, the flexible mast was restrained by cables, permitting modeling of SCOLE as a rigid-body. Starting from an arbitrary state, SCOLE was maneuvered to a specified terminal state using rigid-body minimum-time control law. In the second case, the so called single step optimal control (SSOC) theory is applied to suppress vibrations of the flexible mast mounted as a cantilever beam. Based on the SSOC theory, two parameter optimization algorithms were developed.

Shenhar, J.↗
NASA NTRS1993

A homotopy algorithm for digital optimal projection control GASD-HADOC

The linear-quadratic-gaussian (LQG) compensator was developed to facilitate the design of control laws for multi-input, multi-output (MIMO) systems. The compensator is computed by solving two algebraic equations for which standard closed-loop solutions exist. Unfortunately, the minimal dimension of an LQG compensator is almost always equal to the dimension of the plant and can thus often violate practical implementation constraints on controller order. This deficiency is especially highlighted when considering control-design for high-order systems such as flexible space structures. This deficiency motivated the development of techniques that enable the design of optimal controllers whose dimension is less than that of the design plant. A homotopy approach based on the optimal projection equations that characterize the necessary conditions for optimal reduced-order control. Homotopy algorithms have global convergence properties and hence do not require that the initializing reduced-order controller be close to the optimal reduced-order controller to guarantee convergence. However, the homotopy algorithm previously developed for solving the optimal projection equations has sublinear convergence properties and the convergence slows at higher authority levels and may fail. A new homotopy algorithm for synthesizing optimal reduced-order controllers for discrete-time systems is described. Unlike the previous homotopy approach, the new algorithm is a gradient-based, parameter optimization formulation and was implemented in MATLAB. The results reported may offer the foundation for a reliable approach to optimal, reduced-order controller design.

Collins, Emmanuel G., Jr.↗
DOE OSTI2020

Land Model Testbed: Accelerating Development, Benchmarking and Analysis of Land Surface Models

A Land Model Testbed (LMT), designed to provide a computational framework for systematically assessing model fidelity and supporting rapid development of complex multiscale models, offers a general-purpose workflow for conducting large ensemble simulations of multiple land surface models, post-processing large volumes of model output, and evaluating model results. It leverages existing tools for launching model simulations and the International Land Model Benchmarking (ILAMB) package for assessing model fidelity through comparison with best-available observational datasets. Increased complexity and proliferation of uncertain parameters in process representations in land surface models has driven the need for frequent and intensive testing and evaluating of models to quantify uncertainties and optimize parameters such that results are consistent with observations. The LMT described here meets these needs by providing tools to run thousands of ensemble simulations simultaneously and post-process their output files, by automating execution of an enhanced version of ILAMB with site-specific benchmarks and multivariate functional relationships, and by offering ensemble diagnostics and a customizable dashboard for displaying model performance metrics and associated graphics. We envision the LMT capabilities will serve as a foundational computational resource for a proposed user facility focused on terrestrial multiscale model--data integration.

Sreepathi, Sarat↗
NASA NTRS1990

TROJID: A portable software package for upper-stage trajectory optimization

Performance optimization for upper-stage exoatmospheric vehicles often is performed within the framework of a full capability trajectory simulation package requiring either a large mainframe computer or powerful work-station. Since these software packages tend to include capabilities providing for high-fidelity boost and reentry simulations, the programs usually are quite large and not very portable. The program TROJID is an attempt to provide an environment for the optimization of upper-stage trajectories within a small package capable of being run on a standard desktop microcomputer. Utilizing a state-of-the-art nonlinear programming algorithm and a trajectory simulator implementing impulsive burns and an analytic coast phase propagator, TROJID is capable of producing trajectories for optimal multi-burn upper-stage orbit transfers. The package has been designed to allow full generality in definition of both the trajectory simulator and the parameter optimization problem.

Hammes, Steven M.↗
NASA NTRS2004

OTIS 3.2 Software Released

Trajectory, mission, and vehicle engineers concern themselves with finding the best way for an object to get from one place to another. These engineers rely upon special software to assist them in this. For a number of years, many engineers have used the OTIS program for this assistance. With OTIS, an engineer can fully optimize trajectories for airplanes, launch vehicles like the space shuttle, interplanetary spacecraft, and orbital transfer vehicles. OTIS provides four modes of operation, with each mode providing successively stronger optimization capability. The most powerful mode uses a mathematical method called implicit integration to solve what engineers and mathematicians call the optimal control problem. OTIS 3.2, which was developed at the NASA Glenn Research Center, is the latest release of this industry workhorse and features new capabilities for parameter optimization and mission design. OTIS stands for Optimal Control by Implicit Simulation, and it is implicit integration that makes OTIS so powerful at solving trajectory optimization problems. Why is this so important? The optimization process not only determines how to get from point A to point B, but it can also determine how to do this with the least amount of propellant, with the lightest starting weight, or in the fastest time possible while avoiding certain obstacles along the way. There are numerous conditions that engineers can use to define optimal, or best. OTIS provides a framework for defining the starting and ending points of the trajectory (point A and point B), the constraints on the trajectory (requirements like "avoid these regions where obstacles occur"), and what is being optimized (e.g., minimize propellant). The implicit integration method can find solutions to very complicated problems when there is not a lot of information available about what the optimal trajectory might be. The method was first developed for solving two-point boundary value problems and was adapted for use in OTIS. Implicit integration usually allows OTIS to find solutions to problems much faster than programs that use explicit integration and parametric methods. Consequently, OTIS is best suited to solving very complicated and highly constrained problems.

Riehl, John P.↗
DOE OSTI2021

Follow the power—pathways to steady-state tokamak reactors

Abstract Diagramming pathways of dimensionless power is a potent method for extrapolating between operating points on present-day tokamaks and future burning plasma devices. The heat transport power, current drive power, H-mode threshold power and other ‘plasma physics’ powers can be expressed in dimensionally correct (or normalized) form as Pa 3/4 , where P is the power and a is the plasma minor radius, with the relative gyroradius ( ρ * ) dependence ranging from gyro-Bohm-like for transport ( P a 3 / 4 ∝ ρ * − 3 / 2 ), Bohm-like for current drive ( P a 3 / 4 ∝ ρ * − 5 / 2 ) and worse than Bohm-like for H-mode threshold ( P a 3 / 4 ∝ ρ * − 3 ). The D–T fusion power cannot be normalized in the same fashion since it is governed by nuclear physics, but at fixed B T it scales like P a 3 / 4 ∝ ρ * − 9 / 2 . Other ‘mixed physics’ powers can be incorporated into the dimensionless power framework by holding B T fixed in the same manner. Diagramming these dimensionless powers vs ρ * shows how the pathway to a steady-state reactor can be optimized relative to various operational boundaries. Using a steady-state hybrid discharge with β N = 3.2 from DIII-D as the starting point, a multi-parameter optimization finds an attractive pathway to steady-state operation on ITER using 76 MW of current drive power (fusion gain of Q fus = 8), along with a pathway to Q fus = 20 in a JET-sized steady-state reactor with B T = 10 T.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1977

Image data-processing system for solar astronomy

The paper describes an image data processing system (IDAPS), its hardware/software configuration, and interactive and batch modes of operation for the analysis of the Skylab/Apollo Telescope Mount S056 X-Ray Telescope experiment data. Interactive IDAPS is primarily designed to provide on-line interactive user control of image processing operations for image familiarization, sequence and parameter optimization, and selective feature extraction and analysis. Batch IDAPS follows the normal conventions of card control and data input and output, and is best suited where the desired parameters and sequence of operations are known and when long image-processing times are required. Particular attention is given to the way in which this system has been used in solar astronomy and other investigations. Some recent results obtained by means of IDAPS are presented.

Wilson, R. M.↗
NASA NTRS1994

Large-Scale Numerical Simulations of Human Motion

This paper examines the feasibility of using massively-parallel and vector-processing supercomputers to solve large-scale optimal control problems for human movement. Specifically, we compare the computational expense of determining the optimal controls for the single support phase of walking using a conventional serial machine (a Silicon Graphics Personal Iris 4D25 workstation), a MIMD parallel machine (an Intel iPSC/860 comprising 128 processors), and a parallel-vector-processing machine (a Cray Y-MP 8/864). With the human body modeled as a 14 degree-of-freedom linkage actuated by 46 musculotendinous units, computation of the optimal controls for walking could take up to 3 months of CPU time on the Iris. Both the Cray Y-MP and the Intel iPSC/860 are able to reduce this time to practical levels. The optimal control solution for walking can be found with about 77 hours of CPU time on the Cray, and with about 88 hours of CPU time on the Intel. Although the overall speeds of the Cray and the Intel were found to be similar, the unique capabilities of each machine are best suited to different parts of the optimal control algorithm used. The Intel performed best in the calculation of the derivatives of the performance criterion and the constraints. In contrast, the Cray performed best during parameter optimization of the controls. These results suggest that the ideal computer architecture for solving very large-scale optimal control problems is a hybrid system in which a vector-processing machine is integrated into the communication network of a MIMD parallel machine.

Anderson, Frank C.↗