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 919 records · Page 51

The Simons Observatory: Modeling Optical Systematics in the Large Aperture Telescope

We present geometrical and physical optics simulation results for the Simons Observatory Large Aperture Telescope. This work was developed as part of the general design process for the telescope; allowing us to evaluate the impact of various design choices on performance metrics and potential systematic effects. The primary goal of the simulations was to evaluate the final design of the reflectors and the cold optics which are now being built. We describe non-sequential ray tracing used to inform the design of the cold optics, including absorbers internal to each optics tube. We discuss ray tracing simulations of the telescope structure that allow us to determine geometries that minimize detector loading and mitigate spurious near-field effects that have not been resolved by the internal baffling. We also describe physical optics simulations, performed over a range of frequencies and field locations, that produce estimates of monochromatic far field beam patterns which in turn are used to gauge general optical performance. Finally, we describe simulations that shed light on beam sidelobes from panel gap diffraction.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Recent Hydrodynamic Modeling Enhancements in OpenFAST: Preprint

To improve the hydrodynamic modeling capabilities of the offshore wind design and modeling tool OpenFAST, a restructuring of the HydroDyn hydrodynamics module was undertaken with several new features implemented. The generation of the wave field is now separated from HydroDyn into a new module called SeaState. Unlike previous versions of HydroDyn, which precomputed the wave kinematics at the undisplaced positions of the hydrodynamic nodes, the new SeaState module computes the complete time history of the wave field at the vertices of a user-defined wave grid. During the simulation, the wave kinematics at any point within the grid can be efficiently interpolated, allowing the strip-theory wave loads to be evaluated based on the wave kinematics at the displaced structure position. A phase correction was also implemented for the potential-flow wave-exciting loads to account for large structure displacements. Several wave-stretching methods were implemented for the strip-theory solution along with a wave-load redistribution method that ensures smooth variation of the nodal loads as the nodes of a discrete hydrodynamic mesh enter and exit the water. The load smoothing is included to avoid the excitation of unphysical high-frequency structural vibration during hydroelastic simulations. The MacCamy-Fuchs diffraction correction for the strip-theory solution and the capability to insert a constrained NewWave into extreme stochastic sea states with directional spreading were also added. The present article documents the formulations of the new improvements to HydroDyn with example applications and numerical results.

constrained wave↗

Probing the Magnetic Field Structure in Sgr A* on Black Hole Horizon Scales with Polarized Radiative Transfer Simulations

Magnetic fields are believed to drive accretion and relativistic jets in black hole accretion systems, but the magnetic field structure that controls these phenomena remains uncertain. We perform general relativistic (GR) polarized radiative transfer of time-dependent three-dimensional GR magnetohydrodynamical simulations to model thermal synchrotron emission from the Galactic Center source Sagittarius A* (Sgr A*). We compare our results to new polarimetry measurements by the Event Horizon Telescope (EHT) and show how polarization in the visibility (Fourier) domain distinguishes and constrains accretion flow models with different magnetic field structures. These include models with small-scale fields in disks driven by the magnetorotational instability as well as models with large-scale ordered fields in magnetically arrested disks. We also consider different electron temperature and jet mass-loading prescriptions that control the brightness of the disk, funnel-wall jet, and Blandford–Znajek-driven funnel jet. Our comparisons between the simulations and observations favor models with ordered magnetic fields near the black hole event horizon in Sgr A*, though both disk- and jet-dominated emission can satisfactorily explain most of the current EHT data. We also discuss how the black hole shadow can be filled-in by jet emission or mimicked by the absence of funnel jet emission. We show that stronger model constraints should be possible with upcoming circular polarization and higher frequency (349 GHz) measurements.

79 ASTRONOMY AND ASTROPHYSICS↗

Direct numerical simulation of mechanical response in synthetic additively manufactured microstructures

