Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Finite element model”

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 991 records · Page 55

Rapid Modeling, Assembly and Simulation in Design Optimization

A new capability for design is reviewed. This capability provides for rapid assembly of detail finite element models early in the design process where costs are most effectively impacted. This creates an engineering environment which enables comprehensive analysis and design optimization early in the design process. Graphical interactive computing makes it possible for the engineer to interact with the design while performing comprehensive design studies. This rapid assembly capability is enabled by the use of Interface Technology, to couple independently created models which can be archived and made accessible to the designer. Results are presented to demonstrate the capability.

Housner, Jerry↗

A Finite-Element Approach for Modeling Inviscid and Viscous Compressible Flows using Prismatic Grids

The Galerkin finite-element method is used to solve the Euler and Navier-Stokes equations on prismatic meshes. It is shown that the prismatic grid is advantageous for correctly and efficiently capturing the boundary layers in high Reynolds number flows. It can be captured accurately because of the ability to cluster grid points normal to the body. The efficiency derives from the implicit treatment of the normal direction. To treat the normal direction implicitly, a semi-implicit Runge-Kutta time stepping scheme is developed. The semi-implicit algorithm is validated on simple geometries for inviscid and viscous flows and its convergence history is compared to that of the explicit Runge-Kutta scheme. The semi-implicit scheme is shown to be a factor of 3 to 4 faster in terms of CPU time to convergence.

Pandya, S. A.↗

Coupled Structural, Thermal, Phase-Change and Electromagnetic Analysis for Superconductors

Described are the theoretical development and computer implementation of reliable and efficient methods for the analysis of coupled mechanical problems that involve the interaction of mechanical, thermal, phase-change and electromagnetic subproblems. The focus application has been the modeling of superconductivity and associated quantum-state phase-change phenomena. In support of this objective the work has addressed the following issues: (1) development of variational principles for finite elements, (2) finite element modeling of the electromagnetic problem, (3) coupling of thermal and mechanical effects, and (4) computer implementation and solution of the superconductivity transition problem. The main accomplishments have been: (1) the development of the theory of parametrized and gauged variational principles, (2) the application of those principled to the construction of electromagnetic, thermal and mechanical finite elements, and (3) the coupling of electromagnetic finite elements with thermal and superconducting effects, and (4) the first detailed finite element simulations of bulk superconductors, in particular the Meissner effect and the nature of the normal conducting boundary layer. The theoretical development is described in two volumes. This volume, Volume 1, describes mostly formulations for specific problems. Volume 2 describes generalization of those formulations.

Felippa, C. A.↗

Interior Noise Predictions in the Preliminary Design of the Large Civil Tiltrotor (LCTR2)

A prediction scheme was established to compute sound pressure levels in the interior of a simplified cabin model of the second generation Large Civil Tiltrotor (LCTR2) during cruise conditions, while being excited by turbulent boundary layer flow over the fuselage, or by tiltrotor blade loading and thickness noise. Finite element models of the cabin structure, interior acoustic space, and acoustically absorbent (poro-elastic) materials in the fuselage were generated and combined into a coupled structural-acoustic model. Fluctuating power spectral densities were computed according to the Efimtsov turbulent boundary layer excitation model. Noise associated with the tiltrotor blades was predicted in the time domain as fluctuating surface pressures and converted to power spectral densities at the fuselage skin finite element nodes. A hybrid finite element (FE) approach was used to compute the low frequency acoustic cabin response over the frequency range 6-141 Hz with a 1 Hz bandwidth, and the Statistical Energy Analysis (SEA) approach was used to predict the interior noise for the 125-8000 Hz one-third octave bands.

Grosveld, Ferdinand W.↗

Analysis of thermal and mechanical properties with inventory level of the molten salt storage tank in central receiver concentrating solar power plants

