Search NASA⌕ Search

SEARCH · Search NASA

Results for “level sets”

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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

A coupled fluid-mechanical workflow to simulate the directed energy deposition additive manufacturing process

Simulation of additive manufacturing processes can provide essential insight into material behavior, residual stress, and ultimately, the performance of additively manufactured parts. In this work, we describe a new simulation based workflow utilizing both solid mechanics and fluid mechanics based formulations within the finite element software package SIERRA (Sierra Solid Mechanics Team in Sierra/Solid Mechanics 4.52 User’s Guide SAND2019-2715. Technical report, Sandia National Laboratories, 2011) to enable integrated simulations of directed energy deposition (DED) additive manufacturing processes. In this methodology, a high-fidelity fluid mechanics based model of additive manufacturing is employed as the first step in a simulation workflow. This fluid model uses a level set field to track the location of the boundary between the solid material and background gas and precisely predicts temperatures and material deposition shapes from additive manufacturing process parameters. Next, the resulting deposition shape and temperature field from the fluid model are then mapped into a solid mechanics formulation to provide a more accurate surface topology for radiation and convection boundary conditions and a prescribed temperature field. Solid mechanics simulations are then conducted to predict the evolution of material stresses and microstructure within a part. By combining thermal history and deposition shape from fluid mechanics with residual stress and material property evolutions from solid mechanics, additional fidelity and precision are incorporated into additive manufacturing process simulations providing new insight into complex DED builds.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Spread of Values of a Cantor-Type Fractal Continuous Nonmonotone Function

By using the Q{sub 5}{sup ∗}-representation of numbers[0,1]∋x=β{sub α{sub 1(x)1}}+∑{sub k=2}{sup ∞}(β{sub α{sub k(x)k}}∏{sub j=1}{sup k−1}q{sub α{sub j(x)j}})=Δ{sub α{sub 1(x)α{sub 2(x)…α{sub k(x)…}{sup Q{sub 5}{sup ∗}}}}} determined by the quinary alphabet A{sub 5} ≡ {0, 1, 2, 3, 4} and an infinite stochastic matrix ‖q{sub ik}‖, i ∈ A{sub 5}, k ∈ N, with positive elements (q{sub 0k} + q{sub 1k} + q{sub 2k} + q{sub 3k} + q{sub 4k} = 1) such that ∏{sub k=1}{sup ∞}max{sub i}{q_i_k}=0 and β{sub 0k} = 0, β{sub i + 1, k} = β{sub ik} + q{sub ik}, i= (0,4)-bar , we define a continuous Cantor-type function by the equalityf(Δ{sub α{sub 1…α{sub k…}{sup Q{sub 5}{sup ∗}}}})=δ{sub α{sub 1(x)1}}+∑{sub k=2}{sup ∞}(δ{sub α{sub k(x)k}}∑{sub j=1}{sup k−1}g{sub α{sub j(x)j}})≡Δ{sub α{sub 1(x)…α{sub k(x)…}{sup G}}}, where δ{sub 0n} = 0, δ{sub 1n}=((2+ε{sub n})/4), δ{sub 2n}=(2/4)=δ{sub 3n}, and δ{sub 4n}=((2−ε{sub n})/4), i.e., δ{sub i + 1, n} = δ{sub in} + g{sub in}, n ∈ N, and (ε{sub n}) is a given sequence of real numbers such that 0 ≤ ε{sub n} ≤ 1. We prove that this function is well defined and continuous. Moreover, it does not have intervals of monotonicity, except the intervals where it is constant. A criterion of bounded variation of the function is also established. We are especially interested in the problem of level sets of the function and in the topological and metric properties of the images of Cantor-type sets.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

The Importance of Viscous and Interfacial Forces in the Hydrodynamics of the Top-Submerged-Lance Furnace

The purpose of this work is to focus on the hydrodynamics of a Top-Submerged-Lance (TSL) smelting furnace, understanding how liquid properties and operational parameters act on key factors of a TSL process, such as splashing, mixing, mass transfer area, and bubble development. A deep knowledge of all those aspects is needed since they all influence the smelting reaction rates; hence the efficiency of the reactor. The characterization and scaling of the TSL gas injection are commonly based on the modified Froude number, the ratio of dynamic and gravitational forces. Detailed literature research reveals a potential weakness of this approach, since it does not consider the effects of viscosity and surface tension. To investigate this question an extensive parametric study was performed applying computational fluid dynamics to cold and non-reactive flows, which provided a broad overview of the physics of the flow. The analysis was performed on fluid dynamic properties (liquid density, liquid viscosity, surface tension) and operational variables (gas volume flow, lance immersion depth). The coupled Level Set—Volume of Fluid model, available in the commercial solver ANSYS Fluent{sup ®}, was used to resolve the gas–liquid interface in the multiphase flow. The results of the work underscore the significance of the viscous and interfacial forces for gas injection in smelting slags, confirming the incompleteness of applying only the Froude number to describe such flows.

