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 901 records · Page 50

Detection of ratcheting in finite element calculations

The distinction between a ratcheting and non-ratcheting response is critical for many high temperature design methods. Non-ratcheting is generally considered safe deformation remain bounded over the lifetime of the component while ratcheting is undesirable. As a particular example, the elastic perfectly plastic (EPP) design methods described in recent ASME Section III, Division 5 code cases require a designer to distinguish ratcheting from non-ratcheting for finite element analyses using a relatively simple, elastic perfectly-plastic constitutive response. However, it can be quite difficult to distinguish these two deformation regimes using finite element (FE) analysis particularly in the case where the actual ratcheting strain is small. In practice FE analysis of structures that are analytically in either the plastic shakedown or ratcheting regimes will result in small, cycle to -cycle accumulated strains characteristic of ratcheting. Distinguishing false ratcheting the structure is actually in the plastic shakedown regime from true ratcheting can be challenging. We describe the characteristics of nonlinear FE analysis that cause these false ratcheting strains and describe practical methods for distinguishing a ratcheting from a non-ratcheting response.

Messner, M. C.↗

Application of Nonlinear Seismic Soil-Structure Interaction Analysis for Identification of Seismic Margins at Nuclear Power Plants

Seismic probabilistic risk assessment (SPRA) methods and approaches at nuclear power plants (NPP) were first developed in the 1970s and aspects of them have matured over time as they were applied and incrementally improved. SPRA provides information on risk and risk insights and allows for some accounting for uncertainty and variability. As a result, SPRA is now used as an important basis for risk-informed decision making for both new and operating NPPs in the US and in an increasing number of countries globally. SPRAs are intended to provide best estimates of the various combinations of structural and equipment failures that can lead to a seismic induced core damage event. However, in some instances the current SPRA approach contains large uncertainties, and potentially masks other important events (for instance, it was not the seismic motions that caused the Fukushima core melt events, but the tsunami ingress into the facility). INL has an advanced SPRA research and development (R&D) activity that will identify areas in the calculation process that contain significant uncertainties. One current area of focus is the use of nonlinear soil-structure interaction (NLSSI) analysis methods to accurately capture: 1) nonlinear soil behavior and 2) gapping and sliding between the NPP and soil. The goal of this study is to compare numerical NLSSI analysis results with recorded earthquake ground motions at Fukushima Daichii (Great Tohuku Earthquake) and evaluate the sources of nonlinearity contributing to the observed reduction in peak acceleration. Comparisons are made using recorded data in the free-field (soil column with no structural influence) and recorded data on the NPP basemat (in-structure response). Results presented in this study should identify areas of focus for future R&D activities with the goal of minimizing uncertainty in SPRA calculations. This is not a validation activity since there are too many sources of uncertainty that a numerical analysis would need to consider (variability in soil material properties, structural material properties, etc.). Rather the report will determine if the NLSSI calculations are following similar trends observed in the recorded data (i.e. reductions in maximum acceleration between the free-field and basemat) Numerical NLSSI results presented show maximum accelerations between the free field and basemat were reduced the EW and NS directions. The maximum acceleration in the UD direction increased slightly. The largest reduction in maximum accelerations between the modeled free-field and the NPP basemat resulted in nearly 50% reduction. The observation in reduction of numerical maximum accelerations in the EW and NS directions follows the observed trend in the recorded data. The maximum reductions observed in these NLSSI studies were due to soil nonlinearities, not gapping and sliding (although additional R&D is needed to develop an appropriate approach to model gapping and sliding). This exploratory study highlights the need for additional R&D on developing: (i) improved modeling of soil nonlinearities (soil constitutive models that appropriately capture cyclic soil behavior), (ii) improved modeling of gapping and sliding at the soil-structure interface (to appropriately capture the dissipation of energy at this interface), and (iii) experimental laboratory test data to calibrate the items (i) and (ii).

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Finite element analysis of steady and transiently moving/rolling nonlinear viscoelastic structure. Part 1: Theory

