Search NASA⌕ Search

SEARCH · Search NASA

Results for “fully 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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
NASA NTRS1990

Preconditioned conjugate gradient methods for the compressible Navier-Stokes equations

The compressible Navier-Stokes equations are solved for a variety of two-dimensional inviscid and viscous problems by preconditioned conjugate gradient-like algorithms. Roe's flux difference splitting technique is used to discretize the inviscid fluxes. The viscous terms are discretized by using central differences. An algebraic turbulence model is also incorporated. The system of linear equations which arises out of the linearization of a fully implicit scheme is solved iteratively by the well known methods of GMRES (Generalized Minimum Residual technique) and Chebyschev iteration. Incomplete LU factorization and block diagonal factorization are used as preconditioners. The resulting algorithm is competitive with the best current schemes, but has wide applications in parallel computing and unstructured mesh computations.

Venkatakrishnan, V.↗
NASA NTRS1994

(abstract) Observation by Clementine of a Hugh Opposition Surge on the Moon at Very Small Solar Phase Angles

The Clementine mission enabled the first quantative observations of the Moon at very small solar phase angles. It is well established from Earth-based observations that the Moon exhibits a non-linear increase in brightness as its face becomes fully illuminated to a terrestrial observer. Because the models of the opposition effect are sensitive indicators of surfical compaction state and particle size, observations at small solar phase angles are important to obtain. A recent model for optical coherent backscatter seeks to explain a narrow opposition spike at very small phase angles (< 1 degree) seen on several icy satellites , including Europa, Icarus 90, and Oberon . Over 90 images of the Moon's surface at phase angles less than 0.5(deg) (the minimum phase angle observable from Earth due to the Moon's angular size) were obtained by the Clementine spacecraft. Our analysis of these images shows the moon exhibits a surge in brightness of approximately 20% below 0.25(deg), comparable to the values observed on icy satellites. No color dependence was detected in the lunar phase curve below 2 degrees.

partical size icy satellites Europa Icarus 90 Ober↗
DOE OSTI2023

Deep-Learning-Based Koopman Modeling for Online Control Synthesis of Nonlinear Power System Transient Dynamics

Power system stability and control have become more challenging due to the increasing uncertainty associated with renewable generation. Here, the performance of conventional control is highly driven by the physics-based offline-developed dynamic models that can deviate from the actual system characteristics under different operating conditions and/or configurations. Data-driven approaches based on online measurements can be a better solution to addressing these issues by capturing real-time operation conditions. This article describes a novel fully data-driven probabilistic framework to derive a linear representation of postcontingency grid dynamics and online prescribe control based on the derived model to enhance transient stability. The complex nonlinear power system dynamics is approximated by a linear model by using multiple neural network modules that infer distributions of the observations and introducing a Koopman layer to sample possible Koopman linear models from the inferred distributions. The trained model features linearity that can be easily incorporated into the existing linear control design paradigm and ease the controller design process. The effectiveness of Koopman-based control designs is validated through comparative case studies, which demonstrate increased prediction accuracy and control performance when applied to a power system with heterogeneous generator dynamics.

24 POWER TRANSMISSION AND DISTRIBUTION↗
NASA NTRS1998

On Wave Processes in the Solar Atmosphere

