Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Nonlinear models”

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 577 records · Page 32
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Role of gas composition in weakened nonlinear standing wave excitation and improved plasma radial uniformity in very-high-frequency asymmetric capacitive Ar/CF 4 discharges

Abstract The higher harmonics generated by nonlinear sheath motion would enhance the standing wave effect, and thus lead to center-peaked plasma density profile in very-high-frequency (VHF) capacitive discharges. In this work, a nonlinear transmission line (NTL) model introduced in Zhouet al(2021Plasma Sources Sci. Technol.30125017) has been extended, with radial transport of various particles and nonlinear sheath motion into account, to investigate the effects of CF 4 fraction α on the nonlinear standing wave excitation and plasma radial uniformity in VHF (60 MHz) capacitively coupled Ar/CF 4 plasmas at 3 Pa. The results indicate that for pure Ar discharges (i.e. α = 0%), the nonlinearly excited harmonics with short wavelength significantly enhance the electron power absorption at the radial center, resulting in a pronounced central-high plasma density profile. As α increases, the high-order harmonics are gradually damped due to the increase of resistance, as well as the longer wavelength caused by thicker sheath thickness. Thus, the radial profile of the electron absorbed power density shifts from center-peak to edge-high. Besides, at the radial center, the electron density and Ar + ion density decrease with α , CF 3 + ion density shows an increasing trend, while F − ion density initially rises and then decreases. Moreover, the density profiles of all the species become more uniform at higher CF 4 fraction, due to the suppressed nonlinear standing wave excitation and the longer wavelength of the nonlinear harmonics.

Physics↗
DOE OSTI2021

Current profile and normalized beta control via feedback linearization and Lyapunov techniques

Abstract Simultaneous control of the current profile and normalized plasma beta is an essential control problem in the development of advanced tokamak scenarios. However, this control problem is especially challenging due to the nonlinear nature of the current, heat, and particle transport dynamics, as well as the difficulty to understand and accurately model such processes. In this work, a nonlinear, robust, model-based controller for the simultaneous regulation of the current profile and normalized beta has been designed using feedback linearization and Lyapunov redesign techniques. Feedback linearization avoids approximate linearization of the plasma dynamics, retaining the original physics content of the model. Moreover, the use of Lyapunov redesign techniques makes the controller robust against the uncertainties arising during the modeling process. The controller’s performance in the presence of unknown dynamics is tested in nonlinear, one-dimensional simulations using the Control Oriented Transport SIMulator (COTSIM) code, which employs plasma models that are significantly more complex than those employed for control synthesis.

Physics↗
DOE OSTI2023

False Indications of Dose-Response Nonlinearity in Large Epidemiologic Cancer Radiation Cohort Studies; A Simulation Exercise

This study explores the likely prevalence of false indications of dose-response nonlinearity in large epidemiologic cancer radiation cohort studies (A-bomb survivors, INWORKS, Techa River). Reasons: Increasing numbers of tests of nonlinearity are being made in studies. Hypothesized nonlinear dose-response models have been justified to policy makers by analyses that rely in part on isolated findings that could be statistical fluctuations. After removing dose nonlinearity (linearization) by adjusting person-years of observation at each dose category, indications of nonlinearity, necessarily false, were counted in 5,000 randomized replications of six datasets. The average frequency of any false positive for five indicators of nonlinearity tested against a linear null was roughly 25% in Monte Carlo simulations per study, consistent with binomial calculations, increasing to ~50% within 6 studies assessed. Comparable frequencies were found using Akaike's information criterion (AIC) for model selection or multi-model averaging. False above-zero threshold doses were found more than 50% of the time, averaging to 0.05 Gy, consistent with findings in the 6 studies. Such bias, uncorrected, could distort meta-analyses of multiple studies, because meta-analyses can incorporate high P value findings. AIC-based correction for the extra threshold parameter lowered these false occurrences to 8 to 19%. Given the simulation rates, the possibility of false positives might be noted when isolated findings of nonlinearity are discussed in a regulatory context. When reporting a threshold dose with a P value > 0.05, it would be informative to note the expected high false prevalence rate due to bias.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2021

