Search NASA⌕ Search

SEARCH · Search NASA

Results for “Multiscale 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 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

3D Phase-Field Simulations of Pattern Formation During Freeze Casting

We present the results of a combined experimental and phase-field modeling study of pattern formation during freeze casting. Those studies use unidirectional freezing of simple binary liquid mixtures of water and sugars (sucrose and trehalose) in a temperature gradient, which suffice to produce hierarchical templated structures similar to those observed in more complex multi-component freeze-cast systems including lamellae, undulated ridges, and more exotic “jellyfish-like” substructures. Multiscale 3D phase-field simulations reproduce remarkably well those structures quantitatively and identify key properties of the ice-water interface that control their formation. They further reveal that lamellae form as a result of a novel symmetry-breaking secondary instability of partially faceted cellular structures and pinpoint additional secondary instability mechanisms giving rise to smaller-scale substructures.

Kaihua Ji↗
NASA NTRS2016

Agile Multi-Scale Decompositions for Automatic Image Registration

In recent works, the first and third authors developed an automatic image registration algorithm based on a multiscale hybrid image decomposition with anisotropic shearlets and isotropic wavelets. This prototype showed strong performance, improving robustness over registration with wavelets alone. However, this method imposed a strict hierarchy on the order in which shearlet and wavelet features were used in the registration process, and also involved an unintegrated mixture of MATLAB and C code. In this paper, we introduce a more agile model for generating features, in which a flexible and user-guided mix of shearlet and wavelet features are computed. Compared to the previous prototype, this method introduces a flexibility to the order in which shearlet and wavelet features are used in the registration process. Moreover, the present algorithm is now fully coded in C, making it more efficient and portable than the MATLAB and C prototype. We demonstrate the versatility and computational efficiency of this approach by performing registration experiments with the fully-integrated C algorithm. In particular, meaningful timing studies can now be performed, to give a concrete analysis of the computational costs of the flexible feature extraction. Examples of synthetically warped and real multi-modal images are analyzed.

Science Data Processing↗
NASA NTRS2016

Agile Multi-Scale Decompositions for Automatic Image Registration

In recent works, the first and third authors developed an automatic image registration algorithm based on a multiscale hybrid image decomposition with anisotropic shearlets and isotropic wavelets. This prototype showed strong performance, improving robustness over registration with wavelets alone. However, this method imposed a strict hierarchy on the order in which shearlet and wavelet features were used in the registration process, and also involved an unintegrated mixture of MATLAB and C code. In this paper, we introduce a more agile model for generating features, in which a flexible and user-guided mix of shearlet and wavelet features are computed. Compared to the previous prototype, this method introduces a flexibility to the order in which shearlet and wavelet features are used in the registration process. Moreover, the present algorithm is now fully coded in C, making it more efficient and portable than the MATLAB and C prototype. We demonstrate the versatility and computational efficiency of this approach by performing registration experiments with the fully-integrated C algorithm. In particular, meaningful timing studies can now be performed, to give a concrete analysis of the computational costs of the flexible feature extraction. Examples of synthetically warped and real multi-modal images are analyzed.

Science Data Processing↗
NASA NTRSDate not supplied

A Convolutional Neural Network for Enhancement of Multi-Scale Localization in Granular Metallic Representative Unit Cells

A convolutional neural network was used to enhance the localization of strain and stress for a generalized method of cells model of a metallic microstructure. Enhanced shear strains, measured in terms of the linear regression coefficients as a function of ground truth strains, were improved from inaccurate and uncorrelated (slope=0.003, Rsq=0.000) to accurate and well correlated (slope=0.890, Rsq=0.882) relative to ground truth (slope=1.0, Rsq=1.0). In applying the convolutional neural network, a convolutional stride of 1.0 (padding=’same’) was only modestly effective while strides of 2 or 3 were more effective yet at higher cost. Additional convolutional layers were generally more expensive than additional dense layers, often with limited benefit. The accuracy of enhanced localized shear strains and stress is expected to yield benefits for damage progression models, especially in the context of hierarchical multi-scale methods where the generalized method of cells is applied at the intermediate scale.

Machine Learning↗
NASA NTRS2018

Direct Observation of Electron Distributions Inside Millisecond Duration Electron Holes

