Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “computer simulation”

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 1,099 records · Page 61

Octopus, a computational framework for exploring light-driven phenomena and quantum dynamics in extended and finite systems

Over the last few years, extraordinary advances in experimental and theoretical tools have allowed us to monitor and control matter at short time and atomic scales with a high degree of precision. An appealing and challenging route toward engineering materials with tailored properties is to find ways to design or selectively manipulate materials, especially at the quantum level. To this end, having a state-of-the-art ab initio computer simulation tool that enables a reliable and accurate simulation of light-induced changes in the physical and chemical properties of complex systems is of utmost importance. The first principles real-space-based Octopus project was born with that idea in mind, i.e., to provide a unique framework that allows us to describe non-equilibrium phenomena in molecular complexes, low dimensional materials, and extended systems by accounting for electronic, ionic, and photon quantum mechanical effects within a generalized time-dependent density functional theory. This article aims to present the new features that have been implemented over the last few years, including technical developments related to performance and massive parallelism. We also describe the major theoretical developments to address ultrafast light-driven processes, such as the new theoretical framework of quantum electrodynamics density-functional formalism for the description of novel light–matter hybrid states. Those advances, and others being released soon as part of the Octopus package, will allow the scientific community to simulate and characterize spatial and time-resolved spectroscopies, ultrafast phenomena in molecules and materials, and new emergent states of matter (quantum electrodynamical-materials).

74 ATOMIC AND MOLECULAR PHYSICS↗

The microscopic mechanism of bulk melting of ice

We study the initial stages of homogeneous melting of a hexagonal ice crystal at coexistence and at moderate superheating. Our trajectory-based computer simulation approach provides a comprehensive picture of the events that lead to melting, from the initial accumulation of 5+7 defects, via the formation of L-D and interstitial-vacancy pairs, to the formation of a liquid nucleus. Of the different types of defects that we observe to be involved in melting, a particular kind of 5+7 type defect (type 5) plays a prominent role as it often forms prior to the formation of the initial liquid nucleus and close to the site where the nucleus forms. Hence, like other solids, ice homogeneously melts via the prior accumulation of defects.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Visual Computing Environment

The Visual Computing Environment (VCE) is a NASA Lewis Research Center project to develop a framework for intercomponent and multidisciplinary computational simulations. Many current engineering analysis codes simulate various aspects of aircraft engine operation. For example, existing computational fluid dynamics (CFD) codes can model the airflow through individual engine components such as the inlet, compressor, combustor, turbine, or nozzle. Currently, these codes are run in isolation, making intercomponent and complete system simulations very difficult to perform. In addition, management and utilization of these engineering codes for coupled component simulations is a complex, laborious task, requiring substantial experience and effort. To facilitate multicomponent aircraft engine analysis, the CFD Research Corporation (CFDRC) is developing the VCE system. This system, which is part of NASA's Numerical Propulsion Simulation System (NPSS) program, can couple various engineering disciplines, such as CFD, structural analysis, and thermal analysis. The objectives of VCE are to (1) develop a visual computing environment for controlling the execution of individual simulation codes that are running in parallel and are distributed on heterogeneous host machines in a networked environment, (2) develop numerical coupling algorithms for interchanging boundary conditions between codes with arbitrary grid matching and different levels of dimensionality, (3) provide a graphical interface for simulation setup and control, and (4) provide tools for online visualization and plotting. VCE was designed to provide a distributed, object-oriented environment. Mechanisms are provided for creating and manipulating objects, such as grids, boundary conditions, and solution data. This environment includes parallel virtual machine (PVM) for distributed processing. Users can interactively select and couple any set of codes that have been modified to run in a parallel distributed fashion on a cluster of heterogeneous workstations. A scripting facility allows users to dictate the sequence of events that make up the particular simulation.

Lawrence, Charles↗

Transient cavities in liquids and the nature of the hydrophobic effect

The size distributions of transient cavities in water and organic liquids, obtained from computer simulations, have provided a new means to analyze the nature of the hydrophobic effect and to evaluate the adequacy of different analytical models of this effect. The poor solubility of non-polar solutes in water is attributed to a low probability of finding in water cavities of atomic and molecular size. It has been shown that water applies more force per unit area of cavity surface than do hydrocarbon liquids. Models that successfully capture the main characteristics of the hydrophobic effect must at least include information about the density and the radial distribution of oxygen atoms in liquid water. One such model, quantitatively accurate for molecular solutes of arbitrary shape, is presented.

NASA Center ARC↗

The use of a simulator in a synthetic aperture radar data processing tradeoff study.

