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 937 records · Page 52

Adaptive methods for nonlinear structural dynamics and crashworthiness analysis

The objective is to describe three research thrusts in crashworthiness analysis: adaptivity; mixed time integration, or subcycling, in which different timesteps are used for different parts of the mesh in explicit methods; and methods for contact-impact which are highly vectorizable. The techniques are being developed to improve the accuracy of calculations, ease-of-use of crashworthiness programs, and the speed of calculations. The latter is still of importance because crashworthiness calculations are often made with models of 20,000 to 50,000 elements using explicit time integration and require on the order of 20 to 100 hours on current supercomputers. The methodologies are briefly reviewed and then some example calculations employing these methods are described. The methods are also of value to other nonlinear transient computations.

Belytschko, Ted↗

Theoretical and Computational Studies on Sonic Boom Propagation and Its Submarine Impact

Sonic boom impact on the environment of populated area and habitat has been a major concern for the design, operation, and program planning of super/hypersonic vehicles as well as sE launch. Recent development in sonic boom studies reviewed has made evident need for amending the near-field analysis with nonlinear calculations, but an unambiguous matching procedure to assure waveform-prediction improvement is still lacking. Another problem receiving much attention recently is the renewed issue with 'transition focus booms' particularly the 'superboom' which occurs during a speed change through a threshold Mach number and gives rise to strong wave-focussing effects; however, its intensity and the extent of the impact area not be established from existing methods. A third aspect of a more recent concern is the potential sonic boom impact on pelagic and coastal environment, of which the methodology for defining impact has yet to be adequately developed. The study addresses these three aspects in the frame work of a wave-field analysis for a stratified atmosphere, employing coordinates fixed to the vehicle in steady horizontal motion.

Cheng, H. K.↗

Explicit Dynamic Impact Analysis of the Building 3525 Packages

The safe transport and handling of radioactive material containers is critical to ensure the protection of personnel, facilities, and the environment. Accidental drops during handling pose a significant hazard, necessitating robust design and analysis to mitigate potential consequences. This study presents an explicit dynamic impact analysis of two container designs, the Baby Sugarman and the Core Conduction Cooldown Test Facility (CCCTF) Tall Boy, used in Building 3525 at Oak Ridge National Laboratory. The analysis employed nonlinear finite element modeling in LS-DYNA to simulate free-fall impacts from a height of 16.5 ft across eight critical orientations. The study evaluated deformation, strain failure, and shielding integrity under these extreme conditions. Results demonstrated that both designs maintain structural integrity, with no breaches or loss of shielding under specified administrative controls. Notably, a maximum shielding loss of 10 lbs was observed during a top-corner drop, highlighting the need to limit drop heights to less than 1 ft for certain orientations. These findings underscore the effectiveness of engineering controls and design measures in preventing containment breaches and ensuring compliance with safety requirements. This work provides a comprehensive methodology for impact analysis, contributing valuable insights to the field of radioactive materials packaging and transportation safety.

Martinez, Oscar [ORNL] (ORCID:0000000181814046)↗

The use of experimental design to find the operating maximum power point of PEM fuel cells

Proton Exchange Membrane (PEM) Fuel Cells are difficult to model due to their complex nonlinear nature. In this paper, the development of a PEM Fuel Cells mathematical model based on the Design of Experiment methodology is described. The Design of Experiment provides a very efficient methodology to obtain a mathematical model for the studied multivariable system with only a few experiments. The obtained results can be used for optimization and control of the PEM Fuel Cells systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Nonlinear Ultrasonic Diagnosis and Prognosis of ASR Damage in Dry Cask Storage

