Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “SIMULATION METHODS”

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 847 records · Page 47

Solar panel thermal cycling testing by solar simulation and infrared radiation methods

For the solar panels of the European Space Agency (ESA) satellites OTS/MAROTS and ECS/MARECS the thermal cycling tests were performed by using solar simulation methods. The performance data of two different solar simulators used and the thermal test results are described. The solar simulation thermal cycling tests for the ECS/MARECS solar panels were carried out with the aid of a rotatable multipanel test rig by which simultaneous testing of three solar panels was possible. As an alternative thermal test method, the capability of an infrared radiation method was studied and infrared simulation tests for the ultralight panel and the INTELSAT 5 solar panels were performed. The setup and the characteristics of the infrared radiation unit using a quartz lamp array of approx. 15 sq and LN2-cooled shutter and the thermal test results are presented. The irradiation uniformity, the solar panel temperature distribution, temperature changing rates for both test methods are compared. Results indicate the infrared simulation is an effective solar panel thermal testing method.

Nuss, H. E.↗

Closing the Certification Gaps in Adaptive Flight Control Software

Over the last five decades, extensive research has been performed to design and develop adaptive control systems for aerospace systems and other applications where the capability to change controller behavior at different operating conditions is highly desirable. Although adaptive flight control has been partially implemented through the use of gain-scheduled control, truly adaptive control systems using learning algorithms and on-line system identification methods have not seen commercial deployment. The reason is that the certification process for adaptive flight control software for use in national air space has not yet been decided. The purpose of this paper is to examine the gaps between the state-of-the-art methodologies used to certify conventional (i.e., non-adaptive) flight control system software and what will likely to be needed to satisfy FAA airworthiness requirements. These gaps include the lack of a certification plan or process guide, the need to develop verification and validation tools and methodologies to analyze adaptive controller stability and convergence, as well as the development of metrics to evaluate adaptive controller performance at off-nominal flight conditions. This paper presents the major certification gap areas, a description of the current state of the verification methodologies, and what further research efforts will likely be needed to close the gaps remaining in current certification practices. It is envisioned that closing the gap will require certain advances in simulation methods, comprehensive methods to determine learning algorithm stability and convergence rates, the development of performance metrics for adaptive controllers, the application of formal software assurance methods, the application of on-line software monitoring tools for adaptive controller health assessment, and the development of a certification case for adaptive system safety of flight.

Jacklin, Stephen A.↗

Analysis of the static magnetic field-dependent optical transmission of Ni nanorod colloidal suspensions

The magnetic field-dependent optical transmission of dilute Ni nanorod aqueous suspensions was investigated. A series of four samples of nanorods were synthesized using the AAO template method and processed to stable colloids. The distributions of their length and diameter were characterized by analysis of TEM images and revealed average diameters of ∼25 nm and different lengths in the range of 60 nm–1100 nm. The collinear magnetic and optical anisotropy was studied by static field-dependent transmission measurements of linearly polarized light parallel and perpendicular to the magnetic field direction. The experimental results were modelled assuming the field-dependent orientation distribution function of a superparamagnetic ensemble for the uniaxial ferromagnetic nanorods in liquid dispersion and extinction cross sections for longitudinal and transversal optical polarization derived from different approaches, including the electrostatic approximation and the separation of variables method, both applied to spheroidal particles, as well as finite element method simulations of spheroids and capped cylindrical particles. The extinction cross sections were compared to reveal the differences associated with the approximations of homogeneous polarization and/or particle shape. The consequences of these approximations for the quantitative analysis of magnetic field-dependent optical transmission measurements were investigated and a reliable protocol derived. Furthermore, the changes in optical cross sections induced by electromagnetic interaction between two nanorods in parallel end-to-end and side-by-side configuration as a function of their separation were studied.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

PRESSURE SHIFT AND GRAVITATIONAL REDSHIFT OF BALMER LINES IN WHITE DWARFS: REDISCUSSION

