Search NASA⌕ Search

SEARCH · Search NASA

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 685 records · Page 38

Prediction of Mechanical Properties of Polymers With Various Force Fields

The effect of force field type on the predicted elastic properties of a polyimide is examined using a multiscale modeling technique. Molecular Dynamics simulations are used to predict the atomic structure and elastic properties of the polymer by subjecting a representative volume element of the material to bulk and shear finite deformations. The elastic properties of the polyimide are determined using three force fields: AMBER, OPLS-AA, and MM3. The predicted values of Young s modulus and shear modulus of the polyimide are compared with experimental values. The results indicate that the mechanical properties of the polyimide predicted with the OPLS-AA force field most closely matched those from experiment. The results also indicate that while the complexity of the force field does not have a significant effect on the accuracy of predicted properties, small differences in the force constants and the functional form of individual terms in the force fields determine the accuracy of the force field in predicting the elastic properties of the polyimide.

Odegard, Gregory M.↗

Astronaut Risk Levels During Crew Module (CM) Land Landing

The NASA Engineering Safety Center (NESC) is investigating the merits of water and land landings for the crew exploration vehicle (CEV). The merits of these two options are being studied in terms of cost and risk to the astronauts, vehicle, support personnel, and general public. The objective of the present work is to determine the astronaut dynamic response index (DRI), which measures injury risks. Risks are determined for a range of vertical and horizontal landing velocities. A structural model of the crew module (CM) is developed and computational simulations are performed using a transient dynamic simulation analysis code (LS-DYNA) to determine acceleration profiles. Landing acceleration profiles are input in a human factors model that determines astronaut risk levels. Details of the modeling approach, the resulting accelerations, and astronaut risk levels are provided.

Lawrence, Charles↗

First Directly Retrieved Global Distribution of Tropospheric Column Ozone from GOME: Comparison with the GEOS-CHEM Model

We present the first directly retrieved global distribution of tropospheric column ozone from Global Ozone Monitoring Experiment (GOME) ultraviolet measurements during December 1996 to November 1997. The retrievals clearly show signals due to convection, biomass burning, stratospheric influence, pollution, and transport. They are capable of capturing the spatiotemporal evolution of tropospheric column ozone in response to regional or short time-scale events such as the 1997-1998 El Nino event and a 10-20 DU change within a few days. The global distribution of tropospheric column ozone displays the well-known wave-1 pattern in the tropics, nearly zonal bands of enhanced tropospheric column ozone of 36-48 DU at 20degS-30degS during the austral spring and at 25degN-45degN during the boreal spring and summer, low tropospheric column ozone of <30 DU uniformly distributed south of 35 S during all seasons, and relatively high tropospheric column ozone of >33 DU at some northern high-latitudes during the spring. Simulation from a chemical transport model corroborates most of the above structures, with small biases of <+/-5 DU and consistent seasonal cycles in most regions, especially in the southern hemisphere. However, significant positive biases of 5-20 DU occur in some northern tropical and subtropical regions such as the Middle East during summer. Comparison of GOME with monthly-averaged Measurement of Ozone and Water Vapor by Airbus in-service Aircraft (MOZAIC) tropospheric column ozone for these regions usually shows good consistency within 1 a standard deviations and retrieval uncertainties. Some biases can be accounted for by inadequate sensitivity to lower tropospheric ozone, the different spatiotemporal sampling and the spatiotemporal variations in tropospheric column ozone.

Liu, Xiong↗

Evolution of the Tropical Cyclone Integrated Data Exchange And Analysis System (TC-IDEAS)

The Tropical Cyclone Integrated Data Exchange and Analysis System (TC-IDEAS) is being jointly developed by the Jet Propulsion Laboratory (JPL) and the Marshall Space Flight Center (MSFC) as part of NASA's Hurricane Science Research Program. The long-term goal is to create a comprehensive tropical cyclone database of satellite and airborne observations, in-situ measurements and model simulations containing parameters that pertain to the thermodynamic and microphysical structure of the storms; the air-sea interaction processes; and the large-scale environment.

