Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “structural modeling and simulation”

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 343 records · Page 19

Dynamic Shape Reconstruction of Three-Dimensional Frame Structures Using the Inverse Finite Element Method

A robust and efficient computational method for reconstructing the three-dimensional displacement field of truss, beam, and frame structures, using measured surface-strain data, is presented. Known as shape sensing , this inverse problem has important implications for real-time actuation and control of smart structures, and for monitoring of structural integrity. The present formulation, based on the inverse Finite Element Method (iFEM), uses a least-squares variational principle involving strain measures of Timoshenko theory for stretching, torsion, bending, and transverse shear. Two inverse-frame finite elements are derived using interdependent interpolations whose interior degrees-of-freedom are condensed out at the element level. In addition, relationships between the order of kinematic-element interpolations and the number of required strain gauges are established. As an example problem, a thin-walled, circular cross-section cantilevered beam subjected to harmonic excitations in the presence of structural damping is modeled using iFEM; where, to simulate strain-gauge values and to provide reference displacements, a high-fidelity MSC/NASTRAN shell finite element model is used. Examples of low and high-frequency dynamic motion are analyzed and the solution accuracy examined with respect to various levels of discretization and the number of strain gauges.

Gherlone, Marco↗

Variability of Disk Emission in Pre-Main-Sequence and Related Stars. IV. Investigating the Structural Changes in the Inner Disk Region of MWC 480

We present five epochs of near-IR observations of the protoplanetary disk around MWC 480 (HD 31648) obtained with the SpeX spectrograph on NASA’s Infrared Telescope Facility between 2007 and 2013, inclusive. Using the measured line fluxes in the Pa β and Br γ lines, we found the mass accretion rates to be (1.26–2.30) x 10(exp −7) Solar Mass/yr and (1.4–2.01) x 10(exp −7) Solar Mass/yr, respectively, but which varied by more than 50% from epoch to epoch. The spectral energy distribution reveals a variability of about 30% between 1.5 and 10 μm during this same period of time. We investigated the variability using of the continuum emission of the disk in using the Monte- Carlo Radiative Transfer Code HOCHUNK3D. We find that varying the height of the inner rim successfully produces a change in the NIR flux but lowers the far-IR emission to levels below all measured fluxes. Because the star exhibits bipolar flows, we utilized a structure that simulates an inner disk wind to model the variability in the near-IR, without producing flux levels in the far-IR that are inconsistent with existing data. For this object, variable near-IR emission due to such an outflow is more consistent with the data than changing the scale height of the inner rim of the disk.

Fernandes, Rachel B.↗

Mature thunderstorm cloud top structure - Three-dimensional numerical simulation versus satellite observations

Preliminary results of four runs with a three-dimensional model of the effects of vertical wind shear on cloud top height/temperature structure and the internal properties of isolate midlatitude thunderstorms are reported. The model is being developed as an aid to analyses of GEO remote sensing satellite data. The grid is a 27 x 27 x 20 mesh with 2 km horizontal resolution and 0.9 vertical resolution. The total grid is 54 km on a side and 18 km deep. A second-order Crowley scheme for advecting momentum is extended with a third-order correction for spatial truncation error, and the earth-relative horizontal surface wind components are decreased to 50 percent of their values at 0.45 km. A temperature increase with height is included, together with an initial impulse consisting of a nonrotating cylindrical weak buoyant updraft 10 km in radius. The results of the runs are discussed in terms of the time variation of the vertical velocity extrema, the effects of strong and weak shear on a storm, the cloud top height, the Lagrangian dynamics of a thermal couplet, and data from a real storm.

Schlesinger, R. E.↗

Development of impact design methods for ceramic gas turbine components

Impact damage prediction methods are being developed to aid in the design of ceramic gas turbine engine components with improved impact resistance. Two impact damage modes were characterized: local, near the impact site, and structural, usually fast fracture away from the impact site. Local damage to Si3N4 impacted by Si3N4 spherical projectiles consists of ring and/or radial cracks around the impact point. In a mechanistic model being developed, impact damage is characterized as microcrack nucleation and propagation. The extent of damage is measured as volume fraction of microcracks. Model capability is demonstrated by simulating late impact tests. Structural failure is caused by tensile stress during impact exceeding material strength. The EPIC3 code was successfully used to predict blade structural failures in different size particle impacts on radial and axial blades.

