Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Geometric adaptivity”

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 415 records · Page 23

Estimation of iodine nutrition and thyroid function status in late-gestation pregnant women in the United States: Development and application of a population-based pregnancy model

Previously, a deterministic biologically-based dose-response (BBDR) pregnancy model was developed to evaluate moderate thyroid axis disturbances with and without thyroid-active chemical exposure in a near-term pregnant woman and fetus. In the current study, the existing BBDR model was adapted to include a wider functional range of iodine nutrition, including more severe iodine deficiency conditions, and to incorporate empirically the effects of homeostatic mechanisms. The extended model was further developed into a population-based model and was constructed using a Monte Carlo-based probabilistic framework. In order to characterize total (T4) and free (fT4) thyroxine levels for a given iodine status at the population-level, the distribution of iodine intake for late-gestation pregnant women in the U.S was reconstructed using various reverse dosimetry methods and available biomonitoring data. The range of median (mean) iodine intake values resulting from three different methods of reverse dosimetry tested was 196.5–219.9 μg of iodine/day (228.2–392.9 μg of iodine/day). There was minimal variation in model-predicted maternal serum T4 and ft4 thyroxine levels from use of the three reconstructed distributions of iodine intake; the range of geometric mean for T4 and fT4, was 138–151.7 nmol/L and 7.9–8.7 pmol/L, respectively. The average value of the ratio of the 97.5th percentile to the 2.5th percentile equaled 3.1 and agreed well with similar estimates from recent observations in third-trimester pregnant women in the U.S. In addition, the reconstructed distributions of iodine intake allowed us to estimate nutrient inadequacy for late-gestation pregnant women in the U.S. via the probability approach. The prevalence of iodine inadequacy for third-trimester pregnant women in the U.S. was estimated to be between 21% and 44%. Taken together, the current work provides an improved tool for evaluating iodine nutritional status and the corresponding thyroid function status in pregnant women in the U.S. This model enables future assessments of the relevant risk of thyroid hormone level perturbations due to exposure to thyroid-active chemicals at the population-level. - Highlights: • A population-based thyroid function model for pregnant women was developed. • The model was used specifically study the late-gestation pregnant women in the U.S. • The prevalence of iodine inadequacy was estimated in the sub-population studied. • Developed model well predicts trimester-specific thyroid hormone reference ranges. • The model can be further used to study thyroid perturbations at a population level.

60 APPLIED LIFE SCIENCES↗

Reinforcement Learning for Anomaly Detection in Nuclear Power Plant Operation and Maintenance

In nuclear power plants (NPPs), timely identification of sensor and human errors is critical to ensure safe and efficient plant operations. Anomaly detection models can be employed for this task. However, traditional anomaly detection approaches may have high dependency on labeled datasets and struggle with adaptability in complex, dynamic environments. Reinforcement learning (RL) has demonstrated significant potential in fault diagnosis and anomaly detection; however, its application to anomaly detection in NPPs remains a relatively underexplored research direction. Hence, to address this gap, in this study, we present a novel physics-informed reinforcement learning model, PIRL-AD: Physics-Informed Reinforcement Learning for Anomaly Detection, that integrates domain knowledge from calorimetric equations into the RL framework for enhanced sensor and human error anomaly detection. We evaluate the performance of PIRL-AD against a non-physics informed RL benchmark and a support vector machine (SVM) on data collected from a forced flow loop testbed. Experimental results suggest that PIRL-AD outperforms other baselines on a range of anomalous datasets that include both sensor and human-induced anomalies across key performance metrics, statistically outperforming the RL and SVM benchmarks with respect to geometric mean (respectively, 92.96% vs. 91.06% vs. 83.01%) and F1-score (respectively, 89.23% vs. 86.98% vs. 77.01%). Furthermore, the findings suggest the potential of physics-integrated reinforcement learning models for enhanced anomaly detection performance in NPPs.

Reinforcement learning↗

SU-E-T-334: Dosimetric Impacts of Organ at Risk (OAR) Contour Errors in Prostate Volumetric Modulated Arc Therapy (VMAT) Planning: A Sensitivity Study