36 MATERIALS SCIENCE↗
DOE OSTI2013

Elucidating geochemical response of shallow heterogeneous aquifers to CO2 leakage using high-performance computing: Implications for monitoring of CO2 sequestration

Predicting and quantifying impacts of potential carbon dioxide (CO2) leakage into shallow aquifers that overlie geologic CO2 storage formations is an important part of developing reliable carbon storage techniques. Leakage of CO2 through fractures, faults or faulty wellbores can reduce groundwater pH, inducing geochemical reactions that release solutes into the groundwater and pose a risk of degrading groundwater quality. In order to help quantify this risk, predictions of metal concentrations are needed during geologic storage of CO2. Here, we present regional-scale reactive transport simulations, at relatively fine-scale, of CO2 leakage into shallow aquifers run on the PFLOTRAN platform using high-performance computing. Multiple realizations of heterogeneous permeability distributions were generated using standard geostatistical methods. Increased statistical anisotropy of the permeability field resulted in more lateral and vertical spreading of the plume of impacted water, leading to increased Pb2+ (lead) concentrations and lower pH at a well down gradient of the CO2 leak. Pb2+ concentrations were higher in simulations where calcite was the source of Pb2+ compared to galena. The low solubility of galena effectively buffered the Pb2+ concentrations as galena reached saturation under reducing conditions along the flow path. In all cases, Pb2+ concentrations remained below the maximum contaminant level set by the EPA. Results from this study, compared to natural variability observed in aquifers, suggest that bicarbonate (HCO3) concentrations may be a better geochemical indicator of a CO2 leak under the conditions simulated here.

carbon sequestration↗
DOE OSTI2017

Numerical comparison of bubbling in a waste glass melter

Radioactive tank waste is scheduled for vitrification at the Hanford Tank Waste Treatment and Immobilization Plant, also known as the Vit Plant, being constructed at the Hanford Site. Testing of the pilot-scale, joule-heated DuraMelter 1200 at the Vitreous State Laboratory has shown that bubbling increases the melt rate, and as a result, melter throughput. Computational fluid dynamics (CFD) models of this pilot-scale waste glass melter were developed to improve our understanding of the processes that occur within the melter to aid in process optimization and troubleshooting of the Vit Plant melters. Unfortunately, model validation is complicated by the difficulty of obtaining suitable experimental data for operational melters attributed to the large-scale, radioactivity and high temperatures at which these melters are operated. This study focuses on confirming the fidelity of the CFD models to accurately reproduce the bubbling behavior. Because of the paucity of experimental data at the resolution required for CFD validation, a code-to-code comparison was used to evaluate two common approaches for simulating flows of immiscible fluids on numerical grids and resolving multiphase interfaces. Here, the volume of fluid and level set methods are used to track the dynamically evolving interfaces between the molten glass and the air bubbles. To aid in the verification of the results of these codes, a comparison of the bubble behavior, growth, and frequency of bubble generation are presented and a grid convergence study is performed for the two approaches. In addition to comparing the numerical results to available experimental data and observations from pilot-scale testing, the results are evaluated in relation to laboratory experiments with bubbles injected into viscous fluids, empirical correlations obtained from the published literature, and the basic laws of fluid dynamics. Furthermore, the results of this work can assist in validating waste glass melter models and provide a better understanding of the flow patterns within the WTP melters.

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

Fuel Property Effects on Knock Propensity and Thermal Efficiency in a Direct-Injection Spark-Ignition Engine

