SEARCH · Engineering Papers

Results for “consistent scheme”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 901 records · Page 50

All sources NASA NTRS DOE OSTI

NASA NTRS1990

The resupply interface mechanism RMS compatibility test

Spacecraft on-orbit servicing consists of exchanging components such as payloads, orbital replacement units (ORUs), and consumables. To accomplish the exchange of consumables, the receiving vehicle must mate to the supplier vehicle. Mating can be accomplished by a variety of docking procedures. However, these docking schemes are mission dependent and can vary from shuttle bay berthing to autonomous rendezvous and docking. Satisfying the many docking conditions will require use of an innovative docking device. The device must provide fluid, electrical, pneumatic and data transfer between vehicles. Also, the proper stiffness must be obtained and sustained between the vehicles. A device to accomplish this, the resupply interface mechanism (RIM), was developed. The RIM is a unique device because it grasps the mating vehicle, draws the two vehicles together, simultaneously mates all connectors, and rigidizes the mating devices. The NASA-Johnson Manipulator Development Facility was used to study how compatible the RIM is to on orbit docking and berthing. The facility contains a shuttle cargo bay mockup with a remote manipulator system (RMS). This RMS is used to prepare crew members for shuttle missions involving spacecraft berthing operations. The MDF proved to be an excellant system for testing the RIM/RMS compatibility. The elements examined during the RIM JSC test were: RIM gross and fine alignment; berthing method sequence; visual cuing aids; utility connections; and RIM overall performance. The results showed that the RIM is a good device for spacecraft berthing operations. Mating was accomplished during every test run and all test operators (crew members) felt that the RIM is an effective device. The purpose of the JSC RIM test and its results are discussed.

Jackson, Stewart W.↗

NASA NTRS1993

Global/local interlaminar stress analysis of a grid-stiffened composite panel

A global/local procedure for the computation of the interlaminar stress components at the skin wrap, skin core, and wrap core interfaces for an advanced concept stiffened panel, is described. The procedure consists of a global model of two dimensional shell elements that is used to design a grid stiffened panel with blade type stiffeners, a local model of three dimensional solid elements that is used to compute interlaminar stress components, and a scheme devised to assign displacement boundary conditions for a local model that are based on displacement and rotation data of a few nodes of the global model. A global panel was designed according to strength, stiffness, and stability criteria associated with the design of traditional aircraft wing panels. Interlaminar normal and shearing stress components, computed via the local model, were found to be well below typical tensile normal and shearing strengths of a graphite epoxy material.

Wiggenraad, J. F. M.↗

NASA NTRS1998

Numerical Experiment with Time and Spatial Accuracy of Navier-Stokes Computation For Helicopter Problems

Helicopter flowfields are highly unsteady, nonlinear and three-dimensional. In forward flight and in hover, the rotor blades interact with the tip vortex and wake sheet developed by either itself or the other blades. This interaction, known as blade-vortex interactions (BVI), results in unsteady loading of the blades and can cause a distinctive acoustic signature. Accurate and cost-effective computational fluid dynamic solutions that capture blade-vortex interactions can help rotor designers and engineers to predict rotor performance and to develop designs for low acoustic signature. Such a predictive method must preserve a blade's shed vortex for several blade revolutions before being dissipated. A number of researchers have explored the requirements for this task. This paper will outline some new capabilities that have been added to the NASA Ames' OVERFLOW code to improve its overall accuracy for both vortex capturing and unsteady flows. To highlight these improvements, a number of case studies will be presented. These case studies consist of free convection of a 2-dimensional vortex, dynamically pitching 2-D airfoil including light-stall, and a full 3-D unsteady viscous solution of a helicopter rotor in forward flight In this study both central and upwind difference schemes are modified to be more accurate. Central difference scheme is chosen for this simulation because the flowfield is not dominated by strong shocks. The feature of shock-vortex interaction in such a flow is less important than the dominant blade-vortex interaction. The scheme is second-order accurate in time and solves the thin-layer Navier-Stokes equations in fully-implicit manner at each time-step. The spatial accuracy is either second and fourth-order central difference or third-order upwind difference using Roe-flux and MUSCLE scheme. This paper will highlight and demonstrate the methods for several sample cases and for a helicopter rotor. Preliminary computations on a rotor were performed by using this method and are in the process of documentation.

