Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “nonlinearity”

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 235 records · Page 13

Nonlinear screening of dust grains and structurization of dusty plasma

A review of theoretical ideas on the physics of structurization instability of a homogeneous dusty plasma, i.e., the formation of zones with elevated and depressed density of dust grains and their arrangement into different structures observed in laboratory plasma under microgravity conditions, is presented. Theoretical models of compact dust structures that can form in the nonlinear stage of structurization instability, as well as models of a system of voids (both surrounding a compact structure and formed in the center of the structure), are discussed. Two types of structures with very different dimensions are possible, namely, those smaller or larger than the characteristic mean free path of ions in the plasma flow. Both of them are characterized by relatively regular distributions of dust grains; however, the first ones usually require external confinement, while the structures of the second type can be self-sustained (which is of particular interest). In this review, they are called dust clusters and self-organized dust structures, respectively. Both types of the structures are characterized by new physical processes that take place only in the presence of the dust component. The role of nonlinearities in the screening of highly charged dust grains that are often observed in modern laboratory experiments turns out to be great, but these nonlinearities have not received adequate study as of yet. Although structurization takes place upon both linear and nonlinear screening, it can be substantially different under laboratory and astrophysical conditions. Studies on the nonlinear screening of large charges in plasma began several decades ago; however, up to now, this effect was usually disregarded when interpreting the processes occurring in laboratory dusty plasma. One of the aims of the present review was to demonstrate the possibility of describing the nonlinear screening of individual grains and take it into account with the help of the basic equations for the equilibrium between plasma components when analyzing equilibrium structures. The effect of plasma screening nonlinearity on both the diffusion processes and the forces of dust drag by plasma fluxes is analyzed. It is shown how self-organized dust structures form in these processes. In the limit of very small dust grain charges, the forces acting on the dusty plasma components and the set of basic equations for stationary dust structures (with allowance for nonlinear screening) take a standard form. New qualitative effects, such as the suppression of diffusion due to ion scattering from dust grains and the formation of structures of different configurations, are described. A detailed comparison with previous results is performed. It is shown that the solution of basic nonlinear equations for dust structures yields new qualitative effects. A number of new effects to be studied in future dusty plasma experiments with the formation of structures in spherical chambers are predicted (it is assumed that diffusion will play a significant role under microgravity conditions). Recent ground-based experiments, as well as experiments carried out onboard the International Space Station, directly confirm the nonlinear character of screening and the significant role played by this nonlinearity in the structurization of dusty plasma. Experiments on the formation of structures consisting of smaller dust grains within structures formed of larger grains are discussed. It is shown that those experiments can be interpreted only using the concept of nonlinear screening.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Flat nonlinear optics with intersubband polaritonic metasurfaces

Nonlinear intersubband polaritonic metasurfaces produce some of the strongest second- and third-order nonlinear optical responses reported for condensed matter systems at infrared frequencies. These metasurfaces are fabricated as two-dimensional arrays of nanoresonators from multi-quantum-well semiconductor heterostructures, designed to produce strong nonlinear responses associated with intersubband transitions. By optimally coupling the optical modes of the nanoresonators to vertically polarized intersubband transitions in semiconductor heterostructures, one can boost the nonlinear response associated with intersubband transitions, make intersubband transitions interact with free-space radiation at normal incidence, and hence produce optically thin flat nonlinear optical elements compatible with free-space optical setups. As a result of the strong nonlinear response in these systems, significant nonlinear conversion efficiencies (>0.1 %) can be attained in deeply subwavelength optical films using modest pumping intensities of only 10–100 kW/cm 2 . Subwavelength metasurface thickness relaxes phase-matching constraints limiting the operation of bulk nonlinear crystals. Furthermore, the amplitude and phase of the nonlinear optical response in intersubband polaritonic metasurfaces can be tailored for a specific pump wavelength and a nonlinear process of interest through the co-optimization of quantum engineering of electron states in semiconductor heterostructures and photonic engineering of the metasurface nanoresonators design. Additionally, an applied voltage can dynamically control the amplitude and phase of the nonlinear optical response at a nanoresonator level. Here, we review the current state of the art in this rapidly expanding field, focusing on nonlinear processes supporting second-harmonic generation, saturable absorption, and optical power limiting.

