Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Solid modeling”

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 685 records · Page 38

Direct imaging of shock wave splitting in diamond at Mbar pressure

Understanding the behavior of matter at extreme pressures of the order of a megabar (Mbar) is essential to gain insight into various physical phenomena at macroscales—the formation of planets, young stars, and the cores of super-Earths, and at microscales—damage to ceramic materials and high-pressure plastic transformation and phase transitions in solids. Under dynamic compression of solids up to Mbar pressures, even a solid with high strength exhibits plastic properties, causing the induced shock wave to split in two: an elastic precursor and a plastic shock wave. This phenomenon is described by theoretical models based on indirect measurements of material response. The advent of x-ray free-electron lasers (XFELs) has made it possible to use their ultrashort pulses for direct observations of the propagation of shock waves in solid materials by the method of phase-contrast radiography. However, there is still a lack of comprehensive data for verification of theoretical models of different solids. Here, we present the results of an experiment in which the evolution of the coupled elastic–plastic wave structure in diamond was directly observed and studied with submicrometer spatial resolution, using the unique capabilities of the x-ray free-electron laser (XFEL). The direct measurements allowed, for the first time, the fitting and validation of the 2D failure model for diamond in the range of several Mbar. Our experimental approach opens new possibilities for the direct verification and construction of equations of state of matter in the ultra-high-stress range, which are relevant to solving a variety of problems in high-energy-density physics.

36 MATERIALS SCIENCE↗

Continuous recovery of phosphoric acid and Rare-Earths containing particles from phosphoric acid sludge using a decanter centrifuge

Recovery of rare earth elements (REEs) from various industrial and natural streams currently draws significant attention in efforts to meet the demands of the manufacturing industry. Among many industrial byproducts and waste streams, phosphoric acid sludge could be one of the most economically feasible resources for the recovery of REEs because solid particles in the sludge contain relatively concentrated REEs, up to 3,000 ppm, while the liquid component of the sludge is valuable phosphoric acid (P 2 O 5 ) that can be recovered and returned to the main product. Due to high viscosity and large solids content (e.g., 30–40 %), however, this byproduct stream requires multistep separation and purification processes. In this study, a single-step process involving a continuous-flow decanter centrifuge (CFDC) was employed to investigate its feasibility for continuous solid/liquid separation from real phosphoric acid sludge. High centrifugal forces generated from up to 1500 G gravity acceleration separate solid particles from the sludge, generating a liquid-rich stream and a solids-rich stream at the exit of the CFDC. A single pass of phosphoric-acid sludge through the CFDC yielded 95 % liquid recovery and 90 % recovery of REEs-containing solids from 20 to 34 wt% solids-containing sludge. Finally, a reduced order model developed for the CFDC operation showed good agreement with experimental data, and preliminary technoeconomic analysis revealed potential process feasibility.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Simulating hindered grain boundary diffusion using the smoothed boundary method

Abstract Grain boundaries can greatly affect the transport properties of polycrystalline materials, particularly when the grain size approaches the nanoscale. While grain boundaries often enhance diffusion by providing a fast pathway for chemical transport, some material systems, such as those of solid oxide fuel cells and battery cathode particles, exhibit the opposite behavior, where grain boundaries act to hinder diffusion. To facilitate the study of systems with hindered grain boundary diffusion, we propose a model that utilizes the smoothed boundary method to simulate the dynamic concentration evolution in polycrystalline systems. The model employs domain parameters with diffuse interfaces to describe the grains, thereby enabling solutions with explicit consideration of their complex geometries. The intrinsic error arising from the diffuse interface approach employed in our proposed model is explored by comparing the results against a sharp interface model for a variety of parameter sets. Finally, two case studies are considered to demonstrate potential applications of the model. First, a nanocrystalline yttria-stabilized zirconia solid oxide fuel cell system is investigated, and the effective diffusivities are extracted from the simulation results and are compared to the values obtained through mean-field approximations. Second, the concentration evolution during lithiation of a polycrystalline battery cathode particle is simulated to demonstrate the method’s capability.

Materials Science↗