CloudSat↗

Modeling Powered Aerodynamics for the Orion Launch Abort Vehicle Aerodynamic Database

Modeling the aerodynamics of the Orion Launch Abort Vehicle (LAV) has presented many technical challenges to the developers of the Orion aerodynamic database. During a launch abort event, the aerodynamic environment around the LAV is very complex as multiple solid rocket plumes interact with each other and the vehicle. It is further complicated by vehicle separation events such as between the LAV and the launch vehicle stack or between the launch abort tower and the crew module. The aerodynamic database for the LAV was developed mainly from wind tunnel tests involving powered jet simulations of the rocket exhaust plumes, supported by computational fluid dynamic simulations. However, limitations in both methods have made it difficult to properly capture the aerodynamics of the LAV in experimental and numerical simulations. These limitations have also influenced decisions regarding the modeling and structure of the aerodynamic database for the LAV and led to compromises and creative solutions. Two database modeling approaches are presented in this paper (incremental aerodynamics and total aerodynamics), with examples showing strengths and weaknesses of each approach. In addition, the unique problems presented to the database developers by the large data space required for modeling a launch abort event illustrate the complexities of working with multi-dimensional data.

Chan, David T.↗

WFIRST CGI: Prism and Polarizer Characterization Modes

As part of its technology demonstration, the WFIRST coronagraph instrument (CGI) will demonstrate point source spectroscopy and polarization measurements of disks. The specific implementation of spectroscopy is a zero-deviation Amici prism and slit to be placed on the planet after high contrast has been achieved by CGI. The polarization optics are a set of Wollaston prisms so that orthogonal polarization states can be measured simultaneously. The CGI spectral characterization mode, being designed and built and Goddard Space Flight Center (GSFC), has a spectral resolution of R50 and is designed to accommodate a 15% bandpass spanning 610-785nm. In order to recover Stokes information, there are two sets of Wollaston prisms clocked 45 degrees with respect to one another with each measurement taken in series. The Wollaston design and optical elements are a contribution by the Japanese Aerospace Exploration Agency, with final alignment and testing being done at GSFC. The spectroscopy mode is designed to target Methane absorption features around 730nm, keeping the spectral resolution as low as possible to improve the signal-to-noise ratio and hence reduce detection time. We highlight the requirements for these modes and address the challenge of on-orbit spectral calibration for a deployable slit in the presence of pointing drifts. Of unique interest is how the observatory error budget couples into good stellar spectrum calibration and subsequent speckle subtraction. We also provide further detail on the optomechanical design, its stability based on thermal and structural predictions, modeled performance, and operations concept. These performance metrics are simulated to demonstrate how a slit located at an arbitrary field point is homed onto the planet and converted to a calibrated spectrum.

Tyler D Groff↗

Modal Test and Analysis Correlation of Wind Tunnel Blades for the National Transonic Facility

NASA has a significant maintenance challenge with the aging wind tunnel infrastructure across the agency . Wind tunnel blades, such as the ones in the National Transonic Facility at NASA Langley Research Center, have been subjected to damage from impacts and fatigue. Accurate analytical models of these blades are necessary to either recertify replacement blades to be built or for redesign of these blades. In this paper an upgraded finite element model of these wind tunnel blades with the capability of representing damage is presented. An eigenvalue analysis of the structural model is performed, and the boundary conditions are adjusted to simulate the natural frequencies obtained by experimental testing of a set of twenty five blades in 1992. Additional numerical studies of the effect of delamination damage on the natural frequencies of the blade are conducted using the tuned finite element model.

wind tunnel↗

Modal Test and Analysis Correlation of Wind Tunnel Blades for the National Transonic Facility