Purpose: To investigate the dosimetric impacts of OAR contour errors in prostate VMAT planning with a sensitivity analysis method. Methods: Ten randomly selected prostate VMAT patients were used to simulate OAR contour error in rectum and bladder. For each OAR 12 even spaced contour control points were selected in every three slices. For each simulation, only one of the control points was simulated to move inwards or outwards by up to 6 voxels and coordinates of adjacent voxels to that moving control point were moved accordingly with a 3D-spline smooth function. An in-house software was used to predict OARs dosimetric endpoints based on geometric relationship between PTV and OAR. In this study, the V75, V70, V65, V60 of rectum and V80, V75, V70, V65 of bladder with and without perturbed contours were calculated and compared. Results: The percentage of OAR dose volume difference around the reference OAR contours were plotted as iso-error lines overlaid on CT images. The significant difference was shown in OAR contours adjacent to the PTV where high dose gradient exists. When one of the 12 points adjacent to the PTV moved up to 6 voxels, maximum errors of 1.99%, 2.46%, 2.89%, and 3.19% were found in V75, V70, V65, and V60 of the rectum, and 1.59%, 1.67%, 1.87%, and 1.90% were found in V80, V75, V70, and V65 of the bladder, respectively. Conclusion: We can quantify the dosimetric impact of OAR contouring error by evaluating the percentage of OAR dose volume changes in prostate VMAT planning. For the adaptive radiation therapy, the iso-error lines can provide the planner a guideline that which portion of contour needed to be checked carefully because dosimetric sensitivity of contour errors are different for different part of contours.

60 APPLIED LIFE SCIENCES↗

Verification Suite for the Application of the Limiting Surface Density Method to Arrays of 9975 Packages

Criticality Safety Engineering at the Savannah River Site’s K-Area Complex seeks to reduce the computational burden of analysis for the K-Area Complex. By applying an adapted Limiting Surface Density hand calculation method, the need for lengthy detailed Monte Carlo models of the storage arrays may be reduced or eliminated. After promising initial results using models of idealized but heterogeneous packages, application of the method to large arrays of highly heterogeneous 9975 shipping packages is attempted. A SCALE KENO-VI model of the 9975 package, with true dimensions and radial and axial heterogeneity is developed. A suite of array arrangements representing both arrays with equal number of packages per side (package-cubic) and arrays with unequal number of packages per side (non-cubic) are evaluated in SCALE. The cubic arrangements are used to derive the 9975-specific material and geometry constants for the LSD method. These constants are then used to implement the LSD method and compare the LSD- predicted critical mass per package to SCALE-predicted critical mass per package. Sensitivities to certain material and geometric nuances are evaluated for inclusion in the verification suite. Finally, the ability of the method to predict neutron multiplication changes induced by fissile mass changes is examined as a segue for evaluating credible abnormal conditions with the LSD method.

42 ENGINEERING↗

Verification suite for the application of the limiting surface density method to arrays of 9975 packages

Criticality Safety Engineering at the Savannah River Site's K-Area Complex seeks to reduce the computational burden of analysis for the K-Area Complex. By applying an adapted Limiting Surface Density hand calculation method, the need for lengthy detailed Monte Carlo models of the storage arrays may be reduced or eliminated. After promising initial results using models of idealized but heterogeneous packages, application of the method to large arrays of highly heterogeneous 9975 shipping packages is attempted. A SCALE KENO-VI model of the 9975 package, with true dimensions and radial and axial heterogeneity is developed. A suite of array arrangements representing both arrays with equal number of packages per side (package-cubic) and arrays with unequal number of packages per side (non-cubic) are evaluated in SCALE. The cubic arrangements are used to derive the 9975-specific material and geometry constants for the LSD method. These constants are then used to implement the LSD method and compare the LSD- predicted critical mass per package to SCALE-predicted critical mass per package. Sensitivities to certain material and geometric nuances are evaluated for inclusion in the verification suite. Finally, the ability of the method to predict neutron multiplication changes induced by fissile mass changes is examined as a segue for evaluating credible abnormal conditions with the LSD method. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

WE-G-BRF-03: A Quasi-Cine CBCT Reconstruction Technique for Real-Time On- Board Target Tracking of Lung Cancer Treatment

