Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “statistical mechanics models”

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

Time-dependent THMC properties and microstructural evolution of damaged rocks in excavation damage zone

Modeling coupled thermo-hydro-mechanical-chemical (THMC) processes in host rocks near high-level nuclear waste (HLW) repositories at various time scales is an extremely challenging task. The current study integrates experimental, theoretical, and numerical methods in assessing the evolution of excavation damage zone (EDZ) over time and its implication on the long-term migration of hazardous species. Argillite and rock salt and are the focus of this study. The first part of the report presents a novel time-dependent directional microcrack damage theory for generic brittle rocks. It features detailed statistical description of the microcracks within a damaged solid, permitting a direct upscaling of microscale processes such as crack growth kinetics, crack closure/opening, sliding friction to explain the macroscopic creep, nonlinear elasticity, shear dilation, and loading-unloading hysteresis. This provides a basic platform for describing the anisotropic mechanical and transport properties of damaged rocks during excavation and subsequent THMC loadings. The model is validated and numerically implemented to Finite Element (FE) package ABAQUS through the user-defined material (UMAT) interface and have demonstrated great potential in resolving the time-dependent and anisotropic evolution of damage in the EDZ. The second part of the report presents a multi-scale experimental effort in characterizing the thermal, hydraulic, and mechanical properties of Mancos shale and Avery Island salt. For the Mancos shale, triaxial compression tests are performed at different confining pressures and temperatures to probe its thermomechanical properties relevant to HLW repositories. The obtained stress-strain data are interpreted using the proposed directional damage theory. Post-test specimens are subjected to gas permeability tests to reveal the correlation between permeability and the degree of microcracking. At microscale, temperature-controlled nanoindentation tests are performed and found a linear correlation between fracture toughness and elastic modulus from 25°C to 300°C. For Avery Island salt, we have designed and manufactured a novel relative-humidity controlled uniaxial creep device. Long-term creep tests at low stresses (< 5 MPa) are performed at different levels of relative humidity (RH). Besides confirming the much higher creep rates as one would expect through extrapolating the high-stress creep data, the results reveal that the steady-state creep rate of rock salt is strongly dependent on the ambient RH, an aspect that is often neglected in the literature. Both behaviors can be attributed to the pressure-solution creep mechanism which dominates at low stress and high RH levels. The third part of the report explores a set of numerical strategies in modeling the THMC behavior of porous geomaterials. A fully implicit, monolithic FE solution that can flexibly interface with different material models and coupling mechanisms for THMC problems is developed and verified through the ABAQUS user-defined element (UEL) interface. The scheme is used to study the THM response of a hypothetical HLW storage site with reference to an existing in-situ heater test. Strategies for integrating the UEL and the microcrack UMAT are suggested. Another numerical endeavor of this study is to implement a higher-order asymptotic homogenization method to account for the heterogeneous porous structures. The same method is then extended to perform microstructure-informed thermo-mechanical modeling of generalized continua. The above outcomes of this project provide a strong thrust towards enhancing the fundamental understanding and modeling capability of the long-term evolution of host rocks in EDZ, thus helping achieve the design goal of 1-million-year isolation of high-level nuclear wastes.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

A coupled vacancy diffusion-dislocation dynamics model for the climb-glide motion of jogged screw dislocations

Here, we develop a novel model to study the climb/glide motion of jogged screw dislocations within the discrete dislocation dynamics (DDD) framework. We present results for the dependence of the climb velocity on the applied stress and on the jog size and distribution statistics. We show that the model predictions are consistent with experimental data in both γ-TiAl and Zircaloy-4. The ranges of the applied stress and jog spacing that determine the dominance of one of three dislocation mechanisms are identified. These are the jog dragging, dipole dragging, and dipole bypass mechanisms, respectively. The overall dislocation motion in the jog dragging regime is composed of glide of screw segments and climb of jogs, controlling the plastic strain and the creep rate, respectively. Based on current simulations and on a detailed examination of the predicted jog heights compared to experiments, we advance the hypothesis that a combination of jog dragging and dipole bypass mechanisms is necessary to reproduce the high creep rate observed in some experiments.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Towards Quantum Computing Phase Diagrams of Gauge Theories with Thermal Pure Quantum States

