Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “performance prediction”

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 559 records · Page 31

Propulsion-airframe interactions predictability

An illustration of the possible magnitude of the transonic performance prediction problem is shown. These results derived from a careful and systematic set of tests run on a jet effects model, an aerodynamic and force model, and a pressure model, evaluated test techniques and wind tunnel effects. The results were then carefully compared with flight test data where significant discrepancies in zero-lift drag were observed. A vigorous review of the test techniques identified a number of potential error sources. Some of these were thought to be aircraft roughness and protuberances, tunnel anomalies, sting/model support corrections, and hot gas effects. In addition, there were questions on model metric splitline location and magnitude of corrections for scale and Reynolds number. Flight test inlet/engine characteristics were thought to differ, also, from those of the calibration settings in the ground facilities.

Smith, R. H.↗

Design of a one-year lifetime, spaceborne superfluid helium dewar

The Infrared Astronomical Satellite is an Explorer Mission and a joint venture of the Netherlands, the United Kingdom, and the United States scheduled for launch into earth orbit in 1981. The cryogenic system is a major part of the satellite; it incorporated many unique and state-of-the-art design features to satisfy the requirements of a one-year orbital lifetime, a focal plane temperature less than 4K, minimal launch weight, and zero-gravity operation. The 60-centimeter diameter telescope is contained within a superfluid helium dewar having a capacity of 540 liters. The telescope aperture cover employs an independent cryogenic system containing a 54-liter supercritical helium tank. The aperture cover, which is ejected two weeks after launch, protects the telescope from contamination and provides the low-temperature background needed to perform focal plane health checks. Design and predicted performance of the cryogenic systems are discussed in detail.

Hopkins, R. A.↗

Lunar base thermal management/power system analysis and design

A compilation of several lunar surface thermal management and power system studies completed under contract and IR&D is presented. The work includes analysis and preliminary design of all major components of an integrated thermal management system, including loads determination, active internal acquisition and transport equipment, external transport systems (active and passive), passive insulation, solar shielding, and a range of lunar surface radiator concepts. Several computer codes were utilized in support of this study, including RADSIM to calculate radiation exchange factors and view factors, RADIATOR (developed in-house) for heat rejection system sizing and performance analysis over a lunar day, SURPWER for power system sizing, and CRYSTORE for cryogenic system performance predictions. Although much of the work was performed in support of lunar rover studies, any or all of the results can be applied to a range of surface applications. Output data include thermal loads summaries, subsystem performance data, mass, and volume estimates (where applicable), integrated and worst-case lunar day radiator size/mass and effective sink temperatures for several concepts (shielded and unshielded), and external transport system performance estimates for both single and two-phase (heat pumped) transport loops. Several advanced radiator concepts are presented, along with brief assessments of possible system benefits and potential drawbacks. System point designs are presented for several cases, executed in support of the contract and IR&D studies, although the parametric nature of the analysis is stressed to illustrate applicability of the analysis procedure to a wide variety of lunar surface systems. The reference configuration(s) derived from the various studies will be presented along with supporting criteria. A preliminary design will also be presented for the reference basing scenario, including qualitative data regarding TPS concerns and issues.

Mcghee, Jerry R.↗

Analysis of the Shuttle Orbiter reinforced carbon-carbon oxidation protection system

Reusable, oxidation-protected reinforced carbon-carbon (RCC) has been successfully flown on all Shuttle Orbiter flights. Thermal testing of the silicon carbide-coated RCC to determine its oxidation characteristics has been performed in convective (plasma Arc-Jet) heating facilities. Surface sealant mass loss was characterized as a function of temperature and pressure. High-temperature testing was performed to develop coating recession correlations for predicting performance at the over-temperature flight conditions associated with abort trajectories. Methods for using these test data to establish multi-mission re-use (i.e., mission life) and single mission limits are presented.

Williams, S. D.↗

Advanced Technology Composite Fuselage-Structural Performance

