Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Sensitivity formulations”

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 145 records · Page 8

Variation of solar cell sensitivity and solar radiation on tilted surfaces

The validity is studied that one of various insolation models used to compute solar radiation incident on tilted surfaces from global data measured on horizontal surfaces. The variation of solar cell sensitivity to solar radiation is determined over a wide range of atmospheric condition. A new model was formulated that reduced the deviations between measured and predicted insolation to less than 3 percent. Evaluation of solar cell sensitivity data indicates small change (2-3 percent) in sensitivity from winter to summer for tilted cells. The feasibility of using such global data as a means for calibrating terrestrial solar cells is discussed.

Klucher, T. M.↗

Computing Sensitivities in Evolutionary Systems: A Real-time Reduced Order Modeling Strategy

We present a new methodology for computing sensitivities in evolutionary systems using a model-driven low-rank approximation. To this end, we formulate a variational principle that seeks to minimize the distance between the time derivative of the reduced approximation and sensitivity dynamics. The first order optimality condition of the variational principle leads to a system of closed form evolution equations for an orthonormal basis and corresponding sensitivity coefficients. This approach allows for the computation of sensitivities with respect to a large number of parameters in an accurate and tractable manner by extracting correlations between different sensitivities on the fly. The presented method requires solving forward evolution equations, sidestepping the restrictions imposed by the forward/backward workflow of adjoint sensitivities. For example, the presented method, unlike the adjoint equation, does not impose any input/output load and can be used in applications in which real-time sensitivities are of interest. We demonstrate the utility of the method for three test cases: (1) computing sensitivity with respect to model parameters in the Rössler system, (2) computing sensitivity with respect to an infinite-dimensional forcing parameter in the chaotic Kuramoto--Sivashinsky equation, and (3) computing sensitivity with respect to reaction parameters for species transport in a turbulent reacting flow.

Reduced order model↗

Exploring the Science Trade Space with the JPL Innovation Foundry A-Team

The JPL Innovation Foundry has established a new approach for exploring, developing, and evaluating early concepts with a group called the Architecture Team (A-Team). The A-Team combines innovative collaborative methods and facilitated sessions with subject matter experts and analysis tools to help mature mission concepts. Science, implementation, and programmatic elements are all considered during an ATeam study. In these studies, Concept Maturity Levels (CML) are used to group methods. These levels include idea generation and capture (CML 1), initial feasibility assessment (CML 2), and trade space exploration (CML 3). Methods used for exploring the science objectives, feasibility, and scope will be described including use of a new technique for understanding the most compelling science, called a Science Return Diagram (SRD). In the process of developing the SRD, gradients in the science trade space are uncovered along with their implications for implementation and mission architecture. Special attention is paid towards developing complete investigations, establishing a series of logical claims that lead to the natural selection of a measurement approach. Over 20 science-focused A-Team studies have used these techniques to help science teams refine their mission objectives, make implementation decisions and reveal the mission concept’s most compelling science. This paper will describe the A-Team process for exploring the mission concept's science trade space and the Science Return Diagram technique.In June of 2011 a new collaborative engineering approach forearly concept formulation began in the JPL InnovationFoundry [1], six months later becoming the “A-Team” [2].Responding to a need for exploring mission architecturelevel trades [3], the A-Team precedes Team X [4,5] in asequence of concurrent engineering teams at JPL that can beused to mature a concept from a “cocktail napkin” level ideato a complete mission point design. The A-Team efficientlyexplores the science, implementation, and programmatictrade space in early concept formulation. Small, facilitatedgroups of experts generate innovative ideas, quantitativelyassess feasibility, and discover key sensitivities in the tradespace through collaborative analysis and use of advancedmethods and tools. The A-Team process builds off theexperience within JPL and other recent approaches to earlyconcept formulation [6] including best practices of the JPLInnovation Foundry, Project Systems Engineering &Formulation Section, Team Eureka and the Rapid MissionArchitecture Team[7].The A-Team is a focal point for innovative formulationapproaches and people within JPL. It relies on a largebackground of study resources, creative thinkers and “greybeard” scrutinizers, advanced tools, and subject matterexperts with both breadth and depth in experience andexpertise that are all available at JPL. The A-Team isdesigned to be a rapid and efficient process takingapproximately 6 weeks (the entire process can be as short asjust a few days or as long as up to three months) and costingthe equivalent of a work-month of a full-time employee orless. Studies begin with detailed planning and client reviewfollowed by study sessions, analysis work, and reporting.The staffing on each study is customized to the study goalsand objectives, and it is addressed early in the A-Teamprocess. Sessions are generally half-day or whole-day eventsand conducted over a series of days with focused agendas thatare moderated by a trained facilitator. Preliminary results andknowledge capture are available within hours of each session,and a final report is generally available two weeks later.One of the biggest challenges facing early conceptdevelopment is understanding the gradient in science returnversus various available mission scenarios and payload options. Often times, major areas of scientific inquiry havealready been prioritized by science groups, including throughthe National Research Council’s Decadal Studies inAstronomy, Planetary, and Earth Science. Yet science teamscontinue to struggle, especially in competitive missionsolicitations, to capture the right amount of scope that’sachievable within the cost constraints of the opportunity.Often the desire to completely and comprehensively study ascience area in just one mission (after all, true missionopportunities are rare) drives teams to take on too much,providing requirements that are unachievable within theresources of the opportunity without inducing unacceptableimplementation risk. Alternatively, science teams can seekto reduce risk by using an established instrument, but havenot thought through the traceability and key aspects of thescience question to justify its use. Both scenarios lead to badassumptions at the beginning of the concept development thatcan then ripple through implementation option choices,potentially preventing what would have been a good scienceinvestigation from being selected.The purpose of this paper is first to provide some additionalbackground and summary of the A-Team process, tools,people, and facilities. We then focus on the A-Teammethodology for overcoming the barriers of defining thescience scope well at the early concept development stage.This includes understanding the science story andtraceability, and then examining the gradient in science returnversus key characteristics of observables, developing theright payload and mission requirement specification throughexamining the science and implementation trade space.

