Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “lower and upper parameter bounds”

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 73 records · Page 4

Distributionally Robust Partially Observable Markov Decision Process with Moment-Based Ambiguity

In this paper, we consider a distributionally robust partially observable Markov decision process (DR-POMDP), where the distribution of the transition-observation probabilities is unknown at the beginning of each decision period, but their realizations can be inferred using side information at the end of each period after an action being taken. We build an ambiguity set of the joint distribution using bounded moments via conic constraints and seek an optimal policy to maximize the worst-case (minimum) reward for any distribution in the set. We show that the value function of DR-POMDP is piecewise linear convex with respect to the belief state and propose a heuristic search value iteration method for obtaining lower and upper bounds of the value function. We conduct numerical studies and demonstrate the computational performance of our approach via testing instances of a dynamic epidemic control problem. Our results show that DR-POMDP can produce more robust policies under misspecified distributions of transition-observation probabilities as compared to POMDP but has less costly solutions than robust POMDP. The DR-POMDP policies are also insensitive to varying parameter in the ambiguity set and to noise added to the true transition-observation probability values obtained at the end of each decision period.

97 MATHEMATICS AND COMPUTING↗

Structural tailoring of laminate properties

In this study, structural synthesis techniques have been applied to tailor the composite laminate properties from the micromechanics to the laminate level for a desirable structural response. It is also demonstrated how to establish upper and lower bounds on design variables based on reliability considerations. The question of reliability arises in fiber composite analysis and design because of the inherent scatter that is observed in the constituent material properties in the experimental data. The symmetric composite laminates subject to mechanical loading are considered and no hygrothermal effects are taken into account in the present research. Application examples illustrate the effects of various ply and laminate parameters on the overall properties and structural response of composite laminates.

Thanedar, P. B.↗

Discontinuous Aligned Carbon Fiber Intermediates for Automotive and Related Applications

This work focused on preferentially aligning discontinuous carbon fibers in wet-laid or air-laid processes. It is well known that aligned fibers provides higher directional strength and stiffness. Discontinuous fibers further allow higher degree of draw and formability as the gaps in the fibers allow for higher material movement. The current processes are limited in their ability to align carbon fibers during processing. The aligned fibers have several benefits - (a) in applications where chopped fibers can replace continuous fibers for targeted strength and stiffness metrics, but at a substantially reduced cost; (b) they can tolerate deeper draws than continuous fiber composites in thermo-stamping and compression molding processes; (c) they can be tailored for pultrusion and unidirectional applications. Although pultrusion is primarily a process that adopts continuous fibers, stitch bonded entangled discontinuous fibers can provide unique intermediates. This is analogous to natural coir fibers which get aligned and entangled to produce ropes/rods for example, (d) they can be processed in cross-ply and multi-directional formats, like composite laminates. In this work Neenah Paper partnered with IACMI, UT and ORNL to evaluate structure-process-property relationships with Zoltek carbon fiber. A few process parameters such as machine speed, weight basis, fiber length, effect of fiber sizing, direction of mat lay-up etc. were investigated. The produced mats were converted to thermoplastic composite laminates using polyamide 6 (PA6, nylon) resin. The specific objective of this project is to produce a wet-laid nonwoven carbon fiber mat with a high degree of unidirectional fiber alignment, using discontinuous carbon fibers. The report provides details about the processing, characterization, and lower-upper bound properties.

36 MATERIALS SCIENCE↗

Structural Analysis and Optimization of a Composite Fan Blade for Future Aircraft Engine

