Search NASA⌕ Search

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.

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

At least 739 records · Page 41
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Quantum-Accurate Multiscale Modeling of Shock Hugoniots, Ramp Compression Paths, Structural and Magnetic Phase Transitions, and Transport Properties in Highly Compressed Metals

Fully characterizing high energy density (HED) phenomena using pulsed power facilities (Z machine) and coherent light sources is possible only with complementary numerical modeling for design, diagnostic development, and data interpretation. The exercise of creating numerical tests, that match experimental conditions, builds critical insight that is crucial for the development of a strong fundamental understanding of the physics behind HED phenomena and for the design of next generation pulsed power facilities. The persistence of electron correlation in HED materials arising from Coulomb interactions and the Pauli exclusion principle is one of the greatest challenges for accurate numerical modeling and has hitherto impeded our ability to model HED phenomena across multiple length and time scales at sufficient accuracy. An exemplar is a ferromagnetic material like iron, while familiar and widely used, we lack a simulation capability to characterize the interplay of structure and magnetic effects that govern material strength, kinetics of phase transitions and other transport properties. Herein we construct and demonstrate the Molecular-Spin Dynamics (MSD) simulation capability for iron from ambient to earth core conditions, all software advances are open source and presently available for broad usage. These methods are multi-scale in nature, direct comparisons between high fidelity density functional theory (DFT) and linear-scaling MSD simulations is done throughout this work, with advancements made to MSD allowing for electronic structure changes being reflected in classical dynamics. Main takeaways for the project include insight into the role of magnetic spins on mechanical properties and thermal conductivity, development of accurate interatomic potentials paired with spin Hamiltonians, and characterization of the high pressure melt boundary that is of critical importance to planetary modeling efforts.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2017

Fast neutrino flavor conversions near the supernova core with realistic flavor-dependent angular distributions

It has been recently pointed out that neutrino fluxes from a supernova can show substantial flavor conversions almost immediately above the core. Using linear stability analyses and numerical solutions of the fully nonlinear equations of motion, we perform a detailed study of these fast conversions , focussing on the region just above the supernova core. We carefully specify the instabilities for evolution in space or time, and find that neutrinos travelling towards the core make fast conversions more generic, i.e., possible for a wider range of flux ratios and angular asymmetries that produce a crossing between the zenith-angle spectra of ν {sub e} and ν-bar {sub e} . Using fluxes and angular distributions predicted by supernova simulations, we find that fast conversions can occur within tens of nanoseconds, only a few meters away from the putative neutrinospheres. If these fast flavor conversions indeed take place, they would have important implications for the supernova explosion mechanism and nucleosynthesis.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Considering coasts: Adapting terrestrial models to characterize coastal wetland ecosystems

The Energy Exascale Earth System Model (E3SM) simulates fully coupled processes and interactions among water, energy, carbon and nutrient cycles. E3SM connects vegetation and soil dynamics through nutrient uptake, plant production, litterfall and decomposition as a function of abiotic parameters (e.g. temperature and moisture). However, E3SM is designed to characterize terrestrial ecosystems and connects land and open ocean systems using a single streamflow transport term, ignoring the complex dynamics of energy, water, carbon, and nutrients in coastal systems. The goals of our project were to: (1) Parameterize a point version of E3SM to capture coastal wetland habitats and (2) Determine marsh community responses to increased temperature and elevated CO 2 . We adapted a version of the E3SM land model, previously configured to represent forested bog hydrology to a coastal ecosystem using datasets from field experiments conducted at the Smithsonian Environmental Research Center's Global Change Research Wetland (GCReW). Tidal forcing in a marsh environment was simulated using a two-column system in which the columns are connected by lateral hydrologic flows. One column simulates interactions between vegetation and soil while a second column simulates variation in water level (both tidal and sea level rise). The updated model captures many aspects of the field experiments, showing that plant community responses to environmental change are non-linear, non-additive and different between plant types. Elevated CO 2 treatments increased C 3 plant biomass more than C 4 (33% vs 17%). Temperature exacerbated CO 2 responses in C 3 plants (0 °C: 26%, 5.1 °C: 56%). We were more successful at characterizing C 3 than C 4 responses and simulating above rather than belowground biomass production. Furthermore, the next steps will include updates to key physiological parameters such as root:shoot carbon allocation and the addition of mechanistic feedbacks between vegetation and biogeochemical processes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Interpretation of multiple solutions in fully iterative GF2 and GW schemes using local analysis of two-particle density matrices

