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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
NASA NTRS1998

Computational Modeling Program

An Integrated Product Team (IPT) has been formed at NASA Ames Research Center which has set objectives to investigate devices and processes suitable for meeting NASA requirements on ultrahigh performance computers, fast and low power devices, and high temperature wide bandgap materials. These devices may ultimately be sub-100nm feature-size. Processes and equipment must meet the stringent demands posed by the fabrication of such small devices. Until now, the reactors for Chemical Vapor Deposition (CVD) and plasma processes have been designed by trial and error procedures. Further, once the reactor is in place, optimum processing parameters are found through expensive and time-consuming experimentation. If reliable models are available that describe processes and the operation of the reactors, that chore would be reduced to a routine task while being a cost-effective option. The goal is to develop such a design tool, validate that tool using available data from current generation processes and reactors, and then use that tool to explore avenues for meeting NASA needs for ultrasmall device fabrication. Under the present grant, ARL/Penn State along with other IPT members has been developing models and computer code to meet IPT goals. Some of the accomplishments achieved during the first year of the grant are described in this report

Govindan, T. R.↗
DOE OSTI2017

Interfacial Structure and Electric Field Probed by in Situ Electrochemical Vibrational Stark Effect Spectroscopy and Computational Modeling

Interfacial electric fields play crucial roles in electrochemistry, catalysis, and solar energy conversion. Understanding of the interfacial electric field effects has been hindered by the lack of a direct spectroscopic method to probe of the interfacial field at the molecular level. Here, we report the characterization of the field and interfacial structure at Au/diisocyanide/aqueous electrolyte interfaces, using a combination of in situ electrochemical vibrational sum frequency generation (SFG) spectroscopy, density functional theory (DFT) calculations, and molecular dynamics (MD) simulations. For 1,4-phenylene diisocyanide (PDI), 4,4'-biphenyl diisocyanide (BPDI), and 4,4"-p-terphenyl diisocyanide (TPDI), our results reveal that the frequency of the gold-bound NC stretch mode of the diisocyanide self-assembled monolayer (SAM) increases linearly with the applied potential, suggesting that SFG can be an in situ probe of the strength of the electric field at electrode/electrolyte interfaces. Using DFT-computed Stark tuning rates of model complexes, the electric field strength at the metal/SAM/electrolyte interfaces is estimated to be 108–109 V/m. The linear dependence of the vibrational frequency (and field) with applied potential is consistent with an electrochemical double-layer structure that consists of a Helmholtz layer in contact with a diffused layer. The Helmholtz layer thickness is approximately the same as the molecular length for PDI, suggesting a well-ordered SAM with negligible electrolyte penetration. For BPDI and TPDI, we found that the Helmholtz layer is thinner than the monolayer of molecular adsorbates, indicating that the electrolyte percolates into the SAM, as shown by molecular dynamics simulations of the Au/PDI/electrolyte interface. The reported analysis demonstrates that a combination of in situ SFG probes and computational modeling provides a powerful approach to elucidate the structure of electrochemical interfaces at the detailed molecular level.

Ge, Aimin↗
NASA NTRS2015

Delineating the Impact of Weightlessness on Human Physiology Using Computational Models

Microgravity environment has profound effects on several important human physiological systems. The impact of weightlessness is usually indirect as mediated by changes in the biological fluid flow and transport and alterations in the deformation and stress fields of the compliant tissues. In this context, Fluid-Structural and Fluid-Solid Interaction models provide a valuable tool in delineating the physical origins of the physiological changes so that systematic countermeasures can be devised to reduce their adverse effects. In this presentation, impact of gravity on three human physiological systems will be considered. The first case involves prediction of cardiac shape change and altered stress distributions in weightlessness. The second, presents a fluid-structural-interaction (FSI) analysis and assessment of the vestibular system and explores the reasons behind the unexpected microgravity caloric stimulation test results performed aboard the Skylab. The last case investigates renal stone development in microgravity and the possible impact of re-entry into partial gravity on the development and transport of nucleating, growing, and agglomerating renal calculi in the nephron. Finally, the need for model validation and verification and application of the FSI models to assess the effects of Artificial Gravity (AG) are also briefly discussed.

Computational Models↗
DOE OSTI2014

Computational modeling and analysis of thermoelectric properties of nanoporous silicon

