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 559 records · Page 31
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Enabling Hyper-Differential Sensitivity Analysis for Ill-Posed Inverse Problems

Inverse problems constrained by partial differential equations (PDEs) play a critical role in model development and calibration. In many applications, there are multiple uncertain parameters in a model that must be estimated. However, high dimensionality of the parameters and computational complexity of the PDE solves make such problems challenging. A common approach is to reduce the dimension by fixing some parameters (which we will call auxiliary parameters) to a best estimate and use techniques from PDE-constrained optimization to estimate the other parameters. In this article, hyper-differential sensitivity analysis (HDSA) is used to assess the sensitivity of the solution of the PDE-constrained optimization problem to changes in the auxiliary parameters. Foundational assumptions for HDSA require satisfaction of the optimality conditions which are not always practically feasible as a result of ill-posedness in the inverse problem. Here we introduce novel theoretical and computational approaches to justify and enable HDSA for ill-posed inverse problems by projecting the sensitivities on likelihood informed subspaces and defining a posteriori updates. Our proposed framework is demonstrated on a nonlinear multiphysics inverse problem motivated by estimation of spatially heterogeneous material properties in the presence of spatially distributed parametric modeling uncertainties.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Cs X Si 15 P 21 ( X = Sn or Pb): Polar Noncentrosymmetric Si–P Frameworks Stabilized by Covalent X –P Bonding

Metal silicon phosphides composed of earth-abundant Si and P tend to exhibit semiconducting properties and adopt diverse crystal structures with relatively small additions of structure-directing elements. The potential of silicon phosphide materials in nonlinear optical applications has been hindered by the inability to systematically produce noncentrosymmetric structures with such a flexible framework. Here, in this work, two isostructural compounds with a novel noncentrosymmetric structure were made possible by the inclusion of elements with stereochemically active lone pairs (Sn 2+ and Pb 2+ ). The structures were determined through single-crystal and synchrotron powder X-ray diffraction. Analysis of chemical bonding in real space through the electron localization function revealed stereochemically active Pb 2+ and Sn 2+ species in a trigonal pyramidal coordination with {Pb/Sn}–P bonds. Such covalent bonding between Pb and P is quite uncommon in extended solids and has been reported in a few rare instances. Band structure calculations and linear optical measurements confirm the semiconducting nature of Cs X Si 15 P 21 ( X = Sn or Pb). The synthesis was optimized to yield high-purity polycrystalline samples. The nonlinear optical properties show promising second-harmonic generation (SHG) coefficients from the Kurtz–Perry method. First-principles calculations of the nonlinear optical properties support the experimentally determined SHG values and provide moderate values of birefringence, suggesting Cs X Si 15 P 21 could be phase-matchable and practical nonlinear optical materials in the mid-IR region.

crystal structure↗
DOE OSTI2021

IMEX Runge-Kutta Parareal for Non-diffusive Equations

Parareal is a widely studied parallel-in-time method that can achieve meaningful speedup on certain problems. However, it is well known that the method typically performs poorly on non-diffusive equations. This paper analyzes linear stability and convergence for IMEX Runge-Kutta Parareal methods on non-diffusive equations. By combining standard linear stability analysis with a simple convergence analysis, we find that certain Parareal configurations can achieve parallel speedup on non-diffusive equations. These stable configurations all posses low iteration counts, large block sizes, and a large number of processors. Numerical examples using the nonlinear Schrodinger equation demonstrate the analytical conclusions.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Reduced models for ETG transport in the tokamak pedestal

This paper reports on the development of reduced models for electron temperature gradient (ETG) driven transport in the pedestal. Model development is enabled by a set of 61 nonlinear gyrokinetic simulations with input parameters taken from pedestals in a broad range of experimental scenarios. The simulation data have been consolidated in a new database for gyrokinetic simulation data, the multiscale gyrokinetic database (MGKDB), facilitating the analysis. The modeling approach may be considered a generalization of the standard quasilinear mixing length procedure. The parameter η, the ratio of the density to temperature gradient scale length, emerges as the key parameter for formulating an effective saturation rule. With a single order-unity fitting coefficient, the model achieves an error of 15%. A similar model for ETG particle flux is also described. We also present simple algebraic expressions for the transport informed by an algorithm for symbolic regression.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Insights on the Influence of the Central-Cut Width of the Box Assembly with Removable Component on Its Dynamical Responses

