Search NASA⌕ Search

SEARCH · Search NASA

Results for “explicit solution”

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 1,261 records · Page 70
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Variational control forces for enhanced sampling of nonequilibrium molecular dynamics simulations

We introduce a variational algorithm to estimate the likelihood of a rare event within a nonequilibrium molecular dynamics simulation through the evaluation of an optimal control force. Optimization of a control force within a chosen basis is made possible by explicit forms for the gradients of a cost function in terms of the susceptibility of driven trajectories to changes in variational parameters. We consider probabilities of time-integrated dynamical observables as characterized by their large deviation functions and find that in many cases, the variational estimate is quantitatively accurate. Additionally, we provide expressions to exactly correct the variational estimate that can be evaluated directly. We benchmark this algorithm against the numerically exact solution of a model of a driven particle in a periodic potential, where the control force can be represented with a complete basis. Additionally, we then demonstrate the utility of the algorithm in a model of repulsive particles on a line, which undergo a dynamical phase transition, resulting in singular changes to the form of the optimal control force. In both systems, we find fast convergence and are able to evaluate large deviation functions with significant increases in statistical efficiency over alternative Monte Carlo approaches.

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

A Refinement-by-Superposition -Method for (curl)- and (div)-Conforming Discretizations

Here, we present refinement-by-superposition (RBS) hp-refinement infrastructure for computational electromagnetics (CEMs), which permits exponential rates of convergence. In contrast to dominant approaches to hp-refinement for continuous Galerkin methods, which rely on explicit constraint equations, the multilevel strategy presented drastically reduces the implementation complexity. Through the RBS methodology, enforcement of continuity occurs by construction, enabling arbitrary levels of refinement with ease, and without the practical (but not theoretical) limitations of constrained-node refinement. We outline the construction of the RBS hp-method for refinement with H (curl)- and H (div)-conforming finite cells. Numerical simulations for the 2-D finite element method (FEM) solution of the Maxwell eigenvalue problem demonstrate the effectiveness of RBS hp-refinement. As an additional goal of this work, we aim to promote the use of mixed-order (low- and high-order) elements in practical CEM applications.

42 ENGINEERING↗
NASA NTRS1991

Accuracy of an unstructured-grid upwind-Euler algorithm for the ONERA M6 wing

Improved algorithms for the solution of the three-dimensional, time-dependent Euler equations are presented for aerodynamic analysis involving unstructured dynamic meshes. The improvements have been developed recently to the spatial and temporal discretizations used by unstructured-grid flow solvers. The spatial discretization involves a flux-split approach that is naturally dissipative and captures shock waves sharply with at most one grid point within the shock structure. The temporal discretization involves either an explicit time-integration scheme using a multistage Runge-Kutta procedure or an implicit time-integration scheme using a Gauss-Seidel relaxation procedure, which is computationally efficient for either steady or unsteady flow problems. With the implicit Gauss-Seidel procedure, very large time steps may be used for rapid convergence to steady state, and the step size for unsteady cases may be selected for temporal accuracy rather than for numerical stability. Steady flow results are presented for both the NACA 0012 airfoil and the Office National d'Etudes et de Recherches Aerospatiales M6 wing to demonstrate applications of the new Euler solvers. The paper presents a description of the Euler solvers along with results and comparisons that assess the capability.

Batina, John T.↗
DOE OSTI2022

An Efficient Bayesian Approach to Learning Droplet Collision Kernels: Proof of Concept Using “Cloudy,” a New n -Moment Bulk Microphysics Scheme

The small-scale microphysical processes governing the formation of precipitation particles cannot be resolved explicitly by cloud resolving and climate models. Instead, they are represented by microphysics schemes that are based on a combination of theoretical knowledge, statistical assumptions, and fitting to data (“tuning”). Historically, tuning was done in an ad hoc fashion, leading to parameter choices that are not explainable or repeatable. Recent work has treated it as an inverse problem that can be solved by Bayesian inference. The posterior distribution of the parameters given the data—the solution of Bayesian inference—is found through computationally expensive sampling methods, which require over $\mathcal{O}$(10 5 ) evaluations of the forward model; this is prohibitive for many models. We present a proof of concept of Bayesian learning applied to a new bulk microphysics scheme named “Cloudy,” using the recently developed Calibrate-Emulate-Sample (CES) algorithm. Cloudy models collision-coalescence and collisional breakup of cloud droplets with an adjustable number of prognostic moments and with easily modifiable assumptions for the cloud droplet mass distribution and the collision kernel. The CES algorithm uses machine learning tools to accelerate Bayesian inference by reducing the number of forward evaluations needed to $\mathcal{O}$(10 2 ). It also exhibits a smoothing effect when forward evaluations are polluted by noise. In a suite of perfect-model experiments, we show that CES enables computationally efficient Bayesian inference of parameters in Cloudy from noisy observations of moments of the droplet mass distribution. In an additional imperfect-model experiment, a collision kernel parameter is successfully learned from output generated by a Lagrangian particle-based microphysics model.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1983

