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 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Overcoming geometric limitations in metallic glasses through stretch blow molding

Bulk metallic glasses (BMGs) exhibit remarkable mechanical properties, such as high strength and elasticity, which is often paired with fracture toughness. Their supercooled liquid region gives rise to plastic-like processing and suggests parts and shapes that can otherwise not be obtained for crystalline metals. However, current processing techniques only allow for limited options in terms of geometry, thicknesses uniformity, and shape complexity. Here we introduce a new processing technique, “stretch blow molding,” to expand the range of possible parts and increase the available geometries that can be fabricated with BMGs. Additionally, a model is derived that allows for the quantification and prediction of stretch blow molding and provides insight into its potential use and limitations. We demonstrate that with stretch blow molding overall strains exceeding 2000% are achievable, compared to the previously reported ~150% of blow molding. With the ability to stretch blow mold shapes that were previously unachievable with any other metal fabrication technique in a fast and economical manner, and the superb properties of BMGs, we look forward to a broad commercial adaptation of this technique.

36 MATERIALS SCIENCE↗
NASA NTRS1989

Geometric transformations for video compression and human teleoperator display

A method for bandwidth-efficient processing of video imagery to be viewed by the teleoperator of a remotely-operated vehicle on which the camera is mounted is described. The method comprises image coding, transmission, and reconstruction. It is assumed that the transmission bandpass is the limiting factor rather than encoding/decoding schemata; that image coding and reconstruction will be done within the general abilities of the NASA/TI Programmable Remapper; and that the ratio of retained local detail to the operator's visual resolution is held constant throughout the large-field image that is seen. Novel features include that the compression and reconstruction address certain characteristics of the human visual system, that two-way communication controls a moving 'fovea' in the transformation, and that resolution varies over the image. Conventional motivations accommodated include the Cartesian raster-scan nature of available imagers and display devices and a need for low bandwidth in the image transmission. Unique image processing hardware, NASA's Programmable Remapper, allows demonstration of the method. Once refined, the technology could be adapted to special purpose imagers and display devices, or otherwise to dedicated image processing hardware.

Juday, Richard D.↗
DOE OSTI2016

SU-F-J-186: Enabling Adaptive IMPT with CBCT-Based Dose Recalculation for H&N and Prostate Cancer Patients

Purpose: To enable adaptive intensity modulated proton therapy for sites sensitive to inter-fractional changes on the basis of accurate CBCT-based proton dose calculations. To this aim two CBCT intensity correction methods are considered: planning CT (pCT) to CBCT DIR and projection correction based on pCT DIR prior. Methods: 3 H&N and 3 prostate cancer patients with CBCT images and corresponding projections were used in this study, in addition to pCT and re-planning CT (rpCT) images (H&N only). A virtual CT (vCT) was generated by pCT to CBCT DIR. In a second approach, the vCT was used as prior for scatter correction of the CBCT projections to yield a CBCTcor image. BEV 2D range maps of SFUD IMPT plans were compared. For the prostate cases, the geometric accuracy of the vCT was also evaluated by contour comparison to physician delineation of the CBCTcor and original CBCT. Results: SFUD dose calculations on vCT and CBCTcor were found to be within 3mm for 97% to 99% of 2D range maps. Median range differences compared to rpCT were below 0.5mm. Analysis showed that the DIR-based vCT approach exhibits inaccuracies in the pelvic region due to the very low soft-tissue contrast in the CBCT. The CBCTcor approach yielded results closer to the original CBCT in terms of DICE coefficients than the vCT (median 0.91 vs 0.81) for targets and OARs. In general, the CBCTcor approach was less affected by inaccuracies of the DIR used during the generation of the vCT prior. Conclusion: Both techniques yield 3D CBCT images with intensities equivalent to diagnostic CT and appear suitable for IMPT dose calculation for most sites. For H&N cases, no considerable differences between the two techniques were found, while improved results of the CBCTcor were observed for pelvic cases due to the reduced sensitivity to registration inaccuracies. Deutsche Forschungsgemeinschaft (MAP); Bundesministerium fur Bildung und Forschung (01IB13001)

