Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Dynamic model reduction”

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 253 records · Page 14

Sensitivity of Carbonaceous Aerosol Properties to the Implementation of a Dynamic Aging Parameterization in the Regional Climate Model RegCM

Freshly emitted soot is hydrophobic, but condensation of secondary aerosols and coagulation with other particles modify its hygroscopic optical properties. This conversion is referred to as “aerosol aging.” Many climate models represent this aging process with a fixed aging time scale, whereas in reality, it is a dynamic process that depends on environmental conditions. Here, we implement a dynamic aging parameterization scheme in the regional climate model RegCM4 in place of the fixed aging timescale of 1.15 days (~27.6 h) and examine its impact on the aerosol life cycle over the Indian subcontinent. The conversion from hydrophobic to hydrophilic aerosol is usually lower than 27.6 h over the entire landmass and lower than 10 h over the polluted Indo-Gangetic Basin (IGB), with seasonal variability. Due to the implementation of the dynamic aging scheme, the column burden and surface mass concentration of carbonaceous aerosols increase during the drier season (December–February) when washout is negligible. The burden is reduced during the wet season (June–September) due to a more efficient washout except over the IGB, where a reduction in precipitation as a result of radiative feedbacks increases the aerosol concentrations. Over the polluted IGB, surface dimming increases due to the dynamic aging scheme, with the top of the atmosphere forcing remaining mostly unchanged. As a result, atmospheric heating increases by at least 1.2 W/m 2 . Our results suggest that climate models should incorporate dynamic aging for a more realistic representation of aerosol simulations, especially in highly polluted regions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Reductions in wind farm main bearing rating lives resulting from wake impingement

This paper studies the impacts of wake impingement on main bearing rating lives predicted during the wind turbine design stage. A computational tool chain was developed to explore and quantify these effects across a wind farm populated by 10 MW wind turbines. Wind field and turbine load modelling was undertaken using the Dynamiks Python package, including application of a dynamic wake meandering model. The ISO 281 basic bearing rating life formulation was subsequently applied in order to evaluate impacts from wake effects. Analyses included a two-turbine parametric analysis, followed by a full wind farm analysis undertaken for the TotalControl 32-turbine reference wind farm, including full wind rose simulations across all operational wind speeds. Site conditions were accounted for using a Weibull wind speed distribution and a range of parametric wind direction rose models. Results indicate that wind farm main bearing rating lives are negatively impacted by the effects of wake impingement, resulting in rating life reductions for the analysed wind farm of the order of 16 % on average and as much as 20 %–25 %, both for the locating main bearing. Despite these high sensitivities, it is important to note that these resultant rating lives (i.e. the predicted lives) still far exceed the standard wind turbine operational lifetimes of 20–30 years. Wake impacts were also found to be asymmetrically related to the side on which the rotor is impinged, suggesting that, for the main bearing, there may be a “better” side for wake impingement to occur. Rating life sensitivities to wind rose shape were also observed. While these findings must be interpreted with due consideration for the various methodological limitations present, they provide compelling evidence that wake effects at the wind farm level should necessarily be included when undertaking main bearing operational load modelling, rating life assessment, or other load-related analyses.

17 WIND ENERGY↗

Enhanced Feedstock Characterization and Modeling to Facilitate Optimal Preprocessing and Deconstruction of Corn Stover (Final Report)

This project addresses the challenge of processing corn stover by fractionating this biomass feedstock to both streamline processing and generate new potential co-products. Additionally, the project developed new field-deployable analytical tools that can be coupled with empirical models that were used to predict feedstock properties and processing performance. The overall scope of this project was: (1) identify conditions for optimal corn stover fractionation using a two- stage physical fractionation, (2) assess how physical fractionation impacts properties, partitioning of biomass, and response to processing, (3) further adapt, develop, and validate several advanced characterization tools for assessing biomass properties that can be linked to processing behavior, and (4) develop and validate predictive models based on measurements that can be performed “in the field” or “at the biorefinery gate” to predict feedstock processing behavior (preprocessing and deconstruction). The first objective employed pre-separation processing (size reduction) which was next subjected to enhanced separations to yield fractions enriched or depleted in select compositional components or properties. For the second objective, fractions were screened for their response to post-separation processing (pretreatment and enzymatic hydrolysis). Detailed characterization profiles were developed and dynamic image analysis to assess distribution of particle size and morphology. For the final objective, we utilized these tools to develop empirical models to assess the relative abundance of tissue type in order to assess fractionation efficacy and to predict fraction performance during pretreatment and enzymatic hydrolysis.

