SEARCH · Search NASA
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.
The massively parallel processor is used to perform Monte Carlo simulations for the two dimensional XY model on lattices of sizes up to 128 x 128. A parallel random number generator was constructed, finite size effects were studied, and run times were compared with those on a CRAY X-MP supercomputer.
Lattice filters have properties which make them attractive for the application to the identification of the parameters of flexible structures. Here, the problem of both order and parameter estimation is addressed, where the numerical difficulties of an embedded lattice filter are avoided by applying a periodic lattice filter. With this formulation, it is quite natural to estimate system order by applying Akaike's (1976) information criterion to each channel of the periodic lattice filter. Several examples including nonstationary cases are simulated to demonstrate these facts.
In order to identify the optimum size of macrocrown ether for threading, we first investigated the size and shape of simple crown ethers in the melt at 373 K, and their extent of threading with PEO in the melt using coarse-grained Monte Carlo simulations on the 2nnd (second nearest neighbor diamond) lattice, which is a high coordination lattice whose coarse-grained chains can be reverse mapped into fully atomistic models in continuous space.
We are developing a multiscale framework in computational modeling for the ultra high temperature ceramics (UHTC) ZrB2 and HfB2. These materials are characterized by high melting point, good strength, and reasonable oxidation resistance. They are candidate materials for a number of applications in extreme environments including sharp leading edges of hypersonic aircraft. In particular, we used a combination of ab initio methods, atomistic simulations and continuum computations to obtain insights into fundamental properties of these materials. Ab initio methods were used to compute basic structural, mechanical and thermal properties. From these results, a database was constructed to fit a Tersoff style interatomic potential suitable for atomistic simulations. These potentials were used to evaluate the lattice thermal conductivity of single crystals and the thermal resistance of simple grain boundaries. Finite element method (FEM) computations using atomistic results as inputs were performed with meshes constructed on SEM images thereby modeling the realistic microstructure. These continuum computations showed the reduction in thermal conductivity due to the grain boundary network.
A finite element model of cardiac conduction was used to simulate two experimental protocols: 1) fibrillation threshold measurements and 2) clinical electrophysiologic (EP) testing procedures. The model consisted of a cylindrical lattice whose properties were determined by four parameters: element length, conduction velocity, mean refractory period, and standard deviation of refractory periods. Different stimulation patterns were applied to the lattice under a given set of lattice parameter values and the response of the model was observed through a simulated electrocardiogram. The studies confirm that the model can account for observations made in experimental fibrillation threshold measurements and in clinical EP testing protocols.
A specialized, microstructural lattice model, termed MCFET for combined Monte Carlo Finite Element Technique, was developed which simulates microstructural evolution in material systems where modulated phases occur and the directionality of the modulation is influenced by internal and external stresses. In this approach, the microstructure is discretized onto a fine lattice. Each element in the lattice is labelled in accordance with its microstructural identity. Diffusion of material at elevated temperatures is simulated by allowing exchanges of neighboring elements if the exchange lowers the total energy of the system. A Monte Carlo approach is used to select the exchange site while the change in energy associated with stress fields is computed using a finite element technique. The MCFET analysis was validated by comparing this approach with a closed form, analytical method for stress assisted, shape changes of a single particle in an infinite matrix. Sample MCFET analytical for multiparticle problems were also run and in general the resulting microstructural changes associated with the application of an external stress are similar to that observed in Ni-Al-Cr alloys at elevated temperature.
Two types of ferromagnets, pure iron and steel with varying geometry and microstructural properties, were prepared for the present study of magnetoacoustic emission (MAE). The purpose was to separate the effects of structural property variations from changes caused by differences in the sample geometry. The position and shape of the leading MAE sub-burst and its variation among the samples are explained by magnetic anisotropy and the results of numerical simulations which utilized the concept of self-organized criticality (SOC). The amplitude and duration of the second sub-burst, which previously was thought to occur as a result of a complicated interaction between non-180 deg domain walls and lattice defects, can easily be explained by the results of our simulation.
The Bozzolo, Ferrante and Smith (BFS) method for alloys is applied to the study of the defect structure of B2 FeAI alloys. First-principles Linear Muffin Tin Orbital calculations are used to determine the input parameters to the BFS method used in this work. The calculations successfully determine the phase field of the B2 structure, as well as the dependence with composition of the lattice parameter. Finally, the method is used to perform 'static' simulations where instead of determining the ground state configuration of the alloy with a certain concentration of vacancies, a large number of candidate ordered structures are studied and compared, in order to determine not only the lowest energy configurations but other possible metastable states as well. The results provide a description of the defect structure consistent with available experimental data. The simplicity of the BFS method also allows for a simple explanation of some of the essential features found in the concentration dependence of the heat of formation, lattice parameter and the defect structure.
This paper presents a specialized microstructural lattice model, MCFET (Monte Carlo finite element technique), which simulates microstructural evolution in materials in which strain energy has an important role in determining morphology. The model is capable of accounting for externally applied stress, surface tension, misfit, elastic inhomogeneity, elastic anisotropy, and arbitrary temperatures. The MCFET analysis was found to compare well with the results of analytical calculations of the equilibrium morphologies of isolated particles in an infinite matrix.
The evolution of domain walls in the early universe is studied via 2-D computer simulation. The walls are initially configured on a triangular lattice and then released from the lattice, their evolution driven by wall curvature and by the universal expansion. The walls attain an average velocity of about 0.3c and their surface area per volume (as measured in comoving coordinates) goes down with a slope of -1 with respect to conformal time, regardless of whether the universe is matter or radiation dominated. The additional influence of vacuum pressure causes the energy density to fall away from this slope and steepen, thus allowing a situation in which domain walls can constitute a significant portion of the energy density of the universe without provoking an unacceptably large perturbation upon the microwave background.
This paper discusses the implementation of an adaptive acoustic echo canceler for a hands-free cellular phone operating on a fading channel. The adaptive lattice structure, which is particularly known for faster convergence relative to the conventional tapped-delay-line (TDL) structure, is used in the initialization stage. After convergence, the lattice coefficients are converted into the coefficients for the TDL structure which can accommodate a larger number of taps in real-time operation due to its computational simplicity. The conversion method of the TDL coefficients from the lattice coefficients is derived and the DSP56001 assembly code for the lattice and TDL structure is included, as well as simulation results and the schematic diagram for the hardware implementation.
This paper builds on the computational aeroelastic results published previously and generated in support of the second Aeroelastic Prediction Workshop for the NASA Benchmark Supercritical Wing configuration. The computational results are obtained using FUN3D, an unstructured flow solver developed at the NASA Langley Research Center. Flutter analyses are performed using various aerodynamic models including linear doublet lattice, Euler solutions, Reynolds averaged Navier-Stokes and Delayed Detached Eddy Simulations. The analyses are performed across the transonic Mach range for a range in angle of attack from 0° to 5°, with focus on identifying the transonic dips in the flutter onset boundaries.
A method is presented for decoupling modal amplitudes generated by recursive least squares lattice filters. The proposed method, which involves transformation of basis functions obtained from lattice filters, is illustrated using a free-free beam simulation. The decoupled mode shape functions are required for the implementation of an adaptive control scheme for large space systems.
The effect of binding energy on the growth of clusters on a two-dimensional square lattice was studied by a means of a Monte Carlo simulation. The results indicate that for binding energies higher than 0.5kT, clusters grow relative to the initial random distribution, while for binding energies lower than 0.5 kT, clusters decompose. It was also found that at low binding energies the equilibrium percentage of atoms that are bound in clusters is linearly dependent on the binding energy. The cluster density behaved in a similar manner. This similarity reflects the fact that at low binding energies the growth proceeds mainly by the formation of new clusters rather than by an increase in the size of existing ones.
Full-scale simulations of a Gulfstream G-III aircraft, performed in support of the NASA Acoustic Research Measurements flights, are presented to complement results discussed in earlier studies. The flow solver employed in those studies, Dassault Systèmes’ lattice Boltzmann PowerFLOW®, was also used during this investigation to conduct time-dependent simulations of the entire aircraft in landing configuration with a fully dressed landing gear. The high-fidelity simulations, performed at a Mach number of 0.23 and a Reynolds number of 10.5 × 106 based on mean aerodynamic chord, captured all relevant airframe noise sources. The computations were used to assess the aeroacoustic performance of the main landing gear, with and without noise reduction fairings installed, of a G-III aircraft equipped with Adaptive Compliant Trailing Edge technology and conventional Fowler flaps. To facilitate comparison of predicted noise signatures with effective perceived noise levels obtained from flight test measurements, the “as-flown” nose landing gear geometry, missing in our earlier studies, was added to the simulated G-III aircraft configurations. The high fidelity, synthetic data were post-processed using a Ffowcs-Williams and Hawkings integral approach to estimate farfield acoustic behavior, with pressures on the model solid surface or pressure and velocity components on a permeable surface enveloping the acoustic near field used as input.
Pyrolitic boron nitride was exposed to 310 equivalent sun hours of ultraviolet radiation in a space environment simulator with the objective to evaluate its applicability as a pigment for a thermal control coating and to identify radiation damage using X-ray diffraction techniques. Lattice parameter comparisons show a definite increase in lattice imperfections in the crystal structure resulting from the ultraviolet irradiation. This sensitivity to radiation damage makes pyrolitic boron nitride unsuitable as a pigment for thermal control coating.
We present an architecture of distributed, independent control agents designed to work with the Computer Aided System Engineering and Analysis (CASE/A) simulation tool. CASE/A simulates behavior of Environmental Control and Life Support Systems (ECLSS). We describe a lattice of agents capable of distributed sensing and overcoming certain sensor and effector failures. We address how the architecture can achieve the coordinating functions of a hierarchical command structure while maintaining the robustness and flexibility of independent agents. These agents work between the time steps of the CASE/A simulation tool to arrive at command decisions based on the state variables maintained by CASE/A. Control is evaluated according to both effectiveness (e.g., how well temperature was maintained) and resource utilization (the amount of power and materials used).
A new approach to the prediction of helicopter rotor performance using a free wake analysis was developed. This new method uses a relaxation process that does not suffer from the convergence problems associated with previous time marching simulations. This wake relaxation procedure was coupled to a vortex-lattice, lifting surface loads analysis to produce a novel, self contained performance prediction code: EHPIC (Evaluation of Helicopter Performance using Influence Coefficients). The major technical features of the EHPIC code are described and a substantial amount of background information on the capabilities and proper operation of the code is supplied. Sample problems were undertaken to demonstrate the robustness and flexibility of the basic approach. Also, a performance correlation study was carried out to establish the breadth of applicability of the code, with very favorable results.