Molten salt thermal energy storage (TES) tanks ensure steady power output of concentrating solar power (CSP) plants; however, recent tank failures have highlighted the need for further analysis. Current studies primarily focus on analyzing the molten salt flow, heat transfer, and thermal efficiency. Additionally, research on the latest tank structures is limited and lacks newest experimental validation. This study measures temperature and molten salt inventory levels in the high-temperature tank at a 50 MW central receiver CSP plant, connected to the power grid in 2019. A multi-physics model was developed to evaluate thermal and mechanical properties of TES tanks by combining computational fluid dynamics and finite element modeling using real plant data. Heat loss, temperature, displacement, and stress distribution of the tank at different inventory levels were investigated. Results show that ambient air velocity near the tank roof reaches 2.14 m/s, much higher than 0.2 m/s near the wall. The temperatures of inventory fluid and tank are close, varying slightly at different levels due to thermal conduction and radiation. Because the heat loss strongly depends on temperature, the total tank loss remains nearly constant across inventory levels. Larger temperature gradients and thermal stresses are primarily localized along the tank floor edge and the air-salt interface. Notably, the maximum thermal stress at the tank edge is three times higher than that at the interface. The magnitude of total stress changes by less than 5 MPa with and without thermal load, indicating that high temperatures exert only a minor impact on tank stress. In contrast, thermal load significantly affects tank deformation, particularly at the roof edge, where values exceed 150 mm. Despite the large variation in molten salt levels, tank wall temperatures and displacements present a minor change, suggesting a weak correlation with inventory levels. In conclusion, the findings obtained in this study provide important insights on the TES tank that could be used to optimize tank design and operation strategies.

14 SOLAR ENERGY↗

Active learning with multifidelity modeling for efficient rare event simulation

Here, while multifidelity modeling provides a cost-effective way to conduct uncertainty quantification with computationally expensive models, much greater efficiency can be achieved by adaptively deciding the number of required high-fidelity (HF) simulations, depending on the type and complexity of the problem and the desired accuracy in the results. We propose a framework for active learning with multifidelity modeling emphasizing the efficient estimation of rare events. Our framework works by fusing a low-fidelity (LF) prediction with an HF-inferred correction, filtering the corrected LF prediction to decide whether to call the high-fidelity model, and for enhanced subsequent accuracy, adapting the correction for the LF prediction after every HF model call. The framework does not make any assumptions as to the LF model type or its correlations with the HF model. In addition, for improved robustness when estimating smaller failure probabilities, we propose using dynamic active learning functions that decide when to call the HF model. We demonstrate our framework using several academic case studies (including some high-dimensional problems) and two finite element model case studies: estimating Navier-Stokes velocities using the Stokes approximation and estimating stresses in a transversely isotropic model subjected to displacements via a coarsely meshed isotropic model. Across these case studies, not only did the proposed framework estimate the failure probabilities accurately, but compared with either Monte Carlo or a standard variance reduction method, it also required only a small fraction of the calls to the HF model.

42 ENGINEERING↗

Transient/structural analysis of a combustor under explosive loads

The 8-Foot High Temperature Tunnel (HTT) at NASA Langley Research Center is a combustion-driven blow-down wind tunnel. A major potential failure mode that was considered during the combustor redesign was the possibility of a deflagration and/or detonation in the combustor. If a main burner flame-out were to occur, then unburned fuel gases could accumulate and, if reignited, an explosion could occur. An analysis has been performed to determine the safe operating limits of the combustor under transient explosive loads. The failure criteria was defined and the failure mechanisms were determined for both peak pressures and differential pressure loadings. An overview of the gas dynamics analysis was given. A finite element model was constructed to evaluate 13 transient load cases. The sensitivity of the structure to the frequency content of the transient loading was assessed. In addition, two closed form dynamic analyses were conducted to verify the finite element analysis. It was determined that the differential pressure load or thrust load was the critical load mechanism and that the nozzle is the weak link in the combustor system.

Gregory, Peyton B.↗

Improvements in Nonconservative Force Modelling for TOPEX/POSEIDON

