Search NASA⌕ Search

SEARCH · Search NASA

Results for “shape modeling”

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 505 records · Page 28

Effects of Subscale Size and Shape on Global Energy Dissipation in a Multiscale Model of a Fiber-Reinforced Composite Exhibiting Post-Peak Strain Softening Using Abaqus and FEAMAC

A mesh objective crack band model is implemented in the generalized method of cells (GMC) micromechanics model to predict failure of a composite repeating unit cell (RUC). The micromechanics calculations are achieved using the MAC/GMC core engine within the ImMAC suite of micromechanics codes, developed at the NASA Glenn Research Center. The microscale RUC is linked to a macroscale Abaqus/Standard finite element model using the FEAMAC multiscale framework (included in the ImMAC suite). The effects of the relationship between the characteristic length of the finite element and the size of the microscale RUC on the total energy dissipation of the multiscale model are investigated. A simple 2-D composite square subjected to uniaxial tension is used to demonstrate the effects of scaling the dimensions of the RUC such that the length of the sides of the RUC are equal to the characteristic length of the finite element. These results are compared to simulations where the size of the RUC is fixed, independent of the element size. Simulations are carried out for a variety of mesh densities and element shapes, including square and triangular. Results indicate that a consistent size and shape must be used to yield preserve energy dissipation across the scales.

Pineda, Evan, J.↗

Wake Velocimetry of a Sphere-Cone Model in a Mach 10 Air Freestream using FLEET

Velocity measurements in the wake behind a sphere-cone-shaped vehicle wind tunnel model in a Mach 10 hypersonic air flow using femtosecond laser electronic excitation tagging (FLEET) are reported in this work. The FLEET measurements were performed at 1 kHz using a femtosecond laser centered at 800 nm and an ultrahigh-speed intensified camera system. The results reported here are for a single freestream condition corresponding to approximately Re ∞ /L = 3.6∙10 6 m -1 . The performance of the velocimetry instrument was evaluated in the quiescent test section at conditions relevant to the wake measurements. For velocity measurements in the wake, the FLEET measurement volume was translated to several radial locations from the model centerline in the vertical direction and at a streamwise location corresponding to approximately 12 mm from the payload. Reverse flow with mean velocities ranging from -12 m/s to -48 m/s were observed in the recirculating region of the wake, whereas hypersonic velocities in the range of 1079 m/s to 1183 m/s were observed in the free shear layer. A velocity uncertainty methodology is also outlined and applied for the reported measurements. The velocity data reported in this work is expected to play a significant role in the validation of computational codes modeling the hard-to-predict separated hypersonic wake.

FLEET↗

Wake Velocimetry of a Sphere-Cone Model in a Mach 10 Air Freestream using FLEET

Velocity measurements in the wake behind a sphere-cone-shaped vehicle wind tunnel model in a Mach 10 hypersonic air flow using femtosecond laser electronic excitation tagging (FLEET) are reported in this work. The FLEET measurements were performed at 1 kHz using a femtosecond laser centered at 800 nm and an ultrahigh-speed intensified camera system. The results reported here are for a single freestream condition corresponding to approximately Re ∞ /L = 3.6∙10 6 m -1 . The performance of the velocimetry instrument was evaluated in the quiescent test section at conditions relevant to the wake measurements. For velocity measurements in the wake, the FLEET measurement volume was translated to several radial locations from the model centerline in the vertical direction and at a streamwise location corresponding to approximately 12 mm from the payload. Reverse flow with mean velocities ranging from -12 m/s to -48 m/s were observed in the recirculating region of the wake, whereas hypersonic velocities in the range of 1079 m/s to 1183 m/s were observed in the free shear layer. A velocity uncertainty methodology is also outlined and applied for the reported measurements. The velocity data reported in this work is expected to play a significant role in the validation of computational codes modeling the hard-to-predict separated hypersonic wake.

FLEET↗

Simulation of Stochastic Mud-Crack Damage Formation in an Environmental Barrier Coating

