Search NASA⌕ Search

SEARCH · Search NASA

Results for “Numerical Relativity”

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,405 records · Page 78

Analysis of supersonic combustion flow fields with embedded subsonic regions

The viscous characteristic analysis for supersonic chemically reacting flows was extended to include provisions for analyzing embedded subsonic regions. The numerical method developed to analyze this mixed subsonic-supersonic flow fields is described. The boundary conditions are discussed related to the supersonic-subsonic and subsonic-supersonic transition, as well as a heuristic description of several other numerical schemes for analyzing this problem. An analysis of shock waves generated either by pressure mismatch between the injected fluid and surrounding flow or by chemical heat release is also described.

Dash, S.↗

General Potential Theory of Arbitrary Wing Sections

The problem of determining the two dimensional potential flow around wing sections of any shape is examined. The problem is condensed into the compact form of an integral equation capable of yielding numerical solutions by a direct process. An attempt is made to analyze and coordinate the results of earlier studies relating to properties of wing sections. The existing approximate theory of thin wing sections and the Joukowski theory with its numerous generalizations are reduced to special cases of the general theory of arbitrary sections, permitting a clearer perspective of the entire field. The method which permits the determination of the velocity at any point of an arbitrary section and the associated lift and moments is described. The method is also discussed in terms for developing new shapes of preassigned aerodynamical properties.

Theodorsen, T.↗

The axisymmetric jet counterflow problem

A numerical method for computing the flow field of an underexpanded axisymmetric supersonic jet flowing counter to a uniform supersonic stream is described. The strong shock waves and shear layers which occur in the flow are treated as normal and tangential discontinuities bounding regions of inviscid flow. Within the regions, the unsteady gasdynamic equations are solved in integral form on a moving computational mesh. As the numerical flow field evolves, the mesh is deformed according to the jump relations applied at the discontinuities, thus determining their positions as part of the solution. Also presented are computed results obtained for a wide range of jet total pressure. The standoff distances of the bow shock and jet shock are in good agreement with those measured experimentally, thus demonstrating the accuracy of the numerical method.

Schiff, L. B.↗

Numerical modeling of the energy storage and release in solar flares

This paper reports on investigation of the photospheric magnetic field-line footpoint motion (usually referred to as shear motion) and magnetic flux emerging from below the surface in relation to energy storage in a solar flare. These causality relationships are demonstrated by using numerical magnetohydrodynamic simulations. From these results, one may conclude that the energy stored in solar flares is in the form of currents. The dynamic process through which these currents reach a critical value is discussed as well as how these currents lead to energy release, such as the explosive events of solar flares.

Wu, S. T.↗

A moment-based kinetic theory model for polydisperse gas–particle flows

Starting from a generalized population balance equation and the Boltzmann–Enskog collision model for hard spheres, a kinetic theory model for polydisperse gas–particle flows is introduced. Here, polydispersity results from spherical particles with the same material density but different diameters. The particle size distribution (PSD) of the particles is allowed to evolve in space and time due to physical processes such as mixing. In order to treat a continuous PSD, the particle-phase model is formulated in terms of the moments of the PSD, and velocity moments conditioned on the particle size. Velocity moments up to second order are included, resulting in transport equations for the mass, momentum and granular temperature, all conditioned on the particle size. In the numerical implementation, the PSD is represented using quadrature-based moment methods (QBMM). With QBMM, a continuous PSD can be treated using a relatively small number of moments as compared to class or sectional methods. In this work, a realizable numerical algorithm for solving the moment system in a finite-volume code is proposed, valid for dilute systems wherein frictional forces are negligible. The ability of the proposed model to describe polydisperse gas–particle systems is demonstrated using cluster-induced turbulence and riser flow.

42 ENGINEERING↗

Factorization of Binary Matrices: Rank Relations, Uniqueness and Model Selection of Boolean Decomposition

The application of binary matrices are numerous. Representing a matrix as a mixture of a small collection of latent vectors via low-rank decomposition is often seen as an advantageous method to interpret and analyze data. In this work, we examine the factorizations of binary matrices using standard arithmetic (real and nonnegative) and logical operations (Boolean and $\mathbb{Z}$ 2 ). We examine the relationships between the different ranks, and discuss when factorization is unique. In particular, we characterize when a Boolean factorization X = W$\land$H has a unique W, a unique H (for a fixed W), and when both W and H are unique, given a rank constraint. We introduce a method for robust Boolean model selection, called BMFk, and show on numerical examples that BMFk not only accurately determines the correct number of Boolean latent features but reconstruct the pre-determined factors accurately.