This investigation focuses on the dynamical effects caused by varying the central-cut width within the Box Assembly with Removable Component (BARC) system. The central-cut widths included in this study are a 0.5″ cut, a 0.25″ cut, a thin 0.1″ cut, and a structure that did not have a cut at all. Finite element analysis was conducted to determine the mode shapes and natural frequencies of each of the BARC structures. Structural dynamics experiments were run to examine the effects of the central-cut width on the dynamical responses and nonlinear characteristics of the BARC system. Free vibration testing with an impact hammer was carried out to excite the system and extract the dominant frequencies and directions of the significant responses. A pseudorandom vibration test that allows for the qualitative determination of any nonlinear behavior within the system was performed. This type of behavior can include nonlinear softening, nonlinear hardening, and the most common, nonlinear damping due to the presence of several bolted-joint connections and the possible activation of geometric and inertia nonlinearities. To quantitatively investigate the impacts of the central-cut width on the dynamics of the system, swept sinusoidal testing was conducted. It is determined that almost all systems with central cuts demonstrate the presence of nonlinear softening, but at times, nonlinear hardening trends are seen, particularly in the 0.1″ cut and no-cut systems when testing harmonically. Each of the central-cut systems displays nonlinear damping, with the amount of damping generally increasing as the central cut decreases in size. The effect of the central cut of the BARC system on the mode-switching ability of the system is negligible; however, mode switching takes place when comparing the central-cut configurations to the no-cut one. These results show the significance of accurately measuring the central-cut width and how geometric uncertainty may change the dynamical responses and nonlinear properties of the system.

Padilla, Christopher (ORCID:000900033446732X)↗
DOE OSTI2026

Exploring the Origins of Anti-Ambipolarity in BBL Polymer: Links to Redox Chemistry, Electronic Structure, and Structural Dynamics

We examine the intrinsic physical-chemical properties of the conjugated ladder-type polymer poly(benzimidazobenzophenanthroline) (BBL) in response to electron transfer. We aim at explaining the origin of the anti-ambipolar behavior behind the observed BBL nonlinear response associated with specific device architectures. To elucidate this point, we use theory and computation based on first principles, including density functional theory optimizations, ab initio molecular dynamics, time-dependent DFT, and Marcus-theory analysis. Our results reveal that this redox response is not simply monotonic but follows an alternating odd/even pattern in which gap narrowing and reopening occur sequentially before near-gapless behavior emerges at high charging. Converging theoretical evidence in this work demonstrates that bell shaped conductivity in BBL originates in its fundamental electronic structure and supramolecular organization.

FOS: Physical sciences↗
DOE OSTI2020

Robust PCA-Deep Belief Network Surrogate Model for Distribution System Topology Identification with DERs

With the expansion of distribution networks and increased penetration of distributed energy resources (DERs), it is becoming increasingly important to obtain accurate distribution network topology in real-time. In this paper, a robust principal component analysis coupled deep belief network (PCA-DBN) surrogate model is proposed for distribution system topology identification. It integrates the benefits of robust feature extraction from PCA to deal with data quality issues and filter out noise, and the strength of DBN in capturing the nonlinear relationship between voltage amplitudes and the binary states of switchable connections. This also significantly reduces the DBN training complexity without loss of accuracy. It is shown that the widely used standard deviation of voltage drop and the voltage covariance matrix features yield less accuracy as compared to that of the voltage amplitudes in presence of high penetration of DERs and ZIP loads. Comparison results with other alternatives, such as the random forest (RF), multi-output regression (MOR) and the traditional DBN methods demonstrate that the proposed method can achieve a much higher topology identification accuracy while maintaining robustness to missing data and measurement noise under various penetration levels of DERs.

