Search NASA⌕ Search

SEARCH · Search NASA

Results for “shape modeling”

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,081 records · Page 60

Trajectory Design Employing Convex Optimization for Landing on Irregularly Shaped Asteroids

Mission proposals that land spacecraft on asteroids are becoming increasingly popular. However, in order to have a successful mission the spacecraft must reliably and softly land at the intended landing site with pinpoint precision. The problem under investigation is how to design a propellant optimal powered descent trajectory that can be quickly computed onboard the spacecraft, without interaction from the ground control. The propellant optimal control problem in this work is to determine the optimal finite thrust vector to land the spacecraft at a specified location, in the presence of a highly nonlinear gravity field, subject to various mission and operational constraints. The proposed solution uses convex optimization, a gravity model with higher fidelity than Newtonian, and an iterative solution process for a fixed final time problem. In addition, a second optimization method is wrapped around the convex optimization problem to determine the optimal flight time that yields the lowest propellant usage over all flight times. Gravity models designed for irregularly shaped asteroids are investigated. Success of the algorithm is demonstrated by designing powered descent trajectories for the elongated binary asteroid Castalia.

Pinson, Robin M.↗

Integral boundary-layer models for turbulent separated flows

Two-dimensional, integral boundary-layer methods of the entrainment and dissipation type, with closure relationships which account for turbulent separated flows through the velocity and turbulence models, are developed. Kinematic shape factor relations used in the two methods are evaluated by comparison with recently obtained experimental data for separated flows. Prescribed-displacement-thickness computations are made with both methods and compared to experiment for a series of separated flows, including low speed and transonic shock-induced separation. Effects of upstream history on the turbulent shear stress development are shown to be important for the accurate prediction of the separated flows considered and the experimental trends are reasonably modeled by a rate equation derived from the turbulent kinetic energy equation.

Thomas, J. L.↗

Parallel Infrastructure Modeling and Inversion Module for E4D

