Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “computational mathematics”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Kibble-Zurek mechanism for nonequilibrium phase transitions in driven systems with quenched disorder

Abstract We show that the Kibble-Zurek mechanism applies to nonequilibrium phase transitions found in driven assemblies of superconducting vortices and colloids moving over quenched disorder where a transition occurs from a plastic disordered flowing state to a moving anisotropic crystal. We measure the density of topological defects as a function of quench rate through the nonequilibrium phase transition, and find that on the ordered side of the transition, the topological defect density ρ d scales as a power law, $${\rho }_{{{{{{{{\rm{d}}}}}}}}}\propto 1/{t}_{{{{{{{{\rm{q}}}}}}}}}^{\beta }$$ ρ d ∝ 1 / t q β , where t q is the quench time duration, consistent with the Kibble-Zurek mechanism. We show that scaling with the same exponent holds for varied strengths of quenched disorder and that the exponents fall in the directed percolation universality class. Our results suggest that the Kibble-Zurek mechanism can be applied to the broader class of systems that exhibit absorbing phase transitions.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

A framework for data-driven solution and parameter estimation of PDEs using conditional generative adversarial networks

We employ and adapt the image-to-image translation concept based on conditional generative adversarial networks (cGAN) for learning a forward and an inverse solution operator of partial differential equations (PDEs). We focus on steady-state solutions of coupled hydromechanical processes in heterogeneous porous media and present the parameterization of the spatially heterogeneous coefficients, which is exceedingly difficult using standard reduced-order modeling techniques. We show that our framework provides a speed-up of at least 2,000 times compared to a finite-element solver and achieves a relative root-mean-square error (r.m.s.e.) of less than 2% for forward modeling. For inverse modeling, the framework estimates the heterogeneous coefficients, given an input of pressure and/or displacement fields, with a relative r.m.s.e. of less than 7%, even for cases where the input data are incomplete and contaminated by noise. The framework also provides a speed-up of 120,000 times compared to a Gaussian prior-based inverse modeling approach while also delivering more accurate results.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Dynamical commensuration effect in a two-dimensional Yukawa solid modulated by periodic substrates

Transverse depinning dynamics of a periodic-square-substrate modulated two-dimensional dusty plasma solid driven by a constant force in the longitudinal direction are investigated using Langevin dynamical simulations. When the commensuration ratio (the number ratio of particles to substrate potential wells) is increased, there is a nonmonotonic variation of the critical transverse depinning force, and the local maxima and minima of the critical transverse depinning force precisely correspond to the dynamical commensurate and incommensurate conditions, respectively. Here, the dynamical commensuration effect is also clearly visible in the stable one-dimensional channel particle trajectories and the highly ordered structure, while both the particle trajectories and the structure are more disordered under the incommensurate conditions. The nonmonotonic variation of the critical transverse depinning force is attributed to the stability of the lattice structure at specific commensuration ratios.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Behavior of solitary waves of coupled nonlinear Schrödinger equations subjected to complex external periodic potentials with odd-$\mathcal{PT}$ symmetry

In this work, we discuss the response of both moving and trapped solitary wave solutions of a two-component nonlinear Schrödinger system in 1 + 1 dimensions to an odd-$\mathcal{PT}$ external periodic complex potential. The dynamical behavior of perturbed solitary waves is explored by conducting numerical simulations of the nonlinear system and using a collective coordinate variational approximation. We present case examples corresponding to choices of parameter values and initial conditions involved therein. The results of the collective coordinate approximation are compared against numerical simulations where we observe qualitatively good agreement between the two. Unlike the case for a single-component solitary wave in a complex periodic $\mathcal{PT}$-symmetric potential, the collective coordinate equations do not have a small oscillation regime, and initially the height of the two components changes in opposite directions often causing instability. We find that the dynamic stability criteria we have used in the one-component case are a good indicator for the onset of dynamic instabilities in the present setup.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Nonlinear dimension reduction for surrogate modeling using gradient information

We introduce a method for the nonlinear dimension reduction of a high-dimensional function $u:{\mathbb{R}}^d\rightarrow{\mathbb{R}}$, $d\gg 1$. Our objective is to identify a nonlinear feature map $g:{\mathbb{R}}^d\rightarrow{\mathbb{R}}^m$, with a prescribed intermediate dimension $m\ll d$, so that $u$ can be well approximated by $f\circ g$ for some profile function $f:{\mathbb{R}}^m\rightarrow{\mathbb{R}}$. We propose to build the feature map by aligning the Jacobian $\nabla g$ with the gradient $\nabla u$, and we theoretically analyze the properties of the resulting $g$. Once $g$ is built, we construct $f$ by solving a gradient-enhanced least squares problem. Our practical algorithm uses a sample $\{{\textbf{x}}^{(i)},u({\textbf{x}}^{(i)}),\nabla u({\textbf{x}}^{(i)})\}_{i=1}^N$ and builds both $g$ and $f$ on adaptive downward-closed polynomial spaces, using cross validation to avoid overfitting. We numerically evaluate the performance of our algorithm across different benchmarks, and explore the impact of the intermediate dimension $m$. We show that building a nonlinear feature map $g$ can permit more accurate approximation of $u$ than a linear $g$, for the same input data set.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Transition from susceptible-infected to susceptible-infected-recovered dynamics in a susceptible-cleric-zombie-recovered active matter model

