Search NASA⌕ Search

SEARCH · Search NASA

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 1,279 records · Page 71
All sourcesNASA NTRSDOE OSTI
NASA NTRS1997

The Stability of Radiatively Cooling Jets I. Linear Analysis

The results of a spatial stability analysis of a two-dimensional slab jet, in which optically thin radiative cooling is dynamically important, are presented. We study both magnetized and unmagnetized jets at external Mach numbers of 5 and 20. We model the cooling rate by using two different cooling curves: one appropriate to interstellar gas, and the other to photoionized gas of reduced metallicity. Thus, our results will be applicable to both protostellar (Herbig-Haro) jets and optical jets from active galactic nuclei. We present analytical solutions to the dispersion relations in useful limits and solve the dispersion relations numerically over a broad range of perturbation frequencies. We find that the growth rates and wavelengths of the unstable Kelvin-Helmholtz (K-H) modes are significantly different from the adiabatic limit, and that the form of the cooling function strongly affects the results. In particular, if the cooling curve is a steep function of temperature in the neighborhood of the equilibrium state, then the growth of K-H modes is reduced relative to the adiabatic jet. On the other hand, if the cooling curve is a shallow function of temperature, then the growth of K-H modes can be enhanced relative to the adiabatic jet by the increase in cooling relative to heating in overdense regions. Inclusion of a dynamically important magnetic field does not strongly modify the important differences between an adiabatic jet and a cooling jet, provided the jet is highly supermagnetosonic and not magnetic pressure-dominated. In the latter case, the unstable modes behave more like the transmagnetosonic magnetic pressure-dominated adiabatic limit. We also plot fluid displacement surfaces associated with the various waves in a cooling jet in order to predict the structures that might arise in the nonlinear regime. This analysis predicts that low-frequency surface waves and the lowest order body modes will be the most effective at producing observable features in the jet.

Hardee, Philip E.↗
NASA NTRS2003

Evaluation and Analysis of F-16XL Wind Tunnel Data from Dynamic Tests

A series of wind tunnel tests were conducted in the NASA Langley Research Center as part of an ongoing effort to develop and test mathematical models for aircraft rigid-body aerodynamics in nonlinear unsteady flight regimes. Analysis of measurement accuracy, especially for nonlinear dynamic systems that may exhibit complicated behaviors, is an essential component of this ongoing effort. In this paper, tools for harmonic analysis of dynamic data and assessing measurement accuracy are presented. A linear aerodynamic model is assumed that is appropriate for conventional forced-oscillation experiments, although more general models can be used with these tools. Application of the tools to experimental data is demonstrated and results indicate the levels of uncertainty in output measurements that can arise from experimental setup, calibration procedures, mechanical limitations, and input errors.

Kim, Sungwan↗
NASA NTRS2004

Evaluation and Analysis of F-16XL Wind Tunnel Data From Static and Dynamic Tests

A series of wind tunnel tests were conducted in the NASA Langley Research Center as part of an ongoing effort to develop and test mathematical models for aircraft rigid-body aerodynamics in nonlinear unsteady flight regimes. Analysis of measurement accuracy, especially for nonlinear dynamic systems that may exhibit complicated behaviors, is an essential component of this ongoing effort. In this report, tools for harmonic analysis of dynamic data and assessing measurement accuracy are presented. A linear aerodynamic model is assumed that is appropriate for conventional forced-oscillation experiments, although more general models can be used with these tools. Application of the tools to experimental data is demonstrated and results indicate the levels of uncertainty in output measurements that can arise from experimental setup, calibration procedures, mechanical limitations, and input errors.

Kim, Sungwan↗
NASA NTRS1982

Three-dimensional finite-element analysis of layered composite plates

