Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “computational modeling and simulations”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Quantifying the Impact of Vertical Resolution on the Representation of Marine Boundary Layer Physics for Global-Scale Models

Abstract While GCM horizontal resolution has received the majority of scale improvements in recent years, ample evidence suggests that a model’s vertical resolution exerts a strong control on its ability to accurately simulate the physics of the marine boundary layer. Here we show that, regardless of parameter tuning, the ability of a single-column model (SCM) to simulate the subtropical marine boundary layer improves when its vertical resolution is improved. We introduce a novel objective tuning technique to optimize the parameters of an SCM against profiles of temperature and moisture and their turbulent fluxes, horizontal winds, cloud water, and rainwater from large-eddy simulations (LES). We use this method to identify optimal parameters for simulating marine stratocumulus and shallow cumulus. The novel tuning method utilizes an objective performance metric that accounts for the uncertainty in the LES output, including the covariability between model variables. Optimization is performed independently for different vertical grid spacings and value of time step, ranging from coarse scales often used in current global models (120 m, 180 s) to fine scales often used in parameterization development and large-eddy simulations (10 m, 15 s). Uncertainty-weighted disagreement between the SCM and LES decreases by a factor of ∼5 when vertical grid spacing is improved from 120 to 10 m, with time step reductions being of secondary importance. Model performance is shown to converge at a vertical grid spacing of 20 m, with further refinements to 10 m leading to little further improvement. Significance Statement In successive generations of computer models that simulate Earth’s atmosphere, improvements have been mainly accomplished by reducing the horizontal sizes of discretized grid boxes, while the vertical grid spacing has seen comparatively lesser refinements. Here we advocate for additional attention to be paid to the number of vertical layers in these models, especially in the model layers closest to Earth’s surface where climatologically important marine stratocumulus and shallow cumulus clouds reside. Our experiments show that the ability of a one-dimensional model to represent physical processes important to these clouds is strongly dependent on the model’s vertical grid spacing.

Meteorology & Atmospheric Sciences↗
DOE OSTI2023

Computational model development and validation of fuel dispersal phenomena

As the US nuclear industry is renewing efforts to extend the rod average burnup limits in existing pressurized water reactors, it becomes imperative to investigate the safety implications associated with high burnup fuel fragmentation and the potential dispersal of fragmented fuel into the reactor coolant system. Additionally, as there is a growing interest in transitioning to new nuclear fuel designs to accommodate higher fuel burnup levels, it becomes essential to conduct studies on fuel behavior during design basis accident scenarios, particularly during a postulated loss-of-coolant accident, in order to ensure the safe operation of light-water reactors. The present study develops a computational model to simulate the complex three-phase flow of the fuel dispersal phenomena that occur following a breach in the cladding, when the stream of solid fuel particles and carrier fission gas stream may interact with the surrounding liquid medium. Furthermore, the proposed framework considers the solid phase as discrete Lagrangian particles while the interpenetrating gas–liquid continuum is modeled using a Eulerian framework. The developed simulation framework is validated with experimental results performed in a separate-effect test facility by comparing predicted particle settlement locations with the experimental observations reported in the literature. Key highlights of this study include the interpretation of the high-pressure boundary conditions and the three-phase flow coupling strategy. Also discussed are the gas–liquid hydrodynamics and the behavior of particles as they are transported through the gas and liquid phases.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

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

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

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Diffraction properties of cylindrically bent KAP crystals in energy range of 2.3–7.5 keV using synchrotron radiation

Verification of physics models and computer simulations are heavily reliant upon the accuracy of experimental measurements. Calibration of instrument responses becomes an important step to achieve this goal. This paper presents systematic studies of bent potassium acid phthalate (KAP) crystals using Lawrence Berkeley National Laboratories, Advanced Light Source, beamline 9.3.1 in the energy range of 2.3 to 7.5 keV. A set of KAP crystals, gradually bent from flat up to a 50.8 mm cylindrical curvature. The measured integrated reflectivity for this set of KAP crystals shows good agreement with the X-ray Oriented Program (XOP) calculations when adjusting the Debye–Waller temperature factor and using the multilamellar model in the calculations. Significant differences in rocking curve profiles were observed between experimental measurements and theory. A forward convolution model and software code were developed to include experimental parameters, allowing the investigation of the difference between measurements and calculations. After considering the experimental parameters, good agreements were obtained for the rocking curve profiles for all bending radii with a unique set of parameters. Our results show that XOP can be a useful and reliable tool to predict performance of cylindrically bent KAP crystals in this energy range.

