Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Model fitting”

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 289 records · Page 16

VLBI OBSERVATION OF MICROQUASAR CYG X-3 DURING AN X-RAY STATE TRANSITION FROM SOFT TO HARD IN THE 2007 MAY-JUNE FLARE

We present a radio observation of microquasar Cyg X-3 during an X-ray state transition from ultrasoft to hard state in the 2007 May-June flare using the VLBI Exploration of Radio Astrometry at 22 GHz. During the transition, a short-lived mini-flare of {approx}< 3 hr was detected prior to the major flare. In such a transition, a jet ejection is believed to occur, but there have been no direct observations to support it. An analysis of Gaussian fits to the observed visibility amplitudes shows a time variation of the source axis, or a structural change, during the mini-flare. Our model fits, together with other multiwavelength observations in the radio, soft, and hard X-rays, and the shock-in-jet models for other flaring activities at GHz wavebands, suggest a high possibility of synchrotron flares during the mini-flare, indicative of a predominant contribution from jet activity. Therefore, the mini-flare with an associated structural change is indicative of a jet ejection event in the state transition from ultrasoft to hard state.

79 ASTRONOMY AND ASTROPHYSICS↗

A Novel Eco-Biosorbent for Decontamination of Hazardous Dye from Aqueous Medium

A novel multi-component biosorbent consisted of pine, oak, hornbeam and fir sawdust biomasses was used to eliminate a common hazardous dye (Malachite green, MG) from aqueous solution. The influences of different parameters like solution pH, biosorbent amount, dye concentration and contact time on the removal process were thoroughly studied to evaluate optimum biosorption conditions. The pseudo-first-order, pseudo-second-order, Elovich, logistic and intra-particle diffusion models were used to describe the biosorption kinetics. Logistic model represented the kinetic experimental data well. The experimental equilibrium biosorption data were analyzed by Freundlich, Langmuir and Dubinin–Radushkevich isotherms. Freundlich model fitted better with the experimental data compared to others. In addition, the mean free energy obtained from Dubinin–Radushkevich model and the standard Gibbs free energy change indicated that the nature of biosorption was feasible, spontaneous and physical type. Overall, the novel eco-biosorbent could be used as a promising biosorbent to remove such unsafe dye molecules from aqueous environments.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Age-specific case data reveal varying dengue transmission intensity in US states and territories

Dengue viruses (DENV) are endemic in the US territories of Puerto Rico, American Samoa, and the US Virgin Islands, with focal outbreaks also reported in the states of Florida and Hawaii. However, little is known about the intensity of dengue virus transmission over time and how dengue viruses have shaped the level of immunity in these populations, despite the importance of understanding how and why levels of immunity against dengue may change over time. These changes need to be considered when responding to future outbreaks and enacting dengue management strategies, such as guiding vaccine deployment. We used catalytic models fitted to case surveillance data stratified by age from the ArboNET national arboviral surveillance system to reconstruct the history of recent dengue virus transmission in Puerto Rico, American Samoa, US Virgin Islands, Florida, Hawaii, and Guam. We estimated average annual transmission intensity (i.e., force of infection) of DENV between 2010 and 2019 and the level of seroprevalence by age group in each population. We compared models and found that assuming all reported cases are secondary infections generally fit the surveillance data better than assuming all cases are primary infections. Using the secondary case model, we found that force of infection was highly heterogeneous between jurisdictions and over time within jurisdictions, ranging from 0.00008 (95% CrI: 0.00002–0.0004) in Florida to 0.08 (95% CrI: 0.044–0.14) in American Samoa during the 2010–2019 period. For early 2020, we estimated that seropositivity in 10 year-olds ranged from 0.09% (0.02%–0.54%) in Florida to 56.3% (43.7%–69.3%) in American Samoa. In the absence of serological data, age-specific case notification data collected through routine surveillance combined with mathematical modeling are powerful tools to monitor arbovirus circulation, estimate the level of population immunity, and design dengue management strategies.

60 APPLIED LIFE SCIENCES↗

Characterizing the Roles of Biogeochemical Cycling and Ocean Circulation in Regulating Marine Copper Distributions

