Search NASA⌕ Search

SEARCH · Search NASA

Results for “Generalized Linear 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 397 records · Page 22

Three-dimensional structures of equatorial waves and the resulting super-rotation in the atmosphere of a tidally locked hot Jupiter

Three-dimensional (3D) equatorial trapped waves excited by stellar isolation and the resulting equatorial super-rotating jet in a vertical stratified atmosphere of a tidally locked hot Jupiter are investigated. Taking the hot Jupiter HD 189733b as a fiducial example, we analytically solve linear equations subject to stationary stellar heating with a uniform zonal-mean flow included. We also extract wave information in the final equilibrium state of the atmosphere from our radiative hydrodynamical simulation for HD 189733b. Our analytic wave solutions are able to qualitatively explain the 3D simulation results. Apart from previous wave studies, investigating the vertical structure of waves allows us to explore new wave features such as the wavefronts tilts related to the Rossby-wave resonance as well as dispersive equatorial waves. We also attempt to apply our linear wave analysis to explain some numerical features associated with the equatorial jet development seen in the general circulation model by Showman and Polvani. During the spin-up phase of the equatorial jet, the acceleration of the jet can be in principle boosted by the Rossby-wave resonance. However, we also find that as the jet speed increases, the Rossby-wave structure shifts eastward, while the Kelvin-wave structure remains approximately stationary, leading to the decline of the acceleration rate. Our analytic model of jet evolution implies that there exists only one stable equilibrium state of the atmosphere, possibly implying that the final state of the atmosphere is independent of initial conditions in the linear regime. Limitations of our linear model and future improvements are also discussed.

79 ASTRONOMY AND ASTROPHYSICS↗

Classification of constraints and degrees of freedom for quadratic discrete actions