flat optics↗

Differences between direct current and alternating current capacitance nonlinearities in high-k dielectrics and their relation to hopping conduction

Capacitance nonlinearities were studied in atomic layer deposited HfO{sub 2} films using two types of signals: a pure ac voltage of large magnitude (ac nonlinearities) and a small ac voltage superimposed to a large dc voltage (dc nonlinearities). In theory, ac and dc nonlinearities should be of the same order of magnitude. However, in practice, ac nonlinearities are found to be an order of magnitude higher than dc nonlinearities. Besides capacitance nonlinearities, hopping conduction is studied using low-frequency impedance measurements and is discussed through the correlated barrier hopping model. The link between hopping and nonlinearity is established. The ac nonlinearities are ascribed to the polarization of isolated defect pairs, while dc nonlinearities are attributed to electrode polarization which originates from defect percolation paths. Both the ac and dc capacitance nonlinearities display an exponential variation with voltage, which results from field-induced lowering of the hopping barrier energy.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Nonlinear physics of shear Alfvén waves

Shear Alfvén waves (SAW) play fundamental roles in thermonuclear plasmas of fusion interest, since they are readily excited by energetic particles in the MeV range as well as by the thermal plasma components. Thus, understanding fluctuation induced transport in burning plasmas requires understanding nonlinear SAW physics. There exist two possible routes to nonlinear SAW physics: (i) wave-wave interactions and the resultant spectral energy transfer; (ii) nonlinear wave-particle interactions of SAW instabilities with energetic particles. Within the first route, it is advantageous to understand and describe nonlinear processes in term of proximity of the system to the Alfvénic state, where wave-wave interactions are minimized due to the cancellation of Reynolds and Maxwell stresses. Here, various wave-wave nonlinear dynamics are elucidated in terms of how they break the Alfvénic state. In particular, we discuss the qualitative and quantitative modification of the SAW parametric decay process due to finite ion compressibility and finite ion Larmor radius. We also show that toroidal geometry plays a crucial role in the nonlinear excitation of zonal structures by Alfvén eigenmodes. Within the second route, the coherent nonlinear dynamics of structures in the energetic particle phase space, by which secular resonant particle transport can occur on meso- and macro-scales, must be addressed and understood. These 'nonlinear equilibria' or 'phase-space zonal structures' dynamically evolve on characteristic (fluctuation induced) turbulent transport time scales, which are generally of the same order of the nonlinear time scale of the underlying fluctuations. In this work, we introduce the general structure of nonlinear Schrödinger equations with complex integro-differential nonlinear terms, which govern these physical processes. To elucidate all these aspects, theoretical analyses are presented together with numerical simulation results.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Nonlinear X-ray Optics

Novel X-ray free-electron lasers (XFELs) are capable of producing radiation with unprecedented properties. Specifically, their ultrashort coherent X-ray pulses can generate extreme peak intensities where ordinary rules of light-matter interaction may no longer apply, and nonlinear processes start to become important. XFELs have made it possible to observe some of the most important coherent non-sequential nonlinearities, including X-ray-optical sum frequency generation (SFG), X-ray second harmonic generation (XSHG), two-photon absorption (TPA) and nonlinear two-photon Compton scattering (2PCS). In particular, the 2PCS experiment has led to unexpected results, namely the observation of a substantial anomalous red shift in the energy of photons generated by nonlinear two-photon X-ray Compton scattering in beryllium. The energy shift is in addition to the predicted nonlinear Compton shift and was neither expected from extrapolations from linear X-ray interactions nor from nonlinear effects at optical wavelengths. Our results suggest a novel nonlinear scattering mechanism that can only be observed at high fields at X-ray wavelengths, where the photon energy is a significant fraction of the electron rest mass (and thus a significant momentum transfer to the electron occurs during the scattering) and where the interaction is not dominated by the quiver motion of the electron in the field (ponderomotive potential) as is usually the case at longer (optical) wavelengths. The goal of this research project was the investigation of fundamental nonlinear X-ray-matter interactions. As these processes are typically weak, part of this research is the exploration of more efficient nonlinear X-ray processes by studying novel mechanisms including solid-state effects, the assistance of optical radiation and novel phase-matching techniques. These experiments require extremely large X-ray field strengths, such as produced by XFELs. Nonlinear X-ray effects can lead to important applications in many research fields. For example, it could lead to instantaneous plasma diagnostics for materials in extreme conditions, as a method for combining atomic-scale structural sensitivity with chemical specificity, or as a probe of electron dynamics in solids on the attosecond (1 as = 10 -18 s) time- and Ångstrom (10 -10 m) length scales. In addition, a detailed understanding of nonlinear X-ray matter interactions is crucial for understanding and extending the limits of structural determination at high intensity such as in serial femto-crystallography. The results could have a profound impact on advanced XFEL light sources such as the LCLS-II.