Electrical resistivity tomography ERT is a method of imaging the electrical conductivity of the subsurface. Electrical conductivity is a useful metric for understanding the subsurface because it is governed by geomechanical and geochemical properties that drive subsurface systems. ERT works by injecting current into the subsurface across a pair of electrodes, and measuring the corresponding electrical potential response across another pair of electrodes. Many such measurements are strategically taken across an array of electrodes to produce an ERT data set. These data are then processed through a computationally demanding process known as inversion to produce an image of the subsurface conductivity structure that gave rise to the measurements. Data can be inverted to provide 2D images, 3D images, or in the case of time-lapse 3D imaging, 4D images. ERT is generally not well suited for environments with buried electrically conductive infrastructure such as pipes, tanks, or well casings, because these features tend to dominate and degrade ERT images. This reduces or eliminates the utility of ERT imaging where it would otherwise be highly useful for, for example, imaging fluid migration from leaking pipes, imaging soil contamination beneath leaking subsurface tanks, and monitoring contaminant migration in locations with dense network of metal cased monitoring wells. The location and dimension of buried metallic infrastructure is often known. If so, then the effects of the infrastructure can be explicitly modeled within the ERT imaging algorithm, and thereby removed from the corresponding ERT image. However, there are a number of obstacles limiting this application. 1) Metallic infrastructure cannot be accurately modeled with standard codes because of the large contrast in conductivity between the metal and host material. 2) Modeling infrastructure in true dimension requires the computational mesh to be highly refined near the metal inclusions, which increases computational demands. 3) The ERT imaging algorithm requires specialized modifications to accommodate high conductivity inclusions within the computational mesh. The solution to each of these challenges was implemented within E4D (formerly FERM3D), which is a parallel ERT imaging code developed at PNNL (IPID #30249). The infrastructure modeling module implement in E4D uses a method of decoupling the model at the metallic interface(s) boundaries, into several well posed sub-problems (one for each distinct metallic inclusion) that are subsequently solved and recombined to form the global solution. The approach is based on the immersed interface method, with has been applied for similar problems in other fields (e.g. semiconductor industry). Comparisons to analytic solutions have shown the results to be very accurate, addressing item 1 above. The solution is implemented about an unstructured mesh, which enables arbitrary shapes to be efficiently modelled, thereby addressing item 2 above. In addition, the algorithm is written in parallel and shows excellent scalability, which also addresses equation 2 above. Finally, because only the boundaries of metallic inclusions are modeled, there are no high conductivity cells within the modeling mesh, and the problem described by item 3 above is no longer applicable.

Johnson, TimC.↗

Accuracy problems associated with semi-analytical derivatives of static response

The semianalytical method is widely used for calculating derivatives of static response with respect to design variables for structures modeled by finite elements. This paper shows that the method can have serious accuracy problems for shape design variables in structures modeled by beam elements. The problem was first discovered in the analysis of a complex car model which is described. Next, it is shown that the same accuracy problem occurs for the simplest beam structures. Finally, the problem is shown to be associated with a high ratio of rigid body rotation to elastic deformation, and a test is developed to diagnose the severity of the problem in a given structure. Several examples are presented to show the validity of this test.

Barthelemy, Bruno↗

Net Shape Fabricated Low Cost MHK Pass-Through the Hub Turbine Blades with Integrated Health Management Technology

The primary objective of this project is to develop a three-blade MHK rotor with low manufacturing and maintenance costs. The proposed program will design, fabricate and test a novel half-scale low cost, net shape fabricated single piece three-blade MHK rotor with integrated health management technology to demonstrate significant Capital Expenditures (CAPEX) and Operational Expenditures (OPEX) cost reductions due to the novel design and manufacturing process. The proposed project is divided into three major tasks: Task 1: Single Piece Three-blade Kinetic Hydropower System (KHPS) Rotor Full-Scale and Half-Scale Design; Task 2: Composite Manufacturing Trials and Half-Scale Prototype Rotor Fabrication; and Task 3: Material Characterization and Half-Scale Prototype Test and Evaluation. These three tasks include design and analysis of full-scale and half-scale three-blade rotor prototypes using computational fluid dynamics (CFD) and finite-element analysis (FEA), demonstration of a novel half-scale net shape fabrication process, determination of a fatigue threshold composite strain allowable, three-blade rotor mold design, manufacture of half-scale rotor clam shell mold, three-blade rotor test fixture design and fabrication, development of final manufacturing and test plans, manufacture of the half-scale net shape composite single blade and three-blade prototypes, and test and evaluation of the half-scale rotor.

16 TIDAL AND WAVE POWER↗

Grid-Free 2D Plasma Simulations of the Complex Interaction Between the Solar Wind and Small, Near-Earth Asteroids

We present results from a new grid-free 2D plasma simulation code applied to a small, unmagnetized body immersed in the streaming solar wind plasma. The body was purposely modeled as an irregular shape in order to examine photoemission and solar wind plasma flow in high detail on the dayside, night-side, terminator and surface-depressed 'pocket' regions. Our objective is to examine the overall morphology of the various plasma interaction regions that form around a small body like a small near-Earth asteroid (NEA). We find that the object obstructs the solar wind flow and creates a trailing wake region downstream, which involves the interplay between surface charging and ambipolar plasma expansion. Photoemission is modeled as a steady outflow of electrons from illuminated portions of the surface, and under direct illumination the surface forms a non-monotonic or ''double-sheath'' electric potential upstream of the body, which is important for understanding trajectories and equilibria of lofted dust grains in the presence of a complex asteroid geometry. The largest electric fields are found at the terminators, where ambipolar plasma expansion in the body-sized night-side wake merges seamlessly with the thin photoelectric sheath on the dayside. The pocket regions are found to be especially complex, with nearby sunlit regions of positive potential electrically connected to unlit negative potentials and forming adjacent natural electric dipoles. For objects near the surface, we find electrical dissipation times (through collection of local environmental solar wind currents) that vary over at least 5 orders of magnitude: from 39 Micro(s) inside the near-surface photoelectron cloud under direct sunlight to less than 1 s inside the particle-depleted night-side wake and shadowed pocket regions

Asteroids↗

Morphology modification of gold nanoparticles from nanoshell to C-shape: Improved surface enhanced Raman scattering

Morphology modification of nanostructures is of great interest, because it can be used to fabricate nanostructures which are hard to be done using other methods. Different from traditional lithographic technique which is slow and expensive, morphology modification is easy, cheap, and reproducible. In this paper, modification of the optical and morphological properties of a hollow gold nanoshell (HGNS) is achieved by using H{sub 2}O{sub 2} as an oxidizer. The reshaping of these nanostructures has been demonstrated as a consequence of an oxidation process in which HGNSs are dissolved by H{sub 2}O{sub 2} under the acidic conditions provided by HCl. We investigate the oxidation process by a transmission electron microscope and propose a reshaping model involving four different shapes (HGNS, HGNS with hole, gold nanoring, and C-shaped gold nanoparticle) which are corresponding to the oxidation products of HGNSs at different pH values. Besides, the surface enhanced Raman scattering (SERS) activity of each oxidation product has been evaluated by using rhodamine 6G as the Raman active probe. It has been observed that the C-shaped gold nanoparticles which are corresponding to the oxidation products at the minimum pH value have the highest SERS activity and this result can also be interpreted by discrete-dipole approximation simulations. We demonstrate that the morphology modification of HGNSs becomes possible in a controlled manner using wet chemistry and can be used in preparation of gold nanoparticles such as HGNS with hole, gold nanoring, and C-shaped gold nanoparticle with large SERS activity. These nanostructures must have potential use in many plasmonic areas, including sensing, catalysis, and biomedicine.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Tensor Minkowski Functionals: first application to the CMB

Tensor Minkowski Functionals (TMFs) are tensor generalizations of the usual Minkowski Functionals which are scalar quantities. We introduce them here for use in cosmological analysis, in particular to analyze the Cosmic Microwave Background (CMB) radiation. They encapsulate information about the shapes of structures and the orientation of distributions of structures. We focus on one of the TMFs, namely W {sub 2}{sup 1,1}, which is the (1,1) rank tensor generalization of the genus. The ratio of the eigenvalues of the average of W {sub 2}{sup 1,1} over all structures, α, encodes the net orientation of the structures; and the average of the ratios of the eigenvalues of W {sub 2}{sup 1,1} for each structure, β, encodes the net intrinsic anisotropy of the structures. We have developed a code that computes W {sub 2}{sup 1,1}, and from it α and β, for a set of structures on the 2-dimensional Euclidean plane. We use it to compute α and β as functions of chosen threshold levels for simulated Gaussian and isotropic CMB temperature and E mode fields. We obtain the value of α to be one for both temperature and E mode, which means that we recover the statistical isotropy of density fluctuations that we input in the simulations. We find that the standard ΛCDM model predicts a charateristic shape of β for temperature and E mode as a function of the threshold, and the average over thresholds is β∼ 0.62 for temperature and β∼ 0.63 for E mode. Accurate measurements of α and β can be used to test the standard model of cosmology and to search for deviations from it. For this purpose we compute α and β for temperature and E mode data of various data sets from PLANCK mission. We compare the values measured from observed data with those obtained from simulations to which instrument beam and noise characteristics of the 44GHz frequency channel have been added (which are provided as part of the PLANCK data release). We find very good agreement of β and α between all PLANCK temperature data sets with ΛCDM expectations. E mode data show good agreement for β but α for all data sets deviate from ΛCDM predictions higher than 3−σ. It is most likely that the deviations are probing the anisotropy of the noise field and beam characteristics of the detector rather than the true E mode signal since for 44GHz the signal-to-noise ratio is well below one. This will be further investigated after the full PLANCK data becomes publicly available.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Data-driven Modeling of Hemodynamics and its Role on Thrombus Size and Shape in Aortic Dissections

Aortic dissection is a pathology that manifests due to microstructural defects in the aortic wall. Blood enters the damaged wall through an intimal tear, thereby creating a so-called false lumen and exposing the blood to thrombogenic intramural constituents such as collagen. The natural history of this acute vascular injury thus depends, in part, on thrombus formation, maturation, and possible healing within the false lumen. A key question is: Why do some false lumens thrombose completely while others thrombose partially or little at all? An ability to predict the location and extent of thrombus in subjects with dissection could contribute significantly to clinical decision-making, including interventional design. We develop, for the first time, a data-driven particle-continuum model for thrombus formation in a murine model of aortic dissection. In the proposed model, we simulate a final-value problem in lieu of the original initial-value problem with significantly fewer particles that may grow in size upon activation, thus representing the local concentration of blood-borne species. Numerical results confirm that geometry and local hemodynamics play significant roles in the acute progression of thrombus. Despite geometrical differences between murine and human dissections, mouse models can provide considerable insight and have gained popularity owing to their reproducibility. Our results for three classes of geometrically different false lumens show that thrombus forms and extends to a greater extent in regions with lower bulk shear rates. Dense thrombi are less likely to form in high-shear zones and in the presence of strong vortices. Here, the present data-driven study suggests that the proposed model is robust and can be employed to assess thrombus formation in human aortic dissections.

60 APPLIED LIFE SCIENCES↗

Isospin Decomposition of the Photoproduced Sigma pi System Near the Lambda(1405)

Recent experimental results for the reaction \gamma + p \to K^+ + \Sigma + \pi\ from CLAS at Jefferson Lab are discussed. It was found that the mass distributions or "line shapes" of the three charge combinations \Sigma^+ \pi^-, \Sigma^0 \pi^0 and \Sigma^- \pi^+ differ significantly. Our results show that the \Lambda(1405), as the I=0 constituent of the reaction, must be accompanied by an I > 0 component. We discuss phenomenological fits to the data to test the possible forms and magnitudes of these amplitudes. A two-amplitude I=0 fit of Breit-Wigner form to the \Sigma^0\pi^0 channel alone works quite well. The addition of a single I=1 amplitude works fairly well to model all the line shapes simultaneously.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

TEM Approaches for Microstructure-Informed Prediction of Mechanical Properties in Structural Alloys

Predicting the mechanical performance of structural alloys from their evolving microstructure remains a major challenge in materials science, particularly for nuclear structural materials, where irradiation-induced defects span multiple types and length scales and interact through complex mechanisms. The dispersed barrier hardening (DBH) [1] and Friedel–Kroupa–Hirsch (FKH) [2,3] models have been widely used to evaluate the hardening contributions of individual obstacles and to estimate tensile strength from quantified microstructures; however, when multiple size-dependent obstacles coexist and evolve, predicting temperature-dependent tensile strength becomes significantly more complex, and a fully consistent hardening model is still lacking. Transmission electron microscopy (TEM) plays a central role in refining hardening models and enabling predictive assessments of tensile strength evolution by providing quantitative characterization of dislocations, irradiation-induced defects (e.g., dislocation loops and cavities), precipitates, and grain structure (Fig. 1.). These experimentally measured defect densities are incorporated into physically based hardening models with size- and shape- dependent obstacle strengths [4], using root-sum-square superposition for obstacles of comparable strength and linear superposition for dissimilar ones [5]. In addition, recent advances in TEM [6-8], including high-resolution imaging, 4D-STEM strain mapping, EDS/EELS elemental analysis, and flash-polishing-based TEM specimen preparation and extraction-replica methods (Fig. 2), further improve the accuracy of microstructural quantification. By comparison with prior studies as well as our own results, we show that when TEM-derived microstructural information is carefully integrated with physically grounded hardening models, yield strength (or irradiation-induced hardening) measured at room temperature can be predicted with good quantitative agreement across multiple alloy classes. In-situ TEM combined with high-temperature mechanical testing represents an important next step for refining hardening models by directly probing dislocation–obstacle interactions across varying irradiation doses and temperatures [9]. Because the barrier strength factor (α) depends on both temperature and obstacle size, it should not be treated as a constant fitting parameter; rather, it must be explicitly evaluated to achieve physically meaningful predictions of mechanical behaviour at operating temperatures. This presentation therefore discusses why all strengthening contributions (e.g., Peierls stress, solid-solution strengthening, voids, bubbles, dislocation loops, dislocation lines, and grain boundaries) must be considered collectively, why appropriate superposition methods are essential when obstacles possess different barrier strength factors, how hardness measurements can be meaningfully related to tensile properties, and how TEM-derived microstructural information can be systematically incorporated into hardening models. More broadly, it outlines a pathway toward microstructure-informed prediction of mechanical properties and supports the goal of establishing science-based tools for evaluating structural materials in extreme environments [10].

Lin, Yan-Ru [ORNL] (ORCID:0000000339991473)↗

Two-tier model reduction of viscoelastically damped finite element models

Transient simulations of linear viscoelastically damped structures necessitate critical computational resources to directly integrate the full-order finite element model with time-stepping algorithms. Traditional modal reduction techniques are not directly applicable to these systems since viscoelastic materials depend on time and frequency. A more appropriate reduction basis is obtained from the nonlinear, complex eigenvalue problem, whose eigenvectors capture the appropriate kinematics and enable frequency-based mode selection; unfortunately, the computational cost is prohibitive for computing these modes from large-scale engineering models. To address this shortcoming, this work proposes a novel two-tier reduction procedure to reduce the upfront cost of solving the complex, nonlinear eigenvalue problem. The first reduction step reduces the full-order model with real mode shapes linearized about various centering frequencies to capture the kinematics over a full range of viscoelastic material behavior (glassy, rubbery, and glass-transition zones). This tier-one reduction preserves time-temperature superposition and allows the equations to depend parametrically on operating temperature. The second-level reduction then solves the complex, nonlinear eigenmode solutions in the tier-one reduced space about a fixed temperature to further reduce the equations-of-motion. The method is demonstrated on a cantilevered sandwich plate to display its accuracy and efficiency in comparison to full-order model predictions.

97 MATHEMATICS AND COMPUTING↗

On the mechanics of filament winding. I - A generalized model

In this paper the filament winding process was analyzed to relate process variables to structural parameters. An analytical model has been developed for a filament wound preform by determining the necessary geometric and kinematic conditions. An equation of motion of winding was developed to predict layer by layer geometry using the mandrel shape as the initial boundary condition. This model provides the laws of motion of the traverse stroke relative to the mandrel to ensure that the filament laydown is of the predetermined geodesic path over the mandrel, that the winding density is distributed over the mandrel in the prearranged manner, and that the winding build is being formed to reproduce the preassigned shape or contour of the preform. From this analysis, the winding process parameters can be controlled accordingly to produce any given shape. The analytical model was used to produce computer simulation of the winding process. The dynamic simulation of the winding process visualizes and inspects the layer by layer path of the tow and provides process variables required to produce preforms of optimum structural parameters.

Hamouda, Hechmi↗

On the Properties of Cosmological Ionization Fronts

We investigate the properties of cosmological ionization fronts during the epoch of reionization using the Cosmic Reionization on Computers simulations. By analyzing reionization timing maps, we characterize ionization front velocities and curvatures and their dependence on the density structure of the intergalactic medium (IGM). The velocity distribution of ionization fronts in the simulations indicates that while the barrier-crossing analytical model captures the overall shape in high-velocity regions, it fails to reproduce the low-velocity tail, highlighting the non-Gaussian nature of the IGM’s density field. Ionization front velocities are inversely correlated with local density, propagating faster in underdense regions and more slowly in overdense environments. Faster ionization fronts also lead to higher postionization temperatures, reaching a plateau at ∼2 × 10$^{4}$ K for velocities exceeding 3000 km s$^{−1}$. Examining curvature statistics further establishes a connection between ionization front structure and the normalized density contrast ν, with trends in overdense regions aligning well with barrier-crossing model predictions, while deviations appear in underdense environments due to model limitations. These results provide a detailed characterization of ionization front dynamics and their interaction with the underlying density field, bridging small-scale reionization physics with large-scale observables such as the 21 cm signal and the IGM’s thermal history.

79 ASTRONOMY AND ASTROPHYSICS↗

Strengthening and Improving Yield Asymmetry of Magnesium Alloys by Second Phase Particle Refinement Under the Guidance of Integrated Computational Materials Engineering

Improving yield strength and asymmetry is critical to expand applications of magnesium alloys in industry for higher fuel efficiency and lower CO 2 production. Grain refinement is an efficient method for strengthening low symmetry magnesium alloys, achievable by precipitate refinement. This study provides guidance on how precipitate engineering will improve mechanical properties through grain refinement. Precipitate refinement for improving yield strengths and asymmetry is simulated quantitatively by coupling a stochastic second phase grain refinement model and a modified polycrystalline crystal viscoplasticity φ-model. Using the stochastic second phase grain refinement model, grain size is quantitatively determined from the precipitate size and volume fraction. Yield strengths, yield asymmetry, and deformation behavior are calculated from the modified φ-model. If the precipitate shape and size remain constant, grain size decreases with increasing precipitate volume fraction. If the precipitate volume fraction is kept constant, grain size decreases with decreasing precipitate size during precipitate refinement. Yield strengths increase and asymmetry approves to one with decreasing grain size, contributed by increasing precipitate volume fraction or decreasing precipitate size.

36 MATERIALS SCIENCE↗

DES meets Gaia: discovery of strongly lensed quasars from a multiplet search

We report the discovery, spectroscopic confirmation and first lens models of the first two, strongly lensed quasars from a combined search in WISE and Gaia over the DES footprint. Their Einstein radii span a range between ≈2.0 arcsec and ≈0.4 arcsec. Two of these (WGD2038-4008, RA = 20:38:02.65, Dec. = -40:08:14.64; WGD2021-4115, RA = 20:21:39.45, Dec. = -41:15:57.11) also have confirmed deflector redshifts. The four-image lens WGD2038-4008, with source and deflector redshifts s = 0.777 ± 0.001 and zl = 0.230 ± 0.002, respectively, has a deflector with radius R eff ≈ 3.4 arcsec, stellar mass $\log(M_{\star}/M_{\odot}) = 11.64^{+0.20}_{-0.43}$, and shows extended isophotal shape variation. Simple lens models yield Einstein radii $R_{\rm E}=(1.30\pm0.04)^{\prime\prime},$ axis ratio $q=0.75\pm0.1$ (compatible with that of the starlight) and considerable shear-ellipticity degeneracies. The two-image lens WGD2021-4115 (r.a.=20:21:39.45, dec.=--41:15:57.11) has $z_{s}=1.390\pm0.001$ and $z_l = 0.335 \pm 0.002$, and Einstein radius $R_{\rm E} = (1.1\pm0.1)^{\prime\prime},$ but higher-resolution imaging is needed to accurately separate the deflector and faint quasar image. Finally, we also show high-rank candidate doubles selected this way, some of which have been independently identified with different techniques, and discuss a DES+WISE quasar multiplet selection.

79 ASTRONOMY AND ASTROPHYSICS↗

Determination of cloud optical depth with the multiple field of view sun photometer

An approach to determining cloud optical depth is presented which utilizes a multiple-field-of-view sun photometer, along with modeled behavior of transmittance ratios as a function of the field of view, to infer the optical depth. The data used in the analysis were collected during the FIRE Stratocumulus IFO, which took place from July 15 to 17, 1987. Values which were taken within an hour of local solar noon are presented and compared with the values computed for mu sub 0 = 1.0. The zenith angle for the measurements was between 11 and 18 deg. It is found that this method of inferring the optical depth indicates most clouds with optical depths within these limits. The measured curves have a slightly different shape than do the modeled ratios.

Hein, Paul F.↗

The energy spectra of solar flare electrons

A survey of 50 electron energy spectra from .1 to 100 MeV originating from solar flares was made by the combination of data from two spectrometers onboard the International Sun Earth Explorer-3 spacecraft. The observed spectral shapes of flare events can be divided into two classes through the criteria of fit to an acceleration model. This standard two step acceleration model, which fits the spectral shape of the first class of flares, involves an impulsive step that accelerates particles up to 100 keV and a second step that further accelerates these particles up to 100 MeV by a single shock. This fit fails for the second class of flares that can be characterized as having excessively hard spectra above 1 MeV relative to the predictions of the model. Correlations with soft X-ray and meter radio observations imply that the acceleration of the high energy particles in the second class of flares is dominated by the impulsive phase of the flares.

Evenson, P. A.↗