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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,423 records · Page 79

All sources NASA NTRS DOE OSTI

NASA NTRS2004

Invariant Theory for Dispersed Transverse Isotropy: An Efficient Means for Modeling Fiber Splay

Most soft tissues possess an oriented architecture of collagen fiber bundles, conferring both anisotropy and nonlinearity to their elastic behavior. Transverse isotropy has often been assumed for a subset of these tissues that have a single macroscopically-identifiable preferred fiber direction. Micro-structural studies, however, suggest that, in some tissues, collagen fibers are approximately normally distributed about a mean preferred fiber direction. Structural constitutive equations that account for this dispersion of fibers have been shown to capture the mechanical complexity of these tissues quite well. Such descriptions, however, are computationally cumbersome for two-dimensional (2D) fiber distributions, let alone for fully three-dimensional (3D) fiber populations. In this paper, we develop a new constitutive law for such tissues, based on a novel invariant theory for dispersed transverse isotropy. The invariant theory is based on a novel closed-form splay invariant that can easily handle 3D fiber populations, and that only requires a single parameter in the 2D case. The model is polyconvex and fits biaxial data for aortic valve tissue as accurately as the standard structural model. Modification of the fiber stress-strain law requires no re-formulation of the constitutive tangent matrix, making the model flexible for different types of soft tissues. Most importantly, the model is computationally expedient in a finite-element analysis.

Freed, alan D.↗

DOE OSTI2020

Sparse and Random Sampling Techniques for High-Resolution, Full-Field, BSS-Based Structural Dynamics Identification from Video

Video-based techniques for identification of structural dynamics have the advantage that they are very inexpensive to deploy compared to conventional accelerometer or strain gauge techniques. When structural dynamics from video is accomplished using full-field, high-resolution analysis techniques utilizing algorithms on the pixel time series such as principal components analysis and solutions to blind source separation the added benefit of high-resolution, full-field modal identification is achieved. An important property of video of vibrating structures is that it is particularly sparse. Typically video of vibrating structures has a dimensionality consisting of many thousands or even millions of pixels and hundreds to thousands of frames. However the motion of the vibrating structure can be described using only a few mode shapes and their associated time series. As a result, emerging techniques for sparse and random sampling such as compressive sensing should be applicable to performing modal identification on video. This work presents how full-field, high-resolution, structural dynamics identification frameworks can be coupled with compressive sampling. The techniques described in this work are demonstrated to be able to recover mode shapes from experimental video of vibrating structures when 70% to 90% of the frames from a video captured in the conventional manner are removed.

47 OTHER INSTRUMENTATION↗

DOE OSTI2025

Analysis of the linear and nonlinear stability of Alfven eigenmodes and fish-bones in JET DT discharges: mode identification and shear flows generation

The plasma in future nuclear fusion reactors will be heated by neutral beam injectors (NBIs) and high frequency electromagnetic waves as well as fusion born alpha particles. Energetic particles (EPs), with energies up to two orders of magnitude larger than the thermal plasma, can trigger EP driven modes and induce harmful EP losses, reducing the plasma heating efficiency and the economical viability of the reactor. The present study is dedicated to analyze the Alfven Eigenmode (AE) activity in JET D–T discharges, the closest experiment to reactor-like operation performed until now. There, EP driven modes are induced by the combined effect of tangential NBIs and ion cyclotron resonance heating (ICRH) driven EP. Linear and nonlinear simulations are performed with the gyro-fluid FAR3d code to analyze the AE activity observed in the discharge 99896. The linear simulations reproduce the unstable n = 3 to 5 toroidal AEs (TAE) at the inner plasma region observed in the experiment, triggered by highly energetic passing deuterium populations injected by the tangential NBIs, further accelerated by the effect of the ICRH up to 1 MeV. In addition, fish-bones triggered by energetic trapped hydrogen induced by the ICRH are also reproduced. On the other hand, the alpha particles density is too small to destabilize AEs in the experiment. Nonetheless, increasing artificially the alpha density by one order of magnitude, an n = 1 beta induced AE can be destabilized in the inner plasma region. Nonlinear simulations indicate the generation of zonal structures during the AE/fish-bone saturation phase. TAE and fish-bones causes a rather weak increase of the passing D and trapped H EP (around 2%), respectively. Shear flows and zonal currents are generated during the saturation of TAE and fish-bones. Nonlinear simulations performed for D–T and pure deuterium thermal plasma indicate AE/fish-bone activity is weaker and shear flows are less intense in the pure deuterium case, trends consistent with the experimental observations that also indicates a deterioration of the thermal plasma confinement. Therefore, both numerical studies and experimental evidence indicate the generation of shear flows by AE/fish-bones could be connected with an improvement of the thermal plasma confinement.