Machine-Learned Manifold-Based Models for Large Eddy Simulation of Turbulent Combustion

Reduced-order manifold combustion models are commonly used to lower the cost of reacting Large Eddy Simulations (LES) and may be derived either from physical principles as in flamelet models or from data-driven methods like Principal Component Analysis (PCA). In either case, neural networks are increasingly used as part of these models to provide a nonlinear mapping between a small set of pre-defined variables that parameterize the manifold and outputs of interest, such as reaction rates. In this work, we propose a new manifold-based modeling approach that combines the definition of the manifold-parameterizing variables (linear combinations of species), the nonlinear mapping to the outputs, and closure of filtered quantities for LES into the structure of a single neural network. This allows the process used to train the neural network to simultaneously optimize both the functional form of the model and the identities of the inputs to the model. The new approach can flexibly incorporate thermochemical data from any combustion system; if trained on data from 1D flames it can be interpreted as an optimized flamelet model, but it can also be used to learn models from data from more complex configurations. This work presents a priori evaluations of the new approach in both contexts. Evaluation using data from 1D premixed flames demonstrates the physical interpretability of the manifold variables generated by the new approach. Evaluation using data from direct numerical simulations of turbulent flames shows improved predictions relative to either flamelet or PCA-based models in a more complex configuration.

47 OTHER INSTRUMENTATION↗
DOE OSTI2021

Nonlinear elasticity and damping govern ultrafast dynamics in click beetles

Many small animals use springs and latches to overcome the mechanical power output limitations of their muscles. Click beetles use springs and latches to bend their bodies at the thoracic hinge and then unbend extremely quickly, resulting in a clicking motion. When unconstrained, this quick clicking motion results in a jump. While the jumping motion has been studied in depth, the physical mechanisms enabling fast unbending have not. Here, we first identify and quantify the phases of the clicking motion: latching, loading, and energy release. We detail the motion kinematics and investigate the governing dynamics (forces) of the energy release. We use high-speed synchrotron X-ray imaging to observe and analyze the motion of the hinge’s internal structures of four Elater abruptus specimens. We show evidence that soft cuticle in the hinge contributes to the spring mechanism through rapid recoil. Using spectral analysis and nonlinear system identification, we determine the equation of motion and model the beetle as a nonlinear single-degree-of-freedom oscillator. Quadratic damping and snap-through buckling are identified to be the dominant damping and elastic forces, respectively, driving the angular position during the energy release phase. The methods used in this study provide experimental and analytical guidelines for the analysis of extreme motion, starting from motion observation to identifying the forces causing the movement. Here, the tools demonstrated here can be applied to other organisms to enhance our understanding of the energy storage and release strategies small animals use to achieve extreme accelerations repeatedly.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Dark Energy Survey year 3 results: Constraints on cosmological parameters and galaxy-bias models from galaxy clustering and galaxy-galaxy lensing using the redMaGiC sample