Wu, Ming↗
DOE OSTI2021

SISGR: Chemomechanics of Far-From-Equilibrium Interfaces (COFFEI)

Portable, reliable, and deployable devices for energy storage and conversion require fundamental changes in design of solid-state composites comprising ceramics and metals. These materials comprise the electrodes and electrolytes of next-generation solid-oxide fuel cells and all solid-state batteries, forming solid-state functional composites. The advent of solid-state batteries – which replace liquid electrolytes with solid electrolytes capable of lithium ion transport for reliable energy storage in portable batteries – and the increased demand space for all solid-state fuel cells capable of oxygen reduction at intermediate temperatures remain important challenges for improved material stability and decreased system cost. However, little is understood about three fundamental facets of materials that enable such solid-state energy applications. First, how do such materials deform, fracture, or delaminate under operando conditions? Second, how does such mechanical deformation limit or facilitate electronic and ionic transport within and across such material interfaces? Third, how we can predictably design interface-rich composites to engineer both structural and electrochemical stability? This COFFEI Group comprised expertise from Materials Science & Engineering and Nuclear Science & Engineering to integrate unique in situ experiments, simulations, and fabricated interfaces that address these fundamental questions in solid-state interfaces of nanoscale composites that will guide solid-state electrochemistry, transport kinetics, and mechanical deformation for nonstoichiometric materials that enable such applications. In particular, we built on COFFEI’s understanding of chemomechanical coupling among defect concentrations, ionic transport, electron transport, and stored elastic energy that is particularly acute in the far-from-equilibrium conditions typical of energy device applications. By tailoring our focus to solid-state interfaces, we addressed these important issues by (a) developing and applying advanced in situ and ex situ characterization tools to characterize model materials and interfaces synthesized with molecular-level control, under both laboratory-controlled and extreme environments representative of energy device operation; and (b) employing computational modeling and simulation frameworks to predict transport mechanisms, reactivity and stability of these model materials and interfaces under significant chemical strains typical of energy device operation. Recent progress provided insights to additional materials systems and electrochemomechanical fatigue and fracture that were not fully envisioned when the program was initiated. Specifically, in the final three years of COFFEI we pursued two integrated thrusts, with complementary focus. Thrust I focused on failure-resistant electrochemomechanical composites, while Thrust II focused on strain-modulated conductivity and reactivity across interfaces. In contrast to our initial COFFEI focus, these thrusts concentrated wholly on solid-state material interfacial interactions and included greater integration of multiscale visualization including in situ electron microscopy of strained structures/interactions and mesoscale simulations. Successful development of functionally superior and long-lived battery and fuel cell systems and stress adaptable oxides requires a deeper, fundamental understanding of the coupling among the historically important subfields of solid-state electrochemistry, transport kinetics, and mechanical deformation for nonstoichiometric metal oxide electrodes. In this program, the understanding and the application of chemomechanical coupling of defect concentrations, ionic transport, electro-catalytic activity and stored elastic energy, particularly acute in the far-from-equilibrium conditions typical of energy device applications, are being refined and implications for device operation clarified, including for miniaturized solid-state batteries and fuel cells.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Spatially distributed infection increases viral load in a computational model of SARS-CoV-2 lung infection