The integrated Finite Element Analysis–Micromechanics Analysis Code/Ceramics Analysis and Reliability Evaluation of Structures (FEAMAC/CARES) program was used to simulate the formation of mudflat-cracks from thermomechanical loading on a multi-layered Environmental Barrier Coating (EBC) system deposited on a ceramic substrate. FEAMAC/CARES combines MAC/GMC multiscale composite micromechanics code with CARES/Life probabilistic multiaxial failure criteria code and Abaqusfinite element analysis. In this work, step function elastic modulus reduction of randomly damaged finite elements was used to represent discrete cracking events. The use of many small-sized low-aspect-ratio finite elements enabled the depiction of crack boundaries and formation of mudflat patterned damage. Demonstrated examples include finite element models of button–sized disk–shaped 3-D specimen, and a 2-D model of through-the-thickness cross-section. All models were subjected to a progressive cool down from 1300oC to room temperature. Mudflat crack damage in the coating system resulted from the buildup of residual tensile stresses between the individual material constituents from thermal expansion mismatch. A 2-parameter Weibull distribution characterized the coating layer stochastic strength response and the effect of the Weibull modulus on the formation of damage was studied here.

residual tensile stress↗

Stochastic Simulation of Mudcrack Damage Formation in an Environmental Barrier Coating

The FEAMAC/CARES program, which integrates finite element analysis (FEA) with the MAC/GMC (Micromechanics Analysis Code with Generalized Method of Cells) and the CARES/Life (Ceramics Analysis and Reliability Evaluation of Structures / Life Prediction) programs, was used to simulate the formation of mudcracks during the cooling of a multilayered environmental barrier coating (EBC) deposited on a silicon carbide substrate. FEAMAC/CARES combines the MAC/GMC multiscale micromechanics analysis capability (primarily developed for composite materials) with the CARES/Life probabilistic multiaxial failure criteria (developed for brittle ceramic materials) and Abaqus (Dassault Systèmes) FEA. In this report, elastic modulus reduction of randomly damaged finite elements was used to represent discrete cracking events. The use of many small-sized low-aspect-ratio elements enabled the formation of crack boundaries, leading to development of mudcrack-patterned damage. Finite element models of a disk-shaped three-dimensional specimen and a twodimensional model of a through-the-thickness cross section subjected to progressive cooling from 1,300 °C to an ambient temperature of 23 °C were made. Mudcrack damage in the coating resulted from the buildup of residual tensile stresses between the individual material constituents because of thermal expansion mismatches between coating layers and the substrate. A two-parameter Weibull distribution characterized the coating layer stochastic strength response and allowed the effect of the Weibull modulus on the formation of damage and crack segmentation lengths to be studied. The spontaneous initiation of cracking and crack coalescence resulted in progressively smaller mudcrack cells as cooling progressed, consistent with a fractal-behaved fracture pattern. Other failure modes such as delamination, and possibly spallation, could also be reproduced. The physical basis assumed and the heuristic approach employed, which involves a simple stochastic cellular automaton methodology to approximate the crack growth process, are described. The results ultimately show that a selforganizing mudcrack formation can derive from a Weibull distribution that is used to describe the stochastic strength response of the bulk brittle ceramic material layers of an EBC.

Nemeth, Noel N.↗

Effects of Nose Shape and Spray Control Strips on Emergence and Planing Spray of Hydro-ski Models

The emergence and planing-spray characteristics of flat-bottom surfaces representative of hydro-skis with various bow shapes and deflectors were investigated for trims of 12° and 20° and a speed of 30 feet per second. The emergence-spray characteristics of the models with various bow shapes were investigated for depths of submersion of the bow from 0.50 inch to -0.50 inch. The planing-spray characteristics of the models with the various deflectors were determined for a draft corresponding to a length-beam ratio of 4 . The emergence- and planing-spray patterns are shown in the photographs taken from forward of and above the models and from the side of the models. The most favorable emergence spray was obtained with a bow of triangular plan form and sharp profile. The greatest improvement in planing spray was obtained with vertical chine strips and a transverse barrier strip located forward of the wetted planing area.