Despite the importance of millisecond duration spatial structures [chorus wave nonlinearity or time domain structures (TDS)] to plasma dynamics, there have been no direct observations of the generation and interaction of these waves and TDS with electrons at the millisecond timescale required for their understanding. Through superposition of 0.195 ms Magnetospheric Multiscale Satellite electron measurements inside 37 superposed, millisecond duration electron holes, the first observations of electron spectra and pitch angle distributions on a submillisecond timescale have been obtained. They show that keV electrons inside the superposed electron hole are accelerated by several hundred volts and that the spectrum of electrons inside the electron hole contain several maxima and minima that are explained by a model of electron energy changes on entering the holes. We report the first observation of trapped electrons inside the TDS, in accordance with the theoretical requirement that such electrons must be present in order to form the phase space holes. Mechanisms of electron acceleration by electron holes (through perpendicular energy gain as the TDS moves into a converging magnetic field) and scattering (due to the perpendicular electric field) are discussed.

Mozer, F. S.↗
NASA NTRS2018

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↗
NASA NTRS2018

Weibull-Based Stochastic Simulation of 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 (Micromechanics Analysis Code/Generalized Method of Cells) 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 1300 degrees Centigrade 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↗
NASA NTRS2022

Design of synthetic collagens that assemble into supramolecular banded fibers as a functional biomaterial testbed

Collagens are the most abundant proteins of the extracellular matrix, and the hierarchical folding and supramolecular assembly of collagens into banded fibers is essential for mediating cell-matrix interactions and tissue mechanics. Collagen extracted from animal tissues is a valuable commodity, but suffers from safety and purity issues, limiting its biomaterials applications. Synthetic collagen biomaterials could address these issues, but their construction requires molecular-level control of folding and supramolecular assembly into ordered banded fibers, comparable to those of natural collagens. Here, we show an innovative class of banded fiber-forming synthetic collagens that recapitulate the morphology and some biological properties of natural collagens. The synthetic collagens comprise a functional-driver module that is flanked by adhesive modules that effectively promote their supramolecular assembly. Multiscale simulations support a plausible molecular-level mechanism of supramolecular assembly, allowing precise design of banded fiber morphology. We also experimentally demonstrate that synthetic fibers stimulate osteoblast differentiation at levels comparable to natural collagen. This work thus deepens understanding of collagen biology and disease by providing a ready source of safe, functional biomaterials that bridge the current gap between the simplicity of peptide biophysical models and the complexity of in vivo animal systems.

Jinyuan Hu↗
NASA NTRS2009

Multiscale Computer Simulation of Tensile and Compressive Strain in Polymer- Coated Silica Aerogels

While the low thermal conductivities of silica aerogels have made them of interest to the aerospace community as lightweight thermal insulation, the application of conformal polymer coatings to these gels increases their strength significantly, making them potentially useful as structural materials as well. In this work we perform multiscale computer simulations to investigate the tensile and compressive strain behavior of silica and polymer-coated silica aerogels. Aerogels are made up of clusters of interconnected particles of amorphous silica of less than bulk density. We simulate gel nanostructure using a Diffusion Limited Cluster Aggregation (DLCA) procedure, which produces aggregates that exhibit fractal dimensions similar to those observed in real aerogels. We have previously found that model gels obtained via DLCA exhibited stress-strain curves characteristic of the experimentally observed brittle failure. However, the strain energetics near the expected point of failure were not consistent with such failure. This shortcoming may be due to the fact that the DLCA process produces model gels that are lacking in closed-loop substructures, compared with real gels. Our model gels therefore contain an excess of dangling strands, which tend to unravel under tensile strain, producing non-brittle failure. To address this problem, we have incorporated a modification to the DLCA algorithm that specifically produces closed loops in the model gels. We obtain the strain energetics of interparticle connections via atomistic molecular statics, and abstract the collective energy of the atomic bonds into a Morse potential scaled to describe gel particle interactions. Polymer coatings are similarly described. We apply repeated small uniaxial strains to DLCA clusters, and allow relaxation of the center eighty percent of the cluster between strains. The simulations produce energetics and stress-strain curves for looped and nonlooped clusters, for a variety of densities and interaction parameters.

Good, Brian↗
NASA NTRS2011