47 OTHER INSTRUMENTATION↗

Nonlinear interactions in mixing layers and compressible heated round jets

The nonlinear interactions between a fundamental instability mode and both its harmonics and the changing mean flow are studied using the weakly nonlinear stability theory of Stuart and Watson, and numerical solutions of coupled nonlinear partial differential equations. The first part focuses on incompressible cold (or isothermal; constant temperature throughout) mixing layers, and for these, the first and second Landau constants are calculated as functions of wavenumber and Reynolds number. It is found that the dominant contribution to the Landau constants arises from the mean flow changes and not from the higher harmonics. In order to establish the range of validity of the weakly nonlinear theory, the weakly nonlinear and numerical solutions are compared and the limitation of each is discussed. At small amplitudes and at low-to-moderate Reynolds numbers, the two results compare well in describing the saturation of the fundamental, the distortion of the mean flow, and the initial stages of vorticity roll-up. At larger amplitudes, the interaction between the fundamental, second harmonic, and the mean flow is strongly nonlinear and the numerical solution predicts flow oscillations, whereas the weakly nonlinear theory yields saturation. In the second part, the weakly nonlinear theory is extended to heated (or nonisothermal; mean temperature distribution) subsonic round jets where quadratic and cubic nonlinear interactions are present, and the Landau constants also depend on jet temperature ratio, Mach number and azimuthal mode number. Under exponential growth and nonlinear saturation, it is found that heating and compressibility suppress the growth of instability waves, that the first azimuthal mode is the dominant instability mode, and that the weakly nonlinear solution describes the early stages of the roll-up of an axisymmetric shear layer. The receptivity of a typical jet flow to pulse type input disturbance is also studied by solving the initial value problem and then examining the behavior of the long-time solution.

Jarrah, Yousef Mohd↗

Modulational instability in oppositely doubled-doped directed couplers with a non-kerr-like nonlinear refractive-index change

We theoretically investigate modulation instability (MI) in a two-core nonlinear oppositely directed coupler with a negative-index material channel in the presence of a non-kerr-like nonlinear refractive-index change. Using linear stability analysis, we obtain an expression for the instability gain. It is found that the MI in the nonlinear two-core oppositely directed coupler is significantly influenced by the ratio of the forward- to backward-propagating wave power, nonlinearity, power and saturation. In the normal dispersion regime, the septic nonlinearity is found to suppress the MI by reducing both the gain and width of the instability nonlinearity. But in anomalous dispersion regime, the situation is different because septic nonlinearity enhances the MI by increasing both the gain and width of MI. The increase of septic nonlinearity both in normal and anomalous dispersion regimes increases the gain of the instability nonlinearity. Thus, we report new ways to generate and manipulate the MI and solitons in two-core oppositely directed couplers with a particular emphasis on a negative-index material (NIM) channel in presence of septic nonlinearity.

36 MATERIALS SCIENCE↗

Nonlinear Wavefront Control with All-Dielectric Metasurfaces