The phase diagram of strong interactions in nature at finite temperature and chemical potential remains largely theoretically unexplored due to inadequacy of Monte-Carlo–based computational techniques in overcoming a sign problem. Quantum computing offers a sign-problem-free approach, but evaluating thermal expectation values is generally resource intensive on quantum computers. To facilitate thermodynamic studies of gauge theories, we propose a generalization of the thermal-pure-quantum-state formulation of statistical mechanics applied to constrained gauge-theory dynamics, and numerically demonstrate that the phase diagram of a simple low-dimensional gauge theory is robustly determined using this approach, including mapping a chiral phase transition in the model at finite temperature and chemical potential. Quantum algorithms, resource requirements, and algorithmic and hardware error analysis are further discussed to motivate future implementations. Thermal pure quantum states, therefore, may present a suitable candidate for efficient thermal simulations of gauge theories in the era of quantum computing.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Mapping the kinetic evolution of metastable grain boundaries under non-equilibrium processing

The kinetic evolution of a multiplicity of metastable grain boundaries (GBs) under fast driving conditions are studied by atomistic modeling. Assisted with an enhanced statistical analysis, the energetic evolution of GBs over a broad metastability-temperature space is mapped out, wherein two distinct regimes—an ageing regime and a rejuvenating regime—are retrieved with high fidelity. By comparing the results under various conditions (e.g. random perturbations, isothermal annealing, and fast heating/cooling), it is shown that such ageing/rejuvenating mechanism map is universal, irrespective of the actual stimuli used to elicit the metastable GBs. Here, the ageing/rejuvenating phenomena are demonstrated to stem from the energy imbalance of uphill climbing and downhill dropping during sequential transitions in the system's potential energy landscape. Without the necessity of introducing free parameters, such model can reconcile experimentally measured hardness variation of nanocrystalline metals subjected to femto-second laser irradiation, and it therefore provides a novel perspective on achieving a plurality of interfacial states and facilitating previously inaccessible property regimes.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Direct Monte Carlo Simulations of Chemical Equilibrium Composition of Molecular Fluid Mixtures Under Shock Conditions

A new Monte Carlo simulation method is implemented which holds the number and kind of atoms constant rather than the number and kind of molecules for a molecular fluid mixture. Chemical reactions are allowed as a correlated move of atoms such as an interchange of atoms between molecules. Equilibrium chemical composition is determined as an average over the simulation along with equation of state quantities. Results are presented for N 2 + O 2 → 2NO at 30GPa and 3000K. The equilibrium composition is shown to be sensitive to the potentials between unlike species.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Fuel performance analysis of fully-resolved TRISO compact

The TRi-structural ISOtropic (TRISO) fuel multilayered coating structure offers multiple barriers to fission product release, enhancing safety and performance. The heterogeneous nature of TRISO fuel compacts, comprising thousands of randomly distributed coated fuel particles embedded in a graphite matrix, creates intricate stress fields and thermal gradients that cannot be accurately modeled using simplified one-dimensional or homogenized approaches. Consequently, three-dimensional modeling enables the prediction of fuel compact dimensional changes, internal pressure buildup, and fission product transport pathways under diverse irradiation and thermal conditions. This capability facilitates detailed analysis of particle-to-particle interactions, matrix cracking mechanisms, and the statistical distribution of coating failures, which directly impact fuel performance and safety margins. This capability is particularly critical for advanced reactors, such as high-temperature gas-cooled reactors and other Generation IV reactor designs where TRISO fuel operates at elevated temperatures and burn-up levels. This work introduces a novel method to generate an optimized packing of TRISO compacts and a complete 3D mesh with random distribution of TRISO particles, which are discretized into each coating component layer.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Chemical Understanding of Actinide Separations

