Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Solution structure”

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 253 records · Page 14

Impact of aqueous solution pH on network structure of corrosion–induced surface layers of boroaluminosilicate glass

Understanding of the long-term stability of vitrified nuclear waste glass is critical to its safe storage, but knowledge gaps remain regarding the corrosion mechanism, especially with respect to surface alteration layers whose presence changes the observed corrosion rates. Here, we report the aqueous corrosion of a model boroaluminosilicate glass at ground water relevant pHs as investigated by vibrational spectroscopy to understand the network structure of the alteration layer that we hypothesize controls the corrosion rates. Both specular reflection infrared spectroscopy and infrared variable angle spectroscopic ellipsometry were used to elucidate details of the structure of the alteration layers formed at pH 7 and pH 9. After 2312 days, the thickness of the alteration layer is found to be ~2000 nm for corrosion at pH 7 and ~300 nm for corrosion at pH 9. Here, the interpretation of the vibrational spectra suggests that alteration layers formed at pH 7 and pH 9 both exhibit a porous network with "silica-like" structures, but with differences in the bond lengths and bond angles in the glass network. Moreover, the relative abundance of the hydrous species within the porous network appears to be dependent on pH. Additionally, the spectroscopic data suggest that the B and Ca trapped in the alteration layer during the secondary corrosion of the pH 9 pre-corroded sample in pH 3 solution are hydrated within the nano-porous layer and not bound to the silicate network. However, the precise chemical form remains to be determined. Based on these findings, we propose that the combination of the network structure of the alteration layer and the relative abundance of hydrous species within the alteration layer act to control the transport of ionic species dissolved from the reaction front to the aqueous solution outside the glass.

36 MATERIALS SCIENCE↗

Impact of Solution Chemistry on Growth and Structural Features of Mo-Substituted Spinel Iron Oxides

The effect of crystallizing solution chemistry on the chemistry of subsequently as-grown materials was investigated for Mo-substituted iron oxides prepared by thermally activated co-precipitation. In the presence of Mo ions, we find that varying the oxidation state of the iron precursor from Fe(II) to Fe(III) causes a progressive loss of atomic long-range order with the stabilization of 2-4 nm particles for the sample prepared with Fe(III). The oxidation state of the Fe precursor also affects the distribution of Fe and Mo cations within the spinel structure. Increasing the Fe precursor oxidation state gives decreased Fe-ion occupation and increased Mo-ion occupation of tetrahedral sites, as revealed by the extended X-ray absorption fine structure. The stabilization of Mo within tetrahedral sites appears to be unexpected, considering the octahedral preferred coordination number of Mo(VI). The analysis of the atomic structure of the sample prepared with Fe(III) indicates a local ordering of vacancies and that the occupation of tetrahedral sites by Mo induces a contraction of the interatomic distances within the polyhedra as compared to Fe atoms. Moreover, the occupancy of Mo into the thermodynamic site preference of a Mo dopant in Fe 2 O 3 assessed by density functional theory calculations points to a stronger preference for Mo substitution at octahedral sites. Hence, we suggest that the synthetized compound is thermodynamically metastable, that is, kinetically trapped. Such a state is suggested to be a consequence of the tetrahedral site occupation by Mo ions. The population of these sites, known to be reactive sites enabling particle growth, is concomitant with the stabilization of very small particles. Furthermore, we confirmed our hypothesis by using a blank experiment without Mo ions, further supporting the impact of tetrahedral Mo ions on the growth of iron oxide nanoparticles. Our findings provide new insights into the relationships between the Fechemistry of the crystallizing solution and the structural features of the as-grown Mo-substituted Fe-oxide materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Genarris 2.0: A Random Structure Generator for Molecular Crystals

