Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Nonlinear processes”

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 361 records · Page 20

Inclusive curvaturelike framework for describing dissipation: Metriplectic 4-bracket dynamics

An inclusive framework for joined Hamiltonian and dissipative dynamical systems that are thermodynamically consistent, i.e., preserve energy and produce entropy, is given. The dissipative dynamics of the framework is based on the metriplectic 4-bracket , a quantity like the Poisson bracket defined on phase space functions, but unlike the Poisson bracket has four slots with symmetries and properties motivated by Riemannian curvature. Metriplectic 4-bracket dynamics is generated using two generators, the Hamiltonian and the entropy, with the entropy being a Casimir of the Hamiltonian part of the system. The formalism includes known previous binary bracket theories for dissipation or relaxation as special cases. Rich geometrical significance of the formalism and methods for constructing metriplectic 4-brackets are explored. Many examples of both finite and infinite dimensions are given.

Physics↗

HAIMOS Ensemble Forecasts for Intra-day and Day- Ahead GHI, DNI and Ramps

The objective of this research is to develop a hybrid physics-based/data-driven forecast model to improve direct normal and global horizontal irradiance (DNI and GHI) prediction for horizons ranging from 1 to 72 hours. Project objectives also address key gaps in state-of-the-art solar forecasting: accurate probabilistic solar forecasts and the forecasting of large irradiance ramps (ramp onset and magnitude). The proposed model ensembles Numerical Weather Prediction (NWP) forecasts, determinist physics-based algorithms, and new-generation cloud cover products (high-resolution rapid refresh satellite images and Large Eddy Simulations). The result is the Hybrid Adaptive Input Model Objective Selection (HAIMOS) ensemble model. HAIMOS blends state of the art machine learning methodologies with physics-based models for cloud cover and cloud optical depth forecasts. The technical activities followed a two-pronged strategy. First, the preprocessing of data, the selection of inputs to the nonlinear approximators, the type of approximator and objective functions, and post-processing ensembling techniques included in HAIMOS were all optimized adaptively to find the best model for a specific goal (reduce DNI/GHI forecast error, improve the prediction of ramp onset, etc.). Second, a large effort was put in improving cloud identification and the forecast of cloud cover and cloud optical depth. To this end, new-generation cloud parametrization products were developed in this work. These include improved algorithms to assist in cloud identification, cloud classification and cloud parametrization from satellite images – three key factors in the accuracy of 1 to 6-hours irradiance forecasts and prediction of ramp onset. Furthermore, we also included cloud information extracted from high resolution rapid refresh satellite images (GOES-16) and Large Eddy Simulations (LES). LES was used to model the atmosphere in detail over locations of interest and produce cloud optical depth forecasts. Once these data streams were validated, they were used as input data to the HAIMOS forecast. The model was developed using data from several climatologically distinct locations with potential for high solar penetration. In the last year of the project, we conducted a validation campaign according to the guidelines stipulated by the Topic Area 1 project as described in the FOA. This effort brings, for the first time, proven machine-learning methodologies for generating state-of-the-art solar forecasts interweaved with detailed physics-based models for cloud detection, and cloud optical depth forecasts. HAIMOS will generate accurate irradiance probabilistic forecast to assist in reducing solar generation prediction error. Globally optimized solar forecast models are more likely to impact solar energy stakeholders. The goal of this project was to increase the state-of-the-art forecast skill from their present values of 10 to 35%. At the end of the project, we achieved between 30% and 50% forecast skill across a wide range of horizons for both GHI and DNI.

14 SOLAR ENERGY↗

Breaking Beams: A Kinetic Treatment of the Filamentation and Forward Scatter of Laser Beams in Plasma (Abbreviated Final Report)

Precise and predictable propagation and deposition of laser energy in high energy density (HED) plasma is threatened by the breakup of the beam into filaments. Such filaments can be deflected and sprayed but also locally intensify the light, increasing the likelihood of other problematic interactions with the plasma. Adequately modeling, anticipating, and avoiding such processes in HED experiments requires a thorough understanding of the thresholds and nonlinear saturation mechanisms, particularly for the weakly collisional, high temperature plasmas achieved at the National Ignition Facility (NIF).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Breaking Beams: A Kinetic Treatment of the Filamentation and Forward Scatter of Laser Beams in Plasma

