Search NASA⌕ Search

SEARCH · Search NASA

Results for “Linear”

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 289 records · Page 16

Simplifying activations with linear approximations in neural networks

A key step in Neural Networks is activation. Among the different types of activation functions, sigmoid, tanh, and others involve the usage of exponents for calculation. From a hardware perspective, exponential implementation implies the usage of Taylor series or repeated methods involving many addition, multiplication, and division steps, and as a result are power-hungry and consume many clock cycles. We implement a piecewise linear approximation of the sigmoid function as a replacement for standard sigmoid activation libraries. This approach provides a practical alternative by leveraging piecewise segmentation, which simplifies hardware implementation and improves computational efficiency. In this paper, we detail piecewise functions that can be implemented using linear approximations and their implications for overall model accuracy and performance gain. Our results show that for the DenseNet, ResNet, and GoogLeNet architectures, the piecewise linear approximation of the sigmoid function provides faster execution times compared to the standard TensorFlow sigmoid implementation while maintaining comparable accuracy. Specifically, for MNIST with DenseNet, accuracy reaches 99.91% (Piecewise) vs. 99.97% (Base) with up to 1.31x speedup in execution time. For CIFAR-10 with DenseNet, accuracy improves to 98.97% (Piecewise) vs. 99.40% (Base) while achieving 1.24x faster execution. Similarly, for CIFAR-100 with DenseNet, the accuracy is 97.93% (Piecewise) vs. 98.39% (Base), with a 1.18x execution time reduction. These results confirm the proposed method’s capability to efficiently process large-scale datasets and computationally demanding tasks, offering a practical means to accelerate deep learning models, including LSTMs, without compromising accuracy.

Activation function↗

Radiation induced non-linear oscillations in ITER baseline scenario plasmas in DIII-D

Abstract This work shows how the radiation brought about by metals or metal-equivalent radiators such as Kr and Xe produces non-linear dynamics on otherwise stationary β N flattops of DIII-D ITER Baseline Scenario demonstration discharges. The Kr and Xe gases are used to reproduce the radiative loss rates of W in present machines that operate at core temperatures much lower than the expected ITER temperature. Experiments on DIII-D with injection of Kr and Xe, as well as with sources of intrinsic metals reach the range of radiated fraction values expected in the ITER core and experience slow oscillations in temperature and radiated power. In many cases of high radiated fraction, the core temperature decreases enough for the safety factor profile to rise above the 1/1 rational surface, naturally eliminating sawteeth and occasionally producing a persistent helical core. The oscillations can be reproduced by a modified Lotka–Volterra system for temperature and radiated fraction if diffusion and noise are included, which indicates that the interplay between temperature and radiation can be the main cause of the cyclic nature of the system. A new physics based model which includes equations for temperature, density and input power can also reproduce the oscillations observed in the experiments. The present results suggest that the non-linearity of the system can be increased by the inclusion of the inherently non-linear alpha heating term, which is proportional to ∼ n e 2 T i 2 , and obtains oscillations in the model when added to an otherwise more stationary system.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Mountain Basin Controls on the Snow-to-Streamflow Signal: An AIC-Weighted Multiple Linear Regression Framework

