SEARCH · Search NASA

Results for “fully linear model”

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 37 records · Page 2

All sources NASA NTRS DOE OSTI

DOE OSTI2016

Coarse-grained description of cosmic structure from Szekeres models

We show that the full dynamical freedom of the well known Szekeres models allows for the description of elaborated 3-dimensional networks of cold dark matter structures (over-densities and/or density voids) undergoing ''pancake'' collapse. By reducing Einstein's field equations to a set of evolution equations, which themselves reduce in the linear limit to evolution equations for linear perturbations, we determine the dynamics of such structures, with the spatial comoving location of each structure uniquely specified by standard early Universe initial conditions. By means of a representative example we examine in detail the density contrast, the Hubble flow and peculiar velocities of structures that evolved, from linear initial data at the last scattering surface, to fully non-linear 10–20 Mpc scale configurations today. To motivate further research, we provide a qualitative discussion on the connection of Szekeres models with linear perturbations and the pancake collapse of the Zeldovich approximation. This type of structure modelling provides a coarse grained—but fully relativistic non-linear and non-perturbative —description of evolving large scale cosmic structures before their virialisation, and as such it has an enormous potential for applications in cosmological research.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS1984

A numerical simulation of barotropic instability. II Wave-wave interaction

A fully nonlinear numerical model of the point jet barotropic instability is used to test and confirm the hypothesis that the magnitude of the wave vorticity does not exceed the magnitude of the initial shear. This result arises directly from the local conservation of vorticity following a parcel and the fact that unstable waves are principally confined to the region where the zonal mean vorticity can be smoothed by the wave so as to eliminate the instability. Comparisons are made between fully nonlinear and quasi-linear models of the point jet instability and their tracer transport properties. Differences become particularly evident after wave saturation. The most important effect neglected by the wave-mean flow model appears to be the advection of wave vorticity by the most unstable mode. However, as equilibration of the instability proceeds, the globally averaged properties of both models are found to be similar.

Nielsen, J. E.↗

NASA NTRS1993

An efficient model for coupling structural vibrations with acoustic radiation

The scattering of an incident wave by a flexible panel is studied. The panel vibration is governed by the nonlinear plate equations while the loading on the panel, which is the pressure difference across the panel, depends on the reflected and transmitted waves. Two models are used to calculate this structural-acoustic interaction problem. One solves the three dimensional nonlinear Euler equations for the flow-field coupled with the plate equations (the fully coupled model). The second uses the linear wave equation for the acoustic field and expresses the load as a double integral involving the panel oscillation (the decoupled model). The panel oscillation governed by a system of integro-differential equations is solved numerically and the acoustic field is then defined by an explicit formula. Numerical results are obtained using the two models for linear and nonlinear panel vibrations. The predictions given by these two models are in good agreement but the computational time needed for the 'fully coupled model' is 60 times longer than that for 'the decoupled model'.

Frendi, Abdelkader↗

DOE OSTI2019

Gyrokinetic simulations of m=0 mode in sheared flow Z-pinch

Axisymmetric stability properties of sheared flow Z-pinch plasmas are studied by making use of the gyrokinetic approximation in the long-wavelength limit. Numerical simulations are carried out with the high-order finite-volume code COntinuum Gyrokinetic Edge New Technology (COGENT) and are analyzed for the parameters characteristic of the FuZE experiment. Reduction of the linear growth rate with increasing shear is observed, and the results are elucidated by making use of a local dispersion relation analysis. In addition, COGENT simulations are compared with fully kinetic particle-in-cell simulations, and with an ideal magnetohydrodynamics (MHD) model. Good agreement between the gyrokinetic and fully kinetic results for the linear stability is found, with the gyrokinetic model requiring much less computational time due to its ability to step over particle gyroperiod. The ideal MHD model is found to be consistent with the kinetic models in the long-wavelength part of the spectra ($kρ_i$), while failing to adequately predict short-scale ($kρ_i$) stability. Here, $k$ is the axial wavelength vector and $ρ_i$ is the ion gyroradius.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2014

Is the continuous matter creation cosmology an alternative to ΛCDM?