Ziemer, John K.↗

Space shuttle post-entry and landing analysis. Volume 2: Appendices

Four candidate navigation systems for the space shuttle orbiter approach and landing phase are evaluated in detail. These include three conventional navaid systems and a single-station one-way Doppler system. In each case, a Kalman filter is assumed to be mechanized in the onboard computer, blending the navaid data with IMU and altimeter data. Filter state dimensions ranging from 6 to 24 are involved in the candidate systems. Comprehensive truth models with state dimensions ranging from 63 to 82 are formulated and used to generate detailed error budgets and sensitivity curves illustrating the effect of variations in the size of individual error sources on touchdown accuracy. The projected overall performance of each system is shown in the form of time histories of position and velocity error components.

Crawford, B. S.↗

Development of a model and computer code to describe solar grade silicon production processes

Mechanisms for the SiCl4/Na and SiF4/Na reaction systems were examined. Reaction schemes which include 25 elementary reactions were formulated for each system and run to test the sensitivity of the computed concentration and temperature profiles to the values given estimated rate coefficients. It was found that, for SiCl4/Na, the rate of production of free Si is largely mixing-limited for reasonable rate coefficient estimates. For the SiF4/Na system the results indicate that the endothermicities of many of the reactions involved in producing Si from SiF4/Na cause this system to be chemistry-limited rather than mixing-limited.

Gould, R. K.↗

Optimization by decomposition: A step from hierarchic to non-hierarchic systems

A new, non-hierarchic decomposition is formulated for system optimization that uses system analysis, system sensitivity analysis, temporary decoupled optimizations performed in the design subspaces corresponding to the disciplines and subsystems, and a coordination optimization concerned with the redistribution of responsibility for the constraint satisfaction and design trades among the disciplines and subsystems. The approach amounts to a variation of the well-known method of subspace optimization modified so that the analysis of the entire system is eliminated from the subspace optimization and the subspace optimizations may be performed concurrently.

Sobieszczanski-Sobieski, Jaroslaw↗

Optimization by decomposition: A step from hierarchic to non-hierarchic systems

A new, non-hierarchic decomposition is formulated for system optimization that uses system analysis, system sensitivity analysis, temporary decoupled optimizations performed in the design subspaces corresponding to the disciplines and subsystems, and a coordination optimization concerned with the redistribution of responsibility for the constraint satisfaction and design trades among the disciplines and subsystems, and a coordination optimization concerned with the redistribution of responsibility for the constraint satisfaction and design trades among the disciplines and subsystems. The approach amounts to a variation of the well-known method of subspace optimization modified so that the analysis of the entire system is eliminated from the subspace optimization and the subspace optimizations may be performed concurrently.

Sobieszczanski-Sobieski, Jaroslaw↗

Parallel-vector design sensitivity analysis in structural dynamics