Copper (Cu) is a key micronutrient for marine phytoplankton. Its oceanic biogeochemical cycle has elicited considerable attention due to dissolved Cu exhibiting a unique linear profile with depth. Several processes have been proposed for explaining this behavior. Here, in this study, we characterize the relationships between the observed Cu, PO 4 3– , and Si on a global scale. We find that the depth profiles of Cu resemble those of Si more than of PO 4 3– in the global ocean. To understand their relationships, we couple the biogeochemical and internal circulation processes in a model, fitting optimal Cu:PO 4 3– uptake ratios and remineralization length-scales to replicate the marine Cu distributions. The modeling results suggest that Cu uptake needs in the Southern Ocean are substantially higher than those in other oceanic regions. In addition, our modeling results indicate a deep Cu remineralization in the global ocean. We offer an alternative mechanism that relies on biogeochemical cycling and internal circulation to produce the linear depth profiles of dissolved Cu. Our results suggest that diatoms are likely the major phytoplankton dominating oceanic Cu cycling.

54 ENVIRONMENTAL SCIENCES↗

Physics-informed Estimation of the Covariance Matrix for Various Neutron Spectra

A method for estimating covariance matrices which capture the uncertainties in calculated reactor spectra has been developed. This method is based on perturbing the parameters of a physics-based analytic model fitted to a calculated spectrum. The covariance of the perturbed analytic spectra imposes energy-dependent correlations due to the physics of the neutron processes in the reactor, i.e., a fission component, a 1/E down-scatting component, and a thermal Maxwellian component. An analytic model is developed which is shown to produce good fits to several reactor environments. The covariance matrices produced via this method are then used as the prior spectrum in STAYSL least squares spectrum adjustment where it is combined with integral metrics, such as activation measurements, to produce a high-fidelity neutron spectrum characterization. It was concluded that the methodology showed agreeable results for the ACRR free-field spectrum adjustment in STAYSL resulting in a 𝜒 2 value of 2.21 (per degree of freedom), but further work is needed to describe scattering and interface regions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Primary and Low-Strain Creep Models for 9Cr Tempered Martensitic Steels Including the Effects of Irradiation Softening and High-Helium Re-Hardening

Primary and low-strain creep represents a very important integrity challenge to large, complex structures, like fusion reactors. Here, we develop a predictive empirical primary creep model for 9Cr tempered martensitic steels (TMS), relating the applied stress (σ) to strain (ε), time (t) and temperature (T). The most accurate model is based on the applied σ normalized by the steel’s T-dependent ultimate tensile stress (σo), σ/σo(T). The model, fit to 17 heats of 9Cr TMS, yielded a σ root mean square error (RMSE) of ≈±11 MPa. Notably, the model also provides robust predictions for all the other TMS, when calibrated only by the fusion candidate Eurofer97 database. The model was extended to explore two possible effects of neutron irradiation, which produces both displacements per atom (dpa) and helium (He in atomic parts per million, appm) damage. These effects, which have not been previously considered, include: (a) softening, as a function of dpa, at T > ≈400–450 °C, in low-He fission environments (<1 He/dpa); and (b) subsequent re-hardening in high-He (≥10 He/dpa) fusion first-wall environments. The irradiation effect models predict (a) accelerated primary creep due to irradiation softening; and (b) fully arrested creep due to high-He re-hardening.

Alam, Md Ershadul (ORCID:0000000345968026)↗

Bayesian calibration of terrestrial ecosystem models: A study of advanced Markov chain Monte Carlo methods

Calibration of terrestrial ecosystem models is important but challenging. Bayesian inference implemented by Markov chain Monte Carlo (MCMC) sampling provides a comprehensive framework to estimate model parameters and associated uncertainties using their posterior distributions. The effectiveness and efficiency of the method strongly depend on the MCMC algorithm used. In this study, a Differential Evolution Adaptive Metropolis (DREAM) algorithm was used to estimate posterior distributions of 21 parameters for the data assimilation linked ecosystem carbon (DALEC) model using 14 years of daily net ecosystem exchange data collected at the Harvard Forest Environmental Measurement Site eddy-flux tower. The DREAM is a multi-chain method and uses differential evolution technique for chain movement, allowing it to be efficiently applied to high-dimensional problems, and can reliably estimate heavy-tailed and multimodal distributions that are difficult for single-chain schemes using a Gaussian proposal distribution. The results were evaluated against the popular Adaptive Metropolis (AM) scheme. DREAM indicated that two parameters controlling autumn phenology have multiple modes in their posterior distributions while AM only identified one mode. The calibration of DREAM resulted in a better model fit and predictive performance compared to the AM. DREAM provides means for a good exploration of the posterior distributions of model parameters. Lastly, it reduces the risk of false convergence to a local optimum and potentially improves the predictive performance of the calibrated model.

