Surface damage on honeycomb aircraft panels is often measured manually, and is therefore subject to variation based on inspection personnel. Eddy current testing (ECT) is sensitive to variations in probe-to-specimen spacing, or lift-off, and is thus promising for high resolution profiling of surface damage on aluminum panels. Lower frequency testing also allows inspection through the face sheet, an advantage over optical 3D scanning methods. This paper presents results from the ECT inspection of surface damage on an approximately flat aluminum honeycomb aircraft panel, and compares the measurements to those taken using optical 3D scanning technology. An ECT C-Scan of the dented panel surface was obtained by attaching the probe to a robotic scanning apparatus. Data was taken simultaneously at four frequencies of 25, 100, 400 and 1,600 kHz. A reference surface was then defined that approximated the original, undamaged panel surface, which also compensated for the effects of specimen tilt and thermal drift within the ECT instrument. Data was converted to lift-off using height calibration curves developed for each probe frequency. A damage region of 22,550 mm² area with dents ranging in depth from 0.13-1.01 mm was analyzed. The method was accurate at 1,600 kHz to within 0.05 mm (2s) when compared with 231 measurements taken via optical 3D scanning. Testing at 25 kHz revealed a 3.2 mm cell size within the honeycomb core, which was confirmed via destructive evaluation. As a result, ECT demonstrates potential for implementation as a method for rapid in-field aircraft panel surface damage assessment.

In the current industrial revolution, advanced technologies and methods can be effectively utilized for the detection and verification of defects in high-speed steel filament production. This paper introduces an innovative methodology for the precise detection and verification of micro surface defects found in steel filaments through the application of the Eddy current principle. Permanent magnets are employed to generates a magnetic field with high frequency surrounding a coil of sensor positioned at the filament's output end. The sensor's capacity to detect defects is validated through a meticulous rewinding process, followed by a thorough analysis involving scanning electron microscopy (SEM) and energy-dispersive spectroscopy (EDS). Artificial defects were intentionally introduced into a sample, and their amplitudes were monitored to establish a threshold value. The amplitude signal of these created defect was identified at approximately 10% FSH, which corresponds to a crack depth of about 20 µm. In the experimental production of 182 samples covering 38 km, the defect ratio was notably high, standing at 26.37%. These defects appeared randomly along the length of the samples. The verification results underscore the exceptional precision achieved in the detection of micro surface defects within steel filaments. These defects were primarily characterized by longitudinal scratches and inclusions containing physical tungsten carbide.

Electromagnetic eddy current sensors are commonly used to identify and quantify the surface notches of metals. However, the unintentional tilt of eddy current sensors affects results of size profiling, particularly for the depth profiling. In this paper, based on the eddy current thin-skin regime, a revised algorithm has been proposed for the analytical voltage or impedance of a tilted driver–pickup eddy current sensor scanning across a long ideal notch. Considering the resolution of the measurement, the bespoke driver–pickup, also termed as transmitter-receiver (T-R) sensor is designed with a small mean radius of 1 mm. Besides, the T-R sensor is connected to the electromagnetic instrument and controlled by a scanning stage with high spatial travel resolution , with a limit of 0.2 μm and selected as 0.25 mm. Experiments have been out on the voltage imaging of an aluminium sheet with 7 machined long notches of different depths using T-R sensor under different tilt angles. By fitting the measured voltage (both real and imaginary part) with proposed analytical algorithms, the depth profiling of notches is less affected by the tilt angle of sensors. From the results, the depth of notches can be retrieved within a deviation of 10 % for tilt angles up to 60 degrees.

Currently, electromagnetic and radar methods are the most common non-destructive tests for non-destructive rebar location in structures.. Both methods have some advantages, disadvantages and limitations. This article is an attempt to describe possibilities of particular methods with emphasizing their advantages and limitations. The first stage of tests included results from tests conducted on a lintel beam made of autoclaved aerated concrete (AAC) and reinforced with small-diameter rebars and rebars placed close to each other. The second stage consisted in testing nine lightweight concrete specimens, each with three bars having a diameter of 10, 16 and 20 mm at three different space arrangements, to determine the effect of bar diameter and spacing on the measurement accuracy. The reinforcement of lintels and specimens was tested with two different electromagnetic scanners and a ground-penetrating radar (GPR) device. Received results from direct tests were compared with results from non-destructive tests (NDT). NDT tests conducted in such an element were found to perform the correct assessment of the concrete cover, and at the same time to give ambiguous results with reference to shape and number of rebars.

In the field of quality control, the critical challenge of analyzing microdefects in steel filament holds significant importance. This is particularly vital as steel filaments serve as reinforced fibers in the use and applications within various component manufacturing industries. This paper addresses the crucial requirement of identifying and investigating microdefects on steel filaments. Eddy current signals are used for the identification of microdefects, and an in-depth investigation is conducted. The core objective is to establish the relationship between the depth of defects and the signals detected through the eddy current sensing principle. The threshold of eddy current instrument was set at 10%, corresponding to a created depth of 20 µm, to identify defective specimens. A total of 30 defective samples were analyzed, and the phase angles between the experimental and theoretical results were compared. The depths of defects ranged from 20 to 60 µm, with one sample having a depth exceeding 75 µm. The calculated threshold of 10.18% closely aligns with the set threshold of 10%, with a difference of only 1.77%. The resulting root mean square error (RMSE) was found to be 10.53 degrees, equivalent to 3.49 µm for the difference in depth and phase between measured results and estimated results. This underscores the methodology's accuracy and its applicability across diverse manufacturing industries.

it still undergoes corrosion when exposed to corrosive environments. This paper proposes an evaluation method for assessing the corrosion level of SPCC steel specimens using traditional eddy current testing (ECT), along with two different machine learning approaches. The objective is to classify the corroded state of the specimens into two states: a less corroded state (state-1) and a highly corroded state (state-2). Generative and discriminative models were implemented for classification. The generative classifier was based on the Gaussian mixture model (GMM), while the discriminative model was based on the logistic regression model. The models used features based on the absolute maximum of the perturbation magnetic field components at two different frequencies. The performance of the classifiers was evaluated using metrics such as absolute error, accuracy, precision, recall, and F1 score. The results indicate that high classification accuracy can be achieved based on both methods using traditional eddy current testing.

A reliable and efficient rail track defect detection system is essential for maintaining rail track integrity and avoiding safety hazards and financial losses. Eddy current (EC) testing is a non-destructive technique that can be employed for this purpose. The trade-off between spatial resolution and lift-off should be carefully considered in practical applications to distinguish closely spaced cracks such as those caused by rolling contact fatigue (RCF). A multi-channel eddy current sensor array has been developed to detect defects on rails. Based on the sensor scanning data, defect reconstruction along the rails is achieved using an inverse algorithm that includes both direct and iterative approaches. In experimental evaluations, the EC system with the developed sensor is used to measure defects on a standard test piece of rail with a probe lift-off of 4-6 mm. The reconstruction results clearly reveal cracks at various depths and spacings on the test piece.

In the field of railroad safety, effective detection of surface cracks is critical, necessitating reliable, high-speed non-destructive testing (NDT) methods. This study introduces a hybrid Eddy Current Testing (ECT) probe, specifically engineered for railroad inspection, to address the common issue of ’lift-off noise’ due to varying distances between the probe and test material. Unlike traditional ECT methods, this probe integrates transmit and differential receiver coils, aiming to enhance detection sensitivity and minimize lift-off impact. The study optimizes ECT probes employing different driver coils, emphasizing three main objectives: a) quantitatively evaluating each probe using signal- to-noise ratio (SNR) and outlining a real-time data processing algorithm based on SNR methodology; b) exploring the frequency range proximal to the electrical resonance of the receiver coil; c) examining sensitivity variations across varying lift-off distances. The experimental outcomes indicate that the newly designed probe with figure 8-shape driver coil significantly improves sensitivity in detecting surface cracks on railroads. It achieves an impressive SNR exceeding 100 for defects with minimal dimensions of 1 mm in width and depth. Simulation results closely align with experimental findings, validating the investigation of optimal operational frequency and lift-off distance for selected probe performance, determined to be 0.3 MHz and 0.5 mm, respectively. The realization of this project would lead to notable advancements in enhancing railroad safety by improving crack detection efficiency.

Electromagnetic sensing has been used for diverse applications of non-destructive testing, including the surface inspection, measurement of properties, object characterization. However, the measurement accuracy could be significantly influenced by the lift-off between sensors and samples. To address the issue caused by lift-offs, various strategies have been proposed for the permeability measurement of ferromagnetic steels, which mainly involves different sensor designs and signal features (e.g., the zero-crossing feature). In this paper, a single high-frequency scenario for the permeability retrieval is introduced. By combining the signal of two sensing pairs, the retrieval of magnetic permeability is less affected by the lift-off of sensors. Unlike the previous strategy on reducing the lift-off effect (directly taking the phase term out of the integration) using the Dodd-Deeds analytical method, the proposed method is based on a high-frequency linear feature of the phase term. Therefore, this method has the merit of high accuracy and fast processing for the permeability retrieval (a simplified version of Dodd-Deeds analytical formulas after the integration). Experimental measurement has been carried out on the impedance measurement of designed sensors interrogating ferromagnetic dual-phase steels. For sensor lift-offs of up to 10 mm, the error of the permeability retrieval is controlled within 4 % under the optimal frequency.

