Dynamic Pressure Field Research Articles

Experimental and Numerical Investigation of Self-Burial Mechanism of Pipeline with Spoiler under Steady Flow Conditions

Herein, hydraulic model experiments and numerical simulations were performed to understand the self-burial mechanism of subsea pipelines with spoilers under steady flow conditions. First, scour characteristics and self-burial functions according to the spoiler length-to-pipe diameter ratio (S/D) were investigated through hydraulic experiments. Further, the Navier–Stokes solver was verified. The experimental values of the velocity at the bottom of the pipeline with a spoiler and the pressure on the sand foundation where the pipeline rested were represented with the degree of conformity. Scour characteristics of a sand foundation were investigated from the numerical analysis results of the velocity and vorticity surrounding the pipelines with spoilers. The compilation of results from the hydraulic experiment and numerical analysis showed that the projected area increased when a spoiler was attached to the subsea pipes. This consequently increased the velocity of fluid leaving the top and bottom of the pipe, and high vorticity was formed within and above the sand foundation. This aggravated scouring at the pipe base and increased the top and bottom asymmetry of the dynamic pressure field, which developed a downward force on the pipeline. These two primary effects acting simultaneously under steady flow conditions explained the self-burial of pipelines with a spoiler attachment.

A Contactless Stage Based on Near-Field Acoustic Levitation for Object Handling and Positioning—Concept, Design, Modeling, and Experiments

A technology demonstrating a stage capable of positioning an object in the X – Y plane and performing rotation on a thin layer of air without contact is presented. By creating ultrasonic vibrations in an annulus structure, a thin air layer develops a dynamic acoustic pressure field able to levitate a planar object and to control its position. Motion is obtained by creating traveling pressure waves within the acoustic field giving rise to directional forces acting on the levitated object. The novel configuration is based on three independent acoustic levitation motors. Each capable of levitating and positioning the object simultaneously. Closed-loop control of the object's position is made possible by incorporating several sensors. In particular, the motion and rotation of the mirror-like reflective surface such as the silicon wafers were realized with optical sensors based on a laser mouse. The necessary forces for manipulating the acoustic vibrations and, therefore, the pressure fields are generated by an ${H_\infty }$ controller. The controller is robust enough to handle the strong coupling of the system and the model uncertainties. A laboratory prototype in which the concept was developed is described and validated by a series of experiments. This paper demonstrates the ability of the acoustic levitation systems to manipulate objects without contact.

Observation-based modelling of magnetised coronal mass ejections with EUHFORIA

Context. Coronal mass ejections (CMEs) are the primary source of strong space weather disturbances at Earth. Their geo-effectiveness is largely determined by their dynamic pressure and internal magnetic fields, for which reliable predictions at Earth are not possible with traditional cone CME models.Aims. We study two well-observed Earth-directed CMEs using the EUropean Heliospheric FORecasting Information Asset (EUHFORIA) model, testing for the first time the predictive capabilities of a linear force-free spheromak CME model initialised using parameters derived from remote-sensing observations.Methods. Using observation-based CME input parameters, we performed magnetohydrodynamic simulations of the events with EUHFORIA, using the cone and spheromak CME models.Results. Simulations show that spheromak CMEs propagate faster than cone CMEs when initialised with the same kinematic parameters. We interpret these differences as the result of different Lorentz forces acting within cone and spheromak CMEs, which lead to different CME expansions in the heliosphere. Such discrepancies can be mitigated by initialising spheromak CMEs with a reduced speed corresponding to the radial speed only. Results at Earth provide evidence that the spheromak model improves the predictions ofB(Bz) by up to 12–60 (22–40) percentage points compared to a cone model. Considering virtual spacecraft located within ±10° around Earth,B(Bz) predictions reach 45–70% (58–78%) of the observed peak values. The spheromak model shows inaccurate predictions of the magnetic field parameters at Earth for CMEs propagating away from the Sun-Earth line.Conclusions. The spheromak model successfully predicts the CME properties and arrival time in the case of strictly Earth-directed events, while modelling CMEs propagating away from the Sun-Earth line requires extra care due to limitations related to the assumed spherical shape. The spatial variability of modelling results and the typical uncertainties in the reconstructed CME direction advocate the need to consider predictions at Earth and at virtual spacecraft located around it.

