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

Enhancing the Range and Reliability of the Spacer Layer Imaging Method

The spacer layer imaging method (SLIM) is widely used to measure the thickness of additive and lubricant films, in lubricant development and evaluation, and for fundamental research into elastohydrodynamic lubrication and tribofilm formation mechanisms. The film thickness measurement, as implemented on several popular tribometers, provides powerful, non-destructive in-situ mapping of film topography with nanometre-scale height sensitivity. However, the results can be highly sensitive to experimental procedure, machine condition, and image analysis, in some cases reporting unphysical film thickness trends. The prevailing image analysis techniques make it challenging to interrogate these errors, often hiding their multivariate nonlinear behaviour from the user by spatial averaging. Herein, several common ‘silent errors’ in the SLIM measurement, including colour matching to incorrect fringe orders, and colour drift due to the optical properties of the system or film itself, are discussed, with examples. A robust suite of novel a priori and a posteriori methods to address these issues, and to improve the accuracy and reliability of the measurement, are also presented, including a novel, computationally inexpensive circle-finding algorithm for automated image processing. In combination, these methods allow reliable mapping of films up to at least 800 nm in thickness, representing a significant milestone for the utility of SLIM applied to elastohydrodynamic contact.

EHL film geometry↗
DOE OSTI2022

Nonlinear harmonics coupled by parallel wave propagations in a time-dependent plasma flow

In a time-dependent flow, nonlinear harmonics can be excited by coupling between linear waves and flow-induced harmonic waves. Examining the dispersion relations and selection rules for the coupling, we investigate nonlinearly coupled harmonics for waves propagating along the magnetic field line in a magnetized plasma, as well as waves in an unmagnetized plasma. The coupled harmonics in a plasma flow are described by analytic dispersion relations and selection rules. This nonlinear coupling is corroborated by the particle-in-cell simulation. Here, the coupled-harmonics model describes a mechanism for the excitation of nonlinear harmonics from linear waves in a time-dependent flow. The spectral analysis of the dispersion relation provides a useful way to evaluate the spatiotemporal behavior of a plasma flow.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Convergence of Weak-SINDy Surrogate Models

In this paper, we give an in-depth error analysis for surrogate models generated by a variant of the Sparse Identification of Nonlinear Dynamics (SINDy) method. We start with an overview of a variety of nonlinear system identification techniques, namely SINDy, weak-SINDy, and the occupation kernel method. Under the assumption that the dynamics are a finite linear combination of a set of basis functions, these methods establish a linear system to recover coefficients. We illuminate the structural similarities between these techniques and establish a projection property for the weak-SINDy technique. Following the overview, we analyze the error of surrogate models generated by a simplified version of weak-SINDy. In particular, under the assumption of boundedness of a composition operator given by the solution, we show that (i) the surrogate dynamics converges towards the true dynamics and (ii) the solution of the surrogate model is reasonably close to the true solution. Finally, as an application, we discuss the use of a combination of weak-SINDy surrogate modeling and proper orthogonal decomposition (POD) to build a surrogate model for partial differential equations (PDEs).

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Xyce(™) Parallel Electronic Simulator v.7.5

The Xyce Parallel Electronic Simulator simulates electronic circuit behavior in DC, AC, HB, MPDE and transient mode using standard analog (DAE) and/or device (PDE) device models including several age and radiation aware devices. It supports a variety of computing platforms (both serial and parallel) computers. Lastly, it uses a variety of modern solution algorithms dynamic parallel load-balancing and iterative solvers.! ! Xyce is primarily used to simulate the voltage and current behavior of a circuit network (a network of electronic devices connected via a conductive network). As a tool, it is mainly used for the design and analysis of electronic circuits.! ! Kirchoff's conservation laws are enforced over a network using modified nodal analysis. This results in a set of differential algebraic equations (DAEs). The resulting nonlinear problem is solved iteratively using a fully coupled Newton method, which in turn results in a linear system that is solved by either a standard sparse-direct solver or iteratively using Trilinos linear solver packages, also developed at Sandia National Laboratories.

Source record↗
DOE OSTI2021

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

Tensile-stress effect on ferroelectric Barkhausen noise

Here, this study examines the effect of tensile stress on the ferroelectric properties of Pb(Zr 0.4 Ti 0.6 )O 3 thin film, with a focus on Barkhausen noise, observed for the first time under such conditions. Tensile stress significantly alters domain wall motions, affecting Barkhausen noise more than average polarization. Frequency analysis identifies grain boundaries as primary pinning sites, consistent across stress levels. A nonlinear relationship between stress, domain wall mobility, and polarization is found, where increased stress initially enhances pinning and polarization changes, but this effect diminishes at higher stress levels, indicating a shift in behavior.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Nature and Scalings of Density Fluctuations of Compressible Magnetohydrodynamic Turbulence with Applications to the Solar Wind