deep belief network↗
DOE OSTI2023

Design and Analysis of a Floating-Wind Shallow-Water Mooring System Featuring Polymer Springs

This paper presents a mooring system design featuring polymer springs for the VolturnUS-S 15-MW reference floating wind turbine in site conditions for the New York Bight at a 50-m water depth. Polymer springs have a nonlinear stress-strain curve that allows a stiffer response at low loads and a more flexible response at higher loads, potentially reducing peak mooring line tensions. The mooring dynamics model MoorDyn has been extended to model springs with nonlinear tension-strain curves from a user-inputted look-up table. This MoorDyn modeling advancement is verified against OrcaFlex simulation results. Using MoorDyn's updated capabilities, a spring-equipped catenary mooring system is designed for the 15-MW floating system, along with a baseline catenary mooring system that does not use polymer springs. The floating wind turbine simulator OpenFAST is used to simulate the mooring systems in design-driving load cases to show the effect of polymer springs on key dynamic behaviours. The results show that the spring-equipped design reduces peak tensions by up to 60%, whereas the turbine offsets stay within a maximum of 7.2 m, which is still a reasonable offset limit for cable considerations. The reduction in peak tensions results in a significant decrease in damage equivalent loads - on the order of 50% for upwind lines in fully loaded conditions. These results show that polymer springs can effectively reduce peak tensions and fatigue loads in mooring systems at shallow water depths.

floating offshore wind↗
DOE OSTI2024

Communication Lower Bounds and Optimal Algorithms for Multiple Tensor-Times-Matrix Computation

Multiple tensor-times-matrix (Multi-TTM) is a key computation in algorithms for computing and operating with the Tucker tensor decomposition, which is frequently used in multidimensional data analysis. Here, we establish communication lower bounds that determine how much data movement is required (under mild conditions) to perform the Multi-TTM computation in parallel. The crux of the proof relies on analytically solving a constrained, nonlinear optimization problem. We also present a parallel algorithm to perform this computation that organizes the processors into a logical grid with twice as many modes as the input tensor. We show that, with correct choices of grid dimensions, the communication cost of the algorithm attains the lower bounds and is therefore communication optimal. Finally, we show that our algorithm can significantly reduce communication compared to the straightforward approach of expressing the computation as a sequence of tensor-times-matrix operations when the input and output tensors vary greatly in size.

HBL-inequalities↗
DOE OSTI2022

Artificial Intelligence/Machine Learning Technologies for Advanced Reactors (Workshop Summary Report)

