Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “linearity”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Linear delta systems with additional degrees of freedom and related methods

A linear delta system includes a frame, rails secured to the frame, linear actuators, each linear actuator coupled to a respective rail and configured to translate along a longitudinal length of the respective rail, pairs of parallel rods each operably coupled to a respective linear actuator, a platform coupled to the pairs of parallel rods, structure configured to movable couple an object to the platform; and at least one degree of freedom imparting assembly including a profiled rod extending in a direction parallel to the rails and a drive unit configured to rotate the profiled rod, wherein the at least one degree of freedom imparting assembly is configured to impart a degree of freedom to the object.

Crawford, Anthony L.↗
DOE OSTI2022

On the effect of scalar flux weighing of linearly anisotropic scattering matrices in few-group transport calculations

The majority of the codes available for homogenized group constant generation for deterministic transport calculations apply the approximation of scalar flux weighting during energy group condensation of higher-order anisotropic scattering matrices. In this paper, we discuss the effect of scalar flux weighting of the linearly anisotropic scattering matrices in the frame of S P{sub 3} and S{sub 12} calculations performed for a two-dimensional VVER-440 reactor benchmark. To compare group constants generated for two neutron energy groups, an infinite pin cell was homogenized with Serpent 2 and ERANOS ECCO. Serpent 2 applies scalar flux, while ERANOS ECCO performs current weighting of the linearly anisotropic scattering matrices during energy group condensation. For analyzing the effect of the various weighting options, three simple reactor models were built assuming different core sizes using standard rectangular assemblies with 15*15 fuel pins. Diffusion, SP{sub 3} and S{sub 12} calculations were performed for the 3 models using group constants generated with Serpent 2 and ERANOS ECCO. The effect of scalar flux weighting of linearly anisotropic scattering matrices is shown by comparing the decrease in reactivity due to the decreased reactor size, as well as the assembly power distribution to reference results obtained with Serpent 2 Monte Carlo calculations. Neglecting higher than linearly anisotropic scattering and indirect application of diffusion coefficients in higher-order transport calculations is advised if angular flux-moment spectra weighted higher-order scattering matrices cannot be generated. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Power Analysis of an ePump Applied to the Linear Functions of an Agricultural Planter

Like many other industries, the agricultural industry has recently ex-perienced pressure to reduce vehicle emissions while improving productivity. Electric actuation is perceived as a viable solution to replace or augment hydrau-lic and mechanical actuation. However, electrification presents challenges with regards to linear functions, where hydraulic actuators have advantages in terms of compactness, tolerance to contamination and resistance to shocks. A combined electro-hydraulic actuation architecture can leverage the benefits of both electric and hydraulic actuation, while reducing the drawbacks of both approaches. This work investigates the potential of a centralized electric driven pump (ePump) system powering the pressure-controlled linear functions of an agricul-tural planter, with the goal of improving the operating point of the main supply pump. In this application, the rotary functions are hydraulically actuated, alt-hough such a solution could be applied also to electric rotary actuation. This is accomplished by setting the ePump to boost the pressure supplied by the tractor to the level required by the linear functions. An accumulator is used to stabilize the flow requirements of the linear functions, and control the pressure supplied to the actuators. Two control schemes are proposed for the regulation of the ac-cumulator pressure, one favoring an efficient operating point for the ePump, the other favoring stable steady state operation. A simulation model of the baseline system and proposed system is developed and validated using experimental data from a full-scale machine. Using the vali-dated simulation, both control architectures are then evaluated for improvement in system power consumption and dynamic requirements on the ePump to assess their effectiveness. Both systems demonstrate significant improvement in power consumption over the baseline system, with the best solution improving effi-ciency by 64%.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2022

Non‐Linear Optical Properties of the ( RE ) 3 CuGeS 7 Family of Compounds

Abstract Non‐linear optical materials must possess a balanced combination of laser‐induced damage threshold (LDT) and second‐harmonic generation (SHG) and be phase matchable. In our previous work, chiral and polar La 3 CuGeS 7 was identified as a promising non‐linear optical material. Herein, we report the optimization of non‐linear optical properties through replacement of La with smaller lanthanides. It is determined that Gd 3 CuGeS 7 exhibits the best combination of SHG (1.6× AgGaS 2 at 88–105 μm particle size) and LDT (3× AgGaS 2 , 89 MW/cm 2 ) and is phase matchable. Based on changes in metal‐sulfur bond lengths and angles, we further propose structural optimization through solid‐solution formation and doping.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

