Search NASA⌕ Search

SEARCH · Search NASA

Results for “FREE BOUNDARY”

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 OSTI2024

Observed Variability in Convective Cell Characteristics and Near-Storm Environments across the Sea- and Bay-Breeze Fronts in Southeast Texas

Abstract During the DOE Atmospheric Radiation Measurement (ARM) Tracking Aerosol Convection Interactions Experiment (TRACER) IOP spanning June–September 2022, two fixed ARM sites and a mobile team concurrently sampled the airmass heterogeneity across sea- and bay-breeze fronts around the greater Houston metropolitan region. Here, we quantify the spatiotemporal variability between maritime (coastal/bay side of breeze fronts) and continental (inland side of breeze fronts) air masses over 15 IOP days characterized by strong sea-breeze forcing. We analyze environmental profile data from 177 radiosondes and use S- and C-band radar data to track and quantify the variability in attributes of more than 2300 shallow and transitioning cells across different air masses. The composite analysis of environmental profiles indicates that during the early afternoon, the sea-breeze maritime air mass exhibits lower convective available potential energy (CAPE) than the bay-breeze maritime air mass. As the sea breeze advances inland with time, CAPE within the maritime air mass exceeds that of the continental air mass to the north of the breeze fronts. In general, maritime cells have a larger mean composite reflectivity and cell widths than continental cells; however, the response varies between shallow and transitioning cells. Mean composite 20-dB Z echo-top heights, however, are similar across air masses for both shallow and transitioning cells. The continental and maritime inflow air mass for transitioning cells has significantly different mean values for mixed-layer entrainment CAPE, lifted condensation level, level of free condensation, boundary layer depth, and diluted equilibrium level. For shallow cells, only total precipitable water shows a significant difference. Significance Statement The greater Houston metropolitan area is a natural laboratory for understanding the individual impacts of background meteorology and aerosols on convective clouds. Due to its proximity to the Gulf Coast and Galveston Bay, the Houston region experiences a diurnal precipitation cycle in the summer, driven by convection triggered from sea- and bay-breeze fronts. These fronts act as a boundary between air masses with distinct thermodynamic and environmental characteristics. Convergence along these fronts and interactions between storm outflow and the fronts facilitate convection initiation in different mesoscale air masses. This study quantifies the heterogeneity among these air masses while investigating their influence on cloud microphysics. We find that the effect of airmass heterogeneity is more pronounced for the bulk microphysical properties in shallow clouds.

Sharma, Milind↗
DOE OSTI2020

Spectrum of embrittling potencies and relation to properties of symmetric-tilt grain boundaries

Grain boundaries in metallic materials can exist in a wide range of stable and metastable structures. In addition, the properties of a grain boundary may be altered through solute segregation. In this work, we present a formulation that combines the spectrum of embrittling potencies associated with solute segregation with site-occupancy statistics. As a prototype problem, we illustrate the relation between segregation and embrittlement in the case of S segregation to grain boundaries in Ni. To obtain a population of site segregation energies, we perform molecular statics calculations on 378 different symmetric-tilt grain boundaries and their free surface equivalents, using an embedded-atom method interatomic potential developed specifically for studying embrittlement. Our results show that it is important to consider both the energies associated with embrittlement and the probability of occupancy to describe the general embrittling nature of a grain boundary. When analyzed in isolation, certain grain boundaries show large embrittling potencies; however, that effect is diminished when the probability of S segregation to that grain boundary is considered within a polycrystal. We propose a new quantity, the embrittling estimator, which not only categorizes grain boundaries as embrittling or strengthening, but also considers site occupancy probabilities, so that the embrittlement behavior of grain boundaries within a network of grain boundaries can be compared. Finally, we examine the relationship between embrittlement behavior and innate grain boundary properties, such as the free volume, and find statistical evidence that the complex nature of embrittlement cannot be explained by linear correlations with excess volumes or energies. Ultimately, this combined approach provides a theoretical tool to assist grain boundary engineering of metastable alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Radiation Tolerance of Nanoporous Gadolinium Titanate