A key question in SARS-CoV-2 infection is why viral loads and patient outcomes vary dramatically across individuals. Because spatial-temporal dynamics of viral spread and immune response are challenging to study in vivo, we developed Spatial Immune Model of Coronavirus (SIMCoV), a scalable computational model that simulates hundreds of millions of lung cells, including respiratory epithelial cells and T cells. SIMCoV replicates viral growth dynamics observed in patients and shows how spatially dispersed infections can lead to increased viral loads. The model also shows how the timing and strength of the T cell response can affect viral persistence, oscillations, and control. By incorporating spatial interactions, SIMCoV provides a parsimonious explanation for the dramatically different viral load trajectories among patients by varying only the number of initial sites of infection and the magnitude and timing of the T cell immune response. When the branching airway structure of the lung is explicitly represented, we find that virus spreads faster than in a 2D layer of epithelial cells, but much more slowly than in an undifferentiated 3D grid or in a well-mixed differential equation model. These results illustrate how realistic, spatially explicit computational models can improve understanding of within-host dynamics of SARS-CoV-2 infection.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Platinum nanoparticle compression: Combining in situ TEM and atomistic modeling

The mechanical behavior of nanoparticles governs their performance and stability in many applications. However, the small sizes of technologically relevant nanoparticles, with diameters in the range of 10 nm or less, significantly complicate experimental examination. These small nanoparticles are difficult to manipulate onto commercial test platforms and deform at loads that are below the typical noise floor of the testing instruments. Here, we synthesized small platinum nanoparticles directly onto a mechanical tester and used a modified nanomanipulator to enhance load resolution to the nanonewton scale. We demonstrated the in situ compression of an 11.5-nm platinum nanoparticle with simultaneous high-resolution measurements of load and particle morphology. Molecular dynamics simulations were performed on similarly sized particles to achieve complementary measurements of load and morphology, along with atomic resolution of dislocations. The experimental and simulation results revealed comparable values for the critical resolved shear stress for failure, 1.28 and 1.15 GPa, respectively. Altogether, this investigation demonstrated the promise of, and some initial results from, the combination of atomistic simulations and in situ experiments with an unprecedented combination of high spatial resolution and high load resolution to understand the behavior of metal nanoparticles under compression.

42 ENGINEERING↗
DOE OSTI2022

Sensitivity Analysis of an Electrospray Dehumidification System

Electrically enhanced condensation achieved by an electrospray system can enhance dehumidification in HVAC applications if successfully scaled up. The first step towards scale-up is identifying critical design and operating parameters that significantly impact dehumidification for the electrospray system. Sensitivity analysis is one approach for identifying critical parameters. This paper presents a methodology that can efficiently perform quantitative sensitivity analysis on a computationally expensive high fidelity simulation model, which is a computational fluid dynamics (CFD) model of the electrospray system. The methodology employs surrogate modeling and Sobol sensitivity analysis. We estimate the main and total effect Sobol sensitivity indices of four operating parameters of the electrospray system. The results reveal that the most important parameter is the size of the sprayed droplets, followed by the volumetric flow rate of the bulk air and injected water in the system. The water injection velocity of the spray was found to have little to no impact on dehumidification.

Young, David↗
DOE OSTI2020

A generalized approximate control variate framework for multifidelity uncertainty quantification

We describe and analyze a variance reduction approach for Monte Carlo (MC) sampling that accelerates the estimation of statistics of computationally expensive simulation models using an ensemble of models with lower cost. These lower cost models — which are typically lower fidelity with unknown statistics — are used to reduce the variance in statistical estimators relative to a MC estimator with equivalent cost. We derive the conditions under which our proposed approximate control variate framework recovers existing multifidelity variance reduction schemes as special cases. We demonstrate that existing recursive/nested strategies are suboptimal because they use the additional low-fidelity models only to efficiently estimate the unknown mean of the first low-fidelity model. As a result, they cannot achieve variance reduction beyond that of a control variate estimator that uses a single low-fidelity model with known mean. However, there often exists about an order-of-magnitude gap between the maximum achievable variance reduction using all low-fidelity models and that achieved by a single low-fidelity model with known mean. We show that our proposed approach can exploit this gap to achieve greater variance reduction by using non-recursive sampling schemes. The proposed strategy reduces the total cost of accurately estimating statistics, especially in cases where only low-fidelity simulation models are accessible for additional evaluations. Several analytic examples and an example with a hyperbolic PDE describing elastic wave propagation in heterogeneous media are used to illustrate the main features of the methodology.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