Thermodynamic evaluation and modeling of the CeF 3 -LaF 3 -LiF-NaF system

The thermodynamic properties of the LiF-NaF-CeF 3 -LaF 3 pseudo-quaternary system have been evaluated, and models for the relevant phases have been revised or newly developed. The solid solubility of NaF in LaF 3 and CeF 3 was described in the model based on reported phase-equilibria observations and results from density functional theory calculations. The Gibbs energy functions for the NaCeF 4 and NaLaF 4 phases have been updated with recently reported thermochemical measurements. The Gibbs energy functions and model relations for the LiF-NaF-LaF 3 and LiF-NaF-CeF 3 pseudo-ternary systems were successfully optimized based on thermochemical properties and phase equilibria observations.

CALPHAD↗

NRAP-Open-IAM: FutureGen2 Component Models

This report describes the development and testing of three component models for NRAP-Open-IAM, the National Risk Assessment Partnership’s open-source integrated assessment model. The FutureGen2 Lookup Table Reservoir component model is based on interpolation of data from a set of lookup tables. The lookup tables contain pressures and saturations predicted by multiphase flow simulations performed with a reservoir simulator. The FutureGen2 Above Zone Monitoring Interval (AZMI) component is a surrogate model that can be used to estimate the impact that carbon dioxide (CO 2 ) and brine leaks from the CO 2 storage reservoir at the FutureGen 2.0 site might have had on overlying aquifers or monitoring units were a leak to occur. The model estimates the size of “impact plumes” according to five metrics: pH, Total Dissolved Solids (TDS), pressure, dissolved CO 2 and temperature. The FutureGen2 Aquifer component is similar, but is limited to four metrics: pH, Total Dissolved Solids (TDS), pressure, and dissolved CO 2 . The input parameters for each model are the same, but the Aquifer component is applicable to depths between 100 m and 700 m and the AZMI component is applicable from depths between 700 m and 1050 m.

42 ENGINEERING↗

Modelling a compressible packed bed flow-through washing and deacetylation reactor for corn stover pretreatment

The deacetylation process is a key step in the lignocellulosic biomass pretreatment process that currently is done in batch reactors. A flow-through diffuser type reactor is a promising approach to make the process more efficient and economically favorable. A washing step before deacetylation can further reduce chemical usage and increase the process efficiency. However, the highly compacted biomass particles can cause significant pressure-drop for fluid flow through the bed, which in the worst case can cause self-enhancing pressure-drop leading to reactor blocking. Detailed understanding of the fundamentals of flow through a compressible packed bed is necessary to scale up the washing and deacetylation process. This work closely examines the classical Kozeny Carman equation and proposes modifications to accurately predict the pressure-drop at different compaction states. Here, a compressible packed bed model is further developed based on the interaction between solid mechanics and fluid dynamics. Experimental measurements of water flow through a packed bed of corn stover show good agreement with the developed model. Corn stover washing process scale up simulations show that smaller reactor diameter, higher friction lining material, lower process throughput, and higher process temperature is favorable in stabilizing the bed from extensive compression when the bed height is fixed.

09 BIOMASS FUELS↗

Caustic Neutralization and Precipitation of Acidic Dissolved Simulated Stainless Steel–Clad Plutonium and Plutonium/Uranium Nuclear Fuel