Eddy current testing can be used to interrogate steels but it is hampered by the lift-off distance of the sensor. Previously, the lift-off point of intersection (LOI) feature has been found for the pulsed eddy current (PEC) testing. In this paper, a lift-off invariant inductance (LII) feature is proposed for the multi-frequency eddy current (MEC) testing, which merely targets the ferromagnetic steels. That is, at a certain working frequency, the measured inductance signal is found nearly immune to the lift-off distance of the sensor. Such working frequency and inductance are termed as the lift-off invariant frequency (LIF) and LII. Through simulations and experimental measurements of different steels under the multi-frequency manner, the LII has been verified to be merely related to the sensor parameters and independent of different steels. By referring to the LIF of the test piece and using an iterative inverse solver, one of the steel properties (either the electrical conductivity or magnetic permeability) can be reconstructed with a high accuracy.

In many advanced industrial applications, the thickness is a critical index, especially for metallic coatings. However, the variance of lift-off spacing between sensors and test pieces affects the measured voltage or impedance, which leads to unreliable results from the sensor. Massive research works have been proposed to address the lift-off issue, but few of them applies to the thickness measurement of planar metallic films with finite-size circular (disk) geometry. Previously, a peak-frequency feature from the swept-frequency inductance was used to compensate the measurement error caused by lift-offs, which was based on the slow-changing rate of impedance phase term in the Dodd-Deeds formulas. However, the phase of measured impedance is nearly invariant merely on a limited range of sample thicknesses and working frequencies. Besides, the frequency sweeping is time-consuming, where a recalibration is needed for different sensor setups applied to the online real-time measurement. In this paper, a single-frequency algorithm has been proposed, which is embedded in the measurement instrument for the online real-time retrieval of thickness. Owing to the single-frequency measurement strategy, the proposed method does not need to recalibrate for different sensor setups. The thickness retrieval is based on a triple-coil sensor (with one transmitter and two receivers). The thickness of metallic disk foils is retrieved from the measured electrical resistance of two transmitter-receiver sensing pairs. Experiments on materials of different electrical conductivities (from direct current), thicknesses and planar sizes (radii) have been carried out to verify the proposed method. The error for the thickness retrieval of conductive disk foils is controlled within 5 % for lift-offs up to 5 mm.

Zr-2.5%Nb Pressure tubes (PT) house uranium fuel bundles in the fuel channels of CANDU® nuclear reactors. Preventing a failure mode caused by contact of the PT with an outer calandria tube (CT) is performed by inspection using eddy current (EC) testing and ultrasonic testing (UT) to measure the PT-CT gap. EC gap measurements are particularly sensitive to circumferential variation of the PT’s electrical resistivity, due to microstructural variations. A full-factorial experiment was performed to examine the statistical significance of variations in EC test parameters and manufacturing conditions on average circumferential electrical resistivity of as-manufactured PTs. It was found that 79% of the variance in the data could be attributed to variations caused by any of the test factors or combinations of test factors. The parameters that accounted for the majority of the variance were; 1) heat treatment (HT); 2) HT and EC frequency; 3) probe inner or outer surface placement; and 4) EC frequency. Measurements of circumferential resistivity showed up to ± 2.3% variation from the average of either surface. HT caused average PT resistivity to decrease at a rate of 1.53±0.08 (μΩ⋅cm)/log⁡(hr) and 1.1±0.4 (μΩ⋅cm)/log⁡(hr) for inner and outer PT surfaces, respectively. Results are correlated with literature reported differences in average β_Zr ribbon thickness in the axial-transverse cross-section between inner and outer PT surfaces. Results demonstrate potential for EC based resistivity measurements to characterize variations and changes in microstructure of Zr-2.5%Nb PT material.

In order to solve the problem of large eddy current loss and high temperature rise caused by large amount of permanent magnet, a new type of combined magnetic pole permanent magnet motor is proposed in this paper. The sinusoidal distribution subdomain model of non-uniform air gap rotor is established by using Laplace equation, and the analytical expression of eddy current loss of rotor in uniform air-gap and non-uniform air-gap is derived. The effect of rotor eccentricity on eddy current loss is obtained. According to the characteristics of distributed winding of non-uniform air gap combined pole permanent magnet motor, equivalent treatment is performed to obtain the equivalent thermal conductivity, so as to establish equivalent thermal network model of motor, calculate the temperature value of each component of the motor, and verify the correctness of the thermal network model through the method of magnetothermal bidirectional coupling. Finally, an experimental platform is set up to carry out temperature rise experiments on the two prototypes respectively. Experimental results show that the non-uniform air-gap rotor structure can effectively reduce the rotor eddy current loss and motor temperature rise, and also verify the accuracy of the analytical model calculation results.

Eddy current sensors can be used to test the characteristics and measure the parameters of the conductive samples. As the main obstacle of the multi-frequency eddy current sensor, the lift-off distance affects the effectiveness and accuracy of the measurement. In this paper, a material-independent algorithm has been proposed for the restoration of the lift-off distance when using the multi-frequency eddy current sensor, which is based on the approximation under the thin-skin effect. Experiment testing on the performance of the proposed method is presented. Results show that from the dual-frequency inductance, the lift-off distance could be restored with a maximum error of 0.24 mm for the distance up to 12 mm. Besides, the derived lift-off distance is used for the inversion of the magnetic permeability. Based on a lift-off insensitive inductance (LII) feature, the magnetic permeability of steels can be inversed in an iterative manner, with an error of less than 0.6 % for the lift-off distance up to 12 mm.

As the number of end-of-life vehicles continues to grow, the field of separation technology for broken copper and aluminium parts of end-of-life vehicles has emerged as a significant area of research. As a key technology for the separation of non-magnetic non-ferrous metals, such as copper and aluminium, eddy current separation can achieve efficient separation of non-magnetic non-ferrous metals according to their differing conductivity and density. Nevertheless, further research is required in the separation of large-scale broken copper and aluminium parts of end-of-life vehicles. This paper presents a novel magnetic roller model comprising circumferential and axial periodic arrangements of permanent magnets. A single permanent magnet was constructed as an equivalent current model, and the magnetic scalar potential equations were solved using simulation tools. This allowed the distribution of the magnetic field outside the magnetic roller to be studied. The results demonstrate that the novel permanent magnet array magnetic roller can enhance the magnetic field intensity and action distance of the external magnetic field, which is more conducive to the separation of large-scale broken copper and aluminium parts of scrap cars. The findings of this research provide a theoretical foundation for the separation of scrap car broken copper and aluminium parts.

Defect detection in ferromagnetic substrates is often hampered by non-magnetic coating thickness variation when using conventional eddy current testing technique. The lift-off distance between the sample and the sensor is one of the main obstacles for the thickness measurement of non-magnetic coatings on ferromagnetic substrates when using the eddy current testing technique. Based on the eddy current thin-skin effect and the lift-off insensitive inductance (LII), a simplified iterative algorithm is proposed for reducing the lift-off variation effect using a multi-frequency sensor. Compared to the previous techniques on compensating the lift-off error (e.g., the lift-off point of intersection) while retrieving the thickness, the simplified inductance algorithms avoid the computation burden of integration, which are used as embedded algorithms for the online retrieval of lift-offs via each frequency channel. The LII is determined by the dimension and geometry of the sensor, thus eliminating the need for empirical calibration. The method is validated by means of experimental measurements of the inductance of coatings with different materials and thicknesses on ferrous substrates (dual-phase alloy). The error of the calculated coating thickness has been controlled to within 3 % for an extended lift-off range of up to 10 mm.

For the electromagnetic eddy current testing, various methods have been proposed for reducing the lift-off error on the measurement of samples. In this paper, instead of eliminating the measurement error caused by the lift-off effect, an algorithm has been proposed to directly measure the lift-off distance between the sensor and non-magnetic conductive plates. The algorithm is based on a sample-independent inductance (SII) feature. That is, under high working frequencies, the inductance is found sensitive to the lift-off distance and independent of the test piece under an optimal single high working frequency (43.87 kHz). Furthermore, the predicted lift-off distance is used for the thickness prediction of the non-magnetic conductive samples using an iterative method. Considering the eddy current skin depth, the thickness prediction is operated under a single lower frequency (0.20 kHz). As the inductance has different sensitivities to the lift-off and thickness, the prediction error of the sample thickness is different from that of the lift-off distance. From the experiments on three different nonmagnetic samples – aluminium, copper, and brass, the maximum prediction error of the lift-off distance and sample thickness is 1.1 mm and 5.42 % respectively at the lift-off of 12.0 mm.

