SEARCH · Search NASA

Results for “energy function”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 199 records · Page 11

All sources NASA NTRS DOE OSTI

DOE OSTI2019

Fluctuations and phase transitions of uniaxial and biaxial liquid crystals using a theoretically informed Monte Carlo and a Landau free energy density

In this work, we explore fluctuations during phase transitions of uniaxial and biaxial liquid crystals using a phenomenological free energy functional. We rely on a continuum-level description of the liquid crystal ordering with a tensorial parameter and a temperature dependent Landau polynomial expansion of the tensor's invariants. The free energy functional, over a three-dimensional periodic domain, is integrated with a Gaussian quadrature and minimized with a theoretically informed Monte Carlo method. We reconstruct analytical phase diagrams, following Landau and Doi 's notations, to verify that the free energy relaxation reaches the global minimum. Importantly, our relaxation method is able to follow the thermodynamic behavior provided by other non-phenomenological approaches; we predict the first order character of the isotropic-nematic transition, and we identify the uniaxial-biaxial transition as second order. Finally, we use a finite-size scaling method, using the nematic susceptibility, to calculate the transition temperatures for 4-Cyano-4'-pentylbiphenyl (5CB) and N-(4-methoxybenzylidene)-4-butylaniline (MBBA). Our results show good agreement with experimental values, thereby validating our minimization method. Our approach is an alternative towards the relaxation of temperature dependent continuum-level free energy functionals, in any geometry, and can incorporate complicated elastic and surface energy densities.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

NASA NTRS1990

A unified perspective on robot control - The energy Lyapunov function approach

A unified framework for the stability analysis of robot tracking control is presented. By using an energy-motivated Lyapunov function candidate, the closed-loop stability is shown for a large family of control laws sharing a common structure of proportional and derivative feedback and a model-based feedforward. The feedforward can be zero, partial or complete linearized dynamics, partial or complete nonlinear dynamics, or linearized or nonlinear dynamics with parameter adaptation. As result, the dichotomous approaches to the robot control problem based on the open-loop linearization and nonlinear Lyapunov analysis are both included in this treatment. Furthermore, quantitative estimates of the trade-offs between different schemes in terms of the tracking performance, steady state error, domain of convergence, realtime computation load and required a prior model information are derived.

Wen, John T.↗

DOE OSTI2014

Isospin-invariant Skyrme energy-density-functional approach with axial symmetry

Density functional theory (DFT) is the microscopic tool of choice to describe properties of nuclei over the entire nuclear landscape, with a focus on medium-mass and heavy complex systems. Modern energy density functionals (EDFs) often offer a level of accuracy typical of phenomenological approaches based on parameters locally fitted to the data. It is clear, however, that in order to achieve high quality of predictions to guide spectroscopic studies, current functionals must be improved, especially in the isospin channel. In this respect, experimental studies of short-lived nuclei far from stability offer a unique test of isospin aspects of the many-body theory. In this study, we develop the isospin-invariant Skyrme-EDF method by considering local densities in all possible isospin channels and proton-neutron ($p-n$) mixing terms as mandated by the isospin symmetry. The EDF employed has the most general form that depends quadratically on the isoscalar and isovector densities. We test and benchmark the resulting $p-n$ EDF approach, and study the general properties of the new scheme by means of the cranking in the isospin space. We extend the existing axial DFT solver HFBTHO to the case of isospin-invariant EDF approach with all possible $p-n$ mixing terms. Explicit expressions have been derived for all the densities and potentials that appear in the isospin representation. In practical tests, we consider the Skyrme EDF SkM* and, as a first application, concentrate on Hartree-Fock aspects of the problem, i.e., pairing has been disregarded. Calculations have been performed for the ($A = 78, T ≃ 11), (A = 40, T ≃ 8$), and ($A = 48, T ≃ 4$) isobaric analog chains. Isospin structure of self-consistent $p-n$ mixed solutions has been investigated with and without the Coulomb interaction, which is the sole source of isospin symmetry breaking in our approach. The extended axial HFBTHO solver has been benchmarked against the symmetry-unrestricted HFODD code for deformed and spherical states. We developed and tested a general isospin-invariant Skyrme-EDF framework. The new approach permits spin-isospin densities that may give rise to hitherto unexplored modes in the excitation spectrum. The new formalism has been tested in the Hartree-Fock limit. A systematic comparison between HFODD and HFBTHO results show a maximum deviation of about 10 keV on the total binding energy for deformed nuclei when the Coulomb term is included. Without this term, the results of both solvers agree down to a ~10 eV level.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2024

Impact of tensor couplings with scalar mixing on covariant energy density functionals

