Search NASA⌕ Search

SEARCH · Search NASA

Results for “lattice simulations”

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 145 records · Page 8

An analytical study of the effects of jets located more than one jet diameter above a wing at subsonic speeds

A procedure has been developed for calculating the effects of blowing two jets over a swept tapered wing at low subsonic speeds. The algorithm used is based on a vortex-lattice representation of the wing lifting surface and a line sink-source distribution to simulate the effects of the jet exhaust on the wing lift and drag. The method is limited to those cases in which the jet exhaust does not intersect or wash the wing. The predictions of this relatively simple procedure are in remarkably good agreement with experimentally measured interference lift and interference induced drag.

Putnam, L. E.↗

Lateral Migration and Rotational Motion of Elliptic Particles in Planar Poiseuille Flow

Simulations of elliptic particulate suspensions in the planar Poiseuille flow are performed by using the lattice Boltzmann equation. Effects of the multi-particle on the lateral migration and rotational motion of both neutrally and non-neutrally buoyant elliptic particles are investigated. Low and intermediate total particle volume fraction f(sub a) = 13%, 15%, and 40% are considered in this work.

Qi, Dewei↗

Monte Carlo Simulation of Endlinking Oligomers

This report describes initial efforts to model the endlinking reaction of phenylethynyl-terminated oligomers. Several different molecular weights were simulated using the Bond Fluctuation Monte Carlo technique on a 20 x 20 x 20 unit lattice with periodic boundary conditions. After a monodisperse "melt" was equilibrated, chain ends were linked whenever they came within the allowed bond distance. Ends remained reactive throughout, so that multiple links were permitted. Even under these very liberal crosslinking assumptions, geometrical factors limited the degree of crosslinking. Average crosslink functionalities were 2.3 to 2.6; surprisingly, they did not depend strongly on the chain length. These results agreed well with the degrees of crosslinking inferred from experiment in a cured phenylethynyl-terminated polyimide oligomer.

Hinkley, Jeffrey A.↗

Aeroelastic Modeling and CFD Simulation of Wind-Tunnel Scale Aspect Ratio 13.5 Common Research Model

This paper presents an aeroelastic modeling and simulation study of an aspect ratio 13.5 wind-tunnel scale Common Research Model (CRM) with distributed flaps. A vortex-lattice VSPAERO model of the CRM model is developed. A transonic small disturbance/integral boundary layer correction method is implemented in the VSPAERO model to account for the transonic and viscous flow effects. The structural deformation of the CRM model is calculated using a NASTRAN equivalent beam model. The VSPAERO model is coupled to the NASTRAN equivalent beam model to provide a rapid aero-structural analysis. A validation of the VSPAERO aeroelastic model is conducted by comparing the results to FUN3D CFD aeroelastic simulation results. An aerodynamic database is generated using the developed VSPAERO aeroelastic model for the real-time drag optimization and maneuver load alleviation study of the wind-tunnel scale CRM model.

Aeroelastic Modeling↗

Simulation studies for surfaces and materials strength

From intermolecular force studies, it is now known that the overall non-additive contribution to the lattice enegy is positive so that analysis based on only pairwise additivity suggests a shallower intermolecular potential than the true value. Two body contributions alone are also known to be categorically unable to even qualitatively describe some configurations of molecular clusters in the gas phase or the general relaxation and reconstruction of fcc crystal surfaces. In addition, the many-body contribution was shown to play a key role in the stability of certain crystal structures. In these recent analyses, a relatively simple potential energy function (PEF), comprising only a two-body Mie-type potential plus a three-body Axilrod-Teller-type potential, was found to be extremely effective. This same parametric PEF is applied to describe the bulk stability and surface energy for the diamond cubic structure. To test the stability condition, the FCC, BCC, diamond cubic, graphite and beta-tin structures were considered.

Halicioglu, T.↗

A survey of the state-of-the-art and focused research in range systems

In this one-year renewal of NASA Contract No. 2-304, basic research, development, and implementation in the areas of modern estimation algorithms and digital communication systems have been performed. In the first area, basic study on the conversion of general classes of practical signal processing algorithms into systolic array algorithms is considered, producing four publications. Also studied were the finite word length effects and convergence rates of lattice algorithms, producing two publications. In the second area of study, the use of efficient importance sampling simulation technique for the evaluation of digital communication system performances were studied, producing two publications.

Kung, Yao↗

Studies in turbulence

Various papers on turbulence are presented. Individual topics addressed include: modeling the dissipation rate in rotating turbulent flows, mapping closures for turbulent mixing and reaction, understanding turbulence in vortex dynamics, models for the structure and dynamics of near-wall turbulence, complexity of turbulence near a wall, proper orthogonal decomposition, propagating structures in wall-bounded turbulence flows. Also discussed are: constitutive relation in compressible turbulence, compressible turbulence and shock waves, direct simulation of compressible turbulence in a shear flow, structural genesis in wall-bounded turbulence flows, vortex lattice structure of turbulent shear slows, etiology of shear layer vortices, trilinear coordinates in fluid mechanics.

