Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “source code”

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,567 records · Page 87
All sourcesNASA NTRSDOE OSTI
NASA NTRS1991

ASRM radiation and flowfield prediction status

Existing and proposed methods for the prediction of plume radiation are discussed in terms of their application to the NASA Advanced Solid Rocket Motor (ASRM) and Space Shuttle Main Engine (SSME) projects. Extrapolations of the Solid Rocket Motor (SRM) are discussed with respect to preliminary predictions of the primary and secondary radiation environments. The methodology for radiation and initial plume property predictions are set forth, including a new code for scattering media and independent secondary source models based on flight data. The Monte Carlo code employs a reverse-evaluation approach which traces rays back to their point of absorption in the plume. The SRM sea-level plume model is modified to account for the increased radiation in the ASRM plume due to the ASRM's propellant chemistry. The ASRM cycle-1 environment predictions are shown to identify a potential reason for the shutdown spike identified with pre-SRM staging.

Reardon, J. E.↗
DOE OSTI2014

Benchmarking kinetic calculations of resistive wall mode stability

Validating the calculations of kinetic resistive wall mode (RWM) stability is important for confidently predicting RWM stable operating regions in ITER and other high performance tokamaks for disruption avoidance. Benchmarking the calculations of the Magnetohydrodynamic Resistive Spectrum—Kinetic (MARS-K) [Y. Liu et al., Phys. Plasmas 15, 112503 (2008)], Modification to Ideal Stability by Kinetic effects (MISK) [B. Hu et al., Phys. Plasmas 12, 057301 (2005)], and Perturbed Equilibrium Nonambipolar Transport (PENT) [N. Logan et al., Phys. Plasmas 20, 122507 (2013)] codes for two Solov'ev analytical equilibria and a projected ITER equilibrium has demonstrated good agreement between the codes. The important particle frequencies, the frequency resonance energy integral in which they are used, the marginally stable eigenfunctions, perturbed Lagrangians, and fluid growth rates are all generally consistent between the codes. The most important kinetic effect at low rotation is the resonance between the mode rotation and the trapped thermal particle's precession drift, and MARS-K, MISK, and PENT show good agreement in this term. The different ways the rational surface contribution was treated historically in the codes is identified as a source of disagreement in the bounce and transit resonance terms at higher plasma rotation. Calculations from all of the codes support the present understanding that RWM stability can be increased by kinetic effects at low rotation through precession drift resonance and at high rotation by bounce and transit resonances, while intermediate rotation can remain susceptible to instability. The applicability of benchmarked kinetic stability calculations to experimental results is demonstrated by the prediction of MISK calculations of near marginal growth rates for experimental marginal stability points from the National Spherical Torus Experiment (NSTX) [M. Ono et al., Nucl. Fusion 40, 557 (2000)].

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

The MolSSI QCArchive project: An open-source platform to compute, organize, and share quantum chemistry data

The Molecular Sciences Software Institute's (MolSSI) Quantum Chemistry Archive (QCArchive) project is an umbrella name that covers both a central server hosted by MolSSI for community data and the Python-based software infrastructure that powers automated computation and storage of quantum chemistry (QC) results. The MolSSI-hosted central server provides the computational molecular sciences community a location to freely access tens of millions of QC computations for machine learning, methodology assessment, force-field fitting, and more through a Python interface. Facile, user-friendly mining of the centrally archived quantum chemical data also can be achieved through web applications found at the website. The software infrastructure can be used as a standalone platform to compute, structure, and distribute hundreds of millions of QC computations for individuals or groups of researchers at any scale. The QCArchiveInfrastructure is open-source (BSD-3C), code repositories can be found at github, and releases can be downloaded via PyPI and Conda. This article is categorized under: Electronic Structure Theory > Ab Initio Electronic Structure Methods Software > Quantum Chemistry Data Science > Computer Algorithms and Programming

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

Short-term electricity load forecasting: Application-driven evaluation of machine learning models across spatial and temporal scales

