Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “nonlinear analysis”

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

Dakota, A Multilevel Parallel Object-Oriented Framework for Design Optimization, Parameter Estimation, Uncertainty Quantification, and Sensitivity Analysis: Theory Manual (V.6.15)

The Dakota toolkit provides a flexible and extensible interface between simulation codes and iterative analysis methods. Dakota contains algorithms for optimization with gradient and nongradient-based methods; uncertainty quantification with sampling, reliability, and stochastic expansion methods; parameter estimation with nonlinear least squares methods; and sensitivity/variance analysis with design of experiments and parameter study methods. These capabilities may be used on their own or as components within advanced strategies such as surrogate-based optimization, mixed integer nonlinear programming, or optimization under uncertainty. By employing object-oriented design to implement abstractions of the key components required for iterative systems analyses, the Dakota toolkit provides a flexible and extensible problem-solving environment for design and performance analysis of computational models on high performance computers. This report serves as a theoretical manual for selected algorithms implemented within the Dakota software. It is not intended as a comprehensive theoretical treatment, since a number of existing texts cover general optimization theory, statistical analysis, and other introductory topics. Rather, this manual is intended to summarize a set of Dakota-related research publications in the areas of surrogate-based optimization, uncertainty quantification, and optimization under uncertainty that provide the foundation for many of Dakota's iterative analysis capabilities.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Dakota, A Multilevel Parallel Object-Oriented Framework for Design Optimization, Parameter Estimation, Uncertainty Quantification, and Sensitivity Analysis (V.6.16 User's Manual)

The Dakota toolkit provides a flexible and extensible interface between simulation codes and iterative analysis methods. Dakota contains algorithms for optimization with gradient and nongradient-based methods; uncertainty quantification with sampling, reliability, and stochastic expansion methods; parameter estimation with nonlinear least squares methods; and sensitivity/variance analysis with design of experiments and parameter study methods. These capabilities may be used on their own or as components within advanced strategies such as surrogate-based optimization, mixed integer nonlinear programming, or optimization under uncertainty. By employing object-oriented design to implement abstractions of the key components required for iterative systems analyses, the Dakota toolkit provides a flexible and extensible problem-solving environment for design and performance analysis of computational models on high performance computers. This report serves as a user's manual for the Dakota software and provides capability overviews and procedures for software execution, as well as a variety of example studies.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Dakota, A Multilevel Parallel Object-Oriented Framework for Design Optimization, Parameter Estimation, Uncertainty Quantification, and Sensitivity Analysis: Version 6.16 Theory Manual

The Dakota toolkit provides a flexible and extensible interface between simulation codes and iterative analysis methods. Dakota contains algorithms for optimization with gradient and nongradient-based methods; uncertainty quantification with sampling, reliability, and stochastic expansion methods; parameter estimation with nonlinear least squares methods; and sensitivity/variance analysis with design of experiments and parameter study methods. These capabilities may be used on their own or as components within advanced strategies such as surrogate-based optimization, mixed integer nonlinear programming, or optimization under uncertainty. By employing object-oriented design to implement abstractions of the key components required for iterative systems analyses, the Dakota toolkit provides a flexible and extensible problem-solving environment for design and performance analysis of computational models on high performance computers. This report serves as a theoretical manual for selected algorithms implemented within the Dakota software. It is not intended as a comprehensive theoretical treatment, since a number of existing texts cover general optimization theory, statistical analysis, and other introductory topics. Rather, this manual is intended to summarize a set of Dakota-related research publications in the areas of surrogate-based optimization, uncertainty quantification, and optimization under uncertainty that provide the foundation for many of Dakota’s iterative analysis capabilities.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Time-resolved biphase signatures of quadratic nonlinearity observed in coupled Alfvén eigenmodes on the DIII-D tokamak