A regression-based analysis quantifies how basin characteristics modulate the snow-to-streamflow signal. First, we use the ERA5-Land reanalysis gridded product (European Centre for Medium Range Weather Forecasts reanalysis 5 -Land component) for 4,655 hydrologic unit code - 10 (HUC10) mountain basins across the western United States (US) for water years 1987–2024. Linear regressions are performed for peak snow water equivalent (SWE) and annual streamflow for each mountain basin. Models use ordinary least squares in Python’s statsmodels package. After which, an Akaike Information Criterion (AIC)–weighted ensemble multiple linear regression (MLR) framework with 47 watershed traits is used to predict the linear regression coefficient of determination (r-squared) defining the ability of peak SWE to predict annual streamflow across all mountain basin. Predictor sets are constrained to avoid multicollinearity by excluding models with variance inflation factors (VIF) greater than 5. Mountain basin traits included in the MLR include seasonal climate, topography, vegetation type and structure, and bedrock geology. Accepted models are considered if their AIC is within 2.0 of the model with the minimum AIC, or best model. To compare predictor influence across acceptable models, we computed standardized regression coefficients. To evaluate structural redundancy among models, we constructed binary inclusion vectors for each acceptable model, denoting whether a predictor was present (1) or absent (0). Core predictor variables are defined as occurring in at least 67% of the acceptable models. For this regional analysis, only one model was found acceptable, with higher snow-to-streamflow translation (higher r-squared) occurring in colder mountain basins with higher relative winter precipitation, more snow accumulation and a lower fraction of annual precipitation that falls in the spring and summer. The second component of the data package uses previously published, high-resolution output from an integrated hydrological model of the East River watershed using the U.S. Geological Survey Groundwater and Surface water Flow model (GSFLOW, doi:10.15485/1998576). East River MLR expands upon the approach described above to explore the response of five streamflow metrics—annual streamflow, runoff efficiency, 7-day minimum flow, low-flow duration, and non-perennial stream fraction to snow system indicators including peak SWE, snow-covered area, snow disappearance date, and the fraction of basin area characterized by low-to-no snow, as well as seasonal precipitation and temperature, and annual hydrologic variables representing soil moisture, evapotranspiration (ET), the partitioning of incoming precipitation to evapotranspiration (ET/P), groundwater storage, and groundwater inflow to streams. MLR was done on all water years (P0: 1987-2024) and for each period as determined in the split analysis using pooled regression techniques (P1: 1987-2011 and P2: 2012-2024) to evaluate shifting predictor variable emphasis on streamflow generation. Results indicate that since 2012, peak SWE has lost statistical strength in its prediction of annual streamflow and runoff efficiency, and the indirect influence of spring temperature has emerged as critically important. Low-flow metrics remain largely influenced by soil moisture, vegetation water use and groundwater inflows with summer precipitation becoming a direct influence on minimum summer flow. Together, these data and Python-based analysis tools provide a framework for identifying the key watershed characteristics that control how streamflow responds to snow from year to year. The package also helps quantify uncertainty in statistical models and assess how snow–streamflow relationships vary across regions and over time. This dataset contains comma-separated values files (.csv), text files (.txt), python code files (.py), figure files (.png), and shapefiles (.cpg, .dbf, .prj, .sbn, .sbx, .shp, .xml). Further details on file contents and MLR execution can be found in the readme file and the FLMD files. Work was supported by the Watershed Function Science Focus Area at Lawrence Berkeley National Laboratory funded by the US Department of Energy, Office of Science, Biological and Environmental Research under Contract No. DE-AC02-05CH11231.

54 ENVIRONMENTAL SCIENCES↗

Linear and nonlinear ARMA model parameter estimation using an artificial neural network

This paper addresses parametric system identification of linear and nonlinear dynamic systems by analysis of the input and output signals. Specifically, we investigate the relationship between estimation of the system using a feedforward neural network model and estimation of the system by use of linear and nonlinear autoregressive moving-average (ARMA) models. By utilizing a neural network model incorporating a polynomial activation function, we show the equivalence of the artificial neural network to the linear and nonlinear ARMA models. We compare the parameterization of the estimated system using the neural network and ARMA approaches by utilizing data generated by means of computer simulations. Specifically, we show that the parameters of a simulated ARMA system can be obtained from the neural network analysis of the simulated data or by conventional least squares ARMA analysis. The feasibility of applying neural networks with polynomial activation functions to the analysis of experimental data is explored by application to measurements of heart rate (HR) and instantaneous lung volume (ILV) fluctuations.

NASA Discipline Regulatory Physiology↗

A Linear Programming Approach to Routing Control in Networks of Constrained Nonlinear Positive Systems with Concave Flow Rates