97 MATHEMATICS AND COMPUTING↗

Never the Two Shall Mix: Robust Indel Markers to Ensure the Fidelity of Two Pivotal and Closely-Related Accessions of Brachypodium distachyon

Brachypodium distachyon is an established model for monocotyledonous plants. Numerous markers intended for gene discovery and population genetics have been designed. However to date, very few indel markers with larger and easily scored length polymorphism differences, that distinguish between the two morphologically similar and highly utilized B. distachyon accessions, Bd21, the reference genome accession, and Bd21-3, the transformation-optimal accession, are publically available. In this study, 22 indel markers were designed and utilized to produce length polymorphism differences of 150 bp or more, for easy discrimination between Bd21 and Bd21-3. When tested on four other B. distachyon accessions, one case of multiallelism was observed. It was also shown that the markers could be used to determine homozygosity and heterozygosity at specific loci in a Bd21 x Bd3-1 F2 population. The work done in this study allows researchers to maintain the fidelity of Bd21 and Bd21-3 stocks for both transgenic and nontransgenic studies. It also provides markers that can be utilized in conjunction with others already available for further research on population genetics, gene discovery and gene characterization, all of which are necessary for the relevance of B. distachyon as a model species.

59 BASIC BIOLOGICAL SCIENCES↗

What Controls the Entrainment Rate of Dry Buoyant Thermals with Varying Initial Aspect Ratio?

Here, this study uses theory and numerical simulations to analyze the nondimensional spreading rate a (change in radius with height) of buoyant thermals as they rise and entrain surrounding environmental fluid. A focus is on how a varies with initial thermal aspect ratio A r , defined as height divided by width of the initial buoyancy perturbation. An analytic equation for thermal ascent rate w t that depends on a is derived from the thermal-volume-averaged momentum bud get equation. The thermal top height when w t is maximum, defining a critical height z c , is inversely proportional to a. The height z c also corresponds to the thermal top height when buoyant fluid along the thermal’s vertical axis is fully replaced by entrained nonbuoyant environmental fluid rising from below the thermal. The time scale for this process is controlled by the vertical velocity of parcels rising upward through the thermal’s core. This parcel vertical velocity is approximated from Hill’s analytic spherical vortex, yielding an analytic inverse relation between a and A r . Physically, this a–A r relation is connected to changes in circulation as A r is modified. Numerical simulations of thermals with A r varied from 0.5 to 2 give a values close to the analytic theoretical relation, with a factor of ~3 decrease in a as A r is increased from 0.5 to 2. The theory also explains why a of initially spherical thermals from past laboratory and modeling studies is about 0.15. Overall, this study provides a theoretical underpinning for understanding the entrainment behavior of thermals, relevant to buoyantly driven atmospheric flows.

54 ENVIRONMENTAL SCIENCES↗

Gas Flow Tightly Coupled to Elastoplastic Geomechanics for Tight- and Shale-Gas Reservoirs: Material Failure and Enhanced Permeability

We investigate coupled flow and geomechanics in gas production from extremely low permeability reservoirs such as tight and shale gas reservoirs, using dynamic porosity and permeability during numerical simulation. In particular, we take the intrinsic permeability as a step function of the status of material failure, and the permeability is updated every time step. We consider gas reservoirs with the vertical and horizontal primary fractures, employing the single and dynamic double porosity (dual continuum) models. We modify the multiple porosity constitutive relations for modeling the double porous continua for flow and geomechanics. The numerical results indicate that production of gas causes redistribution of the effective stress fields, increasing the effective shear stress and resulting in plasticity. Shear failure occurs not only near the fracture tips but also away from the primary fractures, which indicates generation of secondary fractures. These secondary fractures increase the permeability significantly, and change the flow pattern, which in turn causes a change in distribution of geomechanical variables. From various numerical tests, we find that shear failure is enhanced by a large pressure drop at the production well, high Biot's coefficient, low frictional and dilation angles. Smaller spacing between the horizontal wells also contributes to faster secondary fracturing. When the dynamic double porosity model is used, we observe a faster evolution of the enhanced permeability areas than that obtained from the single porosity model, mainly due to a higher permeability of the fractures in the double porosity model. These complicated physics for stress sensitive reservoirs cannot properly be captured by the uncoupled or flow-only simulation, and thus tightly coupled flow and geomechanical models are highly recommended to accurately describe the reservoir behavior during gas production in tight and shale gas reservoirs and to smartly design production scenarios.

