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 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
NASA NTRS2004

On Efficient Multigrid Methods for Materials Processing Flows with Small Particles

Multiscale modeling of materials requires simulations of multiple levels of structural hierarchy. The computational efficiency of numerical methods becomes a critical factor for simulating large physical systems with highly desperate length scales. Multigrid methods are known for their superior efficiency in representing/resolving different levels of physical details. The efficiency is achieved by employing interactively different discretizations on different scales (grids). To assist optimization of manufacturing conditions for materials processing with numerous particles (e.g., dispersion of particles, controlling flow viscosity and clusters), a new multigrid algorithm has been developed for a case of multiscale modeling of flows with small particles that have various length scales. The optimal efficiency of the algorithm is crucial for accurate predictions of the effect of processing conditions (e.g., pressure and velocity gradients) on the local flow fields that control the formation of various microstructures or clusters.

Thomas, James↗
NASA NTRS2013

Jovian Plasmas Torus Interaction with Europa. Plasma Wake Structure and Effect of Inductive Magnetic Field: 3D Hybrid Kinetic Simulation

The hybrid kinetic model supports comprehensive simulation of the interaction between different spatial and energetic elements of the Europa moon-magnetosphere system with respect to a variable upstream magnetic field and flux or density distributions of plasma and energetic ions, electrons, and neutral atoms. This capability is critical for improving the interpretation of the existing Europa flyby measurements from the Galileo Orbiter mission, and for planning flyby and orbital measurements (including the surface and atmospheric compositions) for future missions. The simulations are based on recent models of the atmosphere of Europa (Cassidy et al., 2007; Shematovich et al., 2005). In contrast to previous approaches with MHD simulations, the hybrid model allows us to fully take into account the finite gyroradius effect and electron pressure, and to correctly estimate the ion velocity distribution and the fluxes along the magnetic field (assuming an initial Maxwellian velocity distribution for upstream background ions). Photoionization, electron-impact ionization, charge exchange and collisions between the ions and neutrals are also included in our model. We consider the models with Oþ þ and Sþ þ background plasma, and various betas for background ions and electrons, and pickup electrons. The majority of O2 atmosphere is thermal with an extended non-thermal population (Cassidy et al., 2007). In this paper, we discuss two tasks: (1) the plasma wake structure dependence on the parameters of the upstream plasma and Europa's atmosphere (model I, cases (a) and (b) with a homogeneous Jovian magnetosphere field, an inductive magnetic dipole and high oceanic shell conductivity); and (2) estimation of the possible effect of an induced magnetic field arising from oceanic shell conductivity. This effect was estimated based on the difference between the observed and modeled magnetic fields (model II, case (c) with an inhomogeneous Jovian magnetosphere field, an inductive magnetic dipole and low oceanic shell conductivity).

Jovian magnetosphere↗
NASA NTRS2004

A Method for Incorporating Changing Structural Characteristics Due to Propellant Mass Usage in a Launch Vehicle Ascent Simulation

Launch vehicles consume large quantities of propellant quickly, causing the mass properties and structural dynamics of the vehicle to change dramatically. Currently, structural load assessments account for this change with a large collection of structural models representing various propellant fill levels. This creates a large database of models complicating the delivery of reduced models and requiring extensive work for model changes. Presented here is a method to account for these mass changes in a more efficient manner. The method allows for the subtraction of propellant mass as the propellant is used in the simulation. This subtraction is done in the modal domain of the vehicle generalized model. Additional computation required is primarily for constructing the used propellant mass matrix from an initial propellant model and further matrix multiplications and subtractions. An additional eigenvalue solution is required to uncouple the new equations of motion; however, this is a much simplier calculation starting from a system that is already substantially uncoupled. The method was successfully tested in a simulation of Saturn V loads. Results from the method are compared to results from separate structural models for several propellant levels, showing excellent agreement. Further development to encompass more complicated propellant models, including slosh dynamics, is possible.

McGhee, D. S.↗
NASA NTRS1999

Development of GENOA Progressive Failure Parallel Processing Software Systems