Defect sinks play a crucial role in reducing radiation-induced damage accumulation, but their effectiveness varies. This study directly compares the efficacy of grain boundaries and pores (free surfaces) in gadolinium titanate (Gd 2 Ti 2 O 7 ) using in situ ion irradiation transmission electron microscopy (TEM). A unique sample configuration is developed, where each of the three distinct regions, each containing only one primary sink of interest, is in the same TEM lamella for simultaneous irradiation. Regular assessments of crystallinity in each region via nanobeam electron diffraction show that, at both room temperature and 600°C, free surface (pore) defect sinks are more effective than grain boundaries in capturing radiation-induced defects and delaying the onset of amorphization. These findings have significant implications for further microstructural engineering and the design of radiation-resistant microstructures.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Weak boundary conditions for Lagrangian shock hydrodynamics: A high-order finite element implementation on curved boundaries

Here, we propose a new Nitsche-type approach for weak enforcement of normal velocity boundary conditions for a Lagrangian discretization of the compressible shock-hydrodynamics equations using high-order finite elements on curved boundaries. Specifically, the variational formulation is appropriately modified to enforce free-slip wall boundary conditions, without perturbing the structure of the function spaces used to represent the solution, with a considerable simplification with respect to traditional approaches. Total energy is conserved and the resulting mass matrices are constant in time. The robustness and accuracy of the proposed method are validated with an extensive set of tests involving nontrivial curved boundaries.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Origin and regulation of oxygen redox instability in high-voltage battery cathodes

Oxygen redox at high-voltage has emerged as a transformative paradigm for high-energy battery cathodes such as layered transition-metal oxides by offering extra capacity beyond conventional transition-metal redox. However, these cathodes suffer from voltage hysteresis, voltage fade, and capacity drop upon cycling; single-crystalline cathodes have recently shown some improvements but these challenges still remain. We reveal the fundamental origin of oxygen redox instability to be originated from the domain boundaries that are present in single-crystalline cathode particles. By synthesizing single-crystalline cathodes free of domain boundaries, we show that the elimination of domain boundaries enhances the reversible lattice oxygen redox while inhibiting the irreversible oxygen release. Further, this leads to significantly suppressed structural degradation and improved mechanical integrity during battery cycling and abuse heating. The robust oxygen redox enabled through domain boundary control provides practical opportunities towards high-energy, long-cycling, and safe batteries.

25 ENERGY STORAGE↗
DOE OSTI2023

Separation of K + and Bi 3+ displacements in a Pb-free, monoclinic piezoelectric at the morphotropic phase boundary

The best piezoelectric properties of any perovskite oxide known are found in the solid solution of the relaxor Pb (Mg 1/3 Nb 2/3 )O 3 and ferroelectric PbTiO 3 . Despite its impressive properties, this system has limited analogy. We present the compositional exploration of the Pb-free analogue (1-x)(K 1/2 Bi 1/2 )(Mg 1/3 Nb 2/3 )O 3 -x(K 1/2 Bi 1/2 )TiO 3 (KBMN-KBT). We locate the morphotropic phase boundary between x = 0.86 and 0.88 changing from Cm to Pm symmetry and the optimally performing composition at x = 0.88. We report a piezoelectric figure of merit (d 33 *) of 192 pm V —1 from strain measurements. Diffraction methods reveal disordered displacements of K + and Bi 3+ which persist from the KBMN endmember through multiple changes in symmetry. Rearrangement of the Bi 3+ displacements along the uncommon [011] c direction drives the physical response. Ferroelectric, dielectric, and piezoresponse force microscopy are used to study the progression of physical properties through the MPB and attribute the mechanism to a polarization rotation. Taking account for local, short-range, and average structural features yield a balanced perspective on the structure and properties of this system, isolating the driving force within this system to the Bi 3+ bonding configuration. This work yields a strong analogy to the Pb-based analogue, and provides strategies for further optimization.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Effect of different grain boundary diffusion alloys on magnetic properties of Dy-free sintered NdFeB magnet

Dy-free sintered magnets were fabricated by blending Neo powder with different grain boundary diffusion alloy powders. Cu, CeAl and CeAlCu have a negative effect on H cj of magnets, while PrAlCu and AlCuGa have a positive effect. The PrAlCu-added magnet achieves the best magnetic properties among all magnets with the additions of different diffusion alloys. With increasing PrAlCu from 0-10 wt.%, H cj of the magnets gradually increases from the original 14.5 kOe to 19.2 kOe. The magnet with 7.5 wt.% PrAlCu obtains a H cj of 18 kOe and (BH) max of 39.1 MGOe. It is found that Pr and Cu in PrAlCu alloy is mainly distributed at grain boundary and triple junctions, leading to a reduced coupling among grains, thus an enhanced H cj . Here, the grain boundary engineering by adding an appropriate alloy is an effective method to improve H cj of Dy-free NdFeB magnets.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Impact of Lateral Groundwater Flow and Subsurface Lower Boundary Conditions on Atmospheric Boundary Layer Development over Complex Terrain