Results are presented for an investigation of the three-dimensional, geometrically nonlinear, finite-element analysis of the bending of laminated anisotropic composite plates. The individual laminae are treated as homogeneous, transversely isotropic, and linearly elastic. A fully three-dimensional isoparametric finite element with eight modes (i.e., linear element) and 24 degrees of freedom (three displacement components per node) is used. The numerical results obtained using this linear analysis are compared with the exact solutions given in Pagano (1969, 1970). It is found that the results of the linear analysis converge to the exact solution as the mesh is refined.

Putcha, N. S.↗
DOE OSTI2014

Characterizing friction in sliding isolation bearings

Summary The force–displacement behavior of the Friction Pendulum™ (FP) bearing is a function of the coefficient of sliding friction, axial load on the bearing and effective radius of the sliding surface. The coefficient of friction varies during the course of an earthquake with sliding velocity, axial pressure and temperature at the sliding surface. The velocity and axial pressure on the bearing depend on the response of the superstructure to the earthquake shaking. The temperature at an instant in time during earthquake shaking is a function of the histories of the coefficient of friction, sliding velocity and axial pressure, and the travel path of the slider on the sliding surface. A unified framework accommodating the complex interdependence of the coefficient of friction, sliding velocity, axial pressure and temperature is presented for implementation in nonlinear response‐history analysis. Expressions to define the relationship between the coefficient of friction and sliding velocity, axial pressure, and temperature are proposed, based on available experimental data. Response‐history analyses are performed on FP bearings with a range of geometrical and liner mechanical properties and static axial pressure. Friction is described using five different models that consider the dependence of the coefficient of friction on axial pressure, sliding velocity and temperature. Frictional heating is the most important factor that influences the maximum displacement of the isolation system and floor spectral demands if the static axial pressure is high. Isolation system displacements are not significantly affected by considerations of the influence of axial pressure and velocity on the coefficient of friction. Copyright © 2014 John Wiley & Sons, Ltd.

Kumar, Manish↗
DOE OSTI2024

Proximal Galerkin: A Structure-Preserving Finite Element Method for Pointwise Bound Constraints

The proximal Galerkin finite element method is a high-order, low iteration complexity, nonlinear numerical method that preserves the geometric and algebraic structure of pointwise bound constraints in infinite-dimensional function spaces. This paper introduces the proximal Galerkin method and applies it to solve free boundary problems, enforce discrete maximum principles, and develop a scalable, mesh-independent algorithm for optimal design with pointwise bound constraints. This paper also introduces the latent variable proximal point (LVPP) algorithm, from which the proximal Galerkin method derives. When analyzing the classical obstacle problem, we discover that the underlying variational inequality can be replaced by a sequence of second-order partial differential equations (PDEs) that are readily discretized and solved with, e.g., the proximal Galerkin method. Throughout this work, we arrive at several contributions that may be of independent interest. These include (1) a semilinear PDE we refer to as the entropic Poisson equation; (2) an algebraic/geometric connection between high-order positivity-preserving discretizations and certain infinite-dimensional Lie groups; and (3) a gradient-based, bound-preserving algorithm for two-field, density-based topology optimization. The complete proximal Galerkin methodology combines ideas from nonlinear programming, functional analysis, tropical algebra, and differential geometry and can potentially lead to new synergies among these areas as well as within variational and numerical analysis. Open-source implementations of our methods accompany this work to facilitate reproduction and broader adoption.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Simultaneous Decomposition and Dewetting of Nanoscale Alloys: A Comparison of Experiment and Theory

In this work, we consider the coupled process of phase separation and dewetting of metal alloys of nanoscale thickness deposited on solid substrates. The experiments involve applying nanosecond laser pulses that melt the Ag40Ni60 alloy films in two setups: either on thin supporting membranes or on bulk substrates. These two setups allow for extracting both temporal and spatial scales on which the considered processes occur. The theoretical model involves a longwave version of the Cahn–Hilliard formulation used to describe spinodal decomposition, coupled with an asymptotically consistent longwave-based description of dewetting that occurs due to destabilizing interactions between the alloy and the substrate, modeled using the disjoining pressure approach. Careful modeling, combined with linear stability analysis and fully nonlinear simulations, leads to results consistent with the experiments. In particular, we find that the two instability mechanisms occur concurrently, with the phase separation occurring faster and on shorter temporal scales. The modeling results show a strong influence of the temperature dependence of relevant material properties, implying that such a dependence is crucial for the understanding of the experimental findings. The agreement between theory and experiment suggests the utility of the proposed theoretical approach in helping to develop further experiments directed toward formation of metallic alloy nanoparticles of desired properties.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Morphology and Transport of Multivalent Cation-Exchanged Ionomer Membranes Using Perfluorosulfonic Acid–Ce Z+ as a Model System