60 APPLIED LIFE SCIENCES↗
DOE OSTI2024

Algorithm 1049: The Delaunay Density Diagnostic

Accurate approximation of a real-valued function depends on two aspects of the available data: the density of inputs within the domain of interest and the variation of the outputs over that domain. There are few methods for assessing whether the density of inputs is sufficient to identify the relevant variations in outputs—i.e., the “geometric scale” of the function—despite the fact that sampling density is closely tied to the success or failure of an approximation method. In this article, we introduce a general purpose, computational approach to detecting the geometric scale of real-valued functions over a fixed domain using a deterministic interpolation technique from computational geometry. The algorithm is intended to work on scalar data in moderate dimensions (2–10). Our algorithm is based on the observation that a sequence of piecewise linear interpolants will converge to a continuous function at a quadratic rate (in L 2 norm) if and only if the data are sampled densely enough to distinguish the feature from noise (assuming sufficiently regular sampling). We present numerical experiments demonstrating how our method can identify feature scale, estimate uncertainty in feature scale, and assess the sampling density for fixed (i.e., static) datasets of input–output pairs. Finally, we include analytical results in support of our numerical findings and have released lightweight code that can be adapted for use in a variety of data science settings.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

High-speed 3D shape measurement using the optimized composite fringe patterns and stereo-assisted structured light system

In this study, we propose a high-speed 3D shape measurement technique based on the optimized composite fringe patterns and stereo-assisted structured light system. Stereo phase unwrapping, as a new-fashioned method for absolute phase retrieval based on the multi-view geometric constraints, can eliminate the phase ambiguities and obtain a continuous phase map without projecting any additional patterns. However, in order to ensure the stability of phase unwrapping, the period of fringe is generally around 20, which limits the accuracy of 3D measurement. To solve this problem, we develop an optimized method for designing the composite pattern, in which the speckle pattern is embedded into the conventional 4-step phase-shifting fringe patterns without compromising the fringe modulation, and thus the phase measurement accuracy. We also present a simple and effective evaluation criterion for the correlation quality of the designed speckle pattern in order to improve the matching accuracy significantly. When the embedded speckle pattern is demodulated, the periodic ambiguities in the wrapped phase can be eliminated by combining the adaptive window image correlation with geometry constraint. Finally, some mismatched regions are further corrected based on the proposed regional diffusion compensation technique (RDC). These proposed techniques constitute a complete computational framework that allows to effectively recover an accurate, unambiguous, and distortion-free 3D point cloud with only 4 projected patterns. Experimental results verify that our method can achieve high-speed, high-accuracy, robust 3D shape measurement with dense (64-period) fringe patterns at 5000 frames per second.

42 ENGINEERING↗
DOE OSTI2012

A Bayesian nonrigid registration method to enhance intraoperative target definition in image-guided prostate procedures through uncertainty characterization

