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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

251 records · Page 14

All sources NASA NTRS DOE OSTI

NASA NTRS1994

A distributed analysis and visualization system for model and observational data

Software was developed with NASA support to aid in the analysis and display of the massive amounts of data generated from satellites, observational field programs, and from model simulations. This software was developed in the context of the PATHFINDER (Probing ATmospHeric Flows in an Interactive and Distributed EnviRonment) Project. The overall aim of this project is to create a flexible, modular, and distributed environment for data handling, modeling simulations, data analysis, and visualization of atmospheric and fluid flows. Software completed with NASA support includes GEMPAK analysis, data handling, and display modules for which collaborators at NASA had primary responsibility, and prototype software modules for three-dimensional interactive and distributed control and display as well as data handling, for which NSCA was responsible. Overall process control was handled through a scientific and visualization application builder from Silicon Graphics known as the Iris Explorer. In addition, the GEMPAK related work (GEMVIS) was also ported to the Advanced Visualization System (AVS) application builder. Many modules were developed to enhance those already available in Iris Explorer including HDF file support, improved visualization and display, simple lattice math, and the handling of metadata through development of a new grid datatype. Complete source and runtime binaries along with on-line documentation is available via the World Wide Web at: http://redrock.ncsa.uiuc.edu/ PATHFINDER/pathre12/top/top.html.

Wilhelmson, Robert B.↗

NASA NTRS2005

A Comparative Study Using CFD to Predict Iced Airfoil Aerodynamics

WIND, Fluent, and PowerFLOW were used to predict the lift, drag, and moment coefficients of a business-jet airfoil with a rime ice (rough and jagged, but no protruding horns) and with a glaze ice (rough and jagged end has two or more protruding horns) for angles of attack from zero to and after stall. The performance of the following turbulence models were examined by comparing predictions with available experimental data. Spalart-Allmaras (S-A), RNG k-epsilon, shear-stress transport, v(sup 2)-f, and a differential Reynolds stress model with and without non-equilibrium wall functions. For steady RANS simulations, WIND and FLUENT were found to give nearly identical results if the grid about the iced airfoil, the turbulence model, and the order of accuracy of the numerical schemes used are the same. The use of wall functions was found to be acceptable for the rime ice configuration and the flow conditions examined. For rime ice, the S-A model was found to predict accurately until near the stall angle. For glaze ice, the CFD predictions were much less satisfactory for all turbulence models and codes investigated because of the large separated region produced by the horns. For unsteady RANS, WIND and FLUENT did not provide better results. PowerFLOW, based on the Lattice Boltzmann method, gave excellent results for the lift coefficient at and near stall for the rime ice, where the flow is inherently unsteady.

Chi, x.↗

NASA NTRS1980

An experimental investigation of three dimensional low speed minimum interference wind tunnel for high lift wings

As a means to achieve a minimum interference correction wind tunnel, a partially actively controlled test section was experimentally examined. A jet flapped wing with 0.91 m (36 in) span and R = 4.05 was used as a model to create moderately high lift coefficients. The partially controlled test section was simulated using an insert, a rectangular box 0.96 x 1.44 m (3.14 x 4.71 ft) open on both ends in the direction of the tunnel air flow, placed in the University of Washington Aeronautical Laboratories (UWAL) 2.44 x 3.66 m (8 x 12 ft) wind tunnel. A tail located three chords behind the wing was used to measure the downwash at the tail region. The experimental data indicates that, within the range of momentum coefficient examined, it appears to be unnecessary to actively control all four sides of the test section walls in order to achieve the near interference free flow field environment in a small wind tunnel. The remaining wall interference can be satisfactorily corrected by the vortex lattice method.

Shindo, S.↗

NASA NTRS1983

Adaptive control of a flexible beam using least square lattice filters