Boeing is studying the technologies associated with the application of composite materials to commercial transport fuselage structure under the NASA-sponsored contracts for Advanced Technology Composite Aircraft Structures (ATCAS) and Materials Development Omnibus Contract (MDOC). This report addresses the program activities related to structural performance of the selected concepts, including both the design development and subsequent detailed evaluation. Design criteria were developed to ensure compliance with regulatory requirements and typical company objectives. Accurate analysis methods were selected and/or developed where practical, and conservative approaches were used where significant approximations were necessary. Design sizing activities supported subsequent development by providing representative design configurations for structural evaluation and by identifying the critical performance issues. Significant program efforts were directed towards assessing structural performance predictive capability. The structural database collected to perform this assessment was intimately linked to the manufacturing scale-up activities to ensure inclusion of manufacturing-induced performance traits. Mechanical tests were conducted to support the development and critical evaluation of analysis methods addressing internal loads, stability, ultimate strength, attachment and splice strength, and damage tolerance. Unresolved aspects of these performance issues were identified as part of the assessments, providing direction for future development.

Walker, T. H.↗

Progress in the Fabrication and Testing of Telescope Mirrors for the James Webb Space Telescope

The telescope of the James Webb Space Telescope (JWST) is an f/20, three mirror anastigmat design, passively cooled (40K) in an L2 orbit. The design provides diffraction limited performance (Strehl > or = 0.8) at .=2 m. To fit within the launch vehicle envelope (Arianne V), the 6.6 meter primary mirror and the secondary mirror support structure are folded for launch, then deployed and aligned in space. The primary mirror is composed of 18 individual, 1.3 meter (flat flat) hexagonal segments, each adjustable in seven degrees of freedom (six rigid body + radius of curvature) provided by a set of high precision actuators. The actuated secondary mirror (approx.0.74m) is similarly positioned in six degrees of rigid body motion. The approx..70x.51m, fixed tertiary and approx. 0.17m, flat fine steering mirror complete the telescope mirror complement. The telescope is supported by a composite structure optimized for performance at cryogenic temperatures. All telescope mirrors are made of Be with substantial light-weighting (21 kg for each 1.3M primary segment). Additional Be mounting and supporting structure for the high precision (approx.10nm steps) actuators are attached to the primary segments and secondary mirror. All mirrors undergo a process of thermal stabilization to reduce stress. An extensive series of interferometric measurements guide each step of the polishing process. Final polishing must account for any deformation between the ambient temperature of polishing and the cryogenic, operational temperature. This is accomplished by producing highly precise, cryo deformation target maps of each surface which are incorporated into the final polishing cycle. All flight mirrors have now completed polishing, coating with protected Au and final cryo testing, and the telescope is on track to meet all system requirements. We here review the measured performance of the component mirrors and the predicted performance of the flight telescope.

Bowers, Charles↗

GPS Navigation Above 76,000 km for the MMS Mission

NASA's MMS mission, launched in March of 2015,consists of a controlled formation of four spin-stabilized spacecraft in similar highly elliptic orbits reaching apogee at radial distances of 12and 25 Earth radii in the first and second phases of the mission. Navigation for MMS is achieved independently onboard each spacecraft by processing GPS observables using NASA GSFC's Navigator GPS receiver and the Goddard Enhanced Onboard Navigation System (GEONS) extended Kalman filter software. To our knowledge, MMS constitutes, by far, the highest-altitude operational use of GPS to date and represents the culmination of over a decade of high-altitude GPS navigation research and development at NASA GSFC. In this paper we will briefly describe past and ongoing high-altitude GPS research efforts at NASA GSFC and elsewhere, provide details on the design of the MMS GPS navigation system, and present on-orbit performance data. We extrapolate these results to predict performance in the Phase 2b mission orbit, and conclude with a discussion of the implications of the MMS results for future high-altitude GPS navigation, which we believe to be broad and far-reaching.

MMS↗

Characterization of Internal Insulation Thermal Performance

A test campaign was performed to build an SRM-insulation-material thermal-performance database using MSFC’s Solid Fuel Torch that could be leveraged to enhance understanding of the relationship between ablation rates and internal-environment parameters and improve insulation performance predictions. Seven hot-fire tests were performed along with data analysis that revealed the thermal performance differences among the tested materials. Functional Data Analysis was employed to create a model of material ablation rate as a function of six covariates. The relative importance of several environmental parameters was indicated, and the gaps in the operational space covered by the database were identified.

Ablation↗

On the Representativity of Electrode Microstructure Parameters and Their Electrochemical Response for Lithium Ion Batteries

