Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “multilevel 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 235 records · Page 13

Multilevel measurements of surface temperature over undulating terrain planted to barley

A ground and aircraft program was conducted to extend ground based methods for measuring soil moisture and crop water stress to aircraft and satellite altitudes. A 260ha agricultural field in California was used over the 1977-78 growing season. For cloud free days ground based temperature measurements over bare soil were related to soil moisture content. Water stress resulted from too much water, not from lack of it, as was expected. A theoretical examination of the canopy air temperature difference as affected by vapor pressure deficit and net radiation was developed. This analysis shows why surface temperatures delineate crop water stress under conditions of low humidity, but not under high humidity conditions. Multilevel temperatures acquired from the ground, low and high altitude aircraft, and the Heat Capacity Mapping Mission (HCMM) spacecraft were compared for two day and one night overpasses. The U-2 and low altitude temperatures were within 0.5 C. The HCMM data were analyzed using both the pre- and post-launch calibrations, with the former being considerably closer in agreement with the aircraft data than the latter.

Reginato, R. J.↗

Structural optimization by generalized, multilevel decomposition

The developments toward a general multilevel optimization capability and results for a three-level structural optimization are described. The method partitions a structure into a number of substructuring levels where each substructure corresponds to a subsystem in the general case of an engineering system. The method is illustrated by a portal framework that decomposes into individual beams. Each beam is a box that can be further decomposed into stiffened plates. Substructuring for this example spans three different levels: (1) the bottom level of finite elements representing the plates; (2) an intermediate level of beams treated as substructures; and (3) the top level for the assembled structure. The three-level case is now considered to be qualitatively complete.

Sobieszczanski-Sobieski, J.↗

Probabilistic Structural Analysis Methods (PSAM) for select space propulsion system components

The fourth year of technical developments on the Numerical Evaluation of Stochastic Structures Under Stress (NESSUS) system for Probabilistic Structural Analysis Methods is summarized. The effort focused on the continued expansion of the Probabilistic Finite Element Method (PFEM) code, the implementation of the Probabilistic Boundary Element Method (PBEM), and the implementation of the Probabilistic Approximate Methods (PAppM) code. The principal focus for the PFEM code is the addition of a multilevel structural dynamics capability. The strategy includes probabilistic loads, treatment of material, geometry uncertainty, and full probabilistic variables. Enhancements are included for the Fast Probability Integration (FPI) algorithms and the addition of Monte Carlo simulation as an alternate. Work on the expert system and boundary element developments continues. The enhanced capability in the computer codes is validated by applications to a turbine blade and to an oxidizer duct.

Source record↗

The solution of radiative transfer problems in molecular bands without the LTE assumption by accelerated lambda iteration methods

