Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Local Fourier analysis”

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 109 records · Page 6

X-point effects on the ideal MHD modes in tokamaks in the description of dual-poloidal-region safety factor

The flux coordinates with dual-region safety factor (q) in the poloidal direction are developed in this work. The X-point effects on the ideal MHD modes in tokamaks are then analyzed using this coordinate system. Since the X-point effects mainly affect the edge region, the modes localized at the tokamak edge are particularly examined. Two types of modes are studied. The first is related to the conventional peeling or peeling-ballooning modes. The mode existence aligned with the local magnetic field in the poloidally core region as observed experimentally is confirmed. The X points are shown to contribute to a stabilizing effect for the conventionally treated modes with the surface-averaged q and with the tokamak edge portion truncated. The other is the axisymmetric modes localized in the vicinity of X points, which can affect the cross-field-line transport near the X points. The existence of axisymmetric modes points to the possibility of applying a toroidally axisymmetric resonant magnetic perturbation (RMP) in the X-point area for mitigating the edge localized modes, which can be an alternative to the current RMP design. The dual q description also has important implications for the existing non-axisymmetric RMP concept. It helps to understand why the RMP suppression of edge localized modes is difficult to achieve in the double-null tokamak configurations and points to the possibility of further improving the current RMP concept by considering the alignment to the local q.

Fourier analysis↗

Hybrid methods for rotordynamic analysis

Effective procedures are presented for the response analysis of the Space Shuttle Main Engine turbopumps under transient loading conditions. Of particular concern is the determination of the nonlinear response of the systems to rotor imbalance in presence of bearing clearances. The proposed procedures take advantage of the nonlinearities involved being localized at only a few rotor/housing coupling joints. The methods include those based on integral formulations for the incremental solutions involving the transition matrices of the rotor and housing. Alternatively, a convolutional representation of the housing displacements at the coupling points is proposed which would allow performing the transient analysis on a reduced model of the housing. The integral approach is applied to small dynamical models to demonstrate the efficiency of the approach. For purposes of assessing the numerical integration results for the nonlinear rotor/housing systems, a numerical harmonic balance procedure is developed to enable determining all possible harmonic, subharmonic, and nonperiodic solutions of the systems. A brief account of the Fourier approach is presented as applied to a two degree of freedon rotor-support system.

Noah, Sherif T.↗

Trace Gas Measurements on Mars and Earth Using Optical Parametric Generation

Trace gases and their isotopic ratios in planetary atmospheres offer important but subtle clues as to the origins of a planet's atmosphere, hydrology, geology, and potential for biology. An orbiting laser remote sensing instrument is capable of measuring trace gases on a global scale with unprecedented accuracy, and higher spatial resolution that can be obtained by passive instruments. We have developed an active sensing instrument for the remote measurement of trace gases in planetary atmospheres (including Earth). The technique uses widely tunable, seeded optical parametric generation (OPG) to measure methane, CO2, water vapor, and other trace gases in the near and mid-infrared spectral regions. Methane is a strong greenhouse gas on Earth and it is also a potential biogenic marker on Mars and other planets. Methane in the Earth's atmosphere survives for a shorter time than CO2 but its impact on climate change can be larger than CO2. Methane levels have remained relatively constant over the last decade around 1.78 parts per million (ppm) but recent observations indicate that methane levels may be on the rise. Increasing methane concentrations may trigger a positive feedback loop and a subsequent runaway greenhouse effect, where increasing temperatures result in increasing methane levels. The NRC Decadal Survey recognized the importance of global observations of greenhouse gases and called for simultaneous CH4, CO, and CO2 measurements but also underlined the technological limitations for these observations. For Mars, methane measurements are of great interest because of its potential as a strong biogenic marker. A remote sensing instrument that can measure day and night over all seasons and latitudes can identify and localize sources of biogenic gas plumes produced by subsurface chemistry or biology, and aid in the search for extra-terrestrial life. It can identify the dynamics of methane generation over time and latitude and identify future lander mission sites for more detailed in-situ analysis. In this paper we report on remote sensing measurements of methane using a high peak power, widely tunable optical parametric generator (OPG) operating at 3.3 micron and 1.65 micron. The OPG is pumped by a passively q-switched single frequency laser (3ns, 5KHz, 50uJ) and seeded by a diode laser. The spectral width of both signal and idler of seeded OPG is nearly Fourier transform limited. The output of seeded OPG is single frequency with high spectral purity and is widely tunable. Both 1650 nm and 3300 nm can be generated with a conversion efficiency of more than 30%. We have demonstrated detection of methane at 3274 nm and 1650 nm in a cell and also performed open path atmospheric measurements of methane at the same wavelengths. Finally, we were able to demonstrate simultaneous detection of methane at 3270.4 nm and CO2 at 1578.2 nm. In this paper we will discuss the OPG performance and atmospheric open path measurement results.

