Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “nonlinear structural 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,315 records · Page 73
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

Shell stability analysis in a computer aided engineering (CAE) environment

The development of 'DISDECO', the Delft Interactive Shell DEsign COde is described. The purpose of this project is to make the accumulated theoretical, numerical and practical knowledge of the last 25 years or so readily accessible to users interested in the analysis of buckling sensitive structures. With this open ended, hierarchical, interactive computer code the user can access from his workstation successively programs of increasing complexity. The computational modules currently operational in DISDECO provide the prospective user with facilities to calculate the critical buckling loads of stiffened anisotropic shells under combined loading, to investigate the effects the various types of boundary conditions will have on the critical load, and to get a complete picture of the degrading effects the different shapes of possible initial imperfections might cause, all in one interactive session. Once a design is finalized, its collapse load can be verified by running a large refined model remotely from behind the workstation with one of the current generation 2-dimensional codes, with advanced capabilities to handle both geometric and material nonlinearities.

Arbocz, J.↗
NASA NTRS1996

Reliability, Risk and Cost Trade-Offs for Composite Designs

Risk and cost trade-offs have been simulated using a probabilistic method. The probabilistic method accounts for all naturally-occurring uncertainties including those in constituent material properties, fabrication variables, structure geometry and loading conditions. The probability density function of first buckling load for a set of uncertain variables is computed. The probabilistic sensitivity factors of uncertain variables to the first buckling load is calculated. The reliability-based cost for a composite fuselage panel is defined and minimized with respect to requisite design parameters. The optimization is achieved by solving a system of nonlinear algebraic equations whose coefficients are functions of probabilistic sensitivity factors. With optimum design parameters such as the mean and coefficient of variation (representing range of scatter) of uncertain variables, the most efficient and economical manufacturing procedure can be selected. In this paper, optimum values of the requisite design parameters for a predetermined cost due to failure occurrence are computationally determined. The results for the fuselage panel analysis show that the higher the cost due to failure occurrence, the smaller the optimum coefficient of variation of fiber modulus (design parameter) in longitudinal direction.

Shiao, Michael C.↗
NASA NTRS2004

Scatter in Carbon/Silicon Carbide (C/SiC) Composites Quantified

Carbon-fiber-reinforced silicon carbide matrix (C/SiC) composites processed by chemical vapor infiltration are candidate materials for aerospace thermal structures. Carbon fibers can retain properties at very high temperatures, but they are known to have poor oxidation resistance in adverse, high-temperature environments. Nevertheless, the combination of CVI-SiC matrix with higher stiffness and oxidation resistance, the interfacial coating, and additional surface-seal coating provides the necessary protection to the carbon fibers, and makes the material viable for high-temperature space applications operating under harsh environments. Furthermore, C/SiC composites, like other ceramic matrix composites (CMCs), exhibit graceful non-catastrophic failure because of various inherent energy dissipating mechanisms. The material exhibits nonlinearity in deformation even at very low stress levels. This is the result of the severe matrix microcracking present in the as processed composite because of large differences between the coefficients of thermal expansion of the fiber and the matrix. Utilization of these advanced composites in next generation space vehicles will require innovative structural configurations, updated materials, and refined analyses. Structural safety issues for these vehicles are in direct competition with performance and cost. One would have to quantify the uncertainties associated with the design using formal probabilistic methods. Specifically four fundamental aspects on which analyses are based-- (1) loading conditions, (2) material behavior, (3) geometrical configurations, and (4) structural connections between the composite components and baseline structure--are stochastic in nature. A direct way to formally account for uncertainties is to develop probabilistic structural analysis methods where all participating variables are described by appropriate probability density functions. The present work, however, focuses on analyzing the stochastic material behavior of these advanced composites using formal probabilistic analysis methods. Often, some of the desirable property characteristics that allow composites to offer advantages over conventional structural materials (like tailoring of composite properties) and the complexity are in fact responsible for their greater statistical variability and the requirements for more characterization tests. Composite properties are anisotropic as well, having different properties in different directions. This means that characterization of a property such as stiffness--which will vary greatly depending on the orientation of the fiber relative to the direction of the testing--must be repeated for several different directions and loading conditions. The fabrication process for composites also introduces statistical variations in properties and geometry. A composite part is produced in a number of steps, each of which introduces statistical variability. The matrix is usually produced from a combination of raw materials; and the fiber, which has its own set of properties, is often coated or surface treated, introducing yet another source of variability.