We report the detection of nonstationary quadratic coupling between toroidicity-induced Alfvén eigenmodes (TAEs) on sub-millisecond time scales. Identification of phase coherency between multiple TAEs and nonlinearly generated modes is facilitated by wavelet-based bicoherence analysis of time-series from inductive coils, taken from a DIII-D discharge heated by neutral beam injection (NBI). Characterization of nonlinear three-wave interaction is inferred by stationary local bispectrum phase (biphase) and confirmed via bandpass filtering. Biphase dynamics associated with prominent bispectral features are well-resolved in time and consistent with transient quadratic coupling. Onset and duration of nonlinearity are correlated with enhanced amplitude of participating TAEs; coincident changes in amplitude are observed for modes at difference frequency |f TAE,1 - f TAE,2 |.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Subharmonic Shapiro steps in depinning dynamics of a two-dimensional solid dusty plasma modulated by one-dimensional nonlinear deformed periodic substrates

Langevin dynamical simulations are performed to investigate the depinning dynamics of a two-dimensional (2D) solid dusty plasma, which is modulated by one-dimensional (1D) nonlinear deformed periodic substrates, and also driven by the combination of the DC and AC forces. As the DC driving force increases gradually, pronounced subharmonic and harmonic Shapiro steps are discovered under the combined modulation of the deformed substrate and the external AC drive. Here, these observed subharmonic and harmonic Shapiro steps are attributed to the dynamic mode locking. The data analysis results indicate that the nonlinear deformed substrate strongly influences these subharmonic Shapiro steps, which can be accurately diagnosed using the fraction of sixfold coordinated particles. Furthermore, the diagnostic of the kinetic temperature clearly indicates the difference between harmonic and subharmonic Shapiro steps, i.e., the particle motion is slightly less synchronized at subharmonic Shapiro steps, caused by the unstable locations of the deformed substrate.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Scaling Behavior of Thermally Driven Fractures in Deep Low‐Permeability Formations: A Plane Strain Model With 1‐D Heat Conduction

Abstract Injection of cold fluids through/into deep formations may cause significant cooling, thermal stress, and possible thermal fracturing. In this study, the thermal fracturing of low‐permeability formations under one‐dimensional heat conduction was investigated using a plane strain model. Dimensionless governing equations, with dimensionless fracture length , aperture , spacing , time , and effective confining stress , were derived. Solution of single thermal fracture was derived analytically, while solution of multiple fractures with constant (or dynamic) spacing were obtained using the displacement discontinuity method (and stability analysis). For single fracture, increases nonlinearly with and then transitions to scaling law , indicating that late‐time fracture length increases linearly with the square root of cooling time. For constantly spaced fractures, deviates from the single‐fracture solution at a later for a larger , showing slower propagation under inter‐fracture stress interaction . For dynamically spaced fractures, fracture arrest induced by stress interaction was determined by the stability analysis; the fully transient solution provides evolution of dimensionless fracture length, spacing, aperture, and pattern; a similar scaling law, with , obtained shows the effect of both stress interaction and fracture arrest. The solution and scaling law provide fast predictions for all reservoir and cooling conditions using (single) model parameter . Application to a geothermal site with demonstrates that thermal fractures reach 0.67, 6.25, and 78.00 m in length, 0.49, 2.30, and 13.00 m in spacing, and 0.43, 2.09, and 12.19 mm in aperture at 1, 100, and 10,000 days.

Chen, Bin↗
DOE OSTI2025

Grain2Mesh: Mesh Generation for Grain-Scale Nonlinear Elasticity Modeling

The nonlinear hysteretic behavior of rocks under cyclic loading is a crucial area of study in geomechanics. The macroscopic response of a variety of materials has been found to be contingent upon the behavior of the micro-scale structure. This project aims to develop a functional and maintainable software package for generating a multi-phase numerical mesh and accompanying simulation files for finite element modeling used in computational mechanics solvers. Meshes generated from images often lack key preprocessing that reduces noise and prevents mesh element distortion that can increase computational cost. By incorporating user feedback throughout, grain2mesh ensures a high-fidelity mesh that can be used to model grain-scale interactions such as shearing, crack propagation, and interfacial material contrast. Scientific applications of this software include material fracturing, stress-strain analysis for natural and engineered materials, and nonlinear meso-scale analysis.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