The recent pioneering campaigns conducted by the Lead Radius Experiment (PREX) and the Calcium Radius Experiment (CREX) Collaborations have uncovered major deficiencies in the theoretical description of some fundamental properties of atomic nuclei. Following a recent refinement to the isovector sector of covariant energy density functionals we present here additional improvements to the functional by including both tensor couplings and an isoscalar-isovector mixing term in the scalar sector. Motivated by the distinct surface properties of calcium and lead, we expect that the tensor terms that generate derivative couplings will help break the linear correlation between the neutron skin thickness of these two nuclei. Moreover, the addition of these new terms mitigates most of the problems identified by Reed et al. in describing the properties of both finite nuclei and neutron stars. While significant progress has been made in reconciling the PREX-CREX results without compromising other observables, the final resolution awaits the completion of a proper calibration for this new class of functionals. As a result, we expect that powerful reduced basis methods used recently to create efficient emulators will be essential to accomplish this task.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2015

A new Skyrme energy density functional for a better description of spin-isospin resonances

A correct determination of the isospin and spin-isospin properties of the nuclear effective interaction should lead to an accurate description of the Gamow-Teller resonance (GT), the Spin Dipole Resonance (SDR), the Giant Dipole Resonance (GDR) or the Antianalog Giant Dipole Resonance (AGDR), among others. A new Skyrme energy density functional named SAMi is introduced with the aim of going a step forward in setting the bases for a more precise description of spin-isospin resonances [1, 2]. In addition, we will discuss some new features of our analysis on the AGDR in {sup 208}Pb [3] as compared with available experimental data on this resonance [4, 5, 6], and on the GDR [7]. Such study, guided by a simple yet physical pocket formula, has been developed by employing the so called SAMi-J family of systematically varied interactions. This set of interactions is compatible with experimental data for values of the symmetry energy at saturation J and slope parameter L falling in the ranges 31−33 MeV and 75−95 MeV, respectively.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2016

New angles on energy correlation functions

Jet substructure observables, designed to identify specific features within jets, play an essential role at the Large Hadron Collider (LHC), both for searching for signals beyond the Standard Model and for testing QCD in extreme phase space regions. In this paper, we systematically study the structure of infrared and collinear safe substructure observables, defining a generalization of the energy correlation functions to probe n-particle correlations within a jet. These generalized correlators provide a flexible basis for constructing new substructure observables optimized for specific purposes. Focusing on three major targets of the jet substructure community — boosted top tagging, boosted W/Z/H tagging, and quark/gluon discrimination — we use power-counting techniques to identify three new series of powerful discriminants: M i , N i , and U i . The Mi series is designed for use on groomed jets, providing a novel example of observables with improved discrimination power after the removal of soft radiation. The N i series behave parametrically like the N -subjettiness ratio observables, but are defined without respect to subjet axes, exhibiting improved behavior in the unresolved limit. Finally, the U i series improves quark/gluon discrimination by using higher-point correlators to simultaneously probe multiple emissions within a jet. Taken together, these observables broaden the scope for jet substructure studies at the LHC.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2022

Time-resolved electron energy distribution functions at the substrate during a HiPIMS discharge with cathode voltage reversal

Abstract The time-dependent plasma properties of a high-power impulse magnetron sputtering plasma are investigated which include a positive ‘kick’ pulse on the sputtering target 2 μ s after the main negative pulse, this reversing the voltage on the cathode. At a substrate 15 cm distant from the magnetron, the time-dependent electron energy distribution function (EEDF), plasma potential, potential commute time and plasma diffusion properties are measured using a single Langmuir probe. Results show that the positive pulse on the target expels plasma and raises the plasma potential across the chamber on the order of 1 to 2 μ s, which is the time scale of the electron diffusion. The EEDF at the substrate fits a Druyvesteyn distribution during the main negative pulse rising slightly in average energy over time. The distribution is still Druyvesteyn and at the very start of the positive pulse, but then loses the higher energy electrons and drops in average electron energy as the positive pulse progresses. A Boltzmann equation solver, BOLSIG+, was used to predict the EEDF at the substrate during the positive pulse and it agrees best with the measurements assuming a value of 0.2 Td for the E / N (electric field/gas number density).

Physics↗

DOE OSTI2024

Phase space distribution functions and energy distributions of dark matter particles in haloes

ABSTRACT For a spherical dark matter halo with isotropic velocity distribution, the phase space distribution function (DF), the energy distribution, and the density profile form a set of self-consistent description of its equilibrium state, and knowing one is sufficient to determine the other two. The Navarro–Frenk–White density profile (NFW profile) is known to be a good approximation to the spherically averaged density distribution in simulated haloes. The DARKexp energy distribution is also known to compare well with the simulated energy distribution. We present a quantitative assessment of the NFW and DARKexp fits to the simulated DF and energy distribution for a wide range of haloes in a dark-matter-only simulation from the IllustrisTNG Project. As expected, we find that the NFW fits work well except at low energy when the density at small radii deviates from the NFW profile. Further, the NFW and DARKexp fits have comparable accuracy in the region where both fit well, but the DARKexp fits are better at low energy because they require matching of the central gravitational potential. We also find an approximate relation between the energy scale parametrizing the DARKexp energy distribution and that defined by the characteristic density and radius of the NFW profile. This relation may be linked to the relaxation process during halo formation.

