Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “shape approximation”

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 631 records · Page 35
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Manipulating ultrafast even-order nonlinear chiral responses of L-tryptophan by polarization pulse shaping

Molecular chirality has long been monitored in the frequency domain in the ultraviolet, visible, and infrared regimes. Recently developed time-domain approaches can detect time-dependent chiral dynamics by enhancing intrinsically weak chiral signals. Even-order nonlinear signals in chiral molecules have gained attention thanks to their existence in the electric dipole approximation, without relying on the weaker higher-order multipole interactions. We illustrate the optimization of temporal polarization pulse-shaping in various frequency ranges (infrared/optical and optical/X ray) to enhance chiral nonlinear signals. These signals can be recast as an overlap integral of matter and field pseudoscalars which contain the relevant chiral information. Simulations are carried out for second- and fourth-order nonlinear spectroscopies in L-tryptophan.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Learning Optimal Aerodynamic Designs

This project created a framework for efficient, accurate, and scalable deep neural network representations of design optimization problem solutions. The inputs to these DNN representations are the vector of design requirement parameters, the outputs are the optimal design variables, and the goal is to learn the map from inputs to outputs (i.e., inverse design). The team addressed the problem of the optimal shape design of aerodynamic lifting surfaces—in particular aircraft wings—using a Reynolds-Average Navier Stokes model to govern the CFD-based aerodynamic shape optimization. The inverse design map for such problems is very complex and high-dimensional, involving inputs and outputs on the order of 1000s. To approximate this inverse design map, the team developed algorithms to construct parsimonious DNN architectures, which automatically identify low-dimensional manifolds in which design requirements affect optimal shape parameters, and trained these architectures with multifidelity optimization methods. The resulting methodology accurately and automatically designs optimal aerodynamic lifting surfaces with very high accuracy (99%) at interactive speeds, of the order of milliseconds, resulting in factors of one million or more speedup relative to CFD-based design optimization.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1993

Temperature and abundances in the Jovian auroral stratosphere. 1: Ethane as a probe of the millibar region

We report infrared heterodyne spectroscopy (lambda/delta lambda is approximately 10(exp 6)) of C2H6 emission at 11.9 microns from the northern Jovian auroral region, in observations conducted over December 2-7, 1989. Accurately measured line shapes provide information on C2H6 abundance as well as temperature and permit retrieval of the source pressure region. Enhanced emission was observed in the longitude range approximately 150-180 deg at approximately 60 deg north latitude, approximately corresponding to the CH4 7.8 micron hot spot and the region of brightest UV aurora. Significant brightness variations were observed in the hot spot emissions on a time scale of approximately 20 hours. Analysis of the brightest hot spot spectra indicates C2H6 mole fractions of approximately (6.3-6.8) x 10(exp -6) at temperatures of approximately 182-184 K at 1 mbar, compared to mole fractions of (3.8 +/- 1.4) x 10(exp -6) averaged over spectra outside the hot spot at a temperature of approximately 172 K at the same pressure. Fixing the mole fraction to the lower limit retrieved in the quiescent (non-hot spot) region allows the temperature at 1 mbar to be as high as approximately 200 K within the hot spot. These results provide upper limits to the temperature increase near the source of the C2H6 thermal infrared emission. Combined with results from similar measurements of ethylene emission probing the approximately 10-microbar region (Kostiuk et al., this issue), altitude information on the thermal structure of the Jovian auroral stratosphere can be obtained for the first time.

Livengood, Timothy A.↗
NASA NTRSDate not supplied

Additional Large-Drop Ice Accretion Test Results for a Large Scale Swept Wing Section from January 2022