The Stark-induced shift and asymmetry, the so-called pressure shift (PS) of H{sub α} and H{sub β} Balmer lines in spectra of DA white dwarfs (WDs), have been examined in detail as masking effects in measurements of the gravitational redshift in WDs. The results are compared with our earlier ones from a quarter of a century ago. In these earlier papers, the standard, symmetrical Stark line profiles, as a dominant constituent of the Balmer line profiles but shifted as a whole by the PS effect, were applied to all spectrally active layers of the WD atmosphere. At present, in each of the WD layers, the Stark line profiles (especially of H{sub β}) are inherently asymmetrical and shifted due to the effects of strong inhomogeneity of the perturbing fields in plasma. To calculate the Stark line profiles in successive layers of the WD atmosphere we used the modified Full Computer Simulation Method, able to take adequately into account the complexity of local elementary quantum processes in plasma. In the case of the H{sub α} line, the present value of Stark-induced shift of the synthetic H{sub α} line profile is about half the previous one and it is negligible in comparison with the gravitational redshift. In the case of the H{sub β} line, the present value of Stark-induced shift of the synthetic H{sub β} line profile is about twice the previous one. The source of this extra shift is the asymmetry of H{sub β} peaks.

79 ASTRONOMY AND ASTROPHYSICS↗

Ultrasonic fingerprint sensor using a piezoelectric micromachined ultrasonic transducer array integrated with complementary metal oxide semiconductor electronics

This paper presents an ultrasonic fingerprint sensor based on a 24 × 8 array of 22 MHz piezoelectric micromachined ultrasonic transducers (PMUTs) with 100 μm pitch, fully integrated with 180 nm complementary metal oxide semiconductor (CMOS) circuitry through eutectic wafer bonding. Each PMUT is directly bonded to a dedicated CMOS receive amplifier, minimizing electrical parasitics and eliminating the need for through-silicon vias. The array frequency response and vibration mode-shape were characterized using laser Doppler vibrometry and verified via finite element method simulation. The array's acoustic output was measured using a hydrophone to be ∼14 kPa with a 28 V input, in reasonable agreement with predication from analytical calculation. Pulse-echo imaging of a 1D steel grating is demonstrated using electronic scanning of a 20 × 8 sub-array, resulting in 300 mV maximum received amplitude and 5:1 contrast ratio. Because the small size of this array limits the maximum image size, mechanical scanning was used to image a 2D polydimethylsiloxane fingerprint phantom (10 mm × 8 mm) at a 1.2 mm distance from the array.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Nonlinear absorption of short intense laser pulse in multispecies plasma

In the present paper, the detailed investigation concerning the effect of inclusion of heavy negative ions into the finite background plasma on the laser absorption has been carried out by employing particle-in-cell simulation method. For this purpose, in this configuration, the laser energy absorption relying on the nonlinear phenomena such as phase-mixing, wave-breaking, and scattering has been studied in the Raman-Brillouin regime. It is shown that the inclusion of heavy negative ions suppresses the scattering while increases the phase-mixing time. Moreover, it is illustrated that this inclusion can increase the laser absorption in finite plasma environment, after saturation. The obtained results are expected to be relevant to the experiments on the mass spectrometry with laser desorption techniques as well as on the laser-plasma interaction with application to particles acceleration.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Numerical simulation of amplification of frequency-modulated radiation in a gas amplifier THL-100 laser system

Numerical simulation methods have been used to study the effect of frequency modulation of input laser radiation with a pulse duration of 250 ps on the gain characteristics in a gas XeF (C – A) amplifier of the hybrid THL-100 laser system at a pump energy of 270 – 400 J. Using the Wigner distribution function, local radiation spectra are calculated at various points of the input laser beam with positive and negative chirps. The effect of the pump energy of the amplifier on the energy, intensity, and duration of the output laser radiation pulses is investigated. It is shown that with an insignificant change in energy, the pulse duration of the output radiation turns out to be longer with a negative chirp of the input radiation. (lasers)

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Molecular-based analysis of nanoparticle solvation: Classical density functional approach