Monte Carlo Simulations of the Formation Flying Dynamics for the Magnetospheric Multiscale (MMS) Mission

The MMS mission is an ambitious space physics mission that will fly 4 spacecraft in a tetrahedron formation in a series of highly elliptical orbits in order to study magnetic reconnection in the Earth's magnetosphere. The mission design is comprised of a combination of deterministic orbit adjust and random maintenance maneuvers distributed over the 2.5 year mission life. Formal verification of the requirements is achieved by analysis through the use of the End-to-End (ETE) code, which is a modular simulation of the maneuver operations over the entire mission duration. Error models for navigation accuracy (knowledge) and maneuver execution (control) are incorporated to realistically simulate the possible maneuver scenarios that might be realized These error models, coupled with the complex formation flying physics, lead to non-trivial effects that must be taken into account by the ETE automation. Using the ETE code, the MMS Flight Dynamics team was able to demonstrate that the current mission design satisfies the mission requirements.

Schiff, Conrad↗
NASA NTRS2010

Multiscale Simulations of Magnetic Island Coalescence

We describe a new interactive parallel Adaptive Mesh Refinement (AMR) framework written in the Python programming language. This new framework, PyAMR, hides the details of parallel AMR data structures and algorithms (e.g., domain decomposition, grid partition, and inter-process communication), allowing the user to focus on the development of algorithms for advancing the solution of a systems of partial differential equations on a single uniform mesh. We demonstrate the use of PyAMR by simulating the pairwise coalescence of magnetic islands using the resistive Hall MHD equations. Techniques for coupling different physics models on different levels of the AMR grid hierarchy are discussed.

Dorelli, John C.↗
NASA NTRSDate not supplied

Multiscale Simulation of Deployable Composite Structures

In this paper, a multiscale simulation method for analyzing deployable composite structures is presented. Effective shell properties of the composites are obtained based on Mechanics of Structure Genome (MSG) homogenization, and then implemented into a user-subroutine UGENS for structural simulation with shell elements in Abaqus. The column bending test(CBT) of a flat thin flexure and lenticular composite boom in a simplified deployer structure are studied for demonstration. A viscoelastic material model with direct integration is adopted in this paper. The CBT simulation shows good agreement with experiments during relaxation, while errors are observed when comparing residual deformation. It is shown that this CBT model can be calibrated to CBT test results. For the lenticular boom analysis, the complete process of flattening, coiling, stowage, deployment and recovery is simulated with the viscoelastic shell model.

Finite element analysis↗
NASA NTRS1994

Canopy reflectance modelling of semiarid vegetation

Three different types of remote sensing algorithms for estimating vegetation amount and other land surface biophysical parameters were tested for semiarid environments. These included statistical linear models, the Li-Strahler geometric-optical canopy model, and linear spectral mixture analysis. The two study areas were the National Science Foundation's Jornada Long Term Ecological Research site near Las Cruces, NM, in the northern Chihuahuan desert, and the HAPEX-Sahel site near Niamey, Niger, in West Africa, comprising semiarid rangeland and subtropical crop land. The statistical approach (simple and multiple regression) resulted in high correlations between SPOT satellite spectral reflectance and shrub and grass cover, although these correlations varied with the spatial scale of aggregation of the measurements. The Li-Strahler model produced estimated of shrub size and density for both study sites with large standard errors. In the Jornada, the estimates were accurate enough to be useful for characterizing structural differences among three shrub strata. In Niger, the range of shrub cover and size in short-fallow shrublands is so low that the necessity of spatially distributed estimation of shrub size and density is questionable. Spectral mixture analysis of multiscale, multitemporal, multispectral radiometer data and imagery for Niger showed a positive relationship between fractions of spectral endmembers and surface parameters of interest including soil cover, vegetation cover, and leaf area index.

Franklin, Janet↗
NASA NTRS2020

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.↗
NASA NTRS2008

Multiscale Computer Simulation of Failure in Aerogels