Alkali-silica reaction (ASR) is a deleterious chemical process that may occur in cement-based materials such as mortars and concretes, where the hydroxyl ions in the highly alkaline pore solution attack the siloxane groups in the siliceous minerals in the aggregates. The reaction produces a cross-linked alkali-silica gel. The ASR gel swells in the presence of water. Expansion of the gel results in cracking when the swelling-induced stress exceeds the fracture toughness of the concrete. As the ASR continues, cracks may grow and eventually coalesce, which results in reduced service life and a decrease safety of concrete structures. Since concrete is widely used as a critical structural component in dry cask storage of used nuclear fuels, ASR damage poses a significant threat to the sustainability of long term dry cask storage systems. Therefore, techniques for effectively detecting, managing and mitigating ASR damage are needed. Currently, there are no nondestructive methods to accurately detect ASR damage in existing concrete structures. The only current way of accurately assessing ASR damage is to drill a core from an existing structure, and conduct microscopy on this drilled cylindrical core. Clearly, such a practice is not applicable to dry cask storage systems. To meet these needs, this research is aimed at developing (1) a suite of nonlinear ultrasonic quantitative nondestructive evaluation (QNDE) techniques to characterize ASR damage, and (2) a physics-based model for ASR damage evolution using the QNDE data. Outcomes of this research will provide a nondestructive diagnostic tool to evaluate the extent of the ASR damage, and a prognostic tool to estimate the future reliability and safety of the concrete structures in dry cask storage systems

42 ENGINEERING↗

Analysis of simultaneous aerosol and ocean glint retrieval using multi-angle observations

Since early 2000, NASA's Multi-angle Imaging SpectroRadiometer (MISR) instrument has been performing remote sensing retrievals of aerosol optical properties from the polar-orbiting Terra spacecraft. A noteworthy aspect of MISR observations over the ocean is that, for much of the Earth, some of the multi-angle views have contributions from solar reflection by the ocean surface (glint, or glitter), while others do not. Aerosol retrieval algorithms often discard these glint-influenced observations because they can overwhelm the signal and are difficult to predict without knowledge of the (wind-speed-driven) ocean surface roughness. However, theoretical studies have shown that multi-angle observations of a location at geometries with and without reflected sun glint can be a rich source of information, sufficient to support simultaneous retrieval of both the aerosol state and the wind speed at the ocean surface. We are in the early stages of creating such an algorithm. In this paper, we describe our assessment of the appropriate level of parameterization for simultaneous aerosol and ocean surface property retrievals using sun glint. For this purpose, we use generalized nonlinear retrieval analysis (GENRA), an information content assessment (ICA) technique employing Bayesian inference, and simulations from the Ahmad–Fraser iterative radiative transfer code. We find that four parameters are suitable: aerosol optical depth (τ), particle size distribution (expressed as the fine mode fraction f of small particles in a bimodal size distribution), surface wind speed (w), and relative humidity (r, to define the aerosol water content and complex refractive index). None of these parameters define ocean optical properties, as we found that the aerosol state could be retrieved with the nine MISR near-infrared views alone, where the ocean body is strongly absorbing in the open ocean. We also found that retrieval capability varies with observation geometry and that as τ increases so does the ability to determine aerosol intensive optical properties (r and f, while it decreases for w). Increases in w decrease the ability to determine the true value of that parameter but have minimal impact on retrieval of aerosol properties. We explored the benefit of excluding the two most extreme MISR view angles for which radiative transfer with the plane-parallel approximation is less certain, but we found no advantage in doing so. Finally, the impact of treating wind speed as a scalar parameter, rather than as a two-parameter directional wind, was tested. While the simpler scalar model does contribute to overall aerosol uncertainty, it is not sufficiently large to justify the addition of another dimension to parameter space. An algorithm designed upon these principles is in development. It will be used to perform an atmospheric correction with MISR for coincident ocean color (OC) observations by the Moderate Resolution Imaging Spectroradiometer (MODIS) instrument, also on the NASA Terra spacecraft. Unlike MISR, MODIS is a single-view-angle instrument, but it has a more complete set of spectral channels ideal for determination of optical ocean properties. The atmospheric correction of MODIS OC data can therefore benefit from MISR aerosol retrievals. Furthermore, higher-spatial-resolution data from coincident MISR observations may also improve glint screening.

Kirk Knobelspiesse↗

Fluid moments of the nonlinear Landau collision operator

