Search NASA⌕ Search

SEARCH · Search NASA

Results for “force method”

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,693 records · Page 94

A Model for Space Shuttle Orbiter Tire Side Forces Based on NASA Landing Systems Research Aircraft Test Results

Forces generated by the Space Shuttle orbiter tire under varying vertical load, slip angle, speed, and surface conditions were measured using the Landing System Research Aircraft (LSRA). Resulting data were used to calculate a mathematical model for predicting tire forces in orbiter simulations. Tire side and drag forces experienced by an orbiter tire are cataloged as a function of vertical load and slip angle. The mathematical model is compared to existing tire force models for the Space Shuttle orbiter. This report describes the LSRA and a typical test sequence. Testing methods, data reduction, and error analysis are presented. The LSRA testing was conducted on concrete and lakebed runways at the Edwards Air Force Flight Test Center and on concrete runways at the Kennedy Space Center (KSC). Wet runway tire force tests were performed on test strips made at the KSC using different surfacing techniques. Data were corrected for ply steer forces and conicity.

Carter, John F.↗

Wind suppression by X-rays in Cygnus X-3

Context: The radiatively driven wind of the primary star in wind-fed X-ray binaries can be suppressed by the X-ray irradiation of the compact secondary star. This causes feedback between the wind and the X-ray luminosity of the compact star. Aims: We aim to estimate how the wind velocity on the face-on side of the donor star depends on the spectral state of the high-mass X-ray binary Cygnus X-3. Methods: We modeled the supersonic part of the wind by computing the line force (force multiplier) with the Castor, Abbott & Klein formalism and XSTAR physics and by solving the mass conservation and momentum balance equations. We computed the line force locally in the wind considering the radiation fields from both the donor and the compact star in each spectral state. We solved the wind equations at different orbital angles from the line joining the stars and took the effect of wind clumping into account. Wind-induced accretion luminosities were estimated using the Bondi-Hoyle-Lyttleton formalism and computed wind velocities at the compact star. We compared them to those obtained from observations. Results: We found that the ionization potentials of the ions contributing the most to the line force fall in the extreme-UV region(100–230 Å). If the flux in this region is high, the line force is weak, and consequently, the wind velocity is low. We found a correlation between the luminosities estimated from the observations for each spectral state of Cyg X-3 and the computed accretion luminosities assuming moderate wind clumping and a low mass of the compact star. For high wind clumping, this correlation disappears. We compared the XSTAR method used here with the comoving frame method and found that they agree reasonably well with each other. Conclusions. We show that soft X-rays in the extreme-UV region from the compact star penetrate the wind from the donor star and diminish the line force and consequently the wind velocity on the face-on side. This increases the computed accretion luminosities qualitatively in a similar manner as observed in the spectral evolution of Cyg X-3 for a moderate clumping volume filling factor and a compact star mass of a few (2–3) solar masses.

O. Vilhu↗

Variational Coupled Loads Analysis using the Hybrid Parametric Variation Method

Time-domain coupled loads analysis (CLA)is used to determine the response of a launch vehicle and payload system to transient forces, such as liftoff, engine ignitions and shutdowns, jettison events, and atmospheric flight loads, such as buffet. CLA, using Hurty/Craig-Bampton (HCB)component models, is the accepted method for the establishment of design-level loads for launch systems. However, uncertainty in the component models flows into uncertainty in predicted system results. Uncertainty in the structural responses during launch is a significant concern because small variations in launch vehicle and payload mode shapes and their interactions can result in significant variations in system loads. Uncertainty quantification (UQ)is used to determine statistical bounds on prediction accuracy based on model uncertainty. In this paper uncertainty is treated at the HCB component-model level. In an effort to account for model uncertainties and statistically bound their effect on CLA predictions, this work combines CLA with UQ in a process termed variational coupled loads analysis (VCLA). The modeling of uncertainty using a parametric approach, in which input parameters are represented by random variables, is common, but its major drawback is the resulting uncertainty is limited to the form of the nominal model. Uncertainty in model form is one of the biggest contributors to uncertainty in complex built-up structures. Model-form uncertainty can be represented using a nonparametric approach based on random matrix theory (RMT). In this work, UQ is performed using the hybrid parametric variation (HPV)method, which combines parametric with nonparametric uncertainty at the HCB component model level. The HPV method requires the selection of dispersion values for the HCB fixed-interface (FI)eigenvalues, and the HCB mass and stiffness matrices. The dispersions are based upon component test-analysis modal correlation results. During VCLA, random component models are assembled into an ensemble of random systems using a Monte Carlo (MC)approach. CLA is applied to each of the ensemble members to produce an ensemble of system-level responses for statistical analysis. The proposed methodology is demonstrated through its application to a buffet loads analysis of NASA’s Space Launch System (SLS)during the transonic regime fifty seconds after liftoff. Core stage (CS)section shears and moments are recovered, and statistics are computed.

