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,531 records · Page 85
All sourcesNASA NTRSDOE OSTI
DOE OSTI2013

Instability of evaporation fronts in the interstellar medium

The neutral component of the interstellar medium is segregated into the cold neutral medium (CNM) and warm neutral medium (WNM) as a result of thermal instability. It was found that a plane-parallel CNM-WNM evaporation interface, across which the CNM undergoes thermal expansion, is linearly unstable to corrugational disturbances, in complete analogy with the Darrieus-Landau instability (DLI) of terrestrial flames. We perform a full linear stability analysis as well as nonlinear hydrodynamic simulations of the DLI of such evaporation fronts in the presence of thermal conduction. We find that the DLI is suppressed at short length scales by conduction. The length and time scales of the fastest growing mode are inversely proportional to the evaporation flow speed of the CNM and its square, respectively. In the nonlinear stage, the DLI saturates to a steady state where the front deforms to a finger-like shape protruding toward the WNM, without generating turbulence. The evaporation rate at nonlinear saturation is larger than the initial plane-parallel value by a factor of ∼2.4 when the equilibrium thermal pressure is 1800 k {sub B} cm{sup –3} K. The degrees of front deformation and evaporation-rate enhancement at nonlinear saturation are determined primarily by the density ratio between the CNM and WNM. We demonstrate that the Field length in the thermally unstable medium should be resolved by at least four grid points to obtain reliable numerical outcomes involving thermal instability.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Fourier imaging of non-linear structure formation

We perform a Fourier space decomposition of the dynamics of non-linear cosmological structure formation in ΛCDM models. From N -body simulations involving only cold dark matter we calculate 3-dimensional non-linear density, velocity divergence and vorticity Fourier realizations, and use these to calculate the fully non-linear mode coupling integrals in the corresponding fluid equations. Our approach allows for a reconstruction of the amount of mode coupling between any two wavenumbers as a function of redshift. With our Fourier decomposition method we identify the transfer of power from larger to smaller scales, the stable clustering regime, the scale where vorticity becomes important, and the suppression of the non-linear divergence power spectrum as compared to linear theory. Our results can be used to improve and calibrate semi-analytical structure formation models.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2011

SDO/AIA OBSERVATION OF KELVIN-HELMHOLTZ INSTABILITY IN THE SOLAR CORONA

We present observations of the formation, propagation, and decay of vortex-shaped features in coronal images from the Solar Dynamics Observatory associated with an eruption starting at about 2:30 UT on 2010 April 8. The series of vortices were formed along the interface between an erupting (dimming) region and the surrounding corona. They ranged in size from several to 10 arcsec and traveled along the interface at 6-14 km s{sup -1}. The features were clearly visible in six out of the seven different EUV wave bands of the Atmospheric Imaging Assembly. Based on the structure, formation, propagation, and decay of these features, we identified the event as the first observation of the Kelvin-Helmholtz instability (KHI) in the corona in EUV. The interpretation is supported by linear analysis and by a nonlinear 2.5-dimensional magnetohydrodynamic model of KHI. We conclude that the instability is driven by the velocity shear between the erupting and closed magnetic field of the coronal mass ejection. The shear-flow-driven instability can play an important role in energy transfer processes in coronal plasma.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Gravitational memory and compact extra dimensions

Here we develop a general formalism for treating radiative degrees of freedom near $\mathscr{I}^+$ in theories with an arbitrary Ricci-flat internal space. These radiative modes are encoded in a generalized news tensor which decomposes into gravitational, electromagnetic, and scalar components. We find a preferred gauge which simplifies the asymptotic analysis of the full nonlinear Einstein equations and makes the asymptotic symmetry group transparent. This asymptotic symmetry group extends the Bondi–Metzner–Sachs (BMS) group to include angle-dependent isometries of the internal space. We apply this formalism to study memory effects, which are expected to be observed in future experiments, that arise from bursts of higher-dimensional gravitational radiation. We outline how measurements made by gravitational wave observatories might probe properties of the compact extra dimensions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2010

Characterization of etch pits found on a large-grain bulk niobium superconducting radio-frequency resonant cavity

The performance of superconducting radio-frequency (SRF) resonant cavities made of bulk niobium is limited by nonlinear localized effects. Surface analysis of regions of higher power dissipation is thus of intense interest. Such areas (referred to as “hotspots”) were identified in a large-grain single-cell cavity that had been buffered-chemical polished and dissected for examination by high resolution electron microscopy, electron backscattered diffraction microscopy (EBSD), and optical microscopy. Pits with clearly discernible crystal facets were observed in both “hotspot” and “coldspot” specimens. The pits were found in-grain, at bicrystal boundaries, and on tricrystal junctions. They are interpreted as etch pits induced by crystal defects (e.g. dislocations). All coldspots examined had a qualitatively lower density of etch pits or relatively smooth tricrystal boundary junctions. EBSD mapping revealed the crystal orientation surrounding the pits. Locations with high pit density are correlated with higher mean values of the local average misorientation angle distributions, indicating a higher geometrically necessary dislocation content. In addition, a survey of the samples by energy dispersive x-ray analysis did not show any significant contamination of the samples’ surface. In conclusion, the local magnetic field enhancement produced by the sharp-edge features observed on the samples is not sufficient to explain the observed degradation of the cavity quality factor, which starts at peak surface magnetic field as low as 20 mT.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Nonlinear normal modes in electrodynamic systems: A nonperturbative approach