Perfluorosulfonic acids (PFSAs) are commonly used as solid polymer electrolyte membranes (PEMs) in electrochemical energy devices, where they are vulnerable to attack by radical species during operation. A popular strategy to combat this problem is to introduce radical scavengers like cerium (Ce) ions that neutralize these species before they attack the PFSA. Such cation doping creates a multi-ion system, in which understanding the mechanisms of cation solvation and transport becomes important for the effective design and utilization of PFSA–cation systems. Ce ions also provide a representative model system for multication-exchanged ionomers in electrochemical systems. In this study, hydration and conductivity measurements, along with X-ray fluorescence and scattering, are employed to elucidate how Ce ion exchange alters PFSA’s ionic solvation, as well as nano- and mesoscale morphologies, which ultimately control its ion transport properties. A molecular transport model is used to deconvolute the impact of Ce ions on the local solvation structure of water in the membrane from mesoscale changes of the transport pathways. The combined experimental and theoretical analysis reveals a nonlinear decrease in conductivity driven by cation solvation at the molecular level and morphological changes at longer length scales. Migration–diffusion coupling, its nonlinear dependence on ion exchange and hydration, and its overall implications for ionomer performance are also discussed. Finally, these findings have the potential to be translated into other mixed cation–ionomer systems for a wide range of energy and environmental devices.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Dynamics from noisy data with extreme timing uncertainty

Imperfect knowledge of the times at which ‘snapshots’ of a system are recorded degrades our ability to recover dynamical information, and can scramble the sequence of events. In X-ray free-electron lasers, for example, the uncertainty—the so-called timing jitter—between the arrival of an optical trigger (‘pump’) pulse and a probing X-ray pulse can exceed the length of the X-ray pulse by up to two orders of magnitude, marring the otherwise precise time-resolution capabilities of this class of instruments. The widespread notion that little dynamical information is available on timescales shorter than the timing uncertainty has led to various hardware schemes to reduce timing uncertainty. These schemes are expensive, tend to be specific to one experimental approach and cannot be used when the record was created under ill-defined or uncontrolled conditions such as during geological events. Here we present a data-analytical approach, based on singular-value decomposition and nonlinear Laplacian spectral analysis, that can recover the history and dynamics of a system from a dense collection of noisy snapshots spanning a sufficiently large multiple of the timing uncertainty. In this work, the power of the algorithm is demonstrated by extracting the underlying dynamics on the few-femtosecond timescale from noisy experimental X-ray free-electron laser data recorded with 300-femtosecond timing uncertainty. Using a noisy dataset from a pump-probe experiment on the Coulomb explosion of nitrogen molecules, our analysis reveals vibrational wave-packets consisting of components with periods as short as 15 femtoseconds, as well as more rapid changes, which have yet to be fully explored. Our approach can potentially be applied whenever dynamical or historical information is tainted by timing uncertainty.

47 OTHER INSTRUMENTATION↗
DOE OSTI2013

Substituent effects in a series of 1,7-C[subscript 60](R[subscript F])[subscript 2] compounds (R[subscript F] = CF[subscript 3], C[subscript 2]F[subscript 5], n-C[subscrip 3]F[subscript 7], i-C[subscript 3]F[subscript 7], n-C[subscript 4]F[subscript 9], s-C[subscript 4]F[subscript 9], n-C[subscript 8]F[subscript 17]): electron affinities, reduction potentials and E(LUMO) values are not always correlated