Genarris is an open source Python package for generating random molecular crystal structures with physical constraints for seeding crystal structure prediction algorithms and training machine learning models. Here we present a new version of the code, containing several major improvements. A MPI-based parallelization scheme has been implemented, which facilitates the seamless sequential execution of user-defined workflows. A new method for estimating the unit cell volume based on the single molecule structure has been developed using a machine-learned model trained on experimental structures. A new algorithm has been implemented for generating crystal structures with molecules occupying special Wyckoff positions. A new hierarchical structure check procedure has been developed to detect unphysical close contacts efficiently and accurately. New intermolecular distance settings have been implemented for strong hydrogen bonds. To demonstrate these new features, we study two specific cases: benzene and glycine. Genarris finds the experimental structures of the two polymorphs of benzene and the three polymorphs of glycine. Program summary Program Title: Genarris 2.0 Program Files doi: http://dx.doi.org/10.17632/grx6mz4pjn.1 Licensing provisions: BSD-3 Clause Programming language: Python, C External routines/libraries: Spglib, ASE, pymatgen, SciPy, mpi4py, scikit-learn, PyTorch, FHI-aims. Nature of problem: Molecular crystal structure prediction. Solution method: Genarris 2.0 generates molecular crystal structures over the 230 space groups, on general and special Wyckoff positions, using physical constraints. Down-sampling of the generated structures may be performed subsequently, based on molecular crystal packing descriptors and an unsupervised machine learning algorithm. Lastly, ab initio structure relaxation may be performed for the final pool. Depending on the user-defined workflow implemented, Genarris may be used to generate diverse molecular crystal datasets to seed evolutionary algorithms or to train machine learning algorithms or as a standalone crystal structure prediction method. Restrictions: For crystal structure generation, the molecule of interest must be semi-rigid with no bond rotational degrees of freedom. Unusual features: Genarris 2.0 is a highly distributed program, making use of MPI for Python parallelization. The user has the ability to design and implement workflows by executing a user-defined list of procedures. Genarris 2.0 offers new features including a machine learning model for estimating the molecular volume in the solid state from the single molecule structure, structure generation in special Wyckoff positions of space groups, hierarchical structure checks including rigorous treatment of non-orthogonal structures, and clustering and down-selection workflows combining first principles simulations with machine learning. (C) 2020 Elsevier B.V. All rights reserved.

Crystal structure prediction↗

Sensitivity analysis of the optimal solution obtained from the structural damage indentification process

An identification procedure proposed by Shen and Taylor to determine the crack characteristics (location and size of the crack) from dynamic measurements is tested. This procedure was based on minimization of either the 'mean-square' measure of difference between measurement data (natural frequencies and mode shapes) and the corresponding predictions obtained from the computational model. The procedure is tested for simulated damage in the form of symmetric cracks in a simply-supported Bernoulli-Euler beam. The sensitivity of the solution of damage identification problems to the values of parameters that characterize damage is discussed. A sensitivity formula is derived.

Shen, Mo-How H.↗

Atomistic structure of the SARS-CoV-2 pseudoknot in solution from SAXS-driven molecular dynamics

SARS-CoV-2 depends on -1 programmed ribosomal frameshifting (-1 PRF) to express proteins essential for its replication. The RNA pseudoknot stimulating -1 PRF is thus an attractive drug target. However, the structural models of this pseudoknot obtained from cryo-EM and crystallography differ in some important features, leaving the pseudoknot structure unclear. We measured the solution structure of the pseudoknot using small-angle X-ray scattering (SAXS). The measured profile did not agree with profiles computed from the previously solved structures. Beginning with each of these solved structures, we used the SAXS data to direct all atom molecular dynamics (MD) simulations to improve the agreement in profiles. In all cases, this refinement resulted in a bent conformation that more closely resembled the cryo-EM structures than the crystal structure. Applying the same approach to a point mutant abolishing -1 PRF revealed a notably more bent structure with reoriented helices. This work clarifies the dynamic structures of the SARS-CoV-2 pseudoknot in solution.

60 APPLIED LIFE SCIENCES↗

Real-space solution to the electronic structure problem for nearly a million electrons

We report a Kohn–Sham density functional theory calculation of a system with more than 200 000 atoms and 800 000 electrons using a real-space high-order finite-difference method to investigate the electronic structure of large spherical silicon nanoclusters. Our system of choice was a 20 nm large spherical nanocluster with 202 617 silicon atoms and 13 836 hydrogen atoms used to passivate the dangling surface bonds. To speed up the convergence of the eigenspace, we utilized Chebyshev-filtered subspace iteration, and for sparse matrix–vector multiplications, we used blockwise Hilbert space-filling curves, implemented in the PARSEC code. For this calculation, we also replaced our orthonormalization + Rayleigh–Ritz step with a generalized eigenvalue problem step. We utilized all of the 8192 nodes (458 752 processors) on the Frontera machine at the Texas Advanced Computing Center. We achieved two Chebyshev-filtered subspace iterations, yielding a good approximation of the electronic density of states. Our work pushes the limits on the capabilities of the current electronic structure solvers to nearly 106 electrons and demonstrates the potential of the real-space approach to efficiently parallelize large calculations on modern high-performance computing platforms.

Chemistry↗

A family of solution algorithms for nonlinear structural analysis based on relaxation equations