Murthy, Pappu L. N.↗
DOE OSTI2022

Reduced order models for thermal radiative transfer problems based on moment equations and data-driven approximations of the Eddington tensor

Here a new group of structure and asymptotic preserving reduced-order models (ROMs) for multidimensional nonlinear thermal radiative transfer (TRT) problems is presented. They are formulated by means of the nonlinear projective approach and data compression techniques. The nonlinear projection is applied to the Boltzmann transport equation (BTE) to derive a hierarchy of low-order moment equations. Approximation of the Eddington tensor that provides exact closure for the system of moment equations is found with projection-based data-driven methodologies. These include the (i) proper orthogonal decomposition (POD), (ii) dynamic mode decomposition (DMD) and (iii) a variant of the DMD. A parameterization is derived for this ROM for the temperature of radiation incoming to the problem domain (the radiation drive temperature). This parameterization is informed from results of a dimensionless study of the TRT problem. Analysis of the ROMs is performed on the classical Fleck-Cummings TRT multigroup test problem in 2D geometry with a radiation-driven Marshak wave. Numerical results are presented to demonstrate the performance of these ROMs for the simulation of evolving radiation and heat waves. Results show these models to be sufficiently accurate for practical computations with rather low-rank representations of the Eddington tensor. As the rank of the approximation is increased, the errors of solutions generated by the ROMs gradually decreases.

42 ENGINEERING↗
DOE OSTI2025

A machine-learning approach to measure 3D sample properties from 2D Transmission Electron Microscopy images

Transmission Electron Microscopy (TEM) is a powerful tool for the characterization of materials at the nanoscale; however, its inherent two-dimensional (2D) nature poses significant challenges to accurately measure three-dimensional (3D) properties. We introduce a supervised machine-learning model that predicts 3D structural information, such as sample thickness and curvature, from a series of conventional 2D TEM images. The model, a U-Net convolutional neural network, is trained on a large synthetic dataset generated from dynamical diffraction simulations that model TEM’s complex, nonlinear image formation, accounting for sample thickness and curvature. This physically realistic framework enables exploration of a broad parameter space impractical to sample experimentally. We demonstrate that the trained model has accurate predictions for experimental single-crystal silicon samples, achieving performance comparable to established measurement techniques. This work highlights the critical role of robust, simulation-based training in overcoming the limitations of real-world imaging artifacts and inconsistent sample geometries. By integrating machine learning with numerical simulations, we offer an efficient and scalable framework for quantitative TEM analysis, paving the way for more sophisticated 3D characterization of complex materials.

Dynamical diffraction↗
NASA NTRS2003

Use of GTE, Mozaic, Sonex, and UARS-MLS Data in Understanding Tropospheric Processes Critical to Model Development

During the first year we focused on the analysis of data collected on over 7600 commercial aircraft flights (the MOZAIC program). The aim was to further our understanding of the fundamental dynamical processes that drive mesoscale phenomena in the upper troposphere and lower stratosphere, and their effects on the advection of passive scalars. Through these studies we made the following findings. 2001]: We derived the Kolmogorov equation for the third-order velocity structure function on an f-plane. We showed how the sign of the function yields the direction of the energy cascade. The remarkable linearity of the measured off diagonal third-order structure function was studied. We suggested that the Coriolis term, which appears explicitly in this equation, may be crucial in understanding the observed kinetic energy spectra at scales larger than 100 km, instead of the nonlinear advection term as previously assumed. Also, we showed that decreases with latitude for mid- to high-latitudes, and tabulated the values. Ozone concentration structure functions were calculated by restricting the data points to approximately isentropic pairings. In this way we were able to make comparisons with scalar advection theories. We found that, at even the largest scales, there was no evidence for the simple scaling predicted for smooth advection/diffusion.

Newell, Reginald E.↗
DOE OSTI2023

Performance analysis of functionally graded multifunctional piezoelectric energy harvesting microgyroscopes

