SEARCH · Engineering Papers
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.
It is hypothesized that mathematical models of the programming process can be formulated to gauge the sensitivities of that process to various given parameters, and that such models can be calibrated on an empirical basis and used as guides toward maximizing productivity, documentation quality, and programming reliability. Three oversimplified models are presented as illustrations.
The problem of structural optimization of helicopter rotor blades with integrated dynamic and aerodynamic design considerations is addressed. Results of recent optimization work on rotor blades for minimum weight with constraints on multiple coupled natural flap-lag frequencies, blade autorotational inertia and centrifugal stress has been reviewed. A strategy has been defined for the ongoing activities in the integrated dynamic/aerodynamic optimization of rotor blades. As a first step, the integrated dynamic/airload optimization problem has been formulated. To calculate system sensitivity derivatives necessary for the optimization recently developed, Global Sensitivity Equations (GSE) are being investigated. A need for multiple objective functions for the integrated optimization problem has been demonstrated and various techniques for solving the multiple objective function optimization are being investigated. The method called the Global Criteria Approach has been applied to a test problem with the blade in vacuum and the blade weight and the centrifugal stress as the multiple objectives. The results indicate that the method is quite effective in solving optimization problems with conflicting objective functions.
The problem of structural optimization of helicopter rotor blades with integrated dynamic and aerodynamic design considerations is addressed. Results of recent optimization work on rotor blades for minimum weight with constraints on multiple coupled natural flap-lag frequencies, blade autorotational inertia and centrifugal stress has been reviewed. A strategy has been defined for the ongoing activities in the integrated dynamic/aerodynamic optimization of rotor blades. As a first step, the integrated dynamic/airload optimization problem has been formulated. To calculate system sensitivity derivatives necessary for the optimization recently developed, Global Sensitivity Equations (GSE) are being investigated. A need for multiple objective functions for the integrated optimization problem has been demonstrated and various techniques for solving the multiple objective function optimization are being investigated. The method called the Global Criteria Approach has been applied to a test problem with the blade in vacuum and the blade weight and the centrifugal stress as the multiple objectives. The results indicate that the method is quite effective in solving optimization problems with conflicting objective functions.
Dynamic stability constraints are included in a computer program that simultaneously synthesizes a structure and its control system. Two measures of stability, the real part of the system complex eigenvalues and the damping ratio, are examined. The procedure for calculating the sensitivities of the two measures of stability to changes in the structure and its control system is explained. The sensitivities are used to formulate an approximate problem that is solved at each design iteration. The effects of structural damping and noncollated controllers on the synthesis process are discussed.
A gradient-based shape optimization methodology, that is intended for practical three-dimensional aerodynamic applications, has been developed. It is based on the quasi-analytical sensitivities. The flow analysis is rendered by a fully implicit, finite volume formulation of the Euler equations.The aerodynamic sensitivity equation is solved using the alternating-direction-implicit (ADI) algorithm for memory efficiency. A flexible wing geometry model, that is based on surface parameterization and platform schedules, is utilized. The present methodology and its components have been tested via several comparisons. Initially, the flow analysis for for a wing is compared with those obtained using an unfactored, preconditioned conjugate gradient approach (PCG), and an extensively validated CFD code. Then, the sensitivities computed with the present method have been compared with those obtained using the finite-difference and the PCG approaches. Effects of grid refinement and convergence tolerance on the analysis and shape optimization have been explored. Finally the new procedure has been demonstrated in the design of a cranked arrow wing at Mach 2.4. Despite the expected increase in the computational time, the results indicate that shape optimization, which require large numbers of grid points can be resolved with a gradient-based approach.
In order to re-examine trace gas transport in the middle atmosphere for the period May 1991 until April 1995, a "reanalysis" is being performed using an up-to-date version of the DAO's "GEOS" assimilation system. The Reanalysis for Stratospheric Trace-gas Studies (ReSTS) is intended to provide state-of-the-art estimates of the atmosphere during a period when the Upper Atmospheric Research Satellite provided a high density of trace-gas observations, and when the aerosol loading from the eruption of Mount Pinatubo contaminated the lower stratosphere, at the same time performing a natural tracer transport experiment. This study will present the first results from ReSTS, focussing on the improvements over the meteorological analyses produced by the then-operational GEOS-1 data assimilation system; emphasis will be placed on the improved representations of physical processes between GEOS-1 and the current GEOS-4 systems, highlighting the transport properties of the datasets. Alongside the production of a comprehensive atmospheric dataset, important components of ReSTS include performing sensitivity studies to the formulation of the assimilation system (including the representation of physical processes in the GCM, such as feedbacks between ozone/aerosols and meteorology) and to the inclusion of additional data types (including limb-sounding temperature data alongside the TOVS observations). Impacts of some of these factors on the analyzed meteorology and transport will be discussed. Of particular interest are attempts to determine the relative importance of various steps in the assimilation process to the quality of the final analyses.
Mesoscale convective systems (MCSs) are clusters of thunderstorms that are important in Earth's water and energy cycle. Additionally, they are responsible for extreme events such as large hail, strong winds, and extreme precipitation. Automated object-based analyses that track MCSs have become popular since they allow us to identify and follow MCSs over their entire life cycle in a Lagrangian framework. This rise in popularity was accompanied by an increasing number of MCS tracking algorithms, however, little is known about how sensitive analyses are concerning the MCS tracker formulation. Here, we assess differences between six MCS tracking algorithms on South American MCS characteristics and evaluate MCSs in kilometer-scale simulations with observational-based MCSs over 3 years. All trackers are run with a common set of MCS classification criteria to isolate tracker formulation differences. The tracker formulation substantially impacts MCS characteristics such as frequency, size, duration, and contribution to total precipitation. The evaluation of simulated MCS characteristics is less sensitive to the tracker formulation and all trackers agree that the model can capture MCS characteristics well across different South American climate zones. Dominant sources of uncertainty are the segmentation of cloud systems in space and time and the treatment of how MCSs are linked in time. Our results highlight that comparing MCS analyses that use different tracking algorithms is challenging. We provide general guidelines on how MCS characteristics compare between trackers to facilitate a more robust assessment of MCS statistics in future studies.
A new approach to integrated structure/control law design based on multilevel optimization is presented. This new approach is applicable to aircraft and spacecraft and allows for the independent design of the structure and control law. Integration of the designs is achieved through use of an upper level coordination problem formulation within the multilevel optimization framework. The method requires the use of structure and control law design sensitivity information. A general multilevel structure/control law design problem formulation is given, and the use of Linear Quadratic Gaussian (LQG) control law design and design sensitivity methods within the formulation is illustrated. Results of three simple integrated structure/control law design examples are presented. These results show the capability of structure and control law design tradeoffs to improve controlled system performance within the multilevel approach.
A new approach to integrated structure/control law design based on multilevel optimization is presented. This new approach is applicable to aircraft and spacecraft and allows for the independent design of the structure and control law. Integration of the designs is achieved through use of an upper level coordination problem formulation within the multilevel optimization framework. The method requires the use of structure and control law design sensitivity information. A general multilevel structure/control law design problem formulation is given, and the use of Linear Quadratic Gaussian (LQG) control law design and design sensitivity methods within the formulation is illustrated. Results of three simple integrated structure/control law design examples are presented. These results show the capability of structure and control law design tradeoffs to improve controlled system performance within the multilevel approach.
The modal correlation problem is formulated as a constrained optimization problem for validation of finite element models (FEM's). For large-scale structural applications, a pragmatic procedure for substructuring, model verification, and system integration is described to achieve effective modal correlation. The space station substructure FEM's are reduced using Lanczos vectors and integrated into a system FEM using Craig-Bampton component modal synthesis. The optimization code is interfaced with MSC/NASTRAN to solve the problem of modal test/analysis correlation; that is, the problem of validating FEM's for launch and on-orbit coupled loads analysis against experimentally observed frequencies and mode shapes. An iterative perturbation algorithm is derived and implemented to update nonlinear sensitivity (derivatives of eigenvalues and eigenvectors) during optimizer iterations, which reduced the number of finite element analyses.
The modal correlation problem is formulated as a constrained optimization problem for validation of finite element models (FEM's). For large-scale structural applications, a pragmatic procedure for substructuring, model verification, and system integration is described to achieve effective modal correlations. The space station substructure FEM's are reduced using Lanczos vectors and integrated into a system FEM using Craig-Bampton component modal synthesis. The optimization code is interfaced with MSC/NASTRAN to solve the problem of modal test/analysis correlation; that is, the problem of validating FEM's for launch and on-orbit coupled loads analysis against experimentally observed frequencies and mode shapes. An iterative perturbation algorithm is derived and implemented to update nonlinear sensitivity (derivatives of eigenvalues and eigenvectors) during optimizer iterations, which reduced the number of finite element analyses.
Two methods that exploit the availability of sensitivity derivatives are successfully employed to predict uncertainty propagation through Computational Fluid Dynamics (CFD) code for an inviscid airfoil problem. An approximate statistical second-moment method and a Sensitivity Derivative Enhanced Monte Carlo (SDEMC) method are successfully demonstrated on a two-dimensional problem. First- and second-order sensitivity derivatives of code output with respect to code input are obtained through an efficient incremental iterative approach. Given uncertainties in statistically independent, random, normally distributed flow parameters (input variables); these sensitivity derivatives enable one to formulate first- and second-order Taylor Series approximations for the mean and variance of CFD output quantities. Additionally, incorporation of the first-order sensitivity derivatives into the data reduction phase of a conventional Monte Carlo (MC) simulation allows for improved accuracy in determining the first moment of the CFD output. Both methods are compared to results generated using a conventional MC method. The methods that exploit the availability of sensitivity derivatives are found to be valid when considering small deviations from input mean values.
Inherent structural and functional properties of biotechnology-derived therapeutic biologics make them susceptible to light- and temperature-induced degradation and consequently can influence their quality. Photosensitivity of therapeutic proteins continues to be examined, but the commonalities and trends of storage conditions and information about light and temperature sensitivity among currently licensed therapeutic proteins has not been previously surveyed. Using a comprehensive and relational database approach, we conducted a scientific survey of all licensed biotechnology-derived drug products with the goal of providing evidence-based information about recommended storage conditions of formulations sorted by light- and temperature-related attributes as described for each product at licensure. We report the prevalence of indications for light and temperature sensitivity in formulations categorized by their presentation type, number of doses, container type, dosage form and active molecule type. We also report the storage temperature range across formulations and diluents for reconstitution and dilution. Formulations with excipients that potentially facilitate light-induced and thermal degradation were also noted. The result of our analysis indicates that light and temperature sensitivity are prevalent across therapeutic protein formulations. However, when a formulation is reconstituted or diluted, both light and temperature sensitivity are less clear. In addition, light and temperature sensitivity are more well defined in liquid formulations than lyophilized powder formulations, and more well defined in products manufactured in autoinjectors, prefilled-syringes, and pens than products in vials. Overall, our report provides a data-driven summary of storage conditions among therapeutic protein formulations to support the development of future biologic drug products.
Conathane EN-7 (referred to as EN-7) has been used for decades to pot electrical connectors, providing mechanical support for solder joints in cables. Unfortunately, the EN-7 formulation contains a suspect carcinogen and chemical sensitizer, toluene diisocyanate (TDI). Because of this, various groups have been formulating replacement materials, but all have come up short in final properties or in processing. We propose Arathane 5753 HVB as a replacement for EN-7. The properties compare very well with EN-7 and the processing has both advantages and disadvantages over EN-7 as discussed in detail below.
The finite-volume and finite-difference implementations of high-order accurate essentially nonoscillatory shock-capturing schemes are discussed and compared. Results obtained with fourth-order accurate algorithms based on both formulations are examined for accuracy, sensitivity to grid irregularities, resolution of waves that are oblique to the mesh, and computational efficiency. Some algorithm modifications that may be required for a given application are suggested. Conclusions that pertain to the relative merits of both formulations are drawn, and some circumstances for which each might be useful are noted.
The finite-volume and finite-difference implementations of high-order accurate essentially non-oscillatory shock-capturing schemes are discussed and compared. Results obtained with fourth-order accurate algorithms based on both formulations are examined for accuracy, sensitivity to grid irregularities, resolution of waves that are oblique to the mesh, and computational efficiency. Some algorithm modifications that may be required for a given application are suggested. Conclusions that pertain to the relative merits of both formulations are drawn, and some circumstances for which each might be useful are noted.
Quality assurance methods for graphite epoxy/prepregs were developed. Liquid chromatography, differential scanning calorimetry, and gel permeation chromatography were investigated. These methods were applied to a second prepreg system. The resin matrix formulation was correlated with mechanical properties. Dynamic mechanical analysis and fracture toughness methods were investigated. The chromatography and calorimetry techniques were all successfully developed as quality assurance methods for graphite epoxy prepregs. The liquid chromatography method was the most sensitive to changes in resin formulation. The were also successfully applied to the second prepreg system.
We present a new capability of the ice sheet model SICOPOLIS that enables flexible adjoint code generation via source transformation using the open-source algorithmic differentiation (AD) tool OpenAD.The adjoint code enables efficient calculation of the sensitivities of a scalar-valued objective function or quantity of interest (QoI) to a range of important, often spatially varying and uncertain model input variables, including initial and boundary conditions, as well as model parameters.Compared to earlier work on the adjoint code generation of SICOPOLIS, our work makes several important advances:(i) it is embedded within the up-to-date trunk of the SICOPOLIS repository – accounting for 1.5 decades of code development and improvements – and is readily available to the wider community;(ii) the AD tool used, OpenAD, is an open-source tool;(iii) the adjoint code developed is applicable to both Greenland and Antarctica, including grounded ice as well as floating ice shelves, with an extended choice of thermodynamical representations.A number of code refactorization steps were required. They are discussed in detail in an Appendix as they hold lessons for the application of AD to legacy codes at large.As an example application, we examine the sensitivity of the total Antarctic Ice Sheet volume to changes in initial ice thickness, austral summer precipitation, and basal and surface temperatures across the ice sheet.Simulations of Antarctica with floating ice shelves show that over 100 years of simulation the sensitivity of total ice sheet volume to the initial ice thickness and precipitation is almost uniformly positive, while the sensitivities to surface and basal temperature are almost uniformly negative. Sensitivity to austral summer precipitation is largest on floating ice shelves from Queen Maud to Queen Mary Land. The largest sensitivity to initial ice thickness is at outlet glaciers around Antarctica. Comparison between total ice sheet volume sensitivities to surface and basal temperature shows that surface temperature sensitivities are higher broadly across the floating ice shelves, while basal temperature sensitivities are highest at the grounding lines of floating ice shelves and outlet glaciers. A uniformly perturbed region of East Antarctica reveals that, among the four control variables tested here, total ice sheet volume is the most sensitive to variations in austral summer precipitation as formulated in SICOPOLIS.Comparison between adjoint- and finite-difference-derived sensitivities shows good agreement, lending confidence that the AD tool is producing correct adjoint code.The new modeling infrastructure is freely available at http://www.sicopolis.net (last access: 2 April 2020) under the development trunk.