Prediction of roll hydrodynamic characteristics for a ship hull section with a high-order finite volume solver

This work predicts hydrodynamic characteristics for a rolling ship hull section. Flows about the mid-section of a KVLCC2 under harmonic excited roll motions are simulated with solutions of Navier-Stokes (NS) equations filtered by a large eddy simulation (LES) with a finite volume (FV) solver. This solver applies a fourth-order compact scheme to discrete spatial terms, a fourth-order Runge-Kutta scheme to advance equations, a momentum interpolation method (MIM) to avoid pressure and velocity decoupling, a volume of fluid (VOF) method to capture free surface, and a ghost cell immersed boundary method (GCIBM) to simulate rolling ship. Firstly, verification is implemented to confirm the convergency. Validation is followed by comparison with experimental data. Then, systematic assessments are made on the effects of roll amplitude, frequency and bilge keels by analyzing the tendency of roll moment, damping and added mass moment of inertia as well as flow structures. Wave radiation and viscous effects are elucidated through flow structures. As a result, viscous effect is dominant. Roll amplitude, frequency and bilge keels influence roll hydrodynamic characteristics substantially. Besides, causes of vortex generation, shedding and movement, variation of dynamic pressure field and mechanism of bilge keel to mitigate roll motion are disclosed.

Development of a PVDF Sensor Array for Measurement of the Dynamic Pressure Field of the Blade Tip in an Axial Flow Compressor.

Tip clearance flow in axial flow compressor is unavoidable and responsible for pressure losses and noise generation and influences the stability of the compressor. However, necessary flow measurement in the blade tip region is a great challenge due to the small gap width as well as the structure limitation. In this paper, a polyvinylidene fluoride (PVDF) piezoelectric-film sensor array is developed to capture the dynamic pressure field over the blade tip in an axial flow compressor. The PVDF sensor array with 40 evenly distributed sensing points is fabricated directly on a 30 μm thick aluminum-metalized polarized PVDF film through photolithography. Dynamic calibration of the sensor is accomplished using acoustic source as excitation and a microphone as a reference. The test pressure range is up to 3.5 kPa and the sampling frequency is 20 kHz. The sensor presents a high signal-to-noise ratio and good consistency with the reference microphone. Sensitivity, frequency response, linearity, hysteresis, repeatability as well as the influence of temperature are also investigated through the calibration apparatus. The calibration gives credence to the relevance and reliability of this sensor for the application in dynamic pressure field measurement. The sensor is then applied to an actual measurement in a compressor. The output of the PVDF sensor array is also compared with the results of common pressure transducers, and the features of the dynamic pressure filed are discussed. The results indicate that the PVDF sensor array is capable of the dynamic pressure field measurement over the blade tip, and superior to the conventional approaches in installation, spatial resolution, frequency response, and cost. These advantages indicate its potential broad application in pressure measurement, especially for the complex spatial surface or thin-walled structure, such as the blade surface and the thin casing wall of the compressor.

Ungani Oil Field, Canning Basin - Evaluation of a Dolomite Reservoir

Core supported study of the heterogeneous Ungani dolomite reservoir architecture is driving development drilling and upgrades to field resource estimates. Vuggy connected and macro non-connected pore space was directly measured over a 70m continuous core using 3D structural analysis of CT-scans. Plug density measurements indicate non-connected or sub-140 micron resolution contribution of around 1% to 2.5% (pu) for the tight matrix, but all remaining porosity potentially contributes to oil production. The high resolution core porosity data is vertically repositioned and upscaled to calibrate neutron-density and sonic petrophysically derived porosities which are inadequate to resolve productive zones using conventional reservoir cut-offs. Conditioned resistivity image data correlated exceptionally with directly measured connected porosities. Reservoir properties were extrapolated to all wells across the Ungani field giving field net/gross estimates of up to 63% and porosities over 30% pu in some vuggy and brecciaed zones. The heterogeneity and prolific nature of the uppermost 17m of reservoir had not been previously recognised due to poor log data coverage and access at the casing points. Recent re-analysis of this section at Ungani-3 with Chemostrat ICP-OES-MS analysis of ditch cuttings was instrumental in proposing additional drilling to re-target this zone. Mineralogy analysis is used to calculate rock grain densities and help calibrate neutron-density derived porosity logs over the Ungani field. Up-scaled core porosity correlates well with density and sonic porosity logs. Resistivity logs adjusted for minerology can also be used to predict porosity and support the use of resistivity image logs to identify vuggy zones and estimate porosity at a higher resolution than conventional logging tools. A field static model was populated with three facies distributed over vertical zones according to the distribution encountered in the core porosity analysis and well logs, and iteratively matched to the dynamic pressure data and field production history which exhibits field scale multi-Darcy horizontal permeability and protection from vertical water cut. Further drilling and downhole artificial lift is planned to extend field production rates to 3000 bbls/day. Increased confidence in this regionally developed reservoir is supporting further exploration of undrilled prospects in this immature and under-explored trend

