Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Internal energy”

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Developing a Composite Vacuum Insulated Panel (VIP) Insulation/Vinyl Siding Composite Technology for Retrofitting Residential Walls

This paper reports on a project that has attempted to seize on an opportunity to take a giant leap forward in commercializing a technology that would address a key opportunity to achieve energy savings in the existing residential housing stock. It is typically very difficult to increase the R-value of the walls in existing residential housing as doing so from the inside is disruptive and there have not been good solutions to increase thermal envelope performance from the exterior through adding continuous insulation due to the additional work required to fit existing architectural features to the new wall thickness and to reclad the façade. A high thermal resistivity or R-value per inch product that can be used in these applications offers an excellent solution.This paper summarizes a three-year research project to a produce vacuum insulated panel (VIP)/vinyl siding that would have sufficient R-value to meet the continuous insulation requirements of the 2018 International Energy Conservation Code (IECC, 2018) in all climate zones and has a much thinner profile that will facilitate its application to existing residential homes without the need for expensive retrimming of the architectural details. The paper will supply information on the development of the technology, its thermal modelling, testing, voice of user sessions to solicit input from architects, designers, manufacturers, contractors, and installers, and a techno-economic analysis to gauge its competitiveness with existing product options.

Desjarlais, Andre Omer↗
DOE OSTI2022

Developing a Composite Vacuum Insulated Panel (VIP) Insulation/Vinyl Siding Composite Technology for Retrofitting Residential Walls

This paper reports on a project that has attempted to seize on an opportunity to take a giant leap forward in commercializing a technology that would address a key opportunity to achieve energy savings in the existing residential housing stock. It is typically very difficult to increase the R-value of the walls in existing residential housing as doing so from the inside is disruptive and there have not been good solutions to increase thermal envelope performance from the exterior through adding continuous insulation due to the additional work required to fit existing architectural features to the new wall thickness and to reclad the facade. A high thermal resistivity or R-value per inch product that can be used in these applications offers an excellent solution. This paper summarizes a three-year research project to a produce vacuum insulated panel (VIP)/vinyl siding that would have sufficient R-value to meet the continuous insulation requirements of the 2018 International Energy Conservation Code (IECC, 2018) in all climate zones and has a much thinner profile that will facilitate its application to existing residential homes without the need for expensive retrimming of the architectural details. The paper will supply information on the development of the technology, its thermal modelling, testing, voice of user sessions to solicit input from architects, designers, manufacturers, contractors, and installers, and a techno-economic analysis to gauge its competitiveness with existing product options.

Desjarlais, Andre↗
DOE OSTI2026

Thermodynamics and its prediction and CALPHAD modeling: Review, state of the art, and perspectives

