<?xml version="1.0" encoding="utf-8" standalone="yes"?><rss version="2.0" xmlns:atom="http://www.w3.org/2005/Atom"><channel><title>Physics-Based Modeling | Yonghui Liu</title><link>https://lyh360121.github.io/tags/physics-based-modeling/</link><atom:link href="https://lyh360121.github.io/tags/physics-based-modeling/index.xml" rel="self" type="application/rss+xml"/><description>Physics-Based Modeling</description><generator>HugoBlox Kit (https://hugoblox.com)</generator><language>en-us</language><lastBuildDate>Sun, 01 May 2022 00:00:00 +0000</lastBuildDate><image><url>https://lyh360121.github.io/media/sharing.png</url><title>Physics-Based Modeling</title><link>https://lyh360121.github.io/tags/physics-based-modeling/</link></image><item><title>Analytical modeling of resonant converters in frequency domain for wireless power transfer: Continuous current mode (CCM) operation</title><link>https://lyh360121.github.io/publications/2022-zeng-analytical-modeling-of/</link><pubDate>Sun, 01 May 2022 00:00:00 +0000</pubDate><guid>https://lyh360121.github.io/publications/2022-zeng-analytical-modeling-of/</guid><description>&lt;p&gt;Wireless power transfer (WPT) systems move energy across an air gap using a &lt;strong&gt;resonant converter&lt;/strong&gt;: an H-bridge inverter drives a resonant tank (inductors and capacitors) that couples through two loosely-coupled coils, and the received AC is turned back into DC by a &lt;strong&gt;diode-bridge rectifier&lt;/strong&gt;. Getting the steady-state behavior of this converter right — the current and voltage waveforms it settles into — is what lets engineers size components and predict output power before building hardware. The problem is that the diode-bridge rectifier is genuinely non-linear (it switches on and off depending on the polarity of the incoming current), so the waveforms it produces are distorted, not pure sine waves. The industry-standard shortcut, &lt;strong&gt;Fundamental Harmonic Approximation (FHA)&lt;/strong&gt;, ignores that distortion and only tracks the fundamental-frequency component. It is fast, but the paper shows it breaks down once the converter is driven away from its resonant point — and that emerging, higher-power-density WPT designs, which deliberately inject extra harmonics, push it into exactly that regime.&lt;/p&gt;
&lt;p&gt;Prior frequency-domain models that do account for higher harmonics beyond FHA exist, but each one leans on a simplifying assumption to stay solvable: some fix the source driving the resonant network to an idealized sine wave, while others allow an arbitrary source but force the rectifier&amp;rsquo;s DC-side output to be a fixed voltage, a square wave, a pre-specified voltage/current source, or — at best — a purely inductive load. Those assumptions break down for more general, realistic WPT loads, such as the paper&amp;rsquo;s own test hardware, where an output LC filter sits between the rectifier and a resistive load — a combination none of the earlier fixed-output assumptions can fully capture. This paper&amp;rsquo;s contribution is a model that keeps arbitrary input and output configurations while cutting the model down to a &lt;strong&gt;single unknown variable&lt;/strong&gt; — an &amp;ldquo;offset angle&amp;rdquo; describing where the rectifier&amp;rsquo;s input current crosses zero relative to the inverter&amp;rsquo;s output voltage — solved through nothing more than matrix operations (no integration or differentiation needed).&lt;/p&gt;
&lt;p&gt;
&lt;figure id="figure-fig-1-from-the-paper-overall-schematic-of-the-conventional-inductive-wireless-charger-showing-the-h-bridge-inverter-resonant-network-and-diode-bridge-rectifier-that-the-model-captures-end-to-end"&gt;
&lt;div class="flex justify-center "&gt;
&lt;div class="w-full" &gt;
&lt;img alt="Circuit topology of the wireless power transfer system modeled in the paper: an H-bridge inverter driving a resonant network (Cp, Lkp, magnetizing inductance LM, coupled coils N1:N2, Lks, Cs) that feeds a diode-bridge rectifier and output filter/load."
srcset="https://lyh360121.github.io/publications/2022-zeng-analytical-modeling-of/fig-circuit-architecture_hu_733d7f5dc4f8e69.webp 320w, https://lyh360121.github.io/publications/2022-zeng-analytical-modeling-of/fig-circuit-architecture_hu_8e3517799c9abd5f.webp 480w, https://lyh360121.github.io/publications/2022-zeng-analytical-modeling-of/fig-circuit-architecture_hu_fad34934be74eae2.webp 760w"
sizes="(max-width: 480px) 100vw, (max-width: 768px) 90vw, (max-width: 1024px) 80vw, 760px"
src="https://lyh360121.github.io/publications/2022-zeng-analytical-modeling-of/fig-circuit-architecture_hu_733d7f5dc4f8e69.webp"
width="760"
height="233"
loading="lazy" data-zoomable /&gt;&lt;/div&gt;
&lt;/div&gt;&lt;figcaption&gt;
Fig. 1 from the paper: overall schematic of the conventional inductive wireless charger, showing the H-bridge inverter, resonant network, and diode-bridge rectifier that the model captures end to end.
&lt;/figcaption&gt;&lt;/figure&gt;
&lt;/p&gt;
&lt;p&gt;The method works in three steps. First, the non-linear diode-bridge rectifier is captured with a harmonic (Fourier-domain) admittance/impedance matrix that relates every harmonic of the AC-side current and voltage to the DC-side output — including a later refinement that adds the diodes&amp;rsquo; own forward-voltage drop, which matters most for low-output-voltage or high-frequency designs. Second, the passive resonant network (whatever combination of L&amp;rsquo;s and C&amp;rsquo;s the designer chose) is folded in as an ordinary linear two-port network. Combining the two collapses the whole system into one equation with a single unknown: the offset angle. That angle is found by an iterative one-dimensional search — starting from an FHA-based initial guess and refined with a parabolic-interpolation search — which the paper shows converges in under ten iterations, versus gradient-descent methods that can fail outright because the underlying error function isn&amp;rsquo;t smooth.&lt;/p&gt;
&lt;p&gt;To validate the model, the authors built two hardware WPT prototypes — one with series-series (SS) compensation, one with LCC compensation, both based on a dual-coil Qi-style transmitter/receiver — and compared the model&amp;rsquo;s predicted waveforms against PSIM time-domain simulation and bench measurements.&lt;/p&gt;
&lt;p&gt;
&lt;figure id="figure-fig-10-from-the-paper-i_ac-and-v_h-waveforms-from-experiment-time-domain-simulation-the-proposed-calculation-and-fha-for-the-ss-topology-n--50-v_in--1152-v-r_o--10-ω"&gt;
&lt;div class="flex justify-center "&gt;
&lt;div class="w-full" &gt;
&lt;img alt="Measured, simulated, and calculated AC current (i_ac) and H-bridge voltage (v_H) waveforms for the SS-compensated prototype, comparing the proposed method against the conventional FHA calculation."
srcset="https://lyh360121.github.io/publications/2022-zeng-analytical-modeling-of/fig-waveform-validation_hu_a62c96e330ff7f1e.webp 320w, https://lyh360121.github.io/publications/2022-zeng-analytical-modeling-of/fig-waveform-validation_hu_e201b221ded4d41c.webp 480w, https://lyh360121.github.io/publications/2022-zeng-analytical-modeling-of/fig-waveform-validation_hu_dc2d846e7cecae0.webp 760w"
sizes="(max-width: 480px) 100vw, (max-width: 768px) 90vw, (max-width: 1024px) 80vw, 760px"
src="https://lyh360121.github.io/publications/2022-zeng-analytical-modeling-of/fig-waveform-validation_hu_a62c96e330ff7f1e.webp"
width="760"
height="291"
loading="lazy" data-zoomable /&gt;&lt;/div&gt;
&lt;/div&gt;&lt;figcaption&gt;
Fig. 10 from the paper: i_ac and v_H waveforms from experiment, time-domain simulation, the proposed calculation, and FHA, for the SS topology (N = 50, V_in = 11.52 V, R_o = 10 Ω).
&lt;/figcaption&gt;&lt;/figure&gt;
&lt;/p&gt;
&lt;p&gt;The proposed model tracks the measured and simulated waveforms closely, including the current distortion that FHA smooths away — FHA&amp;rsquo;s calculated peak current comes out roughly 16% above the experimental value. On output-power estimation, FHA is off by about 18% for the SS prototype at full load (about 15% at half load), while the proposed method lands within roughly 1-2% of the measured output power at both load points (the paper&amp;rsquo;s headline &amp;ldquo;10% improvement&amp;rdquo; figure, from the abstract). The match is noticeably looser for the LCC prototype — about 11% output-power error at full load, versus about 3% at half load — which the authors attribute mainly to the small-signal LCR-meter measurement of the coils (taken at only 500 mV excitation) not fully capturing the coupled inductor&amp;rsquo;s real, non-linear behavior: its effective AC resistance turns out to be roughly five times higher at the third harmonic than at the fundamental frequency, and its behavior also shifts with current level, rather than reflecting a limitation of the frequency-domain method itself. The authors are explicit that the model is currently restricted to continuous current mode (CCM) — it does not yet cover the discontinuous-current regime that shows up under light loads — and that resonant networks with richer harmonic content (like LCC) need a higher-order harmonic matrix, which increases computation cost and is left for future work.&lt;/p&gt;</description></item><item><title>Unified Model of Sediment Transport Threshold and Rate Across Weak and Intense Subaqueous Bedload, Windblown Sand, and Windblown Snow</title><link>https://lyh360121.github.io/publications/2021-pahtz-unified-model-of/</link><pubDate>Thu, 01 Apr 2021 00:00:00 +0000</pubDate><guid>https://lyh360121.github.io/publications/2021-pahtz-unified-model-of/</guid><description>&lt;p&gt;Grains get pushed along a riverbed, blown across a desert dune, or carried over a snowfield by very different physical processes — yet in each case, a grain that isn&amp;rsquo;t fully suspended in the fluid moves the same basic way: it hops, rolls, or slides along the bed in a sequence of collisions. This is called nonsuspended sediment transport (NST), and it shapes sedimentary landscapes on Earth and other planetary bodies alike. In reality NST is messy: turbulence, an uneven bed surface, and grain-to-grain variability make individual grain trajectories highly erratic. This paper asks whether, despite that mess, one single deterministic model can still predict &lt;em&gt;both&lt;/em&gt; when transport starts (the threshold) and &lt;em&gt;how much&lt;/em&gt; material moves once it does (the rate) — across fluvial settings (grains driven by a stream of liquid) and aeolian settings (grains driven by wind) at once.&lt;/p&gt;
&lt;p&gt;The model&amp;rsquo;s key idea is to represent the entire population of grains — including ones that mostly roll or slide rather than hop — by a single, repeating (periodic) saltation trajectory: a grain leaves the bed, follows a ballistic-like path shaped by the mean turbulent flow profile above the bed, and comes back down, where an empirical &amp;ldquo;rebound law&amp;rdquo; (fit to grain-bed collision experiments and simulations) sets its average rebound velocity as a function of how it hit the bed. The transport threshold is then defined as the weakest flow for which such a self-sustaining periodic trajectory can still exist, with an extra criterion checking whether a grain has enough energy to roll itself out of the most stable pockets of the bed surface.&lt;/p&gt;
&lt;p&gt;
&lt;figure id="figure-the-model-reduces-the-entire-fluctuating-grain-motion-to-a-single-periodic-bouncing-trajectory-driven-by-the-mean-flow-profile-above-the-bed-figure-1"&gt;
&lt;div class="flex justify-center "&gt;
&lt;div class="w-full" &gt;
&lt;img alt="Sketch of the model&amp;rsquo;s core mechanism: a grain bounces in a periodic saltation trajectory just above a granular bed, driven by the mean streamwise flow velocity profile u_x(z); the bed&amp;rsquo;s virtual zero level and crest level are marked below."
srcset="https://lyh360121.github.io/publications/2021-pahtz-unified-model-of/fig-saltation-mechanism_hu_edefbb4515461276.webp 320w, https://lyh360121.github.io/publications/2021-pahtz-unified-model-of/fig-saltation-mechanism_hu_8d4f2cca1f7230ef.webp 480w, https://lyh360121.github.io/publications/2021-pahtz-unified-model-of/fig-saltation-mechanism_hu_7c86099655d1bd7a.webp 760w"
sizes="(max-width: 480px) 100vw, (max-width: 768px) 90vw, (max-width: 1024px) 80vw, 760px"
src="https://lyh360121.github.io/publications/2021-pahtz-unified-model-of/fig-saltation-mechanism_hu_edefbb4515461276.webp"
width="760"
height="630"
loading="lazy" data-zoomable /&gt;&lt;/div&gt;
&lt;/div&gt;&lt;figcaption&gt;
The model reduces the entire fluctuating grain motion to a single periodic bouncing trajectory driven by the mean flow profile above the bed (Figure 1).
&lt;/figcaption&gt;&lt;/figure&gt;
&lt;/p&gt;
&lt;p&gt;What makes the result notable is that none of the model&amp;rsquo;s parameters were tuned to the data it is then tested against — they were fixed beforehand from independent experiments and simulations (grain-bed rebound measurements, earlier discrete-element-method simulations of transport, and separate studies of bed yielding). With those parameters locked in, the model reproduces laboratory measurements and field data of both the transport threshold and the transport rate to within about a factor of two, across density ratios between grain and fluid spanning roughly 2.65 (mineral grains in water) up to around 2000 (mineral grains or snow in air) — that is, across subaqueous bedload, windblown sand, and windblown snow together.&lt;/p&gt;
&lt;p&gt;
&lt;figure id="figure-left-predicted-transport-threshold-curves-versus-galileo-number-compared-with-measured-threshold-data-right-predicted-versus-measured-dimensionless-transport-rate-both-panels-combine-fluvial-mineral-aeolian-mineral-and-aeolian-snow-data-sets-figure-5c-d"&gt;
&lt;div class="flex justify-center "&gt;
&lt;div class="w-full" &gt;
&lt;img alt="Two log-log plots: on the left, the model&amp;rsquo;s predicted transport-threshold curves (solid lines, one per density ratio) plotted against experimental threshold data points across a range of Galileo numbers; on the right, the model&amp;rsquo;s predicted dimensionless transport rate plotted directly against the measured transport rate, with points clustering close to the 1:1 line — both panels showing separate data series for fluvial transport of minerals, aeolian transport of minerals, and aeolian transport of snow."
srcset="https://lyh360121.github.io/publications/2021-pahtz-unified-model-of/fig-model-vs-data_hu_356605e3e161a2c.webp 320w, https://lyh360121.github.io/publications/2021-pahtz-unified-model-of/fig-model-vs-data_hu_9ce2d7cbb2ea63d0.webp 480w, https://lyh360121.github.io/publications/2021-pahtz-unified-model-of/fig-model-vs-data_hu_c0c4aacec1f44819.webp 760w"
sizes="(max-width: 480px) 100vw, (max-width: 768px) 90vw, (max-width: 1024px) 80vw, 760px"
src="https://lyh360121.github.io/publications/2021-pahtz-unified-model-of/fig-model-vs-data_hu_356605e3e161a2c.webp"
width="760"
height="252"
loading="lazy" data-zoomable /&gt;&lt;/div&gt;
&lt;/div&gt;&lt;figcaption&gt;
Left: predicted transport-threshold curves versus Galileo number compared with measured threshold data. Right: predicted versus measured dimensionless transport rate. Both panels combine fluvial mineral, aeolian mineral, and aeolian snow data sets (Figure 5c-d).
&lt;/figcaption&gt;&lt;/figure&gt;
&lt;/p&gt;
&lt;p&gt;A striking side finding concerns cohesion. Windblown sand and snow grains are known to stick together via cohesive bonds, yet the model captures their behavior without ever explicitly representing cohesion. This lines up with independent DEM simulations showing that equilibrium aeolian transport is nearly unaffected by how strong those cohesive bonds are — a fairly counterintuitive result that this model helps explain: in the model&amp;rsquo;s picture, cohesion only matters while a grain is in contact with the bed, and for grains that are mostly saltating (as opposed to rolling), those contacts are so brief that cohesion barely gets a chance to act. The model also yields a simple, physically grounded way to tell apart &amp;ldquo;bedload&amp;rdquo; (where a significant fraction of grains roll or slide in sustained contact with the bed) from &amp;ldquo;saltation&amp;rdquo; (where that fraction is negligible) — consistent with an independent criterion from the authors&amp;rsquo; earlier work — and it further separates viscous from turbulent transport regimes based on how the grain&amp;rsquo;s hop height compares to the flow&amp;rsquo;s viscous sublayer thickness.&lt;/p&gt;</description></item></channel></rss>