We constrain cosmological parameters and galaxy-bias parameters using the combination of galaxy clustering and galaxy-galaxy lensing measurements from the Dark Energy Survey (DES) year-3 data. We describe our modeling framework and choice of scales analyzed, validating their robustness to theoretical uncertainties in small-scale clustering by analyzing simulated data. Using a linear galaxy-bias model and redMaGiC galaxy sample, we obtain 10% constraints on the matter density of the Universe. Here, we also implement a nonlinear galaxy-bias model to probe smaller scales that includes parametrization based on hybrid perturbation theory and find that it leads to a 17% gain in cosmological constraining power. We perform robustness tests of our methodology pipeline and demonstrate stability of the constraints to changes in the theory model. Using the redMaGiC galaxy sample as foreground lens galaxies and adopting the best-fitting cosmological parameters from DES year-1 data, we find the galaxy clustering and galaxy-galaxy lensing measurements to exhibit significant signals akin to decorrelation between galaxies and mass on large scales, which is not expected in any current models. This likely systematic measurement error biases our constraints on galaxy bias and the S 8 parameter. We find that a scale-, redshift- and sky-area-independent phenomenological decorrelation parameter can effectively capture this inconsistency between the galaxy clustering and galaxy-galaxy lensing. We trace the source of this correlation to a color-dependent photometric issue and minimize its impact on our result by changing the selection criteria of redMaGiC galaxies. Using this new sample, our constraints on the S 8 parameter are consistent with previous studies and we find a small shift in the Ω m constraints compared to the fiducial redMaGiC sample. We infer the constraints on the mean host-halo mass of the redMaGiC galaxies in this new sample from the large-scale bias constraints, finding the galaxies occupy halos of mass approximately 1.6 × 10 13 M ⊙ /h.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Comparison of Machine Learning-Based Predictive Models of the Nutrient Loads Delivered from the Mississippi/Atchafalaya River Basin to the Gulf of Mexico

Predicting nutrient loads is essential to understanding and managing one of the environmental issues faced by the northern Gulf of Mexico hypoxic zone, which poses a severe threat to the Gulf’s healthy ecosystem and economy. The development of hypoxia in the Gulf of Mexico is strongly associated with the eutrophication process initiated by excessive nutrient loads. Due to the complexities in the excessive nutrient loads to the Gulf of Mexico, it is challenging to understand and predict the underlying temporal variation of nutrient loads. The study was aimed at identifying an optimal predictive machine learning model to capture and predict nonlinear behavior of the nutrient loads delivered from the Mississippi/Atchafalaya River Basin (MARB) to the Gulf of Mexico. For this purpose, monthly nutrient loads (N and P) in tons were collected from US Geological Survey (USGS) monitoring station 07373420 from 1980 to 2020. Machine learning models—including autoregressive integrated moving average (ARIMA), gaussian process regression (GPR), single-layer multilayer perceptron (MLP), and a long short-term memory (LSTM) with the single hidden layer—were developed to predict the monthly nutrient loads, and model performances were evaluated by standard assessment metrics—Root Mean Square Error (RMSE) and Correlation Coefficient (R). The residuals of predictive models were examined by the Durbin–Watson statistic. The results showed that MLP and LSTM persistently achieved better accuracy in predicting monthly TN and TP loads compared to GPR and ARIMA. In addition, GPR models achieved slightly better test RMSE score than ARIMA models while their correlation coefficients are much lower than ARIMA models. Moreover, MLP performed slightly better than LSTM in predicting monthly TP loads while LSTM slightly outperformed for TN loads. Furthermore, it was found that the optimizer and number of inputs didn’t show effects on the LSTM performance while they exhibited impacts on MLP outcomes. This study explores the capability of machine learning models to accurately predict nonlinearly fluctuating nutrient loads delivered to the Gulf of Mexico. Further efforts focus on improving the accuracy of forecasting using hybrid models which combine several machine learning models with superior predictive performance for nutrient fluxes throughout the MARB.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Hyperelastic nature of the Hoek–Brown criterion

In this article, we propose a nonlinear elasto-plastic model, for which a specific class of hyperbolic elasticity arises as a straight consequence of the yield criterion invariance on the plasticity level. We superimpose this nonlinear elastic (or hyperelastic) behavior with plasticity obeying the associated flow rule. Interestingly, we find that a linear yield criterion on the thermodynamical force associated with plasticity results in a quadratic yield criterion in the stress space. This suggests a specific hyperelastic connection between Mohr–Coulomb and Hoek–Brown (or alternatively between Drucker–Prager and Pan–Hudson) yield criteria. We compare the elasto-plastic responses of standard tests for the Drucker–Prager yield criterion using either linear or the suggested hyperbolic elasticity. Notably, the nonlinear case stands out due to dilatancy saturation observed during cyclic loading in the triaxial compression test. We conclude this study with structural finite element simulations that clearly demonstrate the numerical applicability of the proposed model.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

