Search NASA⌕ Search

SEARCH · Search NASA

Results for “Microscale 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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
NASA NTRS1982

Fluid flow in crystallization processes

Investigations of the interactions of heat, mass, and momentum transport in crystal growth from the melt are described. Studies of the mall-scale floating zone process and on a prototype of the vertical Bridgman growth system were emphasized. In both systems detailed numerical calculations are used to dissect the interplay between fluid convection and dopant segregation. These calculations are based on finite element techniques that make feasible the complete solution of solidification problems which include convection. Other numerical methods were developed for solving thermal models of crystal growth processes with melt/gas menisci and for simulating the microscale instabilities in solidification interfaces.

Brown, R. A.↗
NASA NTRS1998

A Kinematically Consistent Two-Point Correlation Function

A simple kinematically consistent expression for the longitudinal two-point correlation function related to both the integral length scale and the Taylor microscale is obtained. On the inner scale, in a region of width inversely proportional to the turbulent Reynolds number, the function has the appropriate curvature at the origin. The expression for two-point correlation is related to the nonlinear cascade rate, or dissipation epsilon, a quantity that is carried as part of a typical single-point turbulence closure simulation. Constructing an expression for the two-point correlation whose curvature at the origin is the Taylor microscale incorporates one of the fundamental quantities characterizing turbulence, epsilon, into a model for the two-point correlation function. The integral of the function also gives, as is required, an outer integral length scale of the turbulence independent of viscosity. The proposed expression is obtained by kinematic arguments; the intention is to produce a practically applicable expression in terms of simple elementary functions that allow an analytical evaluation, by asymptotic methods, of diverse functionals relevant to single-point turbulence closures. Using the expression devised an example of the asymptotic method by which functionals of the two-point correlation can be evaluated is given.

Ristorcelli, J. R.↗
NASA NTRS2015

Evaluation of the Anisotropic Radiative Conductivity of a Low-Density Carbon Fiber Material from Realistic Microscale Imaging

The radiative heat transfer inside a low-density carbon fiber insulator is analyzed using a three-dimensional direct simulation model. A robust procedure is presented for the numerical calculation of the geometric configuration factor to compute the radiative energy exchange processes among the small discretized surface areas of the fibrous material. The methodology is applied to a polygonal mesh of a fibrous insulator obtained from three-dimensional microscale imaging of the real material. The anisotropic values of the radiative conductivity are calculated for that geometry. The results yield both directional and thermal dependence of the radiative conductivity.

Evaluation↗
NASA NTRSDate not supplied

Overview of Ablative TPS Modeling at NASA Ames

Over the past decade, NASA has invested in efforts to build predictive thermal protection system (TPS) material models from the micro-scale to the macro-scale. To complement the mission design cycle process and reduce the need for extensive testing, NASA is developing modeling and simulation tools that enable characterizing material properties and response to hot plasma experienced during atmospheric entry. Traditional material response and ablation modeling tools, such as the heritage code FIAT, and its multidimensional siblings, TITAN and 3dFIAT, are being complemented with newly developed software such as Icarus and PATO. Both of these programs are three-dimensional, finite-volume solvers that use unstructured meshes and 21st century programming paradigms to allow for efficient parallel simulations. FIAT and Icarus are also used for TPS sizing purposes. Today, these traditional tools are being supplemented with computational materials models at the atomistic level. The scales of interest range from computational chemistry (Density Functional Theory [DFT]), to atomistic simulations (Molecular Dynamics [MD]), to the microscale with the Porous Microstructure Analysis (PuMA) software that was recently awarded the 2022 NASA Software of the Year award. Finally, thermo-structural modeling is also of interest to the TPS Materials branch and done using commercial tools such as MSC MARC, MENTAT, NASTRAN and PATRAN. The present talk will also link the use of these computational tools to current NASA missions and projects associated with challenging and complex vehicles entries/reentries.

materials modeling↗
NASA NTRS2019

Computational Micromechanics Model for the Analysis of Fiber Kinking in Unidirectional FRP