This work provides a broad overview of the separation processes used to isolate actinides and the experimentally and computationally determined chemical characteristics that define those separations. The redox chemistry of the actinides plays a pivotal role in both aqueous and pyrochemical processing separations. The near-overlapping energies of the 6d and 5f orbitals in the light actinides allow for facile adjustment of actinide oxidation states, which is used in many established separation methods. In contrast, the stable, generally 3+ oxidation states of the mid- and heavy actinides can make it difficult to separate them from the similarly lanthanides(III). In aqueous separations, the tendency of the actinides to form anionic and neutral aqueous complexes with a variety of complexants (especially soft donors) is used to achieve high separation factors between chemically similar elements in both solid–liquid separations and liquid–liquid extraction. This selectivity can be further tuned through the use of specialized organic or solid-phase ligands. Pyroprocessing separations utilize the unique redox behavior of the actinides to adjust their distribution between a molten salt electrolyte and either a solid electrode or molten metal phase. Atomic-level insights into the mechanisms underlying actinide separation processes, with the ultimate goal of predicting separation behavior, can be provided by electronic structure and statistical mechanical-based calculation methods.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Quantifying the effect of tow architecture variability on the performance of biaxially braided composite tubes

Targeted experiments are performed with stereoscopic digital image correlation to quantify the effect of irregular braid structure on the local surface mechanical response of carbon-fiber reinforced epoxy composite tubes. Virtual models comparing both nonuniform and uniform tow structure of the outer braid of each specimen are integrated into numerical simulations similar to experiments. Experimental and numerical results agree well. Statistical tests are used to support correlations between as-manufactured braid variation ranges to local increases in stress up to 18% compared to an ideal, uniform structure in the elastic regime. Results support improved manufacturing and modeling efforts of defect sensitive braided composites.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Learning turbulent flows with generative models for super resolution and sparse flow reconstruction

Neural operators are promising surrogates for dynamical systems but when trained with standard L 2 losses they tend to oversmooth fine-scale turbulent structures. Here, we show that combining operator learning with generative modeling overcomes this limitation. We consider three practical turbulent-flow challenges where conventional neural operators fail: spatio-temporal super-resolution, forecasting, and sparse flow reconstruction. For Schlieren jet super-resolution, an adversarially trained neural operator (adv-NO) reduces the energy-spectrum error by 15 × while preserving sharp gradients at neural operator-like inference cost. For 3D homogeneous isotropic turbulence, adv-NO trained on only 160 timesteps from a single trajectory forecasts accurately for five eddy-turnover times and offers 114 × wall-clock speed-up at inference than the baseline diffusion-based forecasters, enabling near-real-time rollouts. For reconstructing cylinder wake flows from highly sparse Particle Tracking Velocimetry-like inputs, a conditional generative model infers full 3D velocity and pressure fields with correct phase alignment and statistics. These advances enable accurate reconstruction and forecasting at low compute cost, bringing near-real-time analysis and control within reach in experimental and computational fluid mechanics.

Fluid dynamics↗
DOE OSTI2020

Hydration Mechanisms in Nanoparticle Interaction and Surface Energetics

Work on the project advanced molecular level understanding, prediction and control of nanoscale hydration in salt solutions under electric stimuli or ionic patterning. Permeation of nanoporous electrodes at preset voltage underlies the function of ultracapacitors. Transitory regulation of wetting in nanoporous media by electric field spans an array of applications in materials, energy storage, and separation sciences. The physically related modulation of nanoparticle solubility by surface charges can significantly extend the range of the nanomaterial applications, improve processing techniques, and potentially alleviate environmental concerns. Pore/solution equilibria, and activated kinetics of liquid gating by nanoconfined electrolytes are challenging problems at the forefront of experimental and theoretical research. Addressing these problems from a molecular perspective required the development of state-of-the-art simulation algorithms in statistical mechanics to capture complex processes in open systems under electric control. Parallel studies of wetting and dispersibility of polar and ionizing particles aim to uncover predictive relations between electrowetting, chemical functionalization, and geometry of nanomaterial particles. By nonequilibrium dynamic modeling, we elucidated electrolyte flow in nanochanels and associated electrokinetic energy conversion. Research on the program provided training opportunities for the next generation of scientists in computational chemistry. Insights, and methods from the project implicate broad segment of researchers in energy and nanosciences, materials and surface chemistry.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

2D Ising Model for Adsorption–induced Enantiopurification of Racemates