A family of hierarchical algorithms for nonlinear structural equations are presented. The algorithms are based on the Davidenko-Branin type homotopy and shown to yield consistent hierarchical perturbation equations. The algorithms appear to be particularly suitable to problems involving bifurcation and limit point calculations. An important by-product of the algorithms is that it provides a systematic and economical means for computing the stepsize at each iteration stage when a Newton-like method is employed to solve the systems of equations. Some sample problems are provided to illustrate the characteristics of the algorithms.

Park, K. C.↗

Element-by-element Solution Procedures for Nonlinear Structural Analysis

Element-by-element approximate factorization procedures are proposed for solving the large finite element equation systems which arise in nonlinear structural mechanics. Architectural and data base advantages of the present algorithms over traditional direct elimination schemes are noted. Results of calculations suggest considerable potential for the methods described.

Hughes, T. J. R.↗

Determining structural performance

An overview is given of the methods and concepts developed to enhance and predict structural dynamic characteristics of advanced aeropropulsion systems. Aeroelasticity, Vibration Control, Dynamic Systems, and Computational Structural Methods are four disciplines that make up the research program at NASA/Lewis Research Center. The Aeroelasticity program develops analytical and experimental methods to minimize flutter and forced vibration of aerospace propulsion systems. Both frequency domain and time domain methods have been developed for applications on the turbofan, turbopump, and advanced turboprop. To improve life and performance, the Vibration Control program conceives, analyzes, develops, and demonstrates new methods to control vibrations in aerospace systems. Active and passive vibration control is accomplished with electromagnetic dampers, magnetic bearings, and piezoelectric crystals to control rotor vibrations. The Dynamic Systems program analyzes and verifies the dynamics of interacting systems, as well as develops concepts and methods for high-temperature dynamic seals. The Computational Structural Methods program uses computer science to improve solutions of structural problems.

Ernst, Michael A.↗

Looking for O(N) Navier-Stokes solutions on non-structured meshes

Multigrid methods are good candidates for the resolution of the system arising in numerical fluid dynamics. However, the question is to know if those algorithms which are efficient for the Poisson equation on structured meshes will still apply well to the Euler and Navier-Stokes equations on unstructured meshes. The study of elliptic problems leads us to define the conditions where a full multigrid strategy has O(N) complexity. The aim of this paper is to build a comparison between the elliptic theory and practical CFD problems. First, as an introduction, we will recall some basic definitions and theorems applied to a model problem. The goal of this section is to point out the different properties that we need to produce an FMG algorithm with O(N) complexity. Then, we will show how we can apply this theory to the fluid dynamics equations such as Euler and Navier-Stokes equations. At last, we present some results which are 2nd-order accurate and some explanations about the behavior of the FMG process.

Morano, Eric↗

Looking for O(N) Navier-Stokes solutions on non-structured meshes

Multigrid methods are good candidates for the resolution of the system arising in Numerical Fluid Dynamics. However, the question is to know if those algorithms which are efficient for the Poissan equation on structured meshes will still apply well to the Euler and Navier-Stokes equations on unstructured meshes. The study of elliptic problems leads us to define the conditions where a Full Multigrid strategy has O(N) complexity. The aim of this paper is to build a comparison between the elliptic theory and practical CFD problems. First, as an introduction, we will recall some basic definitions and theorems applied to a model problem. The goal of this section is to point out the different properties that we need to produce an FMG algorithm with O(N) complexity. Then, we will show how we can apply this theory to the fluid dynamics equations such as Euler and Navier-Stokes equations. At last, we present some results which are 2nd-order accurate and some explanations about the behavior of the FMG process.

Morano, Eric↗

Structural Design Exploration of an Electric Powered Multi-Propulsor Wing Configuration

Advancements in aircraft electric propulsion may enable an expanded operational envelope for electrically powered vehicles compared to their internal combustion engine counterparts. High aspect ratio wings provide additional lift and drag reduction for a proposed multi-propulsor design, however, the challenge is to reduce the weight of wing structures while maintaining adequate structural and aeroelastic margins. Design exploration using a conventional design-and-build philosophy coupled with a finite element method (FEM)-based design of experiments (DOE) strategy are presented to examine high aspect ratio wing structures that have spanwise distributed electric motors. Multiple leading-edge-mounted engine masses presented a challenge to design a wing within acceptable limits for dynamic and aeroelastic stability. Because the first four primary bending eigenmodes of the proposed wing structure are very sensitive to outboard motor placement, safety-of-flight requirements drove the need for multiple spars, rib attachments, and outboard structural reinforcements in the design. Global aeroelasticity became an increasingly important design constraint during the on-going design process, with outboard motor pod flutter ultimately becoming a primary design constraint. Designers successively generated models to examine stress, dynamics, and aeroelasticity concurrently. This research specifically addressed satisfying multi-disciplinary design criteria to generate fluid-structure interaction solution sets, and produced high aspect ratio primary structure designs for the NASA Scalable Convergent Electric Propulsion Technology and Operations Research (SCEPTOR) project in the Aeronautic Research Mission Directorate at NASA. In this paper, a dynamics-driven, quasi-inverse design methodology is presented to address aerodynamic performance goals and structural challenges encountered for the SCEPTOR demonstrator vehicle. These results are compared with a traditional computer aided design based approach.