PETSc TSAdjoint: A Discrete Adjoint ODE Solver for First-Order and Second-Order Sensitivity Analysis

Here, we present a new software system PETSc TSAdjoint for first-order and second order adjoint sensitivity analysis of time-dependent nonlinear differential equations. The derivative calculation in PETSc TSAdjoint is essentially a high-level algorithmic differentiation process. The adjoint models are derived by differentiating the timestepping algorithms and implementing them based on the parallel infrastructure in PETSc. Full differentiation of the library code, including MPI routines, is avoided, and users do not need to derive their own adjoint models for their specific applications. PETSc TSAdjoint can compute the first-order derivative, that is, the gradient of a scalar functional, and the Hessian-vector product, which carries second-order derivative information, while requiring minimal input (a few callbacks) from the users. The adjoint model employs optimal checkpointing schemes in a manner that is transparent to users. Finally, usability, efficiency, and scalability are demonstrated through examples from a variety of applications.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Tree-ring evidence marks year 2022 as the driest spring season in nearly four centuries in the Western Himalayas

The Hindukush-Karakoram-Himalayan region is a crucial freshwater source for billions across South Asia, yet its climate remains poorly understood due to limited long-term records. Winter and spring precipitation govern snow accumulation and downstream water availability in dry months particularly across the Western Himalayas (WH), but recent decades show intensifying droughts with unclear long-term context. Here, we have reconstructed a nearly four-century-long spring i.e. February to May (FMAM) precipitation for the Lahaul region of the (WH), an area dominated by the Western Disturbances. This record was developed using moisture-sensitive Cedrus deodara (Deodar) tree-rings from three high-elevation sites. A regional composite tree-ring-width chronology, developed through a Nested Principal Component Analysis and modeled with a nonlinear Generalized Additive Model (GAM) that explains 71 % of the variance during the calibration period. We identified the last two decades as the most precipitation deficit phase and the year 2022 showing the driest FMAM on record. The observed rise in the FMAM dry episodes post 1999 CE in our reconstruction, corresponds to the meteorological records. This recent drying is linked to a northward shift of the subtropical westerly jet and reduced moisture transport, both associated with unusual sea surface temperature patterns in the tropical Indian Ocean and the Western Pacific Ocean. Our results provide compelling evidence of long-term hydroclimatic instability in the WH and emphasize the value of tree-ring records in extending precipitation histories beyond the instrumental observations. Such reconstructions can be benchmarks to validate high-resolution climate models and formulate adaptation policies to mitigate future risks.

Cedrus deodara↗
DOE OSTI2026

Strain-Modified Raman Responses in Monolayer MoS 2 Nanobubbles Resolved at 5 nm

The formation of nanoscale bubbles is an unavoidable consequence during the transfer of two-dimensional materials onto target substrates, driven by van der Waals interactions at the interface. While often viewed as imperfections, these nanoscale bubbles have garnered considerable scientific interest due to the substantial in-plane strain gradients they induce, which in turn give rise to a variety of intriguing optoelectronic effects, particularly in semiconducting transition metal dichalcogenides. Determining and analyzing the strain distribution within nanobubbles at the nanoscale is crucial for advancing our understanding of these underlying strain-induced effects. Here, we present a high-resolution scanning tunneling microscopy-based tip-enhanced Raman spectroscopic investigation of localized nanoscale strain distribution within the nanobubbles formed between monolayer MoS 2 and Au interface. By employing cryogenic temperature (78 K), we successfully differentiate the nanoscale Raman signatures between nanobubble edge and pristine MoS 2 . We verify a maximum tensile strain of ∽1.15–1.34% at the nanobubble edge, which gradually diminishes toward the center, yielding a cross-sectional strain profile consistent with a doughnut-shaped distribution. Furthermore, we report to achieve ∽5 nm spatial resolution in probing such edge-localized strain within the nanobubble. In addition, comparative average strain analysis of such MoS 2 nanobubbles is conducted via geometric mechanistic analysis such as membrane and nonlinear plate theories, providing key insight into the geometric nature near the bubble edge. Our work provide fundamental information about strain-induced nanoscale chemical understanding of 2D materials on the nanometer scale, paving the way for practical applications of nanobubbles in strain-engineered optoelectronic devices.

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