This paper presents an indirect adaptive control scheme for the control of flexible structures using recursive least square lattice filters. The identification scheme uses lattice filters which provide an on-line estimate of the number of modes, mode shapes and modal amplitudes. These modes are coupled and a transformation to decouple them in order to obtain the natural modes is presented. The decoupled modal amplitude time series are then used in an equation error identification scheme to identify the model parameters in an autoregressive moving average (ARMA) form. The control is based on modal pole placement scheme with the objective of vibration suppression. The control gains are calculated based on the identified ARMA parameters. Before using the identified parameters for control, detailed testing and validation procedures are carried out on the identified parameters. The full adaptive control scheme is demonstrated using the simulation for the 12 foot free-free beam apparatus at NASA Langley Research Center.

Sundararajan, N.↗

NASA NTRS1985

Numerical simulation of comet nuclei. I - Water-ice comets

A one-dimensional numerical model of pure water-ice cometary nuclei is presented, and the influence of the nuclear interior as a heat reservoir on the behavior of the nuclear surface is examined. It is shown that a number of effects, including the thermal inertia due to heat stored in the core and the release of latent heat, which goes entirely into heating the adjacent layers or into sublimation on passing through a phase transition from amorphous to crystalline ice, can help to explain such characteristics as the asymmetrical lightcurve of Comet Halley. Results are given for the cases of Comet Schwassmann-Wachmann 1 and Comet Encke. Consideration is also given to the insulating effect of an evolving dust mantle. The role of this mantle in determining the surface temperature of the ice core is studied as a function of the mass fraction of the dust in the ice-dust mixture and the thermal conductivity of the nucleus. The loose-lattice model of Mendis and Brin (1977) indicates that both high dust to ice ratios and high-core conductivities inhibit mantle blowoff.

Herman, G.↗

NASA NTRS1994

Equilibrium and Dynamics Properties of Poly(oxyethylene) Melts and Related Poly(alkylethers) from Simulations and Ab Initio Calculations

Molecular dynamics simulations of POE melts have been performed utilizing a potential force field parameterized to reproduce conformer energies and rotational energy barriers in dimethoxyethane as determined from ab initio electronic structure calculations. Chain conformations and dimensions of POE from the simulations were found to be in good agreement with predictions of a rotational isomeric state (RIS) model based upon the ab initio conformational. energies. The melt chains were found to be somewhat extended relative to chains at theta conditions. This effect will be discussed in light of neutron scattering experiments which indicate that POE chains are extended in the melt relative to theta solutions. The conformational characteristics of POE chains will also be compared with those of other poly (alkylethers), namely poly(oxymethylene), poly(oxytrimethylene) and poly(oxytetramethylene). Local conformational dynamics were found to be more rapid than in polymethylene. Calculated C-H vector correlation times were found to be in reasonable agreement with experimental values from C-13 NMR spin-lattice relaxation times. The influence of ionic salts on local conformations and dynamics will also be discussed.

Smith, Grant D.↗

NASA NTRSDate not supplied

Multi-point Jig Twist Optimization of Mach 0.745 Transonic Truss-Braced Wing Aircraft and High-Fidelity CFD Validation

This paper presents a multi-point jig twist optimization study of the Mach 0.745 Transonic Truss-Braced Wing (TTBW) aircraft using an in-house developed aero-structural analysis solver VSPAERO coupled to BEAM3D. A vortex-lattice model of the TTBW model is developed, and a transonic small disturbance /integral boundary layer correction method is implemented in the VSPAERO model to account for transonic and viscous flow effects. A correction method for the wing-strut interference aerodynamics is developed and applied to the VSPAERO solver. Also, a structural dynamic finite-element model of the TTBW aircraft is developed. This finite-element model includes the geometric nonlinear effect due to the tension in the struts which causes a deflection-dependent nonlinear stiffness. The VSPAERO model is coupled to the corresponding finite-element model to provide a rapid aero-structural analysis. The VSPAERO model has been validated against wind tunnel test data and FUN3D CFD simulation data. Three different flight conditions corresponding to Mach 0.745 at 40,000 ft are selected for the TTBW aircraft jig twist optimization to reduce the drag coefficients at the design and off-design cruise lift coefficients. The preliminary optimization results show that the TTBW aircraft with the optimized jig twist achieves a drag reduction of about 9 counts, 7 counts, and8 counts corresponding to the lift coefficients 0.68, 0.73, and 0.78, respectively. A high-fidelity CFD solver FUN3D is used to validate the jig twist optimization.