A Data-Driven Method for Modeling Creep-Fatigue Stress- Strain Behavior Using Neural ODEs

In this paper, we introduce a data-driven machine learning approach for modeling one-dimensional stress–strain behavior under cyclic loading, utilizing experimental data from the nickel-based Alloy 617. The study employs uniaxial creep–fatigue test data acquired under various loading histories and compares two distinct neural network-based ODE models. The first model, known as the black-box model, comprehensively describes the strain–stress relationship using a Neural ODE equation. To interpret this black-box model, we apply the Sparse Identification of Nonlinear Dynamical Systems (SINDy) technique, transforming the black-box model into an equation-based model using symbolic regression. The second model, the Neural flow rule model, incorporates Hooke’s Law for the linear elastic component, with the nonlinear part characterized by a Neural ODE. Both models are trained with experimental data to accurately reflect the observed stress–strain behavior. We conduct a detailed comparison with the standard Chaboche model, which includes three back stresses. Our results demonstrate that the neural network-based ODE models precisely capture the experimental creep–fatigue mechanical behavior, exceeding the standard Chaboche model’s accuracy. Furthermore, an interpretable model derived from the black-box neural ODE model through symbolic regression achieves accuracy comparable to the Chaboche model, enhancing its interpretability. The results highlight the potential of neural network-based ODE models to depict complex creep–fatigue behavior, eliminating the necessity for experts to define a specific, material-focused model form.

creep-fatigue↗
DOE OSTI2021

Pronghorn: A Multidimensional Coarse Mesh Application for Advanced Reactor Thermal-Hydraulics

This paper presents an overview of Pronghorn, a multiscale thermal-hydraulic (T/H) application developed by Idaho National Laboratory and the University of California, Berkeley. Pronghorn, built on the open-source finite element Multiphysics Object-Oriented Simulation Environment (MOOSE), leverages state-of-the-art physical models, numerical methods, and nonlinear solvers to deliver fast-running advanced reactor T/H simulation capabilities within a modern software engineering environment. This work summarizes the physical models, multiphysics and multiscale coupling, and numerical discretization in Pronghorn with emphasis on our initial target application to pebble bed reactors (PBRs). A diverse set of applications are shown to depressurized natural circulation in the SANA experiments, forced convection in the Pebble Bed Modular Reactor, three-dimensional (3-D)/one-dimensional coupling of Pronghorn and RELAP-7 systems T/H for loop analysis in the High Temperature Reactor Power Module, and forced convection in the Mark-1 Pebble Bed Fluoride-Salt-Cooled High-Temperature Reactor. A multiphysics coupling of Pronghorn, RELAP-7, and Griffin deterministic neutronics for a gas-cooled PBR demonstrates the capability of the MOOSE framework for reactor design calculations. These applications highlight the verification and validation underlying Pronghorn’s software development while emphasizing features that improve upon capabilities offered by legacy tools in areas such as 3-D unstructured meshing, physics modeling, and multiphysics coupling.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Towards Smart Grids Enhanced Situation Awareness: A Bi-Level Quasi-Static State Estimation Model

Smart Grid situational awareness is provided by Energy Management Systems. A core process of these systems is State Estimation. The great majority of state estimators model the Smart Grid through a set of nonlinear algebraic equations, named the measurement model. Problem formulation considers the Gauss solution. Several model improvements have been presented regarding the Gauss solution, aiming between others to provide measurement noise robustness to the state estimation process. While considerable effort has been focused on such developments, state estimation is still constrained by the implicit modelling error, and thus inevitably vulnerable to cyber-threats. In this work, a state estimation bi-level formal model is presented towards Smart Grids enhanced situational awareness considering the concepts of synthetic measurements and innovation. Comparative test results with the state-of-the-art on the IEEE 14-bus system are presented highlighting improved situational awareness to bad data. Easy-to-implement model, without hard-to-derive parameters, built-on the classic weighted least squares solution, highlight potential aspects for real-life implementation.