09 BIOMASS FUELS↗

Numerical Model Development of a Variable-Geometry Attenuator Wave Energy Converter

Because the wave energy industry is still in its infancy, an optimal design for wave energy converters (WECs) has yet to be established; more work is needed to explore various cost-reduction pathways. The primary cost-reduction pathway considered for this work is the optimization of the geometric profile on an attenuator WEC to maximize power production while, at the same time, minimizing capital expenditures through the use of variable-geometry modules. In this investigation, the variable-geometry modules consist of inflatable bags placed on either side of a base central steel cylinder that would be inflated in low-moderate sea states to maximize power capture and then deflated in moderate-extreme sea states to minimize wave loading. The numerical model and simulation of the attenuator WEC were developed and completed using WEC-Sim, which is an open-source code that is appropriate for use in evaluating the dynamic response of the different WEC models in operational seas. The power production estimates were obtained from the Wave Energy Prize (WEP) sea states, which are representative of U. S. deployment sites, to calculate the average climate capture width that is used in the WEP ACE calculation. Preliminary capital expenditure costs were obtained assuming the base central steel cylinder mass was equal to the fluid displaced mass, minus the mass of the variable-geometry bags. The additional weight required to offset the additional buoyancy from the variable-geometry bags was assumed to come from the addition of seawater ballast. The variable-geometry attenuator model was found to have a similar power capture efficiency as a fixed-body model, but is expected to have a lower characteristic capital expenditure given its more streamlined profile, which demonstrates that variable-geometry modules may provide a realistic cost-reduction pathway to help design a more cost-competitive WEC.