As we transition towards a decarbonized economy, the integration of variable renewable energy resources and new demands (e.g., electric vehicles, heat pumps) into the electricity grid places unprecedented pressure on grid operators to effectively anticipate and manage peak load. In this context, machine learning algorithms are proving to be indispensable for accurate short-term load forecasting, a crucial task to address these challenges. This study benchmarks 6 machine learning algorithms, including three neural networks and three tree-based algorithms, across various levels of spatial aggregation and time horizons (1, 4, 8, 24, and 48 h). The central contribution of this work is the comparison and analysis of load forecasting models not only based on statistical metrics, but also based on a novel error metric, which evaluates the cost implications of forecast errors for power system stakeholders. Results show that tree-based models outperform neural networks, based on statistical metrics, and yield less skewed error distributions for most spatial scales. However, through the lens of the novel error metric, neural networks are the more competitive choice, especially for forecast horizons that exceed 8 h. The study concludes with actionable recommendations to grid operators and highlights the need for the development of error metrics that link forecasting accuracy to operational costs. To promote transparency and open science, the datasets and Python code are open-sourced via a supplementary repository.

Houben, Nikolaus↗
DOE OSTI2023

Computer-controlled finishing via dynamically constraint position-velocity-time scheduler

In a Computer Numerical Controlled (CNC) finishing process, the target material removal from an optical surface is guided by the convolution between the influence function of a machine tool and its dwell time at certain points over the surface. To reduce dynamics stressing and increase machining efficiency, the dwell time must be converted to varying velocities, which are the actual inputs to the machine tool controller. Conventionally, the conversion assumed constant acceleration and relied on linear motion interpolation, which caused discontinuities in velocities. This unsmooth motion affects the material removal distribution, and, thus, the accuracy of the finished surface shape. Many modern CNC machines support the smoother, cubic-polynomial interpolated Position-Velocity-Time (PVT) motion mode; however, the conventional scheduler may fail to provide suitable velocities for the PVT. Here in this study answers this challenge by proposing a novel PVT-based velocity scheduler that achieves smooth motion while considering CNC dynamic limits. Firstly, the principle of the PVT is explained, and the PVT-based velocity scheduler is formulated. Secondly, a quadratic programming is used to optimize the velocities by imposing the CNC dynamic constraints and the C 1 continuities (zeroth and first derivatives are continuous) simultaneously. Thirdly, the smoothness and accuracy of the scheduled velocities are studied on different kinds of tool paths via simulation. Finally, a sub-0.3 nm level surface finishing experiment using ion beam figuring is demonstrated to verify the feasibility of the proposed method. The PVT-based scheduler and simulator code is open-sourced.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Active learning for accelerated design of layered materials

Hetero-structures made from vertically stacked monolayers of transition metal dichalcogenides hold great potential for optoelectronic and thermoelectric devices. Discovery of the optimal layered material for specific applications necessitates the estimation of key material properties, such as electronic band structure and thermal transport coefficients. However, screening of material properties via brute force ab initio calculations of the entire material structure space exceeds the limits of current computing resources. Moreover, the functional dependence of material properties on the structures is often complicated, making simplistic statistical procedures for prediction difficult to employ without large amounts of data collection. Here, we present a Gaussian process regression model, which predicts material properties of an input hetero-structure, as well as an active learning model based on Bayesian optimization, which can efficiently discover the optimal hetero-structure using a minimal number of ab initio calculations. The electronic band gap, conduction/valence band dispersions, and thermoelectric performance are used as representative material properties for prediction and optimization. The Materials Project platform is used for electronic structure computation, while the BoltzTraP code is used to compute thermoelectric properties. Bayesian optimization is shown to significantly reduce the computational cost of discovering the optimal structure when compared with finding an optimal structure by building a regression model to predict material properties. The models can be used for predictions with respect to any material property and our software, including data preparation code based on the Python Materials Genomics (PyMatGen) library as well as python-based machine learning code, is available open source.

36 MATERIALS SCIENCE↗
DOE OSTI2013

Edible oil structures at low and intermediate concentrations. II. Ultra-small angle X-ray scattering of in situ tristearin solids in triolein