This grant was awarded by NASA/MSFC to The University of Alabama in Huntsville (UAH) to investigate the physical processes responsible for heating and wind acceleration in the solar atmosphere, and to construct theoretical, self-consistent and time-dependent solar wind models based on the momentum deposition by finite amplitude and nonlinear Alfven waves. In summary, there are three main goals of the proposed research: (1) Calculate the wave energy spectra and wave energy fluxes carried by magnetic non- magnetic waves. (2) Find out which mechanism dominates in supplying the wave energy to different parts of the solar atmosphere. (3) Use the results obtained in (1) and (2) to construct theoretical, self-consistent and time- dependent models of the solar wind. We have completed the first goal by calculating the amount of non-radiative energy generated in the solar convection zone as acoustic waves and as magnetic tube waves. To calculate the amount of wave energy carried by acoustic waves, we have used the Lighthill-Stein theory for sound generation modified by Musielak, Rosner, Stein & Ulmschneider (1994). The acoustic wave energy fluxes for stars located in different regions of the Hertzsprung-Russell (H-R) diagram have also been computed. The wave energy fluxes carried by longitudinal and transverse waves along magnetic flux tubes have been calculated by using both analytical and numerical methods. Our analytical approach is based a theory developed by Musielak, Rosner & Ulmschnelder and Musielak, Rosner, Gall & Ulmschneider, which allows computing the wave energy fluxes for linear tube waves. A numerical approach has been developed by Huang, Musielak & Ulmschneider and Ulmschneider & Musielak to compute the energy fluxes for nonlinear tube waves. Both methods have been used to calculate the wave energy fluxes for stars located in different regions of the HR diagram (Musielak, Rosner & Ulmschneider 1998; Ulmschneider, Musielak & Fawzy 1998). Having obtained the wave energy fluxes for acoustic and magnetic tube waves, we have investigated the behavior of these waves in the solar and stellar atmospheres. The results of our extensive studies have been published in many papers and presented at numerous scientific meetings. In these studies we have investigated different aspects of propagation of acoustic and magnetic waves, the efficiency of energy transfer along magnetic structures in the solar atmosphere, and behavior of Alfven waves in stgeady and expanding solar and stellar atmospheres. Recently, we have used some of these results to construct first purely theoretical, two component and time-dependent models of solar and stellar chromospheres. Finally, to address the third goal, we have constructed first fully theoretical, self-consistent and time dependent wind models based on the momentum deposition by non-linear Alfven waves. The full set of single-fluid MHD equations with the background flow has been solved by using a modified version of the ZEUS MHD code. The constructed wind models are radially symmetric with the magnetic field decreasing radially and the initial outflow is described by the standard Parker wind solution. In contrast to previous studies, no assumptions regarding wave linearity, wave damping, and wave-flow interaction are made; the models thus naturally account for the backreaction of the wind on the waves as well as for the nonlinear interaction between different types of MHD waves. The models have been used to explain the origin of fast speed streams in solar coronal holes. The obtained results clearly demonstrate that the momentum deposition by Alfven waves in the solar wind can be sufficient to explain the origin of fast stream components of the solar wind. The range of wave amplitudes required to obtain the desired results seems to be in good agreement with recent observations.

Musielak, Z. E.↗
DOE OSTI2013

Multi-channel conduction in redox-based resistive switch modelled using quantum point contact theory

A simple analytic model for the electron transport through filamentary-type structures in Si-rich silica (SiO{sub x})-based resistive switches is proposed. The model is based on a mesoscopic description and is able to account for the linear and nonlinear components of conductance that arise from both fully and partially formed conductive channels spanning the dielectric film. Channels are represented by arrays of identical scatterers whose number and quantum transmission properties determine the current magnitude in the low and high resistance states. We show that the proposed model not only reproduces the experimental current-voltage (I-V) characteristics but also the normalized differential conductance (dln(I)/dln(V)-V) curves of devices under test.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Electron cyclotron drift instability and anomalous transport: two-fluid moment theory and modeling

In the presence of a strong electric field perpendicular to the magnetic field, the electron cross-field (E × B) flow relative to the unmagnetized ions can cause the so-called electron cyclotron drift instability (ECDI) due to resonances of the ion acoustic mode and the electron cyclotron harmonics. This occurs in, for example, collisionless shock ramps in space, and in E × B discharge devices such as Hall thrusters. A prominent feature of ECDI is its capability to induce an electron flow parallel to the background >E field at a speed greatly exceeding predictions by classical collision theory. Such anomalous transport is important due to its role in particle thermalization at space shocks, and in causing plasma flows towards the walls of E × B devices, leading to unfavorable erosion and performance degradation, etc. The development of ECDI and anomalous transport is often considered requiring a fully kinetic treatment. In this work, however, we demonstrate that a reduced variant of this instability, and more importantly, the associated anomalous transport, can be treated self-consistently in a collisionless two-fluid framework without any adjustable collision parameter. By treating both electron and ion species on an equal footing, the free energy due to the inter-species velocity shear allows the growth of an anomalous electron flow parallel to the background E field. We will first present linear analyses of the instability in the two-fluid five- and ten-moment models, and compare them against the fully-kinetic theory. At low temperatures, the two-fluid models predict the fastest-growing mode in good agreement with the kinetic result. Also, by including more (> = 10) moments, secondary (and possibly higher) unstable branches can be recovered. The dependence of the instability on ion-to-electron mass ratio, plasma temperature, and background B field strength is also thoroughly explored. We then carry out direct numerical simulations of the cross-field setup using the five-moment model. The development of the instability, as well as the anomalous transport, is confirmed and in excellent agreement with theoretical predictions. The force balance properties are also studied using the five-moment simulation data. Here this work casts new insights into the nature of ECDI and the associated anomalous transport and demonstrates the potential of the two-fluid moment model in efficient modeling of E × B plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Bad Metallic Transport in a Modified Hubbard Model

