Use of Aria to simulate laser weld pool dynamics for neutron generator production.

This report documents the results for the FY07 ASC Integrated Codes Level 2 Milestone number 2354. The description for this milestone is, 'Demonstrate level set free surface tracking capabilities in ARIA to simulate the dynamics of the formation and time evolution of a weld pool in laser welding applications for neutron generator production'. The specialized boundary conditions and material properties for the laser welding application were implemented and verified by comparison with existing, two-dimensional applications. Analyses of stationary spot welds and traveling line welds were performed and the accuracy of the three-dimensional (3D) level set algorithm is assessed by comparison with 3D moving mesh calculations

Aria 1.5 : user manual.

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 the incompressible Navier-Stokes equations, energy transport equation, species transport equations, nonlinear elastic solid mechanics, and electrostatics as well as generalized scalar, vector and tensor transport equations. Additionally, Aria includes support for arbitrary Lagrangian-Eulerian (ALE) and level set based free and moving boundary tracking. Coupled physics problems are solved in several ways including fully-coupled Newton's method with analytic or numerical sensitivities, fully-coupled Newton-Krylov methods, fully-coupled Picard's method, 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 Aria's more advanced capabilities. Aria is based on the Sierra Framework

Large deformation solid-fluid interaction via a level set approach.

Solidification and blood flow seemingly have little in common, but each involves a fluid in contact with a deformable solid. In these systems, the solid-fluid interface moves as the solid advects and deforms, often traversing the entire domain of interest. Currently, these problems cannot be simulated without innumerable expensive remeshing steps, mesh manipulations or decoupling the solid and fluid motion. Despite the wealth of progress recently made in mechanics modeling, this glaring inadequacy persists. We propose a new technique that tracks the interface implicitly and circumvents the need for remeshing and remapping the solution onto the new mesh. The solid-fluid boundary is tracked with a level set algorithm that changes the equation type dynamically depending on the phases present. This novel approach to coupled mechanics problems promises to give accurate stresses, displacements and velocities in both phases, simultaneously

Solution-verified reliability analysis and design of bistable MEMS using error estimation and adaptivity.

This report documents the results for an FY06 ASC Algorithms Level 2 milestone combining error estimation and adaptivity, uncertainty quantification, and probabilistic design capabilities applied to the analysis and design of bistable MEMS. Through the use of error estimation and adaptive mesh refinement, solution verification can be performed in an automated and parameter-adaptive manner. The resulting uncertainty analysis and probabilistic design studies are shown to be more accurate, efficient, reliable, and convenient

Modeling injection molding of net-shape active ceramic components.