Abstract The solar wind is a magnetized and turbulent plasma. Its turbulence is often dominated by Alfvénic fluctuations and often deemed as nearly incompressible far away from the Sun, as shown by in situ measurements near 1 au. However, for solar wind closer to the Sun, the plasma β decreases (often lower than unity) while the turbulent Mach number M t increases (can approach unity, e.g., transonic fluctuations). These conditions could produce significantly more compressible effects, characterized by enhanced density fluctuations, as seen by several space missions. In this paper, a series of 3D MHD simulations of turbulence are carried out to understand the properties of compressible turbulence, particularly the generation of density fluctuations. We find that, over a broad range of parameter space in plasma β , cross helicity, and polytropic index, the turbulent density fluctuations scale linearly as a function of M t , with the scaling coefficients showing weak dependence on parameters. Furthermore, through detailed spatiotemporal analysis, we show that the density fluctuations are dominated by low-frequency nonlinear structures, rather than compressible MHD eigenwaves. These results could be important for understanding how compressible turbulence contributes to solar wind heating near the Sun.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Experimental Limit on Nonlinear State-Dependent Terms in Quantum Theory

Linear time evolution is one of the fundamental postulates of quantum theory. Past theoretical attempts to introduce nonlinearity into quantum evolution have violated causality. However, a recent theory has introduced nonlinear state-dependent terms in quantum field theory, preserving causality [D. E. Kaplan and S. Rajendran, Phys. Rev. D 105, 055002 (2022)]. We report the results of an experiment that searches for such terms. Our approach, inspired by the Everett many-worlds interpretation of quantum theory, correlates a binary macroscopic classical voltage with the outcome of a projective measurement of a quantum bit, prepared in a coherent superposition state. Measurement results are recorded in a bit string, which is used to control a voltage switch. Presence of a nonzero voltage reading in cases of no applied voltage is the experimental signature of a nonlinear state-dependent shift of the electromagnetic field operator. We implement blinded measurement and data analysis with three control bit strings. Control of systematic effects is realized by producing one of the control bit strings with a classical random-bit generator. The other two bit strings are generated by measurements performed on a superconducting qubit in an IBM Quantum processor and on a N 15 nuclear spin in a nitrogen-vacancy center in diamond. Our measurements find no evidence for electromagnetic quantum state-dependent nonlinearity. We set a bound on the parameter that quantifies this nonlinearity | ε γ | < 4.7 × 10 - 11 , at 90% confidence level.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2023

Growth and Saturation of the Electron Drift Instability in a Crossed Field Plasma

The linear growth and nonlinear energy transfer of the electron drift instability (EDI) are experimentally measured in the plume of a low-temperature, Hall effect discharge. A frequency-based bispectral analysis technique applied to fast ion density fluctuation measurements shows a growth rate function that is qualitatively similar to predictions from the linear instability dispersion relation, but an order of magnitude smaller. Calculation of the nonlinear transfer function indicates multiple three-wave interactions between high-frequency resonances of the instability in addition to an inverse energy cascade toward lower-frequency modes. Furthermore, these results are discussed in the context of recent theoretical, numerical, and experimental efforts on the EDI in Hall effect discharges and how the EDI may impact anomalous cross field transport.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Hybrid Bismuth Halide with Rich Polymorphism and Second Harmonic Generation Response

Hybrid structures have emerged as a promising class of optical materials, due to their ability to couple the robustness of inorganic and the tunability of organic compounds. However, their application in nonlinear optics (NLO) remains limited, largely due to underexplored factors that drive the formation of noncentrosymmetric structures and NLO property characterization. In this work, we explore the formation, structural and temperature polymorphism, and optical properties of the (Et 3 NH) 3 Bi 2 Br 9 composition, which crystallizes as either noncentrosymmetric or centrosymmetric polymorph. Structural analysis showed that the alignment of [Bi 2 Br 9 ] 3– units dictates the symmetry of phases, as well as the nonlinear optical properties, with the triclinic polymorph exhibiting a second harmonic generation (SHG) response both in visible and IR regions (1.29 × KH 2 PO 4 and 0.08 × AgGaS 2 ). Thermal analysis reveals polymorphic phase transitions and low melting points, making them melt-processable and ionic liquid candidates.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Theoretical analysis of the saturation phase of the 1/1 energetic-ion-driven resistive interchange mode