A computational micromechanics (CMM) model is developed to analyze fiber kinking, which is a failure mechanism that takes place in fiber-reinforced composites when they are loaded under longitudinal compression. The CMM model consists of a single AS4 carbon fiber with an initial misalignment embedded in an 8552 polymer matrix. The deformation of the model is governed by periodic boundary conditions (PBC). The relatively simple CMM model enables the evaluation of the role played by initial misalignment of the fiber, shear yielding of the matrix and fiber-matrix debonding. A novel microscale experimental technique devoted to the characterization of the longitudinal compressive strength of the fibers, X(sup f)(sub c) , was developed. By exercising the model and comparing it with several existing models in the literature, the nonlinear shear response of the composite lamina was shown to play a fundamental role not only in the prediction of the compressive strength X(sub c), but also during the post-peak regime in terms of residual stress σ(sub r) and fiber rotation φ. Finally, the influence of the fiber-matrix interface damage (not considered in most other fiber kinking models) on the fiber kinking phenomenon was assessed through a parametric study.

Polymer-matrix composites (PMC)↗
NASA NTRS2017

Flow-Tube Reactor Experiments on the High Temperature Oxidation of Carbon Weaves

Under entry conditions carbon weaves used in thermal protection systems (TPS) decompose via oxidation. Modeling this phenomenon is challenging due to the different regimes encountered along a flight trajectory. Approaches using equilibrium chemistry may lead to over-estimated mass loss and recession at certain conditions. Concurrently, there is a shortcoming of experimental data on carbon weaves to enable development of improved models. In this work, a flow-tube test facility was used to measure the oxidation of carbon weaves at temperatures up to 1500 K. The material tested was the 3D carbon weave used for the heat shield of the NASA Adaptive Deployable Entry and Placement Technology, ADEPT. Oxidation was characterized by quantifying decomposition gases (CO and CO2), by mass measurements, and by microscale surface analysis. The current set of measurements contributes to the development of finite rate chemistry models for carbon fabrics used in woven TPS materials.

Woven TPS↗
NASA NTRS2015

Integrated Modeling of Aerosol, Cloud, Precipitation and Land Processes at Satellite-Resolved Scales

With support from NASA's Modeling and Analysis Program, we have recently developed the NASA Unified-Weather Research and Forecasting model (NU-WRF). NU-WRF is an observation-driven integrated modeling system that represents aerosol, cloud, precipitation and land processes at satelliteresolved scales. "Satellite-resolved" scales (roughly 1e25 km), bridge the continuum between local (microscale), regional (mesoscale) and global (synoptic) processes. NU-WRF is a superset of the National Center for Atmospheric Research (NCAR) Advanced Research WRF (ARW) dynamical core model, achieved by fully integrating the GSFC Land Information System (LIS, already coupled to WRF), the WRF/ Chem enabled version of the Goddard Chemistry Aerosols Radiation Transport (GOCART) model, the Goddard Satellite Data Simulation Unit (G-SDSU), and custom boundary/initial condition preprocessors into a single software release, with source code available by agreement with NASA/GSFC. Full coupling between aerosol, cloud, precipitation and land processes is critical for predicting local and regional water and energy cycles.

Satellite-resolved↗
NASA NTRS2004

Cryogenic Calcite: A Morphologic and Isotopic Analog to the ALH84001 Carbonates

Martian meteorite ALH84001 carbonates preserve large and variable microscale isotopic compositions, which in some way reflect their formation environment. These measurements show large variations (>20%) in the carbon and oxygen isotopic compositions of the carbonates on a 10-20 micron scale that are correlated with chemical composition. However, the utilization of these data sets for interpreting the formation conditions of the carbonates is complex due to lack of suitable terrestrial analogs and the difficulty of modeling under non-equilibrium conditions. Thus, the mechanisms and processes are largely unknown that create and preserve large microscale isotopic variations in carbonate minerals. Experimental tests of the possible environments and mechanisms that lead to large microscale isotopic variations can help address these concerns. One possible mechanism for creating large carbon isotopic variations in carbonates involves the freezing of water. Carbonates precipitate during extensive CO2 degassing that occurs during the freezing process as the fluid s decreasing volume drives CO2 out. This rapid CO2 degassing results in a kinetic isotopic fractionation where the CO2 gas has a much lighter isotopic composition causing an enrichment of 13C in the remaining dissolved bicarbonate. This study seeks to determine the suitability of cryogenically formed carbonates as analogs to ALH84001 carbonates. Specifically, our objective is to determine how accurately models using equilibrium fractionation factors approximate the isotopic compositions of cryogenically precipitated carbonates. This includes determining the accuracy of applying equilibrium fractionation factors during a kinetic process, and determining how isotopic variations in the fluid are preserved in microscale variations in the precipitated carbonates.