A solution framework for linear PDE-constrained mixed-integer problems

Abstract We present a general numerical solution method for control problems with state variables defined by a linear PDE over a finite set of binary or continuous control variables. We show empirically that a naive approach that applies a numerical discretization scheme to the PDEs to derive constraints for a mixed-integer linear program (MILP) leads to systems that are too large to be solved with state-of-the-art solvers for MILPs, especially if we desire an accurate approximation of the state variables. Our framework comprises two techniques to mitigate the rise of computation times with increasing discretization level: First, the linear system is solved for a basis of the control space in a preprocessing step. Second, certain constraints are just imposed on demand via the IBM ILOG CPLEX feature of a lazy constraint callback. These techniques are compared with an approach where the relations obtained by the discretization of the continuous constraints are directly included in the MILP. We demonstrate our approach on two examples: modeling of the spread of wildfire and the mitigation of water contamination. In both examples the computational results demonstrate that the solution time is significantly reduced by our methods. In particular, the dependence of the computation time on the size of the spatial discretization of the PDE is significantly reduced.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Linear in temperature resistivity in the limit of zero temperature from the time reparameterization soft mode

The most puzzling aspect of the ‘strange metal’ behavior of correlated electron compounds is that the linear in temperature resistivity often extends down to low temperatures, lower than natural microscopic energy scales. We consider recently proposed deconfined critical points (or phases) in models of electrons in large dimension lattices with random nearest-neighbor exchange interactions. The criticality is in the class of Sachdev–Ye–Kitaev models, and exhibits a time reparameterization soft mode representing gravity in dual holographic theories. We compute the low temperature resistivity in a large M limit of models with SU(M) spin symmetry, and find that the dominant temperature dependence arises from this soft mode. The resistivity is linear in temperature down to zero temperature at the critical point, with a co-efficient universally proportional to the product of the residual resistivity and the co-efficient of the linear in temperature specific heat. We argue that the time reparameterization soft mode offers a promising and generic mechanism for resolving the strange metal puzzle.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

The transformation of diamond to graphite: Experiments reveal the presence of an intermediate linear carbon phase

Natural diamonds that have been partially replaced by graphite have been observed to occur in natural rocks. While the graphite-to-diamond phase transition has been extensively studied the opposite of this (diamond to graphite) remains poorly understood. We performed high-pressure and temperature hydrous and anhydrous experiments up to 1.0 GPa and 1300 °C using Amplex premium virgin synthetic diamonds (20–40 μm size) as the starting material mixed with Mg(OH) 2 as a source of H 2 O for the hydrous experiments. The experiments revealed that the diamond-to-graphite transformation at P = 1.0 GPa and T = 1300 °C was triggered by the presence of H 2 O and was accomplished through a three-stage process. Stage 1: diamond reacts with a supercritical H 2 O producing an intermediate 200–500 nm size “globular carbon” phase. This phase is a linear carbon chain; i.e. a polyyne or carbyne. Stage 2: the linear carbon chains are unstable and highly reactive, and they decompose by zigzagging and cross-linking to form sp 2 -hybridized structures. Stage 3: normal, disordered, and onion-like graphite is produced by the decomposition of the sp-hybridized carbon chains which are re-organized into sp 2 bonds. Our experiments show that there is no direct transformation from sp 3 C-bonds into sp 2 C-bonds. Our hydrous high-pressure and high-temperature experiments show that the diamond-to-graphite transformation requires an intermediate metastable phase of a linear hydrocarbon. This process also provides a simple mechanism for the substitution of other elements into the graphite structure (e.g. H, S, O).

36 MATERIALS SCIENCE↗
DOE OSTI2023

GPU-resident sparse direct linear solvers for alternating current optimal power flow analysis