NASA has a significant maintenance challenge with the aging wind tunnel infrastructure across the agency . Wind tunnel blades, such as the ones in the National Transonic Facility at NASA Langley Research Center, have been subjected to damage from impacts and fatigue. Accurate analytical models of these blades are necessary to either recertify replacement blades to be built or for redesign of these blades. In this paper an upgraded finite element model of these wind tunnel blades with the capability of representing damage is presented. An eigenvalue analysis of the structural model is performed, and the boundary conditions are adjusted to simulate the natural frequencies obtained by experimental testing of a set of twenty five blades in 1992. Additional numerical studies of the effect of delamination damage on the natural frequencies of the blade are conducted using the tuned finite element model.

wind tunnel↗

Effects of Kinetic Processes in Shaping Io's Global Plasma Environment: A 3D Hybrid Model

The global dynamics of the ionized and neutral gases in the environment of Io plays an important role in the interaction of Jupiter s corotating magnetospheric plasma with Io. Stationary simulations of this problem have already been done using the magnetohydrodynamics (MHD) and the electrodynamics approaches. One of the major results of recent simplified two-fluid model simulations [Saur, J., Neubauer, F.M., Strobel, D.F., Summers, M.E., 2002. J. Geophys. Res. 107 (SMP5), 1-18] was the production of the structure of the double-peak in the magnetic field signature of the Io flyby. These could not be explained before by standard MHD models. In this paper, we present a hybrid simulation for Io with kinetic ions and fluid electrons. This method employs a fluid description for electrons and neutrals, whereas for ions a particle approach is used. We also take into account charge-exchange and photoionization processes and solve self-consistently for electric and magnetic fields. Our model may provide a much more accurate description for the ion dynamics than previous approaches and allows us to account for the realistic anisotropic ion velocity distribution that cannot be done in fluid simulations with isotropic temperatures. The first results of such a simulation of the dynamics of ions in Io s environment are discussed in this paper. Comparison with the Galileo IO flyby results shows that this approach provides an accurate physical basis for the interaction and can therefore naturally reproduce all the observed salient features.

Lipatov, Alexander S.↗

Development of an Effective Finite-rate Oxidation Model for NuSil-coated Charred Carbon Preform Ablators

Recently, a detailed effective finite-rate surface chemistry model was developed for the oxidation of FiberForm [1] using the molecular beam experimental data of Poovathingal et al., [2]. FiberForm is the major building block of the thermal protection system (TPS) material Phenolic Impregnated Carbon Ablator (PICA), commonly used by NASA. The surface chemistry model consists of detailed surface reaction mechanisms such as adsorption, desorption, and several types of Langmuir-Hinshelwood (LH) reactions to characterize the oxygen-carbon interactions at the surface. This model provides excellent agreement with the experimental data for oxidation product compositions and corresponding translational energy distributions. Further, an effective oxidation model was constructed that captures the equivalent interaction of oxygen inside the microstructure (including multiple surface collisions) via a modified reactivity for a smooth wall boundary condition. This enables the use of this model directly into Computational Fluid Dynamics (CFD) codes and Material Response (MR) codes to accurately simulate the gas–surface interactions within FiberForm without using the detailed micro-structure. However, this effective model is valid only for the virgin FiberForm. As the TPS material undergoes ablation, the preform carbon burns and turns into char. In addition, PICA is also coated with a protective silicon coating called NuSil for the purpose of mitigating the spread of phenolic dust, and limit contamination during clean room operations. Fig. 1 shows a X-ray microtomography image of a charred NuSil-coated PICA. The NuSil layer introduces species containing silicon into the product mixture. Further, the reactivity of the carbon within the char layer is different from the virgin FiberForm as shown in Fig. 2. Recently a new set of molecular beam experiments were performed on this NuSil-coated TPS material [3] with the same type of oxygen beam used in the previous experiments. Using the latest experimental data, the previously developed effective model will be extended to account for the charred carbon as well as the NuSil coating. Finally, this new effective model with three phases – preform carbon, char, and NuSil; will be compared and validated against the experimental product compositions.

K Swaminathan Gopalan↗

Can We Model Forest Demography Globally? Benchmarking of State-of-the-Art Demographic DGVMs