Purpose: To develop a quasi-cine CBCT reconstruction technique that uses extremely-small angle (∼3°) projections to generate real-time high-quality lung CBCT images. Method: 4D-CBCT is obtained at the beginning and used as prior images. This study uses extremely-small angle (∼3°) on-board projections acquired at a single respiratory phase to reconstruct the CBCT image at this phase. An adaptive constrained free-form deformation (ACFD) method is developed to deform the prior 4D-CBCT volume at the same phase to reconstruct the new CBCT. Quasi-cine CBCT images are obtained by continuously reconstructing CBCT images at subsequent phases every 3° angle (∼0.5s). Note that the prior 4D-CBCT images are dynamically updated using the latest CBCT images. The 4D digital extended-cardiac-torso (XCAT) phantom was used to evaluate the efficacy of ACFD. A lung patient was simulated with a tumor baseline shift of 2mm along superior-inferior (SI) direction after every respiratory cycle for 5 cycles. Limited-angle projections were simulated for each cycle. The 4D-CBCT reconstructed by these projections were compared with the ground-truth generated in XCAT.Volume-percentage-difference (VPD) and center-of-mass-shift (COMS) were calculated between the reconstructed and the ground-truth tumors to evaluate their geometric differences.The ACFD was also compared to a principal-component-analysis based motion-modeling (MM) method. Results: Using orthogonal-view 3° projections, the VPD/COMS values for tumor baseline shifts of 2mm, 4mm, 6mm, 8mm, 10mm were 11.0%/0.3mm, 25.3%/2.7mm, 22.4%/2.9mm, 49.5%/5.4mm, 77.2%/8.1mm for the MM method, and 2.9%/0.7mm, 3.9%/0.8mm, 6.2%/1mm, 7.9%/1.2mm, 10.1%/1.1mm for the ACFD method. Using orthogonal-view 0° projections (1 projection only), the ACFD method yielded VPD/COMS results of 5.0%/0.9mm, 10.5%/1.2mm, 15.1%/1.4mm, 20.9%/1.6mm and 24.8%/1.6mm. Using single-view instead of orthogonal-view projections yielded less accurate results for ACFD. Conclusion: The ACFD method accurately reconstructs snapshot CBCT images using orthogonal-view 3° projections. It has a great potential to provide real-time quasi-cine CBCT images for verification in lung radiation therapy. The research is supported by grant from Varian Medical Systems.

60 APPLIED LIFE SCIENCES↗

Multifractal characterizations of nonstationarity and intermittency in geophysical fields: Observed, retrieved, or simulated

Geophysical data rarely show any smoothness at any scale, and this often makes comparison with theoretical model output difficult. However, highly fluctuating signals and fractal structures are typical of open dissipative systems with nonlinear dynamics, the focus of most geophysical research. High levels of variability are excited over a large range of scales by the combined actions of external forcing and internal instability. At very small scales we expect geophysical fields to be smooth, but these are rarely resolved with available instrumentation or simulation tools; nondifferentiable and even discontinuous models are therefore in order. We need methods of statistically analyzing geophysical data, whether measured in situ, remotely sensed or even generated by a computer model, that are adapted to these characteristics. An important preliminary task is to define statistically stationary features in generally nonstationary signals. We first discuss a simple criterion for stationarity in finite data streams that exhibit power law energy spectra and then, guided by developments in turbulence studies, we advocate the use of two ways of analyzing the scale dependence of statistical information: singular measures and qth order structure functions. In nonstationary situations, the approach based on singular measures seeks power law behavior in integrals over all possible scales of a nonnegative stationary field derived from the data, leading to a characterization of the intermittency in this (gradient-related) field. In contrast, the approach based on structure functions uses the signal itself, seeking power laws for the statistical moments of absolute increments over arbitrarily large scales, leading to a characterization of the prevailing nonstationarity in both quantitative and qualitative terms. We explain graphically, step by step, both multifractal statistics which are largely complementary to each other. The geometrical manifestations of nonstationarity and intermittency, 'roughness' and 'sparseness', respectively, are illustrated and the associated analytical (differentiability and continuity) properties are discussed. As an example, the two techniques are applied to a series of recent measurements of liquid water distributions inside marine stratocumulus decks; these are found to be multifractal over scales ranging from approximately 60 m to approximately 60 km. Finally, we define the 'mean multifractal plane' and show it to be a simple yet comprehensive tool with many applications including data intercomparison, (dynamical or stochastic) model and retrieval validations.