A series of seven structurally-similar compounds with different pairs of R{sub F} groups were prepared, characterized spectroscopically, and studied by electrochemical methods (cyclic and square-wave voltammetry), low-temperature anion photoelectron spectroscopy, and DFT calculations (five of the compounds are reported here for the first time). This is the first time that a set of seven R{sub F} groups have been compared with respect to their relative effects on E{sub 1/2}(0/-), electron affinity (EA), and the DFT-calculated LUMO energy. The compounds, 1,7-C{sub 60}(R{sub F}){sub 2} (R{sub F} = CF{sub 3}, C{sub 2}F{sub 5}, i-C{sub 3}F{sub 7}, n-C{sub 3}F{sub 7}, s-C{sub 4}F{sub 9}, n-C{sub 4}F{sub 9} and n-C{sub 8}F{sub 21}), were found to have statistically different electron affinities (EA), at the {+-}10 meV level of uncertainty, but virtually identical first reduction potentials, at the {+-}10 mV level of uncertainty. The lack of a correlation between EA and E{sub 1/2}(0/-), and between E(LUMO) and E{sub 1/2}(0/-), for such similar compounds is unprecedented and suggests that explanations for differences in figures of merit for materials and/or devices that are based on equating easily measurable E{sub 1/2}(0/-) values with EAs or E(LUMO) values should be viewed with caution. The solubilities of the seven compounds in toluene varied by nearly a factor of six, but in an unpredictable way, with the C{sub 2}F{sub 5} and s-C{sub 4}F{sub 9} compounds being the most soluble and the i-C{sub 3}F{sub 7} compound being the least soluble. The effects of the different R{sub F} groups on EAs, E(LUMO) values, and solubilities should help fluorine chemists choose the right R{sub F} group to design new materials with improved morphological, electronic, optical, and/or magnetic properties.

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

Role of subdominant stable modes in plasma microturbulence

In gyrokinetic simulations, thousands of degrees of freedom are available to contribute to the fluctuation spectrum. For wavevectors with a single linear instability, the unstable eigenmode accounts for only one of these degrees of freedom. Little has been known about the role of the remaining fluctuations in the turbulent dynamics. In this paper, these fluctuations are characterized as modes in mode decompositions of gyrokinetic distribution functions from nonlinear simulations. This analysis reveals the excitation of a hierarchy of damped modes at the same perpendicular scales as the driving instabilities. Two effects of these subdominant modes are described: First, these damped modes define a potent energy sink, creating a situation where energy drive and energy dissipation peak at the same perpendicular scales. Second, damped modes with tearing parity (even parity about the outboard midplane for A fluctuations) are driven to significant amplitudes and facilitate the development of magnetic stochasticity in electromagnetic gyrokinetic simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

Superconducting cuprate heterostructures for hot electron bolometers

Transport properties of the resistive state of quasi-two dimensional superconducting heterostructures containing ultrathin La{sub 2−x}Sr{sub x}CuO{sub 4} layers synthesized using molecular beam epitaxy are studied. The electron transport exhibits strong deviation from Ohm's law, δV∼γI{sup 3}, with a coefficient γ(T) that correlates with the temperature variation of the resistivity dρ/dT. Close to the normal state, analysis of the nonlinear behavior in terms of electron heating yields an electron-phonon thermal conductance per unit area g{sub e−ph}≈1 W/K cm{sup 2} at T = 20 K, one-two orders of magnitude smaller than in typical superconductors. This makes superconducting LaSrCuO heterostructures to be attractive candidate for the next generation of hot electron bolometers with greatly improved sensitivity.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Decoding the spectroscopic features and time scales of aqueous proton defects