54 ENVIRONMENTAL SCIENCES↗

[Re] Drivers of evapotranspiration from boreal wildfires

Computational reproducibility is a difficult challenge across science. I attempted to use R 3.6.1 to reproduce linear model fits, done originally using v2.6.0 for a 2009 paper on the drivers of large-scale forest evapotranspiration after wildfire. Model outputs were largely identical, aside from minor formatting changes, except for one–out of 12 total– regression in which the median residual value changed very slightly (in the sixth decimal place). I suggest that this essentially successful reproducibility is due to the relative simplicity of the script, its use of only base R functions, and R’s historically conservative approach to breaking changes.

54 ENVIRONMENTAL SCIENCES↗

THE ROLE OF THE ACCRETION DISK, DUST, AND JETS IN THE IR EMISSION OF LOW-LUMINOSITY ACTIVE GALACTIC NUCLEI

We use recent high-resolution infrared (IR; 1-20 μm) photometry to examine the origin of the IR emission in low-luminosity active galactic nuclei (LLAGN). The data are compared with published model fits that describe the spectral energy distribution (SED) of LLAGN in terms of an advection-dominated accretion flow, truncated thin accretion disk, and jet. The truncated disk in these models is usually not luminous enough to explain the observed IR emission, and in all cases its spectral shape is much narrower than the broad IR peaks in the data. Synchrotron radiation from the jet appears to be important in very radio-loud nuclei, but the detection of strong silicate emission features in many objects indicates that dust must also contribute. We investigate this point by fitting the IR SED of NGC 3998 using dusty torus and optically thin (τ{sub mid-IR} ∼ 1) dust shell models. While more detailed modeling is necessary, these initial results suggest that dust may account for the nuclear mid-IR emission of many LLAGN.

79 ASTRONOMY AND ASTROPHYSICS↗

A simplified regression building thermal modelling method for detached two-floor house in U.S.

The residential building sector accounts for approximately 37% of total U.S. electricity consumption. Within the residential building sector, heating and cooling is the main target for peak load shifting/reduction since it is the largest contributor to peak demand. In fact, the flexibility of residential HVAC loads can provide continuous variation of demand to provide grid services by varying their demand over a baseline. The performance of HVAC load control to provide grid services relies heavily on the accuracy of indoor air temperature or cooling/heating demand predictions and therefore the quality of building model. Besides forward models, popular building models are data-driven models which can be broken down into two categories: gray-box model, e.g. Resistance-Capacitance (RC) model and black-box model. RC model, also called lumped capacitance or network model, which is constituted with electrical analogue pattern with resistance (R) and capacitance (C). In general, RC models require considerable computation burden and long periods of data to train limited number of model coefficients. Black-box models have gained increasing interest due to their capability in analyzing large-scale data and flexibility in practical applications. But, the data-mining based (machine learning algorithms/techniques based) models tend to have invisible model structures which poses a problem when trying to use the model for optimal control or model predictive control of the HVAC system. Hence, there is a continuing need for efficient online system identification techniques, which can provide explicit parameters for the model. Traditional regression models fit well for this specific purpose. This paper presents an innovative way to predict average indoor temperature in separate floors of typical detached residential house. A rolling horizon linear regression model, which includes online adaptive correction component, is proposed to predict the temperature difference between downstairs and upstairs. A RC model is used to predict the overall mean indoor air temperature. Since the adaptive algorithm needs to be implemented online, a less computation-demanding polynomial fitting algorithm is adopted. This kind of fitting problem can be cast as linear regression problem with multiple variables, parameters of which can be efficiently obtained by well-known gradient descent method.The validation is conducted by comparing the predicted results with the results from data-mining based models as well as measured data from a real typical detached two-floor house. The results show that the developed method has satisfactory performance in predicting the building indoor temperature in 1st and 2nd floors.

Cui, Borui↗

Constrained partitioning of autotrophic and heterotrophic respiration reduces model uncertainties of forest ecosystem carbon fluxes but not stocks