Proper statistical mechanics understanding of nanoparticle solvation processes requires an accurate description of the molecular structure of the solvent. Achieving this goal with standard molecular dynamics (MD) simulation methods is challenging due to large length scales. An alternative approach to this problem can be formulated using classical density functional theory (cDFT), where a full configurational description of the positions of all the atoms is replaced by collective atomic site densities in the molecule. Using an example of the negatively charged silica-like system in an aqueous polar environment represented by a two-site water model, we demonstrate here that cDFT can reproduce MD data at a fraction of the computational cost. An important implication of this result is the ability to understand how the solvent molecular features may affect the system’s properties at the macroscopic scale. A concrete example highlighted in this work is the analysis of nanoparticle interactions with sizes of up to 100 nm in diameter.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The complex chemical Langevin equation

The chemical Langevin equation (CLE) is a popular simulation method to probe the stochastic dynamics of chemical systems. The CLE’s main disadvantage is its break down in finite time due to the problem of evaluating square roots of negative quantities whenever the molecule numbers become sufficiently small. We show that this issue is not a numerical integration problem, rather in many systems it is intrinsic to all representations of the CLE. Various methods of correcting the CLE have been proposed which avoid its break down. We show that these methods introduce undesirable artefacts in the CLE’s predictions. In particular, for unimolecular systems, these correction methods lead to CLE predictions for the mean concentrations and variance of fluctuations which disagree with those of the chemical master equation. We show that, by extending the domain of the CLE to complex space, break down is eliminated, and the CLE’s accuracy for unimolecular systems is restored. Although the molecule numbers are generally complex, we show that the “complex CLE” predicts real-valued quantities for the mean concentrations, the moments of intrinsic noise, power spectra, and first passage times, hence admitting a physical interpretation. It is also shown to provide a more accurate approximation of the chemical master equation of simple biochemical circuits involving bimolecular reactions than the various corrected forms of the real-valued CLE, the linear-noise approximation and a commonly used two moment-closure approximation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A molecular dynamics study of chloride binding by the cryptand SC24

The capture of chloride from water by the tetraprotonated form of the spherical macrotricyclic molecule SC24 was studied using molecular dynamics simulation methods. This model ionophore represents a broad class of molecules which remove ions from water. Two binding sites for the chloride were found, one inside and one outside the ligand. These sites are separated by a potential energy barrier of approximately 20 kcal mol-1. The major contribution to this barrier comes from dehydration of the chloride. The large, unfavorable dehydration effect is compensated for by an increase in electrostatic attraction between the oppositely charged chloride and cryptand, and by energetically favorable rearrangements of water structure. Additional assistance in crossing the barrier and completing the dehydration of the ion is provided by the shift of three positively charged hydrogen atoms of the cryptand towards the chloride. This structural rigidity is partially responsible for its selectivity.

NASA Discipline Exobiology↗

Laplace Transform–Based Quantum Eigenvalue Transformation via Linear Combination of Hamiltonian Simulation

Eigenvalue transformations, which include solving time-dependent differential equations as a special case, have a wide range of applications in scientific and engineering computation. While quantum algorithms for singular value transformations are well studied, eigenvalue transformations are distinct, especially for nonnormal matrices. Here, we propose an efficient quantum algorithm for performing a class of eigenvalue transformations that can be expressed as a certain type of matrix Laplace transformation. This allows us to significantly extend the recently developed linear combination of Hamiltonian simulation method [D. An, J.-P. Liu, and L. Lin, Phys. Rev. Lett., 131 (2023), 150603; D. An, A. M. Childs, and L. Lin, Commun. Math. Phys. 407, 19 (2026)] to represent a wider class of eigenvalue transformations, such as powers of the matrix inverse, 𝐴 −𝑘 , and the exponential of the matrix inverse, 𝑒 −𝐴 −1 . The latter can be interpreted as the solution of a mass-matrix differential equation of the form form 𝐴⁢𝑢′⁡⁡(𝑡) =−𝑢⁡(𝑡). We demonstrate that our eigenvalue transformation approach can solve this problem without explicitly inverting 𝐴, thereby reducing the computational complexity.