Acid solutions demonstrate a variety of complex structural and dynamical features arising from the presence of multiple interacting reactive proton defects and counterions. Yet, disentangling the transient structural motifs of proton defects in the water hydrogen bond network and the mechanisms for their interconversion remains a formidable challenge. In this work, we use simulations treating the quantum nature of both the electrons and nuclei to show how the experimentally observed spectroscopic features and relaxation time scales can be elucidated using a physically transparent coordinate that encodes the overall asymmetry of the solvation environment of the proton defect. We demonstrate that this coordinate can be used both to discriminate the extremities of the features observed in the linear vibrational spectrum and to explain the molecular motions that give rise to the interconversion time scales observed in recent nonlinear experiments. This analysis provides a unified condensed-phase picture of the proton structure and dynamics that, at its extrema, encompasses proton sharing and spectroscopic features resembling the limiting Eigen [H 3 O(H 2 O) 3 ]+ and Zundel [H(H 2 O) 2 ]+ gas-phase structures, while also describing the rich variety of interconverting environments in the liquid phase.

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

Fluid and kinetic nonlinearities of near-acoustic plasma waves

Large amplitude, near-acoustic plasma waves have a temperature-dependent harmonic content and nonlinear frequency shifts. On cold plasmas, experiments show that fluid nonlinearities dominate, and the observed harmonic content is well predicted by a recent fluid analysis. In contrast, the measured nonlinear frequency shift is smaller than this theory prediction, suggesting the presence of other nonlinear effects neglected in this analysis. When the plasma is heated, kinetic effects become important, and we observe an increase in both the harmonic content and nonlinear frequency shift. Quantitative experiments on the stability of these large amplitude, near-acoustic plasma waves strongly support the idea that trapped particles reduce the amplitude threshold for instability below fluid models. Furthermore, at low temperatures, the broad characteristics of the parametric instability are determined by the frequency detuning of the pump and daughter wave and the wave-wave coupling strength, surprisingly consistent with simple cold fluid, three-wave theories. However, at higher temperatures, the pump wave becomes unstable at half the threshold pump wave amplitude, with similar exponential growth rates as for a cold plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2012

FIELD TOPOLOGY ANALYSIS OF A LONG-LASTING CORONAL SIGMOID

We present the first field topology analysis based on nonlinear force-free field (NLFFF) models of a long-lasting coronal sigmoid observed in 2007 February with the X-Ray Telescope on Hinode. The NLFFF models are built with the flux rope insertion method and give the three-dimensional coronal magnetic field as constrained by observed coronal loop structures and photospheric magnetograms. Based on these models, we have computed horizontal maps of the current and the squashing factor Q for 25 different heights in the corona for all six days of the evolution of the region. We use the squashing factor to quantify the degree of change of the field line linkage and to identify prominent quasi-separatrix layers (QSLs). We discuss the major properties of these QSL maps and devise a way to pick out important QSLs since our calculation cannot reach high values of Q. The complexity in the QSL maps reflects the high degree of fragmentation of the photospheric field. We find main QSLs and current concentrations that outline the flux rope cavity and that become characteristically S-shaped during the evolution of the sigmoid. We note that, although intermittent bald patches exist along the length of the sigmoid during its whole evolution, the flux rope remains stable for several days. However, shortly after the topology of the field exhibits hyperbolic flux tubes (HFT) on February 7 and February 12 the sigmoid loses equilibrium and produces two B-class flares and associated coronal mass ejections (CMEs). The location of the most elevated part of the HFT in our model coincides with the inferred locations of the two flares. Therefore, we suggest that the presence of an HFT in a coronal magnetic configuration may be an indication that the system is ready to erupt. We offer a scenario in which magnetic reconnection at the HFT drives the system toward the marginally stable state. Once this state is reached, loss of equilibrium occurs via the torus instability, producing a CME.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2012

GAS CONDENSATION IN THE GALACTIC HALO