It was recognized prior to the launch of TOPEX/POSEIDON, that the most important source of orbit error other than the gravity field, was due to nonconservative force modelling. Accordingly, an intensive effort was undertaken to study the nonconservative forces acting on the spacecraft using detailed finite element modelling (Antreasian, 1992; Antreasian and Rosborough, 1992). However, this detailed modelling was not suitable for orbit determination, and a simplified eight plate "box-wing" model was developed that took into account the aggregate effect of the various materials and associated thermal properties of each spacecraft surface. The a priori model was later tuned post launch with actual tracking data [Nerem et al., 1994; Marshall and Luthcke, 1994; Marshall et al., 1995]. More recently, Kubitschek (1997] developed a newer box-wing model for TOPEX/POSEIDON, which included updated material properties, accounted for a solar array deflection, and modelled solar array warping due to thermal effects. We have used this updated model as a basis to retune the macromodel for TOPEX/POSEIDON, and report on preliminary results using at least 36 cycles (one year) of SLR and DORIS data in 1993.

Lemoine, Frank G.↗

Modeling of Triangular Lattice Space Structures with Curved Battens

Techniques for simulating an assembly process of lattice structures with curved battens were developed. The shape of the curved battens, the tension in the diagonals, and the compression in the battens were predicted for the assembled model. To be able to perform the assembly simulation, a cable-pulley element was implemented, and geometrically nonlinear finite element analyses were performed. Three types of finite element models were created from assembled lattice structures for studying the effects of design and modeling variations on the load carrying capability. Discrepancies in the predictions from these models were discussed. The effects of diagonal constraint failure were also studied.

Chen, Tzikang↗

A Methodology for Mapping Launch Vehicle Buffet Loads

Buffet loads represent the primary source of high frequency loading for launch vehicles during the ascent portion of flight. Currently, experimental techniques establish the nature of buffeting using a rigid scale model of the vehicle. The buffet forcing functions resulting from such tests are then applied to reduced finite-element models of the full-scale vehicle to determine the response and consequent loading. This paper discusses the techniques required to translate model-derived, empirical buffet forcing functions into responses for the full-scale launch vehicle, as used to determine the buffet loading for NASA's Ares I launch vehicle.

Schwarz, Jordan B.↗

Near Earth Asteroid Scout Thrust and Torque Model

The Near Earth Asteroid (NEA) Scout is a solar sail mission whose objective is to scout at least one Near Earth Asteroid in preparation for manned missions to asteroids. NEA Scout will use a solar sail as the primary means of propulsion. Thus it is important for mission planning to accurately characterize the thrust of the sail. Additionally, the solar sail creates a relatively large solar disturbance torque that must be mitigated. For early mission design studies a flat plate model of the solar sail with a fixed center of pressure was adequate, but as mission concepts and the sail design matured, greater fidelity was required. Here we discuss the progress to a three-dimensional sail model that includes the effects of tension and thermal deformation that has been derived from a large structural Finite Element Model (FEM) developed by the Langley Research Center. We have found that the deformed sail membrane affects torque relatively much more than thrust. We have also found that other than uncertainty over the precise shape, the effect of small (approximately millimeter scale) wrinkles on the diffusivity of the sail is the leading remaining source of uncertainty. We demonstrate that millimeter-scale wrinkles can be modeled analytically as a change in the fraction of specular reflection. Finally we discuss the implications of these results for the NEA Scout mission.

Heaton, Andrew↗

Near Earth Asteroid Scout Solar Sail Thrust and Torque Model