We consider control design for positive compartmental systems in which each compartment's outflow rate is described by a concave function of the amount of material in the compartment.We address the problem of determining the routing of material between compartments to satisfy time-varying state constraints while ensuring that material reaches its intended destination over a finite time horizon. We give sufficient conditions for the existence of a time-varying state-dependent routing strategy which ensures that the closed-loop system satisfies basic network properties of positivity, conservation and interconnection while ensuring that capacity constraints are satisfied, when possible, or adjusted if a solution cannot be found. These conditions are formulated as a linear programming problem. Instances of this linear programming problem can be solved iteratively to generate a solution to the finite horizon routing problem. Results are given for the application of this control design method to an example problem. Key words: linear programming; control of networks; positive systems; controller constraints and structure.

Positive Systems↗

Recent Updates to the GEOS-5 Linear Model

Global Modeling and Assimilation Office (GMAO) is close to having a working 4DVAR system and has developed a linearized version of GEOS-5.This talk outlines a series of improvements made to the linearized dynamics, physics and trajectory.Of particular interest is the development of linearized cloud microphysics, which provides the framework for 'all-sky' data assimilation.

Adjoint↗

A Linear-Elasticity Solver for Higher-Order Space-Time Mesh Deformation

A linear-elasticity approach is presented for the generation of meshes appropriate for a higher-order space-time discontinuous finite-element method. The equations of linear-elasticity are discretized using a higher-order, spatially-continuous, finite-element method. Given an initial finite-element mesh, and a specified boundary displacement, we solve for the mesh displacements to obtain a higher-order curvilinear mesh. Alternatively, for moving-domain problems we use the linear-elasticity approach to solve for a temporally discontinuous mesh velocity on each time-slab and recover a continuous mesh deformation by integrating the velocity. The applicability of this methodology is presented for several benchmark test cases.

higher-order↗

Linear Mixed-Effects Models for Human-in-the-Loop Tracking Experiment Data

Linear mixed-effects models provide several benefits over more traditional statistical inference tests that are particularly useful for most human-in-the-loop tracking experiment data. However, surprisingly, mixed models are virtually not used for the analysis of tracking experiment data. This paper uses of linear mixed-effects models to analyze combined tracking data from two previous human-in-the-loop roll tracking experiments that compared control behavior metrics collected in both a research aircraft and a motion-base simulator. In the experiments, pilots' behavior under 10 different motion configurations with varying motion filter gains and break frequencies was evaluated and compared to that in the real aircraft. The linear mixed-effects model analysis on the combined dataset confirmed the main statistical outcomes of the individual experiments. The main benefits of mixed models for this type of data were demonstrated by successfully combining data from two experiments that used different experimental conditions and of which one had an additional apparatus and the other a missing participant. Finally, the mixed-model analysis was able to explicitly test for scientifically relevant statistical differences in the dependent measures between the aircraft and simulator, as well as both experiments.

manual control↗

Linear Wave Motion from Concentrated Harmonic Sources in Blasius Flow

The motion of individual linear instability waves in shear flows is well described by existing theoretical and numerical methods. However, naturally occuring sources produce coherent wave motions with broadband spanwise wavenumber and frequency spectra, and the different spectral components interact both linearly and nonlinearily. This paper describes a series of calculations for the parameters of three different experiments using locally parallel linear stability theory (LST), the Parabolized Stability Equations (PSE), and Direct Numerical Simulation (DNS). The calculations illustrate the strengths and weaknesses of the different methods, the extent to which the methods agree or disagree, and, finally the extent to which agreement with the measurements can be attained, given that the experiments also have their own difficulties.

linear↗

Photoneutron Production Using an Electron Linear Accelerator for Applications in Neutron Imaging

Photoneutron production is possible using an electron linear accelerator and a target capable of generating photonuclear reactions. A short pulse neutron source can be useful for neutron imaging dynamic experiments. This study is aimed at the feasibility of photoneutron production using a 20 MeV electron linear accelerator and tungsten and depleted uranium targets of various thicknesses. MCNP6 (Monte Carlo N-Particle) code will be used to develop a computational model to estimate total neutron yield, and this will be verified at the Idaho State University’s Accelerator Center. After verification of the neutron yield and energy spectra, an additional MCNP6 model will be developed to analyze the neutron imaging processes. This study will potentially prove it is possible to conduct multi-mode imaging experiments on the anticipated Scorpius electron linear accelerator at the Nevada National Security Sites.

