Search NASA⌕ Search

SEARCH · Search NASA

Results for “force method”

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 793 records · Page 44
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

The Role of SnO2 Processing on Ionic Distribution in Double-Cation-Double Halide Perovskites

Moving toward a future of efficient, accessible, and less carbon-reliant energy devices has been at the forefront of energy research innovations for the past 30 years. Metal-halide perovskite (MHP) thin films have gained significant attention due to their flexibility of device applications and tunable capabilities for improving power conversion efficiency. Serving as a gateway to optimize device performance, consideration must be given to chemical synthesis processing techniques. Therefore, how does common substrate processing techniques influence the behavior of MHP phenomena such as ion migration and strain? Here, we demonstrate how a hybrid approach of chemical bath deposition (CBD) and nanoparticle SnO2 substrate processing significantly improves the performance of (FAPbI3)0.97(MAPbBr3)0.03 by reducing micro-strain in the SnO2 lattice, allowing distribution of K+ from K-Cl treatment of substrates to passivate defects formed at the interface and produce higher current in light and dark environments. X-ray diffraction reveals differences in lattice strain behavior with respect to SnO2 substrate processing methods. Through use of conductive atomic force microscopy (c-AFM), conductivity is measured spatially with MHP morphology, showing higher generation of current in both light and dark conditions for films with hybrid processing. Additionally, time-of-flight secondary ionization mass spectrometry (ToF-SIMS) observed the distribution of K+ at the perovskite/SnO2 interface, indicating K+ passivation of defects to improve the power conversion efficiency (PCE) and device stability. We show how understanding the role of ion distribution at the SnO2 and perovskite interface can help reduce the creating of defects and promote a more efficient MHP device.

conductive atomic force microscopy↗
NASA NTRS1977

Non-linear analysis of vibrations of irregular plates

A numerical perturbation method is used to investigate the forced vibrations of irregular plates. Nonlinear terms associated with the midplane stretching are retained in the analysis. The numerical part of the method involves the use of linear, finite element techniques to determine the free oscillation mode shapes and frequencies and to obtain the linear midplane stress resultants caused by the midplane stretching. Representing the solution as an expansion in terms of these linear mode shapes, these modes and the resultants are used to determine the equations governing the time-dependent coefficients of this expansion. These equations are solved by using the method of multiple scales. Specific solutions are given for the main-resonant vibrations of an elliptical plate in the presence of internal resonances. The results indicate that modes other than the driven mode can be drawn into the steady state response. Though the excitation is composed of a single harmonic, the response may not be periodic. Moreover, the particular types of responses that can occur are highly dependent on the mode being excited and are sensitive to small geometrical changes.

Lobitz, D. W.↗
NASA NTRS2014

Use of Machine Learning Techniques for Iidentification of Robust Teleconnections to East African Rainfall Variability in Observations and Models

Providing advance warning of East African rainfall variations is a particular focus of several groups including those participating in the Famine Early Warming Systems Network. Both seasonal and long-term model projections of climate variability are being used to examine the societal impacts of hydrometeorological variability on seasonal to interannual and longer time scales. The NASA / USAID SERVIR project, which leverages satellite and modeling-based resources for environmental decision making in developing nations, is focusing on the evaluation of both seasonal and climate model projections to develop downscaled scenarios for using in impact modeling. The utility of these projections is reliant on the ability of current models to capture the embedded relationships between East African rainfall and evolving forcing within the coupled ocean-atmosphere-land climate system. Previous studies have posited relationships between variations in El Niño, the Walker circulation, Pacific decadal variability (PDV), and anthropogenic forcing. This study applies machine learning methods (e.g. clustering, probabilistic graphical model, nonlinear PCA) to observational datasets in an attempt to expose the importance of local and remote forcing mechanisms of East African rainfall variability. The ability of the NASA Goddard Earth Observing System (GEOS5) coupled model to capture the associated relationships will be evaluated using Coupled Model Intercomparison Project Phase 5 (CMIP5) simulations.

Roberts, J. Brent↗
DOE OSTI2017

Crack propagation in silica from reactive classical molecular dynamics simulations

