Engineering Papers⌕ Search

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.

At least 469 records · Page 26
All sourcesNASA NTRSDOE OSTI
NASA NTRS1983

Computer modeling of batteries from non-linear circuit elements

A simple non-linear circuit model for battery behavior is given. It is based on time-dependent features of the well-known PIN change storage diode, whose behavior is described by equations similar to those associated with electrochemical cells. The circuit simulation computer program ADVICE was used to predict non-linear response from a topological description of the battery analog built from advice components. By a reasonable choice of one set of parameters, the circuit accurately simulates a wide spectrum of measured non-linear battery responses to within a few millivolts.

Waaben, S.↗
NASA NTRS1990

Assessment of computational models for multilayered composite shells

A review is made of the different approaches used for modeling multilayered composite shells. Discussion focuses on different approaches for developing two-dimensional shear deformation theories; classification of two-dimensional theories based on introducing plausible displacement, strain and/or stress assumptions in the thickness direction; first-order shear deformation theories based on linear displacement assumptions in the thickness coordinate; and efficient computational strategies for anisotropic composite shells. Extensive numerical results are presented showing the effects of variation in the lamination and geometric parameters of simply supported composite cylinders on the accuracy of the static and vibrational responses predicted by eight different modeling approaches (based on two-dimensional shear deformation theories).

Noor, Ahmed K.↗
DOE OSTI2019

Computational modeling of curcumin-based fluorescent probe molecules

In recent years, a series of curcumin analogs have been designed as fluorescent probes for detecting and imaging Aβ peptide aggregates and reactive oxygen species (ROS) in Alzheimer’s disease (AD) brains. In order to gain a better understanding of the photophysical properties of these probe molecules, a systematical computational investigation was performed using the time-dependent density functional theory (TDDFT) calculations. Computed absorption and emission wavelengths well reproduced the spectral shifts among the curcumin analogs. In particular, for a recently proposed pair of probe molecules, CRANAD-5 and CRANAD-61, for sensing ROS in preclinical studies of AD brains, their emission wavelength difference was found to arise from a delocalization of the lowest unoccupied molecular orbital of CRANAD-61 from the curcuminoid backbone to the oxalate moiety. Altogether, this study reaffirms the value of employing TDDFT calculations to assist the design of new curcumin-based fluorescence probes for AD research.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
NASA NTRS1978

JPL Energy Consumption Program (ECP) documentation: A computer model simulating heating, cooling and energy loads in buildings

The engineering manual provides a complete companion documentation about the structure of the main program and subroutines, the preparation of input data, the interpretation of output results, access and use of the program, and the detailed description of all the analytic, logical expressions and flow charts used in computations and program structure. A numerical example is provided and solved completely to show the sequence of computations followed. The program is carefully structured to reduce both user's time and costs without sacrificing accuracy. The user would expect a cost of CPU time of approximately $5.00 per building zone excluding printing costs. The accuracy, on the other hand, measured by deviation of simulated consumption from watt-hour meter readings, was found by many simulation tests not to exceed + or - 10 percent margin.

Lansing, F. L.↗
NASA NTRS1980

NASCAP modelling computations on large optics spacecraft in geosynchronous substorm environments

Satellites in geosynchronous orbits have been found to be charged to significant negative voltages during encounters with geomagnetic substorms. When satellite surfaces are charged, there is a probability of enhanced contamination from charged particles attracted back to the satellite by electrostatic forces. This could be particularly disturbing to large satellites using sensitive optical systems. In this study the NASA Charging Analyzer Program (NASCAP) is used to evaluate qualitatively the possibility of such enhanced contamination on a conceptual version of a large satellite. The evaluation is made by computing surface voltages on the satellite due to encounters with substorm environments and then computing charged-particle trajectories in the electric fields around the satellite. Particular attention is paid to the possibility of contaminants reaching a mirror surface inside a dielectric tube because this mirror represents a shielded optical surface in the satellite model used. Deposition of low energy charged particles from other parts of the spacecraft onto the mirror was found to be possible in the assumed moderate substorm environment condition. In the assumed severe substorm environment condition, however, voltage build up on the inside and edges of the dielectric tube in which the mirror is located prevents contaminants from reaching the mirror surface.

Stevens, N. J.↗
DOE OSTI2020

