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,495 records · Page 83
All sourcesNASA NTRSDOE OSTI
NASA NTRS1994

Discrete crack growth analysis methodology for through cracks in pressurized fuselage structures

A methodology for simulating the growth of long through cracks in the skin of pressurized aircraft fuselage structures is described. Crack trajectories are allowed to be arbitrary and are computed as part of the simulation. The interaction between the mechanical loads acting on the superstructure and the local structural response near the crack tips is accounted for by employing a hierarchical modeling strategy. The structural response for each cracked configuration is obtained using a geometrically nonlinear shell finite element analysis procedure. Four stress intensity factors, two for membrane behavior and two for bending using Kirchhoff plate theory, are computed using an extension of the modified crack closure integral method. Crack trajectories are determined by applying the maximum tangential stress criterion. Crack growth results in localized mesh deletion, and the deletion regions are remeshed automatically using a newly developed all-quadrilateral meshing algorithm. The effectiveness of the methodology and its applicability to performing practical analyses of realistic structures is demonstrated by simulating curvilinear crack growth in a fuselage panel that is representative of a typical narrow-body aircraft. The predicted crack trajectory and fatigue life compare well with measurements of these same quantities from a full-scale pressurized panel test.

Potyondy, David O.↗
NASA NTRS1999

Potential Flow Interactions With Directional Solidification

The effect of convective melt motion on the growth of morphological instabilities in crystal growth has been the focus of many studies in the past decade. While most of the efforts have been directed towards investigating the linear stability aspects, relatively little attention has been devoted to experimental and numerical studies. In a pure morphological case, when there is no flow, morphological changes in the solid-liquid interface are governed by heat conduction and solute distribution. Under the influence of a convective motion, both heat and solute are redistributed, thereby affecting the intrinsic morphological phenomenon. The overall effect of the convective motion could be either stabilizing or destabilizing. Recent investigations have predicted stabilization by a flow parallel to the interface. In the case of non-parallel flows, e.g., stagnation point flow, Brattkus and Davis have found a new flow-induced morphological instability that occurs at long wavelengths and also consists of waves propagating against the flow. Other studies have addressed the nonlinear aspects (Konstantinos and Brown, Wollkind and Segel)). In contrast to the earlier studies, our present investigation focuses on the effects of the potential flow fields typically encountered in Hele-Shaw cells. Such a Hele-Shaw cell can simulate a gravity-free environment in the sense that buoyancy-driven convection is largely suppressed, and hence negligible. Our interest lies both in analyzing the linear stability of the solidification process in the presence of potential flow fields, as well as in performing high-accuracy nonlinear simulations. Linear stability analysis can be performed for the flow configuration mentioned above. It is observed that a parallel potential flow is stabilizing and gives rise to waves traveling downstream. We have built a highly accurate numerical scheme which is validated at small amplitudes by comparing with the analytically predicted results for the pure morphological case. We have been able to observe nonlinear effects at larger times. Preliminary results for the case when flow is imposed also provide good validation at small amplitudes.

Buddhavarapu, Sudhir S.↗
NASA NTRS1999

Progressive Failure Analysis Methodology for Laminated Composite Structures

A progressive failure analysis method has been developed for predicting the failure of laminated composite structures under geometrically nonlinear deformations. The progressive failure analysis uses C(exp 1) shell elements based on classical lamination theory to calculate the in-plane stresses. Several failure criteria, including the maximum strain criterion, Hashin's criterion, and Christensen's criterion, are used to predict the failure mechanisms and several options are available to degrade the material properties after failures. The progressive failure analysis method is implemented in the COMET finite element analysis code and can predict the damage and response of laminated composite structures from initial loading to final failure. The different failure criteria and material degradation methods are compared and assessed by performing analyses of several laminated composite structures. Results from the progressive failure method indicate good correlation with the existing test data except in structural applications where interlaminar stresses are important which may cause failure mechanisms such as debonding or delaminations.

Sleight, David W.↗
DOE OSTI2015

Transient and chaotic low-energy transfers in a system with bistable nonlinearity