A linked-scale coupled model of mass erosion and redistribution in plasma-exposed micro-foam surfaces

Surface evolution due to exposure to harsh environments is of importance in many scientific and technological applications. In plasma-exposed materials, the surface receives charged particles from the plasma, leading to a series of processes that drive the system far from equilibrium and may lead to the severe deterioration of the surface properties. Although surface morphological changes are driven by atomic collisions taking place over picoseconds and nanometers, these processes result in mass loss and redistribution of matter over much larger length and time scales. This necessitates a multi-scale approach capable of capturing the range of processes linking primary atomic collision events with engineering-level surface geometry changes. In this paper, we develop a computational model to simulate the morphological evolution and effective erosion rate of micro-architected tungsten foams during low-energy plasma ion bombardment. Furthermore, the model acts on several length scales, with the energy and angular dependence of the sputtering yield of flat tungsten surfaces determined using the SRIM code based on the binary collision approximation. This information is introduced into a low-fluence, short-term Monte Carlo raytracing model, and further into a high-fluence, long-term particle transport model. In the latter, material particles representing billions of atoms are described in a digitized 3-D representation of the foam structure from X-ray tomography data, and are sputtered off and redeposited using the atomistic information. We show that the redeposition of sputtered atoms leads to partial self-healing in the bottom layers with a sharp reduction in the sputtering coefficient of low density foams, roughly 25% of the solid W value. This is in qualitative agreement with recent experiments on low porosity W structures. At high fluence, the foam structure degrades considerably as there are fewer ligaments available to recapture sputtered atoms and, consequently, the sputtering rate increases again.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Two-layer cathode architecture for high-energy density and high-power density solid state batteries

Solid state batteries with high-energy density and high-power density require the development of thick and energy dense cathodes. Structured cathode electrodes with a double-layer configuration were enabled using a freeze tape casting technique. A bottom dense layer was utilized to enhance the energy density whereas a top porous layer with vertically aligned walls was utilized to enhance the power density. The porous structure of the power layer was achieved by ice templating this layer on top of the densified energy layer of the cathode. This configuration was found to better utilize the active material of the cathodes and was optimized using numerical simulation and computer modeling. Cells with Li metal anode and LiNi 0.6 Mn 0.2 Co 0.2 O 2 (NMC622) at approximately 5 and 20 mg/cm 2 were cycled at 70 °C at different C-rates. Poly(ethylene oxide) (PEO) with lithium bis-trifluoromethanesulfonimide (LiTFSI) was used for the catholyte and the solid-state electrolyte. The structured cathodes exhibited more than double capacity values as well as better Coulombic efficiency compared to non-structured (single-layer) thick cathodes. In conclusion, synchrotron X-ray tomography and scanning electron microscopy were used to characterize the microstructure of the cathodes.

25 ENERGY STORAGE↗
DOE OSTI2025

Temperature sensitivity of the equilibrium neutronics and accident analysis of the HTR-10