Credible soil moisture redistribution schemes are essential to meteorological models, as lower boundary moisture influences the balance of surface turbulent fluxes and atmospheric boundary layer (ABL) development. While land surface models (LSMs) have vastly improved in their hydrologic representation, several commonly held assumptions, such as free-draining lower boundary, one-dimensional moisture flux, and lack of groundwater representation, can bias the terrestrial water balance. This study explores the impact of LSM hydrology representation on ABL development in the Weather Research and Forecasting (WRF) meteorological model. The results of summertime WRF simulations with Noah LSM, characterized by 2-m-thick soil and one-dimensional flow, are shown for a domain in the Colorado Rocky Mountain headwaters region. A reference WRF simulation is compared to 1) the same model with soil moisture initialized by the hydrologic model ParFlow; 2) a deep, free-draining simulation; and 3) WRF coupled to ParFlow, a three-dimensional, integrated groundwater-surface water model. Results show that both lateral transport of groundwater and the rate of drainage from the lower soil layer can weaken or reverse the coupling strength between evaporative fraction and ABL over a 5-month summer period. The resulting shifts in low-level moist convection in river valleys and thermally driven airflows yield strengthened anabatic upslope winds and perturbations to regional precipitation.

Forrester, Mary M.↗
DOE OSTI2023

Development of Dy-free Nd-Fe-B-based sintered magnet through grain boundary engineering using Pr-Cu alloys

For this study, Dy-free sintered magnets were fabricated by mixing Nd-Fe-B powder with various amounts of Pr-Cu powder followed by the common procedures for making Nd-Fe-B based sintered magnets. With the Pr-Cu addition increasing from 0 to 10 wt.%, the obtained magnets’ H cj increases from the original 14.5 kOe to 18.6 kOe. The highest (BH) max achieved was 35.0 MGOe with 7.5 wt.% PrCu. The distribution of the Pr and Cu elements was primarily at grain boundary and triple junctions, leading to a reduced coupling among grains, thus an enhanced H cj .

36 MATERIALS SCIENCE↗
DOE OSTI2021

In Situ TEM Tensile Testing of Bicrystals with Tailored Misorientation Angles

Grain boundaries have complex structural features that influence strength, ductility and fracture in metals and alloys. Grain boundary misorientation angle has been identified as a key parameter that controls their mechanical behavior, but the effect of misorientation angle has been challenging to isolate in polycrystalline materials. Here, we describe the use of bicrystal Au thin films made using a rapid melt growth process to study deformation at a single grain boundary. Tensile testing is performed on bicrystals with different misorientation angles using in situ TEM, as well as on a single crystalline sample. Plastic deformation is initiated through dislocation nucleation from free surfaces. Grain boundary sliding is not observed, and failure occurs away from the grain boundary through plastic collapse in all cases. The failure behavior in these nanoscale bicrystals does not appear to depend on the misorientation angle or grain boundary energy but instead has a more complex dependence on sample surface structure and dislocation activity.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Computing the shape gradient of stellarator coil complexity with respect to the plasma boundary