John R McGehee↗

Structural Modeling of an Osteoblast Subjected to Hypergravity Loading

Osteoblasts, the bone-forming cells, are sensitive to altered mechanical loads such as those induced by fluid shear forces and cyclic uniaxial tension and compression. The cytoskeleton is thought to play an important role in a cell's response to mechanical loading and to coordinate cell shape changes. While some theoretical calculations indicate that a cell should not be able to respond to relatively small changes in gravity given the magnitude of forces exerted by intracellular components, osteoblasts have been shown to respond to microgravity. In this study, we investigated whether an osteoblast can sense a constant hypergravity acceleration of 10G. We asked if this load changes the shape of an osteoblast, and if engineering modeling accurately predicts the shape change. We modeled a cell by representing it as a series of axisymmetric shell structures. To test this model empirically, we performed experiments with cultured osteoblasts to measure changes in cell shape resulting from hypergravity loading.

Globus, Ruth K.↗

DEM simulation of the compression of crushable sand: does the initial particle shape matter?

Advances in DEM modeling, combined with high-resolution X-ray tomography, opened the way for computer models based on virtual replicas of the particles which preserve nearly all facets of their geometry. This leads to simulation advantages, but also high computational costs. Here we tackle a question stemming from this trend: how accurate should particle models be to ensure accuracy? We address this question for the case of the compression of crushable sand. LS-DEM was used to generate three models of Ottawa sand (exact replicas, ellipsoids, and spheres) from digital images of its grains. Compression-induced crushing was simulated for all sets by tracking evolving size and shape distribution. The results confirm that exact replicas provide the closest match of the measurements. However, intermediate degrees of rendering (e.g. ellipsoids preserving volume and aspect ratio of the real grains) led to satisfactory results only marginally different from those of exact replicas. In conclusion, these findings provide an example of the protocols that may be followed to identify the optimal degree of particle approximation which should be regarded as mandatory to achieve a conscious, sustainable use of computational resources.

58 GEOSCIENCES↗

Independent pixel and Monte Carlo estimates of stratocumulus albedo

Monte Carlo radiative transfer methods are employed here to estimate the plane-parallel albedo bias for marine stratocumulus clouds. This is the bias in estimates of the mesoscale-average albedo, which arises from the assumption that cloud liquid water is uniformly distributed. The authors compare such estimates with those based on a more realistic distribution generated from a fractal model of marine stratocumulus clouds belonging to the class of 'bounded cascade' models. In this model the cloud top and base are fixed, so that all variations in cloud shape are ignored. The model generates random variations in liquid water along a single horizontal direction, forming fractal cloud streets while conserving the total liquid water in the cloud field. The model reproduces the mean, variance, and skewness of the vertically integrated cloud liquid water, as well as its observed wavenumber spectrum, which is approximately a power law. The Monte Carlo method keeps track of the three-dimensional paths solar photons take through the cloud field, using a vectorized implementation of a direct technique. The simplifications in the cloud field studied here allow the computations to be accelerated. The Monte Carlo results are compared to those of the independent pixel approximation, which neglects net horizontal photon transport. Differences between the Monte Carlo and independent pixel estimates of the mesoscale-average albedo are on the order of 1% for conservative scattering, while the plane-parallel bias itself is an order of magnitude larger. As cloud absorption increases, the independent pixel approximation agrees even more closely with the Monte Carlo estimates. This result holds for a wide range of sun angles and aspect ratios. Thus, horizontal photon transport can be safely neglected in estimates of the area-average flux for such cloud models. This result relies on the rapid falloff of the wavenumber spectrum of stratocumulus, which ensures that the smaller-scale variability, where the radiative transfer is more three-dimensional, contributes less to the plane-parallel albedo bias than the larger scales, which are more variable. The lack of significant three-dimensional effects also relies on the assumption of a relatively simple geometry. Even with these assumptions, the independent pixel approximation is accurate only for fluxes averaged over large horizontal areas, many photon mean free paths in diameter, and not for local radiance values, which depend strongly on the interaction between neighboring cloud elements.