Pebble-Bed High-Temperature Gas-cooled Reactors (PB-HTGR) are moderated by the graphite in the fuel pebbles and the graphite reflector surrounding the pebble-bed. Because graphite is by far the most abundant material in PB-HTGRs and the primary moderator, accurate modeling of the graphite material, including density, impurities, and temperatures, is crucial for accurate computational modeling and simulation of these reactors. While main characteristics of the graphite components are often known, the local temperature is less well known and often averaged over all components. Here, this work studies the impact of considering accurate temperature profiles in the graphite material on the generation of a small PB-HTGR model at the state of equilibrium operation and on short-term accident progression. The fuel compositions for the PB-HTGR were determined using a jump-in equilibrium modeling method, the Axial Radial Zone Equilibrium Modeling (AR-ZEM) method. In contrast to previous work, the AR-ZEM method was used considering thermal-hydraulic feedback from the MELCOR code to determine temperatures of the fuel pebbles and the surrounding graphite reflector. The consideration of an axial and radial temperature profile in the core and reflector, as opposed to uniform material temperatures, had an impact of almost 1,300 pcm on the equilibrium core eigenvalue and caused significant differences in the discharged plutonium fuel inventory with up to 4.9% and 11.0% for Pu-239 and Pu-242, respectively. To assess the impact on short-term accident progression, two Anticipated Transient Without SCRAM (ATWS) events, a Pressurized Loss of Forced Coolant (PLOFC) and a Control Rod Withdrawal (CRW) with loss of flow, were simulated with MELCOR. The use of temperature profiles in the equilibrium core models did not reveal a significant impact on the temperature, power, or reactivity responses during the transients. In conclusion, a need for consideration of accurate temperature profiles, in particular for the graphite reflector, was found for the generation of equilibrium PB-HTGRs core models using jump-in methods, but detailed temperature profiles may not be necessary when performing conservative transient analysis.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2025

Modeling High Current Pulsed Discharge in AA Battery Cathodes: The Effect of Localized Charging during Rest

During high current operation, substantial heterogeneity develops within battery cathodes, particularly when their thickness is large. Heterogeneity relaxation during subsequent rest is important for understanding battery performance under pulsed conditions. Localized charge balancing phenomena within batteries at zero net current are not well understood and merit investigation. In this work, the heterogeneity within cathodes of commercial alkaline Zn–MnO 2 batteries is measured during discharge and monitored during rest using energy dispersive X-ray diffraction (EDXRD). Significant gradients in protonation form during discharge and partially relax under rest. It is demonstrated that the proton gradient relaxation is through local redox activity at zero net current, where local (de)protonation works to redistribute charge across the cathode thickness. To support this redox-based relaxation, a fundamental kinetic study on prismatic MnO 2 cathodes is conducted to determine an appropriate model to describe both discharge and charge kinetics of MnO 2 . These kinetics are incorporated into a computational model to simulate the proton gradient formation and partial relaxation under identical discharge conditions as the operando EDXRD experiments. Model and experimental data are found to be in excellent agreement, correctly predicting localized charge balancing at rest.

Batteries↗
DOE OSTI2021

Contribution of the Two Types of Ekman Pumping Induced Eddy Heat Flux to the Total Vertical Eddy Heat Flux

Based on eddy-rich (0.1° horizontal resolution for ocean) Community Earth System Model simulations, we compute the Ekman pumping induced vertical eddy heat flux using two formulations: the classical and the Stern-Ekman pumping. The two mechanisms yield similar patterns in the eddy-induced vertical heat flux, showing strong upward heat flux in ocean frontal regions such as the western boundary current regions. Using the Gulf Stream Extension region as an example, we estimate (1) the Stern-Ekman pumping induced eddy heat flux is 23% and 12% of the total vertical eddy heat flux at 50 m depth during summer and winter, respectively; (2) the classical Ekman pumping accounts for 82% and 88% of the Stern-Ekman pumping induced eddy heat flux at 50 m during summer and winter, respectively. The second finding indicates eddy current feedback on wind stress is the primary cause for the upward eddy heat flux generated by Ekman pumping.

58 GEOSCIENCES↗
DOE OSTI2026

Explainable machine learning for incipient anomaly detection in compact molten salt heat exchanger with overlapping feature distributions