To reduce costs and hazardous wastes associated with the production of lead-based active ceramic components, an injection molding process is being investigated to replace the current machining process. Here, lead zirconate titanate (PZT) ceramic particles are suspended in a thermoplastic resin and are injected into a mold and allowed to cool. The part is then bisque fired and sintered to complete the densification process. To help design this new process we use a finite element model to describe the injection molding of the ceramic paste. Flow solutions are obtained using a coupled, finite-element based, Newton-Raphson numerical method based on the GOMA/ARIA suite of Sandia flow solvers. The evolution of the free surface is solved with an advanced level set algorithm. This approach incorporates novel methods for representing surface tension and wetting forces that affect the evolution of the free surface. Thermal, rheological, and wetting properties of the PZT paste are measured for use as input to the model. The viscosity of the PZT is highly dependent both on temperature and shear rate. One challenge in modeling the injection process is coming up with appropriate constitutive equations that capture relevant phenomenology without being too computationally complex. For this reason we model the material as a Carreau fluid and a WLF temperature dependence. Two-dimensional (2D) modeling is performed to explore the effects of the shear in isothermal conditions. Results indicate that very low viscosity regions exist near walls and that these results look similar in terms of meniscus shape and fill times to a simple Newtonian constitutive equation at the shear-thinned viscosity for the paste. These results allow us to pick a representative viscosity to use in fully three-dimensional (3D) simulation, which because of numerical complexities are restricted to using a Newtonian constitutive equation. Further 2D modeling at nonisothermal conditions shows that the choice of representative Newtonian viscosity is dependent on the amount of heating of the initially room temperature mold. An early 3D transient model shows that the initial design of the distributor is sub-optimal. However, these simulations take several months to run on 4 processors of an HP workstation using a preconditioner/solver combination of ILUT/GMRES with fill factors of 3 and PSPG stabilization. Therefore, several modifications to the distributor geometry and orientations of the vents and molds have been investigated using much faster 3D steady-state simulations. The pressure distribution for these steady-state calculations is examined for three different distributor designs to see if this can indicate which geometry has the superior design. The second modification, with a longer distributor, is shown to have flatter, more monotonic isobars perpendicular to the flow direction indicating a better filling process. The effects of the distributor modifications, as well as effects of the mold orientation, have also been examined with laboratory experiments in which the flow of a viscous Newtonian oil entering transparent molds is recorded visually. Here, the flow front is flatter and voids are reduced for the second geometry compared to the original geometry. A horizontal orientation, as opposed to the planned vertical orientation, results in fewer voids. Recently, the Navier-Stokes equations have been stabilized with the Dohrman-Bochev PSPP stabilization method, allowing us to calculate transient 3D simulations with computational times on the order of days instead of months. Validation simulations are performed and compared to the experiments. Many of the trends of the experiments are captured by the level set modeling, though quantitative agreement is lacking mainly due to the high value of the gas phase viscosity necessary for numerical stability, though physically unrealistic. More correct trends are predicted for the vertical model than the horizontal model, which is serendipitous as the actual mold is held in a vertical geometry. The full, transient mold filling calculations indicate that the flow front is flatter and voids may be reduced for the second geometry compared to the original geometry. The validated model is used to predict mold filling for the actual process with the material properties for the PZT paste, the original distributor geometry, and the mold in a vertical orientation. This calculation shows that voids may be trapped at the four corners of the mold opposite the distributor

Observation of Scaling Violations in Scaled Momentum Distributions at HERA

Charged particle production has been measured in deep inelastic scattering (DIS) events over a large range of $x$ and $Q^2$ using the ZEUS detector. The evolution of the scaled momentum, $x_p$, with $Q^2,$ in the range 10 to 1280 $GeV^2$, has been investigated in the current fragmentation region of the Breit frame. The results show clear evidence, in a single experiment, for scaling violations in scaled momenta as a function of $Q^2$.Comment: 21 pages including 4 figures, to be published in Physics Letters B. Two references adde

Computational thermal, chemical, fluid, and solid mechanics for geosystems management.

This document summarizes research performed under the SNL LDRD entitled - Computational Mechanics for Geosystems Management to Support the Energy and Natural Resources Mission. The main accomplishment was development of a foundational SNL capability for computational thermal, chemical, fluid, and solid mechanics analysis of geosystems. The code was developed within the SNL Sierra software system. This report summarizes the capabilities of the simulation code and the supporting research and development conducted under this LDRD. The main goal of this project was the development of a foundational capability for coupled thermal, hydrological, mechanical, chemical (THMC) simulation of heterogeneous geosystems utilizing massively parallel processing. To solve these complex issues, this project integrated research in numerical mathematics and algorithms for chemically reactive multiphase systems with computer science research in adaptive coupled solution control and framework architecture. This report summarizes and demonstrates the capabilities that were developed together with the supporting research underlying the models. Key accomplishments are: (1) General capability for modeling nonisothermal, multiphase, multicomponent flow in heterogeneous porous geologic materials; (2) General capability to model multiphase reactive transport of species in heterogeneous porous media; (3) Constitutive models for describing real, general geomaterials under multiphase conditions utilizing laboratory data; (4) General capability to couple nonisothermal reactive flow with geomechanics (THMC); (5) Phase behavior thermodynamics for the CO2-H2O-NaCl system. General implementation enables modeling of other fluid mixtures. Adaptive look-up tables enable thermodynamic capability to other simulators; (6) Capability for statistical modeling of heterogeneity in geologic materials; and (7) Simulator utilizes unstructured grids on parallel processing computers

Wetting and free surface flow modeling for potting and encapsulation.

