Tuesday, February 22, 2022

Introduction to Polar Diagrams and Optimum VMC

A polar diagram is a graphic presentation of a sailboat's sailing performance in various wind conditions. The same information is also often presented in the form of a table or spreadsheet, which is how the data are originally collected or computed, which in turn is then plotted on a polar diagram as shown below. Programs that compute this data are called velocity prediction programs (VPP).



This is a polar diagram from a specific Farr40 as it appears graphically in LuckGrib above, which includes multiple polar diagrams in its routing function, plus the option to import custom versions. 

Below are the same data in table format. Here we learn that this Farr40 sailing in a true wind speed (TWS) of 20 kts at a true wind angle (TWA) of 135º should in principle have a boat speed of 11.85 kts. We can read this from the diagram, but get more precise readings from the table. As discussed below, these are theoretical or target speeds based on yacht design. It is up to the navigator to determine if these are realistic values or do they need to be modified to match actual performance.


A convenient combination of this information can be presented in a navigation program that does a digital interpolation of the table and diagram simultaneously, as shown below in this example from qtVlm, another large source of polar diagrams.


In this qtVlm display, we set the true wind speed (to any value, does not need to be in the original table data), and then putting the cursor on the diagram creates the blue line which can be moved to the desired TWA as the corresponding boat speed and other data are shown on the right.  This app also has a table view, but it is mainly intended for editing the polar, since it can double interpolate as shown here (i.e., any TWS or TWA) to get specific values.

We see here also extra output that is common to polar diagram presentation, namely the prediction of the best performance directly upwind and directly downwind, as illustrate below.


This display shows what your best TWA is when your destination is directly upwind or directly downwind. For this Farr40, in 10 kts of TWS, the best upwind TWA is 41º at which time your boat speed should read 6.9 kts, and your actual progress up wind, called velocity made good (VMG) would be 5.24 kts. You would tack back and forth till you reached the layline to the mark.

The same is almost true when sailing downwind, but there is a difference. Sailing upwind, by definition you would always be sailing the optimum upwind TWA for that TWS.  Unless there are environmental or tactical reasons to do otherwise, that is the way you best progress upwind—and any target that is not on the layline is by definition upwind.  In other words, you would always choose to sail upwind at the TWA that is right where the red curve meets the green curve on this polar diagram presentation, because that is best speed upwind.  In the real world, one tack might be faster than the other, which affects this decision, but we are not considering such perturbations for now.

Sailing downwind, if the target is indeed dead downwind, then you would again sail the TWA recommended where the red meets the green, jibing back and forth across that line. For example, if the true wind direction (TWD) were 070, blowing toward 250, and the rhumb line course to your target was 250, then that is what you would do for best performance—again assuming speeds were equal on both jibes and you did not have other reasons to favor one jibe over the other.

With this Farr40 in 10 kts of true wind this means we would want to sail at TWA = 150, which we achieve at a heading of 220 (30º to the left of the direction the wind blows toward) where we expect a knotmeter speed = 7.6, and a projected speed downwind of VMG = 6.6, which we assume will be the same on both jibes. 

In practice, however,  the target is not likely to be dead downwind, so the job is not to maximize the speed downwind (VMG), but to maximize the speed in the direction of the desired course, called velocity made course (VMC). It is the projection of the speed you are sailing onto the direction of the desired course, as shown below.


Speed over ground (SOG) and course over ground (COG) are direct outputs from the GPS. The course to next mark (CNM) is computed from the boat's present position and the Lat-Lon of the destination.

We can go over a couple numerical examples of finding VMC from the polar diagram, but in practice this VMC is usually found by experimentation under sail, because all instrument systems can display the VMC being accomplished, as determined by the SOG and COG from the GPS and the known heading to your desired course, CNM. It is a simple, accurate computation that does not depend on anything difficult to measure. 

