Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “test vector”

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 1,459 records · Page 81

High Voltage EEE Parts for EMA/EHA Applications on Manned Launch Vehicles

The objective of this paper is an assessment of high voltage electronic components required for high horsepower electric thrust vector control (TVC) systems for human spaceflight launch critical application. The scope consists of creating of a database of available Grade 1 electrical, electronic and electromechanical (EEE) parts suited to this application, a qualification path for potential non-Grade 1 EEE parts that could be used in these designs, and pathfinder testing to validate aspects of the proposed qualification plan. Advances in the state of the art in high power electric power systems enable high horsepower electric actuators, such as the electromechnical actuator (EMA) and the electro-hydrostatic actuator (EHA), to be used in launch vehicle TVC systems, dramaticly reducing weight, complexity and operating costs. Designs typically use high voltage insulated gate bipolar transistors (HV-IGBT). However, no Grade 1 HV-IGBT exists and it is unlikely that market factors alone will produce such high quality parts. Furthermore, the perception of risk, the lack of qualification methodoloy, the absence of manned space flight heritage and other barriers impede the adoption of commercial grade parts onto the critical path. The method of approach is to identify high voltage electronic component types and key parameters for parts currently used in high horsepower EMA/EHA applications, to search for higher quality substitutes and custom manufacturers, to create a database for these parts, and then to explore ways to qualify these parts for use in human spaceflight launch critical application, including grossly derating and possibly treating hybrid parts as modules. This effort is ongoing, but results thus far include identification of over 60 HV-IGBT from four manufacturers, including some with a high reliability process flow. Voltage ranges for HV-IGBT have been identified, as has screening tests used to characterize HV-IGBT. BSI BS ISO 21350 Space systems Off-the-shelf item utilization, developed from Marshall Work Instruction MWI8060.1 OFF-THE-SHELF HARDWARE UTILIZATION IN FLIGHT HARDWARE DEVELOPMENTwas found to provide guidance for including commercial off-the-shelf (COTS) hardware for use in critical applications.

Griffin, Trent↗

Simplified Applied Loads of a Transonic Truss-Braced Wing Derived from Computational Fluid Dynamics

The transonic truss-braced wing (TTBW) concept is under consideration by the NASA Advanced Air Transport Technology (AATT) project due to its potential to reduce fuel consumption of transport aircraft. To properly evaluate the TTBW concept, detailed structural modeling and subsequent test planning is underway, which requires a simplified set of discrete applied loads that simulate the distributed aerodynamic pressures. Given the unique architecture of the TTBW and its intended operation in the transonic regime, computational fluid dynamics (CFD) is the ideal approach to employ to ensure the aerodynamics of this configuration are sufficiently captured. To convert aerodynamic results into input for structural models, application of a data transfer tool, Discrete Data Transfer Between Dissimilar Meshes (DDTBDM), is presented to resolve aerodynamic pressures from CFD into single load vectors at hardpoints along both the wing and the strut. Evaluation of mapping the aerodynamic forces to hardpoints at rib and spar intersections in a structural model is performed, and the results are also compared with a linear static aeroelastic analysis. The structural responses indicate that distributing forces along only the wing and strut forward spars is insufficient for simulating the distributed aerodynamic loads. Adding moments along the forward spars improved the global response of the wing, but the point moments created unrealistic local rotations within the structure. When distributing loads along both the forward and aft spars of the wing and strut, the global response was less sensitive to the inclusion of point moments, which indicated the arrangement of point forces alone was sufficient input for structural design work and testing and provided a more continuous loading profile.

Finite Element↗

Simplified Applied Loads of a Transonic Truss-Braced Wing Derived from Computational Fluid Dynamic