The Near Earth Asteroid (NEA) Scout is a solar sail mission whose objective is to scout at least one Near Earth Asteroid to help prepare for human missions to Near Earth Asteroids. NEA Scout will launch as a secondary payload on the first SLS-Orion mission. NEA Scout will perform a small trim maneuver shortly after deploy from the spent SLS upper stage using a cold gas propulsion system, but from that point on will depend entirely on the solar sail for thrust. As such, it is important to accurately characterize the thrust of the sail in order to achieve mission success. Additionally, the solar sail creates a relatively large solar disturbance torque that must be mitigated. For early mission design studies a flat plate model of the solar sail with a fixed center of pressure was adequate, but as mission concepts and the sail design matured, greater fidelity was required. Here we discuss the progress to a three-dimensional sail model that includes the effects of tension and thermal deformation that has been derived from a large structural Finite Element Model (FEM) developed by the Langley Research Center. We have found that the deformed sail membrane affects torque relatively much more than thrust; a flat plate model could potentially model thrust well enough to close mission design studies, but a three-dimensional solar sail is essential to control system design. The three-dimensional solar sail model revealed that thermal deformations of unshielded booms would create unacceptably large solar disturbance torques. The original large FEM model was used in control and mission simulations, but was resulted in simulations with prohibitive run times. This led us to adapt the Generalized Sail Model (GSM) of Rios-Reyes. A design reference sail model has been baselined for NEA Scout and has been used to design the mission and control system for the sailcraft. Additionally, since NEA Scout uses reaction wheels for attitude pointing and control, the solar torque model is essentially to successfully design the NEA Scout momentum management control system. We have also updated the estimate of diffusivity used for the aluminized sail material based on optical testing of wrinkled sail material. The model presented here represents the current state of the art of NASA's ability to model solar sail thrust and torque.

Heaton, Andy↗

A higher-order finite element reactive transport model for unstructured and fractured grids

Abstract This work presents a new reactive transport framework that combines a powerful geochemistry engine with advanced numerical methods for flow and transport in subsurface fractured porous media. Specifically, the PhreeqcRM interface (developed by the USGS) is used to take advantage of a large library of equilibrium and kinetic aqueous and fluid-rock reactions, which has been validated by numerous experiments and benchmark studies. Fluid flow is modeled by the Mixed Hybrid Finite Element (FE) method, which provides smooth velocity fields even in highly heterogenous formations with discrete fractures. A multilinear Discontinuous Galerkin FE method is used to solve the multicomponent transport problem. This method is locally mass conserving and its second order convergence significantly reduces numerical dispersion. In terms of thermodynamics, the aqueous phase is considered as a compressible fluid and its properties are derived from a Cubic Plus Association (CPA) equation of state. The new simulator is validated against several benchmark problems (involving, e.g., Fickian and Nernst-Planck diffusion, isotope fractionation, advection-dispersion transport, and rock-fluid reactions) before demonstrating the expanded capabilities offered by the underlying FE foundation, such as high computational efficiency, parallelizability, low numerical dispersion, unstructured 3D gridding, and discrete fraction modeling.

58 GEOSCIENCES↗

Reduction in Computational Cost of Progressive Failure Analysis of Composite Structures

Designing aircraft structures requires efficient modeling approaches to iterate on multiple structural configurations to achieve an optimal design. Typically, damage tolerance is not considered at the design stage because of the high computational cost in its implementation within a finite element modeling approach. Therefore, analytical or empirical approaches are often used to size critical damage-tolerant structures once an optimal design is determined. In this study, the Progressive Release eXplicit Virtual Crack-Closure Technique (PRX-VCCT) is assessed for its capability to cost-effectively evaluate skin stringer separation of a blade-stiffened panel that is subjected to seven point bend loads. An initial verification study was performed to evaluate PRX VCCT to accurately simulate skin-stringer separation with respect to existing cohesive element approaches. Furthermore, the influence of element size, ranging from 0.10 in. to 0.40 in., on the total computational time using the PRX-VCCT is investigated. The results indicate that the PRX-VCCT can be used to accurately simulate skin stringer separation using large element lengths (0.40 in.). Additionally, a significant reduction in the computational time to simulate skin stringer separation is observed using the PRX-VCCT. Large-scale progressive damage analysis using PRX-VCCT can be implemented early in the design cycle of composite structures without requiring a global-to-local modeling approach.

Composites↗

STRUCTURAL MODELING TO SUPPORT POST-YIELD ACCEPTANCE CRITERIA FOR SPENT NUCLEAR FUEL CLADDING