Mechanistic insight into the process of crack growth can be obtained through molecular dynamics (MD) simulations. In this investigation of fracture propagation, a slit crack was introduced into an atomistic amorphous silica model and mode I stress was applied through far–field loading until the crack propagates. Atomic displacements and forces and an Irving–Kirkwood method with a Lagrangian kernel estimator were used to calculate the J –integral of classical fracture mechanics around the crack tip. The resulting fracture toughness (K IC ), 0.76 ± 0.16 MPa√m, agrees with experimental values. In addition, the stress fields and dissipation energies around the slit crack indicate the development of an inelastic region ~30Å in diameter. This is one of the first reports of KIC values obtained from up–scaled atomic–level energies and stresses through the J –integral. In conclusion, the application of the ReaxFF classical MD force field in this study provides the basis for future research into crack growth in multicomponent oxides in a variety of environmental conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Optimized Electromagnetic Fields Levitate Bulk Metallic Materials

The bulk metallic electromagnetic levitation technique has a promising application prospect in condensed physics and materials science. In order to resolve the problem of weak levitation force and the small size of samples, an optimized solution based on finite element analysis (FEA) is developed to seek the best electromagnetic field distribution in the levitation zone. The influence of material properties on the distribution of magnetic flux density, surface pressure, and the Lorentz force were studied by the FEA method with various sample sizes of titanium and copper balls. The levitation capability of cylindrical coil system which consists of a single-layer coil with an inner diameter of 30 mm and the levitated object with a diameter of 20 mm was analyzed. To keep a balance between the turn number and increment efficiency of force, the best number of 7 turns is confirmed. Then, two types of optimized coils were predicted by numerical simulations, and the levitation force is improved drastically. The double-layer coil displays an outstanding levitation capability owing to the strong electromagnetic field generated by multilayer turns. In addition, 5 to 50 g aluminum, titanium, copper, and nickel are levitated by the above three kinds of coils and the levitation capacities are validated from the recorded data of the minimum operating current that can levitate the samples. The experimental results show good agreements with the simulation results, and the minimum current for optimized coils is reduced efficiently.

36 MATERIALS SCIENCE↗
NASA NTRS2004

Turbomachine Sealing and Secondary Flows: Review of Rotordynamics Issues in Inherently Unsteady Flow Systems With Small Clearances - Part 2

Today's computational methods enable the determination of forces in complex systems, but without field validation data, or feedback, there is a high risk of failure when the design envelope is challenged. The data of Childs and Bently and field data reported in NASA Conference Proceedings serve as sources of design information for the development of these computational codes. Over time all turbomachines degrade and instabilities often develop, requiring responsible, accurate, turbomachine diagnostics with proper decisions to prevent failures. Tam et al. (numerical) and Bently and Muszynska (analytical) models corroborate and implicate that destabilizing factors are related through increases in the fluid-force average circumferential velocity. The stability threshold can be controlled by external swirl and swirl brakes and increases in radial fluid film stiffness (e.g., hydrostatic and ambient pressures) to enhance rotor stability. Also cited are drum rotor self-excited oscillations, where the classic fix is to add a split or severed damper ring or cylindrical damper drum, and the Benkert-Wachter work that engendered swirl brake concepts. For a smooth-operating, reliable, long-lived machine, designers must pay very close attention to sealing dynamics and diagnostic methods. Correcting the seals enabled the space shuttle main engine high-pressure fuel turbopump (SSME HPFTP) to operate successfully.

Hendricks, R. C.↗
NASA NTRSDate not supplied

Evaluation of Preliminary Buffet Forcing Function Development based on Wind-Tunnel Tests of Geometrically-similar Models

The current state-of-the-art for launch vehicle buffet load estimation involves wind tunnel testing of highly instrumented rigid models and data processing of experimental unsteady pressure measurements to estimate the fluctuating aerodynamic loads. The model and instrumentation requirements to accomplish this type of testing are significant. The rigid model, in particular, must be built with relatively high geometric fidelity to approximate similar aerodynamic phenomena that will occur during the full-scale vehicle flight. Due to this requirement, several iterative wind tunnel tests are often required during the design maturation of the vehicle. For preliminary buffet load estimation, however, data from previous launch vehicle wind tunnel tests of geometrically similar models can be used to provide initial buffet environment estimates before a high geometric fidelity wind tunnel test has been conducted. In this paper, two components of this preliminary buffet estimation process are described: (1) assembly of buffet forcing functions from multiple data sources and (2) scaling of these buffet forcing functions to geometrically-similar components of a target vehicle configuration. The impacts of these estimation methods on the full-vehicle buffet forcing functions, including their longitudinal distribution along the vehicle, frequency-content at individual stations, and station-to-station coherence, are investigated. Considerations for using these methods to extrapolate wind tunnel data to geometrically similar models are also discussed.