AE↗

NASA NTRS1996

A General Interface Method for Aeroelastic Analysis of Aircraft

The aeroelastic analysis of an aircraft requires an accurate and efficient procedure to couple aerodynamics and structures. The procedure needs an interface method to bridge the gap between the aerodynamic and structural models in order to transform loads and displacements. Such an interface method is described in this report. This interface method transforms loads computed by any aerodynamic code to a structural finite element (FE) model and converts the displacements from the FE model to the aerodynamic model. The approach is based on FE technology in which virtual work is employed to transform the aerodynamic pressures into FE nodal forces. The displacements at the FE nodes are then converted back to aerodynamic grid points on the aircraft surface through the reciprocal theorem in structural engineering. The method allows both high and crude fidelities of both models and does not require an intermediate modeling. In addition, the method performs the conversion of loads and displacements directly between individual aerodynamic grid point and its corresponding structural finite element and, hence, is very efficient for large aircraft models. This report also describes the application of this aero-structure interface method to a simple wing and an MD-90 wing. The results show that the aeroelastic effect is very important. For the simple wing, both linear and nonlinear approaches are used. In the linear approach, the deformation of the structural model is considered small, and the loads from the deformed aerodynamic model are applied to the original geometry of the structure. In the nonlinear approach, the geometry of the structure and its stiffness matrix are updated in every iteration and the increments of loads from the previous iteration are applied to the new structural geometry in order to compute the displacement increments. Additional studies to apply the aero-structure interaction procedure to more complicated geometry will be conducted in the second phase of the present contract.

Tzong, T.↗

DOE OSTI2026

Effects of external magnetic fields on the multi-mode Rayleigh–Taylor instability

The magneto-Rayleigh–Taylor instability (mRTI) is a key process in inertial confinement fusion and is thought to be widespread in the interstellar medium, where it can concentrate plasma into discrete structures. We present resistive MHD simulations of the nonlinear evolution of multi-mode mRTI in both two and three dimensions, examining the effects of uniform external magnetic fields oriented either parallel or perpendicular to the initial interface. In both 2D and 3D, weak parallel fields enhance mixing-zone growth, whereas stronger fields suppress it. For perpendicular fields, growth is initially inhibited but becomes enhanced at later times. These behaviors arise from magnetic tension, which modifies flow anisotropy, buoyancy, drag, and vortex dynamics. The interplay of these mechanisms governs the distinct ways in which magnetic fields influence mRTI evolution.

Flow instabilities↗

DOE OSTI2022

Anisotropic yield models for lattice unit cell structures exploiting orthotropic symmetry

We report numerical homogenization enables efficient computational analysis and design of multiscale structures made of micro-architected materials including lattice unit cells. However, predicting yield is nontrivial because it requires accurate and efficient predictions of the maximum stress inside the homogenized unit cells. To address this challenge, we develop a macroscale anisotropic yield function for micro-architected materials. The yield function depends on the three-dimensional macroscale stress state and the parameters describing a family of micro-architectures, such as the radii of the struts in a lattice unit cell. To ensure accuracy, we determine the maximum stress using three-dimensional continuum finite-element analysis. To ensure efficiency, we construct surrogate models from the aforementioned high-fidelity results for yield prediction. To reduce simulation costs and surrogate modeling complexity, we leverage orthotropic symmetry commonly found in lattice unit cells. In this paper, we provide a thorough presentation of group representations and the systematic procedure to exploit orthotropic domain symmetry in homogenization and surrogate modeling. We illustrate the use of linear homogenization results to predict yield in specific unit cells without further simulations. We furthermore show that surrogate modeling presents a viable option for anisotropic yield prediction in a continuously parametrized family of micro-architectures. More specifically, despite the dimensionality and degree of nonlinearity in the maximum micro von Mises stress within the unit cell, the surrogate models can predict it with less than 5% error at least 90% of the time. Moreover, the largest under-prediction error, which is more critical than the over-prediction error, is typically less than 10%.