Forests in Dynamic Global Vegetation Models (DGVMs) have historically been simulated as area-averaged plant functional types in each gridcell instead of representing communities of trees of different sizes and ages (demography). Just as the behaviour of a tree differs according to its ontogeny, so the behaviour of forests is known to differ depending on their demography. Accurately simulating demography is therefore key in order to address questions on afforestation and management strategies, as well as assessments of resilience of forests to disturbances such as drought and fire or diversity changes after a disturbance. Ultimately, demography determines the overall forest biomass in natural forests and is a key arbiter of growth and mortality rates. DGVMs are now able to simulate size and age structure of the trees in forests. However, these models have so far not been benchmarked alongside each other. We evaluate 6 DGVMs (BiomeE, CABLE-POP, FATES, LPJ-GUESS, JULES-RED, ORCHIDEE) against observations on regrowth dynamics as well as natural forests at boreal, temperate and tropical sites. We examine whether the models capture observed regrowth dynamics after disturbance, well-known stand size structure and well-established processes such as self-thinning. We outline the planned route forward towards a standardised international benchmarking framework for demographic DGVMs.

Dynamic Global Vegetation Models↗

A Study of the Kinematic and Dynamic Processes associated with Mesoscale Snowbands in a Midwest United States Snowstorm on 5 February 2020

Mesoscale snowbands are frequently responsible for dumping large amounts of snow over a relatively small region and pose a considerable forecasting challenge. The present study focuses on the mesoscale snowbands observed on 5 February 2020 over the Midwest United States. In order to improve the forecasting of mesoscale snowbands, it is important to understand the processes that organize these features. To this end, we are investigating the kinematic and dynamic processes that shape the mesoscale snowbands using a combination of high-altitude airborne radar, in-situ aircraft measurements, surface observations, and model simulations. This presentation will primarily focus on the three-dimensional structure of the mesoscale snowbands in terms of EXRAD reflectivity and horizontal winds. The horizontal winds are retrieved from the EXRAD scanning-beam Doppler velocities using a VAD technique. These winds, in combination with virtual potential temperature taken from the HRRR operational analysis, are used to compute the bulk Richardson number in order to investigate any contributions to banding features that might be tied to Kelvin-Helmholtz instability. Our analysis does not support the presence of Kelvin-Helmholtz instability. That said, low, but not subcritical, Richardson numbers were found near the top of some banding features in this event and it is possible that the model analysis is smoothing out the virtual potential temperature fields. Finally, we present arguments for using a flight module designed to sample the along-band structures with particular emphasis on observing generating cells and low-level microphysics.

Charles N. Helms↗

Accommodating a hexagonal ζ- phase Mn 2 N film on a cubic MgO (001) substrate

The thin films of c -plane dominated hexagonal ζ- phase Mn 2 N were successfully grown on a cubic MgO (001) substrate directly using plasma-assisted molecular beam epitaxy. The surface was comprehensively studied through experimental and theoretical approaches. Reflection high energy electron diffraction revealed two pseudo-cubic domains along [100] MgO and [110] MgO , and one hexagonal domain, which is 30° apart from the [100] MgO direction. Scanning tunneling microscopy was used to image and resolve the high-quality surface, revealing a distorted hexagonal surface structure. Furthermore, atomic resolution of the hexagonal domain with a 2 × 2 reconstructed surface is presented. For c-plane surface, theoretical investigations were carried out using first principles studies to determine the surface formation energy for various reconstructed surfaces. The theoretical study shows that the nitrogen terminated 2 × 2 structure with a manganese adatom on the surface as the most stable reconstruction that reproduces the experimental results. A corresponding simulated scanning tunneling microscopy model is also presented, providing strong support for the experimentally observed in-plane lattice structures. Furthermore, the manganese: nitrogen stoichiometry within the bulk and surface is in good agreement with the expected ratio of 2:1.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Parallelization and Visual Analysis of Multidimensional Fields: Application to Ozone Production, Destruction, and Transport in Three Dimensions