Simulated dissolved stainless steel (SS) clad Pu and Pu/U nuclear fuel in HNO 3 was neutralized to a free hydroxide (OH - ) concentration of 0.6 M. A thermal neutron poison, Gd, was added to the simulants at concentrations of either ~3 - 6 g/L or ~37 - 38 g/L. The supernate Pu concentration the day of neutralization ranged from 0.48 to 8.75 mg/L. The supernate Pu concentration of a simplified simulant neutralized to 0.6 M OH - above precipitated solids containing Pu was demonstrated to decrease over 18 days. A significant portion of precipitated Pu was found to be insoluble in 8 M HNO 3 at ambient temperature, but essentially quantitative Pu dissolution was achieved in 11.5 M HNO 3 /0.1 M KF at 100 °C. The difficulty in dissolving the Pu precipitate is believed to be due to the formation of refractory PuO 2 •xH 2 O during the neutralization process. Initial Gd concentrations of ~37 – 38 g/L were found to result in a greater Al precipitation when neutralized to 0.6 M OH - than initial Gd concentrations of ~3 – 6 g/L. Physical properties of the resultant slurries were measured and used to calculate limiting flowrates and slurry velocities by gravity only in transfer piping between the Savannah River Site’s H-Canyon Facility and the Concentration, Storage, and Transfer Facility (CSTF). These results were compared to calculated deposition velocities to predict if solids would settle during the transfer. The Newtonian model was found to be reasonable for each diluted slurry evaluated. Deposition velocities of Pu containing slurries are lower than nuclear fuel slurries primarily composed of U due to the high density of Pu solids. In conclusion, dilution of slurries reduces the margin between the slurry and deposition velocities due to the reduction in viscosity because higher viscous forces on the particles promote maintained suspension.

Actinide Neutralization↗

Retaining Systems Engineering Model Meaning Through Transformation: Demo 2

Digital engineering strategies typically assume that digital engineering models interoperate seamlessly across the multiple different engineering modeling software applications involved, such as model- based systems engineering (MBSE), mechanical computer-aided design (MCAD), electrical computer-aided design (ECAD), and other engineering modeling applications. The presumption is that the data schema in these modeling software applications are structured in the familiar flat- tabular schema like any other software application. Engineering domain-specific applications (e.g., systems, mechanical, electrical, simulation) are typically designed to solve domain-specific problems, necessarily excluding explicit representations of non-domain information to help the engineer focus on the domain problems (system definition, design, simulation). Such exclusions become problematic in inter-domain information exchange. The obvious assumptions of one domain might not be so obvious to experts in another domain. Ambiguity in domain-specific language can erode the ability to enable different domain modeling applications to interoperate, unless the underlying language is understood and used as the basis for translation from one application to another. The engineering modeling software application industry has struggled for decades to enable these applications to interoperate. Industry standards have been developed, but they have not unified the industry. Why is this? The authors assert that the industry has relied on traditional database integration methods. The basic issue prohibiting successful application integration then is that traditional database-driven integration does not consider the distinct languages of each domain. An engineering models meaning is expressed through the underlying language of that engineering domain. In essence, traditional integration methods do not retain the semantic context (meaning) of the model. The basis of this research stems from the widely held assumption that systems engineering models are (or can be) structured according to the underlying semantic ontology of the model. This assumption can be imagined from two thoughts. 1) Digital systems engineering models are often represented using graph theory (the graph of a complex systems model can contain millions of nodes and edges). When examining the nodes one at a time and following the outbound edges of each node one by one, one can end up with rudimentary statements about the model (i.e., node A relates to node B), as in a semantic graph. 2) Likewise, from the study of natural languages, a sentence can be structured into unambiguous triples of subject-predicate-object within formal and highly expressive semantic ontologies. The rudimentary statements about a systems model discerned with graph theory closely mimic the triples used in the ontologies that try to structure natural languages. In other words, a systems models semantic graph can be (or is) structured into an ontology. Additionally, it is well established in industry that through natural language processing (NLP), which provides the means to create language structures, that computers can interpret ontological graphs. Therefore, the authors hypothesized that if the integrity of the underlying semantic structure of a systems model is retained, the contextual meaning of the model is retained. By structuring system models into the triples of the underlying ontology during the transformation from one MBSE application to another, the authors have provided a proof of the concept that the meaning of a system model can be retained during transformation. The authors assert that this is the missing ingredient in effective systems model-to-model interoperability. ACKNOWLEDGEMENTS The authors would like to thank the FY19 Model Interoperability team members who provided a solid foundation for the FY20 team to leverage: John McCloud, for the work he did to guide us toward the right use of technology that will appropriately discover and manipulate ontologies. Carlos Tafoya, for the work he did to develop an application programming interface (API)/Adapter that would export ontology-based data from GENESYS. Peter Chandler, for the work he did to architect our overall integration solution, with an eye toward the future that would influence a large-scale federated production-level systems engineering digital model ecosystem.

42 ENGINEERING↗