Numata, Kenji↗

Scale-Dependent Bedform Migration and Deformation in the Physical and Spectral Domains

Multi-scale bathymetries observed in laboratory channels and natural rivers have posed several challenges to the description of bedform geometry and kinematics. To better quantify scale-dependent migration velocity, high resolution spatio-temporal bed evolution data are analyzed using three independent methods: (i) a bedform tracking method that identifies individual bedforms in longitudinal bed elevation profiles, (ii) a cross-correlation based, image matching technique in consecutive bathymetry scans similar to Particle Image Velocimetry, and (iii) two-dimensional frequency-wavenumber spectra marking a dispersion relation between bedform period and length, allowing estimation of Fourier-based, scale-dependent migration velocities. The comparative studies show that the spectral and bedform tracking methods present small differences in the migration velocity of medium-large bedforms. However, significant deviations occur in the small bedforms, mostly due to bedform sheltering or amplified exposure to the flow, which also contribute to their enhanced scale-dependent deformation. The analysis also shows that some form of scale separation is needed to extract the large and slowly migrating bedforms as they are masked by secondary features when several bedform orders coexist. Bathymetric Image Velocimetry reasonably estimates the averaged migration velocity and exhibits potential to study relations between spanwise bed surface velocity and local bedform kinematics.

58 GEOSCIENCES↗

Multi-harmonic Rutherford island theory

Rutherford island theory, which governs the nonlinear evolution of tearing modes in tokamak plasmas, is generalized to take into account situations in which the conventional one-harmonic approximation is not valid. The analysis incorporates non-inductive currents driven by radio frequency (RF) electromagnetic waves injected into the plasma. A multi-harmonic tearing mode dispersion relation is derived that takes the form of a nonlinear inhomogeneous matrix eigenvalue problem. The dispersion relation is solved in the so-called two-harmonic approximation, in which only the principal Fourier harmonic of the perturbed magnetic flux and its first overtone are included in the calculation. In the absence of RF current drive, the nonlinear behavior of a tearing mode predicted in the two-harmonic approximation does not differ substantially from that predicted in the one-harmonic approximation. On the other hand, RF current drive that is sufficiently localized in the vicinity of the O-points of the mode's magnetic island chain is capable of triggering bifurcations of the O-points (which is impossible in the one-harmonic approximation). However, the current drive is incapable of triggering bifurcations of the island X-points. This finding is significant because Bardóczi and Evans recently observed bifurcations of magnetic island chain O-points in the presence of RF current drive in the DIII-D tokamak but did not observe bifurcations of the X-points. Finally, the changes in the topology of the magnetic island flux-surfaces induced by RF current drive are found to facilitate the stabilization of the tearing mode.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

The periodic variation of spokes in Saturn's rings

The discovery of a periodic variation in spoke activity in Saturn's rings from the analysis of Voyager images is reported. A Fourier power spectrum was computed using a data set generated by quantifying spoke activity observed on the morning (western) half of the rings in Voyager images spanning 5.5 days (about 12 Saturn rotations). A peak approximately 14 sigma above the noise in the frequency domain was found at a period of 621 + or - 22 min. Within the formal error, this value is consistent with the 639.4-min rotation period of Saturn's magnetic field. Maximum spoke activity is most likely to be observed on the morning half of the rings when a particular magnetic field sector coincides with this area. This magnetic sector contains the region which is aligned with local noon at the time of the emission of the SKR (Saturn kilometric radiation). These results suggest that the fundamental period of spoke variation is that of Saturn's magnetic field, and that spoke activity is associated with the region of the field which gives rise to the SKR. Passage of this region through Saturn's shadow may play a significant role in the creation and/or rejuvenation of spokes.