Precise and predictable propagation and deposition of laser energy in high energy density (HED) plasma is threatened by the breakup of the beam into filaments. Such filaments can be deflected and sprayed but also locally intensify the light, increasing the likelihood of other problematic interactions with the plasma. Adequately modeling, anticipating, and avoiding such processes in HED experiments requires a thorough understanding of the thresholds and nonlinear saturation mechanisms, particularly for the weakly collisional, high temperature plasmas achieved at the National Ignition Facility (NIF).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Neural Network Enhanced RKPM for Electrochemical-Mechanical Coupled Damage Modeling of Energy Storage Materials

Energy storage materials undergo significant charge cycling, which makes understanding their reliability and durability fundamental in predicting performance and service life. Strong electrochemical-mechanical coupling and highly anisotropic material properties contribute to the formation and propagation of micro-cracking, largely along material interfaces and grain boundaries. For Li-ion batteries, for example, lithium moving between electrodes during charging and discharging process causes expansion and contraction of grains, and the strongly anisotropic and nonlinearly [Li]-dependent grain material properties can cause grains to expand into and contract away from each other, leading to chemo-mechanical cracking. In the first part of this work, a RKPM based computational framework for solving the coupled solid-phase lithium conservation with Fickian diffusion and the lithium concentration dependent anisotropic mechanical problem subjected to a highly nonlinear Butler-Volmer boundary condition is introduced. The choice of RKPM completeness conditions for lithium concentration and mechanical deformation fields, and the variational consistency condition for the domain integration of the coupled problem is first determined. In the second part of this work, a neural network-enhanced reproducing kernel particle method (NN-RKPM) [1] is leveraged to accurately capture damage and crack propagation throughout the material, by learning the location, orientation, and sharpness of discontinuity while allowing for a coarser nodal distribution than that is necessary for capturing sharp solution transitions using traditional mesh-based methods. NN-RKPM is used to inform how crack opening and closure in turn affect the coupled chemical equations and material microstructure.

damage modeling↗

Block-structured, equal-workload, multi-grid-nesting interface for the Boussinesq wave model FUNWAVE-TVD (Total Variation Diminishing)

Abstract. We describe the development of a block-structured, equal-CPU-load (central processing unit), multi-grid-nesting interface for the Boussinesq wave model FUNWAVE-TVD (Fully Nonlinear Boussinesq Wave Model with Total Variation Diminishing Solver). The new model framework does not interfere with the core solver, and thus the core program, FUNWAVE-TVD, is still a standalone model used for a single grid. The nesting interface manages the time sequencing and two-way nesting processes between the parent grid and child grid with grid refinement in a hierarchical manner. Workload balance in the MPI-based (message passing interface) parallelization is handled by an equal-load scheme. A strategy of shared array allocation is applied for data management that allows for a large number of nested grids without creating additional memory allocations. Four model tests are conducted to verify the nesting algorithm with assessments of model accuracy and the robustness in the application in modeling transoceanic tsunamis and coastal effects.

Choi, Young-Kwang↗

Gamma-radiation-induced negative nonlinear absorption in quartz glass

The effect of ionizing radiation on the nonlinear optical properties of materials needs to be considered for applications that employ high-power laser pulses in extreme radiation environments. The Heraeus Infrasil302 optical quartz glass was irradiated by gamma rays, and its nonlinear absorption was measured using the open-aperture Z-scan technique with nanosecond laser pulses. Irradiation results in a negative nonlinear absorption coefficient whose magnitude increases with irradiation dose. Nonlinear absorption coefficients of -0.78, -1.14, -1.29, and -1.88~ × ~10 -13 m/W are measured for gamma doses of 600 krad, 1.2 Mrad, 3.6 Mrad, and 10 Mrad, respectively. Thermal annealing restores the optical properties of quartz-glass to their unirradiated values. It is suggested that metallic impurities introduced by the manufacturing process have a role in inducing saturated absorption in irradiated quartz glass, which manifests as negative nonlinear absorption coefficient, and is demonstrated by comparing Infrasil302 quartz results to equivalent irradiated measurements for Spectrosil2000 quartz.

07 ISOTOPE AND RADIATION SOURCES↗

Online Monitoring System for Concrete Structures Affected by Alkali-Silica Reaction