The low-energy dynamics of a two-dof system composed of a grounded linear oscillator coupled to a lightweight mass by means of a spring with both cubic nonlinear and negative linear components is investigated. The mechanisms leading to intense energy exchanges between the linear oscillator, excited by a low-energy impulse, and the nonlinear attachment are addressed. For lightly damped systems, it is shown that two main mechanisms arise: Aperiodic alternating in-well and cross-well oscillations of the nonlinear attachment, and secondary nonlinear beats occurring once the dynamics evolves solely in-well. The description of the former dissipative phenomenon is provided in a two-dimensional projection of the phase space, where transitions between in-well and cross-well oscillations are associated with sequences of crossings across a pseudo-separatrix. Whereas the second mechanism is described in terms of secondary limiting phase trajectories of the nonlinear attachment under certain resonance conditions. The analytical treatment of the two aformentioned low-energy transfer mechanisms relies on the reduction of the nonlinear dynamics and consequent analysis of the reduced dynamics by asymptotic techniques. Direct numerical simulations fully validate our analytical predictions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1989

Advanced applications of BEM to inelastic analysis of solids

Newly evolved boundary element solution algorithms are utilized for solving a range of elastoplastic and thermoplastic problems. The first algorithm is the typical iterative procedure in which the unknown boundary solution and the initial stress (or strain) rates are found together in an incremental iterative fashion. The second one is a variable stiffness type approach in which the incremental boundary solution is obtained in a direct (noniterative) manner. The third procedure differs from the two prior algorithms in that volume integration is not needed to incorporate the nonlinear effects in the analysis. Instead, initial stress rates are introduced in the boundary element system by way of particular integrals. The methodology is used in a general purpose, multiregion system that employs quadratic isoparametric shape functions to model the geometry and field variables of the body and can admit up to 15 substructured areas of different material properties.

Banerjee, P. K.↗
DOE OSTI2012

Bounce and cyclic cosmology in extended nonlinear massive gravity

We investigate non-singular bounce and cyclic cosmological evolutions in a universe governed by the extended nonlinear massive gravity, in which the graviton mass is promoted to a scalar-field potential. The extra freedom of the theory can lead to certain energy conditions violations and drive cyclicity with two different mechanisms: either with a suitably chosen scalar-field potential under a given Stückelberg-scalar function, or with a suitably chosen Stückelberg-scalar function under a given scalar-field potential. Our analysis shows that extended nonlinear massive gravity can alter significantly the evolution of the universe at both early and late times.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2011

TIME-FREQUENCY ANALYSIS OF THE SUPERORBITAL MODULATION OF THE X-RAY BINARY SMC X-1 USING THE HILBERT-HUANG TRANSFORM

The high-mass X-ray binary SMC X-1 exhibits a superorbital modulation with a dramatically varying period ranging between {approx}40 days and {approx}60 days. This research studies the time-frequency properties of the superorbital modulation of SMC X-1 based on the observations made by the All-Sky Monitor (ASM) onboard the Rossi X-ray Timing Explorer (RXTE). We analyzed the entire ASM database collected since 1996. The Hilbert-Huang transform (HHT), developed for non-stationary and nonlinear time-series analysis, was adopted to derive the instantaneous superorbital frequency. The resultant Hilbert spectrum is consistent with the dynamic power spectrum as it shows more detailed information in both the time and frequency domains. The RXTE observations show that the superorbital modulation period was mostly between {approx}50 days and {approx}65 days, whereas it changed to {approx}45 days around MJD 50,800 and MJD 54,000. Our analysis further indicates that the instantaneous frequency changed to a timescale of hundreds of days between {approx}MJD 51,500 and {approx}MJD 53,500. Based on the instantaneous phase defined by HHT, we folded the ASM light curve to derive a superorbital profile, from which an asymmetric feature and a low state with barely any X-ray emissions (lasting for {approx}0.3 cycles) were observed. We also calculated the correlation between the mean period and the amplitude of the superorbital modulation. The result is similar to the recently discovered relationship between the superorbital cycle length and the mean X-ray flux for Her X-1.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2011

A Novel Approach to Anharmonicity for a Wealth of Applications in Nonlinear Science Technologies