Moore, James B.↗

Excess chemical potential of small solutes across water--membrane and water--hexane interfaces

The excess chemical potentials of five small, structurally related solutes, CH4, CH3F, CH2F2, CHF3, and CF4, across the water-glycerol 1-monooleate bilayer and water-hexane interfaces were calculated at 300, 310, and 340 K using the particle insertion method. The excess chemical potentials of nonpolar molecules (CH4 and CF4) decrease monotonically or nearly monotonically from water to a nonpolar phase. In contrast, for molecules that possess permanent dipole moments (CH3F, CH2F, and CHF3), the excess chemical potentials exhibit an interfacial minimum that arises from superposition of two monotonically and oppositely changing contributions: electrostatic and nonelectrostatic. The nonelectrostatic term, dominated by the reversible work of creating a cavity that accommodates the solute, decreases, whereas the electrostatic term increases across the interface from water to the membrane interior. In water, the dependence of this term on the dipole moment is accurately described by second order perturbation theory. To achieve the same accuracy at the interface, third order terms must also be included. In the interfacial region, the molecular structure of the solvent influences both the excess chemical potential and solute orientations. The excess chemical potential across the interface increases with temperature, but this effect is rather small. Our analysis indicates that a broad range of small, moderately polar molecules should be surface active at the water-membrane and water-oil interfaces. The biological and medical significance of this result, especially in relation to the mechanism of anesthetic action, is discussed.

NASA Center ARC↗

Analysis and Sizing for Transient Thermal Heating of Insulated Aerospace Vehicle Structures

An analytical solution was derived for the transient response of an insulated structure subjected to a simplified heat pulse. The solution is solely a function of two nondimensional parameters. Simpler functions of these two parameters were developed to approximate the maximum structural temperature over a wide range of parameter values. Techniques were developed to choose constant, effective thermal properties to represent the relevant temperature and pressure-dependent properties for the insulator and structure. A technique was also developed to map a time-varying surface temperature history to an equivalent square heat pulse. Equations were also developed for the minimum mass required to maintain the inner, unheated surface below a specified temperature. In the course of the derivation, two figures of merit were identified. Required insulation masses calculated using the approximate equation were shown to typically agree with finite element results within 10%-20% over the relevant range of parameters studied.

Blosser, Max L.↗

Influence of Kosmotrope and Chaotrope Salts on Water Structural Relaxation

The structural relaxation in water solutions of kosmotrope (structure maker) and chaotrope (structure breaker) salts, namely sodium chloride, potassium chloride, and cesium chloride, were studied through quasielastic neutron scattering measurements. We found that the collective dynamics relaxation time at the structure factor peak obtained using heavy water solutions shows a distinctively different behavior in the kosmotrope as opposed to the chaotrope solutions, increasing with the salt concentration in the former and decreasing in the latter. In both cases the trends are proportional to the concentration dependence of the relative viscosity of the solutions. Furthermore, these results indicate that kosmotropes and chaotropes influence the solution’s viscosity by impacting in opposite ways the hydrogen bond network of water, strengthening it in one case and softening it in the other.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Analytical Solutions for Seismic Fluid-Structure Interaction of Head-Supported Cylindrical Tanks

Fluid-structure interaction (FSI) of earthquake shaking is a design consideration in the civil, nuclear, chemical, and mechanical industries. Preliminary sizing of structures and components in these industries will often use analytical solutions because there is insufficient information to warrant the analysis of numerical models. Here, this study extends prior analytical work on seismic FSI of base-supported tanks in the 1980s to accommodate head-supported tanks with geometries similar to those proposed for new-build, liquid (fluid)-filled, advanced reactors. Analytical solutions are derived for flexible, head-supported, cylindrical tanks subjected to small-amplitude, unidirectional, horizontal seismic motion, and address frequencies of lateral motions of the tank, hydrodynamic pressures, frequencies of waves, wave heights, and reactions at the head support. Solutions are presented for a range of tank dimensions and could be used for a) preliminary design of head-supported cylindrical tanks and their supporting structures with idealized seismic input, and b) verification of numerical models to be used for final seismic analysis, design, and probabilistic risk assessment.

42 ENGINEERING↗