Mechanisms for the spontaneous transformation of achiral chemical systems into states of enantiomeric purity have important ramifications in modern pharmacology and potential relevance to the origins of homochirality in life on Earth. Such mechanisms for enantiopurification are needed for production of chiral pharmaceuticals and other bioactive compounds. Previously proposed chemical mechanisms leading from achiral systems to near homochirality are initiated by a symmetry–breaking step resulting in a minor excess of one enantiomer via statistical fluctuations in enantiomer concentrations. Subsequent irreversible processes then amplify the majority enantiomer concentration while simultaneously suppressing minority enantiomer production. Herein, equilibrium adsorption of amino acid enantiomer mixtures onto chiral and achiral surfaces reveals amplification of surface enantiomeric excess relative to the gas phase; i. e. enantiopurification of chiral adsorbates by adsorption. Furthermore, this adsorption–induced amplification of enantiomeric excess is shown to be well–describe by the 2D Ising model. More importantly, the 2D–Ising model predicts formation of homochiral monolayers from adsorption of racemic mixtures or prochiral molecules on achiral surfaces; i.e. enantiopurification with no apparent chiral driving force.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Atomistic Mechanics of Torn Back Folded Edges of Triangular Voids in Monolayer WS 2

Triangular nanovoids in 2D materials transition metal dichalcogenides have vertex points that cause stress concentration and lead to sharp crack propagation and failure. Here, the atomistic mechanics of back folding around triangular nanovoids in monolayer WS 2 sheets is examined. Combining atomic-resolution images from annular dark-field scanning transmission electron microscopy with reactive molecular modelling, it is revealed that the folding edge formation has statistical preferences under geometric conditions based on the orientation mismatch. It is further investigated how loading directions and strong interlayer friction, interplay with WS 2 lattice's crack preference, govern the deformation and fracture pattern around folding edges. Finally, these results provide fundamental insights into the combination of fracture and folding in flexible monolayer crystals and the resultant Moiré lattices.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Treating random sequential addition via the replica method

While many physical processes are non-equilibrium in nature, the theory and modeling of such phenomena lag behind theoretical treatments of equilibrium systems. The diversity of powerful theoretical tools available to describe equilibrium systems has inspired strategies that map non-equilibrium systems onto equivalent equilibrium analogs so that interrogation with standard statistical mechanical approaches is possible. In this work, we revisit the mapping from the non-equilibrium random sequential addition process onto an equilibrium multi-component mixture via the replica method, allowing for theoretical predictions of non-equilibrium structural quantities. We validate the above approach by comparing the theoretical predictions to numerical simulations of random sequential addition.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Residential Building Energy Efficiency Field Studies: Low-Rise Multifamily