Leakage and Dynamic Force Coefficients for Two Labyrinth Gas Seals: Teeth-on-Stator and Interlocking Teeth Configurations. A Computational Fluid Dynamics Approach to Their Performance

Labyrinth gas seals (LSs) commonly used in turbomachines reduce secondary flow leakage. Conventional see-through labyrinth seal designs include either all teeth-on-stator (TOS) or all teeth-on-rotor (TOR). Experience shows that an interlocking labyrinth seal (ILS), with teeth on both stator and rotor, reduces gas leakage by up to 30% compared to the conventional see-through designs. However, field data for ILS rotordynamic characteristics are still vague and scarce in the literature. This work presents flow predictions for an ILS and a TOS LS, both seals share identical design features, namely radial clearance Cr = 0.2 mm, rotor diameter D = 150 mm, tooth pitch Li = 3.75 mm, and tooth height B = 3 mm. Air enters the seal at supply pressure Pin = 3.8, 6.9 bar (absolute) and temperature of 25 °C. The ratio of gas exit pressure to supply pressure ranges from 0.5 to 0.8, and the rotor speed is fixed at 10 krpm (surface speed of 79 m/s). The analysis implements a computational fluid dynamics (CFD) method with a multi-frequency-orbit rotor whirl model. The CFD predicted mass flow rate for the ILS is ∼ 21% lower than that of the TOS LS, thus making the ILS a more efficient choice. Integration of the dynamic pressure fields in the seal cavities, obtained for excitation frequency (ω) ranging from 12% to 168% of rotor speed (sub and super synchronous whirl), allows an accurate estimation of the seal dynamic force coefficients. For all the considered operating conditions, at low frequency range, the TOS LS shows a negative direct stiffness (K < 0), frequency independent; whereas the ILS has K > 0 that increases with both frequency and supply pressure. For both seals, the magnitude of K decreases when the exit pressure/inlet pressure ratio increases. On the other hand, the cross-coupled stiffness (k) from both seals is frequency dependent, its magnitude increases with gas supply pressure, and k for the ILS is more sensitive to a change in the exit/inlet pressure ratio. Notably, k turns negative for subsynchronous frequencies below rotor speed (Ω) for both the TOS LS and the ILS. The direct damping (C) for the TOS LS remains constant for ω > ½ Ω and has a larger magnitude than the damping for the ILS over the frequency range up to 1.5 Ω. An increase in exit/inlet pressure ratio decreases the direct damping for both seals. The effective damping coefficient, Ceff = (C-k/ω), whenever positive aids to damp vibrations, whereas Ceff < 0 is a potential source for an instability. For frequencies ω/Ω < 1.3, Ceff for the TOS LS is higher in magnitude than that for the ILS. From a rotordynamics point of view, the ILS is not a sound selection albeit it reduces leakage. Comparison of the CFD predicted force coefficients against those from a bulk flow model demonstrates that the later simple model delivers poor results, often contradictory and largely indifferent to the type of seal, ILS or TOS LS. In addition, CFD model predictions are benchmarked against experimental dynamic force coefficients for two TOS LSs published by Ertas et al. (2012, “Rotordynamic Force Coefficients for Three Types of Annular Gas Seals With Inlet Preswirl and High Differential Pressure Ratio,” ASME J. Eng. Gas Turbines Power, 134(4), pp. 04250301–04250312) and Vannini et al. (2014, “Labyrinth Seal and Pocket Damper Seal High Pressure Rotordynamic Test Data,” ASME J. Eng. Gas Turbines Power, 136(2), pp. 022501–022509.)