This report addresses the structural analysis and optimization of a composite fan blade sized for a large aircraft engine. An existing baseline solid metallic fan blade was used as a starting point to develop a hybrid honeycomb sandwich construction with a polymer matrix composite face sheet and honeycomb aluminum core replacing the original baseline solid metallic fan model made of titanium. The focus of this work is to design the sandwich composite blade with the optimum number of plies for the face sheet that will withstand the combined pressure and centrifugal loads while the constraints are satisfied and the baseline aerodynamic and geometric parameters are maintained. To satisfy the requirements, a sandwich construction for the blade is proposed with composite face sheets and a weak core made of honeycomb aluminum material. For aerodynamic considerations, the thickness of the core is optimized whereas the overall blade thickness is held fixed so as to not alter the original airfoil geometry. Weight is taken as the objective function to be minimized by varying the core thickness of the blade within specified upper and lower bounds. Constraints are imposed on radial displacement limitations and ply failure strength. From the optimum design, the minimum number of plies, which will not fail, is back-calculated. The ply lay-up of the blade is adjusted from the calculated number of plies and final structural analysis is performed. Analyses were carried out by utilizing the OpenMDAO Framework, developed at NASA Glenn Research Center combining optimization with structural assessment.

Coroneos, Rula M.↗

Two-dimensional radiative equilibrium of a gray medium in a plane layer bounded by gray nonisothermal walls

Radiative equilibrium temperature and surface heat flux distributions are calculated for an absorbing-emitting gray medium in an infinite plane layer bounded by gray diffuse walls with arbitrary temperature distributions. Superposition is used to obtain the solution for the differential approximation, which yields good accuracy for the optically thick medium. To also obtain accurate results for optically thin and intermediate regimes, the differential approximation is subsequently improved by a number of geometrical parameters, which are derived from the exact expression for the intensity. As an example, the case of constant temperature at the upper wall and a temperature step at the lower wall without heat generation in the medium is presented. Comparison with other available results shows excellent agreement.

Modest, M. F.↗

Estimation of Permanent Magnet Synchronous Motor Parameters using a Linear Unknown Input Interval Observer

This paper develops and evaluate a technique to estimate the speed and unknown load torque of a permanent magnet synchronous motor (PMSM). An interval unknown input observer (UIO) for linear time-invariant (LTI) systems is designed and utilized to estimate the PMSM parameters. Interval UIOs are advantageous as they allow for the estimation of both states and unknown inputs, while accounting for varying uncertainty sources, such as model and measurement uncertainties inherent in dynamic systems and processes. First, the nonlinear PMSM model is linearized to transform it to a suitable form for application of the linear interval UIO. Then, the interval UIO is applied to jointly estimate the motor speed and the unknown load torque disturbance. Assuming that the measurement noise and disturbances are bounded, lower and upper bounds are first computed for the unmeasured state (motor speed) and then for the unknown input (load torque). The proposed approach and its limitations are demonstrated for the nonlinear PMSM model derived from its equivalent electrical circuit.

Elinirina I Robinson↗

Continued research on selected parameters to minimize community annoyance from airport noise

A mathematical model of the annoyance created at an airport by aircraft operations is developed. The model incorporates population distribution considerations around an airport and the annoyance caused by aircraft noise. The objective function of this model corresponds to seeking to minimize total population annoyance created by all aircraft operations in a 24-hour period. Several factors are included in this model as constraint bounded. Demand for flight services is incorporated by including lower bounds on the number of operations by type of aircraft, runway used and time period. Also upper bounds on the number of operations for each runway are included. The mathematical model as formulated is recognized as corresponding to a nonlinear integer mathematical programming problem.

Frair, L.↗

Gear optimization

The use of formal numerical optimization methods for the design of gears is investigated. To achieve this, computer codes were developed for the analysis of spur gears and spiral bevel gears. These codes calculate the life, dynamic load, bending strength, surface durability, gear weight and size, and various geometric parameters. It is necessary to calculate all such important responses because they all represent competing requirements in the design process. The codes developed here were written in subroutine form and coupled to the COPES/ADS general purpose optimization program. This code allows the user to define the optimization problem at the time of program execution. Typical design variables include face width, number of teeth and diametral pitch. The user is free to choose any calculated response as the design objective to minimize or maximize and may impose lower and upper bounds on any calculated responses. Typical examples include life maximization with limits on dynamic load, stress, weight, etc. or minimization of weight subject to limits on life, dynamic load, etc. The research codes were written in modular form for easy expansion and so that they could be combined to create a multiple reduction optimization capability in future.

Vanderplaats, G. N.↗