Functionally graded materials (FGMs) are composite materials with varying material properties in one or more directions. These materials possess distinct characteristics compared to their constituent components. The ability to control material distribution and compositions in FGMs offers improved constraints over the natural frequencies of the system, making them highly suitable for energy harvesting applications. In this study, we focus on a piezoelectric FGM composed of Platinum (Pt) and Lead Titanate Zirconate (PZT). By utilizing FGMs for energy harvesting, we simplify the system from a multilayer structure to a single-layer system, thereby increasing power density by reducing volume. In this work, a power law distribution is used to model the material variation throughout the thickness of the single-layered beam. The governing equations of motion and boundary conditions for FG energy harvesting microgyroscope are derived using Euler-Bernoulli beam theory, the constitutive piezoelectric principle, and the extended Hamilton's principle. To simulate the nonlinear motion of the FG energy harvester, we employ the differential quadrature method (DQM). Subsequently, an investigation is carried out to examine the effects of various factors including material distribution, platinum percentage, base rotation, DC voltage, and electrical load resistance on the energy harvesting microgyroscope with functionally graded materials. Here, the findings suggest that functionally graded energy harvesting gyroscopes can be adjusted to obtain effective energy harvesting and sensing capabilities by carefully selecting the material distribution, input voltage, and electrical load resistance.

42 ENGINEERING↗
DOE OSTI2026

Evolution of the Antarctic Ice Sheet from 2000–2300 and beyond: model sensitivity and uncertainty analysis using MPAS-Albany Land Ice

We present a description of the Antarctic Ice Sheet model configuration submitted to the ISMIP6-Antarctica-2300 experiment using the MPAS-Albany Land Ice model, along with three new sets of simulations: (1) a set of extended simulations to 2500 for three forced experiments and to 2775 for the control experiment; (2) a sensitivity analysis of our model configuration to parameters controlling basal sliding and sub-shelf melt, and to model structural choices including the choice of the energy and stress balances; and (3) a 72-member ensemble run on graphics processing units (GPUs) and analysis of variance to determine the primary sources of uncertainty in our ice-sheet model projections. Our extended simulations predict rapid retreat beginning after 2300 for SSP1-2.6 forcing and after 2500 for present-day (control) forcing, primarily in the Amundsen Sea Embayment. We find that varying the sub-shelf melt parameter between the 5th to 95th percentile values for a mean-Antarctic calibration target results in an up to ∼ ± 40 % change in sea-level contribution relative to our baseline simulations that used the median value. Using a linear basal sliding law reduces sea-level contribution by 51 %–73 % relative to our baseline nonlinear sliding law with an exponent of 1/5. When using basal sliding law exponents of 1/3 and 1/10, the overall difference from our baseline simulations at 2300 is on the order of 10 %. The Amundsen Sea Embayment region displays a strongly non-linear dependence of mass loss on the sliding law exponent, with no discernible relationship between the sliding law exponent and the mass loss by 2300, while the sectors feeding the Ross and Filchner-Ronne ice shelves exhibit more mass loss with a more-plastic sliding law. Our model fidelity sensitivity experiments reveal a 9 %–31 % increase in sea-level contribution when using a depth-integrated stress balance approximation relative to our three-dimensional solver, while using a fixed-in-time temperature field increases sea-level contribution by 14 %–88 % relative to two thermomechanically coupled configurations. Our 72-member ensemble and analysis of variance show that the uncertainty in long-term projections is dominated by the choice of Earth system model forcing and the presence or absence of hydrofracture forcing, rather than uncertainty in sliding and sub-shelf melt parameters.

58 GEOSCIENCES↗
NASA NTRS1993

Buoyancy-driven flows of a radiatively participating fluid in a vertical cylinder heated from below

