Search NASA⌕ Search

SEARCH · Search NASA

Results for “nonlinear experimental dynamics”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 577 records · Page 32
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Dynamics of coupled thermoacoustic modes: Noise and frequency spacing effects

Practical combustion systems consist of multiple thermoacoustic modes that are linearly stable or unstable at a wide range of frequencies. When coupled, the modal interactions are a strong function of the frequency spacing between them – having very significant implications on interpretation of linear stability or instability of modes from experimental data. Frequency spacing between the modes particularly influences the stability and existence of potential limit cycle solutions. In this paper, we extend earlier deterministic studies to include the effects of noise, studying how deterministic dynamics change with the addition of noise and the simultaneous impact of frequency spacing (i.e., closely, or widely spaced). Specifically, this paper characterizes three key effects of noise, which are (1) distribution of limit cycle amplitudes (“diffusion”), (2) shift in the average limit cycle amplitudes (“noise induced drift”), and (3) introduction of new bifurcations that are not present in the deterministic system. Key conclusions from this work are that in highly noisy practical environments, a deterministic/ low noise analysis will not be sufficient to explain even critical qualitative features of the system – such as existence of stable limit cycles or their amplitudes. Furthermore, the nature of noise-induced effects on nonlinear modal interactions can be fundamentally different for degenerate, closely spaced, and widely spaced modes.

42 ENGINEERING↗
NASA NTRS1986

Anisotropic constitutive model for nickel base single crystal alloys: Development and finite element implementation

A tool for the mechanical analysis of nickel base single crystal superalloys, specifically Rene N4, used in gas turbine engine components is developed. This is achieved by a rate dependent anisotropic constitutive model implemented in a nonlinear three dimensional finite element code. The constitutive model is developed from metallurigical concepts utilizing a crystallographic approach. A non Schmid's law formulation is used to model the tension/compression asymmetry and orientation dependence in octahedral slip. Schmid's law is a good approximation to the inelastic response of the material in cube slip. The constitutive equations model the tensile behavior, creep response, and strain rate sensitivity of these alloys. Methods for deriving the material constants from standard tests are presented. The finite element implementation utilizes an initial strain method and twenty noded isoparametric solid elements. The ability to model piecewise linear load histories is included in the finite element code. The constitutive equations are accurately and economically integrated using a second order Adams-Moulton predictor-corrector method with a dynamic time incrementing procedure. Computed results from the finite element code are compared with experimental data for tensile, creep and cyclic tests at 760 deg C. The strain rate sensitivity and stress relaxation capabilities of the model are evaluated.

Dame, L. T.↗
DOE OSTI2025

Enhanced shear stabilization of turbulence in NSTX

In studying a particular non-stationary NSTX L-mode plasma, we observed unexpectedly high levels of flux—first with quasilinear (TGLF) modeling and subsequently with nonlinear gyrokinetic simulations. Upon more detailed analysis, a novel confinement regime was discovered in which a modest increase in E x B shear (beyond baseline experimental estimates) rapidly reduced turbulent transport to levels consistent with power balance. This modest increase is plausible given the errors inherent to the estimation of shearing rates, and the added complexity of the non-stationary (time-dependent) power balance. Remarkably, an additional small increase in shear yields the familiar ion-neoclassical transport level with what appears to be the onset of high-k electron transport only. Although analyses using the TGLF-SAT2 model successfully capture numerous parametric dependencies of this plasma, TGLF does not reproduce the rapid E x B stabilization seen in CGYRO. We believe the results presented should help to better characterize the nonlinear physics of spherical tokamak confinement regimes, provide useful ST datasets for reduced model development, and motivate more accurate experimental diagnosis of E x B shearing rates.

Atomic and molecular collisions↗
NASA NTRS2018

kL-based Linear and Nonlinear Two-Equation Turbulence Models

The development and implementation of kL-based Reynolds Average Navier-Stokes (RANS) two-equation turbulence models are reported herein. The kL is based on Abdol-Hamid's closure and Menter's modification to Rotta's two-equation model. Rotta shows that a reliable transport equation can be formed from the turbulent length scale L, and the turbulent kinetic energy k. Rotta's kL equation is well suited for term-by-term modeling and displays useful features compared to other scale formulation. One of the important differences is the inclusion of higher order velocity derivatives in the source terms of the scale equation. This can enhance the ability of RANS solvers to simulate unsteady flows in URANS mode. The present report documents the formulation of two model levels of turbulence models as implemented in the computational fluid dynamics FUN3D code. The levels are the two-equation linear k-kL and the two-equation algebraic Reynolds stress model (ARSM). Free shear, separated and corner flow cases are documented and compared with experimental, and other turbulence model data. The results show generally very good comparisons with experimental data. The results from this formulation are similar or better than results using the SST two-equation turbulence model. ARSM shows great promise with similar level of computational resources as basic two-equation turbulence models.