A computer simulation of a synthetic aperture radar system is described. The objective in developing the simulator was to make it possible to assess the imaging performance of various configurations against a variety of target types. The simulation involves targets, the radar system itself, and the processor; several analysis and output options are provided. Targets may be internally generated or come from external sources, or externally produced (e.g., real) radar signals can be used. A number of parameters of the radar system can be varied. The characteristics of the processor can be adjusted to a considerable extent. Statistical analyses may be performed. Several output options are available, including imagery.

Lybanon, M.↗

Global simulations of the three-dimensional magnetosphere

Global three-dimensional computer simulations of the magnetosphere using a particle MHD code, reproduce the steady-state Dungey magnetospheric topology in three dimensions. The formation of a compression zone downstream of the tail neutral line that is probably bounded by wake shocks is observed. This compression zone changes its cross-section with distance downstream.

Leboeuf, J. N.↗

Simulation of the electron acoustic instability in the polar cusp

A computer simulation of the beam-driven electron acoustic instability in the polar cusp is used to investigate the temporal evolution of typical cusp electron distribution functions in self-consistently generated wave fields. The simulation results are compared with linear, second-order, and nonlinear theory. The behavior of the instability over a range of plasma parameters, such as beam speed and the cold electron density, is shown. The saturation mechanism of the instability is examined, and the relevance of the results for observations of cusp hiss are discussed.

Lin, C. S.↗

An investigation of the information propagation and entropy transport aspects of Stirling machine numerical simulation

Aspects of the information propagation modeling behavior of integral machine computer simulation programs are investigated in terms of a transmission line. In particular, the effects of pressure-linking and temporal integration algorithms on the amplitude ratio and phase angle predictions are compared against experimental and closed-form analytic data. It is concluded that the discretized, first order conservation balances may not be adequate for modeling information propagation effects at characteristic numbers less than about 24. An entropy transport equation suitable for generalized use in Stirling machine simulation is developed. The equation is evaluated by including it in a simulation of an incompressible oscillating flow apparatus designed to demonstrate the effect of flow oscillations on the enhancement of thermal diffusion. Numerical false diffusion is found to be a major factor inhibiting validation of the simulation predictions with experimental and closed-form analytic data. A generalized false diffusion correction algorithm is developed which allows the numerical results to match their analytic counterparts. Under these conditions, the simulation yields entropy predictions which satisfy Clausius' inequality.

Goldberg, Louis F.↗

Synthesis of Contributed Simulations for OREX Test Cases

A synthesis is presented of the computer simulations of the flow over the Orbital Reentry Vehicle (ORV) at the 92.8 km and 63.6 km Earth altitude trajectory points that were discussed at the First Europe-U.S. High Speed Flow Field Database Workshop Part 2, Napoli, Italy, November 1997. For the materials used on the surface of ORV, the non-catalytic wall condition is appropriate at 92.8 km and the finite-rate catalytic wall condition at 63.6 km. Additional simulations are required for establishing the independency of the discussed results from numerics. The proper modeling of natural phenomena needs further sensitivity studies. The uncertainties of inferred flight data are lacking for a proper evaluation of the presented results.

Mehta, Unmeel B.↗

Crowding-Induced Spatial Organization of Gene Expression in Cell-Sized Vesicles

Cell-free protein synthesis is an important tool for studying gene expression and harnessing it for applications. In cells, gene expression is regulated in part by the spatial organization of transcription and translation. Unfortunately, current cell-free approaches are unable to control the organization of molecular components needed for gene expression, which limits the ability to probe and utilize its effects. In this work, we show, using complementary computational and experimental approaches, that macromolecular crowding can be used to control the spatial organization and translational efficiency of gene expression in cell-sized vesicles. Computer simulations and imaging experiments reveal that, as crowding is increased, DNA plasmids become localized at the inner surface of vesicles. Ribosomes, in contrast, remain uniformly distributed, demonstrating that crowding can be used to differentially organize components of gene expression. We further carried out cell-free protein synthesis reactions in cell-sized vesicles and quantified mRNA and protein abundance. At sufficiently high levels of crowding, we observed localization of mRNA near vesicle surfaces, a decrease in translational efficiency and protein abundance, and anomalous scaling of protein abundance as a function of vesicle size. These results are consistent with high levels of crowding causing altered spatial organization and slower diffusion. Our work demonstrates a straightforward way to control the organization of gene expression in cell-sized vesicles and provides insight into the spatial regulation of gene expression in cells.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Theory of hydrophobicity: transient cavities in molecular liquids

Observation of the size distribution of transient cavities in computer simulations of water, n-hexane, and n-dodecane under benchtop conditions shows that the sizes of cavities are more sharply defined in liquid water but the most-probable-size cavities are about the same size in each of these liquids. The calculated solvent atomic density in contact with these cavities shows that water applies more force per unit area of cavity surface than do the hydrocarbon liquids. This contact density, or "squeezing" force, reaches a maximum near cavity diameters of 2.4 angstroms. The results for liquid water are compared to the predictions of simple theories and, in addition, to results for a reference simple liquid. The numerical data for water at a range of temperatures are analyzed to extract a surface free energy contribution to the work of formation of atomic-size cavities. Comparison with the liquid-vapor interfacial tensions of the model liquids studied here indicates that the surface free energies extracted for atomic-size cavities cannot be accurately identified with the macroscopic surface tensions of the systems.