Integrating renewable resources within the transmission grid at a wide scale poses significant challenges for economic dispatch as it requires analysis with more optimization parameters, constraints, and sources of uncertainty. This motivates the investigation of more efficient computational methods, especially those for solving the underlying linear systems, which typically take more than half of the overall computation time. In this paper, we present our work on sparse linear solvers that take advantage of hardware accelerators, such as graphical processing units (GPUs), and improve the overall performance when used within economic dispatch computations. We treat the problems as sparse, which allows for faster execution but also makes the implementation of numerical methods more challenging. We present the first GPU-native sparse direct solver that can execute on both AMD and NVIDIA GPUs. We demonstrate significant performance improvements when using high-performance linear solvers within alternating current optimal power flow (ACOPF) analysis. Furthermore, we demonstrate the feasibility of getting significant performance improvements by executing the entire computation on GPU-based hardware. Finally, we identify outstanding research issues and opportunities for even better utilization of heterogeneous systems, including those equipped with GPUs.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2025

Leveraging explainable AI to characterize floating-point exceptions in linear solvers

Linear solver packages are central to many scientific, engineering, and machine learning applications. When floating-point exceptions occur in these solvers, e.g., division by zero or overflow, numerical results are compromised and become unreliable. Existing static and dynamic analysis tools can detect such exceptions, but they do not explain why the exceptions occur in terms of the solver inputs. Here, we present a study to characterize the inputs that cause numerical exceptions in linear solver packages. Our approach uses explainable AI (XAI) to find the most relevant characteristics of input matrices that explain the occurrence of exceptions in the solvers. Since training data in this domain is scarce, we perform extensive data gathering and data augmentation to obtain exception-inducing inputs. Our approach uses a repair strategy on the features blamed by XAI to validate that such features indeed explain the exceptions. We compare the LIME and SHAP XAI techniques using a dozen matrix features with three classifiers. We evaluate the approach on three widely used linear solver packages and find that some input characteristics can explain the occurrence of exceptions 100% of the time, in specific solvers and preconditioners.

Explainable AI↗
DOE OSTI2025

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↗
DOE OSTI2021

Linearized modeling and optimization of shared mooring systems

Shared mooring systems, where adjacent platforms are tethered directly to each other, can reduce anchor quantities and mooring line lengths in a floating array but also introduce new modeling and design challenges. This paper presents and demonstrates a first-order approach to modeling and designing shared mooring systems that simplifies these challenges. We formulate a general, linearized model for the force-displacement response of shared mooring systems, including inter-platform couplings. Using this linearization, we realize significant simplifications to the shared mooring system design problem and propose a corresponding design optimization approach. Finally, we demonstrate the complete approach on a variety of shared-mooring floating wind farm layouts with generic design assumptions and constraints. Some fundamental observations about shared mooring systems can be made from the results. Polygonal array layouts with perpendicular anchor positions appear to maximize shared mooring system efficiency. Also, arrays with greater levels of sharing tend to exhibit larger deviations in offsets between platforms. Overall, the results show how the linearized approach can be applied to efficiently evaluate and systematically optimize preliminary shared mooring system designs.

17 WIND ENERGY↗
DOE OSTI2021

Distributed and communication-efficient solutions to linear equations with special sparse structure

In this paper we report two distributed and communication-efficient algorithms based on the multi-agent system are proposed to solve a system of linear equations with the Laplacian sparse system matrix. One algorithm is based on the gradient descent method in optimization. In this algorithm, the agents only share partial information instead of all of their collective state vectors to save significant communication. The other algorithm is obtained by approximating Newton’s method for a faster convergence rate. Although it requires twice as much communication as the first one, it is still communication-efficient given the low dimension of the information shared among agents. The convergence at a linear rate is proved for both algorithms, and a comprehensive comparison of their convergence rate, communication burden, and computation costs is also performed. The proposed algorithms can be applied to various systems to solve those problems that can be modeled as a system of linear equations with a Laplacian sparse system matrix. Simulation results with the electric power system illustrate their effectiveness.

42 ENGINEERING↗
DOE OSTI2022

Unsymmetric Pentacene- and Pentacenequinone-Fused Porphyrins: Understanding the Effect of Cross- and Linear-Conjugation