Lithium-ion battery electrochemical models require an accurate description of the electrodes microstructures to be predictive that can be achieved through nanoscale imaging. Such observations are however limited by their field of view (FOV), as they provide only a subset of the whole electrode volume that does not necessarily represent the whole electrode microstructure heterogeneity, and therefore can bias the microstructure analysis. A microstructure scale electrochemical model was used to investigate lithium plating onset, material non-uniform utilization, and in-plane heterogeneities for an NMC-graphite full cell. To evaluate the representativeness, and thus relevance, of these model predictions, a coupled representativity analysis has been performed on the microstructure parameters and, in a novel way, on the full cell electrochemical response. Electrode microstructure parameters representativeness has been first quantified using the representative volume element (RVE) methodology. The RVE major flaw is that ultimately it can only conclude if a FOV contains representative subvolumes of the FOV, but not if the FOV itself is representative of the electrode volume. Analysis can conclude negatively ('FOV is not representative'), but not positively ('FOV is representative'). One major contribution of this work was to quantify the convergence of the RVE size with the FOV, to actually investigate the FOV representativeness and thus partly remedy this intrinsic limitation. The analysis determined that performing a standard RVE calculation, without exploring its FOV convergence, is likely to strongly underestimate the actual RVE size. The new RVE methodology has been automated in the NREL open-source Microstructure Analysis Toolbox (MATBOX) and is available to the battery community. Representativeness of microstructure parameters is however only an intermediate step, as the end-results of an electrochemical model are performances predictions. Indeed, what is the practical consequence of a given deviation for a microstructure parameter? The microstructure parameter deviation propagations to the 3D microstructure scale electrochemical response have been then quantified for different charge rates. This defines a threshold for the microstructure parameters FOV for a desired maximum deviation of the electrochemical response. Such deviation propagation analysis is analogous to error propagation analysis and is necessary to determine the relevance of microstructure scale model predictions for macroscale predictions. Electrochemical model shows cell representative section areas are increasing with C-rate, due to higher in-plane heterogeneities, indicating larger FOVs are required specifically for fast charge modeling. Therefore, we introduced the novel concept of electrochemical RVE (eRVE) that is a function of the operating conditions (thus defined as a dynamic RVE), with an increasing dependence with the C-rate. Representativity analysis of the investigated cell determined a FOV of 144.4 x 54.4 m2 is large enough to establish a convergence on the representative section areas for low to intermediate C-rate (=2.5C), but not large enough to conclude for higher rates. This work aims to emphasize the importance of representativity analysis for LIB electrode microstructures, as it is required to estimate the error, and thus the relevance, of microstructure parameters intended to be used in macroscale models. The methodology and results can help researchers to select the relevant imaging and associated FOV required to provide accurate enough microstructure parameters.

ADVANCED PROPULSION SYSTEMS↗

Computational Analysis of the Large Scale Low-Boom Supersonic Inlet

This presentation describes two computational fluid dynamic (CFD) analyses done in support of a supersonic inlet test performed at NASA Glenn Research Center in the fall of 2010. The large-scale-low-boom supersonic inlet was designed for a small supersonic aircraft that would cruise at a Mach number of 1.6. It uses an axisymmetric, external compression spike to reduce the Mach number to 0.65 at the fan face. The inlet was tested in the 8x6 supersonic wind tunnel at NASA GRC using conventional pressure probes, pressure sensitive paint, and high-speed schlieren. Two CFD analyses of the inlet were performed before the test, and compared to the experimental data afterwards. Both analyses used the WIND-US code. First, an axisymmetric analysis of the inlet, diffuser, cold pipe, and mass flow plug was performed to predict the performance of the entire system in the wind tunnel. Then a 3-D analysis of the inlet with all its interior struts was performed to predict details of the flow field and effects of angle of attack. Test results showed that the inlet had excellent performance, with a peak total pressure recovery of 96 percent, and a buzz point far outside the engine operating range. The computations agreed very well with the data, with predicted recoveries within 0.3 - 0.5 points of the measurements.

Chima, ROdrick V.↗

A hybrid method for prediction of propeller performance