Coil complexity is a critical consideration in stellarator design. The traditional two-step optimization approach, in which the plasma boundary is optimized for physics properties and the coils are subsequently optimized to be consistent with this boundary, can result in plasma shapes which cannot be produced with sufficiently simple coils. To address this challenge, we propose a method to incorporate considerations of coil complexity in the optimization of the plasma boundary. Coil complexity metrics are computed from the current potential solution obtained with the REGCOIL code (Landreman, Nucl. Fusion , vol. 57, 2017, 046003). While such metrics have previously been included in derivative-free fixed-boundary optimization (Drevlak et al. , Nucl. Fusion , vol. 59, 2018, 016010), we compute the local sensitivity of these metrics with respect to perturbations of the plasma boundary using the shape gradient (Landreman & Paul, Nucl. Fusion , vol. 58, 2018, 076023). We extend REGCOIL to compute derivatives of these metrics with respect to parameters describing the plasma boundary. In keeping with previous research on winding surface optimization (Paul et al. , Nucl. Fusion , vol. 58, 2018, 076015), the shape derivatives are computed with a discrete adjoint method. In contrast with the previous work, derivatives are computed with respect to the plasma surface parameters rather than the winding surface parameters. To further reduce the resolution required to compute the shape gradient, we present a more efficient representation of the plasma surface which uses a single Fourier series to describe the radial distance from a coordinate axis and a spectrally condensed poloidal angle. This representation is advantageous over the standard cylindrical representation used in the VMEC code (Hirshman & Whitson, Phys. Fluids , vol. 26, 1983, pp. 3553–3568), as it provides a uniquely defined poloidal angle, eliminating a null space in the optimization of the plasma surface. In comparison with previous spectral condensation methods (Hirshman & Breslau, Phys. Plasmas , vol. 5, 1998, p. 2664), the modified poloidal angle is obtained algebraically rather than through the solution of a nonlinear optimization problem. The resulting shape gradient highlights features of the plasma boundary that are consistent with simple coils and can be used to couple coil and fixed-boundary optimization.

Physics↗
DOE OSTI2023

Non‐Linear Kinetics of The Lithium Metal Anode on Li 6 PS 5 Cl at High Current Density: Dendrite Growth and the Role of Lithium Microstructure on Creep

Abstract Interfacial instability, viz., pore formation in the lithium metal anode (LMA) during discharge leading to high impedance, current focusing induced solid–electrolyte (SE) fracture during charging, and formation/behaviour of the solid–electrolyte interphase (SEI), at the anode, is one of the major hurdles in the development of solid‐state batteries (SSBs). Also, understanding cell polarization behaviour at high current density is critical to achieving the goal of fast‐charging battery and electric vehicle. Herein, via in situ electrochemical scanning electron microscopy (SEM) measurements, performed with freshly deposited lithium microelectrodes on transgranularly fractured fresh Li6PS5Cl (LPSCl), the LiǀLPSCl interface kinetics are investigated beyond the linear regime. Even at relatively small overvoltages of a few mV, the LiǀLPSCl interface shows non‐linear kinetics. The interface kinetics possibly involve multiple rate‐limiting processes, i.e., ion transport across the SEI and SE|SEI interfaces, as well as charge transfer across the LiǀSEI interface. The total polarization resistance R P of the microelectrode interface is determined to be ≈ 0.8 Ω cm 2 . It is further shown that the nanocrystalline lithium microstructure can lead to a stable LiǀSE interface via Coble creep along with uniform stripping. Also, spatially resolved lithium deposition, i.e., at grain surface flaws, grain boundaries, and flaw‐free surfaces, indicates exceptionally high mechanical endurance of flaw‐free surfaces toward cathodic load (>150 mA cm −2 ). This highlights the prominent role of surface defects in dendrite growth.

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

Microstructure and strengthening of Al-6Ce-3Ni-0.7Fe (wt%) alloy manufactured by laser powder-bed fusion

In this study, an Al-6Ce-3Ni-0.7Fe (wt%) alloy was fabricated via laser powder-bed fusion and its microstructure, thermal stability, tensile properties, and creep properties are investigated for two rapid-solidification rates with different eutectic spacings. For both faster- and slower-cooled states, the as-fabricated alloy mostly shows elongated grains with very fine eutectic networks (~40 nm lamellar width) comprising a high-volume fraction (~18 %) of intermetallic phases. Faster-cooled samples, however, display a less continuous eutectic network with finer spacing (~100 nm), leading to higher Orowan strengthening and therefore superior microhardness and yield stress up to ~350 °C. Upon aging (300 - 450 °C), micron-size Al 9 (Ni,Fe) 2 needle-like precipitates form within grains, and the Al 11 Ce 3 eutectic network spheroidizes, while retaining a submicron width and spacing, resulting in 33–37 % drop in microhardness after 144 h aging at 400 °C. Limited tensile ductility (~6 %) and creep ductility (~1-2 %) are measured for both faster- and slower-cooled alloys due to an inhomogeneous microstructure, where cavitation preferentially initiates at melt-pool boundaries, precipitate-free zones, and/or denuded zones, eventually leading to local fracture. Denuded zones, induced by stress, form due to diffusional flow with stress-dependent orientations, and are identified as microstructurally weak regions leading to strain localization. The present alloy does not show a significant difference in 300 °C creep resistance between: (i) tensile and compressive loading, (ii) faster- and slower-cooled samples, and (iii) continuous and spheroidized eutectics. Creep resistance at 300 °C is comparable to that of a Ce-richer Al-10.5Ce-3.1Ni-1.2Mn (wt%) alloy, despite the formation of denuded zones and needle-like Al 9 (Ni,Fe) 2 precipitates.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Towards defect engineering in hexagonal MoS 2 nanosheets for tuning hydrogen evolution and nitrogen reduction reactions

