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.

SourceResearch connectionFrom yearThrough year

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

At least 613 records · Page 34

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Co-Optimization Scheme for Hybrid Electric Vehicles Powertrain and Exhaust Emission Control System Using Future Speed Prediction

Hybrid electric vehicles (HEVs) have been an effective solution for improved vehicle fuel efficiency and reduced emission pollution. In this paper, a co-optimization scheme is proposed to optimize fuel efficiency for HEVs. Here, the proposed optimization scheme uses obtainable future speed prediction as the basis to optimally tune control parameters for the existing powertrain control system. The ramp-up time of the catalyst temperature to reach its light-off level in the exhaust emission system is also considered as an additional optimization constraint to reduce emission. The Toyota Prius Hybrid Simulink model which is an integrated model for a powertrain and exhaust emission system is validated using real data from several real driving cycle scenarios. Then, to simplify the formulation of the proposed algorithm, the model for the optimization for powertrain and exhaust emission systems is represented by a set of equivalent neural network (NN) models learned using the data generated from the well-validated Toyota Prius Hybrid Simulink model. Using NN models, a co-optimization algorithm is established that provides an optimal tuning of some fuel-sensitive powertrain control parameters using future speed prediction, leading to a novel co-optimization algorithm, achieving on average a further 9.22% fuel savings for the Toyota Prius Hybrid Simulink model.

33 ADVANCED PROPULSION SYSTEMS↗

DOE OSTI2022

Design and execution of a verification, validation, and uncertainty quantification plan for a numerical model of left ventricular flow after LVAD implantation