Niles, P. B.↗
NASA NTRS2002

Deformation, Failure, and Fatigue Life of SiC/Ti-15-3 Laminates Accurately Predicted by MAC/GMC

NASA Glenn Research Center's Micromechanics Analysis Code with Generalized Method of Cells (MAC/GMC) (ref.1) has been extended to enable fully coupled macro-micro deformation, failure, and fatigue life predictions for advanced metal matrix, ceramic matrix, and polymer matrix composites. Because of the multiaxial nature of the code's underlying micromechanics model, GMC--which allows the incorporation of complex local inelastic constitutive models--MAC/GMC finds its most important application in metal matrix composites, like the SiC/Ti-15-3 composite examined here. Furthermore, since GMC predicts the microscale fields within each constituent of the composite material, submodels for local effects such as fiber breakage, interfacial debonding, and matrix fatigue damage can and have been built into MAC/GMC. The present application of MAC/GMC highlights the combination of these features, which has enabled the accurate modeling of the deformation, failure, and life of titanium matrix composites.

Bednarcyk, Brett A.↗
NASA NTRSDate not supplied

Artificial Generation of 2-D Fiber Reinforced Composite Microstructures with Statistically Equivalent Features

Fiber reinforced composites are used widely for their high strength and low weight advantages in various aerospace and automotive applications. While their use may be sought after, modeling of these material requires increasing fidelity at the lower scales to capture accurate material behavior under loading. The first steps in creating statistically equivalent models to real life cases is developing a method of rapid evaluation and artificial microstructure generation. The outlined work is capable of tracking microscale fiber positions and determining regions of localized volume fraction extrema (high and low end). Groupings of high and low volume fraction regions are called clusters and their geometry is used to characterize the microstructure. These cluster features can be evaluated for both artificial models and actual scans, allowing correlation to be established which can ultimately be used to regenerate statistically equivalent models. The results of this work show that if one feature is to be correlated, a model can be generated which matches almost exactly. But once more features are equally taken into account, the regeneration loses accuracy.

micromechanics↗
NASA NTRS2014

Attendee Response to Guiding Questions of the September 2014 DOE Workshop on North Slope of Alaska Priorities

We advocate for several modeling approaches suitable to prove understanding of and capability to model aerosol behavior in nucleating droplets and ice and to model regional cloud systems at the mesoscale. The full complement of data that is ideal for particular study types is increasingly well understood; the ASR ISDAC Science Plan provides an excellent example of listing the data required for each study type targeted by that campaign. Whereas discussion here focused on the microscale and mesoscale, other studies should obviously focus on the regional, seasonal and global scale. When considering ice microphysics alone, laser focus on unconstrained model components could be usefully intensified (e.g., collision-coalescence kernels). Similarly, a focus on boundary layer structure would be deeply valuable to a wide range of model types. The expense of observations motivates efficient planning for observational data sets to support specific study approaches; the concept of closure is often useful (e.g., the surface energy budget).

nucleation↗
NASA NTRSDate not supplied

Modeling the Effective Elasticity of Anisotropic Porous Materials

The development and optimization of composite materials designed for thermal protection of NASA’s spacecraft require understanding their physical response to high-enthalpy environments. To predict their macro-scale properties and behavior, high-fidelity 3D simulations are performed at the microscale on realistic representations of these composites. The digital microstructures are generated either synthetically or through X-ray micro-computed tomography reconstructions. One of the main challenges in the prediction of the structural response of heatshield materials is the computation of the effective elasticity of the fibrous composite, as well as the understanding of the deformation and stresses generated at the microscale. These are driven by the fiber layout within the microstructure and the distribution of the infused matrix. In this effort, the micro-mechanical linear elastic behavior of fibrous ablators is modeled using a numerical method based on the Multi-Point Stress Approximation (MPSA) finite volume scheme, a generalization of the more commonly used Multi-Point Flux Approximation (MPFA) that was presented at the 10th Ablation Workshop. To predict the behavior of fibrous and woven architectures, algorithms that compute the local fiber orientation are used. The implementation of the MPSA was verified using analytical solutions, engineering test cases, and compared against legacy Finite Element Analysis (FEA) software. The stress analysis models were then applied to real geometries used by NASA in thermal protection systems such as fibrous preforms and woven materials and the results were compared to experimental data.