Patrick S Heaney↗
NASA NTRS1990

Use of magnetic sails for advanced exploration missions

The magnetic sail, or magsail, is a field effect device which interacts with the ambient solar wind or interstellar medium over a considerable volume of space to generate drag and lift forces. Two theories describing the method of thrust generation are analyzed and data results are presented. The techniques for maintaining superconductor temperatures in interplanetary space are analyzed and low risk options presented. Comparisons are presented showing mission performance differences between currently proposed spacecraft using chemical and electric propulsion systems, and a Magsail propelled spacecraft capable of generating an average thrust of 250 Newtons at a radius of one A.U. The magsail also provides unique capabilities for interstellar missions, in that at relativistic speeds the magnetic field would ionize and deflect the interstellar medium producing a large drag force. This would make it an ideal brake for decelerating a spacecraft from relativistic speeds and then maneuvering within the target star system.

Andrews, Dana G.↗
NASA NTRS2006

In-Vacuum Photogrammetry of a Ten-Meter Square Solar Sail

Solar sailing is a promising, future in-space propulsion method that uses the small force of reflecting sunlight to accelerate a large, reflective membrane without expendable propellants. One of two solar sail configurations under development by NASA is a striped net approach by L'Garde, Inc. This design uses four inflatably deployed, lightweight booms supporting a network of thin strings onto which four quadrants of ultrathin aluminized membranes are attached. The NASA Langley Research Center (LaRC) provided both experimental and analytical support to L'Garde for validating the structural characteristics of this unique, ultralightweight spacecraft concept. One of LaRC's responsibilities was to develop and apply photogrammetric methods to measure sail shape. The deployed shape provides important information for validating the accuracy of finite-element modeling techniques. Photogrammetry is the science and art of calculating 3D coordinates of targets or other distinguishing features on structures using images. A minimum of two camera views of each target is required for 3D determination, but having four or more camera views is preferable for improved reliability and accuracy. Using retroreflective circular targets typically provides the highest measurement accuracy and automation. References 3 and 4 provide details of photogrammetry technology, and reference 5 discusses previous experiences with photogrammetry for measuring gossamer spacecraft structures such as solar sails. This paper discusses the experimental techniques used to measure a L Garde 10-m solar sail test in vacuum with photogrammetry. The test was conducted at the NASA-Glenn Space Power Facility (SPF) located at Plum Brook Station in Sandusky, Ohio. The SPF is the largest vacuum chamber in the United States, measuring 30 m in diameter by 37 m in height. High vacuum levels (10(exp -6) torr) can be maintained inside the chamber, and cold environments (-195 C) are possible using variable-geometry cryogenic cold walls. This test used a vacuum level of approximately 1 torr (sufficient for structural static/dynamic characterization) and instead of using the cryogenic cold walls, used local LN2 cold plates underneath each of the four cold-rigidizable solar sail booms instead.

Pappa, Richard S.↗
DOE OSTI2018

Functionalisation of Detonation Nanodiamond for Monodispersed, Soluble DNA-Nanodiamond Conjugates Using Mixed Silane Bead-Assisted Sonication Disintegration

Here, nanodiamonds have many attractive properties that make them suitable for a range of biological applications, but their practical use has been limited because nanodiamond conjugates tend to aggregate in solution during or after functionalisation. Here we demonstrate the production of DNA-detonation nanodiamond (DNA-DND) conjugates with high dispersion and solubility using an ultrasonic, mixed-silanization chemistry protocol based on the in situ Bead-Assisted Sonication Disintegration (BASD) silanization method. We use two silanes to achieve these properties: (1) 3-(trihydroxysilyl)propyl methylphosphonate (THPMP); a negatively charged silane that imparts high zeta potential and solubility in solution; and (2) (3-aminopropyl)triethoxysilane (APTES); a commonly used functional silane that contributes an amino group for subsequent bioconjugation. We target these amino groups for covalent conjugation to thiolated, single-stranded DNA oligomers using the heterobifunctional crosslinker sulfosuccinimidyl 4-(N-maleimidomethyl)cyclohexane-1-carboxylate (Sulfo-SMCC). The resulting DNA-DND conjugates are the smallest reported to date, as determined by Dynamic Light Scattering (DLS) and Atomic Force Microscopy (AFM). The functionalisation method we describe is versatile and can be used to produce a wide variety of soluble DND-biomolecule conjugates.