The three-dimensional, spectral transport model used in the current project was first successfully integrated over climatological time scales by Dr. Guang Ping Lou for the simulation of atmospheric N2O using the United Kingdom Meteorological Office (UKMO) 4-dimensional, assimilated wind and temperature data set. A non-parallel, FORTRAN version of this integration using a fairly simple N2O chemistry package containing only photo-chemical reactions was used to verify our initial parallel model results. The integrations reproduced the gross features of the observed stratospheric climatological N2O distributions but also simulated the structure of the stratospheric Antarctic vortex and its evolution. Subsequently, Dr. Thomas Kindler, who produced much of the parallel version of our model, enlarged the N2O model chemistry package to include N2O reactions involving O(D-1) and also introduced assimilated wind data from NASA as well as UKMO. Initially, transport calculations without chemistry were run using Carbon-14 as a non-reactive tracer gas with the result that large differences in the transport properties of the two assimilated wind data sets were apparent from the resultant Carbon-14 distributions. Subsequent calculations for N2O, including its chemistry, with the two input winds data sets with verification from UARS satellite observations have refined the transport differences between the two such that the model's steering capabilities could be used to infer the correct climatological vertical velocity fields required to support the N2O observations. During this process, it was also discovered that both the NASA and the UKMO data contained spurious values in some of the higher frequency wave components, leading to incorrect local transport calculations and ultimately affecting the large scale properties of the model's N2O distributions, particularly at tropical latitudes. Subsequent model runs with wind data that had been filtered to remove some of the high frequency components produced much more realistic N2O distributions. During the past few months, the UKMO wind data base for a complete two-year period was processed into spectral form for model use. This new version of the input transport data base now includes complete temperature fields as well as the necessary wind data. This was done to facilitate advanced chemical calculations in the parallel model which often depend upon temperature. Additional UKMO data is being added as it becomes available.

Schwan, Karsten↗

Evaluating Effectiveness of Modeling Motion System Feedback in the Enhanced Hess Structural Model of the Human Operator

In order to use the Hess Structural Model to predict the need for certain cueing systems, George and Cardullo significantly expanded it by adding motion feedback to the model and incorporating models of the motion system dynamics, motion cueing algorithm and a vestibular system. This paper proposes a methodology to evaluate effectiveness of these innovations by performing a comparison analysis of the model performance with and without the expanded motion feedback. The proposed methodology is composed of two stages. The first stage involves fine-tuning parameters of the original Hess structural model in order to match the actual control behavior recorded during the experiments at NASA Visual Motion Simulator (VMS) facility. The parameter tuning procedure utilizes a new automated parameter identification technique, which was developed at the Man-Machine Systems Lab at SUNY Binghamton. In the second stage of the proposed methodology, an expanded motion feedback is added to the structural model. The resulting performance of the model is then compared to that of the original one. As proposed by Hess, metrics to evaluate the performance of the models include comparison against the crossover models standards imposed on the crossover frequency and phase margin of the overall man-machine system. Preliminary results indicate the advantage of having the model of the motion system and motion cueing incorporated into the model of the human operator. It is also demonstrated that the crossover frequency and the phase margin of the expanded model are well within the limits imposed by the crossover model.

Zaychik, Kirill↗

A Tropospheric Emission Spectrometer HDO/H2O Retrieval Simulator for Climate Models