Cahalan, Robert F.↗

Neural Networks Used to Compare Designed and Measured Time-Average Patterns

Electronic time-average holograms are convenient for comparing the measured vibration modes of fan blades with those calculated by finite-element models. At the NASA Lewis Research Center, neural networks recently were trained to perform what had been a simple visual comparison of the predictions of the design models with the measurements. Finite-element models were used to train neural networks to recognize damage and strain information encoded in subtle changes in the time-average patterns of cantilevers. But the design-grade finite element models were unable to train the neural networks to detect damage in complex blade shapes. The design-model-generated patterns simply did not agree well enough with the measured patterns. Instead, hybrid-training records, with measured time-average patterns as the input and model-generated strain information as the output, were used to effect successful training.

Decker, Arthur J.↗

Effect of Protuberance Shape and Orientation on Space Shuttle Orbiter Boundary-Layer Transition

This document describes an experimental study conducted to examine the effects of protuberances on hypersonic boundary-layer transition. The experiment was conducted in the Langley 20-Inch Mach 6 Tunnel on a series of 0.9%-scale Shuttle Orbiter models. The data were acquired to complement the existing ground-based boundary-layer transition database that was used to develop Version 1.0 of the boundary-layer transition RTF (return-to-flight) tool. The existing ground-based data were all acquired on 0.75%-scale Orbiter models using diamond-shaped ( pizza-box ) trips. The larger model scale facilitated in manufacturing higher fidelity protuberances. The end use of this experimental database will be to develop a technical basis (in the form of a boundary-layer transition correlation) to assess representative protrusion shapes, e.g., gap fillers and protrusions resulting from possible tile repair concepts. The primary objective of this study is to investigate the effects of protuberance-trip location and geometry on Shuttle Orbiter boundary-layer transition. Secondary goals are to assess the effects of gap-filler orientation and other protrusion shapes on boundary-layer transition. Global heat-transfer images using phosphor thermography of the Orbiter windward surface and the corresponding streamwise and spanwise heating distributions were used to infer the state of the boundary layer, i.e., laminar, transitional, or turbulent.

King, RUdolph A.↗

An empirical model of the Venusian outer environment. I - The shape of the dayside solar wind-atmosphere interface. II - The shape and location of the bow shock

Solar wind plasma and magnetic field data and ionospheric data obtained from the Pioneer Venus orbiter are discussed. The variation in the magnetic field pressure within the magnetic barrier is shown to be similar to that expected for the solar wind pressure variations along an obstacle's boundary when a more realistic approximation of the shape of the ionosphere is included. Simultaneous solar wind pressure, ionospheric pressure, and magnetic barrier pressure data reveal that the ionospheric pressure below the ionopause is approximately equal to the solar wind pressure. The magnetic barrier pressure is found to be equal to approximately 2/3 to 3/4 of both the solar wind pressure and the ionospheric pressure. Estimates of the hot plasma pressure contribution to the total pressure in the magnetic barrier vary from 1/4 to 1/3. Just below the ionopause the ionospheric pressure deviations from the mean ionospheric pressure are significant, especially below approximately 400 km; they appear to be indicative of the adjustment of the ionospheric structure to changeing solar wind conditions. A first-order model of the ionopause pressure variations as a function of height and solar-zenith angle is proposed. Venera 9 and Venera 10 bow shock crossings are then analyzed for solar zenith angles in the range from approximately 25 deg to approximately 153 deg. A comparison of Venera 9 and Venera 10 bow shock crossings with those observed by the Pioneer Venus orbiter suggests that both the secular variation and the latitudinal asymmetry may be responsible for the closer shock crossings of the Venera 9 and Venera 10 spacecraft.

Vaisberg, O. L.↗

On the accuracy of shape sensitivity