The matter creation cosmology is revisited, including the evolution of baryons and dark matter particles. The creation process affects only dark matter and not baryons. The dynamics of the ΛCDM model can be reproduced only if two conditions are satisfied: 1) the entropy density production rate and the particle density variation rate are equal and 2) the (negative) pressure associated to the creation process is constant. However, the matter creation model predicts a present dark matter-to-baryon ratio much larger than that observed in massive X-ray clusters of galaxies, representing a potential difficulty for the model. In the linear regime, a fully relativistic treatment indicates that baryons are not affected by the creation process but this is not the case for dark matter. Both components evolve together at early phases but lately the dark matter density contrast decreases since the background tends to a constant value. This behaviour produces a negative growth factor, in disagreement with observations, being a further problem for this cosmology.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2021

Simultaneous Decomposition and Dewetting of Nanoscale Alloys: A Comparison of Experiment and Theory

In this work, we consider the coupled process of phase separation and dewetting of metal alloys of nanoscale thickness deposited on solid substrates. The experiments involve applying nanosecond laser pulses that melt the Ag40Ni60 alloy films in two setups: either on thin supporting membranes or on bulk substrates. These two setups allow for extracting both temporal and spatial scales on which the considered processes occur. The theoretical model involves a longwave version of the Cahn–Hilliard formulation used to describe spinodal decomposition, coupled with an asymptotically consistent longwave-based description of dewetting that occurs due to destabilizing interactions between the alloy and the substrate, modeled using the disjoining pressure approach. Careful modeling, combined with linear stability analysis and fully nonlinear simulations, leads to results consistent with the experiments. In particular, we find that the two instability mechanisms occur concurrently, with the phase separation occurring faster and on shorter temporal scales. The modeling results show a strong influence of the temperature dependence of relevant material properties, implying that such a dependence is crucial for the understanding of the experimental findings. The agreement between theory and experiment suggests the utility of the proposed theoretical approach in helping to develop further experiments directed toward formation of metallic alloy nanoparticles of desired properties.

36 MATERIALS SCIENCE↗

DOE OSTI2016

Linear multispecies gyrokinetic flux tube benchmarks in shaped tokamak plasmas

Verification is the fundamental step that any turbulence simulation code has to be submitted in order to assess the proper implementation of the underlying equations. We have carried out a cross comparison of three flux tube gyrokinetic codes, GENE [F. Jenko et al., Phys. Plasmas 7, 1904 (2000)], GKW [A. G. Peeters et al., Comput. Phys. Commun. 180, 2650 (2009)], and GS2 [W. Dorland et al., Phys. Rev. Lett. 85, 5579 (2000)], focusing our attention on the effect of realistic geometries described by a series of MHD equilibria with increasing shaping complexity. To simplify the effort, the benchmark has been limited to the electrostatic collisionless linear behaviour of the system. A fully gyrokinetic model has been used to describe the dynamics of both ions and electrons. Several tests have been carried out looking at linear stability at ion and electron scales, where for the assumed profiles Ion Temperature Gradient (ITG)/Trapped Electron Modes and Electron Temperature Gradient modes are unstable. The capability of the codes to handle a non-zero ballooning angle has been successfully benchmarked in the ITG regime. Finally, the standard Rosenbluth-Hinton test has been successfully carried out looking at the effect of shaping on Zonal Flows (ZFs) and Geodesic Acoustic Modes (GAMs). Inter-code comparison as well as validation of simulation results against analytical estimates has been accomplished. All the performed tests confirm that plasma elongation strongly stabilizes plasma instabilities as well as leads to a strong increase in ZF residual and GAM damping.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1995

Reverberation mapping by regularized linear inversion

Reverberation mapping of active galactic nucleus (AGN) emission-line regions requires the numerical deconvolution of two time series. We suggest the application of a new method, regularized linear inversion, to the solution of this problem. This method possesses many good features; it imposes no restrictions on the sign of the response function; it can provide clearly defined uncertainty estimates; it involves no guesswork about unmeasured data; it can give a clear indication of when the underlying convolution model is inadequate; and it is computationally very efficient. Using simulated data, we find the minimum S/N and length of the time series in order for this method to work satisfactorily. We also define guidelines for choosing the principal tunable parameter of the method and for interpreting the results. Finally, we reanalyze published data from the 1989 NGC 5548 campaign using this new method and compare the results to those previously obtained by maximum entropy analysis. For some lines we find good agreement, but for others, especially C III lambda(1909) and Si IV lambda(1400), we find significant differences. These can be attributed to the inability of the maximum entropy method to find negative values of the response function, but also illustrate the nonuniqueness of any deconvolution technique. We also find evidence that certain line light curves (e.g., C IV lambda(1549)) cannot be fully described by the simple linear convolution model.