Strongly correlated metals often display anomalous transport, including T -linear resistivity above the Mott-Ioffe-Regel limit. We introduce a tractable microscopic model for bad metals, by restoring in the well-known Hubbard model—with hopping t and on-site repulsion U —a “screened Coulomb” interaction between charge densities that decays exponentially with spatial separation. This interaction lifts the extensive degeneracy in the spectrum of the t = 0 Hubbard model, allowing us to fully characterize the small t electric, thermal, and thermoelectric transport in our strongly correlated model. Throughout the phase diagram we observe T -linear resistivity above the Mott-Ioffe-Regel limit, together with strong violation of the Wiedemann-Franz law and a large thermopower that can undergo sign change.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2019

Wind Turbine Controller to Mitigate Structural Loads on a Floating Wind Turbine Platform

This paper summarizes the control design work that was performed to optimize the controller of a wind turbine on the WindFloat structure. The WindFloat is a semisubmersible floating platform designed to be a support structure for a multimegawatt power-generating wind turbine. A controller developed for a bottom-fixed wind turbine configuration was modified for use when the turbine is mounted on the WindFloat platform. This results in an efficient platform heel resonance mitigation scheme. In addition, several control modules, designed with a coupled linear model, were added to the fixed-bottom baseline controller. The approach was tested in a fully coupled nonlinear aero-hydro-elastic simulation tool in which wind and wave disturbances were modeled. This testing yielded significant improvements in platform global performance and tower-base bending loading.

17 WIND ENERGY↗
DOE OSTI2016

Wind Turbine Controller to Mitigate Structural Loads on a Floating Wind Turbine Platform

This paper summarizes the control design work that was performed to optimize the controller of a wind turbine on the WindFloat structure. The WindFloat is a semi-submersible floating platform designed to be a support structure for a multi-megawatt power-generating wind turbine. A controller developed for a bottom-fixed wind turbine configuration was modified for use when the turbine is mounted on the WindFloat platform. This results in an efficient platform heel resonance mitigation scheme. In addition several control modules, designed with a coupled linear model, were added to the fixed-bottom baseline controller. The approach was tested in a fully coupled nonlinear aero-hydroelastic simulation tool in which wind and wave disturbances were modeled. This testing yielded significant improvements in platform global performance and tower-base-bending loading.

17 WIND ENERGY↗
DOE OSTI2024

CORN (Crop Optimization Realized through Neuralnets)

Traditional models of predicting plant traits are limited because they often rely on linear assumptions that do not fully capture the complexity of biological interactions and DNA-based markers which are static across environments. This project generated a set of RNA-based data from large multi-environment field trials and combined it with advanced machine learning techniques to account for these complex interactions and improve the accuracy of predictions.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2015

Verification of nonlinear particle simulation of radio frequency waves in tokamak

Nonlinear simulation model for radio frequency waves in fusion plasmas has been developed and verified using fully kinetic ion and drift kinetic electron. Ion cyclotron motion in the toroidal geometry is implemented using Boris push in the Boozer coordinates. Linear dispersion relation and nonlinear particle trapping are verified for the lower hybrid wave and ion Bernstein wave (IBW). Parametric decay instability is observed where a large amplitude pump wave decays into an IBW sideband and an ion cyclotron quasimode (ICQM). The ICQM induces an ion perpendicular heating, with a heating rate proportional to the pump wave intensity.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