In a three part series of papers, a generalized finite element analysis scheme is developed to handle the steady and transient response of moving/rolling nonlinear viscoelastic structure. This paper considers the development of the moving/rolling element strategy, including the effects of large deformation kinematics and viscoelasticity modelled by fractional integro-differential operators. To improve the solution strategy, a special hierarchical constraint procedure is developed for the case of steady rolling/translating as well as a transient scheme involving the use of a Grunwaldian representation of the fractional operator. In the second and third parts of the paper, 3-D extensions are developed along with transient contact strategies enabling the handling of impacts with obstructions. Overall, the various developments are benchmarked via comprehensive 2- and 3-D simulations. These are correlated with experimental data to define modelling capabilities.

Padovan, Joe↗

Finite element analysis of steady and transiently moving/rolling nonlinear viscoelastic structure. I - Theory

In a three-part series of papers, a generalized finite element analysis scheme is developed to handle the steady and transient response of moving/rolling nonlinear viscoelastic structure. This paper considers the development of the moving/rolling element strategy, including the effects of large deformation kinematics and viscoelasticity modeled by fractional integrodifferential operators. To improve the solution strategy, a special hierarchical constraint procedure is developed for the case of steady rolling/translating, as well as a transient scheme involving the use of a Grunwaldian representation of the fractional operator.

Padovan, Joe↗

TAP 2: A finite element program for thermal analysis of convectively cooled structures

A finite element computer program (TAP 2) for steady-state and transient thermal analyses of convectively cooled structures is presented. The program has a finite element library of six elements: two conduction/convection elements to model heat transfer in a solid, two convection elements to model heat transfer in a fluid, and two integrated conduction/convection elements to represent combined heat transfer in tubular and plate/fin fluid passages. Nonlinear thermal analysis due to temperature-dependent thermal parameters is performed using the Newton-Raphson iteration method. Transient analyses are performed using an implicit Crank-Nicolson time integration scheme with consistent or lumped capacitance matrices as an option. Program output includes nodal temperatures and element heat fluxes. Pressure drops in fluid passages may be computed as an option. User instructions and sample problems are presented in appendixes.

Thornton, E. A.↗

Monitoring microstructural evolution in irradiated steel with second harmonic generation

Material damage in structural components is driven by microstructural evolution that occurs at low length scales and begins early in component life. In metals, these microstructural features are known to cause measurable changes in the acoustic nonlinearity parameter. Physically, the interaction of a monochromatic ultrasonic wave with microstructural features such as dislocations, precipitates, and vacancies, generates a second harmonic wave that is proportional to the acoustic nonlinearity parameter. These nonlinear ultrasonic techniques thus have the capability to evaluate initial material damage, particularly before crack initiation and propagation occur. This paper discusses how the nonlinear ultrasonic technique of second harmonic generation can be used as a nondestructive evaluation tool to monitor microstructural changes in steel, focusing on characterizing neutron radiation embrittlement in nuclear reactor pressure vessel steels. Current experimental evidence and analytical models linking microstructural evolution with changes in the acoustic nonlinearity parameter are summarized.

36 MATERIALS SCIENCE↗

Metamaterials-based sensor to detect and locate nonlinear elastic sources

In recent years, acoustic metamaterials have attracted increasing scientific interest for very diverse technological applications ranging from sound abatement to ultrasonic imaging, mainly due to their ability to act as band-stop filters. At the same time, the concept of chaotic cavities has been recently proposed as an efficient tool to enhance the quality of nonlinear signal analysis, particularly in the ultrasonic/acoustic case. The goal of the present paper is to merge the two concepts in order to propose a metamaterial-based device that can be used as a natural and selective linear filter for the detection of signals resulting from the propagation of elastic waves in nonlinear materials, e.g., in the presence of damage, and as a detector for the damage itself in time reversal experiments. Numerical simulations demonstrate the feasibility of the approach and the potential of the device in providing improved signal-to-noise ratios and enhanced focusing on the defect locations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Stable Tearing and Buckling Responses of Unstiffened Aluminum Shells with Long Cracks