An important problem in plasma physics is the lack of an accurate and complete description of Coulomb collisions in associated fluid models. To shed light on the problem, this Letter introduces an integral identity involving the multivariate Hermite tensor polynomials and presents a method for computing exact expressions for the fluid moments of the nonlinear Landau collision operator. The proposed methodology provides a systematic and rigorous means of extending the validity of fluid models that have an underlying inverse-square force particle dynamics to arbitrary collisionality and flow.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Time-dependent oral absorption models

The plasma concentration-time profiles following oral administration of drugs are often irregular and cannot be interpreted easily with conventional models based on first- or zero-order absorption kinetics and lag time. Six new models were developed using a time-dependent absorption rate coefficient, ka(t), wherein the time dependency was varied to account for the dynamic processes such as changes in fluid absorption or secretion, in absorption surface area, and in motility with time, in the gastrointestinal tract. In the present study, the plasma concentration profiles of propranolol obtained in human subjects following oral dosing were analyzed using the newly derived models based on mass balance and compared with the conventional models. Nonlinear regression analysis indicated that the conventional compartment model including lag time (CLAG model) could not predict the rapid initial increase in plasma concentration after dosing and the predicted Cmax values were much lower than that observed. On the other hand, all models with the time-dependent absorption rate coefficient, ka(t), were superior to the CLAG model in predicting plasma concentration profiles. Based on Akaike's Information Criterion (AIC), the fluid absorption model without lag time (FA model) exhibited the best overall fit to the data. The two-phase model including lag time, TPLAG model was also found to be a good model judging from the values of sum of squares. This model also described the irregular profiles of plasma concentration with time and frequently predicted Cmax values satisfactorily. A comparison of the absorption rate profiles also suggested that the TPLAG model is better at prediction of irregular absorption kinetics than the FA model. In conclusion, the incorporation of a time-dependent absorption rate coefficient ka(t) allows the prediction of nonlinear absorption characteristics in a more reliable manner.

Non-NASA Center↗

Solving numerically nonlinear systems of balance laws by multivariate sigmoidal functions approximation

A numerical method for solving systems of nonlinear one-dimensional balance laws, based on multivariate sigmoidal functions approximation, is developed. Constructive approximation theorems are first established in both, uniform and L{sup p} norms. A priori as well as a posteriori error estimates are derived for the numerical solutions, when various kinds of sigmoidal functions, such as unit step, logistic and Gompertz functions, are chosen. The residual of the numerical method is also estimated. Numerical examples are given to test the performance of the algorithm, a comparison with the Godunov method is made concerning accuracy and computational cost. Finally, the numerical stability of the method is analyzed.

97 MATHEMATICS AND COMPUTING↗

Demonstration of acoustic monitoring for structural health of microreactors: Through use of neural networks and resonant ultrasound spectroscopy

Nuclear microreactors prioritize modularity and portability and are intended to be a cost-effective technology for non-conventional nuclear markets. As such, the development of microreactors into a safe and feasible solution for energy security applications will necessitate the development of non-destructive technologies to monitor the integrity of inaccessible reactors components during operation. This demonstration applies linear and nonlinear acoustic techniques, in combination with machine learning, to detect and classify mechanical changes (stress and damage) in a test article which are broadly representative of potential operating challenges within a functioning microreactor. All necessary data has been collected for this demonstration, with minor experimental issues identified that can be addressed in follow-on work. Motivated by the expected conditions within a functioning microreactor, we have demonstrated our monitoring techniques on a core-block-like test article using an unstructured excitation source that approximates the noisy acoustic environment expected during reactor operation. At all stress states mechanically applied to the test article, a machine learning model using an artificial neural network was able to classify with 100% accuracy whether a 3D laser vibrometry point measurement was made on an intact or artificially defective test article. Further, model predictions about whether the defect interface was rough or smooth were 95% accurate, indicating the ability of acoustic techniques to recover defect characteristics. Resonant ultrasound spectroscopy (RUS) was also applied to the dataset to provide further quantitative insights about material properties. RUS analysis was ultimately hampered by several minor experimental and data issues, limiting results to certain cases for this demonstration. Last, analysis using nonlinear RUS exhibited sensitivity to changing levels of applied stresses for each intact and defective state. As presented in this demonstration, acoustic monitoring exhibits sensitivity to stress changes, which are of concern due to high thermal gradients expected during startup and operation. Further, our techniques distinguish between measurements made on intact and damaged test articles.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

