Search NASA⌕ Search

SEARCH · Search NASA

Results for “low order model”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 145 records · Page 8

Tools for Assessing Performance (TAP) Project: FY2023 Quarter 2 Report

LANL personnel presented TAP-sponsored work at the 103 rd Annual Meeting of the American Meteorological Society, which was held in Denver, CO in January 2023. There were four different presentations related to TAP-sponsored work that were given by LANL staff and studentsat this meeting: (1) Gloeckler, T., M. Nelson, H. Shah, M. Deshler, L. Wedell, R. Linn, N. Duboc, S. Pol, P. Bieringer, A. Annunzio, B. Martin, and H. Tinnesand, 2023. Testing Validity and Limitations of a Low-Order Diffusive Wake Model. 15 th Symposium on the Urban Environment, January 8-12, 2023, Denver, CO, U.S.A. (2) Nelson, M., H. Shah, T. Gloeckler, M. Deshler, L. Wedell, A. Ortega, R. Linn, N. Duboc, P. Conry, S. Pol, and H. Tinnesand, 2023. The Isolated Building Wake Experiment. Paper 1A.1. 24 th Symposium on Boundary Layers and Turbulence, January 8-12, 2023. Denver, CO, U.S.A. (3) Shah, H., M. Nelson, S. Pol, R. Linn, N. Duboc, and H. Tinnesand, 2023. Detailed Measurement of a Building Wake Under Real-World Atmospheric Conditions. Poster 102. 24 th Symposium on Boundary Layers and Turbulence, January 8-12, 2023, Denver, CO, U.S.A. (4) Wedell, L., M. Nelson, H. Shah, T. Gloeckler, M. Deshler, R. Linn, S. Pol, H. Tinnesand, 2023. The Isolated Building Wake Experiment: Using a Rotating Building to Mulitply the Number of Wake Measurement Locations without Moving the Instruments. Paper 1A.4A. 24 th Symposium on Boundary Layers and Turbulence, January 8-12, 2023, Denver, CO, U.S.A.

54 ENVIRONMENTAL SCIENCES↗

Processing techniques development, volume 3

The author has identified the following significant results. Analysis of the geometric characteristics of the aircraft synthetic aperture radar (SAR) relative to LANDSAT indicated that relatively low order polynominals would model the distortions to subpixel accuracy to bring SAR into registration for good quality imagery. Also the area analyzed was small, about 10 miles square, so this is an additional constraint. For the Air Force/ERIM data, none of the tested methods could achieve subpixel accuracy. Reasons for this is unknown; however, the noisy (high scintillation) nature of the data and attendent unrecognizability of features contribute to this error. It is concluded that the quadratic model would adequately provide distortion modeling for small areas, i.e., 10 to 20 miles square.

Landgrebe, D. A.↗

Radiation budget measurement/model interface

This final report includes research results from the period February, 1981 through November, 1982. Two new results combine to form the final portion of this work. They are the work by Hanna (1982) and Stevens to successfully test and demonstrate a low-order spectral climate model and the work by Ciesielski et al. (1983) to combine and test the new radiation budget results from NIMBUS-7 with earlier satellite measurements. Together, the two related activities set the stage for future research on radiation budget measurement/model interfacing. Such combination of results will lead to new applications of satellite data to climate problems. The objectives of this research under the present contract are therefore satisfied. Additional research reported herein includes the compilation and documentation of the radiation budget data set a Colorado State University and the definition of climate-related experiments suggested after lengthy analysis of the satellite radiation budget experiments.

Vonderhaar, T. H.↗

Propulsive-jet flow field analysis using the three-dimensional Navier-Stokes equations

A three-dimensional Navier-Stokes code has been applied to the analysis of flow fields containing propulsive jets. Specifically, the application was made to a flow field containing a supersonic jet injected at an angle of 90 degrees to a subsonic free stream. Although wind tunnel data were available, the computational results were not readily comparable to the experimental data because of significant differences between the two plume trajectories. Reasons for the differences are suggested in the report and include: (1) incomplete convergence, (2) inadequate grid resolution in the high gradient regions, and (3) use of a low-order turbulence closure model.

Reed, C. L.↗

Bayesian Revision of Residual Detection Power

This paper addresses some issues with quality assessment and quality assurance in response surface modeling experiments executed in wind tunnels. The role of data volume on quality assurance for response surface models is reviewed. Specific wind tunnel response surface modeling experiments are considered for which apparent discrepancies exist between fit quality expectations based on implemented quality assurance tactics, and the actual fit quality achieved in those experiments. These discrepancies are resolved by using Bayesian inference to account for certain imperfections in the assessment methodology. Estimates of the fraction of out-of-tolerance model predictions based on traditional frequentist methods are revised to account for uncertainty in the residual assessment process. The number of sites in the design space for which residuals are out of tolerance is seen to exceed the number of sites where the model actually fails to fit the data. A method is presented to estimate how much of the design space in inadequately modeled by low-order polynomial approximations to the true but unknown underlying response function.