Additive manufacturing (AM) processes for metals can yield as-built microstructures that vary significantly from their cast or wrought counterparts. These microstructural variations can in turn, have profound effects on the properties of a component. Here, a modeling methodology is presented to investigate microstructurally-influenced mechanical response in additively manufactured structures via direct numeral simulation. Three-dimensional, synthetic voxelized microstructures are generated by kinetic Monte Carlo (kMC) additive manufacturing process simulations performed at four scan speeds to create a thin-wall cylindrical geometry notionally constructed using a concentric-pathed directed energy deposition AM process. The kMC simulations utilize a steady-state molten pool geometry that is held constant throughout the study. Resultant microstructures are mapped onto a highly-refined conformal finite-element mesh of a part geometry. A grain-scale anisotropic crystal elasticity model is then used to represent the constitutive response of each grain. Here, the response of the structure subjected to relatively simple load conditions is studied in order to provide understanding of both the influence of AM processing on microstructure as well as the microstructure's influence on the macroscale mechanical response.

42 ENGINEERING↗

Electron diffraction intensity from single crystal silicon in a photoinjector

A method includes simulating diffraction in a transmission geometry of relativistic electron bunches from a crystallographic structure of a crystal thereby simulating diffraction of the relativistic electron bunches into a plurality of Bragg peaks. The method includes selecting a range of angles between a direction of propagation of the relativistic electron bunches and a normal direction of crystal including an angle at which a diffraction portion is maximized. The method includes sequentially accelerating a plurality of physical electron bunches to relativistic energies toward a physical crystal having the crystallographic structure and diffracting the plurality of physical electron bunches off the physical crystal at different angles and measuring the diffraction portion into the respective Bragg peak at the different angles. The method includes selecting a final angle based on the measured diffraction portion into the respective Bragg peak at the different angles and generating a pulse of light.

Graves, William↗

Collaborative Research: Tomographic imaging of laser-plasma structures

The interaction of intense short laser pulses with ionized gases, or plasmas, underlies many applications such as acceleration of elementary particles, production of energy by laser fusion, generation of x-ray and far-infrared “terahertz” pulses for medical and materials probing, remote sensing of explosives and pollutants, and generation of guide stars. Such laser-plasma interactions create tiny electron density structures (analogous to the wake behind a boat) inside the plasma in the shape of waves, bubbles and filaments that move at the speed of light, and evolve as they propagate. Prior to recent work by the PI of this proposal, detailed knowledge of such structures came exclusively from intensive computer simulations. Now “snapshots” of these elusive, light-velocity structures can be taken in the laboratory using dynamic variant of holography, the technique used to produce ID cards and DVDs, and dynamic variant of tomography, the technique used in medicine to image internal bodily organs. These fast visualization techniques are important for understanding, improving and scaling the above-mentioned applications of laser-plasma interactions. In this project, we accomplished three things: 1) We took holographic pictures of a laser-driven plasma-wave in the act of accelerating electrons to high energy, and used computer simulations to understand the pictures. 2) Using results from this experiment to optimize the performance of the accelerator, and the brightness of x-rays that it emits. These x-rays will be useful for medical and materials science applications. 3) We made technical improvements to the holographic technique that enables us to see finer details in the recorded pictures. Four refereed journal papers were published, and two students earned PhDs and moved on to scientific careers in US National Laboratories based on their work under this project.

43 PARTICLE ACCELERATORS↗

Molecular Dynamics Simulation on the Effect of MgO/Al{sub 2}O{sub 3} Ratio on Structure and Properties of Blast Furnace Slag Under Different Basicity Conditions