Multi-point↗

NASA NTRS2010

A Geometric Approach to Modeling Microstructurally Small Fatigue Crack Formation: Simulation and Prediction of Crack Nucleation in AA 7075-T651 - 2

The objective of this paper is to develop further a framework for computationally modeling microstructurally small fatigue crack growth in AA 7075-T651 [1]. The focus is on the nucleation event, when a crack extends from within a second-phase particle into a surrounding grain, since this has been observed to be an initiating mechanism for fatigue crack growth in this alloy. It is hypothesized that nucleation can be predicted by computing a non-local nucleation metric near the crack front. The hypothesis is tested by employing a combination of experimentation and nite element modeling in which various slip-based and energy-based nucleation metrics are tested for validity, where each metric is derived from a continuum crystal plasticity formulation. To investigate each metric, a non-local procedure is developed for the calculation of nucleation metrics in the neighborhood of a crack front. Initially, an idealized baseline model consisting of a single grain containing a semi-ellipsoidal surface particle is studied to investigate the dependence of each nucleation metric on lattice orientation, number of load cycles, and non-local regularization method. This is followed by a comparison of experimental observations and computational results for microstructural models constructed by replicating the observed microstructural geometry near second-phase particles in fatigue specimens. It is found that orientation strongly influences the direction of slip localization and, as a result, in uences the nucleation mechanism. Also, the baseline models, replication models, and past experimental observation consistently suggest that a set of particular grain orientations is most likely to nucleate fatigue cracks. It is found that a continuum crystal plasticity model and a non-local nucleation metric can be used to predict the nucleation event in AA 7075-T651. However, nucleation metric threshold values that correspond to various nucleation governing mechanisms must be calibrated.

Hochhalter, Jake D.↗

NASA NTRS1993

Multidisciplinary aeroelastic analysis of a generic hypersonic vehicle

This paper presents details of a flutter and stability analysis of aerospace structures such as hypersonic vehicles. Both structural and aerodynamic domains are discretized by the common finite element technique. A vibration analysis is first performed by the STARS code employing a block Lanczos solution scheme. This is followed by the generation of a linear aerodynamic grid for subsequent linear flutter analysis within subsonic and supersonic regimes of the flight envelope; the doublet lattice and constant pressure techniques are employed to generate the unsteady aerodynamic forces. Flutter analysis is then performed for several representative flight points. The nonlinear flutter solution is effected by first implementing a CFD solution of the entire vehicle. Thus, a 3-D unstructured grid for the entire flow domain is generated by a moving front technique. A finite element Euler solution is then implemented employing a quasi-implicit as well as an explicit solution scheme. A novel multidisciplinary analysis is next effected that employs modal and aerodynamic data to yield aerodynamic damping characteristics. Such analyses are performed for a number of flight points to yield a large set of pertinent data that define flight flutter characteristics of the vehicle. This paper outlines the finite-element-based integrated analysis procedures in detail, which is followed by the results of numerical analyses of flight flutter simulation.

Gupta, K. K.↗

NASA NTRS2005

MOVPE of GaSb/InGaAsSb Multilayers and Fabrication of Dual Band Photodetectors

