Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “linearized distribution”

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 415 records · Page 23

Birkeland currents in the plasma sheet

A search was conducted for the signatures of Birkeland currents in the Earth's magnetic tail, using observed values of B(sub x) and B(sub y) from large sets of spacecraft data. The data were binned by x and y for -10 greater than x(sub GSM) greater than -35 and absolute value of y(sub GSM) less than or equal to 20 R(sub E) (less than or equal to 30 R(sub E) for x(sub GSM) less than or equal to -25 R(sub E)) and in each bin their distribution in the (B(sub x), B(sub y)) plane was fitted by least squares to a piecewise linear function. That gave average x-y distributions of the flaring angle between B(sub xy) and the x direction, as well as that angle's variation across the thickness of the plasma sheet. Angles obtained in the central plasma sheet differed from those derived near the lobe boundary. That is the expected signature if earthward or tailward Birkeland current sheets are embedded in the plasma sheet, and from this dfiference we derived the dawn-dusk profiles of the tail Birkeland currents for several x(sub GSM) intervals. It was found that (1) the Birkeland currents have the sense of region 1 currents, when mapped to the ionosphere; (2) both the linear current density (kiloamperes/R(sub E)) and the net magnitude of the field-aligned currents decrease rapidly down the tail; (3) the total Birkeland current at x approximately equals -10 R(sub E) equals approximately equals 500-700 kA, which is approx. 30% of the net region 1 current observed at ionospheric altitudes, in agreement with model mapping results; and (4) the B(sub z) and B(sub y) components of the interplanetary magnetic field influence the distribution of Birkeland currents in the tail.

Tsyganenko, Nikolai A.↗

Joint Estimation of Topology and Injection Statistics in Distribution Grids with Missing Nodes

Optimal operation of distribution grid resources relies on accurate estimation of its state and topology. Practical estimation of such quantities is complicated by the limited presence of real-time meters. This article discusses a theoretical framework to jointly estimate the operational topology and statistics of injections in radial distribution grids under limited availability of nodal voltage measurements. In particular, we show that our proposed algorithms are able to provably learn the exact grid topology and injection statistics at all unobserved nodes as long as they are not adjacent. The algorithm design is based on novel ordered trends in voltage magnitude fluctuations at node groups, that are independently of interest for radial physical flow networks. The complexity of the designed algorithms is theoretically analyzed and their performance is validated using both linearized and nonlinear ac power flow samples in test distribution grids.

97 MATHEMATICS AND COMPUTING↗

Inverse sequential procedures for the monitoring of time series

When one or more new values are added to a developing time series, they change its descriptive parameters (mean, variance, trend, coherence). A 'change index (CI)' is developed as a quantitative indicator that the changed parameters remain compatible with the existing 'base' data. CI formulate are derived, in terms of normalized likelihood ratios, for small samples from Poisson, Gaussian, and Chi-Square distributions, and for regression coefficients measuring linear or exponential trends. A substantial parameter change creates a rapid or abrupt CI decrease which persists when the length of the bases is changed. Except for a special Gaussian case, the CI has no simple explicit regions for tests of hypotheses. However, its design ensures that the series sampled need not conform strictly to the distribution form assumed for the parameter estimates. The use of the CI is illustrated with both constructed and observed data samples, processed with a Fortran code 'Sequitor'.

Radok, Uwe↗

Taylor approximation variance reduction for approximation errors in PDE-constrained Bayesian inverse problems

In numerous applications, surrogate models are used as a replacement for accurate parameter-to-observable mappings when solving large-scale inverse problems governed by partial differential equations (PDEs). The surrogate model may be a computationally cheaper alternative to the accurate parameter-to-observable mappings and/or may ignore additional unknowns or sources of uncertainty. The Bayesian approximation error (BAE) approach provides a means to account for the induced uncertainties and approximation errors, i.e. the errors between the accurate parameter-to-observable mapping and the surrogate. The statistics of these errors are, however, in general unknown a priori, and are thus calculated using Monte Carlo sampling. Although the sampling is typically carried out offline, i.e. before considering the data, the process can still represent a computational bottleneck. In this work, we develop a scalable computational approach for reducing the costs associated with the sampling stage of the BAE approach. Specifically, we consider the Taylor expansion of the accurate and surrogate forward models with respect to the uncertain parameter fields either as a control variate for variance reduction or as a means to directly and efficiently approximate the mean and covariance of the approximation errors. We propose efficient methods for evaluating the expressions for the mean and covariance of the Taylor approximations based on linear(-ized) PDE solves. Furthermore, the proposed approach is independent of the dimension of the uncertain parameter, depending instead on the intrinsic dimension of the data, ensuring scalability to high-dimensional problems. The potential benefits of the proposed approach are demonstrated for two high-dimensional inverse problems governed by PDE examples, namely for the estimation of a distributed Robin boundary coefficient in a linear diffusion problem, and for a coefficient estimation problem governed by a nonlinear diffusion problem.