The transonic truss-braced wing (TTBW) concept is under consideration by the NASA Advanced Air Transport Technology (AATT) project due to its potential to reduce fuel consumption of transport aircraft. To properly evaluate the TTBW concept, detailed structural modeling and subsequent test planning is underway, which requires a simplified set of discrete applied loads that simulate the distributed aerodynamic pressures. Given the unique architecture of the TTBW and its intended operation in the transonic regime, computational fluid dynamics (CFD) is the ideal approach to employ to ensure the aerodynamics of this configuration are sufficiently captured. To convert aerodynamic results into input for structural models, application of a data transfer tool, Discrete Data Transfer Between Dissimilar Meshes (DDTBDM), is presented to resolve aerodynamic pressures from CFD into single load vectors at hardpoints along both the wing and the strut. Evaluation of mapping the aerodynamic forces to hardpoints at rib and spar intersections in a structural model is performed, and the results are also compared with a linear static aeroelastic analysis. The structural responses indicate that distributing forces along only the wing and strut forward spars is insufficient for simulating the distributed aerodynamic loads. Adding moments along the forward spars improved the global response of the wing, but the point moments created unrealistic local rotations within the structure. When distributing loads along both the forward and aft spars of the wing and strut, the global response was less sensitive to the inclusion of point moments, which indicated the arrangement of point forces alone was sufficient input for structural design work and testing and provided a more continuous loading profile.

Finite Element↗

Orion Optical Navigation Progress Toward Exploration: Mission 1

Optical navigation of human spacecraft was proposed on Gemini and implemented successfully on Apollo as a means of autonomously operating the vehicle in the event of lost communication with controllers on Earth. It shares a history with the "method of lunar distances" that was used in the 18th century and gained some notoriety after its use by Captain James Cook during his 1768 Pacific voyage of the HMS Endeavor. The Orion emergency return system utilizing optical navigation has matured in design over the last several years, and is currently undergoing the final implementation and test phase in preparation for Exploration Mission 1 (EM-1) in 2019. The software development is being worked as a Government Furnished Equipment (GFE) project delivered as an application within the Core Flight Software of the Orion camera controller module. The mathematical formulation behind the initial ellipse fit in the image processing is detailed in Christian. The non-linear least squares refinement then follows the technique of Mortari as an estimation process of the planetary limb using the sigmoid function. The Orion optical navigation system uses a body fixed camera, a decision that was driven by mass and mechanism constraints. The general concept of operations involves a 2-hour pass once every 24 hours, with passes specifically placed before all maneuvers to supply accurate navigation information to guidance and targeting. The pass lengths are limited by thermal constraints on the vehicle since the OpNav attitude generally deviates from the thermally stable tail-to-sun attitude maintained during the rest of the orbit coast phase. Calibration is scheduled prior to every pass due to the unknown nature of thermal effects on the lens distortion and the mounting platform deformations between the camera and star trackers. The calibration technique is described in detail by Christian, et al. and simultaneously estimates the Brown-Conrady coefficients and the Star Tracker/Camera interlock angles. Accurate attitude information is provided by the star trackers during each pass. Figure 1 shows the various phases of lunar return navigation when the vehicle is in autonomous operation with lost ground communication. The midcourse maneuvers are placed to control the entry interface conditions to the desired corridor for safe landing. The general form of optical navigation on Orion is where still images of the Moon or Earth are processed to find the apparent angular diameter and centroid in the camera focal plane. This raw data is transformed into range and bearing angle measurements using planetary data and precise star tracker inertial attitude. The measurements are then sent to the main flight computer's Kalman filter to update the onboard state vector. The images are, of course, collected over an arc to converge the state and estimate velocity. The same basic technique was used by Apollo to satisfy loss-of-comm, but Apollo used manual crew sightings with a vehicle-integral sextant instead of autonomously processing optical imagery. The software development is past its Critical Design Review, and is progressing through test and certification for human rating. In support of this, a hardware-in-the-loop test rig was developed in the Johnson Space Center Electro-Optics Lab to exercise the OpNav system prior to integrated testing on the Orion vehicle. Figure 2 shows the rig, which the test team has dubbed OCILOT (Orion Camera In the Loop Optical Testbed). Analysis performed to date shows a delivery that satisfies an allowable entry corridor as shown in Figure 3.

Holt, Greg N.↗

Machine Learning–Based Condition Monitoring of a Circulating Water System of a Canadian Nuclear Plant