Krolik, Julian H.↗

DOE OSTI2020

A Bayesian Approach to Regional Decadal Predictability: Sparse Parameter Estimation in High-Dimensional Linear Inverse Models of High-Latitude Sea Surface Temperature Variability

Stochastic reduced models are an important tool in climate systems whose many spatial and temporal scales cannot be fully discretized or underlying physics may not be fully accounted for. One form of reduced model, the linear inverse model (LIM), has been widely used for regional climate predictability studies—typically focusing more on tropical or midlatitude studies. However, most LIM fitting techniques rely on point estimation techniques deriving from fluctuation–dissipation theory. In this methodological study we explore the use of Bayesian inference techniques for LIM parameter estimation of sea surface temperature (SST), to quantify the skillful decadal predictability of Bayesian LIM models at high latitudes. We show that Bayesian methods, when compared to traditional point estimation methods for LIM-type models, provide better calibrated probabilistic skill, while simultaneously providing better point estimates due to the regularization effect of the prior distribution in high-dimensional problems. We compare the effect of several priors, as well as maximum likelihood estimates, on 1) estimating parameter values on a perfect model experiment and 2) producing calibrated 1-yr SST anomaly forecast distributions using a preindustrial control run of the Community Earth System Model (CESM). Finally, we employ a host of probabilistic skill metrics to determine the extent to which an LIM can forecast SST anomalies at high latitudes. We find that the choice of prior distribution has an appreciable impact on estimation outcomes, and priors that emphasize physically relevant properties enhance the model’s ability to capture variability of SST anomalies.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2014

Debye scale turbulence within the electron diffusion layer during magnetic reconnection

During collisionless, anti-parallel magnetic reconnection, the electron diffusion layer is the region of both fieldline breaking and plasma mixing. Due to the in-plane electrostatic fields associated with collisionless reconnection, the inflowing plasmas are accelerated towards the X-line and form counter-streaming beams within the unmagnetized diffusion layer. This configuration is inherently unstable to in-plane electrostatic streaming instabilities provided that there is sufficient scale separation between the Debye length λ{sub D} and the electron skin depth c/ω{sub pe}. This scale separation has hitherto not been well resolved in kinetic simulations. Using both 2D fully kinetic simulations and a simple linear model, we demonstrate that these in-plane streaming instabilities generate Debye scale turbulence within the electron diffusion layer at electron temperatures relevant to magnetic reconnection both in the magnetosphere and in laboratory experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2016

A comparison of WEC control strategies

The operation of Wave Energy Converter (WEC) devices can pose many challenging problems to the Water Power Community. A key research question is how to significantly improve the performance of these WEC devices through improving the control system design. This report summarizes an effort to analyze and improve the performance of WEC through the design and implementation of control systems. Controllers were selected to span the WEC control design space with the aim of building a more comprehensive understanding of different controller capabilities and requirements. To design and evaluate these control strategies, a model scale test-bed WEC was designed for both numerical and experimental testing (see Section 1.1). Seven control strategies have been developed and applied on a numerical model of the selected WEC. This model is capable of performing at a range of levels, spanning from a fully-linear realization to varying levels of nonlinearity. The details of this model and its ongoing development are described in Section 1.2.

30 DIRECT ENERGY CONVERSION↗

DOE OSTI2020

Performance evaluation of surface riding wave energy converter with linear electric generator