60 APPLIED LIFE SCIENCES↗
NASA NTRS1997

Evaluation of heart rate changes: electrocardiographic versus photoplethysmographic methods

The heart rate (HR) variation to forced deep breathing (HRDB) and to the Valsalva maneuver (Valsalva ratio; VR) are the two most widely used tests of cardiovagal function in human subjects. The HR is derived from a continuously running electrocardiographic (ECG) recording. Recently, HR derived from the arterial waveform became available on the Finapres device (FinapHR), but its ability to detect rapid changes in HR remains uncertain. We therefore evaluated HRDB and VR derived from FinapHR using ECG-derived HR (ECGHR) recordings as the standard. We also compared the averaged HR on Finapres (Finapav) with beat-to-beat Finapres (FinapBB) values. Studies were undertaken in 12 subjects with large HR variations: age, 34.5 +/- 9.3 (SD) years; six males and six females. FinapBB values were superimposable upon ECGHR for both HRDB and VR. In contrast, Finapav failed to follow ECGHR for HRDB and followed HRECG with a lag for the VR. To evaluate statistically how closely FinapHR approximated ECGHR, we undertook regression analysis, using mean values for each subject. To compare the two methods, we evaluated the significance of the difference between test and standard values. For HRDB, FinapBB reproducibly recorded HR (R2 = 0.998), and was significantly (p = 0.001) better than Finapav (R2 = 0.616; p < 0.001). For VR, HRBB generated a VR that was not significantly different from the correct values, while HRav generated a value that was slightly but consistently lower than the correct values (p < 0.001). We conclude that FinapHR reliably records HR variations in the beat-to-beat mode for cardiovascular HR tests.

Non-NASA Center↗
NASA NTRS1975

On the prediction of critical flutter conditions from subcritical response data and some related wind-tunnel experience

Methods of interpreting response measurements which could be amenable to flight flutter testing procedures were studied analytically and in the wind tunnel. One suggested scheme, which requires evaluation, is an iterative technique in which derivatives obtained from subcritical response data are used to indicate the approach to flutter. A simplification of this procedure was considered by examining the manner in which a single characteristic of the subcritical response behaves in relation to variations of the density or dynamic pressure in the approach to flutter. The use of this single parameter scheme was examined for random excitation as well as for sinusoidal forcing. The feasibility of the method is illustrated by several examples and the relative merits of random and sinusoidal excitation are discussed.

Houbolt, J. C.↗
NASA NTRSDate not supplied

Correlated Electromagnetic Levitation Actuator: A Reaction Sphere Based Attitude Control System

To address problems experienced by current reaction wheels and control moment gyroscopebased attitude control systems (ACS), researchers at NASA’s Marshall Space Flight Center have begun developing a reaction sphere actuator based on correlated electromagnetic levitation that will be immune to destructive bearing friction, momentum saturation, and gimbal lock. The Correlated Electromagnetic Levitation Actuator (CELA) advances the state of the art of reaction sphere ACSs by employing the concept of correlated magnetics. It is a frictionless, direct-drive reaction sphere that harnesses a unique technology with an array of applications across multiple disciplines. Correlated electromagnets function in a manner that is analogous to a matched filter; the convolution of two signals is peaked at the index representing the greatest match. For CELA, the signals are the patterns of magnetic flux density as a function of position. The magnitude of the convolution equates to an attractive or repulsive force, and these forces can be azimuthal or radial. The development of CELA is based in four distinct disciplines: Advanced Manufacturing, Prototype Development, Electromagnetic Modeling, and Controls. We are developing novel manufacturing techniques required to build arrays of permanent and electromagnet dipoles on curved surfaces. To print the permanent magnetic array, we have developed a probe with pyramidal magnets that will reside on a robotic arm to induce localized magnetic fields on a surface. The probe also includes the ability to erase dipole patterns from a permanent magnet by heating the surface to its Curie temperature. A number of test articles and prototypes have been developed using additive manufacturing methods. These prototypes have included hemispherical motors to test the drive algorithm, and a levitation test bed that demonstrates a magnetic bearing method based on attractive magnetic forces and ratiometric Hall effect sensors. We developed an array of electromagnetic dipoles on a printed circuit board (PCB) with individual H-bridges controlling each coil. This device created various flux density patterns and we measured their magnetic fields using a custom Hall effect 3-D probe and a LabVIEW virtual instrument. These data will serve as a benchmark for characterizing the accuracy of future models. Current work is focused on modeling the magnetic fields of our prototype arrays using COMSOL Finite Element Analysis and verifying the model against our test data. Accurate modeling will allow us to quickly test new patterns of electromagnets and their macro behavior. Eventually, the magnetic field models will be implemented in our controls simulations to facilitate precise control of the reaction sphere. Initial model results agree with field measurements to within 1 G (5% of measured flux density). Currently, we are testing different material properties of the electromagnets and their magnetic fields and thermal effects. These results will be used to refine the design of the electromagnetic dipoles. Our control efforts have centered on developing commutation, levitation, and field pattern shaping hardware in the form of breadboards and PCBs with software running on a local microcontroller. In addition, our partners developed MATLAB Simulink models to demonstrate a PID controller thatmitigates disturbance forces resulting from the interaction of drive and levitation magnetics. Finally, we have designed a three-axis test stand that will be used in future work to demonstrate CELA’s orientation control capability.