Engine knock remains one of the major barriers to further improvement in thermal efficiency of Direct-Injection Spark-Ignition (DISI) engines. While Research Octane Number and Motor Octane Number are often used as standard rating methods for knock resistance of fuels, the impacts of other fuel properties on knock propensity in modern engines such as heat of vaporization (HoV) and laminar flame speed (LFS) require better understanding in order to co-optimize fuels and engine designs to achieve higher thermal efficiency and lower CO2 emission. In the present study, computational fluid dynamics (CFD) is used to model a boosted DISI engine with a focus on knock prediction and fuel property effects. A level-set G-equation model is employed to capture turbulent premixed combustion, and is coupled with a transported Livengood-Wu (L-W) integral approach to predict autoignition in the unburnt region. A criterion associated with the L-W integral is developed to accurately predict knock onset and knock-limited spark-advance. This model is then applied to a sensitivity analysis of HoV and LFS on knock tendency and thermal efficiency. The pressure-temperature trajectory framework is applied and extended to study the fuel effects on auto-ignition process in the engine. An existing efficiency-based merit function, which is derived from experiments for boosted SI engines, is evaluated and improved based on the current CFD results.

DISI↗
DOE OSTI2021

Symmetry-breaking in double gyroid block copolymer films by non-affine distortion

Soft-matter bicontinuous networks find a double gyroid structure from block copolymer (BCP) self-assembly. A gyroid structure composed of dissimilar blocks has proven its potential as a soft crystal capable of tuning structural periodicity and symmetry, of which the lattice dimension is variable with molecular weight of the polymer. Using an asymmetric polystyrene- b -poly(methyl methacrylate) (PS- b -PMMA), in this study we show that the self-assembled gyroid films formed via a solvent vapor annealing (SVA) process undergo unique structural distortion due to directional deformation immediately upon deswelling. During the SVA process with PS- b -PMMA films, transient cylinders developed from the as-cast morphology transform into a cubic gyroid structure in a swollen state. Rapid and spontaneous deswelling processes -the manners in which the films contract along the z-direction while retaining an enlarged lateral dimension of the cubic form -lead to triclinic gyroid structures with z-directional contraction ratios ( C z ) of 2.5 and 2.0, respectively. Our X-ray analysis reveals that the deswelling process of the swollen gyroid films produces a notable symmetry-breaking in non-affine gyroid structure that elicits several forbidden reflections such as {110} and {200} reflections. For further characterization of the symmetry-breaking, we delineate the structural features of noncubic gyroid films by computing electron-density difference maps assisted with X-ray measurements. Level-set approach is accordingly developed to quantitate the structural characteristics of the maps in terms of inversion symmetry-breaking, suggesting its possible application to optical Weyl photonic crystals.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Optimal design of fiber reinforced composite structures and their direct ink write fabrication

The majority of the current structural optimization software does not accommodate manufacturing constraints. Therefore, substantial modifications are imposed upon optimized designs to make them manufacturable and hence nonoptimal. We propose to optimize the design of composite structures that are amenable to Additive Manufacturing (AM). The printing process chosen in our study is based on Direct Ink Writing (DIW) in which short carbon fibers in a thermoset resin are extruded through a moving nozzle to build up a structure. Since the fibers are primarily aligned in the flow direction of the extrudate, the DIW printing trajectory influences the material properties of the composite structure. To accommodate this, the extrudate trajectory follows the contours of parameterized level-set functions. The parametrization allows us to prescribe the material properties and impose many DIW manufacturing constraints such as no-overlap, no-sag, minimum radius of curvature and continuity of the toolpaths. Ultimately, we obtain optimal manufacturable toolpaths that start and finish at a boundary. To minimize the fabrication time, we formulate the linking sequence of the toolpaths as a traveling salesman problem which we solve to obtain the shortest continuous toolpath per layer. Several examples illustrate the optimization procedure. Validation is also performed.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Three-dimensional explicit finite element formulation for shear localization with global tracking of embedded weak discontinuities

Shear localization is a damage mechanism characterized by the concentration of plastic deformation within thin bands of materials, often observed as a precursor to failure in metals. In this study, we present a three-dimensional explicit finite element formulation with embedded weak discontinuities for the treatment of shear localization under dynamic loading conditions. In this formulation, the onset of localization is detected via a material stability analysis suitable for rate-sensitive materials, and the continuity of propagating shear bands is ensured using a global tracking strategy involving a heat-conduction type boundary value problem, solved for a scalar level set function over the global domain of the problem. In addition, we propose a modified quadrature rule that is used to compute the contribution of an individual element to the global finite element arrays, taking into account the position of the embedded shear band within that element.

42 ENGINEERING↗
DOE OSTI2020

A reduced-order shifted boundary method for parametrized incompressible Navier–Stokes equations