Bayesian approximation error↗

Information distribution in distributed microprocessor based flight control systems

This paper presents an optimal control theory that accounts for variable time intervals in the information distribution to control effectors in a distributed microprocessor based flight control system. The theory is developed using a linear process model for the aircraft dynamics and the information distribution process is modeled as a variable time increment process where, at the time that information is supplied to the control effectors, the control effectors know the time of the next information update only in a stochastic sense. An optimal control problem is formulated and solved that provides the control law that minimizes the expected value of a quadratic cost function. An example is presented where the theory is applied to the control of the longitudinal motions of the F8-DFBW aircraft. Theoretical and simulation results indicate that, for the example problem, the optimal cost obtained using a variable time increment Markov information update process where the control effectors know only the past information update intervals and the Markov transition mechanism is almost identical to that obtained using a known uniform information update interval.

Montgomery, R. C.↗

Nonlinear state estimation and feedback control of nonlinear and bilinear distributed parameter systems

This paper presents a theory of nonlinear state observers for nonlinear and bilinear distributed parameter systems. Convergence results are proved for these observers. Linear feedback control derived from such state observers is applied to the distributed parameter system and conditions are presented for closed-loop stability. The emphasis is on finite dimensional state observers and controllers (which can be implemented with on-line computers) and conditions for their successful operation with infinite dimensional distributed parameter systems.

Balas, M. J.↗

Combined linear theory/impact theory method for analysis and design of high speed configurations

Pressure distributions on a wing body at Mach 4.63 are calculated. The combined theory is shown to give improved predictions over either linear theory or impact theory alone. The combined theory is also applied in the inverse design mode to calculate optimum camber slopes at Mach 4.63. Comparisons with optimum camber slopes obtained from unmodified linear theory show large differences. Analysis of the results indicate that the combined theory correctly predicts the effect of thickness on the loading distributions at high Mach numbers, and that finite thickness wings optimized at high Mach numbers using unmodified linear theory will not achieve the minimum drag characteristics for which they are designed.

Brooke, D.↗

A Method of Determining Aerodynamic-Influence Coefficients from Wind-Tunnel Data for Wings at Supersonic Speeds

A method is described for determining aerodynamic-influence coefficients from wind-tunnel data for calculating the steady-state load distribution on a wing with arbitrary angle-of-attack distribution at supersonic speeds. The method combines linearized theory with empirical adjustments in order to give accurate results over a wide range of angles of attack. The experimented data required are pressure distributions measured on a flat wing of the desired planform at the desired Mach number and over the desired range of angles of attack. The method has been tested by applying it to wind-tunnel data measured at Mach numbers of 1.61 and 2.01 on wings of the same planform but of different surface shapes. Influence coefficients adjusted to fit the flat wing gave good predictions of the spanwise and chord-wise distributions of loadings measured on twisted and cambered wings.

Gainer, Patrick A.↗

Continuous analysis of stresses from arbitrary surface loads on a half space

A new form of elemental surface load on a half space is introduced, presuming a quasi-pyramidal variation of load which is doubly linear in each of four rectangular parts of a surface rectangle. Approximations of arbitrary load distributions by sums of such elements are continuous, piecewise linear in two directions and well adaptable. The loads may be normal or tangential. The explicit solutions obtained for all stress and displacement components due to each elemental load involve only elementary functions, are free of the discontinuities which arise with stepwise elements, and are suitable for computing. Some illustrative stress distributions are presented for elemental loads and for multiple pyramidal loads involving both normal and tangential loads. The value of the load continuity in the more complicated analyses of surface cracks is also illustrated.

Bell, J. C.↗

An analysis of the threshold method for measuring area-average rainfall