Gatski, Thomas B.↗

Modeling of the Site Preference in Ternary B2-Ordered Ni-Al-Fe Alloys

The underlying equilibrium structure, site substitution behavior, and lattice parameter of ternary Ni-Fe-Al alloys are determined via Monte Carlo-Metropolis computer simulations and analytical calculations using the BFS method for alloys for the energetics. As a result of the theoretical calculations presented, a simple approach based on the energetics of small atomic clusters is introduced to explain the observed site preference schemes.

Bozzolo, Guillermo H.↗

Lattice Boltzmann and Navier-Stokes Cartesian CFD Approaches for Airframe Noise Predictions

Lattice Boltzmann (LB) and compressible Navier-Stokes (NS) equations based computational fluid dynamics (CFD) approaches are compared for simulating airframe noise. Both LB and NS CFD approaches are implemented within the Launch Ascent and Vehicle Aerodynamics (LAVA) framework. Both schemes utilize the same underlying Cartesian structured mesh paradigm with provision for local adaptive grid refinement and sub-cycling in time. We choose a prototypical massively separated, wake-dominated flow ideally suited for Cartesian-grid based approaches in this study - The partially-dressed, cavity-closed nose landing gear (PDCC-NLG) noise problem from AIAA's Benchmark problems for Airframe Noise Computations (BANC) series of workshops. The relative accuracy and computational efficiency of the two approaches are systematically compared. Detailed comments are made on the potential held by LB to significantly reduce time-to-solution for a desired level of accuracy within the context of modeling airframes noise from first principles.

Airframe Noise Prediction↗

The evaluation of the rolling moments induced by wraparound fins

A possible reason is suggested for the induced rolling moments occurring on wraparound-fin configurations in subsonic flight at zero angle of attack. The subsonic potential flow over the configuration at zero incidence is solved numerically. The body is simulated by a distribution of sources along its axis, and the fins are described by a vortex-lattice method. It is shown that rolling moments can be induced on the antisymmetric fins by the radial flow generated at the base of the configuration, either over the converging separated wake, or over the diverging plume of a rocket motor.

Seginer, A.↗

The evaluation of the rolling moments induced by wraparound fins

A possible reason is suggested for the induced rolling moments occurring on wraparound-fin configurations in subsonic flight at zero angle of attack. The subsonic potential flow over the configuration at zero incidence is solved numerically. The body is simulated by a distribution of sources along its axis, and the fins are described by a vortex-lattice method. It is shown that rolling moments can be induced on the antisymmetric fins by the radial flow generated at the base of the configuration, either over the converging separated wake, or over the diverging plume of rocket motor.

Seginer, A.↗

Computation of rotor aerodynamic loads with a constant vorticity contour free wake model

An analytical method is presented which facilitates the study of isolated rotors with an improved approach to wake simulation. Vortex filaments are simulated along contours of constant sheet strength for the sheet of vorticity resulting from each rotor blade. Curved vortex elements comprise the filaments which can be distorted by the local velocity field. Called the Constant Vorticity Contour wake model, the approach permits the simulation of the blades' wakes corresponding to the full span of the rotor blade. The discretization of the wake of the rotor blade produces spacing and structure that are consistent with the spatial and temporal variations in the loading. A vortex-lattice aerodynamic model of the blade is also included which introduces a finite-element structural model of the blade and consideration of the force and moment trim analysis. Results of the present version of the simulation, called RotorCRAFT, are found to correlate well with H-34 flight-test data.

Quackenbush, Todd R.↗

Self-Consistent Simulation of the Brownian Stage of Dust Growth

It is a widely accepted view that in proto-planetary accretion disks the collision and following sticking of dust particles embedded in the gas eventually leads to the formation of planetesimals (coagulation). For the smallest dust grains, Brownian motion is assumed to be the dominant source of their relative velocities leading to collisions between these dust grains. As the dust grains grow they eventually couple to the turbulent motion of the gas which then drives the coagulation much more efficiently. Many numerical coagulation simulations have been carried out to calculate the fractal dimension of the aggregates, which determines the duration of the ineffective Brownian stage of growth. Predominantly on-lattice and off-lattice methods were used. However, both methods require simplification of the astrophysical conditions. The aggregates found by those methods had a fractal dimension of approximately 2 which is equivalent to a constant, mass-independent friction time. If this value were valid for the conditions in an accretion disk, this would mean that the coagulation process would finally 'freeze out' and the growth of a planetesimal would be impossible within the lifetime of an accretion disk. In order to investigate whether this fractal dimension is model independent, we simulate self-consistently the Brownian stage of the coagulation by an N-particle code. This method has the advantage that no further assumptions about homogeneity of the dust have to be made. In our model, the dust grains are considered as aggregates built up of spheres. The equation of motion of the dust grains is based on the probability density for the diffusive transport within the gas atmosphere. Because of the very low number density of the dust grains, only 2-body-collisions have to be considered. As the Brownian stage of growth is very inefficient, the system is to be simulated over long periods of time. In order to find close particle pairs of the system which are most likely to undergo a collision, we use a particle-in-cell (PIC) method for the early stages of the simulation where the system is still very homogeneous and a tree method later when the particles are more clustered.