Is the Jeffreys' scale a reliable tool for Bayesian model comparison in cosmology?

We are entering an era where progress in cosmology is driven by data, and alternative models will have to be compared and ruled out according to some consistent criterium. The most conservative and widely used approach is Bayesian model comparison. In this paper we explicitly calculate the Bayes factors for all models that are linear with respect to their parameters. We do this in order to test the so called Jeffreys' scale and determine analytically how accurate its predictions are in a simple case where we fully understand and can calculate everything analytically. We also discuss the case of nested models, e.g. one with M{sub 1} and another with M{sub 2} superset of M{sub 1} parameters and we derive analytic expressions for both the Bayes factor and the figure of Merit, defined as the inverse area of the model parameter's confidence contours. With all this machinery and the use of an explicit example we demonstrate that the threshold nature of Jeffreys' scale is not a ''one size fits all'' reliable tool for model comparison and that it may lead to biased conclusions. Furthermore, we discuss the importance of choosing the right basis in the context of models that are linear with respect to their parameters and how that basis affects the parameter estimation and the derived constraints.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Time integration for diffuse interface models for two-phase flow

We propose a variant of the θ-scheme for diffuse interface models for two-phase flow, together with three new linearization techniques for the surface tension. These involve either additional stabilizing force terms, or a fully implicit coupling of the Navier–Stokes and Cahn–Hilliard equation. In the common case that the equations for interface and flow are coupled explicitly, we find a time step restriction which is very different to other two-phase flow models and in particular is independent of the grid size. We also show that the proposed stabilization techniques can lift this time step restriction. Even more pronounced is the performance of the proposed fully implicit scheme which is stable for arbitrarily large time steps. We demonstrate in a Taylor-flow application that this superior coupling between flow and interface equation can decrease the computation time by several orders of magnitude.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Robust scalable initialization for Bayesian variational inference with multi-modal Laplace approximations

Predictive modeling typically relies on Bayesian model calibration to provide uncertainty quantification. Variational inference utilizing fully independent (“mean-field”) Gaussian distributions are often used as approximate probability density functions. This simplification is attractive since the number of variational parameters grows only linearly with the number of unknown model parameters. However, the resulting diagonal covariance structure and unimodal behavior can be too restrictive to provide useful approximations of intractable Bayesian posteriors that exhibit highly non-Gaussian behavior, including multimodality. High-fidelity surrogate posteriors for these problems can be obtained by considering the family of Gaussian mixtures. Gaussian mixtures are capable of capturing multiple modes and approximating any distribution to an arbitrary degree of accuracy, while maintaining some analytical tractability. Unfortunately, variational inference using Gaussian mixtures with full-covariance structures suffers from a quadratic growth in variational parameters with the number of model parameters. The existence of multiple local minima due to strong nonconvex trends in the loss functions often associated with variational inference present additional complications, These challenges motivate the need for robust initialization procedures to improve the performance and computational scalability of variational inference with mixture models. In this work, we propose a method for constructing an initial Gaussian mixture model approximation that can be used to warm-start the iterative solvers for variational inference. The procedure begins with a global optimization stage in model parameter space. In this step, local gradient-based optimization, globalized through multistart, is used to determine a set of local maxima, which we take to approximate the mixture component centers. Around each mode, a local Gaussian approximation is constructed via the Laplace approximation. Finally, the mixture weights are determined through constrained least squares regression. The robustness and scalability of the proposed methodology is demonstrated through application to an ensemble of synthetic tests using high-dimensional, multimodal probability density functions. Here, the practical aspects of the approach are demonstrated with inversion problems in structural dynamics.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

A one-dimensional time-dependent model for studying oxide film growth on metallic surfaces

Corrosion resistance of industrial metals and alloys is largely imparted by the formation of a stable oxide film on the surface that kinetically limits the process of corrosion. Developing a fundamental knowledge of the processes involved in the oxide layer formation and growth is thus important for designing corrosion-resistant alloys. In this work, a model for oxide growth tracking two species (the oxygen and metal vacancies) was presented. The model accounted for species transport and interfacial chemistry. The potential profile was obtained by solving the Poisson equation without needing to invoke the typical linear assumption. Moreover, the fully time-dependent concentration and potential profiles were obtained and the oxide thickness was allowed to evolve by allowing the movement of the domain boundaries at either end. The results were analyzed and the effects of different parameters on the model were discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2018