Continuum-Level Modeling of Li-Ion Battery SEI by Upscaling Atomistically Informed Reaction Mechanisms

Understanding and controlling solid-electrolyte interphase (SEI) formation to stabilize cell performance is a significant challenge for next-generation Li-ion battery technologies. In recent years, computational modeling has become an essential tool in providing fundamental insights into SEI properties and dynamics. However, neither atomistic nor continuum-level approaches alone can capture the complexities of SEI chemistry across all relevant length and time scales. In this work, a continuum-level model is developed that is informed by reaction mechanisms obtained from first-principle calculations. The atomistically informed continuum-level model is used to understand electrolyte degradation, including the decomposition of ethylene carbonate (EC), ethyl methyl carbonate (EMC), and fluoroethylene carbonate (FEC). The model presented here is the most chemically complex continuum-level SEI model in the literature to date. The SEI model is calibrated against experimental irreversible leakage currents and shows qualitative agreement with expected SEI growth trends. The model framework is expected to accelerate fundamental understanding of SEI formation, facilitate mechanism development feedback, and dynamically interact with experimental insights.

continuum-level model↗

Fuel Performance Modeling Internship Final Presentation

This study examines the performance of U-Zr and U-Pd-Zr annular metallic fuel rodlets and details the current status of modeling efforts regarding U-Pu-Zr solid metallic fuel rodlets which were experimentally evaluated as part of the Advanced Fuels Campaign (AFC). The AFC mission is to develop novel fuel technologies and facilitate the implementation of those technologies by industry partners. A key objective is to improve steady-state and transient performance over current fuel types. The experiments of interest in this study included metallic fuel rodlets within HT-9 cladding which were placed in SS-316 capsules and inserted in the Advanced Test Reactor (ATR). Certain mechanical and thermal conditions cannot be directly evaluated through experiments and fuel performance modeling is used to shed light on this evolution over time. In this study, BISON Multiphysics simulations are leveraged to investigate the state of the fuel system throughout and after the experimental conditions.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Fluoride-Cooled High-Temperature Pebble-Bed Reactor Reference Plant Model

In this report we present work performed in Fiscal Year 2022 that demonstrates the modeling and simulation of a fully coupled neutronics thermal hydraulics reference plant model for a fluoride-cooled high-temperature pebble-bed reactor. The multiphysics model is developed on the Nuclear Regulatory Commission’s Comprehensive Reactor Analysis Bundle (BlueCRAB) available on the Idaho National Laboratory’s high-performance computer, which natively and seamlessly couples Griffin, Pronghorn, and the BISON Multiphysics Object-Oriented Simulation Environment based applications. Griffin provides reactor physics capabilities, including depletion to the equilibrium core, k-eigenvalue, adjoint, and transient. The unique direct equilibrium core capability in Griffin is based on a streamline methodology to spatially deplete the pebbles into burnup groups. Pronghorn solves the porous medium equations for the fluid regions and conduction in the solid regions and incorporates a fluidic diode model to simulate the transition from forced to natural convection during accident scenarios. MOOSE modules solves thermal conduction problems for the pebbles and tristructural isotropic in the pebble-bed core, thus providing the fuel and moderator spatial fields for each pebble burnup group. The neutronics feedback relies primarily on fuel, moderator, and reflector temperatures as as well as the FLiBe salt density. Here, we present results for the uncoupled equilibrium core and perform comparisons to equivalent Monte Carlo models. The power distributions and kinetic parameters obtained with Griffin are consistent with those computed with Griffin. We demonstrate a noticeable improvement with the use of discrete ordinates method (SN) transport. The coupled steady-state equilibrium core provides the initial condition for two time-dependent problems: a control rod withdrawal event and an unprotected loss of flow event. In both cases, the reactor design is self-stabilizing and the solutions are consistent with the expected physics. Although this model is prototypical regarding BlueCRAB’s capabilities, its results are consistent with published work by Kairos Power and other research entities. Significant improvements to the model are planned in future work.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

An ALEGRA MHD Spark Model