Environmental Influences on Deep Convective Upscale Growth Rate in Central Argentina From a Convection‐Permitting Simulation

This study uses a convection‐permitting model simulation to describe the environmental conditions under which convective upscale growth occurs in central Argentina, particularly examining environmental parameters when deep convection initially forms that could differentiate the rate of initial upscale growth. Simulated mesoscale convective systems (MCSs) are separated into slow and rapid growth by the rate of spatial growth from convection initiation until reaching the MCS scale. A low‐level jet (LLJ) is found more frequently near the deep convection that experiences rapid growth to an MCS, but its presence alone is not predictive of rapid growth. Using spatially‐averaged parameters, we find that rapid growth to MCSs also occurs in environments that are significantly more thermodynamically favorable with greater low‐level moisture and instability. Fewer significant differences are found in the kinematic environment with only the 0–2 km vertical wind shear magnitude being significantly larger for rapid growth MCSs compared to slow growth MCSs, potentially related to LLJs often peaking near this height. When focusing only on MCSs with the slowest and fastest growth rates, elevated‐layer shear is significantly smaller for very rapid growth MCSs, suggesting elevated‐layer shear may help discriminate between the upper and lower bounds of growth rate. Finally, when upscale growth occurs near the Sierras de Córdoba (SDC) with a LLJ present, rapid growth is also supported by favorable wind shear orientation. However, this does not hold for upscale growth occurring away from the SDC, highlighting the importance of interpreting shear direction relative to the orientation of features initiating deep convection.