A capability consisting of software development and experimental techniques has been developed and is described. The capability is integrated into GENOA-PFA to model polymer matrix composite (PMC) structures. The capability considers the physics and mechanics of composite materials and structure by integration of a hierarchical multilevel macro-scale (lamina, laminate, and structure) and micro scale (fiber, matrix, and interface) simulation analyses. The modeling involves (1) ply layering methodology utilizing FEM elements with through-the-thickness representation, (2) simulation of effects of material defects and conditions (e.g., voids, fiber waviness, and residual stress) on global static and cyclic fatigue strengths, (3) including material nonlinearities (by updating properties periodically) and geometrical nonlinearities (by Lagrangian updating), (4) simulating crack initiation. and growth to failure under static, cyclic, creep, and impact loads. (5) progressive fracture analysis to determine durability and damage tolerance. (6) identifying the percent contribution of various possible composite failure modes involved in critical damage events. and (7) determining sensitivities of failure modes to design parameters (e.g., fiber volume fraction, ply thickness, fiber orientation. and adhesive-bond thickness). GENOA-PFA progressive failure analysis is now ready for use to investigate the effects on structural responses to PMC material degradation from damage induced by static, cyclic (fatigue). creep, and impact loading in 2D/3D PMC structures subjected to hygrothermal environments. Its use will significantly facilitate targeting design parameter changes that will be most effective in reducing the probability of a given failure mode occurring.

Abdi, Frank↗
NASA NTRS2020

Comparison of Surface Fluxes Derived from CYGNSS and Simulated by WRF Model: An MJO Case Study

This study focuses on ocean surface fluxes, mainly the latent heat flux, and their impact on MJO propagation and associated precipitation structures over the Indian Ocean and Maritime Continent. The Coupled-Ocean-Atmosphere-Wave-Sediment Transport (COAWST) model is used to simulate two MJO events during the 2017-2018 season: the December 10 - January 20, 2017 case, which maintained its strong precipitation signal over the Maritime Continent, and the March 1 - 20, 2018 case, which was weaker and did not propagate through the Maritime Continent. Both simulated MJO events show positive biases in surface rainfall compared with GPM IMERG data. During the MJO suppressed phase, the simulations rain more often than the observations. During the active phase, the westward propagating precipitation structures are more organized and much stronger compared with the observations, sometimes forming westward propagating cyclones that weakened the eastward precipitation signals. Two aspects of the surface flux interactions are investigated: the impact of the domain mean surface fluxes, and the impact of storm scale circulations and their interactions with local surface fluxes. Both aspects affect water vapor budget, atmosphere instability and mean flow, through which convection initiation, organization, and propagation are influenced. Model sensitivity tests with different radiation, microphysics, PBL schemes and nudging schemes indicate that in the control simulations, higher SST and surface fluxes, especially during the suppressed period, are the main reason of rainfall overestimation compared with IMERG data. The strong westward propagating signals are caused by both increased atmosphere instability and reduced mean wind shear. Unfortunately, the small differences in mean SST and surface fluxes between different model sensitivity tests are all within the satellite observation error margin, and cannot be directly corroborated by observations. One of the advantages of CYGNSS satellites is that they observe ocean surface wind and heat fluxes underneath strong rainfall events such as the convective systems associated with MJO active phases. Currently we are comparing CYGNSS level 2 surface fluxes retrievals and the model simulations in order to better understand the second aspect of the MJO and surface fluxes interactions, and how this affects MJO strengths and propagations. The interactive atmosphere-ocean-wave model also provides cases that directly comparing satellite observables (the bistatic radar cross section) and the model simulations (through CYGNSS satellite simulator). These discrepancies are more prominent in coupled ocean simulations, mainly due to higher SST and enhanced surface fluxes.

Li, Xiaowen↗
NASA NTRS1978

AAA gunnermodel based on observer theory

The Luenberger observer theory is used to develop a predictive model of a gunner's tracking response in antiaircraft artillery systems. This model is composed of an observer, a feedback controller and a remnant element. An important feature of the model is that the structure is simple, hence a computer simulation requires only a short execution time. A parameter identification program based on the least squares curve fitting method and the Gauss Newton gradient algorithm is developed to determine the parameter values of the gunner model. Thus, a systematic procedure exists for identifying model parameters for a given antiaircraft tracking task. Model predictions of tracking errors are compared with human tracking data obtained from manned simulation experiments. Model predictions are in excellent agreement with the empirical data for several flyby and maneuvering target trajectories.

Kou, R. S.↗
NASA NTRSDate not supplied

Hot Structure Concepts and Materials for Space Exploration Project