Metasurfaces, two-dimensional lattices of nanoscale resonators, offer unique opportunities for functional flat optics and allow the control of the transmission, reflection, and polarization of a wavefront of light. Recently, all-dielectric metasurfaces reached remarkable efficiencies, often matching or out-performing conventional optical elements. The exploitation of the nonlinear optical response of metasurfaces offers a paradigm shift in nonlinear optics, and dielectric nonlinear metasurfaces are expected to enrich subwavelength photonics by enhancing substantially nonlinear response of natural materials combined with the efficient control of the phase of nonlinear waves. Here, we suggest a novel and rather general approach for engineering the wavefront of parametric waves of arbitrary complexity generated by a nonlinear metasurface. We design all-dielectric nonlinear metasurfaces, achieve a highly efficient wavefront control of a third-harmonic field, and demonstrate the generation of nonlinear beams at a designed angle and the generation of nonlinear focusing vortex beams. Lastly, our nonlinear metasurfaces produce phase gradients over a full 0–2π phase range with a 92% diffraction efficiency.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Superposition of elliptic functions as solutions for a large number of nonlinear equations

For a large number of nonlinear equations, both discrete and continuum, we demonstrate a kind of linear superposition. We show that whenever a nonlinear equation admits solutions in terms of both Jacobi elliptic functions cn(x, m) and dn(x, m) with modulus m, then it also admits solutions in terms of their sum as well as difference. We have checked this in the case of several nonlinear equations such as the nonlinear Schrödinger equation, MKdV, a mixed KdV-MKdV system, a mixed quadratic-cubic nonlinear Schrödinger equation, the Ablowitz-Ladik equation, the saturable nonlinear Schrödinger equation, λϕ{sup 4}, the discrete MKdV as well as for several coupled field equations. Further, for a large number of nonlinear equations, we show that whenever a nonlinear equation admits a periodic solution in terms of dn{sup 2}(x, m), it also admits solutions in terms of dn {sup 2}(x,m)±√(m) cn (x,m) dn (x,m), even though cn(x, m)dn(x, m) is not a solution of these nonlinear equations. Finally, we also obtain superposed solutions of various forms for several coupled nonlinear equations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

SNS: A Solution-Based Nonlinear Subspace Method for Time-Dependent Model Order Reduction

Several reduced order models have been successfully developed for nonlinear dynamical systems. To achieve a considerable speed-up, a hyper-reduction step is needed to reduce the computational complexity due to nonlinear terms. Many hyper-reduction techniques require the construction of nonlinear term basis, which introduces a computationally expensive offline phase. A novel way of constructing nonlinear term basis within the hyper-reduction process is introduced. In contrast to the traditional hyper-reduction techniques where the collection of nonlinear term snapshots is required, the SNS method avoids collecting the nonlinear term snapshots. Instead, it uses the solution snapshots that are used for building a solution basis, which enables avoiding an extra data compression of nonlinear term snapshots. As a result, the SNS method provides a more efficient offline strategy than the traditional model order reduction techniques, such as the DEIM, GNAT, and ST-GNAT methods. The SNS method is theoretically justified by the conforming subspace condition and the subspace inclusion relation. It is useful for model order reduction of large-scale nonlinear dynamical problems to reduce the offline cost. It is especially useful for ST-GNAT that has shown promising results, such as a good accuracy with a considerable online speed-up for hyperbolic problems in a recent paper by Choi and Carlberg [SIAM J. Sci. Comput., 41 (2019), pp. A26--A58], because ST-GNAT involves an expensive offline cost related to collecting nonlinear term snapshots. Error analysis for the SNS method is presented. Numerical results support that the accuracy of the solution from the SNS method is comparable to the traditional methods and a considerable speed-up (i.e., a factor of two to a hundred) is achieved in the offline phase.

97 MATHEMATICS AND COMPUTING↗

Nonlinear bleaching, absorption, and scattering of 532-nm-irradiated plasmonic nanoparticles

