Publications

Results 1–25 of 108
Skip to search filters

A conservative, optimization-based semi-lagrangian spectral element method for passive tracer transport

COUPLED PROBLEMS 2015 - Proceedings of the 6th International Conference on Coupled Problems in Science and Engineering

Bochev, Pavel B.; Moe, Scott A.; Peterson, Kara J.; Ridzal, Denis R.

We present a new optimization-based, conservative, and quasi-monotone method for passive tracer transport. The scheme combines high-order spectral element discretization in space with semi-Lagrangian time stepping. Solution of a singly linearly constrained quadratic program with simple bounds enforces conservation and physically motivated solution bounds. The scheme can handle efficiently a large number of passive tracers because the semi-Lagrangian time stepping only needs to evolve the grid points where the primitive variables are stored and allows for larger time steps than a conventional explicit spectral element method. Numerical examples show that the use of optimization to enforce physical properties does not affect significantly the spectral accuracy for smooth solutions. Performance studies reveal the benefits of high-order approximations, including for discontinuous solutions.

More Details

A virtual control coupling approach for problems with non-coincident discrete interfaces

Lecture Notes in Computer Science (including subseries Lecture Notes in Artificial Intelligence and Lecture Notes in Bioinformatics)

Bochev, Pavel B.; Kuberry, Paul A.; Peterson, Kara J.

Independent meshing of subdomains separated by an interface can lead to spatially non-coincident discrete interfaces. We present an optimization-based coupling method for such problems, which does not require a common mesh refinement of the interface, has optimal H1 convergence rates, and passes a patch test. The method minimizes the mismatch of the state and normal stress extensions on discrete interfaces subject to the subdomain equations, while interface “fluxes” provide virtual Neumann controls.

More Details

A virtual control, mesh-free coupling method for non-coincident interfaces

Proceedings of the 6th European Conference on Computational Mechanics: Solids, Structures and Coupled Problems, ECCM 2018 and 7th European Conference on Computational Fluid Dynamics, ECFD 2018

Kuberry, Paul A.; Bochev, Pavel B.; Peterson, Kara J.

We present an optimization approach with two controls for coupling elliptic partial differential equations posed on subdomains sharing an interface that is discretized independently on each subdomain, introducing gaps and overlaps. We use two virtual Neumann controls, one defined on each discrete interface, thereby eliminating the need for a virtual common refinement interface mesh. Global flux conservation is achieved by including the square of the difference of the total flux on each interface in the objective. We use Generalized Moving Least Squares (GMLS) reconstruction to evaluate and compare the subdomain solution and gradients at quadrature points used in the cost functional. The resulting method recovers globally linear solutions and shows optimal L2-norm and H1-norm convergence.

More Details

An overview of magneto-inertial fusion on the Z Machine at Sandia National Laboratories

Yager-Elorriaga, David A.; Gomez, Matthew R.; Ruiz, Daniel E.; Slutz, Stephen A.; Harvey-Thompson, Adam J.; Jennings, Christopher A.; Knapp, Patrick K.; Schmit, Paul S.; Weis, Matthew R.; Awe, Thomas J.; Chandler, Gordon A.; Mangan, Michael M.; Myers, Clayton E.; Fein, Jeffrey R.; Geissel, Matthias G.; Glinsky, Michael E.; Hansen, Stephanie B.; Harding, Eric H.; Lamppa, Derek C.; Webster, Evelyn L.; Rambo, Patrick K.; Robertson, Grafton K.; Savage, Mark E.; Smith, Ian C.; Ampleford, David A.; Beckwith, Kristian B.; Peterson, Kara J.; Porter, John L.; Rochau, G.A.; Sinars, Daniel S.

Abstract not provided.

Arctic Climate Systems Analysis

Ivey, Mark D.; Robinson, David G.; Boslough, Mark B.; Backus, George A.; Peterson, Kara J.; van Bloemen Waanders, Bart G.; Swiler, Laura P.; Desilets, Darin M.; Reinert, Rhonda K.

This study began with a challenge from program area managers at Sandia National Laboratories to technical staff in the energy, climate, and infrastructure security areas: apply a systems-level perspective to existing science and technology program areas in order to determine technology gaps, identify new technical capabilities at Sandia that could be applied to these areas, and identify opportunities for innovation. The Arctic was selected as one of these areas for systems level analyses, and this report documents the results. In this study, an emphasis was placed on the arctic atmosphere since Sandia has been active in atmospheric research in the Arctic since 1997. This study begins with a discussion of the challenges and benefits of analyzing the Arctic as a system. It goes on to discuss current and future needs of the defense, scientific, energy, and intelligence communities for more comprehensive data products related to the Arctic; assess the current state of atmospheric measurement resources available for the Arctic; and explain how the capabilities at Sandia National Laboratories can be used to address the identified technological, data, and modeling needs of the defense, scientific, energy, and intelligence communities for Arctic support.

More Details

Arctic sea ice modeling with the material-point method

Peterson, Kara J.; Bochev, Pavel B.

Arctic sea ice plays an important role in global climate by reflecting solar radiation and insulating the ocean from the atmosphere. Due to feedback effects, the Arctic sea ice cover is changing rapidly. To accurately model this change, high-resolution calculations must incorporate: (1) annual cycle of growth and melt due to radiative forcing; (2) mechanical deformation due to surface winds, ocean currents and Coriolis forces; and (3) localized effects of leads and ridges. We have demonstrated a new mathematical algorithm for solving the sea ice governing equations using the material-point method with an elastic-decohesive constitutive model. An initial comparison with the LANL CICE code indicates that the ice edge is sharper using Materials-Point Method (MPM), but that many of the overall features are similar.

More Details
Results 1–25 of 108
Results 1–25 of 108