Davis, Anthony↗

Sparse-view proton computed tomography using modulated proton beams

Purpose: Proton imaging that uses a modulated proton beam and an intensity detector allows a relatively fast image acquisition compared to the imaging approach based on a trajectory tracking detector. In addition, it requires a relatively simple implementation in a conventional proton therapy equipment. The model of geometric straight ray assumed in conventional computed tomography (CT) image reconstruction is however challenged by multiple-Coulomb scattering and energy straggling in the proton imaging. Radiation dose to the patient is another important issue that has to be taken care of for practical applications. In this work, the authors have investigated iterative image reconstructions after a deconvolution of the sparsely view-sampled data to address these issues in proton CT. Methods: Proton projection images were acquired using the modulated proton beams and the EBT2 film as an intensity detector. Four electron-density cylinders representing normal soft tissues and bone were used as imaged object and scanned at 40 views that are equally separated over 360°. Digitized film images were converted to water-equivalent thickness by use of an empirically derived conversion curve. For improving the image quality, a deconvolution-based image deblurring with an empirically acquired point spread function was employed. They have implemented iterative image reconstruction algorithms such as adaptive steepest descent-projection onto convex sets (ASD-POCS), superiorization method–projection onto convex sets (SM-POCS), superiorization method–expectation maximization (SM-EM), and expectation maximization-total variation minimization (EM-TV). Performance of the four image reconstruction algorithms was analyzed and compared quantitatively via contrast-to-noise ratio (CNR) and root-mean-square-error (RMSE). Results: Objects of higher electron density have been reconstructed more accurately than those of lower density objects. The bone, for example, has been reconstructed within 1% error. EM-based algorithms produced an increased image noise and RMSE as the iteration reaches about 20, while the POCS-based algorithms showed a monotonic convergence with iterations. The ASD-POCS algorithm outperformed the others in terms of CNR, RMSE, and the accuracy of the reconstructed relative stopping power in the region of lung and soft tissues. Conclusions: The four iterative algorithms, i.e., ASD-POCS, SM-POCS, SM-EM, and EM-TV, have been developed and applied for proton CT image reconstruction. Although it still seems that the images need to be improved for practical applications to the treatment planning, proton CT imaging by use of the modulated beams in sparse-view sampling has demonstrated its feasibility.

60 APPLIED LIFE SCIENCES↗

Canopy Spectral Invariants: Application to Classification of Forest Types from Hyperspectral Data - Part 2

Many studies have been conducted to demonstrate the ability of hyperspectral data to discriminate plant dominant species. Most of them have employed the use of empirically based techniques, which are site specific, requires some initial training based on characteristics of known leaf and/or canopy spectra and therefore may not be extendable to operational use or adapted to changing or unknown land cover. In this paper we propose a physically based approach for separation of dominant forest type using hyperspectral data. The radiative transfer theory of canopy spectral invariants underlies the approach, which facilitates parameterization of the canopy reflectance in terms of the leaf spectral scattering and two spectrally invariant and structurally varying variables - recollision and directional escape probabilities. The methodology is based on the idea of retrieving spectrally invariant parameters from hyperspectral data first, and then relating their values to structural characteristics of three-dimensional canopy structure. Theoretical and empirical analyses of ground and airborne data acquired by Airborne Visible/Infrared Imaging Spectrometer (AVIRIS) over two sites in New England, USA, suggest that the canopy spectral invariants convey information about canopy structure at both the macro- and micro-scales. The total escape probability (one minus recollision probability) varies as a power function with the exponent related to the number of nested hierarchical levels present in the pixel. Its base is a geometrical mean of the local total escape probabilities and accounts for the cumulative effect of canopy structure over a wide range of scales. The ratio of the directional to the total escape probability becomes independent of the number of hierarchical levels and is a function of the canopy structure at the macro-scale such as tree spatial distribution, crown shape and size, within-crown foliage density and ground cover. These properties allow for the natural separation of dominant forest classes based on the location of points on the total escape probability vs the ratio log-log plane.