We present a new theory of the anharmonicity of nonlinear oscillations that are exhibited by many physical systems. New physical quantities are introduced that describe the departure from linear or harmonic behavior and as far as extremely anharmonic situations. In order to solve the nonlinear phase equation, the key notion of our theory, which controls the anharmonic behavior, a new and fascinating nonlinear trigonometry is designed. These results provide a general and accurate yet compact description of such signals, by far better than the Fourier description, both quantitatively and qualitatively and will benefit many application fields.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

ESTABLISHING ROCK OUTCROP TIME SERIES AT VARIOUS DEPTHS IN A NONLINEAR SOIL COLUMN

Nonlinear soil-structure interaction (NLSSI) analysis methods are currently being developed at Idaho National Laboratory (INL). Appendix B in ASCE 4 (2016) includes a framework for performing NLSSI analysis. For NLSSI analysis to be widely implemented in industry, it is important to demonstrate that accurate response in a dynamically loaded, nonlinear soil column can be produced. In performing NLSSI analysis, one approach is to consider vertically propagating shear and compressive waves that pass through a nonlinear soil column and into a structure. Using finite element techniques, the structure can be placed in a finite volume of nonlinear soil with a rock outcrop seismic time series applied to the bottom of the nonlinear soil. The necessary depth of the nonlinear soil mesh is dependent on the vertical soil stiffness when the model is in static and dynamic equilibrium. Run time, and therefore size of the meshed soil domain, is an important consideration when developing the NLSSI model. An approach is presented in this paper that establishes compatible rock outcrop seismic time series for reduced nonlinear soil mesh depths. To apply the approach, an example is developed using a reasonable soil column with a reasonable top-of-soil response spectrum. The approach is an iterative method which is applied in two steps. The first step is used to establish a reasonable rock outcrop time series at the highest depth considered. The second step is used to establish compatible rock outcrop time series at shallower depths.

42 ENGINEERING↗
NASA NTRS2000

The Generation of Smooth High Speed Solar Wind From Plume-Interplume Mixing

Ulysses has shown that fast solar wind is extremely smooth, with a variance of less than 5%, in contrast to slow wind with a variance of approximately 30%. Now UVCS has produced the surprising result that the flow speed within coronal holes, the source of fast wind, is not at all smooth. Specifically, Giordano et al. (ApJ, v531, L79-L82, 2000) report that at 1.7 R(sub SUN) the interplume flow speed is typically more than twice the plume flow speed. Other less direct evidence supports this same result, with speeds from less than 300 to over 1000 km/s reported at approximately 5 R(sub SUN). This presents the paradox of how strongly differing plume and interplume flow speeds can exist near the Sun and be absent far from the Sun. The only answer is that plume and interplume material or momentum must be strongly mixed and that the mixing must occur mainly inside 0.3 AU to be consistent with Helios observations of smooth fast wind. Pressure balance structures (PBSs) and He abundance anomalies (Reisenfeld, et al., GRL, v26(13), 1805-1808, 1999) have been identified as interplanetary remnants of plumes, implying momentum mixing is the dominant coronal process. One possible source for plume/interplume momentum mixing is MHD Kelvin-Helmholtz (KH) shear instabilities occurring on the velocity shear interfaces. The velocity shear is a source of free energy and KH fluctuations could, through nonlinear cascade, provide the forcing required for the plasma oscillations (Cranmer, ApJ, v532, 1197-1208, 2000) reported to exist in coronal holes. The physical properties in coronal holes are now sufficiently well known that we can show plume/interplume shear interfaces become unstable to the KH instability at 5 - 10 R(sub SUN). The KH dispersion relation can be used to analyze marginal stability, the most unstable wavelengths, and linear growth rates. Numerical simulations can be used to verify results from the linear analysis and study the nonlinear development of KH modes. Here we will describe the marginal stability criterion for the KH instability, how and where this condition is met in the corona, and the possible character of the resultant fluctuations.

Suess, Steve↗
NASA NTRS1985

Cost Considerations in Nonlinear Finite-Element Computing

Conference paper discusses computational requirements for finiteelement analysis using quasi-linear approach to nonlinear problems. Paper evaluates computational efficiency of different computer architecturtural types in terms of relative cost and computing time.

