skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Deep level-set method for Stefan problems
We propose a level-set approach to characterize the region occupied by the solid in Stefan problems with and without surface tension, based on their recent probabilistic reformulation. The level-set function is parameterized by a feed-forward neural network, whose parameters are trained using the probabilistic formulation of the Stefan growth condition. The algorithm can handle Stefan problems where the liquid is supercooled and can capture surface tension effects through the simulation of particles along the moving boundary together with an efficient approximation of the mean curvature. We demonstrate the effectiveness of the method on a variety of examples with and without radial symmetry.  more » « less
Award ID(s):
2108680
PAR ID:
10519297
Author(s) / Creator(s):
; ;
Publisher / Repository:
Elsevier
Date Published:
Journal Name:
Journal of Computational Physics
Volume:
503
Issue:
C
ISSN:
0021-9991
Page Range / eLocation ID:
112828
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract We present a probabilistic extension of action language $${\cal BC}$$+$$ . Just like $${\cal BC}$$+$$ is defined as a high-level notation of answer set programs for describing transition systems, the proposed language, which we call p $${\cal BC}$$+$$ , is defined as a high-level notation of LP MLN programs—a probabilistic extension of answer set programs. We show how probabilistic reasoning about transition systems, such as prediction, postdiction, and planning problems, as well as probabilistic diagnosis for dynamic domains, can be modeled in p $${\cal BC}$$+$$ and computed using an implementation of LP MLN . 
    more » « less
  2. Assessing the effects of input uncertainty on simulation results for multiphase flows will allow for more robust engineering designs and improved devices. For example, in atomizing jets, surface tension plays a critical role in determining when and how coherent liquid structures break up. Uncertainty in the surface tension coefficient can lead to uncertainty in spray angle, drop size, and velocity distribution. Uncertainty quantification (UQ) determines how input uncertainties affect outputs, and the approach taken can be classified as non-intrusive or intrusive. A classical, non-intrusive approach is the Monte-Carlo scheme, which requires multiple simulation runs using samples from a distribution of inputs. Statistics on output variability are computed from the many simulation outputs. While non-intrusive schemes are straightforward to implement, they can quickly become cost prohibitive, suffer from convergence issues, and have problems with confounding factors, making it difficult to look at uncertainty in multiple variables at once. Alternatively, an intrusive scheme inserts stochastic (uncertain) variables into the governing equations, modifying the mathematics and numerical methods used, but possibly reducing computational cost. In this work, we extend UQ methods developed for single-phase flows to handle gas-liquid multiphase dynamics by developing a stochastic conservative level set approach and a stochastic continuous surface tension method. An oscillating droplet and a 2-D atomizing jet are used to test the method. In these test cases, uncertainty about the surface tension coefficient and initial starting position will be explored, including the impact on breaking/ merging interfaces. 
    more » « less
  3. Abstract In this paper we study a finite‐depth layer of viscous incompressible fluid in dimension , modeled by the Navier‐Stokes equations. The fluid is assumed to be bounded below by a flat rigid surface and above by a free, moving interface. A uniform gravitational field acts perpendicularly to the flat surface, and we consider the cases with and without surface tension acting on the free interface. In addition to these gravity‐capillary effects, we allow for a second force field in the bulk and an external stress tensor on the free interface, both of which are posited to be in traveling wave form, i.e., time‐independent when viewed in a coordinate system moving at a constant velocity parallel to the rigid lower boundary. We prove that, with surface tension in dimension and without surface tension in dimension , for every nontrivial traveling velocity there exists a nonempty open set of force and stress data that give rise to traveling wave solutions. While the existence of inviscid traveling waves is well‐known, to the best of our knowledge this is the first construction of viscous traveling wave solutions. Our proof involves a number of novel analytic ingredients, including: the study of an overdetermined Stokes problem and its underdetermined adjoint problem, a delicate asymptotic development of the symbol for a normal‐stress to normal‐Dirichlet map defined via the Stokes operator, a new scale of specialized anisotropic Sobolev spaces, and the study of a pseudodifferential operator that synthesizes the various operators acting on the free surface functions. © 2022 The Authors.Communications on Pure and Applied Mathematicspublished by Wiley Periodicals LLC. 
    more » « less
  4. Integrated task and motion planning (TAMP) has proven to be a valuable approach to generalizable long-horizon robotic manipulation and navigation problems. However, the typical TAMP problem formulation assumes full observability and deterministic action effects. These assumptions limit the ability of the planner to gather information and make decisions that are risk-aware. We propose a strategy for TAMP with Uncertainty and Risk Awareness (TAMPURA) that is capable of efficiently solving long-horizon planning problems with initial- state and action outcome uncertainty, including problems that require information gathering and avoiding undesirable and irreversible outcomes. Our planner reasons under uncertainty at both the abstract task level and continuous controller level. Given a set of closed-loop goal-conditioned controllers operating in the primitive action space and a description of their preconditions and potential capabilities, we learn a high-level abstraction that can be solved efficiently and then refined to continuous actions for execution. We demonstrate our approach on several robotics problems where uncertainty is a crucial factor and show that reasoning under uncertainty in these problems outperforms previously proposed determinized planning, direct search, and reinforcement learning strategies. Lastly, we demonstrate our planner on two real-world robotics problems using recent ad- vancements in probabilistic perception. 
    more » « less
  5. We develop a consistent adaptive framework in a multilevel collocated grid layout for simulating two-phase flows with adaptive mesh refinement (AMR). The conservative mo-mentum equations and the mass equation are solved in the present consistent framework. This consistent mass and momentum transport treatment greatly improves the accuracy and robustness for simulating two-phase flows with a high density ratio and high Reynolds number. The interface capturing level set method is coupled with the conservative form of the Navier–Stokes equations, and the multilevel reinitialization technique is applied for mass conservation. This adaptive framework allows us to advance all variables level by level using either the subcycling or the non-subcycling method to decouple the data ad-vancement on each level. The accuracy and robustness of the framework are validated by a variety of canonical two-phase flow problems. We demonstrate that the consistent scheme results in a numerically stable solution in flows with high density ratios(up to 106) and high Reynolds numbers(up to 106), while the inconsistent scheme exhibits non-physical fluid behaviors in these tests. Furthermore, it is shown that the subcycling and non-subcycling methods provide consistent results and that both of them can accurately resolve the interfaces of the two-phase flows with surface tension effects. Finally, a 3D breaking wave problem is simulated to show the efficiency and significant speedup of the proposed framework using AMR. 
    more » « less