Ahmad, Jasim↗

NASA NTRS2004

Bayesian Vision for Shape Recovery

We present a new Bayesian vision technique that aims at recovering a shape from two or more noisy observations taken under similar lighting conditions. The shape is parametrized by a piecewise linear height field, textured by a piecewise linear irradiance field, and we assume Gaussian Markovian priors for both shape vertices and irradiance variables. The observation process. also known as rendering, is modeled by a non-affine projection (e.g. perspective projection) followed by a convolution with a piecewise linear point spread function. and contamination by additive Gaussian noise. We assume that the observation parameters are calibrated beforehand. The major novelty of the proposed method consists of marginalizing out the irradiances considered as nuisance parameters, which is achieved by Laplace approximations. This reduces the inference to minimizing an energy that only depends on the shape vertices, and therefore allows an efficient Iterated Conditional Mode (ICM) optimization scheme to be implemented. A Gaussian approximation of the posterior shape density is computed, thus providing estimates both the geometry and its uncertainty. We illustrate the effectiveness of the new method by shape reconstruction results in a 2D case. A 3D version is currently under development and aims at recovering a surface from multiple images, reconstructing the topography by marginalizing out both albedo and shading.

Jalobeanu, Andre↗

DOE OSTI2020

Density-Matrix Based Extended Lagrangian Born–Oppenheimer Molecular Dynamics

Extended Lagrangian Born–Oppenheimer molecular dynamics [ Phys. Rev. Lett. 2008, 100, 123004] is presented for Hartree–Fock theory, where the extended electronic degrees of freedom are represented by a density matrix, including fractional occupation numbers at elevated electronic temperatures. In contrast to regular direct Born–Oppenheimer molecular dynamics simulations, no iterative self-consistent field optimization is required prior to the force evaluations. To sample regions of the potential energy landscape where the gap is small or vanishing, which leads to particular convergence problems in regular direct Born–Oppenheimer molecular dynamics simulations, an adaptive integration scheme for the extended electronic degrees of freedom is presented. The integration scheme is based on a tunable, low-rank approximation of a fourth-order kernel, K, that determines the metric tensor, T ≡ K T K, used in the extended harmonic oscillator of the Lagrangian that generates the dynamics of the electronic degrees of freedom. Here, the formulation and algorithms provide a general guide to implement extended Lagrangian Born–Oppenheimer molecular dynamics for quantum chemistry, density functional theory, and semiempirical methods using a density matrix formalism.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Quandary

Quandary numerically simulates and optimizes the time-evolution of open quantum systems. The underlying dynamics are modelled by Lindblad's master equation, a linear ordinary differential equation (ODE) describing quantum systems interacting with the environment. Quandary solves this ODE numerically by applying a time-stepping integration scheme, and utilizes a gradient-based optimization approach to determine optimal control pulses that drive the quantum system to a desired target state. Two optimization objectives are considered: (a) Unitary gate optimization that finds controls to realize a unitary gate transformation, and (b) optimal reset that aims to drive the quantum system to the ground states. Gradient-based optimization schemes utilizing Petsc's Tao optimization package are applied to generate control pulses that minimize the respective measure. To evaluate the gradient of the objective function, the discrete adjoint method is used while leveraging techniques from Algorithmic Differentiation to produce exact and consistent gradients. To mitigate excessive execution run times, the software can be build together with the XBraid software library which provides a parallelization strategy to distribute the time-evolution of the underlying dynamics onto multiple processor.

Petersson, NilsA.↗

NASA NTRS1973

An analysis of a charring ablator with thermal nonequilibrium, chemical kinetics, and mass transfer

The differential equations governing the transient response of a one-dimensional ablative thermal protection system are presented for thermal nonequilibrium between the pyrolysis gases and the char layer and with finite rate chemical reactions occurring. The system consists of three layers (the char layer, the uncharred layer, and an optical insulation layer) with concentrated heat sinks at the back surface and between the second and third layers. The equations are solved numerically by using a modified implicit finite difference scheme to obtain solutions for the thickness of the charred and uncharred layers, surface recession and pyrolysis rates, solid temperatures, porosity profiles, and profiles of pyrolysis-gas temperature, pressure, composition, and flow rate. Good agreement is obtained between numerical results and exact solutions for a number of simplified cases. The complete numerical analysis is used to obtain solutions for an ablative system subjected to a constant heating environment. Effects of thermal, chemical, and mass transfer processes are shown.