The effect of radiative energy transport on the onset and evolution of natural convective flows is studied in a Rayleigh-Benard system. Steady, axisymmetric flows of a radiatively participating fluid contained in a rigid-walled, vertical cylinder which is heated on the base, cooled on top, and insulated on the side wall are calculated by using the Galerkin FEM. Bifurcation analysis techniques are used to investigate the changes in the flow structure due to internal radiation. The results of this two-parameter study - where the Rayleigh number, Ra and optical thickness, tau, are varied - apply to fluids ranging from opaque to nearly transparent with respect to IR radiation. For any nonopaque fluid, internal radiation eliminates the static state that without radiation, exists for all values of the Rayleigh number. This heat transfer mechanism also destroys a symmetry of the system that relates clockwise and counterclockwise flows. The connectivity between characteristic flow families and the range of Ra where families are stable are found to depend greatly on tau. Results demonstrate the inadequacy of characterizing the behavior of this system using simple notions of radiative transfer in optically thick or thin media; the nonlinear interaction of radiation and flow are far more complicated than these asymptotic limits would imply.

Salinger, A. G.↗
NASA NTRS1987

A Comparative Study of Some Dynamic Stall Models

Three semi-empirical aerodynamic stall models are compared with respect to their lift and moment hysteresis loop prediction, limit cycle behavior, easy implementation, and feasibility in developing the parameters required for stall flutter prediction of advanced turbines. For the comparison of aeroelastic response prediction including stall, a typical section model and a plate structural model are considered. The response analysis includes both plunging and pitching motions of the blades. In model A, a correction to the angle of attack is applied when the angle of attack exceeds the static stall angle. In model B, a synthesis procedure is used for angles of attack above static stall angles and the time history effects are accounted through the Wagner function. In both models the life and moment coefficients for angle of attack below stall are obtained from tabular data for a given Mach number and angle of attack. In model C, referred to an the ONERA model, the life and moment coefficients are given in the form of two differential equations, one for angles below stall, and the other for angles above stall. The parameters of those equations are nonlinear functions of the angle of attack.

Reddy, T. S. R.↗
NASA NTRS1999

A Critical Examination of Wind-Wave Spectral Functional Form

Traditionally, data from random ocean waves are presented in spectral functions. The spectra are the result of Fourier analysis. Fourier spectral analysis has dominated data analysis for, at least, the last hundred years. It has been the standard method for is examining the global amplitude-frequency distributions. Although Fourier transform valid under extremely general conditions, there are some crucial restrictions for the Fourier spectral analysis. The system must be linear, and the data must be stationary- otherwise, the resulting spectrum will make little physical sense. The stationarity requirement is also a common required criterion for most of other available data analysis methods. Nevertheless, few, if any, natural phenomena are linear and stationary. To compound these complications is the imperfection of our probes or numerical schemes the interactions of the imperfect probes even with a perfect linear system can make the final data nonlinear. Furthermore, all the available data are usually of finite duration. Under these conditions, Fourier analysis is of limited use, For lack of alternatives, however, Fourier analysis is still used to process such data. The loose application of Fourier analysis and the insouciant adoption of the stationary and linear assumptions may lead to misleading conclusions. Ocean waves are know to be nonlinear, and the wind system generating the wave field are seldom stationary- As a result, the traditional examination of the spectral form hardly made physical sense. A new method for analyzing nonlinear and nonstationary data has been developed. The key part is the Empirical Mode Decomposition (EMD) method with which any complicated data set can be decomposed into a finite and often small number of Intrinsic Mode Functions (IMF) that serve as the basis of the representation of the data, This decomposition method is adaptive, and, therefore, highly efficient. The IMFs admit well-behaved Hilbert transforms, and yield instantaneous energy and frequency as functions of time that give sharp identifications of imbedded structures. The final presentation of the results is an energy-frequency-time distribution, designated as the Hilbert Spectrum. Among the main conceptual innovations is the introduction of the instantaneous frequencies for complicated data sets, which eliminate the need of spurious harmonics to represent nonlinear and nonstationary signals. The spectral forms of the ocean waves are very different. This new method offers the first physical meaningful view of the wave spectrum. Data from laboratory and field will be presented to illustrate the differences.

Huang, Norden E.↗
NASA NTRS1992

Observations of short large-amplitude magnetic structures at a quasi-parallel shock