Single-pulse irradiation of Au and Ag suspensions of nanospheres and nanodisks with 532-nm 4-ns pulses has identified complex optical nonlinearities while minimizing material damage. For all materials tested, we observe competition between saturable absorption (SA) and reverse SA (RSA), with RSA behavior dominating for intensities above {approx}50 MW/cm{sup 2}. Due to reduced laser damage in single-pulse experiments, the observed intrinsic nonlinear absorption coefficients are the highest reported to date for Au nanoparticles. We find size dependence to the nonlinear absorption enhancement for Au nanoparticles, peaking in magnitude for 80-nm nanospheres and falling off at larger sizes. The nonlinear absorption coefficients for Au and Ag spheres are comparable in magnitude. On the other hand, the nonlinear absorption for Ag disks, when corrected for volume fraction, is several times higher. These trends in nonlinear absorption are correlated to local electric field enhancement through quasi-static mean-field theory. Through variable size aperture measurements, we also separate nonlinear scattering from nonlinear absorption. For all materials tested, we find that nonlinear scattering is highly directional and that its magnitude is comparable to that of nonlinear absorption. These results indicate methods to improve the efficacy of plasmonic nanoparticles as optical limiters in pulsed laser systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Linear and nonlinear dynamic analysis of redundant load path bearingless rotor systems

The goal of this research is to develop the transfer matrix method to treat nonlinear autonomous boundary value problems with multiple branches. The application is the complete nonlinear aeroelastic analysis of multiple-branched rotor blades. Once the development is complete, it can be incorporated into the existing transfer matrix analyses. There are several difficulties to be overcome in reaching this objective. The conventional transfer matrix method is limited in that it is applicable only to linear branch chain-like structures, but consideration of multiple branch modeling is important for bearingless rotors. Also, hingeless and bearingless rotor blade dynamic characteristics (particularly their aeroelasticity problems) are inherently nonlinear. The nonlinear equations of motion and the multiple-branched boundary value problem are treated together using a direct transfer matrix method. First, the formulation is applied to a nonlinear single-branch blade to validate the nonlinear portion of the formulation. The nonlinear system of equations is iteratively solved using a form of Newton-Raphson iteration scheme developed for differential equations of continuous systems. The formulation is then applied to determine the nonlinear steady state trim and aeroelastic stability of a rotor blade in hover with two branches at the root. A comprehensive computer program is developed and is used to obtain numerical results for the (1) free vibration, (2) nonlinearly deformed steady state, (3) free vibration about the nonlinearly deformed steady state, and (4) aeroelastic stability tasks. The numerical results obtained by the present method agree with results from other methods.

Murthy, V. R.↗

Nonlinear Large Deflection Theory with Modified Aeroelastic Lifting Line Aerodynamics for a High Aspect Ratio Flexible Wing

This paper investigates the effect of nonlinear large deflection bending on the aerodynamic performance of a high aspect ratio flexible wing. A set of nonlinear static aeroelastic equations are derived for the large bending deflection of a high aspect ratio wing structure. An analysis is conducted to compare the nonlinear bending theory with the linear bending theory. The results show that the nonlinear bending theory is length-preserving whereas the linear bending theory causes a non-physical effect of lengthening the wing structure under the no axial load condition. A modified lifting line theory is developed to compute the lift and drag coefficients of a wing structure undergoing a large bending deflection. The lift and drag coefficients are more accurately estimated by the nonlinear bending theory due to its length-preserving property. The nonlinear bending theory yields lower lift and span efficiency than the linear bending theory. A coupled aerodynamic-nonlinear finite element model is developed to implement the nonlinear bending theory for a Common Research Model (CRM) flexible wing wind tunnel model to be tested in the University of Washington Aeronautical Laboratory (UWAL). The structural stiffness of the model is designed to give about 10% wing tip deflection which is large enough that could cause the nonlinear deflection effect to become significant. The computational results show that the nonlinear bending theory yields slightly less lift than the linear bending theory for this wind tunnel model. As a result, the linear bending theory is deemed adequate for the CRM wind tunnel model.

Aeroelasticity↗