A workshop on artificial intelligence and machine learning (AI/ML) for advanced reactors (AR) was held October 5-6, 2021. The workshop was to be attended in-person at ANL but COVID restrictions forced the workshop to go virtual. The objectives of the workshop were to identify the most promising AI/ML opportunities for improving advanced reactor design, optimizing plant performance, and enhancing economic competitiveness and to develop an understanding of the scientific, engineering and licensing challenges facing their application. The workshop planning committee included GAIN, EPRI and NEI and members of three national laboratories (ANL, INL, and ORNL). The workshop was attended by more than 200 individuals representing academic and scientific institutions and the nuclear power industry. The definition put forth for an AI/ML system was one that perceives its environment and takes actions that maximize its chance of achieving its goals. In this report AI/ML refers to next generation algorithms that include deep learning, statistical analysis and data analytics and associated scientific computing and their potential application to the design, licensing, operation and maintenance of ARs. These methods typically incorporate models built from process data and may also include data generated by simulations that represent the behavior of a system. The workshop was organized in response to the growing interest in application of AI/ML for improving the economic competitiveness of nuclear energy. Increasingly more resources are being allocated to investigating the benefits of AI/ML methods. The DOE created the Artificial Intelligence & Technology Office to promote their development. And within the Office of Nuclear Energy, resources have been allocated to explore and understand the potential benefits of AI/ML. Additionally, the national laboratories are strategically positioned with DOE computing facilities such as Summit, Perlmutter, Aurora and Frontier that support large-scale simulations, hybrid HPC models with AI surrogates, and the exploration of new types of generative models emerging from multi-model data streams and sources. The workshop was organized with members of the AR community to understand the effort and to identify the level of interest and progress in this emerging technology. The workshop discussions focused on identifying opportunities for AI/ML across diverse areas of the nuclear industry and identifying current scientific and engineering challenges for advanced reactors that might be addressed through transformational uses of AI/ML. Discussion panels focused on four high-interest technical domains for advanced reactors: design, maintenance and operations, energy storage, and materials. The results of those discussions are summarized in this report. This includes opportunities that were identified for exploiting AI techniques and methods to improve the efficacy and efficiency of reactor analysis and to improve the operation and optimization of advanced reactors. Advanced reactor developers expressed an interest in learning more about AI/ML methods and their application. This included understanding whether ML methods can provide an advantage over existing nonlinear data regression methods for collapsing high-fidelity simulation results into faster running models. A consensus emerged that AR advances planned for the next decade will benefit from the use of AI/ML tools. The need exists to understand and model complex systems across length scales and modalities. AI/ML is a tool for discovery that can yield a set of engineering principles for use by nuclear engineers, licensing bodies, and operators to solve problems in plant design, safety analyses, autonomous operation, and predictive maintenance. While AI/ML represents a new set of tools, an awareness by the nuclear community of the full potential is still in the early stages so there is a need to increase awareness. It appears that the wide-spread adoption of AI/ML tools for ARs would be facilitated by future educational workshops that describe foundational methods and capabilities and describe successful applications.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2026

Quantitative analysis of zonal flow influence on turbulent plasmas driven by trapped electrons

The role of zonal flow (ZF) in the turbulence saturation of trapped electron modes (TEM) in magnetically confined plasmas is revisited. Here, in this study, we examine ZF excitation and saturation mechanisms in TEM turbulence using detailed free-energy transfer diagnostics from nonlinear gyrokinetic simulations. When the ion channel becomes subdominant, electrons shift to transferring energy to the zonal component, which results in zonal flows always playing a significant though subdominant role, in the case of the temperature gradient driven TEM. The velocity-space structure of energy transfer reveals that trapped electrons directly couple to ZF. Regarding the saturation physics, the zonal flow advection, drift-wave–drift-wave interactions, and stable modes at the pump waves are compared. When ZF saturation is weak, drift-wave–drift-wave interactions play a larger role in maintaining turbulence saturation. The findings elucidate the roles of nonzonal to zonal flow coupling, stable modes, and nonzonal wave-wave interactions, and provide a basis for improving reduced models of turbulent transport.