Results of a detailed analysis of short large-amplitude magnetic structures (SLAMS) observed at an encounter of the quasi-parallel blow shock by the AMPTE UKS and IRM satellites are presented. Isolated SLAMS, surrounded by solar wind conditions, and embedded SLAMS, which lie within or form the boundary with regions of significant heating and deceleration, are identified. The duration, polarization, and other characteristics of SLAMS are all consistent with their growth directly out of the ULF wave field, including the common occurrence of an attached whistler as found in ULF shocklets. The plasma rest frame propagation speeds and two-spacecraft time delays for all cases show that the SLAMS attempt to propagate upstream against the oncoming flow, but are convected back downstream. The speeds and delays vary systematically with SLAMS amplitude in the way anticipated from nonlinear wave theory, as do their polarization features. Inter-SLAMS regions and boundary regions with solar wind contain hot deflected ions of lesser density than within the SLAMS.

Schwartz, Steven J.↗
NASA NTRS2001

HITEMP Material and Structural Optimization Technology Transfer

The feasibility of adding viscoelasticity and the Generalized Method of Cells (GMC) for micromechanical viscoelastic behavior into the commercial HyperSizer structural analysis and optimization code was investigated. The viscoelasticity methodology was developed in four steps. First, a simplified algorithm was devised to test the iterative time stepping method for simple one-dimensional multiple ply structures. Second, GMC code was made into a callable subroutine and incorporated into the one-dimensional code to test the accuracy and usability of the code. Third, the viscoelastic time-stepping and iterative scheme was incorporated into HyperSizer for homogeneous, isotropic viscoelastic materials. Finally, the GMC was included in a version of HyperSizer. MS Windows executable files implementing each of these steps is delivered with this report, as well as source code. The findings of this research are that both viscoelasticity and GMC are feasible and valuable additions to HyperSizer and that the door is open for more advanced nonlinear capability, such as viscoplasticity.

Collier, Craig S.↗
NASA NTRS1984

Natural convection in steady solidification - Finite element analysis of a two-phase Rayleigh-Benard problem

Galerkin finite-element approximations and Newton's method for solving free boundary problems are combined with computer-implemented techniques from nonlinear perturbation analysis to study solidification problems with natural convection in the melt. The Newton method gives rapid convergence to steady state velocity, temperature and pressure fields and melt-solid interface shapes, and forms the basis for algebraic methods for detecting multiple steady flows and assessing their stability. The power of this combination is demonstrated for a two-phase Rayleigh-Benard problem composed of melt and solid in a veritical cylinder with the thermal boundary conditions arranged so that a static melt with a flat melt-solid interface is always a solution. Multiple cellular flows bifurcating from the static state are detected and followed as Rayleigh number is varied. Changing the boundary conditions to approach those appropriate for the vertical Bridgman solidification system causes imperfections that eliminate the static state. The flow structure in the Bridgman system is related to those for the Rayleigh-Benard system by a continuous evolution of the boundary conditions.

Chang, C. J.↗
DOE OSTI2022

Deconvolving the roles of E × B shear and pedestal structure in the energy confinement quality of super H-mode experiements

Analysis of “super H-mode” experiments on DIII-D has put forward that high plasma toroidal rotation, not high pedestal, plays the essential role in achieving energy confinement quality H 98y2 >>1. Recently, super H-mode experiments with variable input torque have confirmed that high rotation shear discharges have very high levels of H 98y2 (>1.5), independent of the pedestal height, and that high pedestal discharges with low rotation shear have levels of H 98y2 only slightly above 1 (≤1.2). Although some increase in stored energy with higher pedestal occurs, the energy confinement quality mainly depends on the toroidal rotation shear, which varies according to different levels of injected neutral beam torque per particle. Here, quasi-linear gyrofluid modeling achieves a good match of the experiment when including the E×B shear; without including plasma rotation, the modeling predicts a confinement quality consistent with the empirical observation of H 98y2 ~1.2 at low rotation. Nonlinear gyrokinetic transport modeling shows that the effect of E×B turbulence stabilization is far larger than other mechanisms, such as the so-called hot-ion stabilization (T i /T e ) effect. Consistent with these experimental and modeling results are previous simulations of the ITER Baseline Scenario using a super H-mode pedestal solution, which showed the potential to exceed the Q=10 target if the pedestal density could be increased above the Greenwald limit. A close look at these simulations reveals that the predicted energy confinement quality is below 1 even at the highest pedestal pressure. The improvement in Q at higher pedestal density is due to the improved fusion power generation at the higher core density associated with higher pedestal density, not to an improved energy confinement quality.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1993