02 PETROLEUM↗

Validation of Quasi-Invariant Ice Cloud Radiative Quantities with MODIS Satellite-Based Cloud Property Retrievals

Similarity relations applied to ice cloud radiance calculations are theoretically analyzed and numerically validated. If t(1v) and t(1vg) are conserved where t is optical thickness, v the single-scattering albedo, and g the asymmetry factor, it is possible that substantially different phase functions may give rise to similar radiances in both conservative and non-conservative scattering cases, particularly in the case of large optical thicknesses. In addition to theoretical analysis, this study uses operational ice cloud optical thickness retrievals from the Moderate Resolution Imaging Spectroradiometer (MODIS) Level 2 Collection5 (C5) and Collection 6 (C6) cloud property products to verify radiative similarity relations. It is found that, if the MODIS C5 and C6 ice cloud optical thickness values are multiplied by their respective (1wg)factors, the resultant products referred to as the effective optical thicknesses become similar with their ratio values around unity. Furthermore, the ratios of the C5 and C6 ice cloud effective optical thicknesses display an angular variation pattern similar to that of the corresponding ice cloud phase function ratios. The MODIS C5 and C6 values of ice cloud similarity parameter, defined as [(1w)(1(exp. 1/2)wg)]12, also tend to be similar.

Radiative transfer↗

Evaluating the Effectiveness of an Ultrasonic Acoustic Deterrent in Reducing Bat Fatalities at Wind Energy Facilities