In this work, we devised a new WEC (wave energy converter) called SR-WEC (Surface Riding WEC). The SR-WEC consists of two bodies: the outer cylinder with an armature assembly (body #1) and a magnet assembly (body #2) sliding inside the armature. For the SR-WEC, the relative sliding displacement and velocity are caused by gravity acceleration and the outer cylinder's motions, and they lead to electrical power generation. To evaluate its performance, a numerical simulation tool was developed, which solves the fully-coupled floater-mooring-generator dynamics. During the developing stage, the appropriate hydrodynamics model, sliding mechanics model, mooring dynamics model, and LEG (linear electric generator) electro-magnetic model were independently developed and then fully coupled in time domain to account for the cross-coupling interactions among them. Then, the developed simulation tool was verified component by component against various laboratory tests. Subsequently, systematic parametric studies were conducted with several important design parameters under various wave conditions to enhance power generation. After that, the average output power was evaluated in enlarged operational wave conditions. The present SR-WEC is particularly designed to be efficient at low sea states, which is good since they cover the majority of typical annual sea states.

42 ENGINEERING↗

DOE OSTI2018

Scattering forms and the positive geometry of kinematics, color and the worldsheet

The search for a theory of the S-Matrix over the past five decades has revealed surprising geometric structures underlying scattering amplitudes ranging from the string worldsheet to the amplituhedron, but these are all geometries in auxiliary spaces as opposed to the kinematical space where amplitudes actually live. Motivated by recent advances providing a reformulation of the amplituhedron and planar N = 4 SYM amplitudes directly in kinematic space, we propose a novel geometric understanding of amplitudes in more general theories. The key idea is to think of amplitudes not as functions, but rather as differential forms on kinematic space. We explore the resulting picture for a wide range of massless theories in general spacetime dimensions. For the bi-adjoint Φ3 scalar theory, we establish a direct connection between its “scattering form” and a classic polytope — the associahedron — known to mathematicians since the 1960’s. We find an associahedron living naturally in kinematic space, and the tree level amplitude is simply the “canonical form” associated with this “positive geometry”. Fundamental physical properties such as locality and unitarity, as well as novel “soft” limits, are fully determined by the combinatorial geometry of this polytope. Furthermore, the moduli space for the open string worldsheet has also long been recognized as an associahedron. We show that the scattering equations act as a diffeomorphism between the interior of this old “worldsheet associahedron” and the new “kinematic associahedron”, providing a geometric interpretation and simple conceptual derivation of the bi-adjoint CHY formula. We also find “scattering forms” on kinematic space for Yang-Mills theory and the Non-linear Sigma Model, which are dual to the fully color-dressed amplitudes despite having no explicit color factors. This is possible due to a remarkable fact—“Color is Kinematics”— whereby kinematic wedge products in the scattering forms satisfy the same Jacobi relations as color factors. Finally, all our scattering forms are well-defined on the projectivized kinematic space, a property which can be seen to provide a geometric origin for color-kinematics duality.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2015

Past and present cosmic structure in the SDSS DR7 main sample

We present a chrono-cosmography project, aiming at the inference of the four dimensional formation history of the observed large scale structure from its origin to the present epoch. To do so, we perform a full-scale Bayesian analysis of the northern galactic cap of the Sloan Digital Sky Survey (SDSS) Data Release 7 main galaxy sample, relying on a fully probabilistic, physical model of the non-linearly evolved density field. Besides inferring initial conditions from observations, our methodology naturally and accurately reconstructs non-linear features at the present epoch, such as walls and filaments, corresponding to high-order correlation functions generated by late-time structure formation. Our inference framework self-consistently accounts for typical observational systematic and statistical uncertainties such as noise, survey geometry and selection effects. We further account for luminosity dependent galaxy biases and automatic noise calibration within a fully Bayesian approach. As a result, this analysis provides highly-detailed and accurate reconstructions of the present density field on scales larger than ∼ 3 Mpc/h, constrained by SDSS observations. This approach also leads to the first quantitative inference of plausible formation histories of the dynamic large scale structure underlying the observed galaxy distribution. The results described in this work constitute the first full Bayesian non-linear analysis of the cosmic large scale structure with the demonstrated capability of uncertainty quantification. Some of these results will be made publicly available along with this work. The level of detail of inferred results and the high degree of control on observational uncertainties pave the path towards high precision chrono-cosmography, the subject of simultaneously studying the dynamics and the morphology of the inhomogeneous Universe.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

Integrating very-high-resolution UAS data and airborne imaging spectroscopy to map the fractional composition of Arctic plant functional types in Western Alaska

Widespread changes in vegetation cover and composition are driving strong impacts on Arctic ecosystem functioning and global climate feedbacks. An accurate characterization of tundra vegetation composition is required to understand how the Arctic will respond to future climate change. However, quantifying tundra vegetation composition over large areas is challenging as commonly-used satellite observations are too coarse, spatially and spectrally, to differentiate low-lying tundra vegetation types. Recent airborne and spaceborne imaging spectroscopy platforms provide better data to characterize vegetation composition. Yet, our ability to characterize vegetation composition with imaging spectroscopy remains largely unexplored in the Arctic, particularly due to a lack of ground observations needed to train and test classification models. To address this problem, we collected very-high-resolution (VHR, ~5 cm) unoccupied aerial system (UAS) imagery at three low-Arctic tundra sites located on the Seward Peninsula, western Alaska. In this paper, we examine the feasibility of integrating imagery from the UAS and the hyperspectral Airborne Visible/Infrared Imaging Spectrometer, Next Generation (AVIRIS-NG) airborne instrument to map the fractional composition of 12 key Arctic plant functional types (PFTs). To this end, we first mapped the 12 PFTs from our VHR UAS imagery using random forest classification. We then used these UAS-derived PFT maps as ground truth to develop partial least squares regression (PLSR) models to predict the fractional cover (FCover) of each PFT from AVIRIS-NG imagery. Further, we evaluated the performance of our PLSR models using reserved UAS samples, as well as by mapping PFT FCover and dominant PFT for large tundra landscapes. Our results show that 1) Arctic PFTs can be effectively mapped using VHR UAS imagery, with overall accuracy between 86% and 92%, 2) when the UAS mapped PFTs were used to inform PLSR scaling models, the FCover of the 12 PFTs could be effectively estimated from AVIRIS-NG imagery with a mean absolute error (MAE) <0.13, and 3) our PLSR models outperformed traditional, fully constrained least-squares (FCLS) linear mixture analysis and produced high-quality, spatially contiguous PFT FCover and PFT maps that captured vegetation spatial patterns with similar accuracy to those developed from UAS imagery. The developed PLSR models have the potential to be broadly applied for quantifying vegetation composition with AVIRIS-NG images to help monitor tundra vegetation dynamics and improve process-based modeling of tundra ecosystems.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS1988

