Search NASA⌕ Search

SEARCH · Search NASA

Results for “coupling coefficient”

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 1,477 records · Page 82
All sourcesNASA NTRSDOE OSTI
NASA NTRS2012

Assimilation of GIRO Data into a Real-Time IRI

Increasingly accurate and detailed global 3-D specification of the Earth's space plasma environment is required to further understand its intricate organization and behavior. For a long time space physics and aeronomy research has been data starved due to the great variety of natural time scales involved in the plasma phenomenology. We have started developing a new approach to the global ionospheric specification called Real-Time Assimilative Mapping (RTAM). The IRI-RTAM will use data from the Global Ionospheric Radio Observatory (GIRO) to smoothly transform International Reference Ionosphere's (IRI's) background empirical maps of the ionospheric characteristics to match the observations. Such empirical assimilative modeling will provide a high-resolution, global picture of the ionospheric response to various short-term events observed during periods of storm activity or the impact of gravity waves coupling the ionosphere to the lower atmosphere, including timelines of the vertical restructuring of the plasma distribution. It will also contribute to the challenging task of providing a rapid insight into the temporal and spatial space weather development using the real-time GIRO data streams. The new assimilation technique "updates" the IRI electron density distribution while preserving the overall integrity of IRI s typical ionospheric feature representations. The technique adjusts the coefficients of the spherical/diurnal expansions used by the CCIR and URSI-88 model to obtain the global sub-peak electron density distribution. The set of global corrected coefficients can be generated as frequently as every 15 min and easily disseminated using a single real-time RTAM server operated by GIRO.

Galkin, I. A.↗
DOE OSTI2017

Learning from Higgs physics at future Higgs factories

Future Higgs factories can reach impressive precision on Higgs property measurements. In this paper, instead of conventional focus of Higgs precision in certain interaction bases, we explored its sensitivity to new physics models at the electron-positron colliders. In particular, we studied two categories of new physics models, Standard Model (SM) with a real scalar singlet extension, and Two Higgs Double Model (2HDM) as examples of weakly-interacting models, Minimal Composite Higgs Model (MCHM) and three typical patterns of the more general operator counting for strong interacting models as examples of strong dynamics. We performed a global fit to various Higgs search channels to obtain the 95% C.L. constraints on the model parameter space. In the SM with a singlet extension, we obtained the limits on the singlet-doublet mixing angle sin(theta), as well as the more general Wilson coefficients of the induced higher dimensional operators. In the 2HDM, we analyzed tree level effects in tan(beta) vs. cos(beta-alpha) plane, as well as the one-loop contributions from the heavy Higgs bosons in the alignment limit to obtain the constraints on heavy Higgs masses for different types of 2HDM. In strong dynamics models, we obtained lower limits on the strong dynamics scale. In addition, once deviations of Higgs couplings are observed, they can be used to distinguish different models. Here, we also compared the sensitivity of various future Higgs factories, namely Circular Electron Positron Collider (CEPC), Future Circular Collider (FCC)-ee and International Linear Collider (ILC).

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2017

Giant magnetoelectric effects achieved by tuning spin cone symmetry in Y-type hexaferrites

Multiferroics materials, which exhibit coupled magnetic and ferroelectric properties, have attracted tremendous research interest because of their potential in constructing next-generation multifunctional devices. The application of single-phase multiferroics is currently limited by their usually small magnetoelectric effects. Here, we report the realization of giant magnetoelectric effects in a Y-type hexaferrite Ba 0.4 Sr 1.6 Mg 2 Fe 12 O 22 single crystal, which exhibits record-breaking direct and converse magnetoelectric coefficients and a large electric-field-reversed magnetization. We have uncovered the origin of the giant magnetoelectric effects by a systematic study in the Ba 2-x Sr x Mg 2 Fe 12 O 22 family with magnetization, ferroelectricity and neutron diffraction measurements. With the transverse spin cone symmetry restricted to be two-fold, the one-step sharp magnetization reversal is realized and giant magnetoelectric coefficients are achieved. Our study reveals that tuning magnetic symmetry is an effective route to enhance the magnetoelectric effects also in multiferroic hexaferrites.

36 MATERIALS SCIENCE↗
NASA NTRS1987

Absolute flux measurements for swift atoms

While a torsion balance in vacuum can easily measure the momentum transfer from a gas beam impinging on a surface attached to the balance, this measurement depends on the accommodation coefficients of the atoms with the surface and the distribution of the recoil. A torsion balance is described for making absolute flux measurements independent of recoil effects. The torsion balance is a conventional taut suspension wire design and the Young modulus of the wire determines the relationship between the displacement and the applied torque. A compensating magnetic field is applied to maintain zero displacement and provide critical damping. The unique feature is to couple the impinging gas beam to the torsion balance via a Wood's horn, i.e., a thin wall tube with a gradual 90 deg bend. Just as light is trapped in a Wood's horn by specular reflection from the curved surfaces, the gas beam diffuses through the tube. Instead of trapping the beam, the end of the tube is open so that the atoms exit the tube at 90 deg to their original direction. Therefore, all of the forward momentum of the gas beam is transferred to the torsion balance independent of the angle of reflection from the surfaces inside the tube.

Fink, M.↗
DOE OSTI2022

Reduced diffusion and enhanced retention of multiple radionuclides from pore structure characterization of barrier materials for enhanced repository performance

Fluid flow and chemical transport in porous media are the macroscopic consequences of pore structure, which integrates geometry (e.g., pore size and surface area, pore-size distribution) and topology (e.g., pore connectivity). Low-permeability geological media whose pores are poorly interconnected will exhibit the characteristics of anomalous diffusion and sample size-dependent effective porosity, which will strongly impact long-term net diffusion and retention of radionuclides in geological repository settings involving different host rocks and barrier materials. A suite of innovative and complementary experimental approaches is utilized to study the microscopic pore structure and macroscopic fluid flow & chemical transport for a range of host rocks and barrier materials, in addition to standard clay minerals and reference rocks. With a particular focus on quantifying the presence and magnitude of “isolated” pores for a reduced effective porosity in low-permeability geomedia, the integrated methodologies for basic properties and pore structure characterization of these geomedia include X-ray diffraction, thin section petrography, grain size distribution, water immersion porosimetry after vacuum-pulling for full saturation, mercury intrusion porosimetry, nitrogen physisorption, scanning electron microscopy, X-ray computed tomography, and (ultra-)small angle neutron (X-ray) scattering. In addition, custom-designed gas diffusion, tracer recipe involving a range of anionic and cationic chemicals with subsequent analyses by laser ablation and inductively coupled plasma-mass spectrometry, along with batch sorption, column transport, and imbibition tests were conducted for coupled effects of pore structure and chemical retention/transport. From the perspectives of pore structure in conjunction with multiple and complementary approaches to examining a range of sample sizes under different observational scales, we find that the poor pore connectivity is prevalent in low-permeability media (mudstone and crystalline rock) that is related to geological processes (e.g., compaction, diagenesis and thermal maturation). For example, the deep and organic matter-rich mudstones have a much smaller effective porosity than the total porosity (as a result of poor pore connectivity) and associated diffusion coefficient, and the effective porosity & diffusion coefficients are also dependent upon the sample sizes used in the measurement. Similarly, most of the pore space in the shallow mudstone is also controlled by pore-throat diameters in the 5-50 nm range of intergranular pore types from its fine-grained nature, but with an overall good pore connectivity. However, the nm-sized pore space (physically pore-network architecture) and strong sorption capacities (chemical retention from clay minerals) of both shallow and deep mudstones lead to the synergistic retention of cationic radionuclides and their utilities as effective host rocks and barrier materials. Our unique approaches of studying how the micro-scale pore structure affect macro-scale fluid flow, diffusion & retention, and chemical transport produce improved mechanistic understanding, and realistic quantification, of diffusion and retention of typical radionuclides in a range of generic host rocks and barrier materials (clay/shale, salt, crystalline rock, and tuff), with the overall results leading to scientifically-based understanding of enhanced isolation (from both diffusion and retention) of radionuclides and improved confidence on the long-term performance of geological repository to store high-level radioactive wastes. In addition to the training of 25 undergraduates, graduates, and postdocs of UTA, the scientists (organizations) involved in performing this work (e.g., discussion, sample sharing, and operation of SANS and SAXS instruments) include Ed Matteo, Yifeng Wang, and Kristopher Kuhlman (Sandia National Laboratories), Jens Birkholzer, Liange Zheng, Tim Kneafsey, and Sharon Borglin (Lawrence Berkeley National Laboratory), Mavrik Zavarin (Lawrence Livermore National Laboratory), Yukio Tachi and Yuta Fukatsu (Japan Atomic Energy Agency), Mieke de Craen (Euridice, Belgium), Markus Bleuel (NIST), Wei-Ren Chen, Gergely Nagy, Changwoo Do, William Heller, Larry Anovitz, and Kenneth Littrell (ORNL), as well as Jan Illvsky, Ivan Kuzmenko, Ju-Sang Park and Jon Almers (ANL). Key deliverables include a total of 13 peer-reviewed journal articles (nine published and three under review), 23 presentations at scientific conferences (AAPG, AAPG Southwest Section, AGU, Asian Clay Conference, GSA, GSA South-Central Section, IHLRWM, InterPore, International Conference on Chemistry and Migration Behavior of Actinides and Fission Products in the Geosphere, International Conference on Coupled Processes in Fractured Geological Media: Observation, Modeling and Application), and academic institutions (UTA, New Mexico State University; University of Poitiers, France; University of Helsinki, Finland; Uppsala University, Sweden; Istanbul Technical University, Turkey) and other organizations (Andra, France; Posiva Oy, Finland).

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2018

Nonlinear optimization for a tumor invasion PDE model

In this work, we introduce a methodology to approximate unknown parameters that appear on a non-linear reaction–diffusion model of tumor invasion. These equations consider that tumor-induced alteration of micro-environmental pH furnishes a mechanism for cancer invasion. A coupled system reaction–diffusion explaining this model is given by three partial differential equations for the non-dimensional spatial distribution and temporal evolution of the density of normal tissue, the neoplastic tissue growth and the excess concentration of H{sup +} ions. The tumor model parameters have a corresponding biological meaning: the reabsorption rate, the destructive influence of H{sup +} ions in the healthy tissue, the growth rate of tumor tissue and the diffusion coefficient. We propose to solve the direct problem using the Finite Element Method (FEM) and minimize an appropriate functional including both the real data (obtained via in-vitro experiments and fluorescence ratio imaging microscopy) and the numerical solution. The gradient of the functional is computed by the adjoint method.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2018

Structural and Electrochemical Characteristics of Ca-Doped “Flower-like” Li 4 Ti 5 O 12 Motifs as High-Rate Anode Materials for Lithium-Ion Batteries

Doped motifs offer an intriguing structural pathway toward improving conductivity for battery applications. Specifically, Ca-doped, three-dimensional “flower-like” Li 4–x Ca x Ti 5 O 12 (“x” = 0, 0.1, 0.15, and 0.2) micrometer-scale spheres have been successfully prepared for the first time using a simple and reproducible hydrothermal reaction followed by a short calcination process. The products were experimentally characterized by means of X-ray diffraction (XRD), transmission electron microscopy (TEM), scanning electron microscopy (SEM), energy dispersive spectroscopy (EDS) mapping, inductively coupled plasma optical emission spectrometry (ICP-OES), X-ray photoelectron spectroscopy (XPS), cyclic voltammetry (CV), electrochemical impedance spectroscopy (EIS), and galvanostatic charge–discharge testing. Calcium dopant ions were shown to be uniformly distributed within the LTO structure without altering the underlying “flower-like” morphology. The largest lattice expansion and the highest Ti 3+ ratios were noted with XRD and XPS, respectively, whereas increased charge transfer conductivity and decreased Li + -ion diffusion coefficients were displayed in EIS for the Li 4–x Ca x Ti 5 O 12 (“x” = 0.2) sample. The “x” = 0.2 sample yielded a higher rate capability, an excellent reversibility, and a superior cycling stability, delivering 151 and 143 mAh/g under discharge rates of 20C and 40C at cycles 60 and 70, respectively. In addition, a high cycling stability was demonstrated with a capacity retention of 92% after 300 cycles at a very high discharge rate of 20C. In addition, first-principles calculations based on density functional theory (DFT) were conducted with the goal of further elucidating and understanding the nature of the doping mechanism in this study. The DFT calculations not only determined the structure of the Ca-doped Li 4 Ti 5 O 12 , which was found to be in accordance with the experimentally measured XPD pattern, but also yielded valuable insights into the doping-induced effect on both the atomic and electronic structures of Li 4 Ti 5 O 12 .

25 ENERGY STORAGE↗
DOE OSTI2019

Generalised route to effective field theories for quantum systems with local constraints

Some of the exciting phenomena uncovered in strongly correlated systems in recent years – for instance quantum topological order, deconfined quantum criticality and emergent gauge symmetries – appear in systems where the Hilbert space is effectively projected at low energies in a way that imposes local constraints on the original degrees of freedom. Cases in point include spin liquids, valence bond systems, dimer and vertex models. In this work, we use a slave boson description coupled to a large- S path integral formulation to devise a generalised route to obtain effective field theories for such systems. We demonstrate the validity and capability of our approach by studying quantum dimer models and by comparing our results with the existing literature.Field theoretic approaches to date are limited to bipartite lattices, they depend on a gauge-symmetric under-standing of the constraint, and lack generic quantitative predictive power for the coefficients of the terms that appear in the Lagrangians of these systems. Our method overcomes all these shortcomings and we show how the results up to quadratic order compare with the known height description of the square lattice quantum dimer model, as well as with the numerical estimate of the speed of light of the photon excitations on the diamond lattice. Lastly, instanton considerations allow us to infer properties of the finite temperature behaviour in two dimensions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2017

Search for excited states in O 25

Theoretical calculations suggest the presence of low-lying excited states in 25 O. Previous experimental searches by means of proton knockout on 26 F produced no evidence for such excitations. We search for excited states in 25 O using the 24 O(d,p) 25 O reaction. The theoretical analysis of excited states in unbound 25,27 O is based on the configuration interaction approach that accounts for couplings to the scattering continuum. We use invariant-mass spectroscopy to measure neutron-unbound states in 25 O. For the theoretical approach, we use the complex-energy Gamow Shell Model and Density Matrix Renormalization Group method with a finite-range two-body interaction optimized to the bound states and resonances of 23–26 O, assuming a core of 22 O. Here we predict energies, decay widths, and asymptotic normalization coefficients. Our calculations in a large spdf space predict several low-lying excited states in 25 O of positive and negative parity, and we obtain an experimental limit on the relative cross section of a possible J π = 1/2 + state with respect to the ground state of 25 O at σ 1/2+ /σ g.s. = 0.25$^{+1.0}_{-0.25}$. We also discuss how the observation of negative parity states in 25 O could guide the search for the low-lying negative parity states in 27 O. Previous experiments based on the proton knockout of 26 F suffered from the low cross sections for the population of excited states in 25 O because of low spectroscopic factors. In this respect, neutron transfer reactions carry more promise.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Contaminant Transport Modeling for Technology Evaluation and Long-Term Monitoring in the Tims Branch Testbed, SC - 20343

Studies conducted at U.S. DOE sites have shown the presence of heavy metals, radionuclides, and volatile organic compounds in surface water, groundwater, and soil as a result of nuclear activity in the Cold War era. Since the 1990's, innovative cleanup methods have been implemented in the Tims Branch watershed at Savannah River Site (SRS) to limit the contaminant flux to the stream that have reduced the contaminant concentrations to acceptable regulatory levels in the dissolved phase. A tin-based treatment which effectively eliminated all local anthropogenic mercury inputs to this ecosystem resulted in a known step function addition of inert tin oxide particles which now serve as a potential tracer for sedimentation and particle transport processes in the stream. The long-term effectiveness of this and other remediation techniques and the potential for remobilization of adsorbed contaminant in sediment during extreme hydrologic conditions however remains unclear. It is therefore important to understand not only the fate and transport of dissolved contaminants, but also the movement of sediment and the relevant interactions with dissolved contaminant. To narrow this knowledge gap, a study is being conducted using the Tims Branch watershed as a stream-scale ecosystem test-bed to identify the primary transport processes of major contaminants of concern (such as mercury, nickel and uranium) with an emphasis on interactions with sediment transport. This involves the development of a fully distributed hydrologic watershed model of the Tims Branch watershed to predict streamflow under extreme weather conditions, as well as the development of a comprehensive contaminant transport model that can properly account for coupled contaminant and sediment transport. Review of relevant research reports and peer-reviewed journals revealed that aside from advection-dispersion transport of dissolved contaminants, adsorption and desorption with suspended solids and bed sediment also play an important role in the transport of those contaminants of concern. To develop the fully distributed hydrologic watershed model, the MIKE SHE 2-dimensional (2D) land surface/3D groundwater model that simulates surface/subsurface hydrologic processes (such as overland flow, evapotranspiration, and infiltration) was coupled with a 1D streamflow model that accounts for stream water hydraulics (such as hydraulic structures, cross-sections, and network). To model contaminant transport, the MIKE 11 streamflow component was coupled with the MIKE 11 AD module that simulates solute transport through advection and dispersion, and MIKE ECO Lab module that accounts for both sediment transport and interactions with dissolved contaminant. At this stage, the development and optimization of the fully distributed hydrologic model has been completed, achieving satisfactory statistical results between observed and predicted discharge as indicated by a root mean square error (RMSE) of 0.039 cms and a Nash-Sutcliffe efficiency coefficient (NSE) of 0.764. The ongoing development of the contaminant transport model has also yielded realistic results from preliminary tests. Results from this study are a key to evaluating the effectiveness of tin (II)-based mercury treatment of wetlands at the SRS site, and are also relevant to evaluating the potential of using this type of novel remediation technology in other mercury-contaminated stream systems. Knowledge acquired from this research will also support interpretation of historical data on the trends of contaminant concentration distribution in Tims Branch, particularly considering the effect of extreme hydrological events on the stream flow and pollutant transport. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
NASA NTRS1975

Effect of wing planform and canard location and geometry on the longitudinal aerodynamic characteristics of a close-coupled canard wing model at subsonic speeds

A generalized wind-tunnel model with canard and wing planforms typical of highly maneuverable aircraft was tested in the Langley 7- by 10-foot high-speed tunnel at a Mach number of 0.30 to determine the effect of canard location, canard size, wing sweep, and canard strake on canard-wing interference to high angles of attack. The major results of this investigation may be summarized as follows: the high-canard configuration (excluding the canard strake and canard flap), for both the 60 deg and 44 deg swept leading-edge wings, produced the highest maximum lift coefficient and the most linear pitching-moment curves; substantially larger gains in the canard lift and total lift were obtained by adding a strake to the canard located below the wing chord plane rather than by adding a strake to the canard located above the wing chord plane.

Gloss, B. B.↗
DOE OSTI2020

Piezoresponse amplitude and phase quantified for electromechanical characterization

Piezoresponse force microscopy (PFM) is a powerful characterization technique to readily image and manipulate the ferroelectric domains. PFM gives an insight into the strength of local piezoelectric coupling and polarization direction through PFM amplitude and phase, respectively. Converting measured arbitrary units into units of effective piezoelectric constant remains a challenge, and insufficient methods are often used. While most quantification efforts have been spent on quantifying the PFM amplitude signal, little attention has been given to the PFM phase, which is often arbitrarily adjusted to fit expectations. This is problematic when investigating materials with unknown or negative sign of the probed effective electrostrictive coefficient or strong frequency dispersion of electromechanical responses, because assumptions about the PFM phase cannot be reliably made. The PFM phase can, however, provide important information on the polarization orientation and the sign of the effective electrostrictive coefficient probed by PFM. Most notably, the orientation of the PFM hysteresis loop is determined by the PFM phase. Moreover, when presenting PFM data as a combined signal, the resulting response can be artificially lowered or asymmetric if the phase data have not been correctly processed. Here, we explain the PFM amplitude quantification process and demonstrate a path to identify the phase offset required to extract correct meaning from the PFM phase data. We explore different sources of phase offsets including the experimental setup, instrumental contributions, and data analysis. We discuss the physical working principles of PFM and develop a strategy to extract physical meaning from the PFM amplitude and phase.

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

Toward Establishing Uniqueness of Experimentally Determined Transference Numbers

The passage of current through a battery results in the development of concentration gradients in the electrolytic phase. For a fully characterized binary electrolyte, where the conductivity, salt diffusion coefficient, cation transference number, and the thermodynamic factor are known, concentration and potential gradients in the electrolytic phase can be modeled using Newman’s concentrated solution theory. We report two methods for measuring the transference number: the standard method based on electrochemical measurements ( t + , echem 0 ) and electrophoretic NMR ( t + , eNMR 0 ). The electrochemical approach requires combining measurements from multiple experiments; the equations used to determine the cation transference number and the thermodynamic factor are coupled, nonlinear algebraic equations. In the electrophoretic-NMR-based approach, however, the equations used to determine the cation transference number and the thermodynamic factor are decoupled. We find for a liquid electrolyte comprised of a lithium salt dissolved in tetraglyme, the values of the transference numbers obtained by these two methods are distinct. For example, at 30 °C, t + , echem 0 = −1.02 ± 1.11 and t + , eNMR 0 = 0.25 ± 0.04. The corresponding thermodynamic factors are also different. While the magnitude of the predicted concentration gradients based on the two sets of parameters are different, the predicted current-voltage relationships are similar.

Hickson, Darby T. (ORCID:0000000251339755)↗
NASA NTRS1973

The role of eddy turbulence in the development of self-consistent models of the lower and upper thermosphere

Numerical solutions of mutually coupled time dependent equations of continuity, momentum and energy balance are presented to illustrate the effect of eddy turbulence on the neutral composition and temperature of the lower and upper atmosphere. The procedure adopted eliminates the necessity of making ad-hoc assumptions about the conditions at the turbopause level and allows the study of the two regions in a self-consistent manner. From the illustrative examples comprising parametric changes in the eddy diffusion coefficient, the specific roles of eddy turbulence in the development of theoretical models of the thermosphere are discussed.

Chandra, S.↗
NASA NTRS1986

Analytic computation of energy derivatives - Relationships among partial derivatives of a variationally determined function

Formulas are presented relating Taylor series expansion coefficients of three functions of several variables, the energy of the trial wave function (W), the energy computed using the optimized variational wave function (E), and the response function (lambda), under certain conditions. Partial derivatives of lambda are obtained through solution of a recursive system of linear equations, and solution through order n yields derivatives of E through order 2n + 1, extending Puley's application of Wigner's 2n + 1 rule to partial derivatives in couple perturbation theory. An examination of numerical accuracy shows that the usual two-term second derivative formula is less stable than an alternative four-term formula, and that previous claims that energy derivatives are stationary properties of the wave function are fallacious. The results have application to quantum theoretical methods for the computation of derivative properties such as infrared frequencies and intensities.

King, H. F.↗
NASA NTRS2015

Investigation of North American Vegetation Variability under Recent Climate: A Study Using the SSiB4/TRIFFID Biophysical/Dynamic Vegetation Model

Recent studies have shown that current dynamic vegetation models have serious weaknesses in reproducing the observed vegetation dynamics and contribute to bias in climate simulations. This study intends to identify the major factors that underlie the connections between vegetation dynamics and climate variability and investigates vegetation spatial distribution and temporal variability at seasonal to decadal scales over North America (NA) to assess a 2-D biophysical model/dynamic vegetation model's (Simplified Simple Biosphere Model version 4, coupled with the Top-down Representation of Interactive Foliage and Flora Including Dynamics Model (SSiB4/TRIFFID)) ability to simulate these characteristics for the past 60 years (1948 through 2008). Satellite data are employed as constraints for the study and to compare the relationships between vegetation and climate from the observational and the simulation data sets. Trends in NA vegetation over this period are examined. The optimum temperature for photosynthesis, leaf drop threshold temperatures, and competition coefficients in the Lotka-Volterra equation, which describes the population dynamics of species competing for some common resource, have been identified as having major impacts on vegetation spatial distribution and obtaining proper initial vegetation conditions in SSiB4/TRIFFID. The finding that vegetation competition coefficients significantly affect vegetation distribution suggests the importance of including biotic effects in dynamical vegetation modeling. The improved SSiB4/TRIFFID can reproduce the main features of the NA distributions of dominant vegetation types, the vegetation fraction, and leaf area index (LAI), including its seasonal, interannual, and decadal variabilities. The simulated NA LAI also shows a general increasing trend after the 1970s in responding to warming. Both simulation and satellite observations reveal that LAI increased substantially in the southeastern U.S. starting from the 1980s. The effects of the severe drought during 1987-1992 and the last decade in the southwestern U.S. on vegetation are also evident from decreases in the simulated and satellite-derived LAIs. Both simulated and satellite-derived LAIs have the strongest correlations with air temperature at northern middle to high latitudes in spring reflecting the effect of these climatic variables on photosynthesis and phenological processes. Meanwhile, in southwestern dry lands, negative correlations appear due to the heat and moisture stress there during the summer. Furthermore, there are also positive correlations between soil wetness and LAI, which increases from spring to summer. The present study shows both the current improvements and remaining weaknesses in dynamical vegetation models. It also highlights large continental-scale variations that have occurred in NA vegetation over the past six decades and their potential relations to climate. With more observational data availability, more studies with differentmodels and focusing on different regions will be possible and are necessary to achieve comprehensive understanding of the vegetation dynamics and climate interactions.

climate variability↗
NASA NTRS2015

The Effect of Resistance on Rocket Injector Acoustics

Combustion instability, where unsteady heat release couples with acoustic modes, has long been an area of concern in liquid rocket engines. Accurate modeling of the acoustic normal modes of the combustion chamber is important to understanding and preventing combustion instability. This study evaluates the effect of injector resistance on the mode shapes and complex eigen-frequencies of an injector/combustion chamber system by defining a high Mach-flow form of the convective wave equation (see Eq. 1) in COMSOL Multiphysics' Coefficient Form PDE Mathematics Module.

Morgan, C. J.↗
NASA NTRS2015

Solid Rocket Motor Combustion Instability Modeling in COMSOL Multiphysics

Combustion instability modeling of Solid Rocket Motors (SRM) remains a topic of active research. Many rockets display violent fluctuations in pressure, velocity, and temperature originating from the complex interactions between the combustion process, acoustics, and steady-state gas dynamics. Recent advances in defining the energy transport of disturbances within steady flow-fields have been applied by combustion stability modelers to improve the analysis framework [1, 2, 3]. Employing this more accurate global energy balance requires a higher fidelity model of the SRM flow-field and acoustic mode shapes. The current industry standard analysis tool utilizes a one dimensional analysis of the time dependent fluid dynamics along with a quasi-three dimensional propellant grain regression model to determine the SRM ballistics. The code then couples with another application that calculates the eigenvalues of the one dimensional homogenous wave equation. The mean flow parameters and acoustic normal modes are coupled to evaluate the stability theory developed and popularized by Culick [4, 5]. The assumption of a linear, non-dissipative wave in a quiescent fluid remains valid while acoustic amplitudes are small and local gas velocities stay below Mach 0.2. The current study employs the COMSOL multiphysics finite element framework to model the steady flow-field parameters and acoustic normal modes of a generic SRM. The study requires one way coupling of the CFD High Mach Number Flow (HMNF) and mathematics module. The HMNF module evaluates the gas flow inside of a SRM using St. Robert's law to model the solid propellant burn rate, no slip boundary conditions, and the hybrid outflow condition. Results from the HMNF model are verified by comparing the pertinent ballistics parameters with the industry standard code outputs (i.e. pressure drop, thrust, ect.). These results are then used by the coefficient form of the mathematics module to determine the complex eigenvalues of the Acoustic Velocity Potential Equation (AVPE). The mathematics model is truncated at the nozzle sonic line, where a zero flux boundary condition is self-satisfying. The remaining boundaries are modeled with a zero flux boundary condition, assuming zero acoustic absorption on all surfaces. The results of the steady-state CFD and AVPE analyses are used to calculate the linear acoustic growth rate as is defined by Flandro and Jacob [2, 3]. In order to verify the process implemented within COMSOL we first employ the Culick theory and compare the results with the industry standard. After the process is verified, the Flandro/Jacob energy balance theory is employed and results displayed.

Fischbach, Sean R.↗