Retrievals of the isotopic composition of water vapor from the Aura Tropospheric Emission Spectrometer (TES) have unique value in constraining moist processes in climate models. Accurate comparison between simulated and retrieved values requires that model profiles that would be poorly retrieved are excluded, and that an instrument operator be applied to the remaining profiles. Typically, this is done by sampling model output at satellite measurement points and using the quality flags and averaging kernels from individual retrievals at specific places and times. This approach is not reliable when the model meteorological conditions influencing retrieval sensitivity are different from those observed by the instrument at short time scales, which will be the case for free-running climate simulations. In this study, we describe an alternative, categorical approach to applying the instrument operator, implemented within the NASA GISS ModelE general circulation model. Retrieval quality and averaging kernel structure are predicted empirically from model conditions, rather than obtained from collocated satellite observations. This approach can be used for arbitrary model configurations, and requires no agreement between satellite-retrieved and model meteorology at short time scales. To test this approach, nudged simUlations were conducted using both the retrieval-based and categorical operators. Cloud cover, surface temperature and free-tropospheric moisture content were the most important predictors of retrieval quality and averaging kernel structure. There was good agreement between the D fields after applying the retrieval-based and more detailed categorical operators, with increases of up to 30 over the ocean and decreases of up to 40 over land relative to the raw model fields. The categorical operator performed better over the ocean than over land, and requires further refinement for use outside of the tropics. After applying the TES operator, ModelE had D biases of 8 over ocean and 34 over land compared to TES D, which were less than the biases using raw model D fields.

Field, R. D.↗

Comparing computational times for simulations when using PBPK model template and stand-alone implementations of PBPK models

Introduction We previously developed a PBPK model template that consists of a single model “superstructure” with equations and logic found in many physiologically based pharmacokinetic (PBPK) models. Using the template, one can implement PBPK models with different combinations of structures and features. Methods To identify factors that influence computational time required for PBPK model simulations, we conducted timing experiments using various implementations of PBPK models for dichloromethane and chloroform, including template and stand-alone implementations, and simulating four different exposure scenarios. For each experiment, we measured the required computational time and evaluated the impacts of including various model features (e.g., number of output variables calculated) and incorporating various design choices (e.g., different methods for estimating blood concentrations). Results We observed that model implementations that treat body weight and dependent quantities as constant (fixed) parameters can result in a 30% time savings compared with options that treat body weight and dependent quantities as time-varying. We also observed that decreasing the number of state variables by 36% in our PBPK model template led to a decrease of 20–35% in computational time. Other factors, such as the number of output variables, the method for implementing conditional statements, and the method for estimating blood concentrations, did not have large impacts on simulation time. In general, simulations with PBPK model template implementations of models required more time than simulations with stand-alone implementations, but the flexibility and (human) time savings in preparing and reviewing a model implemented using the PBPK model template may justify the increases in computational time requirements. Conclusion Our findings concerning how PBPK model design and implementation decisions impact computational speed can benefit anyone seeking to develop, improve, or apply a PBPK model, with or without the PBPK model template.

Bernstein, Amanda S.↗

Hybrid RANS/LES Simulation of Jet Surface Interaction Noise

NASA’s Commercial Supersonic Technology (CST) project has formulated a technical challenge to design a quiet propulsion system for a low boom supersonic aircraft that meets Federal Aviation Authority’s airport noise regulations with sufficient margin. Several proposed configurations take advantage of shielding from the wing or other air-frame components. Development of carefully validated computational tools are necessary for critically evaluating installation concepts that are currently being proposed to meet the technical challenge. Semi-empirical models that predict the noise reduction potential of arbitrary shielding surfaces are yet to mature. Another key challenge is the systematic assessment of additional noise from the interaction between high speed jet turbulence and a surface in it’s vicinity. As a first step towards predicting noise reduction due to radical installation concepts from first principles, we simulate the noise generated by a high speed turbulent round jet near a simple planar surface. Detailed comparisons are made with a dedicated experiment conducted at NASA’s Glenn Research Center. Sensitivity of far-field noise predictions to grid resolution is systematically documented. A permeable Ffowcs Williams Hawkings (FWH) surface enclosing both the jet and the shielding surface is used to predict far-field noise from the simulated flowfield. Details of the structured overset grids, numerical discretization, and turbulence model are provided. Near-field comparisons to PIV data and far-field comparisons to microphone array measurements are discussed. Excellent agreement for an initial validation study on an isolated free round jet was obtained and the findings were utilized in the jet surface interaction study. The split between shielded and reflected side of the microphone array was captured with good agreement, as well as the peak in the noise spectra due to scattering of turbulent energy into sound by the trailing edge of the surface.

Stich, Gerrit-Daniel↗