The calculation of sensitivity of the response of a structure modeled by finite elements to shape variation is known to be subject to numerical difficulties. The accuracy of a given method is typically measured against the yard stick of finite-difference sensitivity calculation. The present paper demonstrates with a simple example that this approach may be flawed because of discretization errors associated with the finite element mesh. Seven methods for calculating sensitivity derivatives are compared for a two-material beam problem with a moving interface. It is found that as the mesh is refined, displacement sensitivity derivatives converge more slowly than the displacements. Six of the methods agree fairly well, but the adjoint variational surface method provides substantially different results. However, the difference is found to reflect convergence from another direction to the same answer rather than reduced accuracy. Additionally, it is observed that small derivatives are particularly prone to accuracy problems.

Haftka, R. T.↗

Geometric-optical bidirectional reflectance modeling of the discrete crown vegetation canopy - Effect of crown shape and mutual shadowing

Many natural vegetation covers can be regarded as assemblages of plant crowns, on a background, which interact with light as discrete objects. On this basis, geometric optics furnish an approach to the modeling of vegetation canopies that captures the most important features of such growths' bidirectional measurements. Attention is presently given to models which approximate these phenomena and relate the size, shape, and count density of plant crowns to viewing and illumination positions and to crown-background reflectance contrasts.

Li, Xiaowen↗

Optical technique to characterize heavy rain

A technique to measure the size and velocity of raindrops is presented and experimental and numerical results discussed. The technique bases the droplet size on the absolute intensity of the light scattered when the droplets cross a laser-based probe volume. The velocity is obtained from the Doppler frequency. A size range of 0.1 mm to 9 mm was covered with overlapping ranges. Experiments were conducted with simulated rain sprays except that the local LWC was much larger than the values encountered in natural rain. The accuracy of the size measurement was not strongly affected by spray interference, but it showed a clear dependence with the droplet shape. A simple numerical model based on the light reflected from ellipsoids attempts to explain the relationship between droplet shape and absolute scattered light.

Hess, C. F.↗

Quantification of Processing Effects on Filament Wound Pressure Vessels

A computational simulation procedure is described which is designed specifically for the modeling and analysis of filament wound pressure vessels. Cylindrical vessels with spherical or elliptical end caps can be generated automatically. End caps other than spherical or elliptical may be modeled by varying circular sections along the x-axis according to the end cap shape. The finite element model generated is composed of plate type quadrilateral shell elements on the entire vessel surface. This computational procedure can also be used to generate grid, connectivity and material cards (bulk data) for component parts of a larger model. These bulk data are assigned to a user designated file for finite element structural/stress analysis of composite pressure vessels. The procedure accommodates filament wound pressure vessels of all types of shells-of -revolution. It has provisions to readily evaluate initial stresses due to pretension in the winding filaments and residual stresses due to cure temperature.

Aiello, Robert A.↗

A method for the geometrically nonlinear analysis of compressively loaded prismatic composite structures

A method was developed for the geometrically nonlinear analysis of the static response of thin-walled stiffened composite structures loaded in uniaxial or biaxial compression. The method is applicable to arbitrary prismatic configurations composed of linked plate strips, such as stiffened panels and thin-walled columns. The longitudinal ends of the structure are assumed to be simply supported, and geometric shape imperfections can be modeled. The method can predict the nonlinear phenomena of postbuckling strength and imperfection sensitivity which are exhibited by some buckling-dominated structures. The method is computer-based and is semi-analytic in nature, making it computationally economical in comparison to finite element methods. The method uses a perturbation approach based on the use of a series of buckling mode shapes to represent displacement contributions associated with nonlinear response. Displacement contributions which are of second order in the model amplitudes are incorported in addition to the buckling mode shapes. The principle of virtual work is applied using a finite basis of buckling modes, and terms through the third order in the model amplitudes are retained. A set of cubic nonlinear algebraic equations are obtained, from which approximate equilibrium solutions are determined. Buckling mode shapes for the general class of structure are obtained using the VIPASA analysis code within the PASCO stiffened-panel design code. Thus, subject to some additional restrictions in loading and plate anisotropy, structures which can be modeled with respect to buckling behavior by VIPASA can be analyzed with respect to nonlinear response using the new method. Results obtained using the method are compared with both experimental and analytical results in the literature. The configurations investigated include several different unstiffened and blade-stiffening panel configurations, featuring both homogeneous, isotropic materials, and laminated composite material.