Purpose: This study introduces a probabilistic nonrigid registration method for use in image-guided prostate brachytherapy. Intraoperative imaging for prostate procedures, usually transrectal ultrasound (TRUS), is typically inferior to diagnostic-quality imaging of the pelvis such as endorectal magnetic resonance imaging (MRI). MR images contain superior detail of the prostate boundaries and provide substructure features not otherwise visible. Previous efforts to register diagnostic prostate images with the intraoperative coordinate system have been deterministic and did not offer a measure of the registration uncertainty. The authors developed a Bayesian registration method to estimate the posterior distribution on deformations and provide a case-specific measure of the associated registration uncertainty. Methods: The authors adapted a biomechanical-based probabilistic nonrigid method to register diagnostic to intraoperative images by aligning a physician's segmentations of the prostate in the two images. The posterior distribution was characterized with a Markov Chain Monte Carlo method; the maximum a posteriori deformation and the associated uncertainty were estimated from the collection of deformation samples drawn from the posterior distribution. The authors validated the registration method using a dataset created from ten patients with MRI-guided prostate biopsies who had both diagnostic and intraprocedural 3 Tesla MRI scans. The accuracy and precision of the estimated posterior distribution on deformations were evaluated from two predictive distance distributions: between the deformed central zone-peripheral zone (CZ-PZ) interface and the physician-labeled interface, and based on physician-defined landmarks. Geometric margins on the registration of the prostate's peripheral zone were determined from the posterior predictive distance to the CZ-PZ interface separately for the base, mid-gland, and apical regions of the prostate. Results: The authors observed variation in the shape and volume of the segmented prostate in diagnostic and intraprocedural images. The probabilistic method allowed us to convey registration results in terms of posterior distributions, with the dispersion providing a patient-specific estimate of the registration uncertainty. The median of the predictive distance distribution between the deformed prostate boundary and the segmented boundary was Less-Than-Or-Slanted-Equal-To 3 mm (95th percentiles within {+-}4 mm) for all ten patients. The accuracy and precision of the internal deformation was evaluated by comparing the posterior predictive distance distribution for the CZ-PZ interface for each patient, with the median distance ranging from -0.6 to 2.4 mm. Posterior predictive distances between naturally occurring landmarks showed registration errors of Less-Than-Or-Slanted-Equal-To 5 mm in any direction. The uncertainty was not a global measure, but instead was local and varied throughout the registration region. Registration uncertainties were largest in the apical region of the prostate. Conclusions: Using a Bayesian nonrigid registration method, the authors determined the posterior distribution on deformations between diagnostic and intraprocedural MR images and quantified the uncertainty in the registration results. The feasibility of this approach was tested and results were positive. The probabilistic framework allows us to evaluate both patient-specific and location-specific estimates of the uncertainty in the registration result. Although the framework was tested on MR-guided procedures, the preliminary results suggest that it may be applied to TRUS-guided procedures as well, where the addition of diagnostic MR information may have a larger impact on target definition and clinical guidance.

60 APPLIED LIFE SCIENCES↗
NASA NTRS1997

Aerodynamic Shape Sensitivity Analysis and Design Optimization of Complex Configurations Using Unstructured Grids

A three-dimensional unstructured grid approach to aerodynamic shape sensitivity analysis and design optimization has been developed and is extended to model geometrically complex configurations. The advantage of unstructured grids (when compared with a structured-grid approach) is their inherent ability to discretize irregularly shaped domains with greater efficiency and less effort. Hence, this approach is ideally suited for geometrically complex configurations of practical interest. In this work the nonlinear Euler equations are solved using an upwind, cell-centered, finite-volume scheme. The discrete, linearized systems which result from this scheme are solved iteratively by a preconditioned conjugate-gradient-like algorithm known as GMRES for the two-dimensional geometry and a Gauss-Seidel algorithm for the three-dimensional; similar procedures are used to solve the accompanying linear aerodynamic sensitivity equations in incremental iterative form. As shown, this particular form of the sensitivity equation makes large-scale gradient-based aerodynamic optimization possible by taking advantage of memory efficient methods to construct exact Jacobian matrix-vector products. Simple parameterization techniques are utilized for demonstrative purposes. Once the surface has been deformed, the unstructured grid is adapted by considering the mesh as a system of interconnected springs. Grid sensitivities are obtained by differentiating the surface parameterization and the grid adaptation algorithms with ADIFOR (which is an advanced automatic-differentiation software tool). To demonstrate the ability of this procedure to analyze and design complex configurations of practical interest, the sensitivity analysis and shape optimization has been performed for a two-dimensional high-lift multielement airfoil and for a three-dimensional Boeing 747-200 aircraft.

