Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Computer model”

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 829 records · Page 46

Multigrid acceleration and turbulence models for computations of 3D turbulent jets in crossflow

A multigrid method is presented for the calculation of three-dimensional turbulent jets in crossflow. Turbulence closure is achieved with either the standard k-epsilon model or a Reynolds stress model (RSM). Multigrid acceleration enables convergence rates which are far superior to that for a single grid method to be obtained with both turbulence models. With the k-epsilon model the rate approaches that for laminar flow, but with RSM it is somewhat slower. The increased stiffness of the system of equation in the latter may be responsible. Computed results with both turbulence models are compared to experimental data for a pair of opposed jets in crossflow. Both models yield reasonable agreement for the mean flow velocity, but RSM yields better predictions of the Reynolds stresses.

Demuren, A. O.↗

A thermal NO(x) prediction model - Scalar computation module for CFD codes with fluid and kinetic effects

A thermal NO(x) prediction model is developed to interface with a CFD, k-epsilon based code. A converged solution from the CFD code is the input to the postprocessing model for prediction of thermal NO(x). The model uses a decoupled analysis to estimate the equilibrium level of (NO(x))e which is the constant rate limit. This value is used to estimate the flame (NO(x)) and in turn predict the rate of formation at each node using a two-step Zeldovich mechanism. The rate is fixed on the NO(x) production rate plot by estimating the time to reach equilibrium by a differential analysis based on the reaction: O + N2 = NO + N. The rate is integrated in the nonequilibrium time space based on the residence time at each node in the computational domain. The sum of all nodal predictions yields the total NO(x) level.

Mcbeath, Giorgio↗

Computational flow modeling of triply periodic minimal surfaces as feed channel spacers in ultra-high pressure reverse osmosis applications

Triply periodic minimal surfaces (TPMS) are a special class of mathematical surfaces characterized by a high surface area-to-volume ratio. They have generated considerable interest in fields such as acoustics, heat transfer, and membrane-based filtration processes. This study evaluates the performance of four different TPMS designs—Schoen Gyroid, Schoen Crossed Layers of Parallels (CLP), Schoen Transverse Crossed Layers of Parallels (tCLP), and Schwarz-Primitive—when used as feed channel spacers under ultra-high pressure reverse osmosis (UHPRO) conditions, at approximately 200 bar. Our experimentally validated computational fluid dynamics model reveal different flow patterns within the feed channels for each of the four TPMS designs, leading to varying hydrodynamic and permeation properties. Under the simulated UHPRO conditions, the Gyroid and tCLP designs yield up to a 23% increase in average permeate velocity and a 14% reduction in average membrane-surface concentration relative to a non-woven spacer of the same porosity. Furthermore, the enhanced performance comes with an increased feed channel pressure drop, although it only constitutes less than 4% of the operating pressure when extrapolated for a meter-long membrane module. Additionally, the study analyzes the effects of varying inlet velocity and spacer porosity on membrane performance. Overall, this research provides valuable insights into the potential use of TPMS spacers in UHPRO applications.

36 MATERIALS SCIENCE↗

Computing Linear Mathematical Models Of Aircraft

Derivation and Definition of Linear Aircraft Model (LINEAR) computer program provides user with powerful, and flexible, standard, documented, and verified software tool for linearization of mathematical models of aerodynamics of aircraft. Intended for use in software tool to drive linear analysis of stability and design of control laws for aircraft. Capable of both extracting such linearized engine effects as net thrust, torque, and gyroscopic effects, and including these effects in linear model of system. Designed to provide easy selection of state, control, and observation variables used in particular model. Also provides flexibility of allowing alternate formulations of both state and observation equations. Written in FORTRAN.

Duke, Eugene L.↗

A new spherical model for computing the radiation field available for photolysis and heating at twilight

