Search NASA⌕ Search

SEARCH · Search NASA

Results for “parabolic problems”

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.

117 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

A fast Huygens sweeping method for capturing paraxial multi-color optical self-focusing in nematic liquid crystals

We propose a numerically efficient algorithm for simulating the multi-color optical self-focusing phenomena in nematic liquid crystals. The propagation of the nematicon is modeled by a parabolic wave equation coupled with a nonlinear elliptic partial differential equation governing the angle between the crystal and the direction of propagation. Numerically, the paraxial parabolic wave equation is solved by a fast Huygens sweeping method, while the nonlinear elliptic PDE is handled by the alternating direction explicit (ADE) method. The overall algorithm is shown to be numerically efficient for computing high frequency beam propagations. - Highlights: • A simple strategy can already improve the fast Huygens sweeping method for high frequency asymptotic solutions. • The Alternating Direction Explicit method can significant improve the computational efficiency. • Numerical solutions match the Fredericks transition.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Development, Implementation, and Verification of Partially-Ionized Collisional Multifluid Plasma Models in Drekar

We consider the development of multifluid models for partially ionized multispecies plasmas. The models are composed of a standard set of five-moment fluid equations for each species plus a description of electromagnetics. The most general model considered utilizes a full set of fluid equations for each charge state of each atomic species, plus a set of fluid equations for electrons. The fluid equations are coupled through source terms describing electromagnetic coupling, ionization, recombination, charge exchange, and elastic scattering collisions in the low-density coronal limit. The form of each of these source terms is described in detail, and references for required rate coefficients are identified for a diverse range of atomic species. Initial efforts have been made to extend these models to incorporate some higher-density collisional effects, including ionization potential depression and three- body recombination. Some reductions of the general multifluid model are considered. First, a reduced multifluid model is derived which averages over all of the charge states (including neutrals) of each atomic species in the general multifluid model. The resulting model maintains full consistency with the general multifluid model from which it is derived by leveraging a quasi-steady-state collisional ionization equilibrium assumption to recover the ionization fractions required to make use of the general collision models. Further reductions are briefly considered to derive certain components of a single-fluid magnetohydrodynamics (MHD) model. In this case, a generalized Ohm's law is obtained, and the standard MHD resistivity is expressed in terms of the collisional models used in the general multifluid model. A number of numerical considerations required to obtain robust implementations of these multifluid models are discussed. First, an algebraic flux correction (AFC) stabilization approach for a continuous Galerkin finite element discretization of the multifluid system is described in which the characteristic speeds used in the stabilization of the fluid systems are synchronized across all species in the model. It is demonstrated that this synchronization is crucial in order to obtain a robust discretization of the multifluid system. Additionally, several different formulations are considered for describing the electromagnetics portion of the multifluid system using nodal continuous Galerkin finite element discretizations. The formulations considered include a parabolic divergence cleaning method and an implicit projection method for the traditional curl formulation of Maxwell's equations, a purely- hyperbolic potential-based formulation of Maxwell's equations, and a mixed hyperbolic-elliptic potential-based formulation of Maxwell's equations. Some advantages and disadvantages of each formulation are explored to compare solution robustness and the ease of use of each formulation. Numerical results are presented to demonstrate the accuracy and robustness of various components of our implementation. Analytic solutions for a spatially homogeneous damped plasma oscillation are derived in order to verify the implementation of the source terms for electromagnetic coupling and elastic collisions between fluid species. Ionization balance as a function of electron temperature is evaluated for several atomic species of interest by comparing to steady-state calculations using various sets of ionization and recombination rate coefficients. Several test problems in one and two spatial dimensions are used to demonstrate the accuracy and robustness of the discretization and stabilization approach for the fluid components of the multifluid system. This includes standard test problems for electrostatic and electromagnetic shock tubes in the two-fluid and ideal shock-MHD limits, a cylindrical diocotron instability, and the GEM challenge magnetic reconnection problem. A one-dimensional simplified prototype of an argon gas puff configuration as deployed on Sandia's Z-machine is used as a demonstration to exercise the full range of capabilities associated with the general multifluid model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

OPTOM: Optimization of Parabolic Trough - Operations & Maintenance