This project was designed to use thermal video cameras and fatality monitoring to evaluate the effectiveness of an ultrasonic acoustic deterrent (UAD) on bat activity and mortality, respectively. Our goals were to redesign the UAD device and installation infrastructure, determine the placement on wind turbines to optimize safety, compatibility and functionality, and to compare the mortality among the following conditions: Control (deterrents off and turbines feathered up to the manufacturer’s cut-in speed of 3.5 m/s), Deterrent (deterrents on and turbines feathered up to the manufacturer’s cut-in speed of 3.5 m/s), Curtailment (deterrents off and turbines feathered up to 5/ m/s), and combination (deterrents on and turbines feathered up to 5 m/s). The project was divided into a Feasibility Study and Comparative Study. The objectives for the Feasibility Study were to develop an installation strategy, redesign a previous iteration of a UAD to improve performance and weatherization, and test the effectiveness of the deterrents on bat activity. The Feasibility Study was intended to work out potential issues using a relatively small number of devices prior to manufacturing and installing numerous devices for a larger-scale comparative study. The division of the project into these separate studies was based on previous experience and the challenges of assessing the capabilities of an untested UAD. During the initial development and manufacturing of the UAD, NRG Systems decided to use a piezoelectric transducer rather than an electrostatic transducer, which was used by a previous device (i.e. Deaton UAD). NRG Systems conducted lab testing (i.e., IP67 or Ingress Protection) to ensure no water or dust ingress. In addition, shock/drop trials and variations in temperature exposure were conducted as part of the reliability testing. NRG Systems also developed a communications system to allow for continuous performance monitoring of the UADs. For the Feasibility Study, we installed 6 UADs on each of 2 Gamesa G90 2-MW wind turbines (Turbine 14 and 8) at the South Chestnut Wind Energy Facility, Pennsylvania and monitored activity under control and treatment (i.e. Deterrent) conditions using thermal video monitoring. There are two major sources of variation in bat activity (beyond the anticipated treatment effect): 1) environment around the turbines might inherently favor more activity at one than the other; 2) weather conditions on any given night or within season difference (e.g., migration later during the study period) might favor more activity on some nights than on other nights. Because we could only monitor two turbines on any night, we sought to control the potential influences of these two sources by alternating the turbine on which deterrents were activated each night. If there were no loss of data due to technical failures, this design would result in an equal number of deterrent and control nights at each turbine through the monitoring period, balancing the effects of both sources of variation. We compared the time bats spent and the number of events that occurred in overlapping cameras FOV as an indicator of risk, since 80% of the overlapping FOV of the cameras was in the RSA. Equipment failures, majority due to lightning, resulted in only 17 nights with useable data, with unbalanced treatment assignment within turbines and uneven distribution of treatment assignment throughout the observational period. This resulted in a confounding of treatment assignment and seasonal change. Deterrent treatment was measured at Turbine 14 on only 2 of the first 8 usable nights (spanning the period from 8/18-9/17), whereas 6 times on Turbine 8. From 9/18-927, deterrent was on at Turbine 14 on 6 of the remaining 10 nights, and 4 on Turbine 8. We recorded a total of 1,057 bats and observed a reduction in number of events and duration of events at the UAD-activated turbine when it was Turbine 14. When Turbine 8 had the UAD activated, we observed no difference in number or duration of events. Variation between turbines is not unusual and can cause issues when study designs have no true replication, i.e., multiple turbines per treatment. Within-turbine differences suggested a trend for reduced activity when UADs were activated, particularly for Turbine 14. These results may be caused by the overall higher bat activity at Turbine 8 and confounding of treatment assignment and seasonal trends. We mapped 58 bat events in 3-dimensional space (3D), 30 and 28 events during control and treatment conditions, respectively. We observed bats crossing the rotor plane under both control and treatment conditions and observed a total of 40 confirmed or near-collisions (i.e., target close to blade but no visual confirmation of a strike) out of a total of 1,491 medium and high confidence bat observations (880 control, 611 at treatment). Twice as many collisions/possible collisions were observed during control conditions. Given the challenges with the equipment and potential confounding of the data (i.e., different activity levels at the two wind turbines), we were unable to determine whether this initial turbine placement and orientation was optimal. Given no new information on how best to install the devices, we elected to use the same placement and orientation for the comparative study. For the Comparative Study, the objectives were to investigate the relative mortality rates among 4 treatments. We searched the area within 90 m of each turbine daily to recover the highest number of fresh fatalities possible. We were unable to detect a clear reduction in mortality from deterrents alone for any individual species. Surprisingly, mortality rate of the eastern red bat (Lasiurus borealis) was estimated to be 1.3–4.2 times as much when turbines were operating normally and UADs were on than when UADs were off. Reduction in mortality of all bat species combined due to curtailment of turbines was estimated to be between 0%–38%. This effect was nullified when, in addition to curtailment, UADs were on, with 95% confidence interval ranging from a 45% reduction to a 36% increase in mortality. This was likely due to the large proportion of eastern red bats in the total carcass population. Mortality of all low-frequency echolocating bats combined (i.e. hoary bat [L. cinereus], big brown bat [Eptesicus fuscus], silver-haired bat [Lasionycteris noctivagans]) relative to control was lower when curtailed (95% CI: 0%–74%), but the addition of UADs had no detectable effect (95%CI: 13%–79%). The combined treatment reduced mortality in silver-haired bats relative to control by 11%–99%, compared to curtailment (81% reduction–67% increase) or deterrent (82% reduction–67% increase) alone. Because silver-haired bats comprised a large proportion of low-frequency calling bats found during this study, a similar effect was seen for that group. The higher mortality observed for eastern red bats at UAD compared to control could have been caused by several factors, such as the effective range of the UAD, particularly at higher frequencies, behavior, positioning of the devices on the nacelle, or a combination of these. We used 3D thermal videography to compare control and UAD bat behavior from two turbines using a total of 203 3D bat-tracks across 34 nights. We recorded a similar number of bat-tracks between treatment groups, with 51% and 49% for control and UAD, respectively. Due to potential differences in bat behavior around spinning vs stationary turbine blades, we examined UAD effectiveness separately for non-operating turbines (feathered below cut-in speed of 3.5 m/s) and operating (normal operation above wind speed of 3.5 m/s). We found a higher proportion of bat-tracks at operating turbines (82%) compared to non-operating turbines (18%), although this does not account for overall time turbines were operating versus not. At non-operating turbines the UAD appears to be effective at reducing the amount of time, flight length, and number of passes through the rotor plane, compared to control. In addition, we found bats approached turbines similarly between control and treatment turbines, with 61% of control and 63% of UAD bat-tracks originating leeward of the hub. In contrast, at operating turbines, we saw little change in bat behavior in response to UADs. For example, we found an increase in the average duration of bat-tracks between non-operating and operating turbines for UAD but at control turbines we found average duration decreased once turbines became operational. Both control and UAD had a high proportion of bat-tracks that crossed the rotor plane (i.e., collision risk) originate from the windward side when turbines were operational 65% to 92%, respectively. Given that the UAD devices closest to the blades were orientated parallel to the blades, its possible bats were not exposed to the signal until they were close to the turbine blades, as suggested by the slightly higher mean duration within 5 meters of the blades for UAD turbines. Future research should consider concentrating UAD intensity on the areas of risk (i.e. blades) with enough buffer to allow bats to react to the sound before entering the rotor-swept area (RSA). In addition, investigating the potential of installing UAD units windward of the turbine blades (e.g. a hub-mounted UAD), particularly since even under control conditions, a relatively high proportion (65%) of crosses through the blade plane originated windward. 3D thermal videography provided valuable information on future testing strategies (e.g. device placement, UAD orientation) to improve UAD effectiveness when bats are at risk (i.e., operating wind turbines). Across the entire project we experienced issues with the operation and communication with the UADs. Most of the issues occurred during the Feasibility Study and were resolved prior to the Comparability Study. Additional challenges surfaced during the Comparability Study but were remedied immediately and are thought to have little impact on the results. We had logistical constraints at the project that limited our ability use traditional methods in our camera calibration. Several calibrations showed inaccurate scales, which may have been related to inadequate spatial coverage of “points” in the camera calibration volume. Because we were using actual video recordings of bats at a wind turbine, the behavior of bats could have concentrated “points” in specific areas of the turbine (i.e. leeward of nacelle), and limited “points” in other areas, resulting in camera calibration issues. We have plans to address these inconsistencies and improving the software and related methodologies by early 2020.