NASA Center ARC↗

Dataflow computing approach in high-speed digital simulation

New computational tools and methodologies for the digital simulation of continuous systems were explored. Programmability, and cost effective performance in multiprocessor organizations for real time simulation was investigated. Approach is based on functional style languages and data flow computing principles, which allow for the natural representation of parallelism in algorithms and provides a suitable basis for the design of cost effective high performance distributed systems. The objectives of this research are to: (1) perform comparative evaluation of several existing data flow languages and develop an experimental data flow language suitable for real time simulation using multiprocessor systems; (2) investigate the main issues that arise in the architecture and organization of data flow multiprocessors for real time simulation; and (3) develop and apply performance evaluation models in typical applications.

Ercegovac, M. D.↗

Helicopter flight control system design and evaluation for NOE operations using controller inversion techniques

A flight control system design technique is described that allows rapid evaluation of candidate control laws using computer simulated trajectories. The two-step technique first computes the optimum trajectory for the basic aircraft, and then reconstructs the pilot's control displacements necessary to fly this trajectory, given a particular control system architecture. Since the flight control system itself does not add to the acceleration potential of the helicopter, one may evaluate several proposed control laws in terms of the resulting stick inputs for the same maneuver. The method is illustrated through application to implicit and explicit model-following control laws designed for a battlefield helicopter with and without an auxiliary propulsor. It is shown that the technique is a valuable aid for selection of promising control laws for further simulator studies and flight investigations.

Mckillip, Robert M., Jr.↗

Kinetic dissection of individual steps in the poly(C)-directed oligoguanylate synthesis from guanosine 5'-monophosphate 2-methylimidazolide

A kinetic study of oligoguanylate synthesis on a polycytidylate template, poly(C), as a function of the concentration of the activated monomer, guanosine 5'-monophosphate 2-methylimidazolide, 2-MeImpG, is reported. Reactions were run with 0.005-0.045 M 2-MeImpG in the presence of 0.05 M poly(C) at 23 degrees C. The kinetic results are consistent with a reaction scheme (eq 1) that consists of a series of consecutive steps, each step representing the addition of one molecule of 2-MeImpG to the growing oligomer. This scheme allows the calculation of second-order rate constants for every step by analyzing the time-dependent growth of each oligomer. Computer simulations of the course of reaction based on the determined rate constants and eq 1 are in excellent agreement with the product distributions seen in the HPLC profiles. In accord with an earlier study (Fakhrai, H.; Inoue, T.; Orgel, L. E. Tetrahedron 1984, 40, 39), rate constants, ki, for the formation of the tetramer and longer oligomers up to the 16-mer were found to be independent of length and somewhat higher than k3 (formation of trimer), which in turn is much higher than k2 (formation of dimer). The ki (i > or = 4), k3, and k2 values are not true second-order rate constants but vary with monomer concentration. Mechanistic models for the dimerization (Scheme I) and elongation reactions (Scheme II) are proposed that are consistent with our results. These models take into account that the monomer associates with the template in a cooperative manner. Our kinetic analysis allowed the determination of rate constants for the elementary processes of covalent bond formation between two monomers (dimerization) and between an oligomer and a monomer (elongation) on the template. A major conclusion from our study is that bond formation between two monomer units or between a primer and a monomer is assisted by the presence of additional next-neighbor monomer units. This is consistent with recent findings with hairpin oligonucleotides (Wu, T.; Orgel, L. E. J. Am. Chem. Soc. 1992, 114, 317). Our study is the first of its kind that shows the feasibility of a thorough kinetic analysis of a template-directed oligomerization and provides a detailed mechanistic model of these reactions.

Non-NASA Center↗

Kinetics of DSB rejoining and formation of simple chromosome exchange aberrations