Porco, C. C.↗

A time Fourier analysis of zonal averaged ozone heating rates

A time-Fourier analysis is presented for the yearly variation of the zonal averaged ozone heating rates in the middle atmosphere based on a model study. The ozone heating rates are determined by utilizing two-dimensional ozone distributions, the altitude and latitude, and by including the effect of the curved earth's atmosphere. In addition, assumptions are introduced to the yearly variations of the ozone distributions due to the lack of sufficient existing ozone data. Among other results, it is shown that the first harmonic component indicates that the heating rates are completely out of phase between the northern and southern hemispheres. The second Fourier component shows a symmetric pattern with respect to the equator, as well as five distinct local extreme values of the ozone heating rate. The third harmonic component shows a pattern close to that of the first component except in the regions above 70 deg between 45-95 km in both hemispheres.

Wang, P.-H.↗

NASA Tech Briefs, May 2011

Topics covered include: 1) Method to Estimate the Dissolved Air Content in Hydraulic Fluid; 2) Method for Measuring Collimator-Pointing Sensitivity to Temperature Changes; 3) High-Temperature Thermometer Using Cr-Doped GdAlO3 Broadband Luminescence; 4)Metrology Arrangement for Measuring the Positions of Mirrors of a Submillimeter Telescope; 5) On-Wafer S-Parameter Measurements in the 325-508-GHz Band; 6) Reconfigurable Microwave Phase Delay Element for Frequency Reference and Phase-Shifter Applications; 7) High-Speed Isolation Board for Flight Hardware Testing; 8) High-Throughput, Adaptive FFT Architecture for FPGA-Based Spaceborne Data Processors; 9) 3D Orbit Visualization for Earth-Observing Missions; 10) MaROS: Web Visualization of Mars Orbiting and Landed Assets; 11) RAPID: Collaborative Commanding and Monitoring of Lunar Assets; 12) Image Segmentation, Registration, Compression, and Matching; 13) Image Calibration; 14) Rapid ISS Power Availability Simulator; 15) A Method of Strengthening Composite/Metal Joints; 16) Pre-Finishing of SiC for Optical Applications; 17) Optimization of Indium Bump Morphology for Improved Flip Chip Devices; 18) Measuring Moisture Levels in Graphite Epoxy Composite Sandwich Structures; 19) Marshall Convergent Spray Formulation Improvement for High Temperatures; 20) Real-Time Deposition Monitor for Ultrathin Conductive Films; 21) Optimized Li-Ion Electrolytes Containing Triphenyl Phosphate as a Flame-Retardant Additive; 22) Radiation-Resistant Hybrid Lotus Effect for Achieving Photoelectrocatalytic Self-Cleaning Anticontamination Coatings; 23) Improved, Low-Stress Economical Submerged Pipeline; 24) Optical Fiber Array Assemblies for Space Flight on the Lunar Reconnaissance Orbiter; 25) Local Leak Detection and Health Monitoring of Pressurized Tanks; 26) Dielectric Covered Planar Antennas at Submillimeter Wavelengths for Terahertz Imaging; 27) Automated Cryocooler Monitor and Control System; 28) Broadband Achromatic Phase Shifter for a Nulling Interferometer; 29) Super Dwarf Wheat for Growth in Confined Spaces; 30) Fine Guidance Sensing for Coronagraphic Observatories; 31) Single-Antenna Temperature- and Humidity-Sounding Microwave Receiver; 32) Multi-Wavelength, Multi-Beam, and Polarization-Sensitive Laser Transmitter for Surface Mapping; 33) Optical Communications Link to Airborne Transceiver; 34) Ascent Heating Thermal Analysis on Spacecraft Adaptor Fairings; 35) Entanglement in Self-Supervised Dynamics; 36) Prioritized LT Codes; 37) Fast Image Texture Classification Using Decision Trees; 38) Constraint Embedding Technique for Multibody System Dynamics; 39) Improved Systematic Pointing Error Model for the DSN Antennas; 40) Observability and Estimation of Distributed Space Systems via Local Information-Exchange Networks; 41) More-Accurate Model of Flows in Rocket Injectors; 42) In-Orbit Instrument-Pointing Calibration Using the Moon as a Target; 43) Reliability of Ceramic Column Grid Array Interconnect Packages Under Extreme Temperatures; 44) Six Degrees-of-Freedom Ascent Control for Small-Body Touch and Go; and 45) Optical-Path-Difference Linear Mechanism for the Panchromatic Fourier Transform Spectrometer.