Taylor, Arthur C., III↗
NASA NTRS2002

Science opportunity analyzer - a multi-mission tool for planning

For many years the diverse scientific community that supports JPL's wide variety ofinterplanetary space missions has needed a tool in order to plan and develop their experiments. The tool needs to be easily adapted to various mission types and portable to the user community. The Science Opportunity Analyzer, SOA, now in its third year of development, is intended to meet this need. SOA is a java-based application that is designed to enable scientists to identify and analyze opportunities for science observations from spacecraft. It differs from other planning tools in that it does not require an in-depth knowledge of the spacecraft command system or operation modes to begin high level planning. Users can, however, develop increasingly detailed levels of design. SOA consists of six major functions: Opportunity Search, Visualization, Observation Design, Constraint Checking, Data Output and Communications. Opportunity Search is a GUI driven interface to existing search engines that can be used to identify times when a spacecraft is in a specific geometrical relationship with other bodies in the solar system. This function can be used for advanced mission planning as well as for making last minute adjustments to mission sequences in response to trajectory modifications. Visualization is a key aspect of SOA. The user can view observation opportunities in either a 3D representation or as a 2D map projection. The user is given extensive flexibility to customize what is displayed in the view. Observation Design allows the user to orient the spacecraft and visualize the projection of the instrument field of view for that orientation using the same views as Opportunity Search. Constraint Checking is provided to validate various geometrical and physical aspects of an observation design. The user has the ability to easily create custom rules or to use official project-generated flight rules. This capability may also allow scientists to easily impact the cost to science if flight rule changes occur. Data Output generates information based on the spacecraft's trajectory, opportunity search results or based on a created observation. The data can be viewed either in tabular format or as a graph. Finally, SOA is unique in that it is designed to be able to communicate with a variety of existing planning and sequencing tools. From the very beginning SOA was designed with the user in mind. Extensive surveys of the potential user community were conducted in order to develop the software requirements. Throughout the development period, close ties have been maintained with the science community to insure that the tool maintains its user focus. Although development is still in its early stages, SOA is already developing a user community on the Cassini project, which is depending on this tool for their science planning. There are other tools at JPL that do various pieces of what SOA can do; however, there is no other tool which combines all these functions and presents them to the user in such a convenient, cohesive, and easy to use fashion.

SOA science planning mission operations sequence s↗
NASA NTRS2003

Introduction to the Abrupt Wing Stall (AWS) Program

The Abrupt Wing Stall (AWS) Program has addressed the problem of uncommanded, transonic lateral motions, such as wing drop, with experimental, computational, and simulation tools. Background to the establishment of the AWS program is given as well as program objectives. In order to understand the fundamental flow mechanisms that caused the undesirable motions for a pre-production version of the F/A-18E, steady and unsteady flow field details were gathered from dedicated transonic wind-tunnel testing and computational studies. The AWS program has also adapted a free-to-roll (FTR) wind-tunnel testing technique traditionally used for low-speed studies of lateral dynamic stability to the transonic flow regime. This FTR capability was demonstrated first in a proof-of -concept study and then applied to an assessment of four different aircraft configurations. Figures of merit for static testing and for FTR testing have been evaluated for two configurations that demonstrated wing drop susceptibility during full-scale flight conditions (the pre-production F/A-18E and the AV-8B at the extremes of its flight envelope) and two configurations that do not exhibit wing drop (the F/A-18C and the F-16C). Design insights have been obtained from aerodynamic computational studies of the four aircraft configurations and from computations quantifying the impact of the various geometric wing differences between the F/A-18C and the F/A-18E wings. Finally, the AWS program provides guidance for assessing, in the simulator, the impact of experimentally determined lateral activity on flight characteristics before going to flight.

Hall, Robert M.↗
NASA NTRS2003

Introduction to the Abrupt Wing Stall (AWS) Program