Uncertainty Quantification↗

The evolution history of superalloy powders during hot consolidation and plastic deformation

Highlights: • Element segregations in a powder metallurgy superalloy at different states are identified. • The surface energy drives precipitation of oxides, carbides, and γ' particles on powder surface during hot consolidation. • Prior particle boundary structure is distorted and more uniformly refined grains are obtained after hot extrusion. - Abstract: Powder metallurgy superalloys are widely used to produce critical rotating parts in turbine engines. However, the issues like what is the driving force of precipitation on prior particle boundary (PPB) during consolidation and how the PPB structure evolves during plastic deformation still remain to be resolved. In this study, with the characterization of the microstructure evolution from powder, consolidation, to plastic deformation, the segregation of alloying elements on the particle boundaries at elevated temperature was analyzed by equilibrium thermodynamic calculations, and the distortion of PPB structure under large plastic deformation was quantitively simulated by finite element (FE) methods. The results indicated that surface energy serves as the main driving force of PPB precipitation and hot extrusion could distort the PPB structure and yield more uniformly refined grains.

36 MATERIALS SCIENCE↗

Radiative Forcing by Well-Mixed Greenhouse Gases: Estimates from Climate Models in the Intergovernmental Panel on Climate Change (IPCC) Fourth Assessment Report (AR4)

The radiative effects from increased concentrations of well-mixed greenhouse gases (WMGHGs) represent the most significant and best understood anthropogenic forcing of the climate system. The most comprehensive tools for simulating past and future climates influenced by WMGHGs are fully coupled atmosphere-ocean general circulation models (AOGCMs). Because of the importance of WMGHGs as forcing agents it is essential that AOGCMs compute the radiative forcing by these gases as accurately as possible. We present the results of a radiative transfer model intercomparison between the forcings computed by the radiative parameterizations of AOGCMs and by benchmark line-by-line (LBL) codes. The comparison is focused on forcing by CO2, CH4, N2O, CFC-11, CFC-12, and the increased H2O expected in warmer climates. The models included in the intercomparison include several LBL codes and most of the global models submitted to the Intergovernmental Panel on Climate Change (IPCC) Fourth Assessment Report (AR4). In general, the LBL models are in excellent agreement with each other. However, in many cases, there are substantial discrepancies among the AOGCMs and between the AOGCMs and LBL codes. In some cases this is because the AOGCMs neglect particular absorbers, in particular the near-infrared effects of CH4 and N2O, while in others it is due to the methods for modeling the radiative processes. The biases in the AOGCM forcings are generally largest at the surface level. We quantify these differences and discuss the implications for interpreting variations in forcing and response across the multimodel ensemble of AOGCM simulations assembled for the IPCC AR4.

Collins, W. D.↗

Electromagnetic Analysis of ITER Electron Cyclotron Emission Model

The ITER electron cyclotron emission (ECE) diagnostic system is located at diagnostics shielding module 2 (DSM2), equatorial port 9 (EP9), to measure electron temperature profile and electron temperature fluctuations, and also assess nonthermal electron distributions via the oblique view. Therefore, ECE has both radial and oblique views with two hot sources for calibration and two couples of mirrors for two different optical views. It shall receive that the large electromagnetic (EM) loads up to 100-MN/m3 force density due to the eddy currents on the copper mirrors generated by the short transient plasma disruption. The simplified EM analysis model using magnetic field ( B ) data and flux variations (dB/dt) method based on the worst case of plasma disruption, MD_DW_EXP16MS_CATIII, has been developed to calculate the force and moments for each ECE component and support structure for EP9 DSM2 Bay2 and Bay3. The B and dB/dt methods use the constant B and dB/dt assumed for each local bay space, but the real B and dB/dt are basically varied along the radial direction of each bay. The complex global model with the local ECE components and support structure using Maxwell transient has been developed to do EM analysis to compare the results of the B and dB/dt models. The volumetric EM force density of the whole ECE components and support structure can be used for the subsequent structural analysis. In conclusion, combined with the thermal and nuclear loads and seismic and initial loads, the ECE integration analysis can be finalized for the operation case.