In-flight icing is an important consideration that affects aircraft design, performance, certification and safety. Newer regulations combined with increasing demand to reduce fuel burn, emissions and noise are driving a need for improvements in icing simulation capability. To that end, this paper presents the results of additional ice accretion testing conducted in the NASA Icing Research Tunnel in January 2022 with a large swept wing section typical of a modern commercial transport. The model was based upon a section of the Common Research Model wing at the 64% semispan station with a streamwise chord length of 136 in. The test conditions were developed with an icing scaling analysis to generate similar conditions for a small median volumetric diameter (MVD) = 25 μm cloud and a large MVD = 110 μm cloud. A series of tests were conducted over a range of total temperature from -23.8 ˚C to -1.4 ˚C with all other conditions held constant. Another series of tests explored cloud MVD variations from 50 μm to 230 μm while holding constant certain scaling parameters. The variation in ice mass and scanned ice volume across repeat conditions was approximately 50% lower than the uncertainty in the cloud MVD and liquid water content. The measured ice mass and volume calculated from the 3D scans were used to compute the mass to volume ratio that is sometimes referred to as ice density or void fraction. When the ice volume based on the ice shape maximum combined cross section was used to determine this ratio, the resulting values were in the range of 240 to 455 kg/m 3 . This is consistent with analogous values previously reported in the literature. The ice shape mass and volume increased with MVD from 50 to 230 μm at fixed values of the scaling parameters. The ice mass to volume ratio was approximately constant for all of the cases which showed that the accreted mass and volume increased at approximately the same rate. These results demonstrate the significance of cloud MVD on ice shape mass and volume.

Ice accretion, Large drops, Wind tunnel↗
NASA NTRSDate not supplied

Additional Large-Drop Ice Accretion Test Results for a Large Scale Swept Wing Section from January 2022

In-flight icing is an important consideration that affects aircraft design, performance, certification and safety. Newer regulations combined with increasing demand to reduce fuel burn, emissions and noise are driving a need for improvements in icing simulation capability. To that end, this paper presents the results of additional ice accretion testing conducted in the NASA Icing Research Tunnel in January 2022 with a large swept wing section typical of a modern commercial transport. The model was based upon a section of the Common Research Model wing at the 64% semispan station with a streamwise chord length of 136 in. The test conditions were developed with an icing scaling analysis to generate similar conditions for a small median volumetric diameter (MVD) = 25 μm cloud and a large MVD = 110 μm cloud. A series of tests were conducted over a range of total temperature from -23.8 ˚C to -1.4 ˚C with all other conditions held constant. Another series of tests explored cloud MVD variations from 50 μm to 230 μm while holding constant certain scaling parameters. The variation in ice mass and scanned ice volume across repeat conditions was approximately 50% lower than the uncertainty in the cloud MVD and liquid water content. The measured ice mass and volume calculated from the 3D scans were used to compute the mass to volume ratio that is sometimes referred to as ice density or void fraction. When the ice volume based on the ice shape maximum combined cross section was used to determine this ratio, the resulting values were in the range of 240 to 455 kg/m3. This is consistent with analogous values previously reported in the literature. The ice shape mass and volume increased with MVD from 50 to 230 μm at fixed values of the scaling parameters. The ice mass to volume ratio was approximately constant for all of the cases which showed that the accreted mass and volume increased at approximately the same rate. These results demonstrate the significance of cloud MVD on ice shape mass and volume.

Ice accretion↗
DOE OSTI2022

Absolute vs Convective Instabilities and Front Propagation in Lipid Membrane Tubes

Here we analyze the stability of biological membrane tubes, with and without a base flow of lipids. Membrane dynamics are completely specified by two dimensionless numbers: the well-known Föppl-von Kármán number Γ and the recently introduced Scriven-Love number SL, respectively quantifying the base tension and base flow speed. For unstable tubes, the growth rate of a local perturbation depends only on Γ, whereas SL governs the absolute versus convective nature of the instability. Furthermore, nonlinear simulations of unstable tubes reveal an initially localized disturbance result in propagating fronts, which leave a thin atrophied tube in their wake. Depending on the value of Γ, the thin tube is connected to the unperturbed regions via oscillatory or monotonic shape transitions-reminiscent of recent experimental observations on the retraction and atrophy of axons. We elucidate our findings through a weakly nonlinear analysis, which shows membrane dynamics may be approximated by a model of the class of extended Fisher-Kolmogorov equations. Our study sheds light on the pattern selection mechanism in axonal shapes by recognizing the existence of two Lifshitz points, at which the front dynamics undergo steady-to-oscillatory bifurcations.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Spacecraft Geometry Effects on Kinetic Impactor Missions