Abdol-Hamid, K. S.↗
NASA NTRS2006

Experimental and Computational Analysis of Unidirectional Flow Through Stirling Engine Heater Head

A high efficiency Stirling Radioisotope Generator (SRG) is being developed for possible use in long-duration space science missions. NASA s advanced technology goals for next generation Stirling convertors include increasing the Carnot efficiency and percent of Carnot efficiency. To help achieve these goals, a multi-dimensional Computational Fluid Dynamics (CFD) code is being developed to numerically model unsteady fluid flow and heat transfer phenomena of the oscillating working gas inside Stirling convertors. In the absence of transient pressure drop data for the zero mean oscillating multi-dimensional flows present in the Technology Demonstration Convertors on test at NASA Glenn Research Center, unidirectional flow pressure drop test data is used to compare against 2D and 3D computational solutions. This study focuses on tracking pressure drop and mass flow rate data for unidirectional flow though a Stirling heater head using a commercial CFD code (CFD-ACE). The commercial CFD code uses a porous-media model which is dependent on permeability and the inertial coefficient present in the linear and nonlinear terms of the Darcy-Forchheimer equation. Permeability and inertial coefficient were calculated from unidirectional flow test data. CFD simulations of the unidirectional flow test were validated using the porous-media model input parameters which increased simulation accuracy by 14 percent on average.

Wilson, Scott D.↗
NASA NTRS1989

Experimental investigation of the hydrodynamic forces on the shroud of a centrifugal pump impeller

Fluid-induced forces acting on a rotating impeller are known to cause rotor-dynamic problems in turbomachines. The forces generated by leakage flow along the front shroud surface of a centrifugal turbomachine impeller play an important role among these fluid-induced forces. The present research was aimed to gain a better understanding of these shroud forces. An experimental apparatus was designed and constructed to simulate the impeller shroud leakage flow. Hydrodynamic forces and steady and unsteady pressure distributions on the rotating shroud were measured as functions of eccentricity, width of shroud clearance, face seal clearance and shaft rotating speed. The forces measured from the dynamometer and manometers agreed well. The hydrodynamic force matrices were found skew-symmetric and statically unstable. This is qualitatively similar to the result of previous hydrodynamic volute force measurements. Nondimensionalized normal and tangential forces decrease slightly as Reynolds number increases. As the width of the shroud clearance decreases and/or the eccentricity increases, the hydrodynamic forces increase nonlinearly. There was some evidence found that increased front seal clearance could reduce the radial shroud forces and the relative magnitude of the destabilizing tangential force. Subharmonic pressure fluctuations were also observed which may adversely affect the behavior of the rotor system.

Zhuang, Fei↗
NASA NTRS2002

Parametric Identification of Nonlinear Dynamical Systems

In this project, we looked at the application of harmonic balancing as a tool for identifying parameters (HBID) in a nonlinear dynamical systems with chaotic responses. The main idea is to balance the harmonics of periodic orbits extracted from measurements of each coordinate during a chaotic response. The periodic orbits are taken to be approximate solutions to the differential equations that model the system, the form of the differential equations being known, but with unknown parameters to be identified. Below we summarize the main points addressed in this work. The details of the work are attached as drafts of papers, and a thesis, in the appendix. Our study involved the following three parts: (1) Application of the harmonic balance to a simulation case in which the differential equation model has known form for its nonlinear terms, in contrast to a differential equation model which has either power series or interpolating functions to represent the nonlinear terms. We chose a pendulum, which has sinusoidal nonlinearities; (2) Application of the harmonic balance to an experimental system with known nonlinear forms. We chose a double pendulum, for which chaotic response were easily generated. Thus we confronted a two-degree-of-freedom system, which brought forth challenging issues; (3) A study of alternative reconstruction methods. The reconstruction of the phase space is necessary for the extraction of periodic orbits from the chaotic responses, which is needed in this work. Also, characterization of a nonlinear system is done in the reconstructed phase space. Such characterizations are needed to compare models with experiments. Finally, some nonlinear prediction methods can be applied in the reconstructed phase space. We developed two reconstruction methods that may be considered if the common method (method of delays) is not applicable.