Experimental Force Coefficients for Two Sealed Ends Squeeze Film Dampers (Piston Rings and O-Rings): An Assessment of Their Similarities and Differences

Squeeze film dampers (SFDs) in aircraft engines effectively aid to reduce rotor motion amplitudes, in particular when traversing a critical speed, and help to alleviate rotor whirl instabilities. The current work is a long-term endeavor focused on quantifying the dynamic force performance of practical SFDs, exploring novel design damper configurations, and producing physically sound predictive SFD models validated by experimental data. Piston rings (PRs) and O-rings (ORs), commonly used as end seals in SFDs for commercial and military gas turbine engines, respectively, amplify viscous damping in a short physical length and while operating with a modicum of lubricant flow. This paper presents experimental force coefficients (damping and inertia) for two identical geometry SFDs with end seals, one configuration hosts PRs, and the other one ORs. The test rig comprises a stationary journal and bearing cartridge (BC) hosting the SFD and supported on four elastic rods to emulate a squirrel cage. The damper film land length, diameter, and clearance are L = 25.4 mm, D = 5L, and c = 0.373 mm (D/c = 340), respectively. A supply feeds ISO VG 2 oil to the film land at its middle plane through either one hole or three holes, 2.5 mm in diameter, 120 deg apart. In the PRSFD, the lubricant exits through the slit opening at the ring butted ends. The ORs suppress oil leakage; hence, lubricant evacuates through a 1 mm hole at ¼ L near one journal end. The ORs when installed add significant stiffness and damping to the test structure. The ORSFD produces 20% more damping than the PRSFD, whereas both sealed ends SFDs show similar size added mass. For oil supplied at 0.69 bar(g) through a single orifice produces larger damping, 60–80% more than when the damper operates with three oil feedholes. A computational model reproducing the test conditions delivers force coefficients in agreement with the test data. Archival literature calls for measurement of a single pressure signal to estimate SFD reaction forces. For circular centered orbits (CCOs), the dynamic pressure field, in the absence of any geometrical asymmetry or feed/discharge oil condition, “rotates” around the bearing with a speed equal to the whirl frequency. The paper presents force coefficients estimated from (a) measurements of the applied forces and ensuing displacements, and (b) the dynamic pressure recorded at a fixed angular location and “integrated” over the journal surface. The first method delivers a damping coefficient that is large even with lubricant supplied at a low oil supply pressure whereas the inertia coefficient increases steadily with feed pressure. Predictions show good agreement with the test results from measured forces and displacements, in particular the added mass. On the other hand, identified damping and inertia coefficients from dynamic pressures show a marked difference from one pressure sensor to another, and vastly disagreeing with test results from the first method or predictions. The rationale for the discrepancy relies on local distortions in the dynamic pressure fields that show zones of oil vapor cavitation at a near zero absolute pressure and/or with air ingestion producing high frequency spikes from bubble collapsing; both phenomena depend on the magnitude of the oil supply pressure. An increase in lubricant supply pressure suppresses both oil vapor cavitation and air ingestion, which produces an increase of both damping and inertia force coefficients. No prior art compares the performance of a PRSFD vis-à-vis that of an ORSFD. Supplying lubricant with a large enough pressure (flow rate) is crucial to avoid the pervasiveness of air ingestion. Last, the discussion on force coefficients obtained from two distinct methods questions the use of an oversimplifying assumption; the dynamic pressure field is not invariant in a rotating coordinate frame.

Power loss of slipper/swashplate based on elastohydrodynamic Lubrication model in axial piston pump

Energy saving of the hydraulic system is significantly realized by improving the efficiency of the pump. The slipper/swashplate lubricating surface is one of the main energy consumption parts in axial piston pump which is always used in the hydraulic system. Being a hydraulic bearing in axial piston pump, the slipper/swashplate pair exerts both the bearing function and the sealing function under the oscillating working condition. Because of viscous friction and leakage flow, power loss of slipper/swashplate is also one of the main power loss sources in axial piston pump. A precise model of lubricating surfaces, which includes oil film thickness, dynamic pressure field, load carrying ability and energy dissipation, is indispensable to develop more effective lubricating surface morphology. An elastohydrodynamic lubricating oil film model has been established for the sliding surfaces between slipper/swashplate interfaces. The model includes an isothermal fluid model, micro-motion of slipper and pressure deformation of the bounding solid bodies using a partitioned solution scheme. The power loss model between the slipper and swashplate has been developed and solved by numerical analysis method. Numerical results demonstrate that the diameter of the fixed damping orifice in the slipper has an important role in power loss of slipper/swashplate except the axial rotational velocity and operating pressure. Finally, the accuracy of numerical results of power loss is verified by experiments.