This report presents a comprehensive study of developing ultrasonic wave and acoustic emission (AE) techniques for long-term monitoring of alkali-silica reaction (ASR) development in concrete. Small, medium, and full-scale concrete specimens were cast, conditioned, and monitored for various periods, from at least one year to 2.4 years. The reactive concrete specimens contain coarse or fine reactive aggregates to study the effects of different types of reactive aggregates. Confinements were also designed to simulate the 2-dimensional confinement effects of reinforcement in the shield building of nuclear power plants. The concrete specimens were stored in a environmental chamber with high humidity and high temperature to accelerate the ASR development. The ultrasonic monitoring data shows high sensitivity to ASR development and could detect cracking initiation well before visible surface cracks occurred. However, the linear ultrasonic analysis based on wave velocity is strongly affected by temperature variation. Therefore, a nonlinear ultrasonic method was proposed to measure thermally induced nonlinear acoustic responses of concrete (thermal modulation of ultrasonic wave). The measured nonlinear acoustic parameter shows a high correlation with ASR expansion across specimens with different reactive aggregates and confinement conditions. The same conclusion was obtained from nonlinear resonance tests on small concrete prisms. Compared to the linear acoustic methods, the nonlinear tests show high sensitivities to ASR damage from internal microcracking initiation at the early stage to visible cracks at the late stage of ASR. The attributes of the thermal modulation of nonlinear ultrasonic method include high sensitivity, immunity to temperature effects, and strong correlation with ASR expansion, which demonstrate great potentials of the nonlinear ultrasonic method for diagnosis of ASR damage and prediction of concrete deterioration process. Acoustic emission is a passive sensing technique for damage assessment, and access to only one surface is needed even for thick and heavily reinforced elements such as the walls utilized for nuclear shield building. Additionally, relatively few sensors are required to monitor the progression of the damage process. Results indicate that the AE data can be related to the damage rating index, which is a petrography-based means of assessing damage due to ASR in reinforced concrete. Furthermore, AE is capable of detecting ASR damage long before surface cracking is visually noticeable. Boundary conditions play in important role in the progression of ASR damage and differences in boundary conditions are reflected in the AE data. Entropy based data assessment methods provide a means to assess the damage state, and convolutional neural network based data assessment provided a means to assess the damage state in real-time. Results indicate that artificial neural network models may be used as a means to predict volumetric expansion based on AE data.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

A mechanistic interpretation of Nelson curves for PVP failures under high temperature hydrogen attack

As an empirically established design criterion, Nelson curves that relate the service temperature and the allowable hydrogen partial pressure have been developed and utilized for more than sixty years in pressure vessels and piping (PVP) safety design. Despite a relatively clear thermodynamic understanding of the high-temperature-hydrogen-attack (HTHA) problem, the detailed fracture process on the microstructural length scales, however, remains elusive, and a quantitative assessment of the PVP lifetime under HTHA from the available creep fracture dataset is still not possible. This work develops a microstructure-informed and micromechanics-based model by incorporating a synergy between hydrogen transport and intergranular-cavity-based fracture process. Based on the available creep lifetime data of C-0.5Mo steels, we are able to calibrate material constitutive parameters, and then conduct nonlinear finite element simulations that reveal a real-time stress-induced hydrogen diffusional transport along grain boundaries, coupled with a microstructure-explicit failure process, from which Nelson curves can be computed. Such failure analyses allow us to delineate two distinct regimes on the Nelson curves, i.e., dislocation-creep-controlled or grain boundary diffusion-assisted cavity growth. More importantly, we found that a small change of the pipe thickness and applied stresses can significantly shift these lifetime curves. However, these two parameters are usually not provided in Nelson curves, thus limiting their usage in material selection and safety design. In conclusion, this discrepancy can clearly be mitigated by extensive parametric studies from our micromechanical modeling/simulation framework.

36 MATERIALS SCIENCE↗

Why Do Oceanic Nonlinearities Contribute Only Weakly to Extreme El Niño Events?