Feeny, Brian↗
NASA NTRS2000

Kinetics of DSB rejoining and formation of simple chromosome exchange aberrations

PURPOSE: To investigate the role of kinetics in the processing of DNA double strand breaks (DSB), and the formation of simple chromosome exchange aberrations following X-ray exposures to mammalian cells based on an enzymatic approach. METHODS: Using computer simulations based on a biochemical approach, rate-equations that describe the processing of DSB through the formation of a DNA-enzyme complex were formulated. A second model that allows for competition between two processing pathways was also formulated. The formation of simple exchange aberrations was modelled as misrepair during the recombination of single DSB with undamaged DNA. Non-linear coupled differential equations corresponding to biochemical pathways were solved numerically by fitting to experimental data. RESULTS: When mediated by a DSB repair enzyme complex, the processing of single DSB showed a complex behaviour that gives the appearance of fast and slow components of rejoining. This is due to the time-delay caused by the action time of enzymes in biomolecular reactions. It is shown that the kinetic- and dose-responses of simple chromosome exchange aberrations are well described by a recombination model of DSB interacting with undamaged DNA when aberration formation increases with linear dose-dependence. Competition between two or more recombination processes is shown to lead to the formation of simple exchange aberrations with a dose-dependence similar to that of a linear quadratic model. CONCLUSIONS: Using a minimal number of assumptions, the kinetics and dose response observed experimentally for DSB rejoining and the formation of simple chromosome exchange aberrations are shown to be consistent with kinetic models based on enzymatic reaction approaches. A non-linear dose response for simple exchange aberrations is possible in a model of recombination of DNA containing a DSB with undamaged DNA when two or more pathways compete for DSB repair.

DNA Damage↗
DOE OSTI2023

Multi-fidelity physics-informed machine learning for probabilistic damage diagnosis

Machine learning (ML) models are gaining popularity in structural health monitoring (SHM) because of their ability to learn the complex relationship between damage and sensor data. However, the lack of sufficient experimental data for structures with different degrees of damage is a key problem in training ML models for SHM. This problem can be alleviated by using physics-based models to generate the required training data to build physics-informed ML (PIML) models for SHM. However, it takes significant computational effort to perform enough high-fidelity simulations of the diagnostic test. It is thus desirable to know whether the available computational resource budget should be expended on numerous low-fidelity physics simulations, or a small number of high-fidelity simulations, or their combination. In this paper, we investigate this aspect of generating adequate training data for PIML, by constructing multi-fidelity PIML models. We evaluate the performance of several PIML models, trained with different amounts of low-fidelity and high-fidelity data, in locating hidden cracks in concrete structures using a nonlinear dynamics-based diagnosis technique. Here, we find that high-fidelity physics simulations that do not cover the (test and damage) parameter space do not improve the performance of diagnostic PIML models built using data from many low-fidelity physics simulations.

42 ENGINEERING↗
DOE OSTI2025

Detecting isolated resonance curves using fixed frequency voltage control tests

Isolated resonance curves, or isolas, are resonance branches of the harmonically forced system that exist separately from the main nonlinear forced response curve, leading to excessive vibrations. Traditional stepped or swept sine simulations and tests rely on continuation along the frequency parameter, typically resulting in a jump phenomenon along the primary resonance branch, prior to the disconnected isola. The main objective of this research is to propose an approach to identify isolated resonance curves by performing continuation along the input amplitude that initializes the response from a low-amplitude solution in the linear regime. Furthermore, this is achieved with the open-loop fixed frequency voltage control method that continues along the shaker voltage parameter and measures the so-called S-curves, which are theoretically a continuous solution branch that connect to the isola. The methodology is demonstrated on a fixture-wing-pylon assembly with a vibro-impact nonlinearity localized in a pylon subcomponent attachment. Multi-harmonic balance simulations are deployed to compute both the nonlinear forced response curves and S-curves to demonstrate the isola detection strategy on a reduced-order finite element model of the nonlinear system. Swept sine and fixed frequency voltage control tests are then conducted on the physical structure to demonstrate the isola detection experimentally, revealing the existence of the large amplitude vibrations that are undetected in the forces levels and frequencies measured with traditional frequency sweeping.