Hot structures are multifunctional material systems that integrate the functions of thermal protection systems (TPS) and primary structure to increase mass efficiency in aerospace vehicles. The goal of this project was to improve our capabilities for developing lightweight, reusable, and reliable hot structures concepts to enable advanced space exploration missions. Specifically, current hot structure carbon/carbon (C/C) composite material systems have low interlaminar strength, which can lead to material failures at low load levels, and large material property variability in as-fabricated parts. Two parallel paths were followed to address these shortcomings. First, analysis models were developed to identify the conditions under which excessive interlaminar stresses occur and to evaluate the potential of stitching the C/C laminae to improve the interlaminar thermal-structural response. This effort included the development of high-fidelity, fracture-based thermal-structural finite element analysis models to evaluate critical features needed to produce efficient hot structure concepts. The second path focused on developing atomistic and mesoscale simulation tools to model the fabrication process in an effort to reduce the time and cost of production while improving the quality of the resulting parts. Aspects of the process currently being investigated include the infusion of the matrix precursor into the fiber preform, the cure of the precursor molecules into their polymeric form, and the carbonization and graphitization of the polymer to form the carbon matrix material. Results of the 6-month study are presented herein.

thermal↗
NASA NTRS1994

Composite Load Spectra for Select Space Propulsion Structural Components

Generic load models are described with multiple levels of progressive sophistication to simulate the composite (combined) load spectra (CLS) that are induced in space propulsion system components, representative of Space Shuttle Main Engines (SSME), such as transfer ducts, turbine blades and liquid oxygen (LOX) posts. These generic (coupled) models combine the deterministic models for composite load dynamic, acoustic, high-pressure and high rotational speed, etc., load simulation using statistically varying coefficients. These coefficients are then determined using advanced probabilistic simulation methods with and without strategically selected experimental data. The entire simulation process is included in a CLS computer code. Applications of the computer code to various components in conjunction with the PSAM (Probabilistic Structural Analysis Method) to perform probabilistic load evaluation and life prediction evaluations are also described to illustrate the effectiveness of the coupled model approach.

Ho, Hing W.↗
NASA NTRS1992

Pointing-error simulations of the DSS-13 antenna due to wind disturbances

Accurate spacecraft tracking by the NASA Deep Space Network (DSN) antennas must be assured during changing weather conditions. Wind disturbances are the main source of tracking errors. The development of a wind-force model and simulations of wind-induced pointing errors of DSN antennas are presented. The antenna model includes the antenna structure, the elevation and azimuth servos, and the tracking controller. Simulation results show that pointing errors due to wind gusts are of the same order as errors due to static wind pressure and that these errors (similar to those of static wind pressure) satisfy the velocity quadratic law. The presented methodology is used for wind-disturbance estimation and for the design of an antenna controller with wind-disturbance rejection properties.

Gawronski, W.↗
NASA NTRS2020

Structure and Parameter Uncertainty in Centennial Projections of Forest Community Structure and Carbon Cycling

Secondary forest regrowth shapes community succession and biogeochemistry for decades, including in the Upper Great Lakes region. Vegetation models encapsulate our understanding of forest function, and whether models can reproduce multi‐decadal succession patterns is an indication of our ability to predict forest responses to future change. We test the ability of a vegetation model to simulate C cycling and community composition during 100 years of forest regrowth following stand‐replacing disturbance, asking (a) Which processes and parameters are most important to accurately model Upper Midwest forest succession? (b) What is the relative importance of model structure versus parameter values to these predictions? We ran ensembles of the Ecosystem Demography model v2.2 with different representations of processes important to competition for light. We compared the magnitude of structural and parameter uncertainty and assessed which sub‐model–parameter combinations best reproduced observed C fluxes and community composition. On average, our simulations underestimated observed net primary productivity (NPP) and leaf area index (LAI) after 100 years and predicted complete dominance by a single plant functional type (PFT). Out of 4,000 simulations, only nine fell within the observed range of both NPP and LAI, but these predicted unrealistically complete dominance by either early hardwood or pine PFTs. A different set of seven simulations were ecologically plausible but under‐predicted observed NPP and LAI. Parameter uncertainty was large; NPP and LAI ranged from ~0% to >200% of their mean value, and any PFT could become dominant. The two parameters that contributed most to uncertainty in predicted NPP were plant–soil water conductance and growth respiration, both unobservable empirical coefficients. We conclude that (a) parameter uncertainty is more important than structural uncertainty, at least for ED‐2.2 in Upper Midwest forests and (b) simulating both productivity and plant community composition accurately without physically unrealistic parameters remains challenging for demographic vegetation models.

canopy radiative transfer↗
NASA NTRS1986

Simulation of double layers in a model auroral circuit with nonlinear impedance