This paper presents a parallel-vector algorithm for sensitivity calculations in linear structural dynamics. The proposed alternative formulation works efficiently with the reduced system of dynamic equations, since it eliminates the need for expensive and complicated based-vector derivatives, which are required in the conventional reduced system formulation. The relationship between the alternative formulation and the conventional reduced system formulation has been established, and it has been proven analytically that the two approaches are identical when all the mode shapes are included. This paper validates the proposed alternative algorithm through numerical experiments, where only a small number of mode shapes are used. In addition, a modified mode acceleration method is presented, thus not only the displacements but also the velocities and accelerations are shown to be improved.

Zhang, Y.↗

Remote sensing of vegetation canopy photosynthetic and stomatal conductance efficiencies

The problem of remote sensing the canopy photosynthetic and stomatal conductance efficiencies is investigated with the aid of one- and three-dimensional radiative transfer methods coupled to a semi-empirical mechanistic model of leaf photosynthesis and stomatal conductance. Desertlike vegetation is modeled as clumps of leaves randomly distributed on a bright dry soil with partial ground cover. Normalized difference vegetation index (NDVI), canopy photosynthetic (Ep), and stomatal efficiencies (Es) are calculated for various geometrical, optical, and illumination conditions. The contribution of various radiative fluxes to estimates of Ep is evaluated and the magnitude of errors in bulk canopy formulation of problem parameters are quantified. The nature and sensitivity of the relationship between Ep and Es to NDVI is investigated, and an algorithm is proposed for use in operational remote sensing.

Myneni, R. B.↗

Challenges in Modelling of Lightning-Induced Delamination; Effect of Temperature-Dependent Interfacial Properties

Lightning is a major cause of damage in laminated composite aerospace structures during flight. Due to the dielectric nature of Carbon fiber reinforced polymers (CFRPs), the high energy induced by lightning strike transforms into extreme, localized surface temperature accompanied with a high-pressure shockwave resulting in extensive damage. It is crucial to develop a numerical tool capable of predicting the damage induced from a lightning strike to supplement extremely expensive lightning experiments. Delamination is one of the most significant failure modes resulting from a lightning strike. It can be extended well beyond the visible damage zone, and requires sophisticated techniques and equipment to detect. A popular technique used to model delamination is the cohesive zone approach. Since the loading induced from a lightning strike event is assumed to consist of extreme localized heating, the cohesive zone formulation should additionally account for temperature effects. However, the sensitivity to this dependency remains unknown. Therefore, the major focus point of this work is to investigate the importance of this dependency via defining various temperature dependency profiles for the cohesive zone properties, and analyzing the corresponding delamination area. Thus, a detailed numerical model consisting of multidirectional composite plies with temperature-dependent cohesive elements in between is subjected to lightning (excessive amount of heat and pressure) and delamination/damage expansion is studied under specified conditions.

interlaminar damage↗

Predicting Char Yield of High-Temperature Resins

A simulation protocol is developed to predict the char yield of organic resins during high-temperature processing. Such in silico methods can help screen promising new formulations for advanced materials, but previously no chemistry-sensitive technique existed to predict the important experimental value of char yield. The method utilizes a reactive force field (ReaxFF) to model the chemical transformation of precursor monomers into carbonized structures during three processing stages: ramp-up to processing temperatures (~3000 K), pyrolysis, and quenching. Achieving good agreement with experimental char yields requires continuous removal of small byproduct molecules to mimic outgassing, and the application of high pressure to eliminate porosity and encourage graphitization. More than ten different resin chemistries are investigated, including arylacetylenes, cyanate esters, phthalonitriles, and polyimides, representing a diverse group of precursors with respect to initial cyclic content, heteroatoms and reactive groups. The protocol correctly predicts the relative char yield between the investigated chemistries and provides quantitative agreement with experimental values, especially for high char yield resins. The properties of the resins during processing are compared, including outgassing products, morphology of the final chemical configurations, cyclic content and mechanical properties.

polymer simulations, molecular dynamics↗

Predicting Char Yield of High-Temperature Resins

A simulation protocol is developed to predict the char yield of organic resins during high-temperature processing. Simulation protocol in silico methods can help screen promising new formulations for advanced materials, but previously no chemistry-sensitive technique existed to predict the important experimental value of char yield. The method utilizes a reactive force field (ReaxFF) to model the chemical transformation of precursor monomers into carbonized structures during three processing stages: ramp-up to processing temperatures (~3000 K), pyrolysis, and quenching. Achieving good agreement with experimental char yields requires continuous removal of small byproduct molecules to mimic outgassing, and the application of high pressure to eliminate porosity and encourage graphitization. More than ten different resin chemistries are investigated, including arylacetylenes, cyanate esters, phthalonitriles, and polyimides, representing a diverse group of precursors with respect to initial cyclic content, heteroatoms and reactive groups. The protocol correctly predicts the relative char yield between the investigated chemistries and provides quantitative agreement with experimental values, especially for high char yield resins. The properties of the resins during processing are compared, including outgassing products, morphology of the final chemical configurations, cyclic content and mechanical properties.

