Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Weak form”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Electron scattering and neutrino physics

Abstract A thorough understanding of neutrino–nucleus scattering physics is crucial for the successful execution of the entire US neutrino physics program. Neutrino–nucleus interaction constitutes one of the biggest systematic uncertainties in neutrino experiments—both at intermediate energies affecting long-baseline deep underground neutrino experiment, as well as at low energies affecting coherent scattering neutrino program—and could well be the difference between achieving or missing discovery level precision. To this end, electron–nucleus scattering experiments provide vital information to test, assess and validate different nuclear models and event generators intended to test, assess and validate different nuclear models and event generators intended to be used in neutrino experiments. Similarly, for the low-energy neutrino program revolving around the coherent elastic neutrino–nucleus scattering (CEvNS) physics at stopped pion sources, such as at ORNL, the main source of uncertainty in the evaluation of the CEvNS cross section is driven by the underlying nuclear structure, embedded in the weak form factor, of the target nucleus. To this end, parity-violating electron scattering (PVES) experiments, utilizing polarized electron beams, provide vital model-independent information in determining weak form factors. This information is vital in achieving a percent level precision needed to disentangle new physics signals from the standard model expected CEvNS rate. In this white paper, we highlight connections between electron- and neutrino–nucleus scattering physics at energies ranging from 10 s of MeV to a few GeV, review the status of ongoing and planned electron scattering experiments, identify gaps, and lay out a path forward that benefits the neutrino community. We also highlight the systemic challenges with respect to the divide between the nuclear and high-energy physics communities and funding that presents additional hurdles in mobilizing these connections to the benefit of neutrino programs.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

A unified framework of stabilized finite element method for solving the Boltzmann transport equation

This paper presents a unified framework of stabilized finite element method for solving the Boltzmann transport equation. Unlike the traditional Petrov-Galerkin finite element method which modifies the test function to construct the stabilization term, we derive the stabilization methods from the standard Galerkin weak form with Sub-grid scale model. The basic idea of this method is to decompose the unknowns into its numerical solution and residual, with an approximation for the residual and embeds it in the Galerkin weak form to yield a stabilized variational formula. Different approximations of the residual lead to different stabilization methods, all the frequently used stabilized methods, including the Streamline Upwinding Petrov-Galerkin (SUPG) method, Galerkin/Least-Square (GLS) method, and Algebraic Sub-Grid Scale (ASGS) method can be obtained from this framework. The similarities and differences of the different approximations are compared in this paper. The numerical results show that the behaviors of the different methods area similar with the same stabilization parameter, and all these stabilized techniques can obtain a correct and stable solution. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS1997

Numerical Computation of Sensitivities and the Adjoint Approach

We discuss the numerical computation of sensitivities via the adjoint approach in optimization problems governed by differential equations. We focus on the adjoint problem in its weak form. We show how one can avoid some of the problems with the adjoint approach, such as deriving suitable boundary conditions for the adjoint equation. We discuss the convergence of numerical approximations of the costate computed via the weak form of the adjoint problem and show the significance for the discrete adjoint problem.

Lewis, Robert Michael↗
NASA NTRS2002

Development of MLPG and LBIE Methods for Nonlinear Problems of Fracture

The truly meshless local Petrov-Galerkin (MLPG) method holds a great promise in solving boundary value problems, using a local symmetric weak form as a natural approach. In the present paper, in the context of MLPG and the meshless interpolation of a moving least squares (MLS) type, a method which uses primary and secondary nodes in the domain and on the global boundary is introduced, in order to improve the accuracy of solution. The secondary nodes can be placed at any location where one needs to obtain a better resolution. The sub-domains for the shape functions in the MLS approximation are defined only from the primary nodes, and the secondary nodes use the same sub-domains. The shape functions based on the MLS approximation, in an integration domain, have a single type of a rational function, which reduces the difficulty of numerical integration to evaluate the weak form. The present method is very useful in an adaptive calculation, because the secondary nodes can be easily added and/or moved without an additional mesh. The essential boundary conditions can be imposed exactly, and non-convex boundaries can be treated without special techniques. Several numerical examples are presented to illustrate the performance of the present method.

Atluri, Satya N.↗
DOE OSTI2024

Coarse-graining Hamiltonian systems using WSINDy