DeLoach, Richard↗

Stress intensity factor models using mechanics-guided decomposition and symbolic regression

The finite element method can be used to compute accurate stress intensity factors (SIFs) for cracks with complex geometries and boundary conditions. In contrast, handbook solutions act as surrogate SIF models that provide significantly faster evaluation times. However, the development of conventional surrogate SIF models relies on manual development based on low-order parameterizations. This limits surrogate model accuracy and generalizability. Here, in this paper, we develop a framework for the automated development of mechanics-guided handbook SIF solutions by using interpretable machine learning via genetic programming for symbolic regression (GPSR). Formalizing the mechanics-based approach of Raju and Newman, SIF training data is decomposed into multiple subsets. This decomposition enables parallel GPSR model development of subfunctions, each of which accounts for specific geometrical corrections with respect to a known analytical model. Using this mechanics-based approach with GPSR allows for equations to be learned with improved accuracy and reduced complexity relative to the Raju Newman equations while maintaining the inherent interpretability of mathematical expressions. In this paper, we present equations that match the complexity of the Raju Newman equations while having reduced error, as well as equations with similar errors and reduced complexity.

42 ENGINEERING↗

Precision pointing and control of flexible spacecraft

The problem and long term objectives for the precision pointing and control of flexible spacecraft are given. The four basic objectives are stated in terms of two principle tasks. Under Task 1, robust low order controllers, improved structural modeling methods for control applications and identification methods for structural dynamics are being developed. Under Task 2, a lab test experiment for verification of control laws and system identification algorithms is being developed. For Task 1, work has focused on robust low order controller design and some initial considerations for structural modeling in control applications. For Task 2, work has focused on experiment design and fabrication, along with sensor selection and initial digital controller implementation. Conclusions are given.

Bantell, M. H., Jr.↗

Higher Order Modeling In the BEM/FEM Hybrid Formulation

Hybrid formulations using low order curl-conforming bases to represent the total electric field within a finite element region and low order divergence-conforming bases to represent equivalent electric and magnetic currents on the boundary are well known. However, higher-order divergence and curl-conforming bases have been shown to provide significant benefits in convergence rates and accuracy when employed in strictly integral equation and strictly finite element formulations. In this paper, a hybrid electric field formulation employing higher order bases is presented. The paper addresses benefits and issues associated with using higher order divergence-and curl-conforming bases in the hybrid finite element/boundary element electric field formulation. The method of singularity subtraction may be used to compute the self terms of the boundary integral when the bases are of low order. But this method becomes laborious and requires great care when the divergence conforming bases are of higher order. In order to handle these singularities simply and accurately, a generalized Gaussian quadrature method is employed in which the expansion functions account for the singularity. In preliminary tests of the higher order hybrid formulation, the equivalent electric current induced by scattering of a plane wave from a square dielectric cylinder is examined. Accurate results are obtained using only a two-triangle mesh when the current basis is of order 4 or 5. Additional results are presented comparing the error obtained using higher order bases to that obtained using lower order bases when the number of unknowns is approximately equal. Also, convergence rates obtained with higher order bases are compared to those obtained with lower order bases for selected sample problems.

Fink, Patrick W.↗

Compensating for pneumatic distortion in pressure sensing devices

A general numerical technique for obtaining unsteady pressure measurements using conventional pressure sensing technology has been developed. A pneumatic distortion model, based on the Navier-Stokes equations of momentum and continuity, was reduced to a low-order, state-variable model retaining most of the dynamic characteristics of the full model. The reduced-order model is coupled with results from minimum variance estimation theory to develop an algorithm to compensate the effects of pneumatic distortion. Both postflight and real-time algorithms were developed and evaluated using simulated and flight data.

Whitmore, Stephen A.↗

Predictability of weather and climate in a coupled ocean-atmosphere model: A dynamical systems approach

A dynamical systems approach is used to quantify the instantaneous and time-averaged predictability of a low-order moist general circulation model. Specifically, the effects on predictability of incorporating an active ocean circulation, implementing annual solar forcing, and asynchronously coupling the ocean and atmosphere are evaluated. The predictability and structure of the model attractors is compared using the Lyapunov exponents, the local divergence rates, and the correlation, fractal, and Lyapunov dimensions. The Lyapunov exponents measure the average rate of growth of small perturbations on an attractor, while the local divergence rates quantify phase-spatial variations of predictability. These local rates are exploited to efficiently identify and distinguish subtle differences in predictability among attractors. In addition, the predictability of monthly averaged and yearly averaged states is investigated by using attractor reconstruction techniques.

Nese, Jon M.↗

Parametric dynamic mode decomposition for reduced order modeling

