Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “3D 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 1,729 records · Page 96
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Integrated model construction for CO2-EOR monitoring via charged-wellbore casing controlled-source electromagnetics

We present a method for the integrated development of 3D electrical conductivity models for CO2-EOR monitoring with charged wellbore casing controlled-source electromagnetics (CWC-CSEM). The model is constructed through the multiphysics integration of reservoir simulation, seismic, borehole resistivity, and transient EM (TEM) data. The process is performed in two stages. First, a large background conductivity model is constructed from the combination of seismic, borehole, and TEM data. Second, a more detailed and dynamic component of the model is created within the production interval by converting the reservoir simulation parameters to conductivity through Archie’s equation. In this presentation, we demonstrate the complete workflow for building these 3D conductivity site models that can be updated throughout production using the Bell Creek oil field as an example. We then show application of the developed site model to simulating the predicted CWC-CSEM data as one step in the larger problem of reservoir imaging and monitoring of injected CO2 during enhanced oil recovery.

CO2, Sequestration, Monitoring, time-lapse, reserv↗
DOE OSTI2018

Center for Extended Magnetohydrodynamic Modeling

The RPI team efforts focused on providing unstructured mesh technologies to the Center for Extended Magnetohydrodynamic Modeling’s efforts on the development of the M3D-C1 plasma physics code. Specific developments included: Providing the parallel unstructured mesh infrastructure used within the M3D-C1 code; Creating a tool for the generation and subsequent adaptation of the unstructured meshes; A set of procedures to be able to initiate an M3D-C1 simulation in 2D and automatically moving the simulation to a full 3D finite element model when 3D modes begin to appear; Developing a set of physics/mesh based diagonal preconditioners; Supporting the installation of M3D-C1 on various super computer systems; Supplementing the mesh generation procedures to support a model wall region to support the application of an advanced resistive wall model; Providing M3D-C1 with tools to support PIC operations on distributed meshes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Simulating Regional Groundwater Flow in the Vicinity of the Nevada National Security Site, Nye County, Nevada - 16104

The Death Valley Regional Groundwater Flow System (DVRFS) occupies an area of about 70,000 km{sup 2} in south-central Nevada and southeastern California. Interest in this system stems from the need to (1) understand the groundwater flow paths and travel times associated with potential migration of radioactive material from the Nevada National Security Site (NNSS) where underground nuclear testing was conducted from 1956 to 1992; (2) characterize the groundwater system in the vicinity of the proposed high-level radioactive waste repository at Yucca Mountain, Nev.; and (3) address a variety of potential effects from NNSS activities on users down gradient from the NNSS and Yucca Mountain. A single groundwater flow model of the DVRFS has been developed that integrates data and results of two previous regional-scale models and supports Department of Energy (DOE) programs at the NNSS. The numerical model was constrained by organizing hydrogeologic data and interpretations in three-dimensional (3D) geographic information systems and constructing a digital three-dimensional (3D) hydrogeologic framework model (HFM) to represent 27 hydrogeologic units (HGUs) and major structures in the DVRFS region. The 3D hydrogeologic data sets were discretized to 1,500 meters (m) grid cell resolution for input arrays required for the model. The DVRFS was simulated using the USGS 3D groundwater-flow modeling code MODFLOW-2000. Two examples showing how the DVRFS regional model is used by the DOE as part of waste management activities at the NNSS: (1) use of the regional model to provide the boundary conditions for site-scale models, and (2) use of the regional model to assess changes in hydraulic gradient and flow paths in the regional flow system resulting from changes in pumping in various areas of interest. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2020

3D microstructure design of lithium-ion battery electrodes assisted by X-ray nano-computed tomography and modelling

Driving range and fast charge capability of electric vehicles are heavily dependent on the 3D microstructure of lithium-ion batteries (LiBs) and substantial fundamental research is required to optimise electrode design for specific operating conditions. Here we have developed a full microstructure-resolved 3D model using a novel X-ray nano-computed tomography (CT) dual-scan superimposition technique that captures features of the carbon-binder domain. This elucidates how LiB performance is markedly affected by microstructural heterogeneities, particularly under high rate conditions. The elongated shape and wide size distribution of the active particles not only affect the lithium-ion transport but also lead to a heterogeneous current distribution and non-uniform lithiation between particles and along the through-thickness direction. Building on these insights, we propose and compare potential graded-microstructure designs for next-generation battery electrodes. To guide manufacturing of electrode architectures, in-situ X-ray CT is shown to reliably reveal the porosity and tortuosity changes with incremental calendering steps.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2012