An analytical/computational code has been developed to predict the performance of advanced turboprop propellers. The method of solution is based on strip theory and conservation of linear and angular momenta applied to a control volume enclosing the propeller and its far wake. An iterative procedure is used to account for the induced axial and rotational velocities. Sectional aerodynamic coefficients are explicitly accounted for by thin airfoil theory and refinement of the results is obtained using a two-dimensional vortex panel method and Euler equations. The results obtained compare well with wind-tunnel data obtained for advanced propellers, such as the Hamilton Standard SR-2 and SR-7.

Takallu, M. A.↗

Proton Effects and Test Issues for Satellite Designers

Microelectronic and photonic systems in the natural space environment are bombarded by a variety of charged particles including electrons, trapped protons, cosmic rays, and solar particles (protons and other heavy ions). These incident particles cause both ionizing and non-ionizing effects when traversing a device, and the effects can be either transient or permanent. The vast majority of the kinetic energy of an incident proton is lost to ionization, creating the single event effects (SEES) and total ionizing dose (TID) effects. However, the small portion of energy lost in non-ionizing processes causes atoms to be removed from their lattice sites and form permanent electrically active defects in semiconductor materials. These defects, i.e., "displacement damage," can significantly degrade device performance. In general, most of the displacement damage effects in the natural space environment can be attributed to protons since they are plentiful and extremely energetic (and therefore not readily shielded against). For this reason, we consider only proton induced displacement damage in this course. (Nevertheless, we identify solar cells as an important example of a case where both electron and proton damage can be important since only very light shielding is feasible.) The interested reader is encouraged to explore the three previous NSREC and RADECS short courses which also treat displacement damage issues for satellite applications. Part A of this segment of the short course introduces the space environment, proton shielding issues, and requirements specifications for proton-rich environments. In order to exercise the displacement damage analysis tools for on-orbit performance predictions, the requirements document must provide the relevant proton spectra in addition to the usual total ionizing dose-depth curves. Ion-solid interactions and the nature of the displacement damage they generate have been studied extensively for over half a century, yet they still remain a subject of investigation. In this section, a description of the mechanisms by which displacement damage is produced will be followed by a summary of the major consequences for device performance in a space environment. Often the degradation of a device parameter can be characterized by a damage factor (measured in a laboratory using monoenergetic protons) that is simply the change in a particular electrical or optical parameter per unit proton fluence. In addition, we will describe the concept of a non-ionizing energy loss rate (NIEL) which quantifies that portion of the energy lost by an incident ion that goes into displacements. It has been calculated as a function of proton energy, and is analogous to (and has the same units as) the linear energy transfer (LET) for ionizing energy. We will discover that, to first order, the calculated NIEL describes the energy dependence of the measured device damage factors. This observation provides the basis for predicting proton induced device degradation in a space environment based on both the calculated NIEL and relatively few laboratory test measurements. The methodology of such on-orbit device performance predictions will be described, as well as the limitations. Several classes of devices for which displacement damage is a significant (if not the dominant) mode of radiation induced degradation will be presented.

Marshall, Cheryl J.↗

Computational Assessment of a 3-Stage Axial Compressor Which Provides Airflow to the NASA 11- by 11-Foot Transonic Wind Tunnel, Including Design Changes for Increased Performance