In this paper, thermoelectric properties of nanoporous silicon are modeled and studied by using a computational approach. The computational approach combines a quantum non-equilibrium Green's function (NEGF) coupled with the Poisson equation for electrical transport analysis, a phonon Boltzmann transport equation (BTE) for phonon thermal transport analysis and the Wiedemann-Franz law for calculating the electronic thermal conductivity. By solving the NEGF/Poisson equations self-consistently using a finite difference method, the electrical conductivity σ and Seebeck coefficient S of the material are numerically computed. The BTE is solved by using a finite volume method to obtain the phonon thermal conductivity k{sub p} and the Wiedemann-Franz law is used to obtain the electronic thermal conductivity k{sub e}. The figure of merit of nanoporous silicon is calculated by ZT=S{sup 2}σT/(k{sub p}+k{sub e}). The effects of doping density, porosity, temperature, and nanopore size on thermoelectric properties of nanoporous silicon are investigated. It is confirmed that nanoporous silicon has significantly higher thermoelectric energy conversion efficiency than its nonporous counterpart. Specifically, this study shows that, with a n-type doping density of 10{sup 20} cm{sup –3}, a porosity of 36% and nanopore size of 3 nm × 3 nm, the figure of merit ZT can reach 0.32 at 600 K. The results also show that the degradation of electrical conductivity of nanoporous Si due to the inclusion of nanopores is compensated by the large reduction in the phonon thermal conductivity and increase of absolute value of the Seebeck coefficient, resulting in a significantly improved ZT.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS2002

Computational modeling of cardiovascular response to orthostatic stress

The objective of this study is to develop a model of the cardiovascular system capable of simulating the short-term (< or = 5 min) transient and steady-state hemodynamic responses to head-up tilt and lower body negative pressure. The model consists of a closed-loop lumped-parameter representation of the circulation connected to set-point models of the arterial and cardiopulmonary baroreflexes. Model parameters are largely based on literature values. Model verification was performed by comparing the simulation output under baseline conditions and at different levels of orthostatic stress to sets of population-averaged hemodynamic data reported in the literature. On the basis of experimental evidence, we adjusted some model parameters to simulate experimental data. Orthostatic stress simulations are not statistically different from experimental data (two-sided test of significance with Bonferroni adjustment for multiple comparisons). Transient response characteristics of heart rate to tilt also compare well with reported data. A case study is presented on how the model is intended to be used in the future to investigate the effects of post-spaceflight orthostatic intolerance.

Flight Experiment↗
NASA NTRS1982

Computer modeling of thermoelectric generator performance

Features of the DEGRA 2 computer code for simulating the operations of a spacecraft thermoelectric generator are described. The code models the physical processes occurring during operation. Input variables include the thermoelectric couple geometry and composition, the thermoelectric materials' properties, interfaces and insulation in the thermopile, the heat source characteristics, mission trajectory, and generator electrical requirements. Time steps can be specified and sublimation of the leg and hot shoe is accounted for, as are shorts between legs. Calculations are performed for conduction, Peltier, Thomson, and Joule heating, the cold junction can be adjusted for solar radition, and the legs of the thermoelectric couple are segmented to enhance the approximation accuracy. A trial run covering 18 couple modules yielded data with 0.3% accuracy with regard to test data. The model has been successful with selenide materials, SiGe, and SiN4, with output of all critical operational variables.

Chmielewski, A. B.↗
DOE OSTI2024

Computational modeling of phononic pseudocrystal isolators

Methods for the efficient computational prediction of the performance of phononic pseudocrystals (structured materials capable of blocking extraordinary ranges of frequency) in COMSOL and other comparable finite element method codes are set forth. These methods include boundary conditions that make possible halving the size of the computational domain. Also included is an introduction of elastic energy density methods for assessing the extinction of elastic waves within the patterned region.

Swift, Stephen Hales↗
DOE OSTI2015

A UNIFIED COMPUTATIONAL MODEL FOR SOLAR AND STELLAR FLARES

We present a unified computational framework that can be used to describe impulsive flares on the Sun and on dMe stars. The models assume that the flare impulsive phase is caused by a beam of charged particles that is accelerated in the corona and propagates downward depositing energy and momentum along the way. This rapidly heats the lower stellar atmosphere causing it to explosively expand and dramatically brighten. Our models consist of flux tubes that extend from the sub-photosphere into the corona. We simulate how flare-accelerated charged particles propagate down one-dimensional flux tubes and heat the stellar atmosphere using the Fokker–Planck kinetic theory. Detailed radiative transfer is included so that model predictions can be directly compared with observations. The flux of flare-accelerated particles drives return currents which additionally heat the stellar atmosphere. These effects are also included in our models. We examine the impact of the flare-accelerated particle beams on model solar and dMe stellar atmospheres and perform parameter studies varying the injected particle energy spectra. We find the atmospheric response is strongly dependent on the accelerated particle cutoff energy and spectral index.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2015

A Unified Computational Model for Solar and Stellar Flares