Following up on the example above, let's say the desired course was not dead downwind at 250T, but instead it was at 220T, which is 30º to the left of dead downwind.  Our first reaction looking at the polar data shown above might be "Lucky us; our destination is right on the heading that gives  maximum speed downwind."  But things are not quite that simple; on closer look, that is not the optimum heading for our destination.  We have several ways to figure that out directly from the polar plot, independent of what nav program we use. 





Recalling that the polar is aways oriented in the direction of the wind, we can think of looking across it in the direction we want to go, as shown above. The best TWA downwind was 150, meaning with wind at 070 we steer a heading of 220,  which yields  a best VMC = 6.6 kts, but if our destination is actually at the course of 220 (not 250 which is dead downwind) the best TWA = 135, which we achieve at a heading of 205, yielding a best VMC of just over 8 kts.  Again, we do not often need to solve this ahead of time because we can find it experimentally as shown below, but it can help to understand the concept. We can get more precise values from a spreadsheet as shown below, where the "polar speeds" shown were interpolated from the qtVlm function mentioned above.



Most nav apps offer more direct ways to determine the optimum VMC. qtVlm, for example, lets us draw the polar right on the screen that matches the existing wind forecast, or we can mock up a case with its Force-wind option (Grib config/Corrections/Force wind) set for this example to 10 kts from 070. Then we can draw various lines to look at options, as shown below. 