Irreversible wave-mean flow interactions in a mechanistic model of the stratosphere

A wave-1 minor warming is simulated in a mechanistic, global, primitive equation model. The modification of the zonal flow by the rearrangement of potential vorticity on a middle stratopheric isentrope is compared in a fully nonlinear model and in a model with only one wave and the zonal flow (a quasi-linear model). The permanent rearrangement of potential vorticity during the wave episode is more intense and more localized meridionally in the fully nonlinear model, which is able to capture the process of planetary wave breaking in some detail. Additional experiments reveal that the differences between the quasi-linear and nonlinear models persist for a broad range of wave amplitudes, and that the quasi-linear model can qualitatively reproduce the modification of the zonal flow by the wave when the diffusive dissipation of the wave is enhanced. These results are discussed in the context of the theory of barotropic Rossby waves in shear flows, and in comparison with recent numerical simulations of the middle atmosphere.

Robinson, Walter A.↗

DOE OSTI2025

A Performance Portable, Fully Implicit Landau Collision Operator with Batched Linear Solvers

Modern accelerators use hierarchical parallel programming models that enable massive multithreading within a processing element (PE), with multiple PEs per device driven by traditional processes. Batching is a technique for exposing PE-level parallelism in algorithms that have traditionally run on MPI processes or multiple threads within a single process. Opportunities for batching arise in, for example, kinetic discretizations of magnetized plasmas where collisions are advanced in velocity space at each spatial point independently. This paper builds on previous work on a high-performance, fully nonlinear, Landau collision operator by batching the linear solver, as well as batching the spatial point problems and adding new support for multiple grids for multiscale, multispecies problems. An anisotropic relaxation verification test that agrees well with previously published results and analytical models is presented. The performance results from NVIDIA A100 and AMD MI250X nodes are presented with hardware utilization analysis for each architecture. Finally, the entire implicit Landau operator time advance is implemented in Kokkos for performance portability, running entirely on the device and is available in the PETSc numerical library.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2022

Data-driven, multi-moment fluid modeling of Landau damping

Deriving governing equations of complex physical systems based on first principles can be quite challenging when there are certain unknown terms and hidden physical mechanisms in the systems. In this work, we apply a deep learning architecture to learn fluid partial differential equations (PDEs) of a plasma system based on the data acquired from a fully kinetic model. Here, the learned multi-moment fluid PDEs are demonstrated to incorporate kinetic effect such as Landau damping. Based on the learned fluid closure, the data-driven, multi-moment fluid modeling can well reproduce all the physical quantities derived from the fully kinetic model. The calculated damping rate of Landau damping is consistent with both the fully kinetic simulation and the linear theory. The data-driven fluid modeling of PDEs for complex physical systems may be applied to improve the fluid closure and reduce the computational cost of multi-scale modeling of global systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