We consider electromagnetic nonlinear normal modes in cylindrical cavity resonators filled with a nonlinear nondispersive medium. The key feature of the analysis is that exact analytic solutions of the nonlinear field equations are employed to study the mode properties in detail. Based on such a nonperturbative approach, we rigorously prove that the total energy of free nonlinear oscillations in a distributed conservative system, such as that considered in our work, can exactly coincide with the sum of energies of the normal modes of the system. This fact implies that the energy orthogonality property, which has so far been known to hold only for linear oscillations and fields, can also be observed in a nonlinear oscillatory system.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Heavy Ion-Acoustic Solitary Waves and Double Layers in a Multi-Ion Plasma

The formation and propagation of small-amplitude heavy-ion-acoustic (HIA) solitary waves and double layers in an unmagnetized collisionless multicomponent plasma system consisting of superthermal electrons, Boltzmann distributed light ions, and adiabatic positively charged inertial heavy ions are theoretically investigated. The reductive perturbation technique is employed to derive the modified Korteweg–de Vries (mKdV) and standard Gardner (SG) equations. The solitary wave (SW) solution of mKdV and SG equations, as well as double layers (DLs) solution of SG equation, is studied for analysis of higher order nonlinearity. It is found that the plasma system under consideration supports positive and negative potential Gardner solitons, but only positive potential mKdV solitons. In addition, it is shown that, the basic properties of HIA mKdV and Gardner solitons and DLs (viz. polarity, amplitude, width, and phase speed) are incomparably influenced by the adiabaticity effect of heavy ions and the superthermality effect of electrons. The relevance of the present findings to the system of space plasmas, as well as to the system of researchers interest, is specified.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Estimating epistemic uncertainty in soil parameters for nonlinear site response analyses: Introducing the Latin Hypercube Sampling technique

This study quantifies the effects of epistemic uncertainty in soil parameters on nonlinear (NL) site response analysis (SRA) results, validated against the data recorded at a well-instrumented geotechnical downhole array located in Japan. To this end, a one-dimensional soil column model of the Service Hall Array (SHA) near the Kashiwazaki-Kariwa Nuclear Power Plant (KKNPP) is developed using the finite element (FE) program LS-DYNA. The dynamic stress–strain relationship is characterized by a modified two-stage hyperbolic (MTH) NL backbone curve formulation capable of capturing soil behavior at both small- and large-shear strains. The model is then validated against the ground motion recordings to capture the model bias. The uncertainties associated with the shear-wave velocity profile (a small-strain soil property) and soil shear strength (a large-strain soil property) are incorporated in NL SRA to quantify their separate and joint randomization effects on the results. This study proposes using the Latin Hypercube Sampling (LHS) method as an efficient alternative to commonly used methods, such as Standard Monte Carlo (SMC), to account for uncertainty propagation in such reliability analysis. Both low-intensity and design-level records from the recordings at the SHA are applied to study the contribution of the small- and large-strain NL dynamic soil properties. Results from 46,200 NL FE analyses (23,100 per input motion) are presented. Measured and predicted site response, using recorded ground motions at this downhole array, is compared to assess the significance of soil parameter uncertainty on the observed ground motion dispersions. It is demonstrated that increasing the number of soft realizations and implementing higher level earthquake intensity lead to higher ground motion dispersion. Unlike past studies in randomization of Vs profiles with the SMC method, the LHS method is shown to have no significant effect on the predicted median surface response spectra and amplification factors (AFs) for this case study.

Engineering↗
DOE OSTI2012

Center for Gyrokinetic/MHD Hybrid Simulation of Energetic Particle Physics in Toroidal Plasmas (Final Report)