We present a unified computational framework that can be used to describe impulsive flares on the Sun and on dMe stars. The models assume that the flare impulsive phase is caused by a beam of charged particles that is accelerated in the corona and propagates downward depositing energy and momentum along the way. This rapidly heats the lower stellar atmosphere causing it to explosively expand and dramatically brighten. Our models consist of flux tubes that extend from the sub-photosphere into the corona. We simulate how flare-accelerated charged particles propagate down one-dimensional flux tubes and heat the stellar atmosphere using the Fokker-Planck kinetic theory. Detailed radiative transfer is included so that model predictions can be directly compared with observations. The flux of flare-accelerated particles drives return currents which additionally heat the stellar atmosphere. These effects are also included in our models. We examine the impact of the flare-accelerated particle beams on model solar and dMe stellar atmospheres and perform parameter studies varying the injected particle energy spectra. We find the atmospheric response is strongly dependent on the accelerated particle cutoff energy and spectral index.

radiative transfer↗
DOE OSTI2014

Comprehensive computer model for magnetron sputtering. II. Charged particle transport

Discharges for magnetron sputter thin film deposition systems involve complex plasmas that are sensitively dependent on magnetic field configuration and strength, working gas species and pressure, chamber geometry, and discharge power. The authors present a numerical formulation for the general solution of these plasmas as a component of a comprehensive simulation capability for planar magnetron sputtering. This is an extensible, fully three-dimensional model supporting realistic magnetic fields and is self-consistently solvable on a desktop computer. The plasma model features a hybrid approach involving a Monte Carlo treatment of energetic electrons and ions, along with a coupled fluid model for thermalized particles. Validation against a well-known one-dimensional system is presented. Various strategies for improving numerical stability are investigated as is the sensitivity of the solution to various model and process parameters. In particular, the effect of magnetic field, argon gas pressure, and discharge power are studied.

36 MATERIALS SCIENCE↗
DOE OSTI2013

Multi-scale computational model of three-dimensional hemodynamics within a deformable full-body arterial network

In this article, we present a computational multi-scale model of fully three-dimensional and unsteady hemodynamics within the primary large arteries in the human. Computed tomography image data from two different patients were used to reconstruct a nearly complete network of the major arteries from head to foot. A linearized coupled-momentum method for fluid–structure-interaction was used to describe vessel wall deformability and a multi-domain method for outflow boundary condition specification was used to account for the distal circulation. We demonstrated that physiologically realistic results can be obtained from the model by comparing simulated quantities such as regional blood flow, pressure and flow waveforms, and pulse wave velocities to known values in the literature. We also simulated the impact of age-related arterial stiffening on wave propagation phenomena by progressively increasing the stiffness of the central arteries and found that the predicted effects on pressure amplification and pulse wave velocity are in agreement with findings in the clinical literature. This work demonstrates the feasibility of three-dimensional techniques for simulating hemodynamics in a full-body compliant arterial network.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2012

A Multiscale Computational Model of the Response of Swine Epidermis After Acute Irradiation

Radiation exposure from Solar Particle Events can lead to very high skin dose for astronauts on exploration missions outside the protection of the Earth s magnetic field [1]. Assessing the detrimental effects to human skin under such adverse conditions could be predicted by conducting territorial experiments on animal models. In this study we apply a computational approach to simulate the experimental data of the radiation response of swine epidermis, which is closely similar to human epidermis [2]. Incorporating experimentally measured histological and cell kinetic parameters into a multiscale tissue modeling framework, we obtain results of population kinetics and proliferation index comparable to unirradiated and acutely irradiated swine experiments [3]. It is noted the basal cell doubling time is 10 to 16 days in the intact population, but drops to 13.6 hr in the regenerating populations surviving irradiation. This complex 30-fold variation is proposed to be attributed to the shortening of the G1 phase duration. We investigate this radiation induced effect by considering at the sub-cellular level the expression and signaling of TGF-beta, as it is recognized as a key regulatory factor of tissue formation and wound healing [4]. This integrated model will allow us to test the validity of various basic biological rules at the cellular level and sub-cellular mechanisms by qualitatively comparing simulation results with published research, and should lead to a fuller understanding of the pathophysiological effects of ionizing radiation on the skin.

FROM↗
DOE OSTI2025

Computational Modeling of a 3D Printed Recuperator and Subsequent Experimental Loop for Supercritical Carbon Dioxide Cycles