Accurate computation of atmospheric photodissociation and heating rates is needed in photochemical models. These quantities are proportional to the mean intensity of the solar radiation penetrating to various levels in the atmosphere. For large solar zenith angles a solution of the radiative transfer equation valid for a spherical atmosphere is required in order to obtain accurate values of the mean intensity. Such a solution based on a perturbation technique combined with the discrete ordinate method is presented. Mean intensity calculations are carried out for various solar zenith angles. These results are compared with calculations from a plane parallel radiative transfer model in order to assess the importance of using correct geometry around sunrise and sunset. This comparison shows, in agreement with previous investigations, that for solar zenith angles less than 90 deg adequate solutions are obtained for plane parallel geometry as long as spherical geometry is used to compute the direct beam attenuation; but for solar zenith angles greater than 90 deg this pseudospherical plane parallel approximation overstimates the mean intensity.

Dahlback, Arne↗

Description of a mathematical model and computer simulation of separation of the nose cap from the solid rocket booster

A system of equations which models the motion of the Solid Rocket Booster Nose Cap upon separation is described. The computer program which utilizes these equations to generate nose cap trajectories is described in detail. Application of the program to simulate a rocket sled test of the nose cap separation is discussed and the results of the applications are presented. With the information given a user should be able to exercise the computer program with a minimum of effort.

Schwaniger, A. J., Jr.↗

The local wavenumber model for computation of turbulent mixing

We present an overview of the current status in the development of a two-point spectral closure model for turbulent flows, known as the local wavenumber (LWN) model. The model is envisioned as a practical option for applications requiring multi-physics simulations in which statistical hydrodynamics quantities such as Reynolds stresses, turbulent kinetic energy, and measures of mixing such as density-correlations and mix-width evolution, need to be captured with relatively high fidelity. In this review, we present the capabilities of the LWN model since it was first formulated in the early 1990s, for computations of increasing levels of complexity ranging from homogeneous isotropic turbulence, inhomogeneous and anisotropic single-fluid turbulence, to two-species mixing driven by buoyancy forces. The review concludes with a discussion of some of the more theoretical considerations that remain in the development of this model.

42 ENGINEERING↗

Execution environment for intelligent real-time control systems

Modern telerobot control technology requires the integration of symbolic and non-symbolic programming techniques, different models of parallel computations, and various programming paradigms. The Multigraph Architecture, which has been developed for the implementation of intelligent real-time control systems is described. The layered architecture includes specific computational models, integrated execution environment and various high-level tools. A special feature of the architecture is the tight coupling between the symbolic and non-symbolic computations. It supports not only a data interface, but also the integration of the control structures in a parallel computing environment.

Sztipanovits, Janos↗

Fully Three-Dimensional Helical RF Field Effects on TWT Interaction

The phenomenal growth of the satellite communications industry has created a large demand for traveling wave-tubes (TWTs) operating with unprecedented specifications requiring design and production of many novel devices in record time. To achieve this, the TWT industry heavily relies upon computational modeling. However, there is a need for improvement in the TWT industry's computational modeling capabilities as there are often discrepancies between measured TWT data and that predicted by conventional two-dimensional helical TWT interaction codes. This severely limits the analysis and design of novel devices or TWTs with parameters differing from what is conventionally manufactured. In addition, the inaccuracy of current computational tools limits achievable TWT performance because optimized designs require highly accurate models. In particular, collector efficiency, which significantly affects overall TWT efficiency, is difficult to optimize unless extremely accurate spent beam data is predicted by the interaction code.

Kory, Carol L.↗

Performance simulation of the soft gamma-ray concentrator

The soft gamma-ray concentrator is a telescope mission concept utilizing a suitable arrangement of bent multilayer structures of alternating low- and high-density materials. This lens is able to channel gamma-ray photons via total external reflection and concentrate the incident radiation to a point. The channeling technique offers the potential for concentrating gamma rays with focal lengths <10 m and energies >100 keV, beyond the reach of current grazing-incidence hard x-ray mirrors. For the performance estimation of such an instrument, we have developed a flexible set of computer modeling tools to compute the optical properties of multilayer structures, predict the channeling efficiency for a given multilayer configuration, and aid in the optimization of potential gamma-ray concentrator-based telescope designs. This modeling includes the multilayer optical properties calculated by the IMD software, the ray tracing using an IDL code, and the focal plane detector simulation by MEGAlib. We illustrate the potential of this approach by presenting simulated astronomical observations from a balloon-borne platform. The final result, including simulated effective area, instrument sensitivity, and polarization performance, shows that the gamma-ray concentrator will provide greatly increased sensitivity for next-generation soft gamma-ray missions with modest cost and complexity.