43 PARTICLE ACCELERATORS↗

Experimental validation of a Kalman observer using linearized OpenFAST and a fully instrumented 1:70 model

Abstract To enable real‐time monitoring and control strategies for floating offshore wind turbines, accurate information about the state of the system is needed. This paper details the application of a Kalman filter to the UMaine VolturnUS‐S floating wind platform to provide accurate state estimates in real time using minimal system measurements. The midfidelity nonlinear simulation tool OpenFAST was used to generate the underlying linear state‐space model for the Kalman filter. This linear model and its limitations are demonstrated through comparison with experimental data collected on a 1:70 froude‐scaled model of the floating platform and tower. Using a selection of five measurements from the real system, a Kalman filter was developed to provide estimates for the remaining system states and measurements. These estimates were then validated against the experimental values collected from testing of the scale model. Validation of the Kalman filter produced accurate estimates of surge, heave, and tower base bending moment, measurements of which were not available to the Kalman filter. Performance of the Kalman filter was tested and validated over a range of sea conditions from rated wind speed to storm events and demonstrated robustness in the Kalman filter to maintain accuracy across all operating conditions despite significant error in the underlying linear model for extreme conditions.

17 WIND ENERGY↗

Linear Reciprocating Tribometer for In Situ Neutron Reflectometry of Soft Matter

Neutron reflectometry is a technique for measuring structure near planar interfaces that has been previously used to non-destructively characterize the polymer density of hydrated, dilute, and soft materials. Previous investigations have conducted neutron reflectometry measurements of liquids, gels, emulsion, and polymer solutions at rest, in compression, and subject to shear stress. However, correlating structure with tribological properties of soft materials presents significant experimental challenges for prior instruments due to wall slip, sample thickness, and structural heterogeneity (e.g., depth-wise gradients). A linear reciprocating tribometer offers several advantages for in situ neutron reflectometry studies, including uniform velocity profiles, constant shear stress over large regions of interest, and independent control of normal force and sliding velocity during measurements. This work outlines basic considerations for the design of a custom linear reciprocating tribometer that operates in a neutron beamline and includes commissioning measurements. The tribometer is designed to compress soft and hydrated materials against linearly reciprocating silicon disks. The three key design considerations for this tribometer are (1) safety, (2) neutron transmission, and (3) sample positioning. This instrument design will enable in situ studies of soft matter and illuminate the role of interfacial structure on tribological phenomena.

42 ENGINEERING↗

Predicting non-linear stress–strain response of mesostructured cellular materials using supervised autoencoder

Recent breakthroughs in advanced manufacturing capabilities have made it possible to design and print sophisticated topologies of cellular structures using diverse engineering materials such as metals, polymers, and ceramics. In these architectured materials, it is often desirable to tailor the mechanical properties by altering the unit cell topology. This necessitates an in-depth understanding of how the topology of the unit cell structure affects the macroscopic behavior of the material in both the linear and the non-linear regimes encountered under large compression. Here, we have developed a machine learning (ML) approach capable of accelerating the prediction of the stress–strain response of a polymer-based cellular structure under uniaxial confined compression. As part of generating the training data for ML, 60,000 mesostructures were generated using a relatively novel approach based on cellular automata, and their corresponding stress–strain responses were obtained from the finite element simulations. Principal component analysis (PCA) was used to reduce the dimensionality of the stress–strain curves. With only 20 principal components, PCA captured 99.89% of the variance in the stress–strain curves while reducing the dimensionality by 5X. ML using supervised autoencoder was able to successfully speed up the prediction of the non-linear stress–strain response of a unit cell by up to 4600X. The proposed method can serve as an efficient data generation tool and a rapid means for predicting the structure–property relationship through accelerated forward modeling of cellular materials under compaction, in cases where the macroscopic stress–strain response is governed by the unit-cell topology.