Experimentally Validated Computational Modeling of Creep and Creep-Cracking for Nuclear Concrete Structures

In a Nuclear Power Plant, one of the most important components is the concrete nuclear reactor cavity, which serves both a structural and protective function as the biological radiation shield. Given that creep has been identified as a major knowledge gap in the assessment of nuclear structures (NUREG/CR-7153), this work helps to further the understanding of creep behavior of massive concrete containment structures for decades to enable safe and long-term operation of these facilities. This project has developed a robust, experimentally validated model to predict creep in nuclear concrete structures for up 60 years using short-term creep data thereby enabling a longer service life of critical facilities and early detection of structural failure. The work presented in this report is a pairing of computational and experimental methods. For the first time, the time temperature superposition (TTS) principle was successfully used to generate a uniaxial creep compliance master curve to predict mortar creep response for up to 22,500 days (nearly 60 years) at a reference temperature of 20°C. These data were used as input into finite element analysis (FEA) codes that use highly realistic random, 3D concrete microstructures from reconstructed coarse limestone aggregates. Finite element analysis performed provides the ability to quickly upscale mortar viscoelastic behavior to long-term concrete creep/relaxation data. A master creep compliance curve, constructed from the TTS principle, spanning 27 years, was used to validate two and a half decades of simulated concrete creep. Concurrently, three different simulated wall specimens were designed to mimic the behavior of post-tensioned concrete nuclear containment facility vessel walls over time as a result of concrete creep. The specimens were designed with different thicknesses, transverse and longitudinal reinforcement ratios, and level of post-tensioning stress. Each specimen contained various instrumentation to measure internal concrete temperature, concrete strain, and post-tensioning strain hourly for over 3 years. The concrete creep model developed in this project, based on the FEA concrete simulations, was applied to simulate the structural-scale experiments of prestressed concrete walls conducted in this project using the Grizzly code. These models can represent the effects of reinforcing and prestressing. Although there are some discrepancies with the experimental data, the model can predict the overall trends of the creep response in these experiments. One of these experimental models was also applied to an extended time to demonstrate how the findings from this study can be used to predict the behavior of actual structures of interest that have been in service for extended periods of time.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2015

Computational modeling of drug-resistant bacteria. Final report

Initial proposal summary: The evolution of antibiotic-resistant mutants among bacteria (superbugs) is a persistent and growing threat to public health. In many ways, we are engaged in a war with these microorganisms, where the corresponding arms race involves chemical weapons and biological targets. Just as advances in microelectronics, imaging technology and feature recognition software have turned conventional munitions into smart bombs, the long-term objectives of this proposal are to develop highly effective antibiotics using next-generation biomolecular modeling capabilities in tandem with novel subatomic feature detection software. Using model compounds and targets, our design methodology will be validated with correspondingly ultra-high resolution structure-determination methods at premier DOE facilities (single-crystal X-ray diffraction at Argonne National Laboratory, and neutron diffraction at Oak Ridge National Laboratory). The objectives and accomplishments are summarized.

60 APPLIED LIFE SCIENCES↗
NASA NTRS1975

Use of computer modeling to investigate a dynamic interaction problem in the Skylab TACS quad-valve package

A valve opening-response problem encountered during development of a control valve for the Skylab thruster attitude control system (TACS) is described. The problem involved effects of dynamic interaction among valves in the quad-redundant valve package. Also described is a detailed computer simulation of the quad-valve package which was helpful in resolving the problem.

Hesser, R. J.↗
NASA NTRS1987

Development of computer models for correlating data of film cooling of nose cone under hypersonic flow

A 12.5 degree half cone with tangential slot injection at Mach 6.95 was studied to determine the heating rates to the surface of the body near and far downstream of the slot. The cone had a zero degree angle of attack. The heating rates were obtained using a computer program that was developed at NASA-Langley Research Center. The concentration of nitrogen from the slot into the boundary layer was also determined. The ratio of slot to freestream was varied to determine its effect on heating. The numerical heating rates were compared to other correlations obtained from experimental studies as well as theoretical laminar and turbulent results.

Sharpe, L., Jr.↗
DOE OSTI2022

Heuristic Computational Model for Predicting Lignin Solubility in Tailored Organic Solvents