47 OTHER INSTRUMENTATION↗

Anti-Icing Coatings using Ionomer Film Layer Structuring

This research effort examined the application of Nafion polymers in alcohol solvents as an anti-ice surface coating, as a mixture with hydrophilic polymers and freezing point depressant salt systems. Co-soluble systems of Nafion, polymer and salt were applied using dip coating methods to create smooth films for frost observation over a Peltier plate thermal system in ambient laboratory conditions. Cryo-DSC was applied to examine freezing events of the Nafion-surfactant mixtures, but the sensitivity of the measurement was insufficient to determine frost behavior. Collaborations with the Fog Chamber at Sandia-Albuquerque, and in environmental SAXS measurements with CINT-LANL were requested but were not able to be performed under the research duration. Since experimental characterization of these factors is difficult to achieve directly, computational modeling was used to guide the scientific basis for property improvement. Computational modeling was performed to improve understanding of the dynamic association between ionomer side groups and added molecules and deicing salts. The polyacrylic acid in water system was identified at the start of the project as a relevant system for exploring the effect of varying counterions on the properties of fully deprotonated polyacrylic acid (PAA) in the presence of water. Simulations were modeled with four different counterions, two monovalent counterions (K+ and Na+) and two divalent counterions (Ca2+ and Mg2+). The wt% of PAA in these systems was varied from ~10 to 80 wt% PAA for temperatures from 250K to 400K. In the second set of simulations, the interpenetration of water into a dry PAA film was studied for Na+ or Ca2+ counterions for temperatures between 300K and 400K. The result of this project is a sprayable Nafion film composite which resists ice nucleation at -20 °C for periods of greater than three hours. It is composed of Nafion polymer, hydrophilic polyethylene oxide polymer and CaCl2 anti-ice crosslinker. Durability and field performance properties remain to be determined.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Atomic oxygen interactions with protected organic materials on the Long Duration Exposure Facility (LDEF)

The Long Duration Exposure Facility (LDEF) has provided an excellent opportunity to understand the nature of directed atomic oxygen interactions with protected polymers and composites. Although there were relatively few samples of materials with protective coatings on their external surfaces on LDEF which were exposed to a high atomic oxygen fluence, analysis of such samples has enabled an examination of the shape of atomic oxygen undercut cavities at defect sites in the protective coatings. Samples of front-surface aluminized (Kapton) polyimide were inspected by scanning electron microscopy to identify and measure crack defects in the aluminum protective coatings. After chemical removal of the aluminum coating, measurements were also made of the width of the oxidized undercut cavities below the crack defects. The LDEF flight undercut cavity geometries were then compared with Monte Carlo computational model undercut cavity predictions. The comparison of the LDEF results and computational modeling indicates agreement in specific undercut cavity geometries for atomic oxygen reaction probabilities dependent upon the 0.68 to 3.0 power of the energy. However, no single energy dependency was adequate to replicate flight results over a variety of aluminum crack widths.

Banks, Bruce A.↗

Non-Invasive Tension Measurement Devices for Parachute Cordage

The need for lightweight and non-intrusive tension measurements has arisen alongside the development of high-fidelity computer models of textile and fluid dynamics. In order to validate these computer models, data must be gathered in the operational environment without altering the design, construction, or performance of the test article. Current measurement device designs rely on severing a cord and breaking the load path to introduce a load cell. These load cells are very reliable, but introduce an area of high stiffness in the load path, directly affecting the structural response, adding excessive weight, and possibly altering the dynamics of the parachute during a test. To capture the required data for analysis validation without affecting the response of the system, non-invasive measurement devices have been developed and tested by NASA. These tension measurement devices offer minimal impact to the mass, form, fit, and function of the test article, while providing reliable, axial tension measurements for parachute cordage.