Characterization and Analytical Technique↗
NASA NTRS2007

Global Simulation of Electromagnetic Ion Cyclotron Waves

It is very well known that the effects of electromagnetic ion cyclotron (EMIC) waves on ring current (RC) ion and radiation belt (RB) electron dynamics strongly depend on such particle/wave characteristics as the phase-space distribution function, frequency, wave-normal angle, wave energy, and the form of wave spectral energy density. The consequence is that accurate modeling of EMIC waves and RC particles requires robust inclusion of the interdependent dynamics of wave growth/damping, wave propagation, and particles. Such a self-consistent model is being progressively developed by Khazanov et al. This model is based on a system of coupled kinetic equations for the RC and EMIC wave power spectral density along with the ray tracing equations. We will discuss the recent progress in understanding EMIC waves formation mechanisms in the inner magnetosphere. This problem remains unsettled in spite of many years of experimental and theoretical studies. Modern satellite observations by CRRES, Polar and Cluster still do not reveal the whole picture experimentally since they do not stay long enough in the generation region to give a full account of all the spatio-temporal structure of EMIC waves. The complete self-consistent theory taking into account all factors significant for EMIC waves generation remains to be developed. Several mechanisms are discussed with respect to formation of EMIC waves, among them are nonlinear modification of the ionospheric reflection by precipitating energetic protons, modulation of ion-cyclotron instability by long-period (Pc3/4) pulsations, reflection of waves from layers of heavy-ion gyroresonances, and nonlinearities of wave generation process. We show that each of these mechanisms have their attractive features and explains certain part experimental data but any of them, if taken alone, meets some difficulties when compared to observations. We conclude that development of a refined nonlinear theory and further correlated analysis of modern satellite and ground-based data is needed to solve this very intriguing problem.

Khazanov, George V.↗
NASA NTRS2020

Quasi-Static and Dynamic Analysis of Composite Panels: Challenges in Engineering-Level Commercial Software

In engineering applications, where analysis time and efficiency are essential, analytical models and commercially available finite element (FE) material models are the two mainstream methods used to analyze composite crush scenarios. A reliable numerical model addressing the crushing response of composite materials requires a nonlinear progressive damage analysis which includes appropriate failure mechanisms such as matrix cracking, fragmentation, delamination, and fiber breakage. One might be able to address all these interacting failure responses through implementation/validation of a user defined material model or by utilization of a commercially available model, the latter being the most efficient engineering approach. In this paper, the quasi-static and dynamic crushing response of composite specimens, consisting of a C-channel and a corrugated geometry manufactured with carbon fiber/epoxy matrix fabric, are modeled using the progressive composite damage model within ABAQUS. The novel part of this work, which has not been previously addressed in the literature, is the utilization of a detailed parametric study to optimize composite damage model parameters for a given quasi-static crush simulation in ABAQUS. C-channel experiments were used to calibrate the model and then refine/optimize a set of material parameters that best fit the experimental results. These identical parameters were then used to simulate the crush response of the corrugated channel specimens, thereby making it a pure prediction and not a double recalibration. Next, the dynamic response of both panels were predicted, further confirming the predictive capability of the developed numerical model. The effect of parameter variation on the crush response is fully documented for future use in engineering applications utilizing this specific ABAQUS damage model.

Paria Naghipour↗
NASA NTRS2007

Global Simulation of Electromagnetic Ion Cyclotron Waves

It is well known that the effects of electromagnetic ion cyclotron (EMIC) waves on ring current (RC) ion and radiation belt (RB) electron dynamics strongly depend on such particle/wave characteristics as the phase-space distribution function, frequency, wave-normal angle, wave energy, and the form of wave spectral energy density. The consequence is that accurate modeling of EMIC waves and RC particles requires robust inclusion of the interdependent dynamics of wave growth/damping, wave propagation, and particles. Such a self-consistent model is being progressively developed by Khazanov et al. [2002 - 2007]. This model is based on a system of coupled kinetic equations for the RC and EMIC wave power spectral density along with the ray tracing equations. We will discuss the recent progress in understanding EMIC waves formation mechanisms in the inner magnetosphere. This problem remains unsettled in spite of many years of experimental and theoretical studies. Modern satellite observations by CRRES, Polar and Cluster still do not reveal the whole picture experimentally since they do not stay long enough in the generation region to give a full account of all the spatio-temporal structure of EMIC waves. The complete self-consistent theory taking into account all factors significant for EMIC waves generation remains to be developed. Several mechanisms are discussed with respect to formation of EMIC waves, among them are nonlinear modification of the ionospheric reflection by precipitating energetic protons, modulation of ion-cyclotron instability by long-period (Pc3/4) pulsations, reflection of waves from layers of heavy-ion gyroresonances, and nonlinearities of wave generation process. We show that each of these mechanisms have their attractive features and explains certain part experimental data but any of them, if taken alone, meets some difficulties when compared to observations. We conclude that development of a refined nonlinear theory and further correlated analysis of modern satellite and ground-based data is needed to solve this very intriguing problem.