Ultra-small angle X-ray scattering has been used for the first time to elucidate, in situ, the aggregation structure of a model edible oil system. The three-dimensional nano- to micro-structure of tristearin solid particles in triolein solvent was investigated using 5, 10, 15, and 20% solids. Three different sample preparation procedures were investigated: two slow cooling rates of 0.5°/min, case 1 (22 days of storage at room temperature) and case 2 (no storage), and one fast cooling of 30°/min, case 3 (no storage). The length scale investigated, by using the Bonse-Hart camera at beamline ID-15D at the Advanced Photon Source, Argonne National Laboratory, covered the range from 300 Å to 10 μm. The unified fit and the Guinier-Porod models in the Irena software were used to fit the data. The former was used to fit 3 structural levels. Level 1 structures showed that the primary scatterers were essentially 2-dimensional objects for the three cases. The scatterers possessed lateral dimensions between 1000 and 4300 Å. This is consistent with the sizes of crystalline nanoplatelets present which were observed using cryo-TEM. Level 2 structures were aggregates possessing radii of gyration, R{sub g2} between 1800 Å and 12000 Å and fractal dimensions of either D{sub 2}=1 for case 3 or 1.8≤D{sub 2}≤2.1 for case 1 and case 2. D{sub 2} = 1 is consistent with unaggregated 1-dimensional objects. 1.8 ≤ D{sub 2} ≤ 2.1 is consistent with these 1-dimensional objects (below) forming structures characteristic of diffusion or reaction limited cluster-cluster aggregation. Level 3 structures showed that the spatial distribution of the level 2 structures was uniform, on the average, for case 1, with fractal dimension D{sub 3}≈3 while for case 2 and case 3 the fractal dimension was D{sub 3}≈2.2, which suggested that the large-scale distribution had not come to equilibrium. The Guinier-Porod model showed that the structures giving rise to the aggregates with a fractal dimension given by D{sub 2} in the unified fit level 2 model were cylinders described by the parameter s≈1 in the Guinier-Porod model. The size of the base of these cylinders was in agreement with the cryo-TEM observations as well as with the results of the level 1 unified fit model. By estimating the size of the nanoplatelets and understanding the structures formed via their aggregation, it will be possible to engineer novel lipids systems that embody desired functional characteristics.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2015

The CARIBU EBIS control and synchronization system

The Californium Rare Isotope Breeder Upgrade (CARIBU) Electron Beam Ion Source (EBIS) charge breeder has been built and tested. The bases of the CARIBU EBIS electrical system are four voltage platforms on which both DC and pulsed high voltage outputs are controlled. The high voltage output pulses are created with either a combination of a function generator and a high voltage amplifier, or two high voltage DC power supplies and a high voltage solid state switch. Proper synchronization of the pulsed voltages, fundamental to optimizing the charge breeding performance, is achieved with triggering from a digital delay pulse generator. The control system is based on National Instruments realtime controllers and LabVIEW software implementing Functional Global Variables (FGV) to store and access instrument parameters. Fiber optic converters enable network communication and triggering across the platforms.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Synthesis of polycaprolactone-titanium oxide multilayer films by nanosecond laser pulses and electrospinning technique for better implant fabrication

In this study, a combination of electrospinning and laser texturing is introduced as a novel method for increasing the biocompatibility of metal implants. Besides having a rough laser treated surface, the implant benefits from the high porosity and better wettability of an electrospun fibrous structure, which is a more favorable environment for cell proliferation. Titanium samples were patterned using a nanosecond laser beam and were placed as collectors in an electrospinning machine. They were then soaked in simulated body fluid for four weeks. Energy Dispersive X-ray and X-Ray Diffraction results indicate significantly more hydroxyapatite formation on laser treated samples with nanoscale fibers deposited on their surface. This shows that having a laser treated surface underneath the fibrous layer can improve short-term biocompatibility even before degradation of fibers. The thermal conductivity of the electrospun layer, measured using a Hot Disk Transient Plane Source instrument and computer code, was shown to be considerably lower than that of titanium and very close to bone. The presence of this layer can therefore be beneficial in making the implant more compatible to a biological medium. In case of dental implants, it was shown that this layer can act as a thermal barrier while a hot beverage is consumed and it can decrease the temperature rise by about 60%, which avoids any possible damage to newly formed cells during the healing period.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Pushing configuration-interaction to the limit: Towards massively parallel MCSCF calculations