The results of an analytical and experimental study of the nonlinear response of thin, unstiffened, aluminum cylindrical shells with a long longitudinal crack are presented. The shells are analyzed with a nonlinear shell analysis code that accurately accounts for global and local structural response phenomena. Results are presented for internal pressure and for axial compression loads. The effect of initial crack length on the initiation of stable crack growth and unstable crack growth in typical shells subjected to internal pressure loads is predicted using geometrically nonlinear elastic-plastic finite element analyses and the crack-tip-opening angle (CTOA) fracture criterion. The results of these analyses and of the experiments indicate that the pressure required to initiate stable crack growth and unstable crack growth in a shell subjected to internal pressure loads decreases as the initial crack length increases. The effects of crack length on the prebuckling, buckling and postbuckling responses of typical shells subjected to axial compression loads are also described. For this loading condition, the crack length was not allowed to increase as the load was increased. The results of the analyses and of the experiments indicate that the initial buckling load and collapse load for a shell subjected to axial compression loads decrease as the initial crack length increases. Initial buckling causes general instability or collapse of a shell for shorter initial crack lengths. Initial buckling is a stable local response mode for longer initial crack lengths. This stable local buckling response is followed by a stable postbuckling response, which is followed by general or overall instability of the shell.

Starnes, James H., Jr.↗

Stable Tearing and Buckling Responses of Unstiffened Aluminum Shells with Long Cracks

The results of an analytical and experimental study of the nonlinear response of thin, unstiffened, aluminum cylindrical shells with a long longitudinal crack are presented. The shells are analyzed with a nonlinear shell analysis code that accurately accounts for global and local structural response phenomena. Results are presented for internal pressure and for axial compression loads. The effect of initial crack length on the initiation of stable crack growth and unstable crack growth in typical shells subjected to internal pressure loads is predicted using geometrically nonlinear elastic-plastic finite element analyses and the crack-tip-opening angle (CTOA) fracture criterion. The results of these analyses and of the experiments indicate that the pressure required to initiate stable crack growth and unstable crack growth in a shell subjected to internal pressure loads decreases as the initial crack length increases. The effects of crack length on the prebuckling, buckling and postbuckling responses of typical shells subjected to axial compression loads are also described. For this loading condition, the crack length was not allowed to increase as the load was increased. The results of the analyses and of the experiments indicate that the initial buckling load and collapse load for a shell subjected to axial compression loads decrease as the initial crack length increases. Initial buckling causes general instability or collapse of a shell for shorter initial crack lengths. Initial buckling is a stable local response mode for longer initial crack lengths. This stable local buckling response is followed by a stable postbuckling response, which is followed by general or overall instability of the shell.

Starnes, James H., Jr.↗

Buckling and Stable Tearing Responses of Unstiffened Aluminum Shells with Long Cracks

The results of an analytical and experimental study of the nonlinear response of thin, unstiffened, aluminum cylindrical shells with a long longitudinal crack are presented. The shells are analyzed with a nonlinear shell analysis code that accurately accounts for global and local structural response phenomena. Results are presented for internal pressure and for axial compression loads. The effect of initial crack length on the initiation of stable tearing and unstable crack growth in typical shells subjected to internal pressure loads is predicted using geo- metrically nonlinear elastic-plastic finite element analyses. The results of these analyses and of the experiments indicate that the pressure required to initiate stable tearing and unstable tearing in a shell subjected to internal pressure loads decreases as the crack length increases. The effects of crack length on the prebuckling, buckling and postbuckling responses of typical shells subjected to axial compression are also described. For this loading condition, the crack length is held constant. The results of the analyses illustrate the influence of crack length on shell buckling instabilities. The experimental and analytical results correlate well.

Starnes, James H., Jr.↗

A three-dimensional finite-element thermal/mechanical analytical technique for high-performance traveling wave tubes