Laplace transform↗

Mapping and Characterizing the Visual Impacts of the Existing US Wind Turbine Fleet

Visual impacts of wind turbines have been a persistent concern for wind energy development in the United States (US) for decades and remain a major source of project delays and cancellations. Assessments of visual impacts are frequently performed at a local scale for individual projects, but a comprehensive understanding of broader geographic patterns in visual impacts across the US is lacking. This paper presents a visual impact assessment of the existing land-based wind turbine fleet of the contiguous United States (CONUS). The assessment combines geographic information systems and 3D simulation methods to account for key factors driving the visual magnitude of impacts from the installed turbines. The results indicate that, despite the deployment of approximately 70,000 turbines and over 144 gigawatts of land-based wind in the CONUS, the visual impacts are very small when measured as a proportion of land area, population, and sensitive visual resources. Nonetheless, visual impacts are not evenly distributed: people experience a concentrated share in a small number of natural settings, primarily including plains, prairies, and deserts. Finally, we find that although increased density of wind development consistently leads to visual impacts across a greater proportion of land, it does not always lead to impacts to a greater share of the population. These findings suggest that visual impacts from wind energy are generally well-mitigated across the CONUS to date but also highlight the need for a deeper understanding of landscape sensitivity and individual perceptions of wind turbines in the most heavily impacted natural settings.

3D simulation↗

Structural Tunability of Multicompartment Micelles as a Function of Lipophilic–Fluorophilic Block Length Ratio

Structural variation in multicompartment micelles consisting of lipophilic–hydrophilic–fluorophilic (hereafter referred to as BAC) triblock copolymers is investigated using the dissipative particle dynamics simulation method. It is demonstrated from our results that the structure of BAC multicompartment micelles is effectively tuned as a function of the lipophilic–fluorophilic ratio parameter, here termed $R_1$, of the constituent linear triblock copolymers. In particular, a morphological deviation from onion-like ABC micelles arises in BAC micelle systems as $R_1$ increases. The morphologies of BAC micelles with $R_1\ll$1 or $R_1\gg$1 display striking similarities, with the only notable difference being an inversion of the lipophilic and fluorophilic regions. When $R\approx$ 1, segmented worm-like structures with multiple cores are favored in BAC micelle systems. Through this study, it is confirmed that the block length ratio is an effective control parameter to tune the structure of multicompartment micelles.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A particle-in-cell method for the simulation of plasmas based on an unconditionally stable field solver

Here, we propose a new particle-in-cell (PIC) method for the simulation of plasmas based on a recently developed, unconditionally stable solver for the wave equation. This method is not subject to a CFL restriction, limiting the ratio of the time step size to the spatial step size, typical of explicit methods, while maintaining computational cost and code complexity comparable to such explicit schemes. We describe the implementation in one and two dimensions for both electrostatic and electromagnetic cases, and present the results of several standard test problems, showing good agreement with theory with time step sizes much larger than allowed by typical CFL restrictions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Silica Surface States and their Wetting Characteristics

The wetting characteristics of silica (SiO 2 ) surfaces can be described by molecular dynamics (MD) and ab initio simulations, including comparison of silica surfaces (talc (001), siloxanated quartz, tridymite (001) and quartz (001)), some of which have not been considered previously in the literature. Classic MD and ab initio simulation methods have been used to determine the contact angle, interfacial water structure, hydroxylation reaction and hydration energy, the results of which are compared with experimental results reported in the literature. Additionally, wetting of silica surfaces depends on surface polarity and extent of hydroxylation. The non-polar siloxane surfaces are characterized by a contact angle of about 80°, an MD ‘water exclusion zone’ of about 3 Å, a relaxed interfacial water orientation, inertness to hydroxylation and minimal hydration energy. The polar silica surfaces can be wetted by water and have a more ordered interfacial water structure. Silanol groups form at the polar silica surface during hydroxylation reactions, and the calculated hydration energy of -1·2 and -1·6 eV matches the experimental heat of immersion measurements reported in the literature, which correspond to hydrogen (H) bonding with interfacial water. In conclusion, fundamental understanding of silica surfaces is important for understanding flotation phenomena and fluid flow in silica nanopores.