Abstract The DART (Double Asteroid Redirection Test) mission will impact a spacecraft on the secondary (Dimorphos) of the binary asteroid system Didymos in 2022 September, with the goal of altering the orbital period of Dimorphos about Didymos sufficiently to be observed from ground-based observations. Numerical impact modeling is a crucial component in understanding the outcome of the DART experiment, and while many have investigated the effects of target properties, such as material strength and porosity (which remain unknown), an often overlooked factor is the importance of accurately representing the spacecraft itself in such models. Most impact modeling to date has considered simple impactor geometries such as a solid uniform sphere, but in reality the spacecraft is a complex shape full of different components, open spaces, and thin walled structures. At a minimum, a simple solid representation underestimates the surface area of the impact: for a small body such as Dimorphos (approximately 160 m in diameter), the difference between a spacecraft spanning 20 m (including solar arrays) impacting and a sub-1 m idealized shape may be important. In this paper, we compare models impacting high-fidelity models of the spacecraft based on the CAD geometry with various simplified impactors, in order to assess the potential importance of this effect. We find that the difference between the simplest impactor geometries (such as a uniform sphere) and the real spacecraft is measurable, and has an interesting dependence on the material properties of the asteroid itself.

36 MATERIALS SCIENCE↗
NASA NTRS2019

On the Convergence of Numerical Computations for Both Exact and Approximate Solutions for Electromagnetic Scattering by Nonspherical Dielectric Particles

We summarize the size parameter range of the applicability of four lightscattering computational methods for nonspheric dielectric particles. These methods include two exact methods - the extended boundary condition method (EBCM) and the invariant imbedding T-matrix method (II-TM) and two approximate approaches - the physical-geometric optics method (PGOM) and the improved geometric optics method (IGOM). For spheroids, the single-scattering properties computed by EBCM and II-TM agree for size parameters up to 150, and the comparison gives us confidence in using IITM as a benchmark for size parameters up to 150 for other geometries (e.g., hexagonal columns) because the applicability of II-TM with respect to particle shape is generic, as demonstrated in our previous studies involving a complex aggregate. This study demonstrates the convergence of the exact II-TM and approximate PGOM solutions for the complete set of single-scattering properties of a nonspherical shape other than spheroids and circular cylinders with particle sizes of ~48λ, specifically a hexagonal column with a size parameter of length as kL=300, where k=2π/λ and L is the column length. IGOM is also quite accurate except near the exact 180º backscattering direction. This study demonstrates that a synergetic combination of the numerically-exact II-TM and the approximate PGOM can seamlessly cover the entire size parameter range of practical interest. To demonstrate the applicability of the approach, we compute the optical properties of dust particles with a downstream application to the retrieval of dust aerosol optical thickness and effective particle size from polarimetric observations.

Ping Yang↗
DOE OSTI2024

Dynamic Shaping of Grid Response of Multi-Machine Multi-Inverter Systems Through Grid-Forming IBRs

We consider the problem of controlling the frequency response of weakly-coupled multi-machine multi-inverter low-inertia power systems via grid-forming inverter-based resources (IBRs). In contrast to existing methods, our approach relies on dividing the larger system into multiple strongly-coupled subsystems, without ignoring either the underlying network or approximating the subsystem response as an aggregate harmonic mean model. Rather, through a structured clustering and recursive dynamic shaping approach, the frequency response of the overall system to load perturbations is shaped appropriately. We demonstrate the proposed approach for a three-node triangular configuration and a small-scale radial network. Furthermore, previous synchronization analysis for heterogeneous systems requires the machines to satisfy certain proportionality property. In our approach, the effective transfer functions for each cluster can be tuned by the IBRs to satisfy such property, enabling us to apply the shaping control to systems with a wider range of heterogeneous machines.