cyber security↗
DOE OSTI2024

Impact of K-factor on Short Circuit Program Convergence for Inverter-Based Resources during Faults

Classical short-circuit programs that linearize the power network are no longer applicable for inverter based resources (IBRs), necessitating an iterative approach. Phasor domain programs can model the IBRs using an iterative approach considering nonlinear fault responses. In phasor domain models, the IBR can be modeled as a voltage controlled current source (VCCS) in tabular form with positive, negative, and zero sequence information for balanced and unbalance short-circuit faults. In the VCCS modeling of the IBR, positive and negative incremental reactive current, also known as the k-factor, plays an important role in short circuit program convergence. In this work, a few approaches: conventional VCCS modeling with a k-factor of 2, conventional VCCS modeling with a k-factor of 2 with modified pre-fault voltages, VCCS characteristics based on the power flow solution with a k-factor of 2, and VCCS characteristics based on the power flow solution with a k-factor of 1 are investigated for short circuit program convergence under higher IBR penetration. The IEEE 39 bus New England Test System is taken as the test system, and simulations are carried out in PSS®CAPE 15.0.26 simulation software. Simulation results demonstrate that IBR penetration is higher for the VCCS model, which corresponds to the power flow solution with k-factor 1, compared to other approaches.

42 ENGINEERING↗
DOE OSTI2023

Parameter identification for symmetrical Prandtl-Ishlinskii hysteresis model using Gauss-Newton method

The Prandtl-Ishlinskii (PI) hysteresis model has been used to describe nonlinear hysteretic behaviors. In this paper, the PI model is first simplified to accommodate the symmetrical behavior of ferromagnetic hysteresis. A parameter identification method is then developed to extract the parameters of the symmetrical PI model from experimental data. Based on the Gauss-Newton method, regularization technique, and line search method, the model parameters can be retrieved efficiently from the measured data. Numerical examples are presented to demonstrate the performance of the proposed method.

Yan, Su↗
DOE OSTI2025

The role of fast and slow dynamics in nonlinear resonant ultrasound spectroscopy of consolidated granular materials

Abstract Elastic nonlinearity observed in consolidated granular media can be attributed to the combination of slow and fast effects, which give rise to hysteresis and relaxation of both modulus and damping after the sample is perturbed. A consequence is a high level of complexity in the measurements of the sample linear and nonlinear elastic parameters. The results of experiments are dependent on the experimental protocol that is adopted to measure the relevant quantities and it is hard to quantify parameters with accuracy and repeatability. Here we focus on examining Nonlinear Resonant Ultrasound Spectroscopy, showing experimentally the role of slow dynamics in the process and quantifying/discussing its influence on the quantification of nonlinearity. We also propose a model to describe the process, which shows that different contributions to nonlinearity (e.g., classical and hysteretic) could be due to physical features (defects) relaxing with different relaxation times.

Science & Technology - Other Topics↗
DOE OSTI2022

OC6 Phase Ia Definition Document: Validation of Nonlinear Hydrodynamic Loading on the DeepCwind Semisubmersible

The objective of OC6 Phase I is to investigate the persistent under-prediction (about 20% on average) of the structural loads in the OC5-DeepCwind floating semisubmersible wind system, as observed within the previous OC5 Phase II project. The OC5 results indicate that much of this under-prediction originates from the low-frequency response of the system at the surge/pitch natural frequencies, resulting from nonlinear hydrodynamic loading and/or overprediction of hydrodynamic damping. This work package therefore focuses on better understanding this low-frequency response behavior, the applicability of computational-fluid-dynamics (CFD) models and engineering-level hydrodynamic models for predicting the nonlinear hydrodynamic loading, and the level of uncertainty in the measured response characteristics. To address these objectives, two experimental campaigns were developed and performed in 2017 and 2018 by a sub-group of the OC5 project: Floating test campaign: A simplified OC5-DeepCwind semisubmersible model (no turbine, rigid tower) was moored to three taut-spring lines and subjected to motion decay and wave tests (regular and irregular waves), during which motion amplitudes of the model were measured (see Figure 1a). The purpose of this test campaign was to examine the repeatability of the low-frequency response of the system to wave excitation and determine the level of uncertainty in the response behavior.Constrained test campaign: The same model, but without the tower present, was fixed to a carriage and subjected to towing, forced oscillation and wave tests (regular and irregular waves), during which forces exerted upon the model were measured (see Figure 1b). The purpose of this test campaign was to examine the hydrodynamic loading contributions by separating the wave excitation and radiation forces. Both model test campaigns were performed in the concept basin of the Maritime Research Institute Netherlands (MARIN) in Wageningen, the Netherlands and the first test was within the framework of the MARINET2 project. Both experimental campaigns were examined within Phase Ia of the OC6 project. The purpose of this report is to supply the needed information to participants of the OC6 project to build and simulate the experimental campaigns.

