SEARCH · Engineering Papers
Results for “computer modeling”
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Computational Modeling of the Nature and Role of Ga Species for Light Alkane Dehydrogenation Catalyzed by Ga/H-MFI
Ga-exchanged H-MFI zeolites are highly active for the dehydrogenation of light alkanes; however, both the nature of the active gallium species and the associated dehydrogenation mechanism have been difficult to establish. In this study, we examine the activity of Ga species in Ga/H-MFI by calculating the free energy landscapes on which all reactions occur. To this end, we use a hybrid quantum mechanics/molecular mechanics model for all electronic structure calculations. Quantum chemical calculations were carried out with a range-corrected functional and a good representation of dispersive interactions. The molecular mechanics part of our approach captures the long-range effects of Coulombic and dispersive interactions due to atoms in the extended framework. The rate-determining TS (RDTS) is identified by analysis of the free energy landscape for each mechanism, using the energetic span model. Our analysis reveals that, for reduced Ga/H-MFI, univalent and divalent gallium hydrides, [GaH 2 ] + and [GaH] 2+ , respectively, are more active for ethane dehydrogenation in comparison to H + sites and Ga + sites. [GaH] 2+ sites consistently emerge as the most active sites for light alkane dehydrogenation, providing significant enthalpic stabilization to C-H cleavage TSs via alkyl and carbenium dehydrogenation routes. In contrast, carbenium-like C-H cleavage TSs occurring on Brønsted acid sites are enthalpically less favorable due to their limited electronic interactions with framework O atoms and H + sites. Activation enthalpy barriers for dehydrogenation determined using the energetic span analysis are in good agreement with those measured experimentally. Accounting for entropy of activation reveals that constrained TSs become less favorable in free energy with increasing chain length. Here, we also find that an increase in enthalpic favorability of the alkyl mechanism is observed with increasing chain length for the TS responsible for the second C-H cleavage step leading to alkene formation.
Establishing Cost-Effective Computational Models for the Prediction of Lanthanoid Binding in [Ln(NO 3 )] 2+ (with Ln = La to Lu)
Not Available
Dynamics of episodic supershear in the 2023 M7.8 Kahramanmaraş/Pazarcik earthquake, revealed by near-field records and computational modeling
Abstract The 2023 M7.8 Kahramanmaraş/Pazarcik earthquake was larger and more destructive than what had been expected. Here we analyzed nearfield seismic records and developed a dynamic rupture model that reconciles different currently conflicting inversion results and reveals spatially non-uniform propagation speeds in this earthquake, with predominantly supershear speeds observed along the Narli fault and at the southwest (SW) end of the East Anatolian Fault (EAF). The model highlights the critical role of geometric complexity and heterogeneous frictional conditions in facilitating continued propagation and influencing rupture speed. We also constrained the conditions that allowed for the rupture to jump from the Narli fault to EAF and to generate the delayed backpropagating rupture towards the SW. Our findings have important implications for understanding earthquake hazards and guiding future response efforts and demonstrate the value of physics based dynamic modeling fused with near-field data in enhancing our understanding of earthquake mechanisms and improving risk assessment.
Computational modeling predicts the stability of both Pd + and Pd2 + ion-exchanged into H-CHA
Passive NO x adsorbers (PNA) using Pd/zeolites have emerged as a promising solution for the reduction of cold-start emissions from vehicle exhaust. However, the nature of the active sites and the mechanisms underlying NO x adsorption in Pd/zeolites remain a subject of ongoing investigation. In this study, we employ quantum chemical simulations to investigate the structure of Pd species in cation-exchange sites at isolated Al and Al pairs in the 6-ring and 8-ring of the CHA framework before the introduction of NO x . Our calculations show that the speciation of Pd in these exchange sites strongly depends on the precise Al arrangement within the framework, as well as the operating conditions. Ionically dispersed Pd is found to be the most favorable species over a wide range of oxidizing and reducing conditions. Small oligomers of PdO and metallic Pd do not appear to be competitive at either isolated Al or Al pairs. Notably, our calculations show that ion exchange sites other than next–next–nearest neighbor Al pairs in the 6-ring will be preferentially occupied by Pd + instead of Pd 2+ . The stability of Pd + in the zeolite environment is an interesting contrast with its rareness in molecular Pd compounds. Nonetheless, a detailed analysis of the electronic structure shows that predicted Pd oxidation states are consistent with chemical intuition for all complexes investigated in this study. We also discuss the potential ambiguity in Pd characterization provided by typical experimental techniques such as XANES, EXAFS and UV-Vis, and highlight the need for additional EPR spectroscopy studies to further elucidate the initial Pd speciation in zeolites for PNA applications.
Computational models for the berry phase in semiconductor quantum dots
By developing a new model and its finite element implementation, we analyze the Berry phase low-dimensional semiconductor nanostructures, focusing on quantum dots (QDs). In particular, we solve the Schrödinger equation and investigate the evolution of the spin dynamics during the adiabatic transport of the QDs in the 2D plane along circular trajectory. Based on this study, we reveal that the Berry phase is highly sensitive to the Rashba and Dresselhaus spin-orbit lengths.
Computational modeling of electrophotonics nanomaterials: Tunneling in double quantum dots
Single electron localization and tunneling in double quantum dots (DQD) and rings (DQR) and in particular the localized-delocalized states and their spectral distributions are considered in dependence on the geometry of the DQDs (DQRs). The effect of violation of symmetry of DQDs geometry on the tunneling is studied in details. The cases of regular and chaotic geometries are considered. It will be shown that a small violation of symmetry drastically affects localization of electron and that anti-crossing of the levels is the mechanism of tunneling between the localized and delocalized states in DQRs.
Mathematical and computational modeling of the diffraction problems by discrete singularities method
The main objective of this study is reduced the boundary-value problems of scattering and diffraction waves on plane-parallel structures to the singular or hypersingular integral equations. For these cases we use a method of the parametric representations of the integral and pseudo-differential operators. Numerical results of the model scattering problems on periodic and boundary gratings and also on the gratings above a flat screen reflector are presented in this paper.
Dehydration of multilamellar fatty acid membranes: Towards a computational model of the stratum corneum
The level of hydration controls the cohesion between apposed lamellae of saturated free fatty acids found in the lipid matrix of stratum corneum, the outermost layer of mammalian skin. This multilamellar lipid matrix is highly impermeable to water and ions, so that the local hydration shell of its fatty acids may not always be in equilibrium with the acidity and relative humidity, which significantly change over a course of days during skin growth. The homeostasis of the stratum corneum at each moment of its growth likely requires a balance between two factors, which affect in opposite ways the diffusion of hydrophilic species through the stratum corneum: (i) an increase in water order as the lipid lamellae come in closer contact, and (ii) a decrease in water order as the fraction of charged fatty acids is lowered by pH. Herein molecular dynamics simulations are employed to estimate the impact of both effects on water molecules confined between lamellae of fatty acids. Under conditions where membrane undulations are energetically favorable, the charged fatty acids are able to sequester cations around points of contact between lamellae that are fully dehydrated, while essentially maintaining a multilamellar structure for the entire system. This observation suggests that the undulations of the fatty acid lamellae control the diffusion of hydrophilic species through the water phase by altering the positional and rotational order of water molecules in the embedded/occluded “droplets.”
Precise control of polymer coated nanopores by nanoparticle additives: Insights from computational modeling
Polymer-nanoparticle composites are a promising new class of materials for creation of controllable nano-patterned surfaces and nanopores. We use coarse-grained molecular dynamics simulations augmented with analytical theory to study the structural transitions of surface grafted polymer layers (brushes) induced by infiltration of nanoparticles that are attracted to the polymers in the layer. We systematically compare two different polymer brush geometries: one where the polymer chains are grafted to a planar surface and the other where the chains are grafted to the inside of a cylindrical nanochannel. We perform a comprehensive study of the effects of the material parameters such as the polymer chain length, chain grafting density, nanoparticle size, strength of attraction between nanoparticles and polymer monomers, and, in the case of the cylindrically grafted brush, the radius of the cylinder. We find a very general behavioral motif for all geometries and parameter values: the height of the polymer brush is non-monotonic in the nanoparticle concentration in solution. As the nanoparticle concentration increases, the brush height first decreases and after passing through a minimum value begins to increase, resulting in the swelling of the nanoparticle infused brush. These morphological features may be useful for devising tunable “smart” nano-devices whose effective dimensions can be reversibly and precisely adjusted by changing the nanoparticle concentration in solution. The results of approximate Self-Consistent Field Theory (SCFT) calculations, applicable in the regime of strong brush stretching, are compared to the simulation results. The SCFT calculations are found to be qualitatively, even semi-quantitatively, accurate when applied within their intended regime of validity, and provide a useful and efficient tool for modeling such materials.
Uncertainty quantification for computational modelling of laser powder bed fusion
Explore the source record for details and available documents.
Extending and Applying Spartan to Perform Temporal Sensitivity Analyses for Predicting Changes in Influential Biological Pathways in Computational Models
Not Available
Mechanism-Based Energy Regularization in Computational Modeling of Quasibrittle Fracture
Abstract Quasibrittle materials are featured by a strain-softening constitutive behavior under many loading scenarios, which could eventually lead to localization instability. It has long been known that strain localization would result in spurious mesh sensitivity in finite element (FE) simulations. Previous studies have shown that, for the case of fully localized damage, the mesh sensitivity can be mitigated through energy regularization of the material constitutive law. However, depending on the loading configuration and structural geometry, quasibrittle structures could exhibit a complex damage process, which involves both localized and diffused damage patterns at different stages of loading. This study presents a generalized energy regularization method that considers the spatial and temporal evolution of damage pattern. The method introduces a localization parameter, which describes the local damage pattern. The localization parameter governs the energy regularization of the constitutive model, which captures the transition from diffused to localized damage during the failure process. The method is cast into an isotropic damage model, and is further extended to rate-dependent behavior. The energy regularization scheme is directly incorporated into the kinetics of damage growth. The model is applied to simulate static and dynamic failures of ceramic specimens. It is shown that the present model is able to effectively mitigate the spurious mesh sensitivity in FE simulations of both types of failure. The present analysis demonstrates the essential role of mechanism-based energy regularization of constitutive relation in FE simulations of quasibrittle fracture.
Long-term safety assessment of trench-type surface repository at Chernobyl, Ukraine - computer model and comparison with results from simplified models
The Radioactive Waste Disposal Facility (RWDF) Buryakovka was constructed in 1986 as part of the intervention measures after the accident at Chernobyl NPP (ChNPP). Today, the surface repository for solid low and intermediate level waste (LILW) is still being operated but its maximum capacity is nearly reached. Long-existing plans for increasing the capacity of the facility shall be implemented in the framework of the European Commission INSC Programme (Instrument for Nuclear Safety Co-operation). Within the first phase of this project, DBE Technology GmbH prepared a safety analysis report of the facility in its current state (SAR) and a preliminary safety analysis report (PSAR) for a future extended facility based on the planned enlargement. In addition to a detailed mathematical model, also simplified models have been developed to verify results of the former one and enhance confidence in the results. Comparison of the results show that - depending on the boundary conditions - simplifications like modeling the multi trench repository as one generic trench might have very limited influence on the overall results compared to the general uncertainties associated with respective long-term calculations. In addition to their value in regard to verification of more complex models which is important to increase confidence in the overall results, such simplified models can also offer the possibility to carry out time consuming calculations like probabilistic calculations or detailed sensitivity analysis in an economic manner. (authors)
Computational Modeling of Terry Turbine Airflow Testing to Support the Expansion of Operating Band in Beyond Design Basis Conditions.
Abstract not provided.
A Computational Modeling Tool for Welding Repair of Irradiated Materials
Explore the source record for details and available documents.
Computational Modeling of Multi-Pass Rolling Parameters Effect on Resulting Fuel Foil Shape
A focus of the U.S. Department of Energy is to improve production yield and reduce the cost of Low Enriched Uranium (LEU)-molybdenum alloy (U-10Mo) monolithic fuel plates that will be replacing High Enriched Uranium (HEU) oxide dispersion fuels used currently in the United States High Performance Research Reactors (USHPRR). One area of improvement is lowering the high transverse waviness and longitudinal waviness currently present within rolled foils prior to cladding to produce fuel plates. Traditional rolling manufacturing techniques for other metal foils use winders to pull and straighten the foil as it is rolled back and forth to the final thickness. This approach cannot be used to roll thin U-10Mo foils (0.008-0.025” thick) because only small castings can be rolled due to nuclear criticality safety concerns. As a result, the fuel foils are too short (1 m in length) to use traditional winders. Therefore, it is crucial to identify other rolling parameters (i.e., roller friction, axial tension load, roller diameter, and roll pass reduction percent) that might reduce transverse waviness and longitudinal waviness in the rolled fuel foil and develop a high-yield, low-cost multi-pass rolling manufacturing process. This report documents a systematic finite element modeling study to investigate the effects of numerous rolling parameters to reduce resulting transverse waviness and longitudinal waviness in the fuel foil during multi-pass rolling of U-10Mo foils. The rolling of a U-10Mo plate with initial dimensions of 1”x1”x 0.048” is modeled using Abaqus CAE. This rolling is modeled to undergo eight 20% reduction roll passes to a final fuel foil thickness of 0.01”. The elastic-plastic constitutive model of the U-10Mo alloy was input to the fuel foil rolling model. The rollers were modeled as rigid bodies. A comparison of rolling friction coefficients of 0.3 and 0.7 over a wide range of applied axial tension loads were investigated in order to evaluate the effect of using a lubricant during rolling. The effect of roller diameter on the resulting transverse waviness and longitudinal waviness of the fuel foil over a wide range of axial tension loads were also investigated by modeling rollers 7/8” and 3.75” in diameter. The results of this systematic finite element method study will aid manufacturers in producing low transverse waviness and reduced longitudinal waviness in U-10Mo fuel foils.