Song, Jiheon [Hanyang University, Seoul (Korea, Re↗
DOE OSTI2022

Design and Analysis of a Floating-Wind Shallow-Water Mooring System Featuring Polymer Springs: Preprint

In this paper, a mooring system featuring polymer springs is designed for the VolturnUS-S 15 MW reference floating wind turbine in site conditions for the New York Bight at a 50-m water depth. Polymer springs have a nonlinear stress-strain curve that allows a stiffer response at low loads and a more flexible response at higher loads, potentially reducing peak mooring line tensions. The mooring dynamics model MoorDyn has been extended to model springs with nonlinear tension-strain curves from a user-inputted look-up table. This MoorDyn modelling advancement is verified against OrcaFlex simulation results. Using MoorDyn's updated capabilities, a spring-equipped catenary mooring system is designed for the 15 MW floating system, along with a baseline catenary mooring system that does not use polymer springs. The floating wind turbine simulator OpenFAST is used to simulate the mooring systems in design-driving load cases to show the effect of polymer springs on key dynamic behaviours. The results show that the spring-equipped design reduces peak tensions by up to 60%, while the turbine offsets stay within a maximum of 7.2 m, which is still a reasonable offset limit for cable considerations. The reduction in peak tensions results in a significant decrease in damage equivalent loads, on the order of 50% for upwind lines in fully loaded conditions. These results show that polymer springs can effectively reduce peak tensions and fatigue loads in mooring systems at shallow water depths.

floating offshore wind↗
DOE OSTI2020

Eigenfrequency constrained topology optimization of finite strain hyperelastic structures

This paper incorporates hyperelastic materials, nonlinear kinematics, and preloads in eigenfrequency constrained density–based topology optimization. The formulation allows for initial finite deformations and subsequent small harmonic oscillations. The optimization problem is solved by the method of moving asymptotes, and the gradients are calculated using the adjoint method. Both simple and degenerate eigenfrequencies are considered in the sensitivity analysis. A well-posed topology optimization problem is formulated by filtering the volume fraction field. Numerical issues associated with excessive distortion and spurious eigenmodes in void regions are reduced by removing low volume fraction elements. The optimization objective is to maximize stiffness subject to a lower bound on the fundamental eigenfrequency. Numerical examples show that the eigenfrequencies drastically change with the load magnitude, and that the optimization is able to produce designs with the desired fundamental eigenfrequency.

42 ENGINEERING↗
DOE OSTI2021

Quantifying Polaron Mole Fractions and Interpreting Spectral Changes in Molecularly Doped Conjugated Polymers

Molecular doping of conjugated polymers causes bleaching of the neutral absorbance and results in new polaron absorbance transitions in the mid and near infrared. In this report the concentration dependent changes in the spectra for a series of molecularly doped diketopyrrolopyrrole (DPP) co-polymers with a series of ultra-high electron affinity cyanotrimethylenecyclopropane-based dopants is analyzed. With these strong dopants the polaron mole fraction (Θ) reaches saturation. Analysis of the full spectrum enables separation of neutral and polaron signals and quantification of the polaron mole fraction using a simple noninteracting site model. The peak ratios for both neutral and polaron peaks change systematically with increasing polaron mole fraction for all measured polymers. Analysis of the spectral changes indicates that the polaron mole fraction can be quantified to within 5%. While the total change in the absorbance spectrum with increasing polaron mole fraction is linear, the lowest energy polaron peak (P1) grows nonlinearly, which indicates increased polarization/delocalization. Molecular doping of polymers that form either H- or J-aggregates shows systematically different spectral changes in the vibronic peak ratios of the neutral spectra and provides insights into the polymer configuration at undoped sites in the film.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Integration and Optimization of a Waste Heat Driven Organic Rankine Cycle for Power Generation in Wastewater Treatment Plants

The study focuses on achieving energy self-sufficiency in Wastewater Treatment Plants by proposing a comprehensive model for integrating, sizing, and optimizing an Organic Rankine Cycle system. The Organic Rankine Cycle system is designed to utilize waste heat from the gensets at As Samra Wastewater Treatment Plant in Jordan, where it will contribute to the overall electrical energy supply of the plant. Real data from As Samra Wastewater Treatment Plant is used to model and calculate the available waste heat using TRNSYS® software. The Organic Rankine Cycle model is then developed using ASPEN PLUS® software to explore the impact of operational parameters and determine their optimal values for maximizing the plant's energy profile. An economic analysis is conducted to assess the feasibility of the proposed model, considering system components, installation, operation, and maintenance costs. To optimize the Organic Rankine Cycle system, the study employs the Multi-Output Support Vector Regression technique to capture nonlinear relationships between independent variables (fluid type, turbine inlet pressure, turbine inlet temperature, turbine outlet pressure, and mass flow rate) and dependent variables (pump power input, waste heat input, and turbine specific work). The Osprey optimization algorithm is used to address the multi-objective optimization problem, with the proposed Pareto-based Osprey Optimization Algorithm and the Multi-Objective Particle Swarm Optimization technique being employed to evaluate critical performance and economic parameters such as system thermal efficiency, net power output, and the levelized cost of electricity. The results of the optimization strategies indicate that the M-SVR model's prediction accuracy is significantly improved after parameter optimization, with the model returning high R 2 and low Mean Square Error values of 0.991 and 0.00216, respectively. The Pareto-Based Osprey Optimization Algorithm optimizer identifies the best working fluid as Isobutane/Isopentane in a ratio of 66:34, with optimal turbine inlet pressure and temperature of 15 bars and 218 °C, respectively. In conclusion, the Organic Rankine Cycle model at these optimal conditions achieves a cycle efficiency of 19.93% and an Levelized Cost of Electricity value of 0.0353 USD/kWh.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2021

Evidence for a thermally driven charge-density-wave transition in 1T-TaS 2 thin-film devices: Prospects for GHz switching speed

We report on the room-temperature switching of 1T-TaS 2 thin-film charge-density-wave devices, using nanosecond-duration electrical pulsing to construct their time-resolved current–voltage characteristics. The switching action is based upon the nearly commensurate to incommensurate charge-density-wave phase transition in this material, which has a characteristic temperature of 350 K at thermal equilibrium. For sufficiently short pulses, with rise times in the nanosecond range, self-heating of the devices is suppressed, and their current–voltage characteristics are weakly nonlinear and free of hysteresis. This changes as the pulse duration is increased to ~200 ns, where the current develops pronounced hysteresis that evolves nonmonotonically with the pulse duration. By combining the results of our experiments with a numerical analysis of transient heat diffusion in these devices, we clearly reveal the thermal origins of their switching. In spite of this thermal character, our modeling suggests that suitable reduction of the size of these devices should allow their operation at GHz frequencies.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Response of salt intrusion in a tidal estuary to regional climatic forcing

Abstract Salinity distribution in a large tidal estuary is subject to estuarine adjustment under the influences of multiple physical drivers such as freshwater pulses and sea level rise, and is crucial to upstream water quality, aquaculture, and ecosystem functions of the estuary. To better understand the estuarine salinity response to climate change, the unstructured-grid Finite Volume Community Ocean Model was implemented to simulate the salt intrusion in the Delaware Bay Estuary. The model was first validated by multiple observational data sets and subsequently applied in an idealized setting to examine the response of salt front to freshwater pulses in high flow conditions, followed by a long-term drought condition supported by a multi-decadal streamflow drought analysis in the estuary. The model results showed that after the freshwater pulses the salt front location moved further upstream with sea level rise. Under the simulated long-term drought condition, the adjustment timescale of salt intrusion varies nonlinearly with sea level rise. With a significant increase in sea level rise, the adjustment timescale starts to decrease. This shift suggests a transition into a different regime where the estuary becomes more stratified, as indicated by an increasing bulk Simpson number with rising sea levels.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Response Of Salt Intrusion In A Tidal Estuary To Regional Climatic Forcing (FVCOM Output for the baseline scenario)

Salinity distribution in a large tidal estuary is subject to estuarine adjustment under the influences of multiple physical drivers such as freshwater pulses and sea level rise, and is crucial to upstream water quality, aquaculture, and ecosystem functions of the estuary. To better understand the estuarine salinity response to climate change, the unstructured-grid Finite Volume Community Ocean Model was implemented to simulate the salt intrusion in the Delaware Bay Estuary. The model was first validated by multiple observational data sets and subsequently applied in an idealized setting to examine the response of salt front to freshwater pulses in high flow conditions, followed by a long-term drought condition supported by a multi-decadal streamflow drought analysis in the estuary. The model results showed that after the freshwater pulses the salt front location moved further upstream with sea level rise. Under the simulated long-term drought condition, the adjustment timescale of salt intrusion varies nonlinearly with sea level rise. With a significant increase in sea level rise, the adjustment timescale starts to decrease. This shift suggests a transition into a different regime where the estuary becomes more stratified, as indicated by an increasing bulk Simpson number with rising sea levels. These findings have been summarized into a journal manuscript. The data included in this information release request are FVCOM model output associated with the numerical simulations described in the journal manuscript.

Liu, Jinliang↗