The US Department of Energy’s SunShot goals look to reduce the cost of Concentrating Solar Power (CSP) technology to 5¢/kWh for baseload plants. This is about a 50% reduction from current costs. To achieve this cost target, a significant reduction in operation and maintenance (O&M) costs of 40 to 50% is likely needed. Advances are needed in the O&M practices of CSP plants if the technology is to achieve the SunShot cost goals. Digitization of plant performance and O&M data has become a new best practice in the world of renewable energy asset management. Owners and operators of large photovoltaic and wind power plants are working to digitize performance and O&M data at their existing assets, to improve their management of the facilities, to increase performance, reduce O&M costs, and lower the overall life cycle cost of ownership. CSP power plants are behind the curve of other technologies on the digitization of plant information to aid in the plant asset management. This project directly addresses the objective of digitizing the O&M data of the solar field, focusing on three areas: 1) creating a framework for sharing data and information, 2) creating a system for monitoring and managing the maintenance of the solar field collectors, and 3) developing analytic tools to identify issues in the solar field. According to the NREL CSP Best Practices Study, the current practice at many CSP plants is to rely on paper lists, spreadsheets, and email for monitoring and managing problems and maintenance in the solar field. The key element to digitize solar field O&M is the creation of a centralized data archive that all users and systems can interface with. This project developed a centralized relational database framework that allows users and applications to access and share data. Conventional power plants utilize Computerized Maintenance Management Systems (CMMSs) to track the corrective, preventive (scheduled), and predictive maintenance of equipment and subsystems in the power plant. CSP plants use these systems in the power block, but while these systems specialize at tracking maintenance on up to thousands of pieces of equipment, they are not well suited for tracking the tens or hundreds of thousands of components in large commercial CSP or photovoltaic solar fields. In this project we developed a new software application referred to as FieldStatus (TM). This is a specialized database program that is used to track the status of each collector and its components. This application is designed to complement the existing CMMS to enable improved tracking and management of maintenance activities in the solar field. One of the major maintenance tasks for solar fields is maintaining the cleanliness of the mirrors. Although seemingly a relatively straight forward task, it has often proven challenging to maintain high levels of cleanliness in an efficient and cost-effective manner. This project developed new tools and metrics for monitoring and optimizing solar field cleaning resources and overall solar field cleanliness.

14 SOLAR ENERGY↗
DOE OSTI2018

Entropy Signatures of Topological Phase Transitions

We review the behavior of the entropy per particle in various two-dimensional electronic systems. The entropy per particle is an important characteristic of any many-body system that tells how the entropy of the ensemble of electrons changes if one adds one more electron. Recently, it has been demonstrated how the entropy per particle of a two-dimensional electron gas can be extracted from the recharging current dynamics in a planar capacitor geometry. These experiments pave the way to the systematic studies of entropy in various crystal systems including novel two-dimensional crystals such as gapped graphene, germanene, and silicene. Theoretically, the entropy per particle is linked to the temperature derivative of the chemical potential of the electron gas by the Maxwell relation. Using this relation, we calculate the entropy per particle in the vicinity of topological transitions in various two-dimensional electronic systems. We show that the entropy experiences quantized steps at the points of Lifshitz transitions in a two-dimensional electron gas with a parabolic energy spectrum. In contrast, in doubled-gapped Dirac materials, the entropy per particle demonstrates characteristic spikes once the chemical potential passes through the band edges. The transition from a topological to trivial insulator phase in germanene is manifested by the disappearance of a strong zero-energy resonance in the entropy per particle dependence on the chemical potential. We conclude that studies of the entropy per particle shed light on multiple otherwise hidden peculiarities of the electronic band structure of novel two-dimensional crystals.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2019

Investigation of dynamics of nematicons in liquid crystals by extended sinh-Gordon equation expansion method

This paper secures new nematicons in liquid crystals from its governing equation by extended sinh-Gordon equation expansion method with the aid of symbolic computational package Maple. The two laws of nonlinearity, namely Kerr and parabolic laws are studied. As outcomes, some new explicit complex hyperbolic and complex trigonometric function solutions are derived regarding dark, bright, combined dark–bright, singular, combined singular optical, periodic wave, dipole soliton solutions and others. To show the real physical meaning of the derived results, three dimensional graphs with contours and two dimensional graphs are prepared under the choice of proper values of the arbitrary parameters. The derived solutions have been checked back into their corresponding equations with the symbolic computational package Maple. All the combined solutions are found to be new, which may play an important role in the context of liquid crystals. The applied method is found to provide a powerful tool for extracting exact solitary wave solutions and overcome the difficulties of the solitary wave ansatz method.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2018

Dissipative tunnelling by means of scaled trajectories

Dissipative quantum tunnelling through an inverted parabolic barrier is considered in the presence of an electric field. A Schrödinger–Langevin or Kostin quantum–classical transition wave equation is used and applied resulting in a scaled differential equation of motion. A Gaussian wave packet solution to the resulting scaled Kostin nonlinear equation is assumed and compared to the same solution for the scaled linear Caldirola–Kanai equation. The resulting scaled trajectories are obtained at different dynamical regimes and friction cases, showing the gradual decoherence process in this open dynamics. Theoretical results show that the transmission probabilities are always higher in the Kostin approach than in the Caldirola–Kanai approach in the presence or not of an external electric field. This discrepancy should be understood due to the presence of an environment since the corresponding open dynamics should be governed by nonlinear quantum equations, whereas the second approach is issued from an effective Hamiltonian within a linear theory.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