Utku, S.↗
NASA NTRS1995

On the nonlinear stability of the unsteady, viscous flow of an incompressible fluid in a curved pipe

The stability of the flow of an incompressible, viscous fluid through a pipe of circular cross-section curved about a central axis is investigated in a weakly nonlinear regime. A sinusoidal pressure gradient with zero mean is imposed, acting along the pipe. A WKBJ perturbation solution is constructed, taking into account the need for an inner solution in the vicinity of the outer bend, which is obtained by identifying the saddle point of the Taylor number in the complex plane of the cross-sectional angle co-ordinate. The equation governing the nonlinear evolution of the leading order vortex amplitude is thus determined. The stability analysis of this flow to periodic disturbances leads to a partial differential system dependent on three variables, and since the differential operators in this system are periodic in time, Floquet theory may be applied to reduce this system to a coupled infinite system of ordinary differential equations, together with homogeneous uncoupled boundary conditions. The eigenvalues of this system are calculated numerically to predict a critical Taylor number consistent with the analysis of Papageorgiou. A discussion of how nonlinear effects alter the linear stability analysis is also given, and the nature of the instability determined.

Shortis, Trudi A.↗
DOE OSTI2020

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↗
DOE OSTI2017

Light Water Reactor Sustainability Program Interrogation of Alkali-Silica Reaction Degraded Concrete Samples using Acoustic and Thermal Techniques to Support Development of a Structural Health Monitoring Framework

Assessment and management of aging concrete structures in nuclear power plants require a more systematic approach than simple reliance on existing code margins of safety. Health monitoring of concrete structures is performed to understand the current condition of a structure based on heterogeneous measurements and then produce high-confidence actionable information regarding structural integrity. This information can then be used to support operational and maintenance decisions. The proposed diagnosis and prognosis framework consists of four elements: health monitoring, data analytics, uncertainty quantification, and prognosis. The objectives of this ongoing research project focuses on health monitoring and data analytics of concrete slab cast, which is cured using four pockets of reactive aggregates: pure silica, wells, placitas, and spratt. A controlled concrete slab with four aggregates is cast at Vanderbilt University in representative conditions to accelerate degradation due to alkali-silica reaction (ASR). Four nondestructive evaluation (NDE) techniques (thermography, mechanical deformation measurements, nonlinear impact resonance-acoustic spectroscopy, and vibro-acoustic modulation) are used to examine the degraded concrete slab. Heterogeneous data from these techniques are analyzed for assessing ASRrelated degradation in concrete samples. In addition, this report discusses thermography assessment of three cement bricks cured under different conditions to accelerate ASR-related degradation. The application of the digital image correlation technique to study ASR-related degradation on a large concrete sample at the University of Tennessee Knoxville is discussed and initial observations are presented in this report. The results from the data analysis of the NDE techniques on concrete slab cured at Vanderbilt University are discussed in detail in this report. The major outcome of the analysis suggests that the nonlinear impact resonance acoustic spectroscopy is able to detect potential degradation in a concrete slab but is unable to locate or isolate the degradation. Alternatively, the vibro-acoustic technique is able to detect degradation and locate the damage in the concrete slab cast. The outcomes supports development of a structural health monitoring framework.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2019

A data-driven approach to construct a quantitative relationship between microstructural features of fatigue cracks and contact acoustic nonlinearity

This study demonstrates the feasibility of a data-driven approach to construct a quantitative relationship between nonlinear acoustic parameters and microstructural features of contact interfaces. The near-surface nonlinearity is measured using dynamic acousto-elastic testing (DAET) with a surface wave probe, while the microstructural features are extracted from scanning electron microscopy (SEM) images of fatigue cracks. Four aluminum alloy samples, each having a fatigue crack are prepared. Six local nonlinearity parameters are measured at different locations along the crack propagation direction. A total of 40 local measurements are acquired. A principal component analysis (PCA) reveals that all six nonlinearity parameters are correlated and hence can be replaced by one principal component (PC). Fifteen crack micro-geometrical features at each measurement point were extracted from the SEM images. Regression analysis is used to relate the PC of the nonlinearity parameters to the microstructural features at the crack interface. We compare three regression models that take variable selection into account: stepwise multiple linear regression (MLR), stepwise principal component regression (PCR), and least absolute shrinkage and selection operator (LASSO). Despite having different principles, the three predictive models identify two features as the most significant in predicting the interface nonlinearity: the crack aperture (opening) distribution and the distance to the crack tip. The differences between the three models and the physical interpretation of the data-driven predictions are discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2016