Thermodynamics is a science concerning the state of a system, whether it is stable, metastable, or unstable, when interacting with its surroundings. The combined law of thermodynamics derived by Gibbs about 150 years ago laid the foundation of thermodynamics. In Gibbs combined law, the entropy production due to internal processes was not included, and the 2nd law was thus practically removed from the Gibbs combined law, so it is only applicable to systems under equilibrium, thus commonly termed as equilibrium or Gibbs thermodynamics. Gibbs further derived the classical statistical thermodynamics in terms of the probability of configurations in a system in the later 1800's and early 1900's. With the quantum mechanics (QM) developed in 1920's, the QM-based statistical thermodynamics was established and connected to classical statistical thermodynamics at the classical limit as shown by Landau in the 1940's. In 1960's the development of density functional theory (DFT) by Kohn and co-workers enabled the QM prediction of properties of the ground state of a system. On the other hand, the entropy production due to internal processes in non-equilibrium systems was studied separately by Onsager in 1930's and Prigogine and co-workers in the 1950's. In 1960's to 1970's the digitization of thermodynamics was developed by Kaufman in the framework of the CALculation of PHAse Diagrams (CALPHAD) modeling of individual phases with internal degrees of freedom. CALPHAD modeling of thermodynamics and atomic transport properties has enabled computational design of complex materials in the last 50 years. Our recently termed zentropy theory integrates DFT and statistical mechanics through the replacement of the internal energy of each individual configuration by its DFT-predicted free energy. The zentropy theory is capable of accurately predicting the free energy of individual phases, transition temperatures and properties of magnetic and ferroelectric materials with free energies of individual configurations solely from DFT-based calculations and without fitting parameters, and is being tested for other phenomena including superconductivity, quantum criticality, and black holes. Those predictions include the singularity at critical points with divergence of physical properties, negative thermal expansion, and the strongly correlated physics. Furthermore, those individual configurations may thus be considered as the genomic building blocks of individual phases in the spirit of the materials genome®. This has the potential to shift the paradigm of CALPHAD modeling from being heavily dependent on experimental inputs to becoming fully predictive with inputs solely from DFT-based calculations and machine learning models built on those calculations and existing experimental data through newly developed and future open-source tools. Furthermore, through the combined law of thermodynamics including the internal entropy production, it is shown that the kinetic coefficient matrix of independent internal processes is diagonal with respect to the conjugate potentials in the combined law, and the cross phenomena that the phenomenological Onsager flux and reciprocal relationships are due to the dependence of the conjugate potential of a molar quantity on nonconjugate molar quantities and other potentials, which can be predicted by the zentropy theory and CALPHAD modeling.

42 ENGINEERING↗
DOE OSTI2026

Microcanonical Kinetics of Water-Mediated Proton Transfer in 4ABAH + ·(H 2 O) n = 4–6 Clusters (ABA = Aminobenzoic Acid): A Model System for Size-Dependent Relaxation to Ergodic Behavior

Here, we leverage the unique properties of the 4ABAH + · (H 2 O) n clusters (ABA = 4-aminobenzoic acid, n = 4−6) to quantitatively address how a finite, isolated system evolves into an ergodic condition starting from localized arrangements in configuration space. This system adopts two distinct structural isomers in which water molecules cluster around the cationic centers of its two protomers with widely separated positive charge centers. These isomers arise from excess proton attachment to either the acid (O) or amino (N) group on opposite sides of the benzene ring. Both forms are captured and kinetically trapped using cryogenic ion methods and then selectively vibrationally excited through their mutually exclusive IR bands involving NH and OH stretching fundamentals. Because the IR excitation lies below the water binding energy, the system can evolve to explore slow, rare events that lead to the interconversion between the two isomers. The rates of these intracluster reactions are determined by using a pump−probe scheme involving ∼5 ns IR pump and UV probe lasers. The rates occur on the microsecond time scale, leading to steady state populations of the isomers, thus revealing the cluster size-dependent fractionation between the two species at microcanonical equilibrium. The steady state distributions are correlated with the expected trend in the cluster size-dependent reaction energetics, which are in turn consistent with changes in the relative densities of states of the two species. These results thus provide an unusually clear example in which complex, protic-solvent-mediated chemical transformations are captured within a finite system at a precisely determined internal energy.