The Abrupt Wing Stall (AWS) Program has addressed the problem of uncommanded, transonic lateral motions, such as wing drop, with experimental, computational, and simulation tools. Background to the establishment of the AWS program is given as well as program objectives. In order to understand the fundamental flow mechanisms that caused the undesirable motions for a pre-production version of the F/A-18E, steady and unsteady flow field details were gathered from dedicated transonic wind-tunnel testing and computational studies. The AWS program has also adapted a free-to- roll (FTR) wind-tunnel testing technique traditionally used for low-speed studies of lateral dynamic stability to the transonic flow regime. This FTR capability was demonstrated first in a proof-of-concept study and then applied to an assessment of four different aircraft configurations. Figures of merit for static testing and for FTR testing have been evaluated for two configurations that demonstrated wing drop susceptibility during full-scale flight conditions (the pre-production F/A-l8E and the AV-8B at the extremes of its flight envelope) and two configurations that do not exhibit wing drop (the F/A-l8C and the F-l6C). Design insights have been obtained from aerodynamic computational studies of the four aircraft configurations and from computations quantifying the impact of the various geometric wing differences between the F/A-18C and the F/A-l8E wings. Finally, the AWS program provides guidance for assessing, in the simulator, the impact of experimentally determined lateral activity on flight characteristics before going to flight.

Hall, Robert M.↗
DOE OSTI2019

Probability of Pinching Failure in PWR SNF after Long-Term Storage

A probabilistic assessment was performed on pinching failure of PWR SNF rod cladding under accidental SNF cask transport conditions. a set of 3,000 rod pinching scenarios was analyzed, with different geometrical, material, and loading parameters, using a Latin Hypercube Sampling (LHS). The investigated accidents include train derailments leading to cask collisions or cask drops. Pinching failure was related to cladding offset strain, defined as the plastic deformation interval in radial cladding direction normalized with respect to the outer cladding diameter. Finite element (FE) models were used to compute SNF rod deformations and the offset strain demand OSD in each of the assessed pinching scenarios. The numerical models account for the contribution of the degraded fuel pellets to rod pinching stiffness. The offset strain capacity OSC was computed using regression and interpolation models based on ring compression (RCT) test results, and account for hydride-related effects on cladding ductility. The results show a low probability of pinching failure. Hydride-related cladding embrittlement affects the structural integrity of PWR SNF rods after long-term dry storage of up to 300 years. Radial hydride patterns or large amounts of circumferentially oriented hydrides reduce the cladding ductile deformation capability when pinched at low material temperatures. The probability of hydride-related pinching failure of PWR cladding is estimated after long-term dry storage to address: How do PWR SNF rod conditions like hydrogen content, vacuum drying peak cladding hoop stress σ{sub θ}, and cladding material temperature affect cladding ductility? Is hydride-related cladding embrittlement a material degradation mechanism that affects the structural integrity of PWR SNF rods under pinching loading? What are the controlling parameters for deformation demands in PWR SNF rods under pinching loading? Is the pinching load amplitude an important parameter that controls the probability of pinching failure? How do parameters like pellet-cladding gap or fuel pellet degradation affect the pinching loading capacity of PWR SNF rods? Statistical modelling of PWR SNF rod conditions in U.S. inventory: A pseudo rod database was generated based on EIA database. Analysis of train accident scenarios: a probability density function (pdf) of a pinching load distribution P{sub 9m} of a 9-m cask drop is linearly scaled to generate a more representative: distribution P{sub AS}. The new distribution considers train collision velocities, cask drops for on- and off-bridge accidents, and cask impact orientation. LHS generates 3,000 rod pinching scenarios based on nine independent geometrical, material, and loading parameters. Probabilistic FE Analysis: The LHS realizations are used to generate 3,000 FE models. The FE analysis is conducted for P{sub 9m}, and P{sub AS}. The workflow is automatized. The OSC is computed with regression and interpolation models, as a function of the hydrogen content C{sub H}, vacuum drying peak σ{sub θ}, and SNF temperature T{sub FR}: The regression data was obtained in RCTs on PWR cladding. The OSD is derived from the cladding deformations in the FE models. Cladding failure was identified via a comparison of OSC and OSD. A 300 and a 30 years-of-dry-storage duration is investigated: The dry-storage period affects the expected SNF temperature and consequently, the OSC and the OSD. The material temperature LHS population is adapted to account for a SNF temperature variation. The probability of hydride-related, ductile cladding failure during a severe transportation accident is not significantly affected by the pinching load amplitude: The probability is slightly reduced from a range of less than 0.0333% to 0.1667% for P{sub 9m}, to a range of less than 0.0333% to 0.1333% for P{sub AS}, depending on the vacuum drying peak cladding hoop stress modelling assumptions. Further, the SNF dry-storage duration does not significantly affect the probability of ductile pinching failure. Based on available train accidents data, the pinching failure probability for a 5,000 km long transportation route ranges from less than 1.62 E-7 to a value of 6.47 E-7, depending on the vacuum drying peak cladding hoop stress modelling assumptions. The pinching load amplitude is not a controlling parameter for the probability of ductile cladding pinching failure, and does not control the cladding deformations once pellet-cladding contact is established, The pellet-cladding gap width, however, plays a key role. Fuel pellet degradation only affects the deformation capacity, if the pellet stiffness is significantly reduced (more than 90%). Only highly embrittled cladding is likely to fail. The deformation demands on the cladding under pinching loading are small.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2015