17 WIND ENERGY↗

Numerical Integration of Viscoelastic Models

Accurate modeling of viscoelasticity remains an important consideration for a variety of materials (e.g. polymers and inorganic glasses). As such, over the previous decades a substantial body of work has been dedicated to developing appropriate constitutive models for viscoelasticity ranging from initial considerations of linear thermoviscoelasticity to more complex non-linear formulations incorporating fictive temperatures or potential energy clocks including the use of both internal state variable(ISV) and hereditary integral representations. Nonetheless, relatively limited (in comparison to plasticity) attention has been paid to the numerical integration of such schemes. In terms of integral based formulations, Taylor et al. first considered the problem of the integration of a linear viscoelasticity model. That work focused on the integration of the hereditary integrals and demonstrated improved performance of the new scheme with a custom finite element code over an existing finite difference reference. Chambers and Becker, using a free volume based shift factor, also considered the integration of the hereditary integrals and the impact on the problem of a pressurized thick-walled cylinder and developed an adaptive scheme to bound the error. Chambers later developed three-point Gauss and composite integration schemes for the hereditary integrals and noted improved accuracy. With respect to ISV-based schemes, formulations for the non-linear Schapery model have been proposed. However, in those efforts greater attention was paid to convergence of the non-linear solution scheme than impact of numerical integration. Various authors (e.g. Holzapfel and Simo and Hughes) have also studied the use of convolution integrals with differential forms of ISVs for temperature-independent formulations. Regardless, while the "potential energy clock" (PEC) and "simplified potential energy clock"(SPEC) models have been used to study a variety of non-linear responses (e.g.), limited attention has been paid to the numerical performance. As will be discussed later, the "clock" at the center of the formulations includes temperature and complex history dependence making the numerical integration of such a model even more challenging. Thus, in the current work an initial effort towards characterizing the numerical integration of the constitutive model through simplified problems is performed. To that end, in Section 2 the theory of the model is briefly presented while the numerical integration is discussed in Section 3. Results of various studies characterizing the numerical behavior and performance are then given in Section 4. Finally, some concluding remarks and thoughts for follow on works are provided in Section 5.

36 MATERIALS SCIENCE↗

Wave-powered water pump for upwelling in aquaculture: Numerical model and ocean test