Litteken, Douglas A.↗

Thermal Management for Planar Package Power Electronics (CRADA Final Report)

The National Renewable Energy Laboratory (NREL) and John Deere Electronic Solutions (JDES) collaborated to develop and evaluate computer models and simulations related to the thermal performance of semiconductor device packaging for an inverter of an off-road vehicle ("Semiconductor Packaging"). The objective of the research was for NREL and JDES to develop a two or three-dimensional computer-aided design model of the Semiconductor Packaging ("Computer Model") for thermal performance evaluation in a simulation. The project developed computer models of a Semiconductor Package of silicon-carbide semiconductor devices with appropriate thermal dissipation via one or more of the following thermal features: a double-sided planar cooling configuration, an air-cooled configuration, a liquid-cooled configuration for off-road inverter applications in a relevant simulated operating environment. It is noted that surface area and prototype product embodiment consisting of air-cooled and liquid-cooled configurations may and may not be in direct contact with power semiconductor chips. The idea was to explore and develop packaging and thermal management technology for power semiconductors that was most effective in the performance yet had least burden in overall product cost for a given application.

33 ADVANCED PROPULSION SYSTEMS↗

Autonomous hybrid optimization of a SiO 2 plasma etching mechanism

Computational modeling of plasma etching processes at the feature scale relevant to the fabrication of nanometer semiconductor devices is critically dependent on the reaction mechanism representing the physical processes occurring between plasma produced reactant fluxes and the surface, reaction probabilities, yields, rate coefficients, and threshold energies that characterize these processes. The increasing complexity of the structures being fabricated, new materials, and novel gas mixtures increase the complexity of the reaction mechanism used in feature scale models and increase the difficulty in developing the fundamental data required for the mechanism. This challenge is further exacerbated by the fact that acquiring these fundamental data through more complex computational models or experiments is often limited by cost, technical complexity, or inadequate models. In this paper, we discuss a method to automate the selection of fundamental data in a reduced reaction mechanism for feature scale plasma etching of SiO 2 using a fluorocarbon gas mixture by matching predictions of etch profiles to experimental data using a gradient descent (GD)/Nelder–Mead (NM) method hybrid optimization scheme. These methods produce a reaction mechanism that replicates the experimental training data as well as experimental data using related but different etch processes.

36 MATERIALS SCIENCE↗

Modeling of Cardiovascular Response to Weightlessness

It was the hypothesis of this Project that the Simple lack of hydrostatic pressure in microgravity generates several purely physical reactions that underlie and may explain, in part, the cardiovascular response to weightlessness. For instance, hydrostatic pressure within the ventricles of the heart may improve cardiac performance by promoting expansion of ventricular volume during diastole. The lack of hydrostatic pressure in microgravity might, therefore, reduce diastolic filling and cardiac performance. The change in transmural pressure is possible due to the difference in hydrostatic pressure gradients between the blood inside the ventricle and the lung tissue surrounding the ventricle due to their different densities. On the other hand, hydrostatic pressure within the vasculature may reduce cardiac inlet pressures because of the typical location of the heart above the hydrostatic indifference level (the level at which pressure remains constant throughout changes in gravity). Additional physical responses of the body to changing gravitational conditions may influence cardiovascular performance. For instance, fluid shifts from the lower body to the thorax in microgravity may serve to increase central venous pressure (CVP) and boost cardiac output (CO). The concurrent release of gravitational force on the rib cage may tend to increase chest girth and decrease pedcardial pressure, augmenting ventricular filling. The lack of gravity on pulmonary tissue may allow an upward shifting of lung mass, causing a further decrease in pericardial pressure and increased CO. Additional effects include diuresis early in the flight, interstitial fluid shifts, gradual spinal extension and movement of abdominal mass, and redistribution of circulatory impedance because of venous distention in the upper body and the collapse of veins in the lower body. In this project, the cardiovascular responses to changes in intraventricular hydrostatic pressure, in intravascular hydrostatic pressure and, to a limited extent, in extravascular and pedcardial hydrostatic pressure were investigated. A complete hydraulic model of the cardiovascular system was built and flown aboard the NASA KC-135 and a computer model was developed and tested in simulated microgravity. Results obtained with these models have confirmed that a simple lack of hydrostatic pressure within an artificial ventricle causes a decrease in stroke volume. When combined with the acute increase in ventricular pressure associated with the elimination of hydrostatic pressure within the vasculature and the resultant cephalad fluid shift with the models in the upright position, however, stroke volume increased in the models. Imposition of a decreased pedcardial pressure in the computer model and in a simplified hydraulic model increased stroke volume. Physiologic regional fluid shifting was also demonstrated by the models. The unifying parameter characterizing of cardiac response was diastolic ventricular transmural pressure (DVDELTAP) The elimination of intraventricular hydrostatic pressure in O-G decreased DVDELTAP stroke volume, while the elimination of intravascular hydrostatic pressure increased DVDELTAP and stroke volume in the upright posture, but reduced DVDELTAP and stroke volume in the launch posture. The release of gravity on the chest wall and its associated influence on intrathoracic pressure, simulated by a drop in extraventricular pressure4, increased DVDELTAP ans stroke volume.