polymer simulations, molecular dynamics↗

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

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

Robertson, Franklin↗

Sensitivity Analysis of the Static Aeroelastic Response of a Wing

A technique to obtain the sensitivity of the static aeroelastic response of a three dimensional wing model is designed and implemented. The formulation is quite general and accepts any aerodynamic and structural analysis capability. A program to combine the discipline level, or local, sensitivities into global sensitivity derivatives is developed. A variety of representations of the wing pressure field are developed and tested to determine the most accurate and efficient scheme for representing the field outside of the aerodynamic code. Chebyshev polynomials are used to globally fit the pressure field. This approach had some difficulties in representing local variations in the field, so a variety of local interpolation polynomial pressure representations are also implemented. These panel based representations use a constant pressure value, a bilinearly interpolated value. or a biquadraticallv interpolated value. The interpolation polynomial approaches do an excellent job of reducing the numerical problems of the global approach for comparable computational effort. Regardless of the pressure representation used. sensitivity and response results with excellent accuracy have been produced for large integrated quantities such as wing tip deflection and trim angle of attack. The sensitivities of such things as individual generalized displacements have been found with fair accuracy. In general, accuracy is found to be proportional to the relative size of the derivatives to the quantity itself.

Eldred, Lloyd B.↗

On the sensitivity of complex, internally coupled systems

A method is presented for computing sensitivity derivatives with respect to independent (input) variables for complex, internally coupled systems, while avoiding the cost and inaccuracy of finite differencing performed on the entire system analysis. The method entails two alternative algorithms: the first is based on the classical implicit function theorem formulated on residuals of governing equations, and the second develops the system sensitivity equations in a new form using the partial (local) sensitivity derivatives of the output with respect to the input of each part of the system. A few application examples are presented to illustrate the discussion.

Sobieszczanskisobieski, Jaroslaw↗

Sensitivity of complex, internally coupled systems

A method is presented for computing sensitivity derivatives with respect to independent (input) variables for complex, internally coupled systems, while avoiding the cost and inaccuracy of finite differencing performed on the entire system analysis. The method entails two alternative algorithms: the first is based on the classical implicit function theorem formulated on residuals of governing equations, and the second develops the system sensitivity equations in a new form using the partial (local) sensitivity derivatives of the output with respect to the input of each part of the system. A few application examples are presented to illustrate the discussion.

Sobieszczanski-Sobieski, Jaroslaw↗

Sensitivity method for integrated structure/active control law design

The development is described of an integrated structure/active control law design methodology for aeroelastic aircraft applications. A short motivating introduction to aeroservoelasticity is given along with the need for integrated structures/controls design algorithms. Three alternative approaches to development of an integrated design method are briefly discussed with regards to complexity, coordination and tradeoff strategies, and the nature of the resulting solutions. This leads to the formulation of the proposed approach which is based on the concepts of sensitivity of optimum solutions and multi-level decompositions. The concept of sensitivity of optimum is explained in more detail and compared with traditional sensitivity concepts of classical control theory. The analytical sensitivity expressions for the solution of the linear, quadratic cost, Gaussian (LQG) control problem are summarized in terms of the linear regulator solution and the Kalman Filter solution. Numerical results for a state space aeroelastic model of the DAST ARW-II vehicle are given, showing the changes in aircraft responses to variations of a structural parameter, in this case first wing bending natural frequency.

Gilbert, Michael G.↗

Sensitivity analysis of a wing aeroelastic response

A variation of Sobieski's Global Sensitivity Equations (GSE) approach is implemented to obtain the sensitivity of the static aeroelastic response of a three-dimensional wing model. The formulation is quite general and accepts any aerodynamics and structural analysis capability. An interface code is written to convert one analysis's output to the other's input, and visa versa. Local sensitivity derivatives are calculated by either analytic methods or finite difference techniques. A program to combine the local sensitivities, such as the sensitivity of the stiffness matrix or the aerodynamic kernel matrix, into global sensitivity derivatives is developed. The aerodynamic analysis package FAST, using a lifting surface theory, and a structural package, ELAPS, implementing Giles' equivalent plate model are used.

Kapania, Rakesh K.↗