The aim of the present study is to analyze the saturation regime of the energetic-ion-driven resistive interchange mode (EIC) in the LHD plasma. In this work, a set of nonlinear simulations are performed by the FAR3d code that uses a reduced MHD model for the thermal plasma coupled with a gyrofluid model for the energetic particle (EP) species. The hellically trapped EP component is introduced through a modification of the averaged drift velocity operator to include their precessional drift. The nonlinear simulation results show similar 1/1 EIC saturation phases with respect to the experimental observations, reproducing the enhancement of the n/m = 1/1 resistive interchange modes (RIC) amplitude and width as the EP β increases, the EP β threshold for the 1/1 EIC excitation, the further destabilization of the 1/1 EIC as the population of the helically trapped EP increases and the triggering of burst events. The frequency of the 1/1 EIC calculated during the burst event is 9.4 kHz and the 2/2 and 3/3 overtones are destabilized, consistent with the frequency range and the complex mode structure measured in the experiment. In addition, the simulation shows the inward propagation of the 1/1 EIC due to the nonlinear destabilization of the 3/4 and 2/3 energetic particle modes, leading to the partial overlapping between resonances during the burst event. Finally, the analysis of the 1/1 EIC stabilization phase shows the excitation of the 1/1 RIC as soon as the flattening induced by the 1/1 EIC in the pressure profile vanishes, leading to the retrieval of the pressure gradient at the plasma periphery and the overcoming of the RIC stability limit.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Improved bound on nonlinear quantum mechanics using a cryogenic radio frequency experiment

There are strong arguments that quantum mechanics may be nonlinear in its dynamics. A discovery of nonlinearity would hint at a novel understanding of the interplay between gravity and quantum field theory, for example. As such, experiments searching for potential nonlinear effects in the electromagnetic sector are important. Here, in this study, we outline such an experiment, consisting of a stream of random bits (which were generated using Rigetti’s Aspen-M-3 chip) as input to an rf signal generator coupled to a cryogenic detector. Projective measurements of the qubit state, which is originally prepared in an equal superposition, serve as the random binary output of a signal generator. Thereafter, spectral analysis of the rf detector would yield a detectable excess signal predicted to arise from such a nonlinear effect. A comparison between the projective measurements of the quantum bits vs the classical baseline showed no power excess. This sets a new limit on the electromagnetic nonlinearity parameter |ε| ⪅ 1.15 x 10 -12 , at a 90.0% confidence level. This is the most stringent limit on nonlinear quantum mechanics thus far and an improvement by nearly a factor of 50 over the previous experimental limit.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Analytical modeling of the evolution of the nonlinearity parameter of sensitized stainless steel

Austenitic stainless steels are subject to the precipitation of chromium carbides (M 23 C 6 ) during exposure to high temperatures, causing these alloys to be susceptible to intergranular corrosion due to chromium depletion along grain boundaries. The acoustic nonlinearity parameter, β, shows sensitivity to the formation of carbides in these alloys. The Thermo-Calc TC-PRISMA module was used to model the nucleation and growth of grain boundary M 23 C 6 carbides. The model was verified with scanning transmission electron microscopy analysis that allowed measurements of the grain boundary precipitates. The paper introduces a reduced-order model of the acoustic nonlinearity based on the formation of misfit dislocations at the interface of the grain boundary precipitate and matrix to explain the change in β during isothermal aging. Furthermore, a direct relationship between the radius of the M 23 C 6 grain boundary carbides and β was observed and verified with nonlinear ultrasound measurements on 304L and 316L stainless steels.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Observations of modal coupling due to bolted joints in an experimental benchmark structure

The goal of this paper is to present a set of measurements from a benchmark structure containing two bolted joints to support future efforts to predict the damping due to the joints and to model nonlinear coupling between the first two elastic modes. Bolted joints introduce nonlinearities in structures, typically causing a softening in the natural frequency and an increase in damping because of frictional slip between the contact interfaces within the joint. These nonlinearities pose significant challenges when characterizing the response of the structure under a large range of load amplitudes, especially when the modal responses become coupled, causing the effective damping and natural frequency to not only depend on the excitation amplitude of the targeted mode, but also the relative amplitudes of other modes. In this work, two nominally identical benchmark structures, known in some prior works as the S4 beam, are tested to characterize their nonlinear properties for the first two elastic modes. Detailed surface measurements are presented and validated through finite element analysis and reveal distinct contact interactions between the two sets of beams. The free-free test structures are excited with an impact hammer and the transient response is analyzed to extract the damping and frequency backbone curves. A range of impact amplitudes and drive points are used to isolate a single mode or to excite both modes simultaneously. Differences in the nonlinear response correlate with the relative strength of the modes that are excited, allowing one to characterize mode coupling. Each of the beams shows different nonlinear properties for each mode, which is attributed to the different contact pressure distributions between the parts, although the mode coupling relationship is found to be consistent between the two. Lastly, the test data key finding are presented in this paper and the supporting data is available on a public repository for interested researchers.