Dynamic Mode Decomposition (DMD) is a model-order reduction approach, whereby spatial modes of fixed temporal frequencies are extracted from numerical or experimental data sets. The DMD low-rank or reduced operator is typically obtained by singular value decomposition of the temporal data sets. For parameter-dependent models, as found in many multi-query applications such as uncertainty quantification or design optimization, the only parametric DMD technique developed was a stacked approach, with data sets at multiple parameter values were aggregated together, increasing the computational work needed to devise low-rank dynamical reduced-order models. Here in this paper, we present two novel approach to carry out parametric DMD: one based on the interpolation of the reduced-order DMD eigen-pair and the other based on the interpolation of the reduced DMD (Koopman) operator. Numerical results are presented for diffusion-dominated nonlinear dynamical problems, including a multiphysics radiative transfer example. All three parametric DMD approaches are compared.

97 MATHEMATICS AND COMPUTING↗

Spatiotemporal spread of perturbations in power-law models at low temperatures: Exact results for classical out-of-time-order correlators

Here, we present exact results for the classical version of the out-of-time-order commutator (OTOC) for a family of power-law models consisting of N particles in one dimension and confined by an external harmonic potential. These particles are interacting via power-law interaction of the form ∝ $Σ_{i,j=(i≠j)}^{N}$ |x i -x j | -k for every k>1 where x i is the position of the i th particle. We present numerical results for the OTOC for finite N at low temperatures and short enough times so that the system is well approximated by the linearized dynamics around the many-body ground state. In the large-N limit, we compute the ground-state dispersion relation in the absence of external harmonic potential exactly and use it to arrive at analytical results for OTOC. We find excellent agreement between our analytical results and the numerics. We further obtain analytical results in the limit where only linear and leading nonlinear (in momentum) terms in the dispersion relation are included. The resulting OTOC is in agreement with numerics in the vicinity of the edge of the “light cone.” We find remarkably distinct features in OTOC below and above k=3 in terms of going from non-Airy behavior (1 < k < 3) to an Airy universality class (k>3). We present certain additional rich features for the case k=2 that stem from the underlying integrability of the Calogero-Moser model. We present a field theory approach that also assists in understanding certain aspects of OTOC such as the sound speed. Our findings are a step forward towards a more general understanding of the spatiotemporal spread of perturbations in long-range interacting systems.

1-dimensional systems↗

Fingerprinting triangular-lattice antiferromagnet by excitation gaps

CeCd 3 As 3 is a rare-earth triangular-lattice antiferromagnet with large inter-layer separation. Here, our field-dependent heat capacity measurements at dilution fridge temperatures allow us to trace the field-evolution of the spin-excitation gaps throughout the antiferromagnetic and paramagnetic regions. The distinct gap evolution places strong constraints on the microscopic pseudo-spin model, which, in return, yields a close quantitative description of the gap behavior. This analysis provides crucial insights into the nature of the magnetic state of CeCd 3 As 3 , with a certainty regarding its stripe order and low-energy model parameters that sets a compelling paradigm for exploring and understanding the rapidly growing family of the rare-earth-based triangular-lattice systems

36 MATERIALS SCIENCE↗

Bayesian learning of orthogonal embeddings for multi-fidelity Gaussian Processes

Uncertainty propagation in complex engineering systems often poses significant computational challenges related to modeling and quantifying probability distributions of model outputs, as those emerge as the result of various sources of uncertainty that are inherent in the system under investigation. Gaussian Processes regression (GPs) is a robust meta-modeling technique that allows for fast model prediction and exploration of response surfaces. Multi-fidelity variations of GPs further leverage information from cheap and low fidelity model simulations in order to improve their predictive performance on the high fidelity model. In order to cope with the high volume of data required to train GPs in high dimensional design spaces, a common practice is to introduce latent design variables that are typically projections of the original input space to a lower dimensional subspace, and therefore substitute the problem of learning the initial high dimensional mapping, with that of training a GP on a low dimensional space. Here in this paper, we present a Bayesian approach to identify optimal transformations that map the input points to low dimensional latent variables. The \projection" mapping consists of an orthonormal matrix that is considered a priori unknown and needs to be inferred jointly with the GP parameters, conditioned on the available training data. The proposed Bayesian inference scheme relies on a two-step iterative algorithm that samples from the marginal posteriors of the GP parameters and the projection matrix respectively, both using Markov Chain Monte Carlo (MCMC) sampling. In order to take into account the orthogonality constraints imposed on the orthonormal projection matrix, a Geodesic Monte Carlo sampling algorithm is employed, that is suitable for exploiting probability measures on manifolds. We extend the proposed framework to multi-fidelity models using GPs including the scenarios of training multiple outputs together. We validate our framework on three synthetic problems with a known lower-dimensional subspace. The benefits of our proposed framework, are illustrated on the computationally challenging aerodynamic optimization of a last-stage blade for an industrial gas turbine, where we study the effect of an 85-dimensional shape parameterization of a three-dimensional airfoil on two output quantities of interest, specifically on the aerodynamic efficiency and the degree of reaction

