scholarly journals A Three-Dimensional Trajectory Model with Advection Correction for Tropical Cyclones: Algorithm Description and Tests for Accuracy

2019 ◽  
Vol 147 (9) ◽  
pp. 3145-3167 ◽  
Author(s):  
William Miller ◽  
Da-Lin Zhang

Abstract When computing trajectories from model output, gridded winds are often temporally interpolated to a time step shorter than model output intervals to satisfy computational stability constraints. This study investigates whether trajectory accuracy may be improved for tropical cyclone (TC) applications by interpolating the model winds using advection correction (AC) instead of the traditional linear interpolation in time (LI) method. Originally developed for Doppler radar processing, AC algorithms interpolate data in a reference frame that moves with the pattern translation, or advective flow velocity. A previously developed trajectory AC implementation is modified here by extending it to three-dimensional (3D) flows, and the advective flows are defined in cylindrical rather than Cartesian coordinates. This AC algorithm is tested on two model-simulated TC cases, Hurricanes Joaquin (2015) and Wilma (2005). Several variations of the AC algorithm are compared to LI on a sample of 10 201 backward trajectories computed from the modeled 5-min output data, using reference trajectories computed from 1-min output to quantify position errors. Results show that AC of 3D wind vectors using advective flows defined as local gridpoint averages improves the accuracy of most trajectories, with more substantial improvements being found in the inner eyewall where the horizontal flows are dominated by rotating cyclonic wind perturbations. Furthermore, AC eliminates oscillations in vertical velocity along LI backward trajectories run through deep convective updrafts, leading to a ~2.5-km correction in parcel height after 20 min of integration.

Algorithms ◽  
2021 ◽  
Vol 14 (5) ◽  
pp. 129
Author(s):  
Yuan Li ◽  
Ni Zhang ◽  
Yuejiao Gong ◽  
Wentao Mao ◽  
Shiguang Zhang

Compared with continuous elements, discontinuous elements advance in processing the discontinuity of physical variables at corner points and discretized models with complex boundaries. However, the computational accuracy of discontinuous elements is sensitive to the positions of element nodes. To reduce the side effect of the node position on the results, this paper proposes employing partially discontinuous elements to compute the time-domain boundary integral equation of 3D elastodynamics. Using the partially discontinuous element, the nodes located at the corner points will be shrunk into the element, whereas the nodes at the non-corner points remain unchanged. As such, a discrete model that is continuous on surfaces and discontinuous between adjacent surfaces can be generated. First, we present a numerical integration scheme of the partially discontinuous element. For the singular integral, an improved element subdivision method is proposed to reduce the side effect of the time step on the integral accuracy. Then, the effectiveness of the proposed method is verified by two numerical examples. Meanwhile, we study the influence of the positions of the nodes on the stability and accuracy of the computation results by cases. Finally, the recommended value range of the inward shrink ratio of the element nodes is provided.


Author(s):  
Arthur Ecoffet ◽  
Frédéric Poitevin ◽  
Khanh Dao Duc

