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 253 records · Page 14

Transient thermoelasticity and other body force effects in boundary element shape sensitivity analysis

Structural shape design sensitivity analysis (DSA) is presented for objects subjected to gravity, centrifugal, and general thermal loading. DSA employs implicit differentiation of the governing boundary integral equations. Both a surface integral and a particular integral approach are discussed for computation of these effects in response analysis and sensitivity analysis. It is shown that the particular integral approach is capable of computing the sensitivities of objects subjected to nonharmonic temperature distributions. Attention is also given to several techniques for the modeling of temperature distributions in the DSA process. In particular, an approach that eliminates the need for the incorporation of thermal response analysis and sensitivity analysis is described in the overall shape optimization process. A sharp corner formulation is also presented that can account for the jump in the surface traction components at a sharp corner in an interzone boundary.

Kane, James H.↗

Atmospheric Effects of Subsonic Aircraft: Interim Assessment Report of the Advanced Subsonic Technology Program

This first interim assessment of the subsonic assessment (SASS) project attempts to summarize concisely the status of our knowledge concerning the impacts of present and future subsonic aircraft fleets. It also highlights the major areas of scientific uncertainty, through review of existing data bases and model-based sensitivity studies. In view of the need for substantial improvements in both model formulations and experimental databases, this interim assessment cannot provide confident numerical predictions of aviation impacts. However, a number of quantitative estimates are presented, which provide some guidance to policy makers.

Friedl, Randall R.↗

Regime-based TRMM and GV Microphysical Studies at MSFC and UAH

Differences in rain rate between TMI and PR vary systematically with PR Z-profile statistics, whose frequency of occurrence is modified to create seasonal biases in the sub-tropical Southeastern U.S. (and almost certainly elsewhere). Tropical (non-tropical) DSDs in N. Alabama exhibit larger (smaller) D(sub 0), and larger (smaller) N(sub 0) and mu. The formulation process for empirical retrievals of DSD using dual-pol radar is sensitive to D(sub max) assumptions used in the scattering model stage. 4. DSD retrievals from Parsivel disdrometers compare favorably to those of the 2DVD unless rain rates exceed 25 mm/hr and D(sub m) exceeds 2 mm (at which point the Parsivels overestimate D(sub m) and rain rate).

Petersen, Walter A.↗

Optimization of shallow arches against instability using sensitivity derivatives

The author discusses the problem of optimization of shallow frame structures which involve a coupling of axial and bending responses. A shallow arch of a given shape and of given weight is optimized such that its limit point load is maximized. The cross-sectional area, A(x) and the moment of inertia, I(x) of the arch obey the relationship I(x) = rho A(x) sup n, n = 1,2,3 and rho is a specified constant. Analysis of the arch for its limit point calculation involves a geometric nonlinear analysis which is performed using a corotational formulation. The optimization is carried out using a second-order projected Lagrangian algorithm and the sensitivity derivatives of the critical load parameter with respect to the areas of the finite elements of the arch are calculated using implicit differentation. Results are presented for an arch of a specified rise to span ratio under two different loadings and the limitations of the approach for the intermediate rise arches are addressed.

Kamat, Manohar P.↗

Sensitivity of rotor blade vibration characteristics to torsional oscillations.

A theoretical investigation of dynamic response characteristics of helicopter rotor blades in forward flight was carried out with special emphasis on the torsional degrees-of-freedom. The finite element method was applied in the formulation of the coupled equations of motion for flapwise bending and torsion for blades with non-collinear elastic, mass and aerodynamic axes. The sensitivity of blade vibration characteristics with respect to structural, geometric and aerodynamic properties as well as flight conditions was evaluated. Numerical results for sample blades were plotted to show the variation of the coupling between bending and torsional components of the response.

Bratanow, T.↗

Model form and sensitivity analysis of CALPHAD-based nucleation models in b-stabilized Ti alloys