We provide a comprehensive classification of constraints and degrees of freedom for variational discrete systems governed by quadratic actions. This classification is based on the different types of null vectors of the Lagrangian two-form and employs the canonical formalism developed in Dittrich and Höhn [“Constraint analysis for variational discrete systems,” J. Math. Phys. 54, 093505 (2013); e-print http://arxiv.org/abs/arXiv:1303.4294 [math-ph]] and Höhn [“Quantization of systems with temporally varying discretization I: Evolving Hilbert spaces,” J. Math. Phys. 55, 083508 (2014); e-print http://arxiv.org/abs/arXiv:1401.6062 [gr-qc]]. The analysis is carried out in both the classical and quantum theory and applies to systems with both temporally varying or constant discretization. In particular, it is shown explicitly how changes in the discretization, e.g., resulting from canonical coarse graining or refining operations or an evolving background geometry, change the dynamical content of the system. It is demonstrated how, on a temporally varying discretization, constraints, Dirac observables, symmetries, reduced phase spaces, and physical Hilbert spaces become spacetime region dependent. These results are relevant for free field theory on an evolving lattice and linearized discrete gravity models.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

The utilization of satellite data and dynamics in understanding and predicting global weather phenomena

A two layer spectral quasi-geostrophic model is used to simulate the effects of topography on the equilibria, the stability, and the long term evaluation of incipient unstable waves. The flow is forced by latitudinally dependent radiational heating. The nature of the form drag instability of high index equilibria is investigated. The proximity of the equilibrium shear to a resonant value is essential for the instability, provided the equilibrium occurs at a slightly stronger shear than resonance. The properties of the steady Hadley and Rossby required for a thermally forced rotating fluid on a sphere are further explained. An objective parameterization technique is developed for general nonlinear hydrodynamical systems. The typical structure is one in which the rates of change of the dependent variables depend on homogeneous quadratic and linear forms, as well as on inhomogeneous forcing terms. Also documented is a steady, axisymmetric model of the general circulation developed as a basis for climate stability studies. The model includes the effects of heating, rotation, and internal friction, but neglects topography. Included is further research on cloud street phenomena. Orientation angles and horizontal wavelengths of boundary layer rolls and cloud streets are determined from an analysis of a truncated spectral model of three dimensional shallow moist Boussinesq convection in a shearing environment is further explained. Relatively broadly spaced roll clouds have orientations for which the Fourier component of the roll perpendicular shear is nearly zero, but the second corresponds to narrowly spaced rolls having orientations for which the Fourier coefficients of both the perpendicular and the parallel components of the shear are nearly equal.

Shirer, H. N.↗

Computer Model Calibration Based on Image Warping Metrics: An Application for Sea Ice Deformation

Arctic sea ice plays an important role in the global climate. Sea ice models governed by physical equations have been used to simulate the state of the ice including characteristics such as ice thickness, concentration, and motion. More recent models also attempt to capture features such as fractures or leads in the ice. These simulated features can be partially misaligned or misshapen when compared to observational data, whether due to numerical approximation or incomplete physics. In order to make realistic forecasts and improve understanding of the underlying processes, it is necessary to calibrate the numerical model to field data. Traditional calibration methods based on generalized least-square metrics are flawed for linear features such as sea ice cracks. We develop a statistical emulation and calibration framework that accounts for feature misalignment and misshapenness, which involves optimally aligning model output with observed features using cutting-edge image registration techniques. Here, this work can also have application to other physical models which produce coherent structures.

54 ENVIRONMENTAL SCIENCES↗

Mathematical modeling of spinning elastic bodies for modal analysis.

The problem of modal analysis of an elastic appendage on a rotating base is examined to establish the relative advantages of various mathematical models of elastic structures and to extract general inferences concerning the magnitude and character of the influence of spin on the natural frequencies and mode shapes of rotating structures. In realization of the first objective, it is concluded that except for a small class of very special cases the elastic continuum model is devoid of useful results, while for constant nominal spin rate the distributed-mass finite-element model is quite generally tractable, since in the latter case the governing equations are always linear, constant-coefficient, ordinary differential equations. Although with both of these alternatives the details of the formulation generally obscure the essence of the problem and permit very little engineering insight to be gained without extensive computation, this difficulty is not encountered when dealing with simple concentrated mass models.

Likins, P. W.↗

Minimal subspace rotation on the Stiefel manifold for stabilization and enhancement of projection-based reduced order models for the compressible Navier–Stokes equations

For a projection-based reduced order model (ROM) of a fluid flow to be stable and accurate, the dynamics of the truncated subspace must be taken into account. This paper proposes an approach for stabilizing and enhancing projection-based fluid ROMs in which truncated modes are accounted for a priori via a minimal rotation of the projection subspace. Attention is focused on the full non-linear compressible Navier–Stokes equations in specific volume form as a step toward a more general formulation for problems with generic non-linearities. Unlike traditional approaches, no empirical turbulence modeling terms are required, and consistency between the ROM and the Navier–Stokes equation from which the ROM is derived is maintained. Mathematically, the approach is formulated as a trace minimization problem on the Stiefel manifold. As a result, the reproductive as well as predictive capabilities of the method are evaluated on several compressible flow problems, including a problem involving laminar flow over an airfoil with a high angle of attack, and a channel-driven cavity flow problem.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A class of the van Leer-type transport schemes and its application to the moisture transport in a general circulation model

A generalized form of the second-order van Leer transport scheme is derived. Several constraints to the implied subgrid linear distribution are discussed. A very simple positive-definite scheme can be derived directly from the generalized form. A monotonic version of the scheme is applied to the Goddard Laboratory for Atmospheres (GLA) general circulation model (GCM) for the moisture transport calculations, replacing the original fourth-order center-differencing scheme. Comparisons with the original scheme are made in idealized tests as well as in a summer climate simulation using the full GLA GCM. A distinct advantage of the monotonic transport scheme is its ability to transport sharp gradients without producing spurious oscillations and unphysical negative mixing ratio. Within the context of low-resolution climate simulations, the aforementioned characteristics are demonstrated to be very beneficial in regions where cumulus convection is active. The model-produced precipitation pattern using the new transport scheme is more coherently organized both in time and in space, and correlates better with observations. The side effect of the filling algorithm used in conjunction with the original scheme is also discussed, in the context of idealized tests. The major weakness of the proposed transport scheme with a local monotonic constraint is its substantial implicit diffusion at low resolution. Alternative constraints are discussed to counter this problem.

Lin, Shian-Jiann↗

Orbital optimisation in the perfect pairing hierarchy: applications to full-valence calculations on linear polyacenes

Here, we describe the implementation of orbital optimisation for the models in the perfect pairing hierarchy. Orbital optimisation, which is generally necessary to obtain reliable results, is pursued at perfect pairing (PP) and perfect quadruples (PQ) levels of theory for applications on linear polyacenes, which are believed to exhibit strong correlation in the π space. While local minima and σ-π symmetry breaking solutions were found for PP orbitals, no such problems were encountered for PQ orbitals. The PQ orbitals are used for single-point calculations at PP, PQ and perfect hextuples (PH) levels of theory, both only in the π subspace, as well as in the full σπ valence space. It is numerically demonstrated that the inclusion of single excitations is necessary also when optimised orbitals are used. PH is found to yield good agreement with previously published density matrix renormalisation group data in the π space, capturing over 95% of the correlation energy. Full-valence calculations made possible by our novel, efficient code reveal that strong correlations are weaker when larger basis sets or active spaces are employed than in previous calculations. The largest full-valence PH calculations presented correspond to a (192e,192o) problem.

74 ATOMIC AND MOLECULAR PHYSICS↗

Sampling-based Sublinear Low-rank Matrix Arithmetic Framework for Dequantizing Quantum Machine Learning

We present an algorithmic framework for quantum-inspired classical algorithms on close-to-low-rank matrices, generalizing the series of results started by Tang’s breakthrough quantum-inspired algorithm for recommendation systems [STOC’19]. Motivated by quantum linear algebra algorithms and the quantum singular value transformation (SVT) framework of Gilyén et al. [STOC’19], we develop classical algorithms for SVT that run in time independent of input dimension, under suitable quantum-inspired sampling assumptions. Our results give compelling evidence that in the corresponding QRAM data structure input model, quantum SVT does not yield exponential quantum speedups. Since the quantum SVT framework generalizes essentially all known techniques for quantum linear algebra, our results, combined with sampling lemmas from previous work, suffice to generalize all prior results about dequantizing quantum machine learning algorithms. In particular, our classical SVT framework recovers and often improves the dequantization results on recommendation systems, principal component analysis, supervised clustering, support vector machines, low-rank regression, and semidefinite program solving. We also give additional dequantization results on low-rank Hamiltonian simulation and discriminant analysis. Our improvements come from identifying the key feature of the quantum-inspired input model that is at the core of all prior quantum-inspired results: ℓ 2 -norm sampling can approximate matrix products in time independent of their dimension. We reduce all our main results to this fact, making our exposition concise, self-contained, and intuitive.

Computer Science↗

Signal propagation in reversible digital mechanics

Digital mechanics explores information processing through binary, mechanical circuits. This work demonstrates a flexural, mechanical integrated circuit (m-IC) that achieves reversible, non-reciprocal signal propagation through integrated AND logic and memory. Our approach exploits sequential bistable transitions with symmetric energy wells, tunable stiffness, impedance matching, and AND gate non-linearity, to enable signal propagation, repeatability, and reversibility. We present a generalized model of logic kinematics and energetics, validated experimentally, to study energy flows, quantify energetic limits, and identify operating regimes for reversible logic. Macro-scale experiments confirm propagation dynamics, and new fabrication methods extend the architecture to micro-scale devices. By achieving controlled, reversible signal transmission across interconnected logic and memory, this work establishes a scalable platform for robust mechanical computing and adaptive sensing.

Johnson, Hilary A. [Lawrence Livermore National La↗

Over-reflection of slow magnetosonic waves by homogeneous shear flow: Analytical solution

We have analyzed the amplification of slow magnetosonic (or pseudo-Alfvenic) waves (SMW) in incompressible shear flow. As found here, the amplification depends on the component of the wave-vector perpendicular to the direction of the shear flow. Earlier numerical results are consistent with the general analytic solution for the linearized magnetohydrodynamic equations, derived here for the model case of pure homogeneous shear (without Coriolis force). An asymptotically exact analytical formula for the amplification coefficient is derived for the case when the amplification is sufficiently large.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Data-driven predictive probability density function control of fiber length stochastic distribution shaping in refining process

Pulp is the most important raw material in paper in paper industries, whose Fiber Length Stochastic Distribution (FLSD) shaping directly determines the energy consumption and paper quality of the subsequent papermaking processes. However, the mean and variance are insufficient to describe the output FLSD shaping, which display non-Gaussian distribution properties. Therefore, the traditional control, method based on the mean and variance of fiber length is difficult to control the FLSD shaping effectively. In this paper, a novel data-driven predictive probability density function (PDF) control method is proposed for output FLSD shaping in refining process. Primarily, in order to improve the approximation accuracy of the PDF employing RBF neural network, the parameters (center value and width) of each RBF basis functions are tuned via utilizing iterative learning control (ILC) law, and the corresponding estimations of weights law can be obtained. Secondly, considering that the conventional linear model of weights vector has drawbacks of low accuracy and weak generalization ability, random vector functional link networks (RVFLNs) based data driven nonlinear modelling method is employed to characterize the prediction model between the input variables and weights vector. Finally, in order to reduce the randomness of the output FLSD, a minimum entropy control method under mean constraint is employed to design the predictive PDF controller for the output FLSD shaping with the help of the established nonlinear stochastic distribution model. Both simulations and experiments are given to illustrate the effectiveness and practicability of the proposed method.

Probability debsity function, Fiber length distrib↗

Development of control strategies for safe microburst penetration: A progress report

A single-engine, propeller-driven, general-aviation model was incorporated into the nonlinear simulation and into the linear analysis of root loci and frequency response. Full-scale wind tunnel data provided its aerodynamic model, and the thrust model included the airspeed dependent effects of power and propeller efficiency. Also, the parameters of the Jet Transport model were changed to correspond more closely to the Boeing 727. In order to study their effects on steady-state repsonse to vertical wind inputs, altitude and total specific energy (air-relative and inertial) feedback capabilities were added to the nonlinear and linear models. Multiloop system design goals were defined. Attempts were made to develop controllers which achieved these goals.

Psiaki, Mark L.↗

Helicopter mathematical models and control law development for handling qualities research

Progress made in joint NASA/Army research concerning rotorcraft flight-dynamics modeling, design methodologies for rotorcraft flight-control laws, and rotorcraft parameter identification is reviewed. Research into these interactive disciplines is needed to develop the analytical tools necessary to conduct flying qualities investigations using both the ground-based and in-flight simulators, and to permit an efficient means of performing flight test evaluation of rotorcraft flying qualities for specification compliance. The need for the research is particularly acute for rotorcraft because of their mathematical complexity, high order dynamic characteristics, and demanding mission requirements. The research in rotorcraft flight-dynamics modeling is pursued along two general directions: generic nonlinear models and nonlinear models for specific rotorcraft. In addition, linear models are generated that extend their utilization from 1-g flight to high-g maneuvers and expand their frequency range of validity for the design analysis of high-gain flight control systems. A variety of methods ranging from classical frequency-domain approaches to modern time-domain control methodology that are used in the design of rotorcraft flight control laws is reviewed. Also reviewed is a study conducted to investigate the design details associated with high-gain, digital flight control systems for combat rotorcraft. Parameter identification techniques developed for rotorcraft applications are reviewed.

Chen, Robert T. N.↗

Reliability of High-Voltage Tantalum Capacitors. Parts 3 and 4)

Weibull grading test is a powerful technique that allows selection and reliability rating of solid tantalum capacitors for military and space applications. However, inaccuracies in the existing method and non-adequate acceleration factors can result in significant, up to three orders of magnitude, errors in the calculated failure rate of capacitors. This paper analyzes deficiencies of the existing technique and recommends more accurate method of calculations. A physical model presenting failures of tantalum capacitors as time-dependent-dielectric-breakdown is used to determine voltage and temperature acceleration factors and select adequate Weibull grading test conditions. This model is verified by highly accelerated life testing (HALT) at different temperature and voltage conditions for three types of solid chip tantalum capacitors. It is shown that parameters of the model and acceleration factors can be calculated using a general log-linear relationship for the characteristic life with two stress levels.

Teverovsky, Alexander↗

Analysis of Weibull Grading Test for Solid Tantalum Capacitors

Weibull grading test is a powerful technique that allows selection and reliability rating of solid tantalum capacitors for military and space applications. However, inaccuracies in the existing method and non-adequate acceleration factors can result in significant, up to three orders of magnitude, errors in the calculated failure rate of capacitors. This paper analyzes deficiencies of the existing technique and recommends more accurate method of calculations. A physical model presenting failures of tantalum capacitors as time-dependent-dielectric-breakdown is used to determine voltage and temperature acceleration factors and select adequate Weibull grading test conditions. This, model is verified by highly accelerated life testing (HALT) at different temperature and voltage conditions for three types of solid chip tantalum capacitors. It is shown that parameters of the model and acceleration factors can be calculated using a general log-linear relationship for the characteristic life with two stress levels.

Teverovsky, Alexander↗

An asymptotically compatible treatment of traction loading in linearly elastic peridynamic fracture

Meshfree discretizations of state-based peridynamic models are attractive due to their ability to naturally describe fracture of general materials. However, two factors conspire to prevent meshfree discretizations of state-based peridynamics from converging to corresponding local solutions as resolution is increased: quadrature error prevents an accurate prediction of bulk mechanics, and the lack of an explicit boundary representation presents challenges when applying traction loads. Herein, we develop a reformulation of the linear peridynamic solid (LPS) model to address these shortcomings, using improved meshfree quadrature, a reformulation of the nonlocal dilatation, and a consistent handling of the nonlocal traction condition to construct a model with rigorous accuracy guarantees. In particular, these improvements are designed to enforce discrete consistency in the presence of evolving fractures, whose a priori unknown location render consistent treatment difficult. In the absence of fracture, when a corresponding classical continuum mechanics model exists, our improvements provide asymptotically compatible convergence to corresponding local solutions, eliminating surface effects and issues with traction loading which have historically plagued peridynamic discretizations. When fracture occurs, our formulation automatically provides a sharp representation of the fracture surface by breaking bonds, avoiding the loss of mass. We provide rigorous error analysis and demonstrate convergence for a number of benchmarks, including manufactured solutions, free-surface, nonhomogeneous traction loading, and composite material problems. Finally, we validate simulations of brittle fracture against a recent experiment of dynamic crack branching in soda-lime glass, providing evidence that the scheme yields accurate predictions for practical engineering problems.

42 ENGINEERING↗

ISPAN (Interactive Stiffened Panel Analysis): A tool for quick concept evaluation and design trade studies

Interactive Stiffened Panel Analysis (ISPAN) modules, written in FORTRAN, were developed to provide an easy to use tool for creating finite element models of composite material stiffened panels. The modules allow the user to interactively construct, solve and post-process finite element models of four general types of structural panel configurations using only the panel dimensions and properties as input data. Linear, buckling and post-buckling solution capability is provided. This interactive input allows rapid model generation and solution by non finite element users. The results of a parametric study of a blade stiffened panel are presented to demonstrate the usefulness of the ISPAN modules. Also, a non-linear analysis of a test panel was conducted and the results compared to measured data and previous correlation analysis.

Hairr, John W.↗