Due to the presence of non-linear equations, iterative Green’s function methods can result in multiple different solutions even for simple molecular systems. In contrast to the wave-function methods, a detailed and careful analysis of such molecular solutions was not performed before. In this work, we use two-particle density matrices to investigate local spin and charge correlators that quantify the charge resonance and covalent characters of these solutions. When applied within the unrestricted orbital set, spin correlators elucidate the broken symmetry of the solutions, containing necessary information for building effective magnetic Hamiltonians. Based on GW and GF2 calculations of simple molecules and transition metal complexes, we construct Heisenberg Hamiltonians, four-spin-four-center corrections, and biquadratic spin–spin interactions. These Hamiltonian parameterizations are compared to previous wave-function calculations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Reactive Modified Epoxy Resin and Its Miscible Blends Based on Recycled Oligomers from Solvolysis

Chemical depolymerization of fully cured epoxy resin with 20% reactive modifier was successfully performed via a solvent-assisted solvolysis process into low molecular weight recyclable oligomers (RO) at 240 °C in a pressure vessel at 650 psi for 4 h. The thermoset epoxy resin was depolymerized into transparent brown viscous fluid with a higher viscosity than the uncured epoxy resin with approximately 93% yield. Different concentrations of the RO were homogeneously mixed with the pure epoxy resin, and their curing kinetics, viscosity, FTIR, mechanical properties, DMA, and cross-link density were investigated. The curing kinetics of the pure reactive modified epoxy resin (baseline) and its mixtures with RO of different concentrations were investigated under both isothermal and nonisothermal conditions using small amplitude oscillatory shear flow. The elastic and viscous moduli (G′ and G″), complex viscosity (η*), and tan δ values were evaluated at different curing times and temperatures. The G′, G″, and η* increased dramatically, while tan δ decreased strongly by several orders of magnitude at the gel point. The zero-shear viscosity (η 0 ) was determined from the angular frequency dependent on η* based on the Cross model for different blend compositions in the liquid state before curing. The composition dependence of η 0 showed a positive deviation from the linear mixing rule and was well described by the Lecyar model. Here, the apparent activation energy of curing (E a ) was also evaluated according to the Arrhenius equation and was found to be 46 ± 2 kJ/mol regardless of the different contents of RO. For all blends up to 40 wt % RO, only one tan δ peak systematically shifting to lower temperatures with increasing content of RO was observed in the DMA measurements, indicating that the epoxy resin and the RO are miscible with up to 40 wt % RO.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2014

SU-F-BRE-10: Methods to Simulate and Measure the Attenuation for Modeling a Couch Top with Rails for FFF Treatment Delivery On the Varian Edge Linac

Purpose: To measure attenuation for modelling of the KVue Couchtop for 6X and 10X FFF SRS/SBRT treatment Methods: Treatment planning simulation studies were done using 6X FFF beams to estimate the dosimetric impact of KVue couchtops (including the Q-Fix IGRT [carbon fiber] and Calypso [nonconductive Kevlar material]) with a structure model obtained from a research workstation (Eclipse, advanced planning interface (API) v13). Prior to installation on the Varian Edge linac, the couchtop along with (Kevlar) rails were CT scanned with the rails at various positions. An additional scan with the couchtop 15cm above the CT table top was obtained with 20cm solid water to facilitate precised/indexed data acquisition. Measurements for attenuation were obtained for field sizes of 2, 4 and 10 cm{sup 2} at 42 gantry angles including 6 pairs of opposing fields and other angles for oblique delivery where the beams traversed the couchtop and or rails. The delivery was fully automated with xml scripts running in developer mode. The results were then used to determine an accurate structure model for AAA (Eclipse v11) planning of IMRT and RapidArc delivery. Results: The planning simulation relative dose attenuation for oblique entry was not significantly different than the Exact IGRT or BrainLab iBeam couch except that the rails added 6% additional attenuation. The relative attenuation measurements for PA, PA (rails: inner position), oblique, oblique (rails: outer position), oblique (rails: inner position) were: −2.0%, −2.5%, −15.6%, −2.5%, −5.0% for 6X FFF and −1.4%, −1.5%, −12.2%, − 2.5%, −5.0% for 10X FFF with slight decrease in attenuation versus field size. A Couch structure model (with HU values) was developed. Calculation compared to measurement showed good agreement except for oblique (rails: outer position) where differences approached a magnitude of 6%. Conclusion: A model of the couch structures has been developed accounting for attenuation for FFF beams.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2016