The SiO{sub 2}-Al{sub 2}O{sub 3}-CaO-MgO is the basic structural system of blast furnace slag and the composition directly affects the performance of the slag. Molecular dynamics simulations were carried out to analyze the local structure, structural unit, bond angle, transport properties, and enthalpy of the slag with the increase of MgO/Al{sub 2}O{sub 3} mass ratio under different basicity conditions. It was found that the change of MgO/Al{sub 2}O{sub 3} ratio does not affect the short-range order structure, but it will reduce the overall stability of the network structure. Through the analysis of the structural unit of the slag, it was found that the polymerization degree of the system decreases with the increase of MgO/Al{sub 2}O{sub 3} ratio, indicating that the complex network structure of the system is partially depolymerized. Besides, the self-diffusion coefficients of each ion were concluded and the magnitudes were observed to be in the following order: Mg{sup 2+}>Ca{sup 2+}>Al{sup 3+}>O{sup 2-}>Si{sup 4+}. The diffusivity of the slag increases with the increase of MgO/Al{sub 2}O{sub 3} ratio and the corresponding viscosity decreases. And the enthalpy of the slag is increased with the increase of MgO/Al{sub 2}O{sub 3} ratio or basicity, therefore the fuel consumption should also be properly adjusted during the operation of blast furnace to ensure the target temperature.

36 MATERIALS SCIENCE↗

Molecular dynamics simulations of uranyl and plutonyl cations in a task-specific ionic liquid

Ionic liquids (ILs) are a unique class of solvents with potential applications in advanced separation technologies relevant to the nuclear industry. ILs are salts with low melting points and a wide range of tunable physical properties, such as viscosity, hydrophobiciy, conductivity, and liquidus range. ILs have negligible vapor pressure, are often non-flammable, and can have high thermal stability and a wide electrochemical window, making them attractive for use in separations processes relevant to the nuclear industry. Metal salts generally have a low solubility in ILs; however, by incorporating new functional groups into the IL cation or anion that promote complexation with the metal, the solubility can be greatly increased. One such task-specific ionic liquid (TSIL) is 1-carboxy-N, N, N-trimethylglycine bis(trifluoromethylsulfonyl)imide ([Hbet][Tf 2 N]). Water, which is detrimental for electrochemical separations, is a common impurity in ILs and can coordinate with actinyl cations, particularly in ILs containing only weakly coordinating components. Understanding the behavior of actinides in TSIL/water mixtures on a molecular level is vital for designing improved separations processes. Classical molecular dynamics simulations of uranyl(VI) and plutonyl(VI) in 1-ethyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide ([EMIM][Tf 2 N]) with deprotonated Hbet (betaine) and water have been performed to understand the coordination and dynamics of the actinyl cations. We find that betaine is a much stronger ligand than water and prefers to coordinate the metal in a bidentate manner. Potential of mean force simulations yield a relative free energy for betaine coordination of approximately -120 to -90 kJ/mol in mixtures with water. As the amount of betaine coordinated to the actinide increases, the diffusion coefficient of the actinyl cation decreases. Moreover, the betaine ligand is able to bridge between two metal centers, resulting in dimeric complexes with actinide–actinide distances of ~5 Å. Potential of mean force simulations show that these structures are stable, with relative free energies of up to -40 kJ/mol. The crystal structure for [(UO 2 ) 2 (bet) 6 (H 2 O) 2 ][Tf 2 N] 4 shows that the betaine bridges between two uranium atoms to form dimeric complexes similar to those found in our simulations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Lattice Design and Advanced Modeling to Guide the Design of High-Performance Lightweight Structural Materials

Lightweight structural materials are required to increase the mobility of fission batteries. The materials must feature a robust combination of mechanical properties to demonstrate structural resilience. The primary objective of this project is to produce lightweight structural materials whose strength-to-weight ratios exceed those of the current widely used structural materials such as 316L stainless steels (316L SS). To achieve this, advanced modeling and simulation tools were employed to design lattice structures with different lattice parameters and different lattice types. A process was successfully developed for transforming lattice-structures models into Multiphysics Object Oriented Simulation Environment (MOOSE) inputs. Finite element modeling (FEM) was used to simulate the uniaxial tensile testing of the lattice-structured parts to investigate the stress distribution at a given displacement. The preliminary results showed that the lattice-structured sample displayed a lower Young’s modulus in comparison with the solid material and that the unit cell size of the lattice had a minimal effect. The novelty here is to apply up-front modeling to determine the best structure for the application before actually producing the sample. The approach of using modeling as a guiding tool for preliminary material design can significantly save time and cost for material development.