Assessing fragility of a reinforced concrete element to snow avalanches using a non-linear dynamic mass-spring model

Abstract. This paper presents an assessment of the fragility of a reinforced concrete (RC) element subjected to avalanche loads, and more generally to dynamic pressure fields applied orthogonally to a wall, within a reliability framework. In order to obtain accurate numerical results with supportable computation times, a light and efficient Single-Degree-of-Freedom (SDOF) model describing the mechanical response of the RC element is proposed. The model represents its dynamic mechanical response up to failure. Material non-linearity is taken into account by a moment–curvature approach, which describes the overall bending response. The SDOF model is validated under quasi-static and dynamic loading conditions by comparing its results to alternative approaches based on finite element analysis and the yield line theory. Following this, the deterministic SDOF model is embedded within a reliability framework to evaluate the failure probability as a function of the maximal avalanche pressure reached during the loading. Several reliability methods are implemented and compared, suggesting that non-parametric methods provide significant results at a moderate level of computational burden. The sensitivity to material properties, such as tensile and compressive strengths, steel reinforcement ratio, and wall geometry is investigated. The effect of the avalanche loading rate is also underlined and discussed. Finally, the obtained fragility curves are compared with respect to the few proposals available in the snow avalanche engineering field. This approach is systematic and will prove useful in refining formal and practical risk assessments. It could be applied to other similar natural hazards, which induce dynamic pressure fields onto the element at risk (e.g., mudflows, floods) and where potential inertial effects are expected and for which fragility curves are also lacking.

Thermoelasticity Coupling Transient Response of Composite‐Material Direction Pipe

Based on the positive scheme method, the thermal load of the jet flow in an inner composite‐material direction pipe is obtained, and the thermoelasticity coupling transient response is investigated. The positive scheme method with second‐order accuracy is extended for solving the axisymmetric Euler equations, and the supersonic axisymmetric jet flow over a missile afterbody containing jet exhaust is simulated. The correctness of the development for the positive scheme method is verified. With the developed positive scheme method used to simulate the jet flow in the inner direction pipe, the thermal load is obtained. The thermoelasticity coupling finite element model of the composite‐material direction pipe is established, and the stress response under dynamic pressure, unsteady temperature, and coupling state is obtained. Results show that, at the beginning of engine ignition, the effect of dynamic pressure and temperature field on the coupling stress is basically the same, and after that, the contribution of the temperature field to the coupling stress increases, and the thermal stress is the main factor affecting the strength of the composite‐material direction pipe.

Prediction of transmission, reflection and absorption coefficients of periodic structures using a hybrid Wave Based – Finite Element unit cell method

This paper presents a hybrid Wave Based Method – Finite Element unit cell method to predict the absorption, reflection and transmission properties of arbitrary, two-dimensional periodic structures. The planar periodic structure, represented by its unit cell combined with Bloch–Floquet periodicity boundary conditions, is modelled within the Finite Element Method, allowing to represent complex geometries and to include any type of physics. The planar periodic structure is coupled to semi-infinite acoustic domains above and/or below, in which the dynamic pressure field is modelled with the Wave Based Method, applying a wave function set that fulfills the Helmholtz equation and satisfies the Sommerfeld radiation condition and the Bloch–Floquet periodicity conditions inherently. The dynamic fields described within both frameworks are coupled using a direct coupling strategy, accounting for the mutual dynamic interactions via a weighted residual formulation. The method explicitly accounts for the interaction between the unit cell and the surrounding acoustic domain, also accounting for higher order periodic waves. The convergence of the method is analysed and its applicability is shown for a variety of problems, proving it to be a useful tool combining the strengths of two methods.

Modern system of multiphase flow in porous media and its development trend