Source record↗

Investigation of crossed-twin structure formation in magnesium and magnesium alloys

In this work, the effect of alloying addition on the propensity for twin-twin interactions to transform into crossed-twin structures in magnesium alloys is investigated. A full-field elasto-viscoplastic fast Fourier transform (EVP-FFT) framework combined with a discrete twin model and dislocation density-based hardening law for slip strengths is used to calculate the micromechanical fields in the crystalline matrix around the interacting twins. AZ31 and MgLi alloys are selected along with pure Mg to study the influence of plastic anisotropy in connection with alloying elements. These alloys were selected since their plastic anisotropy measure, which is defined as the ratio between the critical resolved shear stress for pyramidal $\langle$c +a$\rangle$ and basal $\langle$a$\rangle$ slip modes, spanned a wide range. To quantify the role of twin thicknesses, we probe a range of impinging twin thicknesses while fixing the recipient twin thickness. The analysis reveals that: (i) the local driving stress for crossed twin structure formation generated from the interaction of the two twins is lower in a low plastically anisotropic material, like a MgLi alloy, than a high plastically anisotropic material like pure Mg and (ii) the critical impinging twin thickness needed to form the crossed twin structure in pure Mg, AZ31 and MgLi alloys is ~0.5, ~0.75 and ~1.5 times the recipient twin thickness. We propose a relationship between the tendency for crossed twin structure formation and the experimentally observed higher ductility in MgLi alloys compared to pure Mg. Finally, one key implication of the findings is that crossed- twin structure formation can be hindered and the ductility of magnesium alloy thereby improved by properly choosing alloying elements that lower the slip strength for pyramidal $\langle$c + a$\rangle$ slip.

36 MATERIALS SCIENCE↗

Advanced infrared astronomy

The CO2 laser heterodyne spectrometer was used at the 3-m IRTF on Mauna Kea to make measurements of Mars during the 1984 opposition. Analysis of the observations of the mesospheric non-thermal emission demonstrated the existence of a warming of the Mars polar mesosphere, similar to the seasonal effects which are well known to occur at the Earth's mesopause. A search for CO2 and NH3 on Comet Halley was done with the new Kitt Peak system as well as with the IRTF heterodyne system. A Lamb-dip absorption cell was designed and constructed. Its use will allow extreme frequency stabilization of the laser local oscillator, which will greatly facilitate measurements of winds and dynamical phenomena. The Lamb-dip cell was used at Kitt Peak to study zonal and meridional winds in the atmosphere of Venus. Water vapor was detected in Comet Halley using Fourier transform spectrometer. The 2.65 micrometer upsilon sub 3 band was seen in emission, confirming non-thermal-equilibrium excitation models for comets. A study was made of the variability of Jovian ethane emission. The average volume mixing ratio of ethane in the Jovian stratosphere was found to be 3x10 to the -6 power, with the greatest variability seen in the auroral regions.

Mumma, M. J.↗

Experimental Assessment of Mammalian Lipidome Complexity Using Multimodal 21 T FTICR Mass Spectrometry Imaging

Herein, we assess the complementarity and complexity of data that can be detected within mammalian lipidome mass spectrometry imaging (MSI) via matrix-assisted laser desorption ionization (MALDI) and nanospray desorption electrospray ionization (nano-DESI). We do so by employing 21 T Fourier transform ion cyclotron resonance mass spectrometry (FTICR-MS) with absorption mode FT processing in both cases, allowing unmatched mass resolving power per unit time (≥613k at m/z 760, 1.536 s transients). While our results demonstrated that molecular coverage and dynamic range capabilities were greater in MALDI analysis, nano-DESI provided superior mass error, and all annotations for both modes had sub-ppm error. Taken together, these experiments highlight the coverage of 1676 lipids and serve as a functional guide for expected lipidome complexity within nano-DESI-MSI and MALDI-MSI. To further assess the lipidome complexity, mass splits (i.e., the difference in mass between neighboring peaks) within single pixels were collated across all pixels from each respective MSI experiment. The spatial localization of these mass splits was powerful in informing whether the observed mass splits were biological or artificial (e.g., matrix related). Mass splits down to 2.4 mDa were observed (i.e., sodium adduct ambiguity) in each experiment, and both modalities highlighted comparable degrees of lipidome complexity. Further, we highlight the persistence of certain mass splits (e.g., 8.9 mDa; double bond ambiguity) independent of ionization biases. Here, we also evaluate the need for ultrahigh mass resolving power for mass splits ≤4.6 mDa (potassium adduct ambiguity) at m/z > 1000, which may only be resolved by advanced FTICR-MS instrumentation.