36 MATERIALS SCIENCE↗

Integrated Simulation Framework for Additively Manufactured Ti-6Al-4V: Melt Pool Dynamics, Microstructure, Solid-State Phase Transformation, and Microelastic Response

In this work, to accelerate the establishment of fundamental understanding of the additive manufacturing (AM) process and its influence on microstructural evolution and related properties, we develop a multiphysics and multiscale modeling framework that integrates: (1) a high-fidelity powder-scale three-dimensional simulation of transient heat transfer and melt flow dynamics, (2) cellular automaton simulation of solidification grain structure and texture, (3) phase-field modeling of precipitation and dissolution of second-phase precipitate during repeated thermal cycles, and (4) microstructure-based micro- and mesoscopic elastic response calculation. Using Ti-6Al-4V as a model system, we demonstrate the application of the integrated framework to simulate complex microstructure evolution during a single-track laser powder bed fusion process and the associated mechanical response. Our modeling framework successfully captures the solidification β grain structure as a function of laser power and scanning speed, α precipitation upon subsequent cooling with different rates, and elastic response of the resulting (α + β) two-phase microstructure. The key features of solidification and second-phase precipitate microstructures, and their dependence on processing parameters, agree well with existing experimental observations. Finally, the established modeling framework is generally applicable to other metallic materials fabricated by AM.

36 MATERIALS SCIENCE↗

Visualization of Solid-State Synthesis for Chalcogenide Na Superionic Conductors by in-situ Neutron Diffraction

Chalcogenide superionic sodium (Na) conductors are great potential as solid electrolytes (SEs) in all-solid-state Na batteries with advantages of high energy density and safety, and cost effectiveness. For solid Na-ion conductors, their crystal structures and ionically conductive properties are strongly influenced by the synthetic approaches and processing parameters. Thus, understanding the synthesis process is essential to control the structures and phases and thereby yields to Na-ion conductors with desirable properties. Thanks to the high-flux and deep-penetrating time-of-flight neutron diffraction (ND), we employed in situ experiments to track real-time structural changes of two chalcogenide SEs (Na 3 SbS 4 and Na 3 SbS 3.5 Se 0.5 ) during the solid-state synthesis. For these two conductors, the ND results reveal a fast one-step reaction for the synthesis and the molten process when heating up, and the recrystallization as well as the cubic-to-tetragonal phase transition up on cooling. Moreover, Se-doping is found to influence the reaction temperatures, lattice parameter and structure stability based on neutron experimental observations and theoretical simulation. This work presents a detailed structural study using in situ neutron diffraction technology for the solid synthesis process of chalcogenide Na-ion conductors, beneficial for the design and synthesis of new solid-state conductors.

25 ENERGY STORAGE↗

Distance‐based reconstruction of protein quaternary structures from inter‐chain contacts