Continuous Structural Transition in Glass-Forming Molten Titanate BaTi 2 O 5

The structure of the model titanate glass former BaTi2O5 has been studied over a wide temperature (T) range in the molten, supercooled, and glassy states under conditions of aerodynamic levitation. Both high-energy X-ray diffraction and Ti K-edge X-ray absorption spectroscopy reveal a continuous structural transition involving reduction of the cation-oxygen (and oxygen-cation) average coordination numbers and bond lengths with increasing T. Ti-0 coordination in the moderately supercooled and equilibrium melt follows a linear trend n(Tio) = 5.4(1)- [3.5(7) x 10(-4)]T [K] (1300 <= T <= 1830 K, T-g = 960 K, T-m = 1660 K). Comparison to the melt-quenched glass implies an increase in partial derivative n(Tio)/partial derivative T at lower T, as T-g is approached from above. Both Ba-0 coordination and bond length also decrease at higher T, and the role of Ba addition is to reduce n(Tio) below its value in pure molten TiO2, which is related to the presence of density maxima in molten BaO-TiO2. Density measurements made by imaging of the levitated melt yielded rho(T) = 4.82(55)- 0.0004(3)T in units of K and g cm(-3). While BaTi2O5 glass likely consists of a fully connected Ti-0 network, free of nonbridging oxygen (NBO) [OTi1 and with at least 13(4)% [OTi3] triclusters, the 1835(40) K equilibrium melt contains at least 10(4)% NBO along with 90(4)% bridging oxygen [OTi2]. The results highlight the fact that glasses can be considered as structural analogues of melts only for those melts deeply supercooled into the glass transition region. The results imply possible fictive T dependence of titanate glass structure, suggesting applications as, e.g., laser written waveguides with large refractive indices and refractive index contrasts. The temperature-dependent structure further implies a super-Arrhenian melt viscosity with consequences for glass manufacture, titanate-rich slags produced in iron smelting, TiO2-bearing magmas, and by analogy silicate melts at high pressures, as found in planetary interiors.

36 MATERIALS SCIENCE↗
DOE OSTI2020

SAM Plug-in Development (Phase I Final Report)

The DOE Office of Nuclear Energy (NE) has created an extensive set of advanced modeling and simulation tools for nuclear engineering analysis. The advanced capabilities of these newer analysis codes require more in-depth training, skills, and knowledge in order to effectively utilize them for the design, analysis, and licensing of advanced nuclear systems and experiments. A high learning curve for inexperienced users may deter organizations from incorporating these tools into their internal processes. This project involved development of a plug-in to the Symbolic Nuclear Analysis Package (SNAP) for the System Analysis Module (SAM) tool. SAM is an advanced system analysis tool for reactor transient analyses being developed at Argonne National Laboratory under the U.S. DOE Office of Nuclear Energy’s Nuclear Energy Advanced Modeling and Simulation (NEAMS) program. SAM utilizes an object-oriented application framework (MOOSE), and its underlying meshing and finite-element library (libMesh) and linear and non-linear solvers (PETSc), to leverage modern advanced software environments and numerical methods. SNAP provides a highly flexible framework for creating, modifying and documenting input for engineering analysis codes such as SAM as well as extensive functionality for submitting, monitoring, and interacting with the codes through an intuitive graphical user interface (GUI). The common user interface provided by SNAP minimizes the learning curve for engineers starting with a new analysis code and provides an intuitive framework for transitioning between different analysis codes. SNAP provides a powerful but intuitive interface to facilitate access to advanced modeling and simulation tools for inexperienced users. Unlike many “form based” GUI’s, SNAP maps each engineering code’s component input to an internal database which manages all component input parameters along with component interconnections. This level of abstraction permits SNAP to support several advanced capabilities such as renodalization, model validation and consistency checks, embedded documentation, model notebook generation, data ownership and reviewer tracking, and variable assignment for inputs to name a few. SNAP includes a built-in Python interpreter and is interfaced to several commercial and open source packages including CPython, MATLAB/OCTAVE, Microsoft Office, Open Office, and SANDIA’s DAKOTA package which provides Uncertainty Quantification analysis through the SNAP plug-ins. Phase I of this project involved development a fully functional basic SAM plug-in to SNAP. This plug-in provides the ability to import existing models, graphically construct, edit and submit models using SNAP’s extensive functionality.

99 GENERAL AND MISCELLANEOUS↗
NASA NTRS2015

A Self-Stabilizing Hybrid Fault-Tolerant Synchronization Protocol

This paper presents a strategy for solving the Byzantine general problem for self-stabilizing a fully connected network from an arbitrary state and in the presence of any number of faults with various severities including any number of arbitrary (Byzantine) faulty nodes. The strategy consists of two parts: first, converting Byzantine faults into symmetric faults, and second, using a proven symmetric-fault tolerant algorithm to solve the general case of the problem. A protocol (algorithm) is also present that tolerates symmetric faults, provided that there are more good nodes than faulty ones. The solution applies to realizable systems, while allowing for differences in the network elements, provided that the number of arbitrary faults is not more than a third of the network size. The only constraint on the behavior of a node is that the interactions with other nodes are restricted to defined links and interfaces. The solution does not rely on assumptions about the initial state of the system and no central clock nor centrally generated signal, pulse, or message is used. Nodes are anonymous, i.e., they do not have unique identities. A mechanical verification of a proposed protocol is also present. A bounded model of the protocol is verified using the Symbolic Model Verifier (SMV). The model checking effort is focused on verifying correctness of the bounded model of the protocol as well as confirming claims of determinism and linear convergence with respect to the self-stabilization period.

Malekpour, Mahyar R.↗
DOE OSTI2020

Pore-scale simulation of gas flow in microscopic permeable media with complex geometries

Lattice Boltzmann method (LBM) is an efficient tool to perform direct numerical simulation of gas flow in micro-size pores. Despite great advances in LBM and its use in describing gas flow in micro-size channels in slip and transitional flow regimes, exploring formulations that are able to capture flow behavior in domains with complex boundary geometries remains a challenging task. As a result, the impact of complexities in pore structures on gas flow may not be fully explored using LBM. In this study, we propose to use the combined bounce-back and Maxwellian diffusive reflection scheme to capture the slip velocity with a new set of slip coefficients. Optimal values of slip coefficients are determined using a gradient based method, where the results from the linearized Boltzmann equation is used as the reference solution. The proposed formulation is further validated against predictions from molecular dynamics simulation in the presence of complex geometries by introducing obstacles of different shapes in straight channels. LBM with the proposed boundary treatment is then used to investigate gas flow in a synthetic microscopic permeable medium. Here, results indicate that gas exhibits different flow configurations as Knudsen number varies, and apparent gas permeability appears to have up to a second-order dependency on the reciprocal mean pressure in slip and transitional flow regimes.

03 NATURAL GAS↗
DOE OSTI2023

Report of METL Mutual Inductance Level Sensor Development for Use in Liquid Metals – FY2023

A robust mutual inductance level sensor (MILS) for use in high temperature liquid metals has been developed at Argonne National Laboratory (ANL). The current design utilizes mineral insulated cables with a 300-series stainless steel sheath, magnesium oxide insulation, and a single copper conductor. Two coils are wrapped on a common core made of a 300-series stainless steel tube. One coil is energized with an alternating current (AC) source, and this electromagnetically couples with the second coil to generate an induced voltage. The voltage is measured to determine the mutual induction between the coils, and the mutual inductance measurement can be used to determine the level of a nearby liquid metal volume. The MILS can be located in an isolating thimble that is fully sealed to a vessel containing liquid metal for ease of maintenance and replacement. When the primary coil is energized with the AC source, the coil not only electromagnetic (EM) couples with the second coil, but also the surrounding liquid metal. The EM coupling with the surrounding liquid metal reduces the EM coupling with the second coil, producing an inverse-linear relationship between liquid metal level and secondary coil voltage. This has all been demonstrated at the Mechanisms Engineering Test Loop (METL) liquid sodium facility at ANL. The most current iteration of the MILS system is the MILS-MK-II. This sensor system has been commissioned in a non-sodium test stand where calibrations were performed using a sodium analog. The calibrations proved to be highly linear and repeatable. The MILS-MK-II has been installed in the METL expansion tank where high temperature sodium tests have been performed. The MILS-MK-II has been calibrated against a known standard at temperature of 300°C, and the calibrations have proved to be highly linear and repeatable. The calibrated MILS-MK-II has been in operation in the METL facility for several thousands of hours at temperatures around 300°C. The experimental data has been used to develop and validate electromagnetic finite element models in COMSOL Multiphysics, and now these models can be used to advance the development of the sensor system. Next steps will include temperature compensation to allow for operation at various temperatures. Additionally, efforts to incorporate a self-calibration methodology are underway.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2015

Controls on terrestrial carbon feedbacks by productivity versus turnover in the CMIP5 Earth System Models

To better understand sources of uncertainty in projections of terrestrial carbon cycle feedbacks, we present an approach to separate the controls on modeled carbon changes. We separate carbon changes into four categories using a linearized, equilibrium approach: those arising from changed inputs (productivity-driven changes), and outputs (turnover-driven changes), of both the live and dead carbon pools. Using Coupled Model Intercomparison Project Phase 5 (CMIP5) simulations for five models, we find that changes to the live pools are primarily explained by productivity-driven changes, with only one model showing large compensating changes to live carbon turnover times. For dead carbon pools, the situation is more complex as all models predict a large reduction in turnover times in response to increases in productivity. This response arises from the common representation of a broad spectrum of decomposition turnover times via a multi-pool approach, in which flux-weighted turnover times are faster than mass-weighted turnover times. This leads to a shift in the distribution of carbon among dead pools in response to changes in inputs, and therefore a transient but long-lived reduction in turnover times. Since this behavior, a reduction in inferred turnover times resulting from an increase in inputs, is superficially similar to priming processes, but occurring without the mechanisms responsible for priming, we call the phenomenon "false priming", and show that it masks much of the intrinsic changes to dead carbon turnover times as a result of changing climate. These patterns hold across the fully coupled, biogeochemically coupled, and radiatively coupled 1 % yr -1 increasing CO 2 experiments. We disaggregate inter-model uncertainty in the globally integrated equilibrium carbon responses to initial turnover times, initial productivity, fractional changes in turnover, and fractional changes in productivity. For both the live and dead carbon pools, inter-model spread in carbon changes arising from initial conditions is dominated by model disagreement on turnover times, whereas inter-model spread in carbon changes from fractional changes to these terms is dominated by model disagreement on changes to productivity in response to both warming and CO 2 fertilization. However, the lack of changing turnover time control on carbon responses, for both live and dead carbon pools, in response to the imposed forcings may arise from a common lack of process representation behind changing turnover times (e.g., allocation and mortality for live carbon; permafrost, microbial dynamics, and mineral stabilization for dead carbon), rather than a true estimate of the importance of these processes.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2010

The Effect of Incident Light Polarization on Vegetation Bidirectional Reflectance Factor

The Laboratory-based Bidirectional Reflectance Factor (BRF) polarization study of vegetation is presented in this paper. The BRF was measured using a short-arc Xenon lamp/monochromator assembly producing an incoherent, tunable light source with a well-defined spectral bandpass at visible and near-infrared wavelengths of interest at 470 nm and 870 nm and coherent light source at 1.656 microns. All vegetation samples were measured using P and S linearly polarized incident light over a range of incident and scatter angles. By comparing these results, we quantitatively examine how the BRF of the samples depends on the polarization of the incident light. The differences are significant, depend strongly on the incident and scatter angles, and can be as high as 120% at 67 deg incident and 470nm. The global nature of Earth's processes requires consistent long-term calibration of all instruments involved in data retrieval. The BRF defines the reflection characteristics of Earth surface. It provides the reflectance of a target in a specific direction as a function of illumination and viewing geometry. The BRF is a function of wavelength and reflects the structural and optical properties of the surface. Various space and airborne radiometric and imaging remote sensing instruments are used in the remote sensing characterization of vegetation canopies and soils, oceans, or especially large pollution sources. The satellite data is validated through comparison with airborne, ground-based and laboratory-based data in an effort to fully understand the vegetation canopy reflectance, The Sun's light is assumed to be unpolarized at the top of the atmosphere; however it becomes polarized to some degree due to atmospheric effects by the time it reaches the vegetation canopy. Although there are numerous atmospheric correction models, laboratory data is needed for model verification and improvement.

Georgiev, Georgi T.↗
DOE OSTI2015

The spherically symmetric gravitational collapse of a clump of solids in a gas

In the subject of planetesimal formation, several mechanisms have been identified that create dense particle clumps in the solar nebula. The present work is concerned with the gravitational collapse of such clumps, idealized as being spherically symmetric. Fully nonlinear simulations using the two-fluid model are carried out (almost) up to the time when a central density singularity forms. We refer to this as the collapse time. The end result of the study is a parametrization of the collapse time, in order that it may be compared with timescales for various disruptive effects to which clumps may be subject in a particular situation. An important effect that determines the collapse time is that as the clump compresses, it also compresses the gas due to drag. This increases gas pressure, which retards particle collapse and can lead to oscillation in the size and density of the clump. In the limit of particles perfectly coupled to the gas, the characteristic ratio of gravitational force to gas pressure becomes relevant and defines a two-phase Jeans parameter, J{sub t}, which is the classical Jeans parameter with the speed of sound replaced by an effective wave speed in the coupled two-fluid medium. The parameter J{sub t} remains useful even away from the perfect coupling limit because it makes the simulation results insensitive to the initial density ratio of particles to gas (Φ{sub 0}) as a separate parameter. A simple ordinary differential equation model is developed. It takes the form of two coupled non-linear oscillators and reproduces key features of the simulations. Finally, a parametric study of the time to collapse is performed and a formula (fit to the simulations) is developed. In the incompressible limit J{sub t}→0, collapse time equals the self-sedimentation time, which is inversely proportional to the Stokes number. As J{sub t} increases, the collapse time decreases with J{sub t} and eventually becomes approximately equal to the dynamical time. Values of collapse time versus clump size are given for a minimum-mass solar nebula. Finally, the timescale of clump erosion due to turbulent strain is estimated.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Exploiting non-linear scales in galaxy–galaxy lensing and galaxy clustering: A forecast for the dark energy survey

ABSTRACT The combination of galaxy–galaxy lensing (GGL) and galaxy clustering is a powerful probe of low-redshift matter clustering, especially if it is extended to the non-linear regime. To this end, we use an N-body and halo occupation distribution (HOD) emulator method to model the redMaGiC sample of colour-selected passive galaxies in the Dark Energy Survey (DES), adding parameters that describe central galaxy incompleteness, galaxy assembly bias, and a scale-independent multiplicative lensing bias Alens. We use this emulator to forecast cosmological constraints attainable from the GGL surface density profile ΔΣ(rp) and the projected galaxy correlation function wp, gg(rp) in the final (Year 6) DES data set over scales $r_p=0.3\!-\!30.0\, h^{-1} \, \mathrm{Mpc}$. For a $3{{\ \rm per\ cent}}$ prior on Alens we forecast precisions of $1.9{{\ \rm per\ cent}}$, $2.0{{\ \rm per\ cent}}$, and $1.9{{\ \rm per\ cent}}$ on Ωm, σ8, and $S_8 \equiv \sigma _8\Omega _m^{0.5}$, marginalized over all halo occupation distribution (HOD) parameters as well as Alens. Adding scales $r_p=0.3\!-\!3.0\, h^{-1} \, \mathrm{Mpc}$ improves the S8 precision by a factor of ∼1.6 relative to a large scale ($3.0\!-\!30.0\, h^{-1} \, \mathrm{Mpc}$) analysis, equivalent to increasing the survey area by a factor of ∼2.6. Sharpening the Alens prior to $1{{\ \rm per\ cent}}$ further improves the S8 precision to $1.1{{\ \rm per\ cent}}$, and it amplifies the gain from including non-linear scales. Our emulator achieves per cent-level accuracy similar to the projected DES statistical uncertainties, demonstrating the feasibility of a fully non-linear analysis. Obtaining precise parameter constraints from multiple galaxy types and from measurements that span linear and non-linear clustering offers many opportunities for internal cross-checks, which can diagnose systematics and demonstrate the robustness of cosmological results.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2016

An Efficient Single Frequency Ho:YLF Laser for IPDA Lidar Applications

A highly efficient, versatile, single frequency 2-micron pulsed laser can be used in a pulsed Differential Absorption Lidar (DIAL) / Integrated Path Differential Absorption (IPDA) instrument to make precise, high-resolution measurements to investigate sources and sinks of CO2. For a direct detection IPDA lidar, the desired 2 m Ho:YLF laser should generate 30-40 mJ pulses at the repetition rate of 100 to 200 Hz, with short pulse length (<100 ns) and better than 2% wall plug efficiency. A Tm fiber laser in-band pumped Ho:YLF laser has been developed to meet this technical challenge. This Ho:YLF laser is designed in a four mirror ring resonator with bow tie configuration, which helps to obtain high beam quality. It is end-pumped by a 40 W linearly polarized Tm fiber laser at 1.94μm. The resonator length is 1.10 meters with output coupler reflectivity at 45%. The laser crystal size is 3 x 3 x 60 mm (w, h, l) with a doping concentration of 0.5% Holmium. The laser beam and pump beam are mode-matched in the active medium. Thus, the pump and laser beams have the same confocal parameters. Mode-matching is also helpful for operating the laser in a single transverse mode. The laser beam waist is slightly less than 0.5 mm at the center of the laser crystal. Based on quasi-four level modeling, pump absorption and saturation depend on laser intensity. Laser amplification and saturation also depend on the pump intensity in the crystal. The laser is injection seeded to obtain the single frequency required by an IPDA lidar measurement. The seed beam is entered into the resonator through an output coupler. The laser is mounted on a water cooled optical bench for stable and reliable operation. The size of the optical bench is 22.16 x 9.20 x 1.25 inches. It is stiffened so that the laser can be operated in any orientation of the optical bench. This packaged Ho:YLF laser is designed for either mobile trailer or airborne platform operation. The engineering prototype Ho:YLF laser has been fully characterized to demonstrate laser performance. Figure 1 shows the laser output power as a function of pump power at different pulse repetition rates from 100 Hz to 333 Hz. The threshold of the laser is less than 14 W. The slope efficiencies are 28%, 40%, 41% and 43% for pulse repetition rates of 100, 200, 250 and 333 Hz, respectively. Maximum power increases with the pulse repetition rate. Output power of 4.2 W, 6 W, 6.7 W, and 7.7 W is achieved for pulse repetition rates of 100, 200, 250 and 333Hz, respectively. This represents the optical conversion efficiency of 16.7%, 22.4%, 23.7%, and 26.5% at these various pulse repletion rates. It is the most efficient and compact Ho:YLF laser demonstrated in the high pulse energy (>20mJ) and moderate pulse repetition rate (100's Hz) operation range. As shown in Figure 1, the maximum pulse energy at 100 Hz is 42 mJ. This is limited due to optical damage. The laser stability is characterized and found to be very stable. A relative pulse energy standard deviation of 2% was measured. The beam quality of the Ho:YLF was measured by a Spiricon infrared laser beam camera. Figure 2 shows the beam profile image of the laser. Both the X-profile of the beam (horizontal direction) and the Y-profile of the beam (vertical direction) are well fitted by a Gaussian profile. The qualitative beam quality measurement shows excellent beam quality in both axis. The M-square value for the laser beam is measured at 1.06 and 1.09 for the x and y axis respectively.-

Yu, J.↗
DOE OSTI2021

A new and alternative look at nonlinear Alfvénic states

The formulation for studying nonlinear Alfvenic states, sustainable in Hall Magnetohydrodynamics (HMHD), becomes considerably simpler and more tractable when circularly polarized Beltrami vectors (the eigenstates of linear HMHD) are used as the basis functions. Nonlinear HMHD is, then, reduced to a rather simple looking set of scalar equations from which a model problem of three interacting Beltrami modes is formulated and analytically solved. The triplet interactions span a variety of familiar nonlinear processes leading to a redistribution as well as periodic exchange of energy. The energy exchange processes (whose strength is measured by an energy exchange/depletion time) will, perhaps, play a dominant role in determining the spectral content of an eventual Alfvenic state. All nonlinearities (sensitive functions of the interacting wave vectors) operate at par, and none is dominant over any substantial region of k-space; their intricate interplay prevents a “universal” picture from emerging; few generalizations on the processes that may, for instance, lead to a turbulent state, are possible. However, the theory can definitely claim: (1) the energy tends to flow from lower to higher k and (2) the higher k z (in the direction of the ambient magnetic field) components of a mode with a given k are depleted/oscillate faster—in some cases much faster. It is noteworthy that the mode coupling is the strongest (with the shortest depletion time) when the participating wave vectors are nearly perpendicular; perhaps, an expected consequence of the curl (cross product) nonlinearities. Numerical simulations will be necessary to help create a fully reliable picture.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Inter-Annual Precipitation Variability Decreases Switchgrass Productivity from Arid to Mesic Environments

Cellulosic biofuels are an important source of renewable biomass within the alternative energy portfolio. Switchgrass (Panicum virgatum L.), a perennial C 4 grass native to North America, is widely studied as a biofuel feedstock for its consistently high yields and minimal input requirements. The influences of precipitation amount and temporal variability on the fertilizer response of switchgrass productivity are not fully understood. Moreover, global climate models predict changes in rainfall patterns towards lower and increasingly variable soil water availability in several productive areas worldwide, which may impact net primary production of biofuel crops. We conducted a meta-analysis of aboveground net primary production of switchgrass from 48 publications encompassing 82 different locations, 11 soil types, 52 switchgrass cultivars, fertilizer inputs between 0 to 896 kg N ha –1 year –1 , and 1 to 6 years of annual productivity measures repeated on the same stand. Productivity of the lowland ecotype doubled with N rates > 131 kg N ha –1 year –1 , but upland ecotype productivity increased only by 50%. Results showed an optimum N rate of 30 to 60 kg N ha –1 year –1 for both ecotypes, after which biomass gain per unit of N added decreased. Growing season precipitation (GSPPT) and inter-annual precipitation variability (inter-PPTvar) affected both ecotypes similarly. Long-term mean annual precipitation (MAP) differentially affected lowland and upland productivity, depending on the N level. Productivity responses to MAP and GSPPT were similar for both upland and lowland ecotypes at none or low N rates. When N increased beyond 60 kg N ha –1 year –1 , lowland cultivars had a greater growth response to MAP than uplands. Productivity increased with increasing GSPPT and MAP and had a positive linear response to MAP ranging from 600 to 1200 mm year –1 . One third of the variability in switchgrass production was accounted for by inter-PPTvar. After accounting for MAP, sites with higher inter-PPTvar had lower switchgrass productivity than sites with lower inter-PPTvar. Increased inter-annual variation in precipitation reduced production of both ecotypes. As a result, predicted changes in the amount and timing of precipitation thus likely will exert greater influence on production of upland than lowland ecotypes of switchgrass.

54 ENVIRONMENTAL SCIENCES↗