A 24 foot diameter 3-stage axial compressor powered by variable-speed induction motors provides the airflow in the closed-return 11- by 11-Foot Transonic Wind Tunnel (11-Foot TWT) Facility at NASA Ames Research Center at Moffett Field, California. The facility is part of the Unitary Plan Wind Tunnel, which was completed in 1955. Since then, upgrades made to the 11-Foot TWT such as flow conditioning devices and instrumentation have increased blockage and pressure loss in the tunnel, somewhat reducing the peak Mach number capability of the test section. Due to erosion effects on the existing aluminum alloy rotor blades, fabrication of new steel rotor blades is planned. This presents an opportunity to increase the Mach number capability of the tunnel by redesigning the compressor for increased pressure ratio. Challenging design constraints exist for any proposed design, demanding the use of the existing driveline, rotor disks, stator vanes, and hub and casing flow paths, so as to minimize cost and installation time. The current effort was undertaken to characterize the performance of the existing compressor design using available design tools and computational fluid dynamics (CFD) codes and subsequently recommend a new compressor design to achieve higher pressure ratio, which directly correlates with increased test section Mach number. The constant cross-sectional area of the compressor leads to highly diffusion factors, which presents a challenge in simulating the existing design. The CFD code APNASA was used to simulate the aerodynamic performance of the existing compressor. The simulations were compared to performance predictions from the HT0300 turbomachinery design and analysis code, and to compressor performance data taken during a 1997 facility test. It was found that the CFD simulations were sensitive to endwall leakages associated with stator buttons, and to a lesser degree, under-stator-platform flow recirculation at the hub. When stator button leakages were modeled, pumping capability increased by over 20 of pressure rise at design point due to a large reduction in aerodynamic blockage at the hub. Incorporating the stator button leakages was crucial to matching test data. Under-stator-platform flow recirculation was thought to be large due to a lack of seals. The effect of this recirculation was assessed with APNASA simulations recirculating 0.5, 1, and 2 of inlet flow about stators 1 and 2, modeled as axisymmetric mass flux boundary conditions on the hub before and after the vanes. The injection of flow ahead of the stators tended to re-energize the boundary layer and reduce hub separations, resulting in about 3 increased stall margin per 1 of inlet flow recirculated. In order to assess the value of the flow recirculation, a mixing plane simulation of the compressor which gridded the under-stator cavities was generated using the ADPAC CFD code. This simulation indicated that about 0.65 of the inlet flow is recirculated around each shrouded stator. This collective information was applied during the redesign of the compressor. A potential design was identified using HT0300 which improved overall pressure ratio by removing pre-swirl into rotor 1, replacing existing NASA 65 series rotors with double circular arc sections, and re-staggering rotors and the existing stators. The performance of the new design predicted by APNASA and HT0300 is compared to the existing design.

Turbomachinery↗

WEC-Sim Modeling of Laminar Scientific's Patented Seesaw Wave Energy Converter: Preprint

Laminar Scientific's patented seesaw wave energy converter was modeled in WEC-Sim to predict performance. The device operates by utilizing ocean surface waves to rotate a truss in pitch about a pivot. The pivot is located at the top of two pylons, which are embedded in the seafloor. The seesaw has a float on either end, and the buoyancy forces from each float cause the system to rise or fall with passing waves. Device performance relies upon seesaw length and ocean wavelength creating an antiphase effect. The seesaw truss has an adjustable length intended to achieve this effect. The operation method enforces a narrow band of wavelengths which induce the largest rotational motion from the device. The hydrodynamic analysis of the device was performed using Capytaine, and the results confirmed that the device operates best in a narrow frequency band. Four float-to-float spacing cases and three pylon radii were examined. The hydrodynamic results indicate a match between the model and the physical expectations for the device, and that varying the pylon radii by 0.1-m increments for three instances creates minimal changes in hydrodynamic properties. Power matrices for three float spacing cases of the device were simulated with Joint North Sea Wave Project spectra waves and optimal power take-off damping in WEC-Sim. The maximum average power production for the 15-m spacing case was 14.1 kW with a 5.0-s peak wave period and 4-m significant wave heights. Plots of capture-width ratios indicated that the device performance was linear and confirmed that the device is optimal in a narrow frequency band. The maximum percentage of the available wave power produced by the 15-m device was approximately 16%. Simulations of the device in regular waves were used to produce plots of average power compared to a ratio of float spacing to wavelength. These plots indicate that the power production is maximized at a ratio of 0.5, and further confirm that the device has a narrow frequency response. The device was simulated at an example field location, where the device produced an annual average power rating of 1.6 kW given an average omnidirectional wave climate of 10.3 kW m-1 and an optimal, linearized power take-off model. While the maximum predicted device performance is reliant upon a narrow band of wave frequencies, the conducted analysis provides an opportunity to improve device design prior to prototyping and testing. Modifying the design to respond to a broader frequency range would improve device performance.

numerical modeling↗

Graph neural networks for CO 2 solubility predictions in Deep Eutectic Solvents