controls↗
NASA NTRSDate not supplied

Correlated Electromagnetic Levitation Actuator: A Reaction Sphere-Based Attitude Control System

To address problems experienced by current reaction wheels and control moment gyroscopebased attitude control systems (ACS), researchers at National Aeronautics and Space Administration’s (NASA’s) Marshall Space Flight Center (MSFC) have begun developing a reaction sphere actuator based on correlated electromagnetic levitation that will be immune to destructive bearing friction, momentum saturation, and gimbal lock. The Correlated Electromagnetic Levitation Actuator (CELA) advances the state of the art of reaction sphere ACSs by employing the concept of correlated magnetics. It is a frictionless, direct-drive reaction sphere that harnesses a unique technology with an array of applications across multiple disciplines. Correlated electromagnets function in a manner that is analogous to a matched filter; the convolution of two signals is peaked at the index representing the greatest match. For CELA, the signals are the patterns of magnetic flux density as a function of position. The magnitude of the convolution equates to an attractive or repulsive force, and these forces can be azimuthal or radial. The development of CELA is based in four distinct disciplines: Advanced Manufacturing, Prototype Development, Electromagnetic Modeling, and Controls. We are developing novel manufacturing techniques required to build arrays of permanent and electromagnet dipoles on curved surfaces. To print the permanent magnetic array, we have developed a probe with pyramidal magnets that will reside on a robotic arm to induce localized magnetic fields on a surface. We have also used high temperature ovens to erase dipole patterns from a permanent magnet by heating the surface to its Curie temperature. A number of test articles and prototypes have been developed using additive manufacturing methods. These prototypes have included hemispherical motors to test the drive algorithm, and a levitation test bed that demonstrates a magnetic bearing method based on attractive magnetic forces and ratiometric Hall effect sensors. We developed an array of electromagnetic dipoles on a printed circuit board (PCB) with individual H-bridges controlling each coil. This device created various flux density patterns and we measured their magnetic fields using a custom Hall effect 3-D probe and a LabVIEW virtual instrument. These data will serve as a benchmark for characterizing the accuracy of future models. Current work is focused on modeling the magnetic fields of our prototype arrays using COMSOL finite element analysis (FEA) and verifying the model against our test data. Accurate modeling will allow us to quickly test new patterns of electromagnets and their macro behavior. Eventually, the magnetic field models will be implemented in our controls simulations to facilitate precise control of the reaction sphere. Initial model results agree with field measurements to within 1 G (5% of measured flux density). Currently, we are testing different material properties of the electromagnets and their magnetic fields and thermal effects. These results will be used to refine the design of the electromagnetic dipoles. Our control efforts have centered on developing commutation, levitation, and field pattern shaping hardware in the form of breadboards and PCBs with software running on a local microcontroller. In addition, our partners developed MATLAB Simulink models to demonstrate a PID controller that mitigates disturbance forces resulting from the interaction of drive and levitation magnetics. Finally, we have designed a three-axis test stand that will be used in future work to demonstrate CELA’s orientation control capability.

reaction sphere↗
DOE OSTI2020

Magnetoelastic sensor for analyzing strain

A magnetoelastic sensor is provided for analyzing strain induced by a force causing mechanical deformation of a region of a substrate. The sensor includes a substrate that has a non-deformable region and a mechanically deformable region, the mechanically deformable region coupled to the non-deformable region by a linking structure. The sensor also has a transformer having a magnetic flux guide. The flux guide includes a magnetoelastic section having a magnetoelastic subsection configured to be strainable by deformation of the mechanically deformable region. In transformer embodiments, one or more drive coil and/or pickup coil has a width greater than the minimal width of the linking structure. Sensor embodiments include numerous configurations. The sensor is useful in methods for analyzing strain and force.