PURPOSE: To investigate the role of kinetics in the processing of DNA double strand breaks (DSB), and the formation of simple chromosome exchange aberrations following X-ray exposures to mammalian cells based on an enzymatic approach. METHODS: Using computer simulations based on a biochemical approach, rate-equations that describe the processing of DSB through the formation of a DNA-enzyme complex were formulated. A second model that allows for competition between two processing pathways was also formulated. The formation of simple exchange aberrations was modelled as misrepair during the recombination of single DSB with undamaged DNA. Non-linear coupled differential equations corresponding to biochemical pathways were solved numerically by fitting to experimental data. RESULTS: When mediated by a DSB repair enzyme complex, the processing of single DSB showed a complex behaviour that gives the appearance of fast and slow components of rejoining. This is due to the time-delay caused by the action time of enzymes in biomolecular reactions. It is shown that the kinetic- and dose-responses of simple chromosome exchange aberrations are well described by a recombination model of DSB interacting with undamaged DNA when aberration formation increases with linear dose-dependence. Competition between two or more recombination processes is shown to lead to the formation of simple exchange aberrations with a dose-dependence similar to that of a linear quadratic model. CONCLUSIONS: Using a minimal number of assumptions, the kinetics and dose response observed experimentally for DSB rejoining and the formation of simple chromosome exchange aberrations are shown to be consistent with kinetic models based on enzymatic reaction approaches. A non-linear dose response for simple exchange aberrations is possible in a model of recombination of DNA containing a DSB with undamaged DNA when two or more pathways compete for DSB repair.

DNA Damage↗

Viewpoint 9--molecular structure of aqueous interfaces

In this review we summarize recent progress in our understanding of the structure of aqueous interfaces emerging from molecular level computer simulations. It is emphasized that the presence of the interface induces specific structural effects which, in turn, influence a wide variety of phenomena occurring near the phase boundaries. At the liquid-vapor interface, the most probable orientations of a water molecule is such that its dipole moment lies parallel to the interface, one O-H bond points toward the vapor and the other O-H bond is directed toward the liquid. The orientational distributions are broad and slightly asymmetric, resulting in an excess dipole moment pointing toward the liquid. These structural preferences persist at interfaces between water and nonpolar liquids, indicating that the interactions between the two liquids in contact are weak. It was found that liquid-liquid interfaces are locally sharp but broadened by capillary waves. One consequence of anisotropic orientations of interfacial water molecules is asymmetric interactions, with respect to the sign of the charge, of ions with the water surface. It was found that even very close to the surface ions retain their hydration shells. New features of aqueous interfaces have been revealed in studies of water-membrane and water-monolayer systems. In particular, water molecules are strongly oriented by the polar head groups of the amphiphilic phase, and they penetrate the hydrophilic head-group region, but not the hydrophobic core. At infinite dilution near interfaces, amphiphilic molecules exhibit behavior different from that in the gas phase or in bulk water. This result sheds new light on the nature of hydrophobic effect in the interfacial regions. The presence of interfaces was also shown to affect both equilibrium and dynamic components of rates of chemical reactions. Applications of continuum models to interfacial problems have been, so far, unsuccessful. This, again, underscores the importance of molecular-level information about interfaces.

Review, Tutorial↗

Flexible rotor dynamics analysis

A digital computer program was developed to analyze the general nonaxisymmetric and nonsynchronous transient and steady-state rotor dynamic performance of a bending- and shear-wise flexible rotor-bearing system under various operating conditions. The effects of rotor material mechanical hysteresis, rotor torsion flexibility, transverse effects of rotor axial and torsional loading and the anisotropic, in-phase and out-of-phase bearing stiffness and damping force and moment coefficients were included in the program to broaden its capability. An optimum solution method was found and incorporated in the computer program. Computer simulation of experimental data was made and qualitative agreements observed. The mathematical formulations, computer program verification, test data simulation, and user instruction was presented and discussed.

Shen, F. A.↗

Combined Experimental and Numerical Simulations of Thermal Barrier Coated Turbine Blades Erosion

A combined experimental and computational study was conducted to investigate the erosion of thermal barrier coated (TBC) blade surfaces by alumina particles ingestion in a single stage turbine. In the experimental investigation, tests of particle surface interactions were performed in specially designed tunnels to determine the erosion rates and particle restitution characteristics under different impact conditions. The experimental results show that the erosion rates increase with increased impingement angle, impact velocity and temperature. In the computational simulations, an Euler-Lagrangian two stage approach is used in obtaining numerical solutions to the three-dimensional compressible Reynolds Averaged Navier-Stokes equations and the particles equations of motion in each blade passage reference frame. User defined functions (UDF) were developed to represent experimentally-based correlations for particle surface interaction models which were employed in the three-dimensional particle trajectory simulations to determine the particle rebound characteristics after each surface impact. The experimentally based erosion UDF model was used to predict the TBC erosion rates on the turbine blade surfaces based on the computed statistical data of the particles impact locations, velocities and angles relative to the blade surface. Computational results are presented for the predicted TBC blade erosion in a single stage commercial APU turbine, for a NASA designed automotive turbine, and for the NASA turbine scaled for modern rotorcraft operating conditions. The erosion patterns in the turbines are discussed for uniform particle ingestion and for particle ingestion concentrated in the inner and outer 5 percent of the stator blade span representing the flow cooling the combustor liner.

Hamed, Awate↗