Abstract Weak form equation learning and surrogate modeling has proven to be computationally efficient and robust to measurement noise in a wide range of applications including ODE, PDE, and SDE discovery, as well as in coarse-graining applications, such as homogenization and mean-field descriptions of interacting particle systems. In this work we extend this coarse-graining capability to the setting of Hamiltonian dynamics which possess approximate symmetries associated with timescale separation. A smooth $$\varepsilon$$ ε -dependent Hamiltonian vector field $$X_\varepsilon$$ X ε possesses an approximate symmetry if the limiting vector field $$X_0=\lim _{\varepsilon \rightarrow 0}X_\varepsilon$$ X 0 = lim ε → 0 X ε possesses an exact symmetry. Such approximate symmetries often lead to the existence of a Hamiltonian system of reduced dimension that may be used to efficiently capture the dynamics of the symmetry-invariant dependent variables. Deriving such reduced systems, or approximating them numerically, is an ongoing challenge. We demonstrate that WSINDy can successfully identify this reduced Hamiltonian system in the presence of large perturbations imparted in the $$\varepsilon >0$$ ε > 0 regime, while remaining robust to extrinsic noise. This is significant in part due to the nontrivial means by which such systems are derived analytically. WSINDy naturally preserves the Hamiltonian structure by restricting to a trial basis of Hamiltonian vector fields. The methodology is computationally efficient, often requiring only a single trajectory to learn the global reduced Hamiltonian, and avoiding forward solves in the learning process. In this way, we argue that weak-form equation learning is particularly well-suited for Hamiltonian coarse-graining. Using nearly-periodic Hamiltonian systems as a prototypical class of systems with approximate symmetries, we show that WSINDy robustly identifies the correct leading-order system, with dimension reduced by at least two, upon observation of the relevant degrees of freedom. While our main contribution is computational, we also provide a contribution to the literature on averaging theory by proving that first-order averaging at the level of vector fields preserves Hamiltonian structure in nearly-periodic Hamiltonian systems. This provides theoretical justification for our approach as WSINDy’s computations occur at the level of Hamiltonian vector fields. We illustrate the efficacy of our proposed method using physically relevant examples, including coupled oscillator dynamics, the Hénon–Heiles system for stellar motion within a galaxy, and the dynamics of charged particles.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Weak-charge form-factor determination at the electron-ion collider

Determining the weak charge form factor, 𝐹 𝑊 ⁡(𝑄 2 ), of nuclei over a continuous range of momentum transfers, 0 ≲ 𝑄 2 ≲ 0.1 GeV 2 , is essential for mapping out the distribution of neutrons in nuclei. The neutron density distribution has significant implications for a broad range of areas, including studies of nuclear structure, neutron stars, and physics beyond the Standard Model. Currently, our knowledge of 𝐹 𝑊 ⁡(𝑄 2 ) comes primarily from fixed target experiments that measure the parity-violating asymmetry in coherent elastic electron-ion scattering. Fixed target experiments, such as CREX and PREX-1,2, have provided high-precision weak charge form factor extractions for the 48 Ca and 208 Pb nuclei, respectively. However, a major limitation of fixed target experiments is that they each provide data only at a single value of 𝑄 2 . With the proposed electron-ion collider (EIC) on the horizon, we explore its potential to impact the determination of the weak charge form factor. While it cannot compete with the precision of fixed target experiments, it can provide data over a wide and continuous range of 𝑄 2 values, and for a wide variety of nuclei. We show that with data corresponding to an integrated luminosity of ℒ ∼ 500/𝐴 fb −1 , where 𝐴 is the nucleus atomic weight, the EIC can significantly impact constraints by lifting degeneracies in theoretical models of the neutron density distribution. Ensuring EIC detector coverage at low 𝑄 2 and large negative pseudorapidities will be essential for such 𝐹 𝑊 ⁡(𝑄 2 ) measurements.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Generalized Boozer coordinates: A natural coordinate system for quasisymmetry

We prove the existence of a straight-field-line coordinate system we call generalized Boozer coordinates. This coordinate system exists for magnetic fields with nested toroidal flux surfaces provided ∮ d l / B   ( j · ∇ ψ ) = 0 , where symbols have their usual meaning, and the integral is taken along closed magnetic field lines. All quasisymmetric fields, regardless of their associated form of equilibria, must satisfy this condition. This coordinate system presents itself as a convenient form to describe general quasisymmetric configurations and their properties. Insight can be gained analytically into the difference between strong and weak forms of quasisymmetry, as well as axisymmetry, and the interaction of quasisymmetry with different forms of equilibria.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Formation of Water Networks on Anionic Perylene

We present infrared photodissociation spectra of hydrated perylene anion clusters with up to four water molecules, as well as electronic structure calculations based on density functional theory. Water molecules form weak hydrogen bonds to the π system of the perylene anion. For clusters with more than one water molecule, water−water hydrogen bonds are formed, which generally appear to be stronger than water−π hydrogen bonds, especially for the trihydrate and tetrahydrate. The resulting water networks exist as water subclusters on the surface of the carbon frame of perylene. Here, we observe temperature-dependent dynamic effects, which highlight large amplitude motions of the water network and the shallowness of the potential energy surfaces governing the structures of these clusters.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1993