Development of 3D full-core ERANOS-2.2/MCNPX-2.7.0 models and neutronic analysis of the BFS-2 zero-power facility

The present paper is addressing the development and validation against experimental data of 3D full-core models of the BFS-2 zero-power fast-reactor using both the deterministic system code ERANOS-2.2 and the stochastic code MCNPX-2.7.0. The model configuration of BFS considered for analysis is the BFS-62-3A benchmark. To extend the - deterministic/stochastic - code-to-code comparison, neutronic parameters, i.e. reactivity, neutron spectrum and reaction rates, were also simulated at the cell level with the Monte Carlo code SERPENT-1.1.7 with two modern data libraries, ENDF-B/VII and JEFF-3.1.1. The BFS-2 critical zero-power facility at the Inst. of Physics and Power Engineering (IPPE) was designed for simulations of the core and shielding of sodium-cooled, fast reactors, for neutron data validation and comparison with experimental results. At the BFS-2 facility, the BFS-62-3A critical benchmark experiment was set-up as a mock-up of the BN-600 reactor, with hybrid MOX fuel and stainless steel reflectors. A UO{sub 2} blanket and a large non-homogeneous stainless-steel reflector surround the core. The lattice is hexagonal of pitch 5.1 cm and metallic dowels are used to keep in central position cylindrical rods made of different types of material (fissile, fertile, blanket, plenum, shielding and absorber). A typical subassembly is formed in piling up various pellets of about 1 cm in height and 4.6 cm in diameter, conferring large heterogeneity in the axial direction. The full-core model development was a complex task due to the large number of subassemblies and the axial subassembly heterogeneity. In ERANOS-2.2, it was necessary to homogenize axially per region the pellets used to form the subassembly. The self-shielded macroscopic cross-sections were calculated using the cell code ECCO in association with JEFF-3.1 and ENDF/B-VI.8 data libraries. The core calculations were performed with broad cross-sections data in 33 neutron energy groups with the solver AVNM in the diffusion approximation, mostly. In MCNPX-2.7.0, a step-by-step approach was used, starting with a model in which the fissile rods were simulated on a homogeneous level, to finally integrate the actual heterogeneous description of the subassemblies. The code-to-code cell analysis performed between ECCO, SERPENT and MCNPX with different modern nuclear data library revealed that the results for the infinite multiplication factor between Monte Carlo and deterministic analysis are in good agreement ({Delta}p < 100 pcm). The differences between the results were observed to be larger for the neutron data libraries, with reactivity differences up to 350 pcm. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

Experimental workflow to estimate model parameters for evaluating long term viscoelastic response of CO2 storage caprocks

Understanding the time-dependent behavior of reservoir and sealing formations is critical to assessing risks associated with geological carbon storage since time-dependent deformation strongly influences mechanical responses of some rock types. Many studies have evaluated the risk of CO2 leakage and induced seismicity by assuming poroelastic rheology in sealing formations. Few have considered viscoelastic or other time-dependent responses, where the existing literature adopts 1D models to represent long-term time-dependent responses. This is primarily because to date, the general form of a reasonable 3D time-dependent model for rocks remains unclear. In this paper, we address this unclear issue by proposing a new workflow to select constitutive modeling parameters to evaluate if a 3D viscoelastic model is reasonable using several-hour-long experimental data and a power-law response to extrapolate to the decades-long time frames of interest in geologic carbon storage. To provide experimental data, we conducted multi-level loading/unloading triaxial relaxation tests with four rock types. The experimental results showed that the maximum load relaxation observed is approximately 49%, with some rock types showing as little as 1.4%. Using a simple linear viscoelastic model, parameters were chosen such that a maximum deviation of 1.5 MPa in axial stress and 7 MPa in radial stress was attained with the extrapolated 30-year data. We found that a reasonable parameter range for the normalized elastic modulus is 0.1~2 for rocks with significant time-dependent responses and 0.01~0.06 for those with small time-dependent responses. No matter how significant time-dependent responses are for rocks considered, our results showed that the relaxation time has a general range of 1~10^10 s, whose time scale can be one or two orders higher than a time frame typically envisioned for CO2 injection projects.