Sasaki, Clayton R. S. [Univ. of Washington, Seattl↗

Comparing Consider-Covariance Analysis with Sigma-Point Consider Filter and Linear-Theory Consider Filter Formulations

Recent literature in applied estimation theory reflects growing interest in the sigma-point (also called unscented ) formulation for optimal sequential state estimation, often describing performance comparisons with extended Kalman filters as applied to specific dynamical problems [c.f. 1, 2, 3]. Favorable attributes of sigma-point filters are described as including a lower expected error for nonlinear even non-differentiable dynamical systems, and a straightforward formulation not requiring derivation or implementation of any partial derivative Jacobian matrices. These attributes are particularly attractive, e.g. in terms of enabling simplified code architecture and streamlined testing, in the formulation of estimators for nonlinear spaceflight mechanics systems, such as filter software onboard deep-space robotic spacecraft. As presented in [4], the Sigma-Point Consider Filter (SPCF) algorithm extends the sigma-point filter algorithm to the problem of consider covariance analysis. Considering parameters in a dynamical system, while estimating its state, provides an upper bound on the estimated state covariance, which is viewed as a conservative approach to designing estimators for problems of general guidance, navigation and control. This is because, whether a parameter in the system model is observable or not, error in the knowledge of the value of a non-estimated parameter will increase the actual uncertainty of the estimated state of the system beyond the level formally indicated by the covariance of an estimator that neglects errors or uncertainty in that parameter. The equations for SPCF covariance evolution are obtained in a fashion similar to the derivation approach taken with standard (i.e. linearized or extended) consider parameterized Kalman filters (c.f. [5]). While in [4] the SPCF and linear-theory consider filter (LTCF) were applied to an illustrative linear dynamics/linear measurement problem, in the present work examines the SPCF as applied to nonlinear sequential consider covariance analysis, i.e. in the presence of nonlinear dynamics and nonlinear measurements. A simple SPCF for orbit determination, exemplifying an algorithm hosted in the guidance, navigation and control (GN&C) computer processor of a hypothetical robotic spacecraft, was implemented, and compared with an identically-parameterized (standard) extended, consider-parameterized Kalman filter. The onboard filtering scenario examined is a hypothetical spacecraft orbit about a small natural body with imperfectly-known mass. The formulations, relative complexities, and performances of the filters are compared and discussed.

Lisano, Michael E.↗

Search for photons above 10 18 eV by simultaneously measuring the atmospheric depth and the muon content of air showers at the Pierre Auger Observatory

The Pierre Auger Observatory is the most sensitive instrument to detect photons with energies above 1 0 17 eV . It measures extensive air showers generated by ultrahigh energy cosmic rays using a hybrid technique that exploits the combination of a fluorescence detector with a ground array of particle detectors. The signatures of a photon-induced air shower are a larger atmospheric depth of the shower maximum ( X max ) and a steeper lateral distribution function, along with a lower number of muons with respect to the bulk of hadron-induced cascades. In this work, a new analysis technique in the energy interval between 1 and 30 EeV ( 1 EeV = 1 0 18 eV ) has been developed by combining the fluorescence detector-based measurement of X max with the specific features of the surface detector signal through a parameter related to the air shower muon content, derived from the universality of the air shower development. No evidence of a statistically significant signal due to photon primaries was found using data collected in about 12 years of operation. Thus, upper bounds to the integral photon flux have been set using a detailed calculation of the detector exposure, in combination with a data-driven background estimation. The derived 95% confidence level upper limits are 0.0403, 0.01113, 0.0035, 0.0023, and 0.0021 km − 2 sr − 1 yr − 1 above 1, 2, 3, 5, and 10 EeV, respectively, leading to the most stringent upper limits on the photon flux in the EeV range. Compared with past results, the upper limits were improved by about 40% for the lowest energy threshold and by a factor 3 above 3 EeV, where no candidates were found and the expected background is negligible. The presented limits can be used to probe the assumptions on chemical composition of ultrahigh energy cosmic rays and allow for the constraint of the mass and lifetime phase space of super-heavy dark matter particles. Published by the American Physical Society 2024

79 ASTRONOMY AND ASTROPHYSICS↗

Deformation and Life Analysis of Composite Flywheel Disk and Multi-disk Systems

In this study an attempt is made to put into perspective the problem of a rotating disk, be it a single disk or a number of concentric disks forming a unit. An analytical model capable of performing an elastic stress analysis for single/multiple, annular/solid, anisotropic/isotropic disk systems, subjected to both pressure surface tractions, body forces (in the form of temperature-changes and rotation fields) and interfacial misfits is derived and discussed. Results of an extensive parametric study are presented to clearly define the key design variables and their associated influence. In general the important parameters were identified as misfit, mean radius, thickness, material property and/or load gradation, and speed; all of which must be simultaneously optimized to achieve the "best" and most reliable design. Also, the important issue of defining proper performance/merit indices (based on the specific stored energy), in the presence of multiaxiality and material anisotropy is addressed. These merit indices are then utilized to discuss the difference between flywheels made from PMC and TMC materials with either an annular or solid geometry. Finally two major aspects of failure analysis, that is the static and cyclic limit (burst) speeds are addressed. In the case of static limit loads, upper, lower, and out-of-plane bounds for disks with constant thickness are presented for both the case of internal pressure loading (as one would see in a hydroburst test) and pure rotation (as in the case of a free spinning disk). The results (interaction diagrams) are displayed graphically in designer friendly format. For the case of fatigue, a representative fatigue/life master curve is illustrated in which the normalized limit speed versus number of applied cycles is given for a cladded TMC disk application.

Arnold, S. M.↗

PMSM Parameter Estimation using a Linear Unknown Input Interval Observer

In this paper, the problem of permanent magnet synchronous motor (PMSM) speed and unknown load torque estimation is addressed. For this purpose, a interval unknown input observer (UIO) for linear time-invariant (LTI) systems is used. First, the PMSM model is linearized in order to make it in a suitable form for the linear interval UIO. Then, the interval UIO is applied to allow the joint estimation of the motor speed and the unknown load torque disturbance. The main advantages of this approach is that it not only allows the joint state and unknown input estimation, but also to take the different uncertainty sources into account. Indeed, taking model and measurement uncertainty into account is crucial. Assuming that the measurement noise and disturbances are bounded, lower and upper bounds are first computed for the unmeasured state (motor speed) and then for the unknown input (load torque). The proposed approach and its limitations are demonstrated with the nonlinear PMSM model derived from its equivalent electrical circuit.

Elinirina I Robinson↗

Mixing and Transition Control Studied

Considerable progress in understanding nonlinear phenomena in both unbounded and wallbounded shear flow transition has been made through the use of a combination of high- Reynolds-number asymptotic and numerical methods. The objective of this continuing work is to fully understand the nonlinear dynamics so that ultimately (1) an effective means of mixing and transition control can be developed and (2) the source terms in the aeroacoustic noise problem can be modeled more accurately. Two important aspects of the work are that (1) the disturbances evolve from strictly linear instability waves on weakly nonparallel mean flows so that the proper upstream conditions are applied in the nonlinear or wave-interaction streamwise region and (2) the asymptotic formulations lead to parabolic problems so that the question of proper out-flow boundary conditions--still a research issue for direct numerical simulations of convectively unstable shear flows--does not arise. Composite expansion techniques are used to obtain solutions that account for both mean-flow-evolution and nonlinear effects. A previously derived theory for the amplitude evolution of a two-dimensional instability wave in an incompressible mixing layer (which is in quantitative agreement with available experimental data for the first nonlinear saturation stage for a plane-jet shear layer, a circular-jet shear layer, and a mixing layer behind a splitter plate) have been extended to include a wave-interaction stage with a three-dimensional subharmonic. The ultimate wave interaction effects can either give rise to explosive growth or an equilibrium solution, both of which are intimately associated with the nonlinear self-interaction of the three dimensional component. The extended theory is being evaluated numerically. In contrast to the mixing-layer situation, earlier comparisons of theoretical predictions based on asymptotic methods and experiments in wall-bounded shear-flow transition have been somewhat lacking in one aspect or another. The current work strongly suggests that the main weakness is the underlying asymptotic representation of the linear "part" of the problem and not the explicit modeling of the nonlinear/wave-interaction effects. Consequently, the long-wave-length/high-Reynolds-number asymptotic limit for the Blasius boundary-layer stability problem was reexamined, and a new dispersion relationship for the instability waves that is uniformly valid for both the upper- and lower branch regions to the required order of approximation was obtained. A comparison with numerical results, obtained by solving the Orr-Sommerfeld stability problem, shows that the asymptotic formula provides surprisingly good results, even for values of the frequency parameter usually encountered in experimental investigations. This is particularly evident in the dynamically important upper-branch region, where much of the nonlinear interactions in transition experiments are believed to take place. The result is important in that it can be used to greatly improve the accuracy of weakly nonlinear critical-layer-based theories, and a consistent nonlinear theory is currently under evaluation.

Source record↗

Provably accurate simulation of gauge theories and bosonic systems

Quantum many-body systems involving bosonic modes or gauge fields have infinite-dimensional local Hilbert spaces which must be truncated to perform simulations of real-time dynamics on classical or quantum computers. To analyze the truncation error, we develop methods for bounding the rate of growth of local quantum numbers such as the occupation number of a mode at a lattice site, or the electric field at a lattice link. Our approach applies to various models of bosons interacting with spins or fermions, and also to both abelian and non-abelian gauge theories. We show that if states in these models are truncated by imposing an upper limit Λ on each local quantum number, and if the initial state has low local quantum numbers, then an error at most ϵ can be achieved by choosing Λ to scale polylogarithmically with ϵ − 1 , an exponential improvement over previous bounds based on energy conservation. For the Hubbard-Holstein model, we numerically compute a bound on Λ that achieves accuracy ϵ , obtaining significantly improved estimates in various parameter regimes. We also establish a criterion for truncating the Hamiltonian with a provable guarantee on the accuracy of time evolution. Building on that result, we formulate quantum algorithms for dynamical simulation of lattice gauge theories and of models with bosonic modes; the gate complexity depends almost linearly on spacetime volume in the former case, and almost quadratically on time in the latter case. We establish a lower bound showing that there are systems involving bosons for which this quadratic scaling with time cannot be improved. By applying our result on the truncation error in time evolution, we also prove that spectrally isolated energy eigenstates can be approximated with accuracy ϵ by truncating local quantum numbers at Λ = polylog ( ϵ − 1 ) .

Tong, Yu↗

Venusian Applications of 3D Convection Modeling

This study models mantle convection on Venus using the 'cubed sphere' code OEDIPUS, which models one-sixth of the planet in spherical geometry. We are attempting to balance internal heating, bottom mantle viscosity, and temperature difference across Venus' mantle, in order to create a realistic model that matches with current planetary observations. We also have begun to run both lower and upper mantle simulations to determine whether layered (as opposed to whole-mantle) convection might produce more efficient heat transfer, as well as to model coronae formation in the upper mantle. Upper mantle simulations are completed using OEDIPUS' Cartesian counterpart, JOCASTA. This summer's central question has been how to define a mantle plume. Traditionally, we have defined a hot plume the region with temperature at or above 40% of the difference between the maximum and horizontally averaged temperature, and a cold plume as the region with 40% of the difference between the minimum and average temperature. For less viscous cases (1020 Pa?s), the plumes generated by that definition lacked vigor, displaying buoyancies 1/100th of those found in previous, higher viscosity simulations (1021 Pa?s). As the mantle plumes with large buoyancy flux are most likely to produce topographic uplift and volcanism, the low viscosity cases' plumes may not produce observable deformation. In an effort to eliminate the smallest plumes, we experimented with different lower bound parameters and temperature percentages.

Venus↗

Probing High-Temperature Amine Chemistry: Is the Reaction NH 3 + NH 2 ⇌ N 2 H 3 + H 2 Important?

The reaction NH 3 + NH 2 ⇌ N 2 H 3 + H 2 (R1) has been identified as a key step to explain experimental results for pyrolysis and oxidation of ammonia. However, no direct experimental or theoretical evidence for the reaction has been reported. In the present work, the reaction was studied by ab initio theory and by reinterpretation of experimental data. We could not locate a transition state for R1 occurring as a direct process, but alternative mechanisms yielded an upper bound to k 1 of 1.5 × 10 13 exp(–58.9 kcal mol –1 /RT) cm 3 mol –1 s –1 over 1000–2500 K, several orders of magnitude lower than values applied in modeling. Consistent with the theoretical work, re-evaluation of NH 3 pyrolysis data supported a very low value of k 1 . However, this finding opens up a novel unresolved issue. Current kinetic models cannot capture the NH 3 oxidation behavior in a number of laminar flow reactor and jet-stirred reactor experiments without adopting an improbably high value for k 1 . Furthermore, important oxidation steps might be underestimated or missing from mechanisms.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Evidence of chaotic pattern in solar flux through a reproducible sequence of period-doubling-type bifurcations

Presented here is a preliminary study of the limits to solar flux intensity prediction, and of whether the general lack of predictability in the solar flux arises from the nonlinear chaotic nature of the Sun's physical activity. Statistical analysis of a chaotic signal can extract only its most gross features, and detailed physical models fail, since even the simplest equations of motion for a nonlinear system can exhibit chaotic behavior. A recent theory by Feigenbaum suggests that nonlinear systems that can be led into chaotic behavior through a sequence of period-doubling bifurcations will exhibit a universal behavior. As the control parameter is increased, the bifurcation points occur in such a way that a proper ratio of these will approach the universal Feigenbaum number. Experimental evidence supporting the applicability of the Feigenbaum scenario to solar flux data is sparse. However, given the hypothesis that the Sun's convection zones are similar to a Rayleigh-Bernard mechanism, we can learn a great deal from the remarkable agreement observed between the prediction by theory (period doubling - a universal route to chaos) and the amplitude decrease of the signal's regular subharmonics. The authors show that period-doubling-type bifurcation is a possible route to a chaotic pattern of solar flux that is distinguishable from the logarithm of its power spectral density. This conclusion is the first positive step toward a reformulation of solar flux by a nonlinear chaotic approach. The ultimate goal of this research is to be able to predict an estimate of the upper and lower bounds for solar flux within its predictable zones. Naturally, it is an important task to identify the time horizons beyond which predictability becomes incompatible with computability.

Ashrafi, S.↗