Here, we investigate a projection-based reduced order model of the steady incompressible Navier–Stokes equations for moderate Reynolds numbers. In particular, we construct an “embedded” reduced basis space, by applying proper orthogonal decomposition to the Shifted Boundary Method, a high-fidelity embedded method recently developed. We focus on the geometrical parametrization through level-set geometries, using a fixed Cartesian background geometry and the associated mesh. This approach avoids both remeshing and the development of a reference domain formulation, as typically done in fitted mesh finite element formulations. Two-dimensional computational examples for one and three parameter dimensions are presented to validate the convergence and the efficacy of the proposed approach.

"embedded" basis↗
DOE OSTI2021

Optimal toolpath design of additive manufactured composite cylindrical structures

Structural design optimization for additive manufacturing is primarily focused on planar layer-by-layer processes and design of composite cylindrical structures does not often accommodate manufacturing constraints. In this study, we propose to optimize the toolpath trajectory of additively manufactured composite cylinders comprised of multiple thin cylindrical annuli. Additionally, our printing process is based on direct ink writing in which short carbon fibers in a thermoset resin are extruded through a moving nozzle to build up a structure. Since the fibers are primarily aligned in the flow direction of the extrudate, the printing trajectory influences the material properties of the composite structure. To accommodate this, we define the toolpaths as contours of parameterized level-set functions. This parameterization allows us to optimize the material properties and impose manufacturing constraints such as no-overlap, no-sag, minimum radius of curvature, and continuity of the toolpaths. Several examples illustrate the optimization procedure.

42 ENGINEERING↗
DOE OSTI2013

Parallel adaptive fluid-structure interaction simulation of explosions impacting on building structures

We pursue a level set approach to couple an Eulerian shock-capturing fluid solver with space-time refinement to an explicit solid dynamics solver for large deformations and fracture. The coupling algorithms considering recursively finer fluid time steps as well as overlapping solver updates are discussed in detail. Our ideas are implemented in the AMROC adaptive fluid solver framework and are used for effective fluid-structure coupling to the general purpose solid dynamics code DYNA3D. Beside simulations verifying the coupled fluid-structure solver and assessing its parallel scalability, the detailed structural analysis of a reinforced concrete column under blast loading and the simulation of a prototypical blast explosion in a realistic multistory building are presented.

Deiterding, Ralf↗
DOE OSTI2017

A broadly-applicable unified closure relation for Taylor bubble rise velocity in pipes with stagnant liquid