Metalorganic vapor phase epitaxy (MOVPE) of GaSb/InGaAsSb multilayer thin films and fabrication of bias-selectable dual band photodetectors are reported. For the dual band photodetectors the short wavelength detector, or the upper p- GaSb/n-GaSb junction photodiode, is placed optically ahead of the long wavelength one, or the lower photodiode. The latter is based on latticed-matched In0.13Ga0.87As0.11Sb0.89 with bandgap near 0.6 eV. Specifically, high quality multilayer thin films are grown sequentially from top to bottom as p+-GaSb/p-GaSb/n-GaSb/n-InGaAsSb/p-InGaAsSb/p-GaSb on undoped p-type GaSb substrate, and as n-GaSb/p-GaSb/p-InGaAsSb/n-InGaAsSb/n-GaSb on Te-doped n-type GaSb substrate respectively. The multilayer thin films are characterized by optical microscope, atomic force microscope (AFM), electron microprobe analyses etc. The photodiode mesa steps are patterned by photolithography with wet chemical etching and the front metallization is carried out by e-beam evaporation with Pd/Ge/Au/Ti/Au to give ohmic contact on both n- and p-type Sb based layer surfaces. Dark I-V measurements show typical diode behavior for both the upper and lower photodiodes. The photoresponsivity measurements indicate that both the upper and lower photodiodes can sense the infrared illumination corresponding to their cutoff wavelengths respectively, comparable with the simulation results. More work is underway to bring the long wavelength band to the medium infrared wavelength region near 4 micrometers.

Xiao, Ye-Gao↗

NASA NTRS2023

HLPW-4: Wall-Modeled LES and Lattice-Boltzmann Technology Focus Group Workshop Summary

TTT↗

NASA NTRS2010

Thermal Properties of Lunar Regolith Simulants

Various high temperature chemical processes have been developed to extract oxygen and metals from lunar regolith. These processes are tested using terrestrial analogues of the regolith. But all practical terrestrial analogs contain H2O and/or OH(-), the presence of which has substantial impact on important system behaviors. We have undertaken studies of lunar regolith simulants to determine the limits of the simulants to validate key components for human survivability during sustained presence on the moon. Differential Thermal Analysis (DTA) yields information on phase transitions and melting temperatures. Themo-Gravimetric Analysis (TGA) with mass spectrometric (MS) determination of evolved gas species yields chemical information on various oxygenated volatiles (water, carbon dioxide, sulfur oxides, nitrogen oxides and phosphorus oxides) and their evolution temperature profiles. The DTA and TGAMS studies included JSC-1A fine, NU-LHT-2M and its proposed feed stocks: anorthosite; dunite; HQ (high quality) glass and the norite from which HQ glass is produced. Fig 1 is a data profile for anorthosite. The DTA (Fig 1a) indicates exothermic transitions at 355 and 490 C and endothermic transitions at 970 and 1235 C. Below the 355 C transition, water (Molecular Weight, MW, 18 in Fig 1c) is lost accounting for approximately 0.1% mass loss due to water removal (Fig 1b). Just above 490 C a second type of water is lost, presumably bound in lattices of secondary minerals. Between 490 and the 970 transition other volatile oxides are lost including those of hydrogen (third water type), carbon (MW = 44), sulfur (MW = 64 and 80), nitrogen (MW 30 and 46) and possibly phosphorus (MW = 79, 95 or 142). Peaks at MW = 35 and 19 may be attributable to loss of chlorine and fluorine respectively. Negative peaks in the NO (MW = 30) and oxygen (MW = 32) MS profiles may indicate the production of NO2 (MW = 46). Because so many compounds are volatilized in this temperature range quantification of the mass loss associated with individual species is difficult. Similar information will be presented for the other materials studied in this investigation.

Street, Kenneth↗

NASA NTRS2010

Thermal Properties of Lunar Regolith Simulants