36 MATERIALS SCIENCE↗

CRPropa 3.1—a low energy extension based on stochastic differential equations

The propagation of charged cosmic rays through the Galactic environment influences all aspects of the observation at Earth. Energy spectrum, composition and arrival directions are changed due to deflections in magnetic fields and interactions with the interstellar medium. Today the transport is simulated with different simulation methods either based on the solution of a transport equation (multi-particle picture) or a solution of an equation of motion (single-particle picture). We developed a new module for the publicly available propagation software CRPropa 3.1, where we implemented an algorithm to solve the transport equation using stochastic differential equations. This technique allows us to use a diffusion tensor which is anisotropic with respect to an arbitrary magnetic background field. The source code of CRPropa is written in C++ with python steering via SWIG which makes it easy to use and computationally fast. In this paper, we present the new low-energy propagation code together with validation procedures that are developed to proof the accuracy of the new implementation. Furthermore, we show first examples of the cosmic ray density evolution, which depends strongly on the ratio of the parallel κ{sub ∥} and perpendicular κ{sub ⊥} diffusion coefficients. This dependency is systematically examined as well the influence of the particle rigidity on the diffusion process.

79 ASTRONOMY AND ASTROPHYSICS↗

Probing the Dependence of the Intergalactic Medium on Large-scale Environment Using the Low-redshift Lyα Forest

We examine the statistics of the low-redshift Lyα forest in a cosmological hydrodynamic simulation using adaptive mesh refinement that has sufficient volume to include distinct large-scale environments. We compare our H i column density distribution of absorbers both with recent work and between two highly refined regions of our simulation: a large-scale overdensity and a large-scale underdensity (on scales of approximately 20 Mpc). We recover the average results presented by Kollmeier et al. using different simulation methods. We further break down these results as a function of environment to examine the detailed dependence of absorber statistics on large-scale density. We find that the slope of the H i column density distribution (CDD) in the range 10{sup 12.5} ≤ N{sub HI}/cm{sup −2} ≤ 10{sup 14.5} depends on environment such that the slope becomes steeper for higher environmental density, and this difference reflects distinct physical conditions of the intergalactic medium (IGM) on these scales. We track this difference to the different temperature structures of filaments in varying environments. Specifically, filaments in the overdensity are hotter, and correspondingly are composed of gas with lower H i fractions than those in underdense environments. Our results highlight that in order to understand the physics driving the H i CDD, we need not only improved accounting of the sources of ionizing UV photons, but also a clearer picture of the physical conditions of the IGM and how it may vary as a function of large-scale environment.

79 ASTRONOMY AND ASTROPHYSICS↗

Compensation of heat load deformations using adaptive optics for the ALS upgrade: a wave optics study

In this paper, a realistic wave optics simulation method has been developed to study how wavefront distortions originating from heat load deformations can be corrected using adaptive X-ray optics. Several planned soft X-ray and tender X-ray insertion-device beamlines in the Advanced Light Source upgrade rely on a common design principle. A flat, first mirror intercepts the white beam; vertical focusing is provided by a variable-line-space monochromator; and horizontal focusing comes from a single, pre-figured, adaptive mirror. A variety of scenarios to cope with thermal distortion in the first mirror are studied by finite-element analysis. The degradation of the intensity distribution at the focal plane is analyzed and the adaptive optics that correct it is modeled. The range of correctable wavefront errors across the operating range of the beamlines is reported in terms of mirror curvature and spatial frequencies. The software developed is a one-dimensional wavefront propagation package made available in the OASYS suite, an adaptable, customizable and efficient beamline modeling platform.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