Song, J.↗

Simulation of Dynamic Crushing Behavior of Simulated Carbon/Epoxy Stanchions Using MAT58 and MAT213 in LS-DYNA

Carbon fiber-reinforced polymer (CFRP) composites have shown promise as a material for structures designed to absorb energy in crush-style impact loading. In complementing the lightweight and tailorable characteristics of CFRP structures, the benefits of impact energy absorption are particularly interesting in aerospace vehicles that require excellent crashworthiness performance. However, simulating the behavior of CFRP structures in dynamic, crushing-style impact loading is challenging because of the many damage and failure modes that are essential to capture in the model. For the present investigation, CFRP stanchions were tested using a crash sled experimental test rig. The stanchions were designed by the Composite Materials Handbook-17 (CMH-17) Crashworthiness Working Group for the purpose of comparing experimental crash sled tests to corresponding simulations of the tests. The stanchions are C-channel shaped and represent a geometry common in the interior of aerospace vehicle structures. Explicit simulations in LSDYNA were performed using a well-established composite material model (MAT58) and a next-generation material model (MAT213). Fully integrated shell elements were utilized instead of 3D solid elements to limit computation time. In all models, each of the sixteen plies were represented with individual layers of shell elements with tiebreak contact between each adjacent ply pair to simulate interlaminar fracture. Simulating the crushing of the stanchions occurred in two phases. First, the material models were calibrated using flat specimens that were manufactured with the same layup as the stanchions. While the original goal was calibration of material-related properties, meshdependent behavior was observed in simulations with either material model, and an unstructured mesh was selected to remediate undesirable mesh-dependent failure modes. Additionally, for both the MAT58 and MAT213 models, it was found that either the crush force or the failure mode could be modeled accurately, but no set of parameters could be identified to attain both results in the same model. Once satisfactory calibration was achieved, the same material parameters were applied to the stanchion crush simulations. The stanchion simulations showed that MAT213 more accurately predicted the experimentally determined crush force, and both material models predicted key aspects of the experimentally observed failure modes.

Composite Materials↗

Ocean Modeling of the North Atlantic

Present modeling of the North Atlantic is inadequate and can be improved in a number of ways. A number of important physical processes are listed in five categories from the viewpoints of how they are treated in isolation, how they are usually represented in present ocean basin models, and how they may be better represented in future models. In the first two categories of vertical boundary processes and internal vertical mixing, parameterizations exist which can easily be incorporated into models and which will have important effects on the simulated structure of the North Atlantic. For the third catagory (mesoscale eddy effects), adequate parameterizations do not exist; but the order of magnitude of the effects is known from observational and process-model studies. A horizontal grid spacing of 100 km or less in required to allow parameterizations with this order of magnitude, as well as to resolve the time-averaged ocean fields. In the fourth category of large scale transports improvements are suggested by way of increased vertical resolution and by the requirement that lateral mixing due to eddies takes place on isopycnal surfaces. Model incorporation of the latter phenomenta is underway. In the fifth category of miscellaneous high-latitude processes, formulations for the treatment of sea ice are available for use. However, the treatment of gravitational instability, which is crucial to deepwater formation in the Atlantic Ocean, will require additional refinements to account for the unresolved physics of chimney formations in the open ocean and buoyant plumes near ocean boundaries.

Albert J Semtner↗

Resin Film Infusion (RFI) Process Modeling for Large Transport Aircraft Wing Structures