Elasticity↗
NASA NTRS2000

Marangoni Effects on Near-Bubble Microscale Transport During Boiling of Binary Fluid Mixtures

In earlier investigations, Marangoni effects were observed to be the dominant mechanism of boiling transport in 2-propanol/water mixtures under reduced gravity conditions. In this investigation we have examined the mechanisms of binary mixture boiling by exploring the transport near a single bubble generated in a binary mixture between a heated surface and cold surface. The temperature field created in the liquid around the bubble produces vaporization over the portion of its interface near the heated surface and condensation over portions of its interface near the cold surface. Experiments were conducted using different mixtures of water and 2-propanol under 1g conditions and under reduced gravity conditions aboard the KC135 aircraft. Since 2-propanol is more volatile than water, there is a lower concentration of 2-propanol near the hot surface and a higher concentration of 2-propanol near the cold plate relative to the bulk quantity. This difference in interface concentration gives rise to strong Marangoni effects that move liquid toward the hot plate in the near bubble region for 2-propanol and water mixtures. In the experiments in this study, the pressure of the test system was maintained at about 5 kPa to achieve the full spectrum of boiling behavior (nucleate boiling, critical heat flux and film boiling) at low temperature and heat flux levels. Heat transfer data and visual documentation of the bubble shape were extracted from the experimental results. In the 1-g experiments at moderate to high heat flux levels, the bubble was observed to grow into a mushroom shape with a larger top portion near the cold plate due to the buoyancy effect. The shape of the bubble was somewhat affected by the cold plate subcooling and the superheat of the heated surface. At low superheat levels for the heated surface, several active nucleation sites were observed, and the vapor stems from them merged to form a larger bubble. The generation rate of vapor is moderate in this regime and the bubble shape is cylindrical in appearance. In some instances, the bubble interface appeared to oscillate. At higher applied heat flux levels, the top of the bubble became larger, apparently to provide more condensing interface area adjacent to the cold plate. Increasing the applied heat flux ultimately led to dry-out of the heated surface, with conditions just prior to dryout corresponding to the maximum heat flux (CHF). A more stable bubble was observed when the system attained the minimum heat flux (for film boiling). In this regime, most of the surface under the bottom of the bubble was dry with nucleate boiling sometimes occuring around the contact perimeter of the bubble at heated surface. Different variations (e.g. gap between two plates, molar concentration of the liquid mixture) of the experiments were examined to determine parametric effects on the boiling process and to determine the best conditions for the KC135 reduced gravity tests. Variation of the gap was found to have a minor impact on the CHF. However, reducing the gap between the hot and cold surface was observed to significantly reduce the minimum heat flux for fixed molar concentration of 2-propanol. In the reduced gravity experiments aboard the KC135 aircraft, the bubble formed in the 6.4 mm gap was generally cylindrical or barrel shaped and it increased its extent laterally as the surface superheat increased. In reduced gravity experiments, dryout of the heated surface under the bubble was observed to occur at a lower superheated temperature than for 1g conditions. Observed features of the boiling process and heat transfer data under reduced gravity will be discussed in detail. The results of the reduced gravity experiments will also be compared to those obtained in comparable 1g experiments. In tandem with the experiments we are also developing a computational model of the transport in the liquid surrounding the bubble during the boiling process. The computational model uses a level set method to model motion of the interface. It will incorporate a macroscale treatment of the transport in the liquid gap between the surfaces and a microscale treatment of transport in the regions between the bubble interface and the solid surfaces. The features of the model will be described in detail. Future research directions suggested by the results to date will also be discussed.

V. Carey↗
NASA NTRS1993

EFEDA - European field experiment in a desertification-threatened area