frequency-shaping control↗
DOE OSTI2024

Dynamic Shaping of Grid Response of Multi-Machine Multi-Inverter Systems Through Grid-Forming IBRs: Preprint

We consider the problem of controlling the frequency response of weakly-coupled multi-machine multi-inverter low-inertia power systems via grid-forming inverter-based resources (IBRs). In contrast to existing methods, our approach relies on dividing the larger system into multiple strongly-coupled subsystems, without ignoring either the underlying network or approximating the subsystem response as an aggregate harmonic mean model. Rather, through a structured clustering and recursive dynamic shaping approach, the frequency response of the overall system to load perturbations is shaped appropriately. We demonstrate the proposed approach for a three-node triangular configuration and a small-scale radial network. Furthermore, for small-scale radial microgrids, we demonstrate the ability of IBRs to tune the effective transfer functions of synchronous machines. This enables us to relax the uniform turbine time-constant assumptions and widen the scope of existing synchronization results for proportionally heterogeneous machines.

frequency-shaping control↗
NASA NTRS1982

Effect of enclosure shape on natural convection velocities

A numerical analysis was performed to compare natural convection velocities in two dimensional enclosures of various shape. The following shapes were investigated: circle, square, horizontal and upright 2 x 1 aspect ratio rectangles, horizontal and upright half circles, diamond. In all cases, the length scale in the various dimensionless parameters, such as Rayleigh number, is defined as the diameter of the equal area circle. Natural convection velocities were calculated for Rayleigh numbers of 1000 and 5000 with the temperature difference taken to be across (1) the maximum horizontal dimension, (2) the median horizontal line (line through centroid) and (3) the horizontal distance such that the temperature gradient is the same for shapes of equal area. For the class of shapes including the square, upright half circle and upright rectangle, the computed velocities were found to agree very closely with that of the equal area circle when the temperature difference is taken to be across the maximum horizontal dimension (condition (a)). The velocities for the horizontal rectangle and half circle were found to be approximately one half that of the equal area circle for the same condition. Better overall agreement among all shapes was obtained by setting the temperature difference across a distance such that the temperature gradients were equal for shapes of equal area.

Robertson, S. J.↗
NASA NTRS1989

Measured and predicted root-mean-square errors in square and triangular antenna mesh facets

Deflection shapes of square and equilateral triangular facets of two tricot-knit, gold plated molybdenum wire mesh antenna materials were measured and compared, on the basis of root mean square (rms) differences, with deflection shapes predicted by linear membrane theory, for several cases of biaxial mesh tension. The two mesh materials contained approximately 10 and 16 holes per linear inch, measured diagonally with respect to the course and wale directions. The deflection measurement system employed a non-contact eddy current proximity probe and an electromagnetic distance sensing probe in conjunction with a precision optical level. Despite experimental uncertainties, rms differences between measured and predicted deflection shapes suggest the following conclusions: that replacing flat antenna facets with facets conforming to parabolically curved structural members yields smaller rms surface error; that potential accuracy gains are greater for equilateral triangular facets than for square facets; and that linear membrane theory can be a useful tool in the design of tricot knit wire mesh antennas.

Fichter, W. B.↗
DOE OSTI2020

An arbitrarily high-order three-dimensional Cartesian-grid method for reconstructing interfaces from volume fraction fields