Combined computational and experimental approaches were used to evaluate defective 2H-MoS 2 nanosheets for their activity and selectivity for hydrogen evolution reaction (HER) and nitrogen reduction reaction (NRR). Density functional theory calculations were used to understand the relationship between HER and NRR activity on the ideal basal MoS 2 plane, seven grain-boundaries, ten single-/few-atom vacancies and anti-sites, and zigzag and armchair edge sites. The results confirm that 2H-MoS 2 should contain several defects with high activity for HER: armchair and zigzag edges, VS vacancy, MoS 2 anti-site, and S-S and Mo-Mo grain boundaries. Considering Gibbs free energy change for all the steps in the NRR mechanism and kinetic barriers for a key NRR step, we have found that activation of perspective NRR selective sites in 2H-MoS 2 , namely V MoS6 and clusters of S-vacancies, would require large overpotential, conditions at which HER dominates. The DFT conclusions are supported by the electrochemical studies of NRR activity and selectivity under aqueous conditions, which show an increase in NRR activity but a decrease in Faradaic efficiency as applied cell potential becomes more negative. The results of this work therefore highlight the challenges in activating natural 2H-MoS 2 for NRR, which would require additional material engineering or reaction condition optimization as a way to suppress HER, decrease the NRR overpotential or preferentially both.

2H-MoS2↗
DOE OSTI2024

Identifying multiple synergistic factors on the susceptibility to stress relaxation cracking in variously heat-treated weldments

The 347H austenitic stainless steel has been widely used for pressure vessels and pipeline (PVP) applications due to its excellent creep and corrosion resistance, which fit ideally to the harsh conditions in petrochemical industries, fossil fuel or nuclear power plants, and modern energy storages. However, a failure mode has been commonly observed with cracks emerging at the heat affected zone (HAZ) of weldments during post-weld heat treatment (PWHT) or under intermediate to high temperature service conditions. This phenomenon is termed as Stress Relaxation Cracking (SRC) since the purpose of PWHT is to relieve the welding-induced residual stress fields, or as Stress Age Cracking (SAC) if failure happens during service. A leading literature explanation of this failure suggests that the residual stress relaxation and the precipitation dissolution and/or re-precipitation occur in the same temperature range, which can lead to locally high strains and thus to crack at the grain boundaries. Since in situ spatial measurements of residual stress fields, microstructural evolution, and failure processes are nearly infeasible, this work recourses to a micromechanical finite element framework that models the high temperature failure as the nucleation and growth of grain boundary cavities, whereas various parameters such as thermomechanical loading history and its evolution, the competition of grain-interior dislocation creep and grain-boundary diffusion in failure lifetime, and microstructural heterogeneities (such as the precipitate free zone near grain boundaries) can be quantitatively incorporated. It can be concluded from these microstructure-explicit simulations that an accurate knowledge of residual stress evolution and a carefully calibrated set of material constitutive parameters are the essential prerequisites for lifetime predictions. The understanding of individual governing factors also leads to a mechanistic interpretation of the observed SRC susceptibility C-curves. In conclusion, these results suggest that the criticality of residual stress evolution, but not the precipitation-induced local strains, be the leading factor for SRC.

347H stainless steel weldments↗
DOE OSTI2026

Gas permeation properties of amorphous zeolitic imidazolate framework membranes made by atomic/molecular layer deposition