As part of an effort to reduce costs and improve quality control in encapsulation and potting processes the Technology Initiative Project ''Defect Free Manufacturing and Assembly'' has completed a computational modeling study of flows representative of those seen in these processes. Flow solutions are obtained using a coupled, finite-element-based, numerical method based on the GOMA/ARIA suite of Sandia flow solvers. The evolution of the free surface is solved with an advanced level set algorithm. This approach incorporates novel methods for representing surface tension and wetting forces that affect the evolution of the free surface. In addition, two commercially available codes, ProCAST and MOLDFLOW, are also used on geometries representing encapsulation processes at the Kansas City Plant. Visual observations of the flow in several geometries are recorded in the laboratory and compared to the models. Wetting properties for the materials in these experiments are measured using a unique flowthrough goniometer

Administration of vitamin D and its metabolites in critically ill adult patients: an updated systematic review with meta-analysis of randomized controlled trials

Background The clinical significance of vitamin D administration in critically ill patients remains inconclusive. The purpose of this systematic review with meta-analysis was to investigate the effect of vitamin D and its metabolites on major clinical outcomes in critically ill patients, including a subgroup analysis based on vitamin D status and route of vitamin D administration. Methods Major databases were searched through February 9, 2022. Randomized controlled trials of adult critically ill patients with an intervention group receiving vitamin D or its metabolites were included. Random-effect meta-analyses were performed to estimate the pooled risk ratio (dichotomized outcomes) or mean difference (continuous outcomes). Risk of bias assessment included the Cochrane tool for assessing risk of bias in randomized trials. Results Sixteen randomized clinical trials with 2449 patients were included. Vitamin D administration was associated with lower overall mortality (16 studies: risk ratio 0.78, 95% confidence interval 0.62–0.97, p = 0.03; I2 = 30%), reduced intensive care unit length of stay (12 studies: mean difference − 3.13 days, 95% CI − 5.36 to − 0.89, n = 1250, p = 0.006; I2 = 70%), and shorter duration of mechanical ventilation (9 studies: mean difference − 5.07 days, 95% CI − 7.42 to − 2.73, n = 572, p < 0.0001; I2 = 54%). Parenteral administration was associated with a greater effect on overall mortality than enteral administration (test of subgroup differences, p = 0.04), whereas studies of parenteral subgroups had lower quality. There were no subgroup differences based on baseline vitamin D levels. Conclusions Vitamin D supplementation in critically ill patients may reduce mortality. Parenteral administration might be associated with a greater impact on mortality. Heterogeneity and assessed certainty among the studies limits the generalizability of the results

Omega-6 sparing effects of parenteral lipid emulsions-an updated systematic review and meta-analysis on clinical outcomes in critically ill patients

Background Parenteral lipid emulsions in critical care are traditionally based on soybean oil (SO) and rich in pro-inflammatory omega-6 fatty acids (FAs). Parenteral nutrition (PN) strategies with the aim of reducing omega-6 FAs may potentially decrease the morbidity and mortality in critically ill patients. Methods A systematic search of MEDLINE, EMBASE, CINAHL and CENTRAL was conducted to identify all randomized controlled trials in critically ill patients published from inception to June 2021, which investigated clinical omega-6 sparing effects. Two independent reviewers extracted bias risk, treatment details, patient characteristics and clinical outcomes. Random effect meta-analysis was performed. Results 1054 studies were identified in our electronic search, 136 trials were assessed for eligibility and 26 trials with 1733 critically ill patients were included. The median methodologic score was 9 out of 14 points (95% confidence interval [CI] 7, 10). Omega-6 FA sparing PN in comparison with traditional lipid emulsions did not decrease overall mortality (20 studies; risk ratio [RR] 0.91; 95% CI 0.76, 1.10; p = 0.34) but hospital length of stay was substantially reduced (6 studies; weighted mean difference [WMD] − 6.88; 95% CI − 11.27, − 2.49; p = 0.002). Among the different lipid emulsions, fish oil (FO) containing PN reduced the length of intensive care (8 studies; WMD − 3.53; 95% CI − 6.16, − 0.90; p = 0.009) and rate of infectious complications (4 studies; RR 0.65; 95% CI 0.44, 0.95; p = 0.03). When FO was administered as a stand-alone medication outside PN, potential mortality benefits were observed compared to standard care. Conclusion Overall, these findings highlight distinctive omega-6 sparing effects attributed to PN. Among the different lipid emulsions, FO in combination with PN or as a stand-alone treatment may have the greatest clinical impact