During June 1991 more than 30 scientific teams worked in Castilla-La Mancha, Spain, studying the energy and water transfer processes between soil, vegetation, and the atmosphere in semiarid conditions within the coordinated European research project EFEDA (European Field Experiment in Desertification-threatened Areas). Measurements were made from the microscale (e.g., measurements on single plants) up to a scale compatible with the grid size of global models. For this purpose three sites were selected 70 km apart and heavily instrumented at a scale in the order of 30 sq km. Aircraft missions, satellite data, and movable equipment were deployed to provide a bridge to the larger scale. This paper gives a description of the experimental design along with some of the preliminary results of this successful experiment.

Bolle, H.-J.↗
NASA NTRS1985

Simulation and modeling of homogeneous, compressed turbulence

Low Reynolds number homogeneous turbulence undergoing low Mach number isotropic and one-dimensional compression was simulated by numerically solving the Navier-Stokes equations. The numerical simulations were performed on a CYBER 205 computer using a 64 x 64 x 64 mesh. A spectral method was used for spatial differencing and the second-order Runge-Kutta method for time advancement. A variety of statistical information was extracted from the computed flow fields. These include three-dimensional energy and dissipation spectra, two-point velocity correlations, one-dimensional energy spectra, turbulent kinetic energy and its dissipation rate, integral length scales, Taylor microscales, and Kolmogorov length scale. Results from the simulated flow fields were used to test one-point closure, two-equation models. A new one-point-closure, three-equation turbulence model which accounts for the effect of compression is proposed. The new model accurately calculates four types of flows (isotropic decay, isotropic compression, one-dimensional compression, and axisymmetric expansion flows) for a wide range of strain rates.

Wu, C. T.↗
NASA NTRS1999

Characteristics of Non-Premixed Turbulent Flames in Microgravity

The momentum of the fuel (and/or air) jet is important in classifying gas-jet diffusion flame behavior. Normal-gravity data on gas-jet flames show that the flame height (non-dimensionalized with respect to an effective diameter) can be correlated to a density weighted Froude number in the buoyancy-dominated limit. In the momentum-dominated limit this non-dimensional flame height asymptotes to a constant value. The momentum-dominated limit under normal gravity conditions is usually obtained for very high injection velocities which in turn results in high values of the injection Reynolds number. This results in a complicated flame structure because of the large number of turbulence scales involved. In order to gain better insight into the structure of these flames it would be useful to reduce the injection Reynolds number while still maintaining turbulent conditions. This can be done in microgravity where momentum-dominated turbulent flames are obtained at much smaller velocities than in normal gravity. In this paper, experimental results on the effects of nozzle diameter and fuel dilution on flame height are discussed. The experimental values are compared with predictions from a numerical procedure utilizing the standard k-epsilon turbulence model. Flame height scaling with nozzle size and dilution is established. Differences between model predictions and measurements are presented. In order to explain these differences, evolutions of turbulent spectra and Taylor microscale along the flame axis are considered.

Hegde, U.↗
NASA NTRS2003

Large- and Small-Scale Ring Current Electrodynamic Coupling

In this talk we will address the two primary issues of ring current (RC) electrodynamic coupling: 1. RC self-consistent magnetosphere-ionosphere coupling that includes calculation of the magnetospheric electric field (large scale electrodynamic coupling); and 2. RC self-consistent coupling with electromagnetic ion cyclotron (EMIC) waves (small scale electrodynamic coupling). Our study will be based on two RC models that we have recently developed in our group. The first model by Khazanov et al. [2002] couples the system of two kinetic equations: one equation which describes the RC ion dynamics and another equation which describes the energy density evolution of EMIC waves. The second model by Khazanov et al. [2003] deals with large scale electrodynamic coupling processes and provides a self-consistent simulation of RC ions and the magnetospheric electric field. There is presently no model that addresses both of these issues simultaneously in a self-consistent calculation. However, the need exists for such a model, because these two processes directly influence each other, with the mesoscale coupling changing the drift paths of the thermal and energetic particle populations in the inner magnetosphere, thereby changing the wave interactions, and the microscale coupling altering the pitch angle distributions and ionospheric conductivities (through increased precipitation), thus changing the field-aligned currents and electric potential structure. The initial thrust of the work will be the development of a combined kinetic model of micro- and meso-scale RC electrodynamic coupling processes and to examine their interactions with each other on a global scale.

Khazanov, G. V.↗