Aerogels have been of interest to the aerospace community primarily for their thermal properties, notably their low thermal conductivities. While such gels are typically fragile, recent advances in the application of conformal polymer layers to these gels has made them potentially useful as lightweight structural materials as well. We have previously performed computer simulations of aerogel thermal conductivity and tensile and compressive failure, with results that are in qualitative, and sometimes quantitative, agreement with experiment. However, recent experiments in our laboratory suggest that gels having similar densities may exhibit substantially different properties. In this work, we extend our original diffusion limited cluster aggregation (DLCA) model for gel structure to incorporate additional variation in DLCA simulation parameters, with the aim of producing DLCA clusters of similar densities that nevertheless have different fractal dimension and secondary particle coordination. We perform particle statics simulations of gel strain on these clusters, and consider the effects of differing DLCA simulation conditions, and the resultant differences in fractal dimension and coordination, on gel strain properties.

Good, Brian S.↗
NASA NTRS1994

Characterizing multiscale variability of zero intermittency in spatial rainfall

In this paper the authors study how zero intermittency in spatial rainfall, as described by the fraction of area covered by rainfall, changes with spatial scale of rainfall measurement or representation. A statistical measure of intermittency that describes the size distribution of 'voids' (nonrainy areas imbedded inside rainy areas) as a function of scale is also introduced. Morphological algorithms are proposed for reconstructing rainfall intermittency at fine scales given the intermittency at coarser scales. These algorithms are envisioned to be useful in hydroclimatological studies where the rainfall spatial variability at the subgrid scale needs to be reconstructed from the results of synoptic- or mesoscale meteorological numerical models. The developed methodologies are demsonstrated and tested using data from a severe springtime midlatitude squall line and a mild midlatitude winter storm monitored by a meteorological radar in Norman, Oklahoma.

Kumar, Praveen↗
NASA NTRS2021

MMS Observations of Energized He+ Pickup Ions at Quasiperpendicular Shocks

Accelerated He+ PUIs, downstream of quasi-perpendicular shocks, are studied as 10a function of the fast mode Mach number (Mf) and shock obliquity (𝜃Bn).We analyze quasi-perpendicular shocks with Mach numbers in the range [1, 7] observed by the Magnetospheric MultiScale (MMS) mission, and compare upstream and downstream He+ velocity distribution functions (VDFs). For each shock event, we characterize the upstream PUI distribution and derive reduced 1D velocity distributions for the selected upstream and downstream intervals. We also compare the upstream to downstream ratio of spectral indices, computed from the He+ perpendicular distributions, to Mf and 𝜃Bn. We find a positive correlation between this spectral index ratio and Mf, which suggests that perpendicular energization of He+ PUIs is enhanced as the shock becomes stronger. These results inform modeling efforts of PUIs and shock acceleration processes, particularly those taking place at the termination shock.

Michael J. Starkey↗
NASA NTRS1998

Arctic Sea Ice Motion from Wavelet Analysis of Satellite Data

Wavelet analysis of DMSP SSM/I (Special Sensor Microwave/Imager) 85 GHz and 37 GHz radiance data, SMMR (Scanning Multichannel Microwave Radiometer) 37 GHz, and NSCAT (NASA Scatterometer) 13.9 GHZ data can be used to obtain daily sea ice drift information for both the northern and southern polar regions. The derived maps of sea ice drift provide both improved spatial coverage over the existing array of Arctic Ocean buoys and better temporal resolution over techniques utilizing data from satellite synthetic aperture radars (SAR). Examples of derived ice-drift maps in the Arctic illustrate large-scale circulation reversals within a period of a couple weeks. Comparisons with ice displacements derived from buoys show good quantitative agreement. NASA Scatterometer (NSCAT) 13.9 GHZ data have been also used for wavelet analysis to derive sea-ice drift. First, the 40' incidence-angle, sigma-zero (surface roughness) daily map of whole Arctic region with 25 km of pixel size from satellite's 600 km swath has been constructed. Then, the similar wavelet transform procedure to SSM/I data can be applied. Various scales of wavelet transform and threshold have been tested. By overlaying , neighbor filtering, and block-averaging the results of multiscale wavelet transforms, the final sea ice drift vectors are much smooth and representative to the sea ice motion. This wavelet analysis procedure is robust and can make a major contribution to the understanding of ice motion over large areas at relatively high temporal resolutions. The results of wavelet analysis of SSM/I and NSCAT images and buoy data can be merged by some data fusion techniques and will help to improve our current knowledge of sea ice drift and related processes through the data assimilation of ocean-ice numerical model.

Liu, Antony K.↗