Sparing Healthy Tissue and Increasing Tumor Dose Using Bayesian Modeling of Geometric Uncertainties for Planning Target Volume Personalization

Objective: To develop a mathematical tool that can update a patient's planning target volume (PTV) partway through a course of radiation therapy to more precisely target the tumor for the remainder of treatment and reduce dose to surrounding healthy tissue. Methods and Materials: Daily on-board imaging was used to collect large datasets of displacements for patients undergoing external beam radiation therapy for solid tumors. Bayesian statistical modeling of these geometric uncertainties was used to optimally trade off between displacement data collected from previously treated patients and the progressively accumulating data from a patient currently partway through treatment, to optimally predict future displacements for that patient. These predictions were used to update the PTV position and margin width for the remainder of treatment, such that the clinical target volume (CTV) was more precisely targeted. Results: Software simulation of dose to CTV and normal tissue for 2 real prostate displacement datasets consisting of 146 and 290 patients treated with a minimum of 30 fractions each showed that re-evaluating the PTV position and margin width after 8 treatment fractions reduced healthy tissue dose by 19% and 17%, respectively, while maintaining CTV dose. Conclusion: Incorporating patient-specific displacement patterns from early in a course of treatment allows PTV adaptation for the remainder of treatment. This substantially reduces the dose to healthy tissues and thus can reduce radiation therapy–induced toxicities, improving patient outcomes.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
NASA NTRS2017

Geometrically Flexible and Efficient Flow Analysis of High Speed Vehicles Via Domain Decomposition, Part 1: Unstructured-Grid Solver for High Speed Flows