Abstract We partitioned the soil carbon dioxide flux ( R s ) into its respective autotrophic and heterotrophic components in a mature temperate‐boreal forest (Howland Forest in Maine, USA). We combined automated chamber measurements of R s with two different partitioning methods: (1) a classic root trenching experiment and (2) a radiocarbon ( 14 C) mass balance approach. With a model‐data fusion approach, we used these data to constrain a parsimonious ecosystem model (FöBAAR), and we investigated differences in modeled C fluxes and pools under both current and future climate scenarios. The trenching experiment indicated that heterotrophic respiration accounted for 53 ± 11% of total R s . In comparison, using the 14 C method, the heterotrophic contribution was 42 ± 9%. For both current and future model runs, incorporating the partitioning data as constraints substantially reduced the uncertainties of autotrophic and heterotrophic respiration fluxes. Moreover, with best fit model parameters, the two partitioning methods yielded fundamentally different estimates of the relative contributions of autotrophic and heterotrophic respiration to total R s , especially at the annual time scale. Surprisingly, however, modeled soil C and biomass C pool size trajectories did not differ significantly between model runs based on the different methods. Instead, model differences in partitioning were compensated for by changes in C allocation, resulting in similar, but still highly uncertain, soil C pool trajectories. Our findings show that incorporating constraints on the partitioning of R s can reduce model uncertainties of fluxes but not pools, and the results are sensitive to the partitioning method used.

Carbone, Mariah S.↗

An Eco-friendly Porous Poly(imide-ether)s for the Efficient Removal of Methylene Blue: Adsorption Kinetics, Isotherm, Thermodynamics and Reuse Performances

Poly(imide-ether)s (PIEs), which is porous in nature, was synthesized using aromatic diamines and dianhydrides via the solution polycondenzation reaction. In the process, PIEs was used as an adsorbent for the removal of basic dye [Methylene blue (MB)] from the aqueous solution. The structure and the thermal stability of PIEs were characterized by FT-IR, {sup 1}H-NMR, XRD and TG analysis, and the adsorption behaviour of the PIEs was confirmed using UV/Visible, FT-IR and SEM with EDX analysis. Batch adsorption experiments were carried out and their parameters (pH 6.95), initial MB concentration (200 mg/L), adsorbent dosage (50 mg), contact time (120 min) and temperature (28 °C) were optimized. The batch adsorption experimental data of the effect of contact time, initial MB concentration and temperature were evaluated using the adsorption kinetics (Pseudo-first order, Pseudo-second order, Elovich and Intra-particle diffusion kinetic), isotherm (Langmuir, Freundlich, Redlich–Peterson and Sips isotherms) and thermodynamic models respectively. In the study, adsorption of MB onto PIEs followed Pseudo-second kinetics and found to be the best fitting model based on the obtained experimental data. The homogeneity and heterogeneity surface nature of PIEs were elucidated with two and three isotherm parameters respectively. The adsorption thermodynamic study indicates that the adsorption of MB onto PIEs was endothermic and spontaneous in nature. The adsorbed PIEs solid waste was utilized to prepare the polymer composite, and its property was characterized. Graphical Abstract: .

54 ENVIRONMENTAL SCIENCES↗

Gas-phase dehydration of tetrahydrofurfuryl alcohol to dihydropyran over γ-Al 2 O 3

Gas-phase dehydration of tetrahydrofurfuryl alcohol (THFA) to 3,4-2H-dihydropyran (DHP) was investigated over solid acid catalysts. A γ-Al 2 O 3 catalyst resulted in 90% DHP yield from THFA. The γ-Al 2 O 3 deactivated due to solid coke formation but was nearly fully re-generable upon a high temperature calcination step (only 2.5% activity loss after 3 rd regeneration). The high catalytic activity and selectivity of γ-Al 2 O 3 was correlated to its high Lewis to Brønsted acid site ratio, as confirmed by ammonia temperature-programmed desorption (NH 3 -TPD) and isopropylamine temperature-programmed desorption (IPA-TPD). According to isotopic-labelling studies, a reaction mechanism was proposed in which THFA initially dehydrates into a carbenium intermediate prior to Wagner-Meerwien rearrangement into DHP. A kinetic model of THFA dehydration over γ-Al 2 O 3 was developed according to kinetic experimental data. The best-fit model implied the rate-determining step was the surface reaction of adsorbed THFA dissociating into adsorbed DHP and water.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Improvement of a mixture experiment model relating the component proportions to the size of nanonized itraconazole particles in extemporary suspensions

A paper by Foglio Bonda et al. published previously in this journal (2016, Vol. 83, pp. 175–183) discussed the use of mixture experiment design and modeling methods to study how the proportions of three components in an extemporaneous oral suspension affected the mean diameter of drug particles (Z ave ). The three components were itraconazole (ITZ), Tween 20 (TW20), and Methocel® E5 (E5). This commentary addresses some errors and other issues in the previous paper, and also discusses an improved model relating proportions of ITZ, TW20, and E5 to Z ave . The improved model contains six of the 10 terms in the full-cubic mixture model, which were selected using a different cross-validation procedure than used in the previous paper. In conclusion, compared to the four-term model presented in the previous paper, the improved model fit the data better, had excellent cross-validation performance, and the predicted Z ave of a validation point was within model uncertainty of the measured value.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Simulations of grain growth in tungsten armor materials under ARC plasma edge operation conditions using an integrated plasma-edge/materials model