Khazanov, G. V.↗
DOE OSTI2022

Quantifying transport and electrocatalytic reaction processes in a gastight rotating cylinder electrode reactor via integration of Computational Fluid Dynamics modeling and experiments

Understanding the complexity of the multiple processes of mass, momentum, charge, and heat transport, and how these affect reaction kinetics at the electrode/electrolyte interface is one of the major challenges in the field of energy and catalysis. The rapid and rational scale-up of electrocatalytic systems to industrial scales require a detailed understanding of nonlinear transport-reaction processes, accessible only through the building of multi-physics models that capture with high fidelity the complexity of real-world devices. The gastight rotating cylinder electrode (RCE) reactor is a promising lab-scale tool that can decouple transport from intrinsic kinetics to generate data for first-principle models useful in the design of industrial, electrochemical reactors. Computational Fluid Dynamics (CFD) studies have previously been used to investigate the bulk flow in RCE reactors for simple corrosion and electroplating processes. However, the quantification of changes in local concentration within the viscous layer where catalysis takes place requires capturing the correct flow conditions inside the hydrodynamic boundary layer near the surface of the electrode. Further, this requires simulations with spatial resolution in the nm and μm scale and temporal resolutions between ms and s scales that are similar to the timescales for reactions on the electrode surface. In this study, experimental electrocatalysis is combined with CFD modeling to elucidate and parameterize the hydrodynamics in a gastight RCE reactor. CFD simulations of the electrochemical ferricyanide reduction reaction under mass transport limited conditions are used to evaluate the validity of the CFD model parameters by comparing calculated dimensionless mass transport descriptors to dimensionless correlations obtained experimentally. Justifications for assumptions and details of the simulation methods used in this study are presented to provide a detailed understanding of the effect that each model parameter has on the ability to accurately simulate electrocatalysis in RCE systems. The simulation methodology reported here is a first step towards the development of multi-scale models for the study of transport dependent electrocatalytic processes, such as the electrochemical transformation of CO 2 to fuels and chemicals.

42 ENGINEERING↗
NASA NTRS2015

Subharmonic Oscillations and Chaos in Dynamic Atomic Force Microscopy

The increasing use of dynamic atomic force microscopy (d-AFM) for nanoscale materials characterization calls for a deeper understanding of the cantilever dynamics influencing scan stability, predictability, and image quality. Model development is critical to such understanding. Renormalization of the equations governing d- AFM provides a simple interpretation of cantilever dynamics as a single spring and mass system with frequency dependent cantilever stiffness and damping parameters. The renormalized model is sufficiently robust to predict the experimentally observed splitting of the free-space cantilever resonance into multiple resonances upon cantilever-sample contact. Central to the model is the representation of the cantilever sample interaction force as a polynomial expansion with coefficients F(sub ij) (i,j = 0, 1, 2) that account for the effective interaction stiffness parameter, the cantilever-to-sample energy transfer, and the amplitude of cantilever oscillation. Application of the Melnikov method to the model equation is shown to predict a homoclinic bifurcation of the Smale horseshoe type leading to a cascade of period doublings with increasing drive displacement amplitude culminating in chaos and loss of image quality. The threshold value of the drive displacement amplitude necessary to initiate subharmonic generation depends on the acoustic drive frequency, the effective damping coefficient, and the nonlinearity of the cantilever-sample interaction force. For parameter values leading to displacement amplitudes below threshold for homoclinic bifurcation other bifurcation scenarios can occur, some of which lead to chaos.

Cantrell, John H.↗
NASA NTRS2016

Development of Semi-Empirical Damping Equation for Baffled Tank with Oblate Spheroidal Dome