Abstract Extreme El Niño events have outsized global impacts and control the El Niño Southern Oscillation (ENSO) warm/cold phases asymmetries. Yet, a consensus regarding the relative contributions of atmospheric and oceanic nonlinearities to their genesis remains elusive. Here, we isolate the contribution of oceanic nonlinearities by conducting paired experiments forced with opposite wind stress anomalies in an oceanic general circulation model, which realistically simulates extreme El Niño events and oceanic nonlinearities thought to contribute to ENSO skewness (Tropical Instability Waves (TIWs), Nonlinear Dynamical Heating (NDH)). Our findings indicate a weak contribution of oceanic nonlinearities to extreme El Niño events in the eastern Pacific, owing to compensatory effects between lateral (NDH and TIWs) and vertical processes. These results hold across different vertical mixing schemes and modifications of the upper‐ocean heat budget mixed layer criterion. Our study reinforces previous research underscoring the pivotal role of atmospheric nonlinearities in shaping extreme El Niño events.

Geology↗

Probabilistic Power Flow Based on a Gaussian Process Emulator

In this letter, a novel Gaussian process emulator is proposed, for the first time, to conduct the probabilistic power-flow calculation. Based on Bayesian inference, a Gaussian process emulator is trained and served as a nonparametric, reduced-order model of the nonlinear power-flow model. This emulator has allowed us to evaluate the time-consuming power-flow solver at the sampled values with a negligible computational cost. The simulations reveal the excellent performance of this method.

42 ENGINEERING↗

Large-Scale Welding Process Simulation by GPU Parallelized Computing

The computational design of industrially relevant welded structures is extremely time consuming due to coupled physics and high nonlinearity. Previously, most welding distortion and residual stress simulations have been limited to small coupons and reduced order (from three-dimensional [3D] to two-dimensional [2D]), or inherent strain approximations were used for large structures. In this current study, an explicit finite element code based on a graphics processing unit was utilized to perform 3D transient thermomechanical simulation of structural components during welding. Laser brazing of aluminum alloy panels as representative of automotive manufacturing scenarios was simulated to predict out-of-plane distortion under different clamping conditions. The predicted deformation pattern and magnitude were validated by laser scanning data of physical assemblies. In addition, the code was used to investigate residual stresses developed during multipass arc welding of a nuclear industry pressurizer surge nozzle and subsequent welding repair where a 3D simulation was necessary. Taking the experimental data as reference, the 3D model predicted better residual stress distribution than a typical 2D asymmetrical model. Stress evolution in welding repair was also presented and discussed in this study. Furthermore, the efficient numerical model made it feasible to use integrated computational welding engineering to simulate welding processes for large-scale structures.

97 MATHEMATICS AND COMPUTING↗

Microbial Models for Simulating Soil Carbon Dynamics: A Review

Abstract Soils store the largest amount of carbon (C) in the biosphere, and the C pool in soil is critical to the global C balance. Numerous microbial models have been developed over the last few decades to represent microbial processes that regulate the responses of soil organic carbon (SOC) to climate change. However, the representation of microbial processes varies, and how microbial processes are incorporated into SOC models has not been well explored. Here, we reviewed 71 microbial models to characterize the microbial processes incorporated into SOC models and analyzed variations in mechanistic complexity. We revealed that (a) four processes (microbial‐mediated decomposition, mineral interaction, microbial necromass recycling, and active and dormant microbial dynamics) are commonly incorporated in microbial models, (b) ∼48% of models simulate only one microbial process (i.e., microbial‐mediated decomposition) and 35% of models simulate two microbial processes: for example, microbial‐mediated decomposition and mineral interaction, (c) more than 80% microbial models use nonlinear equations, such as forward Michaelis‐Menten kinetics, to represent SOC decomposition, (d) the concept of persistence of SOC due to its intrinsic properties has been replaced by organo‐mineral interaction (∼39% of microbial models) that protects SOC from decomposition, and (e) various temperature and moisture modifiers and pH effects have been used to explain the environmental effect on microbial processes. In the future, to realistically incorporate microbial processes into Earth System Models, it is imperative to identify experimental evidence on rate limitation processes and firmly ground model structure on the field and laboratory data.

Environmental Sciences & Ecology↗

The role of collective elasticity on activated structural relaxation, yielding, and steady state flow in hard sphere fluids and colloidal suspensions under strong deformation