Various high temperature chemical processes have been developed to extract oxygen and metals from lunar regolith. These processes are tested using terrestrial analogues of the regolith. But all practical terrestrial analogs contain H2O and/or OH-, the presence of which has substantial impact on important system behaviors. We have undertaken studies of lunar regolith simulants to determine the limits of the simulants to validate key components for human survivability during sustained presence on the moon. Differential Thermal Analysis (DTA) yields information on phase transitions and melting temperatures. Themo-Gravimetric Analysis (TGA) with mass spectrometric (MS) determination of evolved gas species yields chemical information on various oxygenated volatiles (water, carbon dioxide, sulfur oxides, nitrogen oxides and phosphorus oxides) and their evolution temperature profiles. The DTA and TGAMS studies included JSC-1A fine, NU-LHT-2M and its proposed feed stocks: anorthosite; dunite; HQ (high quality) glass and the norite from which HQ glass is produced. Fig 1 is a data profile for anorthosite. The DTA (Fig 1a) indicates exothermic transitions at 355 and 490 C and endothermic transitions at 970 and 1235 C. Below the 355 C transition, water (Molecular Weight, MW, 18 in Fig 1c) is lost accounting for approximately 0.1% mass loss due to water removal (Fig 1b). Just above 490 C a second type of water is lost, presumably bound in lattices of secondary minerals. Between 490 and the 970 transition other volatile oxides are lost including those of hydrogen (third water type), carbon (MW = 44), sulfur (MW = 64 and 80), nitrogen (MW 30 and 46) and possibly phosphorus (MW = 79, 95 or 142). Peaks at MW = 35 and 19 may be attributable to loss of chlorine and fluorine respectively. Negative peaks in the NO (MW = 30) and oxygen (MW = 32) MS profiles may indicate the production of NO2 (MW = 46). Because so many compounds are volatilized in this temperature range quantification of the mass loss associated with individual species is difficult. Similar information will be presented for the other materials studied in this investigation.

Street, Kenneth↗

NASA NTRS1989

The formation of topological defects in phase transitions

It was argued, and fought through numerical work that the results of non-dynamical Monte Carlo computer simulations cannot be applied to describe the formation of topological defects when the correlation length at the Ginzburg temperature is significantly smaller than the horizon size. To test the current hypothesis that infinite strings at formation are essentially described by Brownian walks of size the correlation length at the Ginzburg temperature, fields at the Ginzburg temperature were equilibrated. Infinite structure do not exist in equilibrium for reasonable definitions of the Ginzburg temperature, and horizons must be included in a proper treatment. A phase transition, from small-scale to large-scale string or domain wall structure, is found to occur very close to the Ginzburg temperature, in agreement with recent work. The formation process of domain walls and global strings were investigated through the breaking of initially ordered states. To mimic conditions in the early Universe, cooling times are chosen so that horizons exist in the sample volume when topological structure formation occurs. The classical fields are evolved in real-time by the numerical solution of Langevin equations of motion on a three dimensional spatial lattice. The results indicate that it is possible for most of the string energy to be in small loops, rather than in long strings, at formation.

Hodges, Hardy M.↗

NASA NTRS1998

Nuclear modules for space electric propulsion

Analysis of interplanetary cargo and piloted missions requires calculations of the performances and masses of subsystems to be integrated in a final design. In a preliminary and scoping stage the designer needs to evaluate options iteratively by using fast computer simulations. The Oak Ridge National Laboratory (ORNL) has been involved in the development of models and calculational procedures for the analysis (neutronic and thermal hydraulic) of power sources for nuclear electric propulsion. The nuclear modules will be integrated into the whole simulation of the nuclear electric propulsion system. The vehicles use either a Brayton direct-conversion cycle, using the heated helium from a NERVA-type reactor, or a potassium Rankine cycle, with the working fluid heated on the secondary side of a heat exchanger and lithium on the primary side coming from a fast reactor. Given a set of input conditions, the codes calculate composition. dimensions, volumes, and masses of the core, reflector, control system, pressure vessel, neutron and gamma shields, as well as the thermal hydraulic conditions of the coolant, clad and fuel. Input conditions are power, core life, pressure and temperature of the coolant at the inlet of the core, either the temperature of the coolant at the outlet of the core or the coolant mass flow and the fluences and integrated doses at the cargo area. Using state-of-the-art neutron cross sections and transport codes, a database was created for the neutronic performance of both reactor designs. The free parameters of the models are the moderator/fuel mass ratio for the NERVA reactor and the enrichment and the pitch of the lattice for the fast reactor. Reactivity and energy balance equations are simultaneously solved to find the reactor design. Thermalhydraulic conditions are calculated by solving the one-dimensional versions of the equations of conservation of mass, energy, and momentum with compressible flow.