An integrated model of grain growth deuterium-exposed tungsten polycrystals, consisting of a two-dimensional vertex dynamics model fitted to atomistic data, has been developed to assess the grain growth kinetics of deuterium-exposed polycrystalline tungsten (W). The model tracks the motion of grain boundaries under the effect of driving forces stemming from grain boundary curvature and differential deuterium concentration accumulation. Here, we apply the model to experimentally synthesized tungsten polycrystals under deuterium-saturated conditions relevant to the ARC concept design. The results indicate rapid grain growth kinetics in the near-surface region adjacent to the plasma, where the temperature reaches 1400 K, whereas the microstructure remains stable deeper in the material with the lower temperature of 1000 K. The combined modeling and analysis further reveal that monolithic tungsten produced via conventional fabrication routes is highly susceptible to grain coarsening at temperatures exceeding 1000 K, largely independent of the magnitude of the applied driving force. Moreover, the accumulation of deuterium near grain boundaries has a pronounced inhibitory effect on grain boundary migration. High-angle grain boundaries ( > 50°) contribute more significantly to the overall grain growth process.

36 MATERIALS SCIENCE↗

WE-H-BRA-03: Development of a Model to Include the Evolution of Resistant Tumor Subpopulations Into the Treatment Optimization Process for Schedules Involving Targeted Agents in Chemoradiation Therapy

Purpose: To develop a model that includes the process of resistance development into the treatment optimization process for schedules that include targeted therapies. Further, to validate the approach using clinical data and to apply the model to assess the optimal induction period with targeted agents before curative treatment with chemo-radiation in stage III lung cancer. Methods: Growth of the tumor and its subpopulations is modeled by Gompertzian growth dynamics, resistance induction as a stochastic process. Chemotherapy induced cell kill is modeled by log-cell kill dynamics, targeted agents similarly but restricted to the sensitive population. Radiation induced cell kill is assumed to follow the linear-quadratic model. The validation patient data consist of a cohort of lung cancer patients treated with tyrosine kinase inhibitors that had longitudinal imaging data available. Results: The resistance induction model was successfully validated using clinical trial data from 49 patients treated with targeted agents. The observed recurrence kinetics, with tumors progressing from 1.4–63 months, result in tumor growth equaling a median volume doubling time of 92 days [34–248] and a median fraction of pre-existing resistance of 0.035 [0–0.22], in agreement with previous clinical studies. The model revealed widely varying optimal time points for the use of curative therapy, reaching from ∼1m to >6m depending on the patient’s growth rate and amount of pre-existing resistance. This demonstrates the importance of patient-specific treatment schedules when targeted agents are incorporated into the treatment. Conclusion: We developed a model including evolutionary dynamics of resistant sub-populations with traditional chemotherapy and radiation cell kill models. Fitting to clinical data yielded patient specific growth rates and resistant fraction in agreement with previous studies. Further application of the model demonstrated how proper timing of chemo-radiation could minimize the probability of resistance, increasing tumor control significantly.

60 APPLIED LIFE SCIENCES↗

A PRELIMINARY JUPITER MODEL

In anticipation of new observational results for Jupiter's axial moment of inertia and gravitational zonal harmonic coefficients from the forthcoming Juno orbiter, we present a number of preliminary Jupiter interior models. We combine results from ab initio computer simulations of hydrogen–helium mixtures, including immiscibility calculations, with a new nonperturbative calculation of Jupiter's zonal harmonic coefficients, to derive a self-consistent model for the planet's external gravity and moment of inertia. We assume helium rain modified the interior temperature and composition profiles. Our calculation predicts zonal harmonic values to which measurements can be compared. Although some models fit the observed (pre-Juno) second- and fourth-order zonal harmonics to within their error bars, our preferred reference model predicts a fourth-order zonal harmonic whose absolute value lies above the pre-Juno error bars. This model has a dense core of about 12 Earth masses and a hydrogen–helium-rich envelope with approximately three times solar metallicity.

79 ASTRONOMY AND ASTROPHYSICS↗