A reduced circuit description of the U-shaped potential structure of a discrete auroral arc, consisting of the flank transmission line plus parallel-electric-field region, is used to provide the boundary condition for one-dimensional simulations of the double-layer evolution. The model yields asymptotic scalings of the double-layer potential, as a function of an anomalous transport coefficient alpha and of the perpendicular length scale l(a) of the arc. The arc potential phi(DL) scales approximately linearly with alpha, and for alpha fixed phi (DL) about l(a) to the z power. Using parameters appropriate to the auroral zone acceleration region, potentials of phi (DPL) 10 kV scale to projected ionospheric dimensions of about 1 km, with power flows of the order of magnitude of substorm dissipation rates.

Smith, R. A.↗
NASA NTRS1998

A New Approach to Overcoming Spatial Aliasing in Structural Damage Detection

Aircraft, reusable launch vehicles,unmanned aircraft, and other advanced structures are being built using lightweight composite materials/metals with design safety factors as low as 1.25. These advanced structures operate in uncertain and severe environments and are susceptible to damage such as delamination, fiber/matrix damage, hydrothermal strain in composite materials, and fatigue and cracking in metals. To ensure human safety and load-bearing integrity these structures must be inspected to detect and locate often invisible damage and faults before they become catastrophic. Conventional methods of non-destructive evaluation sometimes miss significant damage and are time consuming and expensive to perform. In contrast, vibrometry or vibration signature techniques are a global method of structural integrity monitoring that potentially can efficiently detect damage on large structures, including damage that is away from sensor locations, and in the interior of structures. However, a barrier problem in damage detection using vibration measurements is the need to measure the vibration response at a large number of points on the structure. Typically, model reduction or expansion procedures such as Guyan reduction or dynamic expansion are attempted to overcome the problem of insufficient measurements or spatial aliasing. These approaches depend on using information from the healthy model, and thus put error into the reduction used to represent the damaged structure. In this paper, a Frequency Response Function technique is used to detect damage to a fixed-free beam. The technique uses measured frequency response functions from the healthy structure as reference data, and then monitors vibration measurements during the life of the structure to detect damage. In an analytical simulation using a finite-element model of a beam, damage was located using only sparse measurements because the technique uses both rotation and translation measurements from the damaged structure.

Schulz, Mark J.↗
NASA NTRS2006

Bearing-Load Modeling and Analysis Study for Mechanically Connected Structures

Bearing-load response for a pin-loaded hole is studied within the context of two-dimensional finite element analyses. Pin-loaded-hole configurations are representative of mechanically connected structures, such as a stiffener fastened to a rib of an isogrid panel, that are idealized as part of a larger structural component. Within this context, the larger structural component may be idealized as a two-dimensional shell finite element model to identify load paths and high stress regions. Finite element modeling and analysis aspects of a pin-loaded hole are considered in the present paper including the use of linear and nonlinear springs to simulate the pin-bearing contact condition. Simulating pin-connected structures within a two-dimensional finite element analysis model using nonlinear spring or gap elements provides an effective way for accurate prediction of the local effective stress state and peak forces.

Knight, Norman F., Jr.↗
NASA NTRS1992

Progress on single barrier varactors for submillimeter wave power generation

Theoretical work on Single Barrier Varactor (SBV) diodes, indicate that the efficiency for a multiplier has a maximum for a considerably smaller capacitance variation than previously thought. The theoretical calculations are performed, both with a simple theoretical model and a complete computer simulation using the method of harmonic balance. Modeling of the SBV is carried out in two steps. First, the semiconductor transport equations are solved simultaneously using a finite difference scheme in one dimension. Secondly, the calculated I-V, and C-V characteristics are input to a multiplier simulator which calculates the optimum impedances, and output powers at the frequencies of interest. Multiple barrier varactors can also be modeled in this way. Several examples on how to design the semiconductor layers to obtain certain characteristics are given. The calculated conversion efficiencies of the modeled structures, in a multiplier circuit, are also presented. Computer simulations for a case study of a 750 GHz multiplier show that InAs diodes perform favorably compared to GaAs diodes. InAs and InGaAs SBV diodes have been fabricated and their current vs. voltage characteristics are presented. In the InAs diode, was the large bandgap semiconductor AlSb used as barrier. The InGaAs diode was grown lattice matched to an InP substrate with InAlAs as a barrier material. The current density is greatly reduced for these two material combinations, compared to that of GaAs/AlGaAs SBV diodes. GaAs based diodes can be biased to higher voltages than InAs diodes.

Nilsen, Svein M.↗
NASA NTRS2013

Uncertainty in Simulating Wheat Yields Under Climate Change