Stress relaxation, 3D time-dependent model, viscoe↗
DOE OSTI2015

Atmospheric heat redistribution and collapse on tidally locked rocky planets

Atmospheric collapse is likely to be of fundamental importance to tidally locked rocky exoplanets, but it remains understudied. Here, general results on the heat transport and stability of tidally locked terrestrial-type atmospheres are reported. First, the problem is modeled with an idealized three-dimensional (3D) general circulation model (GCM) with gray gas radiative transfer. It is shown that over a wide range of parameters that the atmospheric boundary layer, rather than the large-scale circulation, is the key to understanding the planetary energy balance. Through a scaling analysis of the interhemispheric energy transfer, theoretical expressions for the day–night temperature difference and surface wind speed are created that reproduce the GCM results without tuning. Next, the GCM is used with correlated-k radiative transfer to study heat transport for two real gases (CO{sub 2} and CO). For CO{sub 2}, empirical formulae for the collapse pressure as a function of planetary mass and stellar flux are produced, and critical pressures for atmospheric collapse at Earth’s stellar flux are obtained that are around five times higher (0.14 bar) than previous gray gas estimates. These results provide constraints on atmospheric stability that will aid in future interpretations of observations and exoplanet habitability modeling.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

SU-G-IeP2-08: Investigation On Signal Detectability in Volumetric Cone Beam CT Images with Anatomical Background

Purpose: To investigate the slice direction dependent detectability in cone beam CT images with anatomical background. Methods: We generated 3D anatomical background images using breast anatomy model. To generate 3D breast anatomy, we filtered 3D Gaussian noise with a square root of 1/f{sup 3}, and then assigned the attenuation coefficient of glandular (0.8cm{sup −1}) and adipose (0.46 cm{sup −1}) tissues based on voxel values. Projections were acquired by forward projection, and quantum noise was added to the projection data. The projection data were reconstructed by FDK algorithm. We compared the detectability of a 3 mm spherical signal in the image reconstructed from four different backprojection Methods: Hanning weighted ramp filter with linear interpolation (RECON1), Hanning weighted ramp filter with Fourier interpolation (RECON2), ramp filter with linear interpolation (RECON3), and ramp filter with Fourier interpolation (RECON4), respectively. We computed task SNR of the spherical signal in transverse and longitudinal planes using channelized Hotelling observer with Laguerre-Gauss channels. Results: Transverse plane has similar task SNR values for different backprojection methods, while longitudinal plane has a maximum task SNR value in RECON1. For all backprojection methods, longitudinal plane has higher task SNR than transverse plane. Conclusion: In this work, we investigated detectability for different slice direction in cone beam CT images with anatomical background. Longitudinal plane has a higher task SNR than transverse plane, and backprojection with hanning weighted ramp filter with linear interpolation method (i.e., RECON1) produced the highest task SNR among four different backprojection methods. This research was supported by the MSIP (Ministry of Science, ICT and Future Planning), Korea, under the IT Consilience Creative Programs(IITP-2015-R0346-15-1008) supervised by the IITP (Institute for Information & Communications Technology Promotion), Basic Science Research Program through the National Research Foundation of Korea (NRF) funded by the MSIP (2015R1C1A1A01052268) and framework of international cooperation program managed by NRF (NRF-2015K2A1A2067635).

60 APPLIED LIFE SCIENCES↗
DOE OSTI2024

Cellular response of keratinocytes to the entry and accumulation of nanoplastic particles