Left ventricular assist devices (LVADs) are implantable pumps that act as a life support therapy for patients with severe heart failure. Despite improving the survival rate, LVAD therapy can carry major complications. Particularly, the flow distortion introduced by the LVAD in the left ventricle (LV) may induce thrombus formation. While previous works have used numerical models to study the impact of multiple variables in the intra-LV stagnation regions, a comprehensive validation analysis has never been executed. The main goal of this work is to present a model of the LV-LVAD system and to design and follow a verification, validation and uncertainty quantification (VVUQ) plan based on the ASME V&V40 and V&V20 standards to ensure credible predictions. The experiment used to validate the simulation is the SDSU cardiac simulator, a bench mock-up of the cardiovascular system that allows mimicking multiple operation conditions for the heart-LVAD system. The numerical model is based on Alya, the BSC’s in-house platform for numerical modelling. Alya solves the Navier-Stokes equation with an Arbitrary Lagrangian-Eulerian (ALE) formulation in a deformable ventricle and includes pressure-driven valves, a 0D Windkessel model for the arterial output and a LVAD boundary condition modeled through a dynamic pressure-flow performance curve. The designed VVUQ plan involves: (a) a risk analysis and the associated credibility goals; (b) a verification stage to ensure correctness in the numerical solution procedure; (c) a sensitivity analysis to quantify the impact of the inputs on the four quantities of interest (QoIs) (average aortic root flow $Q^{avg}_{Ao}$, maximum aortic root flow $Q^{max}_{Ao}$, average LVAD flow $Q^{avg}_{VAD}$, and maximum LVAD flow $Q^{max}_{VAD}$; (d) an uncertainty quantification using six validation experiments that include extreme operating conditions. Numerical code verification tests ensured correctness of the solution procedure and numerical calculation verification showed a grid convergence index (GCI)95% <3.3%. The total Sobol indices obtained during the sensitivity analysis demonstrated that the ejection fraction, the heart rate, and the pump performance curve coefficients are the most impactful inputs for the analysed QoIs. The Minkowski norm is used as validation metric for the uncertainty quantification. It shows that the midpoint cases have more accurate results when compared to the extreme cases. The total computational cost of the simulations was above 100 [core-years] executed in around three weeks time span in Marenostrum IV supercomputer. This work details a novel numerical model for the LV-LVAD system, that is supported by the design and execution of a VVUQ plan created following recognised international standards. We present a methodology demonstrating that stringent VVUQ according to ASME standards is feasible but computationally expensive.

59 BASIC BIOLOGICAL SCIENCES↗

NASA NTRS1989

Error analysis for semi-analytic displacement derivatives with respect to shape and sizing variables

Sensitivity analysis is fundamental to the solution of structural optimization problems. Consequently, much research has focused on the efficient computation of static displacement derivatives. As originally developed, these methods relied on analytical representations for the derivatives of the structural stiffness matrix (K) with respect to the design variables (b sub i). To extend these methods for use with complex finite element formulations and facilitate their implementation into structural optimization programs using the general finite element method analysis codes, the semi-analytic method was developed. In this method the matrix the derivative of K/the derivative b sub i is approximated by finite difference. Although it is well known that the accuracy of the semi-analytic method is dependent on the finite difference parameter, recent work has suggested that more fundamental inaccuracies exist in the method when used for shape optimization. Another study has argued qualitatively that these errors are related to nonuniform errors in the stiffness matrix derivatives. The accuracy of the semi-analytic method is investigated. A general framework was developed for the error analysis and then it is shown analytically that the errors in the method are entirely accounted for by errors in delta K/delta b sub i. Furthermore, it is demonstrated that acceptable accuracy in the derivatives can be obtained through careful selection of the finite difference parameter.

Fenyes, Peter A.↗

NASA NTRS1990

Free vibrations of thin-walled semicircular graphite-epoxy composite frames

A detailed study is made of the effects of variations in lamination and material parameters of thin walled composite frames on their vibrational characteristics. The structures considered are semicircular thin walled frames with I and J sections. The flanges and webs of the frames are modelled by using 2-D shell and plate finite elements. A mixed formulation is used with the fundamental unknowns consisting of both the generalized displacements and stress resultants in the frames. The frequencies and modes predicted by the 2-D finite element model are compared with those obtained from experiments, as well as with the predictions of a 1-D thin walled beam finite element model. A detailed study is made of the sensitivity of the vibrational response to variations in the fiber orientation, material properties of the individual layers, and boundary conditions.

Noor, Ahmed K.↗

NASA NTRS1992

Underlying modal data issues for detecting damage in truss structures

The modal strain energy distribution and frequency and mode shape sensitivities to structural stiffness changes are studied in order to overcome the current limits of modal data limit measurements in damage location. It is established that damage detection is feasible for members significantly contributing to the strain energy of the measured modes; the modes that are most effective in detecting the damage of critical members are identifiable, and a relationship is formulated between the accuracy of the measured modes and frequencies and the probability of damage detection.

Kashangaki, Thomas A.-L.↗

DOE OSTI2024

Simultaneous shape and topology optimization of inflatable soft robots

Simultaneous shape and topology optimization is used to design pressure-activated inflatable soft robots. The pressure loaded boundary is meshed conformingly and shape optimized, while the morphology of the robot is topology optimized. The design objective is to exert maximum force on an object, i.e. to produce soft “grippers”. The robot’s motion is modeled using nearly incompressible finite deformation hyperelasticity. To ensure stability of the robot, the buckling load factors obtained via linearized buckling analyses are constrained. The finite element method is used to evaluate the optimization cost and constraint functions and the adjoint method is employed to compute their sensitivities. The numerical examples produce pressure-driven soft robots with varying complexity. We also compare our simultaneous optimization results to those obtained via sequential topology and then shape optimization.

42 ENGINEERING↗

DOE OSTI2020

Analysis of the performance of different culvert boundary conditions in 2D shallow flow models

Culverts allow roads to safely traverse small streams or drainage ditches, and their proper design is critical to ensure a safe and reliable transportation network. A correct modelization of these hydraulic structures becomes crucial in the assessment of flood footprints or discharge peak estimation in a risk evaluation plan. The question of how to include culverts comes up frequently when assembling a hydraulic model that requires the presence of as many singular elements as possible. In this work, three different culvert integrations with the surface domain are studied and compared in the context of a 2D shallow water (SW) model. All of them are based on the Federal Highway Administration (FHWA) formulation for the culvert discharge estimation but differ in complexity and in the interaction with the numerical model for surface flow, some of them as internal boundary conditions. Several steady and unsteady validation test cases are presented and the numerical results are compared with the predictions from HEC-RAS 1D and HY-8 software. Furthermore, the culvert area, shape and their sensitivity to the 2D computational mesh is also analyzed.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Negative Index Metamaterial Lens for Subwavelength Microwave Detection

Metamaterials are engineered periodic structures designed to have unique properties not encountered in naturally occurring materials. One such unusual property of metamaterials is the ability to exhibit negative refractive index over a prescribed range of frequencies. A lens made of negative refractive index metamaterials can achieve resolution beyond the diffraction limit. This paper presents the design of a metamaterial lens and its use in far-field microwave imaging for subwavelength defect detection in nondestructive evaluation (NDE). Theoretical formulation and numerical studies of the metamaterial lens design are presented followed by experimental demonstration and characterization of metamaterial behavior. Finally, a microwave homodyne receiver-based system is used in conjunction with the metamaterial lens to develop a far-field microwave NDE sensor system. A subwavelength focal spot of size 0.82λ was obtained. The system is shown to be sensitive to a defect of size 0.17λ × 0.06λ in a Teflon sample. Consecutive positions of the defect with a separation of 0.23λ was resolvable using the proposed system.

47 OTHER INSTRUMENTATION↗

NASA NTRS1990

Free vibrations of thin-walled semicircular graphite-epoxy composite frames

A detailed study is made of the effects of variations in lamination and material parameters of thin walled composite frames on their vibrational characteristics. The structures considered are semicircular thin walled frames with I and J sections. The flanges and webs of the frames are modeled by using 2-D shell and plate finite elements. A mixed formulation is used with the fundamental unknowns consisting of both the generalized displacements and stress resultants in the frame. The frequencies and modes predicted by the 2-D finite element model are compared with those obtained from experiments, as well as with the predictions of a non-dimensional thin walled beam finite element model. A detailed study is made of the sensitivity of the vibrational response to variations in the fiber orientation, material properties of the individual layers, and boundary conditions.

Carden, Huey D.↗

NASA NTRS1990

Free vibrations of thin-walled semicircular graphite-epoxy composite frames

A detailed study is made of the effects of variations in lamination and material parameters of thin walled composite frames on their vibrational characteristics. The structures considered are semicircular thin walled frames with I and J sections. The flanges and webs of the frames are modeled by using 2-D shell and plate finite elements. A mixed formulation is used with the fundamental unknowns consisting of both the generalized displacements and stress resultants in the frame. The frequencies and modes predicted by the 2-D finite element model are compared with those obtained from experiments, as well as with the predictions of a non-dimensional thin walled beam finite element model. A detailed study is made of the sensitivity of the vibrational response to variations in the fiber orientation, material properties of the individual layers, and boundary conditions.

Carden, Huey D.↗

NASA NTRS1993

The evolution equation for the flame surface density in turbulent premixed combustion

The mean reaction rate in flamelet models for turbulent premixed combustion depends on two basic quantities: a mean chemical rate, called the flamelet speed, and the flame surface density. Our previous work had been primarily focused on the problem of the structure and topology of turbulent premixed flames, and it was then determined that the flamelet speed, when space-averaged, is only weakly sensitive to the turbulent flow field. Consequently, the flame surface density is the key quantity that conveys most of the effects of the turbulence on the rate of energy release. In flamelet models, this quantity is obtained via a modeled transport equation called the Sigma-equation. Past theoretical work has produced a rigorous approach that leads to an exact but unclosed formulation for the turbulent Sigma-equation. In the exact Sigma-equation, it appears that the dynamical properties of the flame surface density are determined by a single parameter, namely the turbulent flame stretch. Unfortunately, the turbulent flame stretch as well as the flame surface density is not available from experiments, and, in the absence of experimental data, little is known on the validity of the closure assumptions used in current flamelet models. Direct Numerical Simulation (DNS) is the alternative approach to get basic information on these fundamental quantities. In the present work, three-dimensional DNS of premixed flames in isotropic turbulent flow is used to estimate the different terms appearing in the Sigma-equation. A new methodology is proposed to provide the source and sink terms for the flame surface density, resolved both temporally and spatially throughout the turbulent flame brush. Using this methodology, our objective is to extract the turbulent flame stretch from the DNS data base and then perform extensive comparisons with flamelet models. Thanks to the detailed information produced by the DNS-based analysis, it is expected that this type of comparison will not only underscore the shortcomings of current models, but also suggest ways to improve them.

Trouve, Arnaud↗

NASA NTRS1995

An a priori model for the reduction of nutation observations: KSV(1994.3) nutation series

We discuss the formulation of a new nutation series to be used in the reduction of modern space geodetic data. The motivation for developing such a series is to develop a nutation series that has smaller short period errors than the IAU 1980 nutation series and to provide a series that can be used with techniques such as the Global Positioning System (GPS) that have sensitivity to nutations but can directly separate the effects of nutations from errors in the dynamical force models that effect the satellite orbits. A modern nutation series should allow the errors in the force models for GPS to be better understood. The series is constructed by convolving the Kinoshita and Souchay rigid Earth nutation series with an Earth response function whose parameters are partly based on geophysical models of the Earth and partly estimated from a long series (1979-1993) of very long baseline interferometry (VLBI) estimates of nutation angles. Secular rates of change of the nutation angles to represent corrections to the precession constant and a secular change of the obliquity of the ecliptic are included in the theory. Time dependent amplitudes of the Free Core Nutation (FCN) that is most likely excited by variations in atmospheric pressure are included when the geophysical parameters are estimated. The complex components of the prograde annual nutation are estimated simultaneously with the geophysical parameters because of the large contribution to the nutation from the S(sub 1) atmospheric tide. The weighted root mean square (WRMS) scatter of the nutation angle estimates about this new model are 0.32 mas and the largest correction to the series when the amplitudes of the ten largest nutations are estimated is 0.18 +/- 0.03 mas for the in phase component of the prograde 18. 6 year nutation.

Herring, T. A.↗

DOE OSTI2021

AC Power Flow Based DLMP Calculation and Decomposition Method to Smooth Power Fluctuation of Distributed Renewable Energy Sources

As the penetration of renewable energy sources increases, the growing renewable power variability brings ramping issues to power systems. Meanwhile, the development of distributed energy resources (DERs) makes the distribution systems to provide both energy and ancillary services. To incentivise individual resources and customers to alleviate ramping issues on the demand side, a two-stage distribution locational marginal price (DLMP) calculation and decomposition method is developed to formulate the marginal power ramping price for DERs. In the first stage of the proposed method, a distribution system operator market scheduling model based on AC optimal power flow is designed to estimate the optimal operating point of the distribution system. Subsequently, the voltage and power flow constraints are linearised in stage two to calculate DLMP. Finally, based on the Lagrange function and sensitivity factors, DLMP is decomposed to the marginal costs for active/reactive power, voltage management, power loss and power variability. Case studies demonstrate that the proposed model can effectively smooth the power fluctuation and reduce the ramping flexibility requirements of distribution systems.

AC optimal power flow↗

NASA NTRS1985

Seasonal simulations of the planetary boundary layer and boundary-layer stratocumulus clouds with a general circulation model

The formulation of the planetary boundary layer (PBL) and stratocumulus parametrizations in the UCLA general circulation model (GCM) are briefly summarized, and extensive new results are presented illustrating some aspects of the simulated seasonal changes of the global distributions of PBL depth, stratocumulus cloudiness, cloud-top entrainment instability, the cumulus mass flux, and related fields. Results from three experiments designed to reveal the sensitivity of the GCM results to aspects of the PBL and stratocumulus parametrizations are presented. The GCM results show that the layer cloud instability appears to limit the extent of the marine subtropical stratocumulus regimes, and that instability frequently occurs in association with cumulus convection over land. Cumulus convection acts as a very significant sink of PBL mass throughout the tropics and over the midlatitude continents in winter.

Randall, D. A.↗

NASA NTRS2006

Gain-Scheduled Fault Tolerance Control Under False Identification

An active fault tolerant control (FTC) law is generally sensitive to false identification since the control gain is reconfigured for fault occurrence. In the conventional FTC law design procedure, dynamic variations due to false identification are not considered. In this paper, an FTC synthesis method is developed in order to consider possible variations of closed-loop dynamics under false identification into the control design procedure. An active FTC synthesis problem is formulated into an LMI optimization problem to minimize the upper bound of the induced-L2 norm which can represent the worst-case performance degradation due to false identification. The developed synthesis method is applied for control of the longitudinal motions of FASER (Free-flying Airplane for Subscale Experimental Research). The designed FTC law of the airplane is simulated for pitch angle command tracking under a false identification case.

Shin, Jong-Yeob↗

NASA NTRS2006

Development of a Test Facility for Air Revitalization Technology Evaluation

Lu, Sao-Dung↗

DOE OSTI2023

Exploring pressure-dependent inelastic deformation and failure in bonded granular composites: An energetic materials perspective

In polymer-filled granular composites, damage may develop in mechanical loading prior to material failure. Damage mechanisms such as microcracking or plastic deformation in the binder phase can substantially alter the material’s mesostructure. For energetic materials, such as solid propellants and plastic bonded explosives, these mesostructural changes can have far reaching effects including degraded mechanical properties, potentially increased sensitivity to further insults, and changes in expected performance. Unfortunately, predicting damage is nontrivial due to the complex nature of these composites and the entangled interactions between inelastic mechanisms. In this work, we assess the current literature of experimental knowledge, focusing on the pressure-dependent shear response, and propose a simple simulation framework of bonded particles to study four limiting-case material formulations at both meso- and macro-scales. Further, to construct the four cases, we systematically vary the relative interfacial strength between the polymer binder and granular filler phase and also vary the polymer’s glass transition temperature relative to operating temperature which determines how much the binder can plastically deform. These simulations identify key trends in global mechanical response, such as the emergence of strain hardening or softening regimes with increasing pressure which qualitatively resemble experimental results. By quantifying the activation of different inelastic mechanisms, such as bonds breaking and plastically straining, we identify when each mechanism becomes relevant and provide insight into potential origins for changes in mechanical responses. The locations of broken bonds are also used to define larger, mesoscopic cracks to test various metrics of damage. We primarily focus on triaxial compression, but also test the opposite case of triaxial extension to highlight the impact of Lode angle on mechanical behavior.

36 MATERIALS SCIENCE↗

NASA NTRS1980

Illustration of airfoil shape effect on forward-swept wing divergence

A static aeroelastic analysis is presented of the divergence of untapered wings with conventional and supercritical airfoil sections at sweep angles of zero and -15 deg. One bending and one torsion mode were employed for a uniform rectangular cantilevered beam with the elastic axis at midchord, and calculations were based on a two-dimensional differential equations formulation in the structural coordinate system and in simple strip theory. A minimum divergence speed in the transonic range is obtained which is associated with the rearward shift of the aerodynamic center, and a 17% difference in minimum divergence dynamic pressure is found between a supercritical and a conventional wing. It is noted that although the strip method employed allows the assessment of the sensitivity of airfoil shapes to divergence, three-dimensional transonic aerodynamic methods should be used to predict wing divergence characteristics.

Bland, S. R.↗