The invention relates to the field of oscillatory resonant electromagnetic modes, and in particular to oscillatory resonant electromagnetic modes, with localized slowly evanescent field patterns, for wireless non-radiative energy transfer.
In the early days of electromagnetism, before the electrical-wire grid was deployed, serious interest and effort was devoted towards the development of schemes to transport energy over long distances wirelessly, without any carrier medium. These efforts appear to have met with little, if any, success. Radiative modes of omni-directional antennas, which work very well for information transfer, are not suitable for such energy transfer, because a vast majority of energy is wasted into free space. Directed radiation modes, using lasers or highly-directional antennas, can be efficiently used for energy transfer, even for long distances (transfer distance LTRANS>>LDEV, where LDEV is the characteristic size of the device), but require existence of an uninterruptible line-of-sight and a complicated tracking system in the case of mobile objects.
Rapid development of autonomous electronics of recent years (e.g. laptops, cell-phones, house-hold robots, that all typically rely on chemical energy storage) justifies revisiting investigation of this issue. Today, the existing electrical-wire grid carries energy almost everywhere; even a medium-range wireless non-radiative energy transfer would be quite useful. One scheme currently used for some important applications relies on induction, but it is restricted to very close-range (LTRANS<<LDEV) energy transfers.
According to one aspect of the invention, there is provided an electromagnetic energy transfer device. The electromagnetic energy transfer device includes a first resonator structure receiving energy from an external power supply. The first resonator structure has a first Q-factor. A second resonator structure is positioned distal from the first resonator structure, and supplies useful working power to an external load. The second resonator structure has a second Q-factor. The distance between the two resonators can be larger than the characteristic size of each resonator. Non-radiative energy transfer between the first resonator structure and the second resonator structure is mediated through coupling of their resonant-field evanescent tails.
According to another aspect of the invention, there is provided a method of transferring electromagnetic energy. The method includes providing a first resonator structure receiving energy from an external power supply. The first resonator structure has a first Q-factor. Also, the method includes a second resonator structure being positioned distal from the first resonator structure, and supplying useful working power to an external load. The second resonator structure has a second Q-factor. The distance between the two resonators can be larger than the characteristic size of each resonator. Furthermore, the method includes transferring non-radiative energy between the first resonator structure and the second resonator structure through coupling of their resonant-field evanescent tails.
In another aspect, a method of transferring energy is disclosed including the steps of providing a first resonator structure receiving energy from an external power supply, said first resonator structure having a first resonant frequency ω1, and a first Q-factor Q1, and characteristic size L1. Providing a second resonator structure being positioned distal from said first resonator structure, at closest distance D, said second resonator structure having a second resonant frequency ω2, and a second Q-factor Q2, and characteristic size L2, where the two said frequencies of and ω2 are close to within the narrower of the two resonance widths Γ1, and Γ2, and transferring energy non-radiatively between said first resonator structure and said second resonator structure, said energy transfer being mediated through coupling of their resonant-field evanescent tails, and the rate of energy transfer between said first resonator and said second resonator being denoted by κ, where non-radiative means D is smaller than each of the resonant wavelengths λ1 and λ2, where c is the propagation speed of radiation in the surrounding medium.
Embodiments of the method may include any of the following features. In some embodiments, said resonators have Q1>100 and Q2>100, Q1>200 and Q2>200, Q1>500 and Q2>500, or even Q1>1000 and Q2>1000. In some such embodiments, idsqrtΓ1*Γ2) may be greater than 0.2, greater than 0.5, greater than 1, greater than 2, or even grater than 5. In some such embodiments D/L2 may be greater than 1, greater than 2, greater than 3, greater than 5.
In another aspect, an energy transfer device is disclosed which includes: a first resonator structure receiving energy from an external power supply, said first resonator structure having a first resonant frequency ω1, and a first Q-factor Q1, and characteristic size L1, and a second resonator structure being positioned distal from said first resonator structure, at closest distance D, said second resonator structure having a second resonant frequency ω2, and a second Q-factor Q2, and characteristic size L2.
The two said frequencies of and ω2 are close to within the narrower of the two resonance widths Γ1, and Γ2. The non-radiative energy transfer between said first resonator structure and said second resonator structure is mediated through coupling of their resonant-field evanescent tails, and the rate of energy transfer between said first resonator and said second resonator is denoted by κ. The non-radiative means D is smaller than each of the resonant wavelengths λ1 and λ2, where c is the propagation speed of radiation in the surrounding medium.
Embodiments of the device may include any of the following features. In some embodiments, said resonators have Q1>100 and Q2>100, Q1>200 and Q2>200, Q1>500 and Q2>500, or even Q1>1000 and Q2>1000. In some such embodiments, idsqrt(Γ1*Γ2) may be greater than 0.2, greater than 0.5, greater than 1, greater than 2, or even grater than 5. In some such embodiments D/L2 may be greater than 1, greater than 2, greater than 3, or even greater than 5.
In some embodiments, the resonant field in the device is electromagnetic.
In some embodiments, the first resonator structure includes a dielectric sphere, where the characteristic size L1 is the radius of the sphere.
In some embodiments, the first resonator structure includes a metallic sphere, where the characteristic size L1 is the radius of the sphere.
In some embodiments, the first resonator structure includes a metallodielectric sphere, where the characteristic size L1 is the radius of the sphere.
In some embodiments, the first resonator structure includes a plasmonic sphere, where the characteristic size L1 is the radius of the sphere.
In some embodiments, the first resonator structure includes a polaritonic sphere, where the characteristic size L1 is the radius of the sphere.
In some embodiments, the first resonator structure includes a capacitively-loaded conducting-wire loop, where the characteristic size L1 is the radius of the loop. In some embodiments, the second resonator structure includes a dielectric sphere, where the characteristic size L2 is the radius of the sphere.
In some embodiments, the second resonator structure includes a metallic sphere where the characteristic size L2 is the radius of the sphere.
In some embodiments, the second resonator structure includes a metallodielectric sphere where the characteristic size L2 is the radius of the sphere.
In some embodiments, the second resonator structure includes a plasmonic sphere where the characteristic size L2 is the radius of the sphere.
In some embodiments, the second resonator structure includes a polaritonic sphere where the characteristic size L2 is the radius of the sphere.
In some embodiments, the second resonator structure includes a capacitively-loaded conducting-wire loop where the characteristic size L2 is the radius of the loop.
In some embodiments, the resonant field in the device is acoustic.
It is to be understood that embodiments of the above described methods and devices may include any of the above listed features, alone or in combination.
In contrast to the currently existing schemes, the invention provides the feasibility of using long-lived oscillatory resonant electromagnetic modes, with localized slowly evanescent field patterns, for wireless non-radiative energy transfer. The basis of this technique is that two same-frequency resonant objects tend to couple, while interacting weakly with other off-resonant environmental objects. The purpose of the invention is to quantify this mechanism using specific examples, namely quantitatively address the following questions: up to which distances can such a scheme be efficient and how sensitive is it to external perturbations. Detailed theoretical and numerical analysis show that a mid-range (LTRANS≈few*LDEV) wireless energy-exchange can actually be achieved, while suffering only modest transfer and dissipation of energy into other off-resonant objects.
The omnidirectional but stationary (non-lossy) nature of the near field makes this mechanism suitable for mobile wireless receivers. It could therefore have a variety of possible applications including for example, placing a source connected to the wired electricity network on the ceiling of a factory room, while devices, such as robots, vehicles, computers, or similar, are roaming freely within the room. Other possible applications include electric-engine buses, RFIDs, and perhaps even nano-robots. Similarly, in some embodiments multiple sources can transfer energy to one or more device objects. For example, as explained at in the paragraph bridging pages 4-5 of U.S. Provisional Application No. 60/698,442 to which the present application claims benefit and which is incorporated by reference above, for certain applications having uneven power transfer to the device object as the distance between the device and the source changes, multiple sources can be strategically placed to alleviate the problem, and/or the peak amplitude of the source can be dynamically adjusted.
The range and rate of the inventive wireless energy-transfer scheme are the first subjects of examination, without considering yet energy drainage from the system for use into work. An appropriate analytical framework for modeling the exchange of energy between resonant objects is a weak-coupling approach called “coupled-mode theory”.
where ω1,2 are the individual eigen-frequencies, Γ1,2 are the resonance widths due to the objects' intrinsic (absorption, radiation etc.) losses, κ12,21 are the coupling coefficients, and κ11,22 model the shift in the complex frequency of each object due to the presence of the other.
The approach of Eq. 1 has been shown, on numerous occasions, to provide an excellent description of resonant phenomena for objects of similar complex eigen-frequencies (namely |ω1−ω2|<<|κ12,21| and Γ1≈Γ2), whose resonances are reasonably well defined (namely Γ1,2&Im{κ11,22}<<|κ12,21|) and in the weak coupling limit (namely κ12,21|<<ω1,2). Coincidentally, these requirements also enable optimal operation for energy transfer. Also, Eq. (1) show that the energy exchange can be nearly perfect at exact resonance (ω1=ω2 and Γ1=Γ2), and that the losses are minimal when the “coupling-time” is much shorter than all “loss-times”. Therefore, the invention requires resonant modes of high Q=ω/(2Γ) for low intrinsic-loss rates Γ1,2, and with evanescent tails significantly longer than the characteristic sizes L1 and L2 of the two objects for strong coupling rate |κ12,21| over large distances D, where D is the closest distance between the two objects. This is a regime of operation that has not been studied extensively, since one usually prefers short tails, to minimize interference with nearby devices.
Objects of nearly infinite extent, such as dielectric waveguides, can support guided modes whose evanescent tails are decaying exponentially in the direction away from the object, slowly if tuned close to cutoff, and can have nearly infinite Q. To implement the inventive energy-transfer scheme, such geometries might be suitable for certain applications, but usually finite objects, namely ones that are topologically surrounded everywhere by air, are more appropriate.
Unfortunately, objects of finite extent cannot support electromagnetic states that are exponentially decaying in all directions in air, since in free space: {right arrow over (k)}2=ω2/c2. Because of this, one can show that they cannot support states of infinite Q. However, very long-lived (so-called “high-Q”) states can be found, whose tails display the needed exponential-like decay away from the resonant object over long enough distances before they turn oscillatory (radiative). The limiting surface, where this change in the field behavior happens, is called the “radiation caustic”, and, for the wireless energy-transfer scheme to be based on the near field rather than the far/radiation field, the distance between the coupled objects must be such that one lies within the radiation caustic of the other.
The invention is very general and any type of resonant structure satisfying the above requirements can be used for its implementation. As examples and for definiteness, one can choose to work with two well-known, but quite different electromagnetic resonant systems: dielectric disks and capacitively-loaded conducting-wire loops. Even without optimization, and despite their simplicity, both will be shown to exhibit fairly good performance. Their difference lies mostly in the frequency range of applicability due to practical considerations, for example, in the optical regime dielectrics prevail, since conductive materials are highly lossy.
Consider a 2D dielectric disk cavity of radius r and permittivity ε surrounded by air that supports high-Q whispering-gallery modes, as shown in
The radial modal decay length, which determines the coupling strength κ≡|κ21|=|κ12|, is on the order of the wavelength, therefore, for near-field coupling to take place between cavities whose distance is much larger than their size, one needs subwavelength-sized resonant objects (r<<λ). High-radiation-Q and long-tailed subwavelength resonances can be achieved, when the dielectric permittivity ε is as large as practically possible and the azimuthal field variations (of principal number m) are slow (namely m is small).
One such TE-polarized dielectric-cavity mode, which has the favorable characteristics Qrad=1992 and λ/r=20 using ε=147.7 and m=2, is shown in
With regards to material absorption, typical loss tangents in the microwave (e.g. those listed for the materials above) suggest that Qabs˜ε/Im{ε}˜10000. Combining the effects of radiation and absorption, the above analysis implies that for a properly designed resonant device-object d a value of Qd˜2000 should be achievable. Note though, that the resonant source s will in practice often be immobile, and the restrictions on its allowed geometry and size will typically be much less stringent than the restrictions on the design of the device; therefore, it is reasonable to assume that the radiative losses can be designed to be negligible allowing for Qs˜10000, limited only by absorption.
To calculate now the achievable rate of energy transfer, one can place two of the cavities 20, 22 at distance D between their centers, as shown in
Consider a loop 10 or 12 of N coils of radius r of conducting wire with circular cross-section of radius a surrounded by air, as shown in
For non-radiative coupling one should use the near-field region, whose extent is set roughly by the wavelength λ, therefore the preferable operating regime is that where the loop is small (r<<λ). In this limit, the resistances associated with the two loss channels are respectively Rohm=√{square root over (μOρω/2·)}Nr/a and Rrad=π/6·ηON2 (ωr/c)4, where ρ is the resistivity of the wire material and ηO≈120πΩ is the impedance of free space. The quality factor of such a resonance is then Q=ωL/(Rohm+Rrad) and is highest for some frequency determined by the system parameters: at lower frequencies it is dominated by ohmic loss and at higher frequencies by radiation.
To get a rough estimate in the microwave, one can use one coil (N=1) of copper (ρ=1.69·10−8 Ωm) wire and then for r=1 cm and a=1 mm, appropriate for example for a cell phone, the quality factor peaks to Q=1225 at f=380 MHz, for r=30 cm and a=2 mm for a laptop or a household robot Q=1103 at f=17 MHz, while for r=1 m and a=4 mm (that could be a source loop on a room ceiling) Q=1315 at f=5 MHz. So in general, expected quality factors are Q≈1000-1500 at a λ/r≈50-80, namely suitable for near-field coupling.
The rate for energy transfer between two loops 10 and 12 at distance D between their centers, as shown in
It is important to appreciate the difference between this inductive scheme and the already used close-range inductive schemes for energy transfer in that those schemes are non-resonant. Using coupled-mode theory it is easy to show that, keeping the geometry and the energy stored at the source fixed, the presently proposed resonant-coupling inductive mechanism allows for Q approximately 1000 times more power delivered for work at the device than the traditional non-resonant mechanism, and this is why mid-range energy transfer is now possible. Capacitively-loaded conductive loops are actually being widely used as resonant antennas (for example in cell phones), but those operate in the far-field regime with r/λ˜1, and the radiation Q's are intentionally designed to be small to make the antenna efficient, so they are not appropriate for energy transfer.
Clearly, the success of the inventive resonance-based wireless energy-transfer scheme depends strongly on the robustness of the objects' resonances. Therefore, their sensitivity to the near presence of random non-resonant extraneous objects is another aspect of the proposed scheme that requires analysis. The interaction of an extraneous object with a resonant object can be obtained by a modification of the coupled-mode-theory model in Eq. (1), since the extraneous object either does not have a well-defined resonance or is far-off-resonance, the energy exchange between the resonant and extraneous objects is minimal, so the term κ12 in Eq. (1) can be dropped. The appropriate analytical model for the field amplitude in the resonant object a1(t) becomes:
Namely, the effect of the extraneous object is just a perturbation on the resonance of the resonant object and it is twofold: First, it shifts its resonant frequency through the real part of κ11 thus detuning it from other resonant objects. As shown in
In the first example of resonant objects that have been considered, the class of dielectric disks, small, low-index, low-material-loss or far-away stray objects will induce small scattering and absorption. To examine realistic cases that are more dangerous for reduction in Q, one can therefore place the “test” dielectric disk cavity 40 close to: a) another off-resonance object 42, such as a human being, of large Re{ε}=49 and Im{ε}=16 and of same size but different shape, as shown in
Analytically, for objects that interact with a small perturbation the reduced value of radiation-Q due to scattering could be estimated using the polarization ∫d3r|PX1(r)|2 ∝∫d3r|E1(r)·Re{εX(r)}|2 induced by the resonant cavity 1 inside the extraneous object X=42 or roughened surface X=46. Since in the examined cases either the refractive index or the size of the extraneous objects is large, these first-order perturbation-theory results would not be accurate enough, thus one can only rely on numerical FDTD simulations. The absorption-Q inside these objects can be estimated through Im{κ11}=ω1/2·∫d3r|E1(r)2Im{εX(r)}/∫d3r|E1(r)|2ε(r).
Using these methods, for distances D/r=10, 7, 5, 3 between the cavity and extraneous-object centers one can find that Qrad=1992 is respectively reduced to Qrad=1988, 1258, 702, 226, and that the absorption rate inside the object is Qabs=312530, 86980, 21864, 1662, namely the resonance of the cavity is not detrimentally disturbed from high-index and/or high-loss extraneous objects, unless the (possibly mobile) object comes very close to the cavity. For distances D/r=10, 7, 5, 3, 0 of the cavity to the roughened surface we find respectively Qrad=2101, 2257, 1760, 1110, 572, and Qabs>4000, namely the influence on the initial resonant mode is acceptably low, even in the extreme case when the cavity is embedded on the surface. Note that a close proximity of metallic objects could also significantly scatter the resonant field, but one can assume for simplicity that such objects are not present.
Imagine now a combined system where a resonant source-object s is used to wirelessly transfer energy to a resonant device-object d but there is an off-resonance extraneous-object e present. One can see that the strength of all extrinsic loss mechanisms from e is determined by |Es(re)|2, by the square of the small amplitude of the tails of the resonant source, evaluated at the position re of the extraneous object. In contrast, the coefficient of resonant coupling of energy from the source to the device is determined by the same-order tail amplitude |Es(rd)|, evaluated at the position rd of the device, but this time it is not squared! Therefore, for equal distances of the source to the device and to the extraneous object, the coupling time for energy exchange with the device is much shorter than the time needed for the losses inside the extraneous object to accumulate, especially if the amplitude of the resonant field has an exponential-like decay away from the source. One could actually optimize the performance by designing the system so that the desired coupling is achieved with smaller tails at the source and longer at the device, so that interference to the source from the other objects is minimal.
The above concepts can be verified in the case of dielectric disk cavities by a simulation that combines
In the second example of resonant objects being considered, the conducting-wire loops, the influence of extraneous objects on the resonances is nearly absent. The reason for this is that, in the quasi-static regime of operation (r<<λ) that is being considered, the near field in the air region surrounding the loop is predominantly magnetic, since the electric field is localized inside the capacitor. Therefore, extraneous objects that could interact with this field and act as a perturbation to the resonance are those having significant magnetic properties (magnetic permeability Re{μ}>1 or magnetic loss Im{μ}>0). Since almost all common materials are non-magnetic, they respond to magnetic fields in the same way as free space, and thus will not disturb the resonance of a conducting-wire loop. The only perturbation that is expected to affect these resonances is a close proximity of large metallic structures.
An extremely important implication of the above fact relates to safety considerations for human beings. Humans are also non-magnetic and can sustain strong magnetic fields without undergoing any risk. This is clearly an advantage of this class of resonant systems for many real-world applications. On the other hand, dielectric systems of high (effective) index have the advantages that their efficiencies seem to be higher, judging from the larger achieved values of κ/T, and that they are also applicable to much smaller length-scales, as mentioned before.
Consider now again the combined system of resonant source s and device d in the presence of a human h and a wall, and now let us study the efficiency of this resonance-based energy-transfer scheme, when energy is being drained from the device for use into operational work. One can use the parameters found before: for dielectric disks, absorption-dominated loss at the source Qs˜104, radiation-dominated loss at the device Qd˜103 (which includes scattering from the human and the wall), absorption of the source- and device-energy at the human Qs-h, Qd-h˜104-105 depending on his/her not-very-close distance from the objects, and negligible absorption loss in the wall; for conducting-wire loops, Qs˜Qd˜103, and perturbations from the human and the wall are negligible. With corresponding loss-rates Γ=ω/2Q, distance-dependent coupling κ, and the rate at which working power is extracted Γw, the coupled-mode-theory equation for the device field-amplitude is
Different temporal schemes can be used to extract power from the device and their efficiencies exhibit different dependence on the combined system parameters. Here, one can assume steady state, such that the field amplitude inside the source is maintained constant, namely as(t)=Ase−iωt, so then the field amplitude inside the device is ad(t)=Ade−iωt with Ad=iκ(Γd+Γd-h+Γw)As. Therefore, the power lost at the source is Ps=2Γs|As|2, at the device it is Pd=2Γd|Ad|2, the power absorbed at the human is Ph=2Γs-h|As|2+2Γd-h|Ad|2, and the useful extracted power is Pw=2Γw|Ad|2. From energy conservation, the total power entering the system is Ptotal=Ps+Pd+Ph+Pw. Denote the total loss-rates Γstot=ΓsΓs-h and Γdtot=Γd+Γd-h. Depending on the targeted application, the work-drainage rate should be chosen either Γw=Γdtot to minimize the required energy stored in the resonant objects or Γw=Γdtot+√{square root over (1+κ2/ΓstotΓdtot)}>Γdtot such that the ratio of useful-to-lost powers namely the efficiency ηw=Pw/Ptotal, is maximized for some value of κ. The efficiencies η for the two different choices are shown in
Even better performance should be achievable through optimization of the resonant object designs. Also, by exploiting the earlier mentioned interference effects between the radiation fields of the coupled objects, such as continuous-wave operation at the frequency of the normal mode that has the larger radiation-Q, one could further improve the overall system functionality. Thus the inventive wireless energy-transfer scheme is promising for many modern applications. Although all considerations have been for a static geometry, all the results can be applied directly for the dynamic geometries of mobile objects, since the energy-transfer time κ−1˜1 μs, which is much shorter than any timescale associated with motions of macroscopic objects.
The invention provides a resonance-based scheme for mid-range wireless non-radiative energy transfer. Analyses of very simple implementation geometries provide encouraging performance characteristics for the potential applicability of the proposed mechanism. For example, in the macroscopic world, this scheme could be used to deliver power to robots and/or computers in a factory room, or electric buses on a highway (source-cavity would in this case be a “pipe” running above the highway). In the microscopic world, where much smaller wavelengths would be used and smaller powers are needed, one could use it to implement optical inter-connects for CMOS electronics or else to transfer energy to autonomous nano-objects, without worrying much about the relative alignment between the sources and the devices; energy-transfer distance could be even longer compared to the objects' size, since Im{ε(ω)} of dielectric materials can be much lower at the required optical frequencies than it is at microwave frequencies.
As a venue of future scientific research, different material systems should be investigated for enhanced performance or different range of applicability. For example, it might be possible to significantly improve performance by exploring plasmonic systems. These systems can often have spatial variations of fields on their surface that are much shorter than the free-space wavelength, and it is precisely this feature that enables the required decoupling of the scales: the resonant object can be significantly smaller than the exponential-like tails of its field. Furthermore, one should also investigate using acoustic resonances for applications in which source and device are connected via a common condensed-matter object.
Although the present invention has been shown and described with respect to several preferred embodiments thereof, various changes, omissions and additions to the form and detail thereof, may be made therein, without departing from the spirit and scope of the invention.
This application is a continuation and claims the benefit of priority under 35 USC § 120 to U.S. application Ser. No. 16/851,598, filed on Apr. 17, 2020, which is a continuation of U.S. application Ser. No. 16/184,354, filed on Nov. 8, 2018, now U.S. Pat. No. 10,666,091, which is a continuation of U.S. application Ser. No. 15/793,198, filed on Oct. 25, 2017, now U.S. Pat. No. 10,141,790, which is a continuation of U.S. application Ser. No. 15/083,726, filed on Mar. 29, 2016, now U.S. Pat. No. 9,831,722, which is a continuation of U.S. application Ser. No. 14/629,709, filed on Feb. 24, 2015, now U.S. Pat. No. 9,450,421, which is a continuation of U.S. application Ser. No. 14/302,662, filed on Jun. 12, 2014, now U.S. Pat. No. 9,065,286, which is a continuation of U.S. application Ser. No. 12/639,963, filed on Dec. 16, 2009, now U.S. Pat. No. 8,760,007, which is a continuation of U.S. application Ser. No. 12/553,957, filed on Sep. 3, 2009, which is a continuation of U.S. application Ser. No. 11/481,077 filed on Jul. 5, 2006, now U.S. Pat. No. 7,741,734, which claims priority under 35 USC § 119(e) to U.S. provisional application Ser. No. 60/698,442 filed on Jul. 12, 2005. The contents of the prior applications mentioned above are incorporated herein by reference in their entirety.
This invention was made with government support awarded by the National Science Foundation under Grant No. DMR-0213282. The government has certain rights in this invention.
Number | Date | Country | |
---|---|---|---|
60698442 | Jul 2005 | US |
Number | Date | Country | |
---|---|---|---|
Parent | 16851598 | Apr 2020 | US |
Child | 17462126 | US | |
Parent | 16184354 | Nov 2018 | US |
Child | 16851598 | US | |
Parent | 15793198 | Oct 2017 | US |
Child | 16184354 | US | |
Parent | 15083726 | Mar 2016 | US |
Child | 15793198 | US | |
Parent | 14629709 | Feb 2015 | US |
Child | 15083726 | US | |
Parent | 14302662 | Jun 2014 | US |
Child | 14629709 | US | |
Parent | 12639963 | Dec 2009 | US |
Child | 14302662 | US | |
Parent | 12553957 | Sep 2009 | US |
Child | 12639963 | US | |
Parent | 11481077 | Jul 2006 | US |
Child | 12553957 | US |