42 ENGINEERING↗
DOE OSTI2021

Inelastic peridynamic model for molecular crystal particles

The peridynamic theory of solid mechanics is applied to modeling the deformation and fracture of micrometer-sized particles made of organic crystalline material. A new peridynamic material model is proposed to reproduce the elastic–plastic response, creep, and fracture that are observed in experiments. The model is implemented in a three-dimensional, meshless Lagrangian simulation code. In the small deformation, elastic regime, the model agrees well with classical Hertzian contact analysis for a sphere compressed between rigid plates. Under higher load, material and geometrical nonlinearity is predicted, leading to fracture. Finally, the material parameters for the energetic material CL-20 are evaluated from nanoindentation test data on the cyclic compression and failure of micrometer-sized grains.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Estimation of constituent properties of concrete materials with an artificial neural network based method

Multi-scale models are developed for heterogeneous concrete materials to estimate their macroscopic mechanical properties in terms of micro-structural data. One crucial challenge of those models is the identification of local properties of constituent phases. In this paper, we present an efficient method based on Artificial Neural Networks (ANN). Typical concrete materials are taken as example. A macroscopic analytical strength criterion is established from three steps of nonlinear homogenization procedure. The macroscopic strength of materials is determined as a function of the frictional coefficient and cohesion of solid cement particles at nanometer scale, intra-particle pores, inter-particle pores and aggregates (inclusions). The objective is to identify the nanoscopic frictional coefficient and cohesion of cement particle from measured macroscopic values of uniaxial compression and tensile strengths. For this purpose, a numerical method based on the ANN is developed. With the analytical macroscopic strength criterion, sensitivity studies are first realized to identify the most important micro-structural parameters influencing the macroscopic strength of concrete. A simplified analytical macroscopic strength criterion is then proposed. A large dataset is further constructed through the inversion of the analytical strength criterion by using the aggregates volume fraction, porosity, macroscopic uniaxial tensile and compressive strengths as input variables and the frictional coefficient and cohesion of cement particles as output unknowns. An ANN model containing four hidden layers and 100 neurons in each layer is constructed and trained by using this dataset. Various types of validation of the ANN model are performed. It is found that the proposed ANN based model can effectively predict the frictional coefficient and cohesion of porous cement paste at the microscopic scale with a very good accuracy.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Mathematical Programming Models for Shale Oil & Gas Development: A Review and Perspective

Here, in this paper, we provide a comprehensive review of mathematical programming models for shale oil & gas development, and we offer a perspective on outstanding research opportunities. We distinguish contributions in five major topic areas, namely: (1) development planning, (2) water management, (3) production optimization, (4) supplies, gathering & processing, and (5) life cycle analysis & sustainability. We highlight how various types of mathematical programming models (i.e., linear programs, nonlinear programs, mixed-integer linear programs, mixed-integer nonlinear programs) have been proposed primarily by the Process Systems Engineering community to address the respective decision-making problems, and we highlight instances of successful deployment in industry. Finally, based on a critical assessment of the existing body of work, we identify opportunities for future research across the major topic areas.

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

Three-wave coupling observed between a shear Alfvén wave and a kink-unstable magnetic flux rope

Results from a laboratory experiment are presented in which, for the first time, a shear Alfvén wave is launched using an antenna in a current-carrying plasma column that is tailored to be either stable or unstable to the kink oscillation. As the plasma is driven kink unstable, the frequency power spectrum of the Alfvén wave evolves from a single peak to a peak with multiple sidebands separated by integer multiples of the kink frequency. The main sidebands (one on either side of the launched wave peak in the power spectrum) are analyzed using azimuthal wavenumber matching, perpendicular and parallel wavenumber decomposition, and bispectral time series analysis. The dispersion relation and three-wave matching conditions are satisfied, given each sideband is a propagating Alfvén wave that results from the interaction of the pump Alfvén wave and the co-propagating component of a half-wavelength, standing kink mode. The interaction is shown to generate smaller perpendicular wavelength Alfvén waves that drive energy transport to scales that will approach the dissipation scale of k⊥ρs=1, with k⊥ being the perpendicular wavenumber and ρs being the ion gyroradius at the electron temperature.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