36 MATERIALS SCIENCE↗

An efficient explicit implementation of a near-optimal quantum algorithm for simulating linear dissipative differential equations

We propose an efficient block-encoding technique for the implementation of the Linear Combination of Hamiltonian Simulations (LCHS) for simulating dissipative initial-value problems. This algorithm approximates a target nonunitary operator as a weighted sum of Hamiltonian evolutions, thereby emulating a dissipative problem by mixing various time scales. We introduce an efficient encoding of the LCHS into a quantum circuit based on a simple coordinate transformation that turns the dependence on the summation index into a trigonometric function. Classically, this method is equivalent to the use of a highly accurate Fejér-Clenshaw-Curtis quadrature formula. Quantumly, this significantly simplifies block-encoding of a dissipative problem and allows one to perform an exponential number of Hamiltonian simulations by a single Quantum Signal Processing (QSP) circuit. The resulting LCHS circuit has high success probability and the selector scales logarithmically with the number of terms in the LCHS sum and linearly with time. Careful analysis of error convergence proves that this method is more efficient than other LCHS circuits that have recently appeared in the literature. We verify the quantum circuit and its scaling by simulating it on a digital emulator of fault-tolerant quantum computers and, as a test problem, solve the advection-diffusion equation. The proposed algorithm can be used for simulating a wide class of nonunitary initial-value problems including the Liouville equation with added dissipation and linear embeddings of nonlinear systems, such as the Koopman-von Neumann and Carleman embeddings.