An iterative method based on the use of approximate transfer operators, which was designed initially to solve multilevel NLTE line formation problems in stellar atmospheres, is adapted and applied to the solution of the NLTE molecular band radiative transfer in planetary atmospheres. The matrices to be constructed and inverted are much smaller than those used in the traditional Curtis matrix technique, which makes possible the treatment of more realistic problems using relatively small computers. This technique converges much more rapidly than straightforward iteration between the transfer equation and the equations of statistical equilibrium. A test application of this new technique to the solution of NLTE radiative transfer problems for optically thick and thin bands (the 4.3 micron CO2 band in the Venusian atmosphere and the 4.7 and 2.3 micron CO bands in the earth's atmosphere) is described.

Kutepov, A. A.↗

Algebraic multigrid preconditioners for two-phase flow in porous media with phase transitions [Algebraic multigrid preconditioners for multiphase flow in porous media with phase transitions]

Multiphase flow is a critical process in a wide range of applications, including oil and gas recovery, carbon sequestration, and contaminant remediation. Numerical simulation of multiphase flow requires solving of a large, sparse linear system resulting from the discretization of the partial differential equations modeling the flow. In the case of multiphase multicomponent flow with miscible effect, this is a very challenging task. The problem becomes even more difficult if phase transitions are taken into account. A new approach to handle phase transitions is to formulate the system as a nonlinear complementarity problem (NCP). Unlike in the primary variable switching technique, the set of primary variables in this approach is fixed even when there is phase transition. Not only does this improve the robustness of the nonlinear solver, it opens up the possibility to use multigrid methods to solve the resulting linear system. The disadvantage of the complementarity approach, however, is that when a phase disappears, the linear system has the structure of a saddle point problem and becomes indefinite, and current algebraic multigrid (AMG) algorithms cannot be applied directly. In this study, we explore the effectiveness of a new multilevel strategy, based on the multigrid reduction technique, to deal with problems of this type. We demonstrate the effectiveness of the method through numerical results for the case of two-phase, two-component flow with phase appearance/disappearance. In conclusion, we also show that the strategy is efficient and scales optimally with problem size.

02 PETROLEUM↗

Bayesian Static Parameter Estimation for Partially Observed Diffusions via Multilevel Monte Carlo

In this article we consider static Bayesian parameter estimation for partially observed diffusions that are discretely observed. We work under the assumption that one must resort to discretizing the underlying diffusion process, for instance, using the Euler--Maruyama method. Given this assumption, we show how one can use Markov chain Monte Carlo (MCMC) and particularly particle MCMC [C. Andrieu, A. Doucet, and R. Holenstein, J. R. Stat. Soc. Ser. B Stat. Methodol., 72 (2010), 269--342] to implement a new approximation of the multilevel (ML) Monte Carlo (MC) collapsing sum identity. Our approach comprises constructing an approximate coupling of the posterior density of the joint distribution over parameter and hidden variables at two different discretization levels and then correcting by an importance sampling method. The variance of the weights are independent of the length of the observed data set. The utility of such a method is that, for a prescribed level of mean square error, the cost of this MLMC method is provably less than i.i.d. sampling from the posterior associated to the most precise discretization. However the method here comprises using only known and efficient simulation methodologies. Finally, the theoretical results are illustrated by inference of the parameters of two prototypical processes given noisy partial observations of the process: the first is an Ornstein--Uhlenbeck process and the second is a more general Langevin equation.

97 MATHEMATICS AND COMPUTING↗

Femtosecond Laser-Induced Nano-Joining of Volatile Tellurium Nanotube Memristor

Nanowire/nanotube memristor devices provide great potential for random-access high-density resistance storage. However, fabricating high-quality and stable memristors is still challenging. This paper reports multileveled resistance states of tellurium (Te) nanotube based on the clean-room free femtosecond laser nano-joining method. The temperature for the entire fabrication process was maintained below 190 °C. A femtosecond laser joining technique was used to form nanowire memristor units with enhanced properties. Femtosecond (fs) laser-irradiated silver-tellurium nanotube-silver structures resulted in plasmonic-enhanced optical joining with minimal local thermal effects. This produced a junction between the Te nanotube and the silver film substrate with enhanced electrical contacts. Noticeable changes in memristor behavior were observed after fs laser irradiation. Capacitor-coupled multilevel memristor behavior was observed. Compared to previous metal oxide nanowire-based memristors, the reported Te nanotube memristor system displayed a nearly two-order stronger current response. The research displays that the multileveled resistance state is rewritable with a negative bias.

42 ENGINEERING↗

A multilevel approach to the design of dc-free line codes

A multilevel approach to the design of dc-free line codes is presented. The different levels can be used for different purposes; for example, to control the maximum accumulated charge or to guarantee a certain minimum distance. The advantages of codes designed by this method over similar codes are the improved run-length/accumulated-charge parameters, higher transmission rate, and the systematic nature of the code construction. The multilevel structure allows the redundancy in the signal selection procedure to be allocated efficiently among the different levels. It also allows the use of suboptimal staged decoding procedures that have performance/complexity advantages over maximum-likelihood decoding.

Calderbank, A. R.↗

Multilevel Cybersecurity for Photovoltaic Systems

The motivation behind this project is to protect critical infrastructure in electric power generation pertaining to solar photovoltaic (PV) systems. This growing renewable energy resource is becoming a more vital part of the nation’s energy portfolio, particularly since it has achieved grid-parity to existing generation methods in terms of cost. It is thus vital that steps be taken to ensure the cybersecurity of these assets. The project goal was to devise a multilevel cybersecurity solution to address PV security gaps at the inverter and system levels, and field test the solution under the supervision and review of a US-based solar inverter manufacturer and PV installer/operator. A two-level cyberattack defense approach was formulated whereby the first level, the solar inverter level, hardens individual devices and achieves a deeply cyber-secure inverter. The inverter level security involves a multi-layer defense-in-depth approach for securing the inverter while also providing data for the system level algorithms. The second level, the system level, addresses intrusion detection and restoration involving an ensemble of inverters and relevant systems.

14 SOLAR ENERGY↗

Capacitor Voltage Ripple Reduction of A 4-Level Hybrid Clamped Converter Using Switching State Redundancy for Medium-Voltage Drive

Abstract—A four-level Hybrid Clamped Converter (4LHCC) is a novel promising DC/AC multilevel converter for high-efficiency medium-voltage drives. This paper presents a new control method of the 4L-HCC using the switching state redundancy to reduce the DC-link capacitor voltage ripples at low frequencies (<=30 Hz). The new switching state selection algorithm and flying capacitor balancing method are introduced and analyzed.Both of a3.6 kV 1MW 4L-HCC simulation model and a SiC-based 4L-HCC prototype is established to verify the proposed method. Simulation and experimental result shows that the maximum DC-link capacitor ripple can be reduced effectively from over 50% to less than 5% at low frequencies with the new method. Meanwhile, load voltage and currentdistortion caused by the excessive capacitor rippleisalso mitigatedsignificantly.

hybrid clamped converter, motor drive, voltage bal↗

The art of spacecraft design: A multidisciplinary challenge

Actual design turn-around time has become shorter due to the use of optimization techniques which have been introduced into the design process. It seems that what, how and when to use these optimization techniques may be the key factor for future aircraft engineering operations. Another important aspect of this technique is that complex physical phenomena can be modeled by a simple mathematical equation. The new powerful multilevel methodology reduces time-consuming analysis significantly while maintaining the coupling effects. This simultaneous analysis method stems from the implicit function theorem and system sensitivity derivatives of input variables. Use of the Taylor's series expansion and finite differencing technique for sensitivity derivatives in each discipline makes this approach unique for screening dominant variables from nondominant variables. In this study, the current Computational Fluid Dynamics (CFD) aerodynamic and sensitivity derivative/optimization techniques are applied for a simple cone-type forebody of a high-speed vehicle configuration to understand basic aerodynamic/structure interaction in a hypersonic flight condition.

Abdi, F.↗

Postoperative 3D spine reconstruction by navigating partitioning manifolds

Purpose: The postoperative evaluation of scoliosis patients undergoing corrective treatment is an important task to assess the strategy of the spinal surgery. Using accurate 3D geometric models of the patient’s spine is essential to measure longitudinal changes in the patient’s anatomy. On the other hand, reconstructing the spine in 3D from postoperative radiographs is a challenging problem due to the presence of instrumentation (metallic rods and screws) occluding vertebrae on the spine. Methods: This paper describes the reconstruction problem by searching for the optimal model within a manifold space of articulated spines learned from a training dataset of pathological cases who underwent surgery. The manifold structure is implemented based on a multilevel manifold ensemble to structure the data, incorporating connections between nodes within a single manifold, in addition to connections between different multilevel manifolds, representing subregions with similar characteristics. Results: The reconstruction pipeline was evaluated on x-ray datasets from both preoperative patients and patients with spinal surgery. By comparing the method to ground-truth models, a 3D reconstruction accuracy of 2.24 ± 0.90 mm was obtained from 30 postoperative scoliotic patients, while handling patients with highly deformed spines. Conclusions: This paper illustrates how this manifold model can accurately identify similar spine models by navigating in the low-dimensional space, as well as computing nonlinear charts within local neighborhoods of the embedded space during the testing phase. This technique allows postoperative follow-ups of spinal surgery using personalized 3D spine models and assess surgical strategies for spinal deformities.

61 RADIATION PROTECTION AND DOSIMETRY↗

Diffractive optics fabricated by direct write methods with an electron beam

State-of-the-art diffractive optics are fabricated using e-beam lithography and dry etching techniques to achieve multilevel phase elements with very high diffraction efficiencies. One of the major challenges encountered in fabricating diffractive optics is the small feature size (e.g. for diffractive lenses with small f-number). It is not only the e-beam system which dictates the feature size limitations, but also the alignment systems (mask aligner) and the materials (e-beam and photo resists). In order to allow diffractive optics to be used in new optoelectronic systems, it is necessary not only to fabricate elements with small feature sizes but also to do so in an economical fashion. Since price of a multilevel diffractive optical element is closely related to the e-beam writing time and the number of etching steps, we need to decrease the writing time and etching steps without affecting the quality of the element. To do this one has to utilize the full potentials of the e-beam writing system. In this paper, we will present three diffractive optics fabrication techniques which will reduce the number of process steps, the writing time, and the overall fabrication time for multilevel phase diffractive optics.

Kress, Bernard↗

Hydrothermal synthesis of a novel non-stoichiometric cubic bismuth phosphate (Bi{sub 12}(Bi{sub 4/5–5x}P{sub 5x})O{sub 19.2+5x}) crystal and the coexistence of multilevel Kink-Terrace-Step structures on Kossel (1 0 0) surface

- Highlights: • Pure cubic phase could be synthesized at high alkalinity. • The height of the lowest kinked [0 1 0] steps was about half of the unit cell. • Multilevel Kink-Terrace-Step structures coexisted on Kossel (1 0 0) surface. • Multilevel growth units are PO{sub 4}{sup −} incorporating into the channels of Bi–O frameworks. - Abstract: Nonstoichiometric cubic bismuth phosphate crystals were synthesized and grown by hydrothermal method. Under high alkalinity condition, three dimensional nuclei were generated and combined into cubic plates. The cubic bismuth phosphate crystals were combined by the corresponding cubic plates with various defects and rough surfaces. The cubic plates made (1 0 0) surface of the crystal cleavage face. Additionally, multilevel Kink-Terrace-Step structures on Kossel (1 0 0) surface were characterized through scanning electron microscopy and atomic force microscope techniques, respectively. It revealed that the lowest level Kink-Terrace-Step structure was in nano-scale, and the height of the elemental kinked [0 1 0] steps was about half of the unit cell. The elemental kinked [0 1 0] steps could bunch together and consequently constructed a higher level of Kink-Terrace-Step structures. As a result, multilevel Kink-Terrace-Step structures coexisted on the Kossel (1 0 0) surface of bismuth phosphate crystal, which were owing to the coexistence of multilevel growth units in crystal growth system.

36 MATERIALS SCIENCE↗

An Inrush Current Limiting Method for Grid-Connected Converters Considering Grid Voltage Disturbances

To support the electric power grid, some grid-connected converters are required to ride through abnormal grid conditions, including voltage disturbances. However, at the moment, when a grid voltage disturbance happens, a large inrush current could be induced due to the large difference between the grid and converter output pulsewidth modulation (PWM) voltages. This inrush current could be more severe in the faster switching SiC-based multilevel converter, which only needs a small filter inductance to meet its normal operation requirements, such as harmonics and peak ripple currents. In this article, a PWM mask method is used to limit the inrush current. With this method, the inrush current can be effectively limited to a preset value and can help the converter ride through grid voltage disturbances. Furthermore, the method is validated through experiments conducted on a single-phase SiC-based five-level converter considering low voltage, high voltage, and phase angle disturbances.

42 ENGINEERING↗

Pilot Point Optimization of Mining Boundaries for Lateritic Metal Deposits: Finding the Trade-off Between Dilution and Ore Loss

Geological contacts in lateritic metal deposits (footwall topographies) often delineate the orebody boundaries. Spatial variations seen in such contacts are frequently higher than those for the metal grades of the deposit. Therefore, borehole spacing chosen based on the grade variations cannot adequately capture the geological contact variability. Consequently, models created using such boreholes cause high volumetric uncertainties in the actual and targeted ore extraction volumes, which, in turn, lead to high unplanned dilution and ore losses. In this paper, a method to design optimum ore/mining boundaries for lateritic metal deposits is presented. The proposed approach minimizes the dilution/ore losses and comprises two main steps. First, the uncertainty on the orebody boundary is represented using a set of stochastic realizations generated with a multiple-point statistics algorithm. Then, the optimal orebody boundary is determined using an optimization technique inspired by a model calibration method called Pilot Points. The pilot points represent synthetic elevation values, and they are used to construct smooth mining boundaries using the multilevel B-spline technique. The performance of a generated surface is evaluated using the expected sum of losses in each of the stochastic realizations. The simulated annealing algorithm is used to iteratively determine the pilot point values which minimize the expected losses. The results show a significant reduction in the dilution volume as compared to those obtained from the actual mining operation.

58 GEOSCIENCES↗

TH-E-17A-02: High-Pitch and Sparse-View Helical 4D CT Via Iterative Image Reconstruction Method Based On Tensor Framelet

Purpose: 4D CT is routinely performed during radiation therapy treatment planning of thoracic and abdominal cancers. Compared with the cine mode, the helical mode is advantageous in temporal resolution. However, a low pitch (∼0.1) for 4D CT imaging is often required instead of the standard pitch (∼1) for static imaging, since standard image reconstruction based on analytic method requires the low-pitch scanning in order to satisfy the data sufficient condition when reconstructing each temporal frame individually. In comparison, the flexible iterative method enables the reconstruction of all temporal frames simultaneously, so that the image similarity among frames can be utilized to possibly perform high-pitch and sparse-view helical 4D CT imaging. The purpose of this work is to investigate such an exciting possibility for faster imaging with lower dose. Methods: A key for highpitch and sparse-view helical 4D CT imaging is the simultaneous reconstruction of all temporal frames using the prior that temporal frames are continuous along the temporal direction. In this work, such a prior is regularized through the sparsity transform based on spatiotemporal tensor framelet (TF) as a multilevel and high-order extension of total variation transform. Moreover, GPU-based fast parallel computing of X-ray transform and its adjoint together with split Bregman method is utilized for solving the 4D image reconstruction problem efficiently and accurately. Results: The simulation studies based on 4D NCAT phantoms were performed with various pitches (i.e., 0.1, 0.2, 0.5, and 1) and sparse views (i.e., 400 views per rotation instead of standard >2000 views per rotation), using 3D iterative individual reconstruction method based on 3D TF and 4D iterative simultaneous reconstruction method based on 4D TF respectively. Conclusion: The proposed TF-based simultaneous 4D image reconstruction method enables high-pitch and sparse-view helical 4D CT with lower dose and faster speed.

60 APPLIED LIFE SCIENCES↗

Application of the inverse fast multipole method as a preconditioner in a 3D Helmholtz boundary element method

Here we investigate an efficient preconditioning of iterative methods (such as GMRES) for solving dense linear systems Ax = b that follow from a boundary element method (BEM) for the 3D Helmholtz equation, focusing on the low-frequency regime. While matrix–vector products in GMRES can be accelerated through the low-frequency fast multipole method (LFFMM), the BEM often remains computationally expensive due to the large number of GMRES iterations. We propose the application of the inverse fast multipole method (IFMM) as a preconditioner to accelerate the convergence of GMRES. The IFMM is in essence an approximate direct solver that uses a multilevel hierarchical decomposition and low-rank approximations. The proposed IFMM-based preconditioning has a tunable parameter ε that balances the cost to construct a preconditioner M, which is an approximation of A -1 , and the cost to perform the iterative process by means of M. Namely, using a small (respectively, large) value of ε takes a long (respectively, short) time to construct M, while the number of iterations can be small (respectively, large). A comprehensive set of numerical examples involving various boundary value problems with complicated geometries and mixed boundary conditions is presented to validate the efficiency of the proposed method. We show that the IFMM preconditioner (with a nearly optimal ε of 10 -2 ) clearly outperforms some common preconditioners for the BEM, achieving 1.2–10.8 times speed-up of the computations, in particular when the scale of the underlying scatterer is about five wavelengths or more. In addition, the IFMM preconditioner is capable of solving complicated problems (in a reasonable amount of time) that BD preconditioner can not.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