At Colorado University-Boulder the primary task is to extend our gyrokinetic Particle-in-Cell simulation of tokamak micro-turbulence and transport to the area of energetic particle physics. We have implemented a gyrokinetic ion/massless fluid electron hybrid model in the global δƒ-PIC code GEM, and benchmarked the code with analytic results on the thermal ion radiative damping rate of Toroidal Alfvén Eigenmodes (TAE) and with mode frequency and spatial structure from eigenmode analysis. We also performed nonlinear simulations of both a single-n mode (n is the toroidal mode number) and multiple-n modes, and in the case of single-n, benchmarked the code on the saturation amplitude vs. particle collision rate with analytical theory. Most simulations use the f method for both ions species, but we have explored the full-ƒ method for energetic particles in cases where the burst amplitude of the excited instabilities is large as to cause significant re-distribution or loss of the energetic particles. We used the hybrid model to study the stability of high-n TAEs in ITER. Our simulations show that the most unstable modes in ITER lie in the rage of 10 < n < 20. Thermal ion pressure effect and alpha particles non-perturbative effect are important in determining the mode radial location and stability threshold. The thermal ion Landau damping rate and radiative damping rate from the simulations are compared with analytical estimates. The thermal ion Landau damping is the dominant damping mechanism. Plasma elongation has a strong stabilizing effect on the alpha driven TAEs. The central alpha particle pressure threshold for the most unstable n=15 mode is about β α (0) = 0.7% for the fully shaped ITER equilibrium. We also carried nonlinear simulations of the most unstable n = 15 mode and found that the saturation amplitude for the nominal ITER discharge is too low to cause large redistribution or loss of alpha particles. To include kinetic electron effects in the hybrid model we have studied a kinetic electron closure scheme for the fluid electron model. The most important element of the closure scheme is a complete Ohm’s law for the parallel electric field E || , derived by combining the quasi-neutrality condition, the Ampere’s equation and the v || moment of the gyrokinetic equations. A discretization method for the closure scheme is studied in detail for a three-dimensional shear-less slab plasma. It is found that for long-wavelength shear Alfvén waves the kinetic closure scheme is both more accurate and robust than the previous GEM algorithm using the split-scheme, whereas for the ion-gradient-driven instability the previous algorithm is more efficient. This kinetic electron closure scheme will be implemented in GEM in the future. We have studied the beam driven Reverse Shear Alfven Eigenmodes (RSAE) observed in DIII-D discharge 142111. For this purpose a new scheme for obtaining the electric potential is implemented, i.e., by solving the gyrokinetic moment (GKM) equation, which is essentially the equation for $∂\phi/∂t$ used in GEM’s split-weight scheme, and then integrating in time. Due to charge-neutrality the ExB motions of the equilibrium densities of all species cancel each other and do not cause charge separation if there is no finite Larmor radius effect. The advantage of solving the GKM equation is that this lowest-order cancellation can be made explicit. The GKM approach is found to be more accurate and robust. GEM simulations have reproduced many features of RSAE seen in the experiment, such as frequency chirping and the chirping range. It has been reported by other simulation codes that the shearing direction of the mode structure in the poloidal plane disagrees with observation. We found that the mode structure, including the shearing in the poloidal plane, is in general sensitive to the beam distribution. Using the same beam density profile as in other codes but with a slowing-down distribution in velocity, GEM simulations reproduce the mode shearing direction seen in the experiment. We have carried out extensive nonlinear simulations of RSAE.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Approaches and Recommendations to Address Regulatory Gaps Identified in INL/EXT-15-36945

The information presented in this report completes Task 2.3, identified in Table 1 of 14-CR-04, Mod 3. The title of task 2.3 is, “Develop approaches to address regulatory gaps identified in INL/EXT-15-36945 and prepare recommendations.” The purpose of the task is to identify regulatory gaps related to nonlinear soil structure interaction analysis and seismic isolation of TWR-300 and propose a path to close the gaps and complete the licensing basis.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2016

INSTABILITY OF NON-UNIFORM TOROIDAL MAGNETIC FIELDS IN ACCRETION DISKS

We present a new type of instability that is expected to drive magnetohydrodynamic (MHD) turbulence from a purely toroidal magnetic field in an accretion disk. It is already known that in a differentially rotating system, the uniform toroidal magnetic field is unstable due to magnetorotational instability (MRI) under a non-axisymmetric and vertical perturbation, while it is stable under a purely vertical perturbation. Contrary to the previous study, this paper proposes an unstable mode completely confined to the equatorial plane, driven by the expansive nature of the magnetic pressure gradient force under a non-uniform toroidal field. The basic nature of this growing eigenmode, which we name “magneto-gradient driven instability,” is studied using linear analysis, and the corresponding nonlinear evolution is then investigated using two-dimensional ideal MHD simulations. Although a single localized magnetic field channel alone cannot provide sufficient Maxwell stress to contribute significantly to the angular momentum transport, we find that the mode coupling between neighboring toroidal fields under multiple localized magnetic field channels drastically generates a highly turbulent state and leads to the enhanced transport of angular momentum, which is comparable to the efficiency seen in previous studies on MRIs. This horizontally confined mode may play an important role in the saturation of an MRI through complementray growth with the toroidal MRIs and coupling with magnetic reconnection.

79 ASTRONOMY AND ASTROPHYSICS↗