Dispersed ion structures at the poleward edge of the auroral oval: Low-altitude observations and numerical modeling

We have compared the AUREOL 3 (A3) observations of auroral ion precipitation, particularly ion beams, with the results from the global kinetic model of magnetotail plasma of Ashour-Abdalla et al. (1993). We have identified 101 energetic keV H(+) velocity dispersed precipitating ion structures (VDIS) with fluxes above 10(exp -3) ergs./sq cm./s in the A3 record between the end of 1981 and mid-1984. These beams display a systematic increase in energy with increasing latitude and were observed in a narrow region within less than 1 deg in latitude of the polar cap boundary. The VDIS are the most distinctive feature in the auroral zone of the plasma sheet boundary layer. We report first on a statistical analysis of the possible ralationships between magnetic activity or substorm phase and the VDIS properties. Our particle simulations of the precipitating ions have been extended by using a series of modified versions of the Tsyganenko (1989) magnetic field model and by varying the cross-magnetosphere electric field. In the simulations, plasma from a mantle source is subject to strong nonlinear acceleration, forming beams which flow along the PSBL. Only 3 to 4% of these beams precipitate into the ionosphere to form the VDIS while the majority return to the equatorial plane after mirroring and form the thermalized central plasma sheet. The final energy and the dispersion of the beams in the model depend on the amplitude of the cross-tail electric field. Two unsual observations of low-energy (less than 5 keV) O(+) VDIS, shifted by 4 deg 5 deg in invariant latitude equatorward of H(+) VDIS are analyzed in detail. The sparsity of such O(+) events and the absence of the changes in the flux and frequency of occurrence indicate a solar wind origin for the plasma. Finally, large-scale kinetic modeling, even with its simplifications and assumptions (e.g., static magnetic field, solar wind source), reproduces low-altitude auroral ion features fairly well; it may therefore be presented as an appropriate framework into which data on energization and transport of the hot plasma, obtained in the equatorial plane, could be inserted in the near future.

Bosqued, J. M.↗
NASA NTRS2006

Flutter of High-Speed Civil Transport Flexible Semispan Model: Time-Frequency Analysis

Time/frequency analysis of fluctuations measured by pressure taps and strain gauges in the experimental studies of the flexible semispan model of a high-speed civil transport wing configuration is performed. The interest is in determining the coupling between the aerodynamic loads and structural motions that led to the hard flutter conditions and loss of the model. The results show that, away from the hard flutter point, the aerodynamic loads at all pressure taps near the wing tip and the structural motions contained the same frequency components. On the other hand, in the flow conditions leading to the hard flutter, the frequency content of the pressure fluctuations near the leading and trailing edges varied significantly. This led to contribution to the structural motions over two frequency ranges. The ratio of these ranges was near 2:1, which suggests the possibility of nonlinear structural coupling.

Chabalko, Christopher C.↗
NASA NTRS2011

Enhanced Spectral Anisotropies Near the Proton-Cyclotron Scale: Possible Two-Component Structure in Hall-FLR MHD Turbulence Simulations

Recent analysis of the magnetic correlation function of solar wind fluctuations at 1 AU suggests the existence of two-component structure near the proton-cyclotron scale. Here we use two-and-one-half dimensional and three-dimensional compressible MHD models to look for two-component structure adjacent the proton-cyclotron scale. Our MHD system incorporates both Hall and Finite Larmor Radius (FLR) terms. We find that strong spectral anisotropies appear adjacent the proton-cyclotron scales depending on selections of initial condition and plasma beta. These anisotropies are enhancements on top of related anisotropies that appear in standard MHD turbulence in the presence of a mean magnetic field and are suggestive of one turbulence component along the inertial scales and another component adjacent the dissipative scales. We compute the relative strengths of linear and nonlinear accelerations on the velocity and magnetic fields to gauge the relative influence of terms that drive the system with wave-like (linear) versus turbulent (nonlinear) dynamics.

Ghosh, Sanjoy↗