A new large-scale parallel multiconfigurational self-consistent field (MCSCF) implementation in the open-source NWChem computational chemistry code is presented. The generalized active space approach is used to partition large configuration interaction (CI) vectors and generate a sufficient number of batches that can be distributed to the available cores. Massively parallel CI calculations with large active spaces can be performed. The new parallel MCSCF implementation is tested for the chromium trimer and for an active space of 20 electrons in 20 orbitals, which can now routinely be performed. Unprecedented CI calculations with an active space of 22 electrons in 22 orbitals for the pentacene systems were performed and a single CI iteration calculation with an active space of 24 electrons in 24 orbitals for the chromium tetramer was possible. In conclusion, the chromium tetramer corresponds to a CI expansion of one trillion Slater determinants (914 058 513 424) and is the largest conventional CI calculation attempted up to date.

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

Implications of Voyager 1 observations beyond the heliopause for the local interstellar electron spectrum

Cosmic-ray observations made by the Voyager 1 spacecraft outside the dominant modulating influence of the heliosphere finally allow the comparison of computed galactic spectra with experimental data at lower energies. These computed spectra, based on galactic propagation models, can now be compared with observations at low energies by Voyager 1 and at high energies by the PAMELA space detector at Earth. This improves understanding of basic propagation effects and also provides solar modulation studies with reliable input spectra from 1 MeV to 100 GeV. We set out to reproduce the Voyager 1 electron observations in the energy range of 6-60 MeV, as well as the PAMELA electron spectrum above 10 GeV, using the GALPROP code. By varying the source spectrum and galactic diffusion parameters, specifically the rigidity dependence of spatial diffusion, we find local interstellar spectra that agree with both power-law spectra observed by Voyager 1 beyond the heliopause. The local interstellar spectrum between ∼1 MeV and 100 GeV indicates that it is the combination of two power laws, with E {sup –(1.45} {sup ±} {sup 0.15)} below ∼100 MeV and E {sup –(3.15} {sup ±} {sup 0.05)} above ∼100 MeV. A gradual turn in the spectral shape matching the power laws is found, between 2.0 ± 0.5) GeV and (100 ± 10) MeV. According to our simplified modeling, this transition is caused primarily by galactic propagation effects. We find that the intensity beyond the heliopause at 10 MeV is (350 ± 50) electrons m{sup –2} s{sup –1} sr{sup –1} MeV{sup –1}, decreasing to (50 ± 5) electrons m{sup –2} s{sup –1} sr{sup –1} MeV{sup –1} at 100 MeV.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Development of an open-source segmented blade design tool

As wind turbines continue to grow ever larger to reduce the cost of energy, their blades follow suit, with the largest commercial offshore blades extending past 100 m. Massive blades such as these raise key transportation and manufacturing challenges, especially for land-based turbines. Segmented blades are one solution and are garnering increased industry and research interest. In this work, a detailed mechanical joint model is integrated into the Wind-Plant Integrated System Design and Engineering Model (WISDEM ® ), which will facilitate future segmented blade research and optimization. WISDEM is used to design a wind turbine with 100-m segmented blades. This wind turbine design is compared to other machines with 100-m monolithic blades designed for rail-transportability. The designs are compared in terms of blade mass and cost, turbine capital cost, annual energy production, and levelized cost of energy, with monolithic designs being the lightest and most economical. However, this result may vary by wind plant location. A variety of segmentation joint types exist, and they will inevitably vary in parameters such as cost, spanwise location, and physical characteristics. This work examines the sensitivity of wind turbine design drivers and annual energy production to a variety of the aforementioned parameters, using the open-source wind turbine design codes OpenFAST and WISDEM, finding that joint mass, stiffness, and location can have significant effects on design drivers.

17 WIND ENERGY↗
DOE OSTI2020

ChemEnv : a fast and robust coordination environment identification tool