Difilippo, F. C.↗

NASA NTRS2004

Progressive Response Surfaces

Response surface functions are often used as simple and inexpensive replacements for computationally expensive computer models that simulate the behavior of a complex system over some parameter space. Progressive response surfaces are ones that are built up progressively as global information is added from new sample points in the parameter space. As the response surfaces are globally upgraded based on new information, heuristic indications of the convergence of the response surface approximation to the exact (fitted) function can be inferred. Sampling points can be incrementally added in a structured fashion, or in an unstructured fashion. Whatever the approach, at least in early stages of sampling it is usually desirable to sample the entire parameter space uniformly. At later stages of sampling, depending on the nature of the quantity being resolved, it may be desirable to continue sampling uniformly over the entire parameter space (Progressive response surfaces), or to switch to a focusing/economizing strategy of preferentially sampling certain regions of the parameter space based on information gained in early stages of sampling (Adaptive response surfaces). Here we consider Progressive response surfaces where a balanced indication of global response over the parameter space is desired.We use a variant of Moving Least Squares to fit and interpolate structured and unstructured point sets over the parameter space. On a 2-D test problem we compare response surface accuracy for three incremental sampling methods: Progressive Lattice Sampling; Simple-Random Monte Carlo; and Halton Quasi-Monte-Carlo sequences. We are ultimately after a system for constructing efficiently upgradable response surface approximations with reliable error estimates.

Romero, V. J.↗

NASA NTRSDate not supplied

HLPW-4/GMGW-3: Wall-Modeled LES and Lattice-Boltzmann Technology Focus Group Workshop Summary

A summary of the nine submissions to the Wall-Modeled LES and Lattice-Boltzmann(WMLESLB) Technical Focus Group (TFG) at the 4th High lift Prediction Workshop is provided. The focus of this TFG was to assess the current capabilities of WMLES and Lattice-Boltzmann methods on a complex high-lift configuration across a wide range of angles of attack. Preliminary analysis of the submitted data suggests that>250𝑀spatial degrees of freedom are needed to accurately predict pitching moments at high angles of attack due to large pressure gradients present on the outboard slat and main element for𝛼 >17◦(corrected for free-air).While some scatter is reported in pitching moments at the low-angles of attack (𝛼 <11◦) - likely caused by differences in flap separation possibly due to low Reynolds number effects - excellent agreement is observed between the submissions near the𝐶𝐿,max state. Objective superiority over steady state RANS can be seen in terms of lack of excess outboard separation; a majority of the good quality WMLES and LB submissions predict wedge-shaped separation patterns consistent with the experimental oil flow. Differences in the onset of stall mechanism in the free-air configuration for𝛼 >20◦is reported with two distinct topologies observed. Topology A is characterized by the onset of corner-flow separation which progressively grows to produce a pitch break in free-air with an angle of attack increase from𝛼=20.55◦to𝛼=21.47◦+𝜖where𝜖is a small perturbation (varying between submissions and likely to be within±0.3◦). Topology Bis characterized by boundary layer weakness emanating from the inboard side of the wing-pylon juncture substantially larger than any weakness in the wing-body juncture (due to corner-flow);submissions in this category do not show any tendency for occurrence of a pitch break for the free-air configuration within the vicinity of𝛼=21.47◦. The in-tunnel simulations submitted by3 participants using different discretizations, grids and closure models show excellent agreement with the experiment in terms of a) integrated loads, b) surface flow-topology, and c) mechanism for the onset of inboard stall. Further evidence is provided to demonstrate both qualitative and quantitative superiority of all 3 WMLES submissions over a single steady state RANS submission to the workshop.

TTT↗