TU-FG-BRB-03: Basis Vector Model Based Method for Proton Stopping Power Estimation From Experimental Dual Energy CT Data

Purpose: This work aims at reducing the uncertainty in proton stopping power (SP) estimation by a novel combination of a linear, separable basis vector model (BVM) for stopping power calculation (Med Phys 43:600) and a statistical, model-based dual-energy CT (DECT) image reconstruction algorithm (TMI 35:685). The method was applied to experimental data. Methods: BVM assumes the photon attenuation coefficients, electron densities, and mean excitation energies (I-values) of unknown materials can be approximated by a combination of the corresponding quantities of two reference materials. The DECT projection data for a phantom with 5 different known materials was collected on a Philips Brilliance scanner using two scans at 90 kVp and 140 kVp. The line integral alternating minimization (LIAM) algorithm was used to recover the two BVM coefficient images using the measured source spectra. The proton stopping powers are then estimated from the Bethe-Bloch equation using electron densities and I-values derived from the BVM coefficients. The proton stopping powers and proton ranges for the phantom materials estimated via our BVM based DECT method are compared to ICRU reference values and a post-processing DECT analysis (Yang PMB 55:1343) applied to vendorreconstructed images using the Torikoshi parametric fit model (tPFM). Results: For the phantom materials, the average stopping power estimations for 175 MeV protons derived from our method are within 1% of the ICRU reference values (except for Teflon with a 1.48% error), with an average standard deviation of 0.46% over pixels. The resultant proton ranges agree with the reference values within 2 mm. Conclusion: Our principled DECT iterative reconstruction algorithm, incorporating optimal beam hardening and scatter corrections, in conjunction with a simple linear BVM model, achieves more accurate and robust proton stopping power maps than the post-processing, nonlinear tPFM based DECT analysis applied to conventional reconstructions of low and high energy scans. Funding Support: NIH R01CA 75371; NCI grant R01 CA 149305.

60 APPLIED LIFE SCIENCES↗
NASA NTRS1996

Nonlinear Dynamics of a Foil Bearing Supported Rotor System: Simulation and Analysis

Foil bearings provide noncontacting rotor support through a number of thin metal strips attached around the circumference of a stator and separated from the rotor by a fluid film. The resulting support stiffness is dominated by the characteristics of the foils and is a nonlinear function of the rotor deflection. The present study is concerned with characterizing this nonlinear effect and investigating its influence on rotordynamical behavior. A finite element model is developed for an existing bearing, the force versus deflection relation characterized, and the dynamics of a sample rotor system are studied. Some conclusions are discussed with regard to appropriate ranges of operation for such a system.

Li, Feng↗
NASA NTRS2003

Effect of Forcing Function on Nonlinear Acoustic Standing Waves

Nonlinear acoustic standing waves of high amplitude have been demonstrated by utilizing the effects of resonator shape to prevent the pressure waves from entering saturation. Experimentally, nonlinear acoustic standing waves have been generated by shaking an entire resonating cavity. While this promotes more efficient energy transfer than a piston-driven resonator, it also introduces complicated structural dynamics into the system. Experiments have shown that these dynamics result in resonator forcing functions comprised of a sum of several Fourier modes. However, previous numerical studies of the acoustics generated within the resonator assumed simple sinusoidal waves as the driving force. Using a previously developed numerical code, this paper demonstrates the effects of using a forcing function constructed with a series of harmonic sinusoidal waves on resonating cavities. From these results, a method will be demonstrated which allows the direct numerical analysis of experimentally generated nonlinear acoustic waves in resonators driven by harmonic forcing functions.

Finkheiner, Joshua R.↗