Experimental evidence shows that the area-average rain rate and the fractional area covered by rain rate exceeding a fixed threshold are highly correlated; that is, are highly linearly related. A precise theoretical explanation of this fact is given. The explanation is based on the observation that rain rate has a mixed distribution, one that is a mixture of a discrete distribution and a continuous distribution. Under a homogeneity assumption, the slope of the linear relationship depends only on the continuous part of the distribution and as such is found to be markedly immune to parameter changes. This is illustrated by certain slope surfaces obtained from three specific distributions. The threshold level can be chosen in an optimal way by minimizing a certain distance function defined over the threshold range. In general, the threshold level should be not too far from the mean rain rate conditional on rain. The so-called threshold method advocates measuring rainfall from fractional area exploiting the observed linear relationship of the later with the area average rain rate. The method is potentially useful for the estimation of rainfall from space via satellites.

Kedem, Benjamin↗

Optimum design of structures subject to general periodic loads

A simplified version of Icerman's problem regarding the design of structures subject to a single harmonic load is discussed. The nature of the restrictive conditions that must be placed on the design space in order to ensure an analytic optimum are discussed in detail. Icerman's problem is then extended to include multiple forcing functions with different driving frequencies. And the conditions that now must be placed upon the design space to ensure an analytic optimum are again discussed. An important finding is that all solutions to the optimality condition (analytic stationary design) are local optima, but the global optimum may well be non-analytic. The more general problem of distributing the fixed mass of a linear elastic structure subject to general periodic loads in order to minimize some measure of the steady state deflection is also considered. This response is explicitly expressed in terms of Green's functional and the abstract operators defining the structure. The optimality criterion is derived by differentiating the response with respect to the design parameters. The theory is applicable to finite element as well as distributed parameter models.

Reiss, Robert↗

Incompressible lifting-surface aerodynamics for a rotor-stator combination

Current literature on the three dimensional flow through compressor cascades deals with a row of rotor blades in isolation. Since the distance between the rotor and stator is usually 10 to 20 percent of the blade chord, the aerodynamic interference between them has to be considered for a proper evaluation of the aerothermodynamic performance of the stage. A unified approach to the aerodynamics of the incompressible flow through a stage is presented that uses the lifting surface theory for a compressor cascade of arbitrary camber and thickness distribution. The effects of rotor stator interference are represented as a linear function of the rotor and stator flows separately. The loading distribution on the rotor and stator flows separately. The loading distribution on the rotor and stator blades and the interference factor are determined concurrently through a matrix iteration process.

Ramachandra, S. M.↗

An experimental study of the validity of the heat-field concept for sonic-boom alleviation

An experimental program was carried out in the NASA-Langley 4 ft x 4 ft supersonic pressure tunnel to investigate the validity of the heat-field concept for sonic boom alleviation. The concept involves heating the flow about a supersonic aircraft in such a manner as to obtain an increase in effective aircraft length and yield an effective aircraft shape that will result in a shock-free pressure signature on the ground. First, a basic body-of-revolution representing an SST configuration with its lift equivalence in volume was tested to provide a baseline pressure signature. Second, a model having a 5/2-power area distribution which, according to theory, should yield a linear pressure rise with no front shock wave was tested. Third, the concept of providing the 5/2-power area distribution by using an off-axis slender fin below the basic body was investigated. Then a substantial portion (approximately 40 percent) of the solid fin was replaced by a heat field generated by passing heated nitrogen through the rear of the fin.

Swigart, R. J.↗

Enabling Efficient Sparse Computations using Linear Algebra Aware Compilers

This project developed the LAPIS compiler framework, built on the Multilevel Intermediate Representation (MLIR), to optimize sparse linear algebra operations and support performance portability across diverse architectures. The main innovation of LAPIS is the Kokkos dialect, which allows for lowering codes from a high productivity language to different architectures in an elegant way. The dialect also allows the conversion of lower-level MLIR code to C++ Kokkos code, facilitating the integration of scientific machine learning (SciML) models into applications. To extend LAPIS for distributed memory architectures, a new partition dialect was created to manage the distribution of sparse tensors and express communication patterns for sparse linear algebra operations. This dialect also supports the distributed execution of operators and includes algorithmic optimizations to minimize communication to improve performance. The project also demonstrates that MLIR can enable effective linear algebra-level optimizations, improving performance on different GPUs for both sparse and dense linear algebra kernels. Key applications of LAPIS include sparse linear algebra and graph kernels, TenSQL, a relational database management solution built on GraphBLAS, and the development of subgraph isomorphism and monomorphism kernels, showcasing performance portability. In summary, the LAPIS framework supports productivity, performance, portability, and distributed memory execution, while also enabling linear algebra-level optimizations that are challenging in traditional programming languages, with successful applications ranging from simple sparse linear algebra to complex graph kernels.

97 MATHEMATICS AND COMPUTING↗