Lignin is a random heteropolymer that has been extensively studied as a renewable source of aromatic precursors for high-value chemicals, biofuels, and bioplastics. A key challenge in lignin valorization is the structural and compositional heterogeneity of lignin feedstocks. Solvent-based approaches are commonly used to fractionate lignin to reduce this heterogeneity, but solvent selection can be challenging due to variability in lignin composition. In this work, we developed computational methods to predict good and poor organic solvents as a function of lignin composition. We analyzed 28 different linear pentamer structures, 18 from known libraries and 10 hypothetical polymers, and calculated their activity coefficients in 50 different organic solvents by using the conductor-like screening model for realistic solvents. We used these data to train a regression model that enabled the extensive investigation of the impact of solvent and monolignol compositions on predicted lignin solubility. The exhaustive exploration of solubility trends using model predictions revealed sets of solvents, identified using Kamlet–Taft parameters, that are predicted to promote lignin dissolution regardless of lignin composition. We further identified solvents expected to selectively isolate lignin fractions enriched in certain subunits. Furthermore, these results establish heuristic guidelines for solvent selection that can be used to tailor fractionation processes for lignin feedstocks of distinct composition or to design new processes that isolate fractions with higher proportions of selected subunits.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1992

Reactant conversion in homogeneous turbulence: Mathematical modeling, computational validations and practical applications

Closed form analytical expressions are obtained for predicting the limited rate of reactant conversion in a binary reaction of the type F + rO yields (1 + r) Product in unpremixed homogeneous turbulence. These relations are obtained by means of a single point Probability Density Function (PDF) method based on the Amplitude Mapping Closure. It is demonstrated that with this model, the maximum rate of the reactants' decay can be conveniently expressed in terms of definite integrals of the Parabolic Cylinder Functions. For the cases with complete initial segregation, it is shown that the results agree very closely with those predicted by employing a Beta density of the first kind for an appropriately defined Shvab-Zeldovich scalar variable. With this assumption, the final results can also be expressed in terms of closed form analytical expressions which are based on the Incomplete Beta Functions. With both models, the dependence of the results on the stoichiometric coefficient and the equivalence ratio can be expressed in an explicit manner. For a stoichiometric mixture, the analytical results simplify significantly. In the mapping closure, these results are expressed in terms of simple trigonometric functions. For the Beta density model, they are in the form of Gamma Functions. In all the cases considered, the results are shown to agree well with data generated by Direct Numerical Simulations (DNS). Due to the simplicity of these expressions and because of nice mathematical features of the Parabolic Cylinder and the Incomplete Beta Functions, these models are recommended for estimating the limiting rate of reactant conversion in homogeneous reacting flows. These results also provide useful insights in assessing the extent of validity of turbulence closures in the modeling of unpremixed reacting flows. Some discussions are provided on the extension of the model for treating more complicated reacting systems including realistic kinetics schemes and multi-scalar mixing with finite rate chemical reactions in more complex configurations.

Madnia, C. K.↗
DOE OSTI2013

Comparison of kinetic and extended magnetohydrodynamics computational models for the linear ion temperature gradient instability in slab geometry

We perform linear stability studies of the ion temperature gradient (ITG) instability in unsheared slab geometry using kinetic and extended magnetohydrodynamics (MHD) models, in the regime k{sub ∥}/k{sub ⊥}≪1. The ITG is a parallel (to B) sound wave that may be destabilized by finite ion Larmor radius (FLR) effects in the presence of a gradient in the equilibrium ion temperature. The ITG is stable in both ideal and resistive MHD; for a given temperature scale length L{sub Ti0}, instability requires that either k{sub ⊥}ρ{sub i} or ρ{sub i}/L{sub Ti0} be sufficiently large. Kinetic models capture FLR effects to all orders in either parameter. In the extended MHD model, these effects are captured only to lowest order by means of the Braginskii ion gyro-viscous stress tensor and the ion diamagnetic heat flux. We present the linear electrostatic dispersion relations for the ITG for both kinetic Vlasov and extended MHD (two-fluid) models in the local approximation. In the low frequency fluid regime, these reduce to the same cubic equation for the complex eigenvalue ω=ω{sub r}+iγ. An explicit solution is derived for the growth rate and real frequency in this regime. These are found to depend on a single non-dimensional parameter. We also compute the eigenvalues and the eigenfunctions with the extended MHD code NIMROD, and a hybrid kinetic δf code that assumes six-dimensional Vlasov ions and isothermal fluid electrons, as functions of k{sub ⊥}ρ{sub i} and ρ{sub i}/L{sub Ti0} using a spatially dependent equilibrium. These solutions are compared with each other, and with the predictions of the local kinetic and fluid dispersion relations. Kinetic and fluid calculations agree well at and near the marginal stability point, but diverge as k{sub ⊥}ρ{sub i} or ρ{sub i}/L{sub Ti0} increases. There is good qualitative agreement between the models for the shape of the unstable global eigenfunction for L{sub Ti0}/ρ{sub i}=30 and 20. The results quantify how far fluid calculations can be extended accurately into the kinetic regime. We conclude that for the linear ITG problem in slab geometry with unsheared magnetic field when k{sub ∥}/k{sub ⊥}≪1, the extended MHD model may be a reliable physical model for this problem when ρ{sub i}/L{sub Ti0}<10{sup −2} and k{sub ⊥}ρ{sub i}<0.2.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1997