While modeling a generic pulse transformer, we became interested in the possibility of electric sparks between winding layers in a solid encapsulant. We significantly modified a previously developed ALEGRA MHD model of a generic spark in lexan. The cumulative modifications are significant enough to report here. Possibly the most significant modification was a change in how the simulated spark is initiated from a thin initial channel. The change was from imposing an initial hot temperature to imposing a conductivity floor. The reasons and comparisons of results are included. The second significant change was to replace a fixed current rise rate with an external circuit model. We built a model specifically mimicking the distributed inductance and stray capacitance between the coil turns closest to the modeled spark. Excursions from nominal values examine the sensitivity of resulting behaviors to extreme capacitance and inductance values.

42 ENGINEERING↗

Testing- and Model- Based Optimization of Coal-fired Primary Heater Design for Indirect Supercritical CO 2 Power Cycles (Final Scientific and Technical Report)

The overall objective of this project was to perform the R&D necessary to mitigate the risk associated with the design of a primary heat exchanger for a solid-fired combustion system coupled with an indirect-fired closed-loop Brayton Cycle utilizing supercritical CO 2 . The key technological hurdle was the coupling of a solid-fuel firing system with the primary heater, which poses a singular challenge, which is the management of burner performance and operational conditions in a way to manage heat exchanger tube metal temperatures and temperature ramp rates in the absence of fluid phase change on the inside of the tubes. We designed and built the first ever pseudo power system employing a simple recuperated supercritical CO 2 closed-loop Brayton Cycle coupled to a solid-fuel fired system. Advanced coupled CFD and process modeling were used to design the primary heat exchanger (PHX), which consisted of both radiative and convective sections, to limit tube metal temperatures resulting from the heat release profile of the solid fuel flame near the radiative tubes. The heat exchanger was designed to produce finished CO 2 temperatures of 600 °C a pressure of 20.7 MPa and CO 2 flow of 5.5 kg/s. The constructed PHX was capable of 1.2 MWth heat uptake. During design of the PHX, the modeling showed that most variables influencing flame shape (burner stoichiometric ratio and register velocities and swirl) were not suitable to manage heat flux to the metal surfaces. This is because they substantially increased adiabatic flame temperature through the influence of localized stoichiometric ratio. Excess air and firing rate were the two most powerful variables that could be used to control tube surface temperatures. The coupled system was operated for a total of 407 hours, with the longest continuous run of 248 hours. For 62% of the operational time, the unit was unmanned and in automatic control. The fuels used for the testing included natural gas, two Utah Bituminous coals, woody biomass, and bagasse. During the testing we were able to verify the 1.2 MWth heat uptake and we operated at a finished CO 2 temperature of 607 °C and a pressure of 20.3 MPa simultaneously. The real-time corrosion rate of the Super 304H tube CO 2 surface in the region of the radiative section of the PHX were measured, at an approximate temperature of 550 °C. The two key variables related to corrosion rate are the pressure and flow rate of the CO 2 . A technoeconomic analysis was performed at a scale of 120 MWE. The updated analysis showed that the efficiency of an sCO 2 power producing plant will be related to the pressure drop of the PHX.

01 COAL, LIGNITE, AND PEAT↗

DECOVALEX-2023: Task E Final Report