Abstract Predicting the quaternary structure of protein complex is an important problem. Inter‐chain residue‐residue contact prediction can provide useful information to guide the ab initio reconstruction of quaternary structures. However, few methods have been developed to build quaternary structures from predicted inter‐chain contacts. Here, we develop the first method based on gradient descent optimization (GD) to build quaternary structures of protein dimers utilizing inter‐chain contacts as distance restraints. We evaluate GD on several datasets of homodimers and heterodimers using true/predicted contacts and monomer structures as input. GD consistently performs better than both simulated annealing and Markov Chain Monte Carlo simulation. Starting from an arbitrarily quaternary structure randomly initialized from the tertiary structures of protein chains and using true inter‐chain contacts as input, GD can reconstruct high‐quality structural models for homodimers and heterodimers with average TM‐score ranging from 0.92 to 0.99 and average interface root mean square distance from 0.72 Å to 1.64 Å. On a dataset of 115 homodimers, using predicted inter‐chain contacts as restraints, the average TM‐score of the structural models built by GD is 0.76. For 46% of the homodimers, high‐quality structural models with TM‐score ≥ 0.9 are reconstructed from predicted contacts. There is a strong correlation between the quality of the reconstructed models and the precision and recall of predicted contacts. Only a moderate precision or recall of inter‐chain contact prediction is needed to build good structural models for most homodimers. Moreover, GD improves the quality of quaternary structures predicted by AlphaFold2 on a Critical Assessment of Techniques for Protein Structure Prediction–Critical Assessments of Predictions of Interactions dataset.

59 BASIC BIOLOGICAL SCIENCES↗

Data-Centric Development of Lignin Structure–Solubility Relationships in Deep Eutectic Solvents Using Molecular Simulations

Lignin is a natural source of aromatic chemicals with significant potential as an abundant, renewable feedstock for value-added products. Deep eutectic solvents (DES)–solvents composed of a hydrogen bond donor (HBD) and acceptor (HBA) in varying ratios–have emerged as a highly tunable class of solvents for lignin solubilization. However, the variety of possible DES compositions and limited molecular-scale understanding of lignin solubility makes solvent selection a challenge without laborious trial-and-error experimentation. To address these challenges, we use classical molecular dynamics (MD) simulations to study the interactions of lignin model compounds with various DES–water systems. Quantitative parameters (descriptors) were calculated by postprocessing the MD results and used to train a regression model that predicts experimentally determined solubilities of lignin model compounds. This approach revealed that the most important descriptors of solubility are the system temperature, solute hydrophilicity, and metrics quantifying hydrogen bonding. Maximizing the interactions between solute–HBD (hydrophobic group), water–HBD (hydrophilic group), and water–HBA molecules led to the highest model compound solubility. Our results support a hydrotropic mechanism in which extensive DES–water hydrogen bonding and favorable HBD interactions with the solute promote high solubility. We applied the regression model derived using model compounds to predict the solubility of representative lignin oligomers. The model predicted lignin oligomers’ solubilities in good agreement with experiments, indicating that the simulations of model compounds can be extended to predict the solubility of larger lignin compounds across a range of solvent compositions and temperatures. Furthermore, these findings provide new molecular-scale insight into lignin solubilization mechanisms and a new method for computationally screening potential solvent systems for lignin valorization.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Predicting the preferred morphology of hexagonal boron nitride domain structure on nickel from ReaxFF-based molecular dynamics simulations

An understanding of the nucleation and growth of hexagonal boron nitride (hBN) on nickel substrates is essential to its development as a functional material. In particular, fundamental insights into the formation of the hexagonal lattices with alternating boron (B) and nitrogen (N) atoms could be exploited to control hBN lattice morphologies for targeted applications. Here, the preferred shapes and edge configurations of atomically smooth hBN on Ni(111) were investigated using molecular dynamics (MD) simulations, along with reactive force field (ReaxFF) developed to represent the Ni/B/N system and the lattice-building B–N bond formation. The obtained hBN lattices, from different B : N feed ratios, are able to confirm that hBN domain geometries can indeed be tuned by varying thermodynamic parameters (i.e., chemical potentials of N and B) – a finding that has only been predicted using quantum mechanical theories. Here, we also showed that the nitrogen fed to the system plays a more crucial role in dictating the size of hBN lattices. With an increase of the relative N content, the simulated hBN domain shapes also transition from equilateral triangles to hexagons, again, consistent with the anticipation based on Density Functional Theory (DFT) calculations. Hence, a plausible approach to acquire a desired hBN nanostructure depends on careful control over the synthesis conditions, which now can benefit from reliable molecular simulations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