Kempf, S.↗

System Identification of a Vortex Lattice Aerodynamic Model

The state-space presentation of an aerodynamic vortex model is considered from a classical and system identification perspective. Using an aerodynamic vortex model as a numerical simulator of a wing tunnel experiment, both full state and limited state data or measurements are considered. Two possible approaches for system identification are presented and modal controllability and observability are also considered. The theory then is applied to the system identification of a flow over an aerodynamic delta wing and typical results are presented.

Juang, Jer-Nan↗

Kinetic Monte Carlo Simulations of Rod Eutectics and the Surface Roughening Transition in Binary Alloys

In this paper we report on two related topics: Kinetic Monte Carlo simulations of the steady state growth of rod eutectics from the melt, and a study of the surface roughness of binary alloys. We have implemented a three dimensional kinetic Monte Carlo (kMC) simulation with diffusion by pair exchange only in the liquid phase. Entropies of fusion are first chosen to fit the surface roughness of the pure materials, and the bond energies are derived from the equilibrium phase diagram, by treating the solid and liquid as regular and ideal solutions respectively. A simple cubic lattice oriented in the {100} direction is used. Growth of the rods is initiated from columns of pure B material embedded in an A matrix, arranged in a close packed array with semi-periodic boundary conditions. The simulation cells typically have dimensions of 50 by 87 by 200 unit cells. Steady state growth is compliant with the Jackson-Hunt model. In the kMC simulations, using the spin-one Ising model, growth of each phase is faceted or nonfaceted phases depending on the entropy of fusion. There have been many studies of the surface roughening transition in single component systems, but none for binary alloy systems. The location of the surface roughening transition for the phases of a eutectic alloy determines whether the eutectic morphology will be regular or irregular. We have conducted a study of surface roughness on the spin-one Ising Model with diffusion using kMC. The surface roughness was found to scale with the melting temperature of the alloy as given by the liquidus line on the equilibrium phase diagram. The density of missing lateral bonds at the surface was used as a measure of surface roughness.

Bentz, Daniel N.↗

Inhomogeneous diffusion-limited aggregation

It is demonstrated here that inhomogeneous diffusion-limited aggregation (DLA) model can be used to simulate viscous fingering in a medium with inhomogeneous permeability and homogeneous porosity. The medium consists of a pipe-pore square-lattice network in which all pores have equal volume and the pipes have negligible volume. It is shown that fluctuations in a DLA-based growth process may be tuned by noise reduction, and that fluctuations in the velocity of the moving interface are multiplicative in form.

Selinger, Robin Blumberg↗

Synthesis and Stability of Iron Nanoparticles for Lunar Environment Studies

Simulant of lunar dust is needed when researching the lunar environment. However, unlike the true lunar dust, today s simulants do not contain nanophase iron. Two different processes have been developed to fabricate nanophase iron to be used as part of the lunar dust simulant: (1) Sequentially treating a mixture of ferric chloride, fluorinated carbon, and soda lime glass beads at about 300 C in nitrogen, at room temperature in air, and then at 1050 C in nitrogen. The product includes glass beads that are grey in color, can be attracted by a magnet, and contain alpha-iron nanoparticles (which seem to slowly lose their lattice structure in ambient air during a period of 12 months). This product may have some similarity to the lunar glassy regolith that contains Fe(sup 0). (2) Heating a mixture of carbon black and a lunar simulant (a mixed metal oxide that includes iron oxide) at 1050 C in nitrogen. This process simulates lunar dust reaction to the carbon in a micrometeorite at the time of impact. The product contains a chemically modified simulant that can be attracted by a magnet and has a surface layer whose iron concentration increased during the reaction. The iron was found to be alpha-iron and Fe3O4 nanoparticles, which appear to grow after the fabrication process, but stabilizes after 6 months of ambient air storage.

Hung, Ching-cheh↗

Subsonic steady and unsteady aerodynamic loads on missiles and aircraft

Steady lifting flows over highly swept delta wings at large incidence were studied. After an exhaustive literature review, development of a vortex-lattice method was attempted. To demonstrate the feasibility of using such a method, an existing code was modified. A system of vortex lines to simulate the leading-edge wake was added. The coefficients predicted by the modified code were in good agreement with experimental data.

Source record↗