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 433 records · Page 24

A computer model of the earth's magnetosphere

The interaction of the solar wind with the earth magnetosphere is investigated theoretically by means of three-dimensional MHD simulations, with a focus on the effects of changes in the Bz component of the IMF. A high-resolution (0.5 earth radii) version of the model of Ogino et al. (1986) is employed, and the results are presented in a series of computer-generated maps and diagrams and characterized in detail. Bz of -5 nT is found to be associated with dipolar magnetic-field lines near the earth and very concave lines in the magnetotail, while Bz of +5 nT produces a narrow finger of closed field lines extending into the polar cap. Both IMF orientations have sunward convection near the noon-midnight meridian and region-1-type field-aligned currents on both sides of the plasma-sheet extension.

Ogino, Tatsuki↗

Computational Modeling for Non-equilibrium Shock Tube Flows

Time accurate simulation of non-equilibrium flows inside shock tube facilities presents several challenges from both physical and mathematical aspects. Furthermore, the drastic computational cost makes it non-practical to support real-time experimental test campaign. In this work, we explore other methods for modeling the shock tube prob- lem with the main focus on the post-shock region and the absolute radiation emanating from it. The proposed alternative approach is several orders of magnitude less computa- tionally expansive while still accurate enough with regards to the quantities of interest. Excellent agreement is found with the well-established stagnation-line approach. Comparison with the time-accurate simulation shows good agreement close to the peak values and disagreement of the temperatures relaxation and radiance profiles toward equilibrium, due to shock speed unsteadiness.

Bensassi, Khalil↗

Full Scale 3D Computational Model of the Industrial -Scale Coal Fired Boiler Performance for Temperature Sensor Installation Guidance

Abstract Nearly 30% of the electricity is generated by using coal as the primary fuel in the US. One of the major concerns in coal-fired power plants is the failure of boiler tubes that leads to unscheduled maintenance and has a huge economical and societal impact. High temperature flue gas along with ash pass over the boiler tubes, which over time leads to tube failure. Therefore, developing temperature sensors for harsh environments and install them for temperature sensing and boiler tube lifetime prediction is an urgent need. On the side of sensor development, the location of the sensor installation is important for stable sensing performance and easy calibration. In this study, computational fluid dynamics and heat transfer modeling are adopted to establish a full-scale 3-dimensional model of a coal-fired boiler to investigate the flue gas temperature distribution within the boiler and identify the proper locations for sensor installation. We proposed three criteria to select the temperature sensor installation location: (1) select the boiler tube panel away from the sidewalls, (2) select the boiler tube section closer to the top wall of the boiler; and (3) select the boiler tube on the back of the boiler panel (not directly facing the flue gas flow). In these regions, the flue gas temperature is stable, providing an ideal environment for stable temperature sensing and calibration.

Gupta, Tanuj↗

Comparison of ERBE inferred and model computed clear-sky albedos

Over-ocean clear-sky albedos measured with instruments on the Earth Radiation Budget Satellite (ERBS) are compared with albedos simulated using a radiative transfer model (RTM). The comparison covers the monthly mean albedos for November 1984. The ERBS albedo was calculated with a scene identification algorithm. Techniques used to suppress cloud cover uncertainties are discussed. The plane-parallel delta-Eddington RTM accounted for O3, O2, CO2 and H2O gaseous absorption and background aerosol absorption.

Briegleb, B. P.↗

Electromagnetic and two-photon transition form factors of the pseudoscalar mesons: An algebraic model computation

We compute electromagnetic and two-photon transition form factors of ground-state pseudoscalar mesons: π , K , η c , η b . To this end, we employ an algebraic model based upon the coupled formalism of Schwinger-Dyson and Bethe-Salpeter equations. Within this approach, the dressed quark propagator and the relevant Bethe-Salpeter amplitude encode the internal structure of the corresponding meson. Electromagnetic properties of the meson are probed via the quark-photon interaction. The algebraic model employed by us unifies the treatment of all ground-state pseudoscalar mesons. Its parameters are carefully fitted performing a global analysis of existing experimental data including the knowledge of the charge radii of the mesons studied. We then compute and predict electromagnetic and two-photon transition form factors for a wide range of probing photon momentum-squared which is of direct relevance to the experimental observations carried out thus far or planned at different hadron physics facilities such as the Thomas Jefferson National Accelerator Facility (JLab) and the forthcoming Electron-Ion Collider. We also present comparisons with other theoretical models and approaches and lattice quantum chromodynamics. Published by the American Physical Society 2024