The shear stress generated by wall turbulence is the main cause of wall friction resistance in turbulent flow through pipes. This paper investigates the impact of inserting rest periods (regions of constant Reynolds number) within the pulsating operating cycle of velocity on the fully turbulent flow at large time-averaged Reynolds numbers, using the Large Eddy Simulation (LES) method. The study aims to explore the effect of increasing rest periods within pulsations on drag resistance. The dimensionless shear stress and drag reduction rates during different time periods of rest were analyzed. Numerical simulation results indicate that the pulsating velocity operation mode does not necessarily lead to drag reduction; it may even result in increased resistance. Inserting rest periods within the pulsation cycle can achieve drag reduction effects, with the maximum drag reduction rate reaching 21.8\%. Comparisons with experimental data and Direct Numerical Simulation (DNS) from the literature validate the feasibility of using the LES method for pipe pulsating operation modes.

This paper analyzes and compares models for predicting average magnet losses in interior permanent-magnet motors with fractional-slot concentrated windings due to harmonics in the armature reaction (assuming sinusoidal phase currents). Particularly, loss models adopting different formulations and solutions to the Helmholtz equation to solve for the eddy currents are compared to a simpler model relying on an assumed eddy-current distribution. Boundaries in terms of magnet dimensions and angular frequency are identified (numerically and using an identified approximate analytical expression) to aid the machine designer whether the more simple loss model is applicable or not. The assumption of a uniform flux-density variation (used in the loss models) is also investigated for the case of V-shaped and straight interior permanent magnets. Finally, predicted volumetric loss densities are exemplified for combinations of slot and pole numbers common in automotive applications.

Electromagnetic sensors have been used for inspecting small surface defects of metals. Based on the eddy-current thin-skin regime, a revised algorithm is proposed for a triple-coil drive-pickup eddy-current sensor scanning over long surface crack slots (10 mm) with different rotary angles. The method is validated by the voltage measurement of the designed EC sensor scanning over a benchmark (ferromagnetic) steel with surface defects of different depths and rotary angles. With an additional sensing coil for the designed EC sensor, the defect angle (or orientation) can be measured without spatially and coaxially rotating the excitation coil. By referring to the voltage change (due to the defect) diagram (voltage sum versus voltage different) of two sensing pairs, the rotary angle of the surface crack is retrieved with a maximum residual deviation of 3.5 %.

Previously, various techniques have been proposed for reducing the lift-off effect on the thickness measurement of the non-magnetic films, including the peak-frequency feature and phase feature in the Dodd-Deed analytical formulation. To realise a real-time feedback response on the thickness monitoring, the phase term in the Dodd-Deeds formulation must be taken off the integration. Previous methods were based on the slow change rate of the phase term when compared to the rest of the term – the magnitude term. However, the change rate of the phase term is still considerable for a range of working frequencies. In this paper, a high-frequency feature has been found. That is, the ratio between the imaginary and real part of the phase term is proportional to the integral variable under high frequencies. Based on this proportion relationship, the phase term has been taken out; and a thickness algorithm has been proposed. By combing the measured impedance from the custom-built sensor (three coils), the thickness of the metallic film can be reconstructed. Experiments have been carried out for the verification of the proposed scenario. Results show that the thickness of the metal film can be reconstructed with a small error of less than 2 %, and immune to a reasonable range of lift-offs.

This paper proposes a linear eddy-current feature to determine the radius of a metallic ball in a non-contact manner. An electromagnetic eddy-current sensor with two coils is placed co-axially to the metal ball during measurement. It is well known that the distance between the sensor and test piece (i.e. lift-off) affects eddy-current signals. In this paper, it is found that the peak frequency feature of inductance spectrum is linear to the lift-off spacing between the centre of coil and ball. Besides, the slope of peak frequencies versus lift-offs is linked to the radius of ball. The radius of metallic balls is retrieved from the experimental and embedded analytical result of the slope. Measurements have been carried out on 6 metallic balls with different radii. The radius of the metallic ball can be retrieved with an error of less than 2 %.

Previously, a conductivity invariance phenomena (CIP) has been discovered – at a certain lift-off, the inductance change of the sensor due to a test sample is immune to conductivity variations, i.e. the inductance – lift-off curve passes through a common point at a certain lift-off, termed as conductivity invariance lift-off. However, this conductivity invariance lift-off is fixed for a particular sensor setup, which is not convenient for various sample conditions. In this paper, we propose using two parameters in the coil design – the horizontal and vertical distances between the transmitter and the receiver to control the conductivity invariance lift-off. The relationship between these two parameters and the conductivity invariance lift-off is investigated by simulation and experiments and it has been found that there is an approximate linear relationship between these two parameters and the conductivity invariance lift-off. This is useful for applications where the measurements have restrictions on lift-off, e.g. uneven coating thickness which limits the range of the lift-off of probe during the measurements. Therefore, based on this relationship, it can be easier to adjust the configuration of the probe for a better inspection of the test samples.

In large eddy simulation (LES) of turbulent flows, the most critical dynamical processes to be considered by dynamic subgrid models to account for an average cascade of kinetic energy from the largest to the smallest scales of the flow is not fully clear. Furthermore, evidence of vortex stretching being the primary mechanism of the cascade is not out of the question. In this article, we study some essential statistical characteristics of vortex stretching and its role in dynamic approaches of modeling subgrid-scale turbulence. We have compared the interaction of subgrid stresses with the filtered quantities among four models using invariants of the velocity gradient tensor. This technique is a single unified approach to studying a wide range of length scales in the turbulent flow. In addition, it also provides a rational basis for the statistical characteristics a subgrid model must serve in physical space to ensure an appropriate cascade of kinetic energy. Results indicate that the stretching mechanism extracts energy from the large-scale straining motion and passes it onto small-scale stretched vortices.

In situ aircraft observations in typhoons have been scarce. This paper documents and analyses the aircraft and dropsonde data collected in Super Typhoon Saola (2023) over the northern part of the South China Sea. The wind and turbulence structures of the typhoon are investigated. The turbulence intensities are quantified in terms of turbulent kinetic energy (TKE) and eddy dissipate rate (EDR), and these data are compared with other available estimates of turbulence intensities, such as those based on weather radars, meteorological satellites, and numerical weather predic-tion (NWP) models. It is found that the TKE and EDR are closely correlated, and they are con-sistent with the weather radar/satellite observations as well as NWP-based outputs. Furthermore, the boundary layer inflow, vertical wind profiles, and atmospheric stability are analysed based on the dropsonde observations. The analysed results would advance the understanding of typhoon structures and offer references for the validation of remote sensing observations and NWP models.

Additive Manufacturing (AM) Direct Laser Fabrication (DLF) of Ti-5Al-5V-5Mo-3Cr (Ti5553) is being developed as a method for producing aircraft components. The additive manufacturing process can produce flaws near the surface, such as porosity and material voids, which act as stress raisers, leading to potential component failure. Eddy current testing was investigated to detect flaws on or near the surface of DLF Ti5553 bar samples. For this application, the objective was to develop an eddy current probe capable of detecting flaws 500 µm in diameter, located 1 mm below the component's surface. Two initial sets of coil parameters were chosen: The first, based on successful experiments that demonstrated detection of a near surface flaw in Ti5553 using a transmit-receive array probe, and the second, derived from simulation by Finite Element Method (FEM). An optimized transmit receive coil design, based on the FEM simulations, was constructed. The probe was evaluated on Ti5553 samples containing sub-surface voids of the target size, as well as samples with side-drilled holes and samples with holes drilled from the opposing inspection surface. The probe was able to effectively detect 80% of the sub-surface voids. Limitations included the probe’s inability to detect sub-surface voids near sample edges and a sensitivity to surface roughness, which produces local changes in lift-off. Multifrequency mixing improved signal-to-noise ratio when surface roughness was present on average by 22%. A probe based on that described in this paper could benefit quality assurance of additively manufactured aircraft components.

The three-dimensional structure of Gulf of Mexico’s warm-core rings, detaching from the Loop Current, is investigated using satellite altimetry and a large set of ARGO float profiles. Reconstruction of the Loop Current rings (LCRs) vertical structure from sea surface height observations is made possible by the use of the gravest empirical modes method (GEM). The GEMs are transfer functions that associate a value of temperature and salinity for each variable pair {dynamic height;pressure}, and are computed by estimating an empirical relationship between dynamic height and the vertical thermohaline structure of the ocean. Between 1993 and 2021, 40 LCRs were detected in the altimetry and their three-dimensional thermohaline structure was reconstructed, as well as a number of dynamically relevant variables (geostrophic and cyclogeostrophic velocity, relative vorticity, potential vorticity, available potential energy and kinetic energy density etc.). The structure of a typical LCR was computed by fitting an analytical stream function to the LCRs dynamic height signature and reconstructing its vertical structure with the GEM. The total heat and salt contents and energy of each LCR were computed and their cumulative effect on the Gulf of Mexico’s heat, salt and energy balance is discussed. We show that LCRs have a dramatic impact on these balances and estimate that residual surface heat fluxes of -13 Wm−2 are necessary to compensate their heat input, while the fresh water outflow of the Mississippi river approximately compensates for their salt excess input. An average energy dissipation of O[10−10−10−9] Wkg−1 would be necessary to balance their energy input.