Spent nuclear fuel (SNF) is evaluated for structural failure during storage and transportation scenarios. The U.S. Department of Energy’s Spent Fuel and Waste Science and Technology (SFWST) program has sponsored significant research in quantifying mechanical loads on SNF during storage and transportation scenarios using experimental and modeling methods. The SFWST program has also performed significant research on measuring the mechanical behavior of irradiated SNF as defueled cladding segments and cladding with fuel pellets to measure composite behavior. This paper considers some of the key material data from the Sibling Pin testing and uses structural modeling and analysis methods that have been informed by testing to consider post-yield acceptance criteria for SNF cladding structural analysis. Test data published by Oak Ridge National Laboratory (ORNL) and Pacific Northwest National Laboratory (PNNL) are the foundation for informing the material behavior of the models developed in this study. In particular, four-point bend (4PB) tests of fueled and defueled cladding segments provide significant information about the bending failure mode of SNF. ORNL’s 4PB test data is on fueled cladding segments, so the composite behavior of SNF is demonstrated. This paper describes PNNL’s coincident beam model that was developed to approximate the composite behavior of SNF. This paper also presents PNNL’s structural dynamic finite element models of a cask tip-over scenario, which is predicted to cause the strongest mechanical loads on SNF of all postulated storage and transportation scenarios. SNF bending loads predicted in the cask tip-over scenario and cladding acceptance criteria beyond yield are considered, with justification based on the Sibling Pin test data. ASME Boiler and Pressure Vessel code stress intensity limits are also considered. The ultimate goal of this work is to aid in the justification of structural acceptance criteria for SNF cladding beyond the cladding’s irradiated yield strength for use in structural analysis of all storage and transportation scenarios.

Klymyshyn, Nicholas A.↗

Automated design optimization of supersonic airplane wing structures under dynamic constraints

The problems of the preliminary and first level detail design of supersonic aircraft wings are stated as mathematical programs and solved using automated optimum design techniques. The problem is approached in two phases: the first is a simplified equivalent plate model in which the envelope, planform and structural parameters are varied to produce a design, the second is a finite element model with fixed configuration in which the material distribution is varied. Constraints include flutter, aeroelastically computed stresses and deflections, natural frequency and a variety of geometric limitations.

Fox, R. L.↗

Coordinate transformation and truncation for rotating spacecraft with flexible appendages

The analytical structures of alternative coordinate transformations for the variables that characterize the deformations of finite element models of flexible appendages on rotating spacecraft are examined. Particular emphasis is given to the truncation operations that are essential for efficient numerical simulations of flexible spacecraft, due to the necessarily large number of coordinates employed in a finite element description. A theorem is developed which establishes necessary conditions for the commutativity of the required truncation and inversion operations. The recommended alternative for the rotating elastic appendage is a complex transformation with previously published orthogonality properties, which permit the derivation in this paper of a relatively simple explicit set of transformed and truncated equations of modal vibration in real variables.

Likins, P.↗

Thermal distortion analysis of the space station solar dynamic concentrator

A method was developed to evaluate the thermal distortion of the Space Station Solar Dynamic Concentrator and the effects of thermal distortion on concentrator optical performance. The analytical method includes generating temperature distributions with TRASYS and SINDA models, interfacing the SINDA results with the SINDA-NASTRAN Interface Program (SNIP), calculating thermal distortion with a NASTRAN/PATRAN finite element model, and providing flux distribution maps within the receiver with the ray tracing OFFSET program. Temperature distributions, thermally induced slope errors, and flux distribution maps within the receiver are discussed. Results during a typical orbit indicate that temperatures of the hexagonal panels and triangular facets range between -18 and 99 C (-1 to 210 F), facet rotations are less than 0.2 mrad, and a change in facet radius due to thermal flattening is less than 5 percent. The predicted power loss with thermal distortion effects was less than 0.3 percent. The thermal distortion of the Solar Dynamic concentrator has negligible effect on the flux distribution within the receiver cavity.

Trudell, Jeffrey J.↗