Coordination or local environments have been used to describe, analyze and understand crystal structures for more than a century. Here, a new tool called ChemEnv, which can identify coordination environments in a fast and robust manner, is presented. In contrast to previous tools, the assessment of the coordination environments is not biased by small distortions of the crystal structure. Its robust and fast implementation enables the analysis of large databases of structures. The code is available open source within the pymatgen package and the software can also be used through a web app available on http://crystaltoolkit.org through the Materials Project.

42 ENGINEERING↗
DOE OSTI2014

SU-E-I-66: Radiomics and Image Registration Updates for the Computational Environment for Radiotherapy Research (CERR)

Purpose: To present new tools in CERR for Radiomics, image registration and other software updates and additions. Methods: Radiomics: CERR supports generating 3-D texture metrics based on gray scale co-occurance. Two new ways to calculate texture features were added: (1) Local Texture Averaging: Local texture is calculated around a voxel within the userdefined bounding box. The final texture metrics are the average of local textures for all the voxels. This is useful to detect any local texture patterns within an image. (2) Image Smoothing: A convolution ball of user-defined radius is rolled over an image to smooth out artifacts. The texture metrics are then computed on the smooth image. Image Registration: (1) Support was added to import deformation vector fields as well as non-deformable transformation matrices generated by vendor software and stored in standard DICOM format. (2) Support was added to use image within masks while computing image deformations. CT to MR registration is supported. This registration uses morphological edge information within the images to guide the deformation process. In addition to these features, other noteworthy additions to CERR include (1) Irregularly shaped ROI: This is done by taking intersection between infinitely extended irregular polygons drawn on any of the two views. Such an ROI is more conformal and useful in avoiding any unwanted parts of images that cannot be avoided with the conventional cubic box. The ROI is useful to generate Radiomics metrics. (2) Ability to insert RTDOSE in DICOM format to existing CERR plans. (3) Ability to import multi-frame PET-CT and SPECT-CT while maintaining spatial registration between the two modalities. (4) Ability to compile CERR on Unix-like systems. Results: The new features and updates are available via https://www.github.com/adityaapte/cerr . Conclusion: Features added to CERR increase its utility in Radiomics, Image-Registration and Outcomes modeling.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2014

SU-E-T-212: Influence of the Modulation Index On Daily Quality Assurance in Rapid Arc Treatments

Purpose: At our Institute some measured parameters for daily quality assurance (DQA) of dynamic arc therapy plans showed an unexpected behavior, therefore an investigation of the influence of the magnitude of modulation was conducted. Methods: In our clinical practice all DQAs of dynamic arc therapy plans are measured and analyzed prior to treatments using commercial software. For this study these plans were additionally exported to our in-house software written in MATLAB.The developed software extracted the leaf position, gantry angle, cumulative meterset weight of each control point (CP) and the total number of Monitor Units (MU) of each arc. Based on this information we calculated the leaf travel distance, irradiated segment area, number of MUs and dose rate for each CP. These data allowed us to calculate the modulation indexes (MI) of the plans, applying five different definitions of MI. The results were then correlated to the results of our DQA measurements.To validate the software, additional plans of known MIs were created and analyzed. For confirmation, the calculated parameters were compared to the segmented treatment table (STT) coming from the treatment planning system. Results: All calculated CP-parameters matched the STT by 99% or better. For linac 1, the comparison of the MI evaluation and the DQA results showed a slight tendency: 91.3% failed DQA plans had a MI lower than the average value. For this definition we consider that the lower the MI the higher the modulation. The results of the linac 2 present no significant relevance due to the low sample sizes for each DQA software. Conclusion: Available data and given definitions of the modulation index do not bring conclusive results; one cannot find a clear and distinct correlation with the failure of the DQA. The ongoing analysis with an increased sample size might lead to another conclusion.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2014

SU-E-T-562: Scanned Percent Depth Dose Curve Discrepancy for Photon Beams with Physical Wedge in Place (Varian IX) Using Different Sensitive Volume Ion Chambers

