Lattice Boltzmann Method for Spacecraft Propellant Slosh Simulation
No abstract available
Powers, Joseph↗
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.
Growth morphology with anisotropic surface kinetics
The morphological evolution of crystals growing from an incongruent vapor phase is studied using a Monte Carlo model, and the full range of growth morphologies is recovered. The diffusion in the bulk nutrient and the anisotropy in the interface kinetics are morphologically destabilizing and stabilizing, respectively. For a given set of simulation parameters and lattice symmetries there is a critical size, which scales linearly with the mean free path in the vapor, beyond which a crystal cannot retain its stable, macroscopically faceted growth shape. Surface diffusion stabilizes faceted growth on the shorter scale of the mean surface diffusion length. In simulations with a uniform drift superimposed on the random walk nutrient transport, crystal faces oriented toward the drift show enhanced morphological stability compared to the purely diffusive situation. Rotational drifts with periodic reversal of direction are morphologically stabilizing for all crystal facets.
Xiao, Rong-Fu↗
Lattice Thermal Conductivity of Ultra High Temperature Ceramics (UHTC) ZrB2 and HfB2 from Atomistic Simulations
Ultra high temperature ceramics (UHTC) including ZrB2 and HfB2 are candidate materials for applications in extreme environments because of their high melting point, good mechanical properties and reasonable oxidation resistance. Unlike many ceramics, these materials have high thermal conductivity which can be advantageous, for example, to reduce thermal shock. Recently, we developed Tersoff style interatomic potentials for both ZrB2 and HfB2 appropriate for atomistic simulations. As an application, Green-Kubo molecular dynamics simulations were performed to evaluate the lattice thermal conductivity for single crystals of ZrB2 and HfB2. The atomic mass difference in these binary compounds leads to oscillations in the time correlation function of the heat current. Results at room temperature and at elevated temperatures will be reported.
Lawson, JOhn W.↗
Lattice Thermal Conductivity of Ultra High Temperature Ceramics (UHTC) ZrB2 and HfB2 from Atomistic Simulations
Ultra high temperature ceramics (UHTC) including ZrB2 and HfB2 have a number of properties that make them attractive for applications in extreme environments. One such property is their high thermal conductivity. Computational modeling of these materials will facilitate understanding of fundamental mechanisms, elucidate structure-property relationships, and ultimately accelerate the materials design cycle. Progress in computational modeling of UHTCs however has been limited in part due to the absence of suitable interatomic potentials. Recently, we developed Tersoff style parameterizations of such potentials for both ZrB2 and HfB2 appropriate for atomistic simulations. As an application, Green-Kubo molecular dynamics simulations were performed to evaluate the lattice thermal conductivity for single crystals of ZrB2 and HfB2. The atomic mass difference in these binary compounds leads to oscillations in the time correlation function of the heat current, in contrast to the more typical monotonic decay seen in monoatomic materials such as Silicon, for example. Results at room temperature and at elevated temperatures will be reported.
Lawson, John W.↗
Multijunction upper subcell cascade photovoltaics for space applications
A new class of cascade high-efficiency photovoltaics designed for space-based applications is proposed. The design improves upper subcell performance and avoids electrical and optical losses associated with an intercell ohmic contact. Multijunction upper subcells reduce bulk recombination of photogenerated minority carriers by decreasing the average collection distance, yielding improved spectral response and radiation tolerance. A three-terminal design is employed which circumvents the need for a monolithic intercell contact and, thus, the losses associated with such a contact. Problems related to array interconnection of three-terminal devices may be solved by creating a two-terminal cell from complementary pairs (n-p-n and p-n-p) of three-terminal cells. Simulations of lattice-matched AlGaAs-GaAs and lattice-mismatched AlGaAs-InGaAs cascade cells show that one-sun AM0 efficiencies in excess of 26 and 28 percent, respectively, are possible.
Educato, J. L.↗
Photospheric Current Spikes and Their Possible Association with Flares - Results from an HMI Data Driven Model
A data driven, near photospheric magnetohydrodynamic model predicts spikes in the horizontal current density, and associated resistive heating rate per unit volume Q. The spikes appear as increases by orders of magnitude above background values in neutral line regions (NLRs) of active regions (ARs). The largest spikes typically occur a few hours to a few days prior to M or X flares. The spikes correspond to large vertical derivatives of the horizontal magnetic field. The model takes as input the photospheric magnetic field observed by the Helioseismic & Magnetic Imager (HMI) on the Solar Dynamics Observatory (SDO) satellite. This 2.5 D field is used to determine an analytic expression for a 3 D magnetic field, from which the current density, vector potential, and electric field are computed in every AR pixel for 14 ARs. The field is not assumed to be force-free. The spurious 6, 12, and 24 hour Doppler periods due to SDO orbital motion are filtered out of the time series of the HMI magnetic field for each pixel using a band pass filter. The subset of spikes analyzed at the pixel level are found to occur on HMI and granulation scales of 1 arcsec and 12 minutes. Spikes are found in ARs with and without M or X flares, and outside as well as inside NLRs, but the largest spikes are localized in the NLRs of ARs with M or X flares. The energy to drive the heating associated with the largest current spikes comes from bulk flow kinetic energy, not the electromagnetic field, and the current density is highly non-force free. The results suggest that, in combination with the model, HMI is revealing strong, convection driven, non-force free heating events on granulation scales, and that it is plausible these events are correlated with subsequent M or X flares. More and longer time series need to be analyzed to determine if such a correlation exists. Above an AR dependent threshold value of Q, the number of events N(Q) with heating rates greater than or equal to Q obeys a scale invariant power law distribution for each AR given by N(Q) varies Q(sup -s), where 0.40 less than or equal to S less than or equal to 0.53, with a mean and standard deviation across the 14 ARs of 0.47 and 0.045, showing there is little variation of S from one AR to another. These properties of N(Q) are in close agreement with those of the distribution N(E) for the total energy E of solar flares, determined from observations to be N(E) = constant x E(sup -alpha). From observations of nanoflares in the 0.7 to 4 MK range, and from observations of flares in hard X-rays, it is found that 0.51 less than or equal to alpha less than or equal to 0.57, and 0.4 less than or equal to alpha less than or equal to 0.6, respectively (Crosby et al. 1993, Sol. Phys., 143, 275; Aschwanden & Parnell 2002, ApJ, 572, 1048). Observations also show that, as is found here for the exponent S, there is little variation of alpha with AR (Wheatland 2000, ApJ, 532, 1209), indicating N(E) and N(Q) are largely independent of individual properties of ARs such as area, total magnetic flux, and distribution of current density (i.e. non-potentiality). Therefore the power law scaling of the photospheric heating rate Q computed here on granulation scales is essentially identical to that found for coronal observations of flare energies on scales 1-2 orders of magnitude larger. This suggests the physical mechanisms that cause Q and coronal flares are closely related. It seems likely that Q is the signature of a magnetic reconnection process in an energy range and volume orders of magnitude smaller than those of flares. In this context, at least the larger spikes in Q might be signatures of UV photospheric or lower chromospheric bombs in which plasma is heated to temperatures approximately 10(exp -5) K (Peter et al. 2014, Science 346, 1255726; Judge 2015, ApJ, 808, 116). In addition, lattice based avalanche simulations of flare energy release predict 0.4 less than or equal to alpha less than or equal to 0.5, while analytic, fractal-diffusive self-organized criticality models predict 0.4 less than or equal to alpha less than or equal to 0.67, in excellent agreement with observations, and the results presented here (Aschwanden & Parnell 2002, ApJ, 572, 1048; Aschwanden 2012, A&A, 539, A2; Aschwanden 2013, in "Self Organized Criticality Systems"; Aschwanden et al. 2016, SSR, 198, 47).
Goodman, Michael↗
Solid-state Precipitation in Molecular Dynamics: KMC-MD Hybrid Simulations
A critical drawback of molecular dynamics (MD) simulations when modeling physical properties is the limitation to short timescales. Solid state precipitation is one behavior that is not significantly observable at timescales accessible to MD, and thus has been difficult to simulate at an MD scale. Taking inspiration from Kinetic Monte-Carlo (KMC) methods for simulating diffusion on a rigid lattice, we have developed a hybrid KMC-MD approach to enable simulation of solid-state precipitation during MD simulation implemented in LAMMPS. This method incorporates the atomic neighbor hopping behavior utilized by KMC simulations into a MD simulation with a variable local lattice configuration. By accelerating nearest-neighbor atomic hopping with this technique, solid-state precipitation can be observed at an MD timescale. Utilizing theory based on persistent-embryo methods for liquid-solid precipitation in MD we demonstrate the calculation of critical nucleus size for the gamma’ precipitation out of gamma phase in the Ni-Al system.
Molecular Dynamics↗
Solid-state Precipitation in Molecular Dynamics: kMC-MD Hybrid Simulations
A critical drawback of molecular dynamics (MD) simulations when modeling physical properties is the limitation to short timescales. Solid state precipitation is one behavior that is not significantly observable at timescales accessible to MD, and thus has been difficult to simulate at an MD scale. Taking inspiration from Kinetic Monte-Carlo (KMC) methods for simulating diffusion on a rigid lattice, we have developed a hybrid KMC-MD approach to enable simulation of solid-state precipitation during MD simulation implemented in LAMMPS. This method incorporates the atomic neighbor hopping behavior utilized by KMC simulations into a MD simulation with a variable local lattice configuration. By accelerating nearest-neighbor atomic hopping with this technique, solid-state precipitation can be observed at an MD timescale. Utilizing theory based on persistent-embryo methods for liquid-solid precipitation in MD we demonstrate the calculation of critical nucleus size for the gamma’ precipitation out of gamma phase in the Ni-Al system.
Molecular Dynamics↗
Unsteady Flow Simulation of a Sweeping Jet Actuator Using a Lattice-Boltzmann Method
Active flow control technology is increasingly used in aerospace applications to control flow separation and to improve aerodynamic performance. In this paper, PowerFLOW is used to simulate the flow through a sweeping jet actuator at two different pressure ratios. The lower pressure ratio leads to a high subsonic flow, whereas the high pressure ratio produces a choked flow condition. Comparison of numerical results with experimental data is shown, which includes qualitatively good agreement of pressure histories and spectra. PIV measurements are also available but the simulation overestimates mean and fluctuation quantities outside the actuator. If supply pressure is matched at one point inside the mixing chamber a good qualitative agreement is achieved at all other monitor points.
Duda, B.↗
Determination of Interfacial Energy for Solid-State Precipitation
Precipitation of strengthening phases from solid state is a critical process for many materials. For example, the high-temperature creep resistance of Ni-based superalloys is governed by ’ precipitates. However, the interfacial energy which controls the precipitation is expensive to determine experimentally. For simulations, the atomic nature of precipitation events necessitates using Molecular Dynamics (MD) so relying on fitting to a-priori experimental behavior is not required. However, the solid-state precipitation takes orders of magnitude longer time than is accessible for traditional MD simulation. Taking inspiration from Kinetic Monte-Carlo (KMC) methods for simulating diffusion on a rigid lattice, we have developed a hybrid KMC-MD approach to enable simulation of solid-state precipitation (and determination of the interfacial energy). Implementing these methods in LAMMPS, we calculated the interfacial energy for the ’ precipitation out of a super-saturated phase in Ni-based superalloys. The resulting interfacial energy value for the Ni-Al system is well within the accuracy of experimental values, demonstrating the capability of this method. Enabling the calculation of the interfacial energy without a-priori experimental knowledge of the precipitation unlocks a more efficient ICME-based approach for guiding annealing protocols in Ni-based superalloys and other solid-state precipitation systems.
Precipitation↗
Comparison of Boeing 777 Landing Gear Noise Simulations with Flight Test Data
Acoustic phased microphone array measurements of aircraft flyover noise acquired during the 2005 Quiet Technology Demonstrator II test were used to assess the accuracy of high-fidelity, full-scale simulations of landing gear noise produced by a large civilian aircraft. The simulations, conducted with the lattice Boltzmann solver PowerFLOW®, used a highly accurate digital model of a Boeing 777-300ER aircraft with the nose and main landing gear components replicating the full-scale geometries. The simulations were performed for aircraft parameters that matched those recorded during the flyover test conditions. For benchmarking purposes, several aircraft configurations were simulated: a) nose landing gear deployed with main landing gear and wing high-lift devices stowed, b) nose and main landing gear deployed with wing high-lift devices stowed and c) nose and main landing gear with wing high-lift devices deployed. To facilitate direct comparison with measured data, the simulated data sets were used to generate synthetic pressure records at the same array microphone locations as those used during the flight test. Broadly self-consistent beamforming techniques and procedures were used to process the synthetic pressure records and the measured data. Integration of select regions of the beamform maps containing the nose or main landing gear yielded good agreement between predicted and measured integrated far-field spectra for forward directivity angles where airframe noise is more prominent.
airframe noise↗
Mechanical cosmology: Simulating scalar fluctuations in expanding universes using synthetic mechanical lattices
Explore the source record for details and available documents.
Brendan Rhyno↗
A generalized vortex lattice method for subsonic and supersonic flow applications
If the discrete vortex lattice is considered as an approximation to the surface-distributed vorticity, then the concept of the generalized principal part of an integral yields a residual term to the vorticity-induced velocity field. The proper incorporation of this term to the velocity field generated by the discrete vortex lines renders the present vortex lattice method valid for supersonic flow. Special techniques for simulating nonzero thickness lifting surfaces and fusiform bodies with vortex lattice elements are included. Thickness effects of wing-like components are simulated by a double (biplanar) vortex lattice layer, and fusiform bodies are represented by a vortex grid arranged on a series of concentrical cylindrical surfaces. The analysis of sideslip effects by the subject method is described. Numerical considerations peculiar to the application of these techniques are also discussed. The method has been implemented in a digital computer code. A users manual is included along with a complete FORTRAN compilation, an executed case, and conversion programs for transforming input for the NASA wave drag program.
Miranda, L. R.↗
Equilibrium and kinetic aspects of brittle fracture
The lattice statics method for simulating a brittle crack will be employed to determine: (1) the energy-based criterion for crack stability, (2) at a given stress intensity in the stable region, the kinetic barrier for crack motion, and (3) the assumptions that break down in the Griffith theory when generalized to the lattice model. A simple model is used in the analysis, but this model is general enough to provide a conceptual framework for more realistic models.
Esterling, D. M.↗
Extended applications of the vortex lattice method
The application of the vortex lattice method to problems not usually dealt with by this technique is considered. It is shown that if the discrete vortex lattice is considered as an approximation to surface-distributed vorticity, then the concept of the generalized principal part of an integral yields a residual term to the vortex-induced velocity that renders the vortex lattice method valid for supersonic flow. Special schemes for simulating non-zero thickness lifting surfaces and fusiform bodies with vortex lattice elements are presented. Thickness effects of wing-like components are simulated by a double vortex lattice layer, and fusiform bodies are represented by a vortex grid arranged on a series of concentric cylindrical surfaces. Numerical considerations peculiar to the application of these techniques are briefly discussed.
Miranda, L. R.↗
Aeroacoustic Study of a Subscale Large Civil Transport (STAR) Model – Part 2: Validation of Simulated Results
Aeroacoustic measurements of the 26%-scale, semispan Boeing 777-200 Subsonic Transport Aeroacoustic Research (STAR) model tested in the NASA Ames Research Center 40- by 80-foot wind tunnel were used to ascertain the efficacy of high-fidelity simulations to accurately predict noise from the landing gear of large commercial transports. The simulations, conducted with the lattice Boltzmann solver PowerFLOW®, used a digital replica of the STAR model with or without main landing gear deployed and slats and flaps set to their highest deflection angles to represent aircraft during landing. The computations were performed at a Mach number of 0.21, Reynolds number of 8.2 million based on the model mean aerodynamic chord, and other conditions prevalent during the STAR model test. Measured and computed surface pressures were in very good agreement at most port locations on the model, as were global force coefficients, indicating that the simulations captured the impact of main gear deployment on inboard flap loading. Noise sources produced by the main landing gear and high-lift devices were determined via source localization maps generated with CLEAN from synthetic and experimental data. In general, very good agreement between predicted and measured acoustic source location and relative strength was observed in the maps. Comparisons of far-field noise spectra obtained from the CLEAN deconvolution maps showed remarkable agreement between synthetic and experimental broadband noise at low and medium frequencies. Main landing gear sources for model-scale frequencies above 7,000 Hz could not be resolved with the spatial resolution used during the simulations.
airframe noise↗
High-performance multiprocessor architecture for a 3-D lattice gas model
The lattice gas method has recently emerged as a promising discrete particle simulation method in areas such as fluid dynamics. We present a very high-performance scalable multiprocessor architecture, called ALGE, proposed for the simulation of a realistic 3-D lattice gas model, Henon's 24-bit FCHC isometric model. Each of these VLSI processors is as powerful as a CRAY-2 for this application. ALGE is scalable in the sense that it achieves linear speedup for both fixed and increasing problem sizes with more processors. The core computation of a lattice gas model consists of many repetitions of two alternating phases: particle collision and propagation. Functional decomposition by symmetry group and virtual move are the respective keys to efficient implementation of collision and propagation.
Lee, F.↗
Improved Fermion Hamiltonians for Quantum Simulations
Constructing improved hamiltonians for gauge theories coupled to fermonic matter will be important for improving continuum limit extrapolations of quantum computations. In this talk we will present a formulation for simulating ASQTAD fermions for lattice computation and provide fault tolerant resource costs in terms of primitive group operations. We additionally show that the scaling of energies with respect to the lattice spacing are better than for the unimproved Hamiltonian for toy models.
Erik Joseph Gustafson↗