The ability to solve the equations governing the hypersonic turbulent flow of a real gas on unstructured grids using a spatially-elliptic, 2nd-order accurate, cell-centered, finite-volume method has been recently implemented in the VULCAN-CFD code. This paper describes the key numerical methods and techniques that were found to be required to robustly obtain accurate solutions to hypersonic flows on non-hex-dominant unstructured grids. The methods and techniques described include: an augmented stencil, weighted linear least squares, cell-average gradient method, a robust multidimensional cell-average gradient-limiter process that is consistent with the augmented stencil of the cell-average gradient method and a cell-face gradient method that contains a cell skewness sensitive damping term derived using hyperbolic diffusion based concepts. A data-parallel matrix-based symmetric Gauss-Seidel point-implicit scheme, used to solve the governing equations, is described and shown to be more robust and efficient than a matrix-free alternative. In addition, a y+ adaptive turbulent wall boundary condition methodology is presented. This boundary condition methodology is deigned to automatically switch between a solve-to-the-wall and a wall-matching-function boundary condition based on the local y+ of the 1st cell center off the wall. The aforementioned methods and techniques are then applied to a series of hypersonic and supersonic turbulent flat plate unit tests to examine the efficiency, robustness and convergence behavior of the implicit scheme and to determine the ability of the solve-to-the-wall and y+ adaptive turbulent wall boundary conditions to reproduce the turbulent law-of-the-wall. Finally, the thermally perfect, chemically frozen, Mach 7.8 turbulent flow of air through a scramjet flow-path is computed and compared with experimental data to demonstrate the robustness, accuracy and convergence behavior of the unstructured-grid solver for a realistic 3-D geometry on a non-hex-dominant grid.

White, Jeffery A.↗
DOE OSTI2016

GeMS/GSAOI photometric and astrometric performance in dense stellar fields

Ground-based imagers at 8 m class telescopes assisted by multi-conjugate adaptive optics are primary facilities with which to obtain accurate photometry and proper motions in dense stellar fields. We observed the central region of the globular clusters Liller 1 and NGC 6624 with the Gemini Multi-conjugate Adaptive Optics System (GeMS) feeding the Gemini South Adaptive Optics Imager (GSAOI) currently available at the Gemini South telescope, under different observing conditions. We characterized the stellar point-spread function (PSF) in terms of FWHM, Strehl ratio (SR), and encircled energy (EE), over the field of view (FOV). We found that, for sub-arcsecond seeing at the observed airmass, we can obtain the diffraction-limited PSF (FWHM ≈80 mas), SR ∼ 40%, and EE ≥ 50% with a dispersion around 10% over the FOV of 85″ × 85″, in the K {sub s} band. In the J band the best images provide FWHMs between 60 and 80 mas, SR >10%, and EE>40%. For seeing at the observed airmass exceeding 1″, the performance worsens but it is still possible to perform PSF fitting photometry with 25% EE in J and 40% in K {sub s}. We also computed the geometric distortions of GeMS/GSAOI and we obtained corrected images with an astrometric accuracy of ∼1 mas in a stellar field with high crowding.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

A Spatially Resolved Study of Cold Dust, Molecular Gas, H ii Regions, and Stars in the z = 2.12 Submillimeter Galaxy ALESS67.1

We present detailed studies of a z = 2.12 submillimeter galaxy, ALESS67.1, using sub-arcsecond resolution ALMA, adaptive optics-aided VLT/SINFONI, and Hubble Space Telescope (HST)/CANDELS data to investigate the kinematics and spatial distributions of dust emission (870 μm continuum), {sup 12}CO(J = 3–2), strong optical emission lines, and visible stars. Dynamical modeling of the optical emission lines suggests that ALESS67.1 is not a pure rotating disk but a merger, consistent with the apparent tidal features revealed in the HST imaging. Our sub-arcsecond resolution data set allows us to measure half-light radii for all the tracers, and we find a factor of 4–6 smaller sizes in dust continuum compared to all the other tracers, including {sup 12}CO; also, ultraviolet (UV) and Hα emission are significantly offset from the dust continuum. The spatial mismatch between the UV continuum and the cold dust and gas reservoir supports the explanation that geometrical effects are responsible for the offset of the dusty galaxy on the IRX–β diagram. Using a dynamical method we derive an α{sub CO}=1.8±1.0, consistent with other submillimeter galaxies (SMGs) that also have resolved CO and dust measurements. Assuming a single α{sub CO} value we also derive resolved gas and star formation rate surface densities, and find that the core region of the galaxy (≲5 kpc) follows the trend of mergers on the Schmidt–Kennicutt relationship, whereas the outskirts (≳5 kpc) lie on the locus of normal star-forming galaxies, suggesting different star formation efficiencies within one galaxy. Our results caution against using single size or morphology for different tracers of the star formation activity and gas content of galaxies, and therefore argue the need to use spatially resolved, multi-wavelength observations to interpret the properties of SMGs, and perhaps even for z>1 galaxies in general.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Three-dimensional local ALE-FEM method for fluid flow in domains containing moving boundaries/objects interfaces