Propellant slosh is a potential source of disturbance that can significantly impact the stability of space vehicles. The slosh dynamics are typically represented by a mechanical model of a spring-mass-damper. This mechanical model is then included in the equation of motion of the entire vehicle for Guidance, Navigation and Control analysis. The typical parameters required by the mechanical model include natural frequency of the slosh, slosh mass, slosh mass center location, and the critical damping ratio. A fundamental study has been undertaken at NASA MSFC to understand the fluid damping physics from a ring baffle in the barrel section of a propellant tank. An asymptotic damping equation and CFD blended equation have been derived by NASA MSFC team to complement the popularly used Miles equation at different flow regimes. The new development has found success in providing a nonlinear damping model for the Space Launch System. The purpose of this study is to further extend the semi-empirical damping equations into the oblate spheroidal dome section of the propellant tanks. First, previous experimental data from the spherical baffled tank are collected and analyzed. Several methods of taking the dome curvature effect, including a generalized Miles equation, area projection method, and equalized fill height method, are assessed. CFD simulation is used to shed light on the interaction of vorticity around the baffle with the locally curved wall and liquid-gas interface. The final damping equation will be validated by a recent subscale test with an oblate spheroidal dome conducted at NASA MSFC.

Yang, H. Q.↗
DOE OSTI2023

Network-Wide Traffic Signal Control Using Bilinear System Modeling and Adaptive Optimization

This study proposes a new multi-input multi-output optimal bilinear signal control method in which a bilinear dynamic model approximation is used to capture the nonlinear dynamics of the urban traffic networks. With signal green time splits as the control input and traffic delay changes as the output for each intersections in the network, a bilinear system model was developed, which, on the basis of linear system modeling, takes interactions among traffic delays and signal timing splits into consideration. Based on the bilinear system modeling framework, we conducted two steps in each time interval to derive traffic control strategies: (1) we used the normalized least-squared algorithm to estimate system parameters; and (2) we solved an online optimization problem to obtain the updated traffic control inputs for the signal timing that minimizes future traffic delays. We evaluated the proposed method in a microscopic traffic simulation environment (VISSIM) with a 35-intersection network of Bellevue city in Washington. Two different traffic demand patterns: (1) normal traffic demands; and (2) time-varying traffic demands were simulated to compare the performance of different control strategies. Experimental results show that (1) the proposed bilinear system model can better describe traffic system dynamics than linear-model based methods, such as our previously developed linear-quadratic regulator control; and (2) the proposed method outperforms the state-of-the-art signal control strategies, namely the max-pressure and the self-organizing traffic light control methods. We have also shown that the proposed method is applicable to all other possible network layouts and signal controller phasing structures.

42 ENGINEERING↗
DOE OSTI2023

Superposition principle and nonlinear response in spin glasses

The extended principle of superposition has been a touchstone of spin-glass dynamics for almost 30 years. Here the Uppsala group has demonstrated its validity for the metallic spin glass, CuMn, for magnetic fields H up to 10 Oe at the reduced temperature T r = T/T g = 0.95, where T g is the spin-glass condensation temperature. For H > 10 Oe, they observe a departure from linear response which they ascribe to the development of nonlinear dynamics. The thrust of this paper is to develop a microscopic origin for this behavior by focusing on the time development of the spin-glass correlation length, ξ(t, t w ; H). Here, t is the time after H changes, and t w is the time from the quench for T > T g to the working temperature T until H changes. We connect the growth of ξ(t, t w ; H) to the barrier heights Δ(t w ) that set the dynamics. The effect of H on the magnitude of Δ (t w ) is responsible for affecting differently the two dynamical protocols associated with turning H off (TRM, or thermoremanent magnetization) or on (ZFC, or zero-field-cooled magnetization). This difference is a consequence of nonlinearity based on the effect of H on Δ (t w ). Superposition is preserved if Δ(t w ) is linear in the Hamming distance Hd (proportional to the difference between the self-overlap q EA and the overlap q[Δ(t w )]). However, superposition is violated if Δ(t w ) increases faster than linear in Hd. We have previously shown, through experiment and simulation, that the barriers Δ(t w ) do increase more rapidly than linearly with Hd through the observation that the growth of ξ(t,t w ;H) slows down as ξ(t,t w ; H) increases. In this paper, we display the difference between the zero-field-cooled ξ ZFC (t, t w ; H) and the thermoremanent magnetization ξ TRM (t, t w ; H) correlation lengths as H increases, both experimentally and through numerical simulations, corresponding to the violation of the extended principle of superposition in line with the finding of the Uppsala Group.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