Experimental Investigation of Elastodynamic Nonlinear Response of Dry Intact, Fractured and Saturated Rock

Nonlinear elastodynamic response of fractured rocks carries crucial information on fracture features that can be exploited to forecast flow properties, friction constitutive behavior and poromechanical response. Well-controlled laboratory experiments are designed to measure the nonlinear elastodynamic response of Westerly granite in three states: dry intact, dry fractured and saturated fractured. We study the effect of fracturing and saturation in modifying the elastodynamic response of the rock. Each sample is tested at a normal stress level of 15 MPa. We measure the elastodynamic response of an intact L-shaped sample of Westerly Granite subjected to normal stress oscillations of prescribed amplitudes (0.2–1 MPa) and frequencies (0.1, 1, 10 Hz). Ultrasonic waves transmitted across the sample are used to monitor the evolution of wave velocity before, during and after dynamic stressing. Furthermore, the nonlinearity of the elastodynamic response is measured in terms of: (1) the offset in normalized wave velocity; (2) the amplitude of wave velocity fluctuation during the oscillations; and (3) recovery rate of the wave velocity post-oscillation. We observe that the three nonlinearity parameters show a similar trend. Irrespective of the parameter, the nonlinearity measures higher for sample in dry-intact condition than that for dry-fractured and the saturated-fractured sample exhibits smaller nonlinearity than the dry-fractured sample. As expected, the saturated sample exhibits less nonlinearity than the dry intact and fractured samples due to the presence of interstitial fluid and the resulting increased interface stiffness. Conversely, the dry intact rock shows a higher nonlinearity than the dry fractured. We use numerical simulations to show that the presence of fracture significantly alters the strain distribution across the bulk of the sample and only the contacting asperities are highly strained, thus resulting in a decrease in the measured elastodynamic nonlinearity.

58 GEOSCIENCES↗

Characterization of nonlinear ultrasonic waves behavior while interacting with poor interlayer bonds in large-scale additive manufactured materials

Over the past decades, researchers have developed several nonlinear ultrasonic techniques for quality control of materials commonly used in different applications. Owing to the superior sensitivity of nonlinear ultrasound waves to small defects such as micro-cracks, their applicability in different nondestructive testing (NDT) problems has been investigated in numerous studies. These studies utilize frequency domain analysis to detect the generation of higher harmonics because of the formation of defects in the inspected medium. Frequency domain analysis based on the Fourier transform is a significant approach used in linear systems; however, it may not perform adequately on nonlinear systems. Hence, studies on nonlinear dynamics and physics consider analyzing systems' behavior in the phase-space domain. In contrast to the frequency domain analysis, which can result in information loss, analysis in the phase-space domain retains all the information regarding a system's states. Here, we investigate the nonlinearities induced by poor interlayer bonding in polymer-based additive manufactured parts in the phase-space domain. It is convenient to characterize the nonlinearity in the phase-space domain because it provides a geometrical representation of a system's states. Two types of low quality interlayer bond are considered. The first type is simulated artificially while the second type is manufactured by reducing the bond quality during the printing process. The analysis verified that the received ultrasonic signals exhibit classical nonlinear behavior in the phase-space domain while interacting with simulated poor interlayer bonds. In addition, the results showed that the behavior of ultrasonic waves is amplitude-dependent and evolves into models that have not been previously reported. Furthermore, Largest Lyapunov Exponent (LLE) is used to quantify the behavior of nonlinear ultrasonic waves while interacting with poor interlayer bonds. Using LLE, it was observed that the divergence rate of the phase-space trajectories depends on the amplitude of the excitation. This observation quantitatively proves that nonlinear behavior of ultrasound while interacting with poor interlayer bonds can be amplitude-dependent. The results of both simulated and inherent poor interlayer bond cases showed that LLE can be used as a reliable quantitative damage-sensitive feature to detect and potentially characterize weak bonds, which are difficult to detect using conventional approaches. Additionally, the reported results in the phase-space domain provide a basis for proposing a new mathematical model for ultrasonic waves interacting with poor interlayer bonds.