This investigation completed the verification of a three-dimensional resin transfer molding/resin film infusion (RTM/RFI) process simulation model. The model incorporates resin flow through an anisotropic carbon fiber preform, cure kinetics of the resin, and heat transfer within the preform/tool assembly. The computer model can predict the flow front location, resin pressure distribution, and thermal profiles in the modeled part. The formulation for the flow model is given using the finite element/control volume (FE/CV) technique based on Darcy's Law of creeping flow through a porous media. The FE/CV technique is a numerically efficient method for finding the flow front location and the fluid pressure. The heat transfer model is based on the three-dimensional, transient heat conduction equation, including heat generation. Boundary conditions include specified temperature and convection. The code was designed with a modular approach so the flow and/or the thermal module may be turned on or off as desired. Both models are solved sequentially in a quasi-steady state fashion. A mesh refinement study was completed on a one-element thick model to determine the recommended size of elements that would result in a converged model for a typical RFI analysis. Guidelines are established for checking the convergence of a model, and the recommended element sizes are listed. Several experiments were conducted and computer simulations of the experiments were run to verify the simulation model. Isothermal, non-reacting flow in a T-stiffened section was simulated to verify the flow module. Predicted infiltration times were within 12-20% of measured times. The predicted pressures were approximately 50% of the measured pressures. A study was performed to attempt to explain the difference in pressures. Non-isothermal experiments with a reactive resin were modeled to verify the thermal module and the resin model. Two panels were manufactured using the RFI process. One was a stepped panel and the other was a panel with two 'T' stiffeners. The difference between the predicted infiltration times and the experimental times was 4% to 23%.

Loos, Alfred C.↗

The dependence of coronal hole size on large scale magnetic field strength

The importance of mathematical models of the coronal structure for studies of coronal energetics, to simulate global flows of the solar wind, and to obtain reliable solar terrestrial predictions is discussed. Previous coronal models, including an example of a coronal MHD flow model, are reviewed. The development of a coronal model which is a logical extension of earlier models and which allows a closer relationship to the photospheric magnetic field as it is observed daily is described. The calculations are outlined. The assumptions of the model are: axisymmetric flow with no rotation, resulting in two dimensional flow in a meridional plane; zero viscosity and infinite electrical conductivity; polytropic, single fluid flow; and no momentum addition.

Suess, S. T.↗

Air pollutant transport in a coastal environment. Part 1: Two-dimensional simulations of sea-breeze and mountain effects

Over the southern California coastal region, observations of the vertical distributions of pollutants show that maximum concentrations can occur within temperature inversion layers well above the surface. A mesoscale model is used to study the dynamical phenomena that cause such layers, including sea breezes and mountain flows, and to study the characteristics of air pollutant transport in a coastal environment capped by a temperature inversion. The mathematical and physical structure of the model is described. Two-dimensional simulations corresponding to four configurations of coastal plains and mountains are discussed. The simulations reveal that pollutant transport over a coastal plain is strongly influenced by the topographic configuration, including the height of coastal mountains and their distance from the coastline. Sea breezes induced by land-sea thermal contrasts, as well as upslope winds induced along mountain flanks, both create vertical transport that can lead to the formation of elevated pollution layers. The sea-breeze circulation generates pollution layers by undercutting the mixed layer and lofting pollutants into the stable layer. Heating of mountain slopes acts to vent pollutants above the mountain ridge during the day; during the evening, pollutants can be injected directly into the inversion layer from the decaying upslope flows. In a land-sea configuration with mountains close to the coastline, the sea breeze and heated-mountain flow are strongly coupled. In the afternoon, this interaction can produce upslope flow from which polluted air is detrained into the inversion layer as a return circulation. When the mountains lie farther inland, however, pollutants may be trapped aloft when the mixed layer stabilizes in the late afternoon. As the nocturnal boundary layer forms over the coast in the evening, polluted mixed-layer air is effectively left behind in the inversion layer. In the Los Angeles Basin, the formation mechanism for elevated polluted layers is most similar to our cases with inland mountains.

Lu, Rong↗

Dynamics and control simulation of large space structures

Attention is given to the dynamics and control simulation techniques applicable to large space structures, encompassing the model of the large structure as well as models of its dynamics, its onboard controller, its sensor/actuator systems, and various disturbances to which it may be subjected. Control performance parameters must be computed as functions of time, or of such other design parameters as control gain. A large space antenna is used as an example of the structures in question, illustrating the way in which an efficient dynamics and control simulation structure program can be developed, and how antenna parameters such as reflector surface distortions can be computed, without having to solve a large number of equations for each computation time interval.