A high-speed linear algebra library with automatic parallelism

Parallel or distributed processing is key to getting highest performance workstations. However, designing and implementing efficient parallel algorithms is difficult and error-prone. It is even more difficult to write code that is both portable to and efficient on many different computers. Finally, it is harder still to satisfy the above requirements and include the reliability and ease of use required of commercial software intended for use in a production environment. As a result, the application of parallel processing technology to commercial software has been extremely small even though there are numerous computationally demanding programs that would significantly benefit from application of parallel processing. This paper describes DSSLIB, which is a library of subroutines that perform many of the time-consuming computations in engineering and scientific software. DSSLIB combines the high efficiency and speed of parallel computation with a serial programming model that eliminates many undesirable side-effects of typical parallel code. The result is a simple way to incorporate the power of parallel processing into commercial software without compromising maintainability, reliability, or ease of use. This gives significant advantages over less powerful non-parallel entries in the market.

Boucher, Michael L.↗

Coordinated Optimization and Control of Residential Solid State Power Substations in Electrical Distribution Network

This paper proposes a computationally efficient mixed integer linear programming (MILP) model for the coordinated optimization of solid-state power substations (SSPSs) in a feeder considering the full unbalanced three-phase structure of the distribution grid and its characteristics. The proposed model determines the optimal real and reactive power of each SSPS at the point of common coupling (PCC) that minimizes the operating cost and maximizes the system performance, e.g., voltage regulation and phase balancing. To improve the computational efficiency, an inscribed octagon is introduced to approximate the quadratic capacity constraints of components. Numerical simulation results show the effectiveness of the proposed model and significant improvements in voltage profiles and power imbalance between phases.

Liu, Guodong↗

The ENSO Effect on the Temporal and Spatial Distribution of Global Lightning Activity

The recently reprocessed (1997-2006) OTD/LIS database is used to investigate the global lightning climatology in response to the ENSO cycle. A linear correlation map between lightning anomalies and ENSO (NINO3.4) identifies areas that generally follow patterns similar to precipitation anomalies. We also observed areas where significant lightning/ENSO correlations are found and are not accompanied of significant precipitation/ENSO correlations. An extreme case of the strong decoupling between lightning and precipitation is observed over the Indonesian peninsula (Sumatra) where positive lightning/NINO3.4 correlations are collocated with negative precipitation/NINO3.4 correlations. Evidence of linear relationships between the spatial extent of thunderstorm distribution and the respective NINO3.4 magnitude are presented for different regions on the Earth. Strong coupling is found over areas remote to the main ENSO axis of influence and both during warm and cold ENSO phases. Most of the resulted relationships agree with the tendencies of precipitation related to ENSO empirical maps or documented teleconnection patterns. Over the Australian continent, opposite behavior in terms of thunderstorm activity is noted for warm ENSO phases with NINO3.4 magnitudes with NINO3.4>+l.08 and 0<NqNO3.4<I.08. Finally, we investigate the spatial distribution of areas that consistently portrayed enhanced lightning activity during the main warm/cold (El Nino/La Nina) ENSO episodes of the past decade. The observed patterns show no spatial overlapping and identify areas that in their majority are in agreement with empirical precipitation/ENSO maps. The areas that appear during the warm ENSO phase are found over regions that have been identified as anomalous Hadley circulation ENSO-related patterns. The areas that appear during the cold ENSO phase are found predominantly around the west hemisphere equatorial belt and are in their majority identified by anomalous Walker circulation.

Chronis, Themis G.↗

Structural Verification and Modeling of a Tension Cone Inflatable Aerodynamic Decelerator

Verification analyses were conducted on membrane structures pertaining to a tension cone inflatable aerodynamic decelerator using the analysis code LS-DYNA. The responses of three structures - a cylinder, torus, and tension shell - were compared against linear theory for various loading cases. Stress distribution, buckling behavior, and wrinkling behavior were investigated. In general, agreement between theory and LS-DYNA was very good for all cases investigated. These verification cases exposed the important effects of using a linear elastic liner in membrane structures under compression. Finally, a tension cone wind tunnel test article is modeled in LS-DYNA for which preliminary results are presented. Unlike data from supersonic wind tunnel testing, the segmented tension shell and torus experienced oscillatory behavior when subjected to a steady aerodynamic pressure distribution. This work is presented as a work in progress towards development of a fluid-structures interaction mechanism to investigate aeroelastic behavior of inflatable aerodynamic decelerators.

Tanner, Christopher L.↗