Fluid flow in porous media is the key scientific problem in the development of oil and gas reservoirs. The traditional mechanics of fluid flow in porous media which based on the continuum hypothesis and Darcy′s law plays an important role in developing conventional oil and gas resources. In recent years, unconventional reservoirs are drawing more and more attention all over the world, therefore the development theory and technology, especially the corresponding flow mechanisms have become the hot research issues. The unconventional reservoirs exhibit distinct multiscale characteristics, even with six orders of magnitude difference. In addition, the application of massive multi-stage hydraulic fracturing can induce strong stress interactions. Therefore, the traditional theory of fluid flow in porous media cannot accurately describe the flow characteristics in unconventional reservoirs. In essence, the development of unconventional oil and gas resources involves multiphase fluids (e.g. oil, water and gas) flow in multi-scale porous media with multi-field coupling and various flow patterns. Therefore, the concept of modern system of multiphase flow in porous media is proposed, which means multiphase fluids flowing in multi-scale porous media with multi-field coupling and various flow patterns. The research status and development tendency are reviewed from the aspects of: (1) micro- and nanoscale oil and gas flow simulation; (2) upscaling for reservoir simulation, (3) macroscale flow simulation of unconventional oil and gas reservoirs; (4) simulation of flow in large scale fractured and vuggy carbonate reservoirs and (5) physical simulation of hydrocarbon transport in porous media. More specifically, in nanoscale the density functional theory and molecular simulation method can be used to study the interfacial phenomena to understand the hydrocarbon transport behavior in nanopores and provide key parameters for mesoscale flow simulation. The current study of nanoscale simulation mainly focuses on developing more realistic molecular structure model to represent the heterogeneous shale samples. Microscale simulation methods involve pore network model, lattice Boltzmann method, direct simulation of Navies-Stokes equation, level-set method and smoothed particle hydrodynamics, etc. Digital core and pore network model are the fundamental research platforms. Various methods can be used to reconstruct digital cores with multiscale pore structures and mineral compositions. The complex physicochemical phenomena namely adsorption/desorption, wettability change and boundary effect should be considered in the microscale flow simulations and extensive works have been done in microscale gas flow simulations. The future work on microscale simulation should focus on the multiphase flow mechanisms with multi-field coupling. The multiscale characteristics of unconventional reservoirs indicate the necessity of upscaling process to introduce the microscale flow mechanisms to macroscale. Homogenization theory and volume averaging method are the main upscaling approaches. Current upscaling methods are mostly based on the periodic boundary condition and are unreliable to be used in complex oil and gas reservoirs, which needs further study. In addition, more research needs to be conducted on the upscaling from molecular scale to mesoscale. In macroscale simulations of unconventional oil and gas reservoirs, the fluid-structure interaction should be considered and high efficiency numerical algorithm needs to be established. For large scale fractured and vuggy carbonate reservoirs, the non-Darcy flow characteristics and different flow regimes in vugs and fractures should be taken into account during flow simulation. Physical simulations of hydrocarbon transport in porous media are conducted at two scales: macroscale, nano- and microscale. Macroscale physical simulations aim at monitoring the dynamic saturation and pressure fields change under the realistic reservoir conditions. Nano- and microscale physical simulations are mainly applied to study the fluid transport mechanisms in single pore or throat. In summary, the proposed theory of multiphase fluids flowing in multi-scale porous media with multi-field coupling and various flow patterns can be applied to study the fluid flow problems in unconventional oil and gas industries.

ULF wave analysis and radial diffusion calculation using a global MHD model for the 17 March 2013 and 2015 storms