Higuera-Angulo, I. M. (ORCID:0000000256008875)↗

Understanding Microstructural Evolution During Rapid Heat Treatment of Microalloyed Steels Through Computational Modeling, Advanced Physical Simulation, and Multiscale Characterization Techniques

An AISI 1045 steel modified with vanadium (V) and niobium (Nb) was studied to evaluate microstructural conditioning prior to and throughout a rapid heat treat process. In order to accomplish this, both computational and physical simulation techniques have been employed with the goal of assessing the microstructural evolution in a medium-carbon bar steel during the rapid austenitization and quenching procedures involved in an induction hardening process. The appropriate thermal profiles for induction hardening were obtained through finite element modeling using Flux 2D software. Physical simulations of the induction hardening process were carried out using a Gleeble{sup ®} 3500. Analysis of prior austenite grain size is complemented by observation of nanoscale carbonitride precipitation via transmission electron microscopy, scanning transmission electron microscopy, and high-energy synchrotron small-angle x-ray scattering. Through a combination of characterization techniques, this study presents a deeper understanding of nano- and microstructural changes occurring in a microalloyed steel during an induction hardening process.

36 MATERIALS SCIENCE↗

Formulation and Validation of an Efficient Computational Model for a Dilute, Settling Suspension Undergoing Rotational Mixing

Designing processing equipment for the mixing of settling suspensions is a challenging problem. Achieving low-cost mixing is especially difficult for the application of slowly reacting suspended solids because the cost of impeller power consumption becomes quite high due to the long reaction times (batch mode) or due to large-volume reactors (continuous mode). Further, the usual scale-up metrics for mixing, e.g., constant tip speed and constant power per volume, do not apply well for mixing of suspensions. As an alternative, computational fluid dynamics (CFD) can be useful for analyzing mixing at multiple scales and determining appropriate mixer designs and operating parameters. We developed a mixture model to describe the hydrodynamics of a settling cellulose suspension. The suspension motion is represented as a single velocity field in a computationally efficient Eulerian framework. The solids are represented by a scalar volume-fraction field that undergoes transport due to particle diffusion, settling, fluid advection, and shear stress. A settling model and a viscosity model, both functions of volume fraction, were selected to fit experimental settling and viscosity data, respectively. Simulations were performed with the open-source Nek5000 CFD program, which is based on the high-order spectral-finite-element method. Simulations were performed for the cellulose suspension undergoing mixing in a laboratory-scale vane mixer. The settled-bed heights predicted by the simulations were in semi-quantitative agreement with experimental observations. Further, the simulation results were in quantitative agreement with experimentally obtained torque and mixing-rate data, including a characteristic torque bifurcation. In future work, we plan to couple this CFD model with a reaction-kinetics model for the enzymatic digestion of cellulose, allowing us to predict enzymatic digestion performance for various mixing intensities and novel reactor designs.

09 BIOMASS FUELS↗

Comprehensive silicon solar cell computer modeling

The development of an efficient, comprehensive Si solar cell modeling program that has the capability of simulation accuracy of 5 percent or less is examined. A general investigation of computerized simulation is provided. Computer simulation programs are subdivided into a number of major tasks: (1) analytical method used to represent the physical system; (2) phenomena submodels that comprise the simulation of the system; (3) coding of the analysis and the phenomena submodels; (4) coding scheme that results in efficient use of the CPU so that CPU costs are low; and (5) modularized simulation program with respect to structures that may be analyzed, addition and/or modification of phenomena submodels as new experimental data become available, and the addition of other photovoltaic materials.