Astronomy & Astrophysics↗

NASA NTRS1984

Cross-polar cap potential drop and the energy coupling function

A numerical model is formulated for correlations between the cross polar cap potential, the solar wind, and the AI index and compared with satellite data. The conductivity model was devised from empirically defined relationships between the Hall and Pedersen conductivities and the measured magnitude of the horizontal magnetic disturbance near College, AK. An instantaneous conductivity distribution was generated from data gathered at, e.g., six IMS meridian chain stations. A strong correlation appeared between an energy coupling function (ECF) and the polar cap potential drop. The ECF was derived from the chain station data, implying that potential drop estimates can be made from ground-based magnetometer data.

Ahn, B.-H.↗

DOE OSTI2012

Building a universal nuclear energy density functional (UNEDF)

The long-term vision initiated with UNEDF is to arrive at a comprehensive, quantitative, and unified description of nuclei and their reactions, grounded in the fundamental interactions between the constituent nucleons. We seek to replace current phenomenological models of nuclear structure and reactions with a well-founded microscopic theory that delivers maximum predictive power with well-quantified uncertainties. Specifically, the mission of this project has been three-fold: First, to find an optimal energy density functional (EDF) using all our knowledge of the nucleonic Hamiltonian and basic nuclear properties. Second, to apply the EDF theory and its extensions to validate the functional using all the available relevant nuclear structure and reaction data. Third, to apply the validated theory to properties of interest that cannot be measured, in particular the properties needed for reaction theory.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2023

On the convexity of phase-field fracture formulations: Analytical study and comparison of various degradation functions

Efficient and accurate fracture modeling is of great importance in applications where catastrophic outcomes under extreme scenarios are possible. The phase-field (PF) approach to fracture received significant attention over the past decade, due to its capability to capture complicated fracture patterns (e.g., crack merging and branching). Specifically, crack initiation and propagation are modeled via minimization of the total energy functional, which is regularized with the aid of a phase field. Despite the promising results and modeling capabilities of the PF method in many applications, the solution of fracture problems remains computationally challenging mainly due to the non-convexity of the total energy functional with respect to the combined unknown (phase field and displacement) fields. Understanding the effects of their coupling on convexity is crucial in order to address frequently encountered hurdles in fracture modeling (e.g., inefficient solvers and non-physical crack nucleation). In this paper, we develop convexity criteria for a wide class of PF fracture formulations. For this class of formulations, the second variation of the total energy functional is expressed in terms of Hessian matrices (evaluated at individual material points). Depending on the choice of geometric crack functions and degradation functions, we classify the formulations into three categories and analytically study each one separately. To study the sign of the second variation, we derive inequalities which are satisfied at material points when the Hessian matrix is locally positive semi-definite. These inequalities provide objective criteria for comparing degradation functions. Finally, the applicability of the proposed convexity criteria is demonstrated in the context of a one-dimensional problem, solved using a conventional monolithic solver.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2014

On the role of electron energy distribution function in double frequency heating of electron cyclotron resonance ion source plasmas

Double frequency heating (DFH) is a tool to improve the output of highly charged ions particularly from modern electron cyclotron resonance ion source installations with very high RF-frequencies. In order to gain information on the DFH-mechanism and on the role of the lower injected frequency we have carried out a series of dedicated experiments where we have put emphasis on the creation of a discrete resonance surface also for this lower frequency. Our well-established method of inserting an emissive MD (metal-dielectric) liner into the plasma chamber of the source is used in these experiments as a tool of investigation. In this way, the electron temperature and density for both ECR zones is increased in a controlled manner, allowing conclusions on the role of the change of the electron-energy-distribution function with and without DFH.

43 PARTICLE ACCELERATORS↗

DOE OSTI2014

The role of radiative reabsorption on the electron energy distribution functions in H{sub 2}/He plasma expansion through a tapered nozzle

A collisional-radiative model for the H{sub 2}/He plasma, coupled to a Boltzmann solver for the free electron kinetics is used to investigate the non-equilibrium conditions created in the expansion of an high-temperature plasma flow through a converging-diverging nozzle, starting from the steady state composition at T{sub 0}=10 000 K and p{sub 0}=1 atm in the reservoir. It is shown that the plasma optical thickness plays a major role on the evolution of macroscopic quantities and internal distributions along the nozzle axis. Structured electron energy distribution functions, characterized by long plateaux and humps, are created due to superelastic collisions of cold electrons and electronically excited atomic hydrogen. The magnitudes of the plateaux are orders of magnitude higher in an optically thick plasma compared with a thin plasma, while the electron-electron collisions play a role in smoothing the peaks created by superelastic collisions between cold electrons and H(n>2)

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2012