The susceptible-infected (SI) and susceptible-infected-recovered (SIR) models provide two distinct representations of epidemic evolution, distinguished by whether or not the number of susceptibles always drops to zero at long times. Here, in this study, we introduce a new active matter epidemic model, the “susceptible-cleric-zombie-recovered” (SCZR) model, in which spontaneous recovery is absent but zombies can recover with probability γ via interaction with a cleric. Upon colliding with a zombie, both susceptibles and clerics enter the zombie state with probability β and α, respectively. By changing the initial fraction of clerics or their healing ability rate γ, we can tune the SCZR model between SI dynamics, in which no susceptibles or clerics remain at long times, and SIR dynamics, in which a finite number of clerics and susceptibles survive at long times. The model is relevant to certain real world diseases such as HIV where spontaneous recovery is impossible but where medical interventions by a limited number of caregivers can reduce or eliminate the spread of infection.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2023

Characterizing different motility-induced regimes in active matter with machine learning and noise

Here we examine motility-induced phase separation (MIPS) in two-dimensional run-and-tumble disk systems using both machine learning and noise fluctuation analysis. Our measures suggest that within the MIPS state there are several distinct regimes as a function of density and run time, so that systems with MIPS transitions exhibit an active fluid, an active crystal, and a critical regime. The different regimes can be detected by combining an order parameter extracted from principal component analysis with a cluster stability measurement. The principal component-derived order parameter is maximized in the critical regime, remains low in the active fluid, and has an intermediate value in the active crystal regime. We demonstrate that machine learning can better capture dynamical properties of the MIPS regimes compared to more standard structural measures such as the maximum cluster size. The different regimes can also be characterized via changes in the noise power of the fluctuations in the average speed. In the critical regime, the noise power passes through a maximum and has a broad spectrum with a 1/f 1.6 signature, similar to the noise observed near depinning transitions or for solids undergoing plastic deformation.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Ecological information and approaches needed for risk communication dialogs for acute or chronic environmental crises

Scientists, social scientists, risk communicators, and many others are often thrust into a crisis situation where they need to interact with a range of stakeholders, including governmental personnel (tribal, U.S. federal, state, local), local residents, and other publics, as well as other scientists and other risk communicators in situations where information is incomplete and evolving. This paper provides: (1) an overall framework for thinking about communication during crises, from acute to chronic, and local to widespread, (2) a template for the types of ecological information needed to address public and environmental concerns, and (3) examples to illustrate how this information will aid risk communicators. The main goal is providing an approach to the knowledge needed by communicators to address the challenges of protecting ecological resources during an environmental crisis, or for an on-going, chronic environmental issue. To understand the risk to these ecological resources, it is important to identify the type of event, whether it is acute or chronic (or some combination of these), what receptors are at risk, and what stressors are involved (natural, biological, chemical, radiological). For ecological resources, the key information a communicator needs for a crisis is whether any of the following are present: threatened or endangered species, species of special concern, species groups of concern (e.g., neotropical bird migrants, breeding frogs in vernal ponds, rare plant assemblages), unique or rare habitats, species of commercial and recreational interest, and species/habitats of especial interest for medicinal, cultural, or religious activities. Communication among stakeholders is complicated with respect to risk to ecological receptors because of differences in trust, credibility, empathy, perceptions, world view valuation of the resources, and in many cases, a history of misinformation, disinformation, or no information. Exposure of salmon spawning in the Columbia River to hexavalent chromium from the Hanford Site is used as an example of communication challenges with different stakeholders, including Native Americans with Tribal Treaty rights to the land.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Trust and consequences: Role of community science, perceptions, values, and environmental justice in risk communication

Risk communication is often viewed as imparting information and perhaps as a two-way dialogue. Risk communication inadequacies on the part of both “communicator” and “community members” can lead to adverse consequences and amplify environmental justice disparities. The paper suggests a transformational approach where risk communicators must learn to trust community experts and their knowledge base (and act upon it), where risk information imparted by risk communicators addresses what communities are most concerned about (as well as risk from specific chemicals or radionuclides), and where risk information and assessments address underlying issues and disparities, as well as cultural traditions (among others). Providing risk probabilities is no longer sufficient; western science may not be enough, and community and native scientific knowledge is needed. Risk communication (or information transfer) for environmental risks that are ongoing usually applies to low-income, minority communities—people living in dense inner cities, rural communities, Native American communities—or to people living near a risky facility. Communication within this context requires mutual trust, listening and respect, as well as acceptance of indigenous and community knowledge as equally valuable. Examples are given to illustrate a community perspective.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

The Pele Simulation Suite for Reacting Flows at Exascale

In this work, we present the Pele suite of software tools for compressible and incompressible reacting flows. The Pele suite leverages several different libraries, notably AMReX and SUNDIALS, to achieve performance portability on heterogeneous computing architectures across the supercomputing landscape. The Pele suite is comprised of PeleC, a compressible reacting flow block-structured adaptive mesh refinement solver, PeleLMeX, a low-Mach number reacting flow block-structured adaptive mesh refinement solver, Pele-Physics, a library for transport, thermodynamics, finite rate chemistry, soot, spray and radiation physics. The objective of this paper is (i) to present the code development efforts necessary to achieve highly effective and scalable applications for exascale machines and (ii) to detail the performance results of the Combustion-Pele project applications on Oak Ridge National Laboratory's Frontier. We show good weak and strong scaling results for both PeleC and PeleLMeX up to more than 50 billion cells on more than 4096 Frontier graphics processing unit nodes. We also present a capability demonstration simulation of a dual-fuel pulse compression ignition engine (six adaptive mesh refinement levels, and 60 billion cells or 2.1 trillion degrees of freedom) on Frontier, to date one of the largest simulations performed on the first exascale-class supercomputer.

adaptive mesh refinement↗
DOE OSTI2021

Maximum Size Intersecting Families of Bounded Minimum Positive Co-degree

Let $\mathcal{H}$ be an $r$-uniform hypergraph. The minimum positive co-degree of $\mathcal{H}$, denoted by $\delta_{r-1}^+(\mathcal{H})$, is the minimum $k$ such that if $S$ is an $(r-1)$-set contained in a hyperedge of $\mathcal{H}$, then $S$ is contained in at least $k$ hyperedges of $\mathcal{H}$. For $r\geq k$ fixed and $n$ sufficiently large, we determine the maximum possible size of an intersecting $r$-uniform $n$-vertex hypergraph with minimum positive co-degree $\delta_{r-1}^+(\mathcal{H}) \geq k$ and characterize the unique hypergraph attaining this maximum. This generalizes the Erd\Hos--Ko--Rado theorem which corresponds to the case $k=1$. Furthermore, our proof is based on the delta-system method.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Adaptive Uncertainty Quantification for Stochastic Hyperbolic Conservation Laws

Here, we propose a predictor-corrector adaptive method for the study of hyperbolic partial differential equations (PDEs) under uncertainty. Constructed around the framework of stochastic finite volume (SFV) methods, our approach circumvents sampling schemes or simulation ensembles while also preserving fundamental properties, in particular hyperbolicity of the resulting systems and conservation of the discrete solutions. Furthermore, we augment the existing SFV theory with a priori convergence results for statistical quantities, in particular push-forward densities, which we demonstrate through numerical experiments. By linking refinement indicators to regions of the physical and stochastic spaces, we drive anisotropic refinements of the discretizations, introducing new degrees of freedom where deemed profitable. To illustrate our proposed method, we consider a series of numerical examples for nonlinear hyperbolic PDEs based on Burgers’ and Euler’s equations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Phase separation, edge currents, and Hall effect for active matter with Magnus dynamics

Here we examine run-and-tumble disks in two-dimensional systems where the particles also have a Magnus component to their dynamics. For increased activity, we find that the system forms a motility-induced phase-separated (MIPS) state with chiral edge flow around the clusters, where the direction of the current is correlated with the sign of the Magnus term. The stability of the MIPS state is non-monotonic as a function of increasing Magnus term amplitude, with the MIPS region first extending down to lower activities followed by a break up of MIPS at large Magnus amplitudes into a gel-like state. We examine the dynamics in the presence of quenched disorder and a uniform drive and find that the bulk flow exhibits a drive-dependent Hall angle. This is a result of the side jump effect produced by scattering from the pinning sites and is similar to the behavior found for skyrmions in chiral magnets with quenched disorder.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Superposed hyperbolic kink and pulse solutions of coupled $φ$ 4 , NLS and mKdV equations

Here, in this paper, we obtain novel solutions of a coupled $φ$ 4 , a coupled nonlinear Schrödinger equation and a coupled modified Korteweg de Vries equation which can be re-expressed as a linear superposition of either the sum or the difference of two hyperbolic pulse solutions or the sum of either a two-kink or a kink and an antikink solution. These results demonstrate that the notion of superposed solutions extends to coupled nonlinear equations as well.

97 MATHEMATICS AND COMPUTING↗