Selection of optimal physics-based model and modeling error for a nonlinear fluid-structure interaction system.
Abstract not provided.
Bisaillon, Philippe↗
SEARCH · Engineering Papers
Results for “Structural error”
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.
Perfluoropolyethers: Development of an All-Atom Force Field for Molecular Simulations and Validation with New Experimental Vapor Pressures and Liquid Densities
A force field for perfluoropolyethers (PFPEs) based on the general optimized potentials for liquid simulations all-atom (OPLS-AA) force field has been derived in conjunction with experiments and ab initio quantum mechanical calculations. Vapor pressures and densities of two liquid PFPEs, perfluorodiglyme (CF 3 –O–(CF 2 –CF 2 –O)2–CF 3 ) and perfluorotriglyme (CF 3 –O–(CF 2 –CF 2 –O) 3 –CF 3 ), have been measured experimentally to validate the force field and increase our understanding of the physical properties of PFPEs. Force field parameters build upon those for related molecules (e.g., ethers and perfluoroalkanes) in the OPLS-AA force field, with new parameters introduced for interactions specific to PFPEs. Molecular dynamics simulations using the new force field demonstrate excellent agreement with ab initio calculations at the RHF/6-31G* level for gas-phase torsional energies (<0.5 kcal mol –1 error) and molecular structures for several PFPEs, and also accurately reproduce experimentally determined densities (<0.02 g cm –3 error) and enthalpies of vaporization derived from experimental vapor pressures (<0.3 kcal mol–1). Additional comparisons between experiment and simulation show that polyethers demonstrate a significant decrease in enthalpy of vaporization upon fluorination unlike related molecules (e.g., alkanes and alcohols). Simulation suggests this phenomenon is a result of reduced cohesion in liquid PFPEs due to a reduction in localized associations between backbone oxygen atoms and neighboring molecules.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
SU-F-T-29: The Important of Each Fraction Image-Guided Planning for Postoperative HDR-Brachytherapy in Endometrial Carcinoma
Purpose: Cylindrical applicators are often used for postoperative HDRbrachytherapy in endometrial carcinoma. It has been considered that dosimetric variation between fractions for this treatment is minimal and might not be necessary to perform treatment planning for every fractions. At our institute, it is traditional to perform treatment planning with CT simulation on the first fraction and uses this plan for the rest of treatment. This study was aim to evaluate the errors of critical structure doses between the fractions when simulation and planning were done for first fraction only. Methods: Treatment plans of 10 endometrial carcinoma patients who received postoperative HDR-brachytherapy and underwent CT-simulation for every HDR-fractions at our department were retrospectively reviewed. All of these patients were treated with cylindrical applicator and prescribed dose 15Gy in 3 fractions to 0.5cm from vaginal surface. The treatment plan from the first fraction was used to simulate in second and third CT-simulation. Radiation dose for critical structures in term of Dose-to-2cc (D2cc) were evaluated and compared between planning CT. Results: The D2cc for bladder and rectum were evaluated. For bladder, the mean error of D2cc estimation for second and third fractions was 7.6% (0.1–20.1%, SD=5.7). And the mean error for D2cc of rectum was 8.5% (0.1–29.4%, SD=8.5). Conclusion: The critical structure doses could be significant difference between fractions which may affects treatment outcomes or toxicities. From our data, image-guided brachytherapy at least with CT-Simulation should be done for every treatment fractions.
60 APPLIED LIFE SCIENCES↗
Parallel transport dynamics for mixed quantum states with applications to time-dependent density functional theory
Direct simulation of the von Neumann dynamics for a general (pure or mixed) quantum state can often be expensive. One prominent example is the real-time time-dependent density functional theory (rt-TDDFT), a widely used framework for the first principle description of many-electron dynamics in chemical and materials systems. Practical rt-TDDFT calculations often avoid the direct simulation of the von Neumann equation, and solve instead a set of Schrödinger equations, of which the dynamics is equivalent to that of the von Neumann equation. However, the time step size employed by the Schrödinger dynamics is often much smaller. Here, in order to improve the time step size and the overall efficiency of the simulation, we generalize a recent work of the parallel transport (PT) dynamics for simulating pure states [An, Lin, Multiscale Model. Simul. 18, 612, 2020] to general quantum states. The PT dynamics provides the optimal gauge choice, and can employ a time step size comparable to that of the von Neumann dynamics. Going beyond the linear and near adiabatic regime in previous studies, we find that the error of the PT dynamics can be bounded by certain commutators between Hamiltonians, density matrices, and their derived quantities. Such a commutator structure is not present in the Schrödinger dynamics. We demonstrate that the parallel transport-implicit midpoint (PT-IM) method is a suitable method for simulating the PT dynamics, especially when the spectral radius of the Hamiltonian is large. The commutator structure of the error bound, and numerical results for model rt-TDDFT calculations in both linear and nonlinear regimes, confirm the advantage of the PT dynamics.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Terrain-Influenced Winds and Fire-Fire Interactions in Wildland Fire Simulations [Dissertation]
Ensemble-based approaches to prescribed fire planning cannot be supported by computational fluid dynamics based models like FIRETEC and the Wildland-Urban Interface Fire Dynamics Simulator (WFDS) because they are too computationally expensive and cannot leverage large eddy simulation approaches like CAWFE and WRF-SFIRE because they have too coarse of resolution. QUIC-Fire was developed to fill this gap but it cannot currently address complex terrain, that is typical for instance in the Western United States. This dissertation describes a variety of improvements made to QUIC-Fire and its various incorporated algorithms in an effort to make it a viable tool in simulating wildland and prescribed fires on terrain. The modifications made to QUIC-Fire are described in three chapters. The first chapter describes the extension of the diagnostic wind model QUIC-URB, the wind engine of QUIC-Fire, to a terrain-following coordinate system. The terraininfluenced winds it generates are analyzed and compared. In particular, this chapter presents the mathematical derivation of the wind solver leading to a linear system of equations that are solved through the successive over-relaxation method. The model is validated against a standard test used in previous works (the Askervein Hill) and against a new dataset from measurements in the Socorro Mountains, New Mexico. The terrain-following implementation captures the correct phenomenology for the isolated Askervein Hill, with a wind speed up at the top of the hill. The model agrees well with measurements on the upwind side of the peak, but overestimates speed-up on the downwind side of the hill. This is due to the inability of the model to generate flow separation and wake-eddy dynamics. On a common laptop, the divergence-free wind field is obtained in 6 s, making the solver appealing for coupled fire-atmosphere simulations. The Socorro Mountain is highly complex, with many cliff faces, peaks, and valleys. Although the model captures the magnitude and direction of inlet and outlet areas of the domain, it performs rather poorly in the valley region and in the regions near the steep cliffs. Hence, the model shows good agreement with data in areas of open sloped terrain but lacks in areas where flow separation and thermally driven effects may be present (neither effect is addressed in this work). In the second chapter the implementation of the terrain-following version of QUIC-URB into QUIC-Fire, and the necessary changes needed to include terrain are described. No changes to the underlying fire spread algorithm are made other than what is required to correctly account for the inclusion of terrain. Previously published FIRETEC results that use five different topographies that share the same centerline profile are compared to simulation results from the modified QUIC-Fire that use the same topographies and fuels. QUIC-Fire results show overall similar behaviors in terms of how the topographies affect fire shapes and trends in spread rates. Due to the terrain-following version of QUIC-URB being unable to generate flow separations at the crest of hills, fire spread rates in these regions across all topographies are over-predicted when compared to FIRETEC. Lateral fire growth shows similar trends with FIRETEC between topographies but does not capture the increase in spread due to a diagonal interface between grassland and forested fuel region of the domain. These results suggest that there are three algorithms within QUIC-Fire that could use improvement: how flame tilt angle is accounted for, the incorporation of non-local drag effects, and the inclusion of the wake-eddy parameterizations that are used in QUIC-URB. Lastly, the third chapter describes a modification to the initial guess used for the calculation of the QUIC-URB mass-conserved wind solution during fire simulations. The modification is aimed at improving fire-fire interactions in QUIC-Fire simulations. The modification consists of using the solution from the previous timestep as the starting point for the calculation of the solution for the next timestep. Fire-fire interactions is greatly improved by the change but a new source of error is introduced. Due to how plumes are modelled in QUIC-Fire the new solution contains errors where gaps in the plume structure are present. However, these errors are mostly limited to the upper atmosphere, where they do not affect fire behavior at the surface, and their magnitude isn’t significant enough to discount the amount of new fire phenomenology now captured in QUIC-Fire with the change.
58 GEOSCIENCES↗
Moving beyond the total sea ice extent in gauging model biases
Here, reproducing characteristics of observed sea ice extent remains an important climate modeling challenge. This study describes several approaches to improve how model biases in total sea ice distribution are quantified, and applies them to historically forced simulations contributed to phase 5 of the Coupled Model Intercomparison Project (CMIP5). The quantity of hemispheric total sea ice area, or some measure of its equatorward extent, is often used to evaluate model performance. A new approach is introduced that investigates additional details about the structure of model errors, with an aim to reduce the potential impact of compensating errors when gauging differences between simulated and observed sea ice. Using multiple observational datasets, several new methods are applied to evaluate the climatological spatial distribution and the annual cycle of sea ice cover in 41 CMIP5 models. It is shown that in some models, error compensation can be substantial, for example resulting from too much sea ice in one region and too little in another. Error compensation tends to be larger in models that agree more closely with the observed total sea ice area, which may result from model tuning. The results herein suggest that consideration of only the total hemispheric sea ice area or extent can be misleading when quantitatively comparing how well models agree with observations. Further work is needed to fully develop robust methods to holistically evaluate the ability of models to capture the finescale structure of sea ice characteristics; however, the “sector scale” metric used here aids in reducing the impact of compensating errors in hemispheric integrals.
54 ENVIRONMENTAL SCIENCES↗
Rate constants for fine-structure excitations in O-H collisions with error bars obtained by machine learning
We present an approach using a combination of coupled channel scattering calculations with a machine-learning technique based on Gaussian Process regression to determine the sensitivity of the rate constants for non-adiabatic transitions in inelastic atomic collisions to variations of the underlying adiabatic interaction potentials. Using this approach, we improve the previous computations of the rate constants for the fine-structure transitions in collisions of O({sup 3}P{sub j}) with atomic H. We compute the error bars of the rate constants corresponding to 20% variations of the ab initio potentials and show that this method can be used to determine which of the individual adiabatic potentials are more or less important for the outcome of different fine-structure changing collisions.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Auxiliary-field quantum Monte Carlo calculations of the structural properties of nickel oxide
Auxiliary-field quantum Monte Carlo (AFQMC) has repeatedly shown itself as one of the most accurate quantum many-body methods, capable of simulating both real and model systems. In this work, we investigate the application of AFQMC to realistic strongly correlated materials in periodic Gaussian basis sets. Using nickel oxide (NiO) as an example, we investigate the importance of finite size effects and basis set errors on the structural properties of the correlated solid. We propose benchmark calculations for NiO and compare our results to both experimental measurements and existing theoretical methods.
75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
Supporting the Development of Soft-Error Resilient Message Passing Applications using Simulation
Radiation-induced bit flip faults are of particular concern in extreme-scale high-performance computing systems. This paper presents a simulation-based tool that enables the development of soft-error resilient message passing applications by permitting the investigation of their correctness and performance under various fault conditions. The documented extensions to the Extreme-scale Simulator (xSim) enable the injection of bit flip faults at specific of injection location(s) and fault activation time(s), while supporting a significant degree of configurability of the fault type. Experiments show that the simulation overhead with the new feature is ~2,325% for serial execution and ~1,730% at 128 MPI processes, both with very fine-grain fault injection. Fault injection experiments demonstrate the usefulness of the new feature by injecting bit flips in the input and output matrices of a matrix-matrix multiply application, revealing vulnerability of data structures, masking and error propagation. xSim is the very first simulation-based MPI performance tool that supports both, the injection of process failures and bit flip faults.
Engelmann, Christian↗
AAPM Truth‐based CT (TrueCT) reconstruction grand challenge
Background: This Special Report summarizes the 2022, AAPM grand challenge on Truth-based CT image reconstruction. Purpose: To provide an objective framework for evaluating CT reconstruction methods using virtual imaging resources consisting of a library of simulated CT projection images of a population of human models with various diseases. Methods: Two hundred unique anthropomorphic, computational models were created with varied diseases consisting of 67 emphysema, 67 lung lesions, and 66 liver lesions. The organs were modeled based on clinical CT images of real patients. The emphysematous regions were modeled using segmentations from patient CT cases in the COPDGene Phase I dataset. For the lung and liver lesion cases, 1–6 malignant lesions were created and inserted into the human models, with lesion diameters ranging from 5.6 to 21.9 mm for lung lesions and 3.9 to 14.9 mm for liver lesions. The contrast defined between the liver lesions and liver parenchyma was 82 ± 12 HU, ranging from 50 to 110 HU. Similarly, the contrast between the lung lesions and the lung parenchyma was defined as 781 ± 11 HU, ranging from 725 to 805 HU. For the emphysematous regions, the defined HU values were −950 ± 17 HU ranging from −918 to −979 HU. The developed human models were imaged with a validated CT simulator. The resulting CT sinograms were shared with the participants. The participants reconstructed CT images from the sinograms and sent back their reconstructed images. Further, the reconstructed images were then scored by comparing the results against the corresponding ground truth values. The scores included both task-generic (root mean square error [RMSE] and structural similarity matrix [SSIM]), and task-specific (detectability index [d’] and lesion volume accuracy) metrics. For the cases with multiple lesions, the measured metric was averaged across all the lesions. To combine the metrics with each other, each metric was normalized to a range of 0 to 1 per disease type, with “0” and “1” being the worst and best measured values across all cases of the disease type for all received reconstructions. Results: The True-CT challenge attracted 52 participants, out of which 5 successfully completed the challenge and submitted the requested 200 reconstructions. Across all participants and disease types, SSIM absolute values ranged from 0.22 to 0.90, RMSE from 77.6 to 490.5 HU, d’ from 0.1 to 64.6, and volume accuracy ranged from 1.2 to 753.1 mm3. The overall scores demonstrated that participant “A” had the best performance in all categories, except for the metrics of d’ for lung lesions and RMSE for liver lesions. Participant “A” had an average normalized score of 0.41 ± 0.22, 0.48 ± 0.32, and 0.42 ± 0.33 for the emphysema, lung lesion, and liver lesion cases, respectively. Conclusions: The True-CT challenge successfully enabled objective assessment of CT reconstructions with the unique advantage of access to a diverse population of diseased human models with known ground truth. This study highlights the significant potential of virtual imaging trials in objective assessment of medical imaging technologies.
60 APPLIED LIFE SCIENCES↗
OC6 Phase Ib: Validation of the CFD predictions of difference-frequency wave excitation on a FOWT semisubmersible
During the previous OC5 11 project, state-of-the-art mid-fidelity engineering tools for floating wind systems were found to consistently underpredict the nonlinear, low-frequency responses of semisubmersible offshore wind platforms, leading to substantial errors in the structural loads because of low-frequency surge and pitch resonance. To examine this underprediction, a coordinated investigation with computational fluid dynamics (CFD) simulations and model-basin experiments was carried out. Both investigations involved a fixed and simplified OC5-DeepCwind semisubmersible in bichromatic waves. The wave excitations - especially the nonlinear, difference-frequency excitation - on the structure from the CFD simulations were compared to the experimental measurements for validation, with uncertainty analyses for both the experimental and the CFD results. Furthermore, the wave excitations on each column of the semisubmersible were measured separately in the experiment, allowing the validation of the CFD results to be done on a per-column basis. Overall, the CFD predictions of the difference-frequency excitations agree with the experimental measurements, suggesting the CFD solutions can be used as a reference for tuning and improving the engineering-level tools and can provide a means to better understand the underprediction at low frequencies.
16 TIDAL AND WAVE POWER↗
Impact of Lateral Flow on Surface Water and Energy Budgets over the Southern Great Plains – A Modeling Study
As we see the horizontal grid spacing decrease, treatment of hydrologic processes in land surface models, such as the lateral flow of surface and subsurface flow, need to be explicitly represented. Unlike previous studies that mainly focused on the mountainous regions, in this study the offline WRF-Hydro model is employed to study the impact of lateral flow on soil moisture and energy fluxes over the relatively flat southern Great Plains (SGP). The vast amount of measurements over the SGP provide an unique opportunity to assess the model behavior. In addition, newly developed land surface properties and input forcing are ingested into the model, in an attempt to reduce uncertainties associated with the initial and boundary forcing and help to identify model deficiencies. Our results show that the more realistic inputs (parameters, soil types, forcing) lead to larger underestimation of latent heat flux and dry bias, indicating the existence of model structural uncertainty (embedded errors) in WRF-Hydro that need to be characterized to inform future model development efforts. Including lateral flow processes partly mitigates the model deficiencies in representing hydrologic processes and alleviates the dry bias. In particular, both surface and subsurface lateral flow increase soil moisture mainly over the lower elevations, except that subsurface flow also affects soil moisture over steeper terrains. Additional simulations are performed to assess the effect of routing resolution on model results. When LSM resolution is high, noticeable differences in soil moisture are produced between different routing resolutions especially over steep terrain. Whereas when LSM resolution is coarse, differences between routing resolutions become negligible, especially over flat terrain.
54 ENVIRONMENTAL SCIENCES↗
Fusion of Multiple Models for Improving Gross Primary Production Estimation With Eddy Covariance Data Based on Machine Learning
Abstract Terrestrial gross primary production (GPP) represents the magnitude of CO 2 uptake through vegetation photosynthesis, and is a key variable for carbon cycles between the biosphere and atmosphere. Light use efficiency (LUE) models have been widely used to estimate GPP for its physiological mechanisms and availability of data acquisition and implementation, yet each individual GPP model has exhibited large uncertainties due to input errors and model structure, and further studies of systematic validation, comparison, and fusion of those models with eddy covariance (EC) site data across diverse ecosystem types are still needed in order to further improve GPP estimation. We here compared and fused five GPP models (VPM, EC‐LUE, GOL‐PEM, CHJ, and C‐Fix) across eight ecosystems based on FLUXNET2015 data set using the ensemble methods of Bayesian Model Averaging (BMA), Support Vector Machine (SVM), and Random Forest (RF) separately. Our results showed that for individual models, EC‐LUE gave a better performance to capture interannual variability of GPP than other models, followed by VPM and GLO‐PEM, while CHJ and C‐Fix were more limited in their estimation performance. We found RF and SVM were superior to BMA on merging individual models at various plant functional types (PFTs) and at the scale of individual sites. On the basis of individual models, the fusion methods of BMA, SVM, and RF were examined by a five‐fold cross validation for each ecosystem type, and each method successfully improved the average accuracy of estimation by 8%, 18%, and 19%, respectively.
Environmental Sciences & Ecology↗
Validation of the stochastic inversion algorithm for acoustic travel-time tomography: a large eddy simulation study
Acoustic tomography (AT) is explored as a remote sensing technique to obtain instantaneous snapshots of temperature and velocity fluctuations for wind energy applications. This study integrates Large Eddy Simulation (LES) with the Stochastic Inversion (SI) method to validate the algorithm’s capacity for accurate reconstruction of atmospheric fluctuations. The initial findings demonstrate the efficacy of the method in accurately capturing the predominant flow structures. Normalized L2 error evaluations further inform the algorithm’s precision, with errors accentuated in less sampled peripheral regions. The results underscore the method’s promise as a non-intrusive observational tool, with ongoing development poised to improve its precision and reliability.
17 WIND ENERGY↗
Microscopic theory of onset of decaging and bond-breaking activated dynamics in ultradense fluids with strong short-range attractions
We theoretically study thermally activated “in cage” elementary dynamical processes that precede full structural relaxation in ultra-dense particle liquids interacting via strong short range attractive forces. The analysis is based on a microscopic theory formulated at the particle trajectory level built on the dynamic free energy concept and an explicit treatment of how attractive forces control the formation and lifetime of physical bonds. Mean time scales for bond breaking, the early stage of cage escape, and non-Fickian displacement by a fixed amount are analyzed in the repulsive glass, bonded repulsive (attractive) glass, fluid, and dense gel regimes. The theory predicts a strong length-scale-dependent growth of these time scales with attractive force strength at fixed packing fraction, a much weaker slowing down with density at fixed attraction strength, and a strong decoupling of the shorter bond breaking time with the other two time scales that are controlled mainly by perturbed steric caging. All results are in good accord with simulations, and additional testable predictions are made. In conclusion, the classic statistical mechanical projection approximation of replacing all bare attractive and repulsive forces with a single effective force determined by pair structure incurs major errors for describing processes associated with thermally activated escape from transiently localized states.
36 MATERIALS SCIENCE↗
Structure Factors of Neutron Matter at Finite Temperature
Here, we compute continuum and infinite volume limit extrapolations of the structure factors of neutron matter at finite temperature and density. Using a lattice formulation of leading-order pionless effective field theory, we compute the momentum dependence of the structure factors at finite temperature and at densities beyond the reach of the virial expansion. The Tan contact parameter is computed and the result agrees with the high momentum tail of the vector structure factor. All errors, statistical and systematic, are controlled for. This calculation is a first step towards a model-independent understanding of the linear response of neutron matter at finite temperature.
73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
Optimal Experimental Design for Inverse Problems in the Presence of Observation Correlations
Optimal experimental design (OED) is the general formalism of sensor placement and decisions on the data collection strategy for engineered or natural experiments. This problem is prevalent in many critical fields such as battery design, numerical weather prediction, geosciences, environmental and urban studies. State-of-the-art computational methods for experimental design do not accommodate correlation structure in observational errors produced by many expensive-to-operate devices such as X-ray machines, radars, and satellites. Discarding evident data correlations leads to biased results, higher expenses, and waste of valuable resources. We present a general formulation of the OED formalism for model-constrained large-scale Bayesian linear inverse problems, where measurement errors are generally correlated. The proposed approach utilizes the Hadamard product of matrices to formulate the weighted-likelihood, and is valid for both finite as well as infinite-dimensional Bayesian inverse problems. Extensive numerical experiments are carried out for empirical verification of the proposed approach using an advection-diffusion model, where the objective is to optimally place a small set of sensors, under a limited budget, to predict the concentration of a contaminant in a closed and bounded domain.
Correlated observations↗
Deep machine learning-assisted multiphoton microscopy to reduce light exposure and expedite imaging
Two-photon excitation fluorescence (2PEF) allows imaging of tissue up to about one millimeter in thickness. Typically, reducing fluorescence excitation exposure reduces the quality of the image. However, using deep learning super resolution techniques, these low-resolution images can be converted to high-resolution images. This work explores improving human tissue imaging by applying deep learning to maximize image quality while reducing fluorescence excitation exposure. We analyze two methods: a method based on U-Net, and a patch-based regression method. Both methods are evaluated on a skin dataset and an eye dataset. The eye dataset includes 1200 paired high power and low power images of retinal organoids. The skin dataset contains multiple frames of each sample of human skin. High-resolution images were formed by averaging 70 frames for each sample and low-resolution images were formed by averaging the first 7 and 15 frames for each sample. The skin dataset includes 550 images for each of the resolution levels. We track two measures of performance for the two methods: mean squared error (MSE) and structural similarity index measure (SSIM). For the eye dataset, the patches method achieves an average MSE of 27,611 compared to 146,855 for the U-Net method, and an average SSIM of 0.636 compared to 0.607 for the U-Net method. For the skin dataset, the patches method achieves an average MSE of 3.768 compared to 4.032 for the U-Net method, and an average SSIM of 0.824 compared to 0.783 for the U-Net method. Despite better performance on image quality, the patches method is worse than the U-Net method when comparing the speed of prediction, taking 303 seconds to predict one image compared to less than one second for the U-Net method.
62 RADIOLOGY AND NUCLEAR MEDICINE↗