Sharp, M. Keith↗

Experimental and computational analysis of bending fatigue failure in chopped carbon fiber chip reinforced composites

With a better balance among good mechanical performance, high freedom of design, and low material and manufacturing cost, chopped carbon fiber chip reinforced sheet molding compound (SMC) composites show great potential in different engineering applications. Here in this paper, bending fatigue behaviors of SMC composites considering the heterogeneous fiber orientation distributions have been thoroughly investigated utilizing both experimental and computational methods. First, four-point bending fatigue tests are performed with designed SMC composites, and the local modulus is adopted as a metric to represent the local fiber orientation of two opposing sides. Interestingly, SMC composites with and without large discrepancy in local modulus of opposing sides show different fatigue behaviors. Interrupted tests are conducted to explore the bending fatigue failure mechanism, and the damage processes of valid specimens are also closely examined. We find that the fatigue failure of SMC composites under four-point bending is governed by crack propagation instead of crack initiation. Because of this, the heterogeneous local fiber orientations of both sides of the specimen influence fatigue life. The microstructure of the lower side shows a direct influence while that of the upper side also exhibiting influence which becomes more prominent for high cycle fatigue cases. Furthermore, a hybrid micro–macro computational model is proposed to efficiently study the cyclic bending behavior of SMC composites. The region of interest is reconstructed with a modified random sequential absorption algorithm to conserve all the microstructural details including the heterogeneous fiber orientation, while the rest of the regions are modeled as homogenized macro-scale continua. Combined with a framework to capture the progressive fatigue damage under cyclic bending, the bending fatigue behaviors of SMC composites are accurately captured by the hybrid computational model comparing with our experimental analysis.

36 MATERIALS SCIENCE↗

The dynamics and control of large flexible space structures. Part B: Development of continuum model and computer simulation

The equations of motion of an arbitrary flexible body in orbit were derived. The model includes the effects of gravity with all its higher harmonics. As a specific example, the motion of a long, slender, uniform beam in circular orbit was modelled. The example considers both the inplane and three dimensional motion of the beam in orbit. In the case of planar motion with only flexible vibrations, the pitch motion is not influenced by the elastic motion of the beam. For large values of the square of the ratio of the structural modal frequency to the orbital angular rate the elastic motion was decoupled from the pitch motion. However, for small values of the ratio and small amplitude pitch motion, the elastic motion was governed by a Hill's 3 term equation. Numerical simulation of the equation indicates the possibilities of instability for very low values of the square of the ratio of the modal frequency to the orbit angular rate. Also numerical simulations of the first order nonlinear equations of motion for a long flexible beam in orbit were performed. The effect of varying the initial conditions and the number of modes was demonstrated.

Bainum, P. M.↗