Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Generalized 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 163 records · Page 9

Calculation of standard viscoelastic responses with multiple retardation times through analysis of static force spectroscopy AFM data

Here, we explore the physics of an atomic force microscopy (AFM) cantilever tip interacting with a generalized viscoelastic sample containing an arbitrary number of characteristic times, when the cantilever's base is driven with constant velocity toward the sample. This mode of operation, often called static force spectroscopy (SFS), can be harnessed to thoroughly analyze time-dependent viscoelastic information frequently overlooked in experiments. We generalize the solution of previous authors who have studied the standard linear solid model, and offer a solution applicable to any linear viscoelastic model. This generalization is crucial for the prediction of the model's response over wide ranges of time-scale. As a demonstration, successful predictions of harmonic functions (e.g., loss tangent) over a wide frequency range are obtained through analysis of simulated SFS results. In addition, we show that analysis through the generalized solution and previous expressions is no longer valid when the force does not grow linearly in time, so we also deliver an alternate route for extracting the viscoelastic information, which does not rely on the force linearity assumption. Despite the large amount of theoretical content (included for theoretical rigor's sake), the practical user can also benefit from the new procedures offered and the corresponding explanations.

42 ENGINEERING↗

Comparison of some biased estimation methods (including ordinary subset regression) in the linear model

Ridge, Marquardt's generalized inverse, shrunken, and principal components estimators are discussed in terms of the objectives of point estimation of parameters, estimation of the predictive regression function, and hypothesis testing. It is found that as the normal equations approach singularity, more consideration must be given to estimable functions of the parameters as opposed to estimation of the full parameter vector; that biased estimators all introduce constraints on the parameter space; that adoption of mean squared error as a criterion of goodness should be independent of the degree of singularity; and that ordinary least-squares subset regression is the best overall method.

Sidik, S. M.↗

Development and validation of a general purpose linearization program for rigid aircraft models

A FORTRAN program that provides the user with a powerful and flexible tool for the linearization of aircraft models is discussed. The program LINEAR numerically determines a linear systems model using nonlinear equations of motion and a user-supplied, nonlinear aerodynamic model. The system model determined by LINEAR consists of matrices for both the state and observation equations. The program has been designed to allow easy selection and definition of the state, control, and observation variables to be used in a particular model. Also, included in the report is a comparison of linear and nonlinear models for a high performance aircraft.

Duke, E. L.↗

Development and validation of a general purpose linearization program for rigid aircraft models

This paper discusses a FORTRAN program that provides the user with a powerful and flexible tool for the linearization of aircraft models. The program LINEAR numerically determines a linear systems model using nonlinear equations of motion and a user-supplied, nonlinear aerodynamic model. The system model determined by LINEAR consists of matrices for both the state and observation equations. The program has been designed to allow easy selection and definition of the state, control, and observation variables to be used in a particular model. Also, included in the report is a comparison of linear and nonlinear models for a high-performance aircraft.

Duke, E. L.↗

Quantum effects on Higgs-strahlung events at linear colliders within the general two-Higgs-doublet model

The associated production of neutral Higgs bosons with the Z{sup 0} gauge boson (h{sup 0}Z{sup 0}, H{sup 0}Z{sup 0}) is investigated in the context of the future linear colliders, such as the ILC and CLIC, within the general two-Higgs-doublet model (2HDM). We compute the corresponding cross section for the processes e{sup +}e{sup -{yields}}Z{sup 0}h{sup 0}/Z{sup 0}H{sup 0} at one-loop, including the full set of corrections at O({alpha}{sub ew}{sup 3}) together with the leading O({alpha}{sub ew}{sup 4}) terms, in full consistency with the available theoretical and phenomenological constraints. We find that the wave-function corrections to the external Higgs fields are the dominant source of the quantum effects, which turn out to be large and negative (e.g., {delta}{sigma}/{sigma}{approx}-20%/-60%) in all the {radical}(s) range, and located predominantly in the region around tan{beta}{approx}1 and moderate values of the parameter {lambda}{sub 5} (being {lambda}{sub 5}<0). This behavior can be ultimately traced back to the enhancement potential of the triple Higgs-boson self-couplings, a trademark feature of the 2HDM with no counterpart in the Higgs sector of the minimal supersymmetric standard model. Even under this substantial depletion of the one-loop-corrected signal (which is also highly distinctive with respect to the SM expectation for e{sup +}e{sup -{yields}}Z{sup 0}H), the predicted Higgs-strahlung rates comfortably reach a few tens of femtobarn, which means barely {approx}10{sup 3}-10{sup 4} events per 500 fb{sup -1} of integrated luminosity. Because of their great complementarity, we argue that the combined analysis of the Higgs-strahlung events (h{sup 0}Z{sup 0}, H{sup 0}Z{sup 0}) and the previously computed one-loop Higgs-pair production processes (H{sup 0}A{sup 0}) could be instrumental to probe the structure of the Higgs sector at future linac facilities.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Parallel, iterative solution of sparse linear systems: Models and architectures

A model of a general class of asynchronous, iterative solution methods for linear systems is developed. In the model, the system is solved by creating several cooperating tasks that each compute a portion of the solution vector. A data transfer model predicting both the probability that data must be transferred between two tasks and the amount of data to be transferred is presented. This model is used to derive an execution time model for predicting parallel execution time and an optimal number of tasks given the dimension and sparsity of the coefficient matrix and the costs of computation, synchronization, and communication. The suitability of different parallel architectures for solving randomly sparse linear systems is discussed. Based on the complexity of task scheduling, one parallel architecture, based on a broadcast bus, is presented and analyzed.

Reed, D. A.↗

Modelling intergranular and transgranular micro-cracking in polycrystalline materials

Presented is a grain boundary formulation for intergranular and transgranular micro-cracking in three-dimensional polycrystalline aggregates.. The formulation is based on the displacement and stress boundary integral equations of solid mechanics and it has the advantage of expressing the polycrystalline problem in terms of grain boundary variables only. The individual grains within the polycrystalline morphology are modelled as generally anisotropic linear elastic domains with random spatial orientation. Transgranular micro-cracking is assumed to occur along specific cleavage planes, whose orientation in space within the grains depend upon the crystallographic lattice. Both intergranular and transgranular micro-cracking are modelled using suitably defined cohesive laws, whose parameters characterise the behaviour of the two mechanisms. The algorithm developed to track the inter/transgranular micro-cracking history is presented and discussed. Several numerical tests involving pseudo-3D and fully 3D morphologies are performed and analysed. The presented numerical results show that the developed formulation is capable of tracking the initiation and evolution of both intergranular and transgranular cracking as well as their competition, thus providing a useful tool for the study of damage micro-mechanics.

36 MATERIALS SCIENCE↗

An Analysis of Terrestrial and Aquatic Environmental Controls of Riverine Dissolved Organic Carbon in the Conterminous United States

Riverine carbon cycling is an important, but insufficiently investigated component of the global carbon cycle. Analyses of environmental controls on riverine carbon cycling are critical for improved understanding of mechanisms regulating carbon processing and storage along the terrestrial-aquatic continuum. Here, we compile and analyze riverine dissolved organic carbon (DOC) concentration data from 1402 United States Geological Survey (USGS) gauge stations to examine the spatial variability and environmental controls of DOC concentrations in the United States (U.S.) surface waters. DOC concentrations exhibit high spatial variability, with an average of 6.42 ± 6.47 mg C/ L (Mean ± Standard Deviation). In general, high DOC concentrations occur in the Upper Mississippi River basin and the Southeastern U.S., while low concentrations are mainly distributed in the Western U.S. Single-factor analysis indicates that slope of drainage areas, wetlands, forests, percentage of first-order streams, and instream nutrients (such as nitrogen and phosphorus) pronouncedly influence DOC concentrations, but the explanatory power of each bivariate model is lower than 35%. Analyses based on the general multi-linear regression models suggest DOC concentrations are jointly impacted by multiple factors. Soil properties mainly show positive correlations with DOC concentrations; forest and shrub lands have positive correlations with DOC concentrations, but urban area and croplands demonstrate negative impacts; total instream phosphorus and dam density correlate positively with DOC concentrations. Notably, the relative importance of these environmental controls varies substantially across major U.S. water resource regions. In addition, DOC concentrations and environmental controls also show significant variability from small streams to large rivers, which may be caused by changing carbon sources and removal rates by river orders. In sum, our results reveal that general multi-linear regression analysis of twenty one terrestrial and aquatic environmental factors can partially explain (56%) the DOC concentration variation. In conclusion, this study highlights the complexity of the interactions among these environmental factors in determining DOC concentrations, thus calls for processes-based, non-linear methodologies to constrain uncertainties in riverine DOC cycling.

54 ENVIRONMENTAL SCIENCES↗

Do Remote Camera Arrangements and Image Capture Settings Improve Individual Identification of Golden Eagles?

Individual identification of animals from camera traps has become an important task in wildlife research, but camera deployment methods often do not facilitate this important undertaking. Identification of individual golden eagles (Aquila chrysaetos) is possible using uniquely marked rectrices, but no studies have explored methods to maximize the rate of individual identification from camera images. Furthermore, our objectives were to assess whether different camera heights (1 m vs. 3 m), image capture settings (one image after a 1-min delay vs. burst of 5 images after a 30 sec delay), and arrangements relative to bait (dorsally vs. ventrally aimed) affected views of rectrices on golden eagles and our ability to identify individuals. We conducted our study from 15 December 2016 to 3 March 2017 on the Savannah River Site, South Carolina. First, we developed a scoring system based on views of rectrices and used a linear mixed-effects model to compare image scores among different camera arrangements and image settings. Next, after identifying individual eagles, we used generalized linear mixed-effects models to compare total individual eagle detections, total days an individual was detected, and probability of obtaining an unknown individual identification among camera arrangements and settings. Overall, we scored a total of 27,499 images, with 8,083 providing views of marked rectrices that allowed identification of 18 individual eagles. Average image scores and proportion of images suitable for individual identification were higher from elevated (3 m) camera arrangements than standard arrangements (1 m) across sites. Regardless of camera height, faster frequency of image capture provided more images that could be used to identify individuals and the most trap days per individual. Researchers and managers should consider deploying elevated cameras traps with faster frequency of image capture to improve data quality and potential for analysis of golden eagle populations and trends across the species’ range.

60 APPLIED LIFE SCIENCES↗

A simple nonlinear model for the return to isotropy in turbulence

A quadratic nonlinear generalization of the linear Rotta model for the slow pressure-strain correlation of turbulence is developed. The model is shown to satisfy realizability and to give rise to no stable non-trivial equilibrium solutions for the anisotropy tensor in the case of vanishing mean velocity gradients. The absence of stable non-trivial equilibrium solutions is a necessary condition to ensure that the model predicts a return to isotropy for all relaxational turbulent flows. Both the phase space dynamics and the temporal behavior of the model are examined and compared against experimental data for the return to isotropy problem. It is demonstrated that the quadratic model successfully captures the experimental trends which clearly exhibit nonlinear behavior. Direct comparisons are also made with the predictions of the Rotta model and the Lumley model.

Sarkar, Sutanu↗

A simple nonlinear model for the return to isotropy in turbulence

A quadratic nonlinear generalization of the linear Rotta model for the slow pressure-strain correlation of turbulence is developed. The model is shown to satisfy realizability and to give rise to no stable nontrivial equilibrium solutions for the anisotropy tensor in the case of vanishing mean velocity gradients. The absence of stable nontrivial equilibrium solutions is a necessary condition to ensure that the model predicts a return to isotropy for all relaxational turbulent flows. Both the phase space dynamics and the temporal behavior of the model are examined and compared against experimental data for the return to isotropy problem. It is demonstrated that the quadratic model successfully captures the experimental trends which clearly exhibit nonlinear behavior. Direct comparisons are also made with the predictions of the Rotta model and the Lumley model.

Sarkar, Sutanu↗

Environmental controls on the light use efficiency of terrestrial gross primary production

Abstract Gross primary production (GPP) by terrestrial ecosystems is a key quantity in the global carbon cycle. The instantaneous controls of leaf‐level photosynthesis are well established, but there is still no consensus on the mechanisms by which canopy‐level GPP depends on spatial and temporal variation in the environment. The standard model of photosynthesis provides a robust mechanistic representation for C 3 species; however, additional assumptions are required to “scale up” from leaf to canopy. As a consequence, competing models make inconsistent predictions about how GPP will respond to continuing environmental change. This problem is addressed here by means of an empirical analysis of the light use efficiency (LUE) of GPP inferred from eddy covariance carbon dioxide flux measurements, in situ measurements of photosynthetically active radiation (PAR), and remotely sensed estimates of the fraction of PAR (fAPAR) absorbed by the vegetation canopy. Focusing on LUE allows potential drivers of GPP to be separated from its overriding dependence on light. GPP data from over 100 sites, collated over 20 years and located in a range of biomes and climate zones, were extracted from the FLUXNET2015 database and combined with remotely sensed fAPAR data to estimate daily LUE. Daytime air temperature, vapor pressure deficit, diffuse fraction of solar radiation, and soil moisture were shown to be salient predictors of LUE in a generalized linear mixed‐effects model. The same model design was fitted to site‐based LUE estimates generated by 16 terrestrial ecosystem models. The published models showed wide variation in the shape, the strength, and even the sign of the environmental effects on modeled LUE. These findings highlight important model deficiencies and suggest a need to progress beyond simple “goodness of fit” comparisons of inferred and predicted carbon fluxes toward an approach focused on the functional responses of the underlying dependencies.

54 ENVIRONMENTAL SCIENCES↗

Conservation planning for freshwater–marine carryover effects on Chinook salmon survival

Experiences of migratory species in one habitat may affect their survival in the next habitat, in what is known as carryover effects. These effects are especially relevant for understanding how freshwater experience affects survival in anadromous fishes. Here, we study the carryover effects of juvenile salmon passage through a hydropower system (Snake and Columbia rivers, northwestern United States). To reduce the direct effect of hydrosystem passage on juveniles, some fishes are transported through the hydrosystem in barges, while the others are allowed to migrate in-river. Although hydrosystem survival of transported fishes is greater than that of their run-of-river counterparts, their relative juvenile-to-adult survival (hereafter survival) can be less. We tested for carryover effects using generalized linear mixed effects models of survival with over 1 million tagged Chinook salmon, Oncorhynchus tshawytscha (Walbaum) (Salmonidae), migrating in 1999–2013. Carryover effects were identified with rear-type (wild vs. hatchery), passage-type (run-of-river vs. transported), and freshwater and marine covariates. Importantly, the Pacific Decadal Oscillation (PDO) index characterizing cool/warm (i.e., productive/nonproductive) ocean phases had a strong influence on the relative survival of rear- and passage-types. Specifically, transportation benefited wild Chinook salmon more in cool PDO years, while hatchery counterparts benefited more in warm PDO years. Transportation was detrimental for wild Chinook salmon migrating early in the season, but beneficial for later season migrants. Hatchery counterparts benefited from transportation throughout the season. Altogether, wild fish could benefit from transportation approximately 2 weeks earlier during cool PDO years, with still a benefit to hatchery counterparts. Furthermore, we found some support for hypotheses related to higher survival with increased river flow, high predation in the estuary and plume areas, and faster migration and development-related increased survival with temperature. Thus, pre- and within-season information on local- and broad-scale conditions across habitats can be useful for planning and implementing real-time conservation programs.

54 ENVIRONMENTAL SCIENCES↗

Estimating mortality rates among passerines caught for ringing with mist nets using data from previously ringed birds

Abstract Mist netting is the most commonly used method for catching birds for scientific ringing, but despite decades of use, there have been few attempts to quantify the associated potential risks to the individuals caught. Any incidence of mortality through capture and handling, however low, is of potential ethical concern and may also introduce biases into the data. We estimate the mortality rate associated with capture of previously ringed (recaptured) passerines from the British and Irish Ringing Scheme (c. 1.5 million records) caught using mist nets. The importance of species, age, mass, month, time, previous captures, and an index of predator occurrence was tested using generalized linear mixed‐effects models. The average mortality rate was 0.0011, most of which was reported to occur before the individuals had been extracted from the nets (c. 70% of incidents). Juveniles appeared to be at higher risk and the incidence of predation from mist nets was seasonal, with increased risk during the winter. Species differed in their reported mortality rates with the apparent risk being greatest for Chiffchaff Phylloscopus collybita (0.0029) and Bullfinch Pyrrhula pyrrhula (0.0027). To improve our understanding (and hence minimize risk in future), we recommend collecting more complete data on incidences of mortality, and also injuries; exercising increased care when the species we have identified as being at greater risk are likely to be captured, and ensuring there are robust procedures for the checking of nets (as most reported incidents of mortality occur before handling). We also recommend that all Ringing Schemes should collate and make available data on capture‐related mortality. Overall rates of mortality associated with capture, although, were low and support the use of mist netting as a safe capture technique, without undue bias from mortality, when used by appropriately trained individuals.

Clewley, Gary D.↗

A new approach to estimate fecundity rate from inter‐birth intervals

Abstract Accurate estimates of fecundity rate are key to population assessments and effectively direct conservation efforts. We present a new approach to estimate fecundity rate based on the probability of a female giving birth, conditional on a previous birth t years ago, from which an expected inter‐birth interval ( IBI ) can be estimated. We use generalized linear mixed‐effects models to account for individual and temporal variability and apply the approach to individual reproductive histories of bottlenose dolphins ( Tursiops truncatus ) from the east coast of Scotland. We estimate a fecundity rate of 0.222 (95% CI = 0.218–0.253) and an expected IBI of 4.49 yr (95% CI = 3.94–4.93 yr). We use simulated data samples to show that the approach produces estimates with a minimum bias of <3%. Simulations are also used to investigate the effect of the most common data‐driven biases in the estimates of birth intervals and fecundity rate; we recommend longitudinal studies of at least 10 yr and capture probabilities of at least 0.3 when using this methodology. The approach may be modified to incorporate other parameters of interest and should be applicable to any population with comprehensive data on birth intervals.

Arso Civil, Mònica↗

Estimation of a Self-Consistent Set of Radiobiological Parameters From Hypofractionated Versus Standard Radiation Therapy of Prostate Cancer

Purpose: To determine a self-consistent set of radiobiological parameters in prostate cancer. Methods and Materials: A method to estimate intrinsic radiosensitivity (α), fractionation sensitivity (α/β), repopulation doubling time, number of clonogens, and kick-off time for accelerated repopulation of prostate cancer has been developed. Based on the generalized linear-quadratic model and without assuming the isoeffective hypothesis, the potential applications of the method were investigated using the clinical outcome of biochemical relapse-free survival recently reviewed in the literature. The strengths and limitations of the method, regarding the fitted parameters and 95% confidence intervals (CIs), are also discussed. Results: Our best estimate of α/β is 2.96 Gy (95% CI 2.41-3.53 Gy). The corresponding α value is 0.16 Gy{sup −1} (95% CI 0.14-0.18 Gy{sup −1}), which is compatible with a realistic number of clonogens: 6.5 × 10{sup 6} (95% CI 1.5 × 10{sup 6}-2.1 × 10{sup 7}). The estimated cell doubling time is 5.1 days (95% CI 4.2-7.2 days), very low if compared with that reported in the literature. This corresponds to the dose required to offset the repopulation occurring in 1 day of 0.52 Gy/d (95% CI 0.32-0.68 Gy/d). However, a long kick-off time of 31 days (95% CI 22-41 days) from the start of radiation therapy was found. Conclusion: The proposed analytic/graphic method has allowed the fitting of clinical data, providing a self-consistent set of radiobiological parameters for prostate cancer. With our analysis we confirm a low value for α/β with a correspondingly high value of intrinsic radiosensitivity, a realistic average number of clonogens, a long kick-off time for accelerated repopulation, and a surprisingly fast repopulation that suggests the involvement of subpopulations of specifically tumorigenic stem cells during continuing radiation therapy.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Nonlinear mode coupling and energetics of driven magnetized shear-flow turbulence

To comprehensively understand the saturation of two-dimensional (2D) magnetized Kelvin–Helmholtz-instability-driven turbulence, energy transfer analysis is extended from the traditional interaction between scales to include eigenmode interactions, by using the nonlinear couplings of linear eigenmodes of the ideal instability. While both kinetic and magnetic energies cascade to small scales, a significant fraction of turbulent energy deposited by unstable modes in the fluctuation spectrum is shown to be re-routed to the conjugate-stable modes at the instability scale. They remove energy from the forward cascade at its inception. The remaining cascading energy flux is shown to attenuate exponentially at a small scale, dictated by the large-scale stable modes. Guided by a widely used instability-saturation assumption, a general quasi-linear model of instability is tested by retaining all nonlinear interactions except those that couple to the large-scale stable modes. These complex interactions are analytically removed from the magnetohydrodynamic equations using a novel technique. Observations are an explosive large-scale vortex separation instead of the well-known merger of 2D, a dramatic enhancement in turbulence level and spectral energy fluxes, and a reduced small-scale dissipation length scale. These show the critical role of the stable modes in instability saturation. Possible reduced-order turbulence models are proposed for fusion and astrophysical plasmas, based on eigenmode-expanded energy transfer analyses.

Physics↗

Nonadiabatic transitions in Landau-Zener grids: Integrability and semiclassical theory

We demonstrate that the general model of a linearly time-dependent crossing of two energy bands is integrable. Namely, the Hamiltonian of this model has a quadratically time-dependent commuting operator. We apply this property to four-state Landau-Zener (LZ) models that have previously been used to describe the Landau-Stückelberg interferometry experiments with an electron shuttling between two semiconductor quantum dots. The integrability then leads to simple but nontrivial exact relations for the transition probabilities. In addition, the integrability leads to a semiclassical theory that provides analytical approximation for the transition probabilities in these models for all parameter values. Lastly, the results predict a dynamic phase transition, and show that similarly looking models belong to different topological classes.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