Unsymmetric pentacenequinone-fused (cross-conjugated) and pentacene-fused (linear-conjugated) porphyrins were designed and synthesized. The cross-conjugated (AM 1 –AM 3 ) and linear-conjugated (AM 5 –AM 7 ) porphyrins displayed strikingly different sets of optical and electronic properties, both of which are unusual and nontypical of porphyrins. MCD, DFT, and TDDFT calculations suggest that multiple charge transfer states exist in both π-conjugated systems, which contributes to the complex absorption and MCD spectra of these molecular systems. The general Gouterman’s four-orbital model used to explain porphyrin spectroscopy led to contradicting theoretical and experimental data, and is thus not applicable for these molecular systems. The “2 + 4” and “3 + 3” active spaces have been deduced and have proven effective to interpret the absorption and MCD spectra of the pentacenequinone-fused (cross-conjugated) and pentacene-fused (linear-conjugated) porphyrins, respectively. Spectroelectrochemistry of AM 5 –AM 7 revealed broad and intense IR absorptions in the range of 1500–2500 nm, illustrating the exceptional ability of these pentacene-fused systems to accommodate positive charges. A pronounced metal effect was observed for pentacene-fused porphyrins. While pentacene-fused Ni(II) porphyrin (AM 6 ) demonstrated an abnormal ability to stabilize pentacene with a half-life of >28.3 days, the half-life of the free base and Zn(II) counterparts were normal, similar to those of pentacene analogues. This work provides important and useful information on guiding new material designs.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

X-ray linear dichroic ptychography

Biominerals such as seashells, coral skeletons, bone, and tooth enamel are optically anisotropic crystalline materials with unique nanoscale and microscale organization that translates into exceptional macroscopic mechanical properties, providing inspiration for engineering new and superior biomimetic structures. Using Seriatopora aculeata coral skeleton as a model, here, we experimentally demonstrate X-ray linear dichroic ptychography and map the c-axis orientations of the aragonite (CaCO 3 ) crystals. Linear dichroic phase imaging at the oxygen K-edge energy shows strong polarization-dependent contrast and reveals the presence of both narrow (<35°) and wide (>35°) c-axis angular spread in the coral samples. These X-ray ptychography results are corroborated by four-dimensional (4D) scanning transmission electron microscopy (STEM) on the same samples. Evidence of co-oriented, but disconnected, corallite subdomains indicates jagged crystal boundaries consistent with formation by amorphous nanoparticle attachment. We expect that the combination of X-ray linear dichroic ptychography and 4D STEM could be an important multimodal tool to study nano-crystallites, interfaces, nucleation, and mineral growth of optically anisotropic materials at multiple length scales.

4D scanning transmission electron microscopy↗
DOE OSTI2024

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↗
DOE OSTI2021

A survey of numerical linear algebra methods utilizing mixed-precision arithmetic

The efficient utilization of mixed-precision numerical linear algebra algorithms can offer attractive acceleration to scientific computing applications. Especially with the hardware integration of low-precision special-function units designed for machine learning applications, the traditional numerical algorithms community urgently needs to reconsider the floating point formats used in the distinct operations to efficiently leverage the available compute power. In this study, we provide a comprehensive survey of mixed-precision numerical linear algebra routines, including the underlying concepts, theoretical background, and experimental results for both dense and sparse linear algebra problems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

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↗
DOE OSTI2023

Linear and Nonlinear Optical Properties of Iridium Nanoparticles Grown via Atomic Layer Deposition

Nonlinear optical phenomena enable novel photonic and optoelectronic applications. Especially, metallic nanoparticles and thin films with nonlinear optical properties offer the potential for micro-optical system integration. For this purpose, new nonlinear materials need to be continuously identified, investigated, and utilized for nonlinear optical applications. While noble-metal nanoparticles, nanostructures, and thin films of silver and gold have been widely studied, iridium (Ir) nanoparticles and ultrathin films have not been investigated for nonlinear optical applications yet. Here, we present a combined theoretical and experimental study on the linear and nonlinear optical properties of iridium nanoparticles deposited via atomic layer deposition (ALD). Linear optical constants, such as the effective refractive index and extinction coefficient, were evaluated at different growth stages of nanoparticle formation. Both linear and nonlinear optical properties of these Ir ALD coatings were calculated theoretically using the Maxwell Garnett theory. The third-order susceptibility of iridium nanoparticle samples was experimentally investigated using the z-scan technique. According to the experiment, for an Ir ALD coating with 45 cycles resulting in iridium nanoparticles, the experimentally determined nonlinear third-order susceptibility is about χ Ir (3) = (2.4 – i2.1) × 10 –17 m 2 /V 2 at the fundamental wavelength of 700 nm. The theory fitted to the experimental results predicts a 5 × 10 6 -fold increase around 230 nm. This strong increase is due to the proximity to the Mie resonance of iridium nanoparticles.

36 MATERIALS SCIENCE↗