A coupled ground heat flux–surface energy balance model of evaporation using thermal remote sensing observations

Abstract. One of the major undetermined problems in evaporation (ET) retrieval using thermal infrared remote sensing is the lack of a physically based ground heat flux (G) model and its integration within the surface energy balance (SEB) equation. Here, we present a novel approach based on coupling a thermal inertia (TI)-based mechanistic G model with an analytical surface energy balance model, Surface Temperature Initiated Closure (STIC, version STIC1.2). The coupled model is named STIC-TI. The model is driven by noon–night (13:30 and 01:30 local time) land surface temperature, surface albedo, and a vegetation index from MODIS Aqua in conjunction with a clear-sky net radiation sub-model and ancillary meteorological information. SEB flux estimates from STIC-TI were evaluated with respect to the in situ fluxes from eddy covariance measurements in diverse ecosystems of contrasting aridity in both the Northern Hemisphere and Southern Hemisphere. Sensitivity analysis revealed substantial sensitivity of STIC-TI-derived fluxes due to the land surface temperature uncertainty. An evaluation of noontime G (Gi) estimates showed 12 %–21 % error across six flux tower sites, and a comparison between STIC-TI versus empirical G models also revealed the substantially better performance of the former. While the instantaneous noontime net radiation (RNi) and latent heat flux (LEi) were overestimated (15 % and 25 %), sensible heat flux (Hi) was underestimated (22 %). Overestimation (underestimation) of LEi (Hi) was associated with the overestimation of net available energy (RNi−Gi) and use of unclosed surface energy balance flux measurements in LEi (Hi) validation. The mean percent deviations in Gi and Hi estimates were found to be strongly correlated with satellite day–night view angle difference in parabolic and linear pattern, and a relatively weak correlation was found between day–night view angle difference versus LEi deviation. Findings from this parameter-sparse coupled G–ET model can make a valuable contribution to mapping and monitoring the spatiotemporal variability of ecosystem water stress and evaporation using noon–night thermal infrared observations from future Earth observation satellite missions such as TRISHNA, LSTM, and SBG.

Bhattacharya, Bimal K.↗
DOE OSTI2018

A Study of the von Neumann Stability Analysis of a Semi-implicit Numerical Method Applied Over the Radial Impurity Transport Equation in Tokamak Plasma

The radial impurity transport equation for tokamak plasma is a form of diffusion–convection–reaction equation. The impurity transport equation is solved to determine the distribution of impurity (non-fuel) ion species with different ionization states perpendicular to magnetic surfaces of tokamak plasma. The equation for each charge (ionization) state Z is a non-linear, second-order in space, first-order in time, parabolic partial differential equation coupled to the previous Z − 1 and the next Z + 1 charge states of the impurity species through its reaction term. The number of differential equations to be solved simultaneously is hence determined by the number of ionization states of the impurity species studied. The solution to the set of these coupled equations can be obtained using a semi-implicit numerical method applied on it. The present study describes the application of von Neumann stability analysis over the semi-implicit numerical method applied over the radial impurity transport equation and determines a generic stability criterion for the method. The stability analysis is further illustrated using the geometry of Aditya tokamak installed at the Institute for Plasma Research Gandhinagar, India as an example. The impurity species considered is oxygen (Atomic number = 8). This leads to a set of eight coupled equations for charge states Z = 1 to 8 over which von Neumann analysis is illustrated in present study.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Bound-preserving finite element approximations of the Keller–Segel equations

We report this paper aims to develop numerical approximations of the Keller–Segel equations that mimic at the discrete level the lower bounds and the energy law of the continuous problem. We solve these equations for two unknowns: the organism (or cell) density, which is a positive variable, and the chemoattractant density, which is a non-negative variable. We propose two algorithms, which combine a stabilized finite element method and a semi-implicit time integration. The stabilization consists of a nonlinear artificial diffusion that employs a graph-Laplacian operator and a shock detector that localizes local extrema. As a result, both algorithms turn out to be nonlinear and can generate cell and chemoattractant numerical densities fulfilling lower bounds. However, the first algorithm requires a suitable constraint between the space and time discrete parameters, whereas the second one does not. We design the latter to attain a discrete energy law on acute meshes. We report some numerical experiments to validate the theoretical results on blowup and nonblowup phenomena. In the blowup setting, we identify a locking phenomenon that relates the L ∞ (Ω)-norm to the L 1 (Ω)-norm limiting the growth of the singularity when supported on a macroelement.

97 MATHEMATICS AND COMPUTING↗