A function space approach to state and model error estimation for elliptic systems

An approach is advanced for the concurrent estimation of the state and of the model errors of a system described by elliptic equations. The estimates are obtained by a deterministic least-squares approach that seeks to minimize a quadratic functional of the model errors, or equivalently, to find the vector of smallest norm subject to linear constraints in a suitably defined function space. The minimum norm solution can be obtained by solving either a Fredholm integral equation of the second kind for the case with continuously distributed data or a related matrix equation for the problem with discretely located measurements. Solution of either one of these equations is obtained in a batch-processing mode in which all of the data is processed simultaneously or, in certain restricted geometries, in a spatially scanning mode in which the data is processed recursively. After the methods for computation of the optimal estimates are developed, an analysis of the second-order statistics of the estimates and of the corresponding estimation error is conducted. Based on this analysis, explicit expressions for the mean-square estimation error associated with both the state and model error estimates are then developed.

Rodriguez, G.↗
NASA NTRS1983

A function space approach to state and model error estimation for elliptic systems

An approach is advanced for the concurrent estimation of the state and of the model errors of a system described by elliptic equations. The estimates are obtained by a deterministic least-squares approach that seeks to minimize a quadratic functional of the model errors, or equivalently, to find the vector of smallest norm subject to linear constraints in a suitably defined function space. The minimum norm solution can be obtained by solving either a Fredholm integral equation of the second kind for the case with continuously distributed data or a related matrix equation for the problem with discretely located measurements. Solution of either one of these equations is obtained in a batch-processing mode in which all of the data is processed simultaneously or, in certain restricted geometries, in a spatially scanning mode in which the data is processed recursively. After the methods for computation of the optimal esimates are developed, an analysis of the second-order statistics of the estimates and of the corresponding estimation error is conducted. Based on this analysis, explicit expressions for the mean-square estimation error associated with both the state and model error estimates are then developed. While this paper focuses on theoretical developments, applications arising in the area of large structure static shape determination are contained in a closely related paper (Rodriguez and Scheid, 1982).

Rodriguez, G.↗
NASA NTRS1993

The structure and stability of transonic accretion disks surrounding black holes

Stationary transonic alpha-viscosity models of accretion disks surrounding nonrotating black holes have been investigated. The viscosity is modified such that it vanishes in the supersonic region to ensure its effect does not violate the causality condition. In contrast to previous studies, the viscous stress is taken to be explicitly proportional to the angular velocity gradient and is not assumed to depend solely on the local pressure in the disk. The numerical results reveal that the structure of the innermost regions of the disk are more sensitive to the modified form of the viscosity than to the form of the viscous stress. The critical sonic point is located inside the innermost stable circular orbit of a test particle at 3 Schwarzschild radii. In these solutions, the transition from subsonic to supersonic flow results from pressure effects and not viscous effects. The linear stability of these disks has been examined in the local approximation. It is found that radiative energy transport and viscous stresses in the radial direction can have important effects. As a result, it is shown that the growth rate of the inertial-acoustic mode reaches a maximum at a critical wavelength.

Chen, Xingming↗
NASA NTRS1993

Energy considerations in computational aeroacoustics

A finite-volume multistage time-stepping Euler code is used to investigate the use of CFD algorithms for the direct calculation of acoustics. The 2D compressible inviscid flow about an accelerating or decelerating circular cylinder is used as a model problem. The time evolution of the energy transfer from the cylinder to the fluid, as the cylinder is moved from rest to some nonnegligible velocity, is clearly seen. By examining the temporal and spatial characteristics of the numerical solution, a distinction can be made between the propagating acoustic energy, the convecting energy associated with the entropy change in the fluid, and the energy contained in the local aerodynamic field. Systematic variation of the cylinder acceleration shows that the radiated acoustic energy depends strongly upon the rate of acceleration or deceleration. The computational grid has a large effect on the ratio of acoustic energy to nonphysical entropy associated energy, while the role of the explicit artificial viscosity seems to be of second order. The entropy term was nearly negligible in all cases the cylinder was started slowly.

Brentner, Kenneth S.↗
DOE OSTI2024

Investigation of Ring-heterogeneous Geometry Approximation for Efficient and Practical SFR Analysis