Amorphous metal-organic framework (MOF) membranes are desirable because they may retain some of the molecular sieving properties of their crystalline counterparts while being free of grain boundary defects, which often hinder the consistent achievement of high membrane performance. However, current methods, like melting and compression, for fabricating amorphous MOF membranes involve multi-step processes that require the formation of a crystalline membrane first, that is then amorphized, and therefore, could be challenging to scale. Here, we utilize atomic/molecular layer deposition (ALD/MLD) of diethylzinc (DEZ) and 2-methylimidazole (2mIm) to directly synthesize ultrathin amorphous zeolitic imidazolate framework (aZIF) deposits on γ-alumina-coated α-alumina supports. As the number of ALD/MLD cycles increased from 10 to 300, gas permeances decreased while ideal selectivities increased. Mixture separation factors for C3H6/C3H8, CO2/N2, and H2/C3H8 as high as 4, 37, and 194, respectively, were obtained. At 200 °C and 2.5 bar equimolar feed of H2 and C3H8, a H2/C3H8 mixture separation factor of 185 is obtained with an H2 permeance of ca. 4.15x10-8 mol/m2-s-Pa (124 GPU). Additionally, the membrane achieves a CO2/N2 mixture separation factor of 37 at 25 °C and 1 bar with a CO2 permeance of ca. 3.96x10-8 mol/m2-s-Pa (118 GPU). Considering the vast array of compositionally distinct aZIFs that can be potentially deposited by this approach, an enormously large parameter space for membrane design is emerging to be explored.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Observations of Clouds, Aerosols, Precipitation, and Surface Radiation over the Southern Ocean: An Overview of CAPRICORN, MARCUS, MICRE, and SOCRATES

Weather and climate models are challenged by uncertainties and biases in simulating Southern Ocean (SO) radiative fluxes that trace to a poor understanding of cloud, aerosol, precipitation, and radiative processes, and their interactions. Projects between 2016 and 2018 used in situ probes, radar, lidar, and other instruments to make comprehensive measurements of thermodynamics, surface radiation, cloud, precipitation, aerosol, cloud condensation nuclei (CCN), and ice nucleating particles over the SO cold waters, and in ubiquitous liquid and mixed-phase clouds common to this pristine environment. Data including soundings were collected from the NSF–NCAR G-V aircraft flying north–south gradients south of Tasmania, at Macquarie Island, and on the R/V Investigator and RSV Aurora Australis. Synergistically these data characterize boundary layer and free troposphere environmental properties, and represent the most comprehensive data of this type available south of the oceanic polar front, in the cold sector of SO cyclones, and across seasons. Results show largely pristine environments with numerous small and few large aerosols above cloud, suggesting new particle formation and limited long-range transport from continents, high variability in CCN and cloud droplet concentrations, and ubiquitous supercooled water in thin, multilayered clouds, often with small-scale generating cells near cloud top. These observations demonstrate how cloud properties depend on aerosols while highlighting the importance of dynamics and turbulence that likely drive heterogeneity of cloud phase. Satellite retrievals confirmed low clouds were responsible for radiation biases. Furthermore, the combination of models and observations is examining how aerosols and meteorology couple to control SO water and energy budgets.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

SIERRA Multimechanics Module: Aria User Manual - Version 4.56

Aria is a Galerkin finite element based program for solving coupled-physics problems described by systems of PDEs and is capable of solving nonlinear, implicit, transient and direct-to-steady state problems in two and three dimensions on parallel architectures. The suite of physics currently supported by Aria includes thermal energy transport, species transport, and electrostatics as well as generalized scalar, vector and tensor transport equations. Additionally, Aria includes support for manufacturing process ows via the incompressible Navier-Stokes equations specialized to a low Reynolds number (Re < 1) regime. Enhanced modeling support of manufacturing processing is made possible through use of either arbitrary Lagrangian-Eulerian (ALE) and level set based free and moving boundary tracking in conjunction with quasi-static nonlinear elastic solid mechanics for mesh control. Coupled physics problems are solved in several ways including fully-coupled Newtons method with analytic or numerical sensitivities, fully-coupled Newton-Krylov methods and a loosely-coupled nonlinear iteration about subsets of the system that are solved using combinations of the aforementioned methods. Error estimation, uniform and dynamic h-adaptivity and dynamic load balancing are some of Arias more advanced capabilities.

97 MATHEMATICS AND COMPUTING↗