An adaptive space-time phase field formulation for dynamic fracture of brittle shells based on LR NURBS

Here, we present an adaptive space-time phase field formulation for dynamic fracture of brittle shells. Their deformation is characterized by the Kirchhoff–Love thin shell theory using a curvilinear surface description. All kinematical objects are defined on the shell’s mid-plane. The evolution equation for the phase field is determined by the minimization of an energy functional based on Griffith’s theory of brittle fracture. Membrane and bending contributions to the fracture process are modeled separately and a thickness integration is established for the latter. The coupled system consists of two nonlinear fourth-order PDEs and all quantities are defined on an evolving two-dimensional manifold. Since the weak form requires C 1 -continuity, isogeometric shape functions are used. The mesh is adaptively refined based on the phase field using Locally Refinable (LR) NURBS. Time is discretized based on a generalized-α method using adaptive time-stepping, and the discretized coupled system is solved with a monolithic Newton–Raphson scheme. Finally, the interaction between surface deformation and crack evolution is demonstrated by several numerical examples showing dynamic crack propagation and branching.

42 ENGINEERING↗

Analysis for delamination initiation in postbuckled dropped-ply laminates

The compression strength of dropped-ply, graphite-epoxy laminated plates for the delamination mode of failure is studied by analysis and corroborated with experiments. The nonlinear response of the test specimens is modeled by a geometrically nonlinear finite element analysis. The methodology for predicting delamination is based on a quadratic interlaminar stress criterion evaluated at a characteristic distance from the ply drop-off. The compression strength of specimens exhibiting a linear response is greater than the compression strength of specimens with the same layup exhibiting a geometrically nonlinear response. The analyses for both linear and nonlinear response show that severe interlaminar stress gradients occur in the interfaces at the drop-off because of the thickness/stiffness discontinuity. However, these interlaminar stress distributions are altered in the geometrically nonlinear response such that, with increasing load, their growth at the center of the laminate is retarded while their growth near the unloaded supported edge is increased.

Davila, Carlos G.↗

On the estimation of the attainability set of nonlinear control systems

Analysis of the attainability set and construction of its estimates greatly facilitates the solution of a variety of problems in mathematical control theory. In the paper, the problem of boundedness of the integral funnel of nonlinear controlled system is considered. Some estimates of the attainability sets for a nonlinear controlled system are presented. Theorems on the boundedness of considered integral funnels are proved.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A multiresolution adaptive wavelet method for nonlinear partial differential equations

We report the multiscale complexity of modern problems in computational science and engineering can prohibit the use of traditional numerical methods in multi-dimensional simulations. Therefore, novel algorithms are required in these situations to solve partial differential equations (PDEs) with features evolving on a wide range of spatial and temporal scales. To meet these challenges, we present a multiresolution wavelet algorithm to solve PDEs with significant data compression and explicit error control. We discretize in space by projecting fields and spatial derivative operators onto wavelet basis functions. We provide error estimates for the wavelet representation of fields and their derivatives. Then, our estimates are used to construct a sparse multiresolution discretization which guarantees the prescribed accuracy. Additionally, we embed a predictor-corrector procedure within the temporal integration to dynamically adapt the computational grid and maintain the accuracy of the solution of the PDE as it evolves. We present examples to highlight the accuracy and adaptivity of our approach.

97 MATHEMATICS AND COMPUTING↗

The quiescent H-mode regime for high performance edge localized mode-stable operation in future burning plasmas