In recent years, the U.S. Department of Energy (DOE) has conducted a series of research studies to validate energy efficient building technologies in the field. Much of the work has focused on single-family construction, and some has also addressed commercial energy codes. The work detailed in this DOE-funded study (EE0007616) focuses on low-rise multifamily buildings (three stories or fewer above grade) in various regions of the United States, and reports on how state-level building codes are being implemented, both in terms of observed characteristics and also in terms of estimated energy impacts. Nearly 100 buildings across four states—Illinois, Minnesota, Oregon, and Washington—were sampled, which represent a range of climate types from mild temperature to very cold continental. Both common entry and outdoor entry buildings were included, and a parallel research project evaluated envelope air tightness and current still-evolving air tightness testing methods. Finally, a set of structured interviews of building designers and other relevant professionals was carried to out to gain more insight into this market. To the greatest extent possible, the methodology developed under the project for low-rise multifamily buildings mirrored the approach established by Pacific Northwest National Laboratory (PNNL) for single-family residential buildings (https://www.energy.gov/eere/buildings/downloads/residential-building-energy-code-field-study). This included the general approach to sampling, recruitment, and data collection, as well as data analysis and presentation. The range of permitting dates for the sites encompassed two energy code cycles in most regions. All states in the study had adopted a variation of the International Energy Conservation Code (IECC) for the structure of their state code. The low-rise multifamily occupancy presents a hybrid building type: most of the building’s conditioned floor area was covered by the residential chapter of the code while portions of the building (such as corridors and common spaces) fell under the commercial code chapter. The key items assessed in this work were: Building Shell—exterior wall insulation, ceiling insulation, foundation insulation, windows. Common Areas—HVAC and lighting. Living Units—lighting, ventilation. A few items were not assessed in detail, given their relative paucity in this occupancy type; these included duct leakage, pipe insulation, and hot water circulation controls. Building characteristics were collected via a combination of architectural, mechanical, electrical, and plumbing plan reviews and field inspections, and entered into a spreadsheet-based tool that was later queried to build a database. Data went through quality control both upon arrival and via a later semi-automated review and assurance process. Most of the data are presented graphically so that the reader can quickly assess compliance with the applicable energy codes (both by state and by code year). As a final step, EnergyPlus™ simulations were created for all buildings in the study to estimate both the as-found energy use intensity (EUI) and the energy and CO 2 that could be saved if features that were found to not meet code minimums were brought up to code. The savings estimates were tabulated for each of the four states in the study. The research team found that the single-family approach was largely applicable to low-rise multifamily buildings. This applies to both the data collection and the prototype EUI analysis. Most of the occupied space is living units and falls under residential energy codes, and many characteristics use similar envelope construction and relatively straightforward mechanical systems and lighting. One of the most challenging aspects of this work was to build an effective spreadsheet-based data collection instrument that could allow efficient collection of both building plan and field data. The research team is of the view that other methods could be equally effective if the work is done carefully with diligent quality control. The primary findings for the work center around the thermal envelope and mechanical systems and lighting at the sites: For thermal envelope components, the majority of buildings met or were better than the prescriptive code.This suggests that building designers and builders are aware of code requirements. In some cases, surveyed buildings were designed to qualify for energy efficiency certification programs. These buildings made up at least 20% of sampled buildings in each state. Almost all buildings met mechanical system efficiency requirements (for both living units and common areas). In some cases, sites employed systems that were considerably more efficient than required by the applicable energy code. Dwelling units had a majority of high-efficacy lighting, often in excess of the state’s residential code requirements. While high-efficacy fixtures were also typical in common areas (corridors and stairwells), lighting power densities (LPDs) in these areas were sometimes higher than levels dictated by the applicable part of the state commercial energy code. The simulation models run on a series of low-rise multifamily prototypes, informed by a composite of the field data collected, calculated annual EUIs of between 20 and 50 kBtu/ft2-yr, with the range representing the effects of both building characteristics and building location (climate zone). A detailed process (based on simulations of prototype buildings) was used to estimate the amount of avoided energy use that would occur if 100% adherence to energy codes were attained. The results indicated modest savings are attainable for items such as window thermal performance and common area lighting. The result is overall only a modest potential for additional energy savings, averaging about 10% of EUI.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2021

Machine learning to discover mineral trapping signatures due to CO 2 injection

Mineral trapping is pursued as a geological CO 2 sequestration (GCS) mechanism because it permanently stores CO 2 in solid phases or minerals. However, CO 2 mineral-trapping mechanisms are poorly understood due to (1) lack of sufficient field and laboratory data characterizing these complex processes, and (2) challenges to develop site-specific reactive-transport models coupling fluid flow and geochemical reactions occurring at various temporal (from milliseconds to years) and spatial (from pore (millimeters) to field (kilometers)) scales. Reactive transport with additional complexities such as heterogeneity can make the simulation outputs even more difficult to interpret because of complex nonlinearity and multi-scale interdependencies. Furthermore, the values of model outputs such as concentrations can vary by several orders of magnitude, making it harder to correlate and characterize the impact of the variables via traditional data interpretation techniques such as exploratory data analyses. Recently, machine learning (ML) has shown promise in feature discovery and in highlighting hidden mechanisms that cannot be obtained by existing data-analytics and statistical methods. In this study, we applied an unsupervised ML approach, non-negative matrix factorization with custom -means clustering (NMF) to the data generated by reactive-transport simulations of GCS. The reactive-transport data consisted of 19 attributes, including four physio-chemical variables (pH, porosity, aqueous CO 2 , and sequestered CO 2 ), six chemical species (K + , Na + , HCO, Ca 2+ , Mg 2+ , Fe 2+ ), and four carbonate minerals (calcite, dolomite, siderite, and ankerite), a feldspar mineral (albite), and four clay minerals (illite, clinochlore, kaolinite, and smectite) over a period of 200 years of simulation time. Furthermore, the simulation data used was for Morrow B sandstone at the Farnsworth hydrocarbon unit in Texas. Data are sampled at two locations within the model domain: (1) at the injection well and (2) 200 m west of the injection well. The injection was performed for a period of 10 years. Using NMF, we estimated the temporal interdependencies among the 19 attributes over a span of 200 years. We found that NMF was able to identify four reaction stages and their dominant attributes; these cannot be directly discerned through traditional visualization (e.g., line plots, Pareto analysis, Glyph-based visualization methods) or exploratory data analysis tools of the simulation data. The four stages were: reactions in the injection phase followed by short-, mid-, and long-term reactions. The NMF analysis also revealed that 10 among the 19 attributes are dominant. These dominant attributes for mineral trapping include calcite, dolomite at injection well, siderite at 200 m away from the injection well, clinochlore, kaolinite, Na + , K + , Ca 2+ , Mg 2+ , pH, and aqeuous CO 2 . Finally, at late times (65–200 years), our results showed that calcite plays a major role in mineral trapping with insignificant contribution from siderite, ankerite, and clay minerals. These findings make the proposed unsupervised ML-model attractive for reactive-transport sensing towards real-time GCS monitoring.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Bayesian Physics Informed Spatio-Temporal Network for Streamflow Data Imputation

Reliable reconstruction of incomplete streamflow records is critical for improving hydrological forecasting, flood preparedness, and water resource management. However, large observational gaps and uncertainties in governing physical parameters limit the accuracy of traditional statistical and machinelearning imputation frameworks. To address these challenges, we develop a Bayesian Physics-Informed Spatio-Temporal Network (BPI-STNet) that jointly captures spatial and temporal dependencies while enforcing hydrologic consistency through embedded physical constraints. The framework integrates a GraphSAGE-LSTM architecture to model spatial connectivity across gauges and temporal flow dynamics, coupled with a Bayesian update mechanism to estimate uncertain parameters in a simplified water-balance framework. Unlike conventional physics-informed networks that rely on sampling-based posterior estimation, BPI-STNet derives an analytic solution to the inverse problem, allowing closed-form Bayesian updates of uncertain parameters Λ={α,β,k} using Gaussian priors and likelihoods. Applied to daily observations from the Susquehanna River Basin (1980-2022), BPI-STNet achieves substantial improvements over a purely data-driven RGNN baseline, which reduced RMSE by 23 % and MAE by 9 %, and achieving an average NSE values up to 0.96. The results demonstrate that coupling Bayesian inference with physics-informed learning yields physically consistent, uncertainty-aware reconstructions that preserve the temporal persistence and statistical distribution of observed flows. The proposed framework establishes a generalizable paradigm for data-sparse hydrologic systems where both data fidelity and physical interpretability are essential.

Krishnan Kutty Ambika, Anukesh [ORNL] (ORCID:00000↗
DOE OSTI2022

Hierarchical deep reinforcement learning reveals a modular mechanism of cell movement

Time-lapse images of cells and tissues contain rich information about dynamic cell behaviours, which reflect the underlying processes of proliferation, differentiation and morphogenesis. However, we lack computational tools for effective inference. Here we exploit deep reinforcement learning (DRL) to infer cell–cell interactions and collective cell behaviours in tissue morphogenesis from three-dimensional (3D) time-lapse images. We use hierarchical DRL (HDRL), known for multiscale learning and data efficiency, to examine cell migrations based on images with a ubiquitous nuclear label and simple rules formulated from empirical statistics of the images. When applied to Caenorhabditis elegans embryogenesis, HDRL reveals a multiphase, modular organization of cell movement. Imaging with additional cellular markers confirms the modular organization as a novel migration mechanism, which we term sequential rosettes. Furthermore, HDRL forms a transferable model that successfully differentiates sequential rosettes-based migration from others. Our study demonstrates a powerful approach to infer the underlying biology from time-lapse imaging without prior knowledge.

59 BASIC BIOLOGICAL SCIENCES↗