In this work, we theoretically study the effect of external deformation on activated structural relaxation and aspects of the nonlinear mechanical response of glassy hard sphere fluids in the context of elastically collective nonlinear Langevin equation theory. This microscopic force-based approach describes activated relaxation as a coupled local–nonlocal event involving caging and longer range collective elasticity, with the latter becoming more important and ultimately dominant with increasing packing fraction under equilibrium conditions. The central new question we address is how this physical picture of activated relaxation, and the relative importance of local caging vs collective elasticity physics, depends on external deformation. Theoretical predictions are presented for deformation-induced enhancement of mobility, the onset of relaxation speed up at remarkably low values of stress, strain, or shear rate, apparent power law thinning of the steady state structural relaxation time and viscosity, a non-vanishing activation barrier in the shear thinning regime, an apparent Herschel–Bulkley form of the rate dependence of the steady state shear stress, exponential growth of different measures of a dynamic yield or flow stress with the packing fraction, and reduced fragility and dynamic heterogeneity under deformation. The results are contrasted with experiments and simulations, and qualitative or better agreement is found. An overarching conclusion is that deformation strongly reduces the importance of longer range collective elastic effects relative to the local caging aspect for most, but not all, physical questions, with deformation-dependent fragility and dynamic heterogeneity phenomena being qualitatively sensitive to collective elasticity. Overall, nonlinear rheology is predicted to be a more local problem than quiescent structural relaxation, albeit with deformation-modified activated processes still important.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Sparse regression for plasma physics

Many scientific problems can be formulated as sparse regression, i.e., regression onto a set of parameters when there is a desire or expectation that some of the parameters are exactly zero or do not substantially contribute. This includes many problems in signal and image processing, system identification, optimization, and parameter estimation methods such as Gaussian process regression. Sparsity facilitates exploring high-dimensional spaces while finding parsimonious and interpretable solutions. In the present work, we illustrate some of the important ways in which sparse regression appears in plasma physics and point out recent contributions and remaining challenges to solving these problems in this field. Further, a brief review is provided for the optimization problem and the state-of-the-art solvers, especially for constrained and high-dimensional sparse regression.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Towards a Robust Adaptive Digital Twin for Fusion Applications

The development of a digital twin system for fusion applications is essential for enhancing the prediction, analysis, and optimization of complex plasma processes. Machine learning (ML), particularly deep learning has demonstrated strong capabilities in modeling such highly nonlinear and intricate systems. However, two critical challenges limit the deployment of deep learning-based digital twins: Uncertainty Quantification (UQ) and data drift. UQ is vital for ensuring trustworthy predictions, especially in decision-support scenarios. Additionally, data-driven models are often sensitive to changes in the underlying data distribution, such as shot-to-shot variations in fusion experiments, which can lead to performance degradation over time. To address these challenges, we are developing an uncertainty-aware, adaptive digital twin framework. Our approach incorporates deep learning models enhanced with Gaussian Process approximations for predictive uncertainty estimation, coupled with an online learning mechanism that enables continuous model adaptation to new experimental data. This adaptive capability allows the data driven models to respond effectively to evolving plasma behaviors and equipment conditions. Specifically, to mitigate the effects of shot-to-shot drift, our system updates itself incrementally as new data becomes available, improving both robustness and fidelity. Our vision is to evolve this data driven model into a self-sustaining digital twin system that leverages UQ based feedback to continuously refine itself and potentially support real-time decision making. This presentation will cover a brief background on uncertainty quantification for ML, our ongoing effort on development of UQ capabilities for ML, our data science pipeline from data collection to model development and analysis and online learning framework for modeling coil deflection at DIII-D. I will also briefly touch upon opportunities and challenges in development of digital twin framework.

Sammuli, Brian [General Atomics]↗

Wavefront distortion correction in scanning tunneling microscope image

Here, we report an algorithm to identify and correct distorted wavefronts in atomic resolution scanning tunneling microscope images. This algorithm can be used to correct nonlinear in-plane distortions without prior knowledge of the physical scanning parameters, the characteristics of the piezoelectric actuator, or individual atom positions. The 2D image is first defined as a sum of sinusoidal plane waves, where a nonlinear distortion renders a curve for an otherwise ideal linear wavefront. Using the Fourier transforms of local areas of the image, the algorithm generates a wavefront vector field. The identified wavefronts are subsequently linearized for each plane wave without changing lattice orders, giving rise to distortion corrections. Our algorithm is complementary to conventional post-processing algorithms that require prior detection of real space features, which can also be used to correct nonlinear distortions in 2D images acquired by other microscopy techniques.

47 OTHER INSTRUMENTATION↗