42 ENGINEERING↗

Noise transmission properties and control strategies for composite structures

A study of several component technologies required to apply active control techniques to reduce interior noise in composite aircraft structures is described. The mechanisms of noise transmission in an all composite, large-scale, fuselage model are studied in an experimental program and found similar to mechanisms found in conventional aircraft construction. Two primary conditions of structural acoustic response are found to account for the dominant interior acoustic response. A preliminary study of active noise control in cylinders used piezoceramic actuators as force inputs for a simple aluminum fuselage model. These actuators provided effective control for the same two conditions of noise transmission found in the composite fuselage structure. The use of piezoceramic actuators to apply force inputs overcomes the weight and structural requirements of conventional shaker actuators. Finally, in order to accurately simulate these types of actuators in a cylindrical shell, two analytical models are investigated that apply either in-plane forces or bending moments along the boundaries of a finite patch. It is shown that the bending model may not be as effective as the force model for exciting the low order azimuthal modes that typically dominate the structural acoustic response in these systems. This result will affect the arrangement and distribution of actuators required for effective active control systems.

Silcox, Richard J.↗

A new approach to mixed H2/H infinity controller synthesis using gradient-based parameter optimization methods

In the past few years, the mixed H(sub 2)/H-infinity control problem has been the object of much research interest since it allows the incorporation of robust stability into the LQG framework. The general mixed H(sub 2)/H-infinity design problem has yet to be solved analytically. Numerous schemes have considered upper bounds for the H(sub 2)-performance criterion and/or imposed restrictive constraints on the class of systems under investigation. Furthermore, many modern control applications rely on dynamic models obtained from finite-element analysis and thus involve high-order plant models. Hence the capability to design low-order (fixed-order) controllers is of great importance. In this research a new design method was developed that optimizes the exact H(sub 2)-norm of a certain subsystem subject to robust stability in terms of H-infinity constraints and a minimal number of system assumptions. The derived algorithm is based on a differentiable scalar time-domain penalty function to represent the H-infinity constraints in the overall optimization. The scheme is capable of handling multiple plant conditions and hence multiple performance criteria and H-infinity constraints and incorporates additional constraints such as fixed-order and/or fixed structure controllers. The defined penalty function is applicable to any constraint that is expressible in form of a real symmetric matrix-inequity.

Ly, Uy-Loi↗

The Interplay Between Radiation Pressure and the Photoelectric Instability in Optically Thin Disks of Gas and Dust

In optically thin disks, dust grains are photoelectrically stripped of electrons by starlight, heating nearby gas and possibly creating a dust clumping instability-the photoelectric instability (PeI)-that significantly alters global disk structure. In the current work, we use the Pencil Code to perform the first numerical models of the PeI that include stellar radiation pressure on dust grains in order to explore the parameter regime in which the instability operates. In some models with low gas and dust surface densities, we see a variety of dust structures, including sharp concentric rings. In the most gas- and dust-rich models, nonaxisymmetric clumps, arcs, and spiral arms emerge that represent dust surface density enhancements of factors of ∼5-20. In one high gas surface density model, we include a large, low-order gas viscosity and find that it observably smooths the structures that form in the gas and dust, suggesting that resolved images of a given disk may be useful for deriving constraints on the effective viscosity of its gas. Our models show that radiation pressure does not preclude the formation of complex structure from the PeI, but the qualitative manifestation of the PeI depends strongly on the parameters of the system. The PeI may provide an explanation for unusual disk morphologies, such as the moving blobs of the AU Mic disk, the asymmetric dust distribution of the 49Ceti disk, and the rings and arcs found in the HD 141569A disk.

photoelectric instability (PeI)↗

Theories of multiple equilibria and weather regimes : A critical reexamination. II - Baroclinic two-layer models

Previous results based on low- and intermediate-order truncations of the two-layer model suggest the existence of multiple equilibria and/or multiple weather regimes for the extratropical large-scale flow. The importance of the transient waves in the synoptic scales in organizing the large-scale flow and in the maintenance of weather regimes was emphasized. The result shows that multiple equilibria/weather regimes that are present in lower-order models examined disappear when a sufficient number of modes are kept in the spectral expansion of the solution to the governing partial differential equations. Much of the chaotic behavior of the large-scale flow that is present in intermediate-order models is now found to be spurious. Physical reasons for the drastic modification are offered. A peculiarity in the formulation of most existing two-layer models is noted that also tends to exaggerate the importance of baroclinic processes and increase the degree of unpredictability of the large-scale flow.

Cehelsky, Priscilla↗