Computational Modeling in Plasma Processing for 300 mm Wafers

Migration toward 300 mm wafer size has been initiated recently due to process economics and to meet future demands for integrated circuits. A major issue facing the semiconductor community at this juncture is development of suitable processing equipment, for example, plasma processing reactors that can accomodate 300 mm wafers. In this Invited Talk, scaling of reactors will be discussed with the aid of computational fluid dynamics results. We have undertaken reactor simulations using CFD with reactor geometry, pressure, and precursor flow rates as parameters in a systematic investigation. These simulations provide guidelines for scaling up in reactor design.

Meyyappan, Meyya↗
DOE OSTI2025

Computational Modeling of Graphite Degradation due to Molten Salt Infiltration and Wear

Molten-salt reactors (MSRs) represent a promising next-generation reactor design, with graphite serving as a moderator and/or reflector in several designs. However, due to limited experimental data and operational experience, a technical understanding of the structural integrity of graphite in molten salt environments remains incomplete. This report presents a modeling-based evaluation of graphite degradation in MSR environments, focusing on the effects of salt infiltration in fuel salt-based designs and surface wear in pebble bed reactor designs. The objective of this study is to enhance understanding of the structural integrity challenges posed by these degradation mechanisms and to provide a framework for assessing graphite behavior in MSRs. The first part of the report investigates the phenomenon of molten salt infiltration into graphite. This infiltration occurs when molten salt permeates the interconnected pore structure of the graphite moderator, driven by factors such as pressure differentials and the physical properties of both the salt and graphite. The infiltration process is influenced by characteristics of the pore structure, viscosity of the molten salt, and the interfacial energies between the graphite, salt, and the atmosphere within the graphite pore. Utilizing a coupled multiphysics modeling approach with Grizzly software, the study evaluates the stress induced by internal heat sources due to infiltration, which can lead to structural concerns. This evaluation is crucial for understanding how infiltration affects the mechanical integrity of graphite components in MSRs. The study considers the Molten-Salt Reactor Experiment (MSRE) graphite stringer geometry due to the availability of relevant data. Through detailed finite element analysis, the study examines stress distributions at varying infiltration percentages, revealing that stress levels increase with higher amounts of infiltration. Rare-event simulations, using the parallel subset simulation (PSS) framework, further quantify the failure probabilities under input uncertainties, with a user-specified failure metric. The PSS framework also identifies critical input parameters that significantly affect the stress values, including infiltration amount, thermal conductivity, and power density. Additionally, considering realistic reactor scenarios, the analysis was performed to account for the combined effects of radiation and infiltration, and modeling strategies on how to analyze new reactor designs or new graphite grades are discussed. The second part of the report focuses on wear mechanisms in pebble bed-based MSRs. As graphite fuel pebbles interact with the graphite reflector block, wear can result in material loss and the formation of surface defects, which may act as stress concentrators. A similar multiphysics modeling framework is employed to assess the impact of wear on the structural integrity of graphite components. This study considers a generic fluoride-cooled high-temperature reactor (gFHR) design due to the availability of comprehensive data. Worst-case scenario dimensions of the reflector blocks were analyzed under thermal and radiation conditions. Subsequently, wear in the form of idealized pits and grooves is modeled on the inner surface of the graphite block, with the maximum stress from previous simulations. The simulations show that groove-type defects are more detrimental than pits, leading to higher stress concentrations. Considering worst-case simulation scenarios and experimental wear rates, it was determined that the formation of a surface defect critical enough to affect the stress may not be possible in a gFHR design. Overall, the findings of this research contribute to the development of robust modeling tools for predicting graphite behavior under various operational conditions in MSRs.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS1996