This is the Task E final report for DECOVALEX-2023. Task E is focused on understanding thermal, two-phase hydrological, and mechanical (TH2M) processes, especially related to predicting brine migration in the excavation damaged zone around a heated excavation in salt. Salt is attractive as a disposal medium for radioactive waste because it is self-healing and is essentially impermeable and essentially non-porous in the far field (away from excavations). Investigation of the short-term (days to years) near-field (centimeters to tens of meters) behavior of salt is important for radioactive waste disposal because this early period strongly controls the amount of brine in a salt repository. Brine leads to corrosion of waste forms and waste packages, and possible dissolution of radionuclides with brine transport being a potential transport vector to the accessible environment. The main test case used in Task E is the ongoing Brine Availability Test in Salt (BATS) heater test located underground at the Waste Isolation Pilot Plant (WIPP) near Carlsbad, New Mexico, USA. The Task was divided into a series of Steps. Step 0 was an introduction to processes in salt, that included matching historical unheated brine inflow data from boreholes at WIPP and matching temperature observations during BATS heater test 1a. Step 1 included validation of models against a thermo-poroelastic analytical solution relevant to heated boreholes in salt, and two-phase flow around an excavation in salt. Step 2 required all the individual components covered in steps 0 and 1 to come together to match observed brine inflow behavior during the BATS 1a heater test. There were a range of approaches from the teams, from mechanistic to prescriptive. Given the uncertainties in the problem, some teams used one- or two-dimensional models of the processes, while other teams included more geometrical complexity in three-dimensional models. The key learning points from Task E have been: • Heat conduction through salt typically requires non-linear thermal conductivity (as a function of temperature), but most models do a good job matching observations, given appropriate adjustments to the applied power and some thermocouple locations. • Thermal pressurization requires coupled thermal-hydrological-mechanical (THM) responses that consider the thermal expansion of the fluid and solid phases. • Initialization of two-phase flow models around a borehole or excavation in salt are more realistically represented as “wetting up”, rather than “drying down” (i.e., the initial state after excavation is mostly dry, rather than mostly wet). • The BATS 1a heater test includes a significant release of brine after the end of heating, which requires a large increase in permeability to recreate. Task E has been a great learning experience for all the teams involved, and feedback from the modeling teams has led to changes in the design of follow-on BATS experiments, which are now ongoing underground at WIPP. There was a balance throughout the task between freedom to model phenomena how each team saw fit, and prescriptiveness in problem design to bring the modeling teams closer together to allow attribution of smaller differences between models to different modeling choices. The modeling approaches seem to go through two phases: an early phase of discovery or testing, and a later phase of refinement and improvement. In future modeling efforts, different field data could be used (e.g., BATS 2) and more time should be included in the processes for teams to make multiple model refinement or even significant changes to their conceptual model or setup, based on lessons learned from the modeling exercise.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Frequency-domain hot-wire sensor and 3D model for thermal conductivity measurements of reactive and corrosive materials at high temperatures

High temperature solids and liquids are becoming increasingly important in next-generation energy and manufacturing systems that seek higher efficiencies and lower emissions. Accurate measurements of thermal conductivity at high temperatures are required for the modeling and design of these systems, but commonly employed time-domain measurements can have errors from convection, corrosion, and ambient temperature fluctuations. Here, we describe the development of a frequency-domain hot-wire technique capable of accurately measuring the thermal conductivity of solid and molten compounds from room temperature up to 800 °C. Therefore, by operating in the frequency-domain, we can lock into the harmonic thermal response of the material and reject the influence of ambient temperature fluctuations, and we can keep the probed volume below 1 µl to minimize convection. The design of the microfabricated hot-wire sensor, electrical systems, and insulating wire coating to protect against corrosion is covered in detail. Furthermore, we discuss the development of a full three-dimensional multilayer thermal model that accounts for both radial conduction into the sample and axial conduction along the wire and the effect of wire coatings. The 3D, multilayer model facilitates the measurement of small sample volumes important for material development. A sensitivity analysis and an error propagation calculation of the frequency-domain thermal model are performed to demonstrate what factors are most important for thermal conductivity measurements. Finally, we show thermal conductivity measurements including model data fitting on gas (argon), solid (sulfur), and molten substances over a range of temperatures.

47 OTHER INSTRUMENTATION↗

Investigation of a Hybrid Carbon Capture System at a NGCC Power Plant That Performs Direct Air Capture During Off-Peak Hours

The goal of this project was to investigate the feasibility of coupling NGCC power plants with membrane and solid sorbent carbon capture to achieve net negative carbon emissions. This novel hybrid system has a couple of potential benefits: 1) the carbon capture systems could help the NGCC plant by minimizing its carbon footprint and giving it a new revenue stream from direct air capture, and 2) the NGCC plant could help direct air capture scale up quickly and costeffectively by providing it with electricity, steam and compressed air streams. In order to investigate this hybrid design, we developed a complex, integrated system model consisting of: 1) a NGCC power plant model developed by Dr. Bhattacharyya’s team at WVU, 2) a membrane carbon capture model developed by Dr. Lipscomb’s team at U Toledo, and 3) a solid sorbent carbon capture model developed by Dr. Hornbostel’s team at U Pittsburgh. Mr. Matuszewski’s team at AristoSys LLC integrated these three models into a single, integrated model using FOQUS, an optimization tool developed by the NETL CCSI 2 program. Optimizations were then performed on this integrated model to determine the cost and carbon capture rates under a range of electricity and carbon price scenarios.