42 ENGINEERING↗

NASA NTRS1992

Model of a rotating magnetic cloud

The possibility that magnetic clouds rotate while they propagate antisunward was investigated. Magnetic clouds are modeled as magnetic flux ropes which rotate rigidly about the axis of symmetry. An ideal magnetohydrodynamic model, in which the evolution of the magnetic structure is related to the time evolution of the angular frequency, is developed. A class of 'separable' magnetic fields is employed to reduce the problem to a nonlinear ordinary differential equation for the evolution function, and it is solved numerically. The corresponding effective potential gives rise to two modes of evolution--expansion and oscillation--depending on the energy and on the value of a dimensionless parameter, k. Parameter k depends on the gas pressure, the ratio of the magnetic field components, and the frequency of rotation. There is a critical value of k, k(sub c), above which the oscillatory regime disappears and the flux rope invariably expands, regardless of the energy. Below k(sub c) the energy determines whether the configuration is confined or unbounded. Rotation always helps expansion by lowering the potential barrier. A data example was studied and features which are interpreted as signatures of rotation are presented. The angular speed is comparable to the Alfven speed, and the core of the rotating cloud completes on average one full revolution every three days at 1 AU. The parameter k is calculated from observations, and it is found to be close to, but below, critical. Only three out of the nine clouds examined showed signatures of rotation. Theoretical analysis suggests that close to the Sun rotation effects may play a more important role in the evolution of magnetic clouds than 1 AU.

Farrugia, C. J.↗

DOE OSTI2021

The impact of alkali‐ion intercalation on redox chemistry and mechanical deformations: Case study on intercalation of Li, Na, and K ions into FePO 4 cathode

Abstract Batteries made of charge carriers from Earth‐crust abundant materials (e.g., Na, K, and Mg) have received extensive attention as an alternative to Li‐ion batteries for grid storage. However, a lack of understanding of the behavior of these larger ions in the electrode materials hinders the development of electrode structures suitable for these large ions. In this study, we investigate the impact of alkali ions (Li, Na, and K) on the redox chemistry and mechanical deformations of iron phosphate composite cathodes by using electrochemical techniques and in situ digital image correlation. Na‐ion and Li‐ion intercalation demonstrate a nearly linear correlation between electrochemical strains and the state of charge and discharge. The strain development shows nonlinear dependance on the state of charge and discharge for K ions. Strain rate calculations show that K ion intercalation results in a progressive increase in the strain rate for all cycles. Li and Na intercalation induce nearly constant strain rates with the exception of the first discharge cycle of Na intercalation. When the same amount of ions are inserted into the electrode, the electrode shows the lowest strain generation upon Li intercalation compared to larger alkali ions. Na and K ions induce similar volumetric changes in the electrode when the state of charge and discharge is around 30%. Although the electrode experiences larger absolute strain generation at the end of the discharge cycles upon Na intercalation, strain rates were found to be greater for K ions. Potential‐dependent behaviors also demonstrate more sluggish redox reactions during K intercalation, compared to Li and Na. Our quantitative analysis suggests that the strain rate, rather than the absolute value of strain, is the critical factor in amorphization of the crystalline electrode.

25 ENERGY STORAGE↗

NASA NTRS2020

Feasibility Study of SDAS Instrumentation's Ability to Identify Mobile Launcher (ML)/Crawler-Transporter (CT) Modes During Rollout Operations