For the first time, DIII-D experiments have achieved stationary quiescent H-mode (QH-mode) operation for many energy confinement times at simultaneous ITER-relevant values of beta, confinement, and safety factor, in an ITER-like shape. QH-mode provides excellent energy confinement, even at very low plasma rotation, while operating without edge localized modes (ELMs) and with strong impurity transport via the benign edge harmonic oscillation (EHO). By tailoring the plasma shape to improve the edge stability, the QH-mode operating space has also been extended to densities exceeding 80% of the Greenwald limit, overcoming the long-standing low-density limit of QH-mode operation. In the theory, the density range over which the plasma encounters the kink-peeling boundary widens as the plasma cross-section shaping is increased, thus increasing the QH-mode density threshold. The DIII-D results are in excellent agreement with these predictions, and nonlinear magnetohydrodynamic analysis of reconstructed QH-mode equilibria shows unstable low n kink-peeling modes growing to a saturated level, consistent with the theoretical picture of the EHO. Furthermore, high density operation in the QH-mode regime has opened a path to a new, previously predicted region of parameter space, named “Super H-mode” because it is characterized by very high pedestals that can be more than a factor of two above the peeling-ballooning stability limit for similar ELMing H-mode discharges at the same density.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Analysis of Triangular Rollable and Collapsible (TRAC) Composite Booms under the Effects of Gravity and Twist

Finite element analysis is used to investigate Triangular Rollable and Collapsible (TRAC) composite booms. The TRAC booms are subjected to a suite of loading conditions, including in-plane bending, out-of-plane bending, and axial compression. Simulations containing as-designed boom geometry are performed, as are simulations with geometric deviations associated with the as-built conditions of a sample TRAC boom. Namely, a cross-sectional twist along the length of the boom is incorporated into the analysis to study the effect of the twist on the deformation characteristics of the boom. Twisting in the sample boom was observed after short-term storage in a rolled configuration and may be attributed to differential creep deformation between the inner and outer flanges of the TRAC boom. Simulations are conducted with and without the effects of gravity to understand the significance of gravity on testing conditions. Booms of three lengths are evaluated: 3 meters, 7.5 meters and 30 meters. The mechanical responses of the TRAC booms are discussed with reference to the effects of gravity, boom length and cross-sectional twisting. Gravity is found to have a more significant effect on boom mechanical response as boom length is increased. Cross-sectional twisting is found to have a more significant effect on boom mechanical response in the shorter booms.

High-strain composites↗

Analysis of Triangular Rollable and Collapsible Composite Booms Under the Effects of Gravity and Twist

Finite element analysis is used to investigate Triangular Rollable and Collapsible (TRAC) composite booms. The TRAC booms are subjected to a suite of loading conditions, including in-plane bending, out-of-plane bending, and axial compression. Simulations containing as-designed boom geometry are performed, as are simulations with geometric deviations associated with the as-built conditions of a sample TRAC boom. Namely, a cross-sectional twist along the length of the boom is incorporated into the analysis to study the effect of the twist on the deformation characteristics of the boom. Twisting in the sample boom was observed after short-term storage in a rolled configuration and may be attributed to differential creep deformation between the inner and outer flanges of the TRAC boom. Simulations are conducted with and without the effects of gravity to understand the significance of gravity on testing conditions. Booms of three lengths are evaluated: 3 meters, 7.5 meters and 30 meters. The mechanical responses of the TRAC booms are discussed with reference to the effects of gravity, boom length and cross-sectional twisting. Gravity is found to have a more significant effect on boom mechanical response as boom length is increased. Cross-sectional twisting is found to have a more significant effect on boom mechanical response in the shorter booms.

High-strain composites↗

Pyomo.DOE: An open-source package for model-based design of experiments in Python

Predictive mathematical models are a cornerstone of science and engineering. Yet selecting, calibrating, and validating said science-based models often remains an art in practice. Model-based design of experiments (MBDoE) provides a systematic framework to maximize information gain from experiments while minimizing time and resource costs. But MBDoE remains limited to niche application areas, in part because practitioners must integrate expertise in statistics, computational optimization, and modeling. To help reduce this barrier, we introduce Pyomo.DOE, an open-source package for MBDoE. Pyomo.DOE uses a nonlinear sensitivity analysis code k_aug to quickly approximate the Fisher information matrix and leverages a new stochastic programming abstraction. We demonstrate Pyomo.DOE with the first application of MBDoE to fixed-bed breakthrough experiments, which highlights the power of Pyomo.DOE to quantify the value of experimental modifications a priori for large-scale partial differential-algebraic equation (PDAE) models. Here we also provide a mathematical primer on MBDoE targeted at general chemical engineers.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