Schull, M. A.↗

SU-F-J-223: Patterns of Failure for Laryngeal Cancer Patients Treated with Definitive IMRT: Comparing Two Different Methods for Determining the Origin of Recurrence From Follow-Up PET/CT Scans

Purpose: To determine patterns of failure in laryngeal cancer treated with definitive IMRT by comparing two different methods for identifying the recurrence epicenter on follow-up PET/CT. Methods: We identified 20 patients treated for laryngeal squamous cell carcinoma with definitive IMRT who had loco-regional recurrence diagnosed on PET/CT. Recurrence PET/CT scans were co-registered with the original treatment planning CT using deformable image registration with the VoxAlign deformation engine in MIM Software. Recurrence volumes were delineated on co-registered follow-up scans using a semi-automatic PETedge tool and two separate methods were used to identify the recurrence point of origin: a) Finding the point within the recurrence volume for which the maximum distance to the surface of the surrounding recurrence volume is smaller than for any other point. b) Finding the point within the recurrence volume with the maximum standardized uptake value (SUVmax), without geometric restrictions.For each method the failure pattern was determined as whether the recurrence origin fell within the original high-dose target volumes GTV70, CTV70, PTV70 (receiving 70Gy), intermediate-risk PTV59 (receiving 59.4Gy) or low-risk PTV54 (receiving 54.1Gy), in the original treatment planning CT. Results: 23 primary/nodal recurrences from the 20 patients were analyzed. The three-dimensional distance between the two different origins was on average 10.5mm (std.dev. 10mm). Most recurrences originated in the high-dose target volumes for both methods with 13 (57%) and 11 (48%) in the GTV70 and 20 (87%) and 20 (87%) in the PTV70 for method a) and b), respectively. There was good agreement between the two methods in classifying the origin target volumes with 69% concordance for GTV70, 89% for CTV70 and 100% for PTV70. Conclusion: With strong agreement in patterns of failure between two separate methods for determining recurrence origin, we conclude that most recurrences occurred within the high-dose treatment region, which influences potential risk-adaptive treatment strategies.

60 APPLIED LIFE SCIENCES↗

SU-G-BRA-17: Tracking Multiple Targets with Independent Motion in Real-Time Using a Multi-Leaf Collimator

Purpose: Multiple targets with large intrafraction independent motion are often involved in advanced prostate, lung, abdominal, and head and neck cancer radiotherapy. Current standard of care treats these with the originally planned fields, jeopardizing the treatment outcomes. A real-time multi-leaf collimator (MLC) tracking method has been developed to address this problem for the first time. This study evaluates the geometric uncertainty of the multi-target tracking method. Methods: Four treatment scenarios are simulated based on a prostate IMAT plan to treat a moving prostate target and static pelvic node target: 1) real-time multi-target MLC tracking; 2) real-time prostate-only MLC tracking; 3) correcting for prostate interfraction motion at setup only; and 4) no motion correction. The geometric uncertainty of the treatment is assessed by the sum of the erroneously underexposed target area and overexposed healthy tissue areas for each individual target. Two patient-measured prostate trajectories of average 2 and 5 mm motion magnitude are used for simulations. Results: Real-time multi-target tracking accumulates the least uncertainty overall. As expected, it covers the static nodes similarly well as no motion correction treatment and covers the moving prostate similarly well as the real-time prostate-only tracking. Multi-target tracking reduces >90% of uncertainty for the static nodal target compared to the real-time prostate-only tracking or interfraction motion correction. For prostate target, depending on the motion trajectory which affects the uncertainty due to leaf-fitting, multi-target tracking may or may not perform better than correcting for interfraction prostate motion by shifting patient at setup, but it reduces ∼50% of uncertainty compared to no motion correction. Conclusion: The developed real-time multi-target MLC tracking can adapt for the independently moving targets better than other available treatment adaptations. This will enable PTV margin reduction to minimize health tissue toxicity while remain tumor coverage when treating advanced disease with independently moving targets involved. The authors acknowledge funding support from the Australian NHMRC Australia Fellowship and NHMRC Project Grant No. APP1042375.