Stoll, Frederick↗

Particle Interaction Physics Model Formulation for Plume-Surface Interaction Erosion and Cratering

As part of the Game Changing Development (GCD) Program, funded by NASA’s Space Technology Mission Directorate (STMD), the development of simulation capability for the prediction of extra-terrestrial Plume Surface Interaction (PSI) environments has been undertaken by the Fluid Dynamics Branch at NASA/MSFC. The Predictive Simulation Capability (PSC) Element is focused on creating simulation capability for the reliable and accurate prediction of PSI in Martian (~650 Pa) and Lunar (vacuum) ambient environments. In addition to the predictive simulation capability, the GCD Program also contains a companion Ground Testing Element for development of focused datasets for validation of predictive capability as well as a Flight-focused Instrumentation Element. This paper will present the status of implementing and maturing particle-particle interaction constituent physics models essential in simulating the landing surface granular material flow under PSI effects. This gas-particle multi-phase interaction modeling of plume impingement flow on the extra-terrestrial soil material is performed with the Gas-Granular Flow Solver (GGFS) addressed in a companion paper. The response of regolith particle flow induced by lander PSI requires accurate representation of the regolith granular material fluidic behavior and gas-granular interactions. The lunar regolith, as the extreme example, is poorly sorted with broad particle size distributions and large fines content. It has significant cohesion, due to interlocking particle shapes for the very jagged particles. The combination of particle shape and size distribution has been identified as major drivers in the complex particle flow response and resulting crater shape characteristics of extraterrestrial granular material. Constituent models for spherical particles can be formulated directly from particle kinetics theory. Complex particle shapes can be modeled by gluing together elemental spherical shapes into composite particles, requiring a Discrete Element Model (DEM) particle kinetics modeling approach to extract data and formulate constituent models. Mixture constituent models for poly-disperse mixtures (i.e, containing distribution of particle sizes) have recently been developed. The required non-spherical particle mixture granular material response closure models are then obtained through small-scale unit physics DEM simulations for the range of particle shapes, mixtures and packing densities. The granular material response closure models are then implemented in the Eulerian granular flow formulation. This DEM-based constituent model extraction process and formulation of poly-disperse particle mixtures has been successfully developed by small business and academic partners in the development of the Gas-Granular Flow Solver (GGFS) simulation program simulation framework. The currently implemented capabilities have reached the capability level of modeling bi-disperse, non-spherical particle mixtures is being continuously extended towards computational modeling of full range irregular particle mixtures. Under the GCD project, this technology is being further developed, transferred to NASA analysts, and matured towards application readiness. The predictive simulation capability team under the GCD project has acquired the modeling tools and processes of the DEM based constituent model formulation from the GGFS development team and is developing the capability to replicate the existing process. This is the first important step towards the ability of the NASA team to independently perform such model development in a production setting. Further efforts are underway to migrate the DEM based model simulation process performed with the academic based tools to more capable Open Source, highly parallelized simulation tools for efficient operation on NASA HPC assets. Evaluation of the currently implemented (such as mono-disperse and bi-disperse spherical and irregular shape particle constituent model applications) and continuously evolving full-range particle physics models in the GGFS tool is performed by the NASA team to advance application readiness of the simulations. Application testing for complex PSI erosions and cratering scenarios such as the Apollo LM is performed for axi-symmetric and full 3D simulations to aid the tool developers in achieving practical application readiness for NASA projects. Important validation and application testing will further be performed against experimental data generated under the GCD PSI project experimental component.

Peter A Liever↗