Plastic accumulation in the environment is rapidly increasing, and nanoplastics (NP), byproducts of environmental weathering of bulk plastic waste, pose a significant public health risk. Particles may enter the human body through many possible routes such as ingestion, inhalation, and skin absorption. However, studies on NP penetration and accumulation in human skin are limited. Loss or reduction of the keratinized skin barrier may enhance the skin penetration of NPs. The present study investigated the entry of NPs into a human skin system modeling skin with compromised barrier functions and cellular responses to the intracellular accumulations of NPs. Two in vitro models were employed to simulate human skin lacking keratinized barriers. The first model was an ex vivo human skin culture with the keratinized dermal layer (stratum corneum) removed. The second model was a 3D keratinocyte/dermal fibroblast cell co-culture model with stratified keratinocytes on the top and a monolayer of skin fibroblast cells co-cultured at the bottom. The penetration and accumulation of the NPs in different cell types were observed using fluorescent microscopy, confocal microscopy, and cryogenic electron microscopy (cryo-EM). The cellular responses of keratinocytes and dermal fibroblast cells to stress induced by NPs stress were measured. The genetic regulatory pathway of keratinocytes to the intracellular NPs was identified using transcript analyses and KEGG pathway analysis. The cellular uptake of NPs by skin cells was confirmed by imaging analyses. Transepidermal transport and penetration of NPs through the skin epidermis were observed. According to the gene expression and pathway analyses, an IL-17 signaling pathway was identified as the trigger for cellular responses to internal NP accumulation in the keratinocytes. The transepidermal NPs were also found in co-cultured dermal fibroblast cells and resulted in a large-scale transition from fibroblast cells to myofibroblast cells with enhanced production of α-smooth muscle actin and pro-Collagen Ia. The upregulation of inflammatory factors and cell activation may result in skin inflammation and ultimately trigger immune responses.

36 MATERIALS SCIENCE↗
NASA NTRS2009

Modeling the Strahl in the Solar Wind Using CLUSTER/PEACE Electron Distribution Functions