60 APPLIED LIFE SCIENCES↗

Recent Progress on the Parallel Implementation of Moving-Body Overset Grid Schemes

Viscous calculations about geometrically complex bodies in which there is relative motion between component parts is one of the most computationally demanding problems facing CFD researchers today. This presentation documents results from the first two years of a CHSSI-funded effort within the U.S. Army AFDD to develop scalable dynamic overset grid methods for unsteady viscous calculations with moving-body problems. The first pan of the presentation will focus on results from OVERFLOW-D1, a parallelized moving-body overset grid scheme that employs traditional Chimera methodology. The two processes that dominate the cost of such problems are the flow solution on each component and the intergrid connectivity solution. Parallel implementations of the OVERFLOW flow solver and DCF3D connectivity software are coupled with a proposed two-part static-dynamic load balancing scheme and tested on the IBM SP and Cray T3E multi-processors. The second part of the presentation will cover some recent results from OVERFLOW-D2, a new flow solver that employs Cartesian grids with various levels of refinement, facilitating solution adaption. A study of the parallel performance of the scheme on large distributed- memory multiprocessor computer architectures will be reported.

Wissink, Andrew↗

Six Degrees-of-Freedom Prostate and Lung Tumor Motion Measurements Using Kilovoltage Intrafraction Monitoring

Purpose: Tumor positional uncertainty has been identified as a major issue that deteriorates the efficacy of radiation therapy. Tumor rotational movement, which is not well understood, can result in significant geometric and dosimetric inaccuracies. The objective of this study was to measure 6 degrees-of-freedom (6 DoF) prostate and lung tumor motion, focusing on the more novel rotation, using kilovoltage intrafraction monitoring (KIM). Methods and Materials: Continuous kilovoltage (kV) projections of tumors with gold fiducial markers were acquired during radiation therapy for 267 fractions from 10 prostate cancer patients and immediately before or after radiation therapy for 50 fractions from 3 lung cancer patients. The 6 DoF motion measurements were determined from the individual 3-dimensional (3D) marker positions, after using methods to reject spurious and smooth noisy data, using an iterative closest point algorithm. Results: There were large variations in the magnitude of the tumor rotation among different fractions and patients. Various rotational patterns were observed. The average prostate rotation angles around the left-right (LR), superior-inferior (SI), and anterior-posterior (AP) axes were 1.0 ± 5.0°, 0.6 ± 3.3°, and 0.3 ± 2.0°, respectively. For 35% of the time, the prostate rotated more than 5° about the LR axis, indicating the need for intrafractional adaptation during radiation delivery. For lung patients, the average LR, SI, and AP rotation angles were 0.8 ± 4.2°, −0.8 ± 4.5°, and 1.7 ± 3.1°, respectively. For about 30% of the time, the lung tumors rotated more than 5° around the SI axis. Respiration-induced rotation was detected in 2 of the 3 lung patients. Conclusions: The prostate and lung tumors were found to undergo rotations of more than 5° for about a third of the time. The lung tumor data represent the first 6 DoF tumor motion measured by kV images. The 6 DoF KIM method can enable rotational and translational adaptive radiation therapy and potentially reduce treatment margins.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Blue Rock CropSpanPV - Low Cost Racking for Agricultural Solar Photovoltaics