Computational Modeling of Multi-Phase/Multi-Species Flows with Applications to Liquid Rocket Engines

Accurate prediction of all physical phenomena in a combustion chamber is essential for better understanding of the system performance. Atomization, evaporation, combustion, chemical kinetics, and turbulence are those processes of great importance that need to be well understood. Processes involving the liquid phase in a combustion chamber will be further complicated under supercritical conditions. More advanced and accurate numerical techniques are required to extend our understanding of the above phenomena. A computer program for multi-species/multi-phase flow was developed for NASA/MSFC in 1992. This code, called Liquid Thrust Chamber Performance (LTCP) program takes an Eulerian- Eulerian approach and is based on the Total Variation Diminishing (TVD) technique with Lax-Friedrichs upwind method. Under the NASA/ASEE SFFP the LTCP code was used to predict the performance characteristics of several engines that were of particular interest to NASA. This code was also successful in a combustion detonation study. Converting the program to the PC platform was accomplished which extends usability and makes it available to a wider range of users. The Eulerian formulation of the liquid phase provides a suitable model that can be extended to include combustion modeling under supercritical conditions. The results have been compared against the ones of other codes and available measured data. The algorithm proved to be robust and efficient for problems with stiff source terms.

Navaz, Homayun K.↗
NASA NTRS2006

Computer Modeling of Thermal Convection in Melts to Explain Glass Formation in Low Gravity and on Earth

Experiments conducted up to this time on glass forming melts in the low gravity environment of space show that glasses prepared in low-g are more chemically homogeneous and more resistant to crystallization than the comparable glasses prepared at 1-g on Earth. This result is somewhat surprising and opposite to the accepted concept on glass formation for a melt. A hypothesis based on "shear thinning" of a melt, a decrease in viscosity with increasing shear stress, is proposed as an explanation for the observed low-gravity results. This paper describes detailed simulation procedures to test the role of thermal convection in introducing shear stress in glass forming melts, using a lithium disilcate melt as a model. The simulation system in its idealized version consists of a cylinder that is heated at one end and cooled at the other with gravity acting in a transverse direction to the thermal gradient. The side wall of the cylinder is assumed to be insulating. The governing equations of motion and energy are solved using variable properties for viscosity (Arrehenius and non-Arrehenius behaviors) and density (constant and temperature dependent). Other parametric variables in the calculations include gravity level and gravity vector orientation. The shear stress in the system are then computed as a function of gravity from the calculated values of maximum melt velocity, and its effect on melt viscosity (shear thinning) is predicted. Also included and discussed are the modeling efforts related to other potential convective processes in glass forming melts and their possible effects on melt viscosity.

Ray, Chandra S.↗
DOE OSTI2014

Towards a Computational Model of a Methane Producing Archaeum

Progress towards a complete model of the methanogenic archaeum Methanosarcina acetivorans is reported. We characterized size distribution of the cells using differential interference contrast microscopy, finding them to be ellipsoidal with mean length and width of 2.9 μ m and 2.3 μ m, respectively, when grown on methanol and 30% smaller when grown on acetate. We used the single molecule pull down (SiMPull) technique to measure average copy number of the Mcr complex and ribosomes. A kinetic model for the methanogenesis pathways based on biochemical studies and recent metabolic reconstructions for several related methanogens is presented. In this model, 26 reactions in the methanogenesis pathways are coupled to a cell mass production reaction that updates enzyme concentrations. RNA expression data (RNA-seq) measured for cell cultures grown on acetate and methanol is used to estimate relative protein production per mole of ATP consumed. The model captures the experimentally observed methane production rates for cells growing on methanol and is most sensitive to the number of methyl-coenzyme-M reductase (Mcr) and methyl-tetrahydromethanopterin:coenzyme-M methyltransferase (Mtr) proteins. A draft transcriptional regulation network based on known interactions is proposed which we intend to integrate with the kinetic model to allow dynamic regulation.

59 BASIC BIOLOGICAL SCIENCES↗