High-temperature molten salt-cooled reactors (MSCRs) are a promising next-generation nuclear technology option, offering efficient power conversion and inherent safety features. However, the reliability of these systems depends on the robust operation of heat exchangers (HXs), which are susceptible to failure due to temperature gradients and channel plugging caused by fluid freezing. Conventional monitoring methods, relying on inlet and outlet measurements, lack the spatial resolution needed to detect early-stage faults. We propose a novel design of a compact salt-to-salt matrix-type HX design consisting of interleaved arrays of parallel tubes, with integrated synthetic fiber optic distributed temperature sensing (DTS) to enable localized detection of incipient faults. To evaluate performance of this design, we generate high-fidelity synthetic data using heat transfer computational modeling to simulate channel plugging, and introduce sensor noise for realistic modeling of measurements. The dataset comprises of 97% normal operation and 3% anomaly cases, with each anomaly class representing 1% of the data. These early anomalies result in overlapping temperature profiles between normal and faulty channels, producing a non-separable dataset that challenges traditional classification techniques. We benchmark eight supervised machine learning (ML) models and demonstrate that XGBoost achieves the highest performance. To improve transparency, we develop an explainability framework combining Shapley values and partially ordered sets (POSETs) to quantify and structurally analyze feature importance. This approach identifies both dominant predictors and ambiguous feature relationships, enhancing trust and interpretability. Our results highlight the potential of combining DTS and explainable ML with intelligent feature selection to improve predictive maintenance and ensure operational resilience in advanced nuclear systems.

Prantikos, Konstantinos [Argonne National Laborato↗
DOE OSTI2022

Bibliometric review and recent advances in total scattering pair distribution function analysis: 21 years in retrospect

Global research activities have been driven by the quest to develop and characterize novel materials for technological advancements. The total scattering pair distribution function (TSPDF) is a powerful and versatile characterization technique for examining the structural details of diverse complex materials including liquid, amorphous, disordered crystalline, and nanostructured materials. Thus, it is critical to keep track of research progress, identify research gaps, and future research directions of the application of the TSPDF technique in materials development and discovery. In this work, a bibliometric analysis of literature regarding the TSPDF technique between 2000 and 2021 was conducted using datasets retrieved from the Web of Science database. The research trends based on publication outputs, research subject distribution, co-authorships among institutions, countries/regions, co-citation of referenced sources, and keyword co-occurrence are evaluated and discussed herein. The impact of the TSPDF technique is projected to increase due to its importance in probing emerging functional materials, and the advances in specialized facilities and instrumentation among the scientific communities engaged with it. Finally, current and emerging research hotspots related to TSPDF technique such as catalysis, computer modeling and simulation, pharmaceutics, machine learning, hydrogen storage, battery materials, and layered structured materials are also identified and discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2024

In situ tumor model for longitudinal in silico imaging trials

Abstract Objective.In this article, we introduce a computational model for simulating the growth of breast cancer lesions accounting for the stiffness of surrounding anatomical structures.Approach.In our model, ligaments are classified as the most rigid structures while the softer parts of the breast are occupied by fat and glandular tissues As a result of these variations in tissue elasticity, the rapidly proliferating tumor cells are met with differential resistance. It is found that these cells are likely to circumvent stiffer terrains such as ligaments, instead electing to proliferate preferentially within the more yielding confines of the breast’s soft topography. By manipulating the interstitial tumor pressure in direct proportion to the elastic constants of the tissues surrounding the tumor, this model thus creates the potential for realizing a database of unique lesion morphology sculpted by the distinctive topography of each local anatomical infrastructure. We modeled the growth of simulated lesions within volumes extracted from fatty breast models, developed by Graffet alwith a resolution of 50μm generated with the open-source and readily available Virtual Imaging Clinical Trials for Regulatory Evaluation (VICTRE) imaging pipeline. To visualize and validate the realism of the lesion models, we leveraged the imaging component of the VICTRE pipeline, which replicates the siemens mammomat inspiration mammography system in a digital format. This system was instrumental in generating digital mammogram (DM) images for each breast model containing the simulated lesions.Results.By utilizing the DM images, we were able to effectively illustrate the imaging characteristics of the lesions as they integrated with the anatomical backgrounds. Our research also involved a reader study that compared 25 simulated DM regions of interest (ROIs) with inserted lesions from our models with DM ROIs from the DDSM dataset containing real manifestations of breast cancer. In general the simulation time for the lesions was approximately 2.5 hours, but it varied depending on the lesion’s local environment.Significance.The lesion growth model will facilitate and enhance longitudinal in silico trials investigating the progression of breast cancer.

Engineering↗