New results on the realizability of Reynolds stress turbulence closures

The realizability of Reynolds stress models in homogeneous turbulence is critically assessed from a theoretical standpoint. It is proven that a well known second-order closure formulated by Shih and Lumley using the strong realizability constraints of Schumann is, in fact, not a realizable model. The problem arises from the failure to properly satisfy the necessary positive second time derivative constraint when a principal Reynolds stress vanishes - a fatal flaw that becomes apparent when the non-analytic terms in their model are made single-valued as required on physical grounds. It is furthermore shown that the centrifugal acceleration generated by rotations of the principal axes of the Reynolds stress tensor can make the second derivative singular at the most extreme limits of realizable turbulence. This previously overlooked effect appears to make it impossible to identically satisfy the strong form of realizability in any version of the present generation of second-order closures. On the other hand, models properly formulated to satisfy the weak form of realizability - wherein states of one or two component turbulence are not accessible in finite time are found to be realizable. However, unlike the simpler and more commonly used second order closures, these models can be ill-behaved near the extreme limits of realizable turbulence due to the way that higher-degree nonlinearities are often unnecessarily introduced to satisfy realizability. Illustrative computations of homogeneous shear flows are presented to demonstrate these points which can have important implications for turbulence modeling.

Speziale, Charles G.↗
NASA NTRS1989

Tau approximation techniques for identification of coefficients in parabolic PDE

A variant of the Tau method, called the weak Tau method, is developed on the basis of the weak form of the PDE for use in least-squares parameter estimation; also presented is a suitable abstract convergence framework. The emphasis is on the theoretical framework that allows treatment of the weak Tau method when it is applied to a wide class of inverse problems, including those for diffusion-advection equations, the Fokker-Planck model for population dynamics, and damped beam equations. Extensive numerical testing of the weak Tau method has demonstrated that it compares quite favorably with existing methods.

Banks, H. T.↗
DOE OSTI2023

Precise Measurement of the Neutron Skin Thicknesses of 208Pb and 48Ca

Despite great leaps of development in nuclear physics over the past century, we still do not have a comprehensive understanding of nuclear structure. This is mainly due to the lack of precise knowledge of the neutron distribution inside nucleus. Electromagnetic probes are ineffective in probing neutrons, unlike their charged partners, the protons. Therefore, it is of great importance to constrain the neutron distributions experimentally. Heavy nuclei have more neutrons than protons in order to balance the repulsion between protons. In such neutron-rich nuclei, the extra neutrons are pushed out to the surface by the nuclear symmetry energy, forming the so-called ?neutron skin?. The neutron skin can be probed with a well-established experimental technique ? parity-violating electron scattering (PVES). Using the scattering of longitudinally polarized electrons by an unpolarized target, PREX-II and CREX measure the small parity-violating asymmetry in cross sections. Employing electrons with opposite helicities the weak form factor, the neutron distribution and the neutron skin thickness of the target nucleus are extracted. With excellent beam qualities and dedicated instrumentation at Jefferson Lab, the asymmetry measurements are statistics-limited. We report the results of these two high-precision measurements and their implications on broad topics, from the nuclear structures to the neutron stars.

Zhang, Weibin↗
NASA NTRS2020

Model-Reference Adaptive Control of Distributed Lagrangian Infinite-Dimensional Systems Using Hamilton’s Principle

This paper presents a Hamilton's principle for distributed control of infinite-dimensional systems modeled by a distributed form of the Euler-Lagrange method. The distributed systems are governed by a system of linear partial differential equations in space and time. A generalized potential energy expression is developed that can capture most physical systems including those systems that have no spatial distribution. The Hamilton's principle is applied to derive distributed feedback control methods without resorting to the standard weak-form discretization approach to convert an infinite-dimensional systems to a finite-dimensional systems. It can be shown by the principle of least action that the distributed control synthesized by the Hamilton's principle is a minimum-norm control. A model-reference adaptive control framework is developed for distributed Lagrangian systems in the presence of uncertainty. The theory is demonstrated by an application of adaptive flutter suppression control of a flexible aircraft wing.

Nguyen, Nhan T.↗
DOE OSTI2025

A higher-order finite-element implementation of the nonlinear Fokker–Planck collision operator for charged particle collisions in a low density plasma