Using adaptive mesh refinement (AMR) hydrodynamic simulations of vertically stratified hot halo gas, we examine the conditions under which clouds can form and condense out of the hot halo medium to potentially fuel star formation in the gaseous disk. We find that halo clouds do not develop from linear isobaric perturbations. This is a regime where the cooling time is longer than the Brunt-Vaeisaelae time, confirming previous linear analysis. We extend the analysis into the nonlinear regime by considering mildly or strongly nonlinear perturbations with overdensities up to 100, also varying the initial height, the cloud size, and the metallicity of the gas. Here, the result depends on the ratio of cooling time to the time required to accelerate the cloud to the sound speed (similar to the dynamical time). If the ratio exceeds a critical value near unity, the cloud is accelerated without further cooling and gets disrupted by Kelvin-Helmholtz and/or Rayleigh-Taylor instabilities. If it is less than the critical value, the cloud cools and condenses before disruption. Accreting gas with overdensities of 10-20 is expected to be marginally unstable; the cooling fraction will depend on the metallicity, the size of the incoming cloud, and the distance to the galaxy. Locally enhanced overdensities within cold streams have a higher likelihood of cooling out. Our results have implications on the evolution of clouds seeded by cold accretion that are barely resolved in current cosmological hydrodynamic simulations and absorption line systems detected in galaxy halos.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Global linear and nonlinear gyrokinetic simulations of tearing modes

Abstract To better understand multi-scale interactions between global tearing modes and microturbulence in the Madison Symmetric Torus (MST) reversed-field pinch (RFP), the global gyrokinetic code Gene is modified to describe global tearing mode instability via a shifted Maxwellian (SM) distribution consistent with experimental equilibria. The implementation of the SM is tested and benchmarked by comparisons with different codes and models. Good agreement is obtained in code-code and code-theory comparisons. Linear stability of tearing modes of a non-reversed MST discharge is studied. A collisionality scan is performed to the lowest order unstable modes ( n = 5, n = 6) and shown to behave consistently with theoretical scaling. The nonlinear evolution is simulated, and saturation is found to arise from mode coupling and transfer of energy from the most unstable tearing mode to small-scale stable modes mediated by the m = 2 tearing mode. The work described herein lays the foundation for nonlinear simulation and analysis of the interaction of tearing modes and gyroradius-scale instabilities in RFP plasmas.

Physics↗
DOE OSTI2024

Prediction of pellet mass thresholds for ELM triggering in low-collisionality, ITER-like discharges

Abstract In ITER, pellets are calculated to require more than 8 times the mass than currently planned to reliably trigger edge-localized modes (ELMs). Unmitigated heat flux impulses from ELMs are intolerable in ITER at full power and current. Therefore, ITER operation relies on multiple approaches to control ELM heat fluxes. One method is pellet ELM pacing to instigate small rapid ELMs with low heat flux. Predicting the performance of pellet pacing is critical for ITER, which is expected to operate in a regime with a low-collisionality, peeling-limited pedestal. However, to trigger ELMs the local pressure increase in the expanding pellet cloud pushes the equilibrium over the ballooning stability limit. In this work, linear and nonlinear M3D-C1 simulations are used to predict pellet mass thresholds in DIII-D discharges and ITER scenarios with peeling-limited pedestals. It is found that the distance of the equilibrium’s operational point from the ballooning branch of the pedestal stability boundary strongly changes thresholds. Linear M3D-C1 simulations find a strong dependence of the pellet mass threshold on the poloidal injection location for ITER’s 15 MA, Q = 10 scenario. The required pellet mass at the planned injection locations is 8 to 17 times larger than currently considered. However, such linear simulations do not include pellet ablation physics or time evolution of density and temperature. A new scheme of 2D nonlinear simulations, coupled with linear stability analysis at various steps throughout the nonlinear time evolution, was developed to include such physics and improve on the linear results. These new nonlinear-to-linear simulations confirm previous findings. This result suggests that pellet ELM triggering in ITER could require pellets much larger than those currently planned, which makes ELM-pacing operationally challenging. On the other hand, fueling pellets injected from the high-field side will likely not unintentionally trigger ELMs in an otherwise ELM-stable plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