Current research in high-efficiency, high-performance traveling wave tubes (TWT's) has led to the development of novel thermal/ mechanical computer models for use with helical slow-wave structures. A three-dimensional, finite element computer model and analytical technique used to study the structural integrity and thermal operation of a high-efficiency, diamond-rod, K-band TWT designed for use in advanced space communications systems. This analysis focused on the slow-wave circuit in the radiofrequency section of the TWT, where an inherent localized heating problem existed and where failures were observed during earlier cold compression, or 'coining' fabrication technique that shows great potential for future TWT development efforts. For this analysis, a three-dimensional, finite element model was used along with MARC, a commercially available finite element code, to simulate the fabrication of a diamond-rod TWT. This analysis was conducted by using component and material specifications consistent with actual TWT fabrication and was verified against empirical data. The analysis is nonlinear owing to material plasticity introduced by the forming process and also to geometric nonlinearities presented by the component assembly configuration. The computer model was developed by using the high efficiency, K-band TWT design but is general enough to permit similar analyses to be performed on a wide variety of TWT designs and styles. The results of the TWT operating condition and structural failure mode analysis, as well as a comparison of analytical results to test data are presented.

Bartos, Karen F.↗

A three-dimensional finite-element thermal/mechanical analytical technique for high-performance traveling wave tubes

Current research in high-efficiency, high-performance traveling wave tubes (TWT's) has led to the development of novel thermal/mechanical computer models for use with helical slow-wave structures. A three-dimensional, finite element computer model and analytical technique used to study the structural integrity and thermal operation of a high-efficiency, diamond-rod, K-band TWT designed for use in advanced space communications systems. This analysis focused on the slow-wave circuit in the radiofrequency section of the TWT, where an inherent localized heating problem existed and where failures were observed during earlier cold compression, or 'coining' fabrication technique that shows great potential for future TWT development efforts. For this analysis, a three-dimensional, finite element model was used along with MARC, a commercially available finite element code, to simulate the fabrication of a diamond-rod TWT. This analysis was conducted by using component and material specifications consistent with actual TWT fabrication and was verified against empirical data. The analysis is nonlinear owing to material plasticity introduced by the forming process and also to geometric nonlinearities presented by the component assembly configuration. The computer model was developed by using the high efficiency, K-band TWT design but is general enough to permit similar analyses to be performed on a wide variety of TWT designs and styles. The results of the TWT operating condition and structural failure mode analysis, as well as a comparison of analytical results to test data are presented.

Shalkhauser, Kurt A.↗

Clumpy Langmuir waves in type III radio sources - Comparison of stochastic-growth theory with observations

Detailed comparisons are made between the Langmuir-wave properties predicted by the recently developed stochastic-growth theory of type III sources and those observed by the plasma wave experiment on ISEE 3, after correcting for the main instrumental and selection effects. Analysis of the observed field-strength distribution confirms the theoretically predicted form and implies that wave growth fluctuates both spatially and temporally in sign and magnitude, leading to an extremely clumpy distribution of fields. A cutoff in the field-strength distribution is seen at a few mV/m, corresponding to saturation via nonlinear effects. Analysis of the size distribution of Langmuir clumps yields results in accord with those obtained in earlier work and with the size distribution of ambient density fluctuations in the solar wind. This confirms that the inhomogeneities in the Langmuir growth rate are determined by the density fluctuations and that these fluctuations persist during type III events.

Robinson, P. A.↗

Nonlinear acoustics experimental characterization of microstructure evolution in Inconel 617

Inconel 617 is a candidate material for the intermediate heat exchanger in a very high temperature reactor for the next generation nuclear power plant. This application will require the material to withstand fatigue-ratcheting interaction at temperatures up to 950°C. Therefore nondestructive evaluation and structural health monitoring are important capabilities. Acoustic nonlinearity (which is quantified in terms of a material parameter, the acoustic nonlinearity parameter, β) has been proven to be sensitive to microstructural changes in material. This research develops a robust experimental procedure to track the evolution of damage precursors in laboratory tested Inconel 617 specimens using ultrasonic bulk waves. The results from the acoustic non-linear tests are compared with stereoscope surface damage results. Therefore, the relationship between acoustic nonlinearity and microstructural evaluation can be clearly demonstrated for the specimens tested.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

New Directions in Asymptotically Stable Finite-dimensional Adaptive Control of Linear Distributed Parameter Systems

Distributed Parameter Systems (DPS), such as systems described by partial differential equations, require infinite-dimensional state space descriptions to correctly model their dynamical behavior. However, any adaptive control algorithm must be finite-dimensional in order to be implemented via on-line digital computers. Finite-dimensional adaptive control of linear DPS requires stability analysis of nonlinear, time-varying, infinite-dimensional systems. The structure of nonadaptive finite-dimensional control of linear DPS is summarized as it relates to the existence of limiting systems for adaptive control. Two candidate schemes for finite-dimensional adaptive control of DPS are described and critical issues in infinite-dimensional stability analysis are discussed, in particular, the invariance principle, center manifold theory, and relationships between input-output and internal stability.

Balas, M. J.↗

Analysis of Compressibility Effects and Nonlinear Property Variations in a Supercritical CO2 Mixing Layer

Calculations are performed using the large eddy simulation technique to qualitatively assess the effects of fluid compressibility and thermodynamic nonlinearities on the dynamics of a supercritical flow field. Resulting data are also used to interrogate instantaneous subfilter-scale fields, as modeled by the mixed-dynamic Smagorinsky closure and gradient diffusion model. A three-dimensional wall-resolved calculation of a spatially evolving mixing layer composed solely of carbon dioxide is performed. This is followed by a coarse grid calculation to examine model performance as a function of resolution. Theoretical analysis indicates that the partial derivatives of density with respect to pressure and temperature are important modulators of the temperature and pressure field evolution, respectively. Results pertaining to subfilter velocity-velocity and velocity-temperature fields indicate that the models are most stressed in the regions where compressibility and thermo-fluid nonlinearities dominate the physics.

42 ENGINEERING↗

Testing predictions of electron scale turbulent pedestal transport in two DIII-D ELMy H-modes

In this study, we present linear and nonlinear gyrokinetic analyses in the pedestal region of two DIII-D ELMy H-mode discharges using the CGYRO code. The otherwise matched discharges employ different divertor configurations to investigate the impact of varying recycling and particle source on pedestal profiles. Linear gyrokinetic simulations find electrostatic ion-scale instabilities (ion temperature gradient and trapped electron modes, ITG-TEM) are present just inside the top of the pedestal with growth rates that are enhanced significantly by parallel velocity shear. In the sharp gradient region, E x B shearing rates are comparable or larger than ion scale growth rates, suggesting the suppression of ITG-TEM modes in this region. Instead, the electron temperature profiles are found to be correlated with and just above the electron temperature gradient (ETG) instability thresholds. Using gradients varied within experimental uncertainties, nonlinear electron-scale gyrokinetic simulations predict electron heat fluxes from ETG turbulence, that when added to neoclassical ion thermal transport simulated by NEO, account for 30-60% of the total experimental heat flux. In addition, the neoclassical electron particle flux is found to contribute significantly to the experimental fluxes inferred from SOLPS-ITER analysis. Additional nonlinear gyrokinetic simulations are run varying input gradients to develop a threshold-based reduced model for ETG transport, finding a relatively simple dependence on η e = L ne /L Te . Predictive transport simulations are used to validate this pedestal-specific ETG model, in conjunction with a model for neoclassical particle transport. In both discharges, the predicted electron temperatures are always overpredicted, indicative of the insufficient stiffness in the ETG pedestal model to account for all of the experimental electron thermal transport. In the case of the closed divertor discharge with lower particle source, the predicted electron density is close to the experiment, consistent with the magnitude of neoclassical particle transport in that discharge. However, the density profiles are overpredicted in the open divertor discharge (larger particle source), due to insufficient model transport. In conclusion, the implications for other mechanisms accounting for the remainder of transport in the sharp gradient region in the two discharges are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Comparative analysis of algorithms for compensation of nonlinear distortions in fibre-optic links

We report the results of numerical modelling and comparative analysis of two algorithms for digital signal processing intended for signal demodulation in fibre-optic links (FOLs) operating in a nonlinear mode. The first algorithm is based on the use of nonlinear Schrödinger filters and the second one utilises the ‘reception in general with bit-by-bit decision making’ technique. It is found that the latter algorithm is more noiseproof, compared to the first one (with a SNR gain of about 0.5 dB), in single-channel fibre-optic links, and its operation does not require knowledge of the link parameters; therefore, the second algorithm can also be constructed for wavelength-division multiplexing systems. (paper)

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