Wave-powered upwelling can increase the productivity and survivability of several aquaculture species. This enhancement is due to transporting cold, nutrient-rich ocean water, typically found lower in the water column, to the surface. Macroalgaes, like kelp, exhibit increased growth from these altered conditions. The University of New Hampshire’s (UNH) wave-powered water pump (wave pump) is a point absorber wave energy converter (WEC) that uses ocean waves to create relative motion between a spar buoy and a concentric float which drives an internal pump. A numerical model of the wave pump was developed using WEC-Sim to predict device performance in the ocean. Wave pump performance was evaluated during a five day ocean test near Appledore Island in Maine in March 2023, where volumetric flow rate, relative distance between spar and float, and wave conditions were measured. These data were then used for numerical model validation. The ocean deployment recorded the device’s performance in a variety of sea states, with average significant wave heights up to 0.7 m. The ocean test data were compared to the WEC-Sim numerical model of the device with favorable results. Average values of device stroke period, stroke height, and flow rate agreed between the ocean test and model data to within approximately 16 to 22%. Furthermore, the validated numerical model provides a valuable tool for improving the design and developing a commercial-scale, wave-powered water pump for use in aquaculture.

16 TIDAL AND WAVE POWER↗

Roles of Engineering Correlations in Hypersonic Entry Boundary Layer Transition Prediction

Efforts to design and operate hypersonic entry vehicles are constrained by many considerations that involve all aspects of an entry vehicle system. One of the more significant physical phenomenon that affect entry trajectory and thermal protection system design is the occurrence of boundary layer transition from a laminar to turbulent state. During the Space Shuttle Return To Flight activity following the loss of Columbia and her crew of seven, NASA's entry aerothermodynamics community implemented an engineering correlation based framework for the prediction of boundary layer transition on the Orbiter. The methodology for this implementation relies upon the framework of correlation techniques that have been in use for several decades. What makes the Orbiter boundary layer transition correlation implementation unique is that a statistically significant data set was acquired in multiple ground test facilities, flight data exists to assist in establishing a better correlation and the framework was founded upon state of the art chemical nonequilibrium Navier Stokes flow field simulations. The basic tenets that guided the formulation and implementation of the Orbiter Return To Flight boundary layer transition prediction capability will be reviewed as a recommended format for future empirical correlation efforts. The validity of this approach has since been demonstrated by very favorable comparison of recent entry flight testing performed with the Orbiter Discovery, which will be graphically summarized. These flight data can provide a means to validate discrete protuberance engineering correlation approaches as well as high fidelity prediction methods to higher confidence. The results of these Orbiter engineering and flight test activities only serve to reinforce the essential role that engineering correlations currently exercise in the design and operation of entry vehicles. The framework of information-related to the Orbiter empirical boundary layer transition prediction capability will be utilized to establish a fresh perspective on this role, to illustrate how quantitative statistical evaluations of empirical correlations can and should be used to assess accuracy and to discuss what the authors' perceive as a recent heightened interest in the application of high fidelity numerical modeling of boundary layer transition. Concrete results will also be developed related to empirical boundary layer transition onset correlations. This will include assessment of the discrete protuberance boundary layer transition onset data assembled for the Orbiter configuration during post-Columbia Return To Flight. Assessment of these data will conclude that momentum thickness Reynolds number based correlations have superior coefficients and uncertainty in comparison to roughness height based Reynolds numbers, aka Re(sub k) or Re(sub kk). In addition, linear regression results from roughness height Reynolds number based correlations will be evaluated, leading to a hypothesis that non-continuum effects play a role in the processes associated with incipient boundary layer transition on discrete protuberances.

Campbell, Charles H.↗

On the difference between variational and unitary coupled cluster theories

There have been assertions in the literature that the variational and unitary forms of coupled cluster theory lead to the same energy functional. Numerical evidence from previous authors was inconsistent with this claim, yet the small energy differences found between the two methods and the relatively large number of variational parameters precluded an unequivocal conclusion. Using the Lipkin Hamiltonian, we here present conclusive numerical evidence that the two theories yield different energies. The ambiguities arising from the size of the cluster parameter space are absent in the Lipkin model, particularly when truncating to double excitations. We show that in the symmetry adapted basis under strong correlation, the differences between the variational and unitary models are large, whereas they yield quite similar energies in the weakly correlated regime previously explored. We also provide a qualitative argument rationalizing why these two models cannot be the same. Additionally, we study a generalized non-unitary and non-hermitian variant that contains excitation, de-excitation, and mixed operators with different amplitudes and show that it works best when compared to the traditional, variational, unitary, and extended forms of coupled cluster doubles theories

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Validity of the Relation Between Spontaneous and Stimulated Emissions in Semiconductors