Projections of climate change impacts on crop yields are inherently uncertain1. Uncertainty is often quantified when projecting future greenhouse gas emissions and their influence on climate2. However, multi-model uncertainty analysis of crop responses to climate change is rare because systematic and objective comparisons among process-based crop simulation models1,3 are difficult4. Here we present the largest standardized model intercomparison for climate change impacts so far. We found that individual crop models are able to simulate measured wheat grain yields accurately under a range of environments, particularly if the input information is sufficient. However, simulated climate change impacts vary across models owing to differences in model structures and parameter values. A greater proportion of the uncertainty in climate change impact projections was due to variations among crop models than to variations among downscaled general circulation models. Uncertainties in simulated impacts increased with CO2 concentrations and associated warming. These impact uncertainties can be reduced by improving temperature and CO2 relationships in models and better quantified through use of multi-model ensembles. Less uncertainty in describing how climate change may affect agricultural productivity will aid adaptation strategy development and policymaking.

exhaust gases↗
NASA NTRSDate not supplied

The Impact of All-Sky Hyperspectral Infrared Radiance Assimilation on the Simulation and Forecast of Hurricane Sally in GEOS

Hyperspectral infrared (IR) radiance observations have been one of the major data sources assimilated in the data assimilation system over the last 20 years. However, observations peaking in the lower and mid-troposphere are underutilized in clear-sky radiance data assimilation as the quality control procedure removes a significant portion of cloud-affected observations from the assimilation. To include these cloud-affected observations, one approach is to assimilate them directly under all-sky conditions. The framework for all-sky assimilation of hyperspectral IR radiance observations has been developed in GEOS. The preliminary evaluation of simulated cloud-affected IR observations has been conducted, along with the corresponding sensitivities with respect to all hydrometeors. The symmetric cloud effect has been used as a cloud proxy in the observation error modeling, providing a balanced representation that mitigates the discrepancies between observations and model simulations for IR all-sky assimilation. The observation errors vary for different cloud conditions and have been modeled as cloud amount dependent and inter-channel correlated. Preliminary evaluation of model simulations for IR all-sky assimilation indicates that the simulated brightness temperature tends to have broader structures and lacks small-scale details. Excessive clouds generated by the model are also observed. To address this issue, Hurricane Sally (2020), which formed near the Bahamas, is used as a case study. Several cloud overlap schemes and cloud lookup tables (LUTs), along with other important aspects of all-sky IR radiance assimilation, will be evaluated and tested, aiming to improve the impact on Hurricane Sally's simulation and forecast.

Wei Gu↗
NASA NTRS2011

Long-Term Evolution of the Aerosol Debris Cloud Produced by the 2009 Impact on Jupiter

We present a study of the long-term evolution of the cloud of aerosols produced in the atmosphere of Jupiter by the impact of an object on 19 July 2009. The work is based on images obtained during 5 months from the impact to 31 December 2009 taken in visible continuum wavelengths and from 20 July 2009 to 28 May 2010 taken in near-infrared deep hydrogen-methane absorption bands at 2.1-2.3 micron. The impact cloud expanded zonally from approximately 5000 km (July 19) to 225,000 km (29 October, about 180 deg in longitude), remaining meridionally localized within a latitude band from 53.5 deg S to 61.5 deg S planetographic latitude. During the first two months after its formation the site showed heterogeneous structure with 500-1000 km sized embedded spots. Later the reflectivity of the debris field became more homogeneous due to clump mergers. The cloud was mainly dispersed in longitude by the dominant zonal winds and their meridional shear, during the initial stages, localized motions may have been induced by thermal perturbation caused by the impact's energy deposition. The tracking of individual spots within the impact cloud shows that the westward jet at 56.5 deg S latitude increases its eastward velocity with altitude above the tropopause by 5- 10 m/s. The corresponding vertical wind shear is low, about 1 m/s per scale height in agreement with previous thermal wind estimations. We found evidence for discrete localized meridional motions with speeds of 1-2 m/s. Two numerical models are used to simulate the observed cloud dispersion. One is a pure advection of the aerosols by the winds and their shears. The other uses the EPIC code, a nonlinear calculation of the evolution of the potential vorticity field generated by a heat pulse that simulates the impact. Both models reproduce the observed global structure of the cloud and the dominant zonal dispersion of the aerosols, but not the details of the cloud morphology. The reflectivity of the impact cloud decreased exponentially with a characteristic timescale of 15 days; we can explain this behavior with a radiative transfer model of the cloud optical depth coupled to an advection model of the cloud dispersion by the wind shears. The expected sedimentation time in the stratosphere (altitude levels 5-100 mbar) for the small aerosol particles forming the cloud is 45-200 days, thus aerosols were removed vertically over the long term following their zonal dispersion. No evidence of the cloud was detected 10 months after the impact.

Sanchez-Lavega, A.↗