Accurate prediction of α-phase nucleation and growth in β-stabilized titanium alloys is crucial for designing heat treatments to optimize mechanical properties in additively manufactured lightweight components. Ideally, predictions of nucleation and growth would incorporate both top-down observations of past experimental heat treatments and bottom-up modeling of phase transformations; however, the appropriate method of combining these information sources is not self-evident. Combining top-down and bottom-up information requires a unified form of model that can connect between spatiotemporal scales, as well as sets of fitting parameters that can be identified by each data source. The selection of which parameters to fit to which data source can be made based on expert opinion, or by performing a sensitivity analysis. In solid-solid nucleation, direct observation of the nucleation and growth process is challenging. Most data on the heat treatment-controlled phase transformations are not in-situ. To predict the process and outcome of the nucleation, growth and coarsening of precipitates, theoretical models of the nucleation pathway are used to bridge the gap. Many sources of uncertainty affect the modeling of this nucleation process. It can be influenced by small variations in the thermomechanical processing history, chemical composition, and initial microstructure. If molecular dynamics (MD) simulations are used to determine thermodynamic quantities and inform CALPHAD modeling, additional uncertainty can be introduced and accounted for using Bayesian methods. Top-down uncertainties require additional steps to quantify. The influence of nucleation model form on the sensitivity of predictions to input parameters and physical conditions is the focus of this study. Classical nucleation theory (CNT) allows modeling to formulate the nucleation as homogeneous or, more commonly, heterogeneous. Non-classical nucleation models are also increasingly explored as a means of reconciling top-down and bottom-up data. In this study, the sensitivity of the intragranular nucleation of α in a β-annealed, slow-cooled aging (BASCA) heat treatment of β-stabilized Ti5553 alloy is explored using CNT and both heterogeneous and homogeneous assumptions. The Kampmann-Wagner Numerical model of precipitate nucleation and growth is employed. Using open-source tools (pyCalphad and thermodynamic modeling of TiMo as a surrogate system, a sensitivity analysis is performed to measure variations in key parameters, including chemical driving force, interfacial energy, and diffusivity, as they relate to predictions of precipitate number density. The inclusion of top-down and bottom-up data in selection of nucleation model form is discussed.

Rodriguez Negron, A. M.↗

Numerical analysis of soot emissions from gasoline-ethanol and gasoline-butanol blends under gasoline compression ignition conditions

In the present work, computational fluid dynamics (CFD) simulations of a single-cylinder gasoline compression ignition (GCI) engine were performed to investigate the impact of blending two biofuels, ethanol and n-butanol, with gasoline on the trade-off between combustion phasing and soot emissions under low load conditions. Here, in order to represent market gasoline (RD5-87), a four-component toluene primary reference fuel (TPRF) + ethanol (ETPRF) surrogate (with 20% ethanol by mole; E20) was formulated using a neural network based octane predictor such that the surrogate had the same ethanol content, Research Octane Number (RON) and Octane Sensitivity (S). In addition, a novel skeletal kinetic mechanism for ETPRF and TPRF + n-butanol (BTPRF) blends, incorporating polycyclic aromatic hydrocarbon (PAH) chemistry, was developed. A three-dimensional (3D) engine CFD formulation employing the skeletal mechanism, adaptive mesh refinement (AMR), finite-rate chemistry approach, and hybrid method of moments (HMOM) was adopted to capture the in-cylinder combustion phenomena and soot emissions. The engine CFD model was validated against RD5-87 experimental data for a broad range of start-of-injection (SOI) timings (-21/-27/-36/-45 crank angle degrees (CAD) after top-dead center (aTDC)), with respect to in-cylinder pressure, heat release rate, combustion phasing, and soot emissions. The closed-cycle simulation results were analyzed to elucidate the non-monotonic trend of soot emissions versus SOI timing: SOI-36 > SOI-45 > SOI-21 > SOI-27. Thereafter, the validated CFD model was employed to simulate the combustion of a gasoline-ethanol blend with 45% (by mole) ethanol (E45) and a gasoline-butanol blend with 45% (by mole) n-butanol (B45) under the same operating conditions to study the effects of fuel composition and SOI timing on combustion phasing and soot emissions. The sooting propensity followed the trend: B45 > E20 > E45 at all SOI timings. Overall, it was observed that the autoignition propensity was primarily related to fuel chemistry. On the other hand, sooting propensity showed strong coupling with both fuel chemistry and physical properties, with greater impact of fuel physical properties at advanced SOI timings.

30 DIRECT ENERGY CONVERSION↗

Controllability of Free-piston Stirling Engine/linear Alternator Driving a Dynamic Load

This paper presents the dynamic behavior of a Free-Piston Stirling Engine/linear alternator (FPSE/LA) driving a single-phase fractional horse-power induction motor. The controllability and dynamic stability of the system are discussed by means of sensitivity effects of variations in system parameters, engine controller, operating conditions, and mechanical loading on the induction motor. The approach used expands on a combined mechanical and thermodynamic formulation employed in a previous paper. The application of state-space technique and frequency domain analysis enhances understanding of the dynamic interactions. Engine-alternator parametric sensitivity studies, similar to those of the previous paper, are summarized. Detailed discussions are provided for parametric variations which relate to the engine controller and system operating conditions. The results suggest that the controllability of a FPSE-based power system is enhanced by proper operating conditions and built-in controls.

Kankam, M. David↗

Deep Reinforcement Learning Enabled Physical-Model-Free Two-Timescale Voltage Control Method for Active Distribution Systems

Active distribution networks are being challenged by frequent and rapid voltage violations due to renewable energy integration. Conventional model-based voltage control methods rely on accurate parameters of the distribution networks, which are difficult to achieve in practice. This paper proposes a novel physical-model-free two-timescale voltage control framework for active distribution systems. To achieve fast control of PV inverters, the whole network is first partitioned into several subnetworks using voltage-reactive power sensitivity. Then, the scheduling of PV inverters in the multiple sub-networks is formulated as Markov games and solved by a multi-agent soft actor-critic (MASAC) algorithm, where each subnetwork is modeled as an intelligent agent. All agents are trained in a centralized manner to learn a coordinated strategy while being executed based on only local information for fast response. For the slower time-scale control, OLTCs and switched capacitors are coordinated by a single agent-based SAC algorithm using the global information with considering control behaviors of the inverters. Particularly, the two-level agents are trained concurrently with information exchange according to the reward signal calculated from the data-driven surrogate model. Comparative tests with different benchmark methods on IEEE 33-and 123-bus systems and 342-node low voltage distribution system demonstrate that the proposed method can effectively mitigate the fast voltage violations and achieve systematical coordination of different voltage regulation assets without the knowledge of accurate system model.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Adaptive Recovery Model: Designing Systems for Testing Tracing and Vaccination to Support COVID-19 Recovery Planning.

This report documents a new approach to designing disease control policies that allocate scarce testing, contact tracing, and vaccination resources to better control community transmission of COVID19 or similar diseases. The Adaptive Recovery Model (ARM) combines a deterministic compartmental disease model with a stochastic network disease propagation model to enable us to simulate COVID-19 community spread through the lens of two complementary modeling motifs. ARM contact networks are derived from cell-phone location data that have been anonymized and interpreted as individual arrivals to specic public locations. Modeling disease spread over these networks allows us to identify locations within communities conducive to rapid disease spread. ARM applies this model- and data-derived abstractions of community transmission to evaluate the effectiveness of disease control measures including targeted social distancing, contact tracing, testing and vaccination. The architecture of ARM provides a unique capacity to help decision makers understand how best to deploy scarce testing, tracing and vaccination resources to minimize disease-spread potential in a community. This document details the novel mathematical formulations underlying ARM, presents a dynamical stability analysis of the deterministic model components, a sensitivity analysis of control parameters and network structure, and summarizes a process for deriving contact networks from cell-phone location data. An example use case steps through applying ARM to evaluate three targeted social distancing policies using Bernalillo County, New Mexico as an exemplar test locale. This step-by-step analysis demonstrates how ARM can be used to measure the relative performance of competing public health policies. Initial scenario tests of ARM shows that ARMs design focus on resource utilization rather than simple incidence prediction can provide decision makers with additional quantitative guidance for managing ongoing public health emergencies and planning future responses.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

A calculation of the Rosseland mean opacity of dust grains in primordial solar system nebulae

The Rosseland mean opacity owing to grains was calculated as a function of temperature and density for nebulae having solar elemental abundances. The values of the mean opacity were evaluated with a generalized formulation allowing for anisotropic scattering. The values of the mean opacity do not depend sensitively on the choice of the particle size distribution function, provided that there are few particles having sizes in excess of several tens of microns. The results indicate that thermal convection in primordial nebulae occurs over broader ranges of altitudes at low temperatures than at high temperatures, and for size distributions for which extensive aggregation has not yet occurred.

Pollack, J. B.↗

Numerical parametric studies of spray combustion instability

A coupled numerical algorithm has been developed for studies of combustion instabilities in spray-driven liquid rocket engines. The model couples gas and liquid phase physics using the method of fractional steps. Also introduced is a novel, efficient methodology for accounting for spray formation through direct solution of liquid phase equations. Preliminary parametric studies show marked sensitivity of spray penetration and geometry to droplet diameter, considerations of liquid core, and acoustic interactions. Less sensitivity was shown to the combustion model type although more rigorous (multi-step) formulations may be needed for the differences to become apparent.

Pindera, M. Z.↗

A Flight Dynamics Model for a Multi-Actuated Flexible Rocket Vehicle

A comprehensive set of motion equations for a multi-actuated flight vehicle is presented. The dynamics are derived from a vector approach that generalizes the classical linear perturbation equations for flexible launch vehicles into a coupled three-dimensional model. The effects of nozzle and aerosurface inertial coupling, sloshing propellant, and elasticity are incorporated without restrictions on the position, orientation, or number of model elements. The present formulation is well suited to matrix implementation for large-scale linear stability and sensitivity analysis and is also shown to be extensible to nonlinear time-domain simulation through the application of a special form of Lagrange s equations in quasi-coordinates. The model is validated through frequency-domain response comparison with a high-fidelity planar implementation.

Orr, Jeb S.↗

Adjoint-Based Minimization of X-59 Sonic Boom Noise Via Control Surfaces

A multidisciplinary design optimization methodology to directly minimize the ground-level noise generated by a supersonic aircraft's sonic boom is presented. A gradient-based optimizer is coupled with a Cartesian Euler flow solver and an atmospheric propagation tool to forge a unique design capability. Adjoint formulations for both the flow solver and propagation tool are also coupled to provide sensitivities in an efficient manner. The design method is demonstrated on the X-59 Low-Boom Flight Demonstrator by optimizing control surface deflections to improve ground-level noise while maintaining trimmed cruise flight. Two noise minimization examples and a noise maximization are presented. All optimized designs result in measurable improvement in the ground-level noise objective. Extensive surveys of the design space confirm that the optimization method is effective in finding a local optimum. Moreover, repeated application of the method with varying initial design points also demonstrates the robustness of the method.

ARMD↗

Sonic Boom Ground Noise Minimization Via the Adjoint Method

A multidisciplinary design optimization methodology to directly minimize the ground-level noise generated by the sonic boom of a high altitude supersonic body is presented. A Cartesian Euler flow solver is coupled with an atmospheric propagation tool to create a ground-noise analysis capability for supersonic bodies. Adjoint formulations for both the flow solver and propagation tool are also coupled to compute sensitivities of shape variations in the body in a highly efficient manner. A gradient-based optimizer is then introduced to forge a valuable design capability. Output-based mesh adaptation that is driven directly by ground-level noise is employed to increase accuracy and provide error estimation. The design method is demonstrated first on a simple axisymmetric body with few design variables to evaluate the efficacy of the optimization scheme. Guided by the results of this initial case, the problem is then repeated with somewhat different design variables to further demonstrate the capabilities of the design method. Finally, the method is applied to a real-world problem by optimizing control surface deflection settings of a low-boom aircraft to minimize ground noise while maintaining trimmed, level flight.

ARMD↗

Variation of solar cell sensitivity and solar radiation on tilted surfaces

An empirical study was performed (1) to evaluate the validity of various insolation models used to compute solar radiation incident on tilted surfaces from global data measured on horizontal surfaces and (2) to determine the variation of solar cell sensitivity to solar radiation over a wide range of atmospheric condition. Evaluation of the insolation data indicates that the isotropic sky model of Liu and Jordan underestimates the amount of solar radiation falling on tilted surfaces by as much as 10%. An anisotropic-clear-sky model proposed by Temps and Coulson was also evaluated and found to be deficient under cloudy conditions. A new model, formulated herein, reduced the deviations between measured and predicted insolation to less than 3%. Evaluation of solar cell sensitivity data indicates small change (2-3%) in sensitivity from winter to summer for tilted cells. The feasibility of using such global data as a means for calibrating terrestrial solar cells as done by Treble is discussed.

Klucher, T. M.↗

A New Window into the Baryon Cycle at Cosmic Noon with Line Intensity Mapping: Forecasts for auto- and cross-correlations in [CII]-158$μ$m, HI 21 cm, CO$_{J+1\rightarrow J}$, and H$α$ galaxies

Across the peak of cosmic star formation at $z\sim1-2$, inflow, processing, and feedback drive rapid changes in the spatial distribution and chemical composition of baryons in galaxies and surrounding reservoirs; this baryon cycle can be tomographically mapped by line intensity mapping (LIM) of atomic hydrogen, ionized carbon, and carbon monoxide. We present a simulation-based forecasting framework for detecting auto- and cross-power spectra between spectroscopic surveys of four such tracers at $z\sim0.5-1.7$ mapping the same deep field - TIM, EoRSpec/FYST, MeerKAT, & Euclid. We forward-model 3-D distributions for these tracers from magnetohydrodynamic simulations, directly capturing the two-halo, one-halo, and shot statistics without relying on analytical decompositions. We further detail a signal-to-noise formalism, tailored to LIM surveys with highly anisotropic geometries and Fourier-space coverage. We demonstrate that galaxy cross-correlations will be the dominant discovery channel for current-generation surveys. These instruments will detect the auto-spectra for CO and HI 21 cm and the CO $\times$ 21 cm cross-spectrum at modest S/N $\sim 1-10$, while placing upper limits on the [CII]-158$μ$m signals. [CII], CO, and HI LIM will be $\sim3-30\times$ ($0.5-1.5$ dex) more sensitive to cross-correlation with the Euclid survey, however, than their respective auto-correlations, constraining all three models of line emission at high significance (S/N $\sim 10-40$) within this decade. Finally, we formulate a staged instrumental trajectory with planned or reasonable improvements, including the as-proposed SKA-Mid. We forecast advancing the per-$k$-mode sensitivities of each auto-, galaxy-line, and line-line spectrum by several orders of magnitude, enabling new percent- and sub-percent level constraints on cosmology and the redshift evolution of star formation and the baryon cycle.

Agrawal, Shubh [Pennsylvania U., Dept. Math.]↗

An investigation of using an RQP based method to calculate parameter sensitivity derivatives

Estimation of the sensitivity of problem functions with respect to problem variables forms the basis for many of our modern day algorithms for engineering optimization. The most common application of problem sensitivities has been in the calculation of objective function and constraint partial derivatives for determining search directions and optimality conditions. A second form of sensitivity analysis, parameter sensitivity, has also become an important topic in recent years. By parameter sensitivity, researchers refer to the estimation of changes in the modeling functions and current design point due to small changes in the fixed parameters of the formulation. Methods for calculating these derivatives have been proposed by several authors (Armacost and Fiacco 1974, Sobieski et al 1981, Schmit and Chang 1984, and Vanderplaats and Yoshida 1985). Two drawbacks to estimating parameter sensitivities by current methods have been: (1) the need for second order information about the Lagrangian at the current point, and (2) the estimates assume no change in the active set of constraints. The first of these two problems is addressed here and a new algorithm is proposed that does not require explicit calculation of second order information.

Beltracchi, Todd J.↗