Scalable edge clustering of dynamic graphs via weighted line graphs

Timestamped relational datasets consisting of records (or connections) between pairs of entities are ubiquitous in network science. For applications like peer-to-peer communication, email, various social network interactions, and computer network security, it is useful to organize these records into groups based on how and when they are occurring. Weighted line graphs offer a natural way to model how records are related in such datasets but for large real-world graph topologies, building and utilizing the line graph is prohibitively expensive. Here, we present the framework to cluster the edges of a dynamic graph via the associated line graph that contains two major contributions. The first is a method to work with the line graph implicitly and the second is a distributed scale implementation of an agglomerative hierarchical graph clustering algorithm. We outline a novel hierarchical dynamic graph edge clustering approach that efficiently breaks massive relational datasets into small sets of edges containing events at various timescales. This is in stark contrast to traditional graph clustering algorithms that prioritize highly connected (clique-like) community structures. Our approach relies on constructing a sufficient subgraph of a weighted line graph and applying a hierarchical agglomerative clustering. This approach is related to scalable techniques from spatial clustering, nonlinear-dimension reduction, topological data analysis, and draws particular inspiration from HDBSCAN. As an edge clustering, this method yields an overlapping node clustering. Our algorithm is parallelizable and we demonstrate efficient clustering of a billion-scale, real-world dynamic graph into small edge sets that correlate in topology and time. The entire clustering process for a graph with tens of billions of edges takes just a few minutes of run time on 256 nodes of a distributed compute environment. We argue how the output of the edge clustering is useful for a multitude of data visualization and powerful machine learning tasks, both involving the original massive dynamic graph data and metadata associated with the nodes and edges. Finally, we describe how this approach can be extended to dynamic hypergraphs and dynamic graphs/hypergraphs with unstructured data living on vertices and edges.

Data Analysis↗
DOE OSTI2025

Revisiting the gluon density from the Balitsky-Kovchegov equation with kinematical constraints and large 𝑥 terms

We perform analysis of the small 𝑥 nonlinear evolution equation formulated in momentum space supplemented by higher order terms. The equation is defined in the wide range of transverse momentum and longitudinal momentum fractions extending previous studies performed in Kutak and Kwiecinski [Screening effects in the ultrahigh-energy neutrino interactions, Eur. Phys. J. C 29, 521 (2003).] and in Kutak and Stasto [Unintegrated gluon distribution from modified BK equation, Eur. Phys. J. C 41, 343 (2005).]. The linear part of the equation is motivated by the renormalization group improved small 𝑥 approach, which accounts for resummation of higher orders, and includes collinear splitting functions and kinematical constraints. The solution to the equation is then used to perform the fit to deep inelastic scattering reduced cross section data.

Perturbative QCD↗
DOE OSTI2024

Method for Coupled Electromagnetic and Circuit Simulations to Evaluate Surge Arrester Performance in Protecting Equipment Against E1 HEMP