Abstract Motivation Cryogenic electron microscopy (cryo-EM) offers the unique potential to capture conformational heterogeneity, by solving multiple three-dimensional classes that co-exist within a single cryo-EM image dataset. To investigate the extent and implications of such heterogeneity, we propose to use an optimal-transport-based metric to interpolate barycenters between EM maps and produce morphing trajectories. Results While standard linear interpolation mostly fails to produce realistic transitions, our method yields continuous trajectories that displace densities to morph one map into the other, instead of blending them. Availability and implementation Our method is implemented as a plug-in for ChimeraX called MorphOT, which allows the use of both CPU or GPU resources. The code is publicly available on GitHub (https://github.com/kdd-ubc/MorphOT.git), with documentation containing tutorial and datasets. Supplementary information Supplementary data are available at Bioinformatics online.


2015 ◽  
Vol 54 (3) ◽  
pp. 605-623 ◽  
Author(s):  
Anthony C. Didlake ◽  
Gerald M. Heymsfield ◽  
Lin Tian ◽  
Stephen R. Guimond

AbstractThe coplane analysis technique for mapping the three-dimensional wind field of precipitating systems is applied to the NASA High-Altitude Wind and Rain Airborne Profiler (HIWRAP). HIWRAP is a dual-frequency Doppler radar system with two downward-pointing and conically scanning beams. The coplane technique interpolates radar measurements onto a natural coordinate frame, directly solves for two wind components, and integrates the mass continuity equation to retrieve the unobserved third wind component. This technique is tested using a model simulation of a hurricane and compared with a global optimization retrieval. The coplane method produced lower errors for the cross-track and vertical wind components, while the global optimization method produced lower errors for the along-track wind component. Cross-track and vertical wind errors were dependent upon the accuracy of the estimated boundary condition winds near the surface and at nadir, which were derived by making certain assumptions about the vertical velocity field. The coplane technique was then applied successfully to HIWRAP observations of Hurricane Ingrid (2013). Unlike the global optimization method, the coplane analysis allows for a transparent connection between the radar observations and specific analysis results. With this ability, small-scale features can be analyzed more adequately and erroneous radar measurements can be identified more easily.


2015 ◽  
Vol 784 ◽  
pp. 292-299 ◽  
Author(s):  
Stephan Wulfinghoff ◽  
Marek Fassin ◽  
Stefanie Reese

In this work, two time integration algorithms for the anisotropic damage model proposed by Lemaitre et al. (2000) are compared. Specifically, the standard implicit Euler scheme is compared to an algorithm which implicitly solves the elasto-plastic evolution equations and explicitly computes the damage update. To this end, a three dimensional bending example is solved using the finite element method and the results of the two algorithms are compared for different time step sizes.


1988 ◽  
Vol 141 ◽  
Author(s):  
James P. Lavine ◽  
Gilbert A. Hawkins

AbstractA three-dimensional Monte Carlo computer program has been developed to study the heterogeneous nucleation and growth of oxide precipitates during the thermal treatment of crystalline silicon. In the simulations, oxygen atoms move on a lattice with randomly selected lattice points serving as nucleation sites. The change in free energy that the oxygen cluster would experience in gaining or losing one oxygen atom is used to govern growth or dissolution of the cluster. All the oxygen atoms undergo a jump or a growth decision during each time step of the anneal. The growth and decay kinetics of each nucleation site display interesting fluctuation phenomena. The time dependence of the cluster size generally differs from the expected 3/2 power law due to the fluctuations in oxygen arrival at and incorporation in a precipitate. Competition between growing sites and coarsening are observed.


Sensor Review ◽  
2015 ◽  
Vol 35 (2) ◽  
pp. 157-167 ◽  
Author(s):  
Shengbo Sang ◽  
Ruiyong Zhai ◽  
Wendong Zhang ◽  
Qirui Sun ◽  
Zhaoying Zhou

Purpose – This study aims to design a new low-cost localization platform for estimating the location and orientation of a pedestrian in a building. The micro-electro-mechanical systems (MEMS) sensor error compensation and the algorithm were improved to realize the localization and altitude accuracy. Design/methodology/approach – The platform hardware was designed with common low-performance and inexpensive MEMS sensors, and with a barometric altimeter employed to augment altitude measurement. The inertial navigation system (INS) – extended Kalman filter (EKF) – zero-velocity updating (ZUPT) (INS-EKF-ZUPT [IEZ])-extended methods and pedestrian dead reckoning (PDR) (IEZ + PDR) algorithm were modified and improved with altitude determined by acceleration integration height and pressure altitude. The “AND” logic with acceleration and angular rate data were presented to update the stance phases. Findings – The new platform was tested in real three-dimensional (3D) in-building scenarios, achieved with position errors below 0.5 m for 50-m-long route in corridor and below 0.1 m on stairs. The algorithm is robust enough for both the walking motion and the fast dynamic motion. Originality/value – The paper presents a new self-developed, integrated platform. The IEZ-extended methods, the modified PDR (IEZ + PDR) algorithm and “AND” logic with acceleration and angular rate data can improve the high localization and altitude accuracy. It is a great support for the increasing 3D location demand in indoor cases for universal application with ordinary sensors.


Author(s):  
Nicholas D. Leque ◽  
Ahmet Kahraman

Planet-to-planet load sharing is a major design and manufacturing tolerancing issue in planetary gear sets. Planetary gear sets are advantageous over their countershaft alternatives in many aspects, provided that each planet branch carries a reasonable, preferably equal, share of the torque transmitted. In practice, the load shared among the planets is typically not equal due to the presence of various manufacturing errors. This study aims at enhancing the models for planet load sharing through a three-dimensional formulation of N-planet helical planetary gear sets. Apart from previous models, the proposed model employs a gear mesh load distribution model to capture load and time dependency of the gear meshes iteratively. It includes all three types of manufacturing errors, namely constant errors such as planet pinhole position errors and pinhole diameter errors, constant but assembly dependent errors such as nominal planet tooth thickness errors, planet bore diameter errors, and rotation and assembly dependent errors such as gear eccentricities and run-outs. At the end, the model is used to show combined influence of these errors on planet load sharing to aid designers on how to account for manufacturing tolerances in the design of the gears of a planetary gear set.


2010 ◽  
Vol 3 (5) ◽  
pp. 4459-4495 ◽  
Author(s):  
C. López Carrillo ◽  
D. J. Raymond

Abstract. In this work, we describe an efficient approach for wind retrieval from dual Doppler radar data. The approach produces a gridded field that not only satisfies the observations, but also satisfies the anelastic mass continuity equation. The method is based on the so-called three-dimensional variational approach to the retrieval of wind fields from radar data. The novelty consists in separating the task into steps that reduce the amount of data processed by the global minimization algorithm, while keeping the most relevant information from the radar observations. The method is flexible enough to incorporate observations from several radars, accommodate complex sampling geometries, and readily include dropsonde or sounding observations in the analysis. We demonstrate the usefulness of our method by analyzing a real case with data collected during the TPARC/TCS-08 field campaign.


2014 ◽  
Vol 7 (5) ◽  
pp. 2181-2191 ◽  
Author(s):  
D. Rossi ◽  
A. Maurizi

Abstract. The development and validation of the vertical diffusion module of IL-GLOBO, a Lagrangian transport model coupled online with the Eulerian general circulation model GLOBO, is described. The module simulates the effects of turbulence on particle motion by means of a Lagrangian stochastic model (LSM) consistently with the turbulent diffusion equation used in GLOBO. The implemented LSM integrates particle trajectories, using the native σ-hybrid coordinates of the Eulerian component, and fulfils the well-mixed condition (WMC) in the general case of a variable density profile. The module is validated through a series of 1-D offline numerical experiments by assessing its accuracy in maintaining an initially well-mixed distribution in the vertical. A dynamical time-step selection algorithm with constraints related to the shape of the diffusion coefficient profile is developed and discussed. Finally, the skills of a linear interpolation and a modified Akima spline interpolation method are compared, showing that both satisfy the WMC with significant differences in computational time. A preliminary run of the fully integrated 3-D model confirms the result only for the Akima interpolation scheme while the linear interpolation does not satisfy the WMC with a reasonable choice of the minimum integration time step.


Sign in / Sign up

Export Citation Format

Share Document