A three-dimensional finite element method for the numerical simulations of fluid flow in domains containing moving rigid objects or boundaries is developed. The method falls into the general category of Arbitrary Lagrangian Eulerian methods; it is based on a fixed mesh that is locally adapted in the immediate vicinity of the moving interfaces and reverts to its original shape once the moving interfaces go past the elements. The moving interfaces are defined by separate sets of marker points so that the global mesh is independent of interface movement and the possibility of mesh entanglement is eliminated. The results is a fully robust formulation capable of calculating on domains of complex geometry with moving boundaries or devises that can also have a complex geometry without danger of the mesh becoming unsuitable due to its continuous deformation thus eliminating the need for repeated re-meshing and interpolation. Moreover, the boundary conditions on the interfaces are imposed exactly. This work is intended to support the internal combustion engines simulator KIVA developed at Los Alamos National Laboratories. The model's capabilities are illustrated through application to incompressible flows in different geometrical settings that show the robustness and flexibility of the technique to perform simulations involving moving boundaries in a three-dimensional domain.

16 TIDAL AND WAVE POWER↗
DOE OSTI2016

PHOTO-REVERBERATION MAPPING OF A PROTOPLANETARY ACCRETION DISK AROUND A T TAURI STAR

Theoretical models and spectroscopic observations of newborn stars suggest that protoplantary disks have an inner “wall” at a distance set by the disk interaction with the star. Around T Tauri stars, the size of this disk hole is expected to be on a 0.1 au scale that is unresolved by current adaptive optics imaging, though some model-dependent constraints have been obtained by near-infrared interferometry. Here we report the first measurement of the inner disk wall around a solar-mass young stellar object, YLW 16B in the ρ Ophiuchi star-forming region, by detecting the light-travel time of the variable radiation from the stellar surface to the disk. Consistent time lags were detected on two nights, when the time series in H (1.6 μ m) and K (2.2 μ m) bands were synchronized while the 4.5 μ m emission lagged by 74.5 ± 3.2 s. Considering the nearly edge-on geometry of the disk, the inner rim should be 0.084 au from the protostar on average, with an error of order 0.01 au. This size is likely larger than the range of magnetospheric truncations and consistent with an optically and geometrically thick disk front at the dust sublimation radius at ∼1500 K. The widths of the cross-correlation functions between the data in different wavebands place possible new constraints on the geometry of the disk.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Improvements to the APBS biomolecular solvation software suite: Improvements to the APBS Software Suite

The Adaptive Poisson-Boltzmann Solver (APBS) software was developed to solve the equations of continuum electrostatics for large biomolecular assemblages that has provided impact in the study of a broad range of chemical, biological, and biomedical applications. APBS addresses three key technology challenges for understanding solvation and electrostatics in biomedical applications: accurate and efficient models for biomolecular solvation and electrostatics, robust and scalable software for applying those theories to biomolecular systems, and mechanisms for sharing and analyzing biomolecular electrostatics data in the scientific community. To address new research applications and advancing computational capabilities, we have continually updated APBS and its suite of accompanying software since its release in 2001. In this manuscript, we discuss the models and capabilities that have recently been implemented within the APBS software package including: a Poisson-Boltzmann analytical and a semi-analytical solver, an optimized boundary element solver, a geometry-based geometric flow solvation model, a graph theory based algorithm for determining pKa values, and an improved web-based visualization tool for viewing electrostatics.

Biomedical applications↗