Novikau, I [Lawrence Livermore National Laboratory↗

Climate-eutrophication-anoxia interactions in Late Glacial Soppensee, Switzerland: Forcings, non-linear responses and recovery

Combined effects of climate warming and anthropogenic nutrient loadings lead to lake eutrophication and anoxia globally. Because of chemical feedbacks, lakes under multiple stressors often respond in non-linear ways. However, it remains unclear whether climate change alone can lead to non-linear lake responses in the absence of anthropogenic nutrient disturbances. Here, we investigate the interactions between climate variability, nutrient cycling and trophic state changes, mixing regimes, anoxia and related chemical feedback in a small kettle-hole lake in Switzerland during Late Glacial times (15.2–12.6 cal ka BP), a period known for high-amplitude climate change in pre-anthropogenic times. After its formation during Heinrich Stadial 1 (>15 cal ka BP), Soppensee was oligotrophic and well-mixed. Soppensee became eutrophic and developed anoxia at 14.25 cal ka BP. Phosphorus (P) was released from sediments through the reductive dissolution of Fe-oxyhydroxides, fuelling eutrophication. Eutrophication lagged the Bølling warming (14.65 cal ka BP) by 400 years, suggesting that rising temperatures were not the trigger for eutrophication. Instead, eutrophication responded non-linearly to forest closure (threshold at 76 % arboreal pollen AP), which shielded Soppensee from wind mixing, enhancing lake stratification, anoxia and P release, intensifying eutrophication. These conditions ended during the 200-years cold period of the Aegelsee Oscillation (GI-1d, ca. 14.0 cal ka BP) when the landscape regionally opened (AP<76 %); the lake became well-mixed, oxygenated and P was efficiently sequestered. Throughout the Allerød (13.9–12.8 cal ka BP), enhanced Fe input prompted diagenetic vivianite formation, sequestering P in sediments, naturally remediating lake eutrophication despite closed forests, warm temperatures, lake stratification and anoxia.

Environmental sciences↗

A quasi-linear model of electromagnetic turbulent transport and its application to flux-driven transport predictions for STEP

A quasi-linear reduced transport model is developed from a database of high-β electromagnetic nonlinear gyrokinetic simulations performed with spherical tokamak for energy production (STEP) relevant parameters. The quasi-linear model is fully electromagnetic and accounts for the effect of equilibrium flow shear using a novel approach. Its flux predictions are shown to agree quantitatively with predictions from local nonlinear gyrokinetic simulations across a broad range of STEP-relevant local equilibria. This reduced transport model is implemented in the T3D transport solver that is used to perform the first flux-driven simulations for STEP to account for transport from hybrid kinetic ballooning mode turbulence, which dominates over a wide region of the core plasma. Nonlinear gyrokinetic simulations of the final transport steady state from T3D return turbulent fluxes that are consistent with the reduced model, indicating that the quasi-linear model may also be appropriate for describing the transport steady state. Within the assumption considered here, our simulations support the existence of a transport steady state in STEP with a fusion power comparable to that in the burning flat top of the conceptual design, but do not demonstrate how this state can be accessed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Linear Discriminant Analysis-Based Machine Learning and All-Atom Molecular Dynamics Simulations for Probing Electro-Osmotic Transport in Cationic-Polyelectrolyte-Brush-Grafted Nanochannels

Deciphering the correct mechanisms governing certain phenomena in polyelectrolyte (PE) brush grafted systems, revealed through atomistic simulations, is an extremely challenging problem. In a recent study, our all-atom molecular dynamics (MD) simulations revealed a non-linearly large electroosmotic (EOS) flow (in the presence of an applied electric field) in nanochannels grafted with PMETAC [Poly(2-(methacryloyloxy)ethyl trimethylammonium chloride] brushes. Given the lack of any formal procedure that would have directed us to identify the correct factors responsible for such an occurrence, we needed to spend several months and devote significant analyses to unravel the involved mechanisms. In this paper, we propose a Linear Discriminant Analysis (LDA) based Machine Learning (ML) approach to address this gap. At first, we obtain data on certain basic features from the all-atom MD data. These basic features represent the number of atoms of certain species around one atom of another (or same) species. Here, we obtain such data on basic features for a reference case (case of an EOS flow in PMETAC-brush-grafted nanochannels with a smaller electric field) and a perturbed case (case of an EOS flow in PMETAC-brush-grafted nanochannels with a larger electric field) in bins in which the nanochannel half height has been divided into. These datasets are high-dimensional dataset, to which the LDA is applied. This leads to the projection of the data (between the reference and the perturbed states) in a highly separated form on a 1D line. From such LDA calculations, we are able to identify the relative importance of the different basic features in ensuring this separation of the data (between the reference and the perturbed states) on the 1D line. This relative importance of the different basic features is quantified as “importance scores” for the different features, which in turn tell us what to study and where to study. Such knowledge enables us to rapidly identify the key factors responsible for the non-linearly large EOS transport in PMETAC-brush-grafted nanochannels.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Hydroesterification of Polycyclooctene to Access Linear Ethylene Ethyl Acrylate Copolymers as a Step Toward Polyolefin Functionalization

To advance a strategy of polymer-to-polymer upcycling of waste polyolefin by dehydrogenation then functionalization, we report successful hydroesterification of polycyclooctene (PCOE), an analogue for partially unsaturated polyethylene. Here, we convert PCOE to a linear analog for poly(ethylene-co-ethyl acrylate) (EEA) across a range of ethyl acrylate incorporations (0 to 18 mol % of ethylene units). The ester incorporation was well controlled by reaction time, and the remaining C=C bonds were subsequently hydrogenated. Here, the bulky ethyl acrylate groups did not incorporate into orthorhombic PE crystals, decreasing the crystallinity, crystallite size, and melting temperature with increasing functionalization. Additionally, hydroesterification tuned the dynamic mechanical properties, decreasing both the glass transition temperature and the storage modulus in the rubbery regime with greater functionalization. The linear EEA analogs reported here achieve remarkable extensibility (strain > 4000%) and high toughness, comparable to commercial random and branched EEA. Ultimately, we demonstrate successful conversion of an analogue to dehydrogenated PE to a linear EEA with favorable mechanical properties.

36 MATERIALS SCIENCE↗