Rana, Abhijit [Yale Univ., New Haven, CT (United S↗
DOE OSTI2023

Probabilidad de que los murcielagos sufran barotraumatismo cerca de palas de aerogeneradores en movimiento (Spanish)

In October 2018, the International Energy Agency Wind Task 34 - Working Together to Resolve the Environmental Effects of Wind Energy (WREN) - organized a virtual forum to discuss the likelihood of bats experiencing barotrauma when flying near moving wind turbine blades. The forum included experts in bat biology and physiology, bat and wind turbine interactions, wind turbine technology, and atmospheric sciences. This educational brief summarizes the discussion during the forum and written comments from those who could not attend. Relevant literature was used to provide additional context when needed. Possible explanations regarding the direct cause of bat mortality at operating wind turbines are (1) blunt force trauma caused by turbine blades striking individual bats, often referred to as collisions, and (2) barotrauma, resulting from exposure to pressure changes located near the surface of moving wind turbine blades. While collision- related mortality is easily understood, the mechanism causing barotrauma is more complex. Fast-moving wind turbine blades create regions of high- and low-pressure variations along the blade surfaces. If bats fly within these regions, the rapid change in pressure may cause internal bleeding, damage to lungs or other organs, and damage to the inner ear. However, sufficient data demonstrating barotrauma as a common cause of bat mortality at wind turbines are lacking. Moreover, the pressure variation necessary to cause barotrauma in bats is so close to the surface of turbine blades that there is a higher probability of direct contact with the blades than of solely experiencing barotrauma. Regardless of the cause, bats are interacting with fast-moving wind turbine blades, and these interactions are resulting in fatalities. Resolving this issue will require a better understanding of bat behavior and cost-effective measures to reduce interactions between bats and wind turbines. This is the Spanish translation of NREL/FS-5000-84749, "The Likelihood of Bats Experiencing Barotrauma Near Moving Wind Turbine Blades."

barotrauma↗
DOE OSTI2026

Closure models for the feedback of energetic particles on plasma turbulence

Energetic particles interact with the plasma surrounding them, resonating with certain types of plasma waves to stabilize them while destabilizing others, and changing the character of the background turbulence in ways that have not been fully quantified or understood. Interaction with the turbulent background plasma is key to the acceleration of many types of energetic particles including high-energy cosmic rays, solar energetic particles, and pick-up ions. The acceleration of particles is a process that would ideally be described by a kinetic model, a type of model that follows a probability distribution function (PDF) for all particles in 7-dimensional (x, y, z, v x , v y , v z , t) space. Because of the high dimensionality of a kinetic model, simulations that solve kinetic equations use the largest computational resources currently available, and are yet unable to simulate a realistic number of particles, reach the large scales necessary for astrophysical problems, and use high-precision numerical methods. Two available alternatives to kinetic plasma models have been explored for this problem, with limited success. One is a multi-fluid model produced by a cumulant discarding closure, which evolves coupled equations for the velocity, magnetic field, and internal energy for both the background plasma and the fluid of energetic particles. However, simulations that solve multi-fluid magnetohydrodynamic (MHD) equations are able to include the interaction with energetic particles only in crude ways, typically as an add-on pressure term. The second alternative is to use a hybrid method to couple a fluid description of the background plasma to a kinetic model or a Fokker–Planck model for the energetic particles. These methods are hampered by the physical modeling of the coupling. In this work, we develop a new model, which follows the PDF for all particles; this can be viewed as a step toward physical realism above a multi-fluid MHD model, while also being more computationally efficient than a kinetic model. The equations we develop model both the background plasma and the energetic particles self-consistently. Over the last decade, similar PDF methods have been developed to a high level of sophistication to model reactive flows and turbulent combustion for engineering applications. For treatment of the feedback of the energetic particles on a background plasma, a PDF closure approach should evaluate the mean characteristics, including the density, with better statistical quality than will particle-sampling procedures.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

A New Shape Optimization Approach for Lightweighting Electric Machines Inspired by Additive Manufacturing

Minimizing the mass in electric machines while maintaining superior performance is a new requirement for the advancement of drivetrains used in wind energy and electric mobility. Topology optimization (TO) for lightweighting electric machines using traditional approaches typically explores a restricted design space allowed by standard parametrizable geometry and manufacturing, while advanced methods, such as cell-based density approaches, suffer from a lack of robust manufacturability constraints during the optimization process. To overcome these drawbacks, we explore a grid-independent, boundary optimization where the outer shape of the magnet is parameterized using Bezier curves. We conduct a design of experiments (DOE) to study the effect of different magnet shapes on machine performance by varying the control points on the Bezier curves. A machine-learning-based surrogate model is constructed using the data from the DOE to quantify the relationship between the control points, air-gap torque, and mass. The control points are then optimized to maximize the torque density. The approach is used for minimizing electrical steel mass in the International Energy Agency (IEA) 15-MW radial flux direct-drive wind turbine generator. The new approach to shape optimization resulted in smooth and concise shapes that can be easily additively manufactured with up to a 20-ton reduction in electrical steel mass.

additive manufacturing↗
DOE OSTI2022

A New Shape Optimization Approach for Lightweighting Electric Machines Inspired by Additive Manufacturing: Preprint

Minimizing the mass in electric machines while maintaining superior performance is a new requirement for the advancement of drivetrains used in wind energy and electric mobility. Topology optimization (TO) for lightweighting electric machines using traditional approaches typically explores a restricted design space allowed by standard parametrizable geometry and manufacturing, while advanced methods, such as cell-based density approaches, suffer from a lack of robust manufacturability constraints during the optimization process. To overcome these drawbacks, we explore a grid-independent, boundary TO where the outer shape of the magnet is parameterized using Bezier curves. We conduct a design of experiments (DOE) to study the effect of different magnet shapes on machine performance by varying the control points on the Bezier curves. A machine-learning-based surrogate model is constructed using the data from the DOE to quantify the relationship between the control points, air-gap torque, and mass. The control points are then optimized to maximize the torque density. The approach is used for minimizing electrical steel mass in the International Energy Agency (IEA) 15-MW radial flux direct-drive wind turbine generator. The free-form boundary TO resulted in smooth and concise shapes that can be easily additively manufactured with up to a 20-ton reduction in electrical steel mass.

additive manufacturing↗
DOE OSTI2021

Best Practices Handbook for the Collection and Use of Solar Resource Data for Solar Energy Applications: Third Edition

As the world looks for low-carbon sources of energy, solar power stands out as the single most abundant energy resource on Earth. Harnessing this energy is the challenge for this century. Photovoltaics, solar heating and cooling, and concentrating solar power (CSP) are primary forms of energy applications using sunlight. These solar energy systems use different technologies, collect different fractions of the solar resource, and have different siting requirements and production capabilities. Reliable information about the solar resource is required for every solar energy application. This holds true for small installations on a rooftop as well as for large solar power plants; however, solar resource information is of particular interest for large installations because they require substantial investment, sometimes exceeding 1 billion dollars in construction costs. Before such a project is undertaken, the best possible information about the quality and reliability of the fuel source must be made available. That is, project developers need reliable data about the solar resource available at specific locations, including historic trends with seasonal, daily, hourly, and (preferably) subhourly variability to predict the daily and annual performance of a proposed power plant. Without these data, an accurate financial analysis is not possible. Additionally, with the deployment of large amounts of distributed photovoltaics, there is an urgent need to integrate this source of generation to ensure the reliability and stability of the grid. Forecasting generation from the various sources will allow for larger penetrations of these generation sources because utilities and system operators can then ensure stable grid operations. Developed by the foremost experts in the field who have come together under the umbrella of the International Energy Agency’s Solar Heating and Cooling Task 46, this handbook summarizes state-of-the-art information about all these topics.

14 SOLAR ENERGY↗
DOE OSTI2022

4.2.1.31 Integrated Life Cycle Sustainability Analysis

This project provides the Department of Energy's Bioenergy Technologies Office (BETO) with strategic decision-support for the evaluation of its R&D portfolio by developing, validating, and applying a coherent methodology and consistent model framework to quantify the net effects of an expanding US bioeconomy. The framework fills an analysis gap previously identified by Peer Review and supports a related milestone in BETO's Multi-Year Program Plan. The framework was scoped with inputs from practitioners in academia, national laboratories, and federal agencies. The model is a top-down, economy-wide framework using a coherent methodology to compute environmental and socio-economic metrics. It is purposefully complementary to existing bottom-up, process-based techno-economic and life cycle assessment BETO tools and uses their data as inputs. Presently, the model covers several commercial and near-commercial biofuel routes and an emerging pathway for plastics upcycling. It covers temporal detail across four time-steps and is currently being expanded with a prospective modeling capability. The model has provided analyses for the Third Triennial Report to Congress (RtC3) on the environmental impacts of the Renewable Fuel Standard (RFS2), among others. As part of this project, NREL also provides scientific support to BETO in the International Energy Agency's Technology Collaboration Program on Bioenergy (IEA Bioenergy) Task 45 on Sustainability. Here, NREL evaluates and synthesizes activities that develop, compare, or apply metrics, methods, and tools to quantify sustainability effects of bioeconomy products. NREL also coordinates related national lab involvement and a BETO Working Group on Sustainable Land Management.

bioeconomy↗
DOE OSTI2023

Dual Tunability for Uncatalyzed N‐Alkylation of Primary Amines Enabled by Plasma‐Microdroplet Fusion

Abstract The fusion of non‐thermal plasma with charged microdroplets facilitates catalyst‐free N‐alkylation for a variety of primary amines, without halide salt biproduct generation. Significant reaction enhancement (up to >200×) is observed over microdroplet reactions generated from electrospray. This enhancement for the plasma‐microdroplet system is attributed to the combined effects of energetic collisions and the presence of reactive oxygen species (ROS). The ROS (e.g., O 2 ⋅ − ) act as a proton sink to increase abundance of free neutral amines in the charged microdroplet environment. The effect of ROS on N‐alkylation is confirmed through three unique experiments: (i) utilization of radical scavenging reagent, (ii) characterization of internal energy distribution, and (iii) controls performed without plasma, which lacked reaction acceleration. Establishing plasma discharge in the wake of charged microdroplets as a green synthetic methodology overcomes two major challenges within conventional gas‐phase plasma chemistry, including the lack of selectivity and product scale‐up. Both limitations are overcome here, where dual tunability is achieved by controlling reagent concentration and residence time in the microdroplet environment, affording single or double N‐alkylated products. Products are readily collected yielding milligram quantities in eight hours. These results showcase a novel synthetic strategy that represents a straightforward and sustainable C−N bond‐forming process.

Chemistry↗
DOE OSTI2024

A multiscale anisotropic polymer network model coupled with phase field fracture

The study of polymers has continued to gain substantial attention due to their expanding range of applications, spanning essential engineering fields to emerging domains like stretchable electronics, soft robotics, and implantable sensors. These materials exhibit remarkable properties, primarily stemming from their intricate polymer chain network, which, in turn, increases the complexity of precisely modeling their behavior. Especially for modeling elastomers and their fracture behavior, accurately accounting for the deformations of the polymer chains is vital for predicting the rupture in highly stretched chains. Despite the importance, many robust multiscale continuum frameworks for modeling elastomer fracture tend to simplify network deformations by assuming uniform behavior among chains in all directions. Recognizing this limitation, our study proposes a multiscale fracture model that accounts for the anisotropic nature of elastomer network responses. At the microscale, damage in the chains is assumed to be driven by both the chain's entropy and the internal energy due to molecular bond distortions. In order to bridge the stretching in the chains to the macroscale deformation, we employ the maximal advance path constraint network model, inherently accommodating anisotropic network responses. As a result, chains oriented differently can be predicted to exhibit varying stretch and, consequently, different damage levels. To drive macroscale fracture based on damages in these chains, we utilize the micromorphic regularization theory, which involves the introduction of dual local-global damage variables at the macroscale. The macroscale local damage variable is obtained through the homogenization of the chain damage values, resulting in the prediction of an isotropic material response. The macroscale global damage variable is subjected to nonlocal effects and boundary conditions in a thermodynamically consistent phase field continuum formulation. Moreover, the total dissipation in the system is considered to be mainly due to the breaking of the molecular bonds at the microscale. To validate our model, we employ the double-edge notched tensile test as a benchmark, comparing simulation predictions with existing experimental data. Additionally, to enhance our understanding of the fracturing process, we conduct uniaxial tensile experiments on a square film made up of polydimethylsiloxane (PDMS) rubber embedded with a hole and notches and then compare our simulation predictions with the experimental observations. Furthermore, we visualize the evolution of stretch and damage values in chains oriented along different directions to assess the predictive capacity of the model. In conclusion, the results are also compared with another existing model to evaluate the utility of our model in accurately simulating the fracture behavior of rubber-like materials.

42 ENGINEERING↗
DOE OSTI2024

Review of Technical Photovoltaic Key Performance Indicators and the Importance of Data Quality Routines

Technical key performance indicators (KPIs) are important metrics used to assess and quantitatively summarize various aspects of photovoltaic (PV) systems, including long-term performance, economic viability, and carbon footprint. Herein, a group of experts of the International Energy Agency's Photovoltaic Power Systems Programme Task 13 collect and describ the most important technical KPIs used in the industry. Thereby, a set of best practices for reliably handling PV system data is presented and the impact of data quality and climatic variability on KPI calculation is investigated. Further, the effective use of technical KPIs allows triggering data-driven and informed decisions to optimize PV systems and providing a comprehensive overview of how PV systems operate across different conditions and climates. With the worldwide growth of the PV industry, more companies operate/own PV systems in different regions, where the climatic and seasonal profiles differ. This requires context-aware evaluation of KPIs, or the judicious application of multiple KPIs, to ensure that each asset is evaluated correctly. Beyond that, there is untapped potential in the utilization of KPIs through geospatial mapping and extrapolation of fleet KPIs. This study demonstrates that the uncertainty in KPI estimation is not well understood and depends on data quality, climatic variability, and system configuration.

14 SOLAR ENERGY↗
DOE OSTI2024

OC6 Phase IV: Validation of CFD Models for Stiesdal TetraSpar Floating Offshore Wind Platform

ABSTRACT With only a few floating offshore wind turbine (FOWT) farms deployed anywhere in the world, FOWT technology is still in its infancy, building on a modicum of real‐world experience to advance the nascent industry. To support further development, engineers rely heavily on modeling tools to accurately portray the behavior of these complex systems under realistic environmental conditions. This reliance creates a need for verification and validation of such tools to improve reliability of load and dynamic response prediction and analysis capabilities of FOWT systems. The Offshore Code Comparison Collaboration, Continued with Correlation and unCertainty (OC6) project was created under the framework of the International Energy Agency to address this need and considers a three‐sided verification and validation between engineering level models, computational fluid dynamics (CFD), and experimental results. In this paper, a novel floating offshore wind platform, the Stiesdal TetraSpar, is simulated using CFD under the load conditions defined by Phase IV of the OC6 project. The comparison of these CFD results against the experimental results demonstrated the ability to predict the platform response to waves when imposing the measured wave signals as input. Although validation versus experiment was largely successful, the damping behavior was impacted by uncertainties likely originating from the mooring system and sensor umbilical cable. This extensive comparison effort with multiple CFD practitioners offers insight into best practices to achieve reliable results.

17 WIND ENERGY↗
DOE OSTI2021

Role of scan strategies and heat treatment on grain structure evolution in Fe-Si soft magnetic alloys made by laser-powder bed fusion

A major goal in printing soft magnetic Fe-Si steels using additive manufacturing is to take advantage of the potential for complex geometric designs and site-specific grain control. One major step in the processing of these alloys is understanding how processing parameters might impact how the as-built microstructure responds to annealing (i.e. the annealing response). The impact of scan strategy on the annealing response for thin wall geometries is specifically explored. Two scan strategies were explored for a thin wall geometry that produced a strongly columnar grain structure and equiaxed grain structure. Additionally, samples from both scan strategies annealed at 1200 °C showed a marked difference in annealing response with the more equiaxed sample seeing full recrystallization and grain growth, while the more columnar grain structure saw little change in microstructure. After analysis through characterization techniques and thermal-mechanical simulations Differences in internal energy within the grains were ruled out because calculated GND density values were similar for both samples. The formation of secondary particles was ruled out as a contributing factor due to the type of oxide formations and their size. It was concluded that the contributing factor to the difference in the annealing response were a difference in the resulting grain size and the density of high angle grain boundaries. These two differences were largely attributed to differences in the thermal gradient conditions due to grains preferentially growing in the direction of the steepest thermal gradient.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Advances in benchmarking and round robin testing for PEM water electrolysis: Reference protocol and hardware

While the number of publications in the PEM water electrolysis community increases each year, no common ground concerning reference hardware (test cells and test bench) and testing protocols has been yet established. This would, however, be necessary for the comparability of experimental results. First attempts for such reference hardware and procedures have been made in the framework of the Task 30 Electrolysis within the Technology Collaboration Programme on Advanced Fuel Cells (AFC TCP) of the International Energy Agency (IEA). Since then, improvements of both the test hardware (test cell and components) as well as the measurement protocol were identified, and a revised methodology and key results based on a comprehensive measurement series have been obtained. A detailed protocol for testing commercial reference components with a reference laboratory test cell developed in-house by Fraunhofer ISE is presented. For evaluation of the protocol and the hardware, it was tested at three different institutions at the same time. Impedance spectroscopic and polarization data was acquired and analyzed. The obtained differences in performance were calculated to give the community an expectation window to compare own data to. Finally, the importance of a thorough temperature control and the conditioning phase are demonstrated.

08 HYDROGEN↗
DOE OSTI2023

A multiscale phase field fracture approach based on the non-affine microsphere model for rubber-like materials

Rubber-like materials have a broad scope of applications due to their unique properties like high stretchability and increased toughness. Hence, computational models for simulating their fracture behavior are paramount for designing them against failures. In this study, the phase field fracture approach is integrated with a multiscale polymer model for predicting the fracture behavior in elastomers. At the microscale, damaged polymer chains are modeled to be made up of a number of elastic chain segments pinned together. Using the phase field approach, the damage in the chains is represented using a continuous variable. Both the bond stretch internal energy and the entropic free energy of the chain are assumed to drive the damage, and the advantages of this assumption are expounded. A framework for utilizing the non-affine microsphere model for damaged systems is proposed here by considering the minimization of a hypothetical undamaged free energy, ultimately connecting the chain stretch to the macroscale deformation gradient. At the macroscale, a thermodynamically consistent formulation is derived in which the total dissipation is assumed to be mainly due to the rupture of molecular bonds. Using a monolithic scheme, the proposed model is numerically implemented and the resulting three-dimensional simulation predictions are compared with existing experimental data. The capability of the model to qualitatively predict the propagation of complex crack paths and quantitatively estimate the overall fracture behavior is verified. Additionally, the effect of the length scale parameter on the predicted fracture behavior is studied for an inhomogeneous system.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Dependence of predicted bulk properties of hexagonal hydroxyapatite on exchange-correlation functional

Hydroxyapatite (HA) is the main component of human bones and teeth. HA has also been widely applied in various technological fields due to its unique properties. Reliable computational simulations for this critical material often require input parameters obtained from the first-principles density functional theory (DFT) calculations with appropriate exchange–correlation functionals. Previous DFT calculations are insufficient to verify the reliabilities of various functionals, particularly as they fail to assess predictions for multiple properties of the material. In this paper, we first select 18 different functionals to calculate geometric, elastic, electronic, and thermodynamic properties of hexagonal HA bulk crystal. We find that the results from optB86b-vdW and optB88-vdW functionals with dispersion corrections have the overall best agreement with available experimental data. Then, we choose optB88-vdW functional, as well as PBE functional without dispersion corrections as a comparison, to perform extensive first-principles DFT phonon calculations under the quasiharmonic approximation. Various thermodynamic properties (including phonon contributions to internal energy, entropy, and Helmholtz free energy) and thermal parameters (including thermal expansions of volume, thermal expansion coefficients, heat capacities, isothermal bulk moduli, etc.) versus temperature are consequently obtained. By comparing these quantities, we report that the results from optB88-vdW functional have significantly better agreement with available experimental data than those from PBE functional although the latter has been widely used in previous DFT calculations for HA-based material systems.

36 MATERIALS SCIENCE↗