20 FOSSIL-FUELED POWER PLANTS↗

Three-dimensional traction microscopy accounting for cell-induced matrix degradation

Tractions exerted by cells on the extracellular matrix (ECM) are critical in many important physiological and pathological processes such as embryonic morphogenesis, wound healing, and cancer metastasis. Three-dimensional Traction Microscopy (3DTM) is a tool to quantify cellular tractions by first measuring the displacement field in the ECM in response to these tractions, and then using this measurement to infer tractions. Most applications of 3DTM have assumed that the ECM has spatially-uniform mechanical properties, but cells secrete enzymes that can locally degrade the ECM. In this work, a novel computational method is developed to quantify both cellular tractions and ECM degradation. In particular, the ECM is modeled as a hyperelastic, Neo-Hookean solid, whose material parameters are corrupted by a single degradation parameter. The feasibility of determining both the traction and the degradation parameter is first demonstrated by showing the existence and uniqueness of the solution. An inverse problem is then formulated to determine the nodal values of the traction vector and the degradation parameter, with the objective of minimizing the difference between a predicted and measured displacement field, under the constraint that the predicted displacement field satisfies the equation of equilibrium. The inverse problem is solved by means of a gradient-based optimization approach, and the gradient is computed efficiently using appropriately derived adjoint fields. The computational method is validated in-silico using a geometrically realistic neuronal cell model and synthetic traction and degradation fields. It is found that the method accurately recovers both the traction and degradation fields. Moreover, it is found that neglecting ECM degradation can yield significant errors in traction measurements. Finally, our method can extend the range of context where tractions can be appropriately measured.

42 ENGINEERING↗

Multiphase Processing of the Water-Soluble and Insoluble Phases of Biomass Burning Organic Aerosol

Biomass burning is one of the most significant sources of organic aerosol in the atmosphere. Biomass burning organic aerosol (BBOA) has been observed to undergo liquid– liquid phase separation (LLPS) to give core–shell morphology with the hydrophobic phase encapsulating the hydrophilic phase, potentially impacting the evolution of light-absorbing components, i.e., brown carbon (BrC), through multiphase processes. Here, we demonstrate how multiphase processing differs between the watersoluble (i.e., hydrophilic) and insoluble (i.e., hydrophobic) phases of BBOA in terms of reactive uptake of ozone in a coated-wall flow tube. Effects of relative humidity (RH) and ultraviolet (UV) irradiation were investigated. Experimental timeseries were used to inform simulations using multilayer kinetic modeling. Among non-irradiated thin films, the uptake coefficient was greatest for the water-soluble phase at 75% RH (3 × 10 –5 , corresponding to a diffusion coefficient of BrC, D BrC , of 3 × 10 –9 cm 2 s –1 ) and least for the same phase at 0% RH (1 × 10 –5 , corresponding to D BrC of 1 × 10 –10 cm 2 s –1 ). The uptake coefficient for the water-insoluble phase fell between these two (about 1.5 × 10 –5 ), regardless of RH, and the corresponding D BrC increased only slightly (8 × 10 –10 cm 2 s –1 at 0% RH to 9 × 10 –10 cm 2 s –1 at 75% RH). The uptake coefficients of both phases at 0% RH decreased significantly after UV irradiation, consistent with a transition from viscous liquid to solid and supported by qualitative microscopy observations. Modeling multiphase ozone oxidation of primary BrC components in the atmosphere demonstrated, first, that LLPS may extend the lifetime of water-soluble BBOA encapsulated by water-insoluble species by a factor of 1.5 at moderate to high RH and, also, that UV irradiation may extend the lifetime of both phases by more than a factor of 2.5.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