AbstractThe 17 March 2015 St. Patrick's Day Storm is the largest geomagnetic storm to date of Solar Cycle 24, with a Dst of −223 nT. The magnetopause moved inside geosynchronous orbit under high solar wind dynamic pressure and strong southward interplanetary magnetic field Bz causing loss; however, a subsequent drop in pressure allowed for rapid rebuilding of the radiation belts. The 17 March 2013 storm also shows similar effects on outer zone electrons: first, a rapid dropout due to inward motion of the magnetopause followed by rapid increase in flux above the prestorm level early in the recovery phase and a slow increase over the next 12 days. These phases can be seen in temporal evolution of the electron phase space density measured by the Energetic Particle, Composition, and Thermal Plasma Suite (ECT) instruments on Van Allen Probes. Using the Lyon‐Fedder‐Mobarry global MHD model driven by upstream solar wind measurements, we simulated both St. Patrick's Day 2013 and 2015 events, analyzing Lyon‐Fedder‐Mobarry electric and magnetic fields to calculate radial diffusion coefficients. These coefficients have been implemented in a radial diffusion code, using the measured electron phase space density following the local heating as the initial radial profile and outer boundary condition for subsequent temporal evolution over the next 12 days, beginning 18 March. Agreement with electron phase space density at 1000 MeV/G measured by the MagEIS component of the ECT instrument suite on Van Allen Probes was much improved using radial diffusion coefficients from the MHD simulations relative to coefficients parameterized by a global geomagnetic activity index.

MAVEN observations on a hemispheric asymmetry of precipitating ions toward the Martian upper atmosphere according to the upstream solar wind electric field

AbstractThe Mars Atmosphere and Volatile Evolution (MAVEN) observations show that the global spatial distribution of ions precipitating toward the Martian upper atmosphere has a highly asymmetric pattern relative to the upstream solar wind electric field. MAVEN observations indicate that precipitating planetary heavy ion fluxes measured in the downward solar wind electric field (−E) hemisphere are generally larger than those measured in the upward electric field (+E) hemisphere, as expected from modeling. The −E (+E) hemispheres are defined by the direction of solar wind electric field pointing toward (or away from) the planet. On the other hand, such an asymmetric precipitating pattern relative to the solar wind electric field is less clear around the terminator. Strong precipitating fluxes are sometimes found even in the +E field hemisphere under either strong upstream solar wind dynamic pressure or strong interplanetary magnetic field periods. The results imply that those intense precipitating ion fluxes are observed when the gyroradii of pickup ions are estimated to be relatively small compared with the planetary scale. Therefore, the upstream solar wind parameters are important factors in controlling the global spatial pattern and flux of ions precipitating into the Martian upper atmosphere.

The influences of solar wind pressure and interplanetary magnetic field on global magnetic field and outer radiation belt electrons

AbstractUsing the Van Allen Probe in situ measured magnetic field and electron data, we examine the solar wind dynamic pressure and interplanetary magnetic field (IMF) effects on global magnetic field and outer radiation belt relativistic electrons (≥1.8 MeV). The dynamic pressure enhancements (>2 nPa) cause the dayside magnetic field increase and the nightside magnetic field reduction, whereas the large southward IMFs (Bz‐IMF < −2nT) mainly lead to the decrease of the nightside magnetic field. In the dayside increased magnetic field region (magnetic local time (MLT) ~ 06:00–18:00, and L > 4), the pitch angles of relativistic electrons are mainly pancake distributions with a flux peak around 90° (corresponding anisotropic index A > 0.1), and the higher‐energy electrons have stronger pancake distributions (the larger A), suggesting that the compression‐induced betatron accelerations enhance the dayside pancake distributions. However, in the nighttime decreased magnetic field region (MLT ~ 18:00–06:00, and L ≥ 5), the pitch angles of relativistic electrons become butterfly distributions with two flux peaks around 45° and 135° (A < 0). The spatial range of the nighttime butterfly distributions is almost independent of the relativistic electron energy, but it depends on the magnetic field day‐night asymmetry and the interplanetary conditions. The dynamic pressure enhancements can make the nighttime butterfly distribution extend inward. The large southward IMFs can also lead to the azimuthal expansion of the nighttime butterfly distributions. These variations are consistent with the drift shell splitting and/or magnetopause shadowing effect.

A wave based method to predict the absorption, reflection and transmission coefficient of two-dimensional rigid frame porous structures with periodic inclusions

This paper presents an extension to the Wave Based Method to predict the absorption, reflection and transmission coefficients of a porous material with an embedded periodic set of inclusions. The porous unit cell is described using the Multi-Level methodology and by embedding Bloch-Floquet periodicity conditions in the weighted residual scheme. The dynamic pressure field in the semi-infinite acoustic domains is approximated using a novel wave function set that fulfils the Helmholtz equation, the Bloch-Floquet periodicity conditions and the Sommerfeld radiation condition. The method is meshless and computationally efficient, which makes it well suited for optimisation studies.