Purpose: To investigate and report the discrepancy of scanned percent depth dose (PDD) for photon beams with physical wedge in place when using ion chambers with different sensitive volumes. Methods/Materials: PDD curves of open fields and physical wedged fields (15, 30, 45, and 60 degree wedge) were scanned for photon beams (6MV and 10MV, Varian iX) with field size of 5x5 and 10x10 cm using three common scanning chambers with different sensitive volumes - PTW30013 (0.6cm3), PTW23323 (0.1cm3) and Exradin A16 (0.007cm3). The scanning system software used was OmniPro version 6.2, and the scanning water tank was the Scanditronix Wellhoffer RFA 300.The PDD curves from the three chambers were compared. Results: Scanned PDD curves of the same energy beams for open fields were almost identical between three chambers, but the wedged fields showed non-trivial differences. The largest differences were observed between chamber PTW30013 and Exradin A16. The differences increased as physical wedge angle increased. The differences also increased with depth, and were more pronounced for 6MV beam. Similar patterns were shown for both 5x5 and 10x10 cm field sizes. For open fields, all PDD values agreed with each other within 1% at 10cm depth and within 1.62% at 20 cm depth. For wedged fields, the difference of PDD values between PTW30013 and A16 reached 4.09% at 10cm depth, and 5.97% at 20 cm depth for 6MV with 60 degree physical wedge. Conclusion: We observed a significant difference in scanned PDD curves of photon beams with physical wedge in place obtained when using different sensitive volume ion chambers. The PDD curves scanned with the smallest sensitive volume ion chamber showed significant difference from larger chamber results, beyond 10cm depth. We believe this to be caused by varying response to beam hardening by the wedges.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2016

TH-AB-207A-07: Radiation Dose Simulation for a Newly Proposed Dynamic Bowtie Filters for CT Using Fast Monte Carlo Methods

Purpose: Dynamic bowtie filter is an innovative design capable of modulating the X-ray and balancing the flux in the detectors, and it introduces a new way of patient-specific CT scan optimizations. This study demonstrates the feasibility of performing fast Monte Carlo dose calculation for a type of dynamic bowtie filter for cone-beam CT (Liu et al. 2014 9(7) PloS one) using MIC coprocessors. Methods: The dynamic bowtie filter in question consists of a highly attenuating bowtie component (HB) and a weakly attenuating bowtie (WB). The HB is filled with CeCl3 solution and its surface is defined by a transcendental equation. The WB is an elliptical cylinder filled with air and immersed in the HB. As the scanner rotates, the orientation of WB remains the same with the static patient. In our Monte Carlo simulation, the HB was approximated by 576 boxes. The phantom was a voxelized elliptical cylinder composed of PMMA and surrounded by air (44cm×44cm×40cm, 1000×1000×1 voxels). The dose to the PMMA phantom was tallied with 0.15% statistical uncertainty under 100 kVp source. Two Monte Carlo codes ARCHER and MCNP-6.1 were compared. Both used double-precision. Compiler flags that may trade accuracy for speed were avoided. Results: The wall time of the simulation was 25.4 seconds by ARCHER on a 5110P MIC, 40 seconds on a X5650 CPU, and 523 seconds by the multithreaded MCNP on the same CPU. The high performance of ARCHER is attributed to the parameterized geometry and vectorization of the program hotspots. Conclusion: The dynamic bowtie filter modeled in this study is able to effectively reduce the dynamic range of the detected signals for the photon-counting detectors. With appropriate software optimization methods, the accelerator-based (MIC and GPU) Monte Carlo dose engines have shown good performance and can contribute to patient-specific CT scan optimizations.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2017

Pececillo

Pececillo is a proxy-app for the open source Truchas metal processing code (LA-CC-15-097). It implements many of the physics models used in Truchas: free-surface, incompressible Navier-Stokes fluid dynamics (e.g., water waves); heat transport, material phase change, view factor thermal radiation; species advection-diffusion; quasi-static, elastic/plastic solid mechanics with contact; electomagnetics (Maxwell's equations). The models are simplified versions that retain the fundamental computational complexity of the Truchas models while omitting many non-essential features and modeling capabilities. The purpose is to expose Truchas algorithms in a greatly simplified context where computer science problems related to parallel performance on advanced architectures can be more easily investigated. While Pececillo is capable of performing simulations representative of typical Truchas metal casting, welding, and additive manufacturing simulations, it lacks many of the modeling capabilites needed for real applications.

Jibben, Zechariah↗