The Einstein relation between spontaneous emission and absorption was originally derived for a system consists of a two-state subsystem representing matter and harmonic fields representing radiation. The derivation is based on the detailed balance between these two subsystems under thermal equilibrium. The relationship was later investigated in connection with the interactions between radiation field and solids or semiconductors. The simple derivation dose not hold for semiconductors in general. In certain limiting cases, simple relation was obtained. The validity of this relation is important not only because of its fundamental role connecting two of the most fundamental optical processes in semiconductors, but mostly also because of its wide use as a practical method to measure the optical gain of a semiconductor. The validity of this relation for semiconductors has been an issue of controversial for some time. In this paper we numerically examine the validity of this relationship for several different lineshapes including Lorentzian, Gaussian, Sech, and a convoluted double Lorentzians (CDL). We find out that at relatively low density above transparency level, all first three lineshapes violate the Einstein relation. The relation is approximately valid at high density. At very high density, the validity of the Einstein relation holds well for all three lineshapes. The reason behind this observation is explained. The CDL lineshape has been shown analytically to obey the Einstein relationship previously. We show that for a 2D semiconductor with parabolic bands, the CDL lineshape can be integrated analytically. This analytic lineshape is compared with a simple Lorentzian lineshape.

Ning, Cun-Zheng↗

Unexpected Scaling of Peak Acceleration for a Yielding Mass-Spring System Subjected to a Triangular Base Acceleration Pulse

The use of bounding scenarios is a common practice that greatly simplifies the design and qualification of structures. However, this approach implicitly assumes that the quantities of interest increase monotonically with the input to the structure, which is not necessarily true for nonlinear structures. This paper surveys the literature for observations of nonmonotonic behavior of nonlinear systems and finds such observations in both the earthquake engineering and applied mechanics literature. Numerical simulations of a single degree-of-freedom mass-spring system with an elastic–plastic spring subjected to a triangular base acceleration pulse are then presented, and it is shown that the relative acceleration of this system scales nonmonotonically with the input magnitude in some cases. The equation of motion for this system is solved symbolically and an approximate expression for the relative acceleration is developed, which qualitatively agrees with the nonmonotonic behavior seen in the numerical results. The nonmonotonicity is investigated and found to be a result of dynamics excited by the discontinuous derivative of the base acceleration pulse, the magnitude of which scales nonmonotonically with the input magnitude due to the fact that the first yield of the spring occurs earlier as the input magnitude is increased. The relevance of this finding within the context of defining bounding scenarios is discussed, and it is recommended that modeling be used to perform a survey of the full range of possible inputs prior to defining bounding scenarios.

42 ENGINEERING↗

Soils and sediments host Thermoplasmata archaea encoding novel copper membrane monooxygenases (CuMMOs)

Copper membrane monooxygenases (CuMMOs) play critical roles in the global carbon and nitrogen cycles. Organisms harboring these enzymes perform the first, and rate limiting, step in aerobic oxidation of ammonia, methane, or other simple hydrocarbons. Within archaea, only organisms in the order Nitrososphaerales (Thaumarchaeota) encode CuMMOs, which function exclusively as ammonia monooxygenases. From grassland and hillslope soils and aquifer sediments, we identified 20 genomes from distinct archaeal species encoding divergent CuMMO sequences. These archaea are phylogenetically clustered in a previously unnamed Thermoplasmatota order, herein named the Ca. Angelarchaeales. The CuMMO proteins in Ca. Angelarchaeales are more similar in structure to those in Nitrososphaerales than those of bacteria, and contain all functional residues required for general monooxygenase activity. Ca. Angelarchaeales genomes are significantly enriched in blue copper proteins (BCPs) relative to sibling lineages, including plastocyanin-like electron carriers and divergent nitrite reductase-like (nirk) 2-domain cupredoxin proteins co-located with electron transport machinery. Ca. Angelarchaeales also encode significant capacity for peptide/amino acid uptake and degradation and share numerous electron transport mechanisms with the Nitrososphaerales. Ca. Angelarchaeales are detected at high relative abundance in some of the environments where their genomes originated from. While the exact substrate specificities of the novel CuMMOs identified here have yet to be determined, activity on ammonia is possible given their metabolic and ecological context. The identification of an archaeal CuMMO outside of the Nitrososphaerales significantly expands the known diversity of CuMMO enzymes in archaea and suggests previously unaccounted organisms contribute to critical global nitrogen and/or carbon cycling functions.

59 BASIC BIOLOGICAL SCIENCES↗