Lin, Y. H.↗

GENOA-PFA: Progressive Fracture in Composites Simulated Computationally

GENOA-PFA is a commercial version of the Composite Durability Structural Analysis (CODSTRAN) computer program that simulates the progression of damage ultimately leading to fracture in polymer-matrix-composite (PMC) material structures under various loading and environmental conditions. GENOA-PFA offers several capabilities not available in other programs developed for this purpose, making it preferable for use in analyzing the durability and damage tolerance of complex PMC structures in which the fiber reinforcements occur in two- and three-dimensional weaves and braids. GENOA-PFA implements a progressive-fracture methodology based on the idea that a structure fails when flaws that may initially be small (even microscopic) grow and/or coalesce to a critical dimension where the structure no longer has an adequate safety margin to avoid catastrophic global fracture. Damage is considered to progress through five stages: (1) initiation, (2) growth, (3) accumulation (coalescence of propagating flaws), (4) stable propagation (up to the critical dimension), and (5) unstable or very rapid propagation (beyond the critical dimension) to catastrophic failure. The computational simulation of progressive failure involves formal procedures for identifying the five different stages of damage and for relating the amount of damage at each stage to the overall behavior of the deteriorating structure. In GENOA-PFA, mathematical modeling of the composite physical behavior involves an integration of simulations at multiple, hierarchical scales ranging from the macroscopic (lamina, laminate, and structure) to the microscopic (fiber, matrix, and fiber/matrix interface), as shown in the figure. The code includes algorithms to simulate the progression of damage from various source defects, including (1) through-the-thickness cracks and (2) voids with edge, pocket, internal, or mixed-mode delaminations.

Murthy, Pappu L. N.↗

Simulating the Martian Middle and Upper Atmosphere with the NASA Ames Mars GCM.

The NASA Ames Mars Global Climate Model (MGCM) continues to be utilized to examine the dynamical coupling between the lower atmosphere and middle/upper atmosphere. This coupling impacts the transport of water, dust, and chemistry which augments the thermal structure and circulation. The interactions and coupling within the atmosphere are important for understanding the present-day and past atmosphere and climate. MGCM simulations are compared with observations from recent Mars missions to provide model constraints, and in turn, the simulations can then provide context to the observations. This cycle is important to move us forward in our understanding of planetary atmospheres. We will present on one of the fundamental atmospheric features, the thermal structure, and how the physics in the MGCM can elucidate the observed diurnal and seasonal behavior. Additionally, we will focus on the distribution of chemical species which influences the thermal structure and provides additional information on the coupling between the lower and middle/upper atmosphere. Lastly, we will provide a summary of the model capabilities, as the MGCM continues to be improved upon to properly simulate Mars’ atmosphere from the ground to ~250 km.

Amanda Susanne Brecht↗

Advances and trends in computational structural mechanics

Recent developments in computational structural mechanics are reviewed with reference to computational needs for future structures technology, advances in computational models for material behavior, discrete element technology, assessment and control of numerical simulations of structural response, hybrid analysis, and techniques for large-scale optimization. Research areas in computational structural mechanics which have high potential for meeting future technological needs are identified. These include prediction and analysis of the failure of structural components made of new materials, development of computational strategies and solution methodologies for large-scale structural calculations, and assessment of reliability and adaptive improvement of response predictions.

Noor, A. K.↗

Mars atmospheric dynamics as simulated by the NASA Ames General Circulation Model. II - Transient baroclinic eddies

A large set of experiments performed with the NASA Ames Mars General Circulation Model is analyzed to determine the properties, structure, and dynamics of the simulated transient baroclinic eddies. There is strong transient baroclinic eddy activity in the extratropics of the Northern Hemisphere during the northern autumn, winter, and spring seasons. The eddy activity remains strong for very large dust loadings, though it shifts northward. The eastward propagating eddies are characterized by zonal wavenumbers of 1-4 and periods of about 2-10 days. The properties of the GCM baroclinic eddies in the northern extratropics are compared in detail with analogous properties inferred from Viking Lander meteorology observations.