Here Tthis work describes a newly developed, arbitrarily high-order Cartesian-grid method for reconstructing material interfaces from a volume fraction field. The method begins by identifying all of the grid cells in the volume fraction field that are intersected by the interface and need to be approximated by the reconstruction scheme. Finite-differences are used to calculate the gradient of the volume fraction field and provide an estimate of the surface normal in all of the interfacial grid cells. Groups of connected grid cells are then identified which all have the same dominant component of the normal vector. This grouping by orientation determines the proper dependent variable to use in the surface reconstruction (e.g. for a 2D curve, this step determines if the surface will be approximated by a function of x or y). A cumulative integral over the surface is constructed and fit using b-splines for two-dimensional problems or tensor-product b-splines for three-dimensional problems. This construction allows for the interface to be recovered through application of the second fundamental theorem of calculus. Fitting the cumulative integral with $\mathscr{N}$ th-order b-splines (or tensor-product b-splines) yields an ($\mathscr{N}$-1) th-order convergence rate of the interface shape. Differentiation of the b-spline interface function(s) allows for the high-order approximation of the normal vector and curvature to be obtained directly anywhere along b-spline. Together, the proposed reconstruction technique can achieve arbitrarily high mesh convergence rates. Validation tests are presented with mesh convergence rates ranging from fourth- to tenth-order.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1983

The detection of X rays from Jupiter

X rays in the energy band 0.2-3.0 keV have been detected coming from both polar regions of Jupiter. The observations were made in 1979 and 1981 by using the imaging proportional counter and high resolution imaging detectors on the Einstein X-ray astronomy satellite. The measured flux density of approximately 0.0006/sq cm-sec at earth corresponds to an X ray luminosity of approximately 4 x 10 to the 9th W in the 0.2- to 3.0-keV energy band. The energy spectrum of the X rays is extremely soft and can be characterized by a power law with an exponent of approximately 2.3. Detector energy resolution is insufficient to distinguish a soft line spectrum from a continuum. However, the shape of the response and the observed X ray power indicate that the source of this auroral emission is not electron bremsstrahlung as on the earth, but is most probably line emission from O and S ions with energies between 0.03 and 4.0 MeV/nucleon precipitating from the outer boundary of the Io plasma torus at L approximately 8.

Metzger, A. E.↗
NASA NTRS2007

Iapetus' Geophysics: Rotation Rate, Shape, and Equatorial Ridge

Iapetus has preserved evidence that constrains the modeling of its geophysical history from the time of its accretion until now. The evidence is (a) its present 79.33-day rotation or spin rate, (b) its shape that corresponds to the equilibrium figure for a hydrostatic body rotating with a period of approximately 16 h, and (c) its high, equatorial ridge, which is unique in the Solar System. This paper reports the results of an investigation into the coupling between Iapetus' thermal and orbital evolution for a wide range of conditions including the spatial distributions with time of composition, porosity, short-lived radioactive isotopes (SLRI), and temperature. The thermal model uses conductive heat transfer with temperature-dependent conductivity. Only models with a thick lithosphere and an interior viscosity in the range of about the water ice melting point can explain the observed shape. Short-lived radioactive isotopes provide the heat needed to decrease porosity in Iapetus? early history. This increases thermal conductivity and allows the development of the strong lithosphere that is required to preserve the 16-h rotational shape and the high vertical relief of the topography. Long-lived radioactive isotopes and SLRI raise internal temperatures high enough that significant tidal dissipation can start, and despin Iapetus to synchronous rotation. This occurred several hundred million years after Iapetus formed. The models also constrain the time when Iapetus formed because the successful models are critically dependent upon having just the right amount of heat added by SLRI decay in this early period. The amount of heat available from short-lived radioactivity is not a free parameter but is fixed by the time when Iapetus accreted, by the canonical concentration of Al-26, and, to a lesser extent, by the concentration of Fe-60. The needed amount of heat is available only if Iapetus accreted between 2.5 and 5.0Myr after the formation of the calcium aluminum inclusions as found in meteorites. Models with these features allow us to explain Iapetus? present synchronous rotation, its fossil 16-h shape, and the context within which the equatorial ridge arose.

geophysics↗
NASA NTRS1981

The albedo of snow for partially cloudy skies