Collisions between particles in a low density plasma are described by the Fokker–Planck collision operator. In applications, this nonlinear integro-differential operator is often approximated by linearised or ad-hoc model operators due to computational cost and complexity. In this work, we present an implementation of the nonlinear Fokker–Planck collision operator written in terms of Rosenbluth potentials in the Rosenbluth–MacDonald–Judd (RMJ) form. The Rosenbluth potentials may be obtained either by direct integration or by solving partial differential equations (PDEs) similar to Poisson's equation: we optimise for performance and scalability by using sparse matrices to solve the relevant PDEs. We represent the distribution function using a tensor-product continuous-Galerkin finite-element representation and we derive and describe the implementation of the weak form of the collision operator. We present tests demonstrating a successful implementation using an explicit time integrator and we comment on the speed and accuracy of the operator. Finally, we speculate on the potential for applications in the current and next generation of kinetic plasma models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Numerical integration in the virtual element method with the scaled boundary cubature scheme

Abstract The virtual element method (VEM) is a stabilized Galerkin method on meshes that consist of arbitrary (convex and nonconvex) polygonal and polyhedral elements. A crucial ingredient in the implementation of low‐ and high‐order VEM is the numerical integration of monomials and nonpolynomial functions over such elements. In this article, we apply the recently proposed scaled boundary cubature (SBC) scheme to compute the weak form integrals in various virtual element formulations over polygonal and polyhedral meshes. In doing so, we demonstrate the flexibility of the approach and the accuracy that it delivers on a broad suite of boundary‐value problems in 2D and 3D over polytopes with affine faces as well as on elements with curved boundaries. In addition, the use of the SBC scheme is exemplified in an enriched Poisson formulation of the VEM in which weakly singular functions are required to be integrated. This study establishes the SBC method as a simple, accurate and efficient integration scheme for use in the VEM.

Chin, Eric B.↗
DOE OSTI2025

Inference of phase field fracture models

The phase field approach to modeling fracture uses a diffuse damage field to represent cracks. This representation mollifies singularities that arise in computations with sharp interface models and some of the resultant difficulties in the mathematical and numerical treatment of fracture. Phase field fracture models have proven effective at representing crack propagation, branching, and merging. Specific formulations, beginning with brittle fracture, have also been shown to converge to classical solutions. Extensions to cover the range of material failure, including ductile and cohesive fracture, lead to an array of possible models. There exists a large body of literature focusing on this class of models and on the impact of model form on the predicted crack evolution. However, there have not been systematic studies into how optimal models may be chosen. Here, we take a first step in this direction by developing formal methods for identification of the best parsimonious model of phase field fracture given full-field data on the damage and deformation fields. We consider some of the main models that have been used for the degradation of elastic response due to damage and its propagation. Our approach builds upon Variational System Identification (VSI), a weak form variant of the Sparse Identification of Nonlinear Dynamics (SINDy). Furthermore, in this first communication we focus on synthetically generated data but we also consider central issues associated with the use of experimental full-field data, such as data sparsity and noise.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

A conservative Galerkin solver for the quasilinear diffusion model in magnetized plasmas

We propose a conservative Galerkin scheme for the quasilinear model in three-dimensional momentum space and three-dimensional spectral space, with cylindrical symmetry. We construct an unconditionally conservative weak form and use a discretization that preserves conservation properties independent of the wave emission probability. The discrete operators, combined with a consistent quadrature rule, preserve all the conservation laws rigorously. The proposed scheme is quite general: it works for both relativistic and non-relativistic systems, for both magnetized and unmagnetized plasmas, and even for problems with time-dependent dispersion relations. We represent the particle distribution by continuous basis functions and use discontinuous basis functions for the wave spectral energy density, which enables the application of a positivity-preserving technique. We adopt the marching simplex algorithm, designed initially for computer graphics, for numerical integration on the resonance manifold. Furthermore, the numerical examples with a bump-on-tail initial configuration show how the unstable waves produce strong momentum space diffusion.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

A Particle Method for the Multispecies Landau Equation

Abstract The multispecies Landau collision operator describes the two-particle, small scattering angle or grazing collisions in a plasma made up of different species of particles such as electrons and ions. Recently, a structure preserving deterministic particle method (Carrillo et al. in J. Comput. Phys. 7:100066, 2020) has been developed for the single species spatially homogeneous Landau equation. This method relies on a regularization of the Landau collision operator so that an approximate solution, which is a linear combination of Dirac delta distributions, is well-defined. Based on a weak form of the regularized Landau equation, the time dependent locations of the Dirac delta functions satisfy a system of ordinary differential equations. In this work, we extend this particle method to the multispecies case, and examine its conservation of mass, momentum, and energy, and decay of entropy properties. We show that the equilibrium distribution of the regularized multispecies Landau equation is a Maxwellian distribution, and state a critical condition on the regularization parameters that guarantees a species independent equilibrium temperature. A convergence study comparing an exact multispecies Bobylev-Krook-Wu (BKW) solution to the particle solution shows approximately 2nd order accuracy. Important physical properties such as conservation, decay of entropy, and equilibrium distribution of the particle method are demonstrated with several numerical examples.

Mathematics↗