Clark, R. K.↗

NASA NTRS1974

ATM solar array in-flight performance analysis

The physical and electrical characteristics of the Apollo Telescope Mount (ATM) solar array are described and in-flight performance data are analyzed and compared with predicted results. Two solar cell module configurations were used. Type I module consists of 228 2 x 6 cm solar cells with two cells in parallel and 114 cells in series. Type II modules contain 684 2 x 2 cm cells with six cells in parallel and 114 cells in series. A different interconnection scheme was used for each type. Panels using type II modules with mesh interconnect system performed marginally better than those using type I module with loop interconnect system. The average degradation rate for the ATM array was 8.2% for a 271-day mission.

Thornton, J. P.↗

NASA NTRS1975

Three-dimensional finite element analysis for high velocity impact

A finite element algorithm for solving unsteady, three-dimensional high velocity impact problems is presented. A computer program was developed based on the Eulerian hydroelasto-viscoplastic formulation and the utilization of the theorem of weak solutions. The equations solved consist of conservation of mass, momentum, and energy, equation of state, and appropriate constitutive equations. The solution technique is a time-dependent finite element analysis utilizing three-dimensional isoparametric elements, in conjunction with a generalized two-step time integration scheme. The developed code was demonstrated by solving one-dimensional as well as three-dimensional impact problems for both the inviscid hydrodynamic model and the hydroelasto-viscoplastic model.

Chan, S. T. K.↗

NASA NTRS2008

Single-mode, All-Solid-State Nd:YAG Laser Pumped UV Converter

In this paper, the status of a high-energy, all solid-state Nd:YAG laser pumped nonlinear optics based UV converter development is discussed. The high-energy UV transmitter technology is being developed for ozone sensing applications from space based platforms using differential lidar technique. The goal is to generate greater than 200 mJ/pulse with 10-50 Hz PRF at wavelengths of 308 nm and 320 nm. A diode-pumped, all-solid-state and single longitudinal mode Nd:YAG laser designed to provide conductively cooled operation at 1064 nm has been built and tested. Currently, this pump laser provides an output pulse energy of >1 J/pulse at 50 Hz PRF and a pulsewidth of 22 ns with an electrical-to-optical system efficiency of greater than 7% and a M(sup 2) value of <2. The single frequency UV converter arrangement basically consists of an IR Optical Parametric Oscillator (OPO) and a Sum Frequency Generator (SFG) setups that are pumped by 532 nm wavelength obtained via Second Harmonic Generation (SHG). In this paper, the operation of an inter cavity SFG with CW laser seeding scheme generating 320 nm wavelength is presented. Efforts are underway to improve conversion efficiency of this mJ class UV converter by modifying the spatial beam profile of the pump laser.

Prasad, Narasimha S.↗

NASA NTRSDate not supplied

Wall-Modeled Large Eddy Simulation Method for Unstructured-Grid Navier-Stokes Solvers

This paper reports on the implementation and assessment of a Wall-Modeled Large-Eddy Simulation (WMLES) methodology in an unstructured-grid, node-centered flow solver, FUN3D that is developed and supported at the NASA Langley Research Center. Finite-volume (FV) and finite-element (FE) discretization schemes considered in the study provide formal second-order spatial accuracy. Large-Eddy Simulations (LES) resolve large-scale turbulent-flow features and filter out small-scale effects using the Vreman subgrid-scale model. At solid-wall boundaries, a shear-stress model is employed to provide a proper boundary-flux closure. The nonlinear equations are integrated in time using either an optimized backward difference formula or an implicit multistage Runge-Kutta temporal scheme. The implicit equations at each time step are solved by strong nonlinear iteration schemes. WMLES demonstrations are shown for two high-lift configurations, namely, the McDonnell Douglas 30P30N multielement airfoil and a NASA High-Lift Common Research Model. Results show that the WMLES approaches implemented in the FV and FE discretization methods produce consistent solutions and are capable of capturing key aerodynamic characteristics and flow structures for high-lift configurations at a wide range of angles of attack including maximum-lift conditions. In the 30P30N example, correct trends in the variations of integrated aerodynamic forces and moments, surface pressure distributions, and boundary-layer profiles are captured as the Reynolds number is increased.