17 WIND ENERGY↗
DOE OSTI2025

Modeling kinetic effects of charged vacancies on electromechanical responses of ferroelectrics: Rayleighian approach

Understanding the time-dependent effects of charged vacancies on the electromechanical responses of materials is at the forefront of research for designing materials exhibiting metal-insulator transitions and memristive behavior. A Rayleighian approach is used to develop a model for studying the nonlinear kinetics of the reaction leading to generation of vacancies and electrons via the dissociation of vacancy-electron pairs. Also, diffusion and elastic effects of charged vacancies are considered to model polarization-electric potential and strain-electric potential hysteresis loops. The model captures multiphysics phenomena by introducing couplings among polarization, the electric potential, stress, strain, and concentrations of charged (multivalent) vacancies and electrons (treated as classical negatively charged particles), where the concentrations can vary due to association-dissociation reactions. A derivation of coupled time-dependent equations based on the Rayleighian approach is presented. Three limiting cases of the governing equations are considered, highlighting the effects of (1) nonlinear reaction kinetics on the generation of charged vacancies and electrons, (2) Vegard's law (i.e., the concentration-dependent local strain) on asymmetric strain-electric potential relations, and (3) coupling between a fast component and the slow component of the net polarization on the polarization-electric-field relations. The Rayleighian approach discussed in this work should pave the way for developing a multiscale modeling framework in a thermodynamically consistent manner while capturing multiphysics phenomena in ferroelectric materials. Published by the American Physical Society 2025

Kumar, Rajeev (ORCID:0000000194943488)↗
DOE OSTI2025

Physics-Informed Machine Learning Model for Ceramic Matrix Composite Creep

A physics-informed recurrent neural network (RNN) based surrogate model is developed to emulate the nonlinear, time-dependent constitutive behavior of ceramic matrix composites (CMCs) driven by matrix damage and constituent creep at the microscale. Physics-informed constraints are introduced into the surrogate model through regularization to ground the prediction in physics and improve its predictive capabilities. Training data is generated using the high-fidelity generalized method of cells (HFGMC) approach which calls appropriate creep and damage models for each of the constituents. This coupling permits simulating the nonlinear behavior of CMCs based on constituent response at the microscale along with microstructural features such as fiber and porosity volume fraction and fiber radius. The microscale repeating unit cell is loaded under creep fatigue conditions to replicate the material loading experienced in a turbine engine. Therefore, the RNN-based surrogate model is tasked with predicting, as a function of variable input stress sequence, temperature, and microstructural features, the resulting strain history response while satisfying physical constraints related to creep rate, isochoric inelastic deformation, and strain energy density. The trained surrogate model is shown to effectively match the strain history over quantified distributions of microstructural features and relevant loading regimes and temperatures. Neural network based surrogate models can offer efficient alternatives to running computationally intensive multiscale material models to simulate the nonlinear response of large structural models. Therefore, the presented work provides evidence towards the feasibility of developing, training, and running such models for CMCs with complex microstructures, nonlinear time-dependent material response, and under non-monotonic loading conditions.

ceramic matrix composites↗