The albedo of snow is defined as the ratio of reflected to incident solar energy, and it is an important parameter in the earth's radiation budget analysis and in the study of snowpack's thermal conditions. An approximate model for calculating the incident spectral flux for partially cloudy skies is presented. The input parameters for the calculation are atmospheric precipitable water, turbidity, ozone content, surface pressure, the optical thickness of clouds, and the grain size of snow crystals. The spectral snow reflectance model considers both specular surface reflection and volumetric multiple scattering. The surface reflection is calculated by using a crystal-shape-dependent bidirectional reflectance distribution function; the volumetric multiple scattering is calculated by using a crystal-size-dependent approximate solution in the radiative transfer equation. The model yields spectral and integrated solar flux and snow reflectance as a function of solar elevation and fractional cloud-cover. The illustrative insolation and albedo values were obtained from spectral reflectance and incident flux for representative parameters of Antarctic coastal regions. A simple relationship between grain size and the overcast albedo was obtained. For a set of grain size and shape, the albedo as a function of solar elevation and fractional cloud cover was tabulated.

Choudhury, B. J.↗
DOE OSTI2022

A neural network‐enhanced reproducing kernel particle method for modeling strain localization

Abstract Modeling the localized intensive deformation in a damaged solid requires highly refined discretization for accurate prediction, which significantly increases the computational cost. Although adaptive model refinement can be employed for enhanced effectiveness, it is cumbersome for the traditional mesh‐based methods to perform while modeling the evolving localizations. In this work, neural network‐enhanced reproducing kernel particle method (NN‐RKPM) is proposed, where the location, orientation, and shape of the solution transition near a localization is automatically captured by the NN approximation via a block‐level neural network (NN) optimization. The weights and biases in the blocked parameterization network control the location and orientation of the localization. The designed basic four‐kernel NN block is capable of capturing a triple junction or a quadruple junction topological pattern, while more complicated localization topological patters are captured by the superposition of multiple four‐kernel NN blocks. The standard RK approximation is then utilized to approximate the smooth part of the solution, which permits a much coarser discretization than the high‐resolution discretization needed to capture sharp solution transitions with the conventional methods. A regularization of the NN approximation is additionally introduced for discretization‐independent material responses. The effectiveness of the proposed NN‐RKPM is verified by a series of numerical verifications.

Baek, Jonghyuk↗
DOE OSTI2023

A hybrid 3D/2D field response calculation for liquid argon detectors with PCB based anode plane

Liquid Argon Time Projection Chamber (LArTPC) technology is commonly utilized in neutrino detector designs. It enables detailed reconstruction of neutrino events with high spatial precision and low energy threshold. Its field response (FR) model describes the time-dependent electric currents induced in the anode-plane electrodes when ionization electrons drift nearby. An accurate and precise FR is a crucial input to LArTPC detector simulations and charge reconstruction. Established LArTPC designs have been based on parallel wire planes. It allows accurate and computationally economic two-dimensional (2D) FR models utilizing the translational symmetry along the direction of the wires. Recently, novel LArTPC designs utilize electrodes formed on printed circuit board (PCB) in the shape of strips with through holes. The translational symmetry is no longer a good approximation near the electrodes and a new FR calculation that employs regions with three dimensions (3D) has been developed. Extending the 2D models to 3D would be computationally expensive. Fortuitously, the nature of strips with through holes allows for a computationally economic approach based on the finite-difference method (FDM). In this paper, we present a new software package pochoir that calculates LArTPC field response for these new strip-based anode designs. This package combines 3D calculations in the volume near the electrodes with 2D far-field solutions to achieve fast and precise field response computation. We apply the resulting FR to simulate and reconstruct samples of cosmic-ray muons and 39 Ar decays from a Vertical Drift (VD) detector prototype operated at CERN. We find the difference between real and simulated data within 5%. Current state-of-the-art LArTPC software requires a 2D FR which we provide by averaging over one dimension and estimate that variations lost in this average are smaller than 7%.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