EMIC waves observed at geosynchronous orbit under quiet geomagnetic conditions (Kp ≤ 1)

AbstractWe statistically study the local time distribution of the helium band electromagnetic ion cyclotron (EMIC) waves observed at geosynchronous orbit when geomagnetic activity was low (Kp ≤ 1). In order to identify the geosynchronous EMIC waves, we use high time resolution magnetic field data acquired from GOES 10, 11, and 12 over a 2 year period from 2007 and 2008 and examine the local time distribution of EMIC wave events. Unlike previous studies, which reported high EMIC wave occurrence in the postnoon sector with a peak around 1500–1600 magnetic local time (MLT) during magnetically disturbed times (i.e., storm and/or substorm), we observed that quiet time EMIC waves mostly occur in a region from morning (∼0600 MLT) to afternoon (∼1600 MLT) with a peak around 1100–1200 MLT. To investigate whether the quiet time EMIC wave occurrence has a causal relationship with magnetospheric convection enhancement or solar wind dynamic pressure variations, we performed a superposed epoch analysis of solar wind parameters (solar wind speed, density, dynamic pressure, and interplanetary magnetic field Bz) and geomagnetic indices (AE and SYM‐H). From the superposed epoch analysis we found that solar wind dynamic pressure variation is a more important parameter than AE and SYM‐H for quiet time EMIC wave occurrence.

Periodic variations of oxygen EUV dayglow in the upper atmosphere of Venus: Hisaki/EXCEED observations

AbstractUsing the Extreme Ultraviolet Spectroscope for Exospheric Dynamics (EXCEED) aboard Hisaki and the Solar Extreme Ultraviolet Monitor on the Solar and Heliospheric Observatory, we investigate variations of the extreme ultraviolet (EUV) dayglow brightness for OII 83.4 nm, OI 130.4 nm, and OI 135.6 nm in the Venusian upper atmosphere observed in March–April (period 1), April–May (period 2), and June–July (period 3) in 2014. The result shows that characteristic periodicities exist in the dayglow variations other than the ~27 day solar rotational effect of the solar EUV flux: 1.8, 2.8, 3.1, 4.5, and 9.9 day in period 1; 1.1 day in period 2; and 1.0 and 11 day in period 3. Many of these periodicities are consistent with previous observations and theory. We suggest these periodicities are related to density oscillations of oxygen atoms or photoelectrons in the thermosphere. The cause of these periodicities is still uncertain, but planetary‐scale waves and/or gravity waves propagating from the middle atmosphere, and/or minor periodic variations of the solar EUV radiation flux may play a role. Effects of the solar wind parameters (velocity, dynamic pressure, and interplanetary magnetic field's intensity) on the dayglow variations are also investigated using the Analyser of Space Plasma and Energetic Atoms (ASPERA‐4) and magnetometer aboard Venus Express. Although clear correlation with the dayglow variations is not found, their minor periodicities are similar to the dayglow periodicities. Contribution of the solar wind to the dayglow remains still unknown, but the solar wind parameters might affect the dayglow variations.

Pressure Decimation and Interpolation (PDI) method for a baroclinic non-hydrostatic model

Non-hydrostatic models are computationally expensive in simulating density flows and mass transport problems due to the requirement of sufficient grid resolution to resolve density and flow structures. Numerical tests based on the Non-Hydrostatic Wave Model, NHWAVE (Ma et al., 2012), indicated that up to 70% of the total computational cost may be born by the pressure Poisson solver in cases with high grid resolution in both vertical and horizontal directions. However, recent studies using Poisson solver-based non-hydrostatic models have shown that an accurate prediction of wave dispersion does not require a large number of vertical layers if the dynamic pressure is properly discretized. In this study, we explore the possibility that the solution for the dynamic pressure field may, in general, be decimated to a resolution far coarser than that used in representing velocities and other transported quantities, without sacrificing accuracy of solutions. Following van Reeuwijk (2002), we determine the dynamic pressure field by solving the Poisson equation on a coarser grid and then interpolate the pressure field onto a finer grid used for solving for the remaining dynamic variables. With the Pressure Decimation and Interpolation (PDI) method, computational efficiency is greatly improved. We use three test cases to demonstrate the model’s accuracy and efficiency in modeling density flows.