Lamorte, M. F.↗

Experimental Verification of Computational Models for Laminated Composites

The objective of the research reported herein is to develop a progressive damage methodology capable of predicting the residual strength of continuous fiber-reinforced, laminated, polymer matrix composites with through-penetration damage. The fracture behavior of center-notch tension panels with thin crack-like slits was studied. Since fibers are the major load-carrying constituent in polymer matrix composites, predicting the residual strength of a laminate requires a criterion for fiber fracture. The effects on fiber strain due to other damage mechanisms such as matrix cracking and delaminations must also be modeled. Therefore, the research herein examines the damage mechanisms involved in translaminate fracture and identifies the toughening mechanisms responsible for damage growth resistance in brittle epoxy matrix systems. The mechanics of matrix cracking and fiber fracture are discussed as is the mathematical framework for the progressive damage model developed by the authors. The progressive damage analysis algorithms have been implemented into a general purpose finite element code developed by NASA, the Computational Structural Mechanics Testbed (COMET). Damage growth is numerically simulated and the analytical residual strength predictions are compared to experimental results for a variety of notched panel configurations and materials systems.

Harris, Charles E.↗

Comparison of Analytical Predictions and Experimental Results for a Dual Brayton Power System (Discussion on Test Hardware and Computer Model for a Dual Brayton System)

NASA Glenn Research Center (GRC) contracted Barber-Nichols, Arvada, CO to construct a dual Brayton power conversion system for use as a hardware proof of concept and to validate results from a computational code known as the Closed Cycle System Simulation (CCSS). Initial checkout tests were performed at Barber- Nichols to ready the system for delivery to GRC. This presentation describes the system hardware components and lists the types of checkout tests performed along with a couple issues encountered while conducting the tests. A description of the CCSS model is also presented. The checkout tests did not focus on generating data, therefore, no test data or model analyses are presented.

Johnson, Paul K.↗

Two-Dimensional Computational Model for Wave Rotor Flow Dynamics

A two-dimensional (theta,z) Navier-Stokes solver for multi-port wave rotor flow simulation is described. The finite-volume form of the unsteady thin-layer Navier-Stokes equations are integrated in time on multi-block grids that represent the stationary inlet and outlet ports and the moving rotor passages of the wave rotor. Computed results are compared with three-port wave rotor experimental data. The model is applied to predict the performance of a planned four-port wave rotor experiment. Two-dimensional flow features that reduce machine performance and influence rotor blade and duct wall thermal loads are identified. The performance impact of rounding the inlet port wall, to inhibit separation during passage gradual opening, is assessed.

Welch, Gerard E.↗

User's manual for the REEDM (Rocket Exhaust Effluent Diffusion Model) computer program

The REEDM computer program predicts concentrations, dosages, and depositions downwind from normal and abnormal launches of rocket vehicles at NASA's Kennedy Space Center. The atmospheric dispersion models, cloud-rise models, and other formulas used in the REEDM model are described mathematically Vehicle and source parameters, other pertinent physical properties of the rocket exhaust cloud, and meteorological layering techniques are presented as well as user's instructions for REEDM. Worked example problems are included.

Bjorklund, J. R.↗

Hollow cathodes as electron emitting plasma contactors Theory and computer modeling

Several researchers have suggested using hollow cathodes as plasma contactors for electrodynamic tethers, particularly to prevent the Shuttle Orbiter from charging to large negative potentials. Previous studies have shown that fluid models with anomalous scattering can describe the electron transport in hollow cathode generated plasmas. An improved theory of the hollow cathode plasmas is developed and computational results using the theory are compared with laboratory experiments. Numerical predictions for a hollow cathode plasma source of the type considered for use on the Shuttle are presented, as are three-dimensional NASCAP/LEO calculations of the emitted ion trajectories and the resulting potentials in the vicinity of the Orbiter. The computer calculations show that the hollow cathode plasma source makes vastly superior contact with the ionospheric plasma compared with either an electron gun or passive ion collection by the Orbiter.

Davis, V. A.↗