Taylor bubble velocity in slug flow is a closure relation which significantly affects the prediction of liquid holdup (or void fraction) and pressure gradient in mechanistic models of slug flow for oil and gas pipe applications. In this work, we use a validated Computational Fluid Dynamics (CFD) approach to simulate the motion of Taylor bubbles in pipes; the interface is tracked with a Level-Set method implemented in a commercial code. A large numerical database is generated covering the most ample range of fluid properties and pipe inclination angles explored to date (Eo ∊ [10, 700],Mo ∊ [1 x 10 -6 , 5 x 10 3 , and θ ∊ [0°, 90°]). A unified Taylor bubble rise velocity correlation is extracted from that database. The new correlation predicts the numerical database with 8.6% absolute average relative error and a coefficient of determination R 2 = 0.97 and other available experimental data with 13.0% absolute average relative error and R 2 = 0.84 outperforming existing correlations and models.

Mechanics↗
DOE OSTI2013

Eulerian adaptive finite-difference method for high-velocity impact and penetration problems

Owing to the complex processes involved, faithful prediction of high-velocity impact events demands a simulation method delivering efficient calculations based on comprehensively formulated constitutive models. Such an approach is presented herein, employing a weighted essentially non-oscillatory (WENO) method within an adaptive mesh refinement (AMR) framework for the numerical solution of hyperbolic partial differential equations. Applied widely in computational fluid dynamics, these methods are well suited to the involved locally non-smooth finite deformations, circumventing any requirement for artificial viscosity functions for shock capturing. Application of the methods is facilitated through using a model of solid dynamics based upon hyper-elastic theory comprising kinematic evolution equations for the elastic distortion tensor. The model for finite inelastic deformations is phenomenologically equivalent to Maxwell s model of tangential stress relaxation. Closure relations tailored to the expected high-pressure states are proposed and calibrated for the materials of interest. Sharp interface resolution is achieved by employing level-set functions to track boundary motion, along with a ghost material method to capture the necessary internal boundary conditions for material interactions and stress-free surfaces. The approach is demonstrated for the simulation of high velocity impacts of steel projectiles on aluminium target plates in two and three dimensions.

Barton, Philip T.↗
DOE OSTI2014

Phase field method to optimize dielectric devices for electromagnetic wave propagation

We discuss a phase field method for shape optimization in the context of electromagnetic wave propagation. The proposed method has the same functional capabilities as the level set method for shape optimization. The first advantage of the method is the simplicity of computation, since extra operations such as re-initialization of functions are not required. The second is compatibility with the topology optimization method due to the similar domain representation and the sensitivity analysis. Structural shapes are represented by the phase field function defined in the design domain, and this function is optimized by solving a time-dependent reaction diffusion equation. The artificial double-well potential function used in the equation is derived from sensitivity analysis. We study four types of 2D or 2.5D (axisymmetric) optimization problems. Two are the classical problems of photonic crystal design based on the Bloch theory and photonic crystal wave guide design, and two are the recent topics of designing dielectric left-handed metamaterials and dielectric ring resonators.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Path optimization with limited sensing ability

We propose a computational strategy to find the optimal path for a mobile sensor with limited coverage to traverse a cluttered region. The goal is to find one of the shortest feasible paths to achieve the complete scan of the environment. We pose the problem in the level set framework, and first consider a related question of placing multiple stationary sensors to obtain the full surveillance of the environment. By connecting the stationary locations using the nearest neighbor strategy, we form the initial guess for the path planning problem of the mobile sensor. Then the path is optimized by reducing its length, via solving a system of ordinary differential equations (ODEs), while maintaining the complete scan of the environment. Furthermore, we use intermittent diffusion, which converts the ODEs into stochastic differential equations (SDEs), to find an optimal path whose length is globally minimal. To improve the computation efficiency, we introduce two techniques, one to remove redundant connecting points to reduce the dimension of the system, and the other to deal with the entangled path so the solution can escape the local traps. Numerical examples are shown to illustrate the effectiveness of the proposed method.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Self-consistent field theory simulations of polymers on arbitrary domains

We introduce a framework for simulating the mesoscale self-assembly of block copolymers in arbitrary confined geometries subject to Neumann boundary conditions. We employ a hybrid finite difference/volume approach to discretize the mean-field equations on an irregular domain represented implicitly by a level-set function. The numerical treatment of the Neumann boundary conditions is sharp, i.e. it avoids an artificial smearing in the irregular domain boundary. This strategy enables the study of self-assembly in confined domains and enables the computation of physically meaningful quantities at the domain interface. In addition, we employ adaptive grids encoded with Quad-/Oc-trees in parallel to automatically refine the grid where the statistical fields vary rapidly as well as at the boundary of the confined domain. This approach results in a significant reduction in the number of degrees of freedom and makes the simulations in arbitrary domains using effective boundary conditions computationally efficient in terms of both speed and memory requirement. Finally, in the case of regular periodic domains, where pseudo-spectral approaches are superior to finite differences in terms of CPU time and accuracy, we use the adaptive strategy to store chain propagators, reducing the memory footprint without loss of accuracy in computed physical observables.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2017

A sharp interface method for compressible liquid–vapor flow with phase transition and surface tension

The numerical approximation of non-isothermal liquid–vapor flow within the compressible regime is a difficult task because complex physical effects at the phase interfaces can govern the global flow behavior. We present a sharp interface approach which treats the interface as a shock-wave like discontinuity. Any mixing of fluid phases is avoided by using the flow solver in the bulk regions only, and a ghost-fluid approach close to the interface. The coupling states for the numerical solution in the bulk regions are determined by the solution of local two-phase Riemann problems across the interface. The Riemann solution accounts for the relevant physics by enforcing appropriate jump conditions at the phase boundary. A wide variety of interface effects can be handled in a thermodynamically consistent way. This includes surface tension or mass/energy transfer by phase transition. Moreover, the local normal speed of the interface, which is needed to calculate the time evolution of the interface, is given by the Riemann solution. The interface tracking itself is based on a level-set method. The focus in this paper is the description of the two-phase Riemann solver and its usage within the sharp interface approach. One-dimensional problems are selected to validate the approach. Finally, the three-dimensional simulation of a wobbling droplet and a shock droplet interaction in two dimensions are shown. In both problems phase transition and surface tension determine the global bulk behavior.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