With the need to maintain long-term reliable energy using nuclear power plants, there is an underlying demand to ensure that the maintenance of plant components and systems is also done in an efficient and cost-effective manner. One way to achieve this is by moving from time-based maintenance to condition-based maintenance. The research presented in this paper focuses on applying statistical and machine-learning-based methods to capture anomalies within data for fault detection to further develop into condition monitoring. This paper focuses on system data for a circulating water system (CWS) of a pressurized heavy-water reactor for detecting anomalies. The different methodologies used for detecting and capturing anomalies in the CWS data are matrix profile, density-based spatial clustering of applications with noise (DBSCAN), and support vector machines (SVMs). Matrix profile and DBSCAN are used to distinguish between normal data and anomalous data. This paper presents a hybrid method using DBSCAN and SVM when a portion of the data is used for DBSCAN to generate clusters. This portion of data is then used to train the SVM along with the clusters generated by DBSCAN as output. SVM is then tested on unseen data as a predictive tool, which can work in real time to categorize data points as either normal or anomalous. This paper presents results that show the high accuracies of DBSCAN and SVM in capturing anomalies within the data for a CWS for fault detection. Thus, the maintenance plan would be focused on component condition rather than a time-based schedule by switching to an automated system to identify and predict faults within a CWS.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Crustal Deformation and the Seismic Cycle across the Kodiak Islands, Alaska

The Kodiak Islands are located approximately 130 to 250 km from the Alaska-Aleutian Trench where the Pacific plate is underthrusting the North American plate at a rate of about 57 mm/yr. The southern extent of the 1964 Prince William Sound (${M-w}$ = 9.2) earthquake rupture occurred offshore and beneath the eastern portion of the Kodiak Islands. Here we report GPS results (1993-2001) from northern Kodiak Island that span the transition between the 1964 uplift region along the eastern coast and the region of coseismic subsidence further inland. The horizontal velocity vectors range from 22.9 $\pm$ 2.2 mm/yr at N26.3$\deg$W $\pm$ 2.5$\deg$, about 150 km from the trench, to 5.9 $\pm$ 1.3 mm/yr at N65.9$\deg$W $\pm$ 6.6$\deg$, about 190 km from the trench. Near the northeastern coast of Kodiak the velocity vector above the shallow, locked main thrust zone is between the orientation of PCFC-NOAM plate motion (N22$/deg$W) and the trench-normal (N3O$\deg$W). Further west, our geodetic results suggest the accumulation of shear strain that will be released eventually as left-lateral motion on upper plate faults such as the Kodiak Island fault. These results are consistent with the hypothesis that the difference between the Pacific-North American plate motion and the orientation of the down going slab would lead to 4-8 mm/yr of left-lateral slip. Short-term geodetic uplift rates range from 2 - 14 mm/yr, with the maximum uplift located near the axis of maximum subsidence during the 1964 earthquake. We evaluated alternate interseismic models for Kodiak to test the importance of various mechanisms responsible for crustal deformation rates. These models are based on the plate interface slip history inferred from earlier modeling of coseismic and post-seismic geodetic results. The horizontal (trench perpendicular) and vertical deformation rates across Kodiak are consistent with a model that includes the viscoelastic response to : (1) a downgoing Pacific plate interface that is locked at shallow depths,(2) coseismic slip in the 1964 and (3) interseismic creep below the seismogenic zone. The change in orientation of the horizontal velocity vector occurs down-dip from the locked main thrust zone. In southern Kodiak, the coseismic slip in the 1964 earthquake was smaller than in the northern Kodiak region; yet, the horizontal, interseismic velocities as a function of distance from the trench are comparable to those in northern Kodiak. Based on the earthquake history prior to, and following the 1964 earthquake, we hypothesize that the plate interface in southern Kodiak slips in more frequent large earthquakes than in northern Kodiak.

Sauber, Jeanne↗

It’s all right(-handed neutrinos): a new W ' model for the R D ( &#x2217;<!-- * --> ) anomaly

The measured B-meson semi-leptonic branching ratios R D and R D &#x2217;<!-- * --> have long-standing deviations between theory and experiment. We introduce a model which explains both anomalies through a single interaction by introducing a right-handed neutrino as the missing energy particle. This interaction is mediated by a heavy charged vector boson (W ') which couples only to right-handed quarks and leptons of the Standard Model through the mixing of these particles with new vector-like fermions. Previous W ' models for the R D ( &#x2217;<!-- * --> ) anomaly were strongly constrained from flavor changing neutral currents and direct collider searches for Z' → τ τ resonances. We show that relying on right-handed fermion mixing enables us to avoid these constraints, as well as other severe bounds from electroweak precision tests and neutrino mixing.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