Recent evolution of the offline computing model of the NOvA experiment

The NOvA experiment at Fermilab is a long-baseline neutrino experiment designed to study ν e appearance in a ν μ beam. Over the last few years there has been intense work to streamline the computing infrastructure in preparation for data, which started to flow in from the far detector in Fall 2013. Major accomplishments for this effort include migration to the use of off-site resources through the use of the Open Science Grid and upgrading the file-handling framework from simple disk storage to a tiered system using a comprehensive data management and delivery system to find and access files on either disk or tape storage. NOvA has already produced more than 6.5 million files and more than 1 PB of raw data and Monte Carlo simulation files which are managed under this model. In addition, the current system has demonstrated sustained rates of up to 1 TB/hour of file transfer by the data handling system. NOvA pioneered the use of new tools and this paved the way for their use by other Intensity Frontier experiments at Fermilab. Most importantly, the new framework places the experiment's infrastructure on a firm foundation, and is ready to produce the files needed for first physics.

97 MATHEMATICS AND COMPUTING↗

Computational modelling of Er(3+): Garnet laser materials

The Er(3+) ion has attracted a lot of interest for four reasons: (1) Its (4)I(sub 13/2) yields (4)I(sub 15/2) transition lases in the eyesafe region near 1.5 micron; (2) the (4)I(sub 13/2) transition lases near 2.8 micron, an important wavelength for surgical purposes; (3) it displays surprisingly efficient upconversion with lasing observed at 1.7, 1.2, 0.85, 0.56, 0.55, and 0.47 micron following 1.5 micron pumping; and (4) it has absorption bands at 0.96 and 0.81 micron and thus can be diode pumped. However, properties desirable for upconversion reduce the efficiency of 1.5 and 3 micron laser operation and vice versa. Since all of the processes are influenced by the host via the crystal field induced stark splittings in the Er levels, this project undertook modelling of the host influence on the Er lasinng behavior. While growth and measurement of all ten Er(3+) doped garnets is the surest way of identifying hosts which maximize upconversion (or conversly, 1.5 and 3 micron performance), it is also expensive - costing approximately $10,000/material or approximately $100,000 for the materials computationally investigated here. The calculations were performed using a quantum mechanical point charge model developed by Clyde Morrison at Harry Diamond Laboratories. The programs were used to fit the Er:YAG experimental energy levels so that the crystal field parameters, B(sub nm) could be extracted. From these radial factors, rho (sub n) were determined for Er(3+) in garnets. These, in combination with crystal field components, Anm, available from X-ray data, were used to predict energy levels for Er in the other nine garnet hosts. The levels in Er:YAG were fit with an rms error of 12.2/cm over a 22,000/cm range. Predicted levels for two other garnets for which literature values were available had rms errors of less than 17/cm , showing the calculations to be reliable. Based on resonances between pairs of calculated stark levels, the model predicts GSGG as the best host for 1.5 micron laser operation, GSGG or YSAG as the best host for a 2.8 micron operation, and LuGG as the best host for an upconversion material.

Spangler, Lee H.↗

Free energy of RNA-counterion interactions in a tight-binding model computed by a discrete space mapping

The thermodynamic stability of a folded RNA is intricately tied to the counterions and the free energy of this interaction must be accounted for in any realistic RNA simulations. Extending a tight-binding model published previously, in this paper we investigate the fundamental structure of charges arising from the interaction between small functional RNA molecules and divalent ions such as Mg{sup 2+} that are especially conducive to stabilizing folded conformations. The characteristic nature of these charges is utilized to construct a discretely connected energy landscape that is then traversed via a novel application of a deterministic graph search technique. This search method can be incorporated into larger simulations of small RNA molecules and provides a fast and accurate way to calculate the free energy arising from the interactions between an RNA and divalent counterions. The utility of this algorithm is demonstrated within a fully atomistic Monte Carlo simulation of the P4-P6 domain of the Tetrahymena group I intron, in which it is shown that the counterion-mediated free energy conclusively directs folding into a compact structure.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