Oak Ridge National Laboratory (ORNL), in collaboration with mechanical-thermal energy storage (mTES) provider EarthEn, a US Department of Energy (DOE) Lab-Embedded Entrepreneurship Program (LEEP) recipient at ORNL’s Innovation Crossroads 2023, is utilizing a state-of-the-art patented 3D printing technique to design an additively manufactured (AM) supercritical CO2 (sCO2) recuperator (REC) for EarthEn’s charge/discharge cycle. The AM REC will be printed at ORNL’s Manufacturing Demonstration Facility using Inconel Alloy 718 and tested on a closed-loop, ∼100 kW scale experimental facility that is under construction. The testing will compare the printed design against a commercial-off-the-shelf Printed Circuit Heat Exchanger (PCHE) REC. The design of the sCO2 facility is guided by a Modelica-based system model which is primarily dependent on the open-source TRANSFORM library developed at ORNL and uses the open-source CoolProp library for thermophysical properties of sCO2 via the External Media library. It is envisioned that an iterative process will be followed between the physical loop and the system model wherein the initial experimental data will be used to tune the model, which in turn will be used to guide future loop operation. Simultaneously, the AM REC is being designed using computer-aided design models, and it is also being analyzed for hydraulic and thermomechanical response using commercial computational fluid dynamics software, Simcenter STAR-CCM+, on highperformance computing resources.1

See, Nate [ORNL] (ORCID:0000000178581202)↗
DOE OSTI2018

Efficient Computational Modeling of Human Ventricular Activation and Its Electrocardiographic Representation: A Sensitivity Study

Patient-specific models of the ventricular myocardium, combined with the computational power to run rapid simulations, are approaching the level where they could be used for personalized cardiovascular medicine. A major remaining challenge is determining model parameters from available patient data, especially for models of the Purkinje-myocardial junctions (PMJs): the sites of initial ventricular electrical activation. There are no non-invasive methods for localizing PMJs in patients, and the relationship between the standard clinical ECG and PMJ model parameters is underexplored. Thus, this study aimed to determine the sensitivity of the QRS complex of the ECG to the anatomical location and regional number of PMJs. The QRS complex was simulated using an image-based human torso and biventricular model, and cardiac electrophysiology was simulated using Cardioid. The PMJs were modeled as discrete current injection stimuli, and the location and number of stimuli were varied within initial activation regions based on published experiments. Results indicate that the QRS complex features were most sensitive to the presence or absence of four “seed” stimuli, and adjusting locations of nearby “regional” stimuli provided finer tuning. Decreasing number of regional stimuli by an order of magnitude resulted in virtually no change in the QRS complex. Thus, a minimal 12-stimuli configuration was identified that resulted in physiological excitation, defined by QRS complex feature metrics and ventricular excitation pattern. In conclusion, overall, the sensitivity results suggest that parameterizing PMJ location, rather than number, be given significantly higher priority in future studies creating personalized ventricular models from patient-derived ECGs.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2016

A Theoretical Framework for Calibration in Computer Models: Parametrization, Estimation and Convergence Properties

Calibration parameters in deterministic computer experiments are those attributes that cannot be measured or available in physical experiments. Here, an approach to estimate them by using data from physical experiments and computer simulations. A theoretical framework is given which allows us to study the issues of parameter identifiability and estimation. We define the L 2 -consistency for calibration as a justification for calibration methods. It is shown that a simplified version of the original KO method leads to asymptotically L 2 -inconsistent calibration. This L 2 -inconsistency can be remedied by modifying the original estimation procedure. A novel calibration method, called the L 2 calibration, is proposed and proven to be L 2 -consistent and enjoys optimal convergence rate. Furthermore a numerical example and some mathematical analysis are used to illustrate the source of the L 2 -inconsistency problem.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1978

Energy consumption program: A computer model simulating energy loads in buildings

The JPL energy consumption computer program developed as a useful tool in the on-going building modification studies in the DSN energy conservation project is described. The program simulates building heating and cooling loads and computes thermal and electric energy consumption and cost. The accuracy of computations are not sacrificed, however, since the results lie within + or - 10 percent margin compared to those read from energy meters. The program is carefully structured to reduce both user's time and running cost by asking minimum information from the user and reducing many internal time-consuming computational loops. Many unique features were added to handle two-level electronics control rooms not found in any other program.

Stoller, F. W.↗
NASA NTRS1998

Computer Modeling of the Thermal Conductivity of Cometary Ice

The main objective of this research was to estimate the thermal conductivity of cometry ices from computer simulations of model amorphous ices. This was divided into four specific tasks: (1) Generating samples of amorphous water ices at different microporosities; (2) Comparing the resulting molecular structures of the ices with experimental results, for those densities where data was available; (3) Calculating the thermal conductivities of liquid water and bulk amorphous ices and comparing these results with experimentally determined thermal conductivities; and (4) Investigating how the thermal conductivity of amorphous ice depends upon the microscopic porosity of the samples. The thermal conductivity was found to be only weakly dependent on the microstructure of the amorphous ice. In general, the amorphous ices were found to have thermal conductivities of the same order of magnitude as liquid water. This is in contradiction to recent experimental estimates of the thermal conductivity of amorphous ice, and it is suggested that the extremely low value obtained experimentally is due to larger-scale defects in the ice, such as cracks, but it is not an intrinsic property of the bulk amorphous ice.

Bunch, Theodore E.↗