Search NASA⌕ Search

SEARCH · Search NASA

Results for “Structural simulation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 847 records · Page 47
All sourcesNASA NTRSDOE OSTI
NASA NTRS2016

Stochastic-Strength-Based Damage Simulation Tool for Ceramic Matrix and Polymer Matrix Composite Structures

Stochastic-based, discrete-event progressive damage simulations of ceramic-matrix composite and polymer matrix composite material structures have been enabled through the development of a unique multiscale modeling tool. This effort involves coupling three independently developed software programs: (1) the Micromechanics Analysis Code with Generalized Method of Cells (MAC/GMC), (2) the Ceramics Analysis and Reliability Evaluation of Structures Life Prediction Program (CARES/ Life), and (3) the Abaqus finite element analysis (FEA) program. MAC/GMC contributes multiscale modeling capabilities and micromechanics relations to determine stresses and deformations at the microscale of the composite material repeating unit cell (RUC). CARES/Life contributes statistical multiaxial failure criteria that can be applied to the individual brittle-material constituents of the RUC. Abaqus is used at the global scale to model the overall composite structure. An Abaqus user-defined material (UMAT) interface, referred to here as "FEAMAC/CARES," was developed that enables MAC/GMC and CARES/Life to operate seamlessly with the Abaqus FEA code. For each FEAMAC/CARES simulation trial, the stochastic nature of brittle material strength results in random, discrete damage events, which incrementally progress and lead to ultimate structural failure. This report describes the FEAMAC/CARES methodology and discusses examples that illustrate the performance of the tool. A comprehensive example problem, simulating the progressive damage of laminated ceramic matrix composites under various off-axis loading conditions and including a double notched tensile specimen geometry, is described in a separate report.

composite structures↗
DOE OSTI2020

Patchy particle model of hydrated amorphous silica

Glasses corroded in static aqueous conditions typically reach a slow, residual rate of dissolution as a result of the complex interplay between multiple processes. While intrinsically suited to address this problem, Monte Carlo (MC) models developed to date have relied on lattices to describe the amorphous structure of glass and are thus unable to reproduce the residual rate of glass corrosion. A recently developed MC approach based on amorphous structures derived from molecular dynamics simulations offers a solution [Kerisit and Du J. Non-Cryst. Solids 522 (2019) 119601], but it requires a simple model that can rapidly relax the glass-water interface upon each dissolution/condensation event to retain the ability of MC simulations to reach large spatial and temporal scales. To this end, a patchy particle model of hydrated amorphous silica was developed and evaluated against atomistic simulations in this work. Models of water and amorphous silica were first developed separately, and cross terms were then defined to enable simulations of hydrated amorphous silica structures with varying water content. While the nature of its interaction potential prevents the patchy particle model from reproducing simultaneously the structure and dynamics of the systems of interest with a single set of parameters, it described well the connectivity of hydrated amorphous silica structures and the local coordination geometry of individual species. Therefore, this model opens the door to MC simulations of the residual rate of glass corrosion and offers an alternative to atomistic models to investigate the evolution of silica gels over large spatial and temporal scales.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2025

X-ray and neutron diffraction patterns of the AlCrTiV high-entropy alloy and quaternary Heusler structures

The quaternary alloy AlCrTiV has been proposed as both a lightweight high-entropy alloy and also a functional spin filter material based on the Heusler structure. Experimental investigations to-date, based on X-ray diffraction, offer conflicting interpretations of the structure. Here we simulate diffraction patterns of the various proposed structures to show that neutron diffraction, in particular, can reveal the nature of long-range chemical order and discriminate among distributions of the refractory transition metals. Magnetic contributions to the neutron diffraction are also discussed.