A Local Macroscopic Conservative (LoMaC) Low Rank Tensor Method for the Vlasov Dynamics

Abstract In this paper, we propose a novel Local Macroscopic Conservative (LoMaC) low rank tensor method for simulating the Vlasov-Poisson (VP) system. The LoMaC property refers to the exact local conservation of macroscopic mass, momentum and energy at the discrete level. This is a follow-up work of our previous development of a conservative low rank tensor approach for Vlasov dynamics ( arXiv:2201.10397 ). In that work, we applied a low rank tensor method with a conservative singular value decomposition to the high dimensional VP system to mitigate the curse of dimensionality, while maintaining the local conservation of mass and momentum. However, energy conservation is not guaranteed, which is a critical property to avoid unphysical plasma self-heating or cooling. The new ingredient in the LoMaC low rank tensor algorithm is that we simultaneously evolve the macroscopic conservation laws of mass, momentum and energy using a flux-difference form with kinetic flux vector splitting; then the LoMaC property is realized by projecting the low rank kinetic solution onto a subspace that shares the same macroscopic observables by a conservative orthogonal projection. The algorithm is extended to the high dimensional problems by hierarchical Tuck decomposition of solution tensors and a corresponding conservative projection algorithm. Extensive numerical tests on the VP system are showcased for the algorithm’s efficacy.

Guo, Wei↗

Testing thermal-relic dark matter with a dark-photon mediator

In light of recent DAMIC-M results, we present the status of thermal-relic dark matter χ coupled to a kinetically mixed dark photon A ′ . In the predictive “direct annihilation” regime, m A ′ > m χ , the relic abundance depends on the kinetic mixing parameter, and there is a minimum value compatible with thermal freeze-out. Using only electron- and nuclear-recoil direct-detection results, we find that, for complex scalar dark matter, the direct annihilation regime is now excluded for nearly all values of m χ ; the only exception is the resonant annihilation regime where m A ′ ≈ 2 m χ . Direct annihilation relic targets for other representative models, including Majorana and pseudo-Dirac candidates, remain viable across a wide range of model parameters but will be tested with a combination of dedicated accelerator searches in the near future. In the opposite “secluded annihilation” regime, where m χ > m A ′ , this scenario is excluded by cosmic microwave background measurements for all m χ ≲ 30 GeV . Similar conclusions in both the direct and secluded regimes hold for all anomaly-free vector mediators that couple to the first generation of electrically charged Standard Model particles.