47 OTHER INSTRUMENTATION↗
DOE OSTI2018

Quantifying Tip-Sample Interactions in Vacuum Using Cantilever-Based Sensors: An Analysis

Atomic force microscopy is an analytical characterization method that is able to image a sample’s surface topography at high resolution while simultaneously probing a variety of different sample properties. Such properties include tip-sample interactions, the local measurement of which has gained much popularity in recent years. To this end, either the oscillation frequency or the oscillation amplitude and phase of the vibrating force-sensing cantilever are recorded as a function of tip-sample distance and subsequently converted into quantitative values for the force or interaction potential. Here, we theoretically and experimentally show that the force law obtained from such data acquired under vacuum conditions using the most commonly applied methods may deviate more than previously assumed from the actual interaction when the oscillation amplitude of the probe is of the order of the decay length of the force near the surface, which may result in a non-negligible error if correct absolute values are of importance. Caused by approximations made in the development of the mathematical reconstruction procedures, the related inaccuracies can be effectively suppressed by using oscillation amplitudes sufficiently larger than the decay length. To facilitate efficient data acquisition, we propose a technique that includes modulating the drive amplitude at a constant height from the surface while monitoring the oscillation amplitude and phase. Ultimately, such an amplitude-sweep-based force spectroscopy enables shorter data acquisition times and increased accuracy for quantitative chemical characterization compared to standard approaches that vary the tip-sample distance. An additional advantage is that since no feedback loop is active while executing the amplitude sweep, the force can be consistently recovered deep into the repulsive regime.

36 MATERIALS SCIENCE↗
NASA NTRS2018

Adjustments and Uncertainty Quantification for SLS Aerodynamic Sectional Loads

This paper presents a method for adjusting sectional loads to match target values for integrated force and moment coefficients. In a typical application, the sectional load profile for one flight condition is calculated from Computational Fluid Dynamics (CFD) while the integrated forces and moments are measured in a wind tunnel experiment. These two methods do not generally result in identical predictions, and this leads to an inherent inconsistency between different data products. This paper aims to provide a procedure to remove that inconsistency. A sectional load profile for a launch vehicle splits the rocket into slices along its length and calculates the aerodynamic loading on each slice, which leads to a one-dimensional aerodynamic load profile that is used for structural analysis. Adjusting sectional loads, also known as line loads, is a nontrivial matter due to several consistency constraints. For example, the adjusted sectional normal force profile must be consistent with both the integrated normal force and pitching moment. To avoid such inconsistency issues, this paper presents a method using a Proper Orthogonal Decomposition (POD) to generate basis functions to adjust the sectional load profiles. As a corollary, this correction method enables the creation of an uncertainty quantification for sectional loads that is consistent with the dispersed integrated force and moment database and its uncertainty quantification. Several extensions to this technique, such as applying the method to the surface pressures, are considered.

Quantification↗
DOE OSTI2018

Soft-sphere simulations of a planar shock interaction with a granular bed

Here we consider the problem of shock propagation through a layer of spherical particles. A point particle force model is used to capture the shock-induced aerodynamic force acting upon the particles. The discrete element method (DEM) code liggghts is used to implement the shock-induced force as well as to capture the collisional forces within the system. A volume-fraction-dependent drag correction is applied using Voronoi tessellation to calculate the volume of fluid around each individual particle. A statistically stationary frame is chosen so that spatial and temporal averaging can be performed to calculate ensemble-averaged macroscopic quantities, such as the granular temperature. A parametric study is carried out by varying the coefficient of restitution for three sets of multiphase shock conditions. A self-similar profile is obtained for the granular temperature that is dependent on the coefficient of restitution. A traveling wave structure is observed in the particle concentration downstream of the shock and this instability arises from the volume-fraction-dependent drag force. The intensity of the traveling wave increases significantly as inelastic collisions are introduced. Downstream of the shock, the variance in Voronoi volume fraction is shown to have a strong dependence upon the coefficient of restitution, indicating clustering of particles induced by collisional dissipation. Statistics of the Voronoi volume are computed upstream and downstream of the shock and compared to theoretical results for randomly distributed hard spheres.

42 ENGINEERING↗