Multiscale Simulation of Deployable Composite Structures
Explore the source record for details and available documents.
SEARCH · Engineering Papers
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.
Explore the source record for details and available documents.
The Atlanta Urban Heat Island and Air Quality Project had its genesis in Project ATLANTA (ATlanta Land use Analysis: Temperature and Air quality) that began in 1996. Project ATLANTA examined how high-spatial resolution thermal remote sensing data could be used to derive better measurements of the Urban Heat Island effect over Atlanta. We have explored how these thermal remote sensing, as well as other imaged datasets, can be used to better characterize the urban landscape for improved air quality modeling over the Atlanta area. For the air quality modeling project, the National Land Cover Dataset and the local scale Landpro99 dataset at 30m spatial resolutions have been used to derive land use/land cover characteristics for input into the MM5 mesoscale meteorological model that is one of the foundations for the Community Multiscale Air Quality (CMAQ) model to assess how these data can improve output from CMAQ. Additionally, land use changes to 2030 have been predicted using a Spatial Growth Model (SGM). SGM simulates growth around a region using population, employment and travel demand forecasts. Air quality modeling simulations were conducted using both current and future land cover. Meteorological modeling simulations indicate a 0.5 C increase in daily maximum air temperatures by 2030. Air quality modeling simulations show substantial differences in relative contributions of individual atmospheric pollutant constituents as a result of land cover change. Enhanced boundary layer mixing over the city tends to offset the increase in ozone concentration expected due to higher surface temperatures as a result of urbanization.
A probabilistic method is used to evaluate the structural reliability and risk of select metallic and composite structures. The method is a multiscale, multifunctional and it is based on the most elemental level. A multi-factor interaction model is used to describe the material properties which are subsequently evaluated probabilistically. The metallic structure is a two rotor aircraft engine, while the composite structures consist of laminated plies (multiscale) and the properties of each ply are the multifunctional representation. The structural component is modeled by finite element. The solution method for structural responses is obtained by an updated simulation scheme. The results show that the risk for the two rotor engine is about 0.0001 and the composite built-up structure is also 0.0001.
A probabilistic method is used to evaluate the structural reliability and risk of select metallic and composite structures. The method is a multiscale, multifunctional and it is based on the most elemental level. A multifactor interaction model is used to describe the material properties which are subsequently evaluated probabilistically. The metallic structure is a two rotor aircraft engine, while the composite structures consist of laminated plies (multiscale) and the properties of each ply are the multifunctional representation. The structural component is modeled by finite element. The solution method for structural responses is obtained by an updated simulation scheme. The results show that the risk for the two rotor engine is about 0.0001 and the composite built-up structure is also 0.0001.
Soil dust aerosols are a key component of the climate system, as they interact with short- and long-wave radiation, alter cloud formation processes, affect atmospheric chemistry and play a role in biogeochemical cycles by providing nutrient inputs such as iron and phosphorus. The influence of dust on these processes depends on its physicochemical properties, which, far from being homogeneous, are shaped by its regionally varying mineral composition. The relative amount of minerals in dust depends on the source region and shows a large geographical variability. However, many state-of-the-art Earth system models (ESMs), upon which climate analyses and projections rely, still consider dust mineralogy to be invariant. The explicit representation of minerals in ESMs is more hindered by our limited knowledge of the global soil composition along with the resulting size-resolved airborne mineralogy than by computational constraints. In this work we introduce an explicit mineralogy representation within the state-of-the-art Multiscale Online Nonhydrostatic AtmospheRe CHemistry (MONARCH) model. We review and compare two existing soil mineralogy datasets, which remain a source of uncertainty for dust mineralogy modeling and provide an evaluation of multiannual simulations against available mineralogy observations. Soil mineralogy datasets are based on measurements performed after wet sieving, which breaks the aggregates found in the parent soil. Our model predicts the emitted particle size distribution (PSD) in terms of its constituent minerals based on brittle fragmentation theory (BFT), which reconstructs the emitted mineral aggregates destroyed by wet sieving. Our simulations broadly reproduce the most abundant mineral fractions independently of the soil composition data used. Feldspars and calcite are highly sensitive to the soil mineralogy map, mainly due to the different assumptions made in each soil dataset to extrapolate a handful of soil measurements to arid and semi-arid regions worldwide. For the least abundant or more difficult-to-determine minerals, such as iron oxides, uncertainties in soil mineralogy yield differences in annual mean aerosol mass fractions of up to ∼ 100 %. Although BFT restores coarse aggregates including phyllosilicates that usually break during soil analysis, we still identify an overestimation of coarse quartz mass fractions (above 2 µm in diameter). In a dedicated experiment, we estimate the fraction of dust with undetermined composition as given by a soil map, which makes up ∼ 10 % of the emitted dust mass at the global scale and can be regionally larger. Changes in the underlying soil mineralogy impact our estimates of climate-relevant variables, particularly affecting the regional variability of the single-scattering albedo at solar wavelengths or the total iron deposited over oceans. All in all, this assessment represents a baseline for future model experiments including new mineralogical maps constrained by high-quality spaceborne hyperspectral measurements, such as those arising from the NASA Earth Surface Mineral Dust Source Investigation (EMIT) mission.
This paper presents a comprehensive description and benchmark evaluation of the tropospheric gas-phase chemistry component of the Multiscale Online Nonhydrostatic AtmospheRe CHemistry model (NMMBMONARCH), formerly known as NMMB/BSC-CTM, that can be run on both regional and global domains. Here, we provide an extensive evaluation of a global annual cycle simulation using a variety of background surface stations (EMEP, WDCGG and CASTNET), ozonesondes (WOUDC, CMD and SHADOZ), aircraft data (MOZAIC and several campaigns), and satellite observations (SCIAMACHY and MOPITT).We also include an extensive discussion of our results in comparison to other state-of-the-art models. We note that in this study, we omitted aerosol processes and some natural emissions (lightning and volcano emissions). The model shows a realistic oxidative capacity across the globe. The seasonal cycle for CO is fairly well represented at different locations (correlations around 0.3-0.7 in surface concentrations), although concentrations are underestimated in spring and winter in the Northern Hemisphere, and are overestimated throughout the year at 800 and 500 hPa in the Southern Hemisphere. Nitrogen species are well represented in almost all locations, particularly NO2 in Europe (root mean square error - RMSE - below 5 ppb). The modeled vertical distributions of NOx and HNO3 are in excellent agreement with the observed values and the spatial and seasonal trends of tropospheric NO2 columns correspond well to observations from SCIAMACHY, capturing the highly polluted areas and the biomass burning cycle throughout the year. Over Asia, the model underestimates NOx from March to August, probably due to an underestimation of NOx emissions in the region. Overall, the comparison of the modeled CO and NO2 with MOPITT and SCIAMACHY observations emphasizes the need for more accurate emission rates from anthropogenic and biomass burning sources (i.e., specification of temporal variability).
High resolution observations of dense plasma impulsive structures moving through an ambient background magnetospheric flows were captured by the Magnetospheric Multiscale mission. The observations show particle heating and acceleration, shock-like wave formation, and whistler wave excitation inside the interface between the dense impulsive plasma structures and the ambient plasma. A multiscale hybrid kinetic simulation provides an explanation of the observed wave-particle interactions with the assumption that the dense plasma structures may be represented by plasma clouds which are formed at the magnetopause layer due to reconnection processes. Plain Language Summary: Dense, impulsive plasma structures moving through a background plasma were captured by the NASA Magnetospheric Multiscale mission. The observations show that the dense structures can generate strong perturbations in the electromagnetic field and shock-like waves. Interactions between these electromagnetic waves and the particles results in particle acceleration. 3-D hybrid kinetic modeling (particle description for ions and fluid description for electrons) was used to investigate the plasma physics of the observed structures. It was assumed that the plasma clouds were formed by magnetic field reconnection inside the magnetopause, which is the interface between the solar wind particles and the cold low-density magnetospheric plasma. The work helps us understand the plasma environment at the interface between the Earth and solar wind, near planetary moons, within astrophysical explosions, and possibly at the interface between the solar wind and local interstellar medium.
High strain thin-ply (HS-TPC) technology is being increasingly adopted for high-performanceaerospace applications. Albeit many of these structures such as deployable booms can bemodeled as one-dimensional beam problems, there is a lack of themoviscoelastic beam mod-els to efficiently and accurately simulate HS-TPC. This work will use mechanics of struc-ture genome (MSG) to construct linear thermoviscoelastic beam models that can homogenizethree-dimensional heterogeneous materials made of constituents with time- and temperature-dependent behavior. The formulation derives the transient strain energy based on integralformulation for thermorheologically simple materials subject to finite temperature changeswith the restriction that the strain is small. A lenticular boom is used as a numerical exampleto verify the MSG-based linear thermoviscoelastic beam model against MSG-based shell/platemodel, which has already been validated against experimental data provided by NASA, anddirect numerical simulations performed in a finite element commercial package.
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.
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.
To understand the role of magnetic reconnection in global evolution of magnetosphere and to place spacecraft observations into global context it is essential to perform global simulations with physically motivated model of dissipation that is capable to reproduce reconnection rates predicted by kinetic models. In our efforts to bridge the gap between small scale kinetic modeling and global simulations we introduced an approach that allows to quantify the interaction between large-scale global magnetospheric dynamics and microphysical processes in diffusion regions near reconnection sites. We utilized the high resolution global MHD code BATSRUS and incorporate primary mechanism controlling the dissipation in the vicinity of reconnection sites in terms of kinetic corrections to induction and energy equations. One of the key elements of the multiscale modeling of magnetic reconnection is identification of reconnection sites and boundaries of surrounding diffusion regions where non-MHD corrections are required. Reconnection site search in the equatorial plane implemented in our previous studies is extended to cusp and magnetopause reconnection, as well as for magnetotail reconnection in realistic asymmetric configurations. The role of feedback between the non-ideal effects in diffusion regions and global magnetosphere structure and dynamics will be discussed.
This chapter describes some of the building blocks to ensure a higher level of confidence in the predictability and reliability (PAR) of numerical simulation of multiscale complex nonlinear problems. The focus is on relating PAR of numerical simulations with complex nonlinear phenomena of numerics. To isolate sources of numerical uncertainties, the possible discrepancy between the chosen partial differential equation (PDE) model and the real physics and/or experimental data is set aside. The discussion is restricted to how well numerical schemes can mimic the solution behavior of the underlying PDE model for finite time steps and grid spacings. The situation is complicated by the fact that the available theory for the understanding of nonlinear behavior of numerics is not at a stage to fully analyze the nonlinear Euler and Navier-Stokes equations. The discussion is based on the knowledge gained for nonlinear model problems with known analytical solutions to identify and explain the possible sources and remedies of numerical uncertainties in practical computations. Examples relevant to turbulent flow computations are included.
This talk describes some of the building blocks to ensure a higher level of confidence in the predictability and reliability (PAR) of numerical simulation of multiscale complex nonlinear problems. The focus is on relating PAR of numerical simulations with complex nonlinear phenomena of numerics. To isolate sources of numerical uncertainties, the possible discrepancy between the chosen partial differential equation (PDE) model and the real physics and/or experimental data is set aside. The discussion is restricted to how well numerical schemes can mimic the solution behavior of the underlying PDE model for finite time steps and grid spacings. The situation is complicated by the fact that the available theory for the understanding of nonlinear behavior of numerics is not at a stage to fully analyze the nonlinear Euler and Navier-Stokes equations. The discussion is based on the knowledge gained for nonlinear model problems with known analytical solutions to identify and explain the possible sources and remedies of numerical uncertainties in practical computations. Examples relevant to turbulent flow computations are included.
This work describes some of the procedure to ensure a higher level of confidence in the predictability and reliability (PAR) of numerical simulation of multiscale complex nonlinear problems. The focus is on relating PAR of numerical simulations with complex nonlinear phenomena of numerics. To isolate sources of numerical uncertainties, the possible discrepancy between the chosen partial differential equation (PDE) model and the real physics and/or experimental data is set aside. The discussion is restricted to how well numerical schemes can mimic the solution behavior of the underlying PDE model for finite time steps and grid spacings. The situation is complicated by the fact that the available theory for the understanding of nonlinear behavior of numerics is not at a stage to fully analyze the nonlinear Euler and Navier-Stokes equations. The discussion is based on the knowledge gained for nonlinear model problems with known analytical solutions to identify and explain the possible sources and remedies of numerical uncertainties in practical computations. Examples relevant to turbulent flow computations are included.
This chapter describes some of the building blocks to ensure a higher level of confidence in the predictability and reliability (PAR) of numerical simulation of multiscale complex nonlinear problems. The focus is on relating PAR of numerical simulations with complex nonlinear phenomena of numerics. To isolate sources of numerical uncertainties, the possible discrepancy between the chosen partial differential equation (PDE) model and the real physics and/or experimental data is set aside. The discussion is restricted to how well numerical schemes can mimic the solution behavior of the underlying PDE model for finite time steps and grid spacings. The situation is complicated by the fact that the available theory for the understanding of nonlinear behavior of numerics is not at a stage to fully analyze the nonlinear Euler and Navier-Stokes equations. The discussion is based on the knowledge gained for nonlinear model problems with known analytical solutions to identify and explain the possible sources and remedies of numerical uncertainties in practical computations.
Multiscale progressive failure simulations have been performed for a 3D woven composite considering six length scales, spanning from the woven composite repeating unit cell to the matrix material containing voids. The Multiscale Recursive Micromechanics approach, which enables micromechanics models to call other micromechanics models (or themselves recursively) to consider finer and finer length scales, has been employed. The multiscale model uses both the generalized method of cells and Mori-Tanaka micromechanics theories and considers failure and damage at the lowest length scale. A simple subvolume elimination damage model, as well as the crack band model, are employed for the matrix and their predicted composite responses compared with and without a preexisting binder tow disbond. Results are also compared with experimental data for an AS4 carbon fiber/RTM6 epoxy matrix 3D orthogonal woven PMC, with good agreement in terms of global stiffness and global failure stress.
Modeling atmospheric chemistry at fine resolution globally is computationally expensive; the capability to focus on specific geographic regions using a multi scale grid is desirable. Here, we develop, validate, and demonstrate stretched grids in the GEOS-Chem atmospheric chemistry model in its high-performance implementation (GCHP). These multiscale grids are specified at runtime by four parameters that offer users nimble control of the region that is refined and the resolution of the refinement. We validate the stretched-grid simulation versus global cubed-sphere simulations. We demonstrate the operation and flexibility of stretched-grid simulations with two case studies that compare simulated tropospheric NO2column densities from stretched-grid and cubed-sphere simulations to retrieved column densities from the TROPOspheric Monitoring Instrument (TROPOMI). The first case study uses a stretched grid with a broad refinement covering the contiguous US to produces imulated columns that perform similarly to a C180 (∼50 km) cubed-sphere simulation at less than one-ninth the computational expense. The second case study experiments with a large stretch-factor for a global stretched-grid simulation with a highly localized refinement with∼10 km resolution for California. We find that the refinement improves spatial agreement with TROPOMI columns compared to a C90 cubed-sphere simulation of comparable computational demands. Overall we find that stretched grids in GEOS-Chem are a practical tool for fine resolution regional- or continental-scale simulations of atmospheric chemistry. Stretched grids are available in GEOS-Chem version 13.0.0
Fiber reinforced composite materials, owing to their tailorable thermomechanical and functional properties, allow one to produce a structure that is stronger, stiffer, and lighter than its metal counterpart while performing the same function, yielding a more efficient structure. This not only allows for the improvement of current technologies like aircraft structures, but also enables new technologies like gaseous hydrogen storage for mobility applications, which are otherwise impractical when manufactured using traditional metals due to weight and space restrictions or material embrittlement. However, the use of composites imposes greater design and manufacturing challenges on an engineer, since they are heterogenous, having a distinct structure across multiple length scale, behave generally anisotropically at the structural level and require complex manufacturing and processing methods. Capturing this complex behavior requires detailed numerical simulations, including the modeling of microstructural features like undulations, voids, and fiber alignment. In this paper, multiple repeating unit cells (RUCs), representing filament wound composites, are developed (via a script provided in the Appendix) and analyzed. The refinement of these RUCs is varied, and the analyses are performed using both the Abaqus finite element software and the NASA Multiscale Analysis Tool (NASMAT). A study is undertaken to compare the predicted effective elastic properties of the wound RUC to a laminate representation of the wound RUC, which neglects the undulations. Additionally, two different sets of periodic boundary conditions (PBCs) have been examined. One approximates the real boundary conditions using a standard approach and the other represents the PBCs exactly through the use of an offset. Lastly, a comparison of the local elastic stress fields is made among the models and approaches. Since wound structures are often approximated as laminated structures, it is important to understand the degree to which this assumption is valid, namely by first comparing the elastic constants and local elastic fields. This will provide, on the one hand, information concerning the bulk mechanical behavior and, on the other hand, insights concerning local load distributions and likely damage initiation sites.