CFD; turbulence modeling; High-Lift flow simulatio↗

NASA NTRS2022

LGM Paleoclimate Constraints Inform Cloud Parameterizations and Equilibrium Climate Sensitivity in CESM2

The Community Earth System Model version 2 (CESM2) simulates a high equilibrium climate sensitivity (ECS > 5°C) and a Last Glacial Maximum (LGM) that is substantially colder than proxy temperatures. In this study, we examine the role of cloud parameterizations in simulating the LGM cooling in CESM2. Through substituting different versions of cloud schemes in the atmosphere model, we attribute the excessive LGM cooling to the new CESM2 schemes of cloud microphysics and ice nucleation. Further exploration suggests that removing an inappropriate limiter on cloud ice number (NoNimax) and decreasing the time-step size (substepping) in cloud microphysics largely eliminate the excessive LGM cooling. NoNimax produces a more physically consistent treatment of mixed-phase clouds, which leads to an increase in cloud ice content and a weaker shortwave cloud feedback over mid-to-high latitudes and the Southern Hemisphere subtropics. Microphysical substepping further weakens the shortwave cloud feedback. Based on NoNimax and microphysical substepping, we have developed a paleoclimate-calibrated CESM2 (PaleoCalibr), which simulates well the observed twentieth century warming and spatial characteristics of key cloud and climate variables. PaleoCalibr has a lower ECS (∼4°C) and a 20% weaker aerosol-cloud interaction than CESM2. PaleoCalibr represents a physically more consistent treatment of cloud microphysics than CESM2 and is a valuable tool in climate change studies, especially when a large climate forcing is involved. Our study highlights the unique value of paleoclimate constraints in informing the cloud parameterizations and ultimately the future climate projection.

equilibrium climate sensitivity↗

NASA NTRS2004

Robust Timing Synchronization in Aeronautical Mobile Communication Systems

This work details a study of robust synchronization schemes suitable for satellite to mobile aeronautical applications. A new scheme, the Modified Sliding Window Synchronizer (MSWS), is devised and compared with existing schemes, including the traditional Early-Late Gate Synchronizer (ELGS), the Gardner Zero-Crossing Detector (GZCD), and the Sliding Window Synchronizer (SWS). Performance of the synchronization schemes is evaluated by a set of metrics that indicate performance in digital communications systems. The metrics are convergence time, mean square phase error (or root mean-square phase error), lowest SNR for locking, initial frequency offset performance, midstream frequency offset performance, and system complexity. The performance of the synchronizers is evaluated by means of Matlab simulation models. A simulation platform is devised to model the satellite to mobile aeronautical channel, consisting of a Quadrature Phase Shift Keying modulator, an additive white Gaussian noise channel, and a demodulator front end. Simulation results show that the MSWS provides the most robust performance at the cost of system complexity. The GZCD provides a good tradeoff between robustness and system complexity for communication systems that require high symbol rates or low overall system costs. The ELGS has a high system complexity despite its average performance. Overall, the SWS, originally designed for multi-carrier systems, performs very poorly in single-carrier communications systems. Table 5.1 in Section 5 provides a ranking of each of the synchronization schemes in terms of the metrics set forth in Section 4.1. Details of comparison are given in Section 5. Based on the results presented in Table 5, it is safe to say that the most robust synchronization scheme examined in this work is the high-sample-rate Modified Sliding Window Synchronizer. A close second is its low-sample-rate cousin. The tradeoff between complexity and lowest mean-square phase error determines the rankings of the Gardner Zero-Crossing Detector and both versions of the Early-Late Gate Synchronizer. The least robust models are the high and low-sample-rate Sliding Window Synchronizers. Consequently, the recommended replacement synchronizer for NASA's Advanced Air Transportation Technologies mobile aeronautical communications system is the high-sample-rate Modified Sliding Window Synchronizer. By incorporating this synchronizer into their system, NASA can be assured that their system will be operational in extremely adverse conditions. The quick convergence time of the MSWS should allow the use of high-level protocols. However, if NASA feels that reduced system complexity is the most important aspect of their replacement synchronizer, the Gardner Zero-Crossing Detector would be the best choice.