50 EE - Wind and Water Power Program - Water (EE-4↗

Toward the Prediction and Control of Glass Transition Temperature for Donor–Acceptor Polymers

Semiconducting donor–acceptor (D–A) polymers have attracted considerable attention toward the application of organic electronic and optoelectronic devices. However, a rational design rule for making semiconducting polymers with desired thermal and mechanical properties is currently lacking, which greatly limits the development of new polymers for advanced applications. Here, polydiketopyrrolopyrrole (PDPP)-based D–A polymers with varied alkyl side-chain lengths and backbone moieties are systematically designed, followed by investigating their thermal and thin film mechanical responses. The experimental results show a reduction in both elastic modulus and glass transition temperature (Tg) with increasing side-chain length, which is further verified through coarse-grained molecular dynamics simulations. Informed from experimental results, a mass-per-flexible bond model is developed to capture such observation through a linear correlation between Tg and polymer chain flexibility. Using this model, a wide range of backbone Tg over 80 °C and elastic modulus over 400 MPa can be predicted for PDPP-based polymers. This study highlights the important role of side-chain structure in influencing the thermomechanical performance of conjugated polymers, and provides an effective strategy to design and predict Tg and elastic modulus of future new D–A polymers.

36 MATERIALS SCIENCE↗

A combined theoretical-experimental investigation of thermal transport in low-dose irradiated thorium dioxide

During reactor operation, nuclear fuels are subject to extreme temperature and irradiation conditions which can significantly degrade the fuel's thermal transport properties. The reduction in thermal conductivity of the fuel as a result of irradiation-induced lattice defects is arguably the most important fuel performance metric in regard to reactor efficiency and safety. Because thorium dioxide (ThO2) is suitable as a model system for more complex materials such as UO2 and its mixed oxides, we present a theoretical investigation of thermal conductivity reduction seen in defect-bearing thorium dioxide and compare directly to experimental measurements. Phonon-mediated thermal transport of the fuel is modeled by a solution to the Boltzmann transport equation (BTE) for phonons. A cluster dynamics (CD) model for lattice defect evolution during irradiation predicts defect densities which are used as input to the BTE for modeling phonon-defect scatterings. Phonon scatterings by lattice defects include those from point defects and vacancy clusters and interstitial clusters of various sizes. The CD model is benchmarked against structural defect characterization of irradiated thorium dioxide using electron microscopy. Thermal conductivity predicted by the BTE model is compared to measured values for irradiated thorium dioxide specimens below room temperature to isolate effects of phonon-defect scattering from intrinsic 3-phonon processes, which dominate at higher temperatures. The computed conductivity values are in partial agreement at temperatures close to room temperature while slight deviations are observed at the lowest measured temperatures, suggesting that implemented phonon-defect scattering cross-section expressions may not be adequate for low temperatures. The presented work provides a necessary investigation of the influence of irradiation induced defects on fuel performance and represents a first step toward a full characterization of phonon mediated thermal transport in irradiated materials with complex defect microstructure.

36 MATERIALS SCIENCE↗

A critical review of operational strategies for the management of harmful algal blooms (HABs) in inland reservoirs

Occurrences of freshwater harmful algal blooms (HABs) are increasing on a global scale, largely in part due to increased nutrient input and changing climate patterns. While reservoir management strategies that can influence phytoplankton are known, there is no published guideline or protocol for the management of harmful algal blooms. There is a need to establish what factors are the predominant drivers of blooms, and how common reservoir management strategies specifically influence each factor. The following literature review seeks to establish the benefits and drawbacks of operational management strategies that currently exist. The main focus is altering hydrodynamic conditions (hypolimnetic withdrawals, surface flushing, pulsed inflow, artificial mixing), in order to induce environmental changes within the reservoir itself. Since excess nutrients are one of the biggest contributors to worsening bloom conditions, internal nutrient dynamics and reduction are also discussed. Additionally, we review the predominant seasonal factors (stratification, light, temperature, and wind) that affect likelihood of bloom occurrence and duration. In conclusion, the ultimate objective of this review is to increase understanding of the relationships between HAB drivers and reservoir operations in order to inform the development of data, modeling, and management strategies for the prevention and mitigation of blooms.

54 ENVIRONMENTAL SCIENCES↗

Cation–Ligand Interactions Dictate Salt Partitioning and Diffusivity in Ligand-Functionalized Polymer Membranes

Membranes are an attractive alternative to current thermal separations due to their scalability and energy efficiency in desalinating water. Unfortunately, many of the conventional membrane materials available today are unable to differentiate between ionic solutes, especially alkali cations, compromising their use in ion–ion separations. Inspired by the ion-specific interactions exhibited by biological ion channels, recent research efforts have focused on synthesizing and characterizing new polymeric materials that incorporate ligands into polymer networks to bias solubility and/or diffusivity of one cationic species over another. Despite these efforts, little is known about the influence of incorporating ligands into polymer membranes on solubility and diffusivity of the complexing species. In this study, we first build a qualitative model of salt partitioning, diffusivity, and permeability in generic cation-complexing ligand-functionalized polymer membranes. Next, to validate our model and hypotheses, we perform atomistic molecular dynamics simulations of a 12-crown-4-functionalized membrane in the presence of alkali halide salts at low concentration. Generally, cation complexation enhances cation solubility but decreases diffusivity. Interestingly, the reduction in diffusivity is predicted to be larger than the enhancement in solubility for materials which operate by the mechanisms proposed in our physical picture, ultimately resulting in a reduction in the permeability of the selectively complexing ion.

36 MATERIALS SCIENCE↗

Polariton Control of Molecular Charge Transfer in Perylene Diimide Semiconductors

We report the modulation of molecular charge transfer in a bay-substituted perylene diimide derivative embedded in a planar distributed Bragg reflector microcavity. Angle-resolved reflectance spectra confirm the formation of upper and lower polaritons with clear Rabi splitting, indicating strong coupling between the cavity mode and molecular excitons. Using broadband transient absorption spectroscopy, we compare the ultrafast dynamics of cavity and non-cavity films. While excited-state absorption and stimulated emission pathways remain largely unchanged, kinetic modeling reveals a moderate increase in the charge transfer rate and yield under strong coupling. This enhancement is attributed to a reduction in the effective driving force via the formation of the lower polariton, placing the system deeper into the Marcus inverted regime. Our results demonstrate a promising non-chemical method leveraging cavity quantum electrodynamics to modulate charge separation processes in molecular semiconductors.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Flexible Siting Criteria and Staff Minimization for Micro-Reactors

The economic potential of micro-reactors is vast and underestimated. Commonly-emphasized applications include niche markets such as remote communities, mines and military bases. However, micro-reactors could be used as flexible energy generators also for larger markets, such as mobile and containerized agriculture and manufacturing facilities, district heating, micro-grids for data centers, sea ports, airports and hospitals. The implication is that micro-reactors may have to be deployed also in non-remote locations. Successful implementation of micro-reactors needs a navigable and predictable licensing process, technology-appropriate siting restrictions, risk-informed emergency and safety requirements, and practical operating and maintenance requirements. The primary goal of this project was to develop siting criteria that are tailored to micro-reactors deployable in densely-populated areas, e.g., urban environments. To achieve that goal, we compared the characteristics of the MIT research reactor (MITR) with those of leading micro-reactor concepts (e.g., eVinci, USNC, Aurora), and evaluated whether and how the MITR design basis (e.g., inherent safety features, engineered safety systems, source term, emergency planning and emergency operating procedures) and associated regulations may be applicable to these new micro-reactors as well. What makes MITR a unique analogue in this context is its small power rating (6 MWt) and physical size, mode of operations (24/7 with a somewhat more commercial flavor than typical university reactors), and especially its urban location. Of course significant differences exist, such as mission (power production vs. research) and the reactor design itself. Leveraging the MITR experience, this project was able to generate criteria that will allow micro-reactors to realize their full economic potential as flexible heat and electricity generators for a diverse portfolio of applications in non-remote locations. As such, the outcome of this project might encourage investment in and use of micro-reactors. A second goal of the project was to conceptualize a model of operations for micro-reactors that would minimize the staffing requirements, and thus reduce the cost of electricity and heat generated by these systems. Here too our approach was to systematically review the MITR experience and requirements, as well as survey the innovations in autonomous control technologies and monitoring (e.g., advanced sensors, drones, robotics, AI) that would permit a dramatic reduction in staffing at future micro-reactor installations. The scope of work was expanded after the start date to include also an evaluation of micro-reactor security, using the so-called consequence-based analysis, and the development of a methodology to perform dynamic risk assessment for micro-reactors, using system theory and modeling and simulation.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Extending the Publish/Subscribe Abstraction for High-Performance I/O and Data Management at Extreme Scale

The Adaptable I/O System (ADIOS) represents the culmination of substantial investment in Scientific Data Management, and it has demonstrated success for several important extreme-scale science cases. However, looking towards the exascale and beyond, we see the development of yet more stringent data management requirements that require new abstractions. Therefore, there is an opportunity to attempt to connect the traditional realms of HPC I/O optimization with the Database / Data Management community. As such, in this paper we offer some specific examples from our ongoing work in managing data structures, services, and performance at the extreme scale for scientific computing. Using the publish/subscribe model afforded by ADIOS, we demonstrate a set of services that connect data format, metadata, queries, data reduction, and high-performance delivery. The resulting publish/subscribe framework facilitates connection to on-line workflow systems to enable the dynamic capabilities that will be required for exascale science.

96 KNOWLEDGE MANAGEMENT AND PRESERVATION↗

Effects of social structure and management on risk of disease establishment in wild pigs

Contact heterogeneity among hosts determines invasion and spreading dynamics of infectious disease, thus its characterization is essential for identifying effective disease control strategies. Yet, little is known about the factors shaping contact networks in many 28 wildlife species and how wildlife management actions might affect contact networks. Wild pigs in North America are an invasive, socially-structured species that pose a health concern for domestic swine given their ability to transmit numerous devastating diseases 2 such as African swine fever (ASF). Using proximity loggers and GPS data from 48 wild pigs in Florida and South Carolina, USA, we employed a probabilistic framework to estimate weighted contact networks. We determined the effects of sex, social group, and spatial distribution (monthly home range overlap and distance) on wild pig contact. We also estimated the impacts of management-induced perturbations on contact and inferred their effects on ASF establishment in wild pigs with simulation. Social group membership was the primary factor influencing contacts. Between-group contacts depended primarily on space use characteristics, with fewer contacts among groups separated by >2 km and no contacts among groups >4 km apart within a month. Modeling ASF dynamics on the contact network demonstrated that indirect contacts resulting from baiting (a typical method of attracting wild pigs or game species to a site to enhance recreational hunting) increased the risk of disease establishment by ~33% relative to direct contact. Low-intensity population reduction (<5.9% of the population) had no detectable impact on contact structure but reduced predicted ASF establishment risk relative to no population reduction. Here, we demonstrate an approach for understanding the relative role of spatial, social, and individual-level characteristics in shaping contact networks and predicting their effects on disease establishment risk, thus providing insight for optimizing disease control in spatially- and socially-structured wildlife species.

59 BASIC BIOLOGICAL SCIENCES↗

Strategies for solving the excited-state self-consistent-field problem for highly excited and multiply ionized states

The dynamics of molecules exposed to intense x-ray radiation involve a large number of multiply ionized and highly excited electronic configurations. To model these dynamics a reliable and efficient electronic structure model is imperative. Employing the Hartree-Fock-Slater electronic structure model in combination with the maximum overlap method, we quantify the associated convergence failures when calculating electronic states of carbon monoxide with multiple vacancies in the core and valence levels. Here, we characterize these cases and describe strategies to overcome the convergence problems. The described techniques not only eliminate all convergence issues for CO but also result in a significant reduction of convergence failures for simulations of the x-ray-induced multiple ionization dynamics of the phenol molecule.

74 ATOMIC AND MOLECULAR PHYSICS↗

Co-Active Subspace Methods for the Joint Analysis of Adjacent Computer Models

Active subspace (AS) methods are a valuable tool for understanding the relationship between the inputs and outputs of a Physics simulation. In this article, an elegant generalization of the traditional ASM is developed to assess the co-activity of two computer models. This generalization, which we refer to as a Co-Active Subspace (Co-AS) Method, allows for the joint analysis of two or more computer models allowing for thorough exploration of the alignment (or non-alignment) of the respective gradient spaces. We define co-active directions, co-sensitivity indices, and a scalar “concordance” metric (and complementary “discordance” pseudo-metric) and we demonstrate that these are powerful tools for understanding the behavior of a class of computer models, especially when used to supplement traditional AS analysis. Details for efficient estimation of the Co-AS and an accompanying R package (concordance) are provided. Practical application is demonstrated through analyzing a set of simulated rate stick experiments for PBX 9501, a high explosive, offering insights into complex model dynamics.

97 MATHEMATICS AND COMPUTING↗

Detailed biomass fast pyrolysis kinetics integrated to computational fluid dynamic (CFD) and discrete element modeling framework: Predicting product yields at the bench-scale

Fast pyrolysis is an intricate process due to the variability and anisotropy of lignocellulosic biomass and the complicated chemistry and physics during conversion in a bubbling fluidized bed reactor (BFBR). The complexity of biomass fast pyrolysis lends itself well to computational fluid dynamics (CFD) and discrete element (DEM) analysis, which promises to reduce experimental time and its associated cost. This work investigated switchgrass fast pyrolysis simulated by computational fluid dynamics coupled with a discrete element method to track individual reacting biomass particles throughout a bench-scale BFBR reactor. We accounted for the fast pyrolysis chemistry through a comprehensive reaction scheme with secondary cracking reactions. We performed a three-step reduction for secondary cracking reactions to convert the full cracking scheme into a reduced scheme easily incorporated into our model. We assessed the impact of operational conditions on the steady-state yields of liquid bio-oil, non-condensable gases (NCG), at 550 °C over a range of fluidization numbers (2 – 6 Umf), reported as a ratio to the minimum fluidization velocity (Umf). At steady-state, the volatile bio-oil yield had a range of 49.3–50.4 wt%. Levoglucosan was the primary volatile component present with 21 wt% of the bio-oil while water was the second largest with 20 wt%. The reduction of the secondary reaction schemes did not appreciably affect the overall yields of switchgrass pyrolysis compared to the full secondary scheme.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

CD8 + T cells control SIV infection using both cytolytic effects and non-cytolytic suppression of virus production

Whether CD8 + T lymphocytes control human immunodeficiency virus infection by cytopathic or non-cytopathic mechanisms is not fully understood. Multiple studies highlighted non-cytopathic effects, but one hypothesis is that cytopathic effects of CD8 + T cells occur before viral production. Here, to examine the role of CD8 + T cells prior to virus production, we treated SIVmac251-infected macaques with an integrase inhibitor combined with a CD8-depleting antibody, or with either reagent alone. We analyzed the ensuing viral dynamics using a mathematical model that included infected cells pre- and post- viral DNA integration to compare different immune effector mechanisms. Macaques receiving the integrase inhibitor alone experienced greater viral load decays, reaching lower nadirs on treatment, than those treated also with the CD8 - depleting antibody. Models including CD8 + cell-mediated reduction of viral production (non-cytolytic) were found to best explain the viral profiles across all macaques, in addition an effect in killing infected cells pre-integration (cytolytic) was supported in some of the best models. Our results suggest that CD8 + T cells have both a cytolytic effect on infected cells before viral integration, and a direct, non-cytolytic effect by suppressing viral production.

59 BASIC BIOLOGICAL SCIENCES↗

Efficient learning of power grid voltage control strategies via model-based deep reinforcement learning

Here this article proposes a model-based deep reinforcement learning (DRL) method to design emergency control strategies for short-term voltage stability problems in power systems. Recent advances show promising results for model-free DRL-based methods in power systems control problems. But in power systems applications, these model-free methods have certain issues related to training time (clock time) and sample efficiency; both are critical for making state-of-the-art DRL algorithms practically applicable. DRL-agent learns an optimal policy via a trial-and-error method while interacting with the real-world environment. It is also desirable to minimize the direct interaction of the DRL agent with the real-world power grid due to its safety-critical nature. Additionally, the state-of-the-art DRL-based policies are mostly trained using a physics-based grid simulator where dynamic simulation is computationally intensive, lowering the training efficiency. We propose a novel model-based DRL framework where a deep neural network (DNN)-based dynamic surrogate model (SM), instead of a real-world power grid or physics-based simulation, is utilized within the policy learning framework, making the process faster and more sample efficient. However, having stable training in model-based DRL is challenging because of the complex system dynamics of large-scale power systems. We addressed these issues by incorporating imitation learning to have a warm start in policy learning, reward-shaping, and multi-step loss in surrogate model training. Finally, we achieved 97.5% reduction in samples and 87.7% reduction in training time for an application to the IEEE 300-bus test system.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Calibration and validation of the foundation for a multiphase strength model for tin

In this work, the Common Model of Multi-phase Strength and Equation of State (CMMP) model was applied to tin. Specifically, calibrations of the strength-specific elements of the CMMP foundation were developed with a combination of experiments and theory, and then the model was validated experimentally. The first element of the foundation is a multi-phase analytic treatment of the melt temperature and the shear modulus for the solid phases. These models were parameterized for each phase based on ab initio calculations using the software VASP (Vienna Ab initio Simulations Package) based on density functional theory. The shear modulus model for the ambient phase was validated with ultrasonic sound speed measurements as a function of pressure and temperature. The second element of the foundation is a viscoplastic strength model for the phase, upon which strength for inaccessible higher-pressure phases can be scaled as necessary. The stress–strain response of tin was measured at strain rates of 10 -3 to 3 x 10 3 s -1 and temperatures ranging from 87 to 373 K. The Preston–Tonks–Wallace (PTW) strength model was fit to that data using Bayesian model calibration. For validation, six forward and two reverse Taylor impact experiments were performed at different velocities to measure large plastic deformation of tin at strain rates up to ⁠10 5 s -1 . The PTW model accurately predicted the deformed shapes of the cylinders, with modest discrepancies attributed to the inability of PTW to capture the effects of twinning and dynamic recrystallization. Some material in the simulations of higher velocity Taylor cylinders reached the melting temperature, thus testing the multiphase model because of the presence of a second phase, the liquid. In simulations using a traditional modeling approach, the abrupt reduction of strength upon melt resulted in poor predictions of the deformed shape and non-physical temperatures. With CMMP, the most deformed material points evolved gradually to a mixed solid–liquid but never a fully liquid state, never fully lost strength, stayed at the melt temperature as the latent heat of fusion was absorbed, and predicted the deformed shape well.

36 MATERIALS SCIENCE↗