This paper presents a new approach of eddy current methods for determining the size of the co-axial hole in the metallic circular disk. In recent decades, for the air-cored sensor probe, the impedance change due to the presence of an infinite metal plate can be calculated by the Dodd-Deeds model. However, in practical measurements, the sample cannot match with the condition required - ‘infinite’, thus the Dodd-Deeds model could not be applied to the disk with finite size and certainly not a co-axial hole in the center. In this paper, a dual-constraint analytical method is proposed. That is, the upper and lower limits of the integration are substituted with specific values instead of the original 0 and . Besides, it is found that, once the outer radius of the disk is fixed (i.e. the lower limit of integration is fixed), the upper limit reduces linearly as the size of the coaxial hole increases. Both the FEM simulation and experiments have been carried out to validate this method. The radius of the hole can be estimated based on the dual-constraint integration feature.

The present study investigates the general characteristics of mesoscale eddies in the Global Ocean, as identified and tracked by the TOEddies algorithm (Laxenaire et al., 2018) implemented on a global scale. Applied to satellite observations (AVISO/DUACS) of Absolute Dynamic Topography (ADT), TOEddies provides daily information on eddy dynamical characteristics (e.g., size and intensity) over a 30-year period (1993-2023) and identifies complex eddy-eddy interactions that lead to eddy splitting and merging. By capturing the multiple merging and splitting events that a single eddy can undergo, TOEddies generates a complex network of eddies that challenges the conventional view of a single trajectory associated with a single eddy. Furthermore, it offers an original assessment of interconnected eddy trajectories. A statistical description of the average characteristics of the eddies and their spatial distributions (generation, disappearance and merging/dividing is first provided. A comparative analysis of existing global eddy datasets is then presented. Among the years of observations, several coherent, long-lived mesoscale eddies with lifetimes exceeding 1.5 years stand out. However, we find that only a small fraction of these eddies is comparable between datasets, while their dynamic properties can differ substantially. Furthermore, the detection of eddies from altimetry is combined with 23 years of Argo profile co-localised measurements (2000-2023), allowing us to further investigate their contribution to the global subsurface ocean.

This paper presents a novel theoretical and computational methodology for the generation of an onset turbulent field with prescribed properties in the numerical simulation of an arbitrary viscous flow. The methodology is based on the definition of a suitable distribution of volume force terms in the right–hand side of the Navier–Stokes equations. The distribution is represented by harmonic functions that are randomly variable in time and space. The intensity of the distribution is controlled by a simple PID strategy in order to obtain that the generated turbulent flow matches a prescribed turbulence intensity. A further condition is that a homogeneous isotropic flow is es- tablished downstream of the region where volume force terms are imposed. Although it is general, the proposed methodology is primarily intended for the computational modelling of hydrokinetic turbines in turbulent flows representative of tidal or riverine installations. A first numerical applica- tion is presented by considering the injection of homogeneous and isotropic turbulence with 16% intensity into a uniform unbounded flow. The analysis of statistical properties as auto-correlation, power spectral density, probability density functions, demonstrates that the generated flow tends to achieve satisfactory levels of stationarity and isotropy, whereas the simple control strategy used determines overestimated turbulent intensity levels.

We have found a region in the northwestern Japan Sea near the Subpolar Front where mesoscale eddies regularly form and circulate. The anticyclonic eddy with a strong dynamic feature, vertically extended to the bottom (~3500 m), multilayer core structure and extreme values of temperature and salinity was sampled there in the cruise in May 2004. The altimetry-based eddy tracking showed that it had a life span of nine months with the size reaching 120 km. CTD observational data were used to explore the eddy’s features and vertical structure. The eddy had a multilayer core with different thermohaline characteristics which could be originated by eddy interaction with the surrounding waters during its life. The evolution of the eddy has been studied and analyzed with the help of altimetry-based Lagrangian indicators of water motion. Performing the particle-tracking numerical experiments, we computed the Lagrangian maps that contain information on the origin of water inside the eddy core and its ‘age’ on every day of the eddy’s life. Inspecting daily-computed Lagrangian maps, we documented the essential events in the evolution of the eddy including its formation, splitting, merger, entrainment and detrainment of water, erosion and eventual decay. All these observations have been verified with the infrared satellite images. The new Lagrangian technique has been applied to calculate day by day the content of water of various origin (coastal, subtropical and subarctic) entrained into the eddy core. It has been found that the surface core has been filled mainly with subtropical water originated in the southern flank of the Subpolar Front. The episodes with cardinal changes of ratio of different water types have been recorded and verified with the infrared satellite images and maps of the origin of water inside the eddy core.

This investigation aims to discuss the formation process of eddies and the heat transportation in plasma keyhole arc welding. In order to clarify this issue, the measurement of the convection inside the weld pool, the convection on the weld pool surface, also the temperature distribution on the weld pool surface were carried out. The results showed that two eddies were found in the weld pool, which is controlled mainly through the shear force by the plasma flow acting on the weld pool surface. The magnitude, extent and direction of the shear force are thought to be determined primarily by the variation of keyhole profile. The relative shape and strength of each eddy is largely changed depending on the change of the keyhole profile when nozzle diameter changed. These relative strengths of each eddy are considered to decisively govern the heat transport in the weld pool coinciding with the direction of eddies. A larger eddy near the lower part of the keyhole inside the weld pool was found out in the case of 1.6 mm, meanwhile a upward larger eddy was found out near the upper part of the keyhole inside the weld pool in the case of 2.4 mm.

Deterministic-stochastic empirical mode decomposition (EMD) is used to obtain low-frequency (non-diffusive, i.e., background velocity) and high-frequency (diffusive, i.e., eddies) components from a Lagrangian drifter‘s trajectory. Eddy characteristics are determined from the time series of eddy trajectories from individual Lagrangian drifter such as the eddy radial scale, eddy velocity scale, eddy Rossby number, and eddy-background kinetic energy ratio. A long-term dataset of the SOund Fixing And Ranging float time series obtained near the California coast by the Naval Postgraduate School from 1992 to 2004 at depth between 150 and 600 m (http://www.oc.nps.edu/npsRAFOS/) is used as an example to demonstrate the capability of the deterministic-stochastic EMD.

The characteristics of the eddy mass transport are estimated depending on the values of the parameters of a large-scale flow that forms under the conditions of the shelf seas in the Arctic. For this, the results of numerical simulation of the Kara Sea with a horizontal resolution permitting the development of mesoscale eddies are used. The parameters resulting from numerical experiment are considered as a statistical sample and are analyzed using methods of sensitivity study and clustering of sample elements. Functional dependencies are obtained that are closest to the simulated distributions of quantities. These expressions make it possible, within the framework of large-scale models, to evaluate the characteristics of the cross-isobatic eddy mass transport in the diffusion approximation with a counter-gradient flux. Numerical experiments using the SibCIOM model showed that areas along the Fram branch of the Atlantic waters trajectory in the Arctic as well as the shelf of the East Siberian and Laptev seas with adjacent deep water areas are most sensitive to proposed parameterization of eddy exchanges. Accounting for counter-gradient eddy fluxes turned out to be less important.

Studying the interaction between the upper ocean and the typhoons is crucial to improve our understanding of heat and momentum exchange between the ocean and the atmosphere. In recent years, the upper ocean responses to typhoons have received considerable attention. The sea surface cooling (SSC) process has been repeatedly discussed. In the present work, case studies were examined on five strong and super typhoons that occurred in 2016—LionRock, 1610; Meranti, 1614; Malakas, 1616; Megi, 1617; and Chaba 1618—to search for more evidence and new features of typhoon’s impact on the sea surface environment. The typhoon monitoring data from the Central Meteorological Observatory, the sea surface temperature (SST) data from satellite microwave and infrared remote sensing, and the sea surface height anomaly (SSHA) data from satellite altimeters were used to analyze in detail: the SSC features caused by typhoons, the relationship between the SSC and the typhoon travelling speed, and the variations in cold and warm eddies during typhoon passage. Results showed that: (1) SSC generally occurred during typhoon passage and the degree of SSC was always determined by the strength and the travelling speed of the typhoon, as well as the initial SST. (2) One day before or on the day of typhoon passage, the SSHA slightly increased due to low surface pressure. After the typhoon passed, the SSHA obviously decreased along with the SSC. The pre-existing positive SSHAs, which always represent warm eddies, decreased or disappeared during typhoon passage, whereas negative SSHAs or cold eddies were enhanced. (3) New cold eddies were generated, especially at the turning points of the typhoon path. The presence of warm eddies is suggested to have a strengthening effect on the typhoons.

Lift-offs of the sensor could significantly affect the measurement signal and reconstruction of material properties when using the electromagnetic (inductive) eddy current sensor. Previously, various methods (including novel sensor designs, and features like zero-crossing frequency, lift-off point of intercept) have been used for eliminating the measurement error caused by the lift-off distance effect of the sensor. However, these approaches can only be applied for a small range of lift-off variations. In this paper, a linear relationship has been found between the sensor lift-off and ratio of dual-frequency eddy current signals, particularly under the high working dual frequencies. Based on this linear relationship, the lift-off variation can be reconstructed firstly with a small error of 2.5 % when its actual value is up to 10 mm (10.1 % for 20 mm). The reconstructed lift-off is used to further get the property of the material under a low single frequency. Experiments on different ferrous metals have been carried out for the testing of the reconstruction scheme. Since the inductance is more sensitive to the material property (and less sensitive to the lift-off) under low frequencies, the reconstruction error of the material property is much smaller than that of the lift-off, with 1.4 % under 12 mm (and 4.5 % under 20 mm).

Subgrid scale models for unresolved turbulent fluxes are investigated with a focus on combustion for space propulsion applications. An extension to the gradient model is proposed, introducing a dependency on the local burning regime. The dynamic behavior of the model's coefficients is investigated, and scaling laws are studied. The discussed models are validated using a DNS database of a high-pressure turbulent fuel-rich methane-oxygen diffusion flame. The operating point and turbulence characteristics are selected to resemble those of modern combustors for space propulsion applications to support the future usage of the devised model in this context.

The existence of the hole on a plate affects the calculation of eddy current problems. Consequently, the accuracy for the prediction of the material properties decreases if the effect of the hole is not taken into account. In this paper, a novel analytical method based on the modified integration range is proposed which can address the presence of the hole. Due to the presence of the hole, the conventional Dodd-Deeds analytical solution cannot be used to calculate the inductance change. Therefore, a revised upper integration limit is introduced to replace the original limit -- ∞ when using the co-axially air-core electromagnetic sensor. With the presence of the hole, the magnitude of the received signal reduces, and the peak frequency feature changes. The analytical method is validated by measured and numerical simulation results. It is found that the upper limit is related to the radius of the open hole. With the new technique, the thickness of sample plates with holes can be estimated based on the peak frequency feature.

Swirling coal flames are encountered in swirl burners of boiler furnaces and cyclone coal combustors, and were studied using large-eddy simulation (LES) in reported references. However, the structures of swirling coal flames remain to be further studied. In the present paper LES of swirling coal flames is conducted for two combustors using two gas combustion models. The statistical simulation results are accessed by experiments and the instantaneous simulation results are used to analyze the specific features of flame structures of four cases. It is seen that different inlet flow conditions and different secondary-air flow rate lead to different flame structures

The mean flow and turbulence statistics of the flow through a simplified urban environment, which is an active research area in order to improve the knowledge of turbulent flow in cities, is investigated. This is useful for civil engineering, pedestrian comfort and for health concerns caused by pollutant spreading. In this work, we provide analysis of the turbulence statistics obtained from well-resolved large-eddy simulations (LES). A detailed analysis of this database reveals the impact of the geometry of the urban array on the flow characteristics and provides for a good description of the turbulent features of the flow within a simplified urban environment. The most prominent features of this complex flow include coherent vortical structures such as the so-called arch vortex, the horseshoe vortex and the roof vortex. These structures of the flow have been identified by an analysis of the turbulence statistics. The influence of the geometry of the urban environment (and particularly the street width and the building height) on the overall flow behavior have also been studied. Finally, the well-resolved LES results were compared with the experimental database from Monnier et al. to discuss differences and similarities between the respective urban configurations.

. The impedance change of an induction coil surrounding a metal tube adapter is investigated by the truncated region eigenfunction expansion (TREE) method. Conventional TREE method is inapplicable to this problem as a consequence of the numerical overflow of eigenfunctions of the air-metal multi-subdomain regions. The difficulty is surmounted by a normalization procedure for the numerical eigenfunctions obtained from the finite element method (FEM). Efficient algorithm is devised by the Clenshaw-Curtis quadrature rule for integrals involving the numerical eigenfunctions. Numerical results of the TREE and FEM simulation coincide very well in all cases, and the efficiency of the proposed method is also confirmed.

Wind energy is a rapidly expanding renewable energy technique. Wind farm developers need to understand the interaction between wind farms and the atmospheric flow over complex terrain. Large-eddy simulations provide valuable data for gaining further insight into the impact of rough topography on wind-farm performance. In this research, we investigate the influence of spatial heterogeneity on wind turbine performance. We conducted numerical simulations of a 12×5 wind-turbine array on various rough topographies. First, we evaluated our LES method through mesh convergence analysis, using mean vertical profiles, vertical friction velocity, and resolved and subgrid-scale kinetic energy. Next, we analyze the effects of surface roughness and dispersive stresses on the performance of fully developed large wind farms. Our results demonstrate that the ground roughness element’s flow resistance boosts large wind-farm power production by almost 68% in fully aerodynamic rough surface compared to flat terrain. Dispersive stress analysis revealed that the primary degree of spatial heterogeneity in the wind farm is in the streamwise direction, which is the “wake-occupied” region, and the relative contribution of dispersive shear stress is almost 45% to the overall drag. We also observed that the power performance of the wind farm in complex terrain outperforms the drag. Our study has implications for improving the design of wind turbines and wind farms in complex terrain to increase their efficiency and energy output.

This article examines the wavelet transforms and the machine learning methods in accelerating large eddy simulations of turbulent flows. An overarching challenge of large eddy simulations is to accurately represent the cascade of kinetic energy across the cut-off length scale. Taylor hypothesis suggests that the energy cascade occurs through the process of vortex stretching. However, Richardson hypothesis assumes the self-amplification of the strain field and a hierarchical break down of large eddies lead the energy cascade. The large eddy simulations typically employs the self-amplification of the strain in formulating subgrid models. However, several studies also proposed subgrid models based on vortex stretching. The wavelet-based coherent vortex simulation of turbulence directly accounts for vortex stretching in overall forward scatter of energy, while allowing local backscatter. The wavelet-based large eddy simulation adapts the grid to capture the creation of small-scale eddies, while adopting subgrid models based on the self-amplification of strain. The advancement of artificial intelligence in turbulence modelling is currently evolving around accelerating the numerical simulations of turbulent flow. However, there is a clear connection between the application of wavelet transfroms and neural networks for directly solving the Navier-Stokes equation, indicating some potential benefits of wavelet methods over the neural networks.

Mesoscale eddies, the weather system of the oceans, although being on the scales of O(20-100km), have a disproportionate role in shaping the mean jets such as the separated Gulf Stream in the North Atlantic Ocean, which is on the scale of O(1000km) in the along-jet direction. With the increase in computational power, we are now able to partially resolve the eddies in basin-scale and global ocean simulations, a model resolution often referred to as mesoscale permitting. It is well known, however, that due to grid-scale numerical viscosity, mesoscale permitting simulations have less energetic eddies and consequently weaker eddy feedback onto the mean flow. In this study, we run a quasi-geostrophic model at mesoscale resolving resolution in a double gyre configuration and formulate a deterministic parametrization for the eddy rectification term of potential vorticity (PV), namely, the eddy PV flux divergence. We have moderate success in reproducing the spatial patterns and magnitude of eddy kinetic and potential energy diagnosed from the model. One novel point about our approach is that we account for non-local eddy feedbacks onto the mean flow by solving the eddy PV equation prognostically in addition to the mean flow. In return, we are able to parametrize the variability in total (mean+eddy) PV at each time step instead of solely the mean PV. A closure for the total PV is beneficial as we are able to account for both the mean state and extreme events.

In this study, the generality and prediction accuracy of a generalised series model for large eddy simulation of premixed and non-premixed turbulent combustion is explored. The model is based on Taylor series expansion of the chemical source term in scalar space and implemented into OpenFOAM. The mathematical model does not depend on the combustion regimes, and has the correct limiting behaviour. The numerical error sources are outlined and analysed. The model is first applied to a piloted methane/air non-premixed jet flame (Sandia Flame D). The statistical (time-averaged and RMS) results agree well with experimental measurements, particularly with regard to the mixture fraction, velocity, temperature, and concentrations of major species CH4, CO2, H2O, and O2. The concentrations of intermediates CO and H2 are over-predicted, due to the limitations of the reduced reaction mechanism employed. Then, a Bunsen piloted flame is simulated. Most of the statistical properties of both reactive species and progress variable are well reproduced. The only major discrepancy evident is in the temperature, which is attributed to the experimental uncertainties of temperature fields in the pilot stream. These findings have demonstrated the model`s generality for both premixed and non-premixed combustion simulation, as well as the accuracy of prediction of reactive species distribution.

When undertaking wind assessment around buildings using large eddy simulation (LES), the implementation of the integral length scale at the inlet for inflow generation is controversial, as real atmospheric length scales require huge computational domains. While length scales significantly influence inflow generation in the domain, their effect on the downstream flow field has not yet, been investigated. In this paper, we validate the effectiveness and accuracy of implementing a reduced turbulence integral length scale for inflow generation in LES results at the rooftop of low-rise buildings and develop a technique to estimate the real local length scales using simulation results. We measure the wind locally and calculate the turbulence length scales from the energy spectrum of the wind data and simulation data. According to these results, there is an excellent agreement between the length scale from simulation and measurement when they are scaled with their corresponding freestream/inlet value. These results indicate that a reduced integral length scale can be safely used for LES to provide a reliable prediction of the energy spectrum as well as the length scales around complex geometries. The simulation results were confidently employed to obtain the best location for a wind turbine installation on low-rise buildings.

Numerical simulations of the sub-critical flow past a circular cylinder at the Reynolds number, Re=130,000 are performed using two numerical platforms, commercial, Ansys Fluent and open source, OpenFoam (finite volume method and large eddy simulation based on a differential equation for the sub-grid kinetic energy). An overview of the available experimental data and similar large eddy simulation studies is presented. A detailed analysis of all accumulated data demonstrates satisfactory agreement between them with a dispersion of the order of 10% (for the main integral flow parameters). A detailed comparison of the results obtained by single- and double-precision numerical methods using Ansys Fluent did not reveal any noticeable discrepancies for the integral and local flow parameters, as well as spectral characteristics. It is shown that the behavior of the dynamic system of the fluid dynamic equations computed with single precision is stable by Lyapunov and does not lead to any loss of accuracy. The reconstructed attractors of the dynamic systems in the phase space are limited by an ellipsoid.

High mountain rainforests are vital in the global energy and carbon cycle. Understanding the exchange of energy and carbon plays an important role to reflect responses to climate change. In this study, an eddy-covariance (EC) measurement system installed in the high Andean Mountains of South Ecuador was used. As EC measurements are affected by heterogeneous topography and vegetation height, the main objective was to estimate the effect of the sloped terrain and the forest on the turbulent energy and carbon fluxes considering the energy balance closure (EBC) and the heat storage. The results showed that the performance of the EBC was generally good and estimated to be 79.5%. This could be improved when the heat storage effect was considered. Based on the variability of the residuals in the diel modifications in the imbalances were highlighted. Particularly, during daytime residuals were largest (56.9 W/m² on average) with a clear overestimation. At nighttime mean imbalances were rather weak (6.5 W/m²) and mostly positive, while strongest underestimations developed in the transition period to morning hours (down to -100 W/m²). With respect to the Monin-Obukhov stability parameter ((z-d)/L) and the friction velocity (u*) it was revealed that largest overestimations evolved in weak unstable and very stable conditions associated with large u* values. In contrast, underestimation was related to very unstable conditions. The estimated carbon fluxes were independently modelled with a non-linear regression using a light-response relationship and reached a good performance (R²=0.52). All fluxes were additionally examined in the annual course to estimate whether both, energy and carbon fluxes, resembled the microclimatological conditions of the study site. This unique study demonstrated that EC measurements provide valuable insights into land surface -atmosphere interactions and contribute to our understanding of energy and carbon exchanges. Moreover, the flux data provide an important basis to validate coupled atmosphere ecosystem models.

: In the automotive industry, high-strength plates have been increasingly used to reduce the vehicle body weight due to stringent regulations on fuel economy and safety. Such plates achieve tensile strengths as high as 1.5 GPa owing to the hot-stamping process. However, these plates suffer from scale generation and surface decarburization because of the surface oxidation of the material. Recently, studies on materials subjected to high-frequency induction heating have reported strength improvement owing to the coil shape and cooling method. This study presents the surface temperature changes, mechanical properties, and microstructures of boron steels subjected to longitudinal heating, which uses eddy currents generated from high-frequency induction heating. The surface temperature data were analyzed under different high-frequency induction heating powers (15, 18, 21, 24, 27, and 30 kW) and distances from the specimen (6, 9, and 15 mm). The mechanical properties and microstructures were also analyzed in different high-frequency induction heating coil directions. The correlations between the high-frequency induction heating power and the distance from the specimen with the maximum tensile strength were determined.

Traditional numerical techniques such as sliding mesh, dynamic grid mesh, and others, have many limitations in dealing with flow simulation with large-scale movement of solid boundaries; which is the case for complex-shaped hydraulic turbomachinery such as propellers, pumps and turbines. The immersed boundary (IB) method provides a new approach to solve the above limitations. Therefore, this study proposes an sharp-interface IB method based on the level-set function that is suitable for simulating the flow through turbomachinery with complex geometries. This method is applied to actual three-dimensional numerical simulations of high-Reynolds number propellers using an in-house computational fluid dynamics solver. The results show that the proposed method can provide comparatively accurate predictions of unsteady load coefficients within the propeller flow passage, and capture the correct propeller wake characteristics as well as the interaction between the propeller wake and free surface. This study aims at providing a theoretical basis and engineering reference for the application of the IB method in engineering numerical simulations.

Despite a rapid rise in NBS development in recent years, the methods for evaluating NBS still have certain gaps. We propose an approach based on a combination of remote sensing data and meteorological variables to reconstruct the spatio-temporal variation of net ecosystem exchange from eddy-covariance stations. Lagrangian particle dispersion model was used for upscaling of satellite images and flux towers. We trained data-driven models based on kernel methods separately for each selected land cover class. The results suggest that the proposed approach to quantifying carbon exchange on a medium-to-large scale by blending eddy covariance flux data with moderate resolution satellite and weather data provides a set of key advantages over previously deployed methods: (1) scalability, achieved via the validation design based on a separate set of eddy covariance stations; (2) high spatial and temporal resolution due to use of the Landsat imagery; (3) robust and accurate predictions due to improved data quality control, advanced machine learning techniques, and rigorous validation. The machine learning models yielded high cross-validation results. Overall we present here globally scaled technology for the land sector based on high resolution remote sensing imagery, meteorological variables, and direct carbon flux measurements of eddy covariance flux stations.

Turbulent flows over and within forest canopies have recently been modeled with success using Large Eddy Simulations (LES). Validation exercises against experimental data suggest that models can be applied with a high degree of confidence for many applications, mechanical and physiological plant/atmosphere interaction analysis, seed or pollen dispersal, wildfire spread and firebrand transport, or investigation of causes of eddy-covariance technique bias. Long distances required for shear-induced turbulence to equilibrate, result in the widespread use of cyclic boundary conditions in LES atmospheric boundary layer studies. Vegetation drag dissipates air momentum in the atmosphere, but equilibrium is often achieved through compensatory momentum source, supplied by macro-scale pressure gradient forcing. Unfortunately, both classical Ekman balance or simple spatially-constant pressure gradient techniques for implementing this forcing have major drawbacks in the context of cyclic boundary conditions for the applications listed above. Among them, it is difficult to specify aspects of the mean velocity profile such as a specific desired wind velocity and direction at a reference height. In the present paper, we propose a new technique for capturing the effects of a large-scale pressure gradient force (LSPGF) that can be used at stand scale and enables simulation of realistic and specifiable wind fields. Several variants of this LSPGF are developed and analyzed here and validated against experimental data. Although this LSPGF technique is developed in the context of HIGRAD/FIRETEC wildfire simulations, LSPGF can be used for any LES wind modeling application aimed at generating detailed stand-scale wind fields with resolved turbulence and shear profiles consistent with vegetation structure in the boundary layer.

The eddy current loss should be optimized to be as less as possible for the stability of permanent magnet in high speed permanent magnet synchronous motor (HSPMSM) rotor and ensure the high efficiency and low temperature of the motor. This paper analyzes the eddy current distribution in rotor, with consideration of the conflict of the thickness of sleeve and diameter of the rotor, calculating the eddy current loss (ECL) and the thermal distribution via Separation of variables method for solving Maxwell's equations with analytical hieratical model of ECL constructed. The optimization result of ECL of the HSPMSM whose power and rated speed is 30kw 48000r/min can be got by multi-objective optimization method, combined weighting coefficient method and traversal algorithm based on chaotic local search particle swarm optimization (CLSPSO), utilizing ECL analytical model and other analytical constraints. Related experiment and measurement has been implemented with new approach of loss separation.

Nuclear Magnetic Resonance (NMR) and its various forms are extensively utilized in both research and clinical settings to analyse molecules. NMR data pre-processing, while essential to NMR data analysis, can be uniquely complex. Despite the availability of software tools, understanding these processes can be challenging, complicating the selection of appropriate pre-processing steps. In this review, we elucidate pre-processing steps in the time domain from a mathematical and statistical perspective, explaining direct current offset removal, eddy current correction, shift and linear prediction, weighting, zero filling, and domain transformation using plain language and fundamental mathematical formulas. Our objective is to clarify these processes in simple terms and provide general guidance.

The aim of work is to derive an explicit expression for a function of vertical mixing induced by wind-waves. To this end, in the Navier-Stokes equations, a current is decomposed into four constituents: the mean flow, the wave-orbital motion, the wave-induced turbulent and the background turbulent currents. This decomposition allows separating the wave-induced Reynolds stress, Rw, from the background one, Rb. To make a statistical closure for Rw, the Prandtl approach for the background turbulent fluctuations is used that results in an implicit expression for the wave-induced vertical mixing function, Bv. Expression for Bv is specified based on the author’s results for the eddy viscosity found earlier in the frame of the three-layer concept for a wavy air–sea interface, used for modelling wind-drift currents [1]. Finally, the explicit parameterization for Bv(a, u*, z) is found as a linear function in both the wave amplitude at depth z, a(z), and the friction velocity in the air, u*. The linear dependence of function Bv(a) on the wave amplitude provides the enhanced vertical mixing induced by wind–waves in comparison with function Bv(a) having the cubic dependence found in [2], as far as the wind-wave amplitude a(z) decays exponentially with depth.

The present work is a numerical simulation of velocity and mixture fraction fields in a turbulent non-reaction Propane jet flow issuing into parallel co-flowing air, in isothermal condition which has been experimentally described here: http://www.sandia.gov/TNF/DataArch/ProJet.html. The objective is a better understanding of the flow structure and mixing process, a situation in which there is no chemical interaction and heat-transfer. The two-equation Realizable k-ε eddy viscosity turbulence model has been used to simulate the turbulent flow field on a 2D plane (i.e., on a 5-degree sector of the experimental domain), because Realizable k-ε more accurately predicts the spreading rate of both planar and round jets and presents the best proficiency in comparison with all versions of the k-ε models. Afterward, axial and radial profiles of mean velocities, turbulence energy, mean mixture fraction, the mixture fraction half radius (Lf), and the mass flux diagram have been numerically elicited for grid independent mesh (mesh B) and compared with corresponding experimental data to assess the numerical model. To obtain turbulence kinetic energy, k, due to lack of w' in the experimental data, the assumption of w'=v' seems to be valid. Simulations have demonstrated that mean mixture fraction (at radial profiles at locations x/D=0, 4, 15, 30 and 50) and its half radius, Lf, which characterizes jet width expanse, are prominently well-captured; Moreover, mean vertical velocity fields (in X-direction: Umean) have revealed less accuracy but still conspicuously well-captured. However, mean vertical velocity fields (in Y-direction: Vmean) have disclosed less resolution; Likewise, turbulence kinetic energy, k, have manifested moderate accuracy (between Umean and Vmean). It should be noticed that although numerical results for absolute pressure, p, have been obtained on aforementioned sections, there were no experimental data to compare with them. Thus, the corresponding numerical data has not been demonstrated in this study.

The effects of the coolant pulsation and the film cooling with the plasma aerodynamic actuation (PAA) could be herein explored via large eddy simulations. The electrohydrodynamic force derived from the PAA was solved through the phenomenological plasma model. The Strouhal number of the sinusoidal coolant pulsation and the averaged pulsation blowing ratio were 0.25 and 1.0, respectively. Comprehensive analyses were carried out on the time-averaged flow fields, and the results revealed that the pulsed cooling jet might cause a deeper penetration into the crossflow, and this phenomenon could be remarkably mitigated by the downward force of the PAA. Furthermore, because of the coolant pulsation, the counter rotating vortex pair (CRVP) expanded to lifted off the wall less well, and the PAA lessened the negative lift-off impact and entrapment from the CRVP. Then, the spatial-temporal development of the coherent structures was figured out by the alterations of the centerline temperature, reflecting formation of the intermittent coherent structures rather than hairpin vortices due to the coolant pulsation, and their size and upcast behaviors were reduced by the PAA. Consequently, there was a basic suppression of the turbulent mixing of the cooling fluid with the crossflow. Finally, the three-dimensional streamlines confirmed that the coherent structure dynamic behaviors were significantly regulated by the PAA for alleviating the adverse influences of the coolant pulsation. In summary, through controlling how the dominating coherent structures evolve throughout time and space, the PAA can effectively increase the pulsed film cooling efficiency.

Transport electrification is essential for reducing CO2 emissions, and technologies such as hybrid and range-extended electric vehicles will play a crucial transitional role. Such vehicles employ an internal combustion engine for on-board chemical energy conversion. The Wankel rotary engine should be an excellent candidate for this purpose, offering high power-to-weight ratio, simplicity, compactness, perfect balance, and low cost. Until recently, however, it has not been in production in the automotive market, due, in part, to relatively low combustion efficiency and high fuel consumption and unburnt hydrocarbon emissions, which can be traced to constraints on flame speed, an elongated combustion chamber, and relatively low compression ratios. This work uses large-eddy simulation to study the in-chamber flow in a peripherally ported 225cc Wankel rotary engine, providing insight into these limitations. Flow structures created during the intake phase play a key role in turbulence production but presence of the pinch point inherent to Wankel engine combustion chambers inhibits flame propagation. Two efficiency-enhancement technologies are introduced as disruptive solutions: (i) pre-chamber jet ignition and (ii) the two-stage rotary engine. These concepts overcome the traditional efficiency limitations and show that Wankel rotary engine design can be further enhanced for its role as a range extender in electrified vehicles.

For both conventional and renewable energy conversion processes, computational fluid dynamics (CFD) has been used to address more energy-related challenges in recent decades. Using CFD to investigate vertical-axis wind turbines has become more common in recent years. The main goals of this application have been to more accurately predict the turbine's performance and to comprehend the complicated nature of the complex turbulent flow. The Vertical Axis Wind Turbine (VAWT) simulation for energy-generating applications has several intricate components. One of them is being aware of the chaotic flow that occurs during the first stages of the starting process, which influences overall effectiveness. In this article, the performance of the wind turbine was increased by using a passive flow control approach. The numerical research was carried out using Large Eddy Simulation for four alternative tip speed ratios in both cases, the classic and the optimized, equipped with a vortex trap on the exterior of the blades. The power and torque coefficient variations, as well as the velocity magnitude contours, show that the starting process may begin with a significant improvement in efficiency when flow control is used.

Computational rhinology is a specialized branch of biomechanics leveraging engineering techniques for mathematical modelling and simulation to complement the medical field of rhinology. Computational rhinology has already contributed significantly to advancing our understanding of the nasal function, including airflow patterns, mucosal cooling, particle deposition, and drug delivery, and is foreseen as a crucial element in e.g. the development of virtual surgery as a clinical, patient-specific decision support tool. The current paper delves into the field of computational rhinology from a nasal airflow perspective, highlighting the use of computational fluid dynamics to enhance diagnostics and treatment of breathing disorders. The paper consists of three distinct parts – an introduction to and review of the field of computational rhinology, a review of published literature on in vitro and in silico studies of nasal airflow, and the presentation and analysis of previously unpublished high fidelity CFD simulation data of in silico rhinomanometry. While the two first parts of the paper summarize the current status and challenges in the application of computational tools in rhinology, the last part addresses the gross disagreement commonly observed when comparing in silico and in vivo rhinomanometry results. It is concluded that this discrepancy cannot readily be explained by CFD model deficiencies caused by poor choice of turbulence model, insufficient spatial or temporal resolution, or neglecting transient effects. Hence, alternative explanations such as nasal cavity compliance or drag effects due to nasal hair should be investigated.

An approach based on the OpenFOAM library has been developed to solve a high-speed, multicomponent mixture of a reacting, compressible flow. This work presents comprehensive validation of the newly developed solver called compressibleCentralReactingFoam with different supersonic flows, including shocks, expansion waves, and turbulence-combustion interaction. The comparisons of the simulation results with experimental and computational data confirm the fidelity of this solver for problems involving multicomponent high-speed reactive flows. The gas dynamics of turbulence-chemistry interaction are modeled using a partially-stirred reactor formulation and provide promising results to better understand the complex physics involved in supersonic combustors. A time-scale analysis based on local Damköhler numbers reveals different regimes of turbulent combustion. In the core of the jet flow, the Damköhler number is relatively high indicating that the reaction time scale is smaller than the turbulent mixing time scale. This means that the combustion is controlled by turbulent mixing. In the shear layer, where the heat release rate and the scalar dissipation rate have the highest value, the flame is stabilized due to finite rate chemistry with small Damköhler numbers and a high fraction of fine structure. This solver allows three-dimensional gas dynamic simulation of high-speed multicomponent reactive flows relevant to practical combustion applications.

The water cushion depth of stilling basin with shallow-water cushion is a key factor that affects the flow regime of hydraulic jump in the basin. However, the specific depth at which the water cushion is considered as “shallow” has not be stated clearly for now, and only conceptual description is provided. This paper attempts to specify the best water cushion depth based on the flow regime of hydraulic jump and underflow speed; namely, in case of critical hydraulic jump in the basin, the best water cushion depth is located where the minimum distance to the bottom plate of the stilling basin is 1/5~1/4 of the water cushion depth. The theoretical analysis indicates, at different inclinations of discharge chute (θ) and depth ratios of inlet (m), instead of monotonic change, the Froude number (Fr) at inlet of the stilling basin with shallow-water cushion firstly reduces and then increases as the flow velocity at discharge chute inlet (V) increases; the parameters of inflection point (critical flow velocity and critical Fr) increase as the inclinations of discharge chute (θ) and depth ratios of inlet (m) increase. Such regularity is the theoretical basis for selecting representative study cases. The reliability of the large eddy simulation calculation results are verified by a model test; in the paper, 30 cases including five different Froude numbers and six shallow-water cushion depths are selected, for calculating the hydraulic factors such as flow profile, flow regime and flow velocity in the stilling basin with shallow-water cushion; and the varying pattern between the best depth of stilling basin with shallow-water cushion (depth-to-length ratio) and the inflow Froude number is obtained which indicates that the best depth of stilling basin with shallow-water cushion varies little as the change of the Froude number before reaching the critical Froude number; however, the best depth-to-length ratio of stilling basin with shallow-water cushion increases as the Froude number increases after the critical Froude number is reached. The study results in this paper are of reference significance to design and calculation of the stilling basin with shallow-water cushion.

Evapotranspiration (ET) is a major component of water budget in Bottomland Hardwood Forests (BHF) and is driven by a complex intertwined suite of meteorological variables. The understanding of these interdependencies leading to seasonal variations in ET is crucial in better informing water resource management in the region. We used a structural equation modeling approach to analyze the drivers of ET using Eddy covariance water and heat flux data collected from a BHF, located in the Russel Sage Wildlife Management Area (RSWMA). It consists of mature closed-canopy deciduous hardwood trees with an average canopy height of 27 m. Factor analysis was used to characterize the shared variance among drivers and path analysis was used to quantify independent contributions of individual drivers. In our results, ET and net radiation (Rn) showed similar variability patterns with Vapor Pressure Deficit (VPD) and temperature in Spring, Summer, and Autumn seasons while they differ in Winter season. The path analysis shows that Rn has the strongest influence on ET variations via direct and indirect pathways. In the deciduous forests like BHFs, our results suggest ET is more energy-dependent during growing season (Spring and Summer) and early non-growing season (Autumn) while more temperature-dependent during winter season.

Mathematical modeling is the most important tool for constructing the theory of different kinds of two-phase flows. The review is devoted to the analysis of the advent of the approach of mathematical modeling of two-phase flows, where solid particles act mainly as the dispersed phase. The main problems and features of the study of gas-solid two-phase flows are excluded. The main characteristics of gas flows with solid particles are given and the classification of two-phase flows is developed on their basis. The Lagrangian and Euler approaches to modeling the motion of a dispersed phase (particles) are described. Much attention is paid to the consideration of numerical simulation methods that describe descriptions of turbulent gas flow at different hierarchical levels (RANS, LES, DNS), different levels of description of interphase interactions (one-way coupling (OWC), two-way coupling (TWC) and four-way coupling (FWC)), as well as at different levels of interface resolution (partial-point (PP) and particle-resolved (PR)). Excluded are examples of some studies carried out on the basis of the identified approaches, as well as their exclusion for mathematical modeling of various classes of gas-solid two-phase flows.

This paper reviews how dynamically adaptive wavelet methods can be designed to simulate atmosphere and ocean dynamics in both flat and spherical geometries. We highlight the special features that these models must have in order to be valid for climate modelling applications. These include exact mass conservation and various mimetic properties that ensure the solutions remain physically realistic, even in the under-resolved conditions typical of climate models. Particular attention is paid to the implementation of complex topography in adaptive models. Using \textsc{wavetrisk} as an example, we explain in detail how to build a semi-realistic global atmosphere or ocean model of interest to the geophysical community. We end with a discussion of the challenges that remain to developing a realistic dynamically adaptive atmosphere or ocean climate models. These include scale-aware subgrid scale parameterizations of physical processes, such as clouds. Although we focus on adaptive wavelet methods, many of the topics we discuss are relevant for adaptive mesh refinement (AMR).

The general-purpose Compact Muon Solenoid (CMS) detector at the Large Hadron Collider (LHC) at CERN studies a production of new particles in the proton-proton collisions at the LHC center of mass energy 13.6 TeV. The detector includes the magnet based on the 6 m diameter superconducting solenoid coil operating with the current of 18.164 kA. This current creates a central magnetic flux density of 3.8 T that allows to measure momenta of the produced in collisions charged particles with tracking and muon subdetectors at a high precision. The CMS magnet contains a 10,000 ton flux-return yoke made of construction steel that bends muons in steel blocks magnetized with the solenoid returned magnetic flux. To reconstruct the muon trajectories, and thus, to measure the muon momenta, the drift tube and cathode strip chambers are located between the layers of the steel blocks and serve for this purpose. To describe the distribution of the magnetic flux in the magnet yoke layers, a three-dimensional computer model of the CMS magnet is used. To prove the calculations, the special measurements are performed with the flux loops wound in 22 cross-sections of the flux-return yoke blocks. The measured voltages induced in the flux loops during the CMS magnet ramp ups and downs, as well as during the superconducting coil fast discharges with the 190 s time constant, are integrated over time to obtain the initial magnetic flux densities in the flux loop cross-sections. The measurements obtained during the seven standard ramp downs of the magnet have been analyzed in 2018. From that time three fast discharges are occurred during the standard ramp downs of the magnet. This allows to single out the contributions of the eddy currents, induced in steel, to the flux loop voltages registered during the fast discharges of the coil. Accounting for these contributions to the flux loop measurements during manually triggered fast discharges in 2006 allows to perform the validation of the CMS magnet computer model with a better precision. The technique of the flux loop measurements, and the obtained results are presented and discussed.

This paper attempts to develop a Dynamic Mode Decomposition (DMD) based Reduced Order Model (ROMs) which can quickly but accurately predict the forces and moments experienced by a road vehicle such that they be used by an on-board controller to determine the vehicle's trajectory. DMD can linearize a large dataset of high-dimensional measurements by decomposing them into low-dimensional coherent structures and associated time-dynamics. This ROM then also can be applied to predict the future state of the fluid flow. Existing literature on DMD is limited to low Reynolds number applications. This paper presents DMD analyses of the flow around an idealized road vehicle, called the Ahmed body, at a Reynolds number of 2.7E6. The high dimensional dataset used in this paper was collected from a computational fluid dynamics (CFD) simulation performed using Menter's shear stress transport (SST) turbulence model within the context of Improved Delayed Detached Eddy Simulations (IDDES). The DMD algorithm, as available in literature, was found to suffer nonphysical dampening of the medium-to-high frequency modes. Enhancements to the existing algorithm were explored, and a modified DMD approach is presented in this paper which includes: (a) a requirement of higher sampling rate to obtain higher resolution of data, and (b) a custom filtration process to remove spurious modes. The modified DMD algorithm thus developed was applied to the high Reynolds number, separation dominated flow past the idealized ground vehicle. The effectiveness of the modified algorithm was tested by comparing future predictions of force and moment coefficients as predicted by the DMD-based ROM to the reference CFD simulation data, and found to offer significant improvement.

Forest carbon sequestration is a widely accepted natural climate solution, however, methods to determine net carbon offsets are limited to commercial carbon proxies and CO2 eddy covariance research. Non-CO2 greenhouse gases (GHG) (e.g., CH4, N2O) receive less attention in the context of forests, in part, due to emphasis on CO2 and the operational requirements and cost for three-gas eddy covariance platforms. In this study, Howland forest flux tower (CO2, CH4) and soil flux data (CO2, CH4, N2O), representing net emission reductions, are linked to their respective social costs to estimate commercial revenue if sold as a GHG social cost forest offset product (GHG-SCF). Estimated annual revenue for GHG-SCF products, applicable to realization of a Green New Deal, range from 120,000 covering the site area of 557 acres in 2021, to 12,000,000 for extrapolation to 40,000 acres in 2040, assuming a 3% discount rate. The Howland Forest CO2 flux record for two adjacent towers is compared to California Air Resources Board forest carbon proxy data for compliance sequestration offsets, the only project site where these approaches overlap. Overcrediting, incomplete carbon accounting with annual errors of up to 2,256%, inadequate third-party verification, and limited application to non-CO2 GHG’s are established. In contrast, direct measurement of one or more GHG’s offers new forest products and revenue incentives to restore and conserve forests worldwide.