Xiong, Fu-Qin↗

NASA NTRSDate not supplied

A Holistic DSMC Transport Database for Re-Entry and Ablation Modeling

Hybrid simulation frameworks combining Computational Fluid Dynamics (CFD) and Direct Simulation Monte Carlo (DSMC) are frequently employed to efficiently perform high-fidelity solutions of environments containing combined continuum/rarified flow. The use of DSMC, a stochastic, particle-based method, is necessary for high-Knudsen flow where continuum-based assumptions governing CFD break down. However, the DSMC methodology is generally very computationally inefficient to model the continuum regime. In a CFD/DSMC hybrid approach, obtaining an accurate, high-fidelity solution hinges on the consistent treatment of transport properties and the used thermo-chemical models employed within the two solvers. In principle, in regions where CFD and DSMC are both employed, the same gas mixture under the same conditions should have the same properties, regardless of simulation type. Observed differences should be due to non-equilibrium processes, rather than differences in physical models. While the transport models governing CFD and DSMC simulations are starkly different, they can effectively be linked via their use of reduced Chapman-Enskog collision integrals. In CFD, these integrals are typically stored as fitted polynomial expressions and used to directly compute gas transport properties via mixing rules or the full Chapman-Enskog formulation. In DSMC, they can be used to derive the collision parameters needed for the phenomenological collision cross-section models that govern particle interactions, via a Nelder-Mead optimization scheme. The goal of this work is to provide a unified DSMC transport database encompassing the vast majority of known gas species encountered during atmospheric entry, on Earth or any other Solar body. This goal is largely possible due to recently performed ab-initio quantum chemistry calculations. Combined with other high-fidelity literature sources, the planned database will consist of collision integral data for over 200 neutral and ionized species and over 17000 binary collisions. From these collision integrals, Nelder-Mead optimization is used to compute Variable Soft Sphere (VSS) collision model parameters for DSMC, fitted from 300 K to 20000 K. Initial comparisons of transport properties of relevant equilibrium gas mixtures show great agreement between CFD and DSMC-derived results. The completed database will be able to be readily applied to model binary collisions of any gas mixture containing the included species over the specified temperature range, making it a valuable tool for future planetary probe modeling efforts. An example is shown below. Equilibrium mixture transport properties for a 19-species Titan atmospheric model [4] are computed using both fitted VSS parameters and the original CFD collision integral values. Deviations in computed properties between the two approaches is less than 5% for the entire temperature range.

M R Gosma↗

NASA NTRS2003

Validation of an In-Water, Tower-Shading Correction Scheme

Large offshore structures used for the deployment of optical instruments can significantly perturb the intensity of the light field surrounding the optical measurement point, where different portions of the visible spectrum are subject to different shadowing effects. These effects degrade the quality of the acquired optical data and can reduce the accuracy of several derived quantities, such as those obtained by applying bio-optical algorithms directly to the shadow-perturbed data. As a result, optical remote sensing calibration and validation studies can be impaired if shadowing artifacts are not fully accounted for. In this work, the general in-water shadowing problem is examined for a particular case study. Backward Monte Carlo (MC) radiative transfer computations- performed in a vertically stratified, horizontally inhomogeneous, and realistic ocean-atmosphere system are shown to accurately simulate the shadow-induced relative percent errors affecting the radiance and irradiance data profiles acquired close to an oceanographic tower. Multiparameter optical data processing has provided adequate representation of experimental uncertainties allowing consistent comparison with simulations. The more detailed simulations at the subsurface depth appear to be essentially equivalent to those obtained assuming a simplified ocean-atmosphere system, except in highly stratified waters. MC computations performed in the simplified system can be assumed, therefore, to accurately simulate the optical measurements conducted under more complex sampling conditions (i.e., within waters presenting moderate stratification at most). A previously reported correction scheme, based on the simplified MC simulations, and developed for subsurface shadow-removal processing of in-water optical data taken close to the investigated oceanographic tower, is then validated adequately under most experimental conditions. It appears feasible to generalize the present tower-specific approach to solve other optical sensor shadowing problems pertaining to differently shaped deployment platforms, and also including surrounding structures and instrument casings.