The Space Launch System (SLS) and its Mobile Launcher (ML) will be transported to the launch pad via the Crawler-Transporter (CT) system. Rollout (i.e., transportation) loads produce structural loads on the integrated SLS/Orion Multi-Purpose Crew Vehicle (MPCV) launch vehicle which are of a concern with respect to fatigue. As part of the risk reduction process and in addition to the modal building block test approach that has been adopted by the SLS Program, acceleration data will be obtained during rollout for use in modal parameter estimation. There are several occurrences where the ML/CT will be transported either into the Vertical Assembly Building (VAB) or to the launch pad and back without the SLS stack as part of the Kennedy Space Center (KSC) Exploration Ground Systems (EGS) Integrated Test and Checkout (ITCO). NASA KSC EGS has instrumentation installed on both the ML and CT to record data during rollout, at the launch pad, and during liftoff. The EGS instrumentation on the ML, which includes accelerometers, is referred to as the Sensor Data Acquisition System (SDAS). The EGS instrumentation on the CT, which also includes accelerometers, is referred to as the CT Data Acquisition System (CTDAS). The forces and accelerations applied to the ML and CT during a rollout event will be higher than any of the planned building block modal tests. This can be very beneficial in helping identify nonlinear behavior in the structure. Developing modal parameters from the same test hardware in multiple boundary conditions and under multiple levels of excitation is a key step in developing a well correlated FEM. The purpose of this study was three fold. First, determine the target modes of the ML/CT in its rollout configuration. Second, determine if the test degrees of freedom (DOF) corresponding to the layout of the SDAS/CTDAS accelerometers (i.e. position and orientation) is sufficient to identify the target modes. Third, determine if the Generic Rollout Forcing Functions (GRFF's) is sufficient for identifying the ML/CT target modes accounting for variations in CT speed, modal damping, and sensor/ambient background noise levels. The finding from the first part of this study identified 28 target modes of the ML/CT rollout configuration based upon Modal Effective Mass Fractions (MEFF) and engineering judgement. The finding from the second part of this study showed that the SDAS/CTDAS accelerometers (i.e. position and orientation) are able to identify a sufficient number of the target modes to support model correlation of the ML/CT FEM. The finding from the third part of this study confirms the GRFFs sufficiently excite the ML/CT such that varying quantities of the defined target modes should be able to be extracted when utilizing an Experimental Modal Analysis (EMA) Multi-Input Multi-Output (MIMO) analysis approach. An EMA analysis approach was used because Operational Modal Analysis (OMA) tools were not available and the GRFFs were sufficiently uncorrelated. Two key findings from this third part of the study are that the CT speed does not show a significant impact on the ability to extract the modal parameters and that keeping the ambient background noise observed at each accelerometer location at or below 30 µgrms is essential to the success of this approach.

Winkel, James P.↗

NASA NTRS1996

Analysis and Tests of Reinforced Carbon-Epoxy/Foam-Core Sandwich Panels with Cutouts

The results of a study of a low-cost structurally efficient minimum-gage shear-panel design that can be used in light helicopters are presented. The shear-panel design is based on an integrally stiffened syntactic-foam stabilized-skin with an all-bias-ply tape construction for stabilized-skin concept with an all-bias-ply tape construction for the skins. This sandwich concept is an economical way to increase the panel bending stiffness weight penalty. The panels considered in the study were designed to be buckling resistant up to 100 lbs/in. of shear load and to have an ultimate strength of 300 lbs/in. The panel concept uses unidirectional carbon-epoxy tape on a syntactic adhesive as a stiffener that is co-cured with the skin and is an effective concept for improving panel buckling strength. The panel concept also uses pultruded carbon-epoxy rods embedded in a syntactic adhesive and over-wrapped with a bias-ply carbon-epoxy tape to form a reinforcing beam which is an effective method for redistributing load around rectangular cutout. The buckling strength of the reinforced panels is 83 to 90 percent of the predicted buckling strength based on a linear buckling analysis. The maximum experimental deflection exceeds the maximum deflection predicted by a nonlinear analysis by approximately one panel thickness. The failure strength of the reinforced panels was two and a half to seven times of the buckling strength. This efficient shear-panel design concept exceeds the required ultimate strength requirement of 300 lbs/in by more than 100 percent.

Baker, Donald J.↗

DOE OSTI2021

Quasiclassical Computations of Compton-Scattered Spectra

Quality X-ray sources are crucial to fundamental physics research, medical radiology, humanities research, and materials science. While synchrotron radiation (SR) facilities produce the state-of-the-art emissions with respect to brilliance and frequency tunability, the great expense required to build, maintain, and operate these structures greatly limits their accessibility to researchers. Much of the research conducted at SR facilities, however, may be conducted with inverse Compton sources (ICS). Accelerator-based Compton scattering light sources generate high-energy, high-brilliance emissions. Compton scattering is the process by which a photon scatters o? an electron. ICS offer an affordable, in-lab alternative to SR facilities. Even though SR facilities produce greater intensity emissions, Compton sources provide the same frequency tunability at the intensities suitable for the purposes of many researchers currently fighting for time at SR facilities, i.e., ICS provides intensities suitable for contrast imaging, X-ray fluorescence, X-ray-diffraction, and X-ray spectroscopy. The focus of this work is to create computational models to simulate Compton-scattered spectra. These models have been used to build a theoretical basis for methods of improving the quality of future Compton sources and to preform diagnostic analysis of existing light sources. The theoretical basis of each model is derived from first principles. The numerical methods employed by each model are defined. A full description of the various functionalities of each code will be addressed. Furthermore, an in-depth analysis of spectral bandwidth sources is discussed. The complex physics arising from an extremely high-intensity, nonlinear laser pulse is explored in detail. Methods of frequency modulation of the incident laser, i.e., a method for correcting the nonlinear broadening effects on the scattered spectrum, will also be discussed. The work will conclude with a an exploration of the ongoing research efforts regarding regimes of operation outside of the limits of these current models.

Johnson, Erik↗

DOE OSTI2025

Mechanical Roles of Polysaccharide Assembly and Interactions in Plant Cell Walls

Plants synthesize polysaccharide-based primary cell walls that possess unique microstructures and mechanical properties to accommodate plant growth and provide protection. Here, it remains challenging to assess the role of polysaccharide organization and interactions in the mechanical behavior of primary cell walls owing to their complex microstructure and highly nonlinear mechanical responses. Employing a coarse-grained molecular dynamics model developed for onion epidermal walls, this work explores the conditions under which polysaccharide assembly and interactions might play a significant role in primary cell wall mechanics. Cellulose–cellulose adhesion plays a dominant role in the wall load-bearing capacity, but when cellulose–cellulose adhesion was disrupted computationally, cellulose–xyloglucan adhesion could influence the wall load-bearing capacity. Contrary to the common concept that xyloglucans mechanically tether well-separated cellulose microfibrils, xyloglucans functioned in this case as interfibrillar adhesives capable of transmitting tensile forces between cellulose microfibrils. Our findings may inform design criteria of new materials inspired by plant cell walls.

59 BASIC BIOLOGICAL SCIENCES↗

NASA NTRS2007

Strain Measurement Using FBG on COPV in Stress Rupture Test

White Sands Test Facility (WSTF) was requested to perform ambient temperature hydrostatic pressurization testing of a Space Transportation System (STS) 40-in. Kevlar Composite Overwrapped Pressure Vessel (COPV). The 40-in. vessel was of the same design and approximate age as the STS Main Propulsion System (MPS) and Orbiter Maneuvering System (OMS) vessels. The NASA Engineering Safety Center (NESC) assembled a team of experts and conducted an assessment that involved a review of national Kevlar COPY data. During the review, the STS COPVs were found to be beyond their original certification of ten years. The team observed that the likelihood of STS COPV Stress rupture, a catastrophic burst before leak failure mode, was greater than previously believed. Consequently, a detailed assessment of remaining stress rupture life became necessary. Prior to STS-114, a certification deviation was written for two flights of OV-103 (Discovery) and OV-104 (Atlantis) per rationale that was based on an extensive review of the Lawrence Livermore National Laboratories, COPV data, and revisions to the STS COPV stress levels. In order to obtain flight rationale to extend the certification deviation through the end of the Program, the Orbiter Project Office has directed an interagency COPV team to conduct further testing and analysis to investigate conservatism in the stress rupture model and evaluate material age degradation. Additional analysis of stress rupture life requires understanding the fiber stresses including stress that occurs due to thru-wall composite compression in COPV components. Data must be obtained at both zero gauge pressure (pre-stress) and at the component operating pressure so that this phenomenon can be properly evaluated. The zero gauge pressure stresses are predominantly a result of the autofrettage process used during vessel manufacture. Determining these pre-stresses and the constitutive behavior of the overwrap at pressure will provide necessary information to better predict the remaining life of the STS COPVs. The primary test objective is obtaining data to verify the hypothesis of a radially oriented thru-thickness stress-riser in the COPV composite whose magnitude is a function of the applied pressure and the load history. The anticipated load dependent response follows from the constitutive behavior of the composite overwrap so data to quantify its nonlinear and time dependent response will be sought. The objective of the Fiber Braggs Gratings (FBGs) were to advance the state-of-the-art by developing techniques using FBG sensors that are capable of assessing stress-rupture degradation in Kevlar COPVs in a health monitoring mode (1). Moreover, they sought to answer questions of how embedded sensors affect overall integrity of the structure. And lastly, they sought to provide an important link in the overall stress rupture study that will help close the loop on the COPV fabrication process. NDE inspection methods will be used from start to finish and FBG will be an integral link within the overall chain.

Banks, Curtis↗

DOE OSTI2022

Elastic Solutions to 2D Plane Strain Problems: Nonlinear Contact and Settlement Analysis for Shallow Foundations

The classical Neumann boundary value problem of an isotropic, homogeneous elastic half-plane under plane strain conditions is readdressed as the limiting case of the fully three-dimensional problem. Analytical solutions of the stress and strain tensors are obtained by taking the limit from known three-dimensional solutions. It is shown that the displacement fields for the plane strain problem are not well defined. A small number of simple expressions are developed, which provide a general solution for linearly-varying traction over arbitrary regions on the boundary. A simple, efficient, and rapidly convergent algorithm is developed which uses these solutions as analytic elements and provides a solution approach to the general boundary value problem. The method is verified against known solutions for Hertzian contact between parallel cylinders. Two numerical examples are presented for the analysis of shallow foundation systems. In the first, the boundary conditions are informed by analytical elastoplastic calculations and a strain influence analysis is performed and compared with the Schmertmann method. Subsequently, empirical laboratory contact traction distributions measured by Bauer et al., in both the normal and tangential directions are employed as boundary conditions for an analysis of the underlying stress field.

42 ENGINEERING↗

NASA NTRS1993

On the dynamics of approximating schemes for dissipative nonlinear equations

Since one can rarely write down the analytical solutions to nonlinear dissipative partial differential equations (PDE's), it is important to understand whether, and in what sense, the behavior of approximating schemes to these equations reflects the true dynamics of the original equations. Further, because standard error estimates between approximations of the true solutions coming from spectral methods - finite difference or finite element schemes, for example - and the exact solutions grow exponentially in time, this analysis provides little value in understanding the infinite time behavior of a given approximating scheme. The notion of the global attractor has been useful in quantifying the infinite time behavior of dissipative PDEs, such as the Navier-Stokes equations. Loosely speaking, the global attractor is all that remains of a sufficiently large bounded set in phase space mapped infinitely forward in time under the evolution of the PDE. Though the attractor has been shown to have some nice properties - it is compact, connected, and finite dimensional, for example - it is in general quite complicated. Nevertheless, the global attractor gives a way to understand how the infinite time behavior of approximating schemes such as the ones coming from a finite difference, finite element, or spectral method relates to that of the original PDE. Indeed, one can often show that such approximations also have a global attractor. We therefore only need to understand how the structure of the attractor for the PDE behaves under approximation. This is by no means a trivial task. Several interesting results have been obtained in this direction. However, we will not go into the details. We mention here that approximations generally lose information about the system no matter how accurate they are. There are examples that show certain parts of the attractor may be lost by arbitrary small perturbations of the original equations.

Jones, Donald A.↗

DOE OSTI2024

New systems in MOOSE

The Multiphysics Object-Oriented Simulation Environment (MOOSE) serves as a common library of classes between applications developed for advanced reactor analysis, fusion device engineering, spent fuel cask analysis, geochemistry studies, among other fields. These applications drive the development of the framework to meet their needs. Systems in MOOSE group capabilities that share a common purpose and generally common code. They can be leveraged by all downstream applications, providing extensive code re-use and shared maintenance. They facilitate the discovery by new users of the classes meeting at least partially their needs, and offer the same opportunities for customization as other systems. The addition of a new system to MOOSE opens new ways of solving or discretizing nonlinear problems, of performing distributed postprocessing, and a plethora of other needs. While new systems can be introduced in downstream applications rather than at the framework level, the framework team monitors common needs across the community and often triggers their addition. Documentation, training material, development needs can be centralized, limiting duplicated work across the community. The last three years have seen a large expansion in the capabilities of MOOSE. The supporting role of the framework in the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program has created numerous feature requests to support neutronics, thermal hydraulics, computational fluid dynamics and thermo-mechanics simulations in the Griffin, SAM, Pronghorn and Bison applications respectively. Similarly, laboratory-directed research and development (LDRD) projects in additive manufacturing, high-Reynolds flow simulations, structure optimization also necessitate an expansion of the framework capabilities. This summary reports on the new systems created in MOOSE, their design, their capabilities and some of the relevant interfaces.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2022

304L Can Crush Validation Studies

Accurate prediction of ductile behavior of structural alloys up to and including failure is essential in component or system failure assessment, which is necessary for nuclear weapons alteration and life extensions programs of Sandia National Laboratories. Modeling such behavior requires computational capabilities to robustly capture strong nonlinearities (geometric and material), rate- dependent and temperature-dependent properties, and ductile failure mechanisms. This study's objective is to validate numerical simulations of a high-deformation crush of a stainless steel can. The process consists of identifying a suitable can geometry and loading conditions, conducting the laboratory testing, developing a high-quality Sierra/SM simulation, and then drawing comparisons between model and measurement to assess the fitness of the simulation in regards to material model (plasticity), finite element model construction, and failure model. Following previous material model calibration, a J 2 plasticity model with a microstructural BCJ failure model is employed to model the test specimen made of 304L stainless steel. Simulated results are verified and validated through mesh and mass-scaling convergence studies, parameter sensitivity studies, and a comparison to experimental data. The converged mesh and degree of mass-scaling are the mesh discretization with 140,372 elements, and a mass scaling with a target time increment of 1.0e-6 seconds and time step scale factor of 0.5, respectively. Results from the coupled thermal-mechanical explicit dynamic analysis are comparable to the experimental data. Simulated global force vs displacement (F/D) response predicts key points such as yield, ultimate, and kinks of the experimental F/D response. Furthermore, the final deformed shape of the can and field data predicted from the analysis are similar to that of the deformed can, as measured by 3D optical CMM scans and DIC data from the experiment.

36 MATERIALS SCIENCE↗

NASA NTRS2015

Use of Generalized Fluid System Simulation Program (GFSSP) for Teaching and Performing Senior Design Projects at the Educational Institutions

This paper describes the experience of the authors in using the Generalized Fluid System Simulation Program (GFSSP) in teaching Design of Thermal Systems class at University of Alabama in Huntsville. GFSSP is a finite volume based thermo-fluid system network analysis code, developed at NASA/Marshall Space Flight Center, and is extensively used in NASA, Department of Defense, and aerospace industries for propulsion system design, analysis, and performance evaluation. The educational version of GFSSP is freely available to all US higher education institutions. The main purpose of the paper is to illustrate the utilization of this user-friendly code for the thermal systems design and fluid engineering courses and to encourage the instructors to utilize the code for the class assignments as well as senior design projects. The need for a generalized computer program for thermofluid analysis in a flow network has been felt for a long time in aerospace industries. Designers of thermofluid systems often need to know pressures, temperatures, flow rates, concentrations, and heat transfer rates at different parts of a flow circuit for steady state or transient conditions. Such applications occur in propulsion systems for tank pressurization, internal flow analysis of rocket engine turbopumps, chilldown of cryogenic tanks and transfer lines, and many other applications of gas-liquid systems involving fluid transients and conjugate heat and mass transfer. Computer resource requirements to perform time-dependent, three-dimensional Navier-Stokes computational fluid dynamic (CFD) analysis of such systems are prohibitive and therefore are not practical. Available commercial codes are generally suitable for steady state, single-phase incompressible flow. Because of the proprietary nature of such codes, it is not possible to extend their capability to satisfy the above-mentioned needs. Therefore, the Generalized Fluid System Simulation Program (GFSSP1) has been developed at NASA Marshall Space Flight Center (MSFC) as a general fluid flow system solver capable of handling phase changes, compressibility, mixture thermodynamics and transient operations. It also includes the capability to model external body forces such as gravity and centrifugal effects in a complex flow network. The objectives of GFSSP development are: a) to develop a robust and efficient numerical algorithm to solve a system of equations describing a flow network containing phase changes, mixing, and rotation; and b) to implement the algorithm in a structured, easy-to-use computer program. The analysis of thermofluid dynamics in a complex network requires resolution of the system into fluid nodes and branches, and solid nodes and conductors as shown in Figure 1. Figure 1 shows a schematic and GFSSP flow circuit of a counter-flow heat exchanger. Hot nitrogen gas is flowing through a pipe, colder nitrogen is flowing counter to the hot stream in the annulus pipe and heat transfer occurs through metal tubes. The problem considered is to calculate flowrates and temperature distributions in both streams. GFSSP has a unique data structure, as shown in Figure 2, that allows constructing all possible arrangements of a flow network with no limit on the number of elements. The elements of a flow network are boundary nodes where pressure and temperature are specified, internal nodes where pressure and temperature are calculated, and branches where flowrates are calculated. For conjugate heat transfer problems, there are three additional elements: solid node, ambient node, and conductor. The solid and fluid nodes are connected with solid-fluid conductors. GFSSP solves the conservation equations of mass and energy, and equation of state in internal nodes to calculate pressure, temperature and resident mass. The momentum conservation equation is solved in branches to calculate flowrate. It also solves for energy conservation equations to calculate temperatures of solid nodes. The equations are coupled and nonlinear; therefore, they are solved by an iterative numerical scheme. GFSSP employs a unique numerical scheme known as simultaneous adjustment with successive substitution (SASS), which is a combination of Newton-Raphson and successive substitution methods. The mass and momentum conservation equations and the equation of state are solved by the Newton-Raphson method while the conservation of energy and species are solved by the successive substitution method. GFSSP is linked with two thermodynamic property programs, GASP2 and WASP3 and GASPAK4, that provide thermodynamic and thermophysical properties of selected fluids. Both programs cover a range of pressure and temperature that allows fluid properties to be evaluated for liquid, liquid-vapor (saturation), and vapor region. GASP and WASP provide properties of 12 fluids. GASPAK includes a library of 36 fluids. GFSSP has three major parts. The first part is the graphical user interface (GUI), visual thermofluid analyzer of systems and components (VTASC). VTASC allows users to create a flow circuit by a 'point and click' paradigm. It creates the GFSSP input file after the completion of the model building process. GFSSP's GUI provides the users a platform to build and run their models. It also allows post-processing of results. The network flow circuit is first built using three basic elements: boundary node, internal node, and branch.

Majumdar, A. K.↗