Deep Eutectic Solvents (DESs) are a promising class of solvents for CO 2 capture. DESs are complex mixtures that can be designed to optimize CO solubility and overall capture process efficiency. However, the vast design landscape of DES mixtures makes experimental investigation prohibitive; as such, there is a need for computational models that can quickly and efficiently navigate the design space and inform data collection efforts. In this work, we propose Graph Neural Network (GNN) models for predicting CO 2 solubility for DESs; the GNN leverages a mixture graph representation that captures the molecular structure of the DES components as well as their intermolecular interactions. Here, we compare the GNN framework against alternative architectures (neural networks, graph convolution networks, and random forests) and data representations (molecular fingerprints, sigma profiles, and graphs). We show that the proposed approach offers superior predictive performance; specifically, we show that solubility can be predicted reliably directly from molecular structure (without the need of using sigma profiles as proposed in previous studies). This result is important, as obtaining sigma profiles requires expensive density functional theory computations. We also explored the ability of GNNs to predict solubility for new DES mixtures and operating conditions. We found that the model extrapolates across temperature reliably. However, we also found deficiencies in the ability of the model to predict solubility for DES mixtures, pressures, and molar ratio not included in the training sets; we show that this is due to an inherent lack of chemical diversity in datasets available in the literature. The proposed computational capabilities can thus help navigate the design space of DES and inform data collection efforts. Our models, data, and benchmarks are shared as Python code implemented in Jupyter notebooks.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Regional Assessment for the CO 2 Storage Potential in the Northern Niagaran Pinnacle Reef Trend

The Midwest Regional Carbon Sequestration Partnership (MRCSP) was founded in 2003 as part of the U.S. Department of Energy’s (DOE’s) Regional Carbon Sequestration Partnership initiative. Since its founding, MRCSP has made significant strides toward making CCUS a viable option for states in the region. The public/private consortium, funded through the DOE Regional Carbon Sequestration Initiative, brings together nearly 40 industry partners and 10 states. Battelle, as the project lead, oversees research, development and operations and coordinates activities among the partners. The incremental, phased approach has built a valuable knowledge base for the industry and paved the way for commercial-scale adoption of CCUS technologies. From 2008 to 2020, MRCSP Phase III focused on the development of large-scale injection projects. This report is part of a series of reports prepared under the Midwestern Regional Carbon Sequestration Partnership (MRCSP) Phase III (Development Phase). These reports summarize and detail the findings of the work conducted under the Phase III project. The report describes the regional analysis of the Northern Niagaran Pinnacle Reef Trend (NNPRT) in Northern Michigan as part of the MRCSP programs Regional Assessment goals. The goal of this task was to perform an initial assessment of the geologic storage capacity and injectivity of the Northern Fairway of Michigan’s Niagaran Reef Trend. This was achieved through three main objectives: 1) development of a reef atlas, 2) understanding the regional trends in geology, and 3) estimating CO 2 and CO 2 -EOR resources. Data collected across multiple subtasks were integrated into a mappable database that could be accessed from many of the leading software. The reef atlas was used to map the reefs and associated characteristics across the entire NNPRT. Several methodologies were explored to estimate the CO 2 and CO 2 -EOR resources across the trend which included three scenarios: 1) storage only using fluid substitution and volumetric estimations, 2) CO 2 -EOR which applied proximity analysis concepts and measured performance metrics from the Core Energy reefs to predict performance at all oil reefs, and 3) enhanced storage scenario which combined CO 2 -EOR with fluid substitution to represent maximized storage after completion of CO 2 -EOR.

01 COAL, LIGNITE, AND PEAT↗

Deep Learning Predicts Stress–Strain Relations of Granular Materials Based on Triaxial Testing Data

This study presents an AI-based constitutive modelling framework wherein the prediction model directly learns from triaxial testing data by combining discrete element modelling (DEM) and deep learning. A constitutive learning strategy is proposed based on the generally accepted frame-indifference assumption in constructing material constitutive models. The low-dimensional principal stress-strain sequence pairs, measured from discrete element modelling of triaxial testing, are used to train recurrent neural networks, and then the predicted principal stress sequence is augmented to other high-dimensional or general stress tensor via coordinate transformation. Through detailed hyperparameter investigations, it is found that long short-term memory (LSTM) and gated recurrent unit (GRU) networks have similar prediction performance in constitutive modelling problems, and both satisfactorily predict the stress responses of granular materials subjected to a given unseen strain path. Furthermore, the unique merits and ongoing challenges of data-driven constitutive models for granular materials are discussed.

42 ENGINEERING↗