Hooker, Stanford B.↗

DOE OSTI2025

A Poromechanical Framework for Internal Interactions Induced by Solid Inclusions

The framework of poromechanics is generalized to simulate the multiscale behavior of porous media subjected to internal loadings stemming from the growth of solid inclusions. This generalization is designed to enable the study of anisotropic internal stress generation from solid growth within the pores, while recovering isotropic fluid-induced loading as a particular case. For this purpose, a mathematical strategy to define constitutive tensors in a thermodynamically consistent form is proposed, thus offering new opportunities for determining the poromechanical properties of a porous solid through advanced experimentation or micromechanical models. The framework is specialized by means of established elastic solutions for single pore–matrix interaction, as well as through homogenization schemes considering the interaction among congruent pores. In particular, the second Eshelby solution and the Tanaka–Mori–Benveniste homogenization scheme are used to derive a microporoelastic model. At an elemental scale, the model is tested under mixed control conditions by replicating different scenarios of geomaterial testing. In addition, the model characteristics are outlined with reference to inelastic microscopic loadings replicating chemo-mechanical forcing, such as expansive crystal formation. Through a series of parametric analyses, it is shown that the microstructure of the pores significantly influences the properties of porous media. Most notably, it is shown that the effects of a solid forming within the pores depend in a highly nonlinear fashion on the constitutive characteristics of the inhomogeneities and can therefore not be readily quantified or predicted without models capturing the diverse multiscale interactions among pores, inhomogeneities, and matrix.

chemomechanics↗

DOE OSTI2025

Low-mode nonuniformity in direct-drive ICF implosions due to laser smoothing techniques employed on OMEGA

For successful laser-direct-drive inertial confinement fusion implosions, the laser irradiation must be highly uniform over the target surface. On OMEGA, multiple laser beams are used to illuminate targets quasi-uniformly. High-mode-number nonuniformities due to laser speckle on each individual beam are reduced by splitting each beam into two orthogonal polarizations (i.e., polarization smoothing, or PS) and a range of wavelengths (i.e., smoothing by spectral dispersion) that are dispersed at the target plane. However, cross-beam energy transfer (CBET) is sensitive to both the polarizations and wavelengths of the interacting beams, so the interplay between CBET and the laser-smoothing schemes results in unique intensity variation across each beam profile, which is a systematic source of low-mode drive nonuniformity on OMEGA. Here, we model these effects and find that the predicted ℓ = 1 mode in the laser-absorption distribution is consistent with the systematic core-flow direction that has been determined from the OMEGA implosion database. We also observe good agreement with the measured core-flow directions for two specific sets of implosions (one with PS, the other without PS) when we also account for the measured beam mispointing and the beam power imbalance.

Crossed beam scattering↗

DOE OSTI2022

Non-Blind Deblurring for Fluorescence: A Deformable Latent Space Approach with Kernel Parameterization

We report N\non-blind deblurring (NBD) is a modeling method of the image deblurring problem in computer vision, where the blurring kernel is known or can be externally estimated. In this paper, we attempt to solve a parametric NBD problem, inspired by the simultaneous acquisition of ptychography and fluorescent imaging (FI). Ptychography is an imaging method that favors larger probes, i.e. convolutional kernels, while FI relies on a small probe for high resolution. Also, the kernel can be solved during ptychographic reconstruction. With Ptycho-FI using the same larger kernel, we can perform NBD on the blurred fluorescent images to achieve high-resolution FI, and thus speed up the experiments. To this end, we design a deep latent space deformation network that is directly parameterized by the kernel. The network consists of three components: encoder, deformer, and decoder, where the deformer is specifically meant to rectify the latent space representations of blurred images to a standard latent space, regardless of the kernel. The deformation network is trained with a two-stage training scheme. We conduct extensive experiments to confirm that our parametric model can adapt to drastically different blurring kernels and perform robust deblurring.

77 NANOSCIENCE AND NANOTECHNOLOGY↗