In this work, a total of four solar panel layouts, applied to Agrivoltaics, are investigated. This study presents two hypotheses: 1. Prefabrication - Lower total cost can be achieved by constructing a pre-assembled, pre-wired solar array in a factory, with automation, then rapidly deploying it in the field. 2. Tension Structures - Using tension structures in a way to suspend an overhead solar array also reduces costs when compared to conventional support methods. Designs are presented for frames to hold the panels. Configurations are advanced to connect the frames into arrays. Supports and foundations to hold the arrays in the desired layout are calculated, selected and presented. Cost models are prepared for the four layouts, which include definition of the business enterprise, the manufacturing operation, and the required facility for production at scale. Cost models include, detail material and task/labor take-offs for both manufacturing and installation. Comparisons with conventional system adapted to Agrivoltaics reveal that the two hypotheses do not hold true. A prefabricated array requires additional material that is not needed in the base system and the cost of this additional material is not outweighed by the cost savings of reduced field installation. The cost of wire rope for a tension structure does result in an economy of material, however, the cost of end connections and tensioning devices adds significantly to the cost of the overall tension structure. Foundation loads with tension structures, especially the large wind uplift seen with a fixed tilt solar array, impose significant constraints and high costs. The following conclusions are recommended for further consideration: 1. Best Case Agrivoltaic Configuration - The best value/lowest cost configuration for agrivoltaics is the installation of solar panels mounted on a single axis tracker in either a 2-in-portrait or 2-in landscape arrangement. This is useful for grazing lands and for staple crops which require high light levels. These applications represent the vast majority of the potential agrivoltaic market. Important elements in this configuration are the specific geometric layout coordinated with the various field operations and a control system linked between the farm equipment and single axis tracker. The control system is to monitor location and orientation, then actively tilt the solar panels or brake the field equipment to provide clearance and avoid collisions. 2. Long Span Applications - The use of tension structures for an overhead fixed tilt array is beneficial where long spans are a necessity, especially where the base support in the project location has consolidated rock near the surface or some other solid structure. One application is spanning wide irrigation canals. Another may be as an installation at the top deck of parking structures. 3. Direct Ground Mount Configuration – Having skids at the base of a pre-fabricated array that extend when the array is deployed, could be a successful design. This could be a useful strategy for dense, flat, ground-mount configurations like that employed by a couple of successful commercial operations (5G Maverick and Erthos). Arrays would be orientated north and south, deployed very low on prepared grade. Arrays would be aligned closely next to each other and anchored with small ground screws or simply with ballast. This same product arrangement could be very useful for rapid deployment and set-up of solar arrays for temporary deployments, disaster response and the military operations

14 SOLAR ENERGY↗

Reconnaissance of the HR 8799 exosolar system. II. Astrometry and orbital motion

We present an analysis of the orbital motion of the four substellar objects orbiting HR 8799. Our study relies on the published astrometric history of this system augmented with an epoch obtained with the Project 1640 coronagraph with an integral field spectrograph (IFS) installed at the Palomar Hale telescope. We first focus on the intricacies associated with astrometric estimation using the combination of an extreme adaptive optics system (PALM-3000), a coronagraph, and an IFS. We introduce two new algorithms. The first one retrieves the stellar focal plane position when the star is occulted by a coronagraphic stop. The second one yields precise astrometric and spectrophotometric estimates of faint point sources even when they are initially buried in the speckle noise. The second part of our paper is devoted to studying orbital motion in this system. In order to complement the orbital architectures discussed in the literature, we determine an ensemble of likely Keplerian orbits for HR 8799bcde, using a Bayesian analysis with maximally vague priors regarding the overall configuration of the system. Although the astrometric history is currently too scarce to formally rule out coplanarity, HR 8799d appears to be misaligned with respect to the most likely planes of HR 8799bce orbits. This misalignment is sufficient to question the strictly coplanar assumption made by various authors when identifying a Laplace resonance as a potential architecture. Finally, we establish a high likelihood that HR 8799de have dynamical masses below 13M{sub Jup}, using a loose dynamical survival argument based on geometric close encounters. We illustrate how future dynamical analyses will further constrain dynamical masses in the entire system.

79 ASTRONOMY AND ASTROPHYSICS↗

Finite Element Analysis of Geodesically Stiffened Cylindrical Composite Shells Using a Layerwise Theory

