SEARCH · Engineering Papers
Results for “analytical model”
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.
Oscillation Risks of Grid-Following and Grid-Forming Inverter-Based Resources in Series-Compensated Networks
Here, this paper investigates the dynamic behavior of a grid-connected inverter-based resource (IBR) when connected radially to a series compensated line. Potential interactions between the series compensation and the IBR have been identified for both types: grid-following (GFL) or grid-forming (GFM). The study begins with electromagnetic transient (EMT) simulations to demonstrate stability issues. Subsequently, nonlinear analytical models are formulated in the dq frame, validated against the EMT simulation, and leveraged to assess eigenvalues and participation factors. Influencing factors of the dominant oscillation modes have been identified. The analysis results show that series compensation may make a mode associated with the synchronization unit unstable. Furthermore, customized feedback systems are built for the synchronizing loop. Series compensation can increase the sensitivity of the voltage phase angle towards the synchronizing angle, and introduce phase lag in the real power response towards the synchronizing angle. These factors may cause interactions with the phase-locked loop in GFL-IBR systems and with power-based synchronization in GFM-IBR systems, potentially leading to instability.
Failure Probability Constrained AC Optimal Power Flow
Despite cascading failures being the central cause of blackouts in power transmission systems, existing operational and planning decisions are made largely by ignoring their underlying cascade potential. This paper posits a reliability-aware AC Optimal Power Flow formulation that seeks to design a dispatch point which has a low operator-specified likelihood of triggering a cascade starting from any single component outage. By exploiting a recently developed analytical model of the probability of component failure, our Failure Probability-constrained ACOPF (FP-ACOPF) utilizes the system's expected first failure time as a smoothly tunable and interpretable signature of cascade risk. Here, we use techniques from bilevel optimization and numerical linear algebra to efficiently formulate and solve the FP-ACOPF using off-the-shelf solvers. Extensive simulations on the IEEE 118-bus case show that, when compared to the unconstrained and N-1 security-constrained ACOPF, our probability-constrained dispatch points can significantly lower the probabilities of long severe cascades and of large demand losses, while incurring only minor increases in total generation costs.
Sampling in Long-Screened Wells: Issues, Misconceptions, and Solutions
The issues associated with long-screened wells (LSWs) (and open boreholes) at contaminated sites are well documented in the groundwater literature but are still not fully appreciated in practice. As established in seminal and review papers going back over three decades, the interpretation of sampling results from LSWs is challenging in the presence of vertical hydraulic gradients and borehole flow; furthermore, LSWs allow for vertical redistribution of contamination between aquifer layers. Acknowledgment of these issues has led to the development of new technologies and well designs to enable discrete-zone monitoring (DZM), yet LSWs remain common for many reasons, for example, as multipurpose wells, for geophysical logging, and (or) as legacy installations. Despite the literature on LSWs and despite the adoption of DZM at many sites, the use of LSWs persists and the challenges of interpreting sampling results from LSWs remain. In this issue paper, we provide a conceptual overview of the problems posed by LSWs and review existing literature and past work to improve the interpretation of sampling in LSWs. We draw on experience from previous studies at the Hanford Site in eastern WA, USA, and use synthetic examples to illustrate key concepts and challenges for interpretation. A recently published analytical modeling framework is used to develop illustrative synthetic examples and demonstrate a workflow for building scientific intuition to understand issues around interpreting samples from LSWs, which is critical to effective characterization and groundwater remediation at sites with LSWs.
Shear wave generation by mode conversion in picosecond ultrasonics: Impact of grain orientation and material properties
Abstract Picosecond ultrasonics has seen wide application in the investigation of interfaces in semiconductors and 3D imaging in biological systems. Shear ultrasonic waves are important for elastic constant and grain orientation measurement in ceramic materials. In this study, we investigate the impact of grain orientation and material's elastic properties on the generation efficiency of shear waves by mode conversion at a transducer thin film/substrate interface. The solution of acoustic wave equations suggests that crystal grain orientation has a strong impact on the generated shear wave amplitude. This dependence is found to be closely related to the magnitude of longitudinal component of wave displacements. The applicability of analytical model is validated by the experimental results from time‐domain Brillouin scattering. Moreover, material properties determine acoustic wave amplitudes based on the acoustic mismatch model and, particularly, large elastic anisotropy defined by Zener ratio favors strong shear waves. In the light of this analysis, several recommendations on suitable grain orientations and film materials are made to facilitate shear wave detection.
Verification, Validation, and Uncertainty Quantification in Thermal Hydraulics, Freeman Scholar Lecture (2019)
Engineering problems are generally solved by analytical models or computer codes. These models, in addition to conservation equations, also include many empirical relationships and approximate numerical methods. Each of these components contributes to the uncertainty in the prediction. A systematic approach to judge the applicability of the code to the intended application is needed. It starts from verification of implementation of formulation in the code, identification of important phenomena, finding relevant tests with quantified uncertainty for these phenomena, and validation of the code by comparing predictions with the relevant test data. The relevant tests must address phenomena as expected in the intended application. In case of small size or limited condition tests, the scaling analyses are needed to assess the relevancy of the tests. Finally, a statement of uncertainty in the prediction is needed. Systematic approaches are described to aggregate uncertainties from different components of the code for intended application. Here, verification, validation, and uncertainty quantifications (VVUQs) are briefly described.
Attached Jet Velocity Profiles in Mixing Tanks
Free jets have been studied in detail over much of the last century, but the theory for offset and attached jets remains incomplete. Attached jets differ from free jets in that they lose momentum to nearby surfaces, attenuating their velocities. The velocity profiles of free circular jets are nearly Gaussian, with quantitative mathematical descriptions derived from first principles by Goertler and Tollmien (Rajaratnam, 1976). In contrast, mathematical descriptions of three-dimensional attached jets from circular nozzles remain much less mature. Agelin-Chaab and Tachie (2011) used particle imaging velocimetry of a three-dimensional attached jet to show that the scaled velocity decays with scaled distance from the nozzle with a power law exponent between -1.15 and -1.21, which is larger in magnitude than that of a free jet. However, quantitative analytical expressions for the velocity profiles of attached jets similar to those of free jets remain elusive. This paper addresses this critical gap. Here we evaluate the velocity profiles of three-dimensional offset jets emerging from circular nozzles that that become attached jets. These jets lose momentum due to interactions with nearby surfaces and are important to evaluating flows in mixing vessels and to suspending solids and trapped gases in radioactive waste tanks. Despite the importance of attached jets, prior insight has been purely experimental, limited to overly simplistic analytical models, or restricted to computationally expensive computational fluid dynamics case studies. We compare the expression of Verhoff (1963) to experimental results to find reasonable quantitative agreement. As stated by Agelin-Chaab and Tachie (2011), “detailed velocity measurements of 3D offset jets are rare.” Such remains the case. This study adds to the literature by providing information at two additional Reynolds numbers (1.43 106 and 1.87 106) and evaluating simple but accurate expressions for velocity profiles. These Reynolds numbers are an order of magnitude higher than prior reports. The semi-empirical stream wise velocity profile perpendicular to the surface proposed by Verhoff (1963) is in approximate agreement with these velocity profiles, which is surprising because these attached jets are three-dimensional instead of two-dimensional as evaluated by Verhoff. However, additional work is necessary to fully describe these profiles quantitatively.
Performance of plane wavefront Fizeau interferometers in power spectral density measurements with tilted plane optics
Binary pseudo-random array (BPRA) test samples are useful devices for calibrating the instrument transfer function (ITF) of Fizeau interferometers, interferometric microscopes, and other optical and non-optical surface and wavefront metrology instruments. The intrinsic white noise character of the power spectral density (PSD) function of the BPRA pattern simplifies the extraction of the ITF from the measured PSD. The ITF determined in a dedicated calibration experiment can be used to reconstruct the surface height profile from the measured data, effectively enhancing the instrument’s spatial resolution. For a high confidence reconstruction procedure, a reliable analytical model of the IFT is desirable. Usually, the model accounts for the contributions to the ITF related the imperfections of the instrument’s optical and detector systems. Here, we experimentally demonstrate that in the case of surface height metrology with Fizeau interferometers, the PSD measurements and, therefore, the efficacy of the ITF calibration of the tool, are strongly affected by the instrument data acquisition and processing procedures, as well as by the shape of the optic under test and its alignment with respect to the interferometer.
Extension of the bright high-harmonic photon energy range via nonadiabatic critical phase matching
The concept of critical ionization fraction has been essential for high-harmonic generation, because it dictates the maximum driving laser intensity while preserving the phase matching of harmonics. In this work, we reveal a second, nonadiabatic critical ionization fraction, which substantially extends the phase-matched harmonic energy, arising because of the strong reshaping of the intense laser field in a gas plasma. We validate this understanding through a systematic comparison between experiment and theory for a wide range of laser conditions. In particular, the properties of the high-harmonic spectrum versus the laser intensity undergoes three distinctive scenarios: (i) coincidence with the single-atom cutoff, (ii) strong spectral extension, and (iii) spectral energy saturation. We present an analytical model that predicts the spectral extension and reveals the increasing importance of the nonadiabatic effects for mid-infrared lasers. These findings are important for the development of high-brightness soft x-ray sources for applications in spectroscopy and imaging.
Knots are not for naught: Design, properties, and topology of hierarchical intertwined microarchitected materials
Lightweight and tough engineered materials are often designed with three-dimensional hierarchy and interconnected structural members whose junctions are detrimental to their performance because they serve as stress concentrations for damage accumulation and lower mechanical resilience. We introduce a previously unexplored class of architected materials, whose components are interwoven and contain no junctions, and incorporate micro-knots as building blocks within these hierarchical networks. Tensile experiments, which show close quantitative agreements with an analytical model for overhand knots, reveal that knot topology allows a new regime of deformation capable of shape retention, leading to a ~92% increase in absorbed energy and an up to ~107% increase in failure strain compared to woven structures, along with an up to ~11% increase in specific energy density compared to topologically similar monolithic lattices. Our exploration unlocks knotting and frictional contact to create highly extensible low-density materials with tunable shape reconfiguration and energy absorption capabilities.
Calorimetric wire detector for measurement of atomic hydrogen beams
A calorimetric detector for minimally disruptive measurements of atomic hydrogen beams is described. The calorimeter measures heat released by the recombination of hydrogen atoms into molecules on a thin wire. As a demonstration, the angular distribution of a beam with a peak intensity of ≈ 10 16 atoms/(cm 2 s) is measured by translating the wire across the beam. The data agree well with an analytic model of the beam from the thermal hydrogen atom source. Using the beam shape model, the relative intensity of the beam can be determined to 5% precision or better at any angle.
Hyperdiffusion of dust particles in a turbulent tokamak plasma
The effect of plasma turbulence on the trajectories of dust particles is investigated for the first time. The dynamics of dust particles is computed using the ad-hoc developed Dust Injection Simulator code, using a 3D turbulent plasma background computed with the TOKAM3X code. As a result, the evolution of the particle trajectories is governed by the ion drag force, and the shape of the trajectory is set by the Stokes number $St\propto a_d/n_0$, with $a_d$ the dust radius and $n_0$ the density at the separatrix. The plasma turbulence is observed to scatter the dust particles, exhibiting a hyperdiffusive regime in all cases. The amplitude of the turbulent spread of the trajectories $\Delta r^2$ is shown to depend on the ratio $Ku/St$, with $Ku\propto u_{rms}$ the Kubo number and $u_{rms}$ the fluctuation level of the plasma flow. These results are compared with a simple analytical model, predicting $\Delta r^2\propto (Ku/St)^2t^3$, or $\Delta r^2\propto (u_{rms}n_0/a_d)^2t^3$. As the dust is heated by the plasma fluxes, thermionic emission sets the dust charge, originally negative, to slightly positive values. This results in a substantial reduction of the ion drag force through the suppression of its Coulomb scattering component. The dust grain inertia is then no longer negligible, and drives the transition from a hyperdiffusive regime towards a ballistic one.
Analytical Investigation of the Role of Lateral Mixing in the Evolution of Nonprecipitating Cumulus. Part I: Developing Clouds
Evolution of nonprecipitating cumulus clouds (Cu) at the developing stage under the influence of lateral entrainment and mixing is studied analytically using a minimalistic analytical model. We present a model of an ascending cloud volume (a model of developing Cu) whose structure is determined by the processes of droplet diffusion growth/evaporation and entrainment mixing in the horizontal direction. Spatial and time changes of liquid water content, the adiabatic fraction, droplet concentration, and the mean volume droplet radius are calculated. It is shown that the existence of a nondiluted core in a growing cumulus cloud significantly depends on the cloud width and vertical velocity. While at the updraft velocity of 2 m s −1 the core of a 400-m-wide cloud becomes diluted at distances of a few hundred meters above cloud base, the core of a cloud of 1000-m width remains nondiluted at distances up to 1500 m above cloud base. The explanation of this result is simple: the increase in cloud width and the decrease in the updraft velocity increase the time during which the cloud is diluted due to mixing. Since lateral mixing synchronously decreases both the cloud water content and droplet concentration, the variation of the mean volume droplet radius is low inside the cloud. The approximate quantitative condition for cloud formation in updraft is derived. It is shown that a cloud can arise when its vertical velocity exceeds a critical value. To produce clouds, narrow turbulent plumes should ascend at higher velocity as compared to wider plumes. High humidity of the environment air is favorable for formation of clouds from plumes. The comparison of the obtained results with previously published observational data indicates a reasonable agreement. The results can be useful for parameterization purposes.
Analytical Investigation of the Role of Lateral Mixing in the Evolution of Nonprecipitating Cumulus. Part II: Dissolving Stage
A minimalistic analytical model allowing analysis of the dissolving stage of nonprecipitating convective clouds is proposed. The model takes into account two mechanisms: turbulent mixing with a dry environment and cloud volume settling. The temporal changes in the spatial structure of a cloud and in its immediate environment in the course of cloud dissolving are analyzed. The comparison of the effects of a temperature increase in the course of cloud descent and mixing with dry surrounding air shows that the descent is a dominating factor determining a decrease in the liquid water content (LWC), while mixing has a stronger effect on the cloud shape. Narrowing/broadening of clouds due to lateral mixing with dry air during cloud dissolving is determined by the potential evaporation parameter proportional to the ratio of the saturation deficit in the cloud environment to LWC. An equation for cloud dissolving time is obtained. After a cloud totally dissolves, it leaves behind an area with humidity exceeding that of the environment. The main parameter determining the dissolving time is the downdraft velocity. It should exceed 50 cm s -1 in order to provide reasonable dissolving time. The turbulent intensity, LWC, and humidity of the environment air also have an impact on dissolving time: the lower the LWC and the humidity of environment air, the faster cloud dissolving is. The simple solution presented in this paper can be useful for evaluation of cloud characteristics at the dissolving stage and can be included in procedures of parameterization of cloud cover formed by nonprecipitating or slightly precipitating cumulus clouds (Cu). Values of the environment humidity and temperature, LWC at cloud top, cloud width, vertical velocity of downdraft, and the turbulent coefficient should be parameters of this parameterization.
Theoretical Understanding of the Linear Relationship between Convective Updrafts and Cloud-Base Height for Shallow Cumulus Clouds. Part II: Continental Conditions
This is the Part II of a two-part study that seeks a theoretical understanding of an empirical relationship for shallow cumulus clouds: subcloud updraft velocity covaries linearly with the cloud-base height. This work focuses on continental cumulus clouds that are more strongly forced by surface fluxes and more deviated from equilibrium than those over oceans (Part I). We use a simple analytical model for shallow cumulus that is well tested against a high-resolution (25 m in the horizontal) large-eddy simulation model. Consistent with a conventional idea, we find that surface Bowen ratio is the key variable that regulates the covariability of both parameters: under the same solar insolation, a drier surface allows for stronger buoyancy flux, triggering stronger convection that deepens the subcloud layer. We find that the slope of the Bowen-ratio-regulated relationship between the two parameters (defined as λ) is dependent on both the local time and the stability of the lower free atmosphere. The value of λ decreases with time exponentially from sunrise to early afternoon and linearly from early afternoon to sunset. The value of λ is larger in a more stable atmosphere. In addition, continental λ in the early afternoon more than doubles the oceanic λ. Validation of the theoretical results against ground observations over the Southern Great Plains shows a reasonable agreement. Physical mechanisms underlying the findings are explained from the perspective of different time scales at which updrafts and cloud-base height respond to a surface flux forcing.
Pattern formation on free surfaces via plastic buckling and periodic folding
Spontaneous formation of surface patterns such as folds and creases is widely observed in engineered systems, such as thin films and machined surfaces, as well as in natural systems such as floral petals and leaves. A recognized route for their formation is via surface elastic instabilities. Here we demonstrate an alternate route for their formation —plastic buckling and folding—that occurs when ductile materials are subject to large-strain (>1) deformation. Using in situ imaging of simple-shear deformation of polypropylene, we delineate the mechanics underlying formation of folds and the resulting surface morphology. Here, the fold pattern shows a single wavelength that scales linearly with the size of the plastic zone. An analytical model predicts onset of the plastic buckling instability and characteristics of the folds. Besides obvious applications for controllably generating surface features over large areas, our results suggest a contributing mechanism for origin of surface roughness in solids.
Mechanisms influencing and prediction of tool influence function spots during hemispherical sub-aperture tool polishing on fused silica
Sub-aperture tool polishing of precision optics requires a detailed understanding of the local material removal [tool influence function (TIF)] at the contact spot between the workpiece and tool to achieve high removal determinism and hence precision of the optic relative to the desired/design surface figure. In this study, the mechanisms influencing and the quantitative prediction of the removal rate and shape of TIF spots during polishing of fused silica glass with cerium oxide slurry using a rotating hemispherical pad-foam tool for a wide variety of process conditions (including tool properties, kinematics, and applied displacements) are investigated. The TIF volumetric removal rate can be estimated utilizing the average relative velocity and contact area using a simple analytical model. In addition, stability of the volumetric removal rate for fixed process conditions is shown to be greatly dependent on the pad preparation and amount of tool use (affecting both pad topography and slurry buildup), whose general behavior shows an increase in removal rate followed by stabilization with polishing time. The determination of the TIF removal shape is more complex. An extended version of the Preston removal model is developed to explain a comprehensive set of measured TIF removal shapes to within ∼ <#comment/> 22 % <#comment/> . This model incorporates a number of phenomena impacting the TIF removal shape including: (a) temporal and spatial dependent relative velocity between the workpiece and tool; (b) an elastic mechanics based, as well as hydrodynamic, pressure distribution; (c) a spatially dependent friction coefficient possibly caused by both reduced slurry replenishment in low velocity regions and pad slurry islands (100 µm scale) and porosity (millimeter scale); and (d) a shear-based removal mechanism on the periphery of the contact spot.
Effect of tilt on circular zone plate performance
Fresnel zone plates are frequently used as focusing and imaging optics in x-ray microscopy, as they provide the ease of use of normal incidence optics. We consider here the effects of tilt misalignment on their optical performance, both in the thin optics limit and in the case of zone plates that are sufficiently thick so that volume diffraction effects come into play. Using multislice propagation, we show that simple analytical models describe the tilt sensitivity of thin zone plates and the thickness at which volume diffraction must be considered, and examine numerically the performance of example zone plates for soft x-ray focusing at 0.5 keV and hard x-ray focusing at 10 keV.