Fang, J. [Princeton Plasma Physics Laboratory (PPP↗

Improvement of Numerical Methods in Pronghorn

During the fiscal year of 2021, the finite-volume method (FVM) was deployed in Pronghorn to provide improved efficiency, stability, and accuracy for coarse-mesh, thermal-hydraulics problems. While the main goals of fiscal year 2021 were met, several issues emerged from the early deployment of the finite-volume method in Pronghorn. These issues were: The compressible and incompressible formulations that were implemented are inadequate for many nuclear reactor flow problems; Omission of terms accounting for the porosity and Darcy-Forchheimer body force discontinuities in the Rhie-Chow interpolation lead to oscillations in pressure and velocity at these discontinuities; The FVM lacks a correction for non-orthogonal grids for computing accurate pressure gradient leading to loss of accuracy in regions with skewed elements; The FVM currently uses a monolithic solver. Monolithic solvers have issues dealing with the saddle-point nature of the discretized fluid equations. The result are bad convergence if direct factorization is not used and large memory consumption when direct factorization is used. These four issues are addressed in this report. In particular, we report the completion of the following task: Implementation of a weakly compressible formulation in the MOOSE Navier-Stokes module; Implementation of Moukalled’s method for including body forces in the RCI. Additionally, we identified the need to smooth the porosity using Moukalled’s face-cell smoothing operator; Implementation of a non-orthogonal correction for the Green-Gauss gradient computation; Preliminary implementation of a SIMPLE segregated solver.

97 MATHEMATICS AND COMPUTING↗

Magnetic Properties and Electronic Structure of Co/Zr Multilayers

We investigate the magnetic and electronic properties of Co/Zr multilayers with an emphasis on the interfacial effects. The saturation magnetization, coercivity, and uniaxial anisotropy depend significantly on Co thickness. Spin-wave resonance modes were observed, and the change of resonant field as a function of mode number (n) obeys the so-called n{sup 2} law. The interlayer coupling constant was determined. The electronic properties of the interfacial alloys were studied using the WIEN2k code, employing the density functional theory. Calculations describe satisfactorily the observed magnetic behaviour. The influence of pressure on the magnetic properties of the interfacial Co{sub 1−x}Zr{sub x} alloys was investigated in which we demonstrate that the ferromagnetism is altered by pressure.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Integrating helicoid channels for passive control of fiber alignment in direct-write 3D printing

3D printing of fiber-reinforced composites has been receiving increasing attention as an efficient additive approach enabling lightweight, functional, and high-performance components required for industrial applications. The properties of fiber composites significantly depend on internal microstructures, including fiber orientation, distribution, and degree of alignment. Although multiple strategies have been introduced for controlling the fiber composites microstructures, these strategies are mostly active approaches relying on additional control of mechatronics parts, external forces, or intentionally introduced unstable flows inside narrow nozzles. These methods suffer from higher risk of nozzle clogging and requiring extra parts in 3D printing. Here, we introduce helicoid channels into extrusion system as a passive approach to control the fiber alignment without any additional moving parts in direct-write 3D printing. The helicoid channels automatically guide the composite inks and align fibers before flowing into narrow nozzle space, avoiding clogging and improving printability. The analyses indicate that both helicoid surface to volume ratio and helix angle affect the pre-alignment of fibers, leading to tunable mechanical properties of printed fiber-reinforced composites with increased stiffness and strength up to 77.6% and 47.8%, respectively.

36 MATERIALS SCIENCE↗

Tribological properties of PDA + PTFE coating in oil-lubricated condition

In this work the tribological properties of 45 ± 2 μm-thick polytetrafluoroethylene (PTFE) and polydopamine (PDA) + PTFE coatings in oil-lubricated conditions were investigated under a normal load of 10 N and a linear reciprocating speed of 0.1 m/s. Both the PTFE and the PDA + PTFE coatings were deposited using an in-house spray-coating method. The PDA + PTFE coating was five times more durable compared to the PTFE coating in oil-lubricated conditions. Both PTFE and PDA + PTFE coatings showed lower COF in oil-lubricated conditions than in dry conditions. The nanomechanical properties measured by PeakForce Quantitative Nanomechanical (PFQNM) characterization method showed an increase in Young’s modulus and adhesion force for the PDA + PTFE coatings compared to the PTFE coatings. Chemical analysis showed that the cross-linking between PDA and PTFE polymer chains occurred during the high-temperature sintering procedure. Finally, the higher adhesion of PDA + PTFE coating to the cast iron substrate and the stronger cohesion due to the cross-linking between PDA and PTFE contributed to the higher durability of the PDA + PTFE composite coatings.

36 MATERIALS SCIENCE↗

Monte Carlo simulations for understanding the transport properties of metallic nanowires

We report the conductance calculation of metallic nanowires of different metals under the application of a tensile force. The elongation induced by the tensile force was simulated by means of the canonical Monte Carlo method combined with the embedded atom method to compute the energy of the system. We studied nanowires composed of Ag, Au, Pt and Cu. The conductance was calculated by means of the Landauer formula and the non-equilibrium Green function formalism. As previously verified experimentally and using other simulation techniques, the conductance is quantized. In addition to this, we have found that there is a strong correlation between the magnitude of the conductance and the width of the nanowire at its thinnest point. Conductance histograms were analyzed for all the different metals.

36 MATERIALS SCIENCE↗

Off-Lattice Markov Chain Monte Carlo Simulations of Mechanically Driven Polymers

Here, we develop off-lattice simulations of semiflexible polymer chains subjected to applied mechanical forces by using Markov Chain Monte Carlo. Our approach models the polymer as a chain of fixed length bonds, with configurations updated through adaptive nonlocal Monte Carlo moves. This proposed method enables precise calculation of a polymer’s response to a wide range of mechanical forces, which traditional on-lattice models cannot achieve. Our approach has shown excellent agreement with theoretical predictions of persistence length and end-to-end distance in quiescent states as well as stretching distances under tension. Moreover, our model eliminates the orientational bias present in on-lattice models, which significantly impacts calculations such as the scattering function, a crucial technique for revealing the polymer conformation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Using Neural Network Force Fields to Ascertain the Quality of Ab Initio Simulations of Liquid Water

Accurately simulating the properties of bulk water, despite the apparent simplicity of the molecule, is still a challenge. In order to fully understand and reproduce its complex phase diagram, it is necessary to perform simulations at the ab initio level, including quantum mechanical effects both for electrons and nuclei. This comes at a high computational cost, given that the structural and dynamical properties tend to require long timescales and large simulation cells. In this work, we evaluate the errors that density functional theory (DFT)-based simulations routinely incur into due time- and size-scale limitations. These errors are evaluated using neural-network-trained force fields that are accurate at the level of DFT methods. We compare different exchange and correlation potentials for properties of bulk water that require large timescales. Here we show that structural properties are less dependent on the system size and that dynamical properties such as the diffusion coefficient have a strong dependence on the simulation size and timescale. Our results facilitate comparisons of DFT-based simulation results with experiments and offer a path to discriminate between model and convergence errors in these simulations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Cross-Scale Molecular Analysis of Chemical Heterogeneity in Shale Rocks

The organic and mineralogical heterogeneity in shale at micrometer and nanometer spatial scales contributes to the quality of gas reserves, gas flow mechanisms and gas production. Here, we demonstrate two molecular imaging approaches based on infrared spectroscopy to obtain mineral and kerogen information at these mesoscale spatial resolutions in large-sized shale rock samples. The first method is a modified microscopic attenuated total reflectance measurement that utilizes a large germanium hemisphere combined with a focal plane array detector to rapidly capture chemical images of shale rock surfaces spanning hundreds of micrometers with micrometer spatial resolution. The second method, synchrotron infrared nano-spectroscopy, utilizes a metallic atomic force microscope tip to obtain chemical images of micrometer dimensions but with nanometer spatial resolution. This chemically "deconvoluted" imaging at the nano-pore scale is then used to build a machine learning model to generate a molecular distribution map across scales with a spatial span of 1000 times, which enables high-throughput geochemical characterization in greater details across the nano-pore and micro-grain scales and allows us to identify co-localization of mineral phases with chemically distinct organics and even with gas phase sorbents. Finally, this characterization is fundamental to understand mineral and organic compositions affecting the behavior of shales.

04 OIL SHALES AND TAR SANDS↗

Nonlinear saturation of ballooning modes in stellarators

Ballooning mode saturation is investigated in realistic stellarator configurations using the flux tube approach of Ham et al (2018 Plasma Phys. Control. Fusion 60 075017), Ham et al (2016 Phys. Rev. Lett. 116 235001). The method is adapted to account for the lack of exact force balance in stellarator equilibrium solvers that assume existence of nested flux surfaces. A variational approach for calculating flux tube energy is developed to overcome this force error problem in stellarator numerical equilibria. Saturated (equilibrium) flux tube states that cross 10%–20% of the plasma minor radius are shown to exist for linearly ballooning unstable profiles. It is shown that several features of the displaced flux tube structure in a full nonlinear MHD simulation of Wendelstein 7X are reproduced by our model. Saturated states are found in a compact stellarator equilibrium close but below the marginal ballooning linear instability, i.e. the unperturbed equilibrium is metastable. This suggests that edge-localized-mode-like explosive MHD behavior may be possible in stellarators.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Structure of Offshore Low-Level Jet Turbulence and Implications to Mesoscale-to-Microscale Coupling

This paper explores realistic nonstationary atmospheric boundary layer (ABL) turbulence arising from nonstationarity at the mesoscale, particularly within offshore low-level jets with implications to offshore wind farms, using high-fidelity multiscale large-eddy simulations (LES). To this end, we analyzed the single-point turbulence statistical structure of a North-Atlantic offshore LLJ event simulated using high-resolution LES (AMR-Wind). The nonstationary LLJ is simulated using a mesoscale-to-microscale coupled (MMC) simulation procedure involving data assimilation of mesoscale velocity and temperature data from the Weather Research and Forecasting (WRF) model. Unlike the assimilation of mesoscale velocity data into the LES, the direct assimilation of temperature profiles had a strong impact on turbulence stratification, thereby causing erroneous predictions of turbulence both above and within the jet layer. Various approaches to mitigate this effect have resulted in multiple (four) variants of this MMC strategy. Outcomes from this work clearly show that the turbulence within the low-level jet is a strong function of the MMC approach as the turbulence structure within the low-level jet is dependent on the flux of residual turbulence from outside the jet, which in turn depends on the temperature forcing history. Additionally, the turbulence predicted by all these different methods (as well as the observation data) show similar deviations from equilibrium as evidenced by comparisons with idealized atmospheric turbulence structure obtained using the same numerical method. In general, we observe that the predicted LLJ turbulence tends to differ from canonical ABL turbulence with comparable shear. Particularly, the combination of shear and turbulence observed in such nonstationary low-level turbulence cannot be matched using equilibrium settings and therefore, represents a critical use-case for both testing and leveraging meso–micro coupling strategies.

17 WIND ENERGY↗

Phonon thermal conductance across GaN-AlN interfaces from first principles

The vibrational thermal conductances (G) across GaN-AlN interfaces are computed using a nonequilibrium Green's function formalism in the harmonic limit with bulk and interfacial interatomic force constants (IFCs) fully from density functional theory. Several numerical methods and supercell configurations are employed to examine the sensitivity of G to variances of IFCs. In particular, the effects of supercell size, the enforcement of symmetry constraints, and truncation of IFCs near the interface, and atomic relaxation on phonon transmission and conductance are explored. Our fully first-principles calculations are compared with common approximations and measured G values inferred from thermal conductivity measurements for GaN-AlN superlattices. Here, our calculated value, G~300MWm -2 K -1 , is nearly half that from measurements. This discrepancy is critically analyzed in terms of the physical assumptions of the calculations and the derivation of the experimental values. This work provides guidelines to determine “physically correct” sets of interfacial IFCs from first principles for thermal conductance calculations using minimal computational resources. It also contributes toward developing predictive calculations and a more complete picture of thermal conduction across interfaces, a step toward first-principles multiscale thermal transport.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Uncertainty quantification in breakup reactions

Breakup reactions are one of the favored probes to study loosely bound nuclei, particularly those in the limit of stability forming a halo. In order to interpret such breakup experiments, the continuum discretized coupled channel method is typically used. In this study, the first Bayesian analysis of a breakup reaction model is performed. Herein we use a combination of statistical methods together with a three-body reaction model (the continuum discretized coupled channel method) to quantify the uncertainties on the breakup observables due to the parameters in the effective potential describing the loosely bound projectile of interest. The combination of tools we develop opens the path for a Bayesian analysis of not only breakup processes, but also a wide array of complex processes that require computationally intensive reaction models.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