A Systems Processing Model for Molten Salt Reactors

The chemistry group within the U.S. Department of Energy Molten Salt Reactor Campaign is developing a model for the molten salt reactor system for assessing the suitability of various used fuel processing schemes. In the long run, the developed model is expected to be a tool for better understanding the interplay between the components of the molten salt reactor system including uranium resources, the molten salt reactor, the used fuel processor (e.g., reprocessing plant), and waste management. In this regard, the report proposes a theoretical framework for integrating the various components of the molten salt reactor system. Rather than being buried in implementation details for seemingly heterogeneous components of the reactor system asking for various modeling techniques, this report proposes to approximate and integrate all components with a unified mathematical framework, a discrete-time linear system. An abstracted molten salt reactor system is given as an example to illustrate the effectiveness of the proposed theoretical framework. A fully developed molten salt reactor system model (based on the proposed framework) is expected to be beneficial for various objectives including the transient and equilibrium analysis of fuel cycle scenarios and the identification of potential technology bottlenecks.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
NASA NTRS2007

An Investigation of High-Cycle Fatigue Models for Metallic Structures Exhibiting Snap-Through Response

A study is undertaken to develop a methodology for determining the suitability of various high-cycle fatigue models for metallic structures subjected to combined thermal-acoustic loadings. Two features of this problem differentiate it from the fatigue of structures subject to acoustic loading alone. Potentially large mean stresses associated with the thermally pre- and post-buckled states require models capable of handling those conditions. Snap-through motion between multiple post-buckled equilibrium positions introduces very high alternating stress. The thermal-acoustic time history response of a clamped aluminum beam structure with geometric and material nonlinearities is determined via numerical simulation. A cumulative damage model is employed using a rainflow cycle counting scheme and fatigue estimates are made for 2024-T3 aluminum using various non-zero mean fatigue models, including Walker, Morrow, Morrow with true fracture strength, and MMPDS. A baseline zero-mean model is additionally considered. It is shown that for this material, the Walker model produces the most conservative fatigue estimates when the stress response has a tensile mean introduced by geometric nonlinearity, but remains in the linear elastic range. However, when the loading level is sufficiently high to produce plasticity, the response becomes more fully reversed and the baseline, Morrow, and Morrow with true fracture strength models produce the most conservative fatigue estimates.

Przekop, Adam↗
DOE OSTI2022

A hybrid meshfree discretization to improve the numerical performance of peridynamic models

Efficient and accurate calculation of spatial integrals is of major interest in the numerical implementation of peridynamics (PD). The standard way to perform this calculation is a particle-based approach that discretizes the strong form of the PD governing equation. This approach has rapidly been adopted by the PD community since it offers some advantages. Additionally, it is computationally cheaper than other available schemes, can conveniently handle material separation, and effectively deals with nonlinear PD models. Nevertheless, PD models are still computationally very expensive compared with those based on the classical continuum mechanics theory, particularly for large-scale problems in three dimensions. This results from the nonlocal nature of the PD theory which leads to interactions of each node of a discretized body with multiple surrounding nodes. Here, we propose a new approach to significantly boost the numerical efficiency of PD models. We propose a discretization scheme that employs a simple collocation procedure and is truly meshfree; i.e., it does not depend on any background integration cells. In contrast to the standard scheme, the proposed scheme requires a much smaller set of neighboring nodes (keeping the same physical length scale) to achieve a specific accuracy and is thus computationally more efficient. Our new scheme is applicable to the case of linear PD models and within neighborhoods where the solution can be approximated by smooth basis functions. Therefore, to fully exploit the advantages of both the standard and the proposed schemes, a hybrid discretization is presented that combines both approaches within an adaptive framework. The high performance of the developed framework is illustrated by several numerical examples, including brittle fracture and corrosion problems in two and three dimensions.

42 ENGINEERING↗