Kalliney, Nedjma [Carnegie Mellon Univ., Pittsburg↗
DOE OSTI2020

TDP-43 α-helical structure tunes liquid–liquid phase separation and function

Liquid–liquid phase separation (LLPS) is involved in the formation of membraneless organelles (MLOs) associated with RNA processing. The RNA-binding protein TDP-43 is present in several MLOs, undergoes LLPS, and has been linked to the pathogenesis of amyotrophic lateral sclerosis (ALS). While some ALS-associated mutations in TDP-43 disrupt self-interaction and function, here we show that designed single mutations can enhance TDP-43 assembly and function via modulating helical structure. Using molecular simulation and NMR spectroscopy, we observe large structural changes upon dimerization of TDP-43. Two conserved glycine residues (G335 and G338) are potent inhibitors of helical extension and helix–helix interaction, which are removed in part by variants at these positions, including the ALS-associated G335D. Substitution to helix-enhancing alanine at either of these positions dramatically enhances phase separation in vitro and decreases fluidity of phase-separated TDP-43 reporter compartments in cells. Furthermore, G335A increases TDP-43 splicing function in a minigene assay. Therefore, the TDP-43 helical region serves as a short but uniquely tunable module where application of biophysical principles can precisely control assembly and function in cellular and synthetic biology applications of LLPS.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2024

Chasing the multi-modal plasma response in MAST-U

Abstract Achieving edge localized modes (ELMs) suppression in spherical tokamaks by applying resonant magnetic perturbations (RMPs) has proven challenging. The poloidal spectrum of the applied RMP is a key parameter that has an impact on the capability to mitigate and eventually suppress ELMs. In this work the resistive magnetohydrodynamic code MARS-F (Liu et al 2000 Phys. Plasmas 7 3681) is used to evaluate the possibility of directly measuring the plasma response in MAST-U, and particularly its variation as function of the applied poloidal spectrum, in order to guide the experimental validation of the predicted best RMP configuration for ELM suppression. Toroidal mode number n = 2 RMP is considered to minimize the presence of sidebands, and to avoid the deleterious core coupling of n = 1. Singular Value Decomposition is used to highlight linearly independent structures in the simulated magnetic 3D fields and how those structures can be measured at the wall where the magnetic sensors are located. Alternative ways to measure the multimodal plasma response and how they can be used to infer the best RMP configuration to achieve ELM suppression are also presented, including the plasma displacement and the 3D footprints at the divertor plates.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Compare linear-system solver and preconditioner stacks with emphasis on GPU performance and propose phase-2 NGP solver development pathway

The goal of the ExaWind project is to enable predictive simulations of wind farms comprised of many megawatt-scale turbines situated in complex terrain. Predictive simulations will require computational fluid dynamics (CFD) simulations for which the mesh resolves the geometry of the turbines and captures the rotation and large deflections of blades. Whereas such simulations for a single turbine are arguably petascale class, multi-turbine wind farm simulations will require exascale-class resources. The primary physics codes in the ExaWind project are Nalu-Wind, which is an unstructured-grid solver for the acoustically incompressible Navier-Stokes equations, and OpenFAST, which is a whole-turbine simulation code. The Nalu-Wind model consists of the mass-continuity Poisson-type equation for pressure and a momentum equation for the velocity. For such modeling approaches, simulation times are dominated by linear-system setup and solution for the continuity and momentum systems. For the ExaWind challenge problem, the moving meshes greatly affect overall solver costs as reinitialization of matrices and recomputation of preconditioners is required at every time step. In this report we evaluated GPU-performance baselines for the linear solvers in the Trilinos and hypre solver stacks using two representative Nalu-Wind simulations: an atmospheric boundary layer precursor simulation on a structured mesh, and a fixed-wing simulation using unstructured overset meshes. Both strong-scaling and weak-scaling experiments were conducted on the OLCF supercomputer Summit and similar proxy clusters. We focused on the performance of multi-threaded Gauss-Seidel and two-stage Gauss-Seidel that are extensions of classical Gauss-Seidel; of one-reduce GMRES, a communication-reducing variant of the Krylov GMRES; and algebraic multigrid methods that incorporate the afore-mentioned methods. The team has established that AMG methods are capable of solving linear systems arising from the fixed-wing overset meshes on CPU, a critical intermediate result for ExaWind FY20 Q3 and Q4 milestones. For the fixed-wing strong-scaling study (model with 3M grid-points), the team identified that Nalu-Wind simulations with the new Trilinos and hypre solvers scale to modest GPU counts, maintaining above 70% efficiency up to 6 GPUs. However, there still remain significant bottlenecks to performance: matrix assembly (hypre), AMG setup (hypre and Trilinos) In the weak-scaling experiments (going from 0.4M to 211M gridpoints), it's shown that the solver apply phases are faster on GPUs, but that Nalu-Wind simulation times grow, primarily due to the multigrid-setup process. Finally, based on the report outcomes, we propose a linear solver path-forward for the remainder of the ExaWind project. Near term, the NREL team will continue their work on GPU-based linear-system assembly. They will also investigate how the use of alternatives to the NVIDIA UVM (unified virtual memory) paradigm affects performance. Longer term, the NREL team will evaluate algorithmic performance on other types of accelerators and merge their improvements back to the main hypre repository branch. Near term, the Trilinos team will address performance bottlenecks identified in this milestone, such as implementing a GPU-based segregated momentum solve and reusing matrix graphs across linear-system assembly phases. Longer term, the Trilinos team will do detailed analysis and optimization of multigrid setup.

17 WIND ENERGY↗
DOE OSTI2024

Novel sterol binding domains in bacteria

Sterol lipids are widely present in eukaryotes and play essential roles in signaling and modulating membrane fluidity. Although rare, some bacteria also produce sterols, but their function in bacteria is not known. Moreover, many more species, including pathogens and commensal microbes, acquire or modify sterols from eukaryotic hosts through poorly understood molecular mechanisms. The aerobic methanotroph Methylococcus capsulatus was the first bacterium shown to synthesize sterols, producing a mixture of C-4 methylated sterols that are distinct from those observed in eukaryotes. C-4 methylated sterols are synthesized in the cytosol and localized to the outer membrane, suggesting that a bacterial sterol transport machinery exists. Until now, the identity of such machinery remained a mystery. In this study, we identified three novel proteins that may be the first examples of transporters for bacterial sterol lipids. The proteins, which all belong to well-studied families of bacterial metabolite transporters, are predicted to reside in the inner membrane, periplasm, and outer membrane of M. capsulatus, and may work as a conduit to move modified sterols to the outer membrane. Quantitative analysis of ligand binding revealed their remarkable specificity for 4-methylsterols, and crystallographic structures coupled with docking and molecular dynamics simulations revealed the structural bases for substrate binding by two of the putative transporters. Their striking structural divergence from eukaryotic sterol transporters signals that they form a distinct sterol transport system within the bacterial domain. Finally, bioinformatics revealed the widespread presence of similar transporters in bacterial genomes, including in some pathogens that use host sterol lipids to construct their cell envelopes. The unique folds of these bacterial sterol binding proteins should now guide the discovery of other proteins that handle this essential metabolite.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS1993

Nonlinear dynamic simulation of single- and multi-spool core engines

In this paper a new computational method for accurate simulation of the nonlinear dynamic behavior of single- and multi-spool core engines, turbofan engines, and power generation gas turbine engines is presented. In order to perform the simulation, a modularly structured computer code has been developed which includes individual mathematical modules representing various engine components. The generic structure of the code enables the dynamic simulation of arbitrary engine configurations ranging from single-spool thrust generation to multi-spool thrust/power generation engines under adverse dynamic operating conditions. For precise simulation of turbine and compressor components, row-by-row calculation procedures were implemented that account for the specific turbine and compressor cascade and blade geometry and characteristics. The dynamic behavior of the subject engine is calculated by solving a number of systems of partial differential equations, which describe the unsteady behavior of the individual components. In order to ensure the capability, accuracy, robustness, and reliability of the code, comprehensive critical performance assessment and validation tests were performed. As representatives, three different transient cases with single- and multi-spool thrust and power generation engines were simulated. The transient cases range from operating with a prescribed fuel schedule, to extreme load changes, to generator and turbine shut down.

Schobeiri, T.↗
NASA NTRS1996

Damage progression in Composite Structures

A computational simulation tool is used to evaluate the various stages of damage progression in composite materials during Iosipescu sheat testing. Unidirectional composite specimens with either the major or minor material axis in the load direction are considered. Damage progression characteristics are described for each specimen using two types of boundary conditions. A procedure is outlined regarding the use of computational simulation in composites testing. Iosipescu shear testing using the V-notched beam specimen is a convenient method to measure both shear strength and shear stiffness simultaneously. The evaluation of composite test response can be made more productive and informative via computational simulation of progressive damage and fracture. Computational simulation performs a complete evaluation of laminated composite fracture via assessment of ply and subply level damage/fracture processes.

Minnetyan, Levon↗
DOE OSTI2020

Modeling magnetic confinement of laser-generated plasma in cylindrical geometry leading to disk-shaped structures

Radiation-magnetohydrodynamic simulations were able to reproduce features of the plasma structures observed in recent experiments1, where a laser was used to ablate plasma in a 3 MG magnetic field. The laser ablates the plasma, and localized structures are formed due to the inhibition of heat flow by the magnetic field. The large magnetic pressures cause the plasma to pinch. In an azimuthal fi eld, a disc-shaped plasma is generated. According to simulations, the disc has electron densities that are underdense to the laser, ranging from 10 18 to 10 19 cm -3 , and electron temperatures in the 300 to 1000 eV range during its evolution, similar to experimental data.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Ion acceleration from microstructured targets irradiated by high-intensity picosecond laser pulses

Structures on the front surface of thin foil targets for laser-driven ion acceleration have been proposed to increase the ion source maximum energy and conversion efficiency. While structures have been shown to significantly boost the proton acceleration from pulses of moderate-energy fluence, their performance on tightly focused and high-energy lasers remains unclear. Here, we report the results of laser-driven three-dimensional (3D)-printed microtube targets, focusing on their efficacy for ion acceleration. Using the high-contrast (~10 12) PHELIX laser (150 J, 10 21 W / cm 2 ), we studied the acceleration of ions from 1-μm-thick foils covered with micropillars or microtubes, which we compared with flat foils. The front-surface structures significantly increased the conversion efficiency from laser to light ions, with up to a factor of 5 higher proton number with respect to a flat target, albeit without an increase of the cutoff energy. An optimum diameter was found for the microtube targets. Our findings in this work are supported by a systematic particle-in-cell modeling investigation of ion acceleration using 2D simulations with various structure dimensions. Simulations reproduce the experimental data with good agreement, including the observation of the optimum tube diameter, and reveal that the laser is shuttered by the plasma filling the tubes, explaining why the ion cutoff energy was not increased in this regime.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2019

Numerical Simulation and Experimental Validation of Nondendritic Structure Formation in Magnesium Alloy Under Oscillation and Ultrasonic Vibration

In this study, the formation of nondendritic structures in the primary phase of magnesium alloy solidified under oscillation and ultrasonic vibration was investigated by numerical simulation and experimentally. The growth and motion of a dendrite during solidification was simulated by a combination of the lattice Boltzmann method and the phase-field method. The simulation and experimental results indicated that higher oscillation amplitudes and acoustic streaming made the microstructures change from dendritic to nondendritic in the α-Mg primary phases. A sufficient shear stress and an appropriate flow time in the barrel should be satisfied when a given inclined angle is selected. The effects of the flow and the thermal field on the nucleation behavior, the constitutional undercooling, and the growth morphology are also discussed. It was found that a high shear rate and high turbulence can help homogenize the temperature and the concentration fields and collide and rotate the α-Mg primary phases.

36 MATERIALS SCIENCE↗
NASA NTRS1997

Probabilistic Evaluation of Bolted Joints in Polymer Matrix Composites

Computational methods are described to probabilistically simulate fracture in bolted composite structures. Progressive fracture is simulated via an innovative approach independent of stress intensity factors and fracture toughness. The effect on structure damage of design variable uncertainties is quantified. The Fast Probability Integrator is used to assess the scatter in the composite structure response before and after damage. Sensitivity of the response to design variables is evaluated. The methods are demonstrated for bolted joint polymer matrix composite panels under end loads. The effects of fabrication process are included in the simulation of damage in the bolted panel. The results show that the most effective way to reduce the end displacement at fracture is to control the load and ply thickness.

Chamis, C. C.↗
DOE OSTI2020

Development and Validation of a Simulation Tool to Predict the Combined Structural, Electrical, Electrochemical and Thermal Responses of Automotive Batteries

Battery safety is one of the most important design factors of electrified vehicles (EVs) since battery failure may lead to a catastrophic consequence. To improve the safety of battery systems in EVs, battery behaviors under various abuse conditions should be well understood. One way to examine battery response under extreme conditions is to conduct abuse tests in different scenarios. They can provide first-hand pass/fail information on battery safety, but they are usually expensive and time-consuming, which is inconvenient during design iterations and optimization. On the other hand, computational modeling can predict the battery behaviors in a more efficient and cost-effective way, which becomes an important tool to evaluate battery safety. Modeling battery responses and failure is challenging since it involves multiple coupled physical processes, needs to consider complex material properties and is computationally intensive. This project aims to develop a practical simulation tool to predict the combined structural, electrical, electrochemical, and thermal responses of automotive batteries to crash-induced crush and short circuit and validate it for conditions relevant to automotive crash. Advanced material constitutive models will be utilized to capture the mechanical response of cell components. Moreover, methods will be developed to reduce computational complexity of the model and allow battery safety simulations to extend from cell-level to module/pack level with affordable computational resources.

25 ENERGY STORAGE↗
DOE OSTI2025

MetaHeuristic Feature Selection for Energy Group Optimization and Analysis

Energy discretization is a crucial component of deterministic neutron transport simulations. Metaheuristic (MH) optimizers are effective algorithms to determine group structures that maximize both solution accuracy and computational efficiency. This project establishes a framework for optimizing group structures for PARTISN simulations using the Python library MEALPY. Group structure optimization is formulated as a binary feature selection problem, and results are investigated with permutation and material importance techniques to determine physically relevant energy bounds. We conclude that MH optimizers find group structures that drastically improve flux calculations while preserving k-effective accuracy. Further, we find that individual energy bounds are not necessarily physically relevant, but rather specific energy ranges are.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

A sharp interface Lagrangian-Eulerian method for rigid-body fluid-structure interaction

Herein this paper introduces a sharp interface method to simulate fluid-structure interaction (FSI) involving rigid bodies immersed in viscous incompressible fluids. The capabilities of this methodology are benchmarked using a range of test cases and demonstrated using large-scale models of biomedical FSI. The numerical approach developed herein, which we refer to as an immersed Lagrangian-Eulerian (ILE) method, integrates aspects of partitioned and immersed FSI formulations by solving separate momentum equations for the fluid and solid subdomains, as in a partitioned formulation, while also using non-conforming discretizations of the dynamic fluid and structure regions, as in an immersed formulation. A simple Dirichlet-Neumann coupling scheme is used, in which the motion of the immersed solid is driven by fluid traction forces evaluated along the fluid-structure interface, and the motion of the fluid along that interface is constrained to match the solid velocity and thereby satisfy the no-slip condition. To develop a practical numerical method, we adopt a penalty approach that approximately imposes the no-slip condition along the fluid-structure interface. In the coupling strategy, a separate discretization of the fluid-structure interface is tethered to the volumetric solid mesh via stiff spring-like penalty forces. Our fluid-structure coupling scheme relies on an immersed interface method (IIM) for discrete geometries, which enables the accurate determination of both velocities and stresses along complex internal interfaces. Numerical methods for FSI can suffer from instabilities related to the added mass effect, but computational tests indicate that the methodology introduced here remains stable for selected test cases across a broad range of solid-fluid mass density ratios, including extremely small, nearly equal, equal, and large density ratios. Biomedical FSI demonstration cases include results obtained using this method to simulate the dynamics of a bileaflet mechanical heart valve in a pulse duplicator, and to model transport of blood clots in a patient-averaged anatomical model of the inferior vena cava.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1995

Active Piezoelectric Structures for Tip Clearance Management Assessed

Managing blade tip clearance in turbomachinery stages is critical to developing advanced subsonic propulsion systems. Active casing structures with embedded piezoelectric actuators appear to be a promising solution. They can control static and dynamic tip clearance, compensate for uneven deflections, and accomplish electromechanical coupling at the material level. In addition, they have a compact design. To assess the feasibility of this concept and assist the development of these novel structures, the NASA Lewis Research Center developed in-house computational capabilities for composite structures with piezoelectric actuators and sensors, and subsequently used them to simulate candidate active casing structures. The simulations indicated the potential of active casings to modify the blade tip clearance enough to improve stage efficiency. They also provided valuable design information, such as preliminary actuator configurations (number and location) and the corresponding voltage patterns required to compensate for uneven casing deformations. An active ovalization of a casing with four discrete piezoceramic actuators attached on the outer surface is shown. The center figure shows the predicted radial displacements along the hoop direction that are induced when electrostatic voltage is applied at the piezoceramic actuators. This work, which has demonstrated the capabilities of in-house computational models to analyze and design active casing structures, is expected to contribute toward the development of advanced subsonic engines.

Source record↗
NASA NTRS1994

Design of orbital debris shields for oblique hypervelocity impact

A new impact debris propagation code was written to link CTH simulations of space debris shield perforation to the Lagrangian finite element code DYNA3D, for space structure wall impact simulations. This software (DC3D) simulates debris cloud evolution using a nonlinear elastic-plastic deformable particle dynamics model, and renders computationally tractable the supercomputer simulation of oblique impacts on Whipple shield protected structures. Comparison of three dimensional, oblique impact simulations with experimental data shows good agreement over a range of velocities of interest in the design of orbital debris shielding. Source code developed during this research is provided on the enclosed floppy disk. An abstract based on the work described was submitted to the 1994 Hypervelocity Impact Symposium.

Fahrenthold, Eric P.↗