Building a Universal Nuclear Energy Density Functional

During the period of Dec. 1 2006 – Jun. 30, 2012, the UNEDF collaboration carried out a comprehensive study of all nuclei, based on the most accurate knowledge of the strong nuclear interaction, the most reliable theoretical approaches, the most advanced algorithms, and extensive computational resources, with a view towards scaling to the petaflop platforms and beyond. The long-term vision initiated with UNEDF is to arrive at a comprehensive, quantitative, and unified description of nuclei and their reactions, grounded in the fundamental interactions between the constituent nucleons. We seek to replace current phenomenological models of nuclear structure and reactions with a well-founded microscopic theory that delivers maximum predictive power with well-quantified uncertainties. Specifically, the mission of this project has been three-fold: First, to find an optimal energy density functional (EDF) using all our knowledge of the nucleonic Hamiltonian and basic nuclear properties; Second, to apply the EDF theory and its extensions to validate the functional using all the available relevant nuclear structure and reaction data; Third, to apply the validated theory to properties of interest that cannot be measured, in particular the properties needed for reaction theory.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2012

Building A Universal Nuclear Energy Density Functional (UNEDF)

During the period of Dec. 1 2006 - Jun. 30, 2012, the UNEDF collaboration carried out a comprehensive study of all nuclei, based on the most accurate knowledge of the strong nuclear interaction, the most reliable theoretical approaches, the most advanced algorithms, and extensive computational resources, with a view towards scaling to the petaflop platforms and beyond. The long-term vision initiated with UNEDF is to arrive at a comprehensive, quantitative, and unified description of nuclei and their reactions, grounded in the fundamental interactions between the constituent nucleons. We seek to replace current phenomenological models of nuclear structure and reactions with a well-founded microscopic theory that delivers maximum predictive power with well-quantified uncertainties. Specifically, the mission of this project has been three-fold: first, to find an optimal energy density functional (EDF) using all our knowledge of the nucleonic Hamiltonian and basic nuclear properties; second, to apply the EDF theory and its extensions to validate the functional using all the available relevant nuclear structure and reaction data; third, to apply the validated theory to properties of interest that cannot be measured, in particular the properties needed for reaction theory. The main physics areas of UNEDF, defined at the beginning of the project, were: ab initio structure; ab initio functionals; DFT applications; DFT extensions; reactions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2015

Fluid-solid transition in simple systems using density functional theory

A free energy functional for a crystal which contains both the symmetry-conserved and symmetry-broken parts of the direct pair correlation function has been used to investigate the fluid-solid transition in systems interacting via purely repulsive Weeks-Chandler-Anderson Lennard–Jones potential and the full Lennard–Jones potential. The results found for freezing parameters for the fluid-face centred cubic crystal transition are in very good agreement with simulation results. It is shown that although the contribution made by the symmetry broken part to the grand thermodynamic potential at the freezing point is small compared to that of the symmetry conserving part, its role is crucial in stabilizing the crystalline structure and on values of the freezing parameters.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2016

Mixed Arlequin method for multiscale poromechanics problems: Mixed Arlequin method for multiscale poromechanics problems

An Arlequin poromechanics model is introduced to simulate the hydro-mechanical coupling effects of fluid-infiltrated porous media across different spatial scales within a concurrent computational framework. A two-field poromechanics problem is first recast as the twofold saddle point of an incremental energy functional. We then introduce Lagrange multipliers and compatibility energy functionals to enforce the weak compatibility of hydro-mechanical responses in the overlapped domain. Here, to examine the numerical stability of this hydro-mechanical Arlequin model, we derive a necessary condition for stability, the twofold inf–sup condition for multi-field problems, and establish a modified inf–sup test formulated in the product space of the solution field. We verify the implementation of the Arlequin poromechanics model through benchmark problems covering the entire range of drainage conditions. Finally, through these numerical examples, we demonstrate the performance, robustness, and numerical stability of the Arlequin poromechanics model.

42 ENGINEERING↗

DOE OSTI2018

Optimal Control of the Rigid Inclusion Size in the Problem of Equilibrium of Inhomogeneous Timoshenko Type Plates Containing Cracks

We consider the equilibrium problem for plates with cracks located along rigid inclusions and analyze the dependence of the solution and the derivative of the energy functional on variations of the inclusion size. We establish the solvability of the optimal control problem, where the quality functional is given by the derivative of the energy functional, whereas the control parameter correspond to the inclusion size. A similar analysis is performed for the equilibrium problems for inhomogeneous plates with rigid delaminated inclusions.

97 MATHEMATICS AND COMPUTING↗