Barnes, Jeffrey R.↗

Experimental, physical and numerical modeling of lunar regolith and lunar regolith structures

An extensive series of laboratory strength and deformation experiments have been performed on a lunar regolith simulant. Results of these experiments are compared to results from experiments on real lunar regolith from the Apollo and Luna missions to illustrate the suitability of this material in capturing the engineering properties lunar regolith. In addition, these results are used to calibrate a constitutive model used to describe its stress-strain behavior. This model, in conjunction with numerical analysis techniques, is used to predict the response i.e. material parameters) of lunar simulant under 1/6-g and low confining stress conditions. These tools are also used to predict the displacement response of a lunar soil embankment structure used to cover a first generation human habitat module, which might be used to accommodate the first astronauts revisiting the moon. These predictions are compared to physical models of this structure, which are tested in a geotechnical centrifuge in order to satisfy scaling relationships between prototype and model.

Perkins, Steven W.↗

The use of typed lambda calculus for comprehension and construction of simulation models in the domain of ecology

We are concerned with two important issues in simulation modelling: model comprehension and model construction. Model comprehension is limited because many important choices taken during the modelling process are not documented. This makes it difficult for models to be modified or used by others. A key factor hindering model construction is the vast modelling search space which must be navigated. This is exacerbated by the fact that many modellers are unfamiliar with the terms and concepts catered to by current tools. The root of both problems is the lack of facilities for representing or reasoning about domain concepts in current simulation technology. The basis for our achievements in both of these areas is the development of a language with two distinct levels; one for representing domain information, and the other for representing the simulation model. Of equal importance, is the fact that we make formal connections between these two levels. The domain we are concerned with is ecological modelling. This language, called Elklogic, is based on the typed lambda calculus. Important features include a rich type structure, the use of various higher order functions, and semantics. This enables complex expressions to be constructed from relatively few primitives. The meaning of each expression can be determined in terms of the domain, the simulation model, or the relationship between the two. We describe a novel representation for sets and substructure, and a variety of other general concepts that are especially useful in the ecological domain. We use the type structure in a novel way: for controlling the modelling search space, rather than a proof search space. We facilitate model comprehension by representing modelling decisions that are embodied in the simulation model. We represent the simulation model separately from, but in terms of a domain mode. The explicit links between the two models constitute the modelling decisions. The semantics of Elklogic enables English text to be generated to explain the simulation model in domain terms.

Uschold, Michael↗

Modeling Validation and Its Design Relevance for NTP Fuel Elements

The development of Nuclear Thermal Propulsion engines depends heavily on modeling and simulation of a number of disciplines—neutronic / fluid / thermal / structural are only a start. Further, rockets involve extreme conditions, particularly high material temperatures in the nuclear fuel. This paper considers validation of the underling physical models and codes and how their accuracy influences several design challenges. Several issues are identified, where models, material data, could enable design solutions. In particular, thermal and CTE mismatch stresses are a concern, and stress predictions would benefit from improved mechanical property data at high temperatures. Nuclear inter-element heating causes propellant flow maldistribution and performance risks, but a simple model provides insights. High temperature chemical diffusion should be modeled to understand fuel element mass loss.

Nuclear Propulsion↗

Modeling Validation and Its Design Relevance for NTP Fuel Elements

The development of Nuclear Thermal Propulsion engines depends heavily on modeling and simulation of a number of disciplines—neutronic / fluid / thermal / structural are only a start. Further, rockets involve extreme conditions, particularly high material temperatures in the nuclear fuel. This paper considers validation of the underling physical models and codes and how their accuracy influences several design challenges. Several issues are identified, where models, material data, could enable design solutions. In particular, thermal and CTE mismatch stresses are a concern, and stress predictions would benefit from improved mechanical property data at high temperatures. Nuclear inter-element heating causes propellant flow maldistribution and performance risks, but a simple model provides insights. High temperature chemical diffusion should be modeled to understand fuel element mass loss.

Nuclear Propulsion↗