Krnjaic, Gordan [Fermilab; Chicago U., KICP] (ORCI↗

Sci-Fri PM: Radiation Therapy, Planning, Imaging, and Special Techniques - 04: Assessment of intra-fraction motion during lung SABR VMAT using a custom abdominal compression device

Purpose: Lung SABR patients are treated using Volumetrically Modulated Arc Therapy (VMAT), utilizing 2 arcs with Conebeam CT (CBCT) image-guidance prior to each arc. Intra-fraction imaging can prolong treatment time (up to 20%), and the aim of this study is to determine if it is necessary. Methods: We utilize an in-house abdominal compression device to minimize respiratory motion, 4DCT to define the ITV, a 5 mm PTV margin and a 2–3 mm PRV margin. We treated 23 patients with VMAT, fifteen were treated to 48 Gy in 4 fractions, while eight were treated with up to 60 Gy in 8 fractions. Intrafraction motion was assessed by the translational errors recorded for the second CBCT. Results: There was no significant difference (t-test, p=0.93) in the intra-fraction motion between the patients treated with 4 and 8 fractions, or between the absolute translations in each direction (ANOVA, p=0.17). All 124 intra-fraction CBCT images were analysed and 95% remained localized within the 5 mm PTV margin The mean magnitude of the vector displacement was 1.8 mm. Conclusions: For patients localized with an abdominal compression device, the intrafraction CBCT image may not be necessary, if it is only the tumor coverage that is of concern, as the patients are typically well within the 5 mm PTV margin. On the other hand, if there is a structure with a smaller PRV margin, an intrafraction CBCT is recommended to ensure that the dose limit for the organ at risk is not exceeded.

60 APPLIED LIFE SCIENCES↗

Compressed Sensing with Sparse Corruptions: Fault-Tolerant Sparse Collocation Approximations

The recovery of approximately sparse or compressible coefficients in a polynomial chaos expansion is a common goal in many modern parametric uncertainty quantification (UQ) problems. However, relatively little effort in UQ has been directed toward theoretical and computational strategies for addressing the sparse corruptions problem, where a small number of measurements are highly corrupted. Such a situation has become pertinent today since modern computational frameworks are sufficiently complex with many interdependent components that may introduce hardware and software failures, some of which can be difficult to detect and result in a highly polluted simulation result. In this paper we present a novel compressive sampling--based theoretical analysis for a regularized $\ell^1$ minimization algorithm that aims to recover sparse expansion coefficients in the presence of measurement corruptions. Our recovery results are uniform (the theoretical guarantees hold for all compressible signals and compressible corruptions vectors) and prescribe algorithmic regularization parameters in terms of a user-defined a priori estimate on the ratio of measurements that are believed to be corrupted. We also propose an iteratively reweighted optimization algorithm that automatically refines the value of the regularization parameter and empirically produces superior results. Lastly, our numerical results test our framework on several medium to high dimensional examples of solutions to parameterized differential equations and demonstrate the effectiveness of our approach.

97 MATHEMATICS AND COMPUTING↗

Rapid Spaceborne Mapping of Wildfire Retardant Drops for Active Wildfire Management

Aerial application of fire retardant is a critical tool for managing wildland fire spread. Retardant applications are carefully planned to maximize fire line effectiveness, improve firefighter safety, protect high-value resources and assets, and limit environmental impact. However, topography, wind, visibility, and aircraft orientation can lead to differences between planned drop locations and the actual placement of the retardant. Information on the precise placement and areal extent of the dropped retardant can provide wildland fire managers with key information to (1) adaptively manage event resources, (2) assess the effectiveness of retardant slowing or stopping fire spread, (3) document location in relation to ecologically sensitive areas; and perform or validate cost-accounting for drop services. This study uses Sentinel-2 satellite data and commonly used machine learning classifiers to test an automated approach for detecting and mapping retardant application. We show that a multiclass model (retardant, burned, unburned, and cloud artifact classes) outperforms a single-class retardant model and that image differencing (post-application minus pre-application) outperforms single-image models. Compared to the random forest and support vector machine, the gradient boosting model performed the best with an overall accuracy of 0.88 and an F1 Score of 0.76 for fire retardant, though results were comparable for all three models. Our approach maps the full areal extent of the dropped retardant within minutes of image availability, rather than linear representations currently mapped by aerial GPS surveys. The development of this capability allows for the rapid assessment of retardant effectiveness and documentation of placement in relation to sensitive environments.

54 ENVIRONMENTAL SCIENCES↗

Nonaxisymmetric nozzle technology program - An overview

Nonaxisymmetric nozzles, with in-flight thrust vectoring and/or reversing capabilities, provide the potential for substantial gains in fighter aircraft performance, maneuvering, survivability and ground handling. This paper presents an overview of some of the current programs to develop nonaxisymmetric nozzle technology for future aircraft applications. The initial phase of this program, which consists of a series of subscale model tests, engine/nozzle studies, system integration studies and simulation studies, has been initiated. Successful completion of this program will allow designers of future aircraft to consider, with high confidence, the integration of nonaxisymmetric nozzle concepts to improve aircraft mission performance effectiveness, survivability and cost.

Berrier, B. L.↗

Multivariable control altitude demonstration on the F100 turbofan engine

The control system designed under the Multivariable Control Synthesis (MVCS) program for the F100 turbofan engine is described. The MVCS program, applied the linear quadratic regulator (LQR) synthesis methods in the design of a multivariable engine control system to obtain enhanced performance from cross-coupled controls, maximum use of engine variable geometry, and a systematic design procedure that can be applied efficiently to new engine systems. Basic components of the control system, a reference value generator for deriving a desired equilibrium state and an approximate control vector, a transition model to produce compatible reference point trajectories during gross transients, gain schedules for producing feedback terms appropriate to the flight condition, and integral switching logic to produce acceptable steady-state performance without engine operating limit exceedance are described and the details of the F100 implementation presented. The engine altitude test phase of the MVCS program, and engine responses in a variety of test operating points and power transitions are presented.

Lehtinen, B.↗

Data compression experiments with LANDSAT thematic mapper and Nimbus-7 coastal zone color scanner data

A case study is presented where an image segmentation based compression technique is applied to LANDSAT Thematic Mapper (TM) and Nimbus-7 Coastal Zone Color Scanner (CZCS) data. The compression technique, called Spatially Constrained Clustering (SCC), can be regarded as an adaptive vector quantization approach. The SCC can be applied to either single or multiple spectral bands of image data. The segmented image resulting from SCC is encoded in small rectangular blocks, with the codebook varying from block to block. Lossless compression potential (LDP) of sample TM and CZCS images are evaluated. For the TM test image, the LCP is 2.79. For the CZCS test image the LCP is 1.89, even though when only a cloud-free section of the image is considered the LCP increases to 3.48. Examples of compressed images are shown at several compression ratios ranging from 4 to 15. In the case of TM data, the compressed data are classified using the Bayes' classifier. The results show an improvement in the similarity between the classification results and ground truth when compressed data are used, thus showing that compression is, in fact, a useful first step in the analysis.

Tilton, James C.↗

Long-range sound propagation: A review of some experimental data

Three experimental studies of long range sound propagation carried out or sponsored in the past by NASA are briefly reviewed to provide a partial prospective for some of the analytical studies presented in this symposium. The three studies reviewed cover (1) a unique test of two large rocket engines conducted in such a way as to provide an indication of possible atmospheric scattering loss from a large low-frequency directive sound source, (2) a year-long measurement of low frequency sound propagation which clearly demonstrated the dominant influence of the vertical gradient in the vector sound velocity towards the receiver in defining excess sound attenuation due to refraction, and (3), a series of excess ground attenuation measurements over grass and asphalt surfaces replicated several times under very similar inversion weather conditions.

Sutherland, Louis C.↗

Numerical experiments in geomagnetic modeling

Numerical tests were made, using least squares fitting of a spherical harmonic model, to a selection of Magsat data to determine the practical limits of this technique with modern computers. The resulting (M102189) model, whose coefficients were adjusted up to n = 50, was compared with M07AV6, a previous model which used least squares (on vector data) for coefficients up to n = 29, and Gauss-Legendre quadrature (on Z residuals) to adjust the coefficients up to n = 63. For the new least squares adjustment to n = 50 a condition number of 115 was obtained for the solution matrix, with a resulting precision of 11 significant figures. The M102189 model shows a lower and more Gaussian residual distribution than did M07AV6, though the Gaussian envelope fits to the residual distributions, even for the scalar field, gives "standard deviations' never lower than 6 nT, a factor of three higher than the estimated Magsat observational errors. Ionospheric currents are noted to have a significant effect on the coefficients of the internal potential functions.

Cain, Joseph C.↗

Calibration for Thrust and Airflow Measurements in the CE-22 Advanced Nozzle Test Facility

CE-22 facility procedures and measurements for thrust and airflow calibration obtained with choked-flow ASME nozzles are presented. Six calibration nozzles are used at an inlet total pressure from 20 to 48 psia. Throat areas are from 9.9986 to 39.986 sq. in.. Throat Reynolds number varies from 1.8 to 7.9 million. Nozzle gross thrust coefficient (CFG) uncertainty is 0.25 to 0.75 percent, with smaller uncertainly generally for larger nozzles and higher inlet total pressure. Nozzle discharge coefficient (CDN) uncertainty is 0.15 percent or less for all the data. ASME nozzle calibrations need to be done before and after research model testing to achieve these uncertainties. In addition, facility capability in terms of nozzle pressure ratio (NPR) and nozzle airflow are determined. Nozzle pressure ratio of 50 or more is obtainable at 40 psia for throat areas between 20 and 30 sq. in.. Also presented are results for two of the ASME nozzles vectored at 10deg, a dead-weight check of the vertical (perpendicular to the jet axis) force measurement, a calibration of load cell forces for the effects of facility tank deflection with tank pressure, and the calibration of the metric-break labyrinth seal.

Werner, Roger A.↗