This presentation of the on-screen "speed polar" shows both sides, where as typical presentation is just one side. Some nav apps have a way to think through this and display the optimum VMC without special setup (see qtVlm's pink arrow below), but in all cases we still have to go to the approximate heading we think is right and then try 10 or 15º either side to see what maximizes the VMC. Note that the max values cover a fairly broad range of TWA and the sea state will always have an effect on this.  Below is an illustration of this using the qtVlm simulator.


Beside watching a plot of VMC vs time, qtVlm shows two arrows on the microboard at all times, one showing heading to the course and the other (pink one) is optimum VMC taking into account wind, waves, and currents.  Below is a short video (5:43) showing where the above data came from.



When sailing an ocean route downwind, one of the ongoing concerns is "Are we on the right jibe?"  Keeping track of the optimum VMC is one step to answering that, but as the wind changes and the boat speed changes, we also have to think through what this value would be on the other jibe, and these basics might help with that if you do not have digital tools that display the opposite jibe VMC at all times.


Polar Formats

Polar diagram files require special formats. OpenCPN has an extended discussion of the formats. Luckgrib accepts almost any style at all; qtVlm requires the CSV format to use semicolons, not comma separators. Expedition has their own POL format

The text files for most polars require a CRLF (carriage return, line feed). Exporting an PC excel file with MS DOS txt format adds the proper CRLF at the end of each row, but Mac Excel seems to put just a CR.  If you get an error  message from a program such as "polar does not have sufficient data" it could be do to missing LF. You can check this in a PC by loading the file into Notepad++ and then View/Show Symbols/All Characters, as shown below.



In a Mac you can use terminal command "cat -e" to check this.

macdavid@iMac Downloads % cat -e etap39.csv 
twa\tws,6,8,10,12,14,16,20^M$
0,0,0,0,0,0,0,0^M$
46.5,4.61,0,0,0,0,0,0^M$
45.8,0,5.42,0,0,0,0,0^M$

The ^M$ is the correct CRLF symbol:  ^M is the CR and $ is the LF


Sources of Polars

The programs Expedition, LuckGrib, qtVlm, and others each have an extensive set of polars included. The set within Expedition has been popular amongst racing sailors for many years, but many of these are also found at other sites.

ORC Sailboat Data has a lot and lets you see them graphically, and compare for same boat.

Polar server from qtVlm has a lot of options.

iPolar is a $10 iOS app that creates a very basic starting polar based on boat specs. You can get the needed boat specs for most boats from sailboatdata.com

Get samples for specific yachts using this link to find a country ID.


Nuances to keep in mind

The goal is always to build your polar with actual measurements, but we often have to start out with some stock polar, and the closest one we find might not have data at the low and high ends of the TWA and it might not have data for higher or lower TWS. Below is an example we found that illustrates this. 


We have fairly good TWS coverage, but not good TWA going only from 52 to 150.  This polar is shown below loaded into two different navigation programs that do optimum sailboat routing.




The actual data in the polar are those between the yellow and orange lines, but each program has filled in the missing data using some VPP approximation, heuristic, or simple rule. 

The one on the right has a simple rule. For low wind angles, it is a linear interpolation from last data point in the polar to 0 kts at head to wind (TWA=0). It does not look linear, but that is all due to the polar plot.  On the downwind end, this program assumes the boat speed remains the same at all TWA aft of the last known value, TWA= 150 in this case.  As such, during a routing computation, it will favor dead downwind for optimum VMG, which is not likely right at lower TWS. Note too that this accounts for the apparent kink in the curves. That is not real; it just reflects the intersection of a fixed circle and a curve.

The computed polar on the left uses a heuristic approach, making stronger fall offs in boat speed for lower TWS. At medium and lower TWS either one of these two polar approaches might be closest to the truth, depending on the boat and sails. 

The point being made, however, is both polars include some region of not real data, so after a route is run it would pay to check how much time is spent in no man's land of the polar.  One way to protect yourself from getting routed off to the real no man's land, is to limit the TWA angles in the routing to say, for this boat, no closer than 50 upwind and no farther than say 165 downwind.

But the main job in a case like this is be sure to start collecting real logbook data that directly lets you improve your polar data and especially fill in the missing angle ranges. For a given TWS it would not take long to learn both ends and then just add them to the polar.  We want in the logbook columns for TWA, TWS, boat speed (STW), plus sails set and maybe notes on sea state. Your wind instruments might give you this directly, and if not then you need heading (CTW), STW, AWA, AWS, and from these you compute TWS and TWA. In either case all instruments must be well calibrated: AWS, AWD, and heading sensor, and if heeled you need to know the heel angle. See Modern Marine Weather.

One easy way to fill in the gaps is to use an app like iPolar to get a rough prediction of the curves then scale them to what you think are the right data from your starting polar, and then you will be using some combination of VPP and heuristic to get starting values for the missing data.

There are also sophisticated tools such as Polar Manager that let you then load both the stock polar and the iPolar approximation and then manually average them to some common factor...  which then awaits your confirmation with real data.

I hope to have data shortly to make a demo of creating a working polar for a cruising sailboat, starting from a rough base computed in iPolar.




Saturday, February 19, 2022

ASCAT News — The Good and the Bad.

Scatterometer wind measurements are the truth meter for weather work at sea. There are several sources, but the ASCAT data from EUMETSAT are the gold standard. We can see this data graphically at the STAR website—found easily with a Google search of "ASCAT." We have extensive coverage of this data and how to use it in Modern Marine Weather (MMWX), but the book links are out of date (being updated here!)—in part because the satellites have changed. A is now gone and we use Metop B  and C.

This article has been updated from its origin on May 2021, with new videos at the end.

So the first news is we have the updated links below, there are new data, and it is even better than before. The main bad news is the GRIB format of this data, which offered a super convenient and precise way to test model forecasts is at present no longer available from the Ocean Prediction Center (OPC). It is still listed at the OPC website, which might imply the loss is not permanent. The files look normal and are indeed updated as indicated, but they do not contain any data!

Granted, this wonderful source of digital data could only be conveniently viewed in two popular commercial programs, namely Expedition and LuckGrib—not to mention the free app Panoply, since that one takes extra work. Worse than that, even amongst the many users of these two leading programs, there were not many navigators actually taking advantage of this most important data, which might account for the lingering loss... ie there are not many people complaining. 

(For now, there is no GRIB source for ASCAT, but for those who want it enough to spend some extra work you can download the data in netcdf and then view it in Panoply, but this is not a practical, navigator's solution.)

*** Note added in July 2022: LuckGrib now offers ASCAT once again. They are doing the conversion from netcdf to grib on their own. To our knowledge this is the only grib source of ASCAT data.

Part of the problem in ASCAT popularity is the data are not plug and play like a GFS forecast is, for example. Once we did get the ASCAT grib, life was good, it showed the actual winds on the ocean, not a forecast, but real data as if we are looking at anemometers on a field buoys located every 25 km across the ocean. The problem is the data over some specific point, say 300 nmi in radius, might only be updated 2 or 3 times a day, and those times take some effort to predict. Refer to MMWX for these details.

The good news is, and the motivation of this article is, we have a new, free way to get georeferenced ASCAT data into a navigation program when underway using the program qtVlm. We already had one way described in MMWX, namely look up the file name for the data you want from the STAR site or from the table in MMWX, and then request that image from Saildocs. Once that is in hand, and you check the valid time of the satellite pass (MMWX), you can either just look at it to see the winds and compare to the forecasts, or in Expedition, OpenCPN, or qtVlm  you could manually georeference it and then overlay the model forecast at the corresponding tine right on top of it. In these programs you can store the georeference coordinates and then just periodically update to see if there is new data. The image size you are requesting is 20 to 60m kb. A sample is below.


ASCAT-A descending over Bermuda at 12:33z on May 2, 2021. It is May 2, and not May 3, because this image is from the 22-hr dataset that was most recently updated at 02:07 May 2. Had the small purple time at the bottom been between 00:00z and 02:07z then this would be May 2. (From the original article; A is no longer available.)


These files are always 10º Lat x 15º Lon so they are easy to georeference, plus the grid lines are clearly in view.

Below are samples georeferenced in Expedition and qtVlm with the GFS wind forecast at the corresponding time overlaid on it.


And then zoomed in below, we see the GFS is pretty good since it agrees with the scatterometer wind measurements. This is Metop A, descending. 


The same thing viewed in qtVlm looks like this:





The best comparison of overlaid data like this depends on our choices of transparency and wind barb  colors. This has not been optimized in either of these two programs.

The main new good news from qtVlm is their latest version allows us to enter a live internet link for the image, as opposed to a link to a static file on our hard drive, so that the data are updated automatically when we reload that image. The trick is to create an auto-updating KML file for image and then use that as the target file. This then not only updates automatically, it also brings with it all of the georeferencing information. With your computer linked to your satellite phone, this update should take place just as you do to obtain the latest model forecast.

We have videos online on how to do this.  We make the needed files in Google Earth, which is then another way to look at live, auto updating ASCAT images. The new step is just adding this functionality to a navigation program independent of Google Earth.

We made such KML files  some years ago for the Bermuda Race,  Transpac/Pacific Cup, and for Sydney-Hobart. Our next step will be to update the links and make a convenient place for mariners to download them.  Now that we have lost the GRIB format, these methods become of interest again.

Here is the generic image link using the Bermuda file as the base. Then the colored letters have to be changed for other locations, passes, and satellites.

https://manati.star.nesdis.noaa.gov/ascat_images/cur_25km_META/zooms/WMBds86.png

Change A to B or C for the 3 ASCAT satellites, and change the blue file name to match the location you care about, and then change the d (descending) to a for the ascending pass. 

Thus for each location you care about, there will be 6 files that might be of interest, keeping in mind that A and B are fairly close in time so the earth has not rotated much between them. In rare cases, you might need B having seen A, but not often.  

You can get these files (once you know the name of the ones you want) from Saildocs. For example, send an otherwise blank email to query@saildocs with this as the message:

send https://manati.star.nesdis.noaa.gov/ascat_images/cur_25km_META/zooms/WMBds86.png
send https://manati.star.nesdis.noaa.gov/ascat_images/cur_25km_META/zooms/WMBas86.png

This should get you those two files by return email. Sample file names below.


Click the images above for a better view.

Below are a couple How-to videos:



How to create KML files of ASCAT wind data.


Viewing ASCAT wind data in qtVlm, and one with more timely data.
Links to the Bermuda Race KLM files are in the video description.