Layerwise finite element analyses of geodesically stiffened cylindrical shells are presented. The layerwise laminate theory of Reddy (LWTR) is developed and adapted to circular cylindrical shells. The Ritz variational method is used to develop an analytical approach for studying the buckling of simply supported geodesically stiffened shells with discrete stiffeners. This method utilizes a Lagrange multiplier technique to attach the stiffeners to the shell. The development of the layerwise shells couples a one-dimensional finite element through the thickness with a Navier solution that satisfies the boundary conditions. The buckling results from the Ritz discrete analytical method are compared with smeared buckling results and with NASA Testbed finite element results. The development of layerwise shell and beam finite elements is presented and these elements are used to perform the displacement field, stress, and first-ply failure analyses. The layerwise shell elements are used to model the shell skin and the layerwise beam elements are used to model the stiffeners. This arrangement allows the beam stiffeners to be assembled directly into the global stiffness matrix. A series of analytical studies are made to compare the response of geodesically stiffened shells as a function of loading, shell geometry, shell radii, shell laminate thickness, stiffener height, and geometric nonlinearity. Comparisons of the structural response of geodesically stiffened shells, axial and ring stiffened shells, and unstiffened shells are provided. In addition, interlaminar stress results near the stiffener intersection are presented. First-ply failure analyses for geodesically stiffened shells utilizing the Tsai-Wu failure criterion are presented for a few selected cases.

Gerhard, Craig Steven↗

Online transfer learning strategy for enhancing the scalability and deployment of deep reinforcement learning control in smart buildings

In recent years, advanced control strategies based on Deep Reinforcement Learning (DRL) proved to be effective in optimizing the management of integrated energy systems in buildings, reducing energy costs and improving indoor comfort conditions when compared to traditional reactive controllers. However, the scalability and implementation of DRL controllers are still limited since they require a considerable amount of time before converging to a near-optimal solution. This issue is currently addressed in literature through the offline pre-training of the DRL agent. However this solution results in two main critical issues: (1) the need to develop a building surrogate model to perform the training task, and (2) the need to perform a fine-tuning process over several training episodes to obtain a near-optimal control policy. In this context, this paper introduces an Online Transfer Learning (OTL) strategy that exploits two knowledge-sharing techniques, weight-initialization and imitation learning, to transfer a DRL control policy from a source office building to various target buildings in a simulation environment coupling EnergyPlus and Python. A DRL controller based on discrete Soft Actor–Critic (SAC) is trained on the source building to manage the operation of a cooling system consisting of a chiller and a thermal storage. Several target buildings are defined to benchmark the performance of the OTL strategy with that of a Rule-Based Controller (RBC) and two DRL-based control strategies, deployed in offline and online fashion. The strategy adopted for OTL emulates the real world implementation with a simulation process by implementing the transferred DRL agent for a single episode in the target buildings. Target buildings have the same geometrical features and are served by the same energy system as the source building, but differ in terms of weather conditions, electricity price schedules, occupancy patterns, and building envelope efficiency levels. The results show that the OTL strategy can reduce the cumulated sum of temperature violations on average by 50% and 80% respectively when compared to RBC and online DRL while enhancing the energy system operation with electricity cost savings ranging between 20% and 40%. Furthermore, the OTL agent performs slightly worse than the offline DRL controller but it does not require any modeling effort and can be implemented directly on target buildings emulating a real-world implementation.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Toward Effective Shell Modeling of Wrinkled Thin-Film Membranes Exhibiting Stress Concentrations

Geometrically nonlinear shell finite element analysis has recently been applied to solar-sail membrane problems in order to model the out-of-plane deformations due to structural wrinkling. Whereas certain problems lend themselves to achieving converged nonlinear solutions that compare favorably with experimental observations, solutions to tensioned membranes exhibiting high stress concentrations have been difficult to obtain even with the best nonlinear finite element codes and advanced shell element technology. In this paper, two numerical studies are presented that pave the way to improving the modeling of this class of nonlinear problems. The studies address the issues of mesh refinement and stress-concentration alleviation, and the effects of these modeling strategies on the ability to attain converged nonlinear deformations due to wrinkling. The numerical studies demonstrate that excessive mesh refinement in the regions of stress concentration may be disadvantageous to achieving wrinkled equilibrium states, causing the nonlinear solution to lock in the membrane response mode, while totally discarding the very low-energy bending response that is necessary to cause wrinkling deformation patterns. An element-level, strain-energy density criterion is suggested for facilitating automated, adaptive mesh refinements specifically aimed at the modeling of thin-film membranes undergoing wrinkling deformations.

Tessler, Alexander↗