36 MATERIALS SCIENCE↗

Relationships between nonlinear normal modes and response to random inputs

The ability to model nonlinear structures subject to random excitation is of key importance in designing hypersonic aircraft and other advanced aerospace vehicles. When a structure is linear, superposition can be used to construct its response to a known spectrum in terms of its linear modes. Superposition does not hold for a nonlinear system, but several works have shown that a system's dynamics can still be understood qualitatively in terms of its nonlinear normal modes (NNMs). Here, this work investigates the connection between a structure's undamped nonlinear normal modes and the spectrum of its response to high amplitude random forcing. Two examples are investigated: a spring-mass system and a clamped-clamped beam modeled within a geometrically nonlinear finite element package. In both cases, an intimate connection is observed between the smeared peaks in the response spectrum and the frequency-energy dependence of the nonlinear normal modes. In order to understand the role of coupling between the underlying linear modes, reduced order models with and without modal coupling terms are used to separate the effect of each NNM's backbone from the nonlinear couplings that give rise to internal resonances. In the cases shown here, uncoupled, single-degree-of-freedom nonlinear models are found to predict major features in the response with reasonable accuracy; a highly inexpensive approximation such as this could be useful in design and optimization studies. More importantly, the results show that a reduced order model can be expected to give accurate results only if it is also capable of accurately predicting the frequency-energy dependence of the nonlinear modes that are excited.

42 ENGINEERING↗

Nonlinear Elimination Applied to Radiation Diffusion

We apply a nonlinearly preconditioned, quasi-Newton framework to accelerate the numerical solution of the thermal radiative transfer (TRT) equations. This framework was inspired by the unpublished method that has existed for years in Teton, Lawrence Livermore National Laboratory’s deterministic TRT code. In this paper, we cast this iteration scheme within a formal nonlinear preconditioning framework and compare its performance against other iteration schemes in the framework. With proper choices of iteration controls for the various levels of the solver, we can recover the standard linearized one-step method, a full nonlinear Newton scheme, as well as the method in Teton. In brief, the nonlinear preconditioning TRT scheme formally eliminates the material temperature equation from the nonlinear system in a nonlinear analog of a Schur complement. This nonlinear elimination step involves solving a decoupled nonlinear equation for each spatial degree of freedom and is therefore inexpensive. By applying a quasi-Newton iteration scheme on the new system, we obtain a three-level iteration scheme that is at least as efficient as commonly used TRT schemes. The new method allows full convergence to the nonlinear backward Euler time-discretized system, increasing accuracy and robustness, while using a similar number of linear iterations as the more common linearized one-step methods Eq. (4).

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Calibrating Nonlinear Soil Material Properties for Seismic Analysis Using Soil Material Properties Intended for Linear Analysis

Seismic analysis of nuclear structures is routinely performed using guidance provided in “Seismic Analysis of Safety-Related Nuclear Structures and Commentary (ASCE 4, 1998).” This document, which is currently under revision, provides detailed guidance on linear seismic soil-structure-interaction (SSI) analysis of nuclear structures. To accommodate the linear analysis, soil material properties are typically developed as shear modulus and damping ratio versus cyclic shear strain amplitude. A new Appendix in ASCE 4-2014 (draft) is being added to provide guidance for nonlinear time domain SSI analysis. To accommodate the nonlinear analysis, a more appropriate form of the soil material properties includes shear stress and energy absorbed per cycle versus shear strain. Ideally, nonlinear soil model material properties would be established with soil testing appropriate for the nonlinear constitutive model being used. However, much of the soil testing done for SSI analysis is performed for use with linear analysis techniques. Consequently, a method is described in this paper that uses soil test data intended for linear analysis to develop nonlinear soil material properties. To produce nonlinear material properties that are equivalent to the linear material properties, the linear and nonlinear model hysteresis loops are considered. For equivalent material properties, the shear stress at peak shear strain and energy absorbed per cycle should match when comparing the linear and nonlinear model hysteresis loops. Consequently, nonlinear material properties are selected based on these criteria.

nonlinear↗