We present a study of kinetic properties of the strahl electron velocity distribution functions (VDF's) in the solar wind. These are used to investigate the pitch-angle scattering and stability of the population to interactions with electromagnetic (whistler) fluctuations. The study is based on high time resolution data from the Cluster/PEACE electron spectrometer. Our study focuses on the mechanisms that control and regulate the pitch-angle and stability of strahl electrons in the solar wind; mechanisms that are not yet well understood. Various parameters are investigated such as the strahl-electron density, temperature anisotropy, and electron heat-flux. The goal is to check whether the strahl electrons are constrained by some instability (e.g., the whistler instability), or are maintained by other types of processes. The electron heat-flux and temperature anisotropy are determined by modeling the 3D-VDFs to a spectral spherical harmonic model from which the moments are obtained directly from the spectral coefficients.

Vinas, Adolfo F.↗
NASA NTRS2010

Dst Index in the 2008 GEM Modeling Challenge - Model Performance for Moderate and Strong Magnetic Storms

The GEM 2008 modeling challenge efforts are expanding beyond comparing in-situ measurements in the magnetosphere and ionosphere to include the computation of indices to be compared. The Dst index measures the largest deviations of the horizontal magnetic field at 4 equatorial magnetometers from the quiet-time background field and is commonly used to track the strength of the magnetic disturbance of the magnetosphere during storms. Models can calculate a proxy Dst index in various ways, including using the Dessler-Parker Sckopke relation and the energy of the ring current and Biot-Savart integration of electric currents in the magnetosphere. The GEM modeling challenge investigates 4 space weather events and we compare models available at CCMC against each other and the observed values of Ost. Models used include SWMF/BATSRUS, OpenGGCM, LFM, GUMICS (3D magnetosphere MHD models), Fok-RC, CRCM, RAM-SCB (kinetic drift models of the ring current), WINDMI (magnetosphere-ionosphere electric circuit model), and predictions based on an impulse response function (IRF) model and analytic coupling functions with inputs of solar wind data. In addition to the analysis of model-observation comparisons we look at the way Dst is computed in global magnetosphere models. The default value of Dst computed by the SWMF model is for Bz the Earth's center. In addition to this, we present results obtained at different locations on the Earth's surface. We choose equatorial locations at local noon, dusk (18:00 hours), midnight and dawn (6:00 hours). The different virtual observatory locations reveal the variation around the earth-centered Dst value resulting from the distribution of electric currents in the magnetosphere during different phases of a storm.

Rastaetter, Lutz↗
DOE OSTI2012

Statistical Modeling of the Eye for Multimodal Treatment Planning for External Beam Radiation Therapy of Intraocular Tumors

Purpose: Ocular anatomy and radiation-associated toxicities provide unique challenges for external beam radiation therapy. For treatment planning, precise modeling of organs at risk and tumor volume are crucial. Development of a precise eye model and automatic adaptation of this model to patients' anatomy remain problematic because of organ shape variability. This work introduces the application of a 3-dimensional (3D) statistical shape model as a novel method for precise eye modeling for external beam radiation therapy of intraocular tumors. Methods and Materials: Manual and automatic segmentations were compared for 17 patients, based on head computed tomography (CT) volume scans. A 3D statistical shape model of the cornea, lens, and sclera as well as of the optic disc position was developed. Furthermore, an active shape model was built to enable automatic fitting of the eye model to CT slice stacks. Cross-validation was performed based on leave-one-out tests for all training shapes by measuring dice coefficients and mean segmentation errors between automatic segmentation and manual segmentation by an expert. Results: Cross-validation revealed a dice similarity of 95% {+-} 2% for the sclera and cornea and 91% {+-} 2% for the lens. Overall, mean segmentation error was found to be 0.3 {+-} 0.1 mm. Average segmentation time was 14 {+-} 2 s on a standard personal computer. Conclusions: Our results show that the solution presented outperforms state-of-the-art methods in terms of accuracy, reliability, and robustness. Moreover, the eye model shape as well as its variability is learned from a training set rather than by making shape assumptions (eg, as with the spherical or elliptical model). Therefore, the model appears to be capable of modeling nonspherically and nonelliptically shaped eyes.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2019

3D Moment Tensor Inversion of Underground Chemical Explosions from the Source Physics Experiments

Several physical mechanisms have been proposed to explain the generation of S-waves from underground explosions, such as asymmetries in the source, release of tectonic pre-stress, interactions with the free-surface, spall, and heterogeneities in the Earth. An accurate description of the explosion source processes is an important step towards understanding which of these plausible mechanisms are actively contributing to the generation of S-waves and under what conditions. In this study we explore the application of the seismic moment tensor source to model far-field, low frequency (up to 6 Hz) waveform data of over-buried chemical explosions from the Source Physics Experiment, with a focus on S-wave generation and amplitude predictions. We use an inverse waveform modeling approach to estimate the source properties of the chemical explosions, and compare solutions using different velocity models. 1D and 3D subsurface velocity models are used to characterize wave propagation between the source and receiver. We also performed analysis on wavefield simulations from physic-based explosion source modeling. The analyses show scattering and phase conversion from 3D heterogonies dominate the generation of far-field, S-wave energy observed in data, and that the variability in the recovered deviatoric component of the moment tensor source model are largely a result of inadequately accounting for 3D wave propagation effects in the inversion process.

58 GEOSCIENCES↗
DOE OSTI2021

Three-dimensional paganica fault morphology obtained from hypocenter clustering (L'Aquila 2009 seismic sequence, Central Italy)

In seismic modelling, fault planes are normally assumed to be flat due to the lack of data which can constrain fault morphology. However, incorporating 3D fault morphology is important for modelling several phenomena, for example calculating mainshock induced stress changes. Here we utilize a data-analytical method to unveil the 3D rupture morphology of faults using unsupervised clustering techniques applied to earthquake hypocenters in seismic sequences. We apply this method to the 2009 L'Aquila seismic sequence which involved a M W 6.1 mainshock on April 6th. We use a dataset of about 50,000 relocated events, mostly microearthquakes, reaching magnitude of completeness equal to 0.7. Clustering distinguishes the earthquakes as occurring in three main clusters along with other minor fault segments. We then represent the morphology of the main Paganica fault system (responsible for the largest mainshock) using splines. This method shows promise as a step toward robustly and quickly obtaining 3D rupture morphologies where earthquake sequences have been monitored. The 3D model is presented interactively online, and the processing is presented in an interactive Jupyter Notebook (https://bit.ly/2MnCFdj).

58 GEOSCIENCES↗
DOE OSTI2025

Energy Conversion and Electron Acceleration and Transport in 3D Simulations of Solar Flares

Recent observations and simulations indicate that solar flares undergo extremely complex 3D evolution, making 3D particle transport models essential for understanding electron acceleration and interpreting flare emissions. In this study, we investigate this problem by solving Parker’s transport equation with 3D MHD simulations of solar flares. By examining energy conversion in the 3D system, we evaluate the roles of different acceleration mechanisms, including reconnection current sheet (CS), termination shock (TS), and supra-arcade downflows (SADs). We find that large-amplitude turbulent fluctuations are generated and sustained in the 3D system. The model results demonstrate that a significant number of electrons are accelerated to hundreds of keV and even a few MeV, forming power-law energy spectra. These energetic particles are widely distributed, with concentrations at the TS and in the flare looptop region, consistent with results derived from recent hard X-ray (HXR) and microwave (MW) observations. By selectively turning particle acceleration on or off in specific regions, we find that the CS and SADs effectively accelerate electrons to several hundred keV, while the TS enables further acceleration to MeV. However, no single mechanism can independently account for the significant number of energetic electrons observed. Instead, the mechanisms work synergistically to produce a large population of accelerated electrons. Our model provides spatially and temporally resolved electron distributions in the whole flare region and at the flare footpoints, enabling synthetic HXR and MW emission modeling for comparison with observations. These results offer important insights into electron acceleration and transport in 3D solar flare regions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Mineralized sclerites in the gorgonian coral Leptogorgia chilensis as a natural jamming system

The soft corals (Cnidaria, Octocorallia), a diverse group of colonial marine invertebrates, can reversibly tune their body stiffness in response to external stimuli. This capability is attributed to their dynamic skeletal systems, which consist of thousands of mineralized skeletal elements, called sclerites, embedded within a gel-like matrix that swells/deswells and unjams/jams the sclerites, thus modulating skeletal stiffness. While sclerite morphology is widely used for species identification, its role in the mechanical performance of a soft coral’s skeletal system is largely unknown. Here, we investigated structure-jamming relationships in sclerite-based skeletal architectures using the red gorgonian octocoral Leptogorgia chilensis as a model system. The sclerites of L. chilensis exhibit a shaft-like geometry with two axial branches and two sets of triradiate side branches, which are aligned with the crystallographic symmetry of the constituent magnesium-containing calcite. By combining multiscale three-dimensional (3D) structural characterization, parametric geometrical modeling, 3D printing, mechanical testing, and discrete element simulations, we demonstrate how sclerite geometry achieves a balanced jamming performance in terms of stiffness, weight, strength, and fracture resistance in comparison to alternative geometries parametrically modified from the native sclerites (e.g., changes in the length and number of side branches). Here, we also found that these performance metrics are achieved through the effective interlocking among side and axial branches, which is further enhanced by the fractal-like microscopic spikes on the branch tips. The findings in this natural jamming system offer insights for designing synthetic mechanotunable material architectures for a wide range of applications, from soft robotics to mechanical dampeners.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Formation and evolution of protostellar accretion discs – II. From 3D simulation to a simple semi-analytic model of Class 0/I discs

Here we use a 3D radiative non-ideal magnetohydrodynamic simulation to investigate the formation and evolution of a young protostellar disc from a magnetized pre-stellar core. The simulation covers the first ${\sim }10\, {\rm kyr}$ after protostar formation and shows a massive, weakly magnetized disc with radius that initially grows and then saturates at ${\sim }30\, {\rm au}$. The disc is gravitationally unstable with prominent large-amplitude spiral arms. We use our simulation results and a series of physical arguments to construct a predictive and quantitative physical picture of Class 0/I protostellar disc evolution from several aspects, including (i) the angular-momentum redistribution in the disc, self-regulated by gravitational instability to make most of the disc marginally unstable; (ii) the thermal profile of the disc, well-approximated by a balance between radiative cooling and accretion heating; and (iii) the magnetic-field strength and magnetic-braking rate inside the disc, regulated by non-ideal magnetic diffusion. Using these physical insights, we build a simple 1D semi-analytic model of disc evolution. We show that this 1D model, when coupled to a computationally inexpensive simulation for the evolution of the surrounding pseudo-disc, can be used reliably to predict disc evolution in the Class 0/I phase. The predicted long-term evolution of disc size, which saturates at ${\sim }30\, {\rm au}$ and eventually shrinks, is consistent with a recent observational survey of Class 0/I discs. Such hierarchical modelling of disc evolution circumvents the computational difficulty of tracing disc evolution through Class 0/I phase with direct, numerically converged simulations.

79 ASTRONOMY AND ASTROPHYSICS↗