21T-FTICR-MS↗

Spectrometer technology recommendations

A typical heterodyne remote sensing system contains three major elements: the antenna, the radiometer, and the spectrometer. The radiometer consists of the local oscillator, the mixer, and the intermediate frequency amplifiers. This subsystem performs the function of down converting the high frequency incident thermal emission signal to a lower intermediate frequency. The spectrometer measures the power spectrum of the down-converted signal simultaneously in many contiguous frequency channels. Typical spectrum analysis requirements involve measurement of signal bandwidths of 100 to 1000 MHz with a channel resolution of 0.5 to 10 MHz. Three general approaches are used for spectrometers: (1) filter banks, (2) Acousto-Optic Spectrometers (AOS's), and (3) digital autocorrelators. In contrast to the two frequency domain techniques, an autocorrelator works in the time domain. The autocorrelation function (ACF) of the incoming signal is computed and averaged over the integration time. The averaged ACF is then Fourier transformed to obtain the signal power spectrum. Significant progress was made in the development of sub mm antennas and radiometers. It is now time to begin research in the development of low power spaceborne spectrometers and to reduce their size and weight. The near-term research goal will be to develop a prototype digital autocorrelation spectrometer, using VLSI gate array technology, which will have a small size, low power requirements, and can be used in spacecraft mm and sub mm radiometer systems. The long-range objective of this technology development is to make extremely low power, less than 10 mW/channel, small and stable wideband spectrometers which can be used in future mm and sub mm wavelength space missions such as the Large Deployable Reflector.

Wilson, William J.↗

Quantitative approaches for multiscale structural analysis with atomic resolution electron microscopy

Atomic-resolution imaging with scanning transmission electron microscopy is a powerful tool for characterizing the nanoscale structure of materials, in particular features such as defects, local strains, and symmetry-breaking distortions. In addition to advanced instrumentation, the effectiveness of the technique depends on computational image analysis to extract meaningful features from complex datasets recorded in experiments, which can be complicated by the presence of noise and artifacts, small or overlapping features, and the need to scale analysis over large representative areas. Here, we present image analysis approaches which synergize real and reciprocal space information to efficiently and reliably obtain meaningful structural information with picometer scale precision across hundreds of nanometers of material from atomic-resolution electron microscope images. Damping superstructure peaks in reciprocal space allows symmetry-breaking structural distortions to be disentangled from other sources of inhomogeneity and measured with high precision. Real-space fitting of the wavelike signals resulting from Fourier filtering enables absolute quantification of lattice parameter variations and strain, as well as the uncertainty associated with these measurements. Implementations of these algorithms are made available as an open source python package.

36 MATERIALS SCIENCE↗

White-dwarf variability and the rotation of g-modes

The multiperiodic behavior clearly evident after Fourier analysis of many DA white-dwarf light curves is interpreted as arising from outbursts whose timing is regulated by the rotation and oscillation of the star. A model is proposed on the basis of three main assumptions: (1) a broad array of g-mode oscillations is taking place all the time in DA white-dwarf variables; (2) slow rotation permits most of the g-modes to couple and form a small group of nonlinear modes, each characterized by one of the indices of the spherical harmonic functions involved; and (3) the observed brightenings of the star are produced by upwelling of convection on a large horizontal scale when the antinodes of the nonlinear oscillation patterns periodically rotate into alignment and cause a temporary local enhancement of energy per unit volume. This model is used to match precisely almost every strong periodicity observed in the complex light curves of the DA white dwarfs G207-9, G38-29, G29-38, and HL Tau 76. If the model has been applied correctly, the results indicate that all four of these variables are rotating slowly with periods in the range from 250 to 500 sec.

Wolff, C. L.↗

[Cosmic Microwave Background (CMB) Anisotropies]

One of the main areas of research is the theory of cosmic microwave background (CMB) anisotropies and analysis of CMB data. Using the four year COBE data we were able to improve existing constraints on global shear and vorticity. We found that, in the flat case (which allows for greatest anisotropy), (omega/H)0 less than 10(exp -7), where omega is the vorticity and H is the Hubble constant. This is two orders of magnitude lower than the tightest, previous constraint. We have defined a new set of statistics which quantify the amount of non-Gaussianity in small field cosmic microwave background maps. By looking at the distribution of power around rings in Fourier space, and at the correlations between adjacent rings, one can identify non-Gaussian features which are masked by large scale Gaussian fluctuations. This may be particularly useful for identifying unresolved localized sources and line-like discontinuities. Levin and collaborators devised a method to determine the global geometry of the universe through observations of patterns in the hot and cold spots of the CMB. We have derived properties of the peaks (maxima) of the CMB anisotropies expected in flat and open CDM models. We represent results for angular resolutions ranging from 5 arcmin to 20 arcmin (antenna FWHM), scales that are relevant for the MAP and COBRA/SAMBA space missions and the ground-based interferometer. Results related to galaxy formation and evolution are also discussed.

Silk, Joseph↗

Analytical models for the enhancement of fusion reactivity by turbulence

The reactivity of fusion plasma depends not only on its local density and temperature but also, through a recently identified kinetic effect, on the relative velocities of nearby fluid elements. Turbulence on fine spatial scales, therefore, enhances fusion reactivity. The enhancement is quantified here for general subsonic turbulent flows. Leveraging this effect in the design of inertial confinement fusion experiments could enable substantial energy savings.

Fluid dynamics↗

Time-resolved biphase signatures of quadratic nonlinearity observed in coupled Alfvén eigenmodes on the DIII-D tokamak

We report the detection of nonstationary quadratic coupling between toroidicity-induced Alfvén eigenmodes (TAEs) on sub-millisecond time scales. Identification of phase coherency between multiple TAEs and nonlinearly generated modes is facilitated by wavelet-based bicoherence analysis of time-series from inductive coils, taken from a DIII-D discharge heated by neutral beam injection (NBI). Characterization of nonlinear three-wave interaction is inferred by stationary local bispectrum phase (biphase) and confirmed via bandpass filtering. Biphase dynamics associated with prominent bispectral features are well-resolved in time and consistent with transient quadratic coupling. Onset and duration of nonlinearity are correlated with enhanced amplitude of participating TAEs; coincident changes in amplitude are observed for modes at difference frequency |f TAE,1 - f TAE,2 |.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Parabolized stability equations

The parabolized stability equations (PSE) are a new approach to analyze the streamwise evolution of single or interacting Fourier modes in weakly nonparallel flows such as boundary layers. The concept rests on the decomposition of every mode into a slowly varying amplitude function and a wave function with slowly varying wave number. The neglect of the small second derivatives of the slowly varying functions with respect to the streamwise variable leads to an initial boundary-value problem that can be solved by numerical marching procedures. The PSE approach is valid in convectively unstable flows. The equations for a single mode are closely related to those of the traditional eigenvalue problems for linear stability analysis. However, the PSE approach does not exploit the homogeneity of the problem and, therefore, can be utilized to analyze forced modes and the nonlinear growth and interaction of an initial disturbance field. In contrast to the traditional patching of local solutions, the PSE provide the spatial evolution of modes with proper account for their history. The PSE approach allows studies of secondary instabilities without the constraints of the Floquet analysis and reproduces the established experimental, theoretical, and computational benchmark results on transition up to the breakdown stage. The method matches or exceeds the demonstrated capabilities of current spatial Navier-Stokes solvers at a small fraction of their computational cost. Recent applications include studies on localized or distributed receptivity and prediction of transition in model environments for realistic engineering problems. This report describes the basis, intricacies, and some applications of the PSE methodology.

Herbert, Thorwald↗