Surge arrester behavioral modeling for realistic systems embedded in an E1 high-altitude electromagnetic pulse environment inherently encompasses three interconnected complications: (1) the need to account for signal propagation across two domains, electromagnetics and electrical; (2) the need to include both linear and nonlinear circuit components in the analysis; and (3) the need to understand that the over-current and over-voltage mitigation performance is dependent not only on the properties of the surge arrester and protected load but also on the topology of the overall electrical network. This study presents a framework to address these challenges in a systematic manner to consider the effectiveness of protective measures for a common class of equipment in power generation facilities. Full-wave simulations were carried out to derive circuit-domain (i.e., lumped element–based) equivalent models for the excitation waveform and the physical components of the system. Then, these equivalent models were imported to a circuit solver and combined with a high-frequency surge arrester model to evaluate mitigation performance. The methodology outlined is general enough such that it can be applied for other electromagnetic interference problems that involve E2/E3 HEMP or microwave emissions.

42 ENGINEERING↗
DOE OSTI2022

On the Convergence of Overlapping Schwarz Decomposition for Nonlinear Optimal Control

Here, we study the convergence properties of an overlapping Schwarz decomposition algorithm for solving nonlinear optimal control problems (OCPs). The algorithm decomposes the time domain into a set of overlapping subdomains, and solves all subproblems defined over subdomains in parallel. The convergence is attained by updating primal-dual information at the boundaries of overlapping subdomains. We show that the algorithm exhibits local linear convergence, and that the convergence rate improves exponentially with the overlap size. We also establish global convergence results for a general quadratic programming, which enables the application of the Schwarz scheme inside second-order optimization algorithms (e.g., sequential quadratic programming). The theoretical foundation of our convergence analysis is a sensitivity result of nonlinear OCPs, which we call "exponential decay of sensitivity" (EDS). Intuitively, EDS states that the impact of perturbations at domain boundaries (i.e., initial and terminal time) on the solution decays exponentially as one moves into the domain. Here, we expand a previous analysis available in the literature by showing that EDS holds for both primal and dual solutions of nonlinear OCPs, under uniform second-order sufficient condition, controllability condition, and boundedness condition. We conduct experiments with a quadrotor motion planning problem and a partial differential equations (PDE) control problem to validate our theory, and show that the approach is significantly more efficient than alternating direction method of multipliers and as efficient as the centralized interior-point solver.

42 ENGINEERING↗
DOE OSTI2021

The Relationship between Precipitation and Precipitable Water in CMIP6 Simulations and Implications for Tropical Climatology and Change

It is well documented that over the tropical oceans, column-integrated precipitable water (pw) and precipitation (P) have a nonlinear relationship. In this study moisture budget analysis is used to examine this P–pw relationship in a normalized precipitable water framework. It is shown that the parameters of the nonlinear relationship depend on the vertical structure of moisture convergence. Specifically, the precipitable water values at which precipitation is balanced independently by evaporation versus by moisture convergence define a critical normalized precipitable water, pw nc . This is a measure of convective inhibition that separates tropical precipitation into two regimes: a local evaporation-controlled regime with widespread drizzle and a precipitable water–controlled regime. Most of the 17 CMIP6 historical simulations examined here have higher pw nc compared to ERA5, and more frequently they operate in the drizzle regime. When compared to observations, they overestimate precipitation over the high-evaporation oceanic regions off the equator, thereby producing a ‘‘double ITCZ’’ feature, while underestimating precipitation over the large tropical landmasses and over the climatologically moist oceanic regions near the equator. The responses to warming under the SSP585 scenario are also examined using the normalized precipitable water framework. It is shown that the critical normalized precipitable water value at which evaporation versus moisture convergence balance precipitation decreases as a result of the competing dynamic and thermodynamic responses to warming, resulting in an increase in drizzle and total precipitation. Statistically significant historical trends corresponding to the thermodynamic and dynamic changes are detected in ERA5 and in lowintensity drizzle precipitation in the PERSIANN precipitation dataset.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Demystifying Power System Oscillations - Recent and Ongoing Efforts

Summary of recent and ongoing research efforts in power system oscillation analysis and techniques for locating the source of oscillations. Review of existing technologies, best practice in industry, as well as challenges.

41 EE - Solar Energy Technologies Office (EE-4S)↗