A novel geometry approximation scheme named ring-heterogeneous approximation for intermediate resolution multiphysics analysis of SFRs is presented. This approach particularly targets neutronics coupled with thermo-mechanics feedback, which dominate the passive safety of SFR designs. The ring-heterogeneous approximation casts an array of pins into volume-preserving bands of materials that constitute the pins. It allows to efficiently simulate SFRs with an unstructured mesh and also to explicitly consider the deformation of each material by retaining their heterogeneity. Verification results with single assembly and core problems of the ABTR benchmark demonstrated that the ring-heterogeneous approximation is a reasonable approximation for SFRs. For a 2D full core problem, the ring-heterogeneous approximation significantly reduced the number of elements by at least a factor of 70 from a fully-heterogeneous mesh while almost exactly preserving the reference fully-heterogeneous solutions with the RMS error of only 0.02% in the ring-wise fission reaction rate.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS1979

Boundary conditions for implicit solutions to the compressible Navier-Stokes equations in finite computational domains

The admissible boundary conditions for the unsteady Navier-Stokes equations are insufficient to determine completely the flow variables at grid points on the boundaries of a finite computational domain. A general method is presented in which the boundary conditions are supplemented by conservatively differenced forms of the flow equations. For implicit algorithms, the method provides a fully implicit coupling between boundary and interior points. Implicit boundary-point computation schemes are given for most types of flow boundary encountered in internal and external flow problems. Convergence to steady state is generally faster than with locally explicit boundary-point schemes.

Thomas, P. D.↗
NASA NTRS1993

Development of an explicit multiblock/multigrid flow solver for viscous flows in complex geometries

A new computer program is being developed for doing accurate simulations of compressible viscous flows in complex geometries. The code employs the full compressible Navier-Stokes equations. The eddy viscosity model of Baldwin and Lomax is used to model the effects of turbulence on the flow. A cell centered finite volume discretization is used for all terms in the governing equations. The Advection Upwind Splitting Method (AUSM) is used to compute the inviscid fluxes, while central differencing is used for the diffusive fluxes. A four-stage Runge-Kutta time integration scheme is used to march solutions to steady state, while convergence is enhanced by a multigrid scheme, local time-stepping and implicit residual smoothing. To enable simulations of flows in complex geometries, the code uses composite structured grid systems where all grid lines are continuous at block boundaries (multiblock grids). Example results are shown a flow in a linear cascade, a flow around a circular pin extending between the main walls in a high aspect-ratio channel, and a flow of air in a radial turbine coolant passage.

Steinthorsson, E.↗
NASA NTRS1993

Development of an explicit multiblock/multigrid flow solver for viscous flows in complex geometries

A new computer program is being developed for doing accurate simulations of compressible viscous flows in complex geometries. The code employs the full compressible Navier-Stokes equations. The eddy viscosity model of Baldwin and Lomax is used to model the effects of turbulence on the flow. A cell centered finite volume discretization is used for all terms in the governing equations. The Advection Upwind Splitting Method (AUSM) is used to compute the inviscid fluxes, while central differencing is used for the diffusive fluxes. A four-stage Runge-Kutta time integration scheme is used to march solutions to steady state, while convergence is enhanced by a multigrid scheme, local time-stepping, and implicit residual smoothing. To enable simulations of flows in complex geometries, the code uses composite structured grid systems where all grid lines are continuous at block boundaries (multiblock grids). Example results shown are a flow in a linear cascade, a flow around a circular pin extending between the main walls in a high aspect-ratio channel, and a flow of air in a radial turbine coolant passage.

Steinthorsson, E.↗
DOE OSTI2021

Scale-bridging with the extended/generalized finite element method for linear elastodynamics

This paper presents an extended/generalized finite element method for bridging scales in linear elastodynamics in the absence of scale separation. More precisely, the GFEMgl framework is expanded to enable the numerical solution of multiscale problems through the automated construction of specially-tailored shape functions, thereby enabling high-fidelity finite element modeling on simple, fixed finite element meshes. Furthermore, this introduces time-dependencies in the shape functions in that they are subject to continuous adaptation with time. The temporal aspects of the formulation are investigated by considering the Newmark-β time integration scheme, and the efficacy of mass lumping strategies is explored in an explicit time-stepping scheme. This method is demonstrated on representative wave propagation examples as well as a dynamic fracture problem to assess its accuracy and flexibility.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Efficient Measurement of Length Distribution of 1D Nanoparticles in Solution via Optical Polarimetry

The efficient measurement of the length distribution of nanotubes, nanowires, and other one-dimensional (1D) nanoparticles in solution is important to enable their incorporation into materials and devices and to optimize their processing for properties of interest, such as thermal/electrical conductivity or mechanical strength, in suspensions and composites. We report an electric-field (E-field)-assisted optical-polarimetry technique to measure the length distribution of ensembles of high-aspect-ratio particles in dilute suspension. The degree of alignment of polarizable 1D particles suspended in a fluid under Brownian motion explicitly depends on the E-field strength and the particle length. We show that it is possible to extract the length distribution of 1D nanoparticles suspended in an insulating fluid by applying a range of E-fields and using optical polarimetry to measure the corresponding alignment order parameter. Notably, the method is relatively insensitive to the diameter of the 1D particles, which can be poorly known or vary within a sample. The technique is validated with silver nanowires and carbon nanotubes of known lengths, as well as polymer-depletion-length-separated single-wall carbon nanotube samples with length distributions independently measured with analytical ultracentrifugation. Finally, we demonstrate the ability of the optical-polarimetry technique to quantify changes in the length distribution of ultranarrow, sub-nanometer-diameter single-wall carbon nanotubes under different types and durations of ultrasonication. Within its range of applicability (polarizable 1D nanoparticles in the 0.5 to 15 μm length range, constrained by the voltage stability of the media and the suspended particles), the E-field-assisted optical-polarimetry method is a particularly efficient and accurate method to measure the length distribution of nanowires and nanotubes in suspension.

1D nanoparticles↗
DOE OSTI2024

Maven: a multimodal foundation model for supernova science

Abstract A common setting in astronomy is the availability of a small number of high-quality observations, and larger amounts of either lower-quality observations or synthetic data from simplified models. Time-domain astrophysics is a canonical example of this imbalance, with the number of supernovae observed photometrically outpacing the number observed spectroscopically by multiple orders of magnitude. At the same time, no data-driven models exist to understand these photometric and spectroscopic observables in a common context. Contrastive learning objectives, which have grown in popularity for aligning distinct data modalities in a shared embedding space, provide a potential solution to extract information from these modalities. We present Maven, the first foundation model for supernova science. To construct Maven, we first pre-train our model to align photometry and spectroscopy from 0.5 M synthetic supernovae using a contrastive objective. We then fine-tune the model on 4702 observed supernovae from the Zwicky transient facility. Maven reaches state-of-the-art performance on both classification and redshift estimation, despite the embeddings not being explicitly optimized for these tasks. Through ablation studies, we show that pre-training with synthetic data improves overall performance. In the upcoming era of the Vera C. Rubin observatory, Maven will serve as a valuable tool for leveraging large, unlabeled and multimodal time-domain datasets.

Zhang, Gemma (ORCID:0000000280198082)↗
DOE OSTI2022

Doping-driven resistive collapse of the Mott insulator in a minimal model for VO 2

Here we study the resistive collapse of the Mott insulator state in the dimer Hubbard model. This minimal model has been used to describe the physics of VO 2 , and should be relevant to other strongly correlated materials. It incorporates the physics of correlated dimers and the explicit competition between on-site Coulomb repulsion and magnetic exchange interactions. Our results unveil that between the Mott insulator at half filling and the Fermi liquid metal at high doping there is an intermediate bad metallic phase with exotic features such as a pseudogap, orbital selectivity, and a first-order metal-metal transition. The model is solved within dynamical mean field theory by means of quantum Monte Carlo, which provides the numerically exact solution of the model in the limit of large lattice dimensionality. This model can be considered as a minimal one that captures exotic phenomena associated to the physics of a doped Mott insulator, shading light on their basic physical mechanism.

36 MATERIALS SCIENCE↗
NASA NTRS1981

Transonic swept wings studied by the lifting-line theory

Transonic swept wings are analyzed as a lifting-line problem under a small-disturbance approximation. Basic concepts and principal results of the asymptotic theory are discussed. The study focuses on straight oblique wings and V-shaped swept wings, of which the local centerline curvature can be equated to zero. The three-dimensional (3-D) perturbation of the nonlinear component flow admits a similarity flow structure but requires that all wing sections are generated from a single airfoil profile; the reduced 2-D problems in this case are solved only once for all span stations. Examples of solutions involving high subcritical and slightly supercritical component flows are demonstrated and compared with surface pressure data from 3-D computer codes based on the full-potential equation (FLO 22). Except in the neighborhood of leading edges, where the small-disturbance assumption breaks down, and in the vicinities of wing tips and the symmetry plane, where neither the theory nor the 3-D codes may claim full validity, reasonable agreement is consistently found. The explicit results from the upwash analysis, along with the similarity flow structure, provides a rational approach to the control of 3-D effects in transonic aerodynamic design studies.

Cheng, H. K.↗
NASA NTRS1983

Explicit finite difference predictor and convex corrector with applications to hyperbolic partial differential equations

An explicit finite difference scheme consisting of a predictor and a corrector has been developed and applied to solve some hyperbolic partial differential equations (PDEs). The corrector is a convex-type function which is applied at each time level and at each mesh point. It consists of